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.

A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

Highly Accurate Golden Section Search Algorithms and Fictitious Time Integration Method for Solving Nonlinear Eigenvalue Problems

This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.

Resilience-Oriented Load Restoration Method and Repair Strategies for Regional Integrated Electricity-Natural Gas System

The rising frequency of extreme disaster events seriously threatens the safe and secure operation of the regional integrated electricity-natural gas system(RIENGS).With the growing level of coupling between electric and natural gas systems,it is critical to enhance the load restoration capability of both systems.This paper proposes a coordinated optimization strategy for resilience-enhanced RIENGS load restoration and repair scheduling and transforms it into a mixed integer second-order cone programming(MISOCP)model.The proposed model considers the distribution network reconfiguration and the coordinated repair strategy between the two systems,minimizing the total system load loss cost and repair time.In addition,a bi-directional gas flow model is used to describe the natural gas system,which can provide the RIENGS with more flexibility for load restoration during natural gas system failure.Finally,the effectiveness of the proposed approach is verified by conducting case studies on the test systems RIENGS E13-G7 and RIENGS E123-G20.

U-Net Models for Representing Wind Stress Anomalies over the Tropical Pacific and Their Integrations with an Intermediate Coupled Model for ENSO Studies

El Niño-Southern Oscillation(ENSO)is the strongest interannual climate mode influencing the coupled ocean-atmosphere system in the tropical Pacific,and numerous dynamical and statistical models have been developed to simulate and predict it.In some simplified coupled ocean-atmosphere models,the relationship between sea surface temperature(SST)anomalies and wind stress(τ)anomalies can be constructed by statistical methods,such as singular value decomposition(SVD).In recent years,the applications of artificial intelligence(AI)to climate modeling have shown promising prospects,and the integrations of AI-based models with dynamical models are active areas of research.This study constructs U-Net models for representing the relationship between SSTAs andτanomalies in the tropical Pacific;the UNet-derivedτmodel,denoted asτUNet,is then used to replace the original SVD-basedτmodel of an intermediate coupled model(ICM),forming a newly AI-integrated ICM,referred to as ICM-UNet.The simulation results obtained from ICM-UNet demonstrate their ability to represent the spatiotemporal variability of oceanic and atmospheric anomaly fields in the equatorial Pacific.In the ocean-only case study,theτUNet-derived wind stress anomaly fields are used to force the ocean component of the ICM,the results of which also indicate reasonable simulations of typical ENSO events.These results demonstrate the feasibility of integrating an AI-derived model with a physics-based dynamical model for ENSO modeling studies.Furthermore,the successful integration of the dynamical ocean models with the AI-based atmospheric wind model provides a novel approach to ocean-atmosphere interaction modeling studies.

Research on the development methodology for clinical practice guidelines for organic integration of traditional Chinese medicine and Western medicine

Integrated traditional Chinese medicine(TCM)and Western medicine(WM)is a new medical science grounded in the knowledge bases of both TCM and WM,which then forms a unique modern medical system in China.Integrated TCM and WM has a long history in China,and has made important achievements in the process of clinical diagnosis and treatment.However,the methodological defects in currently published clinical practice guidelines(CPGs)limit its development.The organic integration of TCM and WM is a deeper integration of TCM and WM.To realize the progression of"integration"to"organic integration",a targeted and standardized guideline development methodology is needed.Therefore,the purpose of this study is to establish a standardized development procedure for clinical practice guidelines for the organic integration of TCM and WM to promote the systematic integration of TCM and WM research results into clinical practice guidelines in order to achieve optimal results as the whole is greater than the sum of the parts.

Application of Choquet Integral-Importance-Performance Analysis and TOPSIS Methods in Approaching the Preference Factors of Calligraphy and Seal Engraving Imagery

Classical Chinese characters,presented through calligraphy,seal engraving,or painting,can exhibit different aesthetics and essences of Chinese characters,making them the most important asset of the Chinese people.Calligraphy and seal engraving,as two closely related systems in traditional Chinese art,have developed through the ages.Due to changes in lifestyle and advancements in modern technology,their original functions of daily writing and verification have gradually diminished.Instead,they have increasingly played a significant role in commercial art.This study utilizes the Evaluation Grid Method(EGM)and the Analytic Hierarchy Process(AHP)to research the key preference factors in the application of calligraphy and seal engraving imagery.Different from the traditional 5-point equal interval semantic questionnaire,this study employs a non-equal interval semantic questionnaire with a golden ratio scale,distinguishing the importance ratio of adjacent semantic meanings and highlighting the weighted emphasis on visual aesthetics.Additionally,the study uses Importance-Performance Analysis(IPA)and Technique for Order Preference by Similarity to Ideal Solution(TOPSIS)to obtain the key preference sequence of calligraphy and seal engraving culture.Plus,the Choquet integral comprehensive evaluation is used as a reference for IPA comparison.It is hoped that this study can provide cultural imagery references and research methods,injecting further creativity into industrial design.

An efficient method for train-track-substructure dynamic interaction analysis by implicit-explicit integration and multi-time-step solution

In this work,a method is put forward to obtain the dynamic solution efficiently and accurately for a large-scale train-track-substructure(TTS)system.It is called implicit-explicit integration and multi-time-step solution method(abbreviated as mI-nE-MTS method).The TTS system is divided into train-track subsystem and substruc-ture subsystem.Considering that the root cause of low effi-ciency of obtaining TTS solution lies in solving the alge-braic equation of the substructures,the high-efficient Zhai method,an explicit integration scheme,can be introduced to avoid matrix inversion process.The train-track system is solved by implicitly Park method.Moreover,it is known that the requirement of time step size differs for different sub-systems,integration methods and structural frequency response characteristics.A multi-time-step solution is pro-posed,in which time step size for the train-track subsystem and the substructure subsystem can be arbitrarily chosen once satisfying stability and precision demand,namely the time spent for m implicit integral steps is equal to n explicit integral steps,i.e.,mI=nE as mentioned above.The numeri-cal examples show the accuracy,efficiency,and engineering practicality of the proposed method.

An Integrated Structure Analysis Method of Active Surface Antenna by Using the Simplified Actuator

The main surface of a large reflector antenna is composed of thousands of panels,which are inevitably deformed under natural load,leading to a great deterioration of electrical performance of the antenna.The active surface technique is an effective method to compensate antenna deformation error and has been widely used.The actuator is a complex component,it has not been established in the antenna structure analysis model,which limits the theoretical analysis ability of the active surface technology.To solve this problem,an integrated structure analysis method of active surface antenna by using the simplified actuator is proposed.First,according to the supporting characteristics and adjusting function of the actuator,the complex actuator is simplified a simple structure of support beams,support truss and adjustment beam.Second,the finite element model of the active surface antenna including the simplified actuator is established.Then,the relationship between the adjustment value(load)of adjustment beam and the deformation of the antenna structure is deduced,and the integrated analysis method for realizing the active adjustment of panels is established.Finally,the model and adjustment analysis method of the active surface antenna in this paper is applied to an 8 m antenna,and satisfactory structural analysis results are obtained,which shows the effectiveness and universality of the method,and provides a reference for the modeling and adjustment analysis of the active surface antenna.

Least square method based on Haar wavelet to solve multi-dimensional stochastic Ito-Volterra integral equations

This paper proposes a method combining blue the Haar wavelet and the least square to solve the multi-dimensional stochastic Ito-Volterra integral equation.This approach is to transform stochastic integral equations into a system of algebraic equations.Meanwhile,the error analysis is proven.Finally,the effectiveness of the approach is verified by two numerical examples.

Intelligent Small Sample Defect Detection of Concrete Surface Using Novel Deep Learning Integrating Improved YOLOv5

Dear Editor,This letter presents an intelligent small sample defect detection of concrete surface using novel deep learning integrating the improved YOLOv5 based on the Wasserstein GAN(WGAN)enhancement algorithm.The proposed method is capable of producing top-notch data sets to address the issues of insufficient samples and substandard quality.

Error Control Strategies for Numerical Integrations in Fast Collocation Methods

We propose two error control techniques for numerical integrations in fast multiscale collocation methods for solving Fredholm integral equations of the second kind with weakly singular kernels. Both techniques utilize quadratures for singular integrals using graded points. One has a polynomial order of accuracy if the integrand has a polynomial order of smoothness except at the singular point and the other has exponential order of accuracy if the integrand has an infinite order of smoothness except at the singular point. We estimate the order of convergence and computational complexity of the corresponding approximate solutions of the equation. We prove that the second technique preserves the order of convergence and computational complexity of the original collocation method. Numerical experiments are presented to illustrate the theoretical estimates.

A new method to evaluate integrated production system capacity in oil fields

To improve the design and management of an integrated production system(IPS),a set of mathematical models and workflows are developed for evaluating the capacity of an IPS at steady-state conditions.Combining the conservation laws with applicable multiphase fluid and choke models,these mathematical models are solved to characterize the hydraulics of an integrated system of reservoir,wells,chokes,flowlines,and separator at steady state.The controllable variables such as well count,choke size and separator pressure are adjusted to optimize the performance of the IPs at a specific time.It is found that increasing the well count can increase the bulk flow rate of the production network,but too many wells may increase the manifold pressure,leading to decline of single-well production.Increasing the choke size can improve the capacity of the IPs.The production of the IPs is negatively correlated with the separator pressure.With increasing separator pressure and decreasing choke size,the increment of total fluid production(the capacity of IPS)induced by increasing well count decreases.Validation tests with field examples show a maximum absolute deviation is 1.5%,demonstrating the robustness and validity of the proposed mathematical models and workflows.

Optimal synthesis of heat-integrated distillation configurations using the two-column superstructure

In the realm of the synthesis of heat-integrated distillation configurations,the conventional approach for exploring more heat integration possibilities typically entails the splitting of a single column into a twocolumn configuration.However,this approach frequently necessitates tedious enumeration procedures,resulting in a considerable computational burden.To surmount this formidable challenge,the present study introduces an innovative remedy:The proposition of a superstructure that encompasses both single-column and multiple two-column configurations.Additionally,a simultaneous optimization algorithm is applied to optimize both the process parameters and heat integration structures of the twocolumn configurations.The effectiveness of this approach is demonstrated through a case study focusing on industrial organosilicon separation.The results underscore that the superstructure methodology not only substantially mitigates computational time compared to exhaustive enumeration but also furnishes solutions that exhibit comparable performance.

A comparison study on structure-function relationship of polysaccharides obtained from sea buckthorn berries using different methods:antioxidant and bile acid-binding capacity

In this study,the structural characters,antioxidant activities and bile acid-binding ability of sea buckthorn polysaccharides(HRPs)obtained by the commonly used hot water(HRP-W),pressurized hot water(HRP-H),ultrasonic(HRP-U),acid(HRP-C)and alkali(HRP-A)assisted extraction methods were investigated.The results demonstrated that extraction methods had significant effects on extraction yield,monosaccharide composition,molecular weight,particle size,triple-helical structure,and surface morphology of HRPs except for the major linkage bands.Thermogravimetric analysis showed that HRP-U with filamentous reticular microstructure exhibited better thermal stability.The HRP-A with the lowest molecular weight and highest arabinose content possessed the best antioxidant activities.Moreover,the rheological analysis indicated that HRPs with higher galacturonic acid content and molecular weight showed higher viscosity and stronger crosslinking network(HRP-C,HRP-W and HRP-U),which exhibited stronger bile acid binding capacity.The present findings provide scientific evidence in the preparation technology of sea buckthorn polysaccharides with good antioxidant and bile acid binding capacity which are related to the structure affected by the extraction methods.

Drilling-based measuring method for the c-φ parameter of rock and its field application

The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(RCZ) is built. The formula for an ultimate cutting force is established based on the limit equilibrium principle. The relationship between digital drilling parameters(DDP) and the c-φ parameter(DDP-cφ formula, where c refers to the cohesion and φ refers to the internal friction angle) is derived, and the response of drilling parameters and cutting ratio to the strength parameters is analyzed. The drillingbased measuring method for the c-φ parameter of rock is constructed. The laboratory verification test is then completed, and the difference in results between the drilling test and the compression test is less than 6%. On this basis, in-situ rock drilling tests in a traffic tunnel and a coal mine roadway are carried out, and the strength parameters of the surrounding rock are effectively tested. The average difference ratio of the results is less than 11%, which verifies the effectiveness of the proposed method for obtaining the strength parameters based on digital drilling. This study provides methodological support for field testing of rock strength parameters.

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.

Are the altitudinal patterns of plant diversity derived from field surveys consistent with those from empirical integrated methods?

Field surveys and empirical integrated methods are commonly used in the ecological research to understand the altitudinal pattern of plant diversity of mountains.However,few studies have compared the differences between the two methods on the same scale.Here,we addressed and compared the altitudinal patterns of species richness(SR),phylogenetic diversity(PD),the standardized effect size of phylogenetic diversity(PDses)and mean phylogenetic distance(MPDses)of about 580 angiosperms growing on Mount Kenya from two independent datasets:one is based on our several times field surveys in this mountain and another one is based on empirical data integrated from literatures.We found that the altitudinal diversity patterns of field surveys dataset were consistent with the empirical integrated dataset.Both SR and PD showed hump-shaped patterns along the altitude,and both PDses and MPDses showed monotonically decreasing patterns along the altitude.The ratio of diversity between field surveys dataset and empirical integrated dataset were gradually increase along the altitude.Our research provides new insight for understanding the altitudinal diversity patterns of plants of a tropical mountain.

THE UNIFORM CONVERGENCE OF A DG METHOD FOR A SINGULARLY PERTURBED VOLTERRA INTEGRO-DIFFERENTIAL EQUATION

The purpose of this work is to implement a discontinuous Galerkin(DG)method with a one-sided flux for a singularly perturbed Volterra integro-differential equation(VIDE)with a smooth kernel.First,the regularity property and a decomposition of the exact solution of the singularly perturbed VIDE with the initial condition are provided.Then the existence and uniqueness of the DG solution are proven.Then some appropriate projection-type interpolation operators and their corresponding approximation properties are established.Based on the decomposition of the exact solution and the approximation properties of the projection type interpolants,the DG method achieves the uniform convergence in the L2 norm with respect to the singular perturbation parameter e when the space of polynomials with degree p is used.A numerical experiment validates the theoretical results.Furthermore,an ultra-convergence order 2p+1 at the nodes for the one-sided flux,uniform with respect to the singular perturbation parameter e,is observed numerically.

Material point method simulation of hydro-mechanical behaviour in twophase porous geomaterials: A state-of-the-art review

The material point method(MPM)has been gaining increasing popularity as an appropriate approach to the solution of coupled hydro-mechanical problems involving large deformation.In this paper,we survey the current state-of-the-art in the MPM simulation of hydro-mechanical behaviour in two-phase porous geomaterials.The review covers the recent advances and developments in the MPM and their extensions to capture the coupled hydro-mechanical problems involving large deformations.The focus of this review is aiming at providing a clear picture of what has or has not been developed or implemented for simulating two-phase coupled large deformation problems,which will provide some direct reference for both practitioners and researchers.