We consider the singular Riemann problem for the rectilinear isentropic compressible Euler equations with discontinuous flux,more specifically,for pressureless flow on the left and polytropic flow on the right separat...We consider the singular Riemann problem for the rectilinear isentropic compressible Euler equations with discontinuous flux,more specifically,for pressureless flow on the left and polytropic flow on the right separated by a discontinuity x=x(t).We prove that this problem admits global Radon measure solutions for all kinds of initial data.The over-compressing condition on the discontinuity x=x(t)is not enough to ensure the uniqueness of the solution.However,there is a unique piecewise smooth solution if one proposes a slip condition on the right-side of the curve x=x(t)+0,in addition to the full adhesion condition on its left-side.As an application,we study a free piston problem with the piston in a tube surrounded initially by uniform pressureless flow and a polytropic gas.In particular,we obtain the existence of a piecewise smooth solution for the motion of the piston between a vacuum and a polytropic gas.This indicates that the singular Riemann problem looks like a control problem in the sense that one could adjust the condition on the discontinuity of the flux to obtain the desired flow field.展开更多
We solve the Riemann problems for isentropic compressible Euler equations of polytropic gases in the class of Radon measures,and the solutions admit the concentration of mass.It is found that under the requirement of ...We solve the Riemann problems for isentropic compressible Euler equations of polytropic gases in the class of Radon measures,and the solutions admit the concentration of mass.It is found that under the requirement of satisfying the over-compressing entropy condition:(i)there is a unique delta shock solution,corresponding to the case that has two strong classical Lax shocks;(ii)for the initial data that the classical Riemann solution contains a shock wave and a rarefaction wave,or two shocks with one being weak,there are infinitely many solutions,each consists of a delta shock and a rarefaction wave;(iii)there are no delta shocks for the case that the classical entropy weak solutions consist only of rarefaction waves.These solutions are self-similar.Furthermore,for the generalized Riemann problem with mass concentrated initially at the discontinuous point of initial data,there always exists a unique delta shock for at least a short time.It could be prolonged to a global solution.Not all the solutions are self-similar due to the initial velocity of the concentrated point-mass(particle).Whether the delta shock solutions constructed satisfy the over-compressing entropy condition is clarified.This is the first result on the construction of singular measure solutions to the compressible Euler system of polytropic gases,that is strictly hyperbolic,and whose characteristics are both genuinely nonlinear.We also discuss possible physical interpretations and applications of these new 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.展开更多
A discrete Boltzmann model(DBM) with symmetric velocity discretization is constructed for compressible systems with an adjustable specific heat ratio in the external force field. The proposed two-dimensional(2D) nine-...A discrete Boltzmann model(DBM) with symmetric velocity discretization is constructed for compressible systems with an adjustable specific heat ratio in the external force field. The proposed two-dimensional(2D) nine-velocity scheme has better spatial symmetry and numerical accuracy than the discretized velocity model in literature [Acta Aerodyn. Sin.40 98108(2022)] and owns higher computational efficiency than the one in literature [Phys. Rev. E 99 012142(2019)].In addition, the matrix inversion method is adopted to calculate the discrete equilibrium distribution function and force term, both of which satisfy nine independent kinetic moment relations. Moreover, the DBM could be used to study a few thermodynamic nonequilibrium effects beyond the Euler equations that are recovered from the kinetic model in the hydrodynamic limit via the Chapman–Enskog expansion. Finally, the present method is verified through typical numerical simulations, including the free-falling process, Sod’s shock tube, sound wave, compressible Rayleigh–Taylor instability,and translational motion of a 2D fluid system.展开更多
An additional potential energy distribution function is introduced on the basis of previous D3Q25 model,and the equilibrium distribution function of D3Q25 is obtained by spherical function.A novel three-dimensional(3D...An additional potential energy distribution function is introduced on the basis of previous D3Q25 model,and the equilibrium distribution function of D3Q25 is obtained by spherical function.A novel three-dimensional(3D)shifted lattice model is proposed,therefore a shifted lattice model is introduced into D3Q25.Under the finite volume scheme,several typical compressible calculation examples are used to verify whether the numerical stability of the D3Q25 model can be improved by adding the shifted lattice model.The simulation results show that the numerical stability is indeed improved after adding the shifted lattice model.展开更多
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.展开更多
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.展开更多
A high-order gas kinetic flux solver(GKFS)is presented for simulating inviscid compressible flows.The weighted essentially non-oscillatory(WENO)scheme on a uniform mesh in the finite volume formulation is combined wit...A high-order gas kinetic flux solver(GKFS)is presented for simulating inviscid compressible flows.The weighted essentially non-oscillatory(WENO)scheme on a uniform mesh in the finite volume formulation is combined with the circular function-based GKFS(C-GKFS)to capture more details of the flow fields with fewer grids.Different from most of the current GKFSs,which are constructed based on the Maxwellian distribution function or its equivalent form,the C-GKFS simplifies the Maxwellian distribution function into the circular function,which ensures that the Euler or Navier-Stokes equations can be recovered correctly.This improves the efficiency of the GKFS and reduces its complexity to facilitate the practical application of engineering.Several benchmark cases are simulated,and good agreement can be obtained in comparison with the references,which demonstrates that the high-order C-GKFS can achieve the desired accuracy.展开更多
Fully connected neural networks(FCNNs)have been developed for the closure of subgrid-scale(SGS)stress and SGS heat flux in large-eddy simulations of compressible turbulent channel flow.The FCNNbased SGS model trained ...Fully connected neural networks(FCNNs)have been developed for the closure of subgrid-scale(SGS)stress and SGS heat flux in large-eddy simulations of compressible turbulent channel flow.The FCNNbased SGS model trained using data with Mach number Ma=3.0 and Reynolds number Re=3000 was applied to situations with different Mach numbers and Reynolds numbers.The input variables of the neural network model were the filtered velocity gradients and temperature gradients at a single spatial grid point.The a priori test showed that the FCNN model had a correlation coefficient larger than 0.91 and a relative error smaller than 0.43,with much better reconstructions of SGS unclosed terms than the dynamic Smagorinsky model(DSM).In a posteriori test,the behavior of the FCNN model was marginally better than that of the DSM in predicting the mean velocity profiles,mean temperature profiles,turbulent intensities,total Reynolds stress,total Reynolds heat flux,and mean SGS flux of kinetic energy,and outperformed the Smagorinsky model.展开更多
The numerical solution of compressible flows has become more prevalent than that of incompressible flows.With the help of the artificial compressibility approach,incompressible flows can be solved numerically using th...The numerical solution of compressible flows has become more prevalent than that of incompressible flows.With the help of the artificial compressibility approach,incompressible flows can be solved numerically using the same methods as compressible ones.The artificial compressibility scheme is thus widely used to numerically solve incompressible Navier-Stokes equations.Any numerical method highly depends on its accuracy and speed of convergence.Although the artificial compressibility approach is utilized in several numerical simulations,the effect of the compressibility factor on the accuracy of results and convergence speed has not been investigated for nanofluid flows in previous studies.Therefore,this paper assesses the effect of this factor on the convergence speed and accuracy of results for various types of thermo-flow.To improve the stability and convergence speed of time discretizations,the fifth-order Runge-Kutta method is applied.A computer program has been written in FORTRAN to solve the discretized equations in different Reynolds and Grashof numbers for various grids.The results demonstrate that the artificial compressibility factor has a noticeable effect on the accuracy and convergence rate of the simulation.The optimum artificial compressibility is found to be between 1 and 5.These findings can be utilized to enhance the performance of commercial numerical simulation tools,including ANSYS and COMSOL.展开更多
Inventors:Xiujun Cai,Mingyu Chen,Chen Lu,Yifan Wang,Diyu Huang,Hepan Zhu,Yibin Zhu,Bin Zhang Applicant:Zhejiang University,Zhejiang,China Patent No.:US 10,980,541 B2 Date of Patent:Apr.20,2021 An absorbable and unidir...Inventors:Xiujun Cai,Mingyu Chen,Chen Lu,Yifan Wang,Diyu Huang,Hepan Zhu,Yibin Zhu,Bin Zhang Applicant:Zhejiang University,Zhejiang,China Patent No.:US 10,980,541 B2 Date of Patent:Apr.20,2021 An absorbable and unidirectionally compressible intestinal anastomosis device includes a first through pipe portion and a second through pipe portion including large ends and small ends,respectively;the large ends are used for fixing an intestinal canal,and the small end of the second through pipe portion is inserted into the small end of the first through pipe portion.The two large ends are drawn close to each other to achieve an intestineeintestine anastomosis.The outer sides of circular rings at the ends of the small ends of the first through pipe portion and the second through pipe portion are provided with gradually expanded elastic horn-shaped openings.The device can make two stapled intestinal canals rest relatively,thereby avoiding adjustments required in other assembly modes.展开更多
This paper mainly studies the blowup phenomenon of solutions to the compressible Euler equations with general time-dependent damping for non-isentropic fluids in two and three space dimensions. When the initial data i...This paper mainly studies the blowup phenomenon of solutions to the compressible Euler equations with general time-dependent damping for non-isentropic fluids in two and three space dimensions. When the initial data is assumed to be radially symmetric and the initial density contains vacuum, we obtain that classical solution, especially the density, will blow up on finite time. The results also reveal that damping can really delay the singularity formation.展开更多
High-speed rotor rotation under the low-density condition creates a special low-Reynolds compressible flow around the rotor blade airfoil where the compressibility effect on the laminar separated shear layer occurs. H...High-speed rotor rotation under the low-density condition creates a special low-Reynolds compressible flow around the rotor blade airfoil where the compressibility effect on the laminar separated shear layer occurs. However, the compressibility effect and shock wave generation associated with the increase in the Mach number (M) and the trend change due to their interference have not been clarified. The purpose is to clear the compressibility effect and its impact of shock wave generation on the flow field and aerodynamics. Therefore, we perform a two-dimensional unsteady calculation by Computational fluid dynamics (CFD) analysis using the CLF5605 airfoil used in the Mars helicopter Ingenuity, which succeeded in its first flight on Mars. The calculation conditions are set to the Reynolds number (Re) at 75% rotor span in hovering (Re = 15,400), and the Mach number was varied from incompressible (M = 0.2) to transonic (M = 1.2). The compressible fluid dynamics solver FaSTAR developed by the Japan aerospace exploration agency (JAXA) is used, and calculations are performed under multiple conditions in which the Mach number and angle of attack (α) are swept. The results show that a flow field is similar to that in the Earth’s atmosphere above M = 1.0, such as bow shock at the leading edge, whereas multiple λ-type shock waves are observed over the separated shear layer above α = 3° at M = 0.80. However, no significant difference is found in the C<sub>p</sub> distribution around the airfoil between M = 0.6 and M = 0.8. From the results, it is found that multiple λ-type shock waves have no significant effect on the airfoil surface pressure distribution, the separated shear layer effect is dominant in the surface pressure change and aerodynamic characteristics.展开更多
We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusi...We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusing on discontinuous spectral element semidis-cretizations,we design new controllers for existing methods and for some new embedded Runge-Kutta pairs.We demonstrate the importance of choosing adequate controller parameters and provide a means to obtain these in practice.We compare a wide range of error-control-based methods,along with the common approach in which step size con-trol is based on the Courant-Friedrichs-Lewy(CFL)number.The optimized methods give improved performance and naturally adopt a step size close to the maximum stable CFL number at loose tolerances,while additionally providing control of the temporal error at tighter tolerances.The numerical examples include challenging industrial CFD applications.展开更多
The cavitated bifurcation problem in a solid sphere composed of two compressible hyper-elas-tic materials is examined.The bifurcation solution for the composed sphere under a uniform radial tensileboundary dead-load i...The cavitated bifurcation problem in a solid sphere composed of two compressible hyper-elas-tic materials is examined.The bifurcation solution for the composed sphere under a uniform radial tensileboundary dead-load is obtained.The bifurcation curves and the stress contributions subsequent to the cavita-tion are given.The right and left bifurcation as well as the catastrophe and concentration of stresses are ana-lyzed.The stability of solutions is discussed through an energy comparison.展开更多
A new lattice Boltzmann model for compressible perfect gas is pro-posed.The numerical example shows that it can be used to simulate shock wave andcontact discontinuity.The results are comparable with those obtained by...A new lattice Boltzmann model for compressible perfect gas is pro-posed.The numerical example shows that it can be used to simulate shock wave andcontact discontinuity.The results are comparable with those obtained by traditionalmethods.The ratio of specific heats γmay be chosen according to the requirementof problems.展开更多
By analyzing the characteristics of low Mach number perfect gasflows, a novel Slightly Compressible Model (SCM) for low Mach number perfect gasflows is derived. In view of numerical calculations, this model is proved ...By analyzing the characteristics of low Mach number perfect gasflows, a novel Slightly Compressible Model (SCM) for low Mach number perfect gasflows is derived. In view of numerical calculations, this model is proved very efficient,for it is kept within the p-v frame but does not have to satisfy the time consumingdivergence-free condition in order to get the incompressible Navier-Stokes equationsolutions. Writing the equations in the form of conservation laws, we have derivedthe characteristic systems which are necessary for numerical calculations. A cell-centered finite-volume method with flux difference upwind-biased schemes is used forthe equation solutions and a new Exact Newton Relaxation (ENR) implicit methodis developed. Various computed results are presented to validate the present model.Laminar flow solutions over a circular cylinder with wake developing and vortex shed-ding are presented. Results for inviscid flow over a sphere are compared in excellentagreement with the exact analytic incompressible solution. Three-dimensional viscousflow solutions over sphere and prolate spheroid are also calculated and compared wellwith experiments and other incompressible solutions. Finally, good convergent per-formaces are shown for sphere viscous flows.展开更多
To simulate the transonic atomization jet process in Laval nozzles,to test the law of droplet atomization and distribution,to find a method of supersonic atomization for dust-removing nozzles,and to improve nozzle eff...To simulate the transonic atomization jet process in Laval nozzles,to test the law of droplet atomization and distribution,to find a method of supersonic atomization for dust-removing nozzles,and to improve nozzle efficiency,the finite element method has been used in this study based on the COMSOL computational fluid dynamics module.The study results showed that the process cannot be realized alone under the two-dimensional axisymmetric,three-dimensional and three-dimensional symmetric models,but it can be calculated with the transformation dimension method,which uses the parameter equations generated from the two-dimensional axisymmetric flow field data of the three-dimensional model.The visualization of this complex process,which is difficult to measure and analyze experimentally,was realized in this study.The physical process,macro phenomena and particle distribution of supersonic atomization are analyzed in combination with this simulation.The rationality of the simulation was verified by experiments.A new method for the study of the atomization process and the exploration of its mechanism in a compressible transonic speed flow field based on the Laval nozzle has been provided,and a numerical platform for the study of supersonic atomization dust removal has been established.展开更多
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.展开更多
基金supported by the National Natural Science Foundation of China(11871218,12071298)in part by the Science and Technology Commission of Shanghai Municipality(21JC1402500,22DZ2229014)。
文摘We consider the singular Riemann problem for the rectilinear isentropic compressible Euler equations with discontinuous flux,more specifically,for pressureless flow on the left and polytropic flow on the right separated by a discontinuity x=x(t).We prove that this problem admits global Radon measure solutions for all kinds of initial data.The over-compressing condition on the discontinuity x=x(t)is not enough to ensure the uniqueness of the solution.However,there is a unique piecewise smooth solution if one proposes a slip condition on the right-side of the curve x=x(t)+0,in addition to the full adhesion condition on its left-side.As an application,we study a free piston problem with the piston in a tube surrounded initially by uniform pressureless flow and a polytropic gas.In particular,we obtain the existence of a piecewise smooth solution for the motion of the piston between a vacuum and a polytropic gas.This indicates that the singular Riemann problem looks like a control problem in the sense that one could adjust the condition on the discontinuity of the flux to obtain the desired flow field.
基金supported by the National Natural Science Foundation of China under Grants No.11871218,No.12071298the Science and Technology Commission of Shanghai Municipality under Grant No.18dz2271000.
文摘We solve the Riemann problems for isentropic compressible Euler equations of polytropic gases in the class of Radon measures,and the solutions admit the concentration of mass.It is found that under the requirement of satisfying the over-compressing entropy condition:(i)there is a unique delta shock solution,corresponding to the case that has two strong classical Lax shocks;(ii)for the initial data that the classical Riemann solution contains a shock wave and a rarefaction wave,or two shocks with one being weak,there are infinitely many solutions,each consists of a delta shock and a rarefaction wave;(iii)there are no delta shocks for the case that the classical entropy weak solutions consist only of rarefaction waves.These solutions are self-similar.Furthermore,for the generalized Riemann problem with mass concentrated initially at the discontinuous point of initial data,there always exists a unique delta shock for at least a short time.It could be prolonged to a global solution.Not all the solutions are self-similar due to the initial velocity of the concentrated point-mass(particle).Whether the delta shock solutions constructed satisfy the over-compressing entropy condition is clarified.This is the first result on the construction of singular measure solutions to the compressible Euler system of polytropic gases,that is strictly hyperbolic,and whose characteristics are both genuinely nonlinear.We also discuss possible physical interpretations and applications of these new 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.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 51806116, U2242214, and 11875329)Guangdong Basic and Applied Basic Research Foundation (Grant No. 2022A1515012116)the Natural Science Foundation of Fujian Province, China (Grant Nos. 2021J01652 and 2021J01655)。
文摘A discrete Boltzmann model(DBM) with symmetric velocity discretization is constructed for compressible systems with an adjustable specific heat ratio in the external force field. The proposed two-dimensional(2D) nine-velocity scheme has better spatial symmetry and numerical accuracy than the discretized velocity model in literature [Acta Aerodyn. Sin.40 98108(2022)] and owns higher computational efficiency than the one in literature [Phys. Rev. E 99 012142(2019)].In addition, the matrix inversion method is adopted to calculate the discrete equilibrium distribution function and force term, both of which satisfy nine independent kinetic moment relations. Moreover, the DBM could be used to study a few thermodynamic nonequilibrium effects beyond the Euler equations that are recovered from the kinetic model in the hydrodynamic limit via the Chapman–Enskog expansion. Finally, the present method is verified through typical numerical simulations, including the free-falling process, Sod’s shock tube, sound wave, compressible Rayleigh–Taylor instability,and translational motion of a 2D fluid system.
基金the Youth Program of the National Natural Science Foundation of China(Grant Nos.11972272,12072246,and 12202331)the National Key Project,China(Grant No.GJXM92579)the Natural Science Basic Research Program of Shaanxi Province,China(Program No.2022JQ-028)。
文摘An additional potential energy distribution function is introduced on the basis of previous D3Q25 model,and the equilibrium distribution function of D3Q25 is obtained by spherical function.A novel three-dimensional(3D)shifted lattice model is proposed,therefore a shifted lattice model is introduced into D3Q25.Under the finite volume scheme,several typical compressible calculation examples are used to verify whether the numerical stability of the D3Q25 model can be improved by adding the shifted lattice model.The simulation results show that the numerical stability is indeed improved after adding the shifted lattice model.
基金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.
基金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.
基金Project supported by the National Natural Science Foundation of China(No.12072158)。
文摘A high-order gas kinetic flux solver(GKFS)is presented for simulating inviscid compressible flows.The weighted essentially non-oscillatory(WENO)scheme on a uniform mesh in the finite volume formulation is combined with the circular function-based GKFS(C-GKFS)to capture more details of the flow fields with fewer grids.Different from most of the current GKFSs,which are constructed based on the Maxwellian distribution function or its equivalent form,the C-GKFS simplifies the Maxwellian distribution function into the circular function,which ensures that the Euler or Navier-Stokes equations can be recovered correctly.This improves the efficiency of the GKFS and reduces its complexity to facilitate the practical application of engineering.Several benchmark cases are simulated,and good agreement can be obtained in comparison with the references,which demonstrates that the high-order C-GKFS can achieve the desired accuracy.
基金Financial support provided by the National Natural Science Foundation of China(Grant Nos.11702042 and 91952104)。
文摘Fully connected neural networks(FCNNs)have been developed for the closure of subgrid-scale(SGS)stress and SGS heat flux in large-eddy simulations of compressible turbulent channel flow.The FCNNbased SGS model trained using data with Mach number Ma=3.0 and Reynolds number Re=3000 was applied to situations with different Mach numbers and Reynolds numbers.The input variables of the neural network model were the filtered velocity gradients and temperature gradients at a single spatial grid point.The a priori test showed that the FCNN model had a correlation coefficient larger than 0.91 and a relative error smaller than 0.43,with much better reconstructions of SGS unclosed terms than the dynamic Smagorinsky model(DSM).In a posteriori test,the behavior of the FCNN model was marginally better than that of the DSM in predicting the mean velocity profiles,mean temperature profiles,turbulent intensities,total Reynolds stress,total Reynolds heat flux,and mean SGS flux of kinetic energy,and outperformed the Smagorinsky model.
基金The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through the Large Groups Project under grant number RGP.2/235/43.
文摘The numerical solution of compressible flows has become more prevalent than that of incompressible flows.With the help of the artificial compressibility approach,incompressible flows can be solved numerically using the same methods as compressible ones.The artificial compressibility scheme is thus widely used to numerically solve incompressible Navier-Stokes equations.Any numerical method highly depends on its accuracy and speed of convergence.Although the artificial compressibility approach is utilized in several numerical simulations,the effect of the compressibility factor on the accuracy of results and convergence speed has not been investigated for nanofluid flows in previous studies.Therefore,this paper assesses the effect of this factor on the convergence speed and accuracy of results for various types of thermo-flow.To improve the stability and convergence speed of time discretizations,the fifth-order Runge-Kutta method is applied.A computer program has been written in FORTRAN to solve the discretized equations in different Reynolds and Grashof numbers for various grids.The results demonstrate that the artificial compressibility factor has a noticeable effect on the accuracy and convergence rate of the simulation.The optimum artificial compressibility is found to be between 1 and 5.These findings can be utilized to enhance the performance of commercial numerical simulation tools,including ANSYS and COMSOL.
文摘Inventors:Xiujun Cai,Mingyu Chen,Chen Lu,Yifan Wang,Diyu Huang,Hepan Zhu,Yibin Zhu,Bin Zhang Applicant:Zhejiang University,Zhejiang,China Patent No.:US 10,980,541 B2 Date of Patent:Apr.20,2021 An absorbable and unidirectionally compressible intestinal anastomosis device includes a first through pipe portion and a second through pipe portion including large ends and small ends,respectively;the large ends are used for fixing an intestinal canal,and the small end of the second through pipe portion is inserted into the small end of the first through pipe portion.The two large ends are drawn close to each other to achieve an intestineeintestine anastomosis.The outer sides of circular rings at the ends of the small ends of the first through pipe portion and the second through pipe portion are provided with gradually expanded elastic horn-shaped openings.The device can make two stapled intestinal canals rest relatively,thereby avoiding adjustments required in other assembly modes.
文摘This paper mainly studies the blowup phenomenon of solutions to the compressible Euler equations with general time-dependent damping for non-isentropic fluids in two and three space dimensions. When the initial data is assumed to be radially symmetric and the initial density contains vacuum, we obtain that classical solution, especially the density, will blow up on finite time. The results also reveal that damping can really delay the singularity formation.
文摘High-speed rotor rotation under the low-density condition creates a special low-Reynolds compressible flow around the rotor blade airfoil where the compressibility effect on the laminar separated shear layer occurs. However, the compressibility effect and shock wave generation associated with the increase in the Mach number (M) and the trend change due to their interference have not been clarified. The purpose is to clear the compressibility effect and its impact of shock wave generation on the flow field and aerodynamics. Therefore, we perform a two-dimensional unsteady calculation by Computational fluid dynamics (CFD) analysis using the CLF5605 airfoil used in the Mars helicopter Ingenuity, which succeeded in its first flight on Mars. The calculation conditions are set to the Reynolds number (Re) at 75% rotor span in hovering (Re = 15,400), and the Mach number was varied from incompressible (M = 0.2) to transonic (M = 1.2). The compressible fluid dynamics solver FaSTAR developed by the Japan aerospace exploration agency (JAXA) is used, and calculations are performed under multiple conditions in which the Mach number and angle of attack (α) are swept. The results show that a flow field is similar to that in the Earth’s atmosphere above M = 1.0, such as bow shock at the leading edge, whereas multiple λ-type shock waves are observed over the separated shear layer above α = 3° at M = 0.80. However, no significant difference is found in the C<sub>p</sub> distribution around the airfoil between M = 0.6 and M = 0.8. From the results, it is found that multiple λ-type shock waves have no significant effect on the airfoil surface pressure distribution, the separated shear layer effect is dominant in the surface pressure change and aerodynamic characteristics.
基金Open Access funding enabled and organized by Projekt DEAL.
文摘We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusing on discontinuous spectral element semidis-cretizations,we design new controllers for existing methods and for some new embedded Runge-Kutta pairs.We demonstrate the importance of choosing adequate controller parameters and provide a means to obtain these in practice.We compare a wide range of error-control-based methods,along with the common approach in which step size con-trol is based on the Courant-Friedrichs-Lewy(CFL)number.The optimized methods give improved performance and naturally adopt a step size close to the maximum stable CFL number at loose tolerances,while additionally providing control of the temporal error at tighter tolerances.The numerical examples include challenging industrial CFD applications.
基金the National Natttral Science Foundation of China(No.19802012)
文摘The cavitated bifurcation problem in a solid sphere composed of two compressible hyper-elas-tic materials is examined.The bifurcation solution for the composed sphere under a uniform radial tensileboundary dead-load is obtained.The bifurcation curves and the stress contributions subsequent to the cavita-tion are given.The right and left bifurcation as well as the catastrophe and concentration of stresses are ana-lyzed.The stability of solutions is discussed through an energy comparison.
基金The project supported by the National Natural Science Foundation of China
文摘A new lattice Boltzmann model for compressible perfect gas is pro-posed.The numerical example shows that it can be used to simulate shock wave andcontact discontinuity.The results are comparable with those obtained by traditionalmethods.The ratio of specific heats γmay be chosen according to the requirementof problems.
基金The project supported by the Basic Research on Frontier Problems in Fluid and Aerodynamics in Chinathe National Natural Science Foundation of China (19772069)
文摘By analyzing the characteristics of low Mach number perfect gasflows, a novel Slightly Compressible Model (SCM) for low Mach number perfect gasflows is derived. In view of numerical calculations, this model is proved very efficient,for it is kept within the p-v frame but does not have to satisfy the time consumingdivergence-free condition in order to get the incompressible Navier-Stokes equationsolutions. Writing the equations in the form of conservation laws, we have derivedthe characteristic systems which are necessary for numerical calculations. A cell-centered finite-volume method with flux difference upwind-biased schemes is used forthe equation solutions and a new Exact Newton Relaxation (ENR) implicit methodis developed. Various computed results are presented to validate the present model.Laminar flow solutions over a circular cylinder with wake developing and vortex shed-ding are presented. Results for inviscid flow over a sphere are compared in excellentagreement with the exact analytic incompressible solution. Three-dimensional viscousflow solutions over sphere and prolate spheroid are also calculated and compared wellwith experiments and other incompressible solutions. Finally, good convergent per-formaces are shown for sphere viscous flows.
基金Supported by the National Natural Science Foundation of China (NO: 51704146, 51274116, 51704145).
文摘To simulate the transonic atomization jet process in Laval nozzles,to test the law of droplet atomization and distribution,to find a method of supersonic atomization for dust-removing nozzles,and to improve nozzle efficiency,the finite element method has been used in this study based on the COMSOL computational fluid dynamics module.The study results showed that the process cannot be realized alone under the two-dimensional axisymmetric,three-dimensional and three-dimensional symmetric models,but it can be calculated with the transformation dimension method,which uses the parameter equations generated from the two-dimensional axisymmetric flow field data of the three-dimensional model.The visualization of this complex process,which is difficult to measure and analyze experimentally,was realized in this study.The physical process,macro phenomena and particle distribution of supersonic atomization are analyzed in combination with this simulation.The rationality of the simulation was verified by experiments.A new method for the study of the atomization process and the exploration of its mechanism in a compressible transonic speed flow field based on the Laval nozzle has been provided,and a numerical platform for the study of supersonic atomization dust removal has been established.
基金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.
