A three-dimensional (3D) predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is ...A three-dimensional (3D) predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.展开更多
Flow in tidal rivers periodically propagates upstream or downstream under tidal influence. Hydrodynamic models based on the Saint-Venant equations (the SVN model) are extensively used to model tidal rivers. A force-...Flow in tidal rivers periodically propagates upstream or downstream under tidal influence. Hydrodynamic models based on the Saint-Venant equations (the SVN model) are extensively used to model tidal rivers. A force-corrected term expressed as the combination of flow velocity and the change rate of the tidal fevel was developed to represent tidal effects in the SVN model. A momentum equation incorporating with the corrected term was derived based on Newton's second law. By combing the modified momentum equation with the continuity equation, an improved SVN model for tidal rivers (the ISVN model) was constructed. The simulation of a tidal reach of the Qiantang River shows that the ISVN model performs better than the SVN model. It indicates that the corrected force derived for tidal effects is reasonable; the ISVN model provides an appropriate enhancement of the SVN model for flow simulation of tidal rivers.展开更多
In this paper,we consider the two-dimensional aggregation equation with the shear flow and time-space nonlocal attractive operator.Without the advection,the solution of the aggregation equation may blow up in finite t...In this paper,we consider the two-dimensional aggregation equation with the shear flow and time-space nonlocal attractive operator.Without the advection,the solution of the aggregation equation may blow up in finite time.We show that the shear flow can suppress the blow-up.展开更多
In this paper, the nonreflecting boundary conditions based upon fundamental ideas of the linear analysis are developed for gas dynamic equations, and the modified boundary conditions for Navier-Stokes equations are pr...In this paper, the nonreflecting boundary conditions based upon fundamental ideas of the linear analysis are developed for gas dynamic equations, and the modified boundary conditions for Navier-Stokes equations are proposed as a substitute of the nonreflecting boundary conditions inside boundary layers near rigid walls. These derived boundary conditions are then applied to calculations both for the Euler equations and the Navier-Stokes equations to determine if they can produce acceptable results for the subsonic flows in channels. The numerical results obtained by an implicit second-order upwind difference scheme show the effective- ness and generality of the boundary conditions. Furthermore, the formulae and the analysis performed here may be extended to three dimensional problems.展开更多
The short-range property of interactions between scales in incompressible turbulent flow was examined. Some formulae for the short-range eddy stress were given. A concept of resonant-range interactions between extreme...The short-range property of interactions between scales in incompressible turbulent flow was examined. Some formulae for the short-range eddy stress were given. A concept of resonant-range interactions between extremely contiguous scales was introduced and some formulae for the resonant-range eddy stress were also derived. Multi-scale equations for the incompressible turbulent flows were proposed. Key words turbulence - incompressible flow - interactions between scales - multi-scale equations MSC 2000 76F70展开更多
This paper concerns the sonic-supersonic structures of the transonic crossflow generated by the steady supersonic flow past an infinite cone of arbitrary cross section.Under the conical assumption,the three-dimensiona...This paper concerns the sonic-supersonic structures of the transonic crossflow generated by the steady supersonic flow past an infinite cone of arbitrary cross section.Under the conical assumption,the three-dimensional(3-D)steady Euler equations can be projected onto the unit sphere and the state of fluid can be characterized by the polar and azimuthal angles.Given a segment smooth curve as a conical-sonic line in the polar-azimuthal angle plane,we construct a classical conical-supersonic solution near the curve under some reasonable assumptions.To overcome the difficulty caused by the parabolic degeneracy,we apply the characteristic decomposition technique to transform the Euler equations into a new degenerate hyperbolic system in a partial hodograph plane.The singular terms are isolated from the highly nonlinear complicated system and then can be handled successfully.We establish a smooth local solution to the new system in a suitable weighted metric space and then express the solution in terms of the original variables.展开更多
The accelerated method in solving optimization problems has always been an absorbing topic.Based on the fixedtime(FxT)stability of nonlinear dynamical systems,we provide a unified approach for designing FxT gradient f...The accelerated method in solving optimization problems has always been an absorbing topic.Based on the fixedtime(FxT)stability of nonlinear dynamical systems,we provide a unified approach for designing FxT gradient flows(FxTGFs).First,a general class of nonlinear functions in designing FxTGFs is provided.A unified method for designing first-order FxTGFs is shown under Polyak-Łjasiewicz inequality assumption,a weaker condition than strong convexity.When there exist both bounded and vanishing disturbances in the gradient flow,a specific class of nonsmooth robust FxTGFs with disturbance rejection is presented.Under the strict convexity assumption,Newton-based FxTGFs is given and further extended to solve time-varying optimization.Besides,the proposed FxTGFs are further used for solving equation-constrained optimization.Moreover,an FxT proximal gradient flow with a wide range of parameters is provided for solving nonsmooth composite optimization.To show the effectiveness of various FxTGFs,the static regret analyses for several typical FxTGFs are also provided in detail.Finally,the proposed FxTGFs are applied to solve two network problems,i.e.,the network consensus problem and solving a system linear equations,respectively,from the perspective of optimization.Particularly,by choosing component-wisely sign-preserving functions,these problems can be solved in a distributed way,which extends the existing results.The accelerated convergence and robustness of the proposed FxTGFs are validated in several numerical examples stemming from practical applications.展开更多
For a stochastic differential equation with non-Lipschitz coefficients, we construct, by Euler scheme, a measurable flow of the solution, and we prove the solution is a Markov process.
In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal e...In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.展开更多
The short-range property of interactions between scales in the compressible turbulent flow was examined. An estimation of the short-range scale scope and some formulae for the short-range eddy stress and heat transfe...The short-range property of interactions between scales in the compressible turbulent flow was examined. An estimation of the short-range scale scope and some formulae for the short-range eddy stress and heat transfer etc. were given. A concept of resonant-range interactions between extremely contiguous scales was introduced and some formulae for the resonant-range eddy stress and heat transfer etc. were also given. Multi-scale equations for the compressible turbulent flows were presented. The multi-scale equations are approximately closed and do not contain any empirical constants. The compressibility effects on turbulence are determined by the Farve averaged variables and the nonlinear relationships between the Farve- and physical-averaged variables.展开更多
A local gradient estimate for positive solutions of porous medium equations on complete noncompact Riemannian manifolds under the Ricci flow is derived. Moreover, a global gradient estimate for such equations on compa...A local gradient estimate for positive solutions of porous medium equations on complete noncompact Riemannian manifolds under the Ricci flow is derived. Moreover, a global gradient estimate for such equations on compact Riemannian manifolds is also obtained.展开更多
In this paper, the convergence turbulent flow equations are considered. By rates of solutions to the three-dimensional combining the LP-Lq estimate for the linearized equations and an elaborate energy method, the conv...In this paper, the convergence turbulent flow equations are considered. By rates of solutions to the three-dimensional combining the LP-Lq estimate for the linearized equations and an elaborate energy method, the convergence rates are obtained in various norms for the solution to the equilibrium state in the whole space when the initial perturbation of the equilibrium state is small in the H3-framework. More precisely, the optimal convergence rates of the solutions and their first-order derivatives in the L2-norm are obtained when the LP-norm of the perturbation is bounded for some p ε [1, 6).展开更多
In this paper, a new numerical scheme for solving first-order hyperbolic partial differential equations is proposed and is implemented in the simulation study of macroscopic traffic flow model with constant velocity a...In this paper, a new numerical scheme for solving first-order hyperbolic partial differential equations is proposed and is implemented in the simulation study of macroscopic traffic flow model with constant velocity and linear velocity-density relationship. Macroscopic traffic flow model is first developed by Lighthill Whitham and Richards (LWR) and used to study traffic flow by collective variables such as flow rate, velocity and density. The LWR model is treated as an initial value problem and its numerical simulations are presented using numerical schemes. A variety of numerical schemes are available in literature to solve first order hyperbolic equations. Of these the well-known ones include one-dimensional explicit: Upwind, Downwind, FTCS, and Lax-Friedrichs schemes. Having been studied carefully the space and time mesh sizes, and the patterns of all these schemes, a new scheme has been developed and named as one-dimensional explicit Tolesa numerical scheme. Tolesa numerical scheme is one of the conditionally stable and highest rates of convergence schemes. All the said numerical schemes are applied to solve advection equation pertaining traffic flows. Also the one-dimensional explicit Tolesa numerical scheme is another alternative numerical scheme to solve advection equation and apply to traffic flows model like other well-known one-dimensional explicit schemes. The effect of density of cars on the overall interactions of the vehicles along a given length of the highway and time are investigated. Graphical representations of density profile, velocity profile, flux profile, and in general the fundamental diagrams of vehicles on the highway with different time levels are illustrated. These concepts and results have been arranged systematically in this paper.展开更多
A stable high-order Runge-Kutta discontinuous Galerkin(RKDG) scheme that strictly preserves positivity of the solution is designed to solve the Boltzmann kinetic equation with model collision integrals. Stability is k...A stable high-order Runge-Kutta discontinuous Galerkin(RKDG) scheme that strictly preserves positivity of the solution is designed to solve the Boltzmann kinetic equation with model collision integrals. Stability is kept by accuracy of velocity discretization, conservative calculation of the discrete collision relaxation term, and a limiter. By keeping the time step smaller than the local mean collision time and forcing positivity values of velocity distribution functions on certain points, the limiter can preserve positivity of solutions to the cell average velocity distribution functions. Verification is performed with a normal shock wave at a Mach number 2.05, a hypersonic flow about a two-dimensional(2D) cylinder at Mach numbers 6.0 and 12.0, and an unsteady shock tube flow. The results show that, the scheme is stable and accurate to capture shock structures in steady and unsteady hypersonic rarefied gaseous flows. Compared with two widely used limiters, the current limiter has the advantage of easy implementation and ability of minimizing the influence of accuracy of the original RKDG method.展开更多
Rock avalanches are catastrophic events involving important granular rock masses (>106 m3) and traveling long distances. In exceptional cases, the runout can reach up to tens of kilometers. Even if they are highly ...Rock avalanches are catastrophic events involving important granular rock masses (>106 m3) and traveling long distances. In exceptional cases, the runout can reach up to tens of kilometers. Even if they are highly destructive and uncontrollable events, they give important insights to understand interactions between the displaced masses and landscape conditions. However, those events are not frequent. Therefore, the analogue and numerical modelling gives fundamental inputs to better understand their behavior. The objective of the research is to understand the propagation and spreading of granular mass released at the top of a simple geometry. The flow is unconfined, spreading freely along a 45° slope and deposit on a horizontal surface. The evolution of this analogue rock avalanche was measured from the initiation to its deposition with high speed camera. To simulate the analogue granular flow, a numerical model based on the continuum mechanics approach and the solving of the shallow water equations was used. In this model, the avalanche is described from a eulerian point of view within a continuum framework as single phase of incompressible granular material. The interaction of the flowing layer with the substratum follows a Mohr-Coulomb friction law. Within same initial conditions (slope, volume, basal friction, height of fall and initial velocity), results obtained with the numerical model are similar to those observed in the analogue. In both cases, the runout of the mass is comparable and the size of both deposits matches well. Moreover, both analogue and numerical modeling gave same magnitude of velocities. In this study, we highlighted the importance of the friction on a flowing mass and the influence of the numerical resolution on the propagation. The combination of the fluid dynamic equation with the frictional law enables the self-channelization and the stop of the granular mass.展开更多
This paper presents nonlinear ordinary differential equations (ODES) of the heavier pellets movement for two phase flow, which actually represent a system of equations. The usual methods of solution such as Runge -Kut...This paper presents nonlinear ordinary differential equations (ODES) of the heavier pellets movement for two phase flow, which actually represent a system of equations. The usual methods of solution such as Runge -Kutta method and it's datum results are discussed. This paper solves ODES of general form using variable mesh-length, linearizing the nonlinear terms by finite analysis method, fuilding an iteration sequence, and amending the nonlinear terms by iteration . The conditions of convergent operation of iteration solution is checked. The movement orbit and velocity of the pellets are calculated. Analysis of research results and it's application examples are illustrated.展开更多
Omitting viscosity along flow direction, we have simplified the dimensionless N-Sequations in arbitrary curved coordinate system as the thin layer equations. Using theimplicit approximate-factorization algorithm to so...Omitting viscosity along flow direction, we have simplified the dimensionless N-Sequations in arbitrary curved coordinate system as the thin layer equations. Using theimplicit approximate-factorization algorithm to solve the gas-phase governing equ-ations and the characteristic method to follow the tracks of particles, we then obtainedthe full coupled numerical method of two-phase.transonic, turbulent flow. Here, par- ticle size may be grouped, the subsonic boundary condition at entry of nozzle is ireatedby quasi-characteristic method in reference plane and the algebraic model is used forturbulent flow. These methods are applied in viscous two-phase flow. calculation of ro-cket nozzle and in the prediciton of thrust and specific impulse for solid propellant ro-cket motor. The calculation results are in good agreement with the measurerment va-lues. Moreover, the influences of different particle radius, different particle mass frac-tion and particle size grouped on flow field have been discussed, and the influences of particle two-dimensional radial velosity component and viscosity on specific impulse ofrocket motor have been analysed.The method of this paper possesses the advantage of saving computer time. More important, the effect is more obvious for the calculation of particle size being grouped.展开更多
Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly im...Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly implicit alternating sweeping is implemented in the direction of the third dimension. Very rapid convergence rate is obtained with CFL number reaching the order of 100. The memory resources can be greatly saved too. It is verified that the reflection boundary condition can not be used with flux vector splitting since it will produce too large numerical dissipation. The computed flow fields agree well with experimental results. Only one or two grid points are there within the shock transition zone.展开更多
hi this paper. the author uses the theory of fluid mechanics. dynamics of fluids in Porous media. gas seepage flow in coal seams and combines the tests in the laboratory with the actual coal infusion to have an inves...hi this paper. the author uses the theory of fluid mechanics. dynamics of fluids in Porous media. gas seepage flow in coal seams and combines the tests in the laboratory with the actual coal infusion to have an investigating and study from the theory to the mechanism of coal infusion to wet coal seams. through the analysis to the process of coal infusion the author builds up the mathematical models and has a detailed discussion to the boundary conditions of coal infusion. Because the equation sets to describe coal infusion are non-linear. we have made a simplification to them to use the dimension analysis theory by leading into the non-dimensions of water pressure of coal infusion, seepage flow rate. increment of coal seam moisture and so on Besides the analytic and approximate solutions have also been discussed. At last. we use the scientific research item of the actual coal infusion to illustrate the effects and importance of the theory to direct actual coal infusion and its designs.展开更多
基金supported by the Yunnan Provincial Applied Basic Research Program of China(No. KKSY201207019)
文摘A three-dimensional (3D) predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.
基金supported by the National Key Technologies R&D Program of China for the Eleventh Five-Year Plan Period (Grant No. 2008BAB29B08-02)the Program for the Ministry of Education and State Administration of Foreign Experts Affairs of China (Grant No. B08408)
文摘Flow in tidal rivers periodically propagates upstream or downstream under tidal influence. Hydrodynamic models based on the Saint-Venant equations (the SVN model) are extensively used to model tidal rivers. A force-corrected term expressed as the combination of flow velocity and the change rate of the tidal fevel was developed to represent tidal effects in the SVN model. A momentum equation incorporating with the corrected term was derived based on Newton's second law. By combing the modified momentum equation with the continuity equation, an improved SVN model for tidal rivers (the ISVN model) was constructed. The simulation of a tidal reach of the Qiantang River shows that the ISVN model performs better than the SVN model. It indicates that the corrected force derived for tidal effects is reasonable; the ISVN model provides an appropriate enhancement of the SVN model for flow simulation of tidal rivers.
基金supported by Shanghai Science and Technology Innovation Action Plan(Grant No.21JC1403600)The work of the second author was partially supported by the National Natural Science Foundation of China(Grant No.11831011)Shanghai Science and Technology Innovation Action Plan(Grant No.21JC1403600).
文摘In this paper,we consider the two-dimensional aggregation equation with the shear flow and time-space nonlocal attractive operator.Without the advection,the solution of the aggregation equation may blow up in finite time.We show that the shear flow can suppress the blow-up.
文摘In this paper, the nonreflecting boundary conditions based upon fundamental ideas of the linear analysis are developed for gas dynamic equations, and the modified boundary conditions for Navier-Stokes equations are proposed as a substitute of the nonreflecting boundary conditions inside boundary layers near rigid walls. These derived boundary conditions are then applied to calculations both for the Euler equations and the Navier-Stokes equations to determine if they can produce acceptable results for the subsonic flows in channels. The numerical results obtained by an implicit second-order upwind difference scheme show the effective- ness and generality of the boundary conditions. Furthermore, the formulae and the analysis performed here may be extended to three dimensional problems.
文摘The short-range property of interactions between scales in incompressible turbulent flow was examined. Some formulae for the short-range eddy stress were given. A concept of resonant-range interactions between extremely contiguous scales was introduced and some formulae for the resonant-range eddy stress were also derived. Multi-scale equations for the incompressible turbulent flows were proposed. Key words turbulence - incompressible flow - interactions between scales - multi-scale equations MSC 2000 76F70
基金the two referees for very helpful comments and suggestions to improve the quality of the paper.This work was partially supported by the Natural Science Foundation of Zhejiang province of China(LY21A010017)the National Natural Science Foundation of China(12071106,12171130).
文摘This paper concerns the sonic-supersonic structures of the transonic crossflow generated by the steady supersonic flow past an infinite cone of arbitrary cross section.Under the conical assumption,the three-dimensional(3-D)steady Euler equations can be projected onto the unit sphere and the state of fluid can be characterized by the polar and azimuthal angles.Given a segment smooth curve as a conical-sonic line in the polar-azimuthal angle plane,we construct a classical conical-supersonic solution near the curve under some reasonable assumptions.To overcome the difficulty caused by the parabolic degeneracy,we apply the characteristic decomposition technique to transform the Euler equations into a new degenerate hyperbolic system in a partial hodograph plane.The singular terms are isolated from the highly nonlinear complicated system and then can be handled successfully.We establish a smooth local solution to the new system in a suitable weighted metric space and then express the solution in terms of the original variables.
基金supported by the National Key Research and Development Program of China(2020YFA0714300)the National Natural Science Foundation of China(62003084,62203108,62073079)+3 种基金the Natural Science Foundation of Jiangsu Province of China(BK20200355)the General Joint Fund of the Equipment Advance Research Program of Ministry of Education(8091B022114)Jiangsu Province Excellent Postdoctoral Program(2022ZB131)China Postdoctoral Science Foundation(2022M720720,2023T160105).
文摘The accelerated method in solving optimization problems has always been an absorbing topic.Based on the fixedtime(FxT)stability of nonlinear dynamical systems,we provide a unified approach for designing FxT gradient flows(FxTGFs).First,a general class of nonlinear functions in designing FxTGFs is provided.A unified method for designing first-order FxTGFs is shown under Polyak-Łjasiewicz inequality assumption,a weaker condition than strong convexity.When there exist both bounded and vanishing disturbances in the gradient flow,a specific class of nonsmooth robust FxTGFs with disturbance rejection is presented.Under the strict convexity assumption,Newton-based FxTGFs is given and further extended to solve time-varying optimization.Besides,the proposed FxTGFs are further used for solving equation-constrained optimization.Moreover,an FxT proximal gradient flow with a wide range of parameters is provided for solving nonsmooth composite optimization.To show the effectiveness of various FxTGFs,the static regret analyses for several typical FxTGFs are also provided in detail.Finally,the proposed FxTGFs are applied to solve two network problems,i.e.,the network consensus problem and solving a system linear equations,respectively,from the perspective of optimization.Particularly,by choosing component-wisely sign-preserving functions,these problems can be solved in a distributed way,which extends the existing results.The accelerated convergence and robustness of the proposed FxTGFs are validated in several numerical examples stemming from practical applications.
文摘For a stochastic differential equation with non-Lipschitz coefficients, we construct, by Euler scheme, a measurable flow of the solution, and we prove the solution is a Markov process.
基金supported by the NSF under Grant DMS-1818467Simons Foundation under Grant 961585.
文摘In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.
文摘The short-range property of interactions between scales in the compressible turbulent flow was examined. An estimation of the short-range scale scope and some formulae for the short-range eddy stress and heat transfer etc. were given. A concept of resonant-range interactions between extremely contiguous scales was introduced and some formulae for the resonant-range eddy stress and heat transfer etc. were also given. Multi-scale equations for the compressible turbulent flows were presented. The multi-scale equations are approximately closed and do not contain any empirical constants. The compressibility effects on turbulence are determined by the Farve averaged variables and the nonlinear relationships between the Farve- and physical-averaged variables.
基金Supported by the National Natural Science Foundation of China(11571361)China Scholarship Council
文摘A local gradient estimate for positive solutions of porous medium equations on complete noncompact Riemannian manifolds under the Ricci flow is derived. Moreover, a global gradient estimate for such equations on compact Riemannian manifolds is also obtained.
基金supported by the National Natural Science Foundation of China(Nos.11071057 and 11271052)the Special Fund Project of Mathematical Tian Yuan Fund(No.11226029)
文摘In this paper, the convergence turbulent flow equations are considered. By rates of solutions to the three-dimensional combining the LP-Lq estimate for the linearized equations and an elaborate energy method, the convergence rates are obtained in various norms for the solution to the equilibrium state in the whole space when the initial perturbation of the equilibrium state is small in the H3-framework. More precisely, the optimal convergence rates of the solutions and their first-order derivatives in the L2-norm are obtained when the LP-norm of the perturbation is bounded for some p ε [1, 6).
文摘In this paper, a new numerical scheme for solving first-order hyperbolic partial differential equations is proposed and is implemented in the simulation study of macroscopic traffic flow model with constant velocity and linear velocity-density relationship. Macroscopic traffic flow model is first developed by Lighthill Whitham and Richards (LWR) and used to study traffic flow by collective variables such as flow rate, velocity and density. The LWR model is treated as an initial value problem and its numerical simulations are presented using numerical schemes. A variety of numerical schemes are available in literature to solve first order hyperbolic equations. Of these the well-known ones include one-dimensional explicit: Upwind, Downwind, FTCS, and Lax-Friedrichs schemes. Having been studied carefully the space and time mesh sizes, and the patterns of all these schemes, a new scheme has been developed and named as one-dimensional explicit Tolesa numerical scheme. Tolesa numerical scheme is one of the conditionally stable and highest rates of convergence schemes. All the said numerical schemes are applied to solve advection equation pertaining traffic flows. Also the one-dimensional explicit Tolesa numerical scheme is another alternative numerical scheme to solve advection equation and apply to traffic flows model like other well-known one-dimensional explicit schemes. The effect of density of cars on the overall interactions of the vehicles along a given length of the highway and time are investigated. Graphical representations of density profile, velocity profile, flux profile, and in general the fundamental diagrams of vehicles on the highway with different time levels are illustrated. These concepts and results have been arranged systematically in this paper.
基金Project supported by the National Natural Science Foundation of China(No.11302017)
文摘A stable high-order Runge-Kutta discontinuous Galerkin(RKDG) scheme that strictly preserves positivity of the solution is designed to solve the Boltzmann kinetic equation with model collision integrals. Stability is kept by accuracy of velocity discretization, conservative calculation of the discrete collision relaxation term, and a limiter. By keeping the time step smaller than the local mean collision time and forcing positivity values of velocity distribution functions on certain points, the limiter can preserve positivity of solutions to the cell average velocity distribution functions. Verification is performed with a normal shock wave at a Mach number 2.05, a hypersonic flow about a two-dimensional(2D) cylinder at Mach numbers 6.0 and 12.0, and an unsteady shock tube flow. The results show that, the scheme is stable and accurate to capture shock structures in steady and unsteady hypersonic rarefied gaseous flows. Compared with two widely used limiters, the current limiter has the advantage of easy implementation and ability of minimizing the influence of accuracy of the original RKDG method.
基金Foundation item: Projects(60835005, 90820302) supported by the National Natural Science Foundation of China Project(2007CB311001) supported by the National Basic Research Program of China
文摘Rock avalanches are catastrophic events involving important granular rock masses (>106 m3) and traveling long distances. In exceptional cases, the runout can reach up to tens of kilometers. Even if they are highly destructive and uncontrollable events, they give important insights to understand interactions between the displaced masses and landscape conditions. However, those events are not frequent. Therefore, the analogue and numerical modelling gives fundamental inputs to better understand their behavior. The objective of the research is to understand the propagation and spreading of granular mass released at the top of a simple geometry. The flow is unconfined, spreading freely along a 45° slope and deposit on a horizontal surface. The evolution of this analogue rock avalanche was measured from the initiation to its deposition with high speed camera. To simulate the analogue granular flow, a numerical model based on the continuum mechanics approach and the solving of the shallow water equations was used. In this model, the avalanche is described from a eulerian point of view within a continuum framework as single phase of incompressible granular material. The interaction of the flowing layer with the substratum follows a Mohr-Coulomb friction law. Within same initial conditions (slope, volume, basal friction, height of fall and initial velocity), results obtained with the numerical model are similar to those observed in the analogue. In both cases, the runout of the mass is comparable and the size of both deposits matches well. Moreover, both analogue and numerical modeling gave same magnitude of velocities. In this study, we highlighted the importance of the friction on a flowing mass and the influence of the numerical resolution on the propagation. The combination of the fluid dynamic equation with the frictional law enables the self-channelization and the stop of the granular mass.
文摘This paper presents nonlinear ordinary differential equations (ODES) of the heavier pellets movement for two phase flow, which actually represent a system of equations. The usual methods of solution such as Runge -Kutta method and it's datum results are discussed. This paper solves ODES of general form using variable mesh-length, linearizing the nonlinear terms by finite analysis method, fuilding an iteration sequence, and amending the nonlinear terms by iteration . The conditions of convergent operation of iteration solution is checked. The movement orbit and velocity of the pellets are calculated. Analysis of research results and it's application examples are illustrated.
文摘Omitting viscosity along flow direction, we have simplified the dimensionless N-Sequations in arbitrary curved coordinate system as the thin layer equations. Using theimplicit approximate-factorization algorithm to solve the gas-phase governing equ-ations and the characteristic method to follow the tracks of particles, we then obtainedthe full coupled numerical method of two-phase.transonic, turbulent flow. Here, par- ticle size may be grouped, the subsonic boundary condition at entry of nozzle is ireatedby quasi-characteristic method in reference plane and the algebraic model is used forturbulent flow. These methods are applied in viscous two-phase flow. calculation of ro-cket nozzle and in the prediciton of thrust and specific impulse for solid propellant ro-cket motor. The calculation results are in good agreement with the measurerment va-lues. Moreover, the influences of different particle radius, different particle mass frac-tion and particle size grouped on flow field have been discussed, and the influences of particle two-dimensional radial velosity component and viscosity on specific impulse ofrocket motor have been analysed.The method of this paper possesses the advantage of saving computer time. More important, the effect is more obvious for the calculation of particle size being grouped.
文摘Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly implicit alternating sweeping is implemented in the direction of the third dimension. Very rapid convergence rate is obtained with CFL number reaching the order of 100. The memory resources can be greatly saved too. It is verified that the reflection boundary condition can not be used with flux vector splitting since it will produce too large numerical dissipation. The computed flow fields agree well with experimental results. Only one or two grid points are there within the shock transition zone.
文摘hi this paper. the author uses the theory of fluid mechanics. dynamics of fluids in Porous media. gas seepage flow in coal seams and combines the tests in the laboratory with the actual coal infusion to have an investigating and study from the theory to the mechanism of coal infusion to wet coal seams. through the analysis to the process of coal infusion the author builds up the mathematical models and has a detailed discussion to the boundary conditions of coal infusion. Because the equation sets to describe coal infusion are non-linear. we have made a simplification to them to use the dimension analysis theory by leading into the non-dimensions of water pressure of coal infusion, seepage flow rate. increment of coal seam moisture and so on Besides the analytic and approximate solutions have also been discussed. At last. we use the scientific research item of the actual coal infusion to illustrate the effects and importance of the theory to direct actual coal infusion and its designs.