Search results

This paper investigates the multi‐phase behaviour of droplets injected into a nozzle at two separate wall locations. The physical features of the droplets (rate of mass, density and radius) at each injector location are identical. This system can be described by a two‐phase Eulerian—Eulerian approach that yields classical systems of equations: three for the gaseous phase and three for the dispersed droplet phase. An underlying assumption in the two phase model is that no interaction occurs between droplets. The numerical solution of the model (using the MacCormack scheme) indicates however that the opposite jets do interact to form one jet. This inconsistency is overcome in the current paper by associating the droplets from a given injection location with a separate phase and subsequently solving equations describing a multiphase system (here, three‐phase system). Comparison of numerical predications between the two‐phase and the multiphase model shows significantly different results. In particular the multiphase model shows no jet interaction.

The purpose of this paper is to the study the multiphase bubbles flow motion in a vertical channel with an electroconducting liquid without and under the influence of a magnetic field.

For numerical calculations, the lattice Boltzmann method (LBM) is used, which is based on the kinetic theory for solving fluid mechanics and other physical problems. The phase-field lattice Boltzmann model is developed to simulate the behaviour of multiphase bubble–bubble interaction while rising in the fluid with high density ratios.

The behaviour of the rising bubble flow in a rectangular column of two phases is investigated with the two-dimensional LBM.

The multiphase flow in electroconducting liquids with high ratio of density is studied using the LBM.

The purpose of this study is to propose a novel simulation framework and show that it captures the main effects of the deposition process, such as droplet shape, volume and speed.

In the framework, the time-dependent flow and the fluid-structure interaction between the suspension, the moving piston and the deflection of the jetting head is simulated. The system is modelled as a two-phase system with the surrounding air being one phase and the dense suspension the other. The non-Newtonian suspension is modelled as a mixed single phase with properties determined from material testing. The simulations were performed with two coupled in-house solvers developed at Fraunhofer-Chalmers Centre; IBOFlow, a multiphase flow solver; and LaStFEM, a large strain FEM solver. Experimental deposition was performed with a commercial jet printer and quantitative measurements were made with optical profilometry.

Jetting behaviour was shown to be affected by not only piston motion, fluid rheology and head deformation but also the viscous energy loss in the jetting head nozzle. The simulation results were compared to experimental data obtained from an industrial jetting head and found to match characteristic lengths, speed and volume within ca 10%.

The simulations are based on a rheological description using the Carreau model that does not include a time-dependent relaxation of the fluid. This modelling approach limits the descriptive nature of the deposit after impact on the substrate. The simulation also adopts a continuum approach to the suspension, which will not accurately model the break-off of the droplet filament under the characteristic diameter of the particles in the suspension.

The ability to accurately simulate the deposition of functional materials will enable the efficient development of novel product designs with a minimum of used resources and minimised product development duration.

The ability to accurately simulate the deposition of functional materials will enable the efficient development of novel product designs with a minimum of used resources and therefore an improvement from a sustainability perspective. The ability to plan deposition strategies virtually will also enable a decrease in consumables at manufacturers which will in turn decrease their carbon foot print.

While basic fluid dynamic simulations have been performed to simulate flow through nozzles, the ability to include both fluid-structure interaction and multiphase capability together with a more accurate rheological description of the suspension and with a substrate for surface mount applications has not been published to the knowledge of the authors.

This work aims to investigate the effects of neglecting, modelling or partly resolving turbulent fluctuations of velocity, temperature and concentrations on the predicted turbulence-chemistry interaction in urea-selective non-catalytic reduction (SNCR) systems.

Numerical predictions of the NO conversion efficiency in an industrial urea-SNCR system are compared to experimental data. Reactor models of varying complexity are assessed, ranging from one-dimensional ideal reactor models to state-of-the-art computational fluid dynamics simulations based on the detached-eddy simulation (DES) approach. The models use the same reaction mechanism but differ in the degree to which they resolve the turbulent fluctuations of the gas phase. A methodology for handling of unknown experimental data with regard to providing adequate boundary conditions is also proposed.

One-dimensional reactor models may be useful for a first quick assessment of urea-SNCR system performance. It is critical to account for heat losses, if present, due to the significant sensitivity of the overall process to temperature. The most comprehensive DES setup evaluated is associated with approximately two orders of magnitude higher computational cost than the conventional Reynolds-averaged Navier–Stokes-based simulations. For studies that require a large number of simulations (e.g. optimizations or handling of incomplete experimental data), the less costly approaches may be favored with a tolerable loss of accuracy.

Novel numerical and experimental results are presented to elucidate the role of turbulent fluctuations on the performance of a complex, turbulent, reacting multiphase flow.

A consistent implementation of the general computational framework of unified second-order time accurate integrators via the well-known GSSSS framework in conjunction with the traditional Finite Difference Method is presented to improve the numerical simulations of reactive two-phase flows.

In the present paper, the phase interaction evaluation in the present implementation of the reactive two-phase flows has been derived and implemented to preserve the consistency of the correct time level evaluation during the time integration process for solving the two phase flow dynamics with reactions.

Numerical examples, including the classical Sod shock tube problem and a reactive two-phase flow problem, are exploited to validate the proposed time integration framework and families of algorithms consistently to second order in time accuracy; this is in contrast to the traditional practices which only seem to obtain first-order time accuracy because of the inconsistent time level implementation with respect to the interaction of two phases. The comparisons with the traditional implementation and the advantages of the proposed implementation are given in terms of the improved numerical accuracy in time. The proposed approaches provide a correct numerical simulation implementation to the reactive two-phase flows and can obtain better numerical stability and computational features.

The new algorithmic framework and the consistent time level evaluation extended with the GS4 family encompasses a multitude of past and new schemes and offers a general purpose and unified implementation for fluid dynamics.

This paper aims to propose a three-phase interpenetrating continua model for the numerical simulation of water waves and porous structure interaction.

In contrast with one-fluid formulation or multi-component methods, each phase has its own characteristics, density, velocity, etc., and each point is occupied by all phases. First, the porous structure is modelled as a phase of continua with a penalty force adding on the momentum equation, so the conservation of mass is guaranteed without source terms. Second, the adaptive unstructured mesh modelling with P1DG-P1 elements is used here to decrease the total number of degree of freedom maintaining the same order of accuracy.

Several benchmark problems are used to validate the model, which includes the Darcy flow, classical collapse of water column and water column with a porous structure. The interpenetrating continua model is a suitable approach for water wave and porous structure interaction problem.

The interpenetrating continua model is first applied for the water wave and porous structure interaction problem. First, the structure is modelled as phase of non-viscous fluid with penalty force, so the break of the porous structure, porosity changes can be easily embedded for further complex studies. Second, the mass conservation of fluids is automatically satisfied without special treatment. Finally, adaptive anisotropic mesh in space is employed to reduce the computational cost.

A combination of highly conductive porous media and nanofluids is an efficient way for improving thermal performance of relevant applications. For precisely predicting the flow and thermal transport of nanofluids in porous media, the purpose of this paper is to explore the inter-phase coupling numerical methods.

Based on the lattice Boltzmann (LB) method, this study combines the convective flow, non-equilibrium thermal transport and phase interactions of nanofluids in porous matrix and proposes a new multi-phase LB model. The micro-scale momentum and heat interactions are especially analyzed for nanoparticles, base fluid and solid matrix. A set of three-phase LB equations for the flow/thermal coupling of base fluid, nanoparticles and solid matrix is established.

Distributions of nanoparticles, velocities for nanoparticles and the base fluid, temperatures for three phases and interaction forces are analyzed in detail. Influences of parameters on the nanofluid convection in the porous matrix are examined. Thermal resistance of nanofluid convective transport in porous structures are comprehensively discussed with the models of multi-phases. Results show that the Rayleigh number and the Darcy number have significant influences on the convective characteristics. The result with the three-phase model is mildly larger than that with the local thermal non-equilibrium model.

This paper first creates the multi-phase theoretical model for the complex coupling process of nanofluids in porous structures, which is useful for researchers and technicians in fields of thermal science and computational fluid dynamics.

The purpose of this study is analysis of deformation and vibrations of turbine blades produced by high electrolyte pressure during electrochemical machining.

An experimental setup was designed, experiments were conducted and the obtained results were compared with the finite element results. The deformations were measured according to various flow rates of electrolyte. In finite element calculations, the pressure distribution created by the electrolyte on the blade surface was obtained in the ANSYS® (A finite element analysis software) Fluent software and transferred to the static structural where the deformation analysis was carried out. Three different parameters were examined, namely blade thickness, blade material and electrolyte pressure on blade disk caused by mass flow rate. The deformation results were compared with the gap distances between cathode and anode.

Large deformations were obtained at the free end of the blade and the most curved part of it. The appropriate pressure values for the electrolyte to be used in the production of blisk blades were proposed numerically. It has been determined that high pressure applications are not suitable for gap distance lower than 0.5 mm.

When the literature is examined, it is required that the high speed flow of the electrolyte is desired in order to remove the parts that are separated from the anode from the machining area during electrochemical machining. However, the electrolyte flowing at high speeds causes high pressure in the blisk blades, excessive deformation and vibration of the machined part, and as a result, contact of the anode with the cathode. This study provides important findings for smooth electro chemical machining at high electrolyte flows.

The material point method (MPM)is a particle-based numerical method suitable for solid–liquid simulation and large deformation problems. However, MPM is generally used in solid deformation at present, to develop a multi-physics coupling MPM; the purpose of this study is to extend the MPM to simulate the heat and fluid flow and address the thermal-hydrological (TH) coupling problems.

The porous medium was discretized into two sets of Lagrangian points, and the motion of fluid points follows the Darcy’s law. Two sets of heat transport equations were established for the heat conduction and heat exchange in the pore fluid and solid skeleton. Fractures were considered by adding the porosity gradient term in the governing equations; also a transition function was introduced to smoothen the fracture boundary.

Four cases of heat and fluid flow in porous medium and fractures were presented to verify the feasibility of the proposed method. And the effects of fractures on heat and fluid flow were investigated. Additionally, a case of geothermal extraction was solved and the importance of the interstitial convective heat transfer coefficient was analyzed.

The proposed method extends the conventional MPM, using two sets of material points and two sets of heat transport equations to simulate the heat and fluid flow and address the TH coupling problems, which can be applied in both porous medium and fractures.

The modeling of interphase forces plays a significant role in the numerical simulation of gas–liquid flow in a rotodynamic multiphase pump, which deserves detailed study.

Numerical analysis is conducted to estimate the influence of interphase forces, including drag force, lift force, virtual mass force, wall lubrication force and turbulent dispersion force.

The results show that the magnitude of the interphase forces can be sorted by: drag force > virtual mass force > lift force > turbulent dispersion force > wall lubrication force. The relations between interphase forces and velocity difference of gas–liquid flow and also the interphase forces and gas volume fraction are revealed. The distribution characteristics of interphase forces in the passages from impeller inlet to diffuser outlet are illustrated and analyzed. According to the results, apart from the drag force, the virtual mass force, lift force and turbulent dispersion force are required, whereas wall lubrication force can be neglected for numerical simulation of gas–liquid flow in a rotodynamic multiphase pump. Compared with the conventional numerical method which considers drag force only, the relative errors of predicted pressure rise and efficiency based on the proposed numerical method in account of four major forces can be reduced by 4.95 per cent and 3.00 per cent, respectively.

The numerical analysis reveals the magnitude and distribution of interphase forces inside multiphase pump, which is meaningful for the simulation and design of multiphase pump.