We present a class of arbitrarily high order fully explicit kinetic numerical methods in compressible fluid dynamics,both in time and space,which include the relaxation schemes by Jin and Xin.These methods can use the...We present a class of arbitrarily high order fully explicit kinetic numerical methods in compressible fluid dynamics,both in time and space,which include the relaxation schemes by Jin and Xin.These methods can use the CFL number larger or equal to unity on regular Cartesian meshes for the multi-dimensional case.These kinetic models depend on a small parameter that can be seen as a"Knudsen"number.The method is asymptotic preserving in this Knudsen number.Also,the computational costs of the method are of the same order of a fully explicit scheme.This work is the extension of Abgrall et al.(2022)[3]to multidimensional systems.We have assessed our method on several problems for two-dimensional scalar problems and Euler equations and the scheme has proven to be robust and to achieve the theoretically predicted high order of accuracy on smooth solutions.展开更多
In this paper,we consider the weak solutions of compressible Navier-StokesLandau-Lifshitz-Maxwell(CNSLLM)system for quantum fluids with a linear density dependent viscosity in a 3D torus.By introducing the cold pressu...In this paper,we consider the weak solutions of compressible Navier-StokesLandau-Lifshitz-Maxwell(CNSLLM)system for quantum fluids with a linear density dependent viscosity in a 3D torus.By introducing the cold pressure Pc,we prove the global existence of weak solutions with the pressure P+Pc,where P=Aργwithγ≥1.Our main result extends the one in[13]on the quantum Navier-Stokes equations to the CNSLLM system.展开更多
In this paper, we study the non-isentropic compressible magnetohydrodynamic system with a time periodic external force in R^n. Under the condition that the optimal time decay rates are obtained by spectral analysis, w...In this paper, we study the non-isentropic compressible magnetohydrodynamic system with a time periodic external force in R^n. Under the condition that the optimal time decay rates are obtained by spectral analysis, we show that the existence, uniqueness and time-asymptotic stability of time periodic solutions when the space dimension n 〉 5. Our proof is based on a combination of the energy method and the contraction mapping theorem.展开更多
In this paper, the aim is to establish the local existence of classical solutions for a class of compressible non-Newtonian fluids with vacuum in one-dimensional bounded intervals, under the assumption that the data s...In this paper, the aim is to establish the local existence of classical solutions for a class of compressible non-Newtonian fluids with vacuum in one-dimensional bounded intervals, under the assumption that the data satisfies a natural compatibility condition. For the results, the initial density does not need to be bounded below away from zero.展开更多
The aims of this paper are to discuss global existence and uniqueness of strong solution for a class of isentropic compressible navier-Stokes equations with non-Newtonian in one-dimensional bounded intervals. We prove...The aims of this paper are to discuss global existence and uniqueness of strong solution for a class of isentropic compressible navier-Stokes equations with non-Newtonian in one-dimensional bounded intervals. We prove two global existence results on strong solutions of isentropic compressible Navier-Stokes equations. The first result shows only the existence. And the second one shows the existence and uniqueness result based on the first result, but the uniqueness requires some compatibility condition.展开更多
Recent years the modify ghost fluid method (MGFM) and the real ghost fluid method (RGFM) based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have on...Recent years the modify ghost fluid method (MGFM) and the real ghost fluid method (RGFM) based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.展开更多
This paper investigates the large-time behavior of solutions to an outflow problem for a compressible non-Newtonian fluid in a half space. The main concern is to analyze the phenomena that happens when the compressibl...This paper investigates the large-time behavior of solutions to an outflow problem for a compressible non-Newtonian fluid in a half space. The main concern is to analyze the phenomena that happens when the compressible non-Newtonian fluid blows out through the boundary. Based on the existence of the stationary solution, it is proved that there exists a boundary layer(i.e., the stationary solution) to the outflow problem and the boundary layer is nonlinearly stable under small initial perturbation.展开更多
It is known to all, the spilling of pipeline may cause serious problems, especially when the pipe conveying petroleum, natural gas or other toxic substance. There are countless accidents during past century. Once the ...It is known to all, the spilling of pipeline may cause serious problems, especially when the pipe conveying petroleum, natural gas or other toxic substance. There are countless accidents during past century. Once the spilling occurs, the vibration of the pipe would aggravate spill situation and even result in crack of the pipe. The consequence will be more severe when the fluid inside is compressible. To prevent the detriment of the spilling model is developed by assuming the leakages as orifices or nozzles and a 2-D vertical simply supported pipe is selected to analyze the phenomena of the oscillation. Combining these two models, the oscillation model for the pipe with leakage is set up and the spilling effect is analyzed by numerical method. The amplitude of the pipe oscillation and the normal stress enlarge as the internal velocity increased, while the shear stress changes very little.展开更多
In this article, we consider the blowup criterion for the local strong solution to the compressible fluid-particle interaction model in dimension three with vacuum. We establish a BKM type criterion for possible break...In this article, we consider the blowup criterion for the local strong solution to the compressible fluid-particle interaction model in dimension three with vacuum. We establish a BKM type criterion for possible breakdown of such solutions at critical time in terms of both the L^∞ (0, T; L^6)-norm of the density of particles and the ^L1(0, T; L^∞)-norm of the deformation tensor of velocity gradient.展开更多
A unified numerical scheme for the solutions of the compressible and incompressible Navier-Stokes equations is investigated based on a time-derivative preconditioning algorithm. The primitive variables are pressure, v...A unified numerical scheme for the solutions of the compressible and incompressible Navier-Stokes equations is investigated based on a time-derivative preconditioning algorithm. The primitive variables are pressure, velocities and temperature. The time integration scheme is used in conjunction with a finite volume discretization. The preconditioning is coupled with a high order implicit upwind scheme based on the definition of a Roe's type matrix. Computational capabilities are demonstrated through computations of high Mach number, middle Mach number, very low Mach number, and incompressible flow. It has also been demonstrated that the discontinuous surface in flow field can be captured for the implementation Roe's scheme.展开更多
This paper is devoted to a combined Fourier spectral-finite difference method for solving 3-dimensional, semi-periodic compressible fluid flow problem. The error estimation, as well as the convergence rate, is presented.
Two interface capturing methods are studied for multi fluid flows, governed by the stiffened gas equation of state. The mixture type interface capturing algorithm uses a simple volume fraction model Euler equations wr...Two interface capturing methods are studied for multi fluid flows, governed by the stiffened gas equation of state. The mixture type interface capturing algorithm uses a simple volume fraction model Euler equations written in a quasi conservative form, which is solved by a standard high resolution piecewise parabolic method (PPM) with multi fluid Riemann solver. The level set interface capturing method uses a narrow band ghost fluid method (GFM) with no numerical smearing. Several examples are presented and compared for one and two dimensions, which show the feasibility of the two methods applied to various multi fluid problems.展开更多
This work presents a new application for the Hierarchical Function Expansion Method for the solution of the Navier-Stokes equations for compressible fluids in two dimensions and in high velocity. This method is based ...This work presents a new application for the Hierarchical Function Expansion Method for the solution of the Navier-Stokes equations for compressible fluids in two dimensions and in high velocity. This method is based on the finite elements method using the Petrov-Galerkin formulation, know as SUPG (Streamline Upwind Petrov-Galerkin), applied with the expansion of the variables into hierarchical functions. To test and validate the numerical method proposed as well as the computational program developed simulations are performed for some cases whose theoretical solutions are known. These cases are the following: continuity test, stability and convergence test, temperature step problem, and several oblique shocks. The objective of the last cases is basically to verify the capture of the shock wave by the method developed. The results obtained in the simulations with the proposed method were good both qualitatively and quantitatively when compared with the theoretical solutions. This allows concluding that the objectives of this work are reached.展开更多
The propagation along oz of pulsed sound waves made of sequences of elementary unit pulses U (sin τ) where U is the unit step function and τ = kz -ωt is analyzed using the expansion of U (sin τ) and of the Dirac d...The propagation along oz of pulsed sound waves made of sequences of elementary unit pulses U (sin τ) where U is the unit step function and τ = kz -ωt is analyzed using the expansion of U (sin τ) and of the Dirac distribution δ (sin τ) in terms of τ-nπ where n is an integer. Their properties and how these pulsed sound waves could be generated are discussed.展开更多
In this paper,we discuss the local existence of H^i(i=2,4)solutions for a 1D compressible viscous micropolar fluid model with non-homogeneous temperature boundary.The proof is based on the local existence of solutions...In this paper,we discuss the local existence of H^i(i=2,4)solutions for a 1D compressible viscous micropolar fluid model with non-homogeneous temperature boundary.The proof is based on the local existence of solutions in[1].展开更多
The paper presents the implementation of non-Newtonian fluid properties for compressible multiphase solver in the open source framework OpenFOAM. The transport models for Power Law, Cross Power Law, Casson, Bird-Carre...The paper presents the implementation of non-Newtonian fluid properties for compressible multiphase solver in the open source framework OpenFOAM. The transport models for Power Law, Cross Power Law, Casson, Bird-Carreau and Herschel-Bulkley fluids were included in the thermophysical model library. Appropriate non-Newtonian liquids have been chosen from literature, and pressure driven test simulations are carried out. Therefore, the solver compressibleInterFoam is used to compute air-liquid mixture flows over a backward facing step. A validation of the novel models has been performed by means of a sample-based comparison of the strain rate viscosity relation. The theoretical rheological properties of the selected liquids agree well with the results of the simulated data.展开更多
The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function (DDF) lattice Boltzmann method (LBM) for the compressible Navier-Stokes equations. Firstly, ...The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function (DDF) lattice Boltzmann method (LBM) for the compressible Navier-Stokes equations. Firstly, the method is tested by simulating the low Reynolds number flow at Ma =0. 5,a=0. 0, Re=5 000. Then the simulation of flow around the airfoil is carried out at Ma:0. 5, 0. 85, 1.2; a=-0.05, 1.0, 0.0, respectively. And a better result is obtained by using a local refined grid. It reduces the error produced by the grid at Ma=0. 85. Though the inviscid boundary condition is used to avoid the problem of flow transition to turbulence at high Reynolds numbers, the pressure distribution obtained by the simulation agrees well with that of the experimental results. Thus, it proves the reliability of the method and shows its potential for the compressible flow simulation. The suecessful application to the flow around airfoil lays a foundation of the numerical simulation of turbulence.展开更多
The numerical modelling of the interactions between water waves and floating structures is significant for different areas of the marine sector, especially seakeeping and prediction of wave-induced loads. Seakeeping a...The numerical modelling of the interactions between water waves and floating structures is significant for different areas of the marine sector, especially seakeeping and prediction of wave-induced loads. Seakeeping analysis involving severe flow fluctuations is still quite challenging even for the conventional RANS method. Particle method has been viewed as alternative for such analysis especially those involving deformable boundary, wave breaking and fluid fragmentation around hull shapes. In this paper, the weakly compressible smoothed particle hydrodynamics(WCSPH), a fully Lagrangian particle method, is applied to simulate the symmetric radiation problem for a stationary barge treated as a flexible body. This is carried out by imposing prescribed forced simple harmonic oscillations in heave, pitch and the two-and three-node distortion modes. The resultant,radiation force predictions, namely added mass and fluid damping coefficients, are compared with results from 3-D potential flow boundary element method and 3-D RANS CFD predictions, in order to verify the adopted modelling techniques for WCSPH.WCSPH were found to be in agreement with most results and could predict the fluid actions equally well in most cases.展开更多
We study in this article the compressible heat-conducting Navier-Stokes equations in periodic domain driven by a time-periodic external force. The existence of the strong time-periodic solution is established by a new...We study in this article the compressible heat-conducting Navier-Stokes equations in periodic domain driven by a time-periodic external force. The existence of the strong time-periodic solution is established by a new approach. First, we reformulate the system and consider some decay estimates of the linearized system.Under some smallness and symmetry assumptions on the external force, the existence of the time-periodic solution of the linearized system is then identi?ed as the ?xed point of a Poincare′ map which is obtained by the Tychonoff ?xed point theorem.Although the Tychonoff ?xed point theorem cannot directly ensure the uniqueness,but we could construct a set-valued function, the ?xed point of which is the timeperiodic solution of the original system. At last, the existence of the ?xed point is obtained by the Kakutani ?xed point theorem. In addition, the uniqueness of timeperiodic solution is also studied.展开更多
In this paper,Runge-Kutta Discontinuous Galerkin(RKDG) finite element method is presented to solve the onedimensional inviscid compressible gas dynamic equations in a Lagrangian coordinate.The equations are discreti...In this paper,Runge-Kutta Discontinuous Galerkin(RKDG) finite element method is presented to solve the onedimensional inviscid compressible gas dynamic equations in a Lagrangian coordinate.The equations are discretized by the DG method in space and the temporal discretization is accomplished by the total variation diminishing Runge-Kutta method.A limiter based on the characteristic field decomposition is applied to maintain stability and non-oscillatory property of the RKDG method.For multi-medium fluid simulation,the two cells adjacent to the interface are treated differently from other cells.At first,a linear Riemann solver is applied to calculate the numerical ?ux at the interface.Numerical examples show that there is some oscillation in the vicinity of the interface.Then a nonlinear Riemann solver based on the characteristic formulation of the equation and the discontinuity relations is adopted to calculate the numerical ?ux at the interface,which suppresses the oscillation successfully.Several single-medium and multi-medium fluid examples are given to demonstrate the reliability and efficiency of the algorithm.展开更多
基金funded by the SNF project 200020_204917 entitled"Structure preserving and fast methods for hyperbolic systems of conservation laws".
文摘We present a class of arbitrarily high order fully explicit kinetic numerical methods in compressible fluid dynamics,both in time and space,which include the relaxation schemes by Jin and Xin.These methods can use the CFL number larger or equal to unity on regular Cartesian meshes for the multi-dimensional case.These kinetic models depend on a small parameter that can be seen as a"Knudsen"number.The method is asymptotic preserving in this Knudsen number.Also,the computational costs of the method are of the same order of a fully explicit scheme.This work is the extension of Abgrall et al.(2022)[3]to multidimensional systems.We have assessed our method on several problems for two-dimensional scalar problems and Euler equations and the scheme has proven to be robust and to achieve the theoretically predicted high order of accuracy on smooth solutions.
基金partially supported by the National Natural Sciences Foundation of China(11931010,12061003)。
文摘In this paper,we consider the weak solutions of compressible Navier-StokesLandau-Lifshitz-Maxwell(CNSLLM)system for quantum fluids with a linear density dependent viscosity in a 3D torus.By introducing the cold pressure Pc,we prove the global existence of weak solutions with the pressure P+Pc,where P=Aργwithγ≥1.Our main result extends the one in[13]on the quantum Navier-Stokes equations to the CNSLLM system.
基金Supported by National Natural Science Foundation of China(11271305)
文摘In this paper, we study the non-isentropic compressible magnetohydrodynamic system with a time periodic external force in R^n. Under the condition that the optimal time decay rates are obtained by spectral analysis, we show that the existence, uniqueness and time-asymptotic stability of time periodic solutions when the space dimension n 〉 5. Our proof is based on a combination of the energy method and the contraction mapping theorem.
基金Supported by NSFC(11201371,1331005)Natural Science Foundation of Shaanxi Province(2012JQ020)
文摘In this paper, the aim is to establish the local existence of classical solutions for a class of compressible non-Newtonian fluids with vacuum in one-dimensional bounded intervals, under the assumption that the data satisfies a natural compatibility condition. For the results, the initial density does not need to be bounded below away from zero.
文摘The aims of this paper are to discuss global existence and uniqueness of strong solution for a class of isentropic compressible navier-Stokes equations with non-Newtonian in one-dimensional bounded intervals. We prove two global existence results on strong solutions of isentropic compressible Navier-Stokes equations. The first result shows only the existence. And the second one shows the existence and uniqueness result based on the first result, but the uniqueness requires some compatibility condition.
基金supported by National Science Foundation of China (10576015)
文摘Recent years the modify ghost fluid method (MGFM) and the real ghost fluid method (RGFM) based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.
基金supported by the National Natural Science Foundation of China(11501445)the third author was supported by the National Natural Science Foundation of China(11671319,11331005)
文摘This paper investigates the large-time behavior of solutions to an outflow problem for a compressible non-Newtonian fluid in a half space. The main concern is to analyze the phenomena that happens when the compressible non-Newtonian fluid blows out through the boundary. Based on the existence of the stationary solution, it is proved that there exists a boundary layer(i.e., the stationary solution) to the outflow problem and the boundary layer is nonlinearly stable under small initial perturbation.
基金the support of Thousand Talents Programthe National Natural Science Foundation of China(51479114)special fund for Marine Renewable Energy Project(GHME2014ZC01)
文摘It is known to all, the spilling of pipeline may cause serious problems, especially when the pipe conveying petroleum, natural gas or other toxic substance. There are countless accidents during past century. Once the spilling occurs, the vibration of the pipe would aggravate spill situation and even result in crack of the pipe. The consequence will be more severe when the fluid inside is compressible. To prevent the detriment of the spilling model is developed by assuming the leakages as orifices or nozzles and a 2-D vertical simply supported pipe is selected to analyze the phenomena of the oscillation. Combining these two models, the oscillation model for the pipe with leakage is set up and the spilling effect is analyzed by numerical method. The amplitude of the pipe oscillation and the normal stress enlarge as the internal velocity increased, while the shear stress changes very little.
基金supported by the National Basic Research Program of China(973 Program)(2011CB808002)the National Natural Science Foundation of China(11371152,11128102,11071086,and 11571117)+3 种基金the Natural Science Foundation of Guangdong Province(S2012010010408)the Foundation for Distinguished Young Talents in Higher Education of Guangdong(2015KQNCX095)the Major Foundation of Hanshan Normal University(LZ201403)the Scientific Research Foundation of Graduate School of South China Normal University(2014ssxm04)
文摘In this article, we consider the blowup criterion for the local strong solution to the compressible fluid-particle interaction model in dimension three with vacuum. We establish a BKM type criterion for possible breakdown of such solutions at critical time in terms of both the L^∞ (0, T; L^6)-norm of the density of particles and the ^L1(0, T; L^∞)-norm of the deformation tensor of velocity gradient.
基金Project supported by the National Natural Science Foundation of China(No.50576049) the Foun-dational Scientific Research of National Defence of China(No.A4020060263)Shanghai Leading Academic Discipline Project(No.Y0103)
文摘A unified numerical scheme for the solutions of the compressible and incompressible Navier-Stokes equations is investigated based on a time-derivative preconditioning algorithm. The primitive variables are pressure, velocities and temperature. The time integration scheme is used in conjunction with a finite volume discretization. The preconditioning is coupled with a high order implicit upwind scheme based on the definition of a Roe's type matrix. Computational capabilities are demonstrated through computations of high Mach number, middle Mach number, very low Mach number, and incompressible flow. It has also been demonstrated that the discontinuous surface in flow field can be captured for the implementation Roe's scheme.
文摘This paper is devoted to a combined Fourier spectral-finite difference method for solving 3-dimensional, semi-periodic compressible fluid flow problem. The error estimation, as well as the convergence rate, is presented.
文摘Two interface capturing methods are studied for multi fluid flows, governed by the stiffened gas equation of state. The mixture type interface capturing algorithm uses a simple volume fraction model Euler equations written in a quasi conservative form, which is solved by a standard high resolution piecewise parabolic method (PPM) with multi fluid Riemann solver. The level set interface capturing method uses a narrow band ghost fluid method (GFM) with no numerical smearing. Several examples are presented and compared for one and two dimensions, which show the feasibility of the two methods applied to various multi fluid problems.
文摘This work presents a new application for the Hierarchical Function Expansion Method for the solution of the Navier-Stokes equations for compressible fluids in two dimensions and in high velocity. This method is based on the finite elements method using the Petrov-Galerkin formulation, know as SUPG (Streamline Upwind Petrov-Galerkin), applied with the expansion of the variables into hierarchical functions. To test and validate the numerical method proposed as well as the computational program developed simulations are performed for some cases whose theoretical solutions are known. These cases are the following: continuity test, stability and convergence test, temperature step problem, and several oblique shocks. The objective of the last cases is basically to verify the capture of the shock wave by the method developed. The results obtained in the simulations with the proposed method were good both qualitatively and quantitatively when compared with the theoretical solutions. This allows concluding that the objectives of this work are reached.
文摘The propagation along oz of pulsed sound waves made of sequences of elementary unit pulses U (sin τ) where U is the unit step function and τ = kz -ωt is analyzed using the expansion of U (sin τ) and of the Dirac distribution δ (sin τ) in terms of τ-nπ where n is an integer. Their properties and how these pulsed sound waves could be generated are discussed.
基金Supported by the NNSF of China(11271066)Supported by the grant of Shanghai Education Commission(13ZZ048)
文摘In this paper,we discuss the local existence of H^i(i=2,4)solutions for a 1D compressible viscous micropolar fluid model with non-homogeneous temperature boundary.The proof is based on the local existence of solutions in[1].
文摘The paper presents the implementation of non-Newtonian fluid properties for compressible multiphase solver in the open source framework OpenFOAM. The transport models for Power Law, Cross Power Law, Casson, Bird-Carreau and Herschel-Bulkley fluids were included in the thermophysical model library. Appropriate non-Newtonian liquids have been chosen from literature, and pressure driven test simulations are carried out. Therefore, the solver compressibleInterFoam is used to compute air-liquid mixture flows over a backward facing step. A validation of the novel models has been performed by means of a sample-based comparison of the strain rate viscosity relation. The theoretical rheological properties of the selected liquids agree well with the results of the simulated data.
基金Supported by the Aeronautical Science Foundation of China(20061453020)Foundation for Basic Research of Northwestern Polytechnical University(03)~~
文摘The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function (DDF) lattice Boltzmann method (LBM) for the compressible Navier-Stokes equations. Firstly, the method is tested by simulating the low Reynolds number flow at Ma =0. 5,a=0. 0, Re=5 000. Then the simulation of flow around the airfoil is carried out at Ma:0. 5, 0. 85, 1.2; a=-0.05, 1.0, 0.0, respectively. And a better result is obtained by using a local refined grid. It reduces the error produced by the grid at Ma=0. 85. Though the inviscid boundary condition is used to avoid the problem of flow transition to turbulence at high Reynolds numbers, the pressure distribution obtained by the simulation agrees well with that of the experimental results. Thus, it proves the reliability of the method and shows its potential for the compressible flow simulation. The suecessful application to the flow around airfoil lays a foundation of the numerical simulation of turbulence.
基金funded by the Ministry of Higher Education(MOHE)of Malaysia under the Fundamental Research Grant Scheme(FRGS)No.FRGS17-042-0608
文摘The numerical modelling of the interactions between water waves and floating structures is significant for different areas of the marine sector, especially seakeeping and prediction of wave-induced loads. Seakeeping analysis involving severe flow fluctuations is still quite challenging even for the conventional RANS method. Particle method has been viewed as alternative for such analysis especially those involving deformable boundary, wave breaking and fluid fragmentation around hull shapes. In this paper, the weakly compressible smoothed particle hydrodynamics(WCSPH), a fully Lagrangian particle method, is applied to simulate the symmetric radiation problem for a stationary barge treated as a flexible body. This is carried out by imposing prescribed forced simple harmonic oscillations in heave, pitch and the two-and three-node distortion modes. The resultant,radiation force predictions, namely added mass and fluid damping coefficients, are compared with results from 3-D potential flow boundary element method and 3-D RANS CFD predictions, in order to verify the adopted modelling techniques for WCSPH.WCSPH were found to be in agreement with most results and could predict the fluid actions equally well in most cases.
基金The NSF(20170520047JH) for Young Scientists of Jilin Provincethe Scientific and Technological Project(JJKH20190180KJ) of Jilin Province’s Education Department in Thirteenth Five-Year
文摘We study in this article the compressible heat-conducting Navier-Stokes equations in periodic domain driven by a time-periodic external force. The existence of the strong time-periodic solution is established by a new approach. First, we reformulate the system and consider some decay estimates of the linearized system.Under some smallness and symmetry assumptions on the external force, the existence of the time-periodic solution of the linearized system is then identi?ed as the ?xed point of a Poincare′ map which is obtained by the Tychonoff ?xed point theorem.Although the Tychonoff ?xed point theorem cannot directly ensure the uniqueness,but we could construct a set-valued function, the ?xed point of which is the timeperiodic solution of the original system. At last, the existence of the ?xed point is obtained by the Kakutani ?xed point theorem. In addition, the uniqueness of timeperiodic solution is also studied.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 11261035,11171038,and 10771019)the Science Reaearch Foundation of Institute of Higher Education of Inner Mongolia Autonomous Region,China (Grant No. NJZZ12198)the Natural Science Foundation of Inner Mongolia Autonomous Region,China (Grant No. 2012MS0102)
文摘In this paper,Runge-Kutta Discontinuous Galerkin(RKDG) finite element method is presented to solve the onedimensional inviscid compressible gas dynamic equations in a Lagrangian coordinate.The equations are discretized by the DG method in space and the temporal discretization is accomplished by the total variation diminishing Runge-Kutta method.A limiter based on the characteristic field decomposition is applied to maintain stability and non-oscillatory property of the RKDG method.For multi-medium fluid simulation,the two cells adjacent to the interface are treated differently from other cells.At first,a linear Riemann solver is applied to calculate the numerical ?ux at the interface.Numerical examples show that there is some oscillation in the vicinity of the interface.Then a nonlinear Riemann solver based on the characteristic formulation of the equation and the discontinuity relations is adopted to calculate the numerical ?ux at the interface,which suppresses the oscillation successfully.Several single-medium and multi-medium fluid examples are given to demonstrate the reliability and efficiency of the algorithm.