Conclusions: Our findings suggest that the NMS-E represents a valuable preoperative diagnostic tool for endometriosis, effectively correlating with the disease's severity and surgical outcomes. Incorporating the NMS-E into clinical practice could significantly enhance the management of endometriosis by addressing current diagnostic limitations and guiding surgical planning.

This study aims to precisely characterize the use of magnetoelectric nanoparticles (MENPs) for stimulating peripheral nerves. Numerical methods were employed to quantify the interaction between MENPs and nervous tissue. The influence of MENPs orientation, concentration and distance was assessed in terms of the external potential distribution exerted by the MENPs, the amplification of the exerted MENPs stimulation required to excite the neural fibers and the current injected into the intracellular space. The results highlight the significance of MENPs concentration for stimulation accuracy and efficiency, the impact of MENPs orientation on the electric potential distribution sensed by the nervous tissue, as well as the importance of the MENPs distance over the fibers’ recruitment. Given the considerable variability in the interaction between MENPs and nerves, our research provides a crucial step towards understanding this interaction, offering quantitative support for the application of MENPs in nervous tissue stimulation.

This study examines geothermal energy extraction potential from the basement rock within the Denver-Julesburg (DJ) Basin, focusing on flow performance and heat extraction efficiency of different geothermal well configurations. It specifically compares U-Shaped, V-Shaped, Inclined V-Shaped, and Pipe-in-Pipe configurations against Enhanced Geothermal Systems (EGS) setups. Through numerical modeling, we evaluated the thermal behavior of these systems under various operational scenarios and fracture conditions. Results suggest that while closed-loop systems offer moderate temperature increases, EGS configurations show substantial potential for high-temperature extraction. This underscores the importance of evaluating well configurations in complex geological settings like the DJ Basin. The insights from this study aid in strategic geothermal energy planning and development, marking significant advancements in geothermal technology and setting a foundation for future explorations and optimizations.

Multi-stage hydraulic fracturing with horizontal well technology has been widely used in shale oil reservoir. After hydraulic fracturing, hydraulic fractures and opened beddings are intertwined, which result in a complex fracture network. As well as the migration of multi-phase fluids during fracturing and soaking processes, which lead to complex flowback performance, and brings dif-ficulty to productivity evaluation and flowback strategies optimization. In this paper, taking the Daqing Gulong shale reservoir as an example, a numerical model for flowback period, which considered oil-water-gas three-phase flow and orthogonal fracture network has been established. The characteristics and influencing factors of flowback performance have been deeply studied, and the flowback modes of shale oil are reasonably optimized. Factors such as PTPG (pseu-do-threshold pressure gradient), opened bedding stress sensitivity, opened bedding permeability and matrix permeability have obvious effects on the flowback performance, resulting in signifi-cant variations in production peaks, high yield periods, and decline rates. In addition, three flowback modes are distinguished based on the rate at which the BHP (bottom hole pressure) drops. These correspond to three types of choke modes. Based on the result of numerical model, three flowback modes are optimized and the impacts of different flowback modes on shale oil flowback are clarified. This study is particularly important for understanding the performance of shale oil flowback and guiding the efficient exploitation of shale oil.

The corrugated board packaging industry is increasingly using advanced numerical tools to design and estimate the load capacity of its products. That is why numerical analyzes are becoming a common standard in this branch of manufacturing. Such trend causes either the use of advanced computational models that take into account the full 3D geometry of the flat and wavy layers of corrugated board, or the use of homogenization techniques to simplify the numerical model. The article presents theoretical considerations that extend the numerical homogenization technique already presented in our previous work. The proposed here homogenization procedure also takes into account the creasing and / or perforation of corrugated board, i.e. processes that undoubtedly weaken the stiffness and strength of the corrugated board locally. However, it is not always easy to estimate how exactly these processes affect the bending or torsional stiffness. What is known for sure is that the degradation of stiffness depends, among other things, on the type of cut, its shape, the depth of creasing, as well as their position or direction in relation to the corrugation direction. The method proposed here can be successfully applied to model smeared degradation in a finite element or to define degraded interface stiffnesses on a crease line or a perforation line.

To improve the effect of multi-point stretch forming of sheet metal, it is proposed in this paper to replace fixed ball head with swinging ball head. According to the multi-point dies with different arrangements, this research establishes the finite element models of the following stretch forming, i.e. fixed ball heads with conventional arrangement, swinging ball heads with conventional arrangement, swinging ball heads with declining staggered arrangement, and swinging ball heads with parallel staggered arrangement and then numerical simulation is performed. The simulation results show that by replacing fixed ball head with swinging ball head, the surface indentation of the formed part was effectively suppressed, the stress and tension strain distribution of the formed part was improved and the forming quality was improved; that the thickness of the elastic pad was reduced, the springback was reduced and the forming accuracy was improved; and that when the ball head was applied to multi-point die with staggered arrangement, better forming result was achieved, where the best forming result was achieved in combining the swinging ball heads with the multi-point die with parallel staggered arrangement. The forming experiments were carried out, and the experimental results were consistent with the trend of numerical simulation results, which verified the correctness of the numerical simulation.

Fires are complex multi-physics problems that span wide spatial scale ranges. Capturing this complexity in computationally affordable numerical simulations for process studies and “outer-loop” techniques (e.g., optimization and uncertainty quantification) is a fundamental challenge in reacting flow research. Further complications arise for propagating fires where a priori knowledge of the fire spread rate and direction is typically not available. In such cases, static mesh refinement at all possible fire locations is a computationally inefficient approach to bridging the wide range of spatial scales relevant to fire behavior. In the present study, we address this challenge by incorporating adaptive mesh refinement (AMR) in fireFoam, an OpenFOAM solver for simulations of complex fire phenomena involving pyrolyzing solid surfaces. The AMR functionality in the extended solver, called fireDyMFoam, is load balanced, models gas, solid, and liquid phases, and allows us to dynamically track regions of interest, thus avoiding inefficient over-resolution of areas far from a propagating flame. We demonstrate the AMR capability and computational efficiency for fire spread on vertical panels, showing that the AMR solver reproduces results obtained using much larger statically refined meshes, but at a substantially reduced computational cost. We then leverage the computational efficiency of the AMR solver to demonstrate an optimization framework for fire suppression based on the open-source Dakota toolkit.

The popularity and growth of additive manufacturing(AM) in recent years are due to its ability to provide innovative solutions to complex manufacturing challenges and its potential to transform the way products are designed, developed, and produced. In this work, the numerical simulation geometry of the model extruder from nozzle is established along with the material extrusion process. The high gravity material extrusion system (HG-MEX) properties and characteristics during filament fabrication are evaluated. This work investigates the effects of gravity on the material extrusion process from nozzle extrusion. This work is unique in the incorporation of fluid surface characteristics to represent the performance and contact with the material extrusion in additive manufacturing. The use of simulation in high gravity material extrusion in additive manufacturing is a meaningful and valuable approach to improving the quality and efficiency of the process.

Fiber-metal laminates (FMLs) such as Kevlar reinforced aluminum laminate (ARALL), Carbon reinforced aluminum laminate (CARALL), and Glass reinforced aluminum laminate (GLARE) offer great potential for weight reduction applications in automobile and aerospace construction. In order to investigate the feasibility for utilizing such materials in the form of laminates, sheet hydro-bulging tests are studied under the condition of uniform blank holder force for three-layered aluminum and aluminum-composite laminates using orthogonal carbon and Kevlar as well as glass fiber in the middle. The experimental results validate the finite element results and they exhibited that the forming limit of glass fiber in the middle is the highest among the studied materials, while carbon fiber material performs the worst. Furthermore, the crack modes are different for the three kinds of fiber materials investigated in the research. This study provides fundamental guidance for the selection of multi-layer sheet materials in the future manufacturing field.

This study employs a high-fidelity numerical framework to determine the plastic material flow patterns and temperature distributions that lead to void formation during friction stir welding (FSW), and to relate the void morphologies to the underlying alloy material properties and process conditions. Three aluminum alloys, viz., 6061-T6, 7075-T6, and 5053-H18 were investigated under varying traverse speeds. The choice of aluminum alloys enables investigation of a wide range of thermal and mechanical properties. The numerical simulations were validated using experimental observations of void morphologies in these three alloys. Temperatures, plastic strain rates, and material flow patterns are considered. The key results from this study are: (1) The predicted stir zone and void morphology are in good agreement with the experimental observations, (2) The temperature and plastic strain-rate maps in the steady-state process conditions show a strong dependency on the alloy type and traverse speeds, (3) The material velocity contours provide a good insight into the material flow in the stir zone for the FSW process conditions that result in voids as well as those that do not result in voids. The numerical model and the ensuing parametric studies presented in this work provide a framework for understanding material flow under different process conditions in aluminum alloys, and potentially in other alloys. Furthermore, the utility of the numerical model for making quantitative predictions and investigating different process parameters to reduce void formation is demonstrated.

This paper aims to present a robust computational technique utilizing finite difference schemes for accurately solving time fractional reaction-diffusion models, which are prevalent in chemical and biological phenomena. The time-fractional derivative is treated in the Caputo sense, addressing both linear and nonlinear scenarios. The proposed schemes were rigorously evaluated for stability and convergence. Additionally, the effectiveness of the developed schemes was validated through various linear and nonlinear models, including the Allen-Cahn equation, the KPP-Fisher equation, and the Complex Ginzburg-Landau oscillatory problem. These models were tested in one, two, and three-dimensional spaces to investigate the diverse patterns and dynamics that emerge. Comprehensive numerical results were provided, showcasing different cases of the fractional order parameter, highlighting the schemes' versatility and reliability in capturing complex behaviors in fractional reaction-diffusion dynamics.

In this paper we consider constructions of first derivative approximations using the generating function. The weights of the approximations contain the powers of a parameter whose modulus is less than one. The values of the initial weights are determined, and the convergence and order of the approximations are proved. The paper discusses applications of approximations of first derivative for numerical solution of ordinary and partial differential equations and proposes an algorithm for fast computation of the numerical solution. Proofs of the convergence and accuracy of the numerical solutions are presented and the performance of the numerical methods considered is compared with the Euler method. The main goal of constructing approximations for integer-order derivatives of this type is their application in deriving high-order approximations for fractional derivatives, whose weights have specific properties. The paper proposes the construction of an approximation for the fractional derivative and its application for numerically solving fractional differential equations. The theoretical results for the accuracy and order of the numerical methods are confirmed by the experimental results presented in the paper.

Human number understanding is thought to rely on the analogue number system (ANS), working according to Weber’s law. We propose an alternative account, suggesting that symbolic mathematical knowledge is based on a discrete semantic system (DSS), a representation that stores values in a semantic network, similar to the mental lexicon or to a conceptual network. Here, focusing on the phenomena of numerical distance and size effects in comparison tasks, first we discuss how a DSS model could explain these numerical effects. Second, we demonstrate that DSS model can give quantitatively as appropriate a description of the effects as the ANS model. Finally, we show that symbolic numerical size effect is mainly influenced by the frequency of the symbols, and not by the ratios of their values. This last result suggests that numerical distance and size effects cannot be caused by the ANS, while the DSS model might be the alternative approach that can explain the frequency-based size effect.

This study is dedicated to understanding and simulating the mechanical properties of concrete, with a specific focus on the mesoscopic scale. Investigating concrete at this level involves examining its composition as a heterogeneous amalgamation comprising mortar, aggregates, and the interfacial transition zone (ITZ). Such an approach provides a detailed examination of the material's behavior and attributes. Numerical models are utilized, leveraging the finite element method (FEM), to examine the behavior of concrete. The study employs MATLAB programming to develop three-dimensional models, subsequently subjected to FEM analysis. Various mesoscopic Representative Volume Elements (RVEs) are formulated, encompassing varied shapes and dimensions. The MATLAB framework generates input files for numerical FEM simulations, replicating compression strength tests. As complexity increases with the inclusion of the ITZ, prismatic RVEs are developed. The proposed mesoscopic model establishes a foundational framework for a numerical simulation methodology tailored to laboratory compression tests. It provides detailed insights into concrete behavior, elucidating deformation, and fracture mechanisms. Although not a complete substitute for experimental methods, these models offer a cost-effective and expeditious alternative, pinpointing vulnerable areas and exploring the implications of additional materials on concrete behavior. Experimental data and virtual tests pave the way for mitigating carbon footprint and improving concrete sustainability.

Microscale laser dynamic forming, as a novel high-speed microforming technique, can overcome the shortcomings of traditional microforming methods. However, in practical applications, laser dynamic microforming technology is often affected by the rebound behavior of the workpiece surface, limiting the further improvement of processing quality and efficiency. This paper aims to reduce the rebound effect in laser dynamic microforming by using multi-pulse laser shock loading. The forming results of workpieces under different laser energies and laser impact numbers were studied using experimental and numerical simulation methods. Numerical simulations of the forming results after multiple laser shocks were conducted using ANSYS/LS-DYNA software. These numerical simulation results were then experimentally validated and compared. Surface morphology and microstructure of the workpieces were characterized using a confocal microscope and scanning electron microscope, and energy dispersive spectroscopy (EDS) was used to analyze the chemical element content changes in the collision regions at the bottoms of the workpieces after multi-pulse loading, revealing the collision behavior patterns during the forming process. Finally, the forming laws of workpieces under multiple laser shocks were summarized.

The ground consists of many layers of soil with different properties. The propagation speed and path of seismic waves are affected by different soil layers. It is necessary to understand that layered soil exhibits different dynamic behaviors and responses under the action of seismic waves. This study utilized weathered soil and silica sand as materials to create multi-layered soil conditions with varying degrees of compaction. By conducting 1 g shaking table test on multi-layered soil, the interaction and influence between different soil layers under different earthquake conditions were observed. The approach of numerical analysis aims to validate the experimental results and provide an in-depth understanding of the dynamic behaviors of multi-layered soil grounds during seismic events. The results of this study show that the amplification effect of seismic waves is related to the number of soil layers and the degree of compaction of the soil layers. This indicates that the multi-layered soil ground and the behavior of the soil layers play an important role in the propagation and impact of seismic waves, and this amplification effect is of great significance in actual seismic disaster risk assessment and design.

The United Arab Emirates (UAE) is significantly dependent on desalinated water and groundwater resource, which is expensive and highly energy intensive. Despite the scarce water resource, only 54% of the recycled water was reused in 2015. In this study, an “Oasis” complex comprised of Sustainable Farming Compartments (SFCs) was proposed for reusing treated wastewater to decrease the ambient temperature of the SFC via an evaporative cooling system. A prototype SFC with half the original scale (width = 1.8 m, depth = 1.5 m, front height = 1.2 m back height = 0.9 m) was designed, built, and tested in an environmentally controlled laboratory and field site to evaluate the feasibility and effectiveness of the SFC under the climatic conditions in Abu Dhabi. Based on the experimental results, the temperature drops obtained from the SFC in the laboratory and field site were 5 ̊C at initial relative humidity of 60% and 7- 15 ̊C at initial relative humidity of 50%, respectively. An energy simulation using dynamic numerical simulations was performed in comparison to the results of the experiment. The energy-based dynamic simulation shows good agreement with the experimental results. The total power consumption of the SFC system was approximately three and a half times lower than that of an electrical air conditioner.

The problems of formation and propagation of surface seismic waves in vertically layered and horizontally inhomogeneous media are studied. A mathematical model has been developed for the propagation of surface seismic waves in an inhomogeneous medium based on the use of the fundamental principles of wave dynamics and Fourier transforms. The existence of various types of groups of these waves is proved, and their main properties are studied. It has been established that the reflection coefficients of different groups of a dispersive wave on models of half-spaces (a layered, inhomogeneous medium) with a sharp velocity differentiation are characterized by the same dependence on the angle of incidence as in the absence of dispersion. In the region of the critical angle, which depends on the period of the incident wave, a local decrease in the intensity of wave reflection is observed, which is associated with the origin of the head dispersive wave. It has been found that the reflection coefficients of Rayleigh waves from low-velocity layers, thin (d / λ2 = 0.205) and thick (d / λ2 = 0.82) are approximately 0.2 ± 0.04 and practically do not change with increasing angle of incidence. On the basis of theoretical studies, it was revealed that reflected surface waves carry information both about the elastic properties of the reflecting medium and about the dimensions of the reflecting horizontal inhomogeneity. An experimental dispersion curve of the group velocity of Rayleigh waves on the Central Asia-Caucasus path has been obtained. The seismic stress state of the "subgrade-base" system was studied. Based on the analysis of the results of numerical experiments, it has been established that in order to obtain a complete picture of the seismically stressed state of the subgrade, it is necessary to take into account the contribution to it of each component of the accelerogram of a real earthquake.

Small and micro hydropower represents an attractive solution for electricity generation, with low cost and low environmental impact. The pump-as-turbine (PAT) approach has promise in this application owing to its low purchase and maintenance costs. In this paper, a new method to predict the inverse characteristic of industrial centrifugal pumps is presented. This method is based on results of simulations performed with commercial three-dimensional CFD software. Model results have been first validated in pumping mode using data supplied by pump manufacturers. Then, results have been compared to experimental data for a pump running in reverse condition. Experimentation has been performed on a dedicated test bench installed in the Department of Civil Construction and Environmental Engineering of the University of Naples Federico II. Three different pumps, with different specific speeds, have been analyzed. Using the model results, the inverse characteristic and the best efficiency point have been evaluated. Finally, results of this methodology have been compared to prediction methods available in the literature.

A case of low level windshear at Beijing Daxing International Airport is studied in the present pa-per. It is related to the synoptic scale/mesoscale changes of the airflow pattern over northern China. The shear was captured well by the Doppler LIDAR at the airport and timely alerts had been issued. This study considers the technical feasibility of forecasting the shear feature using a high resolution numerical weather prediction model with a lead time of several hours. It turns out that the shear feature could be forecast well, and is generally consistent with the actual observations. However, there are still discrepancies between reality and the simulation, such as the strength of the shear and persistency of the feature. More cases of windshear at Daxing airport would be studied in the fu-ture in order to understand the nature of windshear over there.

Frontal Polymerization (FP) is recognized as a polymer synthesis method capable of generating a highly localized and self-propagating exothermic front, considered a more rapid and energy-efficient manufacturing scheme for polymer-based fiber-reinforced composite materials. This study conducts a numerical investigation based on the Finite Element Method into the initiation and propagation of frontal polymerization of acrylamide. Utilizing Differential Scanning Calorimetry (DSC) experiments, the study phenomenologically tests the curing kinetics reactions of acrylamide (AM) in deep eutectic solvent (DES)-based polymerizable systems, thereby establishing a curing kinetics model for the autocatalytic reaction of polyacrylamide. By coupling the heat conduction diffusion equation of the frontal polymerization reaction, the Finite Element Method is employed to solve the reaction-diffusion model, aiming to explore the evolution of temperature and degree of curing during the manufacturing process. The results demonstrate the reliability and accuracy of the nth-order autocatalytic model in the study of curing kinetics for DES-based synthesis, indicating that the front propagation speed of the AM-based FP reaction is influenced by the activation temperature, duration of activation temperature, and initial temperature of the solution. Experimental validation of the polyacrylamide manufacturing through frontal polymerization corroborates the accuracy and reliability of the model predictions, providing a reference for the numerical model study and temperature control in the synthesis of acrylamide composite hydrogels via FP.

By installing FBG sensors on the geogrids, smart geogrids can both reinforce and monitor the stability for geogrid-reinforced slopes. In this paper, a geogrid-reinforced sand slope model test was conducted in the laboratory and fiber Bragg grating (FBG) sensing technology is used to measure the strain distribution of the geogrid. Based on the model test, the performance of the reinforced soil slope is simulated by finite element software Midas-GTS, and the stability of the reinforced soil slope is analyzed by strength reduction method. The relationship between the geogrid strain and safety factor is set up. The results indicate that the measured strain and calculated results agree very well. The geogrid strain measured by FBG sensor can be applied to evaluate the stability of geogrid-reinforced sand slopes.

Based on the goals of "high reliability, high frequency, rapid launch, and low cost" for space launch site, and drawing inspiration from the design principles of SpaceX's Starship launch test facility, an integrated dual-sided deflector system for convective cooling and thermal protection is presented. The interaction process between the gas jet and liquid water jet and its effect on the flow field environment are thoroughly studied using numerical calculation methods. Furthermore, Regarding the phase-change heat transfer issue in a compressible gas-liquid two-phase flow and the varying distribution of different bubble shapes and sizes at the gas-liquid interface, a modified Lee model is derived. The rationality and accuracy of the phase-change mass-transfer coefficient values can be improved by adjusting the equivalent diameter and shape factor of the bubbles. The research results demonstrate that compared to the classical Lee model, the modified Lee model can achieve a higher numerical accuracy in predicting the heat and mass transfer processes in gas-liquid two-phase flows. Through comparative analysis with traditional dual-sided deflector and conventional cooling system, the integrated dual-sided deflector system for convective cooling and thermal protection exhibits significant performance advantages in gas flow regulation, flow field environment improvement and erosion protection in the near-ground region of the launch site. Particularly, regarding the thermal environment of the dual-sided deflector, it not only achieves effective flow deflection, but also mitigates the degree of erosion caused by the gas jet on the deflector. This conclusion can provide theoretical references for the thermal protection design of commercial launch vehicle system at space launch site.

To identify the exchange coefficients of a nonlinear polynomial tricompartmental general system of α+β , we follow the procedure; firstly, the recommended solution is to introduce an adequate time t*>0 in a determinable manner. That is, after injecting quantity a into the main compartment, wait a moment for the exchange in the polynomial α+β order nonlinear general system to settle, and then compare this compartment to compartment 2 at t*. Secondly, applying the Taylor formula will linearize the system and identify the exchange coefficients. Finally, we will prove that the linearization method is stable.

Our study essentially concerns the dynamic behavior of an SIRS epidemic model in discrete time. Two equilibrium points are obtained; one is disease-free while the other is endemic. We are interested in the endemic fixed point as well as its asymptotic stability. Depending on the parameters which are estimated using the data from US Department of health and SIRS modelling with optimization, two Flip and Transcritical bifurcations appear. We illustrate their diagrams, as well as their bifurcation curves using the method of Carcasses \cite{carcasses1993determination,carcasses1995singularities} for the Flip bifurcation and by an implicit function deduced from such an equation for the Transcritical bifurcation. We use the scanning of the parametric plane to have a global view of the behavior of the model and to highlight the zones of stability of the existing singularities. A superposition of the bifurcation curves with the parametric plane can show the overlap of the curves with the boundaries of the stability domains, which confirms the smooth running of the simulation and its correspondence with the theory, we finish this article with constrained optimal control applied to infection rate and recruitment rate for an SIRS discrete epidemic model. Pontryagin's maximum principle is used to determine these optimal controls. Finally using COVID-19 data in the USA, we obtain results that demonstrate the effectiveness of the proposed control strategy to mitigate the spread of the pandemic.

This paper introduces a novel multi-strategy enhanced dung beetle optimization (MSDBO) algorithm that is designed to address several issues identified in the standard dung beetle optimization algorithm. Specifically, the MSDBO aims to enhance convergence speed, reduce susceptibility to local optima, and increase search accuracy. By incorporating three strategies: tent chaotic mapping for population initialization, the golden sinusoidal strategy for position updating, and the Lévy flight strategy for balancing exploration and exploitation, the standard dung beetle optimization algorithm is enhanced. The MSDBO algorithm is evaluated using twelve benchmark test functions and compared against five state-of-the-art algorithms. The results consistently show that MSDBO exhibits faster convergence speeds and more accurate solutions than the other algorithms across most of the test functions. In addition, MSDBO is also applied to optimize the parameters of a valve plate, including the close angle, cross angle, triangle groove sizes, and wrap angle. The optimization outcomes reveal that MSDBO effectively minimizes pressure ripples in the piston chamber, resulting in reduced flow rate fluctuations and noise emission compared to the initial design. This study highlights the potential of the MSDBO algorithm in tackling complex nonlinear engineering optimization problems.

This study presents a numerical analysis of rock deformation and damage behaviors induced by dynamic loadings, specifically focusing on the fracturing of tight oil reservoirs. An elastoplastic model based on finite element methods was utilized to simulate the effects of varying loading and perforation parameters. Three distinct scenarios were modeled: a single perforation, boundary loading on three sides, and multiple perforations, each providing unique insights into the mechanical responses of the rock material. The analysis revealed that dynamic effects are highly localized around perforations and loaded boundaries across all scenarios. Acceleration magnitudes are captured with rapid attenuation observed as the distance from the perforation increased. This indicated that the highest stress concentrations and deformation occurred close to the points of loading or perforation. Strain rate analysis in the x direction showed significant fluctuations near the perforations, reflecting intense stress redistributions and compressive deformation, especially prominent in the scenario with multiple perforations. Plastic strain, representing irreversible damage, was found to concentrate near perforation tips and loaded boundaries. However, the pattern of plastic strain accumulation varied: in the small domain scenario with boundary loading, plastic strain did not immediately accumulate next to the loading boundary, whereas in scenarios with single and multiple perforations, damage developed from the loading boundaries. These findings highlight the need for precise control in perforation and dynamic loading processes to manage stress and deformation effectively. It is concluded that optimizing perforation patterns and loading strategies can enhance reservoir stimulation efficiency while minimizing unintended damage to the surrounding formation.

Due to the continued deepening of today’s mining projects, there is a need to develop new or improved rockbolts that are strong and yielding and have high energy absorption when exposed to dynamic loading. Impact tests can provide useful information about the dynamic response of rockbolts but most laboratory tests involve high costs in preparation and result validation. Numerical modeling is an alternative that, in addition to complementing laboratory results, can be used to represent the process of deformation and energy absorption of support elements. In this paper the implementation and results obtained from a finite difference (FDM) numerical model are presented. The model functions as a simulation tool to illuminate all the elements that conform the large-scale (1:1) impact test and their behavior and influence on the dynamic response of a threadbar bolt (22 mm nominal diameter). The model was calibrated using published results and based on these results as well as parametric analysis, the response of each component (steel tube, grout, and bolt) could be identified and its behavior in terms of absorbed energy and displacement could be observed. Results show the model can provide important preliminary information to make design decisions about support elements design.

Hemodynamic relevance of differently located stenoses of the internal jugular veins remains undetermined. It particularly concerns nozzle-like strictures in the upper parts of these veins and stenotic jugular valves located at the end of this vein. This study was aimed at understanding flow disturbances caused by such stenoses. Computational fluid dynamics software, Flowsquare+, was used. We constructed 3-dimensional models of venous outflow, comprising two alternative routes: tube representing the internal jugular vein and irregular network representing the vertebral veins. At the beginning of the tube representing the internal jugular vein, differently-shaped and sized short strictures, representing nozzle-like strictures were built in. At the end of this tube, differently-shaped membranes, representing the jugular valve, were built in. With the use of computational fluid dynamics modeling, we studied how these two obstacles influenced the outflow. We found that the most relevant outflow disturbances were evoked by nozzle-like strictures in the upper part of the internal jugular vein that were either small, long or asymmetrically positioned. Very tight stenotic valves and septum-like malformed valve were equally hemodynamically relevant. These findings suggest that both upper and lower strictures of the internal jugular vein can be of clinical significance.

Recently, there has been a growing interest in the production of electricity from renewable energy sources (RES). The RES investment is characterized by uncertainty, which is long-term, costly, depend on feed-in-tariff and support schemes. In this paper, we address the real option valuation (ROV) of a solar power plant investment. The real option framework is investigated. This framework considers the renewable certificate price, furthermore the cost of delay between establishing and operating the solar power plant. The optimal time of launching the project and assess the value of deferred option are discussed. The new three stage numerical methods are constructed, the Lobatto3C-Milstein (L3CM) methods. The numerical methods are integrated with concept of Black-Scholes option pricing theory, and applied in option valuation for solar energy investment with uncertainty. The numerical results of L3CM, finite difference and Monte Carlo methods are compared to show the efficiency of our methods. Our data set refers to the Arab Republic of Egypt.

The control of the fine particle emission from the coal-fired power plants has been a big challenge in china. Either the conventional electrostatic precipitator or the fabric filter has some shortages. The electrostatic precipitator is hard to keep high efficiency for the submicron particles. The fabric filter could keep high efficiency but with very high pressure drop. One solution is to combine them together. Several types of hybrid particle collectors (HPC) combining the electrostatic precipitator and the fabric filter together have been developed and applied in a number of projects in China. In recent years, numerical simulation technology has been widely used in the design of the collectors. In this paper, an unsteady numerical model is used to investigate the particle deposition characteristics in a scale hybrid collector. The fluid field is modeled by using incompressible Navier-Stokes equations with the RNG turbulence equations. The corona discharge is solved by using a finite volume method. The particle charging is modeled by using a filed-diffusing combined model. And an unsteady cake formation model is used to simulate the filtration process.

The electric vehicle (EV) industry challenges battery joining technologies by requiring higher energy density both by mass and volume. Improving the energy density via new battery chemistry would be the holy grail but is seriously hindered and progresses slowly. In the meantime, alternative ways, such as implementing more efficient cell packaging by minimizing the electrical resistance of joints are of primary focus. In this paper we discuss the challenges associated with the electrical characterisation of laser soldered joints in general, and the minimization of their resistive losses, in particular. In order to assess the impact of the joint resistance on the overall resistance of the sample, the alteration in resistance was monitored as a function of voltage probe distance and modelled by finite element simulation. The experimental measurements showed two different regimes: one far from of the joint area and another in its vicinity and within the joints cross section. The presented results confirm the importance of the thickness of the filler material, the effective and the total soldered area, and the area and the position of the voids within the total soldered area in determining the electrical resistance of joints.

The water exit problem, which is ubiquitous in ocean engineering, is one of the significant research topics in the interaction between navigators and water. The study of the water exit problem can help to improve the structure design for marine ships and underwater weapons for better strength and movement status. However, the water exit problem involves complex processes such as gas-liquid-solid three-phase coupling, cavitation, water separation, liquid surface deformation and fragmentation, making it challenging to study. According to the works carried out by many researchers towards this issue, we summarize the recent developments from three aspects, i.e. theoretical research, numerical simulation and experimental results. In theoretical research, the improved von Kaman model and linearized water exit model are introduced. Several classical experimental devices, data acquisition means and cavitation are introduced in experimental development. Three methods, namely, BEM(boundary element method), FVM(finite volume method) and LES(large eddy simulation) are concluded for numerical simulation. The limitations and shortcomings of these three aspects are then analyzed. Finally, an outlook on the future research improvements and developments of the water exit problem is provided.

The majority of existing numerical simulations of the effect of bridge piers on water movement are based on a limited number of bridge piers at a laboratory scale. Furthermore, some 2D numerical simulations for actual bridge projects have deficiencies, including the use of overly large mesh sizes and an inadequate treatment of bridge piers. In this study, we compare three methods of structural mesh encryption, suspension mesh, and nonstructural mesh based on Delft3D, and apply the optimization scheme to a real bridge project. It is demonstrated that optimal results can be achieved by utilising a grid size comparable to the pier diameter (Dp) in the region away from the piers. In the vicinity of the piers, the grid cell size should be no larger than 1/9 Dp. The suspended grid technique (DD Boundary) can yield results consistent with those obtained using a full-area high-resolution grid, provided that the total number of grids can be reduced and the computational time is considerably reduced. In this study, the unstructured mesh (Deflt3D Flexible Mesh) scheme was unable to capture the oscillations of the wake flow behind the bridge piers. However, the application of the optimized scheme in bridge engineering has demonstrated its practical value. The findings of this study on mesh resolution and suspension mesh scheme can be applied to the Delft3D software and are also useful for other numerical simulation work.

This paper focuses on the combination of the method of moments and the wire-grid approximation as an effective computational technique for modeling symmetrical antennas with low computational cost and quite accurate results. The criteria and conditions for the use of wire-grid surface approximation from various sources are presented together with new recommendations for modeling symmetrical antenna structures using the wire-grid approximation. These recommendations are used to calculate the characteristics of biconical and horn antennas at different frequencies. The results obtained using different grid and mesh settings are compared to those obtained analytically. Moreover, the results are compared to those obtained using the finite difference time domain numerical method, as well as the measured ones. All results are shown to be in a good agreement. The used recommendations for building a symmetrical wire-grid of those symmetrical antenna elements provided the most advantageous parameters of the grid and mesh settings and the wire radius, which are able to give a quite accurate results with low computational cost. Additionally, the known equal area rule was modified for a rectangular grid form. The obtained radiation patterns of a conductive plate using both the original rule and the modified one are compared with the electrodynamic analysis results. It is shown that the use of the modified rule is more accurate when using a rectangle grid form.

COVID-19 infection is extremely dangerous for persons with comorbidities such as kidney disease and diabetes. Although many people of all ages and socioeconomic backgrounds have suffered due to this pandemic, those with underlying medical conditions like kidney disease and diabetes are more vulnerable to becoming infected with COVID-19 related to weak immunity. However, the role of vaccination in the co-infection of COVID-19 with diabetics is not available in the literature. Therefore, we examine the unique challenges presented by the co-infection of COVID-19 in individuals with diabetes, focusing on its disease transmission dynamics. A mathematical modelling technique is to be employed in this work to analyse the dynamics of COVID-19 outbreaks with seven compartments, including vaccination with comorbidities such as diabetes. We also investigated analytically by displaying the solutions' existence, boundedness, positivity and sensitivity. After calculating the basic reproduction numbers, the stability analysis of model equilibrium points is also covered. If the reproduction numbers are less than unity, the disease-free equilibrium points of the model appear to be both locally and globally stable. The findings of this study are that as the vaccination rate rises, the incidence of COVID-19 and its co-infections with diabetes decreases. Effective disease treatment strategies may be based partly on the possible impact of vaccination on the co-infection of COVID-19 related to diabetic disease insights.

In this paper we report a new development on the numerical model for aluminum-steam combustion. This model is based on diffusion flame of continuum regime and the thermal equilibrium between the particle and the flow field, which can be used to calculate the aluminum particle combustion model for two phase calculation conditions. The model prediction is in good agreement with the experimental data. A new type of vortex combustor was proposed for the combustion of aluminum and steam, and the mathematical model of the two phase reacting flow with in this combustor was established. The turbulence effects are modeled using the Reynolds Stress Model (RSM) with Linear Pressure-Strain approach, and the Eddy-Dissipation model is used to simulate the gas phase combustion. Aluminum particles are injected into the vortex combustor and form a swirling flow around the chamber and their trajectories are traced using the Discrete Phase Model (DPM). The simulation results show that the vortex combustor can achieve high efficient combustion of aluminum and steam. The influencing factors, such as the eccentric distance of the inlet of aluminum particles, particle size and steam inlet diameter, etc., are studied. The work described in this paper represents an attempt to the design of a vortex combustor in order to increase aluminum combustion efficiency.

The article presents investigations of rail vehicle bogie of Y25 type. The Y25 bogie family is one of the most commonly used freight car bogie designs. In addition to several significant advantages characterising this design, it also has several disadvantages observed since the beginning of more than fifty years of operation in several types of cargo vehicles. One of these defects observed in real systems is "unsatisfactory running stability", particularly at long straight tracks. The article used commercial engineering software VI-Rail to create a model of the gondola car (type 412W Eaos) with two Y25 bogies. The car model was tested empty and loaded (maximum permissible load). Motion along straight and curved tracks with different radii values was analysed. The vehicle velocity was changed from a few m/s to the maximum values for which stable solutions of the model existed. For each route, a non-linear critical velocity was determined, defining the maximum operating velocity of the modelled car. The model solutions were recorded, while just one was selected to present the results - the first wheelset lateral displacement ylw. The opinion about the "imperfect running quality" in curved tracks was confirmed. The possible appearance of self-exciting wheelset vibrations in the modelled car's operating velocity range in the laden state was also observed. The research results on the impact of the bogie suspension parameters change on the vehicle model stability are presented. The crucial parameter in the bogie suspension was indicated. Reducing its value by several percent about the nominal value increases the critical velocity of the car to values higher than the maximum operating velocity of the modelled vehicle.

The objective of this study, which is divided into two parts, is twofold: to address long standing challenges in the sensing of atmospheric turbulence-induced wavefront aberrations in strong scintillation conditions via comparative analysis of several basic scintillation resistant wavefront sensing (SR-WFS) architectures and iterative phase retrieval (IPR) techniques (Part I, this paper), and to develop a framework for the potential integration of SR-WFS techniques into practical closed-loop non-astronomical atmospheric adaptive optics (AO) systems (Part II). In this paper we consider basic SR-WFS mathematical models and phase retrieval algorithms, tradeoffs in sensor design and phase retrieval technique implementation, and methodologies for WFS parameter optimization and performance assessment. The analysis is based on wave-optics numerical simulations imitating realistic turbulence-induced phase aberrations and intensity scintillations, and optical field propagation inside SR-WFSs. Several potential issues important for practical implementation of SR-WFS and IPR techniques such as requirements for phase retrieval computational grid resolution, tolerance with respect to optical element misalignments, and the impact of camera noise and input light non-monochromaticity are also considered. Results demonstrate that major wavefront sensing requirements desirable for AO operation in strong intensity scintillations can potentially be achieved by transitioning to novel SR-WFS architectures based on iterative phase retrieval techniques.

The material properties and structural characteristics of ballistic composites are crucial to their ballistic performance. Investigation for the influence of the material and structural parameters of UHMWPE body armor on anti-fragmentation penetration performance has been made by numerical simulation in this paper, with a numerical modeling simulating the penetration of a 1.1 g FSP into the target plate of UHMWPE body armor, by which the simulation results demonstrate that the failure process of the body armor target plate primarily involves shear damage, interlayer delamination, and tensile damage. Furthermore, the study evaluates the effects of elastic modulus, tensile strength, shear strength, number of layers, and interlayer strength on the ballistic limit velocity of UHMWPE body armor. The findings reveal that the ballistic limit velocity is most sensitive to changes in shear strength, with variation rates ranging from -18.02% to 10.81%, showing an approximate positive correlation, while the elastic modulus has the smallest impact on ballistic limit velocity, with variation rates ranging from -1.8% to 3.6%. Additionally, appropriate interlayer strength can improve the ballistic limit velocity of the body armor to a certain extent. This study provides theoretical methods and recommendations for optimizing anti-penetration performance of UHMWPE body armor.

Promoting the use of solar energy resources has always faced the challenges of instability and supply-demand mismatch. The key to solving these issues is to efficiently store and utilize solar energy resources using high-performance heat storage devices. This study designed a high-performance shell-and-tube phase change thermal storage device and established a numerical model using ANSYS software to summarize the device's dynamic melting law. To verify the accuracy of the numerical simulation, a performance testing platform for the phase change thermal storage device was built to investigate the impact of factors such as inlet water temperature, inlet water flow rate, type of heat storage, and initial temperature of the device, and to reveal the change law of the device's performance. The results show that the inlet water temperature has the most significant impact on the device's heat storage and release performance. When the device's heat storage or release is used for heating, changing the inlet water flow rate has a weak and limited effect on the device's performance. However, when the device's heat release is used to provide domestic hot water, increasing the makeup water temperature and reducing the inlet water flow rate can significantly improve the device's effective heat release. Furthermore, based on the experimental validation of the model's correctness, this study further simulated and studied the impact of different factors on the device's heat storage process to optimize its structural design and provide technical references for the device's actual operation and installation. The results show that the placement of fins has a negligible effect on the performance of the heat storage device while reducing the fin spacing and increasing the fin thickness can significantly improve the melting efficiency of the phase change material (PCM). Additionally, the heat storage characteristics of the device are significantly better in the vertical installation mode than in the horizontal installation mode. This study provides theoretical guidance and technical references for the design and use of phase change thermal storage devices.

A recently published article posited the use of polymer matrix composites (PMCs) as a replacement material for low-carbon steel (LCS) in large power transformer (LPT) tanks [1]. PMCs offer potential advantages such as weight reduction, corrosion prevention, and self-healing of the tank material which may not be achievable with traditional LCS. LCS can arrest ballistic projectiles at relatively low velocities, but at higher velocities, penetration can occur leading to a potential major transformer failure. Numerical simulations of steel projectile impacts upon LCS plates and glass and/or carbon fiber polymer composite laminates showed that with sufficient thickness both could prevent projectile penetration. This prevention was especially efficient when the plates were coated with a high strain rate sensitive material such as polyurea (PU). Our simulations show that the addition of a 3 mm coating of PU to a standard 10 mm thick LCS plate used in LPT tanks could prevent penetration of a 400 m/s steel projectile. For PMC plates, the simulations with the same projectile velocity showed three times faster deceleration with a 3 mm PU coating than without the coating. Thus, two realistic approaches to improving the ballistic protection of LPTs are presented. In addition to their other advantages, this research clearly shows that PMCs with the coating could provide significant ballistic protection for the next generation of LPT tanks.

Wildland fires are responsible for large socio-economic impacts. Fires affect the environment, damage structures, threaten lives, cause health issues, and involve large suppression costs. These impacts can be mitigated via accurate fire spread forecast to inform the incident management team. We show that a fire forecast system based on a numerical weather prediction (NWP) model coupled with a wildland fire behavior model can provide this forecast. This is illustrated with the Chimney Tops II wildland fire responsible for large socio-economic impacts. The system is run at high horizontal resolution (111 m) over the region affected by the fire to provide a fine representation of the terrain and fuel heterogeneities and explicitly resolve atmospheric turbulence. Our findings suggest that one can use the high spatial resolution winds, fire spread and smoke forecast to minimize the adverse impacts of wildland fires.

This study presents the complexity and sensitivity of chaotic system dynamics in the case of the double pendulum. It applied detailed numerical analyses of the double pendulum in multiple computing platforms in order to demonstrate the complexity in behavior of the system of double pendulums. The equations of motion were derived from the Euler-Lagrange formalism, in order to capture the system's dynamics, which is coupled nonlinearly. These were solved numerically using the efficient Runge-Kutta-Fehlberg method, implemented in Python, R, GNU Octave, and Julia, while runtimes and memory usage were extensively benchmarked across these environments. Time series analyses, including the calculation of Shannon entropy and the Kolmogorov-Smirnov test, quantified the system's unpredictability and sensitivity to infinitesimal perturbations of the initial conditions. Phase space diagrams illustrated the intricate trajectories and strange attractors, as further confirmation of the chaotic nature of the double pendulum. All the findings have a clear indication of the importance of accurate measurements of the initial condition in a chaotic system, contributing to an increased understanding of nonlinear dynamics. Future research directions are faster simulations using Numba and GPU computing, stochastic effects, chaotic synchronization, and applications in climate modeling. This work will be useful for understanding chaos theory and efficient computational approaches in complex systems of dynamical nature.

In the recent years, tubular nanostructures have been related with immense advances in various fields of science and technology. Considerable research efforts have been centered on the theoretical prediction and manufacturing of non-carbon nanotubes (NTs), which meet modern requirements for the developing of novel devices and systems. In this context, diatomic inorganic nanotubes formed by atoms of elements from the 13th group of the periodic table (B, Al, Ga, In, Tl) and nitrogen (N), have received much research attention. In the present work, the elastic properties of single-walled boron nitride, aluminium nitride, gallium nitride, indium nitride and thallium nitride nanotubes were numerically evaluated using the nanoscale continuum modelling approach (also called molecular structural mechanics). The elastic properties (rigidities, surface Young´s and shear moduli, and Poisson´s ratio) of nitride nanotubes are discussed with respect to the bond length of the corresponding diatomic hexagonal lattice. The results obtained contribute to a better understanding of the mechanical response of nitride compound-based nanotubes, covering a broad range from the well-studied boron nitride NTs to the hypothetical thallium nitride NTs.

Local polynomial smoother (LPS) is a weighted local least-squares nonparametric method. It provides a local Taylor series fit of the data at any location and can be directly used in a differential equation to provide a numerical scheme. In this article, we introduce this new nonparametric idea based on local polynomial smoother, for acquiring the numerical solution of the Bagley-Torvik fractional-order differential equations. Furthermore, this paper will present a numerical comparison with some methods, such as legendre operational matrix and pseudo-spectral method. The efficiency and accuracy of the LPS method are demonstrated by two numerical examples.

Predictive numerical models in the study of ground-borne vibrations generated by railway sys-tems have traditionally relied on the segmented approach of subsystems. In this approach, loads are individually applied, and the cumulative effect of rolling stock is obtained through superposi-tion. While this method serves to mitigate computational costs, it may not fully capture the com-plex interactions involved in ground-borne vibrations. Recent advancements in computational and software tools have enabled the development of more sophisticated vibrational pollution models. These advanced models encompass the entire dynamic system, from the rolling stock to the terrain, allowing for continuous simulations with a defined time step. Furthermore, the incorporation of computational contact tools between various ele-ments not only ensures temporal accuracy but also extends the analysis to the frequency domain. This paper explores the evolution of predictive numerical models in the context of ground-borne vibrations generated by trains on ballasted tracks. It highlights the shift from segmented models to continuous simulations, demonstrating the progressive advancement towards more comprehen-sive and accurate representations of railway-induced and generated vibrations in the track sur-roundings as an additional way of evaluate sustainability of infrastructures. In particular, this development of a numerical model from its inception aimed to expand upon ex-isting emission - prediction models by enriching the data with frequency-domain information. This approach has been driven by the aspiration to bridge the gap between numerical predictions and measured signals, striving for a higher level of accuracy in our obtained results. This model can predict the impact of a high-speed rail (HSR) vehicle passing on its own track. It possesses sufficient information in both the time and frequency domains to precisely characterize vibrations. This initiative marks a pivotal advancement, enabling us to more effectively capture the intricacies of ground-borne vibrations and their impact on the surrounding environment due to a deeper comprehension of the occurrences in the frequency domain, there's a heightened ca-pability to accurately characterize vibrations.

The geological and hydrological conditions of tunnels near reservoirs are complex, and the impact of water level changes on the stability and leakage assessment of the lining structure is not considered in the current leakage evaluation. In order to construct an evaluation model for leakage level of tunnels near reservoirs, the influence of water level changes on tunnel stability and changes in environmental conditions on the leakage of tunnels were researched. Based on the AHP and extensibility theory, a hierarchical system for leakage assessment, incorporating values from nine indexes representing three aspects: geological conditions, hydrological conditions, and tunnel engineering. Numerical simulation was used to analyze the influence of water level changes. It was found that the water level change index greatly influences the displacement and stress distribution inside the tunnel structure. The leakage evaluation model was applied to the Tiebeishan Tunnel, resulting in a rating of 3, indicating medium-level leakage. Attention should be paid to water leakage in tunnels with changes in reservoir water levels. The leakage evaluation model for tunnels near reservoirs can effectively assess leakage levels under various conditions, providing a reference for safety assessments of tunnel leakage near reservoirs.

Cascading landslides represent a dynamic and hazardous geological phenomenon, char-15 acterized by a sequential chain of slope failures triggered by various factors such as heavy rainfall, 16 seismic activity, or anthropogenic activities. These events can amplify the damage caused by the 17 initial trigger and propagate instability along a slope, often resulting in significant environmental 18 and societal impacts. The Morino-Rendinara cascading landslide, situated in the Ernici mountains 19 along the border of Abruzzo and Lazio region (Italy), serves as a notable example of the complexities 20 and devastating consequences associated with such events. In March 2021, a substantial debris flow 21 event in Morino obstructed the Liri River, marking the latest step in a series of landslide events 22 characterized by a complexity far beyond initial expectations. This study employs a multidiscipli-23 nary approach, combining conventional techniques and advanced technologies, to understand the 24 complexities of the Morino-Rendinara landslide. Field activities, UAV image acquisition, SAR inter-25 ferometry based on SENTINEL-1 images and pixel offset analysis based on high-resolution Google 26 Earth images, offer insights into the geological and hydrogeological setting of the unstable slope, 27 landslide geometry, mechanism and kinematics. To supplement the analysis a specific FEM slope 28 stability analyses is used to reconstruct the deep geometry of the system emphasizing the modula-29 tion action of groundwater flow to the slope stability

Molten salt electrolysis stands as an indispensable technique for the production of rare earth metals, where the judicious distribution of the electric field within the electrolytic cell is paramount for ensuring the seamless operation of the electrolysis process. This study introduces a pioneering approach by integrating the Nernst-Planck equation with the Butler-Volmer equation, thereby devising a holographic electric field model tailored for the electrolysis of the NdF3-LiF molten salt system. The paper presents pioneering experimental measurements concerning the potential distribution and voltage drop characteristics within a 6kA rare earth electrolytic cell, elucidating the interplay between decomposition voltage and current density, which substantiates the model's precision. Our findings unveil that, under operational conditions of 4010 amperes, the 6kA electrolytic cell incurs a total voltage drop of 6.3 volts, delineating a melt voltage drop of 4.43 volts, an ohmic voltage drop of 0.248 volts, and an interfacial voltage drop of 1.622 volts at the electrode-electrolyte junction. The discrepancy in error is confined to within 0.3% of the actual measured voltage of 6.28 volts.The simulation outcomes from the holographic electric field model have been confirmed to precisely and comprehensively reflect the characteristics of the electric field distribution during the rare earth molten salt electrolysis process, offering valuable theoretical insights for practical manufacturing applications.

The research on the excavation mechanical properties of underpass tunnels has already had certain results, but fewer of them consider the effects of dynamic and static loads on the excavation mechanical properties of underground tunnels at the same time, especially the lack of research cases of double-line highway tunnels with angled underpasses of existing railway tunnels. In this paper, based on the tunnel project of the new double-line Shiqian Highway Tunnel passing under the Hurong Railway at an oblique angle, the mechanical parameters of the rock body in the cross-tunneling area under the influence of static and dynamic load coupling are obtained by using the disturbance quantification coefficients and the Hoek-Brown criterion based on the method of over-advanced geological prediction and the investigation of the surrounding rock on the palisade face, and the engineering-scale finite element model is constructed and the mechanical characteristics of the excavation in the tunnel of the double-lane highway tunnel passing under the existing railroad under the influence of the static and dynamic loads are analyzed. The results show that: the new tunnel rock movement law for the top of the arch sinking, bottom plate bulging, side wall to move outward, the height and width of the arch and bottom plate arch with the tunnel excavation and increase, the side wall rock displacement effect is smaller; left and right lines of the tunnel disturbed area of the rule of change is similar; the existing tunnel bottom plate displacement is more than the top plate and the left and right side walls by the excavation of the time step influence is larger. Typical profile displacement is mainly determined by the distance from the excavation surface; the new tunnel in the four types of calculation conditions of the rock disturbance angle is the same; the existing tunnel roof and side wall Von Mises stress extreme occurs in the tunnel structure, the bottom plate there are pressure relief and pressure-bearing areas; based on the displacement control criteria, the recommended upper and lower step blasting method.

This review is devoted to present the history of CARDIOSIM© software simulator platform, which was developed in Italy to simulate the human cardiovascular and respiratory system. The first version of CARDIOSIM© was developed at the Institute of Biomedical Technologies of the National Research Council in Rome. The first platform version published in 1991 ran on PC with disk operating system (MS-DOS) and was developed using the Turbo Basic language. The last version runs on PC with Microsoft Windows 10 operating system; it is implemented in Visual Basic and C++ languages. The platform has a modular structure consisting of seven different general sections, which can be assembled to reproduce different pathophysiological conditions. The software simulator can reproduce the most important circulatory phenomena in terms of pressure and volume relationships. It represents the whole circulation using a lumped-parameter model and enables the simulation of different cardiovascular conditions according to Starling’s law of the heart and a modified time-varying elastance model. Different mechanical ventilatory and circulatory devices have been implemented in the platform including thoracic artificial lung, ECMO, IABP, pulsatile and continuous right and left ventricular assist devices, biventricular pacemaker and biventricular assist devices. CARDIOSIM© is used in clinical and educational environment.

The bentonite slurry trenches are becoming increasingly popular in the excavation process of trenches, especially for diaphragm wall construction. The problem that needs to be addressed is the stability of bentonite slurry trenches. This paper presents a stability analysis of trenches with temporary support of bentonite slurry, with unit weights ranging from 10.5 to 12.0 kN/m3, in the realistic stratum: C3 in the Hue city area. The analysis employs the shear strength reduction technique (SSR) with Mohr-Coulomb materials for the numerical method to evaluate the factor of safety (FS). The finite element method (FEM) software program (RS) and the finite difference method (FDM) FLAC software were applied. Additionally, the limit equilibrium method (LEM) of Bell-Washbourne and three-dimensional (3D) Bishop are used for calculating trench stability. The analysis results show a good agreement between RS2 and FLAC2D and RS3 and FLAC3D. Secondly, upon comparison, it is noted that the safety factors of the 3D software programs (RS3, FLAC3D) are higher than those of the 2D software programs (RS2, FLAC2D), ranging from 52.3 to 63.0% for trench lengths of 6m. However, for trench lengths of 54m, the safety factor values in 2D and 3D are nearly equal. Thirdly, the safety factor of the Bell-Washbourne method (LEM) is lower than the numerical analysis methods (FLAC and RS programs). Using the three-dimensional numerical method appears to be effective for estimating stability.

This study concerns the transport and retention of polydisperse micron-sized (16 ± 6 µm) of microplastics (MPls) in porous media under varying flow rate conditions. Sorption kinetics were modeled using first-order reversible and irreversible kinetic sorption models, with sensitivity analyses providing insights into each sorption parameter's effect. Both numerical modeling and experimental measurements were employed to assess sand filter retention rates. The impact of flow rate on sorption reveals variations in distribution coefficient (Kd), mass transfer coefficient (β), and irreversible sorption rate (K1). Lower flow rates are correlated with higher Kd and β values, indicating an increase in sorption and diminished mass transfer rates. The findings revealed that an increase in Kd resulted in a more gradual sorption process with a decrease in peak concentration, whereas changes in β influenced the rate of sorption and peak concentration to a lesser extent compared to Kd. Lower K1 values are associated with higher peak concentrations and decreased retention efficiency. Retention rates were evaluated by a numerical model and found as 28 ± 1% at a flow rate of 31 ml min⁻¹ and 17 ± 1% at 65 ml min⁻¹. The introduction of MPls into soil environments has been noted to modify transport dynamics into soil. As a result, these alterations effects hydrological characteristics of soil, thereby impacting quality of groundwater and agricultural output. The mean absolute error (MAE) of 6% between modeled and observed retention rates suggests minor discrepancies. This study highlights the importance of examining retention efficiency and the accuracy of numerical models in porous media during MPl transport.

Monitoring and protecting flowing watercourses is a complex and challenging task that requires the collaboration and coordination of various stakeholders such as governments, industries, farmers, consumers and environmental groups. The study of the dynamics of the concentration of polluting factors and especially the concentrations of heavy metals and highlighting a seasonal variation is a necessary element from this point of view. In this article, we present the results of our analyzes carried out in two measurement campaigns executed in 10 monitoring points along the Danube River, between Braila city and Isaccea city in the pre-deltaic area, during the summer season and autumn season 2022. The importance of this area is given by the fact that the Dan-ube Delta is part of the UNESCO heritage and the monitoring of polluting factors is a necessity in the desire to protect this area. The data measured during the July and August 2022 campaign co-vers a wide range of chemical species: Phosphate, CCO, CBO5, NH4+, N-NO2, N-NO3-, N-Total, P-PO4 3-, SO42-, Cl-, phenols, as well as metals with a harmful effect: Aluminum, Arsenic, Cad-mium, Chromium, Iron. The study includes an evaluation based on the statistical approach of the results to highlight the significant correlations and differences identified between the two data sets. Next, to highlight the obtained results, a numerical model was considered using HEC-RAS and ESRI ArcGis applications in a two-dimensional unsteady flow model in order to get the non-homogenous concentrations’distributions in the studied area. These two-dimensional models have been less studied in the specialized literature. In this way, interesting results could be ob-tained and prediction methods regarding the dynamics of metal concentrations could be struc-tured. The data obtained was used for the terrain model from the USGS service and the flows of the Danube and its two tributaries were simulated using the data provided by the national ser-vices. In this work, we present the results obtained for the dynamics of the concentrations of the metals Aluminum, Arsenic, Cadmium, Chromium and Iron and the evaluation of the specific ab-sorption coefficients for the explanation and correlation with the results of the measurements. The study is part of a more complex project in which the spatio-temporal distribution of the polluting factors in the water was evaluated and the habitats in the study area were inventoried - especially those of community interest. the aspects related to the influence of the distribution of polluting factors on the state of health will be the subject of another article.

A comprehensive modeling framework for the thermoforming of polymer matrix woven laminate composite was developed. Two numerical indicators, the slip path length and traction magnitude, have been identified to be positively correlated to matrix smearing and wrinkling defects. The material model has been calibrated with picture-frame experimental results, and the prediction accuracy for intra-ply shear and thickness distribution were examined with measurements on the physically formed parts. Finally, a parametric study was conducted to determine the relationship between various process parameters and the quality of the formed part. For the trapezoidal part, orienting the laminate at 45-degrees to the mold axis reduces the likelihood of matrix smear and wrinkling defects. Although this laminate orientation yielded a greater spatial variation of part thickness, the thickness deviation is lower than that for the 0-degree orientation case. Higher forming rates were found to result in more defects due to increased slip path lengths and traction magnitudes. Optimal forming speeds should balance rapid shaping within the operating temperature window while avoiding tear. It was discovered that shallower molds benefit from faster ramp rates, while deeper molds require slower rates to manage extensive shearing, stretching, and bending. Faster forming rates lead to smaller thickness increases at high intra-ply shear regions, indicating a shift from intra-ply shear to out-of-plane bending due to the visco-plastic effect of the molten laminate, and can negatively impact part quality. Lastly, it was shown that a well-conceived strategy using darts could improve the part quality by reducing the magnitude of the defect indicators.

This article develops duality principles and numerical results for a large class of non-convex variational models. The main results are based on fundamental tools of convex analysis, duality theory and calculus of variations. More specifically the approach is established for a class of non-convex functionals similar as those found in some models in phase transition. Finally, in some sections we present concerning numerical examples and the respective softwares.

When limited area models of the hydrodynamics of the atmosphere and of the ocean are 1 run over an extended time, variability unrelated to external "drivers" emerge - this variability is 2 colloquially named "hydrodynamical noise" or just "noise". This article summarises what we have 3 learned in the past few years about the properties of such noise, and the implications for numerical 4 experimentation and analysis. The presence of this noise can be identified easily in ensembles of 5 numerical simulations, and it turns out that the intensity of the noise is closely linked to the scale- 6 dependent "memory". At the system level, this "memory" term as given by Hasselmann’s Stochastic 7 Climate model plays a key role. In the case of marginal seas, the process of baroclinic instability 8 modulated by tides and the formation of seasonal thermoclines are significant aspects. Some more 9 general aspects are discussed, such as the applicability of the Stochastic Climate model to systems 10 outside of atmospheric and oceanic dynamics, the irreversibility of tipping points and the challenges 11 of detecting changes beyond a noise level and of the attribution of causes of change.

Climate warming is expected to have pronounced effects on Arctic and Subarctic ecosystems, especially regions underlain by discontinuous and relatively warm permafrost. The main goal of this research is to evaluate the vulnerability and dynamics of permafrost under climate warming across the various ecotypes in respect of ecosystem stability, socioeconomic impact, and for better understanding possible future environmental changes. We suggested the new version of the spatially distributed permafrost dynamics model (GIPL2-MPI), which is developed in the Geophysical Institute, University of Alaska Fairbanks. This model is based on the ecosystem approach to simulate the permafrost dynamics, which we are discussing in this paper. We combined ground-based observations and numerical freeze/thaw modeling using climate-ecosystem-permafrost interactions to understand the physical processes and mechanisms controlling permafrost physical state. We predict the changes in permafrost conditions using output from two GCMs (NCAR-CCSM4 and GFDL-CM3) and Five-Model Average Ensemble for the RCP-4.5 and RCP-8.5 scaled down to 1 by 1 km spatial resolution (https://uaf-snap.org/) across entire Alaska. Our result shows that by the end of the current century widespread near-surface permafrost degradation could begin everywhere in Alaska southward of the Brooks Range as well as across some spots at the North Slope Alaska.

In this study, a Traveling-Wave Mach-Zehnder intensity Modulator (TW-MZM) was designed and optimized for six different electro-optical (EO) crystals: Lithium Niobate (LNB), Potassium Niobate (KNB), Lithium Titanate (LTO), Beta Barium Borate (BBO), Cadmium Telluride (CdTe), and Indium phosphide (InP). The performance of each EO crystal, including optical and radio frequency (RF) loss, applied voltage, and modulation bandwidth, was estimated and compared. The results suggested that, in theory, KNB, LTO, BBO, and CdTe have the potential to outperform LNB. However, it should be noted that the loss associated with KNB and LTO is comparable to that of LNB. The findings demonstrated that BBO and CdTe exhibit a modulation bandwidth exceeding 100 GHz and demonstrate the lowest loss among the considered crystals based on the assumed geometry.

The mechanical properties of the materials stored inside agricultural silos influence the loads generated inside them. The different silo calculation methodologies allow predicting these loads under static and dynamic conditions if the mechanical properties of the stored materials are known. It has been known for many years that these properties are highly dependent on the moisture content of the material. However, to date there are not many studies focused on its determination. The goal of this research is the determination of the angle of internal friction, apparent cohesion, dilatancy angle and apparent specific weight of maize subjected to different moisture contents. Direct shear and oedometer assay devices were used for this purpose. In addition, the maize samples used were moistened using a climatic chamber. Moisture contents applied to maize sample ranged from 9.3% to 17.4%. The results obtained for the angle of internal friction and apparent cohesion, as well as those obtained for the apparent specific weight, were similar to those reported in the literature. However, the values obtained for the dilatancy angle of maize as a function of moisture content could not be compared because nothing has been published so far. The values obtained for this parameter overlap with those published for this material under ambient conditions. In addition, no direct relationship between the values obtained for this parameter and the moisture content of the sample could be found. Finally, the apparent specific weight increased with increasing moisture content of the samples, which differed from that published in the literature. The values provided here can be used for the optimization of storage and handling structures for granular agricultural materials using numerical methods.

Tegucigalpa, the capital city of Honduras, has the highest number of landslides recorded in the country. In a previous study, landslides were found to be mainly concentrated in areas with colluvium and residual soils. The present study focused on assessing the physical vulnerability of houses and the risk linked to landslides in these types of soils. As an input for this assessment, the study generated landslide risk vulnerability functions based on empirical data. The study's findings suggest that the estimated Annual Average Losses associated with landslides in this region is USD 7.26 million. Moreover, the study developed a novel approach for the spatial and temporal estimation of landslide hazard probability, which overcomes some of the limitations identified in previous landslide risk assessment studies.

The development of computational acoustics allows simulation of sound generation and propagation in complex environment. In particular, meshfree methods are widely used to solve acoustics problems through arbitrarily distributed field points and approximation smoothness flexibility. As a Lagrangian meshfree method, smoothed particle hydrodynamics (SPH) method reduce the difficulty in solving problems with deformable boundaries, complex topologies, or multiphase medium. The traditional SPH method has been applied in acoustic simulation. This study presents the corrective smoothed particle method (CSPM), which is a combination of SPH kernel estimate and Taylor series expansion. The CSPM is introduced as a Lagrangian approach to improve accuracy in solving acoustic wave equations in the time domain. Moreover, a boundary treatment technique based on the hybrid meshfree and finite difference time domain (FDTD) method is proposed to represent different acoustic boundaries with particles. To model sound propagation in pipes with different boundaries, soft, rigid, and absorbing boundary conditions are built with this technique. Numerical results show that the CSPM algorithm is consistent and demonstrates convergence with exact solutions. Main computational parameters are discussed, and different boundary conditions are validated to be effective for benchmark problems in computational acoustics.

The numerical simulation of entropy generation at mixed natural free and forced convection flow in a square lid-driven cavity filled with an eletrically conducting binary fluids saturated porous media in the presence of a magnetic field is performed in this investigation. Both the top and bottom horizontal walls of the cavity are kept at constant and different temperatures while the left and right vertical walls are adiabatic and insulated. The magnetic field is applied in normal direction to the cavity. The Darcy model, including Brinkman term relative to viscous effects and Forchheimer term due to inertial forces, is used for the momentum equations, and the SIMPLER algorithm, based on the control volume finite-element method approach is used to solve the pressure –velocity coupling. The flow pattern and the heat transfer characteristics inside the cavity are presented in the form of isotherms, streamlines and isentropic lines and average Nusselt number trend. The entropy generation rate is determined for various values of Darcy number (10^-3 ≤ Da ≤ 1) and for a range of Hartmann number (0 ≤ Ha ≤ 10²). This study is limited to Pr = 7 which corresponds to water. It is observed that entropy generation is very affected by the variations of dimensionless control parameters considered.

This research will investigate how heat treatments affect the behaviour of 3D printed gears made from Polyamide 6 (PA6), a thermoplastic polymer with excellent mechanical properties. However, imperfect structures can result from the 3D printing process, which involves rapid cooling of the molten polymer. Gears were manufactured from PA6 by means of a dedicated 3D printer and their behaviour in gears was then analysed by means of various specific test protocols. Some of the gears were subjected to annealing heat treatments to determine material behaviour, which was a key aspect of this study. Finite Element Analysis (FEA) was used to investigate the effect of heat treatments on the mechanical properties of PA6 thermoplastic material, specifically in the area of 3D printed gear mechanisms. The simulations included both static and transient conditions to determine the influence of different heat treatments on the mechanical properties of PA6. Through the use of comprehensive FEA, the research investigates the deformation, stress distribution and strain rates in 3D printed gears and elucidates the correlation between the mechanical properties and the contact pressure on the tooth flanks during gear operation. As conclusion the study mainly uses numerical simulations based on mechanical characteristics determined in our laboratory. The contribution of the study lies in the novelty of highlighting, through the results obtained, the influence of heat treatments on wear and mechanical stress related factors.

The dynamic response of the shed-tunnel structure under rockfall impact attracts much attention in recent years. However, the research on the dynamic response of the shed-tunnel structure under rockfall impact mostly focuses on the rockfall single impact, while in practical engineering, the shed-tunnel structure often encounters the rockfall repeated impacts during its service life. Therefore, this research focuses on exploring the dynamic response characteristics of the shed-tunnel structure under the rockfall repeated impacts. First of all, based on the model test of the shed-tunnel under rockfall impact as a reference, the FEM (Finite Element Method)-SPH (Smoothed Particle Hydrodynamics) coupling numerical calculation model is established based on ANSYS/LS-DYNA finite element code. The numerical simulation of the dynamic response of the shed-tunnel structure under rockfall impact is realized, and the rationality of the model is verified. Then, with this model and the full restart technology of LS-DYNA code, the effects of four factors, e.g. rockfall mass, rockfall impact velocity, rockfall impact angle and rockfall shape, on the impact force and impact depth of the buffer layer, the maximum plastic strain and axial force of the rebar, the shed roof vertical displacement and the shed-tunnel plastic strain are studied. These results show that the cumulative deformation and damage of the shed-tunnel structure caused by the rockfall repeated impacts play an important role in the destruction of the shed-tunnel structure, and the special attention should be paid to this issue.

The shaft hydropower plant (SHPP) is a novel hydraulic concept for low-head hydropower sites with several environmental and operational advantages over conventional layouts. However, the first two projects implementing this concept have shown comparatively high construction costs and project risks. Therefore, further optimization is required to the economic attractiveness and enabling broader market adoption. Initial model tests recommend a square-shaped shaft inlet with a three-sided approach flow for low-loss and fish-friendly inflow conditions. Yet, this design requires significant space for structural implementation and may be unsuitable for use with multiple shafts or as an extension of non-powered dams and weirs. This research paper presents the application of a computational fluid dynamics simulation setup to evaluate the hydraulic per-formance of various design configurations, especially, alternative design layouts with one-sided approach flow without further physical model tests,. The simulation setup is calibrated against observations including head loss and velocity measurements from the physical model tests, and its satisfactory performance enables the analysis of alternative design layouts. This aims to derive the most significant design parameters for achieving the desired hydraulic conditions at the intake. Increasing the flow depth before the intake and enlarging the inlet area have the most significant impact, while increasing the overflow of the front gate has the least significant effect.. The chosen CFD application is deemed suitable for hydraulic design optimization and provides guidance on the key parameters to focus on for tailored site-specific design development.

Modular steel buildings (MSBs) interlocking inter-module connections (IMCs) have piqued researchers’ interest. However, the existing literature has not yet studied the optimisation of its plate thicknesses. This paper, therefore, focuses on optimising interlocking (IMCs) plate thickness in (MSBs), leveraging previous experimental and numerical simulation methodologies. Various numerical models for four MSBs interlocking (IMCs) with plate thicknesses (4 mm, 6 mm, 10 mm, and 12 mm) were developed under compression loading conditions. The study’s novelty lay in its comprehensive parametric analysis, which evaluated the slip phenomenon in (MSBs) connections and introduced a slip prediction model validated against empirical data. Further innovative machine learning approach was utilised to predict slip values based on applied force through a random forest regression model trained using the 'Treebagger' function. Sensitivity analysis and comparisons with alternative methods were utilised to underscore the reliability and applicability of the findings. The results showed that 11.03 mm was the optimal plate thickness for interlocking (IMCs) in modular steel buildings, with up to 8.08% reduced material costs. Increasing plate thickness boosts its deformation resistance by up to 50.75%. ‘TreeBagger’ random forest for anomaly detection and machine learning improves slip prediction up to 7% at higher force levels.

The construction of tunnels excavated by the conventional method in densely populated urban environments requires an adequate characterization of the loads acting on the primary lining during the excavation process to ensure that the ground is deformed and stresses around the tunnel are relieved, simultaneously complying with the failure and serviceability limits of international standards while minimizing damage to nearby structures. In this paper, common lining design criteria are revisited, through the numerical simulation of an instrumented tunnel section which is part of a 4.5 km long metro line currently under construction in Mexico City. Key needs for improvement in current design approaches are identified. The tunnel was instrumented with load cells, extensometers, and topographical references for convergences and divergences. A three-dimensional finite difference model of the instrumented section was developed, and the load transfer mechanisms between the excavated soil and the primary lining was analyzed. Then, the numerical simulation of the contribution of the secondary lining in the overall stability for sustained load was established, along with the expected ground settlements, which can significantly affect nearby structures. Results gathered from this research are key for updating lining design criteria for urban tunnels built in stiff brittle soils.

Space-based gravitational wave detection is extremely sensitive to disturbances. The Keplerian configuration cannot accurately reflect the variations in spacecraft configuration. Planetary gravitational disturbances are one of the main sources. Numerical simulation is an effective method to investigate the impact of perturbation on spacecraft orbits. This study shows that, in the context of the Taiji project, Earth’s gravity is an essential factor in the change in heliocentric formation configuration, contributing to the relative acceleration between spacecrafts in the order of O(10−6)m·s−2. Considering 00:00:00 on October 27, 2032 as the initial orbiting moment, under the influence of Earth’s gravitational perturbation, the maximum relative change in armlengths and variation rates of armlengths for Taiji is 1.6×105km, 32m·s−1 respectively, compared with the unperturbed Keplerian orbit. Additionally, by considering the gravitational perturbations of Venus and Jupiter, the armlength and relative velocity for Taiji are reduced by 16.01% and 17.45%, respectively, compared with when only considering that of Earth. The maximum amplitude of the formation motion indicator changes with the orbit entry time. Results show that the relative velocity increase between the spacecrafts is minimal when the initial orbital moment occurs in July. Moreover, the numerical simulation results are inconsistent when using different ephemerises. The differences between ephemeris DE440 and DE430 are smaller than those between DE440 and DE421.

The fractional derivatives are a generalization the derivatives of integer order and find applications in studying memory processes in various scientific fields. Numerical methods are used to solve and analyze fractional models of real world problems. In this paper, we consider an approximation of the Caputo fractional derivative and its asymptotic formula, whose generating function is the polylogarithmic function. In the paper, we prove the convergence of the approximation and derive an estimate for the error and order. We consider an application of the approximation for the construction of finite difference schemes for numerical solution of the two-term ordinary fractional differential equation and the time-fractional Black-Scholes equation for option pricing. The properties of the approximation are used to prove the convergence of the methods used for numerical solution of the fractional differential equations. The theoretical results for the error and order of the methods are illustrated by the results of numerical experiments.

The annual energy matrix has been changing in the last years because of the necessity of less dependence on fossil fuels, which are running out on the planet. Therefore, a study was started to simulate the most efficient production of H2 through COMSOL Multiphysics®, as it is a gas that in the future will be essential to supply the world's energy necessities, due to its easy to get and insert in piped gas pipes. From this study, we were able to present an optimization of parameters that presents good indications of how to obtain efficient hydrogen production in an idealized reactor.

Tsunami modeling and simulation has changed in the past few years more than it had in decades, especially so with respect to coastal inundation. Among other things, this change is supported by the approaching era of exa-scale computing, whether via GPU or more likely forms of hybrid computing whose presence is growing across the geosciences. For reasons identified across this review, exa-scale computing efforts will impact the on-shore, highly turbulent r\'egime to higher degree than the 2D shallow water equations used to model tsunami propagation in the open ocean. This short review describes the different approaches to tsunami modeling from generation to impact and underlines the limits of each model based on the flow r\'egime. Moreover, from the perspective of a future comprehensive multi-scale modeling infrastructure to simulate a full tsunami, we underline the current challenges associated with this approach and review the few efforts that are currently underway to achieve this goal. A table of existing tsunami software packages is provided along with an open Github repository to allow developers and model users to update the table with additional models as they are published and help with model discoverability.

This paper proposes a lightweight joint with bamboo culm cross connection, which is constructed by screwing one side of the screw into the other side of the wooden plug, its simple construction makes it easy to manufacture and cost-effective. This form connects adjacent bamboo culms to form corner joints, and bolt connects adjacent bamboo culm frames to form edge joints, to explore the mechanical properties of the combination joints formed by corner joints and edge joints, by combining experiments and numerical simulations, the method discovered and explained the relationship and mechanism between the number and position of bolts and the stiffness of edge joints. A calculation model for frame-unit bamboo culm members considering semi-stiffness of joints is proposed, and the conversion equations between linear stiffness of corner and edge joint and their rotational stiffnesses are derived. By establishing the series and parallel relationship between the corner and edge joints in the assembled members, a generalized relational equation between the total stiffness of the assembled joint and the stiffness of the corner and edge joints is derived. A numerical simulation method for frame-unit bamboo culm members with semi-rigid joints is proposed, which can obtain the stiffness calculation formula of edge joints and the average rotational stiffness value of semi-rigid corner joints. A finite element analysis method based on semi-rigid joints is proposed for frame-unit bamboo culm members, and the internal force sketches and calculation methods for bamboo culms and joints are presented, which can be used for the internal force analysis and design of frame-unit bamboo culm lattice structures.

In this paper, low circumferential reciprocating load foot-scale tests were performed on two truncated PHC B 600 130 tubular piles with bearing nodes to characterize the damage process and morphology of the specimens and to investigate the load-carrying performance of the members. The test results reveal that under the action of tensile-bending-shear loading, the bearing concrete in the node area buckles and is damaged, the anchored reinforcement in the node area yields, the constraint is weakened, an articulation point is formed, and the node rotational capacity increases. When the embedment depth increases from 200 mm to 300 mm, the ultimate bearing capacities of the positive and negative nodes increase by 57.60% and 54.60%, respectively. A numerical simulation is used to verify the test results. Considering the three types of piles with truncated nodes, the numerical simulation is used to analyze the node-bearing capacity at different embedment depths. Finally, two preferred node types are proposed as follows: pipe pile core-filled longitudinal reinforcement anchored into the bearing node and pipe pile body longitudinal reinforcement anchored into the bearing node + pipe pile core-filled longitudinal reinforcement anchored into the bearing node, with preferred embedment depths of 350 mm and 200 mm, respectively.

Introduction: Alzheimer's disease (AD) is the first cause of dementia worldwide without any current curative treatment. Facing an increasing prevalence and its associated costs, AD represents a public health challenge. Usual diagnostic methods still rely on extended interviews and paper tests provided by an exterior examiner. We aim to create a novel, quick cognitive-screening tool on a numerical tablet. Methods: This program, built and edited with Unity®, runs on Android® for the Samsung Galaxy Tab S7 FE®. Composed of seven tasks inspired by the Mini-Mental Status Examination and the Montréal Cognitive Assessment, it browses several cognitive functions. The architectural design of this tablet application is distinguished by its multifaceted capabilities, encompassing not only seamless offline functionality but also a mechanism to ensure the singularity of data amalgamated from diverse sites. Additionally, a paramount emphasis is placed on safeguarding the confidentiality of patient information in the healthcare domain. Furthermore, the application empowers individual site managers to access and peruse specific datasets, enhancing their operational efficacy and decision-making processes. We performed a preliminary usability assessment among healthy subjects. Results: 24 healthy patients were included with a final F-SUS score of "excellent". Participants perceived the tool as simple to use and achieved the test in a mean time of 142 seconds. No technical errors occurred. Conclusion: These preliminary results suggest that our new assessment on a numerical tablet allows a short cognitive screening.

A method for solving the turbulence equations embedded in the sigma ocean general ocean circulation model is proposed. Like the general circulation model, the turbulence equations are solved using the splitting method by physical processes. The turbulence equations are split into two main stages describing transport-diffusion and generation-dissipation processes. Parameterization of turbulence in the framework of equations allows, at the generation-dissipation stage, to use both numerical and analytical solutions and to ensure high efficiency of the algorithm. The results of large-scale ocean dynamics simulation taking into account the parameterization of vertical turbulent exchange are considered. Numerical experiments were carried out using k-omega turbulence model embedded to the Institute of Numerical Mathematics Ocean general circulation Model (INMOM). Both the circulation and turbulence models are solved using the splitting method with respect to physical processes. The coupled model is used to simulate the hydrophysical fields of the North Atlantic and Arctic Oceans for 1948--2009. The model has a horizontal resolution of 0.25 degree and 40 sigma-levels along the vertical. The sensitivity of the solution to the changes in mixing parameterization is studied. Experiments demonstrate that taking into account the climatic annual mean buoyancy frequency improves the reproduction of large-scale ocean characteristics. There is a positive effect of Prandtl number variations for reproducing the upper mixed layer depth. The experiments also demonstrate the computational effectiveness of the proposed approach in solving the turbulence equations.

In the present paper, laminar forced convection nanofluid flows in a uniformly heated horizontal tube were revisited by direct numerical simulations. Single and two-phase models were employed with constant and temperature-dependent properties. Comparisons with experimental data showed that the mixture model performs better than the single-phase model in the all cases studied. Temperature-dependent fluid properties also resulted in a better prediction of the thermal field. A particular attention was paid to the grid arrangement. The two-phase model was used then confidently to investigate the influence of the nanoparticle size on the heat and fluid flow with a particular emphasis on the sedimentation process. Four nanoparticle diameters were considered: 10, 42, 100 and 200 nm for both copper-water and alumina/water nanofluids. For the largest diameter dnp = 200 nm, the Cu nanoparticles were more sedimented by around 80 %, while the Al₂O₃ nanoparticles sedimented only by 2.5 %. Besides, it was found that increasing the Reynolds number improved the heat transfer rate, while it decreased the friction factor allowing the nanoparticles to stay more dispersed in the base fluid. The effect of nanoparticle type on the heat transfer coefficient was also investigated for six different water-based nanofluids. Results showed that the Cu-water nanofluid achieved the highest heat transfer coefficient, followed by C, Al₂O₃, CuO, TiO₂, and SiO₂, respectively. All results were presented and discussed for four different values of the concentration in nanoparticles, namely φ = 0, 0.6, 1 and 1.6%. Empirical correlations for the friction coefficient and the average Nusselt number were also provided summarizing all the presented results.

Transition between free-surface and pressurized flows is an crucial phenomenon in many hydraulic systems, including water distribution systems, urban drainage systems, etc. During the transition, the force exerted on the structures changes drastically, thus it is meaningful to simulate this process. However, severe numerical oscillations are widely observed behind filling-bores, causing unphysical pressure variations and even computation failure. In this paper, some oscillation-suppressing approaches are reviewed and evaluated on a benchmark model. Then a new oscillation-suppressing approach is proposed to admit numerical viscosity when the water surface is at proximity of conduct roof which has first order accuracy. This approach adds numerical viscosity when water surface is at the proximity of conduct roof. It can sufficiently suppress numerical oscillations under an acoustic wave speed of 1000m/s and is simple to apply. In comparison with two experiments, the simulation results of this method show good agreement and little numerical oscillations. The results in this paper can help readers to choose an appropriate oscillation-suppressing method to improve the robustness and accuracy of flow regime transition simulations.

The study of the influence of the nanoparticle volume fraction and aspect ratio of microchannels on the fluid flow and heat transfer characteristics of nanofluids in microchannels is important in the optimal design of heat dissipation systems with high heat flux. In this work, the computational fluid dynamics method was adopted to simulate the flow and heat transfer characteristics of two types of water–Al2O3 nanofluids with two different volume fractions and five types of microchannel heat sinks with different aspect ratios. Results showed that increasing the nanoparticle volume fraction reduced the average temperature of the liquid–solid heat transfer surface and thereby improved the heat transfer capacity of the nanofluids. Meanwhile, the increase of the nanoparticle volume fraction led to a considerable increase in the pumping power of the system. Changing the aspect ratio of the microchannel effectively improved the heat transfer capacity of the heat sink. Moreover, increasing the aspect ratio effectively reduced the average temperature of the heating surface of the heat sink without significantly increasing the flow resistance loss. When the aspect ratio exceeded 30, the heat transfer coefficient did not increase with the increase of the aspect ratio. The results of this work may offer guiding significance for the optimal design of high heat flux microchannel heat sinks.

Advanced Finite Element (FE) modeling and simulations were performed on vehicular crashes into a Three-bar Metal Bridge Rail (TMBR). The FE models of a sedan, a pickup truck, and a TMBR section were adopted in the crash simulations subject to Manual for Assessing Safety Hardware (MASH) Test Level 2 (TL-2) and Test Level 3 (TL-3) requirements. The test vehicle models were first validated using full-scale physical crash tests conducted on a two-bar metal bridge using a sedan and a pickup truck with similar overall physical properties and sizes to their respective vehicles used in the simulations. The validated vehicular models were then used to evaluate the crash performance of the TMBR using MASH Evaluation Criteria for structural adequacy, occupant risk, and post-impact trajectory. The TMBR met all MASH TL-2 requirements but failed to meet the MASH TL-3 Criteria H and N requirements when impacted by the sedan. The TMBR was also evaluated under in-service conditions (behind a 1.52-m wide sidewalk) and impacted by the sedan under MASH TL-3 conditions. The simulation results showed that the TMBR behind a sidewalk met all safety requirements except for the occupant impact velocity in the longitudinal direction, which exceeded the MASH limit by 3.93%.

Based on the nonlinear failure criterion and modified tangential technique, the upper bound solutions of supporting pressure on the deep tunnel face were obtained under pore water pressure conditions. The influence of parameters on the supporting pressure and collapse range was investigated according to the unlimited block failure mechanism. It is found that the upper bound solutions of supporting pressure increase with the growth of the nonlinear coefficient and pore water pressure coefficient. The collapse range of the tunnel face scales out with the increase of the nonlinear coefficient and shrinks when increasing the pore water pressure coefficient. Moreover, with the increase of the nonlinear coefficient, the impact strength on supporting pressure and collapse range declines gradually. According to the calculating results, it is manifest that both the pore water pressure and nonlinear criterion factors have negative impacts on the stability of the tunnel face. Thus, more attention should be paid to these parameters to ensure face stability in deep tunnel construction.

This study introduces a new approach to modeling galaxy rotation curves by integrating concepts from Perelman’s functional with modifications to General Relativity. While initial numerical simulations indicate that our P-E Integration model may offer a closer fit to observed data than existing theories, we recognize the preliminary nature of these findings [14]. The results, though promising, should be interpreted cautiously. Further validation across different galactic systems and refinement of the model are necessary to confirm its broader applicability [15]. This work aims to contribute to the understanding of galactic dynamics, with the hope of opening new avenues for future research.

We provide a very general review of the resonant transmission line method for optical fiber problems. The method has been found to work seamlessly for a variety of difficult problems including elliptical and eccentric core fibers as well as “holey”, photonic crystal fibers. A latest version has been shown to offer great versatility with respect to cases of unconventional, inhomogeneous index profiles.

In this paper major machine learning (ML) tools and the most important applications developed elsewhere for numerical weather and climate modeling systems (NWCMS) are reviewed. NWCMSs are briefly introduced. The most important papers published in this field in recent years are reviewed. The advantages and limitations of the ML approach in applications to NWCMS are briefly discussed. Currently, this field is experiencing explosive growth. Several important papers are published every week. Thus, this paper should be considered a simple introduction to the problem.

Groundwater models serve as support tools to among others: assess water resources, evaluate management strategies, design remediation systems and optimize monitoring networks. Thus, the assimilation of information from observations into models is crucial to improve forecasts and reduce uncertainty of their results. As more information is collected routinely due to the use of automatic sensors, data loggers and real time transmission systems; groundwater modelers are becoming increasingly aware of the importance of using sophisticated tools to perform model calibration in combination with sensitivity and uncertainty analysis. Despite their usefulness, available approaches to perform this kind of analyses still present some challenges such as non-unique solution for the parameter estimation problem, high computational burden and a need of a deep understanding of the theoretical basis for the correct interpretation and use of their results, in particular the ones related to uncertainty analysis. We present a brief derivation of the main equations that serve as basis for this kind of analysis. We demonstrate how to use them to estimate parameters, assess the sensitivity and quantify the uncertainty of the model results using an example inspired by a real world setting. We analyze some of the main pitfalls that can occur when performing such kind of analyses and comment on practical approaches to overcome them. We also demonstrate that including groundwater flow estimations, although helpful in constraining the solution of the inverse problem as shown previously, may be difficult to apply in practice and, in some cases, may not provide enough information to significantly constrain the set of potential solutions.

Some cognitive functions, shared by humans and certain animals, were acquired early in the course of phylogeny and, in humans, are operational in their primitive form shortly after birth. This is the case for the quantification of discrete objects. The further phylogenetic evolution of the human brain allows such functions to be reconstructed in a much more sophisticated way during child development. Certain functional characteristics of the brain (plasticity, multiple cognitive processes involved in the same response, interactions and substitution relationships between those processes) provide degrees of freedom that open up the possibility of different pathways of reconstruction. The within- and between-individual variability of these developmental pathways offers an original window on the dynamics of development. Here, I will illustrate this theoretical approach to cognitive development—which can be called “reconstructivist” and “pluralistic”—using children's construction of number as an example.

In ocean engineering, accurate and efficient numerical simulations are crucial. Leveraging axisymmetry vastly enhances the efficiency of a simulation, reducing a three-dimensional scenario to a two-dimensional simulation. Exploiting axisymmetry enables a significant reduction in computational demand while maintaining simulation fidelity. In this paper, axisymmetry is leveraged to perform a large ensemble of simulations, in order to evaluate and maximise the accuracy and efficiency of a Computational Fluid Dynamics (CFD) simulation. In particular, the free decay motion and wave radiation from a heaving semi-submerged sphere are simulated, using the Reynolds Averaged Navier Stokes (RANS) solver, interFOAM, in the opensource finite volume CFD software OpenFOAM. Validated against highly accurate experimental data, extensive parametric studies are conducted, previously limited by computational constraints, which facilitate the refinement of simulation setups. More than 50 iterations of the same heaving sphere simulation are performed, informing efficient trade-offs between computational cost and accuracy across various simulation parameters and mesh configurations. Ultimately, by employing axisymmetry, this research contributes to the development of more accurate and efficient numerical modeling in ocean engineering.

The present study aims to develop an efficient numerical method for computing the diffraction and radiation of water waves with horizontal long cylindrical structures, such as floating breakwaters. A higher-order scheme is used to discretize geometry of the structure as well as the relevant physical quantities. As the kernel of this method, Wehausen’s free-surface Green function is calculated by a newly-developed Gauss-Kronrod adaptive quadrature algorithm after elimination of its Cauchy-type singularities. To improve computational efficiency, a Chebyshev approximation approach is applied to a fast calculation of the Green function that needs evaluation thousands of times. In addition, OpenMP parallel technique is used to the formation of influence coefficient matrix, which significantly reduces CPU time. Finally, computations are performed on wave exciting forces and hydrodynamic coefficients for the long cylindrical structures, either floating or submerged. Comparison with other numerical and analytical methods demonstrates good performance of the present method.

The paper continues the study of efficient algorithms for the computation of zeta functions over the complex plane. We aim to apply the modifications of algorithms to the investigation of underlying fractal structures associated with the Riemann zeta function. We discuss the computational complexity and numerical aspects of the implemented algorithms based on series with binomial-like coefficients.

A multiscale modeling concept for numerical simulation of multiphysics turbulent flow utilizing map-based advection is described. The approach is outlined with emphasis on its theoretical foundations and physical interpretations in order to establish the context for subsequent presentation of the associated numerical algorithms and the results of validation studies. The model formulation is a synthesis of existing methods, modified and extended in order to obtain a qualitatively new capability. The salient feature of the approach is that time advancement of the flow is fully resolved both spatially and temporally, albeit with modeled advancement processes restricted to one spatial dimension. This one-dimensional advancement is the basis of a bottom-up modeling approach in which three-dimensional space is discretized into under-resolved mesh cells, each of which contains an instantiation of the modeled one-dimensional advancement. Filtering is done only to provide inputs to a pressure correction that enforces continuity and to obtain mesh-scale-filtered outputs if desired. The one-dimensional advancement, the pressure correction, and coupling of one-dimensional instantiations using a Lagrangian implementation of mesh-resolved volume fluxes is sufficient to advance the three-dimensional flow without time advancing coarse-grained equations, a feature that motivates the designation of the approach as autonomous microscale evolution (AME). In this sense, the one-dimensional treatment is not a closure because there are no unclosed terms to evaluate. However, the approach is additionally suitable for use as a subgrid-scale closure of existing large-eddy-simulation methods. The potential capabilities and limitations of both of these implementations of the approach are assessed conceptually and with reference to demonstrated capabilities of related methods.

This paper presents a hydrodynamic numerical model for a pitching wave energy converter (WEC). The model uses potential wave theory and is based on Cummins' equation, with nonlinear hydrostatic restoring stiffness and excitation forces based on instantaneous body position and water surface elevation. The numerical model can include non-linear forces, like quadratic drag, power-take-off and other forces that may account for unknown viscous effects observed in experiments. The paper discusses the applicability and limitations of the code and presents the cases where assumptions and simplifications can be made. The goal is to conclude on the simplest, yet accurate, version of the model by evaluating its accuracy using experimental data. The case study for validation is Floating Power Plant's (FPP's) WEC [1]. In the full-scale commercial project, FPP's device consists of a semisubmersible platform hosting a wind turbine (5-8MW) and 4 WECs, each one connected to the platform by a rotation shaft. Due to the configuration of the platform, strong interactions occur between the WECs and the structure, as they are very closely spaced. In order to validate a numerical model able to simulate these hydrodynamic interactions, wave basin experiments with a similar but simplified setup were performed.

This paper investigates the effects that rain drop erosion has on the integrity of wind turbine materials. This research involves analysis of impact models and theoretical equations to simulate the impact phenomena when a raindrop collides with a wind turbine blade at high speed. Impact phenomena that occur when a raindrop impacts with a wind turbine blade were evaluated. Additionally, the use of various types of software was considered to simulate this impact event. The model results of the raindrop impact with respect to a comparison of the governing equations and impact behaviour are addressed. Future work is outlined including adding variables to the erosion modelling algorithms.

A reduced-dimension (RD) method based on the proper orthogonal decomposition (POD) technology and the linearized Crank-Nicolson mixed finite element (CNMFE) scheme for solving the 2D nonlinear Extended Fisher-Kolmogorov (EFK) equation is proposed. The method reduces CPU runtime and error accumulation by reducing the dimension of the unknown CNMFE solution coefficient vectors. For this purpose, The CNMFE scheme of the EFK equation is established and the uniqueness, stability and convergence of the CNMFE solutions are discussed. Subsequently, the matrix-based RDCNMFE scheme is derived by applying the POD method. Furthermore, the uniqueness, stability and error estimates of the linearized RDCNMFE solution are proved. Finally, numerical experiments are carried out to validate the theoretical findings. In addition, we contrast the RDCNMFE method with the CNMFE method, highlighting the advantages of the dimensionality reduction method.

Currently, the use of numerical models for reproducing the landscape evolution of a river basin is part of the day-by-day research activities of fluvial engineers and geomorphologists. However, despite landscape modelling is based on a rather long tradition, and scientists and practitioners are trying to schematize the processes involved in the evolution of a landscape since decades, there is still the need for improving both the knowledge of the physical mechanisms and their numerical coding. The present review focuses on the first aspect, discussing the main components of a landscape evolution model and their more common schematizations, presenting possible open questions to be addressed towards an improvement of the reliability of such kind of models in describing the fluvial geomorphology.

For the polymer electrolyte fuel cell (PEFC), the separator thickness plays an important role to determine the weight, the volume and the costs of the PEFC. In addition, thermal management, i.e. temperature distribution is also important for the PEFC system to obtain higher performance. However, there were few reports investigating the relationship between the temperature distribution and the power generation performance such as current density distribution of PEFC operated at higher temperature (HT-PEFC). The aim of this study is to clarify the impact of separator thickness on the relationship between the temperature distribution and the current density distribution of HT-PEFC. This study has also investigated the impact of separator thickness on the relationship between the temperature distribution and the gases such as H2, O2 and the current density distribution of HT-PEFC numerically using CFD software COMSOL Multiphysics. The separator thicknesses of 2.0 mm, 1.5 mm and 1.0 mm were studied. In the study, the operation temperature and the relative humidity (RH) of supply gas were also changed. As a result, it was revealed that the optimum separator thickness was 2.0 mm to obtain higher power generation of HT-PEFC. Since the heat capacity in the case of the separator thickness of 2.0 mm was the largest among the separators investigated in this study, the dehydration of PEM and catalyst layer was smaller compared with the thinner separator thicknesses. It was also revealed that the impacts of separator thickness on distribution of gases, e.g. O2, H2O, and current density distribution became larger under higher temperature and lower RH conditions. This study also concluded that the tendency of temperature distribution matched those of O2, H2O and current density distributions.

The scope of this work was part of our previous research on Pulsed Power Plasma Stimulation Technique. Electromagnetic fields generated during a pulsed power plasma discharge were simulated using a finite element method in the COMSOL RF module. The field distributions were calculated during and after the pulse. The current paper takes advantage of prior work measuring and modeling the current distribution in the plasma arc. Agreement with laboratory experimental measurements provides support for extension of the model to reservoir scale. The validated model was used to quantify the signal attenuation level in different medium environment.

A simple optimization strategy for the computation of 3D finite-differencing kernels on many-cores architectures is proposed. The 3D finite-differencing computation is split direction-by-direction and exploits two level of parallelism: in-core vectorization and multi-threads shared-memory parallelization. The main application of this method is to accelerate the high-order stencil computations in numerical relativity codes.

The previously developed models for fuel droplet heating and evaporation processes, mainly the Discrete Multi Component Model (DMCM), and Multi-Dimensional Quasi-Discrete Model (MDQDM) are investigated for the aerodynamic combustion simulation. The models have been recently improved, and generalised for a broad range of bio-fossil fuel blends so that the application areas are broadened with increased accuracy. The main distinctive features of these models are that they consider the impacts of species thermal conductivities and diffusivities within the droplets to account for the temperature gradient, transient diffusion of species and recirculation. A formulation of fuel surrogates is made, using the recently introduced model, referred to as ‘’Complex Fuel Surrogate Model (CFSM)’’ and analysing their heating, evaporation, and combustion characteristics. The CFSM is aimed to reduce the full composition of fuel to a much smaller number of components based on their mass fractions, and to formulate fuel surrogates. Such approach has provided a proof of concept with the implementation of the developed model into a commercial CFD code ANSYS-Fluent. A case study is made for the CFD modelling of gas-turbine engine using kerosene fuel surrogate. The surrogate is proposed using the CFSM. The model is implemented into ANSYS-Fluent via a user-defined function to provide the first full simulation of the combustion process. Detailed chemical mechanism is also implemented into ANSYS Chemkin for the combustion study.

A single fracture is the basic unit of fracture medium, and the roughness of fracture wall surface is an important factor influencing hydraulic characteristics of the flow in bedrock fracture. However, effects of the shape and density of roughness elements (various bulges/pits on rough fracture wall surfaces) on water flow in a single rough fracture have not been thoroughly discovered. Thus the water flow in single fracture with different shapes and densities of roughness elements was simulated by using Fluent software in this study. The results show that in wider fractures the flow rate mainly depends on fracture aperture, while in narrow and close fracture medium the surface roughness of fracture wall is the main factor of head loss of seepage; there is a negative power exponential relation between the hydraulic gradient index and the average fracture aperture, i.e. with the increase of fracture aperture, the relative roughness of fracture and the influence weight of hydraulic gradient both decrease; and in symmetrical-uncoupled fractures there is a super-cubic relation between the discharge per unit width and average aperture. Above results would help to deepen the understanding of rough fracture seepage.