Software systems are vulnerable to security breaches as they expand in complexity and functionality.The confidentiality,integrity,and availability of data are gravely threatened by flaws in a system’s design,implemen...Software systems are vulnerable to security breaches as they expand in complexity and functionality.The confidentiality,integrity,and availability of data are gravely threatened by flaws in a system’s design,implementation,or configuration.To guarantee the durability&robustness of the software,vulnerability identification and fixation have become crucial areas of focus for developers,cybersecurity experts and industries.This paper presents a thorough multi-phase mathematical model for efficient patch management and vulnerability detection.To uniquely model these processes,the model incorporated the notion of the learning phenomenon in describing vulnerability fixation using a logistic learning function.Furthermore,the authors have used numerical methods to approximate the solution of the proposed framework where an analytical solution is difficult to attain.The suggested systematic architecture has been demonstrated through statistical analysis using patch datasets,which offers a solid basis for the research conclusions.According to computational research,learning dynamics improves security response and results in more effective vulnerability management.The suggested model offers a systematic approach to proactive vulnerability mitigation and has important uses in risk assessment,software maintenance,and cybersecurity.This study helps create more robust software systems by increasing patch management effectiveness,which benefits developers,cybersecurity experts,and sectors looking to reduce security threats in a growing digital world.展开更多
This article reviews the application and progress of deep learning in efficient numerical computing methods.Deep learning,as an important branch of machine learning,provides new ideas for numerical computation by cons...This article reviews the application and progress of deep learning in efficient numerical computing methods.Deep learning,as an important branch of machine learning,provides new ideas for numerical computation by constructing multi-layer neural networks to simulate the learning process of the human brain.The article explores the application of deep learning in solving partial differential equations,optimizing problems,and data-driven modeling,and analyzes its advantages in computational efficiency,accuracy,and adaptability.At the same time,this article also points out the challenges faced by deep learning numerical computation methods in terms of computational efficiency,interpretability,and generalization ability,and proposes strategies and future development directions for integrating with traditional numerical methods.展开更多
This review paper provides a comprehensive introduction to various numerical methods for the phase-field model used to simulate the phase separation dynamics of diblock copolymer melts.Diblock copolymer systems form c...This review paper provides a comprehensive introduction to various numerical methods for the phase-field model used to simulate the phase separation dynamics of diblock copolymer melts.Diblock copolymer systems form complex structures at the nanometer scale and play a significant role in various applications.The phase-field model,in particular,is essential for describing the formation and evolution of these structures and is widely used as a tool to effectively predict the movement of phase boundaries and the distribution of phases over time.In this paper,we discuss the principles and implementations of various numerical methodologies for this model and analyze the strengths,limitations,stability,accuracy,and computational efficiency of each method.Traditional approaches such as Fourier spectral methods,finite difference methods and alternating direction explicit methods are reviewed,as well as recent advancements such as the invariant energy quadratization method and the scalar auxiliary variable scheme are also presented.In addition,we introduce examples of the phase-field model,which are fingerprint image restoration and 3D printing.These examples demonstrate the extensive applicability of the reviewed methods and models.展开更多
The feasibility of using a problem-dependent method to solve systems of second order ODEs is corroborated by an eigen-based theory and a methodology to develop such a numerical method is constructed.The key steps of t...The feasibility of using a problem-dependent method to solve systems of second order ODEs is corroborated by an eigen-based theory and a methodology to develop such a numerical method is constructed.The key steps of this methodology are to decouple a system of ODEs of second order into a set of uncoupled ODEs of second order;next,an eigen-dependent method is proposed to approximate the solution of each uncoupled ODE of second order.It is vital to transform all eigen-dependent methods to a problem-dependent method to bypass an Eigen analysis.The development of an eigen-dependent method plays a key role in this methodology so that slow eigenmodes can be accurately integrated while there is no instability or excessive amplitude growth in fast eigenmodes.This can explain why a problem-dependent method can simultaneously combine the explicitness of each step and A-stability.Consequently,huge computational efforts can be saved for solving nonlinear stiff problems.A new family of problem-dependent methods is developed in this work so that the feasibility of the proposed methodology can be affirmed.It has almost the same performance as that of the HHT-αmethod.However,it can save more than 99.5%of CPU demand in approximating a solution for a system of 1000 nonlinear second order ODEs.展开更多
In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error...In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.展开更多
To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based sim...To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.展开更多
As a pyrometallurgical process,circulating fluidized bed(CFB) roasting has good potential for application in desulfurization of high-sulfur bauxite.The gas-solid distribution and reaction during CFB roasting of high-s...As a pyrometallurgical process,circulating fluidized bed(CFB) roasting has good potential for application in desulfurization of high-sulfur bauxite.The gas-solid distribution and reaction during CFB roasting of high-sulfur bauxite were simulated using the computational particle fluid dynamics(CPFD) method.The effect of primary air flow velocity on particle velocity,particle volume distribution,furnace temperature distribution and pressure distribution were investigated.Under the condition of the same total flow of natural gas,the impact of the number of inlets on the desulfurization efficiency,atmosphere mass fraction distribution and temperature distribution in the furnace was further investigated.展开更多
After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir develo...After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir development.Conventional numerical reservoir simulation methodologies typically employ static assumptions for model construction,presuming invariant reservoir geological parameters throughout the development process while neglecting the reservoir’s temporal evolution characteristics.Although such simplifications reduce computational complexity,they introduce substantial descriptive inaccuracies.Therefore,this paper proposes a meshless numerical simulation method for reservoirs that considers time-varying characteristics.This method avoids the meshing in traditional numerical simulation methods.From the fluid flow perspective,the reservoir’s computational domain is discretized into a series of connection units.An influence domain with a certain radius centered on the nodes is selected,and one-dimensional connection units are established between the nodes to achieve the characterization of the flow topology structure of the reservoir.In order to reflect the dynamic evolution of the reservoir’s physical properties during the water injection development process,the time-varying characteristics are incorporated into the formula of the seepage characteristic parameters in the meshless calculation.The change relationship of the permeability under different surface fluxes is considered to update the calculated connection conductivity in real time.By combining with the seepage control equation for solution,a time-varying meshless numerical simulation method is formed.The results show that compared with the numerical simulationmethod of the connection elementmethod(CEM)that only considers static parameters,this method has higher simulation accuracy and can better simulate the real migration and distribution of oil and water in the reservoir.Thismethod improves the accuracy of reservoir numerical simulation and the development effect of oilfields,providing a scientific basis for optimizing the water injection strategy,adjusting the production plan,and extending the effective production cycle of the oilfield.展开更多
Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce f...Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce fatigue failure and even cause unpredictable drilling accidents.Therefore,it is important to study the ViV characteristics of deepwater drilling riser and reveal the main controlling factors for ensuring the safe and efficient operation of deepwater drilling engineering.In this paper,the ViV of deepwater drilling riser is numerically simulated in time domain based on the discrete vortex method(DvM).A hydrodynamic analysis model and governing equation of VIV is proposed with considering the effect of riser motion using DVM and slice method,where the governing equation is solved by Runge-Kutta method.Model validation is performed,which verified the correctness and accuracy of the mechanical model and the solution method.On this basis,the influence of the number of control points,current velocity,riser outer diameter,shear flow and top tension on the ViV characteristics of deepwater drilling risers are discussed in detail.The results show that with the increase of current velocity,the vibration amplitude of deepwater drilling riser decreases obviously,while the vibration frequency increases gradually.However,if the outer diameter of riser increases,the vibration amplitude increases,while the vibration frequency decreases gradually.The top tension also has great influence on the VIV of riser.When the top tension is 1.25 G,the VIV is suppressed to a certain extent.This study has guiding significance for optimal design and engineering control of deepwater drilling riser.展开更多
0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation ph...0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation phases of railway projects(Yan et al.,2023;Chen et al.,2022;Fanos and Pradhan,2018).展开更多
The modeling of crack growth in three-dimensional(3D)space poses significant challenges in rock mechanics due to the complex numerical computation involved in simulating crack propagation and interaction in rock mater...The modeling of crack growth in three-dimensional(3D)space poses significant challenges in rock mechanics due to the complex numerical computation involved in simulating crack propagation and interaction in rock materials.In this study,we present a novel approach that introduces a 3D numerical manifold method(3D-NMM)with a geometric kernel to enhance computational efficiency.Specifically,the maximum tensile stress criterion is adopted as a crack growth criterion to achieve strong discontinuous crack growth,and a local crack tracking algorithm and an angle correction technique are incorporated to address minor limitations of the algorithm in a 3D model.The implementation of the program is carried out in Python,using object-oriented programming in two independent modules:a calculation module and a crack module.Furthermore,we propose feasible improvements to enhance the performance of the algorithm.Finally,we demonstrate the feasibility and effectiveness of the enhanced algorithm in the 3D-NMM using four numerical examples.This study establishes the potential of the 3DNMM,combined with the local tracking algorithm,for accurately modeling 3D crack propagation in brittle rock materials.展开更多
The boundary knot method(BKM)is a simple boundary-type meshless method.Due to the use of non-singular general solutions rather than singular fundamental solutions,BKM does not need to consider the artificial boundary....The boundary knot method(BKM)is a simple boundary-type meshless method.Due to the use of non-singular general solutions rather than singular fundamental solutions,BKM does not need to consider the artificial boundary.Therefore,this method has the merits of purely meshless,easy to program,high solution accuracy and so on.In this paper,we investigate the effectiveness of the BKM for solving Helmholtz-type problems under various conditions through a series of novel numerical experiments.The results demonstrate that the BKM is efficient and achieves high computational accuracy for problems with smooth or continuous boundary conditions.However,when applied to discontinuous boundary problems,the method exhibits significant numerical instability,potentially leading to substantial deviations in the computed results.Finally,three potential improvement strategies are proposed to mitigate this limitation.展开更多
In this work, an efficient spectral method is proposed to solve the fourth-order eigenvalue problem in cylinder domain. Firstly, the key point of this method is to decompose the original model into a kind of decoupled...In this work, an efficient spectral method is proposed to solve the fourth-order eigenvalue problem in cylinder domain. Firstly, the key point of this method is to decompose the original model into a kind of decoupled two-dimensional eigenvalue problem by cylindrical coordinate transformation and Fourier series expansion, and deduce the crucial essential pole conditions. Secondly, we define a kind of weighted Sobolev spaces, and establish a suitable variational formula and its discrete form for each two-dimensional eigenvalue problem. Furthermore, we derive the equivalent operator formulas and obtain some prior error estimates of spectral theory of compact operators. More importantly, we further obtained error estimates for approximating eigenvalues and eigenfunctions by using two newly constructed projection operators. Finally,some numerical experiments are performed to validate our theoretical results and algorithm.展开更多
The protection effectiveness of traditional Lightning Strike Protection(LSP)for composite rotor blade of helicopter can be diminished due to the explosion risk in overlapping attachment under lightning strike,so a new...The protection effectiveness of traditional Lightning Strike Protection(LSP)for composite rotor blade of helicopter can be diminished due to the explosion risk in overlapping attachment under lightning strike,so a new protection method based on Air Breakdown and insulating adhesive layer(AB-LSP method)was designed to avoid it.In this study,a numerical method was developed to simulate the electrical breakdown,and verified by experiment results.Based on this method,a Finite Element Model(FEM)was established to investigate the effect of two factors(breakdown strength and initial ablation temperature of adhesive layer)on the LSP effectiveness.The results show that the breakdown strength impacts more to the ablation damage in composite than that of high-temperature resistance.Then,another FEM was established to predict the ablation damage by lightning strike in the AB-LSP method protected composite rotor blade.The mechanisms and potential key parameters(magnitude of lightning current,discharge channel location,adhesive layer thickness,and air gap width)that could affect the protection effectiveness were analyzed.The introduction of air breakdown changes the current conduction path and reduces explosion risk.After rational design,this method can offer effective lightning protection for composite helicopter rotor blade and other composite structures.展开更多
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 a...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.展开更多
To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’...To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.展开更多
Quadratic matrix equations arise in many elds of scienti c computing and engineering applications.In this paper,we consider a class of quadratic matrix equations.Under a certain condition,we rst prove the existence of...Quadratic matrix equations arise in many elds of scienti c computing and engineering applications.In this paper,we consider a class of quadratic matrix equations.Under a certain condition,we rst prove the existence of minimal nonnegative solution for this quadratic matrix equation,and then propose some numerical methods for solving it.Convergence analysis and numerical examples are given to verify the theories and the numerical methods of this paper.展开更多
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 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.展开更多
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 challeng...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.展开更多
Atmospheric models are physical equations based on the ideal gas law. Applied to the atmosphere, this law yields equations for water, vapor (gas), ice, air, humidity, dryness, fire, and heat, thus defining the model o...Atmospheric models are physical equations based on the ideal gas law. Applied to the atmosphere, this law yields equations for water, vapor (gas), ice, air, humidity, dryness, fire, and heat, thus defining the model of key atmospheric parameters. The distribution of these parameters across the entire planet Earth is the origin of the formation of the climatic cycle, which is a normal climatic variation. To do this, the Earth is divided into eight (8) parts according to the number of key parameters to be defined in a physical representation of the model. Following this distribution, numerical models calculate the constants for the formation of water, vapor, ice, dryness, thermal energy (fire), heat, air, and humidity. These models vary in complexity depending on the indirect trigonometric direction and simplicity in the sum of neighboring models. Note that the constants obtained from the equations yield 275.156˚K (2.006˚C) for water, 273.1596˚K (0.00963˚C) for vapor, 273.1633˚K (0.0133˚C) for ice, 0.00365 in/s for atmospheric dryness, 1.996 in<sup>2</sup>/s for humidity, 2.993 in<sup>2</sup>/s for air, 1 J for thermal energy of fire, and 0.9963 J for heat. In summary, this study aims to define the main parameters and natural phenomena contributing to the modification of planetary climate. .展开更多
基金supported by grants received by the first author and third author from the Institute of Eminence,Delhi University,Delhi,India,as part of the Faculty Research Program via Ref.No./IoE/2024-25/12/FRP.
文摘Software systems are vulnerable to security breaches as they expand in complexity and functionality.The confidentiality,integrity,and availability of data are gravely threatened by flaws in a system’s design,implementation,or configuration.To guarantee the durability&robustness of the software,vulnerability identification and fixation have become crucial areas of focus for developers,cybersecurity experts and industries.This paper presents a thorough multi-phase mathematical model for efficient patch management and vulnerability detection.To uniquely model these processes,the model incorporated the notion of the learning phenomenon in describing vulnerability fixation using a logistic learning function.Furthermore,the authors have used numerical methods to approximate the solution of the proposed framework where an analytical solution is difficult to attain.The suggested systematic architecture has been demonstrated through statistical analysis using patch datasets,which offers a solid basis for the research conclusions.According to computational research,learning dynamics improves security response and results in more effective vulnerability management.The suggested model offers a systematic approach to proactive vulnerability mitigation and has important uses in risk assessment,software maintenance,and cybersecurity.This study helps create more robust software systems by increasing patch management effectiveness,which benefits developers,cybersecurity experts,and sectors looking to reduce security threats in a growing digital world.
文摘This article reviews the application and progress of deep learning in efficient numerical computing methods.Deep learning,as an important branch of machine learning,provides new ideas for numerical computation by constructing multi-layer neural networks to simulate the learning process of the human brain.The article explores the application of deep learning in solving partial differential equations,optimizing problems,and data-driven modeling,and analyzes its advantages in computational efficiency,accuracy,and adaptability.At the same time,this article also points out the challenges faced by deep learning numerical computation methods in terms of computational efficiency,interpretability,and generalization ability,and proposes strategies and future development directions for integrating with traditional numerical methods.
文摘This review paper provides a comprehensive introduction to various numerical methods for the phase-field model used to simulate the phase separation dynamics of diblock copolymer melts.Diblock copolymer systems form complex structures at the nanometer scale and play a significant role in various applications.The phase-field model,in particular,is essential for describing the formation and evolution of these structures and is widely used as a tool to effectively predict the movement of phase boundaries and the distribution of phases over time.In this paper,we discuss the principles and implementations of various numerical methodologies for this model and analyze the strengths,limitations,stability,accuracy,and computational efficiency of each method.Traditional approaches such as Fourier spectral methods,finite difference methods and alternating direction explicit methods are reviewed,as well as recent advancements such as the invariant energy quadratization method and the scalar auxiliary variable scheme are also presented.In addition,we introduce examples of the phase-field model,which are fingerprint image restoration and 3D printing.These examples demonstrate the extensive applicability of the reviewed methods and models.
文摘The feasibility of using a problem-dependent method to solve systems of second order ODEs is corroborated by an eigen-based theory and a methodology to develop such a numerical method is constructed.The key steps of this methodology are to decouple a system of ODEs of second order into a set of uncoupled ODEs of second order;next,an eigen-dependent method is proposed to approximate the solution of each uncoupled ODE of second order.It is vital to transform all eigen-dependent methods to a problem-dependent method to bypass an Eigen analysis.The development of an eigen-dependent method plays a key role in this methodology so that slow eigenmodes can be accurately integrated while there is no instability or excessive amplitude growth in fast eigenmodes.This can explain why a problem-dependent method can simultaneously combine the explicitness of each step and A-stability.Consequently,huge computational efforts can be saved for solving nonlinear stiff problems.A new family of problem-dependent methods is developed in this work so that the feasibility of the proposed methodology can be affirmed.It has almost the same performance as that of the HHT-αmethod.However,it can save more than 99.5%of CPU demand in approximating a solution for a system of 1000 nonlinear second order ODEs.
文摘In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52271317 and 52071149)the Fundamental Research Funds for the Central Universities(HUST:2019kfy XJJS007)。
文摘To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.
基金supported by the National Key Research and Development Program of China(2022YFC2904400)Guangxi Science and Technology Major Project(Gui Ke AA23023033)。
文摘As a pyrometallurgical process,circulating fluidized bed(CFB) roasting has good potential for application in desulfurization of high-sulfur bauxite.The gas-solid distribution and reaction during CFB roasting of high-sulfur bauxite were simulated using the computational particle fluid dynamics(CPFD) method.The effect of primary air flow velocity on particle velocity,particle volume distribution,furnace temperature distribution and pressure distribution were investigated.Under the condition of the same total flow of natural gas,the impact of the number of inlets on the desulfurization efficiency,atmosphere mass fraction distribution and temperature distribution in the furnace was further investigated.
基金funded by the 14th Five-Year Plan Major Science and Technology Project of CNOOC project number KJGG2021-0506.
文摘After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir development.Conventional numerical reservoir simulation methodologies typically employ static assumptions for model construction,presuming invariant reservoir geological parameters throughout the development process while neglecting the reservoir’s temporal evolution characteristics.Although such simplifications reduce computational complexity,they introduce substantial descriptive inaccuracies.Therefore,this paper proposes a meshless numerical simulation method for reservoirs that considers time-varying characteristics.This method avoids the meshing in traditional numerical simulation methods.From the fluid flow perspective,the reservoir’s computational domain is discretized into a series of connection units.An influence domain with a certain radius centered on the nodes is selected,and one-dimensional connection units are established between the nodes to achieve the characterization of the flow topology structure of the reservoir.In order to reflect the dynamic evolution of the reservoir’s physical properties during the water injection development process,the time-varying characteristics are incorporated into the formula of the seepage characteristic parameters in the meshless calculation.The change relationship of the permeability under different surface fluxes is considered to update the calculated connection conductivity in real time.By combining with the seepage control equation for solution,a time-varying meshless numerical simulation method is formed.The results show that compared with the numerical simulationmethod of the connection elementmethod(CEM)that only considers static parameters,this method has higher simulation accuracy and can better simulate the real migration and distribution of oil and water in the reservoir.Thismethod improves the accuracy of reservoir numerical simulation and the development effect of oilfields,providing a scientific basis for optimizing the water injection strategy,adjusting the production plan,and extending the effective production cycle of the oilfield.
基金the financial support from National Key R&D Program of China(Grant number:2024YFC2815100)Natural Science Foundation of China(Grant number:52322110)Beijing Nova Program(Grant number:20230484341).
文摘Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce fatigue failure and even cause unpredictable drilling accidents.Therefore,it is important to study the ViV characteristics of deepwater drilling riser and reveal the main controlling factors for ensuring the safe and efficient operation of deepwater drilling engineering.In this paper,the ViV of deepwater drilling riser is numerically simulated in time domain based on the discrete vortex method(DvM).A hydrodynamic analysis model and governing equation of VIV is proposed with considering the effect of riser motion using DVM and slice method,where the governing equation is solved by Runge-Kutta method.Model validation is performed,which verified the correctness and accuracy of the mechanical model and the solution method.On this basis,the influence of the number of control points,current velocity,riser outer diameter,shear flow and top tension on the ViV characteristics of deepwater drilling risers are discussed in detail.The results show that with the increase of current velocity,the vibration amplitude of deepwater drilling riser decreases obviously,while the vibration frequency increases gradually.However,if the outer diameter of riser increases,the vibration amplitude increases,while the vibration frequency decreases gradually.The top tension also has great influence on the VIV of riser.When the top tension is 1.25 G,the VIV is suppressed to a certain extent.This study has guiding significance for optimal design and engineering control of deepwater drilling riser.
基金supported by the Open Research Fund of Key Laboratory of Geological Hazards on Three Gorges Reservoir Area(China Three Gorges University),Ministry of Education(No.2022KDZ03)the Science and Technology Projects of Yunnan Provincial Science and Technology Department(No.202401AT070328)+1 种基金the Young talents project of“Xingdian Talent Support Program”in Yunnan Province(No.YNWR-QNBJ-2020-019)the Fund Project of China Academy of Railway Sciences Co.,Ltd.(No.2021YJ178)。
文摘0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation phases of railway projects(Yan et al.,2023;Chen et al.,2022;Fanos and Pradhan,2018).
基金supported by the National Natural Science Foundation of China(Grant Nos.42172312 and 52211540395)support from the Institut Universitaire de France(IUF).
文摘The modeling of crack growth in three-dimensional(3D)space poses significant challenges in rock mechanics due to the complex numerical computation involved in simulating crack propagation and interaction in rock materials.In this study,we present a novel approach that introduces a 3D numerical manifold method(3D-NMM)with a geometric kernel to enhance computational efficiency.Specifically,the maximum tensile stress criterion is adopted as a crack growth criterion to achieve strong discontinuous crack growth,and a local crack tracking algorithm and an angle correction technique are incorporated to address minor limitations of the algorithm in a 3D model.The implementation of the program is carried out in Python,using object-oriented programming in two independent modules:a calculation module and a crack module.Furthermore,we propose feasible improvements to enhance the performance of the algorithm.Finally,we demonstrate the feasibility and effectiveness of the enhanced algorithm in the 3D-NMM using four numerical examples.This study establishes the potential of the 3DNMM,combined with the local tracking algorithm,for accurately modeling 3D crack propagation in brittle rock materials.
基金Supported by the Key Scientific Research Plan of Colleges and Universities in Henan Province(23B140006)。
文摘The boundary knot method(BKM)is a simple boundary-type meshless method.Due to the use of non-singular general solutions rather than singular fundamental solutions,BKM does not need to consider the artificial boundary.Therefore,this method has the merits of purely meshless,easy to program,high solution accuracy and so on.In this paper,we investigate the effectiveness of the BKM for solving Helmholtz-type problems under various conditions through a series of novel numerical experiments.The results demonstrate that the BKM is efficient and achieves high computational accuracy for problems with smooth or continuous boundary conditions.However,when applied to discontinuous boundary problems,the method exhibits significant numerical instability,potentially leading to substantial deviations in the computed results.Finally,three potential improvement strategies are proposed to mitigate this limitation.
基金Supported by the National Natural Science Foundation of China(Grant No.12261017)the Scientific Research Foundation of Guizhou University of Finance and Economics(Grant No.2022ZCZX077)。
文摘In this work, an efficient spectral method is proposed to solve the fourth-order eigenvalue problem in cylinder domain. Firstly, the key point of this method is to decompose the original model into a kind of decoupled two-dimensional eigenvalue problem by cylindrical coordinate transformation and Fourier series expansion, and deduce the crucial essential pole conditions. Secondly, we define a kind of weighted Sobolev spaces, and establish a suitable variational formula and its discrete form for each two-dimensional eigenvalue problem. Furthermore, we derive the equivalent operator formulas and obtain some prior error estimates of spectral theory of compact operators. More importantly, we further obtained error estimates for approximating eigenvalues and eigenfunctions by using two newly constructed projection operators. Finally,some numerical experiments are performed to validate our theoretical results and algorithm.
文摘The protection effectiveness of traditional Lightning Strike Protection(LSP)for composite rotor blade of helicopter can be diminished due to the explosion risk in overlapping attachment under lightning strike,so a new protection method based on Air Breakdown and insulating adhesive layer(AB-LSP method)was designed to avoid it.In this study,a numerical method was developed to simulate the electrical breakdown,and verified by experiment results.Based on this method,a Finite Element Model(FEM)was established to investigate the effect of two factors(breakdown strength and initial ablation temperature of adhesive layer)on the LSP effectiveness.The results show that the breakdown strength impacts more to the ablation damage in composite than that of high-temperature resistance.Then,another FEM was established to predict the ablation damage by lightning strike in the AB-LSP method protected composite rotor blade.The mechanisms and potential key parameters(magnitude of lightning current,discharge channel location,adhesive layer thickness,and air gap width)that could affect the protection effectiveness were analyzed.The introduction of air breakdown changes the current conduction path and reduces explosion risk.After rational design,this method can offer effective lightning protection for composite helicopter rotor blade and other composite structures.
基金the Postgraduate Research and Practice Innovation Program of Jiangsu Province(Grant No.KYCX22_0621)the National Natural Science Foundation of China(Grant No.52209130)Jiangsu Funding Program for Excellent Postdoctoral Talent.
文摘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.
文摘To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.
基金Supported by the National Natural Science Foundation of China(12001395)the special fund for Science and Technology Innovation Teams of Shanxi Province(202204051002018)+1 种基金Research Project Supported by Shanxi Scholarship Council of China(2022-169)Graduate Education Innovation Project of Taiyuan Normal University(SYYJSYC-2314)。
文摘Quadratic matrix equations arise in many elds of scienti c computing and engineering applications.In this paper,we consider a class of quadratic matrix equations.Under a certain condition,we rst prove the existence of minimal nonnegative solution for this quadratic matrix equation,and then propose some numerical methods for solving it.Convergence analysis and numerical examples are given to verify the theories and the numerical methods of this paper.
基金supported by the National Science Foundation grant DMS-1818998.
文摘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.
基金Deanship of Research and Graduate Studies at King Khalid University for funding this work through large Research Project under Grant Number RGP2/302/45supported by the Deanship of Scientific Research,Vice Presidency forGraduate Studies and Scientific Research,King Faisal University,Saudi Arabia(Grant Number A426).
文摘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.
文摘Atmospheric models are physical equations based on the ideal gas law. Applied to the atmosphere, this law yields equations for water, vapor (gas), ice, air, humidity, dryness, fire, and heat, thus defining the model of key atmospheric parameters. The distribution of these parameters across the entire planet Earth is the origin of the formation of the climatic cycle, which is a normal climatic variation. To do this, the Earth is divided into eight (8) parts according to the number of key parameters to be defined in a physical representation of the model. Following this distribution, numerical models calculate the constants for the formation of water, vapor, ice, dryness, thermal energy (fire), heat, air, and humidity. These models vary in complexity depending on the indirect trigonometric direction and simplicity in the sum of neighboring models. Note that the constants obtained from the equations yield 275.156˚K (2.006˚C) for water, 273.1596˚K (0.00963˚C) for vapor, 273.1633˚K (0.0133˚C) for ice, 0.00365 in/s for atmospheric dryness, 1.996 in<sup>2</sup>/s for humidity, 2.993 in<sup>2</sup>/s for air, 1 J for thermal energy of fire, and 0.9963 J for heat. In summary, this study aims to define the main parameters and natural phenomena contributing to the modification of planetary climate. .
