Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.

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.

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.

Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

Cascading multi-segment rupture process of the 2023 Turkish earthquake doublet on a complex fault system revealed by teleseismic P wave back projection method

In this study,the vertical components of broadband teleseismic P wave data recorded by China Earthquake Network are used to image the rupture processes of the February 6th,2023 Turkish earthquake doublet via back projection analysis.Data in two frequency bands(0.5-2 Hz and 1-3 Hz)are used in the imaging processes.The results show that the rupture of the first event extends about 200 km to the northeast and about 150 km to the southwest,lasting~90 s in total.The southwestern rupture is triggered by the northeastern rupture,demonstrating a sequential bidirectional unilateral rupture pattern.The rupture of the second event extends approximately 80 km in both northeast and west directions,lasting~35 s in total and demonstrates a typical bilateral rupture feature.The cascading ruptures on both sides also reflect the occurrence of selective rupture behaviors on bifurcated faults.In addition,we observe super-shear ruptures on certain fault sections with relatively straight fault structures and sparse aftershocks.

A smoothed assumed enhanced strain method for frictional contact with constant strain elements

This paper presents a finite element framework for imposing frictional contact conditions on embedded fracture faces,implemented by the constant-strain assumed enhanced strain(AES)method,where penalty method is used to impose both non-penetration constraint and Coulomb’s law of friction.The proposed constant-strain AES method for modeling embedded frictional contact can be cast into an integration algorithm similar to those used in the classical plasticity theory,where displacement jump is calculated from the local traction equilibrium at Gauss point,so the method does not introduce any additional global degrees of freedom.Moreover,constant-strain elements are often desirable in practice because they can be easily created automatically for large-scale engineering applications with complicated geometries.As encountered in other enriched finite element methods for frictional contact,the problem of normal contact pressure oscillations is also observed in the constant-strain AES method.Therefore,we developed a strain-smoothing procedure to effectively mitigate the oscillations.We investigated and verified the proposed AES framework through several numerical examples,and illustrated the capability of this method in solving challenging nonlinear frictional contact problems.

A High-Accuracy Curve Boundary Recognition Method Based on the Lattice Boltzmann Method and Immersed Moving Boundary Method

Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Lattice Boltzmann(LBM)and the immersion boundary method based on solid ratio(IMB)have limitations in identifying custom curved boundaries.Meanwhile,IBM based on velocity correction(IBM-VC)suffers from inaccuracies and numerical instability.Therefore,this study introduces a high-accuracy curve boundary recognition method(IMB-CB),which identifies boundary nodes by moving the search box,and corrects the weighting function in LBM by calculating the solid ratio of the boundary nodes,achieving accurate recognition of custom curve boundaries.In addition,curve boundary image and dot methods are utilized to verify IMB-CB.The findings revealed that IMB-CB can accurately identify the boundary,showing an error of less than 1.8%with 500 lattices.Also,the flow in the custom curve boundary and aerodynamic characteristics of the NACA0012 airfoil are calculated and compared to IBM-VC.Results showed that IMB-CB yields lower lift and drag coefficient errors than IBM-VC,with a 1.45%drag coefficient error.In addition,the characteristic curve of IMB-CB is very stable,whereas that of IBM-VC is not.For the moving boundary problem,LBM-IMB-CB with discrete element method(DEM)is capable of accurately simulating the physical phenomena of multi-moving particle flow in complex curved pipelines.This research proposes a new curve boundary recognition method,which can significantly promote the stability and accuracy of fluid-solid interaction simulations and thus has huge applications in engineering.

Optimization of Random Feature Method in the High-Precision Regime

Machine learning has been widely used for solving partial differential equations(PDEs)in recent years,among which the random feature method(RFM)exhibits spectral accuracy and can compete with traditional solvers in terms of both accuracy and efficiency.Potentially,the optimization problem in the RFM is more difficult to solve than those that arise in traditional methods.Unlike the broader machine-learning research,which frequently targets tasks within the low-precision regime,our study focuses on the high-precision regime crucial for solving PDEs.In this work,we study this problem from the following aspects:(i)we analyze the coeffcient matrix that arises in the RFM by studying the distribution of singular values;(ii)we investigate whether the continuous training causes the overfitting issue;(ii)we test direct and iterative methods as well as randomized methods for solving the optimization problem.Based on these results,we find that direct methods are superior to other methods if memory is not an issue,while iterative methods typically have low accuracy and can be improved by preconditioning to some extent.

The Adjustment Analysis Method of the Active Surface Antenna Based on Convolutional Neural Network

Active surface technique is one of the key technologies to ensure the reflector accuracy of the millimeter/submillimeter wave large reflector antenna.The antenna is complex,large-scale,and high-precision equipment,and its active surfaces are affected by various factors that are difficult to comprehensively deal with.In this paper,based on the advantage of the deep learning method that can be improved through data learning,we propose the active adjustment value analysis method of large reflector antenna based on deep learning.This method constructs a neural network model for antenna active adjustment analysis in view of the fact that a large reflector antenna consists of multiple panels spliced together.Based on the constraint that a single actuator has to support multiple panels(usually 4),an autonomously learned neural network emphasis layer module is designed to enhance the adaptability of the active adjustment neural network model.The classical 8-meter antenna is used as a case study,the actuators have a mean adjustment error of 0.00252 mm,and the corresponding antenna surface error is0.00523 mm.This active adjustment result shows the effectiveness of the method in this paper.

Development and Validation of a UPLC Method for the Determination of Docetaxel and Its Related Substances in Pharmaceutical Dosage Forms, an Antineoplastic Agent

A novel, simple, and sensitive Ultra Performance Liquid Chromatography (UPLC) method was developed and validated for the quantification of process-related impurities and degradants, as well as the assay of Docetaxel. The stability-indicating capability of the method was demonstrated through forced degradation studies and a comprehensive mass balance evaluation. Chromatographic separation was achieved using an ACQUITY UPLC BEH C18 column (100 × 2.1 mm, 1.7 µm), with gradient elution. The mobile phase A comprised a mixture of water, methanol, and acetonitrile (500:300:200, v/v/v), while mobile phase B was acetonitrile and water (800:200, v/v). The flow rate was set at 0.4 mL/min, with detection at 232 nm using a photodiode array detector. The method exhibited excellent performance, with a tailing factor of 1.10 for Docetaxel. The method was rigorously validated for precision, accuracy, linearity, LOD, LOQ, ruggedness, specificity, and robustness. Forced degradation studies confirmed the method’s suitability for stability analysis. Stability testing on the drug substance was conducted following ICH guidelines.

Upper and Lower Bounds of the α-Universal Triple I Method for Unified Interval Implications

Theα-universal triple I(α-UTI)method is a recognized scheme in the field of fuzzy reasoning,whichwas proposed by our research group previously.The robustness of fuzzy reasoning determines the quality of reasoning algorithms to a large extent,which is quantified by calculating the disparity between the output of fuzzy reasoning with interference and the output without interference.Therefore,in this study,the interval robustness(embodied as the interval stability)of theα-UTI method is explored in the interval-valued fuzzy environment.To begin with,the stability of theα-UTI method is explored for the case of an individual rule,and the upper and lower bounds of its results are estimated,using four kinds of unified interval implications(including the R-interval implication,the S-interval implication,the QL-interval implication and the interval t-norm implication).Through analysis,it is found that theα-UTI method exhibits good interval stability for an individual rule.Moreover,the stability of theα-UTI method is revealed in the case of multiple rules,and the upper and lower bounds of its outcomes are estimated.The results show that theα-UTI method is stable for multiple rules when four kinds of unified interval implications are used,respectively.Lastly,theα-UTI reasoning chain method is presented,which contains a chain structure with multiple layers.The corresponding solutions and their interval perturbations are investigated.It is found that theα-UTI reasoning chain method is stable in the case of chain reasoning.Two application examples in affective computing are given to verify the stability of theα-UTImethod.In summary,through theoretical proof and example verification,it is found that theα-UTImethod has good interval robustness with four kinds of unified interval implications aiming at the situations of an individual rule,multi-rule and reasoning chain.

A novel method for determining the resistivity of compressed superconducting materials

The measurement of resistivity in a compressed material within a diamond anvil cell presents significant challenges.The high-pressure exper-imental setup makes it difficult to directly measure the size changes induced by pressure in the three crystallographic directions of the sample.In this study,we introduce a novel and effective method that addresses these technical challenges.This method is anticipated to offer a valuable foundation for high-pressure investigations on quantum materials,particularly those with anisotropic layered structures.

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.

A Hierarchical Method for Locating the Interferometric Fringes of Celestial Sources in the Visibility Data

In source detection in the Tianlai project,locating the interferometric fringe in visibility data accurately will influence downstream tasks drastically,such as physical parameter estimation and weak source exploration.Considering that traditional locating methods are time-consuming and supervised methods require a great quantity of expensive labeled data,in this paper,we first investigate characteristics of interferometric fringes in the simulation and real scenario separately,and integrate an almost parameter-free unsupervised clustering method and seeding filling or eraser algorithm to propose a hierarchical plug and play method to improve location accuracy.Then,we apply our method to locate single and multiple sources’interferometric fringes in simulation data.Next,we apply our method to real data taken from the Tianlai radio telescope array.Finally,we compare with unsupervised methods that are state of the art.These results show that our method has robustness in different scenarios and can improve location measurement accuracy effectively.

Research on the Extraction of OLR Anomaly Prior to Ms 7.5 Sand Point,Alaska Earthquake based on IPI Method

Based on the understanding that the seismic fault system is a nonlinear complex system,Rundle(1995)introduced the nonlinear threshold system used in meteorology to analyze the ocean-atmosphere interface and the El Ni?o Southern Oscillation into the study of seismic activity changes,and then proposed the PI method(Rundle et al.,2000a,b).Wu et al.(2011)modified the Pattern Informatics Method named MPI to extract the ionospheric anomaly by using data from DEMETER satellites which is suitable for 1–3 months short-term prediction.

Numerical Investigations on the Resonance Errors of Multiscale Discontinuous Galerkin Methods for One-Dimensional Stationary Schrödinger Equation

In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al.in J Sci Comput 66:321–345,2016;Dong and Wang in J Comput Appl Math 380:1–11,2020)for a one-dimensional stationary Schrödinger equation.Previous work showed that penalty parameters were required to be positive in error analysis,but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes.In this work,by performing extensive numerical experiments,we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods,and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.

The Cumulative Method for Multiplication and Division

This paper provides a method of the process of computation called the cumulative method, it is based upon repeated cumulative process. The cumulative method is being adapted to the purposes of computation, particularly multiplication and division. The operations of multiplication and division are represented by algebraic formulas. An advantage of the method is that the cumulative process can be performed on decimal numbers. The present paper aims to establish a basic and useful formula valid for the two fundamental arithmetic operations of multiplication and division. The new cumulative method proved to be more flexible and made it possible to extend the multiplication and division based on repeated addition/subtraction to decimal numbers.

Efficient Preparation and Electrochemical Properties of BiOCl/Graphite by One-step Solvothermal Method

We used the surface-pretreated graphite paper(Gp)as a carrier and loaded BiOCl with high selectivity to Cl^(-)on its surface by solvothermal method to form BiOCl@Gp electrode.The morphology,structure,and composition of the materials were characterized by scanning electron microscopy and nitrogen adsorption/desorption,and the results showed that the spherical BiOCl particles were uniformly dispersed on the surface of the Gp,forming a mesoporous BiOCl@Gp composite with a specific surface area of 22.82 m^(2)/g and a pore volume of 0.043 cm3/g.Furthermore,cyclic voltammetry and electrochemical impedance spectroscopy were used to test the electrochemical properties of the composites,and the stability of BiOCl and the high conductivity of Gp were synergistic,the BiOCl@Gp exhibited a specific capacitance of 30.2 F·g^(-1) at a current density of 0.5 A·g^(-1),and the selectivity of the BiOCl@Gp materials for Cl^(-)was significantly higher than that of SO_(4)^(2-),NO_(2)^(-),and HCO_(3)^(-).Therefore,BiOCl@Gp composite electrode materials can be used for the selective adsorption of Cl^(-)in wastewater,in order to achieve efficient wastewater recycling.

Analytical Method for the Wave Diffraction of Asymmetrically Arranged Breakwaters

The layout forms of several breakwater structures can be generalized as asymmetrical arrangements in actual engineering.However,the problem of wave diffraction around asymmetrically arranged breakwaters has not been adequately investigated.In this study,we propose an analytical method of wave diffraction for regular waves passing through asymmetrically arranged breakwaters,and we use the Nyström method to obtain the analytical solution numerically.We compared the results of this method with those of previous analytical solutions and with numerical results to demonstrate the validity of our approach.We also provided diffraction coefficient diagrams of breakwaters with different layout forms.Moreover,we described the analytical expression for the problem of diffraction through long-wave incident breakwaters and presented an analysis of the relationship between the diffraction coefficients and the widths of breakwater gates.The analytical method presented in this study contributes to the limited literature on the theory of wave diffraction through asymmetrically arranged breakwaters.

Synthesis of Reviews on Auscultation, Approaches, and Methods for Engineering Structures

Topometric auscultation is used to monitor the durability of structures, measure deformations linked to the structure of a structure or to the movement of the ground over a part of the globe, set up warning systems, etc. It first appeared as a visual method and rapidly evolved through the various techniques used. Some of these techniques using topography are used in several fields (civil engineering, geodesy, topography, mechanics, nuclear engineering, hydraulics, physics, etc.). These topometric techniques have undergone major changes as a result of technological advances, growing needs in the monitoring of movements or deformations, increased requirements and new challenges. The methodology adopted depends on the measuring instrument used, the parameters to be estimated and access to the area to be measured. There are two types of methods: destructive and non-destructive. In addition to the visual method, they can also be classified as mechanical, physico-chemical, dynamometric, electrophysical and geometric. The estimated parameter varies according to the methodology adopted. It can be defined by coordinates, distances, potential, electrical resistance, etc.