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.

STABILIZATION EFFECT OF FRICTIONS FOR TRANSONIC SHOCKS IN STEADY COMPRESSIBLE EULER FLOWS PASSING THREE-DIMENSIONAL DUCTS

Transonic shocks play a pivotal role in designation of supersonic inlets and ramjets.For the three-dimensional steady non-isentropic compressible Euler system with frictions,we constructe a family of transonic shock solutions in rectilinear ducts with square cross-sections.In this article,we are devoted to proving rigorously that a large class of these transonic shock solutions are stable,under multidimensional small perturbations of the upcoming supersonic flows and back pressures at the exits of ducts in suitable function spaces.This manifests that frictions have a stabilization effect on transonic shocks in ducts,in consideration of previous works which shown that transonic shocks in purely steady Euler flows are not stable in such ducts.Except its implications to applications,because frictions lead to a stronger coupling between the elliptic and hyperbolic parts of the three-dimensional steady subsonic Euler system,we develop the framework established in previous works to study more complex and interesting Venttsel problems of nonlocal elliptic equations.

Three-Dimensional Lattice Boltzmann Model for High-Speed Compressible Flows

A highly efficient three-dimensional （31）） Lattice Boltzmann （LB） model for high-speed compressible flows is proposed. This model is developed from the original one by Kataoka and Tsutahara [Phys. Rev. E 69 （2004） 056702]. The convection term is discretized by the Non-oscillatory, containing No free parameters and Dissipative （NND） scheme, which effectively damps oscillations at discontinuities. To be more consistent with the kinetic theory of viscosity and to further improve the numerical stability, an additional dissipation term is introduced. Model parameters are chosen in such a way that the von Neumann stability criterion is satisfied. The new model is validated by well-known benchmarks, （i） Riemann problems, including the problem with Lax shock tube and a newly designed shock tube problem with high Mach number; （ii） reaction of shock wave on droplet or bubble. Good agreements are obtained between LB results and exact ones or previously reported solutions. The model is capable of simulating flows from subsonic to supersonic and capturing jumps resulted from shock waves.

NUMERICAL INVESTIGATION OF THREE-DIMENSIONAL VISCOUS INCOMPRESSIBLE FLOWS IN DIVERGENT CURVED CHANNELS AND TURBULENT MODEL STUDY

In order to make the numerical calculation of viscous flows more convenient for the flows in channel with complicated profile governing equations expressed in the arbitrary curvilinear coordinates were derived by means of Favre density-weighted averaged method, and a turbulent model with effect of curvature modification was also derived. The numerical calculation of laminar and turbulent flown in divergent curved channels was carried out by means of parabolizeil computation method. The calculating results were used to analyze and investigate the aerodynamic performance of talor cascades in compressors preliminarily.

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.

Effect of Weft Binding Structure on Compression Properties of Three-Dimensional Woven Spacer Fabrics and Composites

With the wide use of three-dimensional woven spacer composites(3DWSCs),the market expects greater mechanical properties from this material.By changing the weft fastening method of the traditional I-shape pile yarns,we designed three-dimensional woven spacer fabrics(3DWSFs)and 3DWSCs with the weft V-shape to improve the compression performance of traditional 3DWSFs.The effects of weft binding structures,V-pile densities,and V-shaped angle were investigated in this paper.It is found that the compression resistance of 3DWSFs with the weft V-shape is improved compared to that with the weft I-shape,the fabric height recovery rate is as high as 95.7%,and the average elastic recovery rate is 59.39%.When the interlayer pile yarn density is the same,the weft V-shaped and weft I-shaped 3DWSCs have similar flatwise pressure and edgewise pressure performance.The compression properties of the composite improve as the density of the V-pile yarns increases.The flatwise compression load decreases as the V-shaped angle decreases.When the V-shaped angle is 28°and 42°,the latitudinal V-shaped 3DWSCs perform exceptionally well in terms of anti-compression cushioning.The V-shaped weft binding method offers a novel approach to structural design of 3DWSCs.

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.

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.

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.

UNSTEADY/STEADY NUMERICAL SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON ARTIFICIAL COMPRESSIBILITY

A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.

Three-dimensional numerical simulation of the flow past a circular cylinder based on LES method

The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three dimensional N-S equations were obtained by the finite volume method. The focus of this numerical simulation was to research the characteristics of pressure distribution （drag and litt forces） and vortex tubes at high Reynolds numbers. The results of the calculations showed that the forces at every section in the spanwise direction of the cylinder were symmetrical about the middle section and smaller than the forces calculated in two dimensional cases. Moreover, the flow around the cylinder obviously presents three dimensional 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.

GLOBAL EXISTENCE AND BLOW-UP PHENOMENA OF CLASSICAL SOLUTIONS FOR THE SYSTEM OF COMPRESSIBLE ADIABATIC FLOW THROUGH POROUS MEDIA

By means of maximum principle for nonlinear hyperbolic systems, the results given by HSIAO Ling and D. Serre was improved for Cauchy problem of compressible adiabatic flow through porous media, and a complete result on the global existence and the blow-up phenomena of classical solutions of these systems. These results show that the dissipation is strong enough to preserve the smoothness of ‘small ’ solution.

Sports compression garments improve resting markers of venous return and muscle blood flow in male basketball players

Background:The benefits associated with sports compression garments are thought to be closely related to enhanced blood flow.However,findings are equivocal,possibly due to heterogeneity in the techniques used for measuring blood flow,the garment types used,and the pressures applied.This study combined Doppler ultrasound and near-infrared spectroscopy technologies to provide the first comprehensive assessment of the effects of 3 sports compression garment types on markers of venous return and muscle blood flow at rest.Methods:Resting lower-limb blood flow measures(markers of venous return,muscle blood flow,and muscle oxygenation)of 22 elite,junior,male basketball players(age=17.2±0.9 years,mean±SD)were assessed in 4 separate conditions:no compression(CON),compression tights(TIGHTS),compression shorts(SHORTS),and compression socks(SOCKS).Markers of venous return(cross-sectional area,time-averaged mean and peak blood flow velocity,and venous blood flow)were measured via Doppler ultrasound at the popliteal and common femoral veins.Muscle blood flow and muscle oxygenation were measured in the gastrocnemius medialis and vastus lateralis using near-infrared spectroscopy.Results:Popliteal markers of venous return were higher in TIGHTS compared to CON(p<0.01)and SHORTS(p<0.01),with SOCKS values higher compared with CON(p<0.05).Common femoral vein markers of venous return were higher for all conditions compared to CON(p<0.05),with TIGHTS values also higher compared to SOCKS(p<0.05).Gastrocnemius medialis blood flow was higher for TIGHTS compared to CON(p=0.000),SOCKS(p=0.012),and SHORTS(p=0.000),with SOCKS higher compared to SHORTS(p=0.046).Vastus lateralis blood flow was higher for TIGHTS compared to CON(p=0.028)and SOCKS(p=0.019),with SHORTS also higher compared to CON(p=0.012)and SOCKS(p=0.005).Gastrocnemius medialis oxygenation was higher for TIGHTS compared to CON(p=0.003),SOCKS(p=0.033),and SHORTS(p=0.003),with SOCKS higher compared to CON(p=0.044)and SHORTS(p=0.032).Vastus lateralis oxygenation was higher for TIGHTS compared to CON(p=0.020)and SOCKS(p=0.006).Conclusion:Markers of venous return,muscle blood flow,and muscle oxygenation are increased with sports compression garments.TIGHTS are most effective,potentially because of the larger body area compressed.

Direct numerical simulation of compressible turbulent flows

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.

A NOVEL SLIGHTLY COMPRESSIBLE MODEL FOR LOW MACH NUMBER PERFECT GAS FLOW CALCULATION

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.

A BLOW-UP CRITERION FOR COMPRESSIBLE VISCOUS HEAT-CONDUCTIVE FLOWS

We study an initial boundary value problem for the Navier-Stokes equations of compressible viscous heat-conductive fluids in a 2-D periodic domain or the unit square domain. We establish a blow-up criterion for the local strong solutions in terms of the gradient of the velocity only, which coincides with the famous Beale-Kato-Majda criterion for ideal incompressible flows.

A GHOST FLUID BASED FRONT TRACKING METHOD FOR MULTIMEDIUM COMPRESSIBLE FLOWS

Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.

Numerical simulation of the dimensional transformation of atomization in a supersonic aerodynamic atomization dust-removing nozzle based on transonic speed compressible flow

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.

Three-dimensional Computational Fluid Dynamics Modeling of Two-phase Flow in a Structured Packing Column

Characterizing the complex two-phase hydrodynamics in structured packed columns requires a power- ful modeling tool. The traditional two-dimensional model exhibits limitations when one attempts to model the de- tailed two-phase flow inside the columns. The present paper presents a three-dimensional computational fluid dy- namics （CFD） model to simulate the two-phase flow in a representative unit of the column. The unit consists of an CFD calculations on column packed with Flexipak 1Y were implemented within the volume of fluid （VOF） mathe- matical framework. The CFD model was validated by comparing the calculated thickness of liquid film with the available experimental data. Special attention was given to quantitative analysis of the effects of gravity on the hy- drodynamics. Fluctuations in the liquid mass flow rate and the calculated pressure drop loss were found to be quali- tatively in agreement with the experimental observations.