PERTURBATIONAL FINITE VOLUME METHOD FOR THE SOLUTION OF 2-D NAVIER-STOKES EQUATIONS ON UNSTRUCTURED AND STRUCTURED COLOCATED MESHES

Based on the first-order upwind and second-order central type of finite volume (UFV and CFV) scheme, upwind and central type of perturbation finite volume (UPFV and CPFV) schemes of the Navier-Stokes equations were developed. In PFV method, the mass fluxes of across the cell faces of the control volume (CV) were expanded into power series of the grid spacing and the coefficients of the power series were determined by means of the conservation equation itself. The UPFV and CPFV scheme respectively uses the same nodes and expressions as those of the normal first-order upwind and second-order central scheme, which is apt to programming. The results of numerical experiments about the flow in a lid-driven cavity and the problem of transport of a scalar quantity in a known velocity field show that compared to the first-order UFV and second-order CFV schemes, upwind PFV scheme is higher accuracy and resolution, especially better robustness. The numerical computation to flow in a lid-driven cavity shows that the under-relaxation factor can be arbitrarily selected ranging from (0.3) to (0.8) and convergence perform excellent with Reynolds number variation from 10~2 to 10~4.

Three-dimensional double-rough-walled modeling of fluid flow through self-affine shear fractures

This study proposes a double-rough-walled fracture model to represent the natural geometries of rough fractures.The rough surface is generated using a modified successive random additions(SRA)algorithm and the aperture distribution during shearing is calculated using a mechanistic model.The shear-flow simulations are performed by directly solving the Navier-Stokes(NS)equations.The results show that the double-rough-walled fracture model can improve the accuracy of fluid flow simulations by approximately 14.99%-19.77%,compared with the commonly used single-rough-walled fracture model.The ratio of flow rate to hydraulic gradient increases by one order of magnitude for fluids in a linear flow regime with increment of shear displacement from 2.2 mm to 2.6 mm.By solving the NS equations,the inertial effect is taken into account and the significant eddies are simulated and numerically visualized,which are not easy to be captured in conventional experiments.The anisotropy of fluid flow in the linear regime during shearing is robustly enhanced as the shearing advances;however,it is either increased or decreased for fluids in the nonlinear flow regime,depending on the geometry of shear-induced void spaces between the two rough walls of the fracture.The present study provides a method to represent the real geometry of fractures during shearing and to simulate fluid flow by directly solving the NS equations,which can be potentially utilized in many applications such as heat and mass transfer,contaminant transport,and coupled hydro-thermo-mechanical processes within rock fractures/fracture networks.

Momentum and heat transfer of a special case of the unsteady stagnation-point flow

This paper investigates the unsteady stagnation-point flow and heat transfer over a moving plate with mass transfer,which is also an exact solution to the unsteady Navier-Stokes(NS)equations.The boundary layer energy equation is solved with the closed form solutions for prescribed wall temperature and prescribed wall heat flux conditions.The wall temperature and heat flux have power dependence on both time and spatial distance.The solution domain,the velocity distribution,the flow field,and the temperature distribution in the fluids are studied for different controlling parameters.These parameters include the Prandtl number,the mass transfer parameter at the wall,the wall moving parameter,the time power index,and the spatial power index.It is found that two solution branches exist for certain combinations of the controlling parameters for the flow and heat transfer problems.The heat transfer solutions are given by the confluent hypergeometric function of the first kind,which can be simplified into the incomplete gamma functions for special conditions.The wall heat flux and temperature profiles show very complicated variation behaviors.The wall heat flux can have multiple poles under certain given controlling parameters,and the temperature can have significant oscillations with overshoot and negative values in the boundary layers.The relationship between the number of poles in the wall heat flux and the number of zero-crossing points is identified.The difference in the results of the prescribed wall temperature case and the prescribed wall heat flux case is analyzed.Results given in this paper provide a rare closed form analytical solution to the entire unsteady NS equations,which can be used as a benchmark problem for numerical code validation.

Unsteady magnetohydrodynamic stagnation point flow — closed-form analytical solutions

This paper investigates the unsteady stagnation point flow and heat transfer of magnetohydrodynamic(MHD) fluids over a moving permeable flat surface. The unsteady Navier-Stokes(NS) equations are transformed into a similarity nonlinear ordinary differential equation, and a closed form solution is obtained for the unsteadiness parameter of 2. The boundary layer energy equation is transformed into a similarity equation,and is solved for a constant wall temperature and a time-dependent uniform wall heat flux case. The solution domain, velocity, and temperature profiles are calculated for different combinations of parameters including the Prandtl number, mass transfer parameter, wall moving parameter, and magnetic parameter. Two solution branches are obtained for certain combinations of the controlling parameters, and a stability analysis demonstrates that the lower solution branch is not stable. The present solutions provide an exact solution to the entire unsteady MHD NS equations, which can be used for validating the numerical code of computational fluid dynamics.

Gaussian process hydrodynamics

We present a Gaussian process(GP)approach,called Gaussian process hydrodynamics(GPH)for approximating the solution to the Euler and Navier-Stokes(NS)equations.Similar to smoothed particle hydrodynamics(SPH),GPH is a Lagrangian particle-based approach that involves the tracking of a finite number of particles transported by a flow.However,these particles do not represent mollified particles of matter but carry discrete/partial information about the continuous flow.Closure is achieved by placing a divergence-free GP priorξon the velocity field and conditioning it on the vorticity at the particle locations.Known physics(e.g.,the Richardson cascade and velocityincrement power laws)is incorporated into the GP prior by using physics-informed additive kernels.This is equivalent to expressingξas a sum of independent GPsξl,which we call modes,acting at different scales(each modeξlself-activates to represent the formation of eddies at the corresponding scales).This approach enables a quantitative analysis of the Richardson cascade through the analysis of the activation of these modes,and enables us to analyze coarse-grain turbulence statistically rather than deterministically.Because GPH is formulated by using the vorticity equations,it does not require solving a pressure equation.By enforcing incompressibility and fluid-structure boundary conditions through the selection of a kernel,GPH requires significantly fewer particles than SPH.Because GPH has a natural probabilistic interpretation,the numerical results come with uncertainty estimates,enabling their incorporation into an uncertainty quantification(UQ)pipeline and adding/removing particles(quanta of information)in an adapted manner.The proposed approach is suitable for analysis because it inherits the complexity of state-of-the-art solvers for dense kernel matrices and results in a natural definition of turbulence as information loss.Numerical experiments support the importance of selecting physics-informed kernels and illustrate the major impact of such kernels on the accuracy and stability.Because the proposed approach uses a Bayesian interpretation,it naturally enables data assimilation and predictions and estimations by mixing simulation data and experimental data.

Numerical Simulation and Experimental Investigation of 3 D Separated Flow Field around a Blunt Body

To study the effects of factors such as electromagnetic field and electron temperature on the ion extraction characteristics (extraction time, collision loss ratio) in atomic vapor laser isotope separation (AVLIS), 2 D electron equilibrium model was used to analyze and simulate the ion extraction process. The collisions between particles mainly considered charge exchanges between isotopic ions and atoms, which were treated by using cross section. The results show that the electric field and electron temperature have significant effects on the ion extraction characteristics: the stronger the electric field is, the higher the electron temperature is; the shorter the extraction time is, the less the collision loss ratio is, and moderate constant magnetic field has few effects on them. Key words atomic vapor laser isotope separation (AVLIS); ion extraction; charge exchange cross section; collision loss ratio

HYBRID CARTESIAN GRID/GRIDLESS METHOD FOR CALCULATING VISCOUS FLOWS OVER MULTI-ELEMENT AIRFOILS

A hybrid Cartesian grid/gridless method is developed for calculating viscous flows over multi-element airfoils.The method adopts an unstructured Cartesian grid to cover most areas of the computational domain and leaves only small region adjacent to the aerodynamic bodies to be filled with the cloud of points used in the gridless methods,which results in a better combination of the computational efficiency of the Cartesian grid and the flexibility of the gridless method in handling complex geometries.The clouds of points in the local gridless region are implemented in an anisotropic way according to the features of the thin boundary layer of the viscous flows over the airfoils,and the clouds of points at the vicinity of the interface between the grid and the gridless regions are also controlled by using an adaptive refinement technique during the generation of the unstructured Cartesian grid.An implementation of the resulting hybrid method is presented for solving two-dimensional compressible Navier-Stokes(NS)equations.The simulations of the viscous flows over a RAE2822airfoil or a two-element airfoil are successfully carried out,and the obtained results agree well with the available experimental data.