The operator splitting method is used to deal with the Navier-Stokes equation, in which the physical process described by the equation is decomposed into two processes： a diffusion process and a convection process; a...The operator splitting method is used to deal with the Navier-Stokes equation, in which the physical process described by the equation is decomposed into two processes： a diffusion process and a convection process; and the finite element equation is established. The velocity field in the element is described by the shape function of the isoparametric element with nine nodes and the pressure field is described by the interpolation function of the four nodes at the vertex of the isoparametric element with nine nodes. The subroutine of the element and the integrated finite element code are generated by the Finite Element Program Generator （FEPG） successfully. The numerical simulation about the incompressible viscous liquid flowing over a cylinder is carded out. The solution agrees with the experimental results very well.展开更多
Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton...Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton and the relaxation mechanisms of porosity and saturation（capillary pressure）.So it accurately simulates the numerical attenuation property of the wavefields and is much closer to actual earth media in exploration than the equivalent liquid model and the unsaturated porous medium model on the basis of open system theory.The velocity and attenuation for different wave modes in this medium have been discussed in previous literature but studies of the complete wave-field have not been reported.In our work,wave equations with the relaxation mechanisms of capillary pressure and the porosity are derived.Furthermore,the wavefield and its characteristics are studied using the numerical finite element method.The results show that the slow P3-wave in the non-wetting phase can be observed clearly in the seismic band.The relaxation of capillary pressure and the porosity greatly affect the displacement of the non-wetting phase.More specifically,the displacement decreases with increasing relaxation coefficient.展开更多
In practical engineering,only pressure sensors are allowed to install to detect leakage in most of oil transportation pipelines,while flowmeters are only installed at the toll ports.For incompressible fluid,the leakag...In practical engineering,only pressure sensors are allowed to install to detect leakage in most of oil transportation pipelines,while flowmeters are only installed at the toll ports.For incompressible fluid,the leakage rate and amount cannot be accurately calculated through critical pressure conditions.In this paper,a micro-element body of the pipeline was intercepted for calculation.The relationship between radial displacement and pressure of pipe wall was studied based on the stress-strain equation.Then,the strain response of pipeline volume with pipeline pressure was obtained.The change in volume expansion of pipeline was used to characterize leakage of incompressible fluid.Finally,the calculation model of leakage amount of incompressible fluid was obtained.To verify the above theory,the pipeline expansion model under pressure was established by COMSOL software for simulation.Both simulation results and deduction equations show that the volumetric change has a quadratic parabolic relationship with the change of pipeline pressure.However,the relationship between them can be approximately linear when the pressure change is not too large.In addition,the leakage of incompressible fluid under the pressure of 0 MPa-0.8 MPa was obtained by experiments.The experimental results verify the linear relationship between leakage of incompressible fluid and the change of pipeline pressure.The theoretical and experimental results provide a basis for the calculation of leakage of incompressible fluid in the pipeline.展开更多
In this paper, an efficient multigrid fictitious boundary method （MFBM） coupled with the FEM solver package FEATFLOW was used for the detailed simulation of incompressible viscous flows around one or more moving NAC...In this paper, an efficient multigrid fictitious boundary method （MFBM） coupled with the FEM solver package FEATFLOW was used for the detailed simulation of incompressible viscous flows around one or more moving NACA0012 airfoils. The calculations were carded on a fixed multigrid finite element mesh on which fluid equations were satisfied everywhere, and the airfoils were allowed to move freely through the mesh. The MFBM was employed to treat interactions between the fluid and the airfoils The motion of the airfoils was modeled by Newton-Euler equations. Numerical results of experiments verify that this method provides an efficient way to simulate incompressible viscous flows around moving airfoils.展开更多
The pinch instability for a cylindrical jet of liquid metal passed through by an axial electrical current is investigated. Besides the pinch effect originating from surface tension, the Lorentz force, created by the a...The pinch instability for a cylindrical jet of liquid metal passed through by an axial electrical current is investigated. Besides the pinch effect originating from surface tension, the Lorentz force, created by the axial current density and the corresponding azimuthal magnetic field, causes an electromagnetic pinch effect. This effect has drawn attention in electrical engineering, because it can be used in the construction of liquid metal current limit- ers with self-healing properties. In this paper a simple model is derived using the shallow water approximation： the equations describing the full system are reduced to two one-dimensional evolution equations for the axial velocity and the radius of the jet. A stability analysis for this reduced system is carried out yielding critical current density and the growth rate for the instability. To investigate the nonlinear behaviour of the pinch instability for finite perturbations simulations, the shallow water model are performed.展开更多
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.展开更多
Flow around two rotating side-by-side circular cylinders of equal diameter D is numerically studied at the Rey- nolds number 40〈 Re 〈200 and various rotation rate 8i. The incoming flow is assumed to be two-dimension...Flow around two rotating side-by-side circular cylinders of equal diameter D is numerically studied at the Rey- nolds number 40〈 Re 〈200 and various rotation rate 8i. The incoming flow is assumed to be two-dimensional laminar flow. The governing equations are the incompressible Navier-Stokes equations and solved by the finite volume method （FVM）. The ratio of the center-to-center spacing to the cylinder diameter is T/D=2. The objective of the present work is to investigate the effect of rotational speed and Reynolds number on the stability of the flow. The simulation results are compared with the experimental data and a good agreement is achieved. The stability of the flow is analyzed by using the energy gradient theory, which produces the energy gradient function K to identify the region where the flow is the most prone to be destabilized and the degree of the destabilization. Numerical results reveal that K is the most significant at the separated shear layers of the cylinder pair. With Re in- creases, the length of the wake is shorter and the vortex shedding generally exhibits a symmetrical distribution for θi〈θcrit. It is also shown that the unsteady vortex shedding can be suppressed by rotating the cylinders in the counter-rotating mode.展开更多
In this paper,a two-level search method for searching transfer opportunities between interplanetary halo orbits,exploiting the invariant manifolds of the restricted three-body problem,is proposed.In the method,the fir...In this paper,a two-level search method for searching transfer opportunities between interplanetary halo orbits,exploiting the invariant manifolds of the restricted three-body problem,is proposed.In the method,the first-level search procedure is performed under the conditions of the initial time of escape manifold trajectory of the Sun-Earth halo orbit and the terminal time of capture manifold of the target planet fixed,by solving the optimal two-impulsive heliocentric trajectory to connect the two manifold trajectories.The contour map,helpful to the understanding of the global characteristics of the transfer opportunities,taking the initial time of escape manifold and the terminal time of capture manifold as variables,the optimal velocity increment of the first-level search as objective function,is used for the second-level search.Finally,taking the Earth-Mars and Earth-Venus halo to halo transfers for example,the transfer opportunities in 2015-2017 are searched.The results show the effectiveness of the proposed method and reveal the property of quasi-period of transfer opportunities between interplanetary halo orbits.展开更多
A numerical algorithm using a bilinear or linear finite element and semi-implicit three-step method is presented for the analysis of incompressible viscous fluid problems. The streamline upwind/Petrov-Galerkin (SUPG) ...A numerical algorithm using a bilinear or linear finite element and semi-implicit three-step method is presented for the analysis of incompressible viscous fluid problems. The streamline upwind/Petrov-Galerkin (SUPG) stabilization scheme is used for the formulation of the Navier-Stokes equations. For the spatial discretization, the convection term is treated explicitly, while the viscous term is treated implicitly, and for the temporal discretization, a three-step method is employed. The present method is applied to simulate the lid driven cavity problems with different geometries at low and high Reynolds numbers. The results compared with other numerical experiments are found to be feasible and satisfactory.展开更多
This article discusses computational methods for the numerical simulation of unsteady Bingham visco-plastic flow. These methods are based on time-discretization by operator-splitting and take advantage of a characteri...This article discusses computational methods for the numerical simulation of unsteady Bingham visco-plastic flow. These methods are based on time-discretization by operator-splitting and take advantage of a characterization of the solutions involving some kind of Lagrange multipliers. The full discretization is achieved by combining the above operator-splitting methods with finite element approximations, the advection being treated by a wave-like equation 'equivalent' formulation easier to implement than the method of characteristics or high order upwinding methods. The authors illustrate the methodology discussed in this article with the results of numerical experiments concerning the simulation of wall driven cavity Bingham flow in two dimensions.展开更多
This paper studies the incompressible limit and stability of global strong solutions to the threedimensional full compressible Navier-Stokes equations, where the initial data satisfy the "well-prepared" cond...This paper studies the incompressible limit and stability of global strong solutions to the threedimensional full compressible Navier-Stokes equations, where the initial data satisfy the "well-prepared" conditions and the velocity field and temperature enjoy the slip boundary condition and convective boundary condition, respectively. The uniform estimates with respect to both the Mach number ∈(0, ∈] and time t ∈ [0, ∞) are established by deriving a differential inequality with decay property, where ∈∈(0, 1] is a constant.As the Mach number vanishes, the global solution to full compressible Navier-Stokes equations converges to the one of isentropic incompressible Navier-Stokes equations in t ∈ [0, +∞). Moreover, we prove the exponentially asymptotic stability for the global solutions of both the compressible system and its limiting incompressible system.展开更多
文摘The operator splitting method is used to deal with the Navier-Stokes equation, in which the physical process described by the equation is decomposed into two processes： a diffusion process and a convection process; and the finite element equation is established. The velocity field in the element is described by the shape function of the isoparametric element with nine nodes and the pressure field is described by the interpolation function of the four nodes at the vertex of the isoparametric element with nine nodes. The subroutine of the element and the integrated finite element code are generated by the Finite Element Program Generator （FEPG） successfully. The numerical simulation about the incompressible viscous liquid flowing over a cylinder is carded out. The solution agrees with the experimental results very well.
基金supported by the 973 Program (Grant No.2007CB209505)the National Natural Science Foundation of China (Grant No.40674061,40704019)
文摘Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton and the relaxation mechanisms of porosity and saturation（capillary pressure）.So it accurately simulates the numerical attenuation property of the wavefields and is much closer to actual earth media in exploration than the equivalent liquid model and the unsaturated porous medium model on the basis of open system theory.The velocity and attenuation for different wave modes in this medium have been discussed in previous literature but studies of the complete wave-field have not been reported.In our work,wave equations with the relaxation mechanisms of capillary pressure and the porosity are derived.Furthermore,the wavefield and its characteristics are studied using the numerical finite element method.The results show that the slow P3-wave in the non-wetting phase can be observed clearly in the seismic band.The relaxation of capillary pressure and the porosity greatly affect the displacement of the non-wetting phase.More specifically,the displacement decreases with increasing relaxation coefficient.
文摘In practical engineering,only pressure sensors are allowed to install to detect leakage in most of oil transportation pipelines,while flowmeters are only installed at the toll ports.For incompressible fluid,the leakage rate and amount cannot be accurately calculated through critical pressure conditions.In this paper,a micro-element body of the pipeline was intercepted for calculation.The relationship between radial displacement and pressure of pipe wall was studied based on the stress-strain equation.Then,the strain response of pipeline volume with pipeline pressure was obtained.The change in volume expansion of pipeline was used to characterize leakage of incompressible fluid.Finally,the calculation model of leakage amount of incompressible fluid was obtained.To verify the above theory,the pipeline expansion model under pressure was established by COMSOL software for simulation.Both simulation results and deduction equations show that the volumetric change has a quadratic parabolic relationship with the change of pipeline pressure.However,the relationship between them can be approximately linear when the pressure change is not too large.In addition,the leakage of incompressible fluid under the pressure of 0 MPa-0.8 MPa was obtained by experiments.The experimental results verify the linear relationship between leakage of incompressible fluid and the change of pipeline pressure.The theoretical and experimental results provide a basis for the calculation of leakage of incompressible fluid in the pipeline.
基金Supported by National 863 Plan Project of Ministry of Science and Technology of China under Grant No. 2006AA09Z354National Natural Science Foundation of China under Grant No. 10672101.
文摘In this paper, an efficient multigrid fictitious boundary method （MFBM） coupled with the FEM solver package FEATFLOW was used for the detailed simulation of incompressible viscous flows around one or more moving NACA0012 airfoils. The calculations were carded on a fixed multigrid finite element mesh on which fluid equations were satisfied everywhere, and the airfoils were allowed to move freely through the mesh. The MFBM was employed to treat interactions between the fluid and the airfoils The motion of the airfoils was modeled by Newton-Euler equations. Numerical results of experiments verify that this method provides an efficient way to simulate incompressible viscous flows around moving airfoils.
基金the Deutsche Forschungsgemeinschaft in the French-German DFG-CNRS research program 'Numerische Strmungssimulation-Simulation Numérique d'Ecoulements'National Nataral Science Foundation of China under granted number 10772044
文摘The pinch instability for a cylindrical jet of liquid metal passed through by an axial electrical current is investigated. Besides the pinch effect originating from surface tension, the Lorentz force, created by the axial current density and the corresponding azimuthal magnetic field, causes an electromagnetic pinch effect. This effect has drawn attention in electrical engineering, because it can be used in the construction of liquid metal current limit- ers with self-healing properties. In this paper a simple model is derived using the shallow water approximation： the equations describing the full system are reduced to two one-dimensional evolution equations for the axial velocity and the radius of the jet. A stability analysis for this reduced system is carried out yielding critical current density and the growth rate for the instability. To investigate the nonlinear behaviour of the pinch instability for finite perturbations simulations, the shallow water model are performed.
基金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.
基金supported by National Natural Science Foundation of China(51579224)Zhejiang Province Key Science and Technology Innovation Team Project(2013TD18)Zhejiang Province Science and Technology Plan Project(2017C34007)
文摘Flow around two rotating side-by-side circular cylinders of equal diameter D is numerically studied at the Rey- nolds number 40〈 Re 〈200 and various rotation rate 8i. The incoming flow is assumed to be two-dimensional laminar flow. The governing equations are the incompressible Navier-Stokes equations and solved by the finite volume method （FVM）. The ratio of the center-to-center spacing to the cylinder diameter is T/D=2. The objective of the present work is to investigate the effect of rotational speed and Reynolds number on the stability of the flow. The simulation results are compared with the experimental data and a good agreement is achieved. The stability of the flow is analyzed by using the energy gradient theory, which produces the energy gradient function K to identify the region where the flow is the most prone to be destabilized and the degree of the destabilization. Numerical results reveal that K is the most significant at the separated shear layers of the cylinder pair. With Re in- creases, the length of the wake is shorter and the vortex shedding generally exhibits a symmetrical distribution for θi〈θcrit. It is also shown that the unsteady vortex shedding can be suppressed by rotating the cylinders in the counter-rotating mode.
基金supported by the National Basic Research Program of China ("973" Program) (Grant No. 2012CB720000)the National Natural Science Foundation of China (Grant Nos. 10832004 and 11102020)
文摘In this paper,a two-level search method for searching transfer opportunities between interplanetary halo orbits,exploiting the invariant manifolds of the restricted three-body problem,is proposed.In the method,the first-level search procedure is performed under the conditions of the initial time of escape manifold trajectory of the Sun-Earth halo orbit and the terminal time of capture manifold of the target planet fixed,by solving the optimal two-impulsive heliocentric trajectory to connect the two manifold trajectories.The contour map,helpful to the understanding of the global characteristics of the transfer opportunities,taking the initial time of escape manifold and the terminal time of capture manifold as variables,the optimal velocity increment of the first-level search as objective function,is used for the second-level search.Finally,taking the Earth-Mars and Earth-Venus halo to halo transfers for example,the transfer opportunities in 2015-2017 are searched.The results show the effectiveness of the proposed method and reveal the property of quasi-period of transfer opportunities between interplanetary halo orbits.
基金Project supported by the National Natural Science Foundation of China (No.51078230)the Research Fund for the Doctoral Program of Higher Education of China (No.200802480056)the Key Project of Fund of Science and Technology Development of Shanghai (No.10JC1407900),China
文摘A numerical algorithm using a bilinear or linear finite element and semi-implicit three-step method is presented for the analysis of incompressible viscous fluid problems. The streamline upwind/Petrov-Galerkin (SUPG) stabilization scheme is used for the formulation of the Navier-Stokes equations. For the spatial discretization, the convection term is treated explicitly, while the viscous term is treated implicitly, and for the temporal discretization, a three-step method is employed. The present method is applied to simulate the lid driven cavity problems with different geometries at low and high Reynolds numbers. The results compared with other numerical experiments are found to be feasible and satisfactory.
文摘This article discusses computational methods for the numerical simulation of unsteady Bingham visco-plastic flow. These methods are based on time-discretization by operator-splitting and take advantage of a characterization of the solutions involving some kind of Lagrange multipliers. The full discretization is achieved by combining the above operator-splitting methods with finite element approximations, the advection being treated by a wave-like equation 'equivalent' formulation easier to implement than the method of characteristics or high order upwinding methods. The authors illustrate the methodology discussed in this article with the results of numerical experiments concerning the simulation of wall driven cavity Bingham flow in two dimensions.
基金supported by National Natural Science Foundation of China (Grant No. 11471334)Program for New Century Excellent Talents in University (Grant No. NCET-12-0085)
文摘This paper studies the incompressible limit and stability of global strong solutions to the threedimensional full compressible Navier-Stokes equations, where the initial data satisfy the "well-prepared" conditions and the velocity field and temperature enjoy the slip boundary condition and convective boundary condition, respectively. The uniform estimates with respect to both the Mach number ∈(0, ∈] and time t ∈ [0, ∞) are established by deriving a differential inequality with decay property, where ∈∈(0, 1] is a constant.As the Mach number vanishes, the global solution to full compressible Navier-Stokes equations converges to the one of isentropic incompressible Navier-Stokes equations in t ∈ [0, +∞). Moreover, we prove the exponentially asymptotic stability for the global solutions of both the compressible system and its limiting incompressible system.
