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.展开更多
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.展开更多
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.展开更多
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.展开更多
This paper is concerned with the global well-posedness of the solution to the compressible Navier-Stokes/Allen-Cahn system and its sharp interface limit in one-dimensional space.For the perturbations with small energy...This paper is concerned with the global well-posedness of the solution to the compressible Navier-Stokes/Allen-Cahn system and its sharp interface limit in one-dimensional space.For the perturbations with small energy but possibly large oscillations of rarefaction wave solutions near phase separation,and where the strength of the initial phase field could be arbitrarily large,we prove that the solution of the Cauchy problem exists for all time,and converges to the centered rarefaction wave solution of the corresponding standard two-phase Euler equation as the viscosity and the thickness of the interface tend to zero.The proof is mainly based on a scaling argument and a basic energy method.展开更多
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.展开更多
In this paper,we are concerned with a three-dimensional non-isothermal model for the compressible nematic liquid crystal flows in a periodic domain.Under some smallness and structural assumptions imposed on the time-p...In this paper,we are concerned with a three-dimensional non-isothermal model for the compressible nematic liquid crystal flows in a periodic domain.Under some smallness and structural assumptions imposed on the time-periodic force,we establish the existence of the time-periodic solutions to the system by using a regularized approximation scheme and the topological degree theory.We also prove a uniqueness result via energy estimates.展开更多
The polyurethane foam(PU)compressible layer is a viable solution to the problem of damage to the secondary lining in squeezing tunnels.Nevertheless,the mechanical behaviour of the multi-layer yielding supports has not...The polyurethane foam(PU)compressible layer is a viable solution to the problem of damage to the secondary lining in squeezing tunnels.Nevertheless,the mechanical behaviour of the multi-layer yielding supports has not been thoroughly investigated.To fill this gap,large-scale model tests were conducted in this study.The synergistic load-bearing mechanics were analyzed using the convergenceconfinement method.Two types of multi-layer yielding supports with different thicknesses(2.5 cm,3.75 cm and 5 cm)of PU compressible layers were investigated respectively.Digital image correlation(DIC)analysis and acoustic emission(AE)techniques were used for detecting the deformation fields and damage evolution of the multi-layer yielding supports in real-time.Results indicated that the loaddisplacement relationship of the multi-layer yielding supports could be divided into the crack initiation,crack propagation,strain-hardening,and failure stages.Compared with those of the stiff support,the toughness,deformability and ultimate load of the yielding supports were increased by an average of 225%,61%and 32%,respectively.Additionally,the PU compressible layer is positioned between two primary linings to allow the yielding support to have greater mechanical properties.The analysis of the synergistic bearing effect suggested that the thickness of PU compressible layer and its location significantly affect the mechanical properties of the yielding supports.The use of yielding supports with a compressible layer positioned between the primary and secondary linings is recommended to mitigate the effects of high geo-stress in squeezing tunnels.展开更多
We present a class of arbitrarily high order fully explicit kinetic numerical methods in compressible fluid dynamics,both in time and space,which include the relaxation schemes by Jin and Xin.These methods can use the...We present a class of arbitrarily high order fully explicit kinetic numerical methods in compressible fluid dynamics,both in time and space,which include the relaxation schemes by Jin and Xin.These methods can use the CFL number larger or equal to unity on regular Cartesian meshes for the multi-dimensional case.These kinetic models depend on a small parameter that can be seen as a"Knudsen"number.The method is asymptotic preserving in this Knudsen number.Also,the computational costs of the method are of the same order of a fully explicit scheme.This work is the extension of Abgrall et al.(2022)[3]to multidimensional systems.We have assessed our method on several problems for two-dimensional scalar problems and Euler equations and the scheme has proven to be robust and to achieve the theoretically predicted high order of accuracy on smooth solutions.展开更多
In this paper,we establish some regularity conditions on the density and velocity fields to guarantee the energy conservation of the weak solutions for the three-dimensional compressible nematic liquid crystal flow in...In this paper,we establish some regularity conditions on the density and velocity fields to guarantee the energy conservation of the weak solutions for the three-dimensional compressible nematic liquid crystal flow in the periodic domain.展开更多
This paper presents a high-order discontinuous Galerkin(DG)finite-element method to solve the barotropic version of the conservative symmetric hyperbolic and thermodynamically compatible(SHTC)model of compressible two...This paper presents a high-order discontinuous Galerkin(DG)finite-element method to solve the barotropic version of the conservative symmetric hyperbolic and thermodynamically compatible(SHTC)model of compressible two-phase flow,introduced by Romenski et al.in[59,62],in multiple space dimensions.In the absence of algebraic source terms,the model is endowed with a curl constraint on the relative velocity field.In this paper,the hyperbolicity of the system is studied for the first time in the multidimensional case,showing that the original model is only weakly hyperbolic in multiple space dimensions.To restore the strong hyperbolicity,two different methodologies are used:(i)the explicit symmetrization of the system,which can be achieved by adding terms that contain linear combinations of the curl involution,similar to the Godunov-Powell terms in the MHD equations;(ii)the use of the hyperbolic generalized Lagrangian multiplier(GLM)curl-cleaning approach forwarded.The PDE system is solved using a high-order ADER-DG method with a posteriori subcell finite-volume limiter to deal with shock waves and the steep gradients in the volume fraction commonly appearing in the solutions of this type of model.To illustrate the performance of the method,several different test cases and benchmark problems have been run,showing the high order of the scheme and the good agreement when compared to reference solutions computed with other well-known methods.展开更多
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.展开更多
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.展开更多
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 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.展开更多
基金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.
基金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 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 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 National Natural Science Foundation of China(12361044)supported by the National Natural Science Foundation of China(12171024,11971217,11971020)supported by the Academic and Technical Leaders Training Plan of Jiangxi Province(20212BCJ23027)。
文摘This paper is concerned with the global well-posedness of the solution to the compressible Navier-Stokes/Allen-Cahn system and its sharp interface limit in one-dimensional space.For the perturbations with small energy but possibly large oscillations of rarefaction wave solutions near phase separation,and where the strength of the initial phase field could be arbitrarily large,we prove that the solution of the Cauchy problem exists for all time,and converges to the centered rarefaction wave solution of the corresponding standard two-phase Euler equation as the viscosity and the thickness of the interface tend to zero.The proof is mainly based on a scaling argument and a basic energy method.
基金supported by the 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.
基金partially supported by the Science and Technology Research Program of Chongqing Municipal Education Commission(KJQN202100523,KJQN202000536)the National Natural Science Foundation of China(12001074)+3 种基金the Natural Science Foundation of Chongqing(cstc2020jcyj-msxmX0606)supported by the National Natural Science Foundation of Chongqing(CSTB2023NSCQ-MSX0278)the Science and Technology Research Program of Chongqing Municipal Education Commission(KJZD-K202100503)the Research Project of Chongqing Education Commission(CXQT21014)。
文摘In this paper,we are concerned with a three-dimensional non-isothermal model for the compressible nematic liquid crystal flows in a periodic domain.Under some smallness and structural assumptions imposed on the time-periodic force,we establish the existence of the time-periodic solutions to the system by using a regularized approximation scheme and the topological degree theory.We also prove a uniqueness result via energy estimates.
基金supported by the National Key Research and Development Program of China (Grant No.2021YFB2600800)the National Key Research and Development 451 Program of China (Grant No.2021YFC3100803)the Guangdong Innovative and Entrepreneurial Research Team Program (Grant No.2016ZT06N340).
文摘The polyurethane foam(PU)compressible layer is a viable solution to the problem of damage to the secondary lining in squeezing tunnels.Nevertheless,the mechanical behaviour of the multi-layer yielding supports has not been thoroughly investigated.To fill this gap,large-scale model tests were conducted in this study.The synergistic load-bearing mechanics were analyzed using the convergenceconfinement method.Two types of multi-layer yielding supports with different thicknesses(2.5 cm,3.75 cm and 5 cm)of PU compressible layers were investigated respectively.Digital image correlation(DIC)analysis and acoustic emission(AE)techniques were used for detecting the deformation fields and damage evolution of the multi-layer yielding supports in real-time.Results indicated that the loaddisplacement relationship of the multi-layer yielding supports could be divided into the crack initiation,crack propagation,strain-hardening,and failure stages.Compared with those of the stiff support,the toughness,deformability and ultimate load of the yielding supports were increased by an average of 225%,61%and 32%,respectively.Additionally,the PU compressible layer is positioned between two primary linings to allow the yielding support to have greater mechanical properties.The analysis of the synergistic bearing effect suggested that the thickness of PU compressible layer and its location significantly affect the mechanical properties of the yielding supports.The use of yielding supports with a compressible layer positioned between the primary and secondary linings is recommended to mitigate the effects of high geo-stress in squeezing tunnels.
基金funded by the SNF project 200020_204917 entitled"Structure preserving and fast methods for hyperbolic systems of conservation laws".
文摘We present a class of arbitrarily high order fully explicit kinetic numerical methods in compressible fluid dynamics,both in time and space,which include the relaxation schemes by Jin and Xin.These methods can use the CFL number larger or equal to unity on regular Cartesian meshes for the multi-dimensional case.These kinetic models depend on a small parameter that can be seen as a"Knudsen"number.The method is asymptotic preserving in this Knudsen number.Also,the computational costs of the method are of the same order of a fully explicit scheme.This work is the extension of Abgrall et al.(2022)[3]to multidimensional systems.We have assessed our method on several problems for two-dimensional scalar problems and Euler equations and the scheme has proven to be robust and to achieve the theoretically predicted high order of accuracy on smooth solutions.
基金support by the NSFC(12071391,12231016)the Guangdong Basic and Applied Basic Research Foundation(2022A1515010860)support by the China Postdoctoral Science Foundation(2023M742401)。
文摘In this paper,we establish some regularity conditions on the density and velocity fields to guarantee the energy conservation of the weak solutions for the three-dimensional compressible nematic liquid crystal flow in the periodic domain.
基金Initiative 2018–2027 attributed to DICAM of the University of Trento(grant L.232/2016)the PRIN 2022 project High-order structure-preserving semi-implicit schemes for hyperbolic equations and by the European Union-Next GenerationEU(PNRR,Spoke 7 CN HPC).
文摘This paper presents a high-order discontinuous Galerkin(DG)finite-element method to solve the barotropic version of the conservative symmetric hyperbolic and thermodynamically compatible(SHTC)model of compressible two-phase flow,introduced by Romenski et al.in[59,62],in multiple space dimensions.In the absence of algebraic source terms,the model is endowed with a curl constraint on the relative velocity field.In this paper,the hyperbolicity of the system is studied for the first time in the multidimensional case,showing that the original model is only weakly hyperbolic in multiple space dimensions.To restore the strong hyperbolicity,two different methodologies are used:(i)the explicit symmetrization of the system,which can be achieved by adding terms that contain linear combinations of the curl involution,similar to the Godunov-Powell terms in the MHD equations;(ii)the use of the hyperbolic generalized Lagrangian multiplier(GLM)curl-cleaning approach forwarded.The PDE system is solved using a high-order ADER-DG method with a posteriori subcell finite-volume limiter to deal with shock waves and the steep gradients in the volume fraction commonly appearing in the solutions of this type of model.To illustrate the performance of the method,several different test cases and benchmark problems have been run,showing the high order of the scheme and the good agreement when compared to reference solutions computed with other well-known methods.
基金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 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.
基金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.
基金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.
