Polyurethane foam,when used as a compressible layer in deep soft rock tunnels,offers a feasible solution to reduce the support pressure on the secondary lining.The foam spraying method using sprayed polyurethane mater...Polyurethane foam,when used as a compressible layer in deep soft rock tunnels,offers a feasible solution to reduce the support pressure on the secondary lining.The foam spraying method using sprayed polyurethane material is convenient for engineering applications;however,the compressive behaviour and feasibility of sprayed polyurethane material as a compressible layer remain unclear.To address this gap,this study conducts uniaxial compression tests and scanning electron microscope(SEM)tests to investigate the compressive behaviour of the rigid foams fabricated from a self-developed polyurethane spray material.A peridynamics model for the composite lining with a polyurethane compressible layer is then established.After validating the proposed method by comparison with two tests,a parametric study is carried out to investigate the damage evolution of the composite lining with a polyurethane compressible layer under various combinations of large deformations and compressible layer parameters.The results indicate that the polyurethane compressible layer effectively reduces the radial deformation and damage index of the secondary lining while increasing the damage susceptibility of the primary lining.The thickness of the polyurethane compressible layer significantly influences the prevention effect of large deformation-induced damage to the secondary lining within the density range of 50e100 kg/m^(3).In accordance with the experimental and simulation results,a simple,yet reasonable and convenient approach for determining the key parameters of the polyurethane compressible layer is proposed,along with a classification scheme for the parameters of the polyurethane compressible layer.展开更多
Numerical simulations were conducted on a 10-blade Sevik rotor ingesting wake downstream of two turbulence-generating grids.These simulations were based on implicit large-eddy simulation(ILES)and the boundary data imm...Numerical simulations were conducted on a 10-blade Sevik rotor ingesting wake downstream of two turbulence-generating grids.These simulations were based on implicit large-eddy simulation(ILES)and the boundary data immersion method(BDIM)for compressible flows,which were solved using a fully self-programmed Fortran code.Results show that the predicted thrust spectrum aligns closely with the experimental measurements.In addition,it captures the thrust dipole directivity of the noise around the rotating propeller due to random pressure pulsations on the blades,as well as the flow structures simultaneously.Furthermore,the differences in the statistical characteristics,flow structures,and low-frequency broadband thrust spectra due to different turbulence levels were investigated.This analysis indicates that the interaction between the upstream,which is characterized by a lower turbulence level and a higher turbulent length of scale,and the rotating propeller results in a lower amplitude in force spectra and a slight increase in the scale of tip vortices.展开更多
This paper concerns the Cauchy problem of 3D compressible micropolar fluids in the whole space R^(3). For regular initial data with m0E0 is suitable small, where m0 and E0 represent the upper bound of initial density ...This paper concerns the Cauchy problem of 3D compressible micropolar fluids in the whole space R^(3). For regular initial data with m0E0 is suitable small, where m0 and E0 represent the upper bound of initial density and initial energy, we prove that if ρ0 ∈ Lγ ∩ H3 with γ ∈ (1, 6), then the problem possesses a unique global classical solution on R^(3) × [0, T] with any T ∈ (0, ∞). It’s worth noting that both the vacuum states and possible random largeness of initial energy are allowed.展开更多
In this paper,we develop a fourth-order conservative wavelet-based shock-capturing scheme.The scheme is constructed by combining a wavelet collocation upwind method with the monotonic tangent of hyperbola for interfac...In this paper,we develop a fourth-order conservative wavelet-based shock-capturing scheme.The scheme is constructed by combining a wavelet collocation upwind method with the monotonic tangent of hyperbola for interface capturing(THINC)technique.We employ boundary variation diminishing(BVD)reconstruction to enhance the scheme’s effectiveness in handling shocks.First,we prove that wavelet collocation upwind schemes based on interpolating wavelets can be reformulated into a conservative form within the framework of wavelet theory,forming the foundation of the proposed scheme.The new fourthorder accurate scheme possesses significantly better spectral resolution than the fifth-and even seventh-order WENO-Z(weighted essentially non-oscillatory)schemes over the entire wave-number range.Moreover,the inherent low-pass filtering property of the wavelet bases allows them to filter high-frequency numerical oscillations,endowing the wavelet upwind scheme with robustness and accuracy in solving problems under extreme conditions.Notably,due to the wavelet multiresolution approximation,the proposed scheme possesses a distinctive shape-preserving property absent in the WENO-Z schemes and the fifth-order schemes with BVD reconstruction based on polynomials.Furthermore,compared to the fifth-order scheme with BVD reconstruction based on polynomials—which is significantly superior to the WENO schemes—the proposed scheme further enhances the ability to capture discontinuities.展开更多
The recovery and utilization of ubiquitous low-grade heat are crucial for mitigating the fossil energy crisis.However,uncontrolled spontaneous heat dissipation limits its practical application.Inspired by skeletal mus...The recovery and utilization of ubiquitous low-grade heat are crucial for mitigating the fossil energy crisis.However,uncontrolled spontaneous heat dissipation limits its practical application.Inspired by skeletal muscle thermogenesis,we develop a compressible wood phase change gel with mechano-controlled heat release by infiltrating xylitol gel into wood aerogel.The xylitol gel can store recovered low-grade heat for at least 1 month by leveraging its inherent energy barrier.The hierarchically aligned lamellar structure of wood aerogel facilitates mechanical adaptation,hydrogen bond formation,and energy dissipation between the wood aerogel and the xylitol gel,increasing the compressive strength and toughness of wood phase change gel fivefold compared to xylitol gel.This enhancement effect enables repetitive contact-separation motions between the wood phase change gel and the substrate during radial compression,overcoming the energy barrier and releasing approximately 178.6 J g−1 of heat.As a proof-of-concept,the wood phase change gel serves as the hot side in a thermoelectric generator,providing about 2.13 W m^(−2) of clean electricity by the controlled utilization of recovered solar heat.This study presents a sustainable method to achieve off-grid electricity generation through the controlled utilization of recovered low-grade heat.展开更多
Multiphase flows widely exist in various scientific and engineering fields,and strongly compressible multiphase flows commonly occur in practical applications,which makes them an important part of computational fluid ...Multiphase flows widely exist in various scientific and engineering fields,and strongly compressible multiphase flows commonly occur in practical applications,which makes them an important part of computational fluid dynamics.In this study,an axisymmetric adaptive multiresolution smooth particle hydrodynamics(SPH)model is proposed to solve various strongly compressible multiphase flow problems.In the present model,the governing equations are discretized in cylindrical polar coordinates,and an improved volume adaptive scheme is developed to better solve the problem of excessive volume change in strongly compressible multiphase flows.On this basis,combined with the adaptive particle refinement technique,an adaptive multiresolution scheme is proposed in this study.In addition,the high-order differential operator and diffusion correction term are utilized to improve the accuracy and stability.The effectiveness of the model is verified by testing four typical strongly compressible multiphase flow problems.By comparing the results of adaptive multiresolution SPH with other numerical results or experimental data,we can conclude that the present SPH method effectively models strongly compressible multiphase flows.展开更多
We investigate the long time existence of strong solutions to the initial value problem for the three-dimensional non-isentropic compressible Navier-Stokes-Korteweg system.Under the conditions of slight density and te...We investigate the long time existence of strong solutions to the initial value problem for the three-dimensional non-isentropic compressible Navier-Stokes-Korteweg system.Under the conditions of slight density and temperature variations,we verify that the full compressible Navier-Stokes-Korteweg equations admit a unique strong solution as long as the solution of the limiting system exists,when the Mach number is sufficiently small.Furthermore,we deduce the uniform convergence of strong solutions for the compressible system toward those for the corresponding incompressible system on the time interval in which the solution exists.展开更多
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.展开更多
This article is devoted to the study of the so-called incompressible limit for solutions of the compressible liquid crystals system. We consider the problem in the whole space R^N and a bounded domain of R^N with Diri...This article is devoted to the study of the so-called incompressible limit for solutions of the compressible liquid crystals system. We consider the problem in the whole space R^N and a bounded domain of R^N with Dirichlet boundary conditions. Here we set the number of dimension N = 2 or 3.展开更多
A hybrid central-upwind scheme is proposed. Two sub-schemes, the central difference scheme and the Roets flux difference splitting scheme, are hybridized by means of a binary sensor function. In order to examine the c...A hybrid central-upwind scheme is proposed. Two sub-schemes, the central difference scheme and the Roets flux difference splitting scheme, are hybridized by means of a binary sensor function. In order to examine the capability of the proposed hybrid scheme in computing compressible turbulent flow around a curved surface body, especially the flow involving shock wave, three typical eases are investigated by using detached-eddy simulation technique. Numerical results show good agreements with the experimental measurements. The present hybrid scheme can be applied to simulating the compressible flow around a curved surface body involving shock wave and turbulence.展开更多
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.展开更多
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.展开更多
This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-laye...This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.展开更多
The cavitated bifurcation problem in a solid sphere composed oftwo compressible hyper-elas- tic materials is examined. Thebifurcation solution for the composed sphere under a uniform radialtensile boundary dead-load i...The cavitated bifurcation problem in a solid sphere composed oftwo compressible hyper-elas- tic materials is examined. Thebifurcation solution for the composed sphere under a uniform radialtensile boundary dead-load is obtained. The bifurcation curves andthe stress contributions subsequent to the cavita- tion are given.The right and left bifurcation as well as the catastrophe andconcentration of stresses are ana- lyzed. The stability of solutionsis discussed through an energy comparison.展开更多
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.展开更多
This is a continuation of the article（Comm.Partial Differential Equations 26（2001）965）.In this article,the authors consider the one-dimensional compressible isentropic Navier-Stokes equations with gravitational fo...This is a continuation of the article（Comm.Partial Differential Equations 26（2001）965）.In this article,the authors consider the one-dimensional compressible isentropic Navier-Stokes equations with gravitational force,fixed boundary condition,a general pressure and the density-dependent viscosity coefficient when the viscous gas connects to vacuum state with a jump in density.Precisely,the viscosity coefficient μ is proportional to ρ^θ and 0〈θ〈1/2,where ρ is the density,and the pressure P=P（ρ）is a general pressure.The global existence and the uniqueness of weak solution are proved.展开更多
By analyzing the characteristics of low Mach number perfect gas flows, a novel Slightly Compressible Model (SCM) for low Mach number perect gas flows is derived. In view of numerical calculations, this model is proved...By analyzing the characteristics of low Mach number perfect gas flows, a novel Slightly Compressible Model (SCM) for low Mach number perect gas flows is derived. In view of numerical calculations, this model is proved very efficient, for it is kept within thep-v frame but does not have to satisfy the time consuming divergence-free condition in order to get the incompressible Navier-Stokes equation solution. Writing the equations in the form of conservation laws, we have derived the characteristic systems which are necessary for numerical calculations. A cell-centered finite-volume method with flux difference upwind-biased schemes is used for the equation solutions and a new Exact Newton Relaxation (ENR) implicit method is developed. Various computed results are presented to validate the present model. Laminar flow solutions over a circular cylinder with wake developing and vortex shedding are presented. Results for inviscid flow over a sphere are compared in excellent agreement with the exact analytic incompressible solution. Three-dimensional viscous flow solutions over sphere and prolate spheroid are also calculated and compared well with experiments and other incompressible solutions. Finally, good convergent performances are shown for sphere viscous flows.展开更多
This article is concerned with the time periodic solution to the isentropic compressible Navier-Stokes equations in a periodic domain. Using an approach of parabolic regularization, we first obtain the existence of th...This article is concerned with the time periodic solution to the isentropic compressible Navier-Stokes equations in a periodic domain. Using an approach of parabolic regularization, we first obtain the existence of the time periodic solution to a regularized problem under some smallness and symmetry assumptions on the external force. The result for the original compressible Navier-Stokes equations is then obtained by a limiting process. The uniqueness of the periodic solution is also given.展开更多
In this paper, we consider the global existence of classical solution to the 3-D compressible Navier-Stokes equations with a density-dependent viscosity coefficient λ(ρ)provided that the initial energy is small in s...In this paper, we consider the global existence of classical solution to the 3-D compressible Navier-Stokes equations with a density-dependent viscosity coefficient λ(ρ)provided that the initial energy is small in some sense. In our result, we give a relation between the initial energy and the viscosity coefficient μ, and it shows that the initial energy can be large if the coefficient of the viscosity μ is taken to be large, which implies that large viscosity μ means large solution.展开更多
基金financially supported by the National Key Research and Development Program of China(Grant Nos.2023YFB2604005)the National Key Research and Development 451 Program of China(Grant No.2021YFC3100803)the Yangtze River Water Science Research Joint Fund Key Project of National Natural Science Foundation of China(Grant No.U2340231).
文摘Polyurethane foam,when used as a compressible layer in deep soft rock tunnels,offers a feasible solution to reduce the support pressure on the secondary lining.The foam spraying method using sprayed polyurethane material is convenient for engineering applications;however,the compressive behaviour and feasibility of sprayed polyurethane material as a compressible layer remain unclear.To address this gap,this study conducts uniaxial compression tests and scanning electron microscope(SEM)tests to investigate the compressive behaviour of the rigid foams fabricated from a self-developed polyurethane spray material.A peridynamics model for the composite lining with a polyurethane compressible layer is then established.After validating the proposed method by comparison with two tests,a parametric study is carried out to investigate the damage evolution of the composite lining with a polyurethane compressible layer under various combinations of large deformations and compressible layer parameters.The results indicate that the polyurethane compressible layer effectively reduces the radial deformation and damage index of the secondary lining while increasing the damage susceptibility of the primary lining.The thickness of the polyurethane compressible layer significantly influences the prevention effect of large deformation-induced damage to the secondary lining within the density range of 50e100 kg/m^(3).In accordance with the experimental and simulation results,a simple,yet reasonable and convenient approach for determining the key parameters of the polyurethane compressible layer is proposed,along with a classification scheme for the parameters of the polyurethane compressible layer.
基金Supported by the National Key R&D Program of China(2022YFB3303500).
文摘Numerical simulations were conducted on a 10-blade Sevik rotor ingesting wake downstream of two turbulence-generating grids.These simulations were based on implicit large-eddy simulation(ILES)and the boundary data immersion method(BDIM)for compressible flows,which were solved using a fully self-programmed Fortran code.Results show that the predicted thrust spectrum aligns closely with the experimental measurements.In addition,it captures the thrust dipole directivity of the noise around the rotating propeller due to random pressure pulsations on the blades,as well as the flow structures simultaneously.Furthermore,the differences in the statistical characteristics,flow structures,and low-frequency broadband thrust spectra due to different turbulence levels were investigated.This analysis indicates that the interaction between the upstream,which is characterized by a lower turbulence level and a higher turbulent length of scale,and the rotating propeller results in a lower amplitude in force spectra and a slight increase in the scale of tip vortices.
基金supported by the Natural Science Foundation of Shandong Province of China(ZR2024MA033ZR2021QA049).
文摘This paper concerns the Cauchy problem of 3D compressible micropolar fluids in the whole space R^(3). For regular initial data with m0E0 is suitable small, where m0 and E0 represent the upper bound of initial density and initial energy, we prove that if ρ0 ∈ Lγ ∩ H3 with γ ∈ (1, 6), then the problem possesses a unique global classical solution on R^(3) × [0, T] with any T ∈ (0, ∞). It’s worth noting that both the vacuum states and possible random largeness of initial energy are allowed.
基金supported by the National Natural Science Foundation of China(Grant No.11925204).
文摘In this paper,we develop a fourth-order conservative wavelet-based shock-capturing scheme.The scheme is constructed by combining a wavelet collocation upwind method with the monotonic tangent of hyperbola for interface capturing(THINC)technique.We employ boundary variation diminishing(BVD)reconstruction to enhance the scheme’s effectiveness in handling shocks.First,we prove that wavelet collocation upwind schemes based on interpolating wavelets can be reformulated into a conservative form within the framework of wavelet theory,forming the foundation of the proposed scheme.The new fourthorder accurate scheme possesses significantly better spectral resolution than the fifth-and even seventh-order WENO-Z(weighted essentially non-oscillatory)schemes over the entire wave-number range.Moreover,the inherent low-pass filtering property of the wavelet bases allows them to filter high-frequency numerical oscillations,endowing the wavelet upwind scheme with robustness and accuracy in solving problems under extreme conditions.Notably,due to the wavelet multiresolution approximation,the proposed scheme possesses a distinctive shape-preserving property absent in the WENO-Z schemes and the fifth-order schemes with BVD reconstruction based on polynomials.Furthermore,compared to the fifth-order scheme with BVD reconstruction based on polynomials—which is significantly superior to the WENO schemes—the proposed scheme further enhances the ability to capture discontinuities.
基金supported by the National Key R&D Program of China (2023YFD2201403)the National Natural Science Foundation of China (Grant Nos. 32171693, 32201482)+1 种基金the Heilongjiang Natural Science Foundation Outstanding Youth project (Grant No. YQ2022C002)College Students'Innovative Entrepreneurial Training Plan Program (202410225338)
文摘The recovery and utilization of ubiquitous low-grade heat are crucial for mitigating the fossil energy crisis.However,uncontrolled spontaneous heat dissipation limits its practical application.Inspired by skeletal muscle thermogenesis,we develop a compressible wood phase change gel with mechano-controlled heat release by infiltrating xylitol gel into wood aerogel.The xylitol gel can store recovered low-grade heat for at least 1 month by leveraging its inherent energy barrier.The hierarchically aligned lamellar structure of wood aerogel facilitates mechanical adaptation,hydrogen bond formation,and energy dissipation between the wood aerogel and the xylitol gel,increasing the compressive strength and toughness of wood phase change gel fivefold compared to xylitol gel.This enhancement effect enables repetitive contact-separation motions between the wood phase change gel and the substrate during radial compression,overcoming the energy barrier and releasing approximately 178.6 J g−1 of heat.As a proof-of-concept,the wood phase change gel serves as the hot side in a thermoelectric generator,providing about 2.13 W m^(−2) of clean electricity by the controlled utilization of recovered solar heat.This study presents a sustainable method to achieve off-grid electricity generation through the controlled utilization of recovered low-grade heat.
基金Supported by the Guangxi Natural Science Foundation(Grant No.2021GXNSFBA196008)the Guangxi Science and Technology Development Program(Grant No.GuikeAD22035189).
文摘Multiphase flows widely exist in various scientific and engineering fields,and strongly compressible multiphase flows commonly occur in practical applications,which makes them an important part of computational fluid dynamics.In this study,an axisymmetric adaptive multiresolution smooth particle hydrodynamics(SPH)model is proposed to solve various strongly compressible multiphase flow problems.In the present model,the governing equations are discretized in cylindrical polar coordinates,and an improved volume adaptive scheme is developed to better solve the problem of excessive volume change in strongly compressible multiphase flows.On this basis,combined with the adaptive particle refinement technique,an adaptive multiresolution scheme is proposed in this study.In addition,the high-order differential operator and diffusion correction term are utilized to improve the accuracy and stability.The effectiveness of the model is verified by testing four typical strongly compressible multiphase flow problems.By comparing the results of adaptive multiresolution SPH with other numerical results or experimental data,we can conclude that the present SPH method effectively models strongly compressible multiphase flows.
文摘We investigate the long time existence of strong solutions to the initial value problem for the three-dimensional non-isentropic compressible Navier-Stokes-Korteweg system.Under the conditions of slight density and temperature variations,we verify that the full compressible Navier-Stokes-Korteweg equations admit a unique strong solution as long as the solution of the limiting system exists,when the Mach number is sufficiently small.Furthermore,we deduce the uniform convergence of strong solutions for the compressible system toward those for the corresponding incompressible system on the time interval in which the solution exists.
基金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.
基金supported by National Natural Science Foundation of China-NSAF: 11071043, 11131005
文摘This article is devoted to the study of the so-called incompressible limit for solutions of the compressible liquid crystals system. We consider the problem in the whole space R^N and a bounded domain of R^N with Dirichlet boundary conditions. Here we set the number of dimension N = 2 or 3.
基金Supported by the National Science Foundation for Post-doctoral Scientists of China(20100481141,201104567)the Natural Science Foundation of Jiangsu Province(BK2011723)the Planned Projects for Postdoctoral Research Foundation of Jiangsu Province(0902001C)~~
文摘A hybrid central-upwind scheme is proposed. Two sub-schemes, the central difference scheme and the Roets flux difference splitting scheme, are hybridized by means of a binary sensor function. In order to examine the capability of the proposed hybrid scheme in computing compressible turbulent flow around a curved surface body, especially the flow involving shock wave, three typical eases are investigated by using detached-eddy simulation technique. Numerical results show good agreements with the experimental measurements. The present hybrid scheme can be applied to simulating the compressible flow around a curved surface body involving shock wave and turbulence.
基金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.
基金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.
基金supported by the National Basic Research Program of China(2009CB724100)the National Natural Science Foundation of China(10632050,10872205,11072248).
文摘This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.
基金the National Natttral Science Foundation of China(No.19802012)
文摘The cavitated bifurcation problem in a solid sphere composed oftwo compressible hyper-elas- tic materials is examined. Thebifurcation solution for the composed sphere under a uniform radialtensile boundary dead-load is obtained. The bifurcation curves andthe stress contributions subsequent to the cavita- tion are given.The right and left bifurcation as well as the catastrophe andconcentration of stresses are ana- lyzed. The stability of solutionsis discussed through an energy comparison.
基金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.
基金Program for New Century ExcellentTalents in University(NCET-04-0745)the Key Project of the National Natural Science Foundation of China(10431060)
文摘This is a continuation of the article（Comm.Partial Differential Equations 26（2001）965）.In this article,the authors consider the one-dimensional compressible isentropic Navier-Stokes equations with gravitational force,fixed boundary condition,a general pressure and the density-dependent viscosity coefficient when the viscous gas connects to vacuum state with a jump in density.Precisely,the viscosity coefficient μ is proportional to ρ^θ and 0〈θ〈1/2,where ρ is the density,and the pressure P=P（ρ）is a general pressure.The global existence and the uniqueness of weak solution are proved.
基金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 gas flows, a novel Slightly Compressible Model (SCM) for low Mach number perect gas flows is derived. In view of numerical calculations, this model is proved very efficient, for it is kept within thep-v frame but does not have to satisfy the time consuming divergence-free condition in order to get the incompressible Navier-Stokes equation solution. Writing the equations in the form of conservation laws, we have derived the characteristic systems which are necessary for numerical calculations. A cell-centered finite-volume method with flux difference upwind-biased schemes is used for the equation solutions and a new Exact Newton Relaxation (ENR) implicit method is developed. Various computed results are presented to validate the present model. Laminar flow solutions over a circular cylinder with wake developing and vortex shedding are presented. Results for inviscid flow over a sphere are compared in excellent agreement with the exact analytic incompressible solution. Three-dimensional viscous flow solutions over sphere and prolate spheroid are also calculated and compared well with experiments and other incompressible solutions. Finally, good convergent performances are shown for sphere viscous flows.
基金supported by the Program for New Century Excellent Talents in University of the Ministry of Education(NCET-13-0804)NSFC(11471127)+3 种基金Guangdong Natural Science Funds for Distinguished Young Scholar(2015A030306029)The Excellent Young Teachers Program of Guangdong Province(HS2015007)Pearl River S&T Nova Program of Guangzhou(2013J2200064)supported by the General Research Fund of Hong Kong,City U 104511
文摘This article is concerned with the time periodic solution to the isentropic compressible Navier-Stokes equations in a periodic domain. Using an approach of parabolic regularization, we first obtain the existence of the time periodic solution to a regularized problem under some smallness and symmetry assumptions on the external force. The result for the original compressible Navier-Stokes equations is then obtained by a limiting process. The uniqueness of the periodic solution is also given.
文摘In this paper, we consider the global existence of classical solution to the 3-D compressible Navier-Stokes equations with a density-dependent viscosity coefficient λ(ρ)provided that the initial energy is small in some sense. In our result, we give a relation between the initial energy and the viscosity coefficient μ, and it shows that the initial energy can be large if the coefficient of the viscosity μ is taken to be large, which implies that large viscosity μ means large solution.
