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.

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.

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.

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.

AHermitian C^(2) Differential Reproducing Kernel Interpolation Meshless Method for the 3D Microstructure-Dependent Static Flexural Analysis of Simply Supported and Functionally Graded Microplates

This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C^(2) differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.

Efficient method to calculate the eigenvalues of the Zakharov–Shabat system

A numerical method is proposed to calculate the eigenvalues of the Zakharov–Shabat system based on Chebyshev polynomials. A mapping in the form of tanh(ax) is constructed according to the asymptotic of the potential function for the Zakharov–Shabat eigenvalue problem. The mapping can distribute Chebyshev nodes very well considering the gradient for the potential function. Using Chebyshev polynomials, tanh(ax) mapping, and Chebyshev nodes, the Zakharov–Shabat eigenvalue problem is transformed into a matrix eigenvalue problem. This method has good convergence for the Satsuma–Yajima potential and the convergence rate is faster than the Fourier collocation method. This method is not only suitable for simple potential functions but also converges quickly for a complex Y-shape potential. It can also be further extended to other linear eigenvalue problems.

On the spreading behavior of a droplet on a circular cylinder using the lattice Boltzmann method

The study of a droplet spreading on a circular cylinder under gravity was carried out using the pseudo-potential lattice Boltzmann high-density ratios multiphase model with a non-ideal Peng–Robinson equation of state. The calculation results indicate that the motion of the droplet on the cylinder can be divided into three stages: spreading, sliding, and aggregating.The contact length and contact time of a droplet on a cylindrical surface can be affected by factors such as the wettability gradient of the cylindrical wall, the Bond number, and droplet size. Furthermore, phase diagrams showing the relationship between Bond number, cylinder wall wettability gradient, and contact time as well as maximum contact length for three different droplet sizes are given. A theoretical foundation for additional research into the heat and mass transfer process between the droplet and the cylinder can be established by comprehending the variable rules of maximum contact length and contact time.

Application of the finite analytic numerical method to a flowdependent variational data assimilation

An anisotropic diffusion filter can be used to model a flow-dependent background error covariance matrix,which can be achieved by solving the advection-diffusion equation.Because of the directionality of the advection term,the discrete method needs to be chosen very carefully.The finite analytic method is an alternative scheme to solve the advection-diffusion equation.As a combination of analytical and numerical methods,it not only has high calculation accuracy but also holds the characteristic of the auto upwind.To demonstrate its ability,the one-dimensional steady and unsteady advection-diffusion equation numerical examples are respectively solved by the finite analytic method.The more widely used upwind difference method is used as a control approach.The result indicates that the finite analytic method has higher accuracy than the upwind difference method.For the two-dimensional case,the finite analytic method still has a better performance.In the three-dimensional variational assimilation experiment,the finite analytic method can effectively improve analysis field accuracy,and its effect is significantly better than the upwind difference and the central difference method.Moreover,it is still a more effective solution method in the strong flow region where the advective-diffusion filter performs most prominently.

Stability and Time-Step Constraints of Implicit-Explicit Runge-Kutta Methods for the Linearized Korteweg-de Vries Equation

This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.

A combined method using Lattice Boltzmann Method(LBM)and Finite Volume Method(FVM)to simulate geothermal reservoirs in Enhanced Geothermal System(EGS)

With the development of industrial activities,global warming has accelerated due to excessive emission of CO_(2).Enhanced Geothermal System(EGS)utilizes deep geothermal heat for power generation.Although porous medium theory is commonly employed to model geothermal reservoirs in EGS,Hot Dry Rock(HDR)presents a challenge as it consists of impermeable granite with zero porosity,potentially distorting the physical interpretation.To address this,the Lattice Boltzmann Method(LBM)is employed to simulate CO_(2)flow within geothermal reservoirs and the Finite Volume Method(FVM)to solve the energy conservation equation for temperature distribution.This combined method of LBM and FVM is imple-mented using MATLAB.The results showed that the Reynolds numbers(Re)of 3,000 and 8,000 lead to higher heat extraction rates from geothermal reservoirs.However,higher Re values may accelerate thermal breakthrough,posing challenges to EGS operation.Meanwhile,non-equilibrium of density in fractures becomes more pronounced during the system's life cycle,with non-Darcy's law becoming significant at Re values of 3,000 and 8,000.Density stratification due to buoyancy effects significantly impacts temperature distribution within geothermal reservoirs,with buoyancy effects at Re=100 under gravitational influence being noteworthy.Larger Re values(3,000 and 8,000)induce stronger forced convection,leading to more uniform density distribution.The addition of proppant negatively affects heat transfer performance in geothermal reservoirs,especially in single fractures.Practical engineering considerations should determine the quantity of proppant through detailed numerical simulations.

Simplified quantitative analysis method and its application in the insitu synthesized copper-based azide chips

Copper-based azide(Cu(N_(3))2 or CuN_(3),CA)chips synthesized by in-situ azide reaction and utilized in miniaturized explosive systems has become a hot research topic in recent years.However,the advantages of in-situ synthesis method,including small size and low dosage,bring about difficulties in quantitative analysis and differences in ignition capabilities of CA chips.The aim of present work is to develop a simplified quantitative analysis method for accurate and safe analysis of components in CA chips to evaluate and investigate the corresponding ignition ability.In this work,Cu(N_(3))2 and CuN_(3)components in CA chips were separated through dissolution and distillation by utilizing the difference in solubility and corresponding content was obtained by measuring N_(3)-concentration through spectrophotometry.The spectrophotometry method was optimized by studying influencing factors and the recovery rate of different separation methods was studied,ensuring the accuracy and reproducibility of test results.The optimized method is linear in range from 1.0-25.0 mg/L,with a correlation coefficient R^(2)=0.9998,which meets the requirements of CA chips with a milligram-level content test.Compared with the existing ICP method,component analysis results of CA chips obtained by spectrophotometry are closer to real component content in samples and have satisfactory accuracy.Moreover,as its application in miniaturized explosive systems,the ignition ability of CA chips with different component contents for direct ink writing CL-20 and the corresponding mechanism was studied.This study provided a basis and idea for the design and performance evaluation of CA chips in miniaturized explosive systems.

Inversion method of deflection of the vertical based on SWOT wide-swath altimeter data

The deflection of the vertical(DOV)is the key information in the study of ocean gravity field.However,in most areas,the precision of the prime component of DoV is significantly lower than that of the meridian component.To obtain higher accuracy and resolution of ocean gravity information,researchers have proposed a novel altimeter called the wide-swath altimeter.This altimeter allows for the simultaneous acquisition of high-precision and high-resolution two-dimensional measurements of sea surface height(SSH).In this paper,the Surface Water and Ocean Topography(SWOT)mission with a wide-swath altimeter on board is selected for research.One cycle of swoT sea surface height data is simulated to inverse the DOV in the Arabian Sea(45°E—80°E,0°-30°N),and the inversion results are compared with those of conventional altimeter data.The results demonstrate that the difference between the meridian and prime components derived from the inversion of swoT wide-swath data is minimal,significantly outperforming the inversion results of conventional nadir altimeter data.The advantage of swoT wide-swath altimeter lies in its ability to use the multi-directional geoid slope at any sea surface measurement point to invert the components in the meridian and prime directions.To investigate the impact of this advantage on inversion precision,this paper employs a method to calculate the gradient of the geoid in multiple directions to invert DoV components.The improvement effect of calculating the gradient of the geoid in multiple directions on the precision of DoV component is analyzed.It is found that the accuracy of DoV inversion has significantly improved with the increase of geodetic gradient calculation direction.In addition,the effects of various errors and grid spacing in SwoT wide sea surface height data on the precision of Dov inversion are also analyzed.

Numerical Simulation of the Seismic Response of a Horizontal Storage Tank Based on a SPH-FEM Coupling Method

A coupled numerical calculation method combining smooth particle hydrodynamics(SPH)and the finite element method(FEM)was implemented to investigate the seismic response of horizontal storage tanks.Anumericalmodel of a horizontal storage tank featuring a free liquid surface under seismic action was constructed using the SPH–FEM coupling method.The stored liquid was discretized using SPH particles,while the tank and supports were discretized using the FEM.The interaction between the stored liquid and the tank was simulated by using the meshless particle contact method.Then,the numerical simulation results were compared and analyzed against seismic simulation shaking table test data to validate the method.Subsequently,a series of numerical models,considering different liquid storage volumes and seismic effects,were constructed to obtain time history data of base shear and top center displacement,which revealed the seismic performance of horizontal storage tanks.Numerical simulation results and experimental data showed good agreement,with an error rate of less than 18.85%.And this conformity signifies the rationality of the SPH-FEM coupling method.The base shear and top center displacement values obtained by the coupled SPH-FEM method were only 53.3% to 69.1% of those calculated by the equivalent mass method employed in the current code.As the stored liquid volume increased,the seismic response of the horizontal storage tank exhibited a gradual upward trend,with the seismic response increasing from 73% to 388% for every 35% increase in stored liquid volume.The maximum von Mises stress of the tank and the supports remained below the steel yield strength during the earthquake.The coupled SPH-FEM method holds certain advantages in studying the seismic problems of tanks with complex structural forms,particularly due to the representation of the flow field distribution during earthquakes by involving reservoir fluid participation.