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.展开更多
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.展开更多
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.展开更多
A finite volume method based unstructured grid is presented to solve the two dimensional viscous and incompressible flow.The method is based on the pressure-correction concept and solved by using a semi-staggered grid...A finite volume method based unstructured grid is presented to solve the two dimensional viscous and incompressible flow.The method is based on the pressure-correction concept and solved by using a semi-staggered grid technique.The computational procedure can handle cells of arbitrary shapes,although solutions presented in this paper were only involved with triangular and quadrilateral cells.The pressure or pressure-correction value was stored on the vertex of cells.The mass conservation equation was discretized on the dual cells surrounding the vertex of primary cells,while the velocity components and other scale variables were saved on the central of primary cells.Since the semi-staggered arrangement can’t guarantee a strong coupling relationship between pressure and velocity,thus a weak coupling relationship leads to the oscillations for pressure and velocity.In order to eliminate such an oscillation,a special interpolation scheme was used to construct the pressure-correction equation.Computational results of several viscous flow problems show good agreement with the analytical or numerical results in previous literature.This semi-staggered grid method can be applied to arbitrary shape elements,while it has the most efficiency for triangular cells.展开更多
In order to make the numerical calculation of viscous flows more convenient for the flows in channel with complicated profile governing equations expressed in the arbitrary curvilinear coordinates were derived by mean...In order to make the numerical calculation of viscous flows more convenient for the flows in channel with complicated profile governing equations expressed in the arbitrary curvilinear coordinates were derived by means of Favre density-weighted averaged method, and a turbulent model with effect of curvature modification was also derived. The numerical calculation of laminar and turbulent flown in divergent curved channels was carried out by means of parabolizeil computation method. The calculating results were used to analyze and investigate the aerodynamic performance of talor cascades in compressors preliminarily.展开更多
The possible states in the flow past two identical cylinders in tandem arrangements are investigated. The effect of the gap （L/D = 1.5, 1.75 and 2） between the two cylinders at Reynolds number （Re = 52,639） is tak...The possible states in the flow past two identical cylinders in tandem arrangements are investigated. The effect of the gap （L/D = 1.5, 1.75 and 2） between the two cylinders at Reynolds number （Re = 52,639） is taken into consideration. The presence of three-dimensional flow structures was observed to include notable changes to the response of the flow as result of variation of cylinder separation. A number of planes （z/h = 0.02, 0.25, 0.5 and 0.98） were taken at 20 step times of interval 0.005 s. to cover the details of flow along the cylinders. CFD FLUENT program was used to detect the flow structure. It is observed that the gap between the two cylinders affects the flow regime, i.e., there is no distinct vortex shedding downstream of the first cylinder.展开更多
This paper mainly studies the well-posedness of steady incompressible impinging jet flow problem through a 3D axisymmetric finitely long nozzle.This problem originates from the physical phenomena encountered in practi...This paper mainly studies the well-posedness of steady incompressible impinging jet flow problem through a 3D axisymmetric finitely long nozzle.This problem originates from the physical phenomena encountered in practical engineering fields,such as in short take-off and vertical landing(STOVL)aircraft.Nowadays many intricate phenomena associated with impinging jet flows remain inadequately elucidated,which limits the ability to optimize aircraft design.Given a boundary condition in the inlet,the impinging jet problem is transformed into a Bernoulli-type free boundary problem according to the stream function.Then the variational method is used to study the corresponding variational problem with one parameter,thereby the wellposedness is established.The main conclusion is as follows.For a 3D axisymmetric finitely long nozzle and an infinitely long vertical wall,given an axial velocity in the inlet of nozzle,there exists a unique smooth incom‑pressible impinging jet flow such that the free boundary initiates smoothly at the endpoint of the nozzle and extends to infinity along the vertical wall at far fields.The key point is to investigate the regularity of the corner where the nozzle and the vertical axis intersect.展开更多
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.展开更多
Smoothed particle hydrodynamics （SPH） is a Lagrangian, meshfree particle method and has been widely applied to diffe- rent areas in engineering and science. Since its original extension to modeling free surface flow...Smoothed particle hydrodynamics （SPH） is a Lagrangian, meshfree particle method and has been widely applied to diffe- rent areas in engineering and science. Since its original extension to modeling free surface flows by Monaghan in 1994, SPH has been gradually developed into an attractive approach for modeling viscous incompressible fluid flows. This paper presents an overview on the recent progresses of SPH in modeling viscous incompressible flows in four major aspects which are closely related to the computational accuracy of SPH simulations. The advantages and disadvantages of different SPH particle approximation sche- mes, pressure field solution approaches, solid boundary treatment algorithms and particle adapting algorithms are described and analyzed. Some new perspectives and fuRtre trends in SPH modeling of viscous incompressible flows are discussed.展开更多
A lattice Boltzmann flux solver(LBFS)is presented in this work for simulation of incompressible viscous and inviscid flows.The new solver is based on Chapman-Enskog expansion analysis,which is the bridge to link Navie...A lattice Boltzmann flux solver(LBFS)is presented in this work for simulation of incompressible viscous and inviscid flows.The new solver is based on Chapman-Enskog expansion analysis,which is the bridge to link Navier-Stokes(N-S)equations and lattice Boltzmann equation(LBE).The macroscopic differential equations are discretized by the finite volume method,where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers.The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.LBFS is validated by its application to simulate the viscous decaying vortex flow,the driven cavity flow,the viscous flow past a circular cylinder,and the inviscid flow past a circular cylinder.The obtained numerical results compare very well with available data in the literature,which show that LBFS has the second order of accuracy in space,and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.展开更多
A fast, matrix-free implicit method based on the Lower-Upper Symmetric Gauss-Seidel （LU-SGS） method was developed to solve the three-dimensional incompressible Reynolds-averaged Navier-Stokes equations on multi-bloc...A fast, matrix-free implicit method based on the Lower-Upper Symmetric Gauss-Seidel （LU-SGS） method was developed to solve the three-dimensional incompressible Reynolds-averaged Navier-Stokes equations on multi-block structured grids. The method was applied to the simulations of a variety of flowfields around 3D complex underwater bodies with different appendages. The numerical procedure and results show that the method is efficient, reliable and robust for steady viscous flow simulations.展开更多
In this study,a three-dimensional artificial compressibility solver based on the average-state Harten-Lax-van Leer-Contact(HLLC)[13]type Riemann solution is first proposed and developed to solve the time-dependent inc...In this study,a three-dimensional artificial compressibility solver based on the average-state Harten-Lax-van Leer-Contact(HLLC)[13]type Riemann solution is first proposed and developed to solve the time-dependent incompressible flow equations.To implement unsteady flow calculations,a dual time stepping strategy including the LU decomposition method is used in the pseudo-time iteration and the second-order accurate backward difference is adopted to discretize the unsteady flow term.Also a third-order accurate HLLC numerical flux is derived for approximating the inviscid terms.To verify numerical accuracy,flows over a lid-driven cavity and an oscillating flat plate are chosen as the benchmark tests.In addition,the current solver is extended to solve blood flows in a realistic human aorta measured from MRI(Magnetic Resonance Imaging).The simulation geometry was derived from a three-dimensional reconstruction of a series of two-dimensional slices obtained in vivo.Numerical results demonstrate wall stresses were highly dynamic,but were generally high along the outer wall in the vicinity of the branches and low along the inner wall,particularly in the descending aorta.The maximum wall stress distribution is presented on the aortic arch in the systole.In addition,extensive counter-clockwise secondary flows and three-dimensional helical vortex influenced considerably by the presence of vessel contraction,torsion and the branches were shown in the descending aorta in the late systole and early diastolic cycles.展开更多
The authors study the Rayleigh-Taylor instability for two incompressible immiscible fluids with or without surface tension, evolving with a free interface in the presence of a uniform gravitational field in Eulerian c...The authors study the Rayleigh-Taylor instability for two incompressible immiscible fluids with or without surface tension, evolving with a free interface in the presence of a uniform gravitational field in Eulerian coordinates. To deal with the free surface, instead of using the transformation to Lagrangian coordinates, the perturbed equations in Eulerian coordinates are transformed to an integral form and the two-fluid flow is formulated as a single-fluid flow in a fixed domain, thus offering an alternative approach to deal with the jump conditions at the free interface. First, the linearized problem around the steady state which describes a denser immiscible fluid lying above a light one with a free interface separating the two fluids, both fluids being in(unstable) equilibrium is analyzed. By a general method of studying a family of modes, the smooth(when restricted to each fluid domain) solutions to the linearized problem that grow exponentially fast in time in Sobolev spaces are constructed, thus leading to a global instability result for the linearized problem.Then, by using these pathological solutions, the global instability for the corresponding nonlinear problem in an appropriate sense is demonstrated.展开更多
In this article we discuss a new full discrete scheme for the numerical solution of the Navier-Stokes equations modeling viscous incompressible flow. This scheme consists of nonlinear Galerkin method using mixed fini...In this article we discuss a new full discrete scheme for the numerical solution of the Navier-Stokes equations modeling viscous incompressible flow. This scheme consists of nonlinear Galerkin method using mixed finite elements and Crank-Nicolson method. Next, we provide the second-order convergence accuracy of numerical solution corresponding to this scheme. Compared with the usual Galerkin scheme, this scheme can save a large amount of computational time under the same convergence accuracy. (Author abstract) 8 Refs.展开更多
Fixed-interval smoothing,as one of the most important types of state estimation,has been concerned in many practical problems especially in the analysis of flight test data.However,the existing sequential filters and ...Fixed-interval smoothing,as one of the most important types of state estimation,has been concerned in many practical problems especially in the analysis of flight test data.However,the existing sequential filters and smoothers usually cannot deal with nonlinear or high-dimensional systems well.A state-of-the-art technique is employed in this study to explore the fixed-interval smoothing problem of a conceptual two-dimensional airfoil model in incompressible flow from noisy measurement data.Therein,the governing equations of the airfoil model are assumed to be known or only partially known.A single objective optimization problem is constructed with the classical Runge–Kutta scheme,and then estimations of the system states,the measurement noise and even the unknown parameters are obtained simultaneously through minimizing the objective function.Effectiveness and feasibility of the method are examined under several simulated measurement data corrupted by different measurement noises.All the obtained results indicate that the introduced algorithm is applicable for the airfoil model with cubic or free-play structural nonlinearity and leads to accurate state and parameter estimations.Besides,it is highly robust to Gaussian white and even more complex heavy-tailed measurement noises.It should be emphasized that the employed algorithm is still effective to high-dimensional nonlinear aeroelastic systems.展开更多
This paper focuses on the low-order nonconforming rectangular and quadrilateral finite elements approximation of incompressible flow.Beyond the previous research works,we propose a general strategy to construct the ba...This paper focuses on the low-order nonconforming rectangular and quadrilateral finite elements approximation of incompressible flow.Beyond the previous research works,we propose a general strategy to construct the basis functions.Under several specific constraints,the optimal error estimates are obtained,i.e.,the first order accuracy of the velocities in H1-norm and the pressure in L2-norm,as well as the second order accuracy of the velocities in L2-norm.Besides,we clarify the differences between rectangular and quadrilateral finite element approximation.In addition,we give several examples to verify the validity of our error estimates.展开更多
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展开更多
Fourth-order stream-function methods are proposed for the time dependent, incom- pressible Navier-Stokes and Boussinesq equations. Wide difference stencils are used instead of compact ones and the boundary terms are h...Fourth-order stream-function methods are proposed for the time dependent, incom- pressible Navier-Stokes and Boussinesq equations. Wide difference stencils are used instead of compact ones and the boundary terms are handled by extrapolating the stream-function values inside the computational domain to grid points outside, up to fourth-order in the noslip condition. Formal error analysis is done for a simple model problem, showing that this extrapolation introduces numerical boundary layers at fifth-order in the stream-function. The fourth-order convergence in velocity of the proposed method for the full problem is shown numerically.展开更多
文摘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.
文摘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.
文摘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 Natural Science Foundation of China(No.51676208)the Fundamental Research Funds for the Central Universities(No.18CX07012A and No.19CX05002A)support from the Major Program of the Natural Science Foundation of Shandong Province(No.ZR2019ZD11).
文摘A finite volume method based unstructured grid is presented to solve the two dimensional viscous and incompressible flow.The method is based on the pressure-correction concept and solved by using a semi-staggered grid technique.The computational procedure can handle cells of arbitrary shapes,although solutions presented in this paper were only involved with triangular and quadrilateral cells.The pressure or pressure-correction value was stored on the vertex of cells.The mass conservation equation was discretized on the dual cells surrounding the vertex of primary cells,while the velocity components and other scale variables were saved on the central of primary cells.Since the semi-staggered arrangement can’t guarantee a strong coupling relationship between pressure and velocity,thus a weak coupling relationship leads to the oscillations for pressure and velocity.In order to eliminate such an oscillation,a special interpolation scheme was used to construct the pressure-correction equation.Computational results of several viscous flow problems show good agreement with the analytical or numerical results in previous literature.This semi-staggered grid method can be applied to arbitrary shape elements,while it has the most efficiency for triangular cells.
文摘In order to make the numerical calculation of viscous flows more convenient for the flows in channel with complicated profile governing equations expressed in the arbitrary curvilinear coordinates were derived by means of Favre density-weighted averaged method, and a turbulent model with effect of curvature modification was also derived. The numerical calculation of laminar and turbulent flown in divergent curved channels was carried out by means of parabolizeil computation method. The calculating results were used to analyze and investigate the aerodynamic performance of talor cascades in compressors preliminarily.
文摘The possible states in the flow past two identical cylinders in tandem arrangements are investigated. The effect of the gap （L/D = 1.5, 1.75 and 2） between the two cylinders at Reynolds number （Re = 52,639） is taken into consideration. The presence of three-dimensional flow structures was observed to include notable changes to the response of the flow as result of variation of cylinder separation. A number of planes （z/h = 0.02, 0.25, 0.5 and 0.98） were taken at 20 step times of interval 0.005 s. to cover the details of flow along the cylinders. CFD FLUENT program was used to detect the flow structure. It is observed that the gap between the two cylinders affects the flow regime, i.e., there is no distinct vortex shedding downstream of the first cylinder.
文摘This paper mainly studies the well-posedness of steady incompressible impinging jet flow problem through a 3D axisymmetric finitely long nozzle.This problem originates from the physical phenomena encountered in practical engineering fields,such as in short take-off and vertical landing(STOVL)aircraft.Nowadays many intricate phenomena associated with impinging jet flows remain inadequately elucidated,which limits the ability to optimize aircraft design.Given a boundary condition in the inlet,the impinging jet problem is transformed into a Bernoulli-type free boundary problem according to the stream function.Then the variational method is used to study the corresponding variational problem with one parameter,thereby the wellposedness is established.The main conclusion is as follows.For a 3D axisymmetric finitely long nozzle and an infinitely long vertical wall,given an axial velocity in the inlet of nozzle,there exists a unique smooth incom‑pressible impinging jet flow such that the free boundary initiates smoothly at the endpoint of the nozzle and extends to infinity along the vertical wall at far fields.The key point is to investigate the regularity of the corner where the nozzle and the vertical axis intersect.
基金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.
基金Project supported by the National Natural Science Foun-dation of China(Grant Nos.11172306,U1530110)the Institu-te of Systems Engineering,China Academy of Engineering Physics(Grant No.2013KJZ01)
文摘Smoothed particle hydrodynamics （SPH） is a Lagrangian, meshfree particle method and has been widely applied to diffe- rent areas in engineering and science. Since its original extension to modeling free surface flows by Monaghan in 1994, SPH has been gradually developed into an attractive approach for modeling viscous incompressible fluid flows. This paper presents an overview on the recent progresses of SPH in modeling viscous incompressible flows in four major aspects which are closely related to the computational accuracy of SPH simulations. The advantages and disadvantages of different SPH particle approximation sche- mes, pressure field solution approaches, solid boundary treatment algorithms and particle adapting algorithms are described and analyzed. Some new perspectives and fuRtre trends in SPH modeling of viscous incompressible flows are discussed.
文摘A lattice Boltzmann flux solver(LBFS)is presented in this work for simulation of incompressible viscous and inviscid flows.The new solver is based on Chapman-Enskog expansion analysis,which is the bridge to link Navier-Stokes(N-S)equations and lattice Boltzmann equation(LBE).The macroscopic differential equations are discretized by the finite volume method,where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers.The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.LBFS is validated by its application to simulate the viscous decaying vortex flow,the driven cavity flow,the viscous flow past a circular cylinder,and the inviscid flow past a circular cylinder.The obtained numerical results compare very well with available data in the literature,which show that LBFS has the second order of accuracy in space,and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.
基金Project supported by the National Postdoctor Foundation (Grant No: 2003033308).
文摘A fast, matrix-free implicit method based on the Lower-Upper Symmetric Gauss-Seidel （LU-SGS） method was developed to solve the three-dimensional incompressible Reynolds-averaged Navier-Stokes equations on multi-block structured grids. The method was applied to the simulations of a variety of flowfields around 3D complex underwater bodies with different appendages. The numerical procedure and results show that the method is efficient, reliable and robust for steady viscous flow simulations.
文摘In this study,a three-dimensional artificial compressibility solver based on the average-state Harten-Lax-van Leer-Contact(HLLC)[13]type Riemann solution is first proposed and developed to solve the time-dependent incompressible flow equations.To implement unsteady flow calculations,a dual time stepping strategy including the LU decomposition method is used in the pseudo-time iteration and the second-order accurate backward difference is adopted to discretize the unsteady flow term.Also a third-order accurate HLLC numerical flux is derived for approximating the inviscid terms.To verify numerical accuracy,flows over a lid-driven cavity and an oscillating flat plate are chosen as the benchmark tests.In addition,the current solver is extended to solve blood flows in a realistic human aorta measured from MRI(Magnetic Resonance Imaging).The simulation geometry was derived from a three-dimensional reconstruction of a series of two-dimensional slices obtained in vivo.Numerical results demonstrate wall stresses were highly dynamic,but were generally high along the outer wall in the vicinity of the branches and low along the inner wall,particularly in the descending aorta.The maximum wall stress distribution is presented on the aortic arch in the systole.In addition,extensive counter-clockwise secondary flows and three-dimensional helical vortex influenced considerably by the presence of vessel contraction,torsion and the branches were shown in the descending aorta in the late systole and early diastolic cycles.
基金supported by the National Natural Science Foundation of China(Nos.11101044,11271051,11229101,11301083,11371065,11471134)the Fujian Provincial Natural Science Foundation of China(No.2014J01011)+1 种基金the National Basic Research Program(No.2011CB309705)the Beijing Center for Mathematics and Information Interdisciplinary Sciences
文摘The authors study the Rayleigh-Taylor instability for two incompressible immiscible fluids with or without surface tension, evolving with a free interface in the presence of a uniform gravitational field in Eulerian coordinates. To deal with the free surface, instead of using the transformation to Lagrangian coordinates, the perturbed equations in Eulerian coordinates are transformed to an integral form and the two-fluid flow is formulated as a single-fluid flow in a fixed domain, thus offering an alternative approach to deal with the jump conditions at the free interface. First, the linearized problem around the steady state which describes a denser immiscible fluid lying above a light one with a free interface separating the two fluids, both fluids being in(unstable) equilibrium is analyzed. By a general method of studying a family of modes, the smooth(when restricted to each fluid domain) solutions to the linearized problem that grow exponentially fast in time in Sobolev spaces are constructed, thus leading to a global instability result for the linearized problem.Then, by using these pathological solutions, the global instability for the corresponding nonlinear problem in an appropriate sense is demonstrated.
文摘In this article we discuss a new full discrete scheme for the numerical solution of the Navier-Stokes equations modeling viscous incompressible flow. This scheme consists of nonlinear Galerkin method using mixed finite elements and Crank-Nicolson method. Next, we provide the second-order convergence accuracy of numerical solution corresponding to this scheme. Compared with the usual Galerkin scheme, this scheme can save a large amount of computational time under the same convergence accuracy. (Author abstract) 8 Refs.
基金supported by the National Natural Sciencs Fundation of China(Grants 12072264.11772255)the Fundamental Research Funds for the Central Universities,the National Key Research and Development Program of China(Grant 2018AAA0102201)+2 种基金the Research Funds for Interdisciplinary Subject of Northwestern Polytechnical University,the Shaanxi Project for Distinguished Young Scholars,the Shaanxi Provincial Key R&D Program(Grants 2O2OKW-013.2019TD-010)the Innovation Foundation for Doctor Dissertation of Northwestern Polytechnical University(Grant CX201962)J.K.was sup ported by the Russian Ministry of Science and Education Agreement(Grant 075-15-2020-808).
文摘Fixed-interval smoothing,as one of the most important types of state estimation,has been concerned in many practical problems especially in the analysis of flight test data.However,the existing sequential filters and smoothers usually cannot deal with nonlinear or high-dimensional systems well.A state-of-the-art technique is employed in this study to explore the fixed-interval smoothing problem of a conceptual two-dimensional airfoil model in incompressible flow from noisy measurement data.Therein,the governing equations of the airfoil model are assumed to be known or only partially known.A single objective optimization problem is constructed with the classical Runge–Kutta scheme,and then estimations of the system states,the measurement noise and even the unknown parameters are obtained simultaneously through minimizing the objective function.Effectiveness and feasibility of the method are examined under several simulated measurement data corrupted by different measurement noises.All the obtained results indicate that the introduced algorithm is applicable for the airfoil model with cubic or free-play structural nonlinearity and leads to accurate state and parameter estimations.Besides,it is highly robust to Gaussian white and even more complex heavy-tailed measurement noises.It should be emphasized that the employed algorithm is still effective to high-dimensional nonlinear aeroelastic systems.
基金supported by National Natural Science Foundation of China(GrantNo.11071139)National Basic Research Program of China(Grant No.2011CB309705)Tsinghua University Initiative Scientific Research Program
文摘This paper focuses on the low-order nonconforming rectangular and quadrilateral finite elements approximation of incompressible flow.Beyond the previous research works,we propose a general strategy to construct the basis functions.Under several specific constraints,the optimal error estimates are obtained,i.e.,the first order accuracy of the velocities in H1-norm and the pressure in L2-norm,as well as the second order accuracy of the velocities in L2-norm.Besides,we clarify the differences between rectangular and quadrilateral finite element approximation.In addition,we give several examples to verify the validity of our error estimates.
文摘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
基金funding from NSF under grants DMS-0713670 and ACI-0204932funding from NSERC Canada that supported this work
文摘Fourth-order stream-function methods are proposed for the time dependent, incom- pressible Navier-Stokes and Boussinesq equations. Wide difference stencils are used instead of compact ones and the boundary terms are handled by extrapolating the stream-function values inside the computational domain to grid points outside, up to fourth-order in the noslip condition. Formal error analysis is done for a simple model problem, showing that this extrapolation introduces numerical boundary layers at fifth-order in the stream-function. The fourth-order convergence in velocity of the proposed method for the full problem is shown numerically.
