In this article we will present pure three dimensional analytic solutions for the Navier-Stokes and the continuity equations in Cartesian coordinates. The key idea is the three-dimensional generalization of the well-k...In this article we will present pure three dimensional analytic solutions for the Navier-Stokes and the continuity equations in Cartesian coordinates. The key idea is the three-dimensional generalization of the well-known self-similar Ansatz of Barenblatt. A geometrical interpretation of the Ansatz is given also. The results are the Kummer functions or strongly related. Our final formula is compared with other results obtained from group theoretical approaches.展开更多
In this paper, homotopy analysis method (HAM) and Padé approximant will be considered for finding analytical solution of three-dimensional viscous flow near an infinite rotating disk which is a well-known classic...In this paper, homotopy analysis method (HAM) and Padé approximant will be considered for finding analytical solution of three-dimensional viscous flow near an infinite rotating disk which is a well-known classical problem in fluid mechanics. The solution is compared to the numerical (fourth-order Runge-Kutta) solution and the convergence of the obtained series solution is carefully analyzed. The results illustrate that HAM-Padé is an appropriate method in solving the systems of nonlinear equations.展开更多
This article gives a general model using specific periodic special functions, that is, degenerate elliptic Weierstrass P functions composed with the LambertW function, whose presence in the governing equations through...This article gives a general model using specific periodic special functions, that is, degenerate elliptic Weierstrass P functions composed with the LambertW function, whose presence in the governing equations through the forcing terms simplify the periodic Navier Stokes equations (PNS) at the centers of arbitrary r balls of the 3-Torus. The continuity equation is satisfied together with spatially periodic boundary conditions. The yicomponent forcing terms consist of a function F as part of its expression that is arbitrarily small in an r ball where it is associated with a singular forcing expression both for inviscid and viscous cases. As a result, a significant simplification occurs with a v3(vifor all velocity components) only governing PDE resulting. The extension of three restricted subspaces in each of the principal directions in the Cartesian plane is shown as the Cartesian product ℋ=Jx,t×Jy,t×Jz,t. On each of these subspaces vi,i=1,2,3is continuous and there exists a linear independent subspace associated with the argument of the W function. Here the 3-Torus is built up from each compact segment of length 2R on each of the axes on the 3 principal directions x, y, and z. The form of the scaled velocities for non zero scaled δis related to the definition of the W function such that e−W(ξ)=W(ξ)ξwhere ξdepends on t and proportional to δ→0for infinite time t. The ratio Wξis equal to 1, making the limit δ→0finite and well defined. Considering r balls where the function F=(x−ai)2+(y−bi)2+(z−ci)2−ηset equal to −1e+rwhere r>0. is such that the forcing is singular at every distance r of centres of cubes each containing an r-ball. At the centre of the balls, the forcing is infinite. The main idea is that a system of singular initial value problems with infinite forcing is to be solved for where the velocities are shown to be locally Hölder continuous. It is proven that the limit of these singular problems shifts the finite time blowup time ti∗for first and higher derivatives to t=∞thereby indicating that there is no finite time blowup. Results in the literature can provide a systematic approach to study both large space and time behaviour for singular solutions to the Navier Stokes equations. Among the references, it has been shown that mathematical tools can be applied to study the asymptotic properties of solutions.展开更多
In this paper, the theory of constructing optimal dynamical systems based on weighted residual presented by Wu & Sha is applied to three-dimensional Navier-Stokes equations, and the optimal dynamical system modeli...In this paper, the theory of constructing optimal dynamical systems based on weighted residual presented by Wu & Sha is applied to three-dimensional Navier-Stokes equations, and the optimal dynamical system modeling equations are derived. Then the multiscale global optimization method based on coarse graining analysis is presented, by which a set of approximate global optimal bases is directly obtained from Navier-Stokes equations and the construction of optimal dynamical systems is realized. The optimal bases show good properties, such as showing the physical properties of complex flows and the turbulent vortex structures, being intrinsic to real physical problem and dynamical systems, and having scaling symmetry in mathematics, etc.. In conclusion, using fewer terms of optimal bases will approach the exact solutions of Navier-Stokes equations, and the dynamical systems based on them show the most optimal behavior.展开更多
A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and...A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.展开更多
In the study of the regularity criteria for Leray weak solutions to threedimensional Navier-Stokes equations, two sufficient conditions such that the horizontal velocity u satisfies u∈L2(0,T;BMO(R3)) or u∈L^2/1...In the study of the regularity criteria for Leray weak solutions to threedimensional Navier-Stokes equations, two sufficient conditions such that the horizontal velocity u satisfies u∈L2(0,T;BMO(R3)) or u∈L^2/1+r(0,T;B∞,∞(R3)) for 0 〈 r 〈 1 are considered.展开更多
This paper is devoted to the mixed Legendre spectral-finite element approximation of the three-dimensional, non-periodic, unsteady Navier-Stokes equations. A class of fully discrete schemes are constructed with artifi...This paper is devoted to the mixed Legendre spectral-finite element approximation of the three-dimensional, non-periodic, unsteady Navier-Stokes equations. A class of fully discrete schemes are constructed with artificial compression. The generalized stability and convergence are proved strictly on the assumption that the two-dimensional inf-sup condition of the finite element approximation is satisfied.展开更多
The generalized Navier-Stokes equations with damping are considered.We will show that the generalized Navier-Stokes equations with damping |u|^β-1u have weak solutions for anyβ>1 and 0<α<5/4,and we will us...The generalized Navier-Stokes equations with damping are considered.We will show that the generalized Navier-Stokes equations with damping |u|^β-1u have weak solutions for anyβ>1 and 0<α<5/4,and we will use the Fourier splitting method to prove the L2 decay of weak solutions forβ>2 and 0<α<3/4.展开更多
A global weak solution to the isentropic Navier-Stokes equation with initial data around a constant state in the L^(1)∩BV class was constructed in[1].In the current paper,we will continue to study the uniqueness and ...A global weak solution to the isentropic Navier-Stokes equation with initial data around a constant state in the L^(1)∩BV class was constructed in[1].In the current paper,we will continue to study the uniqueness and regularity of the constructed solution.The key ingredients are the Holder continuity estimates of the heat kernel in both spatial and time variables.With these finer estimates,we obtain higher order regularity of the constructed solution to Navier-Stokes equation,so that all of the derivatives in the equation of conservative form are in the strong sense.Moreover,this regularity also allows us to identify a function space such that the stability of the solutions can be established there,which eventually implies the uniqueness.展开更多
In this paper,we study the controllability of compressible Navier-Stokes equations with density dependent viscosities.For when the shear viscosityμis a positive constant and the bulk viscosityλis a function of the d...In this paper,we study the controllability of compressible Navier-Stokes equations with density dependent viscosities.For when the shear viscosityμis a positive constant and the bulk viscosityλis a function of the density,it is proven that the system is exactly locally controllable to a constant target trajectory by using boundary control functions.展开更多
The analytical solution of the multi-dimensional,time-fractional model of Navier-Stokes equation using the triple and quadruple Elzaki transformdecompositionmethod is presented in this article.The aforesaidmodel is an...The analytical solution of the multi-dimensional,time-fractional model of Navier-Stokes equation using the triple and quadruple Elzaki transformdecompositionmethod is presented in this article.The aforesaidmodel is analyzed by employing Caputo fractional derivative.We deliberated three stimulating examples that correspond to the triple and quadruple Elzaki transform decomposition methods,respectively.The findings illustrate that the established approaches are extremely helpful in obtaining exact and approximate solutions to the problems.The exact and estimated solutions are delineated via numerical simulation.The proposed analysis indicates that the projected configuration is extremely meticulous,highly efficient,and precise in understanding the behavior of complex evolutionary problems of both fractional and integer order that classify affiliated scientific fields and technology.展开更多
We consider the global well-posedness of strong solutions to the Cauchy problem of compressible barotropic Navier-Stokes equations in R^(2). By exploiting the global-in-time estimate to the two-dimensional(2D for shor...We consider the global well-posedness of strong solutions to the Cauchy problem of compressible barotropic Navier-Stokes equations in R^(2). By exploiting the global-in-time estimate to the two-dimensional(2D for short) classical incompressible Navier-Stokes equations and using techniques developed in(SIAM J Math Anal, 2020, 52(2): 1806–1843), we derive the global existence of solutions provided that the initial data satisfies some smallness condition. In particular, the initial velocity with some arbitrary Besov norm of potential part Pu_0 and large high oscillation are allowed in our results. Moreover, we also construct an example with the initial data involving such a smallness condition, but with a norm that is arbitrarily large.展开更多
We investigate the low Mach number limit for the isentropic compressible NavierStokes equations with a revised Maxwell's law(with Galilean invariance) in R^(3). By applying the uniform estimates of the error syste...We investigate the low Mach number limit for the isentropic compressible NavierStokes equations with a revised Maxwell's law(with Galilean invariance) in R^(3). By applying the uniform estimates of the error system, it is proven that the solutions of the isentropic Navier-Stokes equations with a revised Maxwell's law converge to that of the incompressible Navier-Stokes equations as the Mach number tends to zero. Moreover, the convergence rates are also obtained.展开更多
Many applications in fluid mechanics require the numerical solution of sequences of linear systems typically issued from finite element discretization of the Navier-Stokes equations. The resulting matrices then exhibi...Many applications in fluid mechanics require the numerical solution of sequences of linear systems typically issued from finite element discretization of the Navier-Stokes equations. The resulting matrices then exhibit a saddle point structure. To achieve this task, a Newton-based root-finding algorithm is usually employed which in turn necessitates to solve a saddle point system at every Newton iteration. The involved linear systems being large scale and ill-conditioned, effective linear solvers must be implemented. Here, we develop and test several methods for solving the saddle point systems, considering in particular the LU factorization, as direct approach, and the preconditioned generalized minimal residual (ΡGMRES) solver, an iterative approach. We apply the various solvers within the root-finding algorithm for Flow over backward facing step systems. The particularity of Flow over backward facing step system is an interesting case for studying the performance and solution strategy of a turbulence model. In this case, the flow is subjected to a sudden increase of cross-sectional area, resulting in a separation of flow starting at the point of expansion, making the system of differential equations particularly stiff. We assess the performance of the direct and iterative solvers in terms of computational time, numbers of Newton iterations and time steps.展开更多
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.展开更多
This paper is concerned with the global well-posedness of the solution to the compressible Navier-Stokes/Allen-Cahn system and its sharp interface limit in one-dimensional space.For the perturbations with small energy...This paper is concerned with the global well-posedness of the solution to the compressible Navier-Stokes/Allen-Cahn system and its sharp interface limit in one-dimensional space.For the perturbations with small energy but possibly large oscillations of rarefaction wave solutions near phase separation,and where the strength of the initial phase field could be arbitrarily large,we prove that the solution of the Cauchy problem exists for all time,and converges to the centered rarefaction wave solution of the corresponding standard two-phase Euler equation as the viscosity and the thickness of the interface tend to zero.The proof is mainly based on a scaling argument and a basic energy method.展开更多
In this paper,we study the one-dimensional motion of viscous gas near a vacuum,with the gas connecting to a vacuum state with a jump in density.The interface behavior,the pointwise decay rates of the density function ...In this paper,we study the one-dimensional motion of viscous gas near a vacuum,with the gas connecting to a vacuum state with a jump in density.The interface behavior,the pointwise decay rates of the density function and the expanding rates of the interface are obtained with the viscosity coefficientμ(ρ)=ρ^(α)for any 0<α<1;this includes the timeweighted boundedness from below and above.The smoothness of the solution is discussed.Moreover,we construct a class of self-similar classical solutions which exhibit some interesting properties,such as optimal estimates.The present paper extends the results in[Luo T,Xin Z P,Yang T.SIAM J Math Anal,2000,31(6):1175-1191]to the jump boundary conditions case with density-dependent viscosity.展开更多
The preconditioning method is used to solve the low Mach number flow. The space discritisation scheme is the Roe scheme and the DES turbulence model is used. Then, the low Mach number turbulence flow around the NACA00...The preconditioning method is used to solve the low Mach number flow. The space discritisation scheme is the Roe scheme and the DES turbulence model is used. Then, the low Mach number turbulence flow around the NACA0012 airfoil is used to verify the efficiency of the proposed method. Two cases of the low Mach number flows around the multi-element airfoil and the circular cylinder are also used to test the proposed method. Numerical results show that the methods combined the preconditioning method and compressible Navier-Stokes equations are efficient to solve low Mach number flows.展开更多
In this paper,we study a one-dimensional motion of viscous gas near vacuum. We are interested in the case that the gas is in contact with the vacuum at a finite interval. This is a free boundary problem for the one-di...In this paper,we study a one-dimensional motion of viscous gas near vacuum. We are interested in the case that the gas is in contact with the vacuum at a finite interval. This is a free boundary problem for the one-dimensional isentropic Navier-Stokes equations, and the free boundaries are the interfaces separating the gas from vacuum,across which the density changes discontinuosly.Smoothness of the solutions and the uniqueness of the weak solutions are also discussed.The present paper extends results in Luo-Xin-Yang[12] to the jump boundary conditions case.展开更多
The zero dissipation limit of the compressible heat-conducting Navier–Stokes equations in the presence of the shock is investigated. It is shown that when the heat conduction coefficient κ and the viscosity coeffici...The zero dissipation limit of the compressible heat-conducting Navier–Stokes equations in the presence of the shock is investigated. It is shown that when the heat conduction coefficient κ and the viscosity coefficient ε satisfy κ = O(ε), κ/ε≥ c 〉 0, as ε→ 0 (see (1.3)), if the solution of the corresponding Euler equations is piecewise smooth with shock wave satisfying the Lax entropy condition, then there exists a smooth solution to the Navier–Stokes equations, which converges to the piecewise smooth shock solution of the Euler equations away from the shock discontinuity at a rate of ε. The proof is given by a combination of the energy estimates and the matched asymptotic analysis introduced in [3].展开更多
文摘In this article we will present pure three dimensional analytic solutions for the Navier-Stokes and the continuity equations in Cartesian coordinates. The key idea is the three-dimensional generalization of the well-known self-similar Ansatz of Barenblatt. A geometrical interpretation of the Ansatz is given also. The results are the Kummer functions or strongly related. Our final formula is compared with other results obtained from group theoretical approaches.
文摘In this paper, homotopy analysis method (HAM) and Padé approximant will be considered for finding analytical solution of three-dimensional viscous flow near an infinite rotating disk which is a well-known classical problem in fluid mechanics. The solution is compared to the numerical (fourth-order Runge-Kutta) solution and the convergence of the obtained series solution is carefully analyzed. The results illustrate that HAM-Padé is an appropriate method in solving the systems of nonlinear equations.
文摘This article gives a general model using specific periodic special functions, that is, degenerate elliptic Weierstrass P functions composed with the LambertW function, whose presence in the governing equations through the forcing terms simplify the periodic Navier Stokes equations (PNS) at the centers of arbitrary r balls of the 3-Torus. The continuity equation is satisfied together with spatially periodic boundary conditions. The yicomponent forcing terms consist of a function F as part of its expression that is arbitrarily small in an r ball where it is associated with a singular forcing expression both for inviscid and viscous cases. As a result, a significant simplification occurs with a v3(vifor all velocity components) only governing PDE resulting. The extension of three restricted subspaces in each of the principal directions in the Cartesian plane is shown as the Cartesian product ℋ=Jx,t×Jy,t×Jz,t. On each of these subspaces vi,i=1,2,3is continuous and there exists a linear independent subspace associated with the argument of the W function. Here the 3-Torus is built up from each compact segment of length 2R on each of the axes on the 3 principal directions x, y, and z. The form of the scaled velocities for non zero scaled δis related to the definition of the W function such that e−W(ξ)=W(ξ)ξwhere ξdepends on t and proportional to δ→0for infinite time t. The ratio Wξis equal to 1, making the limit δ→0finite and well defined. Considering r balls where the function F=(x−ai)2+(y−bi)2+(z−ci)2−ηset equal to −1e+rwhere r>0. is such that the forcing is singular at every distance r of centres of cubes each containing an r-ball. At the centre of the balls, the forcing is infinite. The main idea is that a system of singular initial value problems with infinite forcing is to be solved for where the velocities are shown to be locally Hölder continuous. It is proven that the limit of these singular problems shifts the finite time blowup time ti∗for first and higher derivatives to t=∞thereby indicating that there is no finite time blowup. Results in the literature can provide a systematic approach to study both large space and time behaviour for singular solutions to the Navier Stokes equations. Among the references, it has been shown that mathematical tools can be applied to study the asymptotic properties of solutions.
基金supported by the National Natural Science Foundation of China(Grant Nos.11372068 and 11572350)the National Basic Research Program of China(Grant No.2014CB744104)
文摘In this paper, the theory of constructing optimal dynamical systems based on weighted residual presented by Wu & Sha is applied to three-dimensional Navier-Stokes equations, and the optimal dynamical system modeling equations are derived. Then the multiscale global optimization method based on coarse graining analysis is presented, by which a set of approximate global optimal bases is directly obtained from Navier-Stokes equations and the construction of optimal dynamical systems is realized. The optimal bases show good properties, such as showing the physical properties of complex flows and the turbulent vortex structures, being intrinsic to real physical problem and dynamical systems, and having scaling symmetry in mathematics, etc.. In conclusion, using fewer terms of optimal bases will approach the exact solutions of Navier-Stokes equations, and the dynamical systems based on them show the most optimal behavior.
文摘A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.
文摘In the study of the regularity criteria for Leray weak solutions to threedimensional Navier-Stokes equations, two sufficient conditions such that the horizontal velocity u satisfies u∈L2(0,T;BMO(R3)) or u∈L^2/1+r(0,T;B∞,∞(R3)) for 0 〈 r 〈 1 are considered.
文摘This paper is devoted to the mixed Legendre spectral-finite element approximation of the three-dimensional, non-periodic, unsteady Navier-Stokes equations. A class of fully discrete schemes are constructed with artificial compression. The generalized stability and convergence are proved strictly on the assumption that the two-dimensional inf-sup condition of the finite element approximation is satisfied.
文摘The generalized Navier-Stokes equations with damping are considered.We will show that the generalized Navier-Stokes equations with damping |u|^β-1u have weak solutions for anyβ>1 and 0<α<5/4,and we will use the Fourier splitting method to prove the L2 decay of weak solutions forβ>2 and 0<α<3/4.
基金partially the National Key R&D Program of China(2022YFA1007300)the NSFC(11901386,12031013)+2 种基金the Strategic Priority Research Program of the Chinese Academy of Sciences(XDA25010403)the NSFC(11801194,11971188)the Hubei Key Laboratory of Engineering Modeling and Scientific Computing。
文摘A global weak solution to the isentropic Navier-Stokes equation with initial data around a constant state in the L^(1)∩BV class was constructed in[1].In the current paper,we will continue to study the uniqueness and regularity of the constructed solution.The key ingredients are the Holder continuity estimates of the heat kernel in both spatial and time variables.With these finer estimates,we obtain higher order regularity of the constructed solution to Navier-Stokes equation,so that all of the derivatives in the equation of conservative form are in the strong sense.Moreover,this regularity also allows us to identify a function space such that the stability of the solutions can be established there,which eventually implies the uniqueness.
基金partially supported by the National Science Foundation of China(11971320,11971496)the National Key R&D Program of China(2020YFA0712500)the Guangdong Basic and Applied Basic Research Foundation(2020A1515010530)。
文摘In this paper,we study the controllability of compressible Navier-Stokes equations with density dependent viscosities.For when the shear viscosityμis a positive constant and the bulk viscosityλis a function of the density,it is proven that the system is exactly locally controllable to a constant target trajectory by using boundary control functions.
基金supported by the Natural Science Foundation of China(GrantNos.61673169,11301127,11701176,11626101,11601485).
文摘The analytical solution of the multi-dimensional,time-fractional model of Navier-Stokes equation using the triple and quadruple Elzaki transformdecompositionmethod is presented in this article.The aforesaidmodel is analyzed by employing Caputo fractional derivative.We deliberated three stimulating examples that correspond to the triple and quadruple Elzaki transform decomposition methods,respectively.The findings illustrate that the established approaches are extremely helpful in obtaining exact and approximate solutions to the problems.The exact and estimated solutions are delineated via numerical simulation.The proposed analysis indicates that the projected configuration is extremely meticulous,highly efficient,and precise in understanding the behavior of complex evolutionary problems of both fractional and integer order that classify affiliated scientific fields and technology.
基金Zhai was partially supported by the Guangdong Provincial Natural Science Foundation (2022A1515011977)the Science and Technology Program of Shenzhen(20200806104726001)+1 种基金Zhong was partially supported by the NNSF of China (11901474, 12071359)the Exceptional Young Talents Project of Chongqing Talent (cstc2021ycjh-bgzxm0153)。
文摘We consider the global well-posedness of strong solutions to the Cauchy problem of compressible barotropic Navier-Stokes equations in R^(2). By exploiting the global-in-time estimate to the two-dimensional(2D for short) classical incompressible Navier-Stokes equations and using techniques developed in(SIAM J Math Anal, 2020, 52(2): 1806–1843), we derive the global existence of solutions provided that the initial data satisfies some smallness condition. In particular, the initial velocity with some arbitrary Besov norm of potential part Pu_0 and large high oscillation are allowed in our results. Moreover, we also construct an example with the initial data involving such a smallness condition, but with a norm that is arbitrarily large.
基金Yuxi HU was supported by the NNSFC (11701556)the Yue Qi Young Scholar ProjectChina University of Mining and Technology (Beijing)。
文摘We investigate the low Mach number limit for the isentropic compressible NavierStokes equations with a revised Maxwell's law(with Galilean invariance) in R^(3). By applying the uniform estimates of the error system, it is proven that the solutions of the isentropic Navier-Stokes equations with a revised Maxwell's law converge to that of the incompressible Navier-Stokes equations as the Mach number tends to zero. Moreover, the convergence rates are also obtained.
文摘Many applications in fluid mechanics require the numerical solution of sequences of linear systems typically issued from finite element discretization of the Navier-Stokes equations. The resulting matrices then exhibit a saddle point structure. To achieve this task, a Newton-based root-finding algorithm is usually employed which in turn necessitates to solve a saddle point system at every Newton iteration. The involved linear systems being large scale and ill-conditioned, effective linear solvers must be implemented. Here, we develop and test several methods for solving the saddle point systems, considering in particular the LU factorization, as direct approach, and the preconditioned generalized minimal residual (ΡGMRES) solver, an iterative approach. We apply the various solvers within the root-finding algorithm for Flow over backward facing step systems. The particularity of Flow over backward facing step system is an interesting case for studying the performance and solution strategy of a turbulence model. In this case, the flow is subjected to a sudden increase of cross-sectional area, resulting in a separation of flow starting at the point of expansion, making the system of differential equations particularly stiff. We assess the performance of the direct and iterative solvers in terms of computational time, numbers of Newton iterations and time steps.
文摘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 National Natural Science Foundation of China(12361044)supported by the National Natural Science Foundation of China(12171024,11971217,11971020)supported by the Academic and Technical Leaders Training Plan of Jiangxi Province(20212BCJ23027)。
文摘This paper is concerned with the global well-posedness of the solution to the compressible Navier-Stokes/Allen-Cahn system and its sharp interface limit in one-dimensional space.For the perturbations with small energy but possibly large oscillations of rarefaction wave solutions near phase separation,and where the strength of the initial phase field could be arbitrarily large,we prove that the solution of the Cauchy problem exists for all time,and converges to the centered rarefaction wave solution of the corresponding standard two-phase Euler equation as the viscosity and the thickness of the interface tend to zero.The proof is mainly based on a scaling argument and a basic energy method.
基金supported by the NSFC(11931013)the GXNSF(2022GXNSFDA035078)。
文摘In this paper,we study the one-dimensional motion of viscous gas near a vacuum,with the gas connecting to a vacuum state with a jump in density.The interface behavior,the pointwise decay rates of the density function and the expanding rates of the interface are obtained with the viscosity coefficientμ(ρ)=ρ^(α)for any 0<α<1;this includes the timeweighted boundedness from below and above.The smoothness of the solution is discussed.Moreover,we construct a class of self-similar classical solutions which exhibit some interesting properties,such as optimal estimates.The present paper extends the results in[Luo T,Xin Z P,Yang T.SIAM J Math Anal,2000,31(6):1175-1191]to the jump boundary conditions case with density-dependent viscosity.
文摘The preconditioning method is used to solve the low Mach number flow. The space discritisation scheme is the Roe scheme and the DES turbulence model is used. Then, the low Mach number turbulence flow around the NACA0012 airfoil is used to verify the efficiency of the proposed method. Two cases of the low Mach number flows around the multi-element airfoil and the circular cylinder are also used to test the proposed method. Numerical results show that the methods combined the preconditioning method and compressible Navier-Stokes equations are efficient to solve low Mach number flows.
文摘In this paper,we study a one-dimensional motion of viscous gas near vacuum. We are interested in the case that the gas is in contact with the vacuum at a finite interval. This is a free boundary problem for the one-dimensional isentropic Navier-Stokes equations, and the free boundaries are the interfaces separating the gas from vacuum,across which the density changes discontinuosly.Smoothness of the solutions and the uniqueness of the weak solutions are also discussed.The present paper extends results in Luo-Xin-Yang[12] to the jump boundary conditions case.
基金the Knowledge Innovation Program of the Chinese Academy of Sciences
文摘The zero dissipation limit of the compressible heat-conducting Navier–Stokes equations in the presence of the shock is investigated. It is shown that when the heat conduction coefficient κ and the viscosity coefficient ε satisfy κ = O(ε), κ/ε≥ c 〉 0, as ε→ 0 (see (1.3)), if the solution of the corresponding Euler equations is piecewise smooth with shock wave satisfying the Lax entropy condition, then there exists a smooth solution to the Navier–Stokes equations, which converges to the piecewise smooth shock solution of the Euler equations away from the shock discontinuity at a rate of ε. The proof is given by a combination of the energy estimates and the matched asymptotic analysis introduced in [3].