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.展开更多
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.展开更多
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.展开更多
Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The t...Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.展开更多
In this investigation,a hybrid approach integrating the IDDES turbulence model and FW-H is employed to forecast the hydroacoustic of the rim driven thruster(RDT)under non-cavitation and uniform flow conditions at vary...In this investigation,a hybrid approach integrating the IDDES turbulence model and FW-H is employed to forecast the hydroacoustic of the rim driven thruster(RDT)under non-cavitation and uniform flow conditions at varying loading conditions(J=0.3 and J=0.6).It is revealed that the quadrupole term contribution in the P-FWH method significantly affects the monopole term in the low-frequency region,while it mainly affects the dipole term in the high-frequency region.Specifically,the overall sound pressure levels(SPL)of the RDT using the P-FWH method are 2.27 dB,10.03 dB,and 16.73 dB at the receiving points from R1 to R3 under the heavy-loaded condition,while they increase by 0.67 dB at R1,and decrease by 14.93 dB at R2,and 22.20 dB at R3,for the light-loaded condition.The study also utilizes the pressure-time derivatives to visualize the numerical noise and to pinpoint the dynamics of the vortex cores,and the optimization of the grid design can significantly reduce the numerical noise.The computational accuracy of the P-FWH method can meet the noise requirements for the preliminary design of rim driven thrusters.展开更多
In order to predict the long-term rutting of asphalt pavement, the effective temperature for pavement rutting is calculated using the numerical simulation method. The transient temperature field of asphalt pavement wa...In order to predict the long-term rutting of asphalt pavement, the effective temperature for pavement rutting is calculated using the numerical simulation method. The transient temperature field of asphalt pavement was simulated based on actual meteorological data of Nanjing. 24-hour rutting development under a transient temperature field was calculated in each month. The rutting depth accumulated under the static temperature field was also estimated and the relationship between constant temperature parameters was analyzed. Then the effective temperature for pavement rutting was determined based on the rutting equivalence principle. The results show that the monthly effective temperature is above 40 t in July and August, while in June and September it ranges from 30 to 40 Rutting development can be ignored when the monthly effective temperature is less than 30 t. The yearly effective temperature for rutting in Nanjing is around 38. 5 t. The long-term rutting prediction model based on the effective temperature can reflect the influences of meteorological factors and traffic time distribution.展开更多
In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step ...In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.展开更多
Built on the integral formulas in Part I,numerical methods are developed for computing velocity potential and streamfunction in a limited domain.When there is no inner boundary（around a data hole） inside the domain,...Built on the integral formulas in Part I,numerical methods are developed for computing velocity potential and streamfunction in a limited domain.When there is no inner boundary（around a data hole） inside the domain,the total solution is the sum of the internally and externally induced parts.For the internally induced part,three numerical schemes（grid-staggering,local-nesting and piecewise continuous integration） are designed to deal with the singularity of the Green＇s function encountered in numerical calculations.For the externally induced part,by setting the velocity potential（or streamfunction） component to zero,the other component of the solution can be computed in two ways：（1） Solve for the density function from its boundary integral equation and then construct the solution from the boundary integral of the density function.（2） Use the Cauchy integral to construct the solution directly.The boundary integral can be discretized on a uniform grid along the boundary.By using local-nesting（or piecewise continuous integration）,the scheme is refined to enhance the discretization accuracy of the boundary integral around each corner point（or along the entire boundary）.When the domain is not free of data holes,the total solution contains a data-hole-induced part,and the Cauchy integral method is extended to construct the externally induced solution with irregular external and internal boundaries.An automated algorithm is designed to facilitate the integrations along the irregular external and internal boundaries.Numerical experiments are performed to evaluate the accuracy and efficiency of each scheme relative to others.展开更多
Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlight...Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.展开更多
One of the critical aspects in mine design is slope stability analysis and the determination of stable slopes. In the Chador- Malu iron ore mine, one of the most important iron ore mines in central Iran, it was consid...One of the critical aspects in mine design is slope stability analysis and the determination of stable slopes. In the Chador- Malu iron ore mine, one of the most important iron ore mines in central Iran, it was considered vital to perform a comprehensive slope stability analysis. At first, we divided the existing rock hosting pit into six zones and a geotechnical map was prepared. Then, the value of MRMR （Mining Rock Mass Rating） was determined for each zone. Owing to the fact that the Chador-Malu iron ore mine is located in a highly tectonic area and the rock mass completely crushed, the Hoek-Brown failure criterion was found suitable to estimate geo-mechanical parameters. After that, the value of cohesion （c） and friction angle （tp） were calculated for different geotechnical zones and relative graphs and equations were derived as a function of slope height. The stability analyses using numerical and limit equilibrium methods showed that some instability problems might occur by increasing the slope height. Therefore, stable slopes for each geotechnical zone and prepared sections were calculated and presented as a function of slope height.展开更多
Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concret...Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concrete pillars (CP) were carried out. During the construction of the second CP, in the space between the remaining pillars, one of the pillars failed leading to the progressive failure of other pillars until 4 000 m 2 of mine had collapsed within a few minutes. In this work, this phenomenon is described by applying both numerical and empirical methods and the respective results are compared. The results of numerical modelling are found to be closer to the actual condition than those of the empirical method. Also, a width-to-height (W/H) ratio less than 1, an inadequate support system and the absence of a detailed program for pillar recovery are shown to be the most important causes of the Domino failure in this mine.展开更多
With the development of modern science and technology, especially computer science, the numerical simulation method has been widely used in material hot-working. Mary achievements have been made in this field by using...With the development of modern science and technology, especially computer science, the numerical simulation method has been widely used in material hot-working. Mary achievements have been made in this field by using the numerical simulation method. The numerical simulation method, especially finite element method fully described in this paper.Applications of the numerical simulation method in material hot-working are also discussed. Finally, the future of the numerical simulation method is outlined.展开更多
Several parameter identification methods of thermal response test were evaluated through numerical and experimental study.A three-dimensional finite-volume numerical model was established under the assumption that the...Several parameter identification methods of thermal response test were evaluated through numerical and experimental study.A three-dimensional finite-volume numerical model was established under the assumption that the soil thermal conductivity had been known in the simulation of thermal response test.The thermal response curve was firstly obtained through numerical calculation.Then,the accuracy of the numerical model was verified with measured data obtained through a thermal response test.Based on the numerical and experimental thermal response curves,the thermal conductivity of the soil was calculated by different parameter identification methods.The calculated results were compared with the assumed value and then the accuracy of these methods was evaluated.Furthermore,the effects of test time,variable data quality,borehole radius,initial ground temperature,and heat injection rate were analyzed.The results show that the method based on cylinder-source model has a low precision and the identified thermal conductivity decreases with an increase in borehole radius.For parameter estimation,the measuring accuracy of the initial temperature of the deep ground soil has greater effect on identified thermal conductivity.展开更多
Model error is one of the key factors restricting the accuracy of numerical weather prediction （NWP）. Considering the continuous evolution of the atmosphere, the observed data （ignoring the measurement error） can ...Model error is one of the key factors restricting the accuracy of numerical weather prediction （NWP）. Considering the continuous evolution of the atmosphere, the observed data （ignoring the measurement error） can be viewed as a series of solutions of an accurate model governing the actual atmosphere. Model error is represented as an unknown term in the accurate model, thus NWP can be considered as an inverse problem to uncover the unknown error term. The inverse problem models can absorb long periods of observed data to generate model error correction procedures. They thus resolve the deficiency and faultiness of the NWP schemes employing only the initial-time data. In this study we construct two inverse problem models to estimate and extrapolate the time-varying and spatial-varying model errors in both the historical and forecast periods by using recent observations and analogue phenomena of the atmosphere. Numerical experiment on Burgers＇ equation has illustrated the substantial forecast improvement using inverse problem algorithms. The proposed inverse problem methods of suppressing NWP errors will be useful in future high accuracy applications of NWP.展开更多
It is believed that the microseismicity induced by mining effect and gas gradient disturbance stress is a precursor to the essential characteristics of roadway unstability. In order to effectively identify and evaluat...It is believed that the microseismicity induced by mining effect and gas gradient disturbance stress is a precursor to the essential characteristics of roadway unstability. In order to effectively identify and evaluate the stability of coal roadways in the process of mine development and extraction, a microseismic monitoring system was deployed for the study of the stress evolution process, damage degree and distribution characteristics in the tailgate and headgate. The mine under study is the 62113 outburst working face of Xin Zhuangzi coalmine in Huainan mining area. The whole process of microfractures initiation,extension, interaction and coalescence mechanisms during the progressive failure processes of the coal rock within the delineated and typical event clusters were investigated by means of a two dimensional realistic failure process analysis code(RFPA2D-Flow). The results show that the microseismic events gradually create different-sized event clusters. The microseismicity of the tailgate is significantly higher than that of the headgate. The study indicates that the greater anomalous stress region matches the area where microfractures continuously develop and finally connect to each other and form a fissure zone.Due to the mine layout and stress concentration, the ruptured area is mainly located on the left shoulder of the tailgate roof. The potential anomalous stress region of the coal roadway obtained by numerical simulation is relatively in good agreement with the trend of spatial macro evolution of coal rock microfractures captured by the microseismic monitoring system. The research results can provide important basis for understanding instability failure mechanism of deep roadway and microseismic activity law in complex geologic conditions, and it ultimately can be used to guide the selection and optimization of reinforcement and protection scheme.展开更多
This paper focuses on the numerical stability of the block θ methods adapted to differential equations with a delay argument. For the block θ methods, an interpolation procedure is introduced which leads to the nume...This paper focuses on the numerical stability of the block θ methods adapted to differential equations with a delay argument. For the block θ methods, an interpolation procedure is introduced which leads to the numerical processes that satisfy an important asymptotic stability condition related to the class of test problems y′(t)=ay(t)+by(t-τ) with a,b∈C, Re(a)<-|b| and τ>0. We prove that the block θ method is GP stable if and only if the method is A stable for ordinary differential equations. Furthermore, it is proved that the P and GP stability are equivalent for the block θ method.展开更多
This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solutio...This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1.展开更多
基金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.
文摘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.
文摘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.
基金Supported by the National Natural Science Foundation of China under Grant No.51975138the High-Tech Ship Scientific Research Project from the Ministry of Industry and Information Technology under Grant No.CJ05N20the National Defense Basic Research Project under Grant No.JCKY2023604C006.
文摘Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.
基金The National Natural Science Foundation of China(Grant No.52201376)the Natural Science Foundation of Hubei Province,China(Grant No.2023AFB683).
文摘In this investigation,a hybrid approach integrating the IDDES turbulence model and FW-H is employed to forecast the hydroacoustic of the rim driven thruster(RDT)under non-cavitation and uniform flow conditions at varying loading conditions(J=0.3 and J=0.6).It is revealed that the quadrupole term contribution in the P-FWH method significantly affects the monopole term in the low-frequency region,while it mainly affects the dipole term in the high-frequency region.Specifically,the overall sound pressure levels(SPL)of the RDT using the P-FWH method are 2.27 dB,10.03 dB,and 16.73 dB at the receiving points from R1 to R3 under the heavy-loaded condition,while they increase by 0.67 dB at R1,and decrease by 14.93 dB at R2,and 22.20 dB at R3,for the light-loaded condition.The study also utilizes the pressure-time derivatives to visualize the numerical noise and to pinpoint the dynamics of the vortex cores,and the optimization of the grid design can significantly reduce the numerical noise.The computational accuracy of the P-FWH method can meet the noise requirements for the preliminary design of rim driven thrusters.
基金The National Natural Science Foundation of China(No.51378121)the Fok Ying Tung Education Foundation(No.141076)the Scientific Innovation Research of College Graduates in Jiangsu Province(No.KYLX_0164)
文摘In order to predict the long-term rutting of asphalt pavement, the effective temperature for pavement rutting is calculated using the numerical simulation method. The transient temperature field of asphalt pavement was simulated based on actual meteorological data of Nanjing. 24-hour rutting development under a transient temperature field was calculated in each month. The rutting depth accumulated under the static temperature field was also estimated and the relationship between constant temperature parameters was analyzed. Then the effective temperature for pavement rutting was determined based on the rutting equivalence principle. The results show that the monthly effective temperature is above 40 t in July and August, while in June and September it ranges from 30 to 40 Rutting development can be ignored when the monthly effective temperature is less than 30 t. The yearly effective temperature for rutting in Nanjing is around 38. 5 t. The long-term rutting prediction model based on the effective temperature can reflect the influences of meteorological factors and traffic time distribution.
基金Project(41630642)supported by the Key Project of National Natural Science Foundation of ChinaProject(51974360)supported by the National Natural Science Foundation of ChinaProject(2018JJ3656)supported by the Natural Science Foundation of Hunan Province,China。
文摘In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.
基金supported by the Office of Naval Research (Grant No.N000141010778) to the University of Oklahomathe National Natural Sciences Foundation of China (Grant Nos. 40930950,41075043,and 4092116037) to the Institute of Atmospheric Physicsprovided by NOAA/Office of Oceanic and Atmospheric Research under NOAA-University of Oklahoma Cooperative Agreement No. (NA17RJ1227),U.S. Department of Commerce
文摘Built on the integral formulas in Part I,numerical methods are developed for computing velocity potential and streamfunction in a limited domain.When there is no inner boundary（around a data hole） inside the domain,the total solution is the sum of the internally and externally induced parts.For the internally induced part,three numerical schemes（grid-staggering,local-nesting and piecewise continuous integration） are designed to deal with the singularity of the Green＇s function encountered in numerical calculations.For the externally induced part,by setting the velocity potential（or streamfunction） component to zero,the other component of the solution can be computed in two ways：（1） Solve for the density function from its boundary integral equation and then construct the solution from the boundary integral of the density function.（2） Use the Cauchy integral to construct the solution directly.The boundary integral can be discretized on a uniform grid along the boundary.By using local-nesting（or piecewise continuous integration）,the scheme is refined to enhance the discretization accuracy of the boundary integral around each corner point（or along the entire boundary）.When the domain is not free of data holes,the total solution contains a data-hole-induced part,and the Cauchy integral method is extended to construct the externally induced solution with irregular external and internal boundaries.An automated algorithm is designed to facilitate the integrations along the irregular external and internal boundaries.Numerical experiments are performed to evaluate the accuracy and efficiency of each scheme relative to others.
文摘Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.
文摘One of the critical aspects in mine design is slope stability analysis and the determination of stable slopes. In the Chador- Malu iron ore mine, one of the most important iron ore mines in central Iran, it was considered vital to perform a comprehensive slope stability analysis. At first, we divided the existing rock hosting pit into six zones and a geotechnical map was prepared. Then, the value of MRMR （Mining Rock Mass Rating） was determined for each zone. Owing to the fact that the Chador-Malu iron ore mine is located in a highly tectonic area and the rock mass completely crushed, the Hoek-Brown failure criterion was found suitable to estimate geo-mechanical parameters. After that, the value of cohesion （c） and friction angle （tp） were calculated for different geotechnical zones and relative graphs and equations were derived as a function of slope height. The stability analyses using numerical and limit equilibrium methods showed that some instability problems might occur by increasing the slope height. Therefore, stable slopes for each geotechnical zone and prepared sections were calculated and presented as a function of slope height.
文摘Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concrete pillars (CP) were carried out. During the construction of the second CP, in the space between the remaining pillars, one of the pillars failed leading to the progressive failure of other pillars until 4 000 m 2 of mine had collapsed within a few minutes. In this work, this phenomenon is described by applying both numerical and empirical methods and the respective results are compared. The results of numerical modelling are found to be closer to the actual condition than those of the empirical method. Also, a width-to-height (W/H) ratio less than 1, an inadequate support system and the absence of a detailed program for pillar recovery are shown to be the most important causes of the Domino failure in this mine.
文摘With the development of modern science and technology, especially computer science, the numerical simulation method has been widely used in material hot-working. Mary achievements have been made in this field by using the numerical simulation method. The numerical simulation method, especially finite element method fully described in this paper.Applications of the numerical simulation method in material hot-working are also discussed. Finally, the future of the numerical simulation method is outlined.
基金Project(xjj20100078) supported by the Fundamental Research Funds for the Central Universities in China
文摘Several parameter identification methods of thermal response test were evaluated through numerical and experimental study.A three-dimensional finite-volume numerical model was established under the assumption that the soil thermal conductivity had been known in the simulation of thermal response test.The thermal response curve was firstly obtained through numerical calculation.Then,the accuracy of the numerical model was verified with measured data obtained through a thermal response test.Based on the numerical and experimental thermal response curves,the thermal conductivity of the soil was calculated by different parameter identification methods.The calculated results were compared with the assumed value and then the accuracy of these methods was evaluated.Furthermore,the effects of test time,variable data quality,borehole radius,initial ground temperature,and heat injection rate were analyzed.The results show that the method based on cylinder-source model has a low precision and the identified thermal conductivity decreases with an increase in borehole radius.For parameter estimation,the measuring accuracy of the initial temperature of the deep ground soil has greater effect on identified thermal conductivity.
基金Project supported by the Special Scientific Research Project for Public Interest(Grant No.GYHY201206009)the Fundamental Research Funds for the Central Universities,China(Grant Nos.lzujbky-2012-13 and lzujbky-2013-11)the National Basic Research Program of China(Grant Nos.2012CB955902 and 2013CB430204)
文摘Model error is one of the key factors restricting the accuracy of numerical weather prediction （NWP）. Considering the continuous evolution of the atmosphere, the observed data （ignoring the measurement error） can be viewed as a series of solutions of an accurate model governing the actual atmosphere. Model error is represented as an unknown term in the accurate model, thus NWP can be considered as an inverse problem to uncover the unknown error term. The inverse problem models can absorb long periods of observed data to generate model error correction procedures. They thus resolve the deficiency and faultiness of the NWP schemes employing only the initial-time data. In this study we construct two inverse problem models to estimate and extrapolate the time-varying and spatial-varying model errors in both the historical and forecast periods by using recent observations and analogue phenomena of the atmosphere. Numerical experiment on Burgers＇ equation has illustrated the substantial forecast improvement using inverse problem algorithms. The proposed inverse problem methods of suppressing NWP errors will be useful in future high accuracy applications of NWP.
基金Financial support for this work, provided by the National Natural Science Foundation of China (Nos.51674189,51304154,and 51327007)the Youth Science and technology new star of Shaanxi Province (No.2016KJXX-37)the Scientific research plan of Shaanxi Education Department (No.16JK1487)
文摘It is believed that the microseismicity induced by mining effect and gas gradient disturbance stress is a precursor to the essential characteristics of roadway unstability. In order to effectively identify and evaluate the stability of coal roadways in the process of mine development and extraction, a microseismic monitoring system was deployed for the study of the stress evolution process, damage degree and distribution characteristics in the tailgate and headgate. The mine under study is the 62113 outburst working face of Xin Zhuangzi coalmine in Huainan mining area. The whole process of microfractures initiation,extension, interaction and coalescence mechanisms during the progressive failure processes of the coal rock within the delineated and typical event clusters were investigated by means of a two dimensional realistic failure process analysis code(RFPA2D-Flow). The results show that the microseismic events gradually create different-sized event clusters. The microseismicity of the tailgate is significantly higher than that of the headgate. The study indicates that the greater anomalous stress region matches the area where microfractures continuously develop and finally connect to each other and form a fissure zone.Due to the mine layout and stress concentration, the ruptured area is mainly located on the left shoulder of the tailgate roof. The potential anomalous stress region of the coal roadway obtained by numerical simulation is relatively in good agreement with the trend of spatial macro evolution of coal rock microfractures captured by the microseismic monitoring system. The research results can provide important basis for understanding instability failure mechanism of deep roadway and microseismic activity law in complex geologic conditions, and it ultimately can be used to guide the selection and optimization of reinforcement and protection scheme.
文摘This paper focuses on the numerical stability of the block θ methods adapted to differential equations with a delay argument. For the block θ methods, an interpolation procedure is introduced which leads to the numerical processes that satisfy an important asymptotic stability condition related to the class of test problems y′(t)=ay(t)+by(t-τ) with a,b∈C, Re(a)<-|b| and τ>0. We prove that the block θ method is GP stable if and only if the method is A stable for ordinary differential equations. Furthermore, it is proved that the P and GP stability are equivalent for the block θ method.
文摘This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1.
