Bound-Preserving Discontinuous Galerkin Methods with Modified Patankar Time Integrations for Chemical Reacting Flows

In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.

An inverse analysis of fluid flow through granular media using differentiable lattice Boltzmann method

This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.

Simulation of gas–liquid two-phase flow in a flow-focusing microchannel with the lattice Boltzmann method

A lattice Boltzmann method for gas–liquid two-phase flow involving non-Newtonian fluids is developed. Bubble formation in a flow-focusing microchannel is simulated by the method. The influences of flow rate ratio, surface tension,wetting properties, and rheological characteristics of the fluid on the two-phase flow are analyzed. The results indicate that the flow pattern transfers from slug flow to dry-plug flow with a sufficiently small capillary number. Due to the presence of three-phase contact lines, the contact angle has a more significant effect on the dry-plug flow pattern than on the slug flow pattern. The deformation of the front and rear meniscus of a bubble in the shear-thinning fluid can be explained by the variation of the capillary number. The reduced viscosity and increased contact angle are beneficial for the drag reduction in a microchannel. It also demonstrates the effectiveness of the current method to simulate the gas–liquid two-phase flow in a microchannel.

Three-dimensional color particle image velocimetry based on a cross-correlation and optical flow method

Rainbow particle image velocimetry(PIV)can restore the three-dimensional velocity field of particles with a single camera;however,it requires a relatively long time to complete the reconstruction.This paper proposes a hybrid algorithm that combines the fast Fourier transform(FFT)based co-correlation algorithm and the Horn–Schunck(HS)optical flow pyramid iterative algorithm to increase the reconstruction speed.The Rankine vortex simulation experiment was performed,in which the particle velocity field was reconstructed using the proposed algorithm and the rainbow PIV method.The average endpoint error and average angular error of the proposed algorithm were roughly the same as those of the rainbow PIV algorithm;nevertheless,the reconstruction time was 20%shorter.Furthermore,the effect of velocity magnitude and particle density on the reconstruction results was analyzed.In the end,the performance of the proposed algorithm was verified using real experimental single-vortex and double-vortex datasets,from which a similar particle velocity field was obtained compared with the rainbow PIV algorithm.The results show that the reconstruction speed of the proposed hybrid algorithm is approximately 25%faster than that of the rainbow PIV algorithm.

Arbitrary High-Order Fully-Decoupled Numerical Schemes for Phase-Field Models of Two-Phase Incompressible Flows

Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.

THE NONLINEAR STABILITY OF PLANE PARALLEL SHEAR FLOWS WITH RESPECT TO TILTED

The nonlinear stability of plane parallel shear flows with respect to tilted perturbations is studied by energy methods.Tilted perturbation refers to the fact that perturbations form an angleθ∈(0,π/2)with the direction of the basic flows.By defining an energy functional,it is proven that plane parallel shear flows are unconditionally nonlinearly exponentially stable for tilted streamwise perturbation when the Reynolds number is below a certain critical value and the boundary conditions are either rigid or stress-free.In the case of stress-free boundaries,by taking advantage of the poloidal-toroidal decomposition of a solenoidal field to define energy functionals,it can be even shown that plane parallel shear flows are unconditionally nonlinearly exponentially stable for all Reynolds numbers,where the tilted perturbation can be either spanwise or streamwise.

Modelling of debris-flow susceptibility and propagation: a case study from Northwest Himalaya

The geological and geographical position of the Northwest Himalayas makes it a vulnerable area for mass movements particularly landslides and debris flows. Mass movements have had a substantial impact on the study area which is extending along Karakorum Highway(KKH) from Besham to Chilas. Intense seismicity, deep gorges, steep terrain and extreme climatic events trigger multiple mountain hazards along the KKH, among which debris flow is recognized as the most destructive geohazard. This study aims to prepare a field-based debris flow inventory map at a regional scale along a 200 km stretch from Besham to Chilas. A total of 117 debris flows were identified in the field, and subsequently, a point-based debris-flow inventory and catchment delineation were performed through Arc GIS analysis. Regional scale debris flow susceptibility and propagation maps were prepared using Weighted Overlay Method(WOM) and Flow-R technique sequentially. Predisposing factors include slope, slope aspect, elevation, Topographic Roughness Index(TRI), Topographic Wetness Index(TWI), stream buffer, distance to faults, lithology rainfall, curvature, and collapsed material layer. The dataset was randomly divided into training data(75%) and validation data(25%). Results were validated through the Receiver Operator Characteristics(ROC) curve. Results show that Area Under the Curve(AUC) using WOM model is 79.2%. Flow-R propagation of debris flow shows that the 13.15%, 22.94%, and 63.91% areas are very high, high, and low susceptible to debris flow respectively. The propagation predicated by Flow-R validates the naturally occurring debris flow propagation as observed in the field surveys. The output of this research will provide valuable input to the decision makers for the site selection, designing of the prevention system, and for the protection of current infrastructure.

Glacial debris flow susceptibility mapping based on combined models in the Parlung Tsangpo Basin,China

Machine learning(ML)-based prediction models for mapping hazard(e.g.,landslide and debris flow)susceptibility have been widely developed in recent research.However,in some specific areas,ML models have limited application because of the uncertainties in identifying negative samples.The Parlung Tsangpo Basin exemplifies a region prone to recurrent glacial debris flows(GDFs)and is characterized by a prominent landform featuring deep gullies.Considering the limitations of the ML model,we developed and compared two combined statistical models(FA-WE and FA-IC)based on factor analysis(FA),weight of evidence(WE),and the information content(IC)method.The final GDF susceptibility maps were generated by selecting 8 most important static factors and considering the influence of precipitation.The results show that the FA-IC model has the best performance.The areas with a very high susceptibility to GDFs are primarily located in the narrow valley section upstream,on both sides of the valley in the middle and downstream of the Parlung Tsangpo River,and in the narrow valley section of each tributary.These areas encompass 86 gullies and are characterized as"narrow and steep".

Volumetric lattice Boltzmann method for pore-scale mass diffusionadvection process in geopolymer porous structures

Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advection process within porous structures is essential for material design.In this study,we present advancements in the volumetric lattice Boltzmann method(VLBM)for modeling and simulating pore-scale diffusion-advection of radioactive isotopes within geopolymer porous structures.These structures are created using the phase field method(PFM)to precisely control pore architectures.In our VLBM approach,we introduce a concentration field of an isotope seamlessly coupled with the velocity field and solve it by the time evolution of its particle population function.To address the computational intensity inherent in the coupled lattice Boltzmann equations for velocity and concentration fields,we implement graphics processing unit(GPU)parallelization.Validation of the developed model involves examining the flow and diffusion fields in porous structures.Remarkably,good agreement is observed for both the velocity field from VLBM and multiphysics object-oriented simulation environment(MOOSE),and the concentration field from VLBM and the finite difference method(FDM).Furthermore,we investigate the effects of background flow,species diffusivity,and porosity on the diffusion-advection behavior by varying the background flow velocity,diffusion coefficient,and pore volume fraction,respectively.Notably,all three parameters exert an influence on the diffusion-advection process.Increased background flow and diffusivity markedly accelerate the process due to increased advection intensity and enhanced diffusion capability,respectively.Conversely,increasing the porosity has a less significant effect,causing a slight slowdown of the diffusion-advection process due to the expanded pore volume.This comprehensive parametric study provides valuable insights into the kinetics of isotope uptake in porous structures,facilitating the development of porous materials for nuclear waste treatment applications.

Fatigue Safety Assessment of Concrete Continuous Rigid Frame Bridge Based on Rain Flow Counting Method and Health Monitoring Data

The fatigue of concrete structures will gradually appear after being subjected to alternating loads for a long time,and the accidents caused by fatigue failure of bridge structures also appear from time to time.Aiming at the problem of degradation of long-span continuous rigid frame bridges due to fatigue and environmental effects,this paper suggests a method to analyze the fatigue degradation mechanism of this type of bridge,which combines long-term in-site monitoring data collected by the health monitoring system(HMS)and fatigue theory.In the paper,the authors mainly carry out the research work in the following aspects:First of all,a long-span continuous rigid frame bridge installed with HMS is used as an example,and a large amount of health monitoring data have been acquired,which can provide efficient information for fatigue in terms of equivalent stress range and cumulative number of stress cycles;next,for calculating the cumulative fatigue damage of the bridge structure,fatigue stress spectrum got by rain flow counting method,S-N curves and damage criteria are used for fatigue damage analysis.Moreover,it was considered a linear accumulation damage through the Palmgren-Miner rule for the counting of stress cycles.The health monitoring data are adopted to obtain fatigue stress data and the rain flow counting method is used to count the amplitude varying fatigue stress.The proposed fatigue reliability approach in the paper can estimate the fatigue damage degree and its evolution law of bridge structures well,and also can help bridge engineers do the assessment of future service duration.

Equivalent Method of Integrated Power Generation System of Wind, Photovoltaic and Energy Storage in Power Flow Calculation and Transient Simulation

针对工程实际开展风光储联合发电系统在潮流计算和机电暂态仿真中的等值方法研究，旨在为大容量风光储联合发电系统的并网仿真分析奠定基础。将潮流计算的等值分为单元机组和集电系统2部分来研究。单元机组等值采用根据不同控制模式选取不同节点类型的方法，针对集电系统等值提出基于损耗不变原则的方法。等值模型和详细模型的算例结果表明，潮流计算等值方法具有较好的精度。在机电暂态仿真动态等值中，基于实际工程计算的最严重工况分析原则，提出运行在满出力点的单机“倍乘”等值模型，为工程计算中的风光储联合发电系统动态等值提供了一种解决方案。

Calculating Method for Influence of Material Flow on Energy Consumption in Steel Manufacturing Process

From the viewpoint of systems energy conservation, the influences of material flow on its energy consumption in a steel manufacturing process is an important subject. The quantitative analysis of the relationship between material flow and the energy intensity is useful to save energy in steel industry. Based on the concept of standard material flow diagram, all possible situations of ferric material flow in steel manufacturing process are analyzed. The expressions of the influence of material flow deviated from standard material flow diagram on energy consumption are put forward.

Analytical investigation of Jeffery-Hamel flow with high magnetic field and nanoparticle by Adomian decomposition method

In this study, the effects of magnetic field and nanoparticle on the Jeffery- Hamel flow are studied using a powerful analytical method called the Adomian decomposition method （ADM）. The traditional Navier-Stokes equation of fluid mechanics and Maxwell＇s electromagnetism governing equations are reduced to nonlinear ordinary differential equations to model the problem. The obtained results are well agreed with that of the Runge-Kutta method. The present plots confirm that the method has high accuracy for different a, Ha, and Re numbers. The flow field inside the divergent channel is studied for various values of Hartmann ：number and angle of channel. The effect of nanoparticle volume fraction in the absence of magnetic field is investigated.

Simulation of bluff body stabilized flows with hybrid RANS and PDF method

The motivation of this study is to investigate the turbulence-chemistry interactions by using probability density function （PDF） method. A consistent hybrid Reynolds Averaged Navier-Stokes （RANS）/PDF method is used to simulate the turbulent non-reacting and reacting flows. The joint fluctuating velocity-frequency-composition PDF equation coupled with the Reynolds averaged density, momentum and energy equations are solved on unstructured meshes by the Lagrangian Monte Carlo （MC） method combined with the finite volume （FV） method. The simulation of the axisymmetric bluff body stabilized non-reacting flow fields is presented in this paper. The calculated length of the recirculation zone is in good agreement with the experimental data. Moreover, the significant change of the flow pattern with the increase of the jet-to-coflow momentum flux ratio is well predicted. In addition, comparisons are made between the joint PDF model and two different Reynolds stress models.

A control volume based finite element method for simulating incompressible two-phase flow in heterogeneous porous media and its application to reservoir engineering

Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.

The Construction Method for Solving Radial Flow Problem through the Homogeneous Reservoir

On the basis of similar structure of solutions of ordinary differential equation (ODE) boundary value problem, the similar construction method was put forward by solving problems of fluid flow in porous media through the homogeneous reservoir. It is indicate that the pressure distribution of dimensionless reservoir and bottom hole in Laplace space, which take on the radial flow, also shows similar structure, and the internal relationship between the above solutions were illustrated in detail.

RESEARCH ON METHOD TO CALCULATE VELOCITIES OF SOLID PHASE AND LIQUID PHASE IN DEBRIS FLOW

Velocities of solid phase and liquid phase in debris flow are one key problem to research on impact and abrasion mechanism of banks and control structures under action of debris flow. Debris flow was simplified as two-phase liquid composed of solid phase with the same diameter particles and liquid phase with the same mechanical features. Assume debris flow was one-dimension two-phase liquid moving to one direction, then general equations of velocities of solid phase and liquid phase were founded in two-phase theory. Methods to calculate average pressures, volume forces and surface forces of debris flow control volume were established. Specially, surface forces were ascertained using Bingham＇s rheology equation of liquid phase and Bagnold＇s testing results about interaction between particles of solid phase. Proportional coefficient of velocities between liquid phase and solid phase was put forward, meanwhile, divergent coefficient between theoretical velocity and real velocity of solid phase was provided too. To state succinctly before, method to calculate velocities of solid phase and liquid phase was obtained through solution to general equations. The method is suitable for both viscous debris flow and thin debris flow. Additionally, velocities every phase can be identified through analyzing deposits in-situ after occurring of debris flow. It is obvious from engineering case the result in the method is consistent to that in real-time field observation.

Simulating high Reynolds number flow in two-dimensional lid-driven cavity by multi-relaxation-time lattice Boltzmann method

By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number （Re） flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.

A simple method for solving unidirectional methane gas flow in coal seam based on similarity solution

The equation used to model the unidirectional flow of methane gas in coal seams is usually formulated as a nonlinear partial differential equation, which needs to be solved numerically with a computer program.Nevertheless, for people without access to the computer program, the conventional numerical method may be inconvenient. Thus, the objective here is to seek some method simpler than the conventional one for solving the flow problem. A commonly used model of the unidirectional methane gas flow is considered, where the methane adsorption is described by the Langmuir isotherm and the free gas is treated as real gas. By introducing the similarity solution, a simple method for solving the flow model is proposed, which can be done on a hand calculator. It is shown by two examples that the gas pressure profile obtained by the proposed method agrees well with the direct numerical solution of the flow model.

MHD flow of upper-convected Maxwell fluid over porous stretching sheet using successive Taylor series linearization method

This paper investigates the magnetohydrodynamic （MHD） boundary layer flow of an incompressible upper-convected Maxwell （UCM） fluid over a porous stretching surface. Similarity transformations are used to reduce the governing partial differential equations into a kind of nonlinear ordinary differential equations. The nonlinear prob- lem is solved by using the successive Taylor series linearization method （STSLM）. The computations for velocity components are carried out for the emerging parameters. The numerical values of the skin friction coefficient are presented and analyzed for various parameters of interest in the problem.