Two kinds of bank erosion mechanisms were analyzed, including fluvial and non-fluvial controlled mechanisms, and mechanical methods of simulating the erosion processes of cohesive, non-cohesive and composite riverbank...Two kinds of bank erosion mechanisms were analyzed, including fluvial and non-fluvial controlled mechanisms, and mechanical methods of simulating the erosion processes of cohesive, non-cohesive and composite riverbanks were improved. Then a two-dimensional numerical model of the channel deformation was developed, consisting of a 2D flow and sediment transport submodel and bank-erosion submodels of different soil riverbanks. In the model, a new technique for updating the bank geometry during the bed evolution was presented, which combines closely two kinds of submodels. The proposed model is capable of not only predicting the processes of flood routing and longitudinal channel deformation in natural rivers, but also simulating the processes of lateral channel deformation, especially the processes of lateral erosion and failure of cohesive, non-cohesive and composite riverbanks.展开更多
Machine learning-assisted methods for rapid and accurate prediction of temperature field,mushy zone,and grain size were proposed for the heating−cooling combined mold(HCCM)horizontal continuous casting of C70250 alloy...Machine learning-assisted methods for rapid and accurate prediction of temperature field,mushy zone,and grain size were proposed for the heating−cooling combined mold(HCCM)horizontal continuous casting of C70250 alloy plates.First,finite element simulations of casting processes were carried out with various parameters to build a dataset.Subsequently,different machine learning algorithms were employed to achieve high precision in predicting temperature fields,mushy zone locations,mushy zone inclination angle,and billet grain size.Finally,the process parameters were quickly optimized using a strategy consisting of random generation,prediction,and screening,allowing the mushy zone to be controlled to the desired target.The optimized parameters are 1234℃for heating mold temperature,47 mm/min for casting speed,and 10 L/min for cooling water flow rate.The optimized mushy zone is located in the middle of the second heat insulation section and has an inclination angle of roughly 7°.展开更多
Rock is geometrically and mechanically multiscale in nature,and the traditional phenomenological laws at the macroscale cannot render a quantitative relationship between microscopic damage of rocks and overall rock st...Rock is geometrically and mechanically multiscale in nature,and the traditional phenomenological laws at the macroscale cannot render a quantitative relationship between microscopic damage of rocks and overall rock structural degradation.This may lead to problems in the evaluation of rock structure stability and safe life.Multiscale numerical modeling is regarded as an effective way to gain insight into factors affecting rock properties from a cross-scale view.This study compiles the history of theoretical developments and numerical techniques related to rock multiscale issues according to different modeling architectures,that is,the homogenization theory,the hierarchical approach,and the concurrent approach.For these approaches,their benefits,drawbacks,and application scope are underlined.Despite the considerable attempts that have been made,some key issues still result in multiple challenges.Therefore,this study points out the perspectives of rock multiscale issues so as to provide a research direction for the future.The review results show that,in addition to numerical techniques,for example,high-performance computing,more attention should be paid to the development of an advanced constitutive model with consideration of fine geometrical descriptions of rock to facilitate solutions to multiscale problems in rock mechanics and rock engineering.展开更多
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.展开更多
A three-dimensional numerical model of laser-arc hybrid plasma for aluminum alloy fillet joints is developed in this study.This mod-el accounts for the geometric complexity of fillet joints,the physical properties of ...A three-dimensional numerical model of laser-arc hybrid plasma for aluminum alloy fillet joints is developed in this study.This mod-el accounts for the geometric complexity of fillet joints,the physical properties of shielding gases with varying He-Ar ratios,and the coupling between arc plasma and laser-induced metal plume.The accuracy of the model is validated using a high-speed camera.The effects of varying He contents in the shielding gas on both the temperature and flow velocity of hybrid plasma,as well as the distribu-tion of laser-induced metal vapor mass,were investigated separately.The maximum temperature and size of arc plasma decrease as the He volume ratio increases,the arc distribution becomes more concentrated,and its flow velocity initially decreases and then sharply increases.At high helium content,both the flow velocity of hybrid plasma and metal vapor are high,the metal vapor is con-centrated on the right side of keyhole,and its flow appears chaotic.The flow state of arc plasma is most stable when the shielding gas consists of 50%He+50%Ar.展开更多
Erosion is an important issue in soil science and is related to many environmental problems,such as soil erosion and sediment transport.Establishing a simulation model suitable for soil erosion prediction is of great ...Erosion is an important issue in soil science and is related to many environmental problems,such as soil erosion and sediment transport.Establishing a simulation model suitable for soil erosion prediction is of great significance not only to accurately predict the process of soil separation by runoff,but also improve the physical model of soil erosion.In this study,we develop a graphic processing unit(GPU)-based numerical model that combines two-dimensional(2D)hydrodynamic and Green-Ampt(G-A)infiltration modelling to simulate soil erosion.A Godunov-type scheme on a uniform and structured square grid is then generated to solve the relevant shallow water equations(SWEs).The highlight of this study is the use of GPU-based acceleration technology to enable numerical models to simulate slope and watershed erosion in an efficient and high-resolution manner.The results show that the hydrodynamic model performs well in simulating soil erosion process.Soil erosion is studied by conducting calculation verification at the slope and basin scales.The first case involves simulating soil erosion process of a slope surface under indoor artificial rainfall conditions from 0 to 1000 s,and there is a good agreement between the simulated values and the measured values for the runoff velocity.The second case is a river basin experiment(Coquet River Basin)that involves watershed erosion.Simulations of the erosion depth change and erosion cumulative amount of the basin during a period of 1-40 h show an elevation difference of erosion at 0.5-3.0 m,especially during the period of 20-30 h.Nine cross sections in the basin are selected for simulation and the results reveal that the depth of erosion change value ranges from-0.86 to-2.79 m and the depth of deposition change value varies from 0.38 to 1.02 m.The findings indicate that the developed GPU-based hydrogeomorphological model can reproduce soil erosion processes.These results are valuable for rainfall runoff and soil erosion predictions on rilled hillslopes and river basins.展开更多
In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-d...In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-dimensional nonlinear time fractional thermal diffusion model.The time Caputo fractional derivative is approximated by using the L2-1formula,the first-order derivative and nonlinear term are discretized by some second-order approximation formulas,and the quadratic finite element is used to approximate the spatial direction.The error accuracy O(h3+t2)is obtained,which is verified by the numerical results.展开更多
The deep lateral growth process of the plateau is a hot topic in earth science.As the growth boundary zone of the northeastern Tibetan Plateau,southern Ningxia geomorphic zone is controlled by four major arcuate fault...The deep lateral growth process of the plateau is a hot topic in earth science.As the growth boundary zone of the northeastern Tibetan Plateau,southern Ningxia geomorphic zone is controlled by four major arcuate fault zones(Fig.1),namely Haiyuan fault(HYF),XiangshanTianjingshan fault(XTSF),Yantongshan fault(YTSF).展开更多
To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional (2D) time- domain transient electromagnetic secondary field of t...To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional (2D) time- domain transient electromagnetic secondary field of the line source based on the DuFort- Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force (Emf), and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the fiat layered model with high precision--the maximum relative error is less than 0.01% between our method and the analytical method--and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.展开更多
This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in r...This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in real seismic data application. Our results reveal that there is an offset limitation for both PP- and PS-waves in elliptical anisotropy fitting, and that PS-waves show a wider applicable offset range and larger observable azimuthal anisotropy than PP-waves. The major axis of the elliptical fit to the amplitudes of the R-component is perpendicular to the fracture strike, which is opposite to that in PP-wave analysis. The azimuthal interval travel time of PS-waves shows a nearly elliptical distribution and the major axis of the fit ellipse is perpendicular to the fracture strike, which is same as that in PP-wave analysis. For data within the applicable offset range, the anisotropic magnitude obtained from amplitude and travel time attributes of PP- and PS-waves exhibits a dependence on fracture density, and the major to minor axis ratio of the fit ellipse may be used to infer the relative distribution of fracture densities.展开更多
In order to reasonably simulate tidal currents around small structures such as piles in a large-scale model domain, a 2-D hydrodynamic integrated model for Bohai Sea is established with the finite element method. The ...In order to reasonably simulate tidal currents around small structures such as piles in a large-scale model domain, a 2-D hydrodynamic integrated model for Bohai Sea is established with the finite element method. The grid can be discretionarily refined as a non-structure triangle or quadrilateral so that piers can be treated as one or several impermeable elements with an area of 20 to 30 km^2 in a model domain over 85 700 km^2. The computational results of tidal levels and horizontal velocities are in good agreement with the field data. Based on the computed results by the model, the layout of an open 105 DWT liquefied natural gas (LNG)terminal in Caofeidian, Bohal Sea is effectively and reasonably optimized. It can be concluded that the model is suitable and reasonable for direct simulation of tidal currents around small structures in projects.展开更多
To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of ...To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.展开更多
Haihua Islands is a large artificial island in Danzhou, Hainan. The construction of Haihua Islands changes the hYdrodynamic environment of Yangpu waters, and further affects its morphological change. Delft3D is used t...Haihua Islands is a large artificial island in Danzhou, Hainan. The construction of Haihua Islands changes the hYdrodynamic environment of Yangpu waters, and further affects its morphological change. Delft3D is used to set up a two dimensional nested hydrodynamic and sediment model for Yangpu waters in this paper, and this paper focuses on simulating the velocity and morphological change due to the construction of Haihua Islands after the verification of the model. The seabed deposition is small because of low suspended sediment concentration and less sand source near Yangpu waters. The bed level erodes in the south area of Xiaochan Reef and the Yangpu channel due to the velocity increase in the area.展开更多
Automation and intelligence have become the primary trends in the design of investment casting processes.However,the design of gating and riser systems still lacks precise quantitative evaluation criteria.Numerical si...Automation and intelligence have become the primary trends in the design of investment casting processes.However,the design of gating and riser systems still lacks precise quantitative evaluation criteria.Numerical simulation plays a significant role in quantitatively evaluating current processes and making targeted improvements,but its limitations lie in the inability to dynamically reflect the formation outcomes of castings under varying process conditions,making real-time adjustments to gating and riser designs challenging.In this study,an automated design model for gating and riser systems based on integrated parametric 3D modeling-simulation framework is proposed,which enhances the flexibility and usability of evaluating the casting process by simulation.Firstly,geometric feature extraction technology is employed to obtain the geometric information of the target casting.Based on this information,an automated design framework for gating and riser systems is established,incorporating multiple structural parameters for real-time process control.Subsequently,the simulation results for various structural parameters are analyzed,and the influence of these parameters on casting formation is thoroughly investigated.Finally,the optimal design scheme is generated and validated through experimental verification.Simulation analysis and experimental results show that using a larger gate neck(24 mm in side length) and external risers promotes a more uniform temperature distribution and a more stable flow state,effectively eliminating shrinkage cavities and enhancing process yield by 15%.展开更多
This paper investigates ruin,capital injection,and dividends for a two-dimensional risk model.The model posits that surplus levels of insurance companies are governed by a perturbed composite Poisson risk model.This m...This paper investigates ruin,capital injection,and dividends for a two-dimensional risk model.The model posits that surplus levels of insurance companies are governed by a perturbed composite Poisson risk model.This model introduces a dependence between the two surplus levels,present in both the associated perturbations and the claims resulting from common shocks.Critical levels of capital injection and dividends are established for each of the two risks.The surplus levels are observed discretely at fixed intervals,guiding decisions on capital injection,dividends,and ruin at these junctures.This study employs a two-dimensional Fourier cosine series expansion method to approximate the finite time expected discounted operating cost until ruin.The ensuing approximation error is also quantified.The validity and accuracy of the method are corroborated through numerical examples.Furthermore,the research delves into the optimal capital allocation problem.展开更多
This review summarizes the scientific and technical progress in atmospheric modeling in China since 2011,including the dynamical core,model physics,data assimilation,ensemble forecasting,and model evaluation strategie...This review summarizes the scientific and technical progress in atmospheric modeling in China since 2011,including the dynamical core,model physics,data assimilation,ensemble forecasting,and model evaluation strategies.In terms of the dynamical core,important efforts have been made in the improvement of the existing model formulations and in exploring new modeling approaches that can better adapt to massively parallel computers and global multiscale modeling.With regard to model physics,various achievements in physical representations have been made,especially a trend toward scale-aware parameterization for accommodating the increase of model resolution.In the field of data assimilation,a 4D-Var system has been developed and is operationally used by the National Meteorological Center of China,and its performance is promising.Furthermore,ensemble forecasting has played a more important role in operational forecast systems and progressed in many fundamental techniques.Model evaluation strategies,including key performance metrics and standardized experimental protocols,have been proposed and widely applied to better understand the strengths and weaknesses of the systems,offering key routes for model improvement.The paper concludes with a concise summary of the status quo and a brief outlook in terms of future development.展开更多
In this paper, a heavy sea fog episode that occurred over the Yellow Sea on 9 March 2005 is investigated. The sea fog patch, with a spatial scale of several hundred kilometers at its mature stage, reduced visibility a...In this paper, a heavy sea fog episode that occurred over the Yellow Sea on 9 March 2005 is investigated. The sea fog patch, with a spatial scale of several hundred kilometers at its mature stage, reduced visibility along the Shandong Peninsula coast to 100 m or much less at some sites. Satellite images, surface observations and soundings at islands and coasts, and analyses from the Japan Meteorology Agency (JMA) axe used to describe and analyze this event. The analysis indicates that this sea fog can be categorized as advection cooling fog. The main features of this sea fog including fog area and its movement axe reasonably reproduced by the Fifth-generation Pennsylvania State University/National Center for Atmospheric Research Mesoscale Model (MM5). Model results suggest that the formation and evolution of this event can be outlined as: (1) southerly warm/moist advection of low-level air resulted in a strong sea-surface-based inversion with a thickness of about 600 m; (2) when the inversion moved from the warmer East Sea to the colder Yellow Sea, a thermal internal boundary layer (TIBL) gradually formed at the base of the inversion while the sea fog grew in response to cooling and moistening by turbulence mixing; (3) the sea fog developed as the TIBL moved northward and (4) strong northerly cold and dry wind destroyed the TIBL and dissipated the sea fog. The principal findings of this study axe that sea fog forms in response to relatively persistent southerly waxm/moist wind and a cold sea surface, and that turbulence mixing by wind shear is the primary mechanism for the cooling and moistening the marine layer. In addition, the study of sensitivity experiments indicates that deterministic numerical modeling offers a promising approach to the prediction of sea fog over the Yellow Sea but it may be more efficient to consider ensemble numerical modeling because of the extreme sensitivity to model input.展开更多
Recent advances in three-dimensional numerical simulations of mantle convection have aided in approximately reproducing continental movement since the Pangea breakup at 200 Ma. These have also led to a better understa...Recent advances in three-dimensional numerical simulations of mantle convection have aided in approximately reproducing continental movement since the Pangea breakup at 200 Ma. These have also led to a better understanding of the thermal and mechanical coupling between mantle convection and surface plate motion and predictions of the configuration of the next supercontinent. The simulations of mantle convection from 200 Ma to the present reveals that the development of large-scale cold mantle downwellings in the North Tethys Ocean at the earlier stage of the Pangea breakup triggered the northward movement of the Indian subcontinent. The model of high temperature anomaly region beneath Pangea resulting from the thermal insulation effect support the breakup of Pangea in the real Earth time scale, as also suggested in previous geological and geodynamic models. However, considering the low radioactive heat generation rate of the depleted upper mantle, the high temperature anomaly region might have been generated by upwelling plumes with contribution of deep subducted TTG(tonalite-trondhjemite-granite) materials enriched in radiogenic elements. Integrating the numerical results of mantle convection from 200 Ma to the present, and from the present to the future, it is considered that the mantle drag force acting on the base of continents may be comparable to the slab pull force, which implies that convection in the shallower part of the mantle is strongly coupled with surface plate motion.展开更多
Stress-induced failure is among the most common causes of instability in Canadian deep underground mines.Open stoping is the most widely practiced underground excavation method in these mines,and creates large stopes ...Stress-induced failure is among the most common causes of instability in Canadian deep underground mines.Open stoping is the most widely practiced underground excavation method in these mines,and creates large stopes which are subjected to stress-induced failure.The probability of failure(POF)depends on many factors,of which the geometry of an open stope is especially important.In this study,a methodology is proposed to assess the effect of stope geometrical parameters on the POF,using numerical modelling.Different ranges for each input parameter are defined according to previous surveys on open stope geometry in a number of Canadian underground mines.A Monte-Carlo simulation technique is combined with the finite difference code FLAC3D,to generate model realizations containing stopes with different geometrical features.The probability of failure(POF)for different categories of stope geometry,is calculated by considering two modes of failure;relaxation-related gravity driven(tensile)failure and rock mass brittle failure.The individual and interactive effects of stope geometrical parameters on the POF,are analyzed using a general multi-level factorial design.Finally,mathematical optimization techniques are employed to estimate the most stable stope conditions,by determining the optimal ranges for each stope’s geometrical parameter.展开更多
Vehicle electrification,an important method for reducing carbon emissions from road transport,has been promoted globally.In this study,we analyze how individuals adapt to this transition in transportation and its subs...Vehicle electrification,an important method for reducing carbon emissions from road transport,has been promoted globally.In this study,we analyze how individuals adapt to this transition in transportation and its subsequent impact on urban structure.Considering the varying travel costs associated with electric and fuel vehicles,we analyze the dynamic choices of households concerning house locations and vehicle types in a two-dimensional monocentric city.A spatial equilibrium is developed to model the interactions between urban density,vehicle age and vehicle type.An agent-based microeconomic residential choice model dynamically coupled with a house rent market is developed to analyze household choices of home locations and vehicle energy types,considering vehicle ages and competition for public charging piles.Key findings from our proposed models show that the proportion of electric vehicles(EVs)peaks at over 50%by the end of the first scrappage period,accompanied by more than a 40%increase in commuting distance and time compared to the scenario with only fuel vehicles.Simulation experiments on a theoretical grid indicate that heterogeneity-induced residential segregation can lead to urban sprawl and congestion.Furthermore,households with EVs tend to be located farther from the city center,and an increase in EV ownership contributes to urban expansion.Our study provides insights into how individuals adapt to EV transitions and the resulting impacts on home locations and land use changes.It offers a novel perspective on the dynamic interactions between EV adoption and urban development.展开更多
基金the National Key Research Program of Basic Science Foundation (Grant No. G19990436) partially by the Science Fund for Creative Research Groups of the National Natural Science Foundation of China (Grant No. 50221903).
文摘Two kinds of bank erosion mechanisms were analyzed, including fluvial and non-fluvial controlled mechanisms, and mechanical methods of simulating the erosion processes of cohesive, non-cohesive and composite riverbanks were improved. Then a two-dimensional numerical model of the channel deformation was developed, consisting of a 2D flow and sediment transport submodel and bank-erosion submodels of different soil riverbanks. In the model, a new technique for updating the bank geometry during the bed evolution was presented, which combines closely two kinds of submodels. The proposed model is capable of not only predicting the processes of flood routing and longitudinal channel deformation in natural rivers, but also simulating the processes of lateral channel deformation, especially the processes of lateral erosion and failure of cohesive, non-cohesive and composite riverbanks.
基金financially supported by the National Key Research and Development Program of China (No. 2023YFB3812601)the National Natural Science Foundation of China (No. 51925401)the Young Elite Scientists Sponsorship Program by CAST, China (No. 2022QNRC001)。
文摘Machine learning-assisted methods for rapid and accurate prediction of temperature field,mushy zone,and grain size were proposed for the heating−cooling combined mold(HCCM)horizontal continuous casting of C70250 alloy plates.First,finite element simulations of casting processes were carried out with various parameters to build a dataset.Subsequently,different machine learning algorithms were employed to achieve high precision in predicting temperature fields,mushy zone locations,mushy zone inclination angle,and billet grain size.Finally,the process parameters were quickly optimized using a strategy consisting of random generation,prediction,and screening,allowing the mushy zone to be controlled to the desired target.The optimized parameters are 1234℃for heating mold temperature,47 mm/min for casting speed,and 10 L/min for cooling water flow rate.The optimized mushy zone is located in the middle of the second heat insulation section and has an inclination angle of roughly 7°.
基金National Natural Science Foundation of China,Grant/Award Numbers:52192691,52192690。
文摘Rock is geometrically and mechanically multiscale in nature,and the traditional phenomenological laws at the macroscale cannot render a quantitative relationship between microscopic damage of rocks and overall rock structural degradation.This may lead to problems in the evaluation of rock structure stability and safe life.Multiscale numerical modeling is regarded as an effective way to gain insight into factors affecting rock properties from a cross-scale view.This study compiles the history of theoretical developments and numerical techniques related to rock multiscale issues according to different modeling architectures,that is,the homogenization theory,the hierarchical approach,and the concurrent approach.For these approaches,their benefits,drawbacks,and application scope are underlined.Despite the considerable attempts that have been made,some key issues still result in multiple challenges.Therefore,this study points out the perspectives of rock multiscale issues so as to provide a research direction for the future.The review results show that,in addition to numerical techniques,for example,high-performance computing,more attention should be paid to the development of an advanced constitutive model with consideration of fine geometrical descriptions of rock to facilitate solutions to multiscale problems in rock mechanics and rock engineering.
基金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.
基金supported by the National Natural Science Foundation of China(Grant No.52375340,51975263,52405366).
文摘A three-dimensional numerical model of laser-arc hybrid plasma for aluminum alloy fillet joints is developed in this study.This mod-el accounts for the geometric complexity of fillet joints,the physical properties of shielding gases with varying He-Ar ratios,and the coupling between arc plasma and laser-induced metal plume.The accuracy of the model is validated using a high-speed camera.The effects of varying He contents in the shielding gas on both the temperature and flow velocity of hybrid plasma,as well as the distribu-tion of laser-induced metal vapor mass,were investigated separately.The maximum temperature and size of arc plasma decrease as the He volume ratio increases,the arc distribution becomes more concentrated,and its flow velocity initially decreases and then sharply increases.At high helium content,both the flow velocity of hybrid plasma and metal vapor are high,the metal vapor is con-centrated on the right side of keyhole,and its flow appears chaotic.The flow state of arc plasma is most stable when the shielding gas consists of 50%He+50%Ar.
基金This research was funded by the National Natural Science Foundation of China(52079106,52009104,51609199)the National Key Research and Development Program of China(2016YFC0402704).
文摘Erosion is an important issue in soil science and is related to many environmental problems,such as soil erosion and sediment transport.Establishing a simulation model suitable for soil erosion prediction is of great significance not only to accurately predict the process of soil separation by runoff,but also improve the physical model of soil erosion.In this study,we develop a graphic processing unit(GPU)-based numerical model that combines two-dimensional(2D)hydrodynamic and Green-Ampt(G-A)infiltration modelling to simulate soil erosion.A Godunov-type scheme on a uniform and structured square grid is then generated to solve the relevant shallow water equations(SWEs).The highlight of this study is the use of GPU-based acceleration technology to enable numerical models to simulate slope and watershed erosion in an efficient and high-resolution manner.The results show that the hydrodynamic model performs well in simulating soil erosion process.Soil erosion is studied by conducting calculation verification at the slope and basin scales.The first case involves simulating soil erosion process of a slope surface under indoor artificial rainfall conditions from 0 to 1000 s,and there is a good agreement between the simulated values and the measured values for the runoff velocity.The second case is a river basin experiment(Coquet River Basin)that involves watershed erosion.Simulations of the erosion depth change and erosion cumulative amount of the basin during a period of 1-40 h show an elevation difference of erosion at 0.5-3.0 m,especially during the period of 20-30 h.Nine cross sections in the basin are selected for simulation and the results reveal that the depth of erosion change value ranges from-0.86 to-2.79 m and the depth of deposition change value varies from 0.38 to 1.02 m.The findings indicate that the developed GPU-based hydrogeomorphological model can reproduce soil erosion processes.These results are valuable for rainfall runoff and soil erosion predictions on rilled hillslopes and river basins.
基金the National Natural Science Fund(11661058,11761053)Natural Science Fund of Inner Mongolia Autonomous Region(2016MS0102,2017MS0107)+1 种基金Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region(NJYT-17-A07)National Undergraduate Innovative Training Project of Inner Mongolia University(201710126026).
文摘In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-dimensional nonlinear time fractional thermal diffusion model.The time Caputo fractional derivative is approximated by using the L2-1formula,the first-order derivative and nonlinear term are discretized by some second-order approximation formulas,and the quadratic finite element is used to approximate the spatial direction.The error accuracy O(h3+t2)is obtained,which is verified by the numerical results.
基金supported by the Geological Survey Project of China Geological Survey(Grant No.DD20190018)
文摘The deep lateral growth process of the plateau is a hot topic in earth science.As the growth boundary zone of the northeastern Tibetan Plateau,southern Ningxia geomorphic zone is controlled by four major arcuate fault zones(Fig.1),namely Haiyuan fault(HYF),XiangshanTianjingshan fault(XTSF),Yantongshan fault(YTSF).
基金supported by the National High Technology Research and Development Program (863 Program)(2009AA06Z108)
文摘To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional (2D) time- domain transient electromagnetic secondary field of the line source based on the DuFort- Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force (Emf), and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the fiat layered model with high precision--the maximum relative error is less than 0.01% between our method and the analytical method--and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.
文摘This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in real seismic data application. Our results reveal that there is an offset limitation for both PP- and PS-waves in elliptical anisotropy fitting, and that PS-waves show a wider applicable offset range and larger observable azimuthal anisotropy than PP-waves. The major axis of the elliptical fit to the amplitudes of the R-component is perpendicular to the fracture strike, which is opposite to that in PP-wave analysis. The azimuthal interval travel time of PS-waves shows a nearly elliptical distribution and the major axis of the fit ellipse is perpendicular to the fracture strike, which is same as that in PP-wave analysis. For data within the applicable offset range, the anisotropic magnitude obtained from amplitude and travel time attributes of PP- and PS-waves exhibits a dependence on fracture density, and the major to minor axis ratio of the fit ellipse may be used to infer the relative distribution of fracture densities.
文摘In order to reasonably simulate tidal currents around small structures such as piles in a large-scale model domain, a 2-D hydrodynamic integrated model for Bohai Sea is established with the finite element method. The grid can be discretionarily refined as a non-structure triangle or quadrilateral so that piers can be treated as one or several impermeable elements with an area of 20 to 30 km^2 in a model domain over 85 700 km^2. The computational results of tidal levels and horizontal velocities are in good agreement with the field data. Based on the computed results by the model, the layout of an open 105 DWT liquefied natural gas (LNG)terminal in Caofeidian, Bohal Sea is effectively and reasonably optimized. It can be concluded that the model is suitable and reasonable for direct simulation of tidal currents around small structures in projects.
文摘To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.
文摘Haihua Islands is a large artificial island in Danzhou, Hainan. The construction of Haihua Islands changes the hYdrodynamic environment of Yangpu waters, and further affects its morphological change. Delft3D is used to set up a two dimensional nested hydrodynamic and sediment model for Yangpu waters in this paper, and this paper focuses on simulating the velocity and morphological change due to the construction of Haihua Islands after the verification of the model. The seabed deposition is small because of low suspended sediment concentration and less sand source near Yangpu waters. The bed level erodes in the south area of Xiaochan Reef and the Yangpu channel due to the velocity increase in the area.
基金financially supported by the National Key Research and Development Program of China (2022YFB3706802)。
文摘Automation and intelligence have become the primary trends in the design of investment casting processes.However,the design of gating and riser systems still lacks precise quantitative evaluation criteria.Numerical simulation plays a significant role in quantitatively evaluating current processes and making targeted improvements,but its limitations lie in the inability to dynamically reflect the formation outcomes of castings under varying process conditions,making real-time adjustments to gating and riser designs challenging.In this study,an automated design model for gating and riser systems based on integrated parametric 3D modeling-simulation framework is proposed,which enhances the flexibility and usability of evaluating the casting process by simulation.Firstly,geometric feature extraction technology is employed to obtain the geometric information of the target casting.Based on this information,an automated design framework for gating and riser systems is established,incorporating multiple structural parameters for real-time process control.Subsequently,the simulation results for various structural parameters are analyzed,and the influence of these parameters on casting formation is thoroughly investigated.Finally,the optimal design scheme is generated and validated through experimental verification.Simulation analysis and experimental results show that using a larger gate neck(24 mm in side length) and external risers promotes a more uniform temperature distribution and a more stable flow state,effectively eliminating shrinkage cavities and enhancing process yield by 15%.
基金supported by the Shihezi University High-Level Talents Research Startup Project(Project No.RCZK202521)the National Natural Science Foundation of China(Grant Nos.12271066,11871121,12171405)+1 种基金the Chongqing Natural Science Foundation Joint Fund for Innovation and Development Project(Project No.CSTB2024NSCQLZX0085)the Chongqing Normal University Foundation(Grant No.23XLB018).
文摘This paper investigates ruin,capital injection,and dividends for a two-dimensional risk model.The model posits that surplus levels of insurance companies are governed by a perturbed composite Poisson risk model.This model introduces a dependence between the two surplus levels,present in both the associated perturbations and the claims resulting from common shocks.Critical levels of capital injection and dividends are established for each of the two risks.The surplus levels are observed discretely at fixed intervals,guiding decisions on capital injection,dividends,and ruin at these junctures.This study employs a two-dimensional Fourier cosine series expansion method to approximate the finite time expected discounted operating cost until ruin.The ensuing approximation error is also quantified.The validity and accuracy of the method are corroborated through numerical examples.Furthermore,the research delves into the optimal capital allocation problem.
基金supported by the National Key R&D Program of China(Grant No.2017YFC1502202)the National Science Foundation of China(Grant Nos.41675075 and41875135)the National Key R&D Program of China(Grant No.2016YFA0602101)
文摘This review summarizes the scientific and technical progress in atmospheric modeling in China since 2011,including the dynamical core,model physics,data assimilation,ensemble forecasting,and model evaluation strategies.In terms of the dynamical core,important efforts have been made in the improvement of the existing model formulations and in exploring new modeling approaches that can better adapt to massively parallel computers and global multiscale modeling.With regard to model physics,various achievements in physical representations have been made,especially a trend toward scale-aware parameterization for accommodating the increase of model resolution.In the field of data assimilation,a 4D-Var system has been developed and is operationally used by the National Meteorological Center of China,and its performance is promising.Furthermore,ensemble forecasting has played a more important role in operational forecast systems and progressed in many fundamental techniques.Model evaluation strategies,including key performance metrics and standardized experimental protocols,have been proposed and widely applied to better understand the strengths and weaknesses of the systems,offering key routes for model improvement.The paper concludes with a concise summary of the status quo and a brief outlook in terms of future development.
文摘In this paper, a heavy sea fog episode that occurred over the Yellow Sea on 9 March 2005 is investigated. The sea fog patch, with a spatial scale of several hundred kilometers at its mature stage, reduced visibility along the Shandong Peninsula coast to 100 m or much less at some sites. Satellite images, surface observations and soundings at islands and coasts, and analyses from the Japan Meteorology Agency (JMA) axe used to describe and analyze this event. The analysis indicates that this sea fog can be categorized as advection cooling fog. The main features of this sea fog including fog area and its movement axe reasonably reproduced by the Fifth-generation Pennsylvania State University/National Center for Atmospheric Research Mesoscale Model (MM5). Model results suggest that the formation and evolution of this event can be outlined as: (1) southerly warm/moist advection of low-level air resulted in a strong sea-surface-based inversion with a thickness of about 600 m; (2) when the inversion moved from the warmer East Sea to the colder Yellow Sea, a thermal internal boundary layer (TIBL) gradually formed at the base of the inversion while the sea fog grew in response to cooling and moistening by turbulence mixing; (3) the sea fog developed as the TIBL moved northward and (4) strong northerly cold and dry wind destroyed the TIBL and dissipated the sea fog. The principal findings of this study axe that sea fog forms in response to relatively persistent southerly waxm/moist wind and a cold sea surface, and that turbulence mixing by wind shear is the primary mechanism for the cooling and moistening the marine layer. In addition, the study of sensitivity experiments indicates that deterministic numerical modeling offers a promising approach to the prediction of sea fog over the Yellow Sea but it may be more efficient to consider ensemble numerical modeling because of the extreme sensitivity to model input.
文摘Recent advances in three-dimensional numerical simulations of mantle convection have aided in approximately reproducing continental movement since the Pangea breakup at 200 Ma. These have also led to a better understanding of the thermal and mechanical coupling between mantle convection and surface plate motion and predictions of the configuration of the next supercontinent. The simulations of mantle convection from 200 Ma to the present reveals that the development of large-scale cold mantle downwellings in the North Tethys Ocean at the earlier stage of the Pangea breakup triggered the northward movement of the Indian subcontinent. The model of high temperature anomaly region beneath Pangea resulting from the thermal insulation effect support the breakup of Pangea in the real Earth time scale, as also suggested in previous geological and geodynamic models. However, considering the low radioactive heat generation rate of the depleted upper mantle, the high temperature anomaly region might have been generated by upwelling plumes with contribution of deep subducted TTG(tonalite-trondhjemite-granite) materials enriched in radiogenic elements. Integrating the numerical results of mantle convection from 200 Ma to the present, and from the present to the future, it is considered that the mantle drag force acting on the base of continents may be comparable to the slab pull force, which implies that convection in the shallower part of the mantle is strongly coupled with surface plate motion.
基金funded by a grant from Natural Sciences and Engineering Research Council of Canada (NSERC)the authors would like to acknowledge the Niobec mine (Saint-Honoré, QuébecQuébec)
文摘Stress-induced failure is among the most common causes of instability in Canadian deep underground mines.Open stoping is the most widely practiced underground excavation method in these mines,and creates large stopes which are subjected to stress-induced failure.The probability of failure(POF)depends on many factors,of which the geometry of an open stope is especially important.In this study,a methodology is proposed to assess the effect of stope geometrical parameters on the POF,using numerical modelling.Different ranges for each input parameter are defined according to previous surveys on open stope geometry in a number of Canadian underground mines.A Monte-Carlo simulation technique is combined with the finite difference code FLAC3D,to generate model realizations containing stopes with different geometrical features.The probability of failure(POF)for different categories of stope geometry,is calculated by considering two modes of failure;relaxation-related gravity driven(tensile)failure and rock mass brittle failure.The individual and interactive effects of stope geometrical parameters on the POF,are analyzed using a general multi-level factorial design.Finally,mathematical optimization techniques are employed to estimate the most stable stope conditions,by determining the optimal ranges for each stope’s geometrical parameter.
基金supported by National Natural Science Foundation of China(72288101,72361137002,and 72101018)the Dutch Research Council(NWO Grant 482.22.01).
文摘Vehicle electrification,an important method for reducing carbon emissions from road transport,has been promoted globally.In this study,we analyze how individuals adapt to this transition in transportation and its subsequent impact on urban structure.Considering the varying travel costs associated with electric and fuel vehicles,we analyze the dynamic choices of households concerning house locations and vehicle types in a two-dimensional monocentric city.A spatial equilibrium is developed to model the interactions between urban density,vehicle age and vehicle type.An agent-based microeconomic residential choice model dynamically coupled with a house rent market is developed to analyze household choices of home locations and vehicle energy types,considering vehicle ages and competition for public charging piles.Key findings from our proposed models show that the proportion of electric vehicles(EVs)peaks at over 50%by the end of the first scrappage period,accompanied by more than a 40%increase in commuting distance and time compared to the scenario with only fuel vehicles.Simulation experiments on a theoretical grid indicate that heterogeneity-induced residential segregation can lead to urban sprawl and congestion.Furthermore,households with EVs tend to be located farther from the city center,and an increase in EV ownership contributes to urban expansion.Our study provides insights into how individuals adapt to EV transitions and the resulting impacts on home locations and land use changes.It offers a novel perspective on the dynamic interactions between EV adoption and urban development.