For solving two-dimensional incompressible flow in the vorticity form by the fourth-order compact finite difference scheme and explicit strong stability preserving temporal discretizations,we show that the simple boun...For solving two-dimensional incompressible flow in the vorticity form by the fourth-order compact finite difference scheme and explicit strong stability preserving temporal discretizations,we show that the simple bound-preserving limiter in Li et al.(SIAM J Numer Anal 56:3308–3345,2018)can enforce the strict bounds of the vorticity,if the velocity field satisfies a discrete divergence free constraint.For reducing oscillations,a modified TVB limiter adapted from Cockburn and Shu(SIAM J Numer Anal 31:607–627,1994)is constructed without affecting the bound-preserving property.This bound-preserving finite difference method can be used for any passive convection equation with a divergence free velocity field.展开更多
Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple the...Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.展开更多
The numerical solution of compressible flows has become more prevalent than that of incompressible flows.With the help of the artificial compressibility approach,incompressible flows can be solved numerically using th...The numerical solution of compressible flows has become more prevalent than that of incompressible flows.With the help of the artificial compressibility approach,incompressible flows can be solved numerically using the same methods as compressible ones.The artificial compressibility scheme is thus widely used to numerically solve incompressible Navier-Stokes equations.Any numerical method highly depends on its accuracy and speed of convergence.Although the artificial compressibility approach is utilized in several numerical simulations,the effect of the compressibility factor on the accuracy of results and convergence speed has not been investigated for nanofluid flows in previous studies.Therefore,this paper assesses the effect of this factor on the convergence speed and accuracy of results for various types of thermo-flow.To improve the stability and convergence speed of time discretizations,the fifth-order Runge-Kutta method is applied.A computer program has been written in FORTRAN to solve the discretized equations in different Reynolds and Grashof numbers for various grids.The results demonstrate that the artificial compressibility factor has a noticeable effect on the accuracy and convergence rate of the simulation.The optimum artificial compressibility is found to be between 1 and 5.These findings can be utilized to enhance the performance of commercial numerical simulation tools,including ANSYS and COMSOL.展开更多
Physics-informed deep learning has drawn tremendous interest in recent years to solve computational physics problems,whose basic concept is to embed physical laws to constrain/inform neural networks,with the need of l...Physics-informed deep learning has drawn tremendous interest in recent years to solve computational physics problems,whose basic concept is to embed physical laws to constrain/inform neural networks,with the need of less data for training a reliable model.This can be achieved by incorporating the residual of physics equations into the loss function.Through minimizing the loss function,the network could approximate the solution.In this paper,we propose a mixed-variable scheme of physics-informed neural network(PINN)for fluid dynamics and apply it to simulate steady and transient laminar flows at low Reynolds numbers.A parametric study indicates that the mixed-variable scheme can improve the PINN trainability and the solution accuracy.The predicted velocity and pressure fields by the proposed PINN approach are also compared with the reference numerical solutions.Simulation results demonstrate great potential of the proposed PINN for fluid flow simulation with a high accuracy.展开更多
Incompressible viscous flows on curved surfaces are considered with respect to the interplay of surface geometry, curvature, and vorticity dynamics. Free flows and cylindrical wakes over a Gaussian bump are numericall...Incompressible viscous flows on curved surfaces are considered with respect to the interplay of surface geometry, curvature, and vorticity dynamics. Free flows and cylindrical wakes over a Gaussian bump are numerically solved using a surface vorticity- stream function formulation. Numerical simulations show that the Gaussian curvature can generate vorticity, and non-uniformity of the Gaussian curvature is the main cause. In the cylindrical wake, the bump dominated by the positive Gaussian curvature can significantly affect the vortex street by forming velocity depression and changing vorticity transport. The results may provide possibilities for manipulating surface flows through local change in the surface geometry.展开更多
In this paper, we propose a novel incompressible finite-difference lattice Boltzmann Equation (FDLBE). Because source terms that reflect the interaction between phases can be accurately described, the new model is s...In this paper, we propose a novel incompressible finite-difference lattice Boltzmann Equation (FDLBE). Because source terms that reflect the interaction between phases can be accurately described, the new model is suitable for simulating two-way coupling incompressible multiphase flow The 2-D particle-laden flow over a backward-facing step is chosen as a test case to validate the present method. Favorable results are obtained and the present scheme is shown to have good prospects in practical applications.展开更多
A combined characteristic-based split algorithm and all adaptive meshing technique for analyzing two-dimensional viscous incompressible flow are presented. Tile method uses the three-node triangular element with equal...A combined characteristic-based split algorithm and all adaptive meshing technique for analyzing two-dimensional viscous incompressible flow are presented. Tile method uses the three-node triangular element with equal-order interpolation functions for all variables of tile velocity components and pressure. The main advantage of the combined nlethod is that it inlproves the sohltion accuracy by coupling an error estinla- tion procedure to an adaptive meshing technique that generates small elements in regions with a large change ill sohmtion gradients, mid at the same time, larger elements in the other regions. The performance of the combined procedure is evaluated by analyzing one test case of the flow past a cylinder, for their transient and steady-state flow behaviors.展开更多
In this paper we derive LPS's criterion for the breakdown of classical solutions to the incompressible nematic liquid crystal flow, a simplified version of Ericksen-Leslie system modeling the hydrodynamic evolution o...In this paper we derive LPS's criterion for the breakdown of classical solutions to the incompressible nematic liquid crystal flow, a simplified version of Ericksen-Leslie system modeling the hydrodynamic evolution of nematic liquid crystals in R^3. We show that if 0 〈 T 〈 +∞ is the maximal time interval for the unique smooth solution u ∈ C^∞([0, T),R^3),then |u|+|△d|∈L^q([0,T],L^p(R^3)),where p and q satisfy the Ladyzhenskaya-Prodi-Serrin's condition:3/p+2/q=1 and p∈(3,+∞].展开更多
The short-range property of interactions between scales in incompressible turbulent flow was examined. Some formulae for the short-range eddy stress were given. A concept of resonant-range interactions between extreme...The short-range property of interactions between scales in incompressible turbulent flow was examined. Some formulae for the short-range eddy stress were given. A concept of resonant-range interactions between extremely contiguous scales was introduced and some formulae for the resonant-range eddy stress were also derived. Multi-scale equations for the incompressible turbulent flows were proposed. Key words turbulence - incompressible flow - interactions between scales - multi-scale equations MSC 2000 76F70展开更多
A method combining the immersed boundary technique and a multi- relaxation-time (MRT) lattice Boltzmann flux solver (LBFS) is presented for numerical simulation of incompressible flows over circular and elliptic c...A method combining the immersed boundary technique and a multi- relaxation-time (MRT) lattice Boltzmann flux solver (LBFS) is presented for numerical simulation of incompressible flows over circular and elliptic cylinders and NACA 0012 Airfoil. The method uses a simple Cartesian mesh to simulate flows past immersed complicated bodies. With the Chapman-Enskog expansion analysis, a transform is performed between the Navier-Stokes and lattice Boltzmann equations (LBEs). The LBFS is used to discretize the macroscopic differential equations with a finite volume method and evaluate the interface fluxes through local reconstruction of the lattice Boltzmann solution. The immersed boundary technique is used to correct the intermediate velocity around the solid boundary to satisfy the no-slip boundary condition. Agreement of simulation results with the data found in the literature shows reliability of the proposed method in simulating laminar flows on a Cartesian mesh.展开更多
This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional (2D) incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fra...This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional (2D) incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling. The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by join- ing the centroid of cells sharing the common vertex. For the temporal integration of the momentum equations, an im- plicit second-order scheme is utilized to enhance the com- putational stability and eliminate the time step limit due to the diffusion term. The momentum equations are discretized by the vertex-centered finite volume method (FVM) and the pressure Poisson equation is solved by the Galerkin finite el- ement method (FEM). The momentum interpolation is used to damp out the spurious pressure wiggles. The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both veloc- ity and pressure. The classic test cases, the lid-driven cavity flow, the skew cavity flow and the backward-facing step flow, show that numerical results are in good agreement with the published benchmark solutions.展开更多
The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation (DNS) of incompressible turbulent reacting channel flows. The instantaneous DNS results show the n...The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation (DNS) of incompressible turbulent reacting channel flows. The instantaneous DNS results show the near-wall strip structures of concentration and temperature fluctuations. The DNS statistical results give the budget of the terms in the correlation equations, showing that the production and dissipation terms are most important. The DNS statistical data are used to validate the closure model in RANS second-order moment (SOM) combustion model. It is found that the simulated diffusion and production terms are in agreement with the DNS data in most flow regions, except in the near-wall region, where the near-wall modification should be made, and the closure model for the dissipation term needs further improvement. The algebraic second-order moment (ASOM) combustion model is well validated by DNS.展开更多
We present a cut-cell method for the simulation of 2D incompressible flows past obstacles.It consists in using the MAC scheme on cartesian grids and imposing Dirchlet boundary conditions for the velocity field on the ...We present a cut-cell method for the simulation of 2D incompressible flows past obstacles.It consists in using the MAC scheme on cartesian grids and imposing Dirchlet boundary conditions for the velocity field on the boundary of solid structures following the Shortley-Weller formulation.In order to ensure local conservation properties,viscous and convecting terms are discretized in a finite volume way.The scheme is second order implicit in time for the linear part,the linear systems are solved by the use of the capacitance matrix method for non-moving obstacles.Numerical results of flows around an impulsively started circular cylinder are presented which confirm the efficiency of the method,for Reynolds numbers 1000 and 3000.An example of flows around a moving rigid body at Reynolds number 800 is also shown,a solver using the PETSc-Library has been prefered in this context to solve the linear systems.展开更多
A unified numerical scheme for the solutions of the compressible and incompressible Navier-Stokes equations is investigated based on a time-derivative preconditioning algorithm. The primitive variables are pressure, v...A unified numerical scheme for the solutions of the compressible and incompressible Navier-Stokes equations is investigated based on a time-derivative preconditioning algorithm. The primitive variables are pressure, velocities and temperature. The time integration scheme is used in conjunction with a finite volume discretization. The preconditioning is coupled with a high order implicit upwind scheme based on the definition of a Roe's type matrix. Computational capabilities are demonstrated through computations of high Mach number, middle Mach number, very low Mach number, and incompressible flow. It has also been demonstrated that the discontinuous surface in flow field can be captured for the implementation Roe's scheme.展开更多
Three exact solutions are obtained for 2-D incompressible potential flows around two moving circles in three cases: (i) expansion (or contraction) of themselves, (ii) approaching (or departing from) each other, (iii) ...Three exact solutions are obtained for 2-D incompressible potential flows around two moving circles in three cases: (i) expansion (or contraction) of themselves, (ii) approaching (or departing from) each other, (iii) moving perpendicularly to the line connecting the centres in opposite directions. Meanwhile, an- other set of two exact solutions is obtained for 2-D incompressible potential flows between two moving eccen- tric circles in two cases: moving parallelly or perpendicularly to the line connecting the centres.展开更多
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.展开更多
A simple method is proposed, for incremental static analysis of a set of inter-colliding particles, simulating 2D flow. Within each step of proposed algorithm, the particles perform small displacements, proportional t...A simple method is proposed, for incremental static analysis of a set of inter-colliding particles, simulating 2D flow. Within each step of proposed algorithm, the particles perform small displacements, proportional to the out-of-balance forces, acting on them. Numerical experiments show that if the liquid is confined within boundaries of a set of inter-communicating vessels, then the proposed method converges to a final equilibrium state. This incremental static analysis approximates dynamic behavior with strong damping and can provide information, as a first approximation to 2D movement of a liquid. In the initial arrangement of particles, a rhombic element is proposed, which assures satisfactory incompressibility of the fluid. Based on the proposed algorithm, a simple and short computer program (a “pocket” program) has been developed, with only about 120 Fortran instructions. This program is first applied to an amount of liquid, contained in a single vessel. A coarse and refined discretization is tried. In final equilibrium state of liquid, the distribution on hydro-static pressure on vessel boundaries, obtained by proposed computational model, is found in satisfactory approximation with corresponding theoretical data. Then, an opening is formed, at the bottom of a vertical boundary of initial vessel, and the liquid is allowed to flow gradually to an adjacent vessel. Almost whole amount of liquid is transferred, from first to second vessel, except of few drops-particles, which remain, in equilibrium, at the bottom of initial vessel. In the final equilibrium state of liquid, in the second vessel, the free surface level of the liquid confirms that the proposed rhombing element assures a satisfactory incompressibility of the fluid.展开更多
Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method(GAHM) is presented for incompressible two-phase flow.A front tracking method(FTM) is used as the skeleto...Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method(GAHM) is presented for incompressible two-phase flow.A front tracking method(FTM) is used as the skeleton of the GAHM for low mass loss and resources.Smooth eulerian level set values are calculated from the FTM interface,and are used for a local interface reconstruction.The reconstruction avoids marker particle redistribution and enables an automatic treatment of interfacial topology change.The cubic Hermit interpolation is employed in all steps of the GAHM to capture subgrid structures within a single spacial cell.The performance of the GAHM is carefully evaluated in a benchmark test.Results show significant improvements of mass loss,clear subgrid structures,highly accurate derivatives(normals and curvatures) and low cost.The GAHM is further coupled with an incompressible multiphase flow solver,Super CE/SE,for more complex and practical applications.The updated solver is evaluated through comparison with an early droplet research.展开更多
Generally the incompressible viscous flow problem is described by the Navier-Stokes equation. Based on the weighted residual method the discrete formulation of element-free Galerkin is inferred in this paper. By the s...Generally the incompressible viscous flow problem is described by the Navier-Stokes equation. Based on the weighted residual method the discrete formulation of element-free Galerkin is inferred in this paper. By the step-bystep computation in the field of time, and adopting the least-square estimation of the-same-order shift, this paper has calculated both velocity and pressure from the decoupling independent equations. Each time fraction Newton-Raphson iterative method is applied for the velocity and pressure. Finally, this paper puts the method into practice of the shear-drive cavity flow, verifying the validity, high accuracy and stability.展开更多
A new scalar projection method presented for simulating incompressible flows with variable density is proposed. It reverses conventional projection algorithm by computing first the irrotational component of the veloci...A new scalar projection method presented for simulating incompressible flows with variable density is proposed. It reverses conventional projection algorithm by computing first the irrotational component of the velocity and then the pressure. The first phase of the projection is purely kinematics. The predicted velocity field is subjected to a discrete Hodge-Helmholtz decomposition. The second phase of upgrade of pressure from the density uses Stokes’ theorem to explicitly compute the pressure. If all or part of the boundary conditions is then fixed on the divergence free physical field, the system required to be solved for the scalar potential of velocity becomes a Poisson equation with constant coefficients fitted with Dirichlet conditions.展开更多
文摘For solving two-dimensional incompressible flow in the vorticity form by the fourth-order compact finite difference scheme and explicit strong stability preserving temporal discretizations,we show that the simple bound-preserving limiter in Li et al.(SIAM J Numer Anal 56:3308–3345,2018)can enforce the strict bounds of the vorticity,if the velocity field satisfies a discrete divergence free constraint.For reducing oscillations,a modified TVB limiter adapted from Cockburn and Shu(SIAM J Numer Anal 31:607–627,1994)is constructed without affecting the bound-preserving property.This bound-preserving finite difference method can be used for any passive convection equation with a divergence free velocity field.
基金supported by the NSFC Grant no.12271492the Natural Science Foundation of Henan Province of China Grant no.222300420550+1 种基金supported by the NSFC Grant no.12271498the National Key R&D Program of China Grant no.2022YFA1005202/2022YFA1005200.
文摘Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.
基金The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through the Large Groups Project under grant number RGP.2/235/43.
文摘The numerical solution of compressible flows has become more prevalent than that of incompressible flows.With the help of the artificial compressibility approach,incompressible flows can be solved numerically using the same methods as compressible ones.The artificial compressibility scheme is thus widely used to numerically solve incompressible Navier-Stokes equations.Any numerical method highly depends on its accuracy and speed of convergence.Although the artificial compressibility approach is utilized in several numerical simulations,the effect of the compressibility factor on the accuracy of results and convergence speed has not been investigated for nanofluid flows in previous studies.Therefore,this paper assesses the effect of this factor on the convergence speed and accuracy of results for various types of thermo-flow.To improve the stability and convergence speed of time discretizations,the fifth-order Runge-Kutta method is applied.A computer program has been written in FORTRAN to solve the discretized equations in different Reynolds and Grashof numbers for various grids.The results demonstrate that the artificial compressibility factor has a noticeable effect on the accuracy and convergence rate of the simulation.The optimum artificial compressibility is found to be between 1 and 5.These findings can be utilized to enhance the performance of commercial numerical simulation tools,including ANSYS and COMSOL.
文摘Physics-informed deep learning has drawn tremendous interest in recent years to solve computational physics problems,whose basic concept is to embed physical laws to constrain/inform neural networks,with the need of less data for training a reliable model.This can be achieved by incorporating the residual of physics equations into the loss function.Through minimizing the loss function,the network could approximate the solution.In this paper,we propose a mixed-variable scheme of physics-informed neural network(PINN)for fluid dynamics and apply it to simulate steady and transient laminar flows at low Reynolds numbers.A parametric study indicates that the mixed-variable scheme can improve the PINN trainability and the solution accuracy.The predicted velocity and pressure fields by the proposed PINN approach are also compared with the reference numerical solutions.Simulation results demonstrate great potential of the proposed PINN for fluid flow simulation with a high accuracy.
基金supported by the National Natural Science Foundation of China(Nos.11472082 and11172069)
文摘Incompressible viscous flows on curved surfaces are considered with respect to the interplay of surface geometry, curvature, and vorticity dynamics. Free flows and cylindrical wakes over a Gaussian bump are numerically solved using a surface vorticity- stream function formulation. Numerical simulations show that the Gaussian curvature can generate vorticity, and non-uniformity of the Gaussian curvature is the main cause. In the cylindrical wake, the bump dominated by the positive Gaussian curvature can significantly affect the vortex street by forming velocity depression and changing vorticity transport. The results may provide possibilities for manipulating surface flows through local change in the surface geometry.
基金The project supported by the National Natural Science Foundation of China(60073044)the State Key Development Programme for Basic Research of China(G1990022207).
文摘In this paper, we propose a novel incompressible finite-difference lattice Boltzmann Equation (FDLBE). Because source terms that reflect the interaction between phases can be accurately described, the new model is suitable for simulating two-way coupling incompressible multiphase flow The 2-D particle-laden flow over a backward-facing step is chosen as a test case to validate the present method. Favorable results are obtained and the present scheme is shown to have good prospects in practical applications.
文摘A combined characteristic-based split algorithm and all adaptive meshing technique for analyzing two-dimensional viscous incompressible flow are presented. Tile method uses the three-node triangular element with equal-order interpolation functions for all variables of tile velocity components and pressure. The main advantage of the combined nlethod is that it inlproves the sohltion accuracy by coupling an error estinla- tion procedure to an adaptive meshing technique that generates small elements in regions with a large change ill sohmtion gradients, mid at the same time, larger elements in the other regions. The performance of the combined procedure is evaluated by analyzing one test case of the flow past a cylinder, for their transient and steady-state flow behaviors.
基金Supported by National Natural Science Foundation of China (10976026, 11271305, 11301439, 11226174)
文摘In this paper we derive LPS's criterion for the breakdown of classical solutions to the incompressible nematic liquid crystal flow, a simplified version of Ericksen-Leslie system modeling the hydrodynamic evolution of nematic liquid crystals in R^3. We show that if 0 〈 T 〈 +∞ is the maximal time interval for the unique smooth solution u ∈ C^∞([0, T),R^3),then |u|+|△d|∈L^q([0,T],L^p(R^3)),where p and q satisfy the Ladyzhenskaya-Prodi-Serrin's condition:3/p+2/q=1 and p∈(3,+∞].
文摘The short-range property of interactions between scales in incompressible turbulent flow was examined. Some formulae for the short-range eddy stress were given. A concept of resonant-range interactions between extremely contiguous scales was introduced and some formulae for the resonant-range eddy stress were also derived. Multi-scale equations for the incompressible turbulent flows were proposed. Key words turbulence - incompressible flow - interactions between scales - multi-scale equations MSC 2000 76F70
文摘A method combining the immersed boundary technique and a multi- relaxation-time (MRT) lattice Boltzmann flux solver (LBFS) is presented for numerical simulation of incompressible flows over circular and elliptic cylinders and NACA 0012 Airfoil. The method uses a simple Cartesian mesh to simulate flows past immersed complicated bodies. With the Chapman-Enskog expansion analysis, a transform is performed between the Navier-Stokes and lattice Boltzmann equations (LBEs). The LBFS is used to discretize the macroscopic differential equations with a finite volume method and evaluate the interface fluxes through local reconstruction of the lattice Boltzmann solution. The immersed boundary technique is used to correct the intermediate velocity around the solid boundary to satisfy the no-slip boundary condition. Agreement of simulation results with the data found in the literature shows reliability of the proposed method in simulating laminar flows on a Cartesian mesh.
基金supported by the Natural Science Foundation of China (11061021)the Program of Higher-level talents of Inner Mongolia University (SPH-IMU,Z200901004)the Scientific Research Projection of Higher Schools of Inner Mongolia(NJ10016,NJ10006)
文摘This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional (2D) incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling. The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by join- ing the centroid of cells sharing the common vertex. For the temporal integration of the momentum equations, an im- plicit second-order scheme is utilized to enhance the com- putational stability and eliminate the time step limit due to the diffusion term. The momentum equations are discretized by the vertex-centered finite volume method (FVM) and the pressure Poisson equation is solved by the Galerkin finite el- ement method (FEM). The momentum interpolation is used to damp out the spurious pressure wiggles. The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both veloc- ity and pressure. The classic test cases, the lid-driven cavity flow, the skew cavity flow and the backward-facing step flow, show that numerical results are in good agreement with the published benchmark solutions.
基金Supported by the National Natural Science Foundation of China (50606026, 50736006).
文摘The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation (DNS) of incompressible turbulent reacting channel flows. The instantaneous DNS results show the near-wall strip structures of concentration and temperature fluctuations. The DNS statistical results give the budget of the terms in the correlation equations, showing that the production and dissipation terms are most important. The DNS statistical data are used to validate the closure model in RANS second-order moment (SOM) combustion model. It is found that the simulated diffusion and production terms are in agreement with the DNS data in most flow regions, except in the near-wall region, where the near-wall modification should be made, and the closure model for the dissipation term needs further improvement. The algebraic second-order moment (ASOM) combustion model is well validated by DNS.
文摘We present a cut-cell method for the simulation of 2D incompressible flows past obstacles.It consists in using the MAC scheme on cartesian grids and imposing Dirchlet boundary conditions for the velocity field on the boundary of solid structures following the Shortley-Weller formulation.In order to ensure local conservation properties,viscous and convecting terms are discretized in a finite volume way.The scheme is second order implicit in time for the linear part,the linear systems are solved by the use of the capacitance matrix method for non-moving obstacles.Numerical results of flows around an impulsively started circular cylinder are presented which confirm the efficiency of the method,for Reynolds numbers 1000 and 3000.An example of flows around a moving rigid body at Reynolds number 800 is also shown,a solver using the PETSc-Library has been prefered in this context to solve the linear systems.
基金Project supported by the National Natural Science Foundation of China(No.50576049) the Foun-dational Scientific Research of National Defence of China(No.A4020060263)Shanghai Leading Academic Discipline Project(No.Y0103)
文摘A unified numerical scheme for the solutions of the compressible and incompressible Navier-Stokes equations is investigated based on a time-derivative preconditioning algorithm. The primitive variables are pressure, velocities and temperature. The time integration scheme is used in conjunction with a finite volume discretization. The preconditioning is coupled with a high order implicit upwind scheme based on the definition of a Roe's type matrix. Computational capabilities are demonstrated through computations of high Mach number, middle Mach number, very low Mach number, and incompressible flow. It has also been demonstrated that the discontinuous surface in flow field can be captured for the implementation Roe's scheme.
文摘Three exact solutions are obtained for 2-D incompressible potential flows around two moving circles in three cases: (i) expansion (or contraction) of themselves, (ii) approaching (or departing from) each other, (iii) moving perpendicularly to the line connecting the centres in opposite directions. Meanwhile, an- other set of two exact solutions is obtained for 2-D incompressible potential flows between two moving eccen- tric circles in two cases: moving parallelly or perpendicularly to the line connecting the centres.
文摘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.
文摘A simple method is proposed, for incremental static analysis of a set of inter-colliding particles, simulating 2D flow. Within each step of proposed algorithm, the particles perform small displacements, proportional to the out-of-balance forces, acting on them. Numerical experiments show that if the liquid is confined within boundaries of a set of inter-communicating vessels, then the proposed method converges to a final equilibrium state. This incremental static analysis approximates dynamic behavior with strong damping and can provide information, as a first approximation to 2D movement of a liquid. In the initial arrangement of particles, a rhombic element is proposed, which assures satisfactory incompressibility of the fluid. Based on the proposed algorithm, a simple and short computer program (a “pocket” program) has been developed, with only about 120 Fortran instructions. This program is first applied to an amount of liquid, contained in a single vessel. A coarse and refined discretization is tried. In final equilibrium state of liquid, the distribution on hydro-static pressure on vessel boundaries, obtained by proposed computational model, is found in satisfactory approximation with corresponding theoretical data. Then, an opening is formed, at the bottom of a vertical boundary of initial vessel, and the liquid is allowed to flow gradually to an adjacent vessel. Almost whole amount of liquid is transferred, from first to second vessel, except of few drops-particles, which remain, in equilibrium, at the bottom of initial vessel. In the final equilibrium state of liquid, in the second vessel, the free surface level of the liquid confirms that the proposed rhombing element assures a satisfactory incompressibility of the fluid.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10972010,11028206,11371069,11372052,11402029,and 11472060)the Science and Technology Development Foundation of China Academy of Engineering Physics(CAEP),China(Grant No.2014B0201030)the Defense Industrial Technology Development Program of China(Grant No.B1520132012)
文摘Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method(GAHM) is presented for incompressible two-phase flow.A front tracking method(FTM) is used as the skeleton of the GAHM for low mass loss and resources.Smooth eulerian level set values are calculated from the FTM interface,and are used for a local interface reconstruction.The reconstruction avoids marker particle redistribution and enables an automatic treatment of interfacial topology change.The cubic Hermit interpolation is employed in all steps of the GAHM to capture subgrid structures within a single spacial cell.The performance of the GAHM is carefully evaluated in a benchmark test.Results show significant improvements of mass loss,clear subgrid structures,highly accurate derivatives(normals and curvatures) and low cost.The GAHM is further coupled with an incompressible multiphase flow solver,Super CE/SE,for more complex and practical applications.The updated solver is evaluated through comparison with an early droplet research.
文摘Generally the incompressible viscous flow problem is described by the Navier-Stokes equation. Based on the weighted residual method the discrete formulation of element-free Galerkin is inferred in this paper. By the step-bystep computation in the field of time, and adopting the least-square estimation of the-same-order shift, this paper has calculated both velocity and pressure from the decoupling independent equations. Each time fraction Newton-Raphson iterative method is applied for the velocity and pressure. Finally, this paper puts the method into practice of the shear-drive cavity flow, verifying the validity, high accuracy and stability.
文摘A new scalar projection method presented for simulating incompressible flows with variable density is proposed. It reverses conventional projection algorithm by computing first the irrotational component of the velocity and then the pressure. The first phase of the projection is purely kinematics. The predicted velocity field is subjected to a discrete Hodge-Helmholtz decomposition. The second phase of upgrade of pressure from the density uses Stokes’ theorem to explicitly compute the pressure. If all or part of the boundary conditions is then fixed on the divergence free physical field, the system required to be solved for the scalar potential of velocity becomes a Poisson equation with constant coefficients fitted with Dirichlet conditions.