Stability Analysis and Performance Evaluation of Additive Mixed-Precision Runge-Kutta Methods

Additive Runge-Kutta methods designed for preserving highly accurate solutions in mixed-precision computation were previously proposed and analyzed.These specially designed methods use reduced precision for the implicit computations and full precision for the explicit computations.In this work,we analyze the stability properties of these methods and their sensitivity to the low-precision rounding errors,and demonstrate their performance in terms of accuracy and efficiency.We develop codes in FORTRAN and Julia to solve nonlinear systems of ODEs and PDEs using the mixed-precision additive Runge-Kutta(MP-ARK)methods.The convergence,accuracy,and runtime of these methods are explored.We show that for a given level of accuracy,suitably chosen MP-ARK methods may provide significant reductions in runtime.

Method of fabricating artificial rock specimens based on extrusion free forming(EFF)3D printing

Three-dimensional(3D)printing technology has been widely used to create artificial rock samples in rock mechanics.While 3D printing can create complex fractures,the material still lacks sufficient similarity to natural rock.Extrusion free forming(EFF)is a 3D printing technique that uses clay as the printing material and cures the specimens through high-temperature sintering.In this study,we attempted to use the EFF technology to fabricate artificial rock specimens.The results show the physico-mechanical properties of the specimens are significantly affected by the sintering temperature,while the nozzle diameter and layer thickness also have a certain impact.The specimens are primarily composed of SiO_(2),with mineral compositions similar to that of natural rocks.The density,uniaxial compressive strength(UCS),elastic modulus,and tensile strength of the printed specimens fall in the range of 1.65–2.54 g/cm3,16.46–50.49 MPa,2.17–13.35 GPa,and 0.82–17.18 MPa,respectively.It is capable of simulating different types of rocks,especially mudstone,sandstone,limestone,and gneiss.However,the simulation of hard rocks with UCS exceeding 50 MPa still requires validation.

Recommendation Method for Contrastive Enhancement of Neighborhood Information

Knowledge graph can assist in improving recommendation performance and is widely applied in various person-alized recommendation domains.However,existing knowledge-aware recommendation methods face challenges such as weak user-item interaction supervisory signals and noise in the knowledge graph.To tackle these issues,this paper proposes a neighbor information contrast-enhanced recommendation method by adding subtle noise to construct contrast views and employing contrastive learning to strengthen supervisory signals and reduce knowledge noise.Specifically,first,this paper adopts heterogeneous propagation and knowledge-aware attention networks to obtain multi-order neighbor embedding of users and items,mining the high-order neighbor informa-tion of users and items.Next,in the neighbor information,this paper introduces weak noise following a uniform distribution to construct neighbor contrast views,effectively reducing the time overhead of view construction.This paper then performs contrastive learning between neighbor views to promote the uniformity of view information,adjusting the neighbor structure,and achieving the goal of reducing the knowledge noise in the knowledge graph.Finally,this paper introduces multi-task learning to mitigate the problem of weak supervisory signals.To validate the effectiveness of our method,experiments are conducted on theMovieLens-1M,MovieLens-20M,Book-Crossing,and Last-FM datasets.The results showthat compared to the best baselines,our method shows significant improvements in AUC and F1.

Dynamic Coupling Analysis of Semisubmersible Platform Float-over Method for Docking Case

In this paper,the multi-body coupled dynamic characteristics of a semisubmersible platform and an HYSY 229 barge were investigated.First,coupled hydrodynamic analysis of the HYSY 229 barge and the semisubmersible platform was performed.Relevant hydrodynamic parameters were obtained using the retardation function method of three-dimensional frequency-domain potential flow theory.The results of the hydrodynamic analysis were highly consistent with the test findings,verifying the accuracy of the multifloating hydrodynamic coupling analysis,and key hydrodynamic parameters were solved for different water depths and the coupling effect.According to the obtained results,the hydrodynamic influence was the largest in shallow waters when the coupling effect was considered.Furthermore,the coupled motion equation combined with viscous damping,fender system,and mooring system was established,and the hydrodynamics,floating body motion,and dynamic response of the fender system were analyzed.Motion analysis revealed good agreement among the surge,sway,and yaw motions of the two floating bodies.However,when the wave period reached 10 s,the motion of the two floating bodies showed severe shock,and a relative motion was also observed.Therefore,excessive constraints should be added between the two floating bodies during construction to ensure construction safety.The numerical analysis and model test results of the semisubmersible platform and HYSY 229 barge at a water depth of 42 m and sea conditions of 0°,45°,and 90° were in good agreement,and the error was less than 5%.The maximum movement of the HYSY 229 barge reached 2.61 m in the sway direction,whereas that of the semisubmersible platform was 2.11 m.During construction,excessive constraints should be added between the two floating bodies to limit their relative movement and ensure construction safety.

Numerical Simulation of Surrounding Rock Deformation and Grouting Reinforcement of Cross-Fault Tunnel under Different Excavation Methods

Tunnel construction is susceptible to accidents such as loosening, deformation, collapse, and water inrush, especiallyunder complex geological conditions like dense fault areas. These accidents can cause instability and damageto the tunnel. As a result, it is essential to conduct research on tunnel construction and grouting reinforcementtechnology in fault fracture zones to address these issues and ensure the safety of tunnel excavation projects. Thisstudy utilized the Xianglushan cross-fault tunnel to conduct a comprehensive analysis on the construction, support,and reinforcement of a tunnel crossing a fault fracture zone using the three-dimensional finite element numericalmethod. The study yielded the following research conclusions: The excavation conditions of the cross-fault tunnelarray were analyzed to determine the optimal construction method for excavation while controlling deformationand stress in the surrounding rock. The middle partition method (CD method) was found to be the most suitable.Additionally, the effects of advanced reinforcement grouting on the cross-fault fracture zone tunnel were studied,and the optimal combination of grouting reinforcement range (140°) and grouting thickness (1m) was determined.The stress and deformation data obtained fromon-site monitoring of the surrounding rock was slightly lower thanthe numerical simulation results. However, the change trend of both sets of data was found to be consistent. Theseresearch findings provide technical analysis and data support for the construction and design of cross-fault tunnels.

Energy Stable Nodal DG Methods for Maxwell’s Equations of Mixed-Order Form in Nonlinear Optical Media

In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.

Research on the design method for uniform wear of shield cutters in sand-pebble strata

During shield tunneling in highly abrasive formations such as sand–pebble strata,nonuniform wear of shield cutters is inevitable due to the different cutting distances.Frequent downtimes and cutter replacements have become major obstacles to long-distance shield driving in sand–pebble strata.Based on the cutter wear characteristics in sand–pebble strata in Beijing,a design methodology for the cutterhead and cutters was established in this study to achieve uniform wear of all cutters by the principle of frictional wear.The applicability of the design method was verified through three-dimensional simulations using the engineering discrete element method.The results show that uniform wear of all cutters on the cutterhead could be achieved by installing different numbers of cutters on each trajectory radius and designing a curved spoke with a certain arch height according to the shield diameter.Under the uniform wear scheme,the cutter wear coefficient is greatly reduced,and the largest shield driving distance is increased by approximately 47%over the engineering scheme.The research results indicate that the problem of nonuniform cutter wear in shield excavation could be overcome,thereby providing guiding significance for theoretical innovation and construction of long-distance shield excavation in highly abrasive strata.

Method of Monitoring Three-dimensional Mining Surface Deformation Based on InSAR

In order to solve the problems of small monitoring range,long time and high cost of existing sedimentation observation methods,based on two-view sentinel No.1 radar images of Guqiao mining area in Huainan City from November 4,2017 to November 28,2017,surface change information was obtained in combination with D-InSAR,and the three-dimensional surface deformation was monitored by two-pass method and single line of sight D-InSAR method.The results show that during the research period of 24 d,the maximum deformation of the mining area reached 71 mm,and the southern subsidence was the most obvious,which was in line with the mining subsidence law.The maximum displacement from the north to the south was about 250 mm,while the maximum displacement from the east to the west was about 80 mm,and the maximum subsidence in the center was 110 mm.It is concluded that D-InSAR technique has a good effect on the inversion of the mining subsidence,and this method is suitable for three-dimensional surface monitoring in areas with similar geological conditions.The monitoring results have certain reference value.

Application of Elzaki Transform Method to Market Volatility Using the Black-Scholes Model

Black-Scholes Model (B-SM) simulates the dynamics of financial market and contains instruments such as options and puts which are major indices requiring solution. B-SM is known to estimate the correct prices of European Stock options and establish the theoretical foundation for Option pricing. Therefore, this paper evaluates the Black-Schole model in simulating the European call in a cash flow in the dependent drift and focuses on obtaining analytic and then approximate solution for the model. The work also examines Fokker Planck Equation (FPE) and extracts the link between FPE and B-SM for non equilibrium systems. The B-SM is then solved via the Elzaki transform method (ETM). The computational procedures were obtained using MAPLE 18 with the solution provided in the form of convergent series.

Optimizing Bucket Elevator Performance through a Blend of Discrete Element Method, Response Surface Methodology, and Firefly Algorithm Approaches

This research introduces a novel approach to enhancing bucket elevator design and operation through the integration of discrete element method(DEM)simulation,design of experiments(DOE),and metaheuristic optimization algorithms.Specifically,the study employs the firefly algorithm(FA),a metaheuristic optimization technique,to optimize bucket elevator parameters for maximizing transport mass and mass flow rate discharge of granular materials under specified working conditions.The experimental methodology involves several key steps:screening experiments to identify significant factors affecting bucket elevator operation,central composite design(CCD)experiments to further explore these factors,and response surface methodology(RSM)to create predictive models for transport mass and mass flow rate discharge.The FA algorithm is then applied to optimize these models,and the results are validated through simulation and empirical experiments.The study validates the optimized parameters through simulation and empirical experiments,comparing results with DEM simulation.The outcomes demonstrate the effectiveness of the FA algorithm in identifying optimal bucket parameters,showcasing less than 10%and 15%deviation for transport mass and mass flow rate discharge,respectively,between predicted and actual values.Overall,this research provides insights into the critical factors influencing bucket elevator operation and offers a systematic methodology for optimizing bucket parameters,contributing to more efficient material handling in various industrial applications.

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.

Analyzing the Form-Finding of a Large-Span Transversely Stiffened Suspended Cable System: A Method Considering Construction Processes

The precise control of the shape of transversely stiffened suspended cable systems is crucial. However, existing form-finding methods primarily rely on iterative calculations that treat loads as fixed known conditions. These methods are inefficient and fail to accurately control shape results. In this study, we propose a form-finding method that analyzes the load response of models under different sag and stress levels, taking into account the construction process. To analyze the system, a structural finite element model was established in ANSYS, and geometric nonlinear analysis was conducted using the Newton-Raphson method. The form-finding analysis results demonstrate that the proposed method achieves precise control of shape, with a maximum shape error ranging from 0.33% to 0.98%. Furthermore, the relationships between loads and tension forces are influenced by the deformed shape of the structures, exhibiting significant geometric nonlinear characteristics. Meanwhile, the load response analysis reveals that the stress level of the self-equilibrium state in the transversely stiffened suspended cable system is primarily governed by strength criteria, while shape is predominantly controlled by stiffness criteria. Importantly, by simulating the initial tensioning process as an initial condition, this method solves for a counterweight that satisfies the requirements and achieves a self-equilibrium state with the desired shape. The shape of the self-equilibrium state is precisely controlled by simulating the construction process. Overall, this work presents a new method for analyzing the form-finding process of large-span transversely stiffened suspended cable system, considering the construction process which was often overlooked in previous studies.

Research trends in methods for controlling macro-micro motion platforms

With ongoing economic,scientific,and technological developments,the electronic devices used in daily lives are developing toward precision and miniaturization,and so the demand for high-precision manufacturing machinery is expanding.The most important piece of equipment in modern high-precision manufacturing is the macro-micro motion platform(M3P),which offers high speed,precision,and efficiency and has macro-micro motion coupling characteristics due to its mechanical design and composition of its driving components.Therefore,the design of the control system is crucial for the overall precision of the platform;conventional proportional–integral–derivative control cannot meet the system requirements,and so M3Ps are the subject of a growing range of modern control strategies.This paper begins by describing the development history of M3Ps,followed by their platform structure and motion control system components,and then in-depth assessments of the macro,micro,and macro-micro control systems.In addition to examining the advantages and disadvantages of current macro-micro motion control,recent technological breakthroughs are noted.Finally,based on existing problems,future directions for M3P control systems are given,and the present conclusions offer guidelines for future work on M3Ps.

The dynamic relaxation form finding method aided with advanced recurrent neural network

How to establish a self‐equilibrium configuration is vital for further kinematics and dynamics analyses of tensegrity mechanism.In this study,for investigating tensegrity form‐finding problems,a concise and efficient dynamic relaxation‐noise tolerant zeroing neural network(DR‐NTZNN)form‐finding algorithm is established through analysing the physical properties of tensegrity structures.In addition,the non‐linear constrained opti-misation problem which transformed from the form‐finding problem is solved by a sequential quadratic programming algorithm.Moreover,the noise may produce in the form‐finding process that includes the round‐off errors which are brought by the approximate matrix and restart point calculating course,disturbance caused by external force and manufacturing error when constructing a tensegrity structure.Hence,for the purpose of suppressing the noise,a noise tolerant zeroing neural network is presented to solve the search direction,which can endow the anti‐noise capability to the form‐finding model and enhance the calculation capability.Besides,the dynamic relaxation method is contributed to seek the nodal coordinates rapidly when the search direction is acquired.The numerical results show the form‐finding model has a huge capability for high‐dimensional free form cable‐strut mechanisms with complicated topology.Eventually,comparing with other existing form‐finding methods,the contrast simulations reveal the excellent anti‐noise performance and calculation capacity of DR‐NTZNN form‐finding algorithm.

Selecting the Technology of Sodium Silicate Modified Poplar with the Highest Performance by Fuzzy Orthogonal Method

Sodium silicate modification can improve the overall performance of wood.The modification process has a great influence on the properties of modified wood.In this study,a new method was introduced to analyze the wood modification process,and the properties of modified wood were studied.Poplar wood was modified with sodium silicate by vacuum-pressure impregnation.After screening using single-factor experiments,an orthogonal experiment was carried out with solution concentration,impregnation time,impregnation pressure,and the cycle times as experimental factors.The modified poplar with the best properties was selected by fuzzy mathematics and characterized by SEM,FT-IR,XRD and TG.The results showed that some lignin and hemicellulose were removed from the wood due to the alkaline action of sodium silicate,and the orderly crystal area of poplar became disorderly,resulting in the reduction of crystallinity of the modified poplar wood.FT-IR analysis showed that sodium silicate was hydrolyzed to form polysilicic acid in wood,and structural analysis revealed the formation of Si-O-Si and Si-O-C,indicating that sodium silicate reacted with fibers on the wood cell wall.TG-DTG curves showed that the final residual mass of modified poplar wood increased from 25%to 67%,and the temperature of the maximum loss rate decreased from 343℃ to 276℃.The heat release and smoke release of modified poplar wood decreased obviously.This kind of material with high strength and fire resistance can be used in the outdoor building and indoor furniture.

Numerical simulation on the seismic performance of retrofitted masonry walls based on the combined finite-discrete element method

Due to the long construction life,improper design methods,brittle material properties and poor construction techniques,most existing masonry structures do not perform well during earthquakes.The retrofitting method using an external steel-meshed mortar layer is widely used to retrofit existing masonry buildings.Assessing the seismic performance of masonry walls reinforced by an external steel-meshed mortar layer reasonably and effectively is a difficult subject in the research field of masonry structures.Based on the combined finite-discrete elements method,the numerical models of retrofitted brick walls with four different masonry mortar strengths by an external mortar layer are established.The shear strength of mortar and the contact between the retrofitted mortar layer and the brick blocks are discussed in detail.The failure patterns and load-displacement curves of the retrofitted brick walls were obtained by applying low cycle reciprocating loads to the numerical model,and the bearing capacity and the failure mechanism of the retrofitted walls were obtained by comparing the failure patterns,ultimate bearing capacity,deformability and other aspects with the tests.This study provides a basis for improving the seismic strengthening design method of masonry structures and helps to better assess the seismic performance of masonry structures after retrofitting.

Local Buckling-Induced Forming Method to Produce Metal Bellows

A novel buckling-induced forming method is proposed to produce metal bellows.The tube billet is firstly treated by local heating and cooling,and the axial loading is applied on both ends of the tube,then the buckling occurs at the designated position and forms a convolution.In this paper,a forming apparatus is designed and developed to produce both discontinuous and continuous bellows of 304 stainless steel,and their characteristics are discussed respectively.Furthermore,the influences of process parameters and geometric parameters on the final convolution profile are deeply studied based on FEM analysis.The results suggest that the steel bellows fabricated by the presented buckling-induced forming method have a uniform shape and no obvious reduction of wall thickness.Meanwhile,the forming force required in the process is quite small.

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.

A Modified Formulation of Singular Boundary Method for Exterior Acoustics

This paper proposes amodified formulation of the singular boundarymethod(SBM)by introducing the combined Helmholtz integral equation formulation(CHIEF)and the self-regularization technique to exterior acoustics.In the SBM,the concept of the origin intensity factor(OIF)is introduced to avoid the singularities of the fundamental solutions.The SBM belongs to the meshless boundary collocation methods.The additional use of the CHIEF scheme and the self-regularization technique in the SBM guarantees the unique solution of the exterior acoustics accurately and efficiently.Consequently,by using the SBM coupled with the CHIEF scheme and the self-regularization technique,the accuracy of the numerical solution can be improved,especially near the corresponding internal characteristic frequencies.Several numerical examples of two-dimensional and threedimensional benchmark examples about exterior acoustics are used to verify the effectiveness and accuracy of the proposed method.The proposed numerical results are compared with the analytical solutions and the solutions obtained by the other numerical 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.