The two-dimensional turbulent evaporating gas-droplet two-phase flows in an afterburner diffusor of turbo-fan jet engines are simulated by the k-ε turbulence model and the particle trajectory model. Comparison of pre...The two-dimensional turbulent evaporating gas-droplet two-phase flows in an afterburner diffusor of turbo-fan jet engines are simulated by the k-ε turbulence model and the particle trajectory model. Comparison of predicted gas velocity and temperature distributions with experimental results for the cases without liquid spray shows pretty good agreement. Gas-droplet two-phase flow predictions give plausible droplet trajectories, fuel-vapor concentration distribution, gas-phase velocity and temperature field in presence of liquid droplets. One run of computation with this method is made for a particular afterburner. The results indicate that the location of the atomizers is not favorable to flame stabilization and combustion efficiency. The proposed numerical modeling can also be adopted for optimization design and performance evaluation of afterburner combustors of turbo-fan jet engines.展开更多
This paper presents a mass and momentum conservative semi-implicit finite volume(FV)scheme for complex non-hydrostatic free surface flows,interacting with moving solid obstacles.A simplified incompressible Baer-Nunzia...This paper presents a mass and momentum conservative semi-implicit finite volume(FV)scheme for complex non-hydrostatic free surface flows,interacting with moving solid obstacles.A simplified incompressible Baer-Nunziato type model is considered for two-phase flows containing a liquid phase,a solid phase,and the surrounding void.According to the so-called diffuse interface approach,the different phases and consequently the void are described by means of a scalar volume fraction function for each phase.In our numerical scheme,the dynamics of the liquid phase and the motion of the solid are decoupled.The solid is assumed to be a moving rigid body,whose motion is prescribed.Only after the advection of the solid volume fraction,the dynamics of the liquid phase is considered.As usual in semi-implicit schemes,we employ staggered Cartesian control volumes and treat the nonlinear convective terms explicitly,while the pressure terms are treated implicitly.The non-conservative products arising in the transport equation for the solid volume fraction are treated by a path-conservative approach.The resulting semi-implicit FV discretization of the mass and momentum equations leads to a mildly nonlinear system for the pressure which can be efficiently solved with a nested Newton-type technique.The time step size is only limited by the velocities of the two phases contained in the domain,and not by the gravity wave speed nor by the stiff algebraic relaxation source term,which requires an implicit discretization.The resulting semi-implicit algorithm is first validated on a set of classical incompressible Navier-Stokes test problems and later also adds a fixed and moving solid phase.展开更多
An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating c...An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating coordinate system, and continuity is conserved by a mass-weighted method to solve the filtered governing equations. In the cur- rent second-order SGS model, the SGS stress is a function of both the resolved strain-rate and rotation-rate tensors, and the model parameters are obtained from the dimensional consistency and the invariants of the strain-rate and the rotation-rate tensors. In the numerical calculation, the finite volume method is used to discretize the governing equations with a staggered grid system. The SIMPLEC algorithm is applied for the solution of the discretized governing equations. Body- fitted coordinates are used to simulate the two-phase flows in complex geometries. Finally the second-order dynamic SGS model is successfully applied to simulate the dense turbu-lent particle-liquid two-phase flows in a centrifugal impeller. The predicted pressure and velocity distributions are in good agreement with experimental results.展开更多
In this paper, a consistent projection-based streamline upwind/pressure stabilizing Petrov-Galerkin （SUPG/PSPG） extended finite element method （XFEM） is presented to model incompressible immiscible two-phase flows...In this paper, a consistent projection-based streamline upwind/pressure stabilizing Petrov-Galerkin （SUPG/PSPG） extended finite element method （XFEM） is presented to model incompressible immiscible two-phase flows. As the application of linear elements in SUPG/PSPG schemes gives rise to inconsistency in stabilization terms due to the inability to regenerate the diffusive term from viscous stresses, the numerical accuracy would deteriorate dramatically. To address this issue, projections of convection and pressure gradient terms are constructed and incorporated into the stabilization formulation in our method. This would substantially recover the consistency and free the practitioner from burdensome computations of most items in the residual. Moreover, the XFEM is employed to consider in a convenient way the fluid properties that have interfacial jumps leading to discontinuities in the velocity and pressure fields as well as the projections. A number of numerical examples are analyzed to demonstrate the complete recovery of consistency, the reproduction of interfacial discontinuities and the ability of the proposed projection-based SUPG/PSPG XFEM to model two-phase flows with open and closed interfaces.展开更多
The dense solid-phase governing equations for two-phase flows are obtained by using the kinetic theory of gas molecules.Assuming that the solid-phase velocity distributions obey the Maxwell equations,the collision ter...The dense solid-phase governing equations for two-phase flows are obtained by using the kinetic theory of gas molecules.Assuming that the solid-phase velocity distributions obey the Maxwell equations,the collision term for particles under dense two-phase flow conditions is also derived. In comparison with the governing equations of a dilute two-phase flow,the solid-particle's governing equations are developed for a dense turbulent solid-liquid flow by adopting some relevant terms from the dilute two-phase governing equations.Based on Cauchy-Helmholtz theorem and Smagorinsky model, a second-order dynamic sub-grid-scale(SGS)model,in which the sub-grid-scale stress is a function of both the strain-rate tensor and the rotation-rate tensor,is proposed to model the two-phase governing equations by applying dimension analyses.Applying the SIMPLEC algorithm and staggering grid system to the two-phase discretized governing equations and employing the slip boundary conditions on the walls,the velocity and pressure fields,and the volumetric concentration are calculated.The simulation results are in a fairly good agreement with experimental data in two operating cases in a conduit with a rectangular cross-section and these comparisons imply that these models are practical.展开更多
Modelling and simulations are conducted on velocity slip and interfacial momentum transfer for supersonic two-pha.se (gas-droplet) flow in the transient section inside and outside a Laval jet(LJ). The initial velocity...Modelling and simulations are conducted on velocity slip and interfacial momentum transfer for supersonic two-pha.se (gas-droplet) flow in the transient section inside and outside a Laval jet(LJ). The initial velocity slip between gas and droplets causes an interfacial momentum transfer flux as high as (2.0-5.0) x 104 Pa. The relaxation time corresponding to this transient process is in the range of 0.015-0.090ms for the two-phase flow formed inside the LJ and less than 0.5ms outside the LJ. It demonstrates the unique performance of this system for application to fast chemical reactions using electrically active media with a lifetime in the order of 1 ms. Through the simulations of the transient processes with initial Mach number Mg from 2.783 to 4.194 at different axial positions inside the LJ, it is found that Mg has the strongest effect on the process. The momentum flux increases as the Mach number decreases. Due to compression by the shock wave at the end of the LJ, the flow pattern becomes two dimensional and viscous outside the LJ. Laser Doppler velocirneter (LDV) measurements of droplet velocities outside the LJ are in reasonably good agreement with the results of the simulation.展开更多
The flow patterns and the void fraction related to a gas-liquid two-phase flow in a small channel are experimentally studied.The test channel is a transparent quartz glass circular channel with an inner diameter of 6....The flow patterns and the void fraction related to a gas-liquid two-phase flow in a small channel are experimentally studied.The test channel is a transparent quartz glass circular channel with an inner diameter of 6.68 mm.The working fluids are air and water and their superficial velocities range from 0.014 to 8.127 m/s and from 0.0238 to 0.556 m/s,respectively.The void fraction is determined using the flow pattern images captured by a high-speed camera,while quick closing valves are used for verification.Four flow patterns are analyzed in experiments:slug flow,bubbly flow,annular flow and stratified flow.For intermittent flows(bubbly flow and slug flow),the cross-sectional void fraction is in a borderline condition while its probability distribution function(PDF)image displays a bimodal structure.For continuous flows(annular flow and stratified flow)the cross-sectional void fraction behaves as a fluctuating continuous curve while the(PDF)image displays a single peak structure.The volumetric void fraction data are also compared with available predictive formulas,and the results show that the agreement is very good.An effort is also provided to improve the so-called Gregory and Scott model using the available data.展开更多
The gas-solid two-phase flous of the precalciner were simulated by different multiphase models,such as mixture model,the Enderium model,including mixture and dispersed,and discrete phase model（DPM）,The results of th...The gas-solid two-phase flous of the precalciner were simulated by different multiphase models,such as mixture model,the Enderium model,including mixture and dispersed,and discrete phase model（DPM）,The results of the different multiphase models were analyzed and compared.showing the rationality of the diffusion and mixture of the cenment raic meals and coal poroder some extent Moreover,the results also shose the rationality of the given inlets parameters of actual process of the precalciner.展开更多
A statistical formalism overcoming some conceptual and practical difficulties arising in existing two-phase flow (2PHF) mathematical modelling has been applied to propose a model for dilute 2PHF turbulent Hows. Phase ...A statistical formalism overcoming some conceptual and practical difficulties arising in existing two-phase flow (2PHF) mathematical modelling has been applied to propose a model for dilute 2PHF turbulent Hows. Phase interaction terms with a clear physical meaning enter the equations and the formalism provides some guidelines for the avoidance of closure assumptions or the rational approximation of these terms. Continuous phase averaged continuity, momentum, turbulent kinetic energy and turbulence dissipation rate equations have been rigorously and systematically obtained in a single step. These equations display a structure similar to that for single-phase flows. It is also assumed that dispersed phase dynamics is well described by a probability density function (pdf) equation and Eulerian continuity, momentum and fluctuating kinetic energy equations for the dispersed phase are deduced. An extension of the standard k-e turbulence model for the continuous phase is used. A gradient transport model is adopted for the dispersed phase fluctuating fluxes of momentum and kinetic energy at the non-colliding, large inertia limit. This model is then used to predict the behaviour of three axisymmetric turbulent jets of air laden with solid particles varying in size and concentration. Qualitative and quantitative numerical predictions compare reasonably well with the three different sets of experimental results, studying the influence of particle size, loading ratio and flow confinement velocity.展开更多
A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carr...A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carrier fluid in a two-phase flow The continuity, the momentum, K and εequations are modeled. In this model,the solid-liquid slip veloeites, the particle-particte interactions and the interactions between two phases are considered,The sandy water pipe turbulent flows are sueeessfuly predicted by this turbulince model.展开更多
We propose a hybrid scheme for computations of incompressible two-phase flows. The incompressible constraint has been replaced by a pressure Poisson-like equation and then the pressure is updated by the modified marke...We propose a hybrid scheme for computations of incompressible two-phase flows. The incompressible constraint has been replaced by a pressure Poisson-like equation and then the pressure is updated by the modified marker and cell method. Meanwhile, the moment equations in the incompressible Navier-Stokes equations are solved by our semidiscrete Hermite central-upwind scheme, and the interface between the two fluids is considered to be continuous and is described implicitly as the 0.5 level set of a smooth function being a smeared out Heaviside function. It is here named the hybrid scheme. Some numerical experiments are successfully carried out, which verify the desired efficiency and accuracy of our hybrid scheme.展开更多
Particle-gas two-phase flows show significantly different behaviors compared to single gas flow through a convergent-divergent nozzle. Non-equilibrium effects, thermal and velocity lag results to the inefficiency of n...Particle-gas two-phase flows show significantly different behaviors compared to single gas flow through a convergent-divergent nozzle. Non-equilibrium effects, thermal and velocity lag results to the inefficiency of nozzle performance. In the present studies, theoretical analysis and numerical simulations were carried out to investigate particle-gas flows in a C-D nozzle. Homogeneous equilibrium model that no lag in velocity and temperature occurs between particles and gas phase was used to derive mass flow rate and sound speed of multiphase flows. Two-phase flows are regarded as isentropic flows that isentropic relations can be used for homogeneous equilibrium model. Discrete phase model （DPM） where interaction with continuous phase and discrete random walk model were considered was used to calculate particle- gas flows. Particle mass loadings were varied to investigate their effects on choking phenomena of particle-gas flows. Mass flow rate and sound speed of mixture flows were theoretically calculated by homogeneous equilibrium model and compared with numerical results. Shock wave structure and particle number density were also obtained to be different at different particle mass loading and operating pressure conditions.展开更多
In the present paper, random-choice method (RCM) and second-order GRP difference method, which are high resolution methods used for pure gas flows with shocks, are extended and employed to study the problem of one-dim...In the present paper, random-choice method (RCM) and second-order GRP difference method, which are high resolution methods used for pure gas flows with shocks, are extended and employed to study the problem of one-dimensional unsteady two-phase flows. The two-phase shock wave and the flow field behind it in a dusty gas shock tube are calculated and the time-dependent change of the flow parameters for the gas and particle phase are obtained. The numerical results indicate that both the two methods can give the relaxation structure of the two-phase shocks with a sharp discontinuous front and that the GRP method has the advantages of less time-consuming and higher accuracy over the RCM method.展开更多
In this paper,we propose a novel probe to identify phases in any two-phase flows where one phase is conductive and the other nonconductive.We can further obtain many parameters such as void fraction,bubble velocity,an...In this paper,we propose a novel probe to identify phases in any two-phase flows where one phase is conductive and the other nonconductive.We can further obtain many parameters such as void fraction,bubble velocity,and interfacial area concentration.Compared with the traditional probe,the novel probe has unique advantages that it is less dependent on water conductance or distance between the electrodes,and that the amplitude is bigger between high and low levels.Theoretical analyses showed that the measurement error became higher when water conductance decreases or distance increases,which is consistent with the theoretical analyses.Experimental results showed that the output signal kept constant with salt content of 0-5% and electrode distance of 0-30 mm in tap water.The level difference was up to 6.4 V,resulting in identifying two phases easily.Time traces of phase identification were completely consistent with the flow structures.展开更多
In this study,an efficient regularized lattice Boltzmann model aimed at solving the consistent and conservative phase-field model is developed.This model is composed of the conservative Allen-Cahn equation,the momentu...In this study,an efficient regularized lattice Boltzmann model aimed at solving the consistent and conservative phase-field model is developed.This model is composed of the conservative Allen-Cahn equation,the momentum equation featuring a modified mass flux,and the associated consistency conditions.Consequently,two distribution functions are introduced within the framework of the regularized lattice Boltzmann model:one dedicated to the conservative Allen-Cahn equation,and the other designed for addressing the fluid dynamics equations.In order to accurately recover the momentum equation and ensure the consistency of mass and momentum transport,a simple force distribution function with a auxiliary term is incorporated into the regularized lattice Boltzmann model.To assess the capabilities of the current regularized lattice Boltzmann model,simulations of various two-phase flow problems with substantial density ratios have been conducted,including layered Poiseuille flow and spinodal decomposition.These simulations demonstrate excellent agreement with previously published numerical results.Additionally,numerical investigations into Rayleigh-Taylor instability indicate that the present regularized lattice Boltzmann model can accurately and stably track interfaces with high precision.展开更多
The comparative study on seven equation models with two different six equations modelfor compressible two-phase flow analysis is proposed. The seven equations model isderived for compressible two-phase flow that is in...The comparative study on seven equation models with two different six equations modelfor compressible two-phase flow analysis is proposed. The seven equations model isderived for compressible two-phase flow that is in the nonconservation form. In thepresent work, two different six equations model are derived for two pressures, two velocities and single temperature with the derivation of the equation of state. The closingequation for one of the six equations model is energy conservation equation while anotherone is closed by entropy balance equation. The partial differential form of governingequations is hyperbolic and written in the conservative form. At this point, the set ofgoverning equations are derived based on the principle of extended thermodynamics.The method of solving single temperature from both six equation models are simple anddirect solution can be obtained. Numerical simulation has been tried using one of the sixequation models for air–water shock tube problems. Explicit fourth order Runge–Kuttascheme is used with Finite Volume Shock Capturing method for solving the governingequations numerically. The pressure, velocity and volume fraction variations are captured along the shock tube length through flow solver. Experimental work is carried outto magnify the initial stage of liquid injection into a gas. The outcome of six equationsmodel for compressible two-phase flow has revealed the multi-phase flow characteristicsthat are similar to the actual conditions.展开更多
We propose a multiple relaxation time entropic realization of a recent twophase flow lattice Boltzmann model[S.A.Hosseini,B.Dorschner,and I.V.Karlin,Journal of Fluid Mechanics 953(2022)].While the original model with ...We propose a multiple relaxation time entropic realization of a recent twophase flow lattice Boltzmann model[S.A.Hosseini,B.Dorschner,and I.V.Karlin,Journal of Fluid Mechanics 953(2022)].While the original model with a single relaxation time allows us to reach large density ratios,it is limited in terms of stability with respect to non-dimensional viscosity and velocity.Herewe showthat the entropic multiple relaxation time model extends the stability limits of the model significantly,which allows us to reach larger Reynolds numbers for a given grid resolution.The thermodynamic properties of the solver,using the Peng-Robinson equation of state,are studied first using simple configurations.Co-existence densities and temperature scaling of both the interface thickness and the surface tension are shown to agree well with theory.The model is then used to simulate the impact of a drop onto a thin liquid film with density and viscosity ratios matching those of water and air both in two and three dimensions.The results are in very good agreement with theoretically predicted scaling laws and experimental data.展开更多
A compressible and multiphase flows solver has been developed for the study of liquid/gas flows involving shock waves and strong expansion waves leading to cavitation.This solver has a structure similar to those of th...A compressible and multiphase flows solver has been developed for the study of liquid/gas flows involving shock waves and strong expansion waves leading to cavitation.This solver has a structure similar to those of the one-fluid Euler solvers,differing from them by the presence of a void ratio transport-equation.The model and the system of equations to be simulated are presented.Results are displayed for shock and expansion tube problems,shock-bubble interaction and underwater explosion.Close agreement with reference solutions,obtained from explicit finite volume approaches,is demonstrated.Different numerical methods are additionally displayed to provide comparable and improved computational efficiency to the model and the system of equations.The overall procedure is therefore very well suited for use in general two-phase fluid flow simulations.展开更多
Deep-sea mineral resource transportation predominantly utilizes hydraulic pipeline methodology.Environmental factors induce vibrations in flexible pipelines,thereby affecting the internal flow characteristics.Therefor...Deep-sea mineral resource transportation predominantly utilizes hydraulic pipeline methodology.Environmental factors induce vibrations in flexible pipelines,thereby affecting the internal flow characteristics.Therefore,real-time monitoring of solid–liquid two-phase flow in pipelines is crucial for system maintenance.This study develops an autoencoder-based deep learning framework to reconstruct three-dimensional solid–liquid two-phase flow within flexible vibrating pipelines utilizing sparse wall information from sensors.Within this framework,separate X-model and F-model with distinct hidden-layer structures are established to reconstruct the coordinates and flow field information on the computational domain grid of the pipeline under traveling wave vibration.Following hyperparameter optimization,the models achieved high reconstruction accuracy,demonstrating R^(2)values of 0.990 and 0.945,respectively.The models’robustness is evaluated across three aspects:vibration parameters,physical fields,and vibration modes,demonstrating good reconstruction performance.Results concerning sensors show that 20 sensors(0.06%of total grids)achieve a balance between accuracy and cost,with superior accuracy obtained when arranged along the full length of the pipe compared to a dense arrangement at the front end.The models exhibited a signal-to-noise ratio tolerance of approximately 27 dB,with reconstruction accuracy being more affected by sensor failures at both ends of the pipeline.展开更多
This work investigated the dynamic behavior of vertical pipes conveying gas-liquid two-phase flow when subjected to external excitations at both ends.Even with minimal excitation amplitude,resonance can occur when the...This work investigated the dynamic behavior of vertical pipes conveying gas-liquid two-phase flow when subjected to external excitations at both ends.Even with minimal excitation amplitude,resonance can occur when the excitation frequency aligns with the natural frequency of the pipe,significantly increasing the degree of operational risk.The governing equation of motion based on the Euler-Bernoulli beam is derived for the relative deflection with stationary simply supported ends,with the effects of the external excitations represented by source terms distributed along the pipe length.The fourth-order partial differential equation is solved via the generalized integral transform technique(GITT),with the solution successfully verified via comparison with results in the literature.A comprehensive analysis of the vibration phenomena and changes in the motion state of the pipe is conducted for three classes of external excitation conditions:same frequency and amplitude(SFSA),same frequency but different amplitudes(SFDA),and different frequencies and amplitudes(DFDA).The numerical results show that with increasing gas volume fraction,the position corresponding to the maximum vibration displacement shifts upward.Compared with conditions without external excitation,the vibration displacement of the pipe conveying two-phase flow under external excitation increases significantly.The frequency of external excitation has a significant effect on the dynamic behavior of a pipe conveying two-phase flow.展开更多
文摘The two-dimensional turbulent evaporating gas-droplet two-phase flows in an afterburner diffusor of turbo-fan jet engines are simulated by the k-ε turbulence model and the particle trajectory model. Comparison of predicted gas velocity and temperature distributions with experimental results for the cases without liquid spray shows pretty good agreement. Gas-droplet two-phase flow predictions give plausible droplet trajectories, fuel-vapor concentration distribution, gas-phase velocity and temperature field in presence of liquid droplets. One run of computation with this method is made for a particular afterburner. The results indicate that the location of the atomizers is not favorable to flame stabilization and combustion efficiency. The proposed numerical modeling can also be adopted for optimization design and performance evaluation of afterburner combustors of turbo-fan jet engines.
基金funded by the Italian Ministry of Education,University and Research(MIUR)in the frame of the Departments of Excellence Initiative 2018-2027 attributed to DICAM of the University of Trento(grant L.232/2016)in the frame of the PRIN 2017 project Innovative numerical methods for evolutionary partial differential equations and applications,the PRIN 2022 project High order structure-preserving semi-implicit schemes for hyperbolic equations.D.is member of INdAM GNCS and was also co-funded by the European Union NextGenerationEU(PNRR,Spoke 7 CN HPC).Views and opinions expressed are however those of the author(s)only and do not necessarily reflect those of the European Union or the European Research Council.Neither the European Union nor the granting authority can be held responsible for them.
文摘This paper presents a mass and momentum conservative semi-implicit finite volume(FV)scheme for complex non-hydrostatic free surface flows,interacting with moving solid obstacles.A simplified incompressible Baer-Nunziato type model is considered for two-phase flows containing a liquid phase,a solid phase,and the surrounding void.According to the so-called diffuse interface approach,the different phases and consequently the void are described by means of a scalar volume fraction function for each phase.In our numerical scheme,the dynamics of the liquid phase and the motion of the solid are decoupled.The solid is assumed to be a moving rigid body,whose motion is prescribed.Only after the advection of the solid volume fraction,the dynamics of the liquid phase is considered.As usual in semi-implicit schemes,we employ staggered Cartesian control volumes and treat the nonlinear convective terms explicitly,while the pressure terms are treated implicitly.The non-conservative products arising in the transport equation for the solid volume fraction are treated by a path-conservative approach.The resulting semi-implicit FV discretization of the mass and momentum equations leads to a mildly nonlinear system for the pressure which can be efficiently solved with a nested Newton-type technique.The time step size is only limited by the velocities of the two phases contained in the domain,and not by the gravity wave speed nor by the stiff algebraic relaxation source term,which requires an implicit discretization.The resulting semi-implicit algorithm is first validated on a set of classical incompressible Navier-Stokes test problems and later also adds a fixed and moving solid phase.
基金the National Natural Science Foundation of China(50779069 and 90510007)the Start-up Scientific Research Foundation of China Agricultural University(2006021)the Beijing Natural Science Foundation(3071002).
文摘An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating coordinate system, and continuity is conserved by a mass-weighted method to solve the filtered governing equations. In the cur- rent second-order SGS model, the SGS stress is a function of both the resolved strain-rate and rotation-rate tensors, and the model parameters are obtained from the dimensional consistency and the invariants of the strain-rate and the rotation-rate tensors. In the numerical calculation, the finite volume method is used to discretize the governing equations with a staggered grid system. The SIMPLEC algorithm is applied for the solution of the discretized governing equations. Body- fitted coordinates are used to simulate the two-phase flows in complex geometries. Finally the second-order dynamic SGS model is successfully applied to simulate the dense turbu-lent particle-liquid two-phase flows in a centrifugal impeller. The predicted pressure and velocity distributions are in good agreement with experimental results.
文摘In this paper, a consistent projection-based streamline upwind/pressure stabilizing Petrov-Galerkin （SUPG/PSPG） extended finite element method （XFEM） is presented to model incompressible immiscible two-phase flows. As the application of linear elements in SUPG/PSPG schemes gives rise to inconsistency in stabilization terms due to the inability to regenerate the diffusive term from viscous stresses, the numerical accuracy would deteriorate dramatically. To address this issue, projections of convection and pressure gradient terms are constructed and incorporated into the stabilization formulation in our method. This would substantially recover the consistency and free the practitioner from burdensome computations of most items in the residual. Moreover, the XFEM is employed to consider in a convenient way the fluid properties that have interfacial jumps leading to discontinuities in the velocity and pressure fields as well as the projections. A number of numerical examples are analyzed to demonstrate the complete recovery of consistency, the reproduction of interfacial discontinuities and the ability of the proposed projection-based SUPG/PSPG XFEM to model two-phase flows with open and closed interfaces.
基金The project supported by the National Natural Science Foundation of China (50176022)
文摘The dense solid-phase governing equations for two-phase flows are obtained by using the kinetic theory of gas molecules.Assuming that the solid-phase velocity distributions obey the Maxwell equations,the collision term for particles under dense two-phase flow conditions is also derived. In comparison with the governing equations of a dilute two-phase flow,the solid-particle's governing equations are developed for a dense turbulent solid-liquid flow by adopting some relevant terms from the dilute two-phase governing equations.Based on Cauchy-Helmholtz theorem and Smagorinsky model, a second-order dynamic sub-grid-scale(SGS)model,in which the sub-grid-scale stress is a function of both the strain-rate tensor and the rotation-rate tensor,is proposed to model the two-phase governing equations by applying dimension analyses.Applying the SIMPLEC algorithm and staggering grid system to the two-phase discretized governing equations and employing the slip boundary conditions on the walls,the velocity and pressure fields,and the volumetric concentration are calculated.The simulation results are in a fairly good agreement with experimental data in two operating cases in a conduit with a rectangular cross-section and these comparisons imply that these models are practical.
基金Supported by the National Natural Science Foundation of China (No. 29876022) and Grant of State Key Laboratory of High Speed Hydrodynamics (No. 2007).
文摘Modelling and simulations are conducted on velocity slip and interfacial momentum transfer for supersonic two-pha.se (gas-droplet) flow in the transient section inside and outside a Laval jet(LJ). The initial velocity slip between gas and droplets causes an interfacial momentum transfer flux as high as (2.0-5.0) x 104 Pa. The relaxation time corresponding to this transient process is in the range of 0.015-0.090ms for the two-phase flow formed inside the LJ and less than 0.5ms outside the LJ. It demonstrates the unique performance of this system for application to fast chemical reactions using electrically active media with a lifetime in the order of 1 ms. Through the simulations of the transient processes with initial Mach number Mg from 2.783 to 4.194 at different axial positions inside the LJ, it is found that Mg has the strongest effect on the process. The momentum flux increases as the Mach number decreases. Due to compression by the shock wave at the end of the LJ, the flow pattern becomes two dimensional and viscous outside the LJ. Laser Doppler velocirneter (LDV) measurements of droplet velocities outside the LJ are in reasonably good agreement with the results of the simulation.
基金This work was supported by the Guangdong Basic and Applied Basic Research Foundation(2019A1515111116)Key R&D Program of Shandong Province(Nos.2019GSF109051,2019GGX101030)+1 种基金Shandong Provincial Postdoctoral Innovation Project(No.201902002)Foundation of Shandong University for Young Scholar’s Future Plans.
文摘The flow patterns and the void fraction related to a gas-liquid two-phase flow in a small channel are experimentally studied.The test channel is a transparent quartz glass circular channel with an inner diameter of 6.68 mm.The working fluids are air and water and their superficial velocities range from 0.014 to 8.127 m/s and from 0.0238 to 0.556 m/s,respectively.The void fraction is determined using the flow pattern images captured by a high-speed camera,while quick closing valves are used for verification.Four flow patterns are analyzed in experiments:slug flow,bubbly flow,annular flow and stratified flow.For intermittent flows(bubbly flow and slug flow),the cross-sectional void fraction is in a borderline condition while its probability distribution function(PDF)image displays a bimodal structure.For continuous flows(annular flow and stratified flow)the cross-sectional void fraction behaves as a fluctuating continuous curve while the(PDF)image displays a single peak structure.The volumetric void fraction data are also compared with available predictive formulas,and the results show that the agreement is very good.An effort is also provided to improve the so-called Gregory and Scott model using the available data.
文摘The gas-solid two-phase flous of the precalciner were simulated by different multiphase models,such as mixture model,the Enderium model,including mixture and dispersed,and discrete phase model（DPM）,The results of the different multiphase models were analyzed and compared.showing the rationality of the diffusion and mixture of the cenment raic meals and coal poroder some extent Moreover,the results also shose the rationality of the given inlets parameters of actual process of the precalciner.
基金Supported by the Spanish CICYTR &D National Programs,under contract PB91-0699.
文摘A statistical formalism overcoming some conceptual and practical difficulties arising in existing two-phase flow (2PHF) mathematical modelling has been applied to propose a model for dilute 2PHF turbulent Hows. Phase interaction terms with a clear physical meaning enter the equations and the formalism provides some guidelines for the avoidance of closure assumptions or the rational approximation of these terms. Continuous phase averaged continuity, momentum, turbulent kinetic energy and turbulence dissipation rate equations have been rigorously and systematically obtained in a single step. These equations display a structure similar to that for single-phase flows. It is also assumed that dispersed phase dynamics is well described by a probability density function (pdf) equation and Eulerian continuity, momentum and fluctuating kinetic energy equations for the dispersed phase are deduced. An extension of the standard k-e turbulence model for the continuous phase is used. A gradient transport model is adopted for the dispersed phase fluctuating fluxes of momentum and kinetic energy at the non-colliding, large inertia limit. This model is then used to predict the behaviour of three axisymmetric turbulent jets of air laden with solid particles varying in size and concentration. Qualitative and quantitative numerical predictions compare reasonably well with the three different sets of experimental results, studying the influence of particle size, loading ratio and flow confinement velocity.
文摘A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carrier fluid in a two-phase flow The continuity, the momentum, K and εequations are modeled. In this model,the solid-liquid slip veloeites, the particle-particte interactions and the interactions between two phases are considered,The sandy water pipe turbulent flows are sueeessfuly predicted by this turbulince model.
文摘We propose a hybrid scheme for computations of incompressible two-phase flows. The incompressible constraint has been replaced by a pressure Poisson-like equation and then the pressure is updated by the modified marker and cell method. Meanwhile, the moment equations in the incompressible Navier-Stokes equations are solved by our semidiscrete Hermite central-upwind scheme, and the interface between the two fluids is considered to be continuous and is described implicitly as the 0.5 level set of a smooth function being a smeared out Heaviside function. It is here named the hybrid scheme. Some numerical experiments are successfully carried out, which verify the desired efficiency and accuracy of our hybrid scheme.
文摘Particle-gas two-phase flows show significantly different behaviors compared to single gas flow through a convergent-divergent nozzle. Non-equilibrium effects, thermal and velocity lag results to the inefficiency of nozzle performance. In the present studies, theoretical analysis and numerical simulations were carried out to investigate particle-gas flows in a C-D nozzle. Homogeneous equilibrium model that no lag in velocity and temperature occurs between particles and gas phase was used to derive mass flow rate and sound speed of multiphase flows. Two-phase flows are regarded as isentropic flows that isentropic relations can be used for homogeneous equilibrium model. Discrete phase model （DPM） where interaction with continuous phase and discrete random walk model were considered was used to calculate particle- gas flows. Particle mass loadings were varied to investigate their effects on choking phenomena of particle-gas flows. Mass flow rate and sound speed of mixture flows were theoretically calculated by homogeneous equilibrium model and compared with numerical results. Shock wave structure and particle number density were also obtained to be different at different particle mass loading and operating pressure conditions.
文摘In the present paper, random-choice method (RCM) and second-order GRP difference method, which are high resolution methods used for pure gas flows with shocks, are extended and employed to study the problem of one-dimensional unsteady two-phase flows. The two-phase shock wave and the flow field behind it in a dusty gas shock tube are calculated and the time-dependent change of the flow parameters for the gas and particle phase are obtained. The numerical results indicate that both the two methods can give the relaxation structure of the two-phase shocks with a sharp discontinuous front and that the GRP method has the advantages of less time-consuming and higher accuracy over the RCM method.
基金Support by the National Natural Science Foundation of China (No. 51006068)
文摘In this paper,we propose a novel probe to identify phases in any two-phase flows where one phase is conductive and the other nonconductive.We can further obtain many parameters such as void fraction,bubble velocity,and interfacial area concentration.Compared with the traditional probe,the novel probe has unique advantages that it is less dependent on water conductance or distance between the electrodes,and that the amplitude is bigger between high and low levels.Theoretical analyses showed that the measurement error became higher when water conductance decreases or distance increases,which is consistent with the theoretical analyses.Experimental results showed that the output signal kept constant with salt content of 0-5% and electrode distance of 0-30 mm in tap water.The level difference was up to 6.4 V,resulting in identifying two phases easily.Time traces of phase identification were completely consistent with the flow structures.
基金supported by the Natural Science Foundation of China(Grant No.12372286)Excellent Youth Foundation of Changsha Scientific Committee(Grant No.kq2306023)+3 种基金the Research Foundation of Education Bureau of Hunan Province,China(Grant No.20A097)supported by the Natural Science Foundation of China(Grant No.12301554)the Natural Science Foundation of Jiangsu Province(Grant No.BK20220368)Natural Science Fund for Colleges and Universities in Jiangsu Province(Grant No.22KJD110003).
文摘In this study,an efficient regularized lattice Boltzmann model aimed at solving the consistent and conservative phase-field model is developed.This model is composed of the conservative Allen-Cahn equation,the momentum equation featuring a modified mass flux,and the associated consistency conditions.Consequently,two distribution functions are introduced within the framework of the regularized lattice Boltzmann model:one dedicated to the conservative Allen-Cahn equation,and the other designed for addressing the fluid dynamics equations.In order to accurately recover the momentum equation and ensure the consistency of mass and momentum transport,a simple force distribution function with a auxiliary term is incorporated into the regularized lattice Boltzmann model.To assess the capabilities of the current regularized lattice Boltzmann model,simulations of various two-phase flow problems with substantial density ratios have been conducted,including layered Poiseuille flow and spinodal decomposition.These simulations demonstrate excellent agreement with previously published numerical results.Additionally,numerical investigations into Rayleigh-Taylor instability indicate that the present regularized lattice Boltzmann model can accurately and stably track interfaces with high precision.
文摘The comparative study on seven equation models with two different six equations modelfor compressible two-phase flow analysis is proposed. The seven equations model isderived for compressible two-phase flow that is in the nonconservation form. In thepresent work, two different six equations model are derived for two pressures, two velocities and single temperature with the derivation of the equation of state. The closingequation for one of the six equations model is energy conservation equation while anotherone is closed by entropy balance equation. The partial differential form of governingequations is hyperbolic and written in the conservative form. At this point, the set ofgoverning equations are derived based on the principle of extended thermodynamics.The method of solving single temperature from both six equation models are simple anddirect solution can be obtained. Numerical simulation has been tried using one of the sixequation models for air–water shock tube problems. Explicit fourth order Runge–Kuttascheme is used with Finite Volume Shock Capturing method for solving the governingequations numerically. The pressure, velocity and volume fraction variations are captured along the shock tube length through flow solver. Experimental work is carried outto magnify the initial stage of liquid injection into a gas. The outcome of six equationsmodel for compressible two-phase flow has revealed the multi-phase flow characteristicsthat are similar to the actual conditions.
基金supported by European Research Council(ERC)Advanced Grant no.834763-PonD(S.A.H,B.D.and I.K.)the Swiss National Science Foundation(SNSF)grant No.200021-172640(S.A.H.)。
文摘We propose a multiple relaxation time entropic realization of a recent twophase flow lattice Boltzmann model[S.A.Hosseini,B.Dorschner,and I.V.Karlin,Journal of Fluid Mechanics 953(2022)].While the original model with a single relaxation time allows us to reach large density ratios,it is limited in terms of stability with respect to non-dimensional viscosity and velocity.Herewe showthat the entropic multiple relaxation time model extends the stability limits of the model significantly,which allows us to reach larger Reynolds numbers for a given grid resolution.The thermodynamic properties of the solver,using the Peng-Robinson equation of state,are studied first using simple configurations.Co-existence densities and temperature scaling of both the interface thickness and the surface tension are shown to agree well with theory.The model is then used to simulate the impact of a drop onto a thin liquid film with density and viscosity ratios matching those of water and air both in two and three dimensions.The results are in very good agreement with theoretically predicted scaling laws and experimental data.
基金The authors gratefully thank K.Tang and A.Beccantini fromthe Commissariata l’Energie Atomique for having provided the numerical solutions computed with their sevenequation model.The second author would like to particularly acknowledge the support provided by the German Jordanian University through the project SEED-SNRE 7-2014.
文摘A compressible and multiphase flows solver has been developed for the study of liquid/gas flows involving shock waves and strong expansion waves leading to cavitation.This solver has a structure similar to those of the one-fluid Euler solvers,differing from them by the presence of a void ratio transport-equation.The model and the system of equations to be simulated are presented.Results are displayed for shock and expansion tube problems,shock-bubble interaction and underwater explosion.Close agreement with reference solutions,obtained from explicit finite volume approaches,is demonstrated.Different numerical methods are additionally displayed to provide comparable and improved computational efficiency to the model and the system of equations.The overall procedure is therefore very well suited for use in general two-phase fluid flow simulations.
基金financial support by the National Natural Science Foundation of China (Nos.52471293 and 12372270)the National Youth Science Foundation of China (Nos.52101322 and 52108375)+3 种基金the Program for Intergovernmental International S&T Cooperation Projects of Shanghai Municipality, China (Nos.24510711100 and 22160710200)The Oceanic Interdisciplinary Program of Shanghai Jiao Tong University (No.SL2022PT101)funded by the Open Fund of the State Key Laboratory of Coastal and Offshore Engineering of Dalian University of Technology (No.LP2415)National Key R&D Program of China (No.2023YFC2811600)
文摘Deep-sea mineral resource transportation predominantly utilizes hydraulic pipeline methodology.Environmental factors induce vibrations in flexible pipelines,thereby affecting the internal flow characteristics.Therefore,real-time monitoring of solid–liquid two-phase flow in pipelines is crucial for system maintenance.This study develops an autoencoder-based deep learning framework to reconstruct three-dimensional solid–liquid two-phase flow within flexible vibrating pipelines utilizing sparse wall information from sensors.Within this framework,separate X-model and F-model with distinct hidden-layer structures are established to reconstruct the coordinates and flow field information on the computational domain grid of the pipeline under traveling wave vibration.Following hyperparameter optimization,the models achieved high reconstruction accuracy,demonstrating R^(2)values of 0.990 and 0.945,respectively.The models’robustness is evaluated across three aspects:vibration parameters,physical fields,and vibration modes,demonstrating good reconstruction performance.Results concerning sensors show that 20 sensors(0.06%of total grids)achieve a balance between accuracy and cost,with superior accuracy obtained when arranged along the full length of the pipe compared to a dense arrangement at the front end.The models exhibited a signal-to-noise ratio tolerance of approximately 27 dB,with reconstruction accuracy being more affected by sensor failures at both ends of the pipeline.
基金financially supported by the Key Research and Development Program of Shandong Province(Grant Nos.2022CXGC020405,2023CXGC010415 and 2025TSGCCZZB0238)the National Natural Science Foundation of China(Grant No.52171288)the financial support from CNPq,FAPERJ,ANP,Embrapii,and China National Petroleum Corporation(CNPC).
文摘This work investigated the dynamic behavior of vertical pipes conveying gas-liquid two-phase flow when subjected to external excitations at both ends.Even with minimal excitation amplitude,resonance can occur when the excitation frequency aligns with the natural frequency of the pipe,significantly increasing the degree of operational risk.The governing equation of motion based on the Euler-Bernoulli beam is derived for the relative deflection with stationary simply supported ends,with the effects of the external excitations represented by source terms distributed along the pipe length.The fourth-order partial differential equation is solved via the generalized integral transform technique(GITT),with the solution successfully verified via comparison with results in the literature.A comprehensive analysis of the vibration phenomena and changes in the motion state of the pipe is conducted for three classes of external excitation conditions:same frequency and amplitude(SFSA),same frequency but different amplitudes(SFDA),and different frequencies and amplitudes(DFDA).The numerical results show that with increasing gas volume fraction,the position corresponding to the maximum vibration displacement shifts upward.Compared with conditions without external excitation,the vibration displacement of the pipe conveying two-phase flow under external excitation increases significantly.The frequency of external excitation has a significant effect on the dynamic behavior of a pipe conveying two-phase flow.
