Implementation of a particle-in-cell method for the energy solver in 3D spherical geodynamic modeling

The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.

Modeling footing resting on anisotropic sand using material point method

Sand typically exhibits anisotropic internal structure which may significantly influence its mechanical behavior. The material point method (MPM) can eliminate mesh distortion and thus is suitable for investigating geotechnical problems with large deformation. In this study, an advanced anisotropic critical state theory (ACST)-based soil model is implemented in MPM to study the response of strip footing resting on anisotropic sand. The capability of the model is verified by simulating several element tests and strip footing tests with different soil densities and fabric bedding plane orientations. For the footing problem with a vertical load, as the fabric bedding plane orientation increases, the bearing capacity decreases and its corresponding settlement increases. The failure pattern becomes asymmetrical when the bedding plane orientation or the loading direction is inclined. A comparison between the simulation results predicted by the anisotropic and isotropic models is made, which demonstrates that neglecting the fabric anisotropy may lead to the overestimation of the bearing capacity.

Uncertainties of landslide susceptibility prediction: Influences of random errors in landslide conditioning factors and errors reduction by low pass filter method

In the existing landslide susceptibility prediction(LSP)models,the influences of random errors in landslide conditioning factors on LSP are not considered,instead the original conditioning factors are directly taken as the model inputs,which brings uncertainties to LSP results.This study aims to reveal the influence rules of the different proportional random errors in conditioning factors on the LSP un-certainties,and further explore a method which can effectively reduce the random errors in conditioning factors.The original conditioning factors are firstly used to construct original factors-based LSP models,and then different random errors of 5%,10%,15% and 20%are added to these original factors for con-structing relevant errors-based LSP models.Secondly,low-pass filter-based LSP models are constructed by eliminating the random errors using low-pass filter method.Thirdly,the Ruijin County of China with 370 landslides and 16 conditioning factors are used as study case.Three typical machine learning models,i.e.multilayer perceptron(MLP),support vector machine(SVM)and random forest(RF),are selected as LSP models.Finally,the LSP uncertainties are discussed and results show that:(1)The low-pass filter can effectively reduce the random errors in conditioning factors to decrease the LSP uncertainties.(2)With the proportions of random errors increasing from 5%to 20%,the LSP uncertainty increases continuously.(3)The original factors-based models are feasible for LSP in the absence of more accurate conditioning factors.(4)The influence degrees of two uncertainty issues,machine learning models and different proportions of random errors,on the LSP modeling are large and basically the same.(5)The Shapley values effectively explain the internal mechanism of machine learning model predicting landslide sus-ceptibility.In conclusion,greater proportion of random errors in conditioning factors results in higher LSP uncertainty,and low-pass filter can effectively reduce these random errors.

Modeling One Dimensional Two-Cell Model with Tumor Interaction Using Krylov Subspace Methods

A brain tumor occurs when abnormal cells grow, sometimes very rapidly, into an abnormal mass of tissue. The tumor can infect normal tissue, so there is an interaction between healthy and infected cell. The aim of this paper is to propose some efficient and accurate numerical methods for the computational solution of one-dimensional continuous basic models for the growth and control of brain tumors. After computing the analytical solution, we construct approximations of the solution to the problem using a standard second order finite difference method for space discretization and the Crank-Nicolson method for time discretization. Then, we investigate the convergence behavior of Conjugate gradient and generalized minimum residual as Krylov subspace methods to solve the tridiagonal toeplitz matrix system derived.

Multi-Objective Optimization of Fused Deposition Modeling for Mechanical Properties of Biopolymer Parts Using the Grey-Taguchi Method

The urgent need to develop customized functional products only possible by 3D printing had realized when faced with the unavailability of medical devices like surgical instruments during the coronavirus-19 disease and the ondemand necessity to perform surgery during space missions.Biopolymers have recently been the most appropriate option for fabricating surgical instruments via 3D printing in terms of cheaper and faster processing.Among all 3D printing techniques,fused deposition modelling(FDM)is a low-cost and more rapid printing technique.This article proposes the fabrication of surgical instruments,namely,forceps and hemostat using the fused deposition modeling(FDM)process.Excellent mechanical properties are the only indicator to judge the quality of the functional parts.The mechanical properties of FDM-processed parts depend on various process parameters.These parameters are layer height,infill pattern,top/bottom pattern,number of top/bottom layers,infill density,flow,number of shells,printing temperature,build plate temperature,printing speed,and fan speed.Tensile strength and modulus of elasticity are chosen as evaluation indexes to ascertain the mechanical properties of polylactic acid(PLA)parts printed by FDM.The experiments have performed through Taguchi’s L27orthogonal array(OA).Variance analysis(ANOVA)ascertains the significance of the process parameters and their percent contributions to the evaluation indexes.Finally,as a multiobjective optimization technique,grey relational analysis(GRA)obtains an optimal set of FDM process parameters to fabricate the best parts with comprehensive mechanical properties.Scanning electron microscopy(SEM)examines the types of defects and strong bonding between rasters.The proposed research ensures the successful fabrication of functional surgical tools with substantial ultimate tensile strength(42.6 MPa)and modulus of elasticity(3274 MPa).

3D forward modeling and response characteristics of low-resistivity overburden of the CFS-PML absorbing boundary for ground-well transient electromagnetic method

This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.

Displacement Field Variable Modeling Method for Heterogeneous Materials in Wind Power Blade Core Plates

In order to study the mechanical properties of the heterogeneous core plate of the wind turbine blade,a modeling method of the core plate based on displacement field variables is proposed.Firstly,the wind turbine blade core plate was modeled according to the theory of modeling heterogeneous material characteristics.Secondly,the three-point bending finite element model of the wind turbine blade core plate was solved by the display dynamic equation to obtain the deformation pattern and force-deformation relationship of the core plate.Finally,the three-point bending static test was conducted to compare with the finite element analysis.The test results show that:the damage form of the wind turbine blade core plate includes elasticity,yield,and failure stages.The main failure modes are plastic deformation,core material collapse,and panel-core delamination.The failure load measured by the test is 1.59 kN,which is basically consistent with the load-displacement result obtained by the simulation,with a difference of only 1.9%,which verifies the validity and reliability of the model.It provides data references for wind turbine blade structure design.

Hydrodynamic Performance of a Newly-Designed Pelagic and Demersal Trawls Using Physical Modeling and Analytical Methods for Cameroonian Industrial Fisheries

This study proposed the newly-designed Pelagic and demersal trawls for the fishing vessels operating in Cameroonian waters in pelagic and demersal fishing grounds. The engineering performances of both trawls were investigated using physical modelling method and analytical method based on the predicted equations. In a flume tank, a series of physical model tests based on Tauti’s law were performed to investigate the hydrodynamic and geometrical performances of both trawls and to assess the applicability of the analytical methods based on predicted equations. The results showed that in model scale, the working towing speed and door spread for the pelagic trawl were 3.5 knots and 1.85 m, respectively, and for the bottom trawl net they were 4.0 knots and 1.8 m. At that speed and door spread, the drag force, net opening height, and wing-end spread of the pelagic model trawl were 36.73 N, 0.89 m, and 0.86 m, respectively, and the swept area was 0.76 m2. Bottom trawl speed and door spread were 30.43 N, 0.38 m, and 0.45 m, respectively, and the swept area was 0.25 m2. The maximum difference between the experimental and analytical results of hydrodynamic performances was less than 56.22% and 41.45%, respectively, for pelagic and bottom trawls, the results of the geometrical performances obtained using predicted equations were close to the experimental results in the flume tank with a maximum relative error less than 12.85%. The newly developed pelagic and bottom trawls had advanced engineering performance for high catch efficiency and selectivity and could be used in commercial fishing operations in Cameroonian waters.

Noise-induced phase transition in the Vicsek model through eigen microstate methodology

This paper presents a comprehensive framework for analyzing phase transitions in collective models such as theVicsek model under various noise types. The Vicsek model, focusing on understanding the collective behaviors of socialanimals, is known due to its discontinuous phase transitions under vector noise. However, its behavior under scalar noiseremains less conclusive. Renowned for its efficacy in the analysis of complex systems under both equilibrium and nonequilibriumstates, the eigen microstate method is employed here for a quantitative examination of the phase transitions inthe Vicsek model under both vector and scalar noises. The study finds that the Vicsek model exhibits discontinuous phasetransitions regardless of noise type. Furthermore, the dichotomy method is utilized to identify the critical points for thesephase transitions. A significant finding is the observed increase in the critical point for discontinuous phase transitions withescalation of population density.

Modularized and Parametric Modeling Technology for Finite Element Simulations of Underground Engineering under Complicated Geological Conditions

The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.

Dynamic Modeling and Experimental Verification of an RPR Type Compliant Paralle Mechanism with Low Orders

Efficiency of calculating a dynamic response is an important point of the compliant mechanism for posture adjustment.Dynamic modeling with low orders of a 2R1T compliant parallel mechanism is studied in the paper.The mechanism with two out-of-plane rotational and one lifting degrees of freedom(DoFs)plays an important role in posture adjustment.Based on elastic beam theory,the stiffness matrix and mass matrix of the beam element are established where the moment of inertia is considered.To improve solving efficiency,a dynamic model with low orders of the mechanism is established based on a modified modal synthesis method.Firstly,each branch of the RPR type mechanism is divided into a substructure.Subsequently,a set of hypothetical modes of each substructure is obtained based on the C-B method.Finally,dynamic equation of the whole mechanism is established by the substructure assembly.A dynamic experiment is conducted to verify the dynamic characteristics of the compliant mechanism.

Novel damage constitutive models and new quantitative identification method for stress thresholds of rocks under uniaxial compression

Four key stress thresholds exist in the compression process of rocks,i.e.,crack closure stress(σ_(cc)),crack initiation stress(σ_(ci)),crack damage stress(σ_(cd))and compressive strength(σ_(c)).The quantitative identifications of the first three stress thresholds are of great significance for characterizing the microcrack growth and damage evolution of rocks under compression.In this paper,a new method based on damage constitutive model is proposed to quantitatively measure the stress thresholds of rocks.Firstly,two different damage constitutive models were constructed based on acoustic emission(AE)counts and Weibull distribution function considering the compaction stages of the rock and the bearing capacity of the damage element.Then,the accumulative AE counts method(ACLM),AE count rate method(CRM)and constitutive model method(CMM)were introduced to determine the stress thresholds of rocks.Finally,the stress thresholds of 9 different rocks were identified by ACLM,CRM,and CMM.The results show that the theoretical stress−strain curves obtained from the two damage constitutive models are in good agreement with that of the experimental data,and the differences between the two damage constitutive models mainly come from the evolutionary differences of the damage variables.The results of the stress thresholds identified by the CMM are in good agreement with those identified by the AE methods,i.e.,ACLM and CRM.Therefore,the proposed CMM can be used to determine the stress thresholds of rocks.

Emergent topological ordered phase for the Ising-XY model revealed by cluster-updating Monte Carlo method

The two-component cold atom systems with anisotropic hopping amplitudes can be phenomenologically described by a two-dimensional Ising-XY coupled model with spatial anisotropy.At low temperatures,theoretical predictions[Phys.Rev.A 72053604(2005)]and[arXiv:0706.1609]indicate the existence of a topological ordered phase characterized by Ising and XY disorder but with 2XY ordering.However,due to ergodic difficulties faced by Monte Carlo methods at low temperatures,this topological phase has not been numerically explored.We propose a linear cluster updating Monte Carlo method,which flips spins without rejection in the anisotropy limit but does not change the energy.Using this scheme and conventional Monte Carlo methods,we succeed in revealing the nature of topological phases with half-vortices and domain walls.In the constructed global phase diagram,Ising and XY-type transitions are very close to each other and differ significantly from the schematic phase diagram reported earlier.We also propose and explore a wide range of quantities,including magnetism,superfluidity,specific heat,susceptibility,and even percolation susceptibility,and obtain consistent and reliable results.Furthermore,we observed first-order transitions characterized by common intersection points in magnetizations for different system sizes,as opposed to the conventional phase transition where Binder cumulants of various sizes share common intersections.The critical exponents of different types of phase transitions are reasonably fitted.The results are useful to help cold atom experiments explore the half-vortex topological phase.

Impacts of Aggregation Methods and Trophospecies Number on the Structure and Function of Marine Food Webs

Aggregation of species with similar ecological properties is one of the effective methods to simplify food web researches.However,species aggregation will affect not only the complexity of modeling process but also the accuracy of models’outputs.Selection of aggregation methods and the number of trophospecies are the keys to study the simplification of food web.In this study,three aggregation methods,including taxonomic aggregation(TA),structural equivalence aggregation(SEA),and self-organizing maps(SOM),were analyzed and compared with the linear inverse model–Markov Chain Monte Carlo(LIM-MCMC)model.Impacts of aggregation methods and trophospecies number on food webs were evaluated based on the robustness and unitless of ecological net-work indices.Results showed that aggregation method of SEA performed better than the other two methods in estimating food web structure and function indices.The effects of aggregation methods were driven by the differences in species aggregation principles,which will alter food web structure and function through the redistribution of energy flow.According to the results of mean absolute percentage error(MAPE)which can be applied to evaluate the accuracy of the model,we found that MAPE in food web indices will increase with the reducing trophospecies number,and MAPE in food web function indices were smaller and more stable than those in food web structure indices.Therefore,trade-off between simplifying food webs and reflecting the status of ecosystem should be con-sidered in food web studies.These findings highlight the importance of aggregation methods and trophospecies number in the analy-sis of food web simplification.This study provided a framework to explore the extent to which food web models are affected by dif-ferent species aggregation,and will provide scientific basis for the construction of food webs.

High-Order Decoupled and Bound Preserving Local Discontinuous Galerkin Methods for a Class of Chemotaxis Models

In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.

Evaluation of soil erosion vulnerability in Hubei Province of China using RUSLE model and combination weighting method

Soil erosion has been recognized as a critical environmental issue worldwide.While previous studies have primarily focused on watershed-scale soil erosion vulnerability from a natural factor perspective,there is a notable gap in understanding the intricate interplay between natural and socio-economic factors,especially in the context of spatial heterogeneity and nonlinear impacts of human-land interactions.To address this,our study evaluates the soil erosion vulnerability at a provincial scale,taking Hubei Province as a case study to explore the combined effects of natural and socio-economic factors.We developed an evaluation index system based on 15 indicators of soil erosion vulnerability:exposure,sensitivity,and adaptability.In addition,the combination weighting method was applied to determine index weights,and the spatial interaction was analyzed using spatial autocorrelation,geographical temporally weighted regression and geographical detector.The results showed an overall decreasing soil erosion intensity in Hubei Province during 2000 and 2020.The soil erosion vulnerability increased before 2000 and then.The areas with high soil erosion vulnerability were mainly confined in the central and southern regions of Hubei Province(Xiantao,Tianmen,Qianjiang and Ezhou)with obvious spatial aggregation that intensified over time.Natural factors(habitat quality index)had negative impacts on soil erosion vulnerability,whereas socio-economic factors(population density)showed substantial spatial variability in their influences.There was a positive correlation between soil erosion vulnerability and erosion intensity,with the correlation coefficients ranging from-0.41 and 0.93.The increase of slope was found to enhance the positive correlation between soil erosion vulnerability and intensity.

Modeling asymmetric fracture mechanics of Mg alloy wire in drawing process

In this study,a numerical analysis was conducted on the ductile fracture of a 2-mm diameter Mg-1Zn-0.5Mn-0.5Sr-0.1Ca alloy wire during drawing.The hexagonally close-packed crystal structure of Mg alloys causes asymmetric fracture behavior,especially in the compression region.The aim of this study is to develop a comprehensive damage model for Mg alloy wire that accurately predicts ductile fracture,with a focus on the compression region.A novel experimental method was introduced to measure the ductile fracture of Mg alloy wires under different stress states.The wire drawing process was simulated using the Generalized Incremental Stress-State dependent damage(GISSMO)Model and the Semi-Analytical Model for Polymers(SAMP)model.The damage model's prediction and the experimental results were found to be in excellent agreement,especially in determining crack initiation.Computational analysis established a safe zone diagram for die angle and reduction ratio,and experimental validation confirmed the feasibility of this approach.The proposed damage model can provide a practical and reliable analysis for optimizing the drawing process of Mg alloy wire.

Numerical Analysis of Bacterial Meningitis Stochastic Delayed Epidemic Model through Computational Methods

Based on theWorld Health Organization(WHO),Meningitis is a severe infection of the meninges,the membranes covering the brain and spinal cord.It is a devastating disease and remains a significant public health challenge.This study investigates a bacterial meningitis model through deterministic and stochastic versions.Four-compartment population dynamics explain the concept,particularly the susceptible population,carrier,infected,and recovered.The model predicts the nonnegative equilibrium points and reproduction number,i.e.,the Meningitis-Free Equilibrium(MFE),and Meningitis-Existing Equilibrium(MEE).For the stochastic version of the existing deterministicmodel,the twomethodologies studied are transition probabilities and non-parametric perturbations.Also,positivity,boundedness,extinction,and disease persistence are studiedrigorouslywiththe helpofwell-known theorems.Standard and nonstandard techniques such as EulerMaruyama,stochastic Euler,stochastic Runge Kutta,and stochastic nonstandard finite difference in the sense of delay have been presented for computational analysis of the stochastic model.Unfortunately,standard methods fail to restore the biological properties of the model,so the stochastic nonstandard finite difference approximation is offered as an efficient,low-cost,and independent of time step size.In addition,the convergence,local,and global stability around the equilibria of the nonstandard computational method is studied by assuming the perturbation effect is zero.The simulations and comparison of the methods are presented to support the theoretical results and for the best visualization of results.

A vector sum analysis method for stability evolution of expansive soil slope considering shear zone damage softening

Slope stability analysis is a classical mechanical problem in geotechnical engineering and engineering geology.It is of great significance to study the stability evolution of expansive soil slopes for engineering construction in expansive soil areas.Most of the existing studies evaluate the slope stability by analyzing the limit equilibrium state of the slope,and the analysis method for the stability evolution considering the damage softening of the shear zone is lacking.In this study,the large deformation shear mechanical behavior of expansive soil was investigated by ring shear test.The damage softening characteristic of expansive soil in the shear zone was analyzed,and a shear damage model reflecting the damage softening behavior of expansive soil was derived based on the damage theory.Finally,by skillfully combining the vector sum method and the shear damage model,an analysis method for the stability evolution of the expansive soil slope considering the shear zone damage softening was proposed.The results show that the shear zone subjected to large displacement shear deformation exhibits an obvious damage softening phenomenon.The damage variable equation based on the logistic function can be well used to describe the shear damage characteristics of expansive soil,and the proposed shear damage model is in good agreement with the ring shear test results.The vector sum method considering the damage softening behavior of the shear zone can be well applied to analyze the stability evolution characteristics of the expansive soil slope.The stability factor of the expansive soil slope decreases with the increase of shear displacement,showing an obvious progressive failure behavior.

Exponentially Convergent Multiscale Finite Element Method

We provide a concise review of the exponentially convergent multiscale finite element method(ExpMsFEM)for efficient model reduction of PDEs in heterogeneous media without scale separation and in high-frequency wave propagation.The ExpMsFEM is built on the non-overlapped domain decomposition in the classical MsFEM while enriching the approximation space systematically to achieve a nearly exponential convergence rate regarding the number of basis functions.Unlike most generalizations of the MsFEM in the literature,the ExpMsFEM does not rely on any partition of unity functions.In general,it is necessary to use function representations dependent on the right-hand side to break the algebraic Kolmogorov n-width barrier to achieve exponential convergence.Indeed,there are online and offline parts in the function representation provided by the ExpMsFEM.The online part depends on the right-hand side locally and can be computed in parallel efficiently.The offline part contains basis functions that are used in the Galerkin method to assemble the stiffness matrix;they are all independent of the right-hand side,so the stiffness matrix can be used repeatedly in multi-query scenarios.