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.展开更多
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.展开更多
The purpose of this study was to investigate the effect of bolt profile on load transfer mechanism of fully grouted bolts in jointed rocks using analytical and numerical methods. Based on the analytical method with de...The purpose of this study was to investigate the effect of bolt profile on load transfer mechanism of fully grouted bolts in jointed rocks using analytical and numerical methods. Based on the analytical method with development of methods, a new model is presented. To validate the analytical model, five different profiles modeled by ANSYS software. The profile of rock bolts T3 and T4with load transfer capacity,respectively 180 and 195 kN in the jointed rocks was selected as the optimum profiles. Finally, the selected profiles were examined in Tabas Coal Mine. FLAC analysis indicates that patterns 6+7 with2 NO flexi bolt 4 m better than other patterns within the faulted zone.展开更多
Statistical distributions are used to model wind speed,and the twoparameters Weibull distribution has proven its effectiveness at characterizing wind speed.Accurate estimation of Weibull parameters,the scale(c)and sha...Statistical distributions are used to model wind speed,and the twoparameters Weibull distribution has proven its effectiveness at characterizing wind speed.Accurate estimation of Weibull parameters,the scale(c)and shape(k),is crucial in describing the actual wind speed data and evaluating the wind energy potential.Therefore,this study compares the most common conventional numerical(CN)estimation methods and the recent intelligent optimization algorithms(IOA)to show how precise estimation of c and k affects the wind energy resource assessments.In addition,this study conducts technical and economic feasibility studies for five sites in the northern part of Saudi Arabia,namely Aljouf,Rafha,Tabuk,Turaif,and Yanbo.Results exhibit that IOAs have better performance in attaining optimal Weibull parameters and provided an adequate description of the observed wind speed data.Also,with six wind turbine technologies rating between 1 and 3MW,the technical and economic assessment results reveal that the CN methods tend to overestimate the energy output and underestimate the cost of energy($/kWh)compared to the assessments by IOAs.The energy cost analyses show that Turaif is the windiest site,with an electricity cost of$0.016906/kWh.The highest wind energy output is obtained with the wind turbine having a rated power of 2.5 MW at all considered sites with electricity costs not exceeding$0.02739/kWh.Finally,the outcomes of this study exhibit the potential of wind energy in Saudi Arabia,and its environmental goals can be acquired by harvesting wind energy.展开更多
In wall-bounded turbulent flow calculations, the past focus has been directed to the modelling of the Reynolds-stress gradients. Not much attention has been paid to the effects of the numerical methods used to calcula...In wall-bounded turbulent flow calculations, the past focus has been directed to the modelling of the Reynolds-stress gradients. Not much attention has been paid to the effects of the numerical methods used to calculate these terms and the modelled equations. Discrepancies between model calculations and measurements are quite often attributed to incorrect modelling, while the suitability and accuracy of the numerical methods used are seldom scrutinized. Instead, alternate near-wall and Reynolds-stress models are proposed to remedy the incorrect turbulent flow calculations. On the other hand, if care is not taken in the numerical treatment of the Reynolds-stress gradient terms, physically unrealistic results and solution instability could occur. Previous studies by the author and his collaborators on the effects of numerical methods have shown that some of the more commonly used numerical methods could enhance numerical stability in the solution procedure but would introduce considerable inaccuracy to the results. The flow cases chosen to demonstrate these inaccuracies are a backstep flow and flow in a square duct, where flow complexities are present. The current investigation attempts to show that the above-mentioned effects of numerical methods could also occur in the calculation of a developing plane channel flow, where flow complexities are absent. In addition, this study shows that the results thus obtained lead to a predicted skin friction coefficient that is influenced more by the numerical method used than by the turbulence model invoked. Together, these results show that numerical treatment of the Reynolds-stress gradients in the equations play an important role, even for a developing plane channel flow.展开更多
Viscous fluid flows contain abundant "physical phenomena and the viscous fluid dynamics is of wide applications in the fields of natural and engineering sciences. After the basic equations of viscousfluiddynamics...Viscous fluid flows contain abundant "physical phenomena and the viscous fluid dynamics is of wide applications in the fields of natural and engineering sciences. After the basic equations of viscousfluiddynamics (i.e., the Navier-Stokes equations) came out, one of the most important contributions to the discipline was the boundary layer (BL) theory and the BL equations presented by Prandtl展开更多
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.展开更多
To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based sim...To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.展开更多
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.展开更多
In the field of discretization-based meshfree/meshless methods,the improvements in the higher-order consistency,stability,and computational efficiency are of great concerns in computational science and numerical solut...In the field of discretization-based meshfree/meshless methods,the improvements in the higher-order consistency,stability,and computational efficiency are of great concerns in computational science and numerical solutions to partial differential equations.Various alternative numerical methods of the finite particle method(FPM)frame have been extended from mathematical theories to numerical applications separately.As a comprehensive numerical scheme,this study suggests a unified resolved program for numerically investigating their accuracy,stability,consistency,computational efficiency,and practical applicability in industrial engineering contexts.The high-order finite particle method(HFPM)and corrected methods based on the multivariate Taylor series expansion are constructed and analyzed to investigate the whole applicability in different benchmarks of computational fluid dynamics.Specifically,four benchmarks are designed purposefully from statical exact solutions to multifaceted hydrodynamic tests,which possess different numerical performances on the particle consistency,numerical discretized forms,particle distributions,and transient time evolutional stabilities.This study offers a numerical reference for the current unified resolved program.展开更多
The modeling of crack growth in three-dimensional(3D)space poses significant challenges in rock mechanics due to the complex numerical computation involved in simulating crack propagation and interaction in rock mater...The modeling of crack growth in three-dimensional(3D)space poses significant challenges in rock mechanics due to the complex numerical computation involved in simulating crack propagation and interaction in rock materials.In this study,we present a novel approach that introduces a 3D numerical manifold method(3D-NMM)with a geometric kernel to enhance computational efficiency.Specifically,the maximum tensile stress criterion is adopted as a crack growth criterion to achieve strong discontinuous crack growth,and a local crack tracking algorithm and an angle correction technique are incorporated to address minor limitations of the algorithm in a 3D model.The implementation of the program is carried out in Python,using object-oriented programming in two independent modules:a calculation module and a crack module.Furthermore,we propose feasible improvements to enhance the performance of the algorithm.Finally,we demonstrate the feasibility and effectiveness of the enhanced algorithm in the 3D-NMM using four numerical examples.This study establishes the potential of the 3DNMM,combined with the local tracking algorithm,for accurately modeling 3D crack propagation in brittle rock materials.展开更多
The best way to check the validity of our theories(models)is by direct comparison with the experiment(observations).However,this process suffers from numerical inaccuracies,which are not frequently studied and often r...The best way to check the validity of our theories(models)is by direct comparison with the experiment(observations).However,this process suffers from numerical inaccuracies,which are not frequently studied and often remain mostly unknown.In this study,we focus on addressing the numerical inaccuracies intrinsic to the process of comparing theory and observations.To achieve this goal,we built four-dimensional(4D)spectral grids for Wolf–Rayet stars(WC and WN spectral classes)and blue supergiants characterized by low metallicity similar to that of the Small Magellanic Cloud.In contrast to lighter(three-dimensional)grids,which rely on a priori assumptions about certain stellar parameters(e.g.,wind velocity)and thus have limited applicability,our 4D grids vary four independent parameters,enabling more flexible and broadly applicable spectral fitting.Utilizing these 4D grids,we developed and validated a fitting approach facilitating direct fits to observed spectra.Through rigorous testing on designated“test”models,we demonstrated that the numerical precision of derived essential stellar parameters,including effective temperature,mass-loss rate,luminosity,and wind velocity,is better than 0.05 dex.Furthermore,we explored the influence of unaccounted factors,including variations in the metal abundances,wind acceleration laws,and clumping,on the precision of the derived parameters.The results indicate that the first two factors have the strongest influence on the numerical accuracy of the derived stellar parameters.Variations in abundances predominantly influenced the mass-loss rate for weak-wind scenarios,while effective temperature and luminosity remained robust.We found that the wind acceleration law influences the numerical uncertainty of the derived wind parameters mostly for models with weak winds.Interestingly,different degrees of clumping demonstrated good precision for spectra with strong winds,contrasting with a decrease in the precision for weak-wind cases.We found also that the accuracy of our approach depends on spectral range and the inclusion of ultraviolet spectral range improves the precision of derived parameters,especially for an object with weak winds.展开更多
This study investigates the technique of variational calculus applied to estimate the slope stability considering the mechanism of planar failure.The critical plane failure surface should be determined because it theo...This study investigates the technique of variational calculus applied to estimate the slope stability considering the mechanism of planar failure.The critical plane failure surface should be determined because it theoretically indicates the most unfavorable plane to be considered when stabilizing a slope to rectify the instability generated by several statistically possible planes.This generates integrals that can be solved by numerical methods,such as the Newton Cotes and the finite differences methods.Additionally,a system of nonlinear equations is obtained and solved.The surface of the critical planar failure is determined by applying the condition of transversality in mobile boundaries,for which various examples are provided.The number of slices is varied in one of the examples,while the surface of the critical planar failure is determined in the others.Results are compared using analytical methods through axis rotations.All the results obtained by considering normal stress,safety factors,and critical planar failure are nearly the same;however,in this research,a study is carried out for“n”number of slices using programming methods.Sub-routines are important because they can be applied in slopes with different geometry,surcharge,interstitial pressure,and pseudo-static load.展开更多
Efficiency and accuracy are two major concerns in numerical solutions of the Poisson-Boltzmann equation for applications in chemistry and biophysics.Recent developments in boundary element methods,interface methods,ad...Efficiency and accuracy are two major concerns in numerical solutions of the Poisson-Boltzmann equation for applications in chemistry and biophysics.Recent developments in boundary element methods,interface methods,adaptive methods,finite element methods,and other approaches for the Poisson-Boltzmann equation as well as related mesh generation techniques are reviewed.We also discussed the challenging problems and possible future work,in particular,for the aim of biophysical applications.展开更多
Numerical predictions are made for Laminar Forced convection heat transfer with and without buoyancy effects for Supercritical Nitrogen flowing over an isothermal horizontal flat plate with a heated surface facing dow...Numerical predictions are made for Laminar Forced convection heat transfer with and without buoyancy effects for Supercritical Nitrogen flowing over an isothermal horizontal flat plate with a heated surface facing downwards.Computations are performed by varying the value ofΔT from5 to 30 K and P_(∞)/P_(cr)ratio from1.1 to 1.5.Variation of all the thermophysical properties of supercritical Nitrogen is considered.The wall temperatures are chosen in such a way that two values of Tw are less than T∗(T*is the temperature at which the fluid has a maximum value of Cp for the given pressure),one value equal to T∗and two values greater than T∗.Three different values of U∞are used to obtain Re∞range of 3.6×10_(4)to 4.74×10^(5)for forced convection without buoyancy effects and Gr_(∞)/Re^(2)_(∞)range of 0.011 to 3.107 for the case where buoyancy effects are predominant.Six different forms of correlations are proposed based on numerical predictions and are compared with actual numerical predictions.It has been found that in all six forms of correlations,the maximum deviations are found to occur in those cases where the pseudocritical temperature TT∗lies between the wall temperature and bulk fluid temperature.展开更多
Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on t...Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.展开更多
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.展开更多
To ensure the safe performance of deep-sea mining vehicles(DSMVs),it is necessary to study the mechanical characteristics of the interaction between the seabed soil and the track plate.The rotation and digging motions...To ensure the safe performance of deep-sea mining vehicles(DSMVs),it is necessary to study the mechanical characteristics of the interaction between the seabed soil and the track plate.The rotation and digging motions of the track plate are important links in the contact between the driving mechanism of the DSMV and seabed soil.In this study,a numerical simulation is conducted using the coupled Eulerian–Lagrangian(CEL)large deformation numerical method to investigate the interaction between the track plate of the DSMV and the seabed soil under two working conditions:rotating condition and digging condition.First,a soil numerical model is established based on the elastoplastic mechanical characterization using the basic physical and mechanical properties of the seabed soil obtained by in situ sampling.Subsequently,the soil disturbance mechanism and the dynamic mechanical response of the track plate under rotating and digging conditions are obtained through the analysis of the sensitivity of the motion parameters,the grouser structure,the layered soil features and the soil heterogeneity.The results indicate that the above parameters remarkably influence the interaction between the DSMV and the seabed soil.Therefore,it is important to consider the rotating and digging motion of the DSMV in practical engineering to develop a detailed optimization design of the track plate.展开更多
As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accura...As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.展开更多
文摘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.
基金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.
文摘The purpose of this study was to investigate the effect of bolt profile on load transfer mechanism of fully grouted bolts in jointed rocks using analytical and numerical methods. Based on the analytical method with development of methods, a new model is presented. To validate the analytical model, five different profiles modeled by ANSYS software. The profile of rock bolts T3 and T4with load transfer capacity,respectively 180 and 195 kN in the jointed rocks was selected as the optimum profiles. Finally, the selected profiles were examined in Tabas Coal Mine. FLAC analysis indicates that patterns 6+7 with2 NO flexi bolt 4 m better than other patterns within the faulted zone.
基金The author extends his appreciation to theDeputyship forResearch&Innovation,Ministry of Education,Saudi Arabia for funding this research work through the Project Number(QUIF-4-3-3-33891)。
文摘Statistical distributions are used to model wind speed,and the twoparameters Weibull distribution has proven its effectiveness at characterizing wind speed.Accurate estimation of Weibull parameters,the scale(c)and shape(k),is crucial in describing the actual wind speed data and evaluating the wind energy potential.Therefore,this study compares the most common conventional numerical(CN)estimation methods and the recent intelligent optimization algorithms(IOA)to show how precise estimation of c and k affects the wind energy resource assessments.In addition,this study conducts technical and economic feasibility studies for five sites in the northern part of Saudi Arabia,namely Aljouf,Rafha,Tabuk,Turaif,and Yanbo.Results exhibit that IOAs have better performance in attaining optimal Weibull parameters and provided an adequate description of the observed wind speed data.Also,with six wind turbine technologies rating between 1 and 3MW,the technical and economic assessment results reveal that the CN methods tend to overestimate the energy output and underestimate the cost of energy($/kWh)compared to the assessments by IOAs.The energy cost analyses show that Turaif is the windiest site,with an electricity cost of$0.016906/kWh.The highest wind energy output is obtained with the wind turbine having a rated power of 2.5 MW at all considered sites with electricity costs not exceeding$0.02739/kWh.Finally,the outcomes of this study exhibit the potential of wind energy in Saudi Arabia,and its environmental goals can be acquired by harvesting wind energy.
文摘In wall-bounded turbulent flow calculations, the past focus has been directed to the modelling of the Reynolds-stress gradients. Not much attention has been paid to the effects of the numerical methods used to calculate these terms and the modelled equations. Discrepancies between model calculations and measurements are quite often attributed to incorrect modelling, while the suitability and accuracy of the numerical methods used are seldom scrutinized. Instead, alternate near-wall and Reynolds-stress models are proposed to remedy the incorrect turbulent flow calculations. On the other hand, if care is not taken in the numerical treatment of the Reynolds-stress gradient terms, physically unrealistic results and solution instability could occur. Previous studies by the author and his collaborators on the effects of numerical methods have shown that some of the more commonly used numerical methods could enhance numerical stability in the solution procedure but would introduce considerable inaccuracy to the results. The flow cases chosen to demonstrate these inaccuracies are a backstep flow and flow in a square duct, where flow complexities are present. The current investigation attempts to show that the above-mentioned effects of numerical methods could also occur in the calculation of a developing plane channel flow, where flow complexities are absent. In addition, this study shows that the results thus obtained lead to a predicted skin friction coefficient that is influenced more by the numerical method used than by the turbulence model invoked. Together, these results show that numerical treatment of the Reynolds-stress gradients in the equations play an important role, even for a developing plane channel flow.
文摘Viscous fluid flows contain abundant "physical phenomena and the viscous fluid dynamics is of wide applications in the fields of natural and engineering sciences. After the basic equations of viscousfluiddynamics (i.e., the Navier-Stokes equations) came out, one of the most important contributions to the discipline was the boundary layer (BL) theory and the BL equations presented by Prandtl
基金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.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52271317 and 52071149)the Fundamental Research Funds for the Central Universities(HUST:2019kfy XJJS007)。
文摘To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.
基金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 the National Natural Science Foundation of China(No.12002290)。
文摘In the field of discretization-based meshfree/meshless methods,the improvements in the higher-order consistency,stability,and computational efficiency are of great concerns in computational science and numerical solutions to partial differential equations.Various alternative numerical methods of the finite particle method(FPM)frame have been extended from mathematical theories to numerical applications separately.As a comprehensive numerical scheme,this study suggests a unified resolved program for numerically investigating their accuracy,stability,consistency,computational efficiency,and practical applicability in industrial engineering contexts.The high-order finite particle method(HFPM)and corrected methods based on the multivariate Taylor series expansion are constructed and analyzed to investigate the whole applicability in different benchmarks of computational fluid dynamics.Specifically,four benchmarks are designed purposefully from statical exact solutions to multifaceted hydrodynamic tests,which possess different numerical performances on the particle consistency,numerical discretized forms,particle distributions,and transient time evolutional stabilities.This study offers a numerical reference for the current unified resolved program.
基金supported by the National Natural Science Foundation of China(Grant Nos.42172312 and 52211540395)support from the Institut Universitaire de France(IUF).
文摘The modeling of crack growth in three-dimensional(3D)space poses significant challenges in rock mechanics due to the complex numerical computation involved in simulating crack propagation and interaction in rock materials.In this study,we present a novel approach that introduces a 3D numerical manifold method(3D-NMM)with a geometric kernel to enhance computational efficiency.Specifically,the maximum tensile stress criterion is adopted as a crack growth criterion to achieve strong discontinuous crack growth,and a local crack tracking algorithm and an angle correction technique are incorporated to address minor limitations of the algorithm in a 3D model.The implementation of the program is carried out in Python,using object-oriented programming in two independent modules:a calculation module and a crack module.Furthermore,we propose feasible improvements to enhance the performance of the algorithm.Finally,we demonstrate the feasibility and effectiveness of the enhanced algorithm in the 3D-NMM using four numerical examples.This study establishes the potential of the 3DNMM,combined with the local tracking algorithm,for accurately modeling 3D crack propagation in brittle rock materials.
文摘The best way to check the validity of our theories(models)is by direct comparison with the experiment(observations).However,this process suffers from numerical inaccuracies,which are not frequently studied and often remain mostly unknown.In this study,we focus on addressing the numerical inaccuracies intrinsic to the process of comparing theory and observations.To achieve this goal,we built four-dimensional(4D)spectral grids for Wolf–Rayet stars(WC and WN spectral classes)and blue supergiants characterized by low metallicity similar to that of the Small Magellanic Cloud.In contrast to lighter(three-dimensional)grids,which rely on a priori assumptions about certain stellar parameters(e.g.,wind velocity)and thus have limited applicability,our 4D grids vary four independent parameters,enabling more flexible and broadly applicable spectral fitting.Utilizing these 4D grids,we developed and validated a fitting approach facilitating direct fits to observed spectra.Through rigorous testing on designated“test”models,we demonstrated that the numerical precision of derived essential stellar parameters,including effective temperature,mass-loss rate,luminosity,and wind velocity,is better than 0.05 dex.Furthermore,we explored the influence of unaccounted factors,including variations in the metal abundances,wind acceleration laws,and clumping,on the precision of the derived parameters.The results indicate that the first two factors have the strongest influence on the numerical accuracy of the derived stellar parameters.Variations in abundances predominantly influenced the mass-loss rate for weak-wind scenarios,while effective temperature and luminosity remained robust.We found that the wind acceleration law influences the numerical uncertainty of the derived wind parameters mostly for models with weak winds.Interestingly,different degrees of clumping demonstrated good precision for spectra with strong winds,contrasting with a decrease in the precision for weak-wind cases.We found also that the accuracy of our approach depends on spectral range and the inclusion of ultraviolet spectral range improves the precision of derived parameters,especially for an object with weak winds.
文摘This study investigates the technique of variational calculus applied to estimate the slope stability considering the mechanism of planar failure.The critical plane failure surface should be determined because it theoretically indicates the most unfavorable plane to be considered when stabilizing a slope to rectify the instability generated by several statistically possible planes.This generates integrals that can be solved by numerical methods,such as the Newton Cotes and the finite differences methods.Additionally,a system of nonlinear equations is obtained and solved.The surface of the critical planar failure is determined by applying the condition of transversality in mobile boundaries,for which various examples are provided.The number of slices is varied in one of the examples,while the surface of the critical planar failure is determined in the others.Results are compared using analytical methods through axis rotations.All the results obtained by considering normal stress,safety factors,and critical planar failure are nearly the same;however,in this research,a study is carried out for“n”number of slices using programming methods.Sub-routines are important because they can be applied in slopes with different geometry,surcharge,interstitial pressure,and pseudo-static load.
基金the NIH,NSF,the Howard Hughes Medical Institute,National Biomedical Computing Resource,the NSF Center for Theoretical Biological Physics,SDSC,the W.M.Keck Foundation,and Accelrys,Inc.Michael Holst was supported in part by NSF Awards 0411723,0511766,and 0225630,and DOE Awards DEFG02-05ER25707 and DE-FG02-04ER25620.
文摘Efficiency and accuracy are two major concerns in numerical solutions of the Poisson-Boltzmann equation for applications in chemistry and biophysics.Recent developments in boundary element methods,interface methods,adaptive methods,finite element methods,and other approaches for the Poisson-Boltzmann equation as well as related mesh generation techniques are reviewed.We also discussed the challenging problems and possible future work,in particular,for the aim of biophysical applications.
文摘Numerical predictions are made for Laminar Forced convection heat transfer with and without buoyancy effects for Supercritical Nitrogen flowing over an isothermal horizontal flat plate with a heated surface facing downwards.Computations are performed by varying the value ofΔT from5 to 30 K and P_(∞)/P_(cr)ratio from1.1 to 1.5.Variation of all the thermophysical properties of supercritical Nitrogen is considered.The wall temperatures are chosen in such a way that two values of Tw are less than T∗(T*is the temperature at which the fluid has a maximum value of Cp for the given pressure),one value equal to T∗and two values greater than T∗.Three different values of U∞are used to obtain Re∞range of 3.6×10_(4)to 4.74×10^(5)for forced convection without buoyancy effects and Gr_(∞)/Re^(2)_(∞)range of 0.011 to 3.107 for the case where buoyancy effects are predominant.Six different forms of correlations are proposed based on numerical predictions and are compared with actual numerical predictions.It has been found that in all six forms of correlations,the maximum deviations are found to occur in those cases where the pseudocritical temperature TT∗lies between the wall temperature and bulk fluid temperature.
文摘Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.
文摘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.
基金supported by the Natural Science Foundation of Hainan Province(Grant No.520LH015)the Fundamental Research Funds for the Central Universities and the Major Projects of Strategic Emerging Industries in Shanghai(Grant No.BH3230001).
文摘To ensure the safe performance of deep-sea mining vehicles(DSMVs),it is necessary to study the mechanical characteristics of the interaction between the seabed soil and the track plate.The rotation and digging motions of the track plate are important links in the contact between the driving mechanism of the DSMV and seabed soil.In this study,a numerical simulation is conducted using the coupled Eulerian–Lagrangian(CEL)large deformation numerical method to investigate the interaction between the track plate of the DSMV and the seabed soil under two working conditions:rotating condition and digging condition.First,a soil numerical model is established based on the elastoplastic mechanical characterization using the basic physical and mechanical properties of the seabed soil obtained by in situ sampling.Subsequently,the soil disturbance mechanism and the dynamic mechanical response of the track plate under rotating and digging conditions are obtained through the analysis of the sensitivity of the motion parameters,the grouser structure,the layered soil features and the soil heterogeneity.The results indicate that the above parameters remarkably influence the interaction between the DSMV and the seabed soil.Therefore,it is important to consider the rotating and digging motion of the DSMV in practical engineering to develop a detailed optimization design of the track plate.
基金supported by the National Natural Science Foundation of China(Grant No.42277165)the Fundamental Research Funds for the Central Universities,China University of Geosciences(Wuhan)(Grant No.CUGCJ1821)the National Overseas Study Fund(Grant No.202106410040).
文摘As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.
