THE RIEMANN PROBLEM FOR ISENTROPIC COMPRESSIBLE EULER EQUATIONS WITH DISCONTINUOUS FLUX

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.

STABILITY OF THE RAREFACTION WAVE IN THE SINGULAR LIMIT OF A SHARP INTERFACE PROBLEM FOR THE COMPRESSIBLE NAVIER-STOKES/ALLEN-CAHN SYSTEM

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.

THE EXISTENCE AND UNIQUENESS OF TIME-PERIODIC SOLUTIONS TO THE NON-ISOTHERMAL MODEL FOR COMPRESSIBLE NEMATIC LIQUID CRYSTALS IN A PERIODIC DOMAIN

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.

An Arbitrarily High Order and Asymptotic Preserving Kinetic Scheme in Compressible Fluid Dynamic

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.

ENERGY CONSERVATION FOR THE WEAK SOLUTIONS TO THE 3D COMPRESSIBLE NEMATIC LIQUID CRYSTAL FLOW

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.

Radon Measure Solutions to Riemann Problems for Isentropic Compressible Euler Equations of Polytropic Gases

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.

Artificial neural network-based subgrid-scale models for LES of compressible turbulent channel flow

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.

GLOBAL WEAK SOLUTIONS OF COMPRESSIBLE NAVIER-STOKES-LANDAU-LIFSHITZ-MAXWELL EQUATIONS FOR QUANTUM FLUIDS IN DIMENSION THREE

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 with symmetric velocity discretization for compressible flow

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.

Application of shifted lattice model to 3D compressible lattice Boltzmann method

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.

ON LOCAL CONTROLLABILITY FOR COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH DENSITY DEPENDENT VISCOSITIES

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.

GLOBAL SOLUTIONS TO THE 2D COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH SOME LARGE INITIAL DATA

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.

THE LOW MACH NUMBER LIMIT FOR ISENTROPIC COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH A REVISED MAXWELL'S LAW

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.

Gas kinetic flux solver based finite volume weighted essentially non-oscillatory scheme for inviscid compressible flows

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.

Impact of Artificial Compressibility on the Numerical Solution of Incompressible Nanofluid Flow

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.

Absorbable and unidirectionally compressible intestine-intestine stapler

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.

Blowup of Solutions to the Non-Isentropic Compressible Euler Equations with Time-Dependent Damping and Vacuum

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.

Numerical Study on Low-Reynolds Compressible Flows around Mars Helicopter Rotor Blade Airfoil

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 Cp 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.

COMPRESSIBLE FLOW SIMULATION AROUND AIRFOIL BASED ON LATTICE BOLTZMANN METHOD

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.

INCOMPRESSIBLE LIMIT OF A COMPRESSIBLE LIQUID CRYSTALS SYSTEM

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.