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.

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 the last section we present a concerning numerical example and the respective software.

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 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.

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.

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.

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.

A novel concept of ion separation Anisotropic Diffusion Mobility Analyser (ADMA) with Individual Ion Counter (IIC) was shown by numerical modelling to be able to separate ions with high resolution and ultra-high sensitivity. An ADMA-IIC device has been developed and tested at sheath flow rates below 40 l/min with a number of chemicals relevant to airport security and military applications. The resolving power of the ADMA was observed to be in the range of 40 to 80. These data were close to the resolution obtained in numerical modelling predictions. It is greater than or comparable with the RP of currently available on-the-market devices for airport security and other applications. It was also confirmed that the ADMA technology enables very fast detection of analytes, at least 100 times faster than ion mobility spectrometers. It was shown that Brownian diffusion can be controlled and even overpowered by generating highly anisotropic diffusion in the ADMA device. The ultimate sensitivity of single-ion detection below the ppq range has been demonstrated.

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.

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.

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.

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.

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.

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.

: This study presents energy and exergy analyses conducted on a Gifford-McMahon cycle refriger-ator. Through detailed numerical simulations, the performance of the refrigerator is thoroughly assessed, with a particular focus on energy and exergy losses. The obtained results are compared with Sankey and Grassman diagrams, providing visual representations of energy and exergy flows within the system. The localization of exergy losses is carefully analysed to identify key areas for optimization aimed at enhancing efficiency and power output. By pinpointing sources of ex-ergy destruction, targeted improvements can be implemented to maximize the overall perfor-mance of the refrigerator. Overall, this research contributes to a deeper understanding of the thermodynamic behaviour of Gifford-McMahon cycle refrigerators and provides valuable insights for optimizing their performance in practical 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.

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.

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.

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.

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.

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.

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.

This study focuses on the optimal deployment problem, and determines the optimal size of a military force to send to the battle field. The decision is optimized, based on an objective function, that considers the cost of deployment, the cost of the time it takes to win the battle, and the costs of killed and wounded soldiers with equipment. The cost of deployment is modeled as an explicit function of the number of deployed troops and the value of a victory with access to a free territory, is modeled as a function of the length of the time it takes to win the battle. The cost of lost troops and equipment, is a function of the size of the reduction of these lives and resources. An objective function, based on these values and costs, is optimized, under different parameter assumptions. The battle dynamics is modeled via the Lanchester differential equation system based on the principles of directed fire. First, the deterministic problem is solved analytically, via derivations and comparative statics analysis. General mathematical results are reported, including the directions of changes of the optimal deployment decisions, under the influence of alternative types of parameter changes. Then, the first order optimum condition from the analytical model, in combination with numerically specified parameter values, is used to determine optimal values of the levels of deployment in different situations. A concrete numerical case, based on documented facts from the Battle of Iwo Jima, during WW II, is analyzed, and the optimal US deployment decisions are determined under different assumptions. The known attrition coefficients of both armies, from USA and Japan, and the initial size of the Japanese force, are parameters. The analysis is also based on some parameters without empirical documentation, that are necessary to include to make optimization possible. These parameter values are motivated in the text. The optimal solutions are found via Newton- Raphson iteration. Finally, a stochastic version of the optimal deployment problem is defined. The attrition parameters are considered as stochastic, before the deployment decisions have been made. The attrition parameters of the two armies have the same expected values as in the deterministic analysis, are independent of each other, have correlation zero, and have relative standard deviations of 20%. All possible deployment decisions, with 5000 units intervals, from 0 to 150000 troops, are investigated, and the optimal decisions are selected. The analytical, and the two numerical, methods, all show that the optimal deployment level is a decreasing function of the marginal deployment cost, an increasing function of the marginal cost of the time to win the battle, an increasing function of the marginal cost of killed and wounded soldiers and lost equipment, an increasing function of the initial size of the opposing army, an increasing function of the efficiency of the soldiers in the opposing army and a decreasing function of the efficiency of the soldiers in the deployed army. With stochastic attrition parameters, the stochastic model also shows that the probability to win the battle is an increasing function of the size of the deployed army. When the optimal deployment level is selected, the probability of a victory is usually less than 100%, since it would be too expensive to guarantee a victory with 100% probability. Some of many results of relevance to the Battle of Iwo Jima, are the following: In the deterministic Case 0 analysis, the optimal US deployment level is 66200, the time to win the battle is 30 days and 14000 US soldiers are killed or wounded. If the marginal cost of the time it takes to win a victory doubles, the optimal deployment increases to 75400, the time to win a victory is reduced to 26 days, and less than 12000 soldiers are killed or wounded. In the stochastic Case 0 analysis, the optimal US deployment level is 65000, the expected time to win the battle is 46 days and almost 25000 US soldiers are expected to be killed or wounded. If the cost per killed or wounded soldier increases from 0 to 5 M $US, the optimal deployment level increases to 75000. Then, the victory is expected to appear after 35 days and 19900 US soldiers are expected to be killed or wounded.

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.

In the process of use and manufacture, carbon fiber foam sandwich structures were often damaged by low-energy impact, resulting in performance degradation. Therefore, it was necessary to study the damage caused by low-speed impact of composite sandwich structures. In this paper, based on Hashin failure criterion, an equivalent finite element model of carbon fiber foam sandwich panel under low velocity impact was established. The model was used to simulate the damage of foam sandwich panel with [±45°/±45°/ (core) /±45°/±45°] ply structure under the impact energy of 10.58J, 21.17J, 31.75J and 42.34J. The simulation results of impact damage depth were compared with the experimental results. The error was less than 10%, which proved the rationality of the impact equivalent model. The model was used to predict and analyze the damage of foam sandwich panels with [±45°/ (core) /±45°], [±45°/ (0°, 90°) / (core) /±45°] and [±45°/ (0°, 90°) (core) / (0°, 90°) /±45°] ply structures under 21.17J impact energy. By comparing and analyzing the damage situation, impact force response time and impact velocity response time, the low energy impact resistance was analyzed. The results showed that increasing the number of ply structure of [±45°] can reduce the impact damage degree and improve the bearing capacity of sandwich panels.

Data obtained by Direct Numerical Simulations (DNS) of pressure-driven turbulent channel flow are studied in the range 180 ≤ Reτ≤ 10,000. Reynolds number effects on the Mean Velocity Profile (MVP) and second order statistics are analyzed with a view of finding logarithmic behavior in the overlap region or even further from the wall, well in the boundary layer’s outer region. Values of Kármán constant for the MVPs and the Townsend-Perry constants for the streamwise and spanwise fluctuation variances are determined for the Reynolds numbers considered. A data-driven model of the MVP, proposed and validated for zero pressure gradient flow over a flat plate, is employed for pressure-driven channel flow by appropriately adjusting Coles’ strength of the wake function parameter, Π. There is excellent agreement between the analytic model predictions of MVP and the DNS-computed MVP as well as of the Reynolds shear stress profile. The skin friction coefficient Cf is calculated analytically. The agreement between the analytical model predictions and the DNS-based computed discrete values of Cf is excellent.

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.

In rock blasting for mining production, stress waves play a major role in rock fracturing along with explosion gases. Better energy distribution improves fragmentation and safety, lowers production costs, increases productivity, and controls ore losses and dilution. Blast outcomes vary significantly with the choice of the explosive and the properties of the rock mass encountered. This study analyzes the effects of rock mass and explosives properties on the blast outcomes through numerical simulation using data from the case study and later validates the simulation results from the field blast fragmentation. The outcomes suggest that, for a given set of rock properties, the choice of explosive has a major influence on the resulting fragmentation. Strong explosives favor large fracture extents in hard rocks, while less strong explosives offer a better distribution of explosive energy and fractures. The presence of rock structures such as rock con-tacts and joints influences the propagation of stress waves and fractures depending on the structures' material properties, intensity and orientations, and the direction and strength of the stress wave. To achieve effective fragmentation, the blast design should mitigate the effect of variability in the rock mass by ensuring adequate energy distribution within the limits of geo-metrical design.

This paper introduces an innovative approach to numerically model Structure-Soil-Structure Interaction (SSSI) by integrating the Boundary Element Method (BEM) and the Finite Element Method (FEM) in a coupled manner. To assess the accuracy of the proposed method, a comparative study is undertaken, comparing its outcomes with those generated by the conventional FEM technique. Alongside accuracy, the computational efficiency aspect is crucial for the analysis of large-scale SSSI problems. Hence, the computational performance of the coupled BEM-FEM method undergoes a thorough examination and is compared with that of the standalone FEM method. The results from these comparisons illustrate the superior capabilities of the proposed method in comparison to the FEM method. The novel approach provides more reliable results compared to traditional FEM methods, serving as a valuable tool for engineers and researchers involved in structural analysis and design.

Alzheimer's disease (AD) is the first cause of dementia worldwide and represents a public health challenge. Current diagnostic methods still rely on extended interviews and paper tests. We aim to create a novel, quick cognitive-screening tool on a numerical tablet. 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 by allowing them to access and peruse specific datasets, thereby enhancing their operational efficacy and decision-making processes. We performed a usability assessment among 24 healthy patients 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, confirming that a short assessment on a numerical tablet is possible.

Due to the ongoing global warming on the Earth, permafrost degradation has been extensively taking place, which poses a substantial threat to civil and industrial facilities and infrastructure elements, as well as to the utilization of natural resources in the Arctic and high-latitude regions. In order to prevent the negative consequences of permafrost thawing under the foundations of constructions, various geophysical techniques for monitoring permafrost have been proposed and applied so far: temperature, electrical, seismic and many others. At the same time, non-stationary electromagnetic methods seem to have significant potential and a number of important practical advantages. We propose a cross-borehole exploration system for the highest localization of target objects in the cryolithozone. A novel mathematical apparatus for three-dimensional modeling of transient electromagnetic (TEM) signals by the vector finite element method has been developed. The original combination of the latter, the Sumudu integral transform and artificial neural networks makes it possible to examine spatially heterogeneous objects of the cryolithozone with a high contrast of geoelectric parameters, significantly reducing computational costs. We have studied the sensitivity of the TEM signals to the boundary between thawed and frozen rocks depending on the inter-borehole distance and spatial orientations of the transmitters and receivers. Numerical simulation results of the TEM monitoring of industrial facilities located on permafrost are considered. The formation of a talik has been shown to significantly manifest itself in the measured electromagnetic responses, which enables timely prevention of industrial disasters and environmental catastrophes.

This is the second in a series of two articles introducing the Gravitas computational general relativity framework, in which we now focus upon the design and capabilities of Gravitas's numerical subsystem, including its ability to perform general 3+1 decompositions of spacetime via the ADM formalism, its support for the definition and construction of arbitrary Cauchy surfaces as initial data, its support for the definition and enforcement of arbitrary gauge and coordinate conditions, its various algorithms for ensuring the satisfaction of the ADM Hamiltonian and momentum constraints, and its unique (totally-unstructured) adaptive refinement algorithms based on hypergraph rewriting via Wolfram model evolution. Particular attention will be paid to the seamless integration between Gravitas's symbolic and numerical subsystems, its ability to configure, run, analyze and visualize complex numerical relativity simulations and their outputs within a single notebook environment, and its capabilities for handling generic curvilinear coordinate systems and spacetimes with general (and often highly non-trivial) topologies using its specialized and highly efficient hypergraph-based numerical algorithms. We also provide illustrations of Gravitas's functionality for the visualization of hypergraph geometries and spacetime embedding diagrams in two and three dimensions, the ability for Gravitas's symbolic and numerical subsystems to be used in concert for the extraction of gravitational wave signals and other crucial simulation data, and Gravitas's in-built library of standard initial data, matter distributions and gauge conditions. We conclude by demonstrating how the numerical subsystem can be used to set up, run, visualize and analyze a standard yet nevertheless reasonably challenging numerical relativity test case: a binary black hole collision and merger within a vacuum spacetime (including the extraction of its outgoing gravitational wave profile).

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.

The loads generated inside agricultural silos under static and dynamic conditions depend on the mechanical properties of the materials stored inside them. The silo calculation methodologies are based on these mechanical properties. Although it is known that moisture content greatly influences the values reached by these mechanical properties, just few studies have been conducted to determine them. The present work determines the angle of internal friction, the apparent cohesion, the dilatancy angle and the apparent specific weight of wheat when subjected to different moisture contents. Direct shear and oedometer assay devices were used. In addition, a climatic chamber was used to moisten the wheat samples used in this work. The moisture content applied to the samples ranged from 7.5% to 15.5%. From the different assays carried out it could be observed that the values of the angle of internal friction and apparent cohesion obtained from direct shear assays, as well as those related to the apparent specific weight obtained from oedometer assays were like those found in the literature. The apparent specific weight increased as moisture content increased, reaching a peak value at a moisture content of 13.4%, followed by a drop at a moisture content of 15.5%. The values here provided can be used in silo load calculations involving numerical methods.

The application of numerical simulation methods for studying vertical-axis wind turbines offers significant advantages over conventional horizontal-axis wind turbines regarding wind energy utilisation, start-up wind speed, and environmental impact. This paper aims to investigate and establish the aerodynamic design performance of asymmetrical (NACA 4412) airfoils used in the application of vertical-axis wind turbines. The computational fluid dynamics (CFD) simulation uses Ansys software based on the Reynolds Averaged Navier-Stokes (RANS) equations. The aerodynamic parameters (pressure and velocity) are evaluated and compared within a range of angles of attack (AoA) 5-12 degrees at a fixed Reynolds number, as is a usual condition for wind turbine application. The lift coefficient (Cl), drag coefficient (Cd), and Cl/Cd ratio are obtained within the range AoA for the asymmetrical airfoils, and the performances have been predicted using the potential flow panel and CFD methods. This study suggests that while Xfoil may not be as precise as computational fluid dynamics (CFD), its predicted lift and drag coefficients generally align with experimental results. These help us to understand the best air-foil selection in a wind turbine application and predict its optimum AoA at which it performs at its best.

Certain configurations of liquefied natural gas refueling stations exhibit a deficiency in managing boil-off gas. Furthermore, the ill-conceived linkage between the submersible pump and the gas storage tank pipeline leads to impeded natural gas transmission. This study employs Computational Fluid Dynamics (CFD) methodology to scrutinize the hydrodynamic attributes of T-type tee and dovetail tee configurations implemented in the pipeline design connecting the submersible pump and storage tank in an LNG filling station across diverse operational scenarios. The T-type tee induces detachment of the primary flow from the inner wall due to inertial forces, resulting in vortex formation and heightened resistance, accompanied by increased energy dissipation. Conversely, the dovetail tee, with its circular wall transition, mitigates vortex generation, reducing the separation zone and consequently minimizing resistance and energy loss. The maximum static differential pressure between the inlet and outlet of the dovetail tee is diminished by 52.52% when compared to that of the T-type tee. In the most noteworthy scenarios, the local resistance coefficient of the dovetail tee is even less than 1/10 of that of the T-type tee. In practical engineering applications, the incorporation of dovetail tees leads to a 17.58% reduction in the required elevation difference. This results in a more uniform flow field and enhanced stability in both flow rate and pressure. Such improvements contribute to heightened engineering efficiency and environmental sustainability, particularly evident in the design of LNG filling stations.

The emerging industry of offshore wind turbines mounted on floating bases has garnered significant attention from both academia and industry. The desire to understand the complex physics of these floating structures has led to the development of numerical and physical modelling techniques. While physical testing has traditionally been employed, there is a growing focus on cost-effective and accurate high-fidelity numerical modelling as a potential alternative or supplement. However, commonly used numerical engineering tools in the offshore industry are considered mid- to low-fidelity and may lack the desired precision for floating offshore wind turbines (FOWTs). Given the complexity of these simulation codes, it is crucial to validate their accuracy. To address this, the International Energy Agency (IEA) Wind Technology Collaboration Programme initiated various research endeavors, including the Offshore Code Comparison Collaboration (OC3), Offshore Code Comparison Collaboration Continuation (OC4), Offshore Code Comparison Collaboration Continuation with Correlation (OC5), and the recent Offshore Code Comparison Collaboration Continued with Correlation and Uncertainty (OC6) projects. This study offers a comprehensive survey of the simulation tools available for FOWTs which were part of OC projects, focusing particularly on horizontal axis wind turbines (HAWTs) and highlighting their capabilities and fundamental theories.

A problem of efficacy of crystal growth methods for crystallization from solution or melt has been investigated. Axial vibrational control (AVC) technique was considered as a perspective method to manage both heat-mass transfer and chemical component composition of the melts in the case of crystallization of complex chemical compounds. Numerical modeling and the search for generalized dependencies made it possible to predict the AVC parameters that provide optimal heat and mass transfer modes for creating flat liquid-solid interfaces, as well as the component composition of dissociated melts of various chemical compounds – Ge, NaNO3, CdTe.

The results of numerical experiments with the chemistry-climate model of the Institute of Numerical Mathematics of the Russian Academy of Sciences - Russian State Hydro-Meteorological University (INM-RSHU CCM) in comparison with the MERRA2 reanalysis data are presented. To assess the impact of sea surface temperature (SST) on stratospheric processes and the ozone layer, numerical experiments were carried out for scenarios corresponding to the El Niño and La Niña Southern Oscillation phases. It is shown that the El Niño phase contributes to the enhancement of heat and mass transfer to the polar stratosphere, which strengthen the Brewer-Dobson circulation and the flux of wave activity that affect the zonal wind. As a result, during El Niño, the probability of sudden stratospheric warming increases and the stratospheric polar vortex (SPV) weakens, re-sulting in an increase in total polar stratospheric ozone. The La Niña phase, on the contrary, is accompanied by a weakening of heat and mass transfer from the troposphere to the stratosphere, as a result of which the stratospheric polar vortex strengthens and the total column ozone decreases. It is shown that the effect associated with the Southern Oscillation is more pronounced in the Northern Hemisphere than in the Southern Hemisphere. Numerical experiments were also carried out for scenarios in which the sea surface temperature was set corresponding to the early 1980s, and the results were compared with calculations in which the SST was set corresponding to the late 2010s. This allowed us to assess the influence of the multiyear trend in the SST increase on stra-tospheric processes and the ozone layer. It was shown that the SST increase is most pronounced in the Arctic, which corresponds to the phenomenon of "Arctic amplification", and therefore, the heat and mass transfer to the stratosphere increases in the Northern Hemisphere, which contributes to some weakening of the SPV and increase of the total column ozone. But, compared to the Southern Oscillation, the influence of the SST trend on the SPV and total column ozone is much weaker and is observed mainly in the Northern Hemisphere.

The study examines the behavior of the Mechanically Stabilized Earth (MSE) wall using available reinforcement materials in Danang, Vietnam. The MSE was reinforced by the self-fabricated galvanized steel grids using CB300V steel with 3 cm ribs. The backfill soil is a sandy clay soil from the local area with a low cohesion. A full-scale model with full instrumentation was installed to investigate the distribution of tensile forces along the reinforcement layers. The highest load that caused the wall to collapse due to internal instability (reinforcement rupture) was 302 kN/m2, which is 15 times the design load of 20 kN/m2. The failure surface within the reinforced soil had a parabolic sliding shape which is similar to the theoretical studies. At the failure load level, the maximum lateral displacement at the top of the wall facing was small (3.9 mm) which is significantly lower than the allowable displacements for the retaining wall. Furthermore, a numerical model using FLAC software was applied to simulate the performance of the MSE wall. The modeling results are in good agreement with the physical model. Thus, the self-fabricated galvanized steel grids could be confidentially used in combination with the local backfill soil for the MSE walls.

This is the second paper in a cycle investigating the exact solution of loop equations in decaying turbulence. We perform numerical simulations of the Euler ensemble, suggested in the previous work, as a solution to the loop equations. We designed novel algorithms for simulation, which take a small amount of computer RAM so that only the CPU time grows linearly with the size of the system and its statistics. This algorithm allows us to simulate systems with billions of discontinuity points on our fractal curve, dual to the decaying turbulence. For the vorticity correlation function, we obtain \textbf{quantum fractal laws} with regime changes between universal discrete values of indexes. The traditional description of turbulence with fractal or multifractal scaling laws does not apply here. The measured conditional probabilities of fluctuating variables are smooth functions of the logarithm of scale, with statistical errors being negligible at our large number of random samples $T = 167,772,160$ with $N = 200,000,000$ points on a fractal curve. The quantum jumps arise from the analytical distribution of scaling variable $X = \frac{\cot^2(\pi p/q)}{q^2}$ where $\frac{p}{q}$ is a random fraction. This distribution is related to the totient summatory function and is a discontinuous function of $X$. In particular, the energy spectrum is a devil's staircase with tilted uneven steps with slopes between $-1$ and $-2.03$. The logarithmic derivative of energy decay $n(t)$ as a function of time is jumping down the stairs of universal levels between $1.50$ and $1.53$. The quantitative verification of this quantization would require more precise experimental data.

This study explores how the properties of heater materials impact nucleate pool boiling using a comprehensive simulation model. The model integrates well-known heat transfer mechanisms, allowing us to assess the effects of two distinct heater materials. Findings indicate that materials with higher thermal conductivity, such as copper, notably improve refrigeration efficiency in nucle- ate boiling. The study provides insights into the relationship between bubble growth, microlayer recovery beneath a bubble, temperature fluctuations, and heater properties. Comparisons between copper and silicon oxide underscore variations in bubble frequency, attributed to differences in bubble growth time, microlayer recovery time, and material-dependent behavior. The influence of neighboring bubbling sites, especially pronounced in silicon oxide due to its low thermal conductivity and diffusivity, was observed. Temperature variations in such materials became highly visible due to their significant inertia in recovery. Simulation results align well with empirical correlations, confirming the model’s versatility in capturing the intricate phenomena of nucleate pool boiling.

The review aims to examine how well numerical indicators of business bankruptcy can predict outcomes. The paper examines ten critical studies that concentrate on statistical models for bankruptcy predictions utilizing PRISMA criteria. The findings highlight the usefulness of numerical indicators in indicating financial hardship, particularly cash flow ratios. In order to connect theoretical knowledge with real-world corporate strategy applications, the study ends by reaffirming the importance of these indicators in strategic decision-making.

We study the tight focusing properties of generic bright and dark HyOPS beams in the presence of primary coma. The role of the polarization singularity index and handedness of the polarization of the HyOPS beams on the focused structure has been discussed. Results have been presented for the total intensity, component intensities, and component phase distributions for left-, and right-handed bright and dark star and lemon types singularities. The presence of primary coma distorted the focal plane intensity distributions for both positive and negative index generic C-points. Coma is known to disturb the circular symmetry of the focal plane intensity distribution. Similarly in tight focusing polarization is known to disturb the symmetry. Therefore, a beam with structured and inhomogeneous polarization distribution tightly focused under the influence of coma is a fit case to study. It is found that the presence of primary coma aberration in the focusing system produces a positional shift of the high-intensity peaks and a reduction of the intensity on one side of the center. As the strength of the primary coma increases, the focal plane intensity distributions shift more and more toward the right from the initial position. Unlike the scalar vortex case, in the case of HyOPS beams, the focal plane intensity distribution undergoes rotation, as the helicity of the HyOPS beams is inverted, in addition to shifting. All the component phase distributions are found to be embedded with phase vortices of charge ±1.

This paper offers a comprehensive examination of the process involved in developing and automating supervised end-to-end machine learning workflows for forecasting and classification purposes. It offers a complete overview of the components (i.e. feature engineering, model selection, etc), principles (i.e. bias-variance decomposition, model complexity, overfitting, model sensitivity to feature assumptions and scaling, output interpretability, etc), models (i.e. neural networks, regression models, etc), methods (i.e. Cross-Validation, data augmentation, etc), metrics (i.e. Mean Squared Error, F1-score, etc) and tools that rule most supervised learning applications with numerical and categorical data, as well as their integration, automation, and deployment. The end goal and contribution of this paper is the education and guidance of the non-AI expert academic community over complete and rigorous machine learning pipelines and data science practices, from problem scoping to design and state-of-the-art automation tools, including basic principles and reasoning in the choice of methods. The paper delves into the critical stages of supervised machine learning workflow development, many of which are often omitted by researchers due to brevity, and covers foundational concepts essential for understanding and optimizing a functional machine learning workflow, thereby offering a holistic view of task-specific application development for applied researchers who are not AI experts.

The problem of intensification of melt crystal growth process has been analyzed using CdTe as an actual material. Numerical simulation of 100 mm diameter CdTe crystal growth by VGF technique has been carried out. The heat-mass transfer was controlled by introducing of low-frequency oscillating baffle into the melt, so called axial vibrational control (AVC) technique. The baffle configuration has been optimized to destroy solid «tails», which were formed near the crucible walls at high cooling rates due to the low thermoconductivity and the corresponding latent heat. Analysis of CdTe homogeneity range showed that during fast crystal cooling Te micro precipitations were formed resulting from decay of oversaturated Cd-riched nonstoichiometric solid solution at Bridgman crystal growth technique. A VGF grown CdTe crystal after full crystallization stays inside the phase field of the high temperature wurtzite polymorph. This makes it possible to go through the polymorph transition without Te micro-precipitating using the advantages of VGF specific feature of very slow cooling.

Nuclear fusion reactors are composed by several complex components whose behaviour may be not certain a priori. This uncertainty may have a significant importance in the evolution of fault transients in the machine, causing unexpected damages to its components. For this reason a suitable method for the uncertainty propagation during those transient is required. The Monte Carlo method would be the reference option, but it is, in most of the cases, not applicable due to the large amount of required simulations to be repeated. In this context the Polynomial Chaos Expansion has been considered as a valuable alternative. It allows to create a surrogate model of the original one in terms of orthogonal polynomials. Then, the uncertainty quantification is performed repeatedly calling this much simpler and faster model. Using the fast current discharge in the Divertor Tokamak Test Toroidal Field (DTT TF) coils as a reference scenario, this method has been applied: the uncertainty on the parameters of the Fast Discharge Unit (FDU) varistor disks is propagated to the simulated electrical and electro-magnetic relevant effects. Eventually two worst case scenarios have been analyzed from the thermal-hydraulic point of view with the 4C code, simulating a fast current discharge as a consequence of a coil quench. Eventually it has been demonstrated that the uncertainty on the inputs (varistor parameters) strongly propagates obtaining a wide range of possible scenarios in case of accidental transients. This result underlines the necessity to take into account and propagate all possible uncertainties in the design of a fusion reactor, according to the Best Estimate Plus Uncertainty approach.

In the case of solid slabs made of reinforced concrete that are usually subjected to bending, large areas of the structure are stressed well below their load-bearing capacity or remain stress-free. Contrary to this are shell structures, which can bridge large spans with little material if designed well. To improve the efficiency of ceiling slabs, we want to utilize the shell load-bearing behaviour on a smaller scale by dissolving the solid interior accordingly. In order to be able to study a wide range of such constructions virtually, a parametric multiobjective simulation environment is to be developed in an ongoing research project, basic analysis approaches of which are presented in this paper. In addition to the basic workflow and the programs used, the material models used, and their calibration are described on the basis of tests on textile reinforced concrete (TRC) samples. Initial tests on two individual TRC served to verify the model.

Hollow and concrete-filled steel tubes (CFSTs) are widely used as columns in many structural systems and local buckling can occur when they are subjected to axial compression loading. The critical regions for the local buckling of the columns are near the ends of the column where the moments are the highest. To strengthen the critical regions of the CFSTs the CFRP provides a simple and effective method which used previously as an enhancement for both of strength and ductility. So, this study suggests investigating the axial load behaviour of CFST strengthened with CFRP using FE method after a verification study and then a various parametric studies were carried out to cover the effect of the CFRP layers, number, confinement length and position. The FE results proved that, one CFRP sheet with thickness equal to 1.2 mm was strengthened the composite column by 8.5%, two CFRP sheets with total thickness equal to 2.4 mm was strengthened the composite column by 23.5%, three CFRP sheets with total thickness equal to 3.6 mm was strengthened the composite column by 35.1% and four CFRP sheets with total thickness equal to 4.8 mm was strengthened the composite column by 44.5%. Moreover, this study proved that by strengthened the composite columns with CFRP sheets with different lengths, the axial load resistance was improved by 8.5%, 4.6%, 0.1% and 0.0% for length percentages 100%, 75%, 50% and 25% respectively.

In the case of solid slabs made of reinforced concrete that are usually subjected to bending, large areas of the structure are stressed well below their load-bearing capacity. Contrary to this are shell structures, which can bridge large spans with little material if designed according to the force flow. To improve the efficiency of ceiling slabs, we want to utilize the shell load-bearing behavior on a smaller scale by resolving the solid interior accordingly. In order to be able to study a wide range of such constructions virtually, a parametric multiobjective simulation environment is being developed in an ongoing research project. The basic analysis approaches which were implemented are presented in this paper. The basic workflow, the used programs and material models and their calibration on the tests on textile-reinforced concrete (TRC) samples are described.

A mathematical model that comprehensively captures the real behavior of riverbed deformation, encompassing all pertinent effects, is developed. The underwater slope reformation process, with the generatrix aligned along the flow velocity in the model, is considered. A numerical model is introduced to calculate the flow involving a deformable bottom, and the model's validation is established through rigorous analysis of experimental findings. The research firmly confirms the suitability of the proposed mathematical and numerical model for describing deformations in uneven and unsteady river flows, including the movement of dredging slots and channel quarries. The model's minimal equation count and reliance on empirical constants demonstrate the model's efficiency. The model's predictions align strongly with experimental data, although optimal values of empirical coefficients vary slightly across different experiments. Hence, there is a call for further investigation to derive more universally applicable closure relationships for the model. The importance of validating the model with reliable field data and its potential extension to accommodate hydraulically diverse soils is emphasized. Such an extension is feasible due to the concentration transfer equation, enabling independent calculations for particle fractions of varying sizes as long as the total particle concentration in the stream remains within reasonable limits. This dedicated research contributes significantly to understanding riverbed deformations and advancing accurate modeling and management of riverine environments.

The research on the impact of excavation of underpass tunnels has already had certain results, but there is a lack of research cases of double-line highway tunnels with oblique angles underpass the existing railway tunnels, especially the method of determining the area of the impact of the dynamic and static loads of the new tunnels and existing trains to be enriched by considering them at the same time. In this paper, based on the tunnel project of the new double-line Shiqian Highway Tunnel with oblique angles underpass the Hurong Railway in Wanshoushan. By constructing a three-dimensional finite element model and simulating the application of static tunnel excavation load, dynamic blasting load and dynamic train operation load of the existing tunnel, the overall displacement of the existing tunnel and the settlement value of the bottom plate of the track surface by the static tunnel excavation load were analyzed. Then, the stresses, vibration speeds and displacements of the tunnel due to the dynamic blasting loads and the dynamic train operation loads were obtained. The results show that: the area of influence of the static loads of the new tunnel excavation on the existing tunnel is divided into three types of perturbations presenting strong, weak and slight; the area affected by the blasting dynamic load of new tunnel excavation is a circular domain with the datum point in the tunnel section as the center; the area affected by the dynamic load of train operation in the existing tunnel is an ellipse with the center of the track surface of the existing tunnel as the datum; there is an anomalous shape in the area affected by the blasting dynamic load of new tunnel excavation and the area affected by the dynamic load of train operation in the existing tunnel.

Plastic pollution of the ocean is currently one of the most serious environmental threat. Analysis of observations in the Arctic ocean shows that microplastic particles have been found in Arctic snow cover, sea ice, water and sediments. Possible pathways and deposition locations of microplastics in the Arctic Ocean can be assessed based on 3D numerical modeling of thermohaline structure and water circulation, sea ice conditions and drift. The problem is that in addition to transport by currents, a number of other processes affect the state of microplastics. These include mechanical fragmentation, freezing into sea ice and release back into the water, biofouling, adsorption of contaminants and other factors. In this study, based on scenario calculations for a 5-year period, we analyzed the possible spread of microplastics brought by Siberian rivers to the Kara Sea shelf. The Lagrangian particle model used daily 3D numerical simulation data to simulate microplastic transport by ocean currents and by sea ice drift. The results of a series of scenario calculations show how the distribution of particles and their subsequent deposition depend on their type (density), size, processes of freezing into the ice and biofouling.

In the southeastern coast of Brazil, the bays of Ilha Grande and Sepetiba are linked by the Ilha Grande Channel, where remarkably vigorous currents have been consistently observed. The sheer intensity of these currents defies a simple explanation through tidal forcing alone, prompting researchers to delve into potential influences such as basin resonance and the baroclinic effect. This study aims to elucidate the role of remote meteorological effects within this complex hydrodynamic system. A numerical approach with a coastal model nested within an ocean model was employed, enabling an in-depth examination of the intricate interplay between meteorological and tidal forcings. The study revealed a significant finding: the lag in signal propagation plays a pivotal role in determining how these signals impact the dynamics of the bays. The astronomical signal exhibits a minimal lag along the coast (1 minute), and leads to water level differences between the sea and the coastline, resulting in the generation of tidal currents at the bay entrances. On the other hand, the remote meteorological signal, with a stronger signal lag along the coast (4.92 hours), leads to the creation of water level difference between the bay’s entrances, inducing significant fluxes along the narrow Ilha Grande Channel.

The first paper of this two-part series paper discussed the development of a novel lightweight 3D wood dowel-reinforced glass epoxy hybrid composite material (3DdrFRP) and its manufacturing procedures. It also experimentally compared the performance of scaled utility poles made from conventional 2D E-glass epoxy and 3DdrFRP materials. In the second part, we developed robust finite element (FE) models in LS-DYNA, calibrated them based on experimental results, and used them to simulate the performance of equivalent 2D and 3D poles, proving the integrity of the numerical model. Additionally, we developed a simplified analytical calculation method for engi-neers to determine the stiffness of 3D-DrFRP poles accurately and quickly.

Recently, more wind turbine systems are being installed in deep waters far from the coast. Several concepts of floating wind turbine systems (FWTSs) are developed, among which the semisubmersible platform, due to its applicability in different water depths, good hydrodynamic performance, and facility in the installation process, constitutes the most explored technology compared to the others. However, a significant obstacle to the industrialization of this technology is a design of a cost-effective FWTS, which can be achieved by optimizing the geometry, size and weight of the floating platform, along with the mooring system. This is only possible by selecting a method capable of accurately analyzing the FWTS coupled hydro-aero-structural dynamics at each design stage. Accordingly, this paper aims at providing a detailed overview of the most common coupled numerical and physical methods, including their basic assumptions, formulations, limitations, and costs, used for analyzing the dynamics of FWTSs, mainly those supported by a semisubmersible, to assist the choice of the most suitable method at each design phase of FWTSs. Finally, the article discusses possible future research to address challenges in modeling FWTSs dynamics that persist to date.

The CMT (Cold Metal Transfer) arc welding process is a convenient instrument for the assembly of thin sheets [1], [2]. More particularly, it permits the decrease of the amount of heat transferred to the parts to-be assembled, thus reducing their subsequent deformations. The present work aims at simulating the CMT welding of thin stainless-steel sheets in order to develop an approach that allows the prediction of temperature fields, as well as final deformations induced by the welding process. For that purpose, instrumented tests and numerical simulations were set aiming at comparing the experiments with simulations. Butt-welding of stainless-steel sheets was considered, with a thickness ranging from 1 to 1.2 mm. The weld seam samples were observed in order to set an equivalent heat source for each configuration. Moreover, the electric current and voltage were surveyed, and temperature measurements were performed using K-type thermocouples. In addition, displacement measurements were obtained using the DIC (Digital Image Correlation) technique. Thermomechanical simulations were then carried out using an equivalent heat source approach, considering the phase changes of elements, i.e. the transitions from solid to liquid and from liquid to solid. These simulations also consider the actual geometry of seams induced by addition of a filler material.

When the formation flying consists of satellites that are independent from each other and forms a virtual synthetic aperture telescope that replaces a single large telescope with all of its optical elements it is required mutual coordination between the satellites with high accuracy in order to achieve good optical observations. In this paper we consider the problem of developing a control system for a tetrahedral satellite formation flying representing a synthetic aperture telescope in the geostationary orbit for high resolution monitoring of fire in the forest area in infrared spectrum. Synthesizing the image with good quality requires to keep the configuration of formation with µm-class accuracy. Study of dynamics of passive uncontrolled motion of formation flying showed a significant deviation of the configuration from the required shape with a high frequency due to the action of gravitational forces. To keep the configuration in the required form with µm accuracy the analysis of efficiency of various controllers was carried out in the process of numerical simulation. The simulation results made it possible to highlight the features of using various approaches to the development of control system for the satellite formation flying representing a synthetic aperture telescope.

In the research we consider a structure-fluid interaction of oscillating piezodevice and acoustic liquid. Such problem requires finding resonance or anti-resonance frequencies of the solid structure to use them in the latter frequency-response analyses step. After finite element discretization of the issue, we need to solve a generalized eigenvalue problem (GEP). The most efficient numerical method for solving the partial GEP is the shift-and-invert Lanczos method. Each step of this algorithm requires solving a large scale system of linear algebraic equations (SLAEs) with special properties. Various ways of ordering grid nodes lead to coefficient matrices with different structure. Several methods are compared in terms of efficiency. The results of the numerical experiments are presented.

In this paper, four lithological parameters, i.e. thermal conductivity, particle size, porosity and saturation, are investigated by combining experimental observations with numerical simulations to study the influence laws and mechanisms of action on the combustion process of methane hydrate sediments. The variation of combustion characteristics parameters such as flame height, effective combustion ratio and dimensionless discharge water mass with the lithological parameters are studied. In addition, the combustion mechanism of lithologic parameters on methane hydrate deposits is revealed. Combining the experimental results and simulation calculations, to optimize the combustion of methane hydrate sediments, it is recommended to use methane hydrate sediment samples with high saturation and low thermal conductivity, while the oxidant concentration and porosity of methane hydrate sediment samples should be increased.

The treatment of critically ill patients remains an evolving and controversial issue. Mechanical circulatory support through a percutaneous approach for the management of cardiogenic shock has taken place in recent years. The combined use of IABP and Impella 2.5 device may have a role to play for this group of patients. A simulation approach may help with a quantitative as-sessment of the achievable degree of assistance. In this paper, we analyse the interaction be-tween Impella 2.5 pump applied as LVAD and IABP using the numerical simulator of the cardi-ovascular system developed in our laboratory. Starting with pathological conditions reproduced using literature data, the simulations were performed setting different rotational speeds for the LVAD and driving the IABP in full mode (1:1) or partial mode (1:2, 1:3 and 1:4). The effects in-duced by drug administration during the assistance were simulated too. The haemodynamic pa-rameters under investigation were: mean aortic pressure, systolic and diastolic aortic pressure, mean pulmonary arterial pressure, mean left and right atrial pressure, cardiac output, cardiac index, left and right ventricular end-systolic volume, left ventricular end-diastolic volume and mean coronary blood flow. The energetic variables considered in the study were: left and right ventricular external work and left and right atrial pressure-volume area. The outcome of our simulations shows that the combined use of IABP and Impella 2.5 achieves adequate support in the acute phase of cardiogenic shock compared to each standalone device. This would allow fur-ther stabilisation and transfer to a transplant centre should escalation of treatment be required.

This paper presents the results of laboratory tests of CFRP (carbon fibre-reinforced polymer) laminates, which allowed the development of numerical material models. The obtained data were used in a further stage to perform numerical simulations of four variants of medical tabletops differing, among other features, in the shape of the cross-section. Maximum deflections and effort in the composite material were analyzed. The final step was to perform a laboratory test for one of the tabletop versions, the results of which confirmed the correctness of the numerical calculations. This work is aimed at both researchers and designers involved in the practical application of fibre-reinforced polymer matrix composites.

The soil-structure interaction represents an important parameter when analyzing the dynamic behavior of bridges, taking into account the integration between soil, foundation, and structure to ensure an adequate structural representation. In this paper, a computational modeling strategy using soil-structure interaction (SSI) methodology was applied to represent the dynamic behavior of a pedestrian bridge (LOP). The computational models were validated using monitoring results from an extensive structural health monitoring campaign. The bridge monitoring program was characterized using embedded sensors and cutting-edge technology for remote monitoring: ground-based radar interferometry and video motion amplification. The instrumentation results were compared to establish the capabilities and advantages of this novel remote monitoring technique. The experimental results obtained from free and forced vibrations in the asset served as a basis for validating the developed computational model and allowed evaluating the asset be-havior in service for different loading conditions. This paper emphasizes the SSI approach to represent the pedestrian bridge dynamic behavior more accurately, allowing computational models to represent the LOP structural responses over the course of its life-cycle.

Folk residence in the cold region of Northeast China is an important part of China's traditional folk residence, with significant geographical representation. In this paper, we take the residential compounds of villages and towns in the cold region of Northeast China as the research object, summarize four typical layout forms of residential compounds of villages and towns in the cold region of Northeast China, simulate the wind environment of compounds of different layout forms in winter by using the computer numerical simulation software Fluent, and then analyze and evaluate the wind environment of compounds by combining with the evaluation standard of wind speed ratio and considering the area of wind shadows, number of vortices and other factors.

This study develops a new definition of fractional derivative that mixes the definitions of fractional derivatives with singular and non-singular kernels. Such developed definition encompasses many types of fractional derivatives, such as the Riemann-Liouville and Caputo fractional derivatives for singular kernel type as well as the Caputo-Fabrizio, the Atangana-Baleanu and the generalized Hattaf fractional derivatives for non-singular kernel type. The associate fractional integral of the new mixed fractional derivative is rigorously introduced. Furthermore, newly numerical scheme is developed to approximate the solutions of a class of fractional differential equations (FDEs) involving the mixed fractional derivative. Finally, an application to computational biology is presented.

A teoria diz que o número interage com os espaços com base no significado da matemática, desta forma é possível obter conhecimento numérico. A teoria mostra a possibilidade de usar espaços relacionados entre si para explicar as razões pelas quais dois números pares formam um par, dois números ímpares formam um par, um número par e um ímpar formam um número ímpar, três números ímpares formam um número ímpar, três números pares formam um número ímpar. um par, um ímpar e dois pares formam um ímpar, dois ímpares e um par formam um par.

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.

An investigation of compressible flow over two-dimensional obstacles is presented in this paper. When shock waves are formed on surfaces due to compressibility and these shocks interfere with the boundary layer, undesirable effects such as friction, separation, and flow instability occur. To weaken the shock and reduce its effects, various methods have been employed in the past, including surface suction, surface blowing, and vortex generators. Using numerical solutions, the present work attempts to reduce the negative effects of the shock generated by the obstacle by creating porosity on its surface. Numerical methods such as the finite volume method and the Navier-Stokes equations as well as the Spalart-Allmaras model are used. By creating suction and blowing and increasing the cross-sectional area of the flow, porous surfaces reduce the effects of shock. Achieving around 50% porosity on the obstacle weakened the shock and reduced the Mach number after it, according to this paper's results. Also, the area under the skin friction diagram has decreased by about 38% as a result of less interference between the shock and the boundary layer and blowing caused by the porous surface.

1960s Europe saw a large number of residential neighborhoods built to house those migrating from the countryside. Today, more than 50 years later, these neighborhoods suffer high levels of functional, social, and technical obsolescence. In response to this, University of Seville developed the Aura Strategy as an intervention methodology to find global solutions to issues in outdated neighborhoods. To provide visibility to this aspect of the Aura strategy, the retrofit proposal presented in this article provides a solution to improve the roofing of buildings in a case study neighborhood (Polígono de San Pablo, Seville) and an analysis of the results in terms of energy and financial savings for local residents. The results show that for a population of roughly 18,000 (in 2018), net savings, including energy and retrofitting costs, ranging from nearly 6.5 to over 8.6 million can be made over the 20-year life cycle span. Likewise, the results obtained on the reduction of thermal loads indicate a 72% decrease in energy consumption, equivalent to a saving of close to 4,500 tons of greenhouse gas emissions for the district and the entire life cycle time period, with the consequent benefits on the impact on air quality and and the fight against climate change. climate change.

This article presents the impact analysis of space harmonics (SH) presence on a single-sided three-phase linear induction motor, along with a comprehensive parametric investigation. The presence of SH often reduces the LIM performances. The electromagnetic phenomena are governed by Maxwell's equations associated with equivalent circuit modelling. The chosen mathematical model uses 2D formulation with magnetic vector potential A. The model implementation is performed with finite element method on Gmsh-GetDP free platform. The electromagnetic thrust is calculated in current excitation case by a numerical model (FEM) with and without term generating SH in order to highlight the effect of SH. The adaptation to voltage supply operation is obtained by equivalent electric circuits (EEC), through calculation of the operational impedance. The main objective of the study consists in a comparative exploration of SH presence using two numerical models (single and multi-harmonics), and an assessment of parameters effect on space harmonic presence and various machine characteristics such as thrust, efficiency, and power factor. The motor’s characteristics (ie: thrust, efficiency and power factor) strongly depend on parameters such as the pole pair number and conductivity. Improving the operation and maximizing the performances of such a machine for a given specification requires the use of optimization algorithms.

The transportation of feed pellets via pneumatic conveying plays an improtant role in feed industry. In recent years, many efforts are made both exprimentally and numerically to gain a better understanding of the underlying mechanism for pneumatic conveying. Though few of those research is directly related to feed pellets, it can be beneficial to draw a full picture of rencent development in pneumatic conveying, with emphasis on aspects extendable to feed transportation. This paper reviews first the major factors concerning feed pellets pneumatic conveying, such as particle breakage, pipe wear and electrostatic. Building on that, new advancement including emerging new technologies in experiments and simulations are discussed, with a focus on the experimental measurements and mathematical models, repectively. The needs for future research and the analogy to feed pellets tranportation is given lastly.

Several factors impact the durability of concrete bridge decks, including traffic loads, fatigue, temperature changes, environmental stress, and maintenance activities. Detecting problems such as corrosion, delamination, or concrete degradation early on can lower maintenance costs. Nondestructive evaluation (NDE) techniques can detect these issues at early stages. Each NDE method, meanwhile, has limitations that reduce the accuracy of the assessment. In this study, multiple NDE technologies were combined with machine learning algorithms to improve the interpretation of half-cell potential (HCP) and electrical resistivity (ER) measurements. A parametric study was performed to analyze the influence of five parameters on HCP and ER measurements, such as degree of saturation, corrosion length, delamination depth, concrete cover, and moisture condition of delamination. The results were obtained through finite element simulations and used to build two machine learning algorithms, a classification algorithm and a regression algorithm, based on the Random Forest methodology. The algorithms were tested using data collected from a bridge deck in the BEAST® facility. Both machine learning algorithms were effective in improving the interpretation of ER and HCP measurements using data from multiple NDE technologies.

Understanding the dynamic interaction between piles and surrounding soil under vehicular impacts is essential for effectively designing and optimizing soil-embedded vehicle barrier systems. The complex behavior of pile-soil systems under impact loading, attributed to the soil’s nonlinear behavior and large deformation experienced by both components, presents significant simulation challenges. Popular computation techniques, such as the Updated Lagrangian Finite Element Method (UL-FEM), encounter difficulties in scenarios marked by large soil deformation, e.g., impacts involving rigid piles. While mesh-free particle and discrete element methods offer another option, their computational demands for field-scale pile-soil impact simulations are considerable. We introduce the erosion method to bridge this gap by integrating UL-FEM with an erosion algorithm designed for simulating large soil deformations during vehicular impacts. Validation against established physical impact tests confirmed the method’s effectiveness for flexible and rigid pile failure mechanisms. Additionally, this method was used to investigate the effects of soil mesh density, soil domain sizes, and boundary conditions on the dynamic impact response of pile-soil systems. Our findings provide guidelines for optimal soil domain size, mesh density, and boundary conditions. This investigation sets the stage for improved, computationally efficient techniques for the pile-soil impact problem, leading to better pile designs for vehicular impacts.

This study focuses on the formulation and analysis of a numerical method to compute singularly perturbed functional differential equations (SP-FDEs) involving large mixed shifts. SP-FDEs are a class of differential equations with both shifts and small perturbation parameter. Such equations arise in various scientific fields to model phenomena with memory effects and perturbed dynamics. In this particular study, the considered SP-FDEs involve large negative shift (delay) and large positive shift (advance). Dealing with large mixed shifts in SP-FDEs gives rise to challenges in terms of stability analysis, numerical methods, and convergence of a solution. In this paper, the influence of the large mixed shifts is treated by choosing a uniform mesh discretized in such a way that the term with the shifts lie on the nodal points. Then, using the Numerov’s finite difference method, the continuous problem is transformed into a finite difference scheme and the influence of the perturbation parameter is handled by introducing a fitting factor in the difference scheme. Stability estimate and convergence analysis of the developed scheme are investigated and proved. Numerical experiments are carried out to demonstrate the validity and applicability of the developed scheme. It is shown that the numerical scheme obtained in this work is parameter-robust of second order convergence rate.

Using numerical software to establish a three-dimensional geological model of shield tunneling through fault fracture zones, considering the synchronous grouting effect at the shield tail, the disturbance of various geotechnical structures in the fault fracture zone under deep burial and high water pressure conditions is studied. According to the comparison between the numerical model analysis results and the on-site inspection results of the entire process of Shiziyang B3 section, synchronous grouting at the shield tail is the main reason for the changes in soil and water pressure in the fault fracture zone area. The soil and water pressure changes strongly during shield tunneling, and the maximum excess pore pressure value occurs when the shield tail passes. When the shield tunneling machine is far away, the soil and water stress shows a downward trend. The impact of synchronous grouting decreases with the increase of tunneling distance. Due to the difference in permeability, the variation law of pore pressure in the fault fracture zone is significantly different from the original rock, leading to a faster increase in pore pressure in the fracture zone area and a wider range of influence.

The steel grid structure with bolted ball joints is gaining popularity in offshore platforms and aquaculture cages, but the corrosion resistance of these joints is a critical safety concern. This article introduces an innovative method to enhance the sealing performance of bolted ball joints, preventing seawater ingress and corrosion. The approach involves creating a sealed surface within the contact gap between the sleeve and connecting components. Circular grooves and sealing washers are added to both ends of the sleeve. A two-dimensional finite element analysis model of the bolted ball joint with the sealing structure was created using SOLIDWORKS and ANSYS. The study analyzes the sealing gasket's contact pressure at various compression levels and evaluates its performance with bubble tests for air tightness. Research results show a linear relationship between contact pressure and compression rate, with a sealing pressure of 5.95MPa achieved at a 14% compression rate. This simple, cost-effective sealing method effectively prevents corrosion from seawater, making it ideal for marine environments.

Bare slope instability is a prevalent concern. The root system of, herbaceous vegetation enhances the shear strength of shallow–slope soil. Indoor experiments were conducted on rootless undisturbed soil (RUS) and undisturbed soil with a root system (USRS) using a triaxial compression apparatus to analyze the slope stability of composite soil with a Tagetes erecta root system. Significance tests and correlation analysis of the factors affecting shear performance were conducted using R software. The slope reinforcement, effect by the plant root system was simulated under 24 working conditions using the MIDAS finite element method. The results revealed the influence of the–root content (RC), moisture content (MC), and stress on the shear strength of USRS, the contribution degree, and the variables’ influences on slope stability. Both RUS and USRS exhibited strain hardening during shearing. The internal friction angle (φ) and cohesion (c) of USRS were negatively and positively correlated with the RC and MC (root burial depth), respectively, and a good fit was obtained for the relationship. The maximum deviatoric stress during shear failure was 1.29 times higher for USRS than for RUS. The RC (root depth) was positively correlated with the slope safety coefficient and the slope of the line under different working conditions, whereas the slope gradient was negatively correlated with the slope safety coefficient. The reinforcement effect by the root system resulted in a 13.2% increase in the safety coefficient and improved stability of slopes with a gradient larger than 1.5%. This article investigated the mechanism of the root–soil system and, the effects of different influencing factors on the shear strength of the soil, and slope stability. The findings provide new insights into shallow slope stability in practical slope protection projects.

In the context of post-harvest pest management in agricultural products, the adoption of modified atmospheres presents an eco-friendly alternative to conventional pesticides. This study focuses on nitrogen gas as a potential agent for insect control in stored commodities, utilizing computational simulations to investigate its penetration and distribution within two common storage configurations: chamber-contained pallets and silos. By employing the convection-diffusion equation, the simulations reveal insights into nitrogen distribution dynamics and duration required for oxygen reduction. The results highlight the influence of boundary conditions, commodity porosity, and convection effects on nitrogen dispersion. In chamber scenarios, boundary condition type significantly affects oxygen reduction time, while in silos, convection and diffusion interact to establish a consistent diffusion layer thickness. This research enhances the understanding of modified atmosphere techniques and their potential for sustainable pest management. It provides practitioners with valuable quantitative insights, paving the way for optimized modified atmosphere strategies in real-world agricultural contexts.

The fatigue delamination damage is one of the most important fatigue failure modes of laminated composite structures. However, the developments of the theoretical framework, mathematical/physical model and numerical simulation of fatigue delamination are still challenging issues, and there is a lack of systematic classification. This article reviews the experimental phenomena of the onset and propagation stage of delamination under fatigue loading, and summarizes the commonly used phenomenological models for laminated composite structures. The research methods, general modeling formulas and development prospects of the phenomenological models are presented in detail. In the light of the analysis of finite element model (FEM) for the laminated composite structures, a number of simulation methods of fatigue delamination crack propagation damage model are carefully classified. Furthermore, within the theoretical framework of strength theory model (SM), fracture mechanics model (FM), damage mechanics model (DM) and hybrid model (HM), the complete procedure, application, capability assessment, advantages and limitations of the models are also discussed in detail. Finally, several valuable suggestions are given for the research interests, and challenging technical problems that need to be overcome on the fatigue delamination prediction of laminated composite structures in the future.

Typhoid fever is an infectious disease that affects humanity worldwide particularly in the lower socio-economic communities where many individuals are exposed to a dirty environment and unclean food. A mathematical model is formulated to analyze the impact of control measures such as vaccination of susceptible humans, treatment of infected humans and sanitation in the different socio-economic communities. The model assumed that the population comprises of two socio-economic classes. The essential dynamical system analysis of our model was appropriately carried out. The impact of the control measures was analyzed, and the optimal control theory was applied on the control model to explore the impact of the different control measures. Numerical simulation of the models and the optimal controls were carried out and results obtained indicates that the overall combination of the control measures eradicates Typhoid fever in the population but the controls are more optimal in the higher socio-economic communities.

In the process of construction and operation of the pumping station, a large number of axial flow pumps pose a threat to the survival of fish in the natural environment. Through experiments, the survival conditions of many kinds of fish in extreme positive and negative pressure environment are explored, and -40 kPa is established as the low pressure damage threshold, which is used as a reference to analyze the probability of fish damage caused by pumps. The immersion boundary-lattice Boltzmann numerical method (LES-IB-LBM) based on large eddy simulation is used to optimize the parameters of the prototype runner of axial flow pump. A new type of runner with leading edge thickness of 11.4mm and blade cutting angle of 18 ° is optimized. This kind of runner has both hydraulic performance and fish-passing ability. The runner of the fish-friendly axial flow pump effectively guides the fish into the low risk area and reduces the risk of friction, shear and impact damage to the fish. the total impact mortality ratio of the runner of the axial flow pump before and after optimization is close to 7:1. The pressure damage analysis shows that the lowest negative pressure of the fish is -59.2 kPa, which is higher than -64.5 kPa of the prototype pump, and the highest positive pressure of the body surface is 11.8 kPa, which is lower than the 43.8 kPa of the prototype pump. Shear damage analysis shows that the fish will not be damaged by shear in the prototype pump and fish-friendly pump.

Experimental methodologies for fatigue lifetime prediction are time-intensive and susceptible to environmental variables. While the cohesive zone model is popular for predicting adhesive fatigue lifetime, entropy-based methods have also shown promise. This study aims to provide an understanding of the durability characteristics of CFRP adhesive joints by incorporating an entropy damage model within the context of the finite element method (FEM), examined the effects of different adhesive layer thicknesses on SLS models. As the adhesive layer increased in thickness, DAMAGE VARIABLES initially rose and then declined, peaking at 0.3mm. This finding provides crucial understanding into the stress behavior at the resin-CFRP interface and the resin's fatigue mechanisms.

The effects of wastewater step-feeding on the performance of Horizontal Subsurface Flow Constructed Wetlands (HSF CWs) are numerically investigated. The purpose is to check if this alternative feeding technique increases the ability of HSF CWs to remove pollutants. Two methodologies are used: Initially, the Tanks-In-series (TIS) methodology, based on Finite Volume Method (FVM), is analyzed by using the volumetric degradation coefficient λ. In this case the operation of CW is similar with a series of Continuous Stirred-Tank Reactors (CSTR), under steady conditions. Then, the Step-Feeding (SF) procedure is presented, when the CW is operated like a Plug Flow Reactor (PFR). The results show that the SF procedure does not improve the performance of HSF CWs on removing Biochemical Oxygen Demand (BOD), operating under Mediterranean conditions.

The oscillation of an elastically mounted rigid cylinder excited by steady fluid flow is investigated. Regarding the nonlinearity of real practical structures like marine risers and the stay cable of a long-span bridge, the dynamic behavior of the circular cylinder is described by double Duffing equations and the aerodynamic force performance of the wake flow is expressed by wake oscillator equation. Unlike previous studies, in the present investigation, attention is focused on coupling the wake oscillator equations, taking into account quadratic terms. Following this approach, the cylinder's combined cross-flow and in-line Vortex-Induced Vibrations (VIV) are modeled more accurately. Empirical coefficients are adjusted through comparison with previous credible experimental studies and the effects of changing coefficients of the quadratic terms on major VIV parameters are investigated. The oscillating amplitude calculated by the present model is close to that of the experiment. The relative error of results obtained by the present coupled model is lower than the previous model. Moreover, the present model successfully predicts the moving trajectories of a circular cylinder under VIV in a figure-of-eight shape.

The investigation of plasma channels induced by focused ultra-short 1030-nm laser pulses in bulk of synthetic HPHT diamond revealed strong dependence of their size and shape on the used numerical aperture of the lens (NA=0.15-0.45). It was shown that at loose focusing conditions, it is possible to significantly increase the length of the plasma channel with a slight increase in pulse power, while tight focusing allows to obtain more compact structures in the same range of used powers. Such dependence paves the way to new possibilities in 3D processing of transparent dielectrics, allowing, for example, to vary the spatial parameters of modified regions without changing the setup, but only by controlling the lens aperture, which seems very promising for industrial applications.

This paper presents an original method of additive manufacturing of cylindrical parts with variable circumference thickness, which allows the control of the deposition of molten material using an algorithm for decomposing the part geometry into volumetric elements with known dimensional configuration. In the absence of a post-processor capable of controlling additive manufacturing on a 5-axis numerical control machine, control of the deposition of molten material is done using parameterized programs, which can control both the feed speeds of the machine tool axes and the specific functions of the printing equipment. Additive manufacturing can make a positive contribution to sustainable development compared to traditional manufacturing technologies, thus making a positive contribution for a sustainable future. The aim of the work is to make it possible to 3D print parts with variable wall thickness using a CNC machining centre. To obtain the variable thickness layer we have implemented an original method of deposition with molten material (FDM), the coordination of the system composed of physical elements respectively programmable elements, is realized through control functions materialized in parameterized part programs, the generated outputs being the variable speed of the machine axes on a circular trajectory, the angular positioning, the filament advance.

Bistatic scattering coefficient have important application of ocean active microwave remote sensing. In this paper, a general expression of the composite scattering coefficient for multiple dielectric targets of one-dimensional fractal sea surface is derived by MOM method. Secondly, the bistatic scattering coefficient is computed, there are two class missile targets and one class ship target of one-dimensional fractal sea surface, and the influence of the water dielectric constant, fractal dimension, wind speed, the target dielectric constant and spatial position is discussed, and the corresponding conclusions are given. The proposed algorithm is able to solve some composite scattering problems, such as “dielectric or conductive object + rough surface”, “dielectric or conductive floating object +rough ”and it has certain generalization and reference.

Bending pipe is a common component of long distance pipeline. Accurately understanding the flow law of hydrate particles in elbow pipe is of great significance for optimizing pipeline design, improving production efficiency of gas transmission pipeline and ensuring pipeline safety. Taking the flow of hydrate particles in a Bending pipe as the object of study, the spiral flow attenuation and hydrate particle deposition under different bending pipe angles, different rate of bending pipe to diameter, different Reynolds number and different torsion were studied by numerical simulation method. The results show that the larger the swirl number, the higher the spiral flow intensity. When the fluid enters the bending pipe, the Angle of the bending pipe is larger, the torsional rate is smaller, and the Reynolds number is larger, the swirl number is larger, the intensity of the swirl flow is stronger, and the swirl number is larger at the same position. However, the larger the bending pipe to diameter ratio, the greater the change of the total swirl number and the weaker the strength at outflow, but the slower the change of swirl number, indicating that the spiral flow attenuation is slower. The larger the bending pipe Angle, the larger the bending pipe diameter rate, the smaller the torsion rate and the larger the Reynolds number, the smaller the deposition rate of the particles after flowing through the pipe. The results show that the spiral flow strength can be maintained in a proper way to reduce the deposition and ensure the safety of the conveying process.

In this study, three new demountable joints consisting of reinforced concrete columns and steel beams were proposed and their seismic performance was investigated using cyclic loading tests. The test results demonstrated that the three demountable RCS joints had good seismic performance. Using the finite element software ABAQUS, the influence of parameters such as beam flange thickness, bolt strength, and connector steel strength on the seismic performance of each joint was analyzed, and the influence of different parameters on the seismic behavior of the joint was determined. The results also showed that the three demountable RCS joints were sensitive to changes in the thickness of the steel beam flange, while the connector steel strength and bolt type had little impact on the joint's capacity. Additionally, shear capacity calculation formulas for the joint core area provided by different codes and researchers were compared with test and finite element results. The results showed that the calculation results of the ASCE guidance method, the Nishiyama method, and the CECS 347-2013 method were higher than the test values, while the calculation results of the Para method were lower because it did not consider the contribution of the cylindrical steel plate to the shear capacity of the joint core area.

The educational and play-related activities of children are proceeded mainly indoors in a kin-dergarten. The exposure of children to high concentrations of PM2.5 and CO2 indoors have leaded to various harmful effects and negatively impact educational outcomes in a kindergarten. Therefore, this study presents a numerical model of CO2 and PM2.5 concentration based on two classrooms in a kindergarten. Using this numerical model, we present different scenarios of oper-ating mechanical ventilation and air purifiers in kindergartens with the aim of optimizing for reducing the concentrations of CO2 and PM2.5. We found that the amount of ventilation required to maintain good air quality, per child, was approximately 20.4 m3/h. However, we also found that as the amount of ventilation increased, so did the concentration of indoor PM2.5; we found that this issue can be resolved through the use of a high-grade filter (i.e., a MERV 13 grade filter with a collection efficiency of 75%). This study provides a scientific basis for reducing PM2.5concentrations in kindergartens, while keeping CO2 levels low.

With more and more transportation tunnels have been and will be constructed in loess areas in Northwest China with high earthquake potential, the overall stability of portal section under earthquakes actions and the related aseismic countermeasures attracted the attentions from both the scholars and engineers, especially the tunnels in the upper slope connecting the high bridges crossing rivers or valleys. In order to study the dynamic response characteristics and damage evolution of steep loess slope with a tunnel under the earthquake actions, a large-scale shaking table tests were performed on steep loess slope with a tunnel. Wenchuan-tangyu (WT) wave and El Centro (El) wave records were applied on the model to investigate the displacement response and acceleration response of loess slope with a tunnel under the horizontal (X) and the combined action of horizontal and vertical (X-Z) seismic loads, respectively. In particular, three-dimensional non-contact optical measurement techniques were used to obtain the slope surface displacements. The results showed that the main deformation pattern of the slope was horizontal movement and settlement when the seismic wave input was in the X and X-Z directions, respectively. However, the X direction seismic wave had a greater impact on the deformation of the slope, and the tunnel portal slope was destroyed under the action of a large horizontal seismic acceleration finally. Slope failure ahead of a tunnel could be divided into four stages, i.e. elastic deformation stage, plastic deformation accumulation stage, local failure stage, and overall failure stage. The peak ground displacement of X direction (PGDX) and peak ground displacement of Z direction (PGDZ) of the slope surface increased with the increasing of the input peak ground acceleration (PGA) when the input wave was same waveform and same direction. The existence of the tunnel had a great influence on the PGA and the PGA amplification factor (PGAAF) of the soil mass surrounding it. This was because the seismic waves encounter a tunnel surface with clear differences in the physical properties of the medium during their propagation in the slope, thereby forming a strong reflection and refraction effect, and the amplitude changed significantly. Numerical simulation results were basically consistent with the experimental results.

It is imperative to examine the performance of wind turbines both in steady and unsteady states due to their significant role in generating renewable energy by utilizing wind power. The computational fluid dynamics method (CFD) is one of the most efficient tools in this field, and the unsteady state has a great deal of importance when predicting turbine performance. Using the CFD method, the hydrodynamic performance of two types of offshore and onshore horizontal-axis wind turbines in the unsteady state was examined. Offshore wind turbines were modeled using the SST turbulence model and onshore wind turbines were modeled using the large eddy simulation (LES). Three-dimensional results of the simulation were compared and validated with experimental results, and good agreement between the two results indicates that the assumptions and method were accurate. Vortices were also studied in relation to the blade tip and their formation and development. Increasing the number of blades results in a reduction in torque, as well as a rise in vortices caused by the tips of the rotor blades. Due to the intensity of turbulence, the distribution of pressure on the surfaces will become irregular as the wind speed increases.

Controlling injection parameters is paramount when it comes to producing high-quality green parts using powder injection molding. This work combines experimental and numerical approaches to study the impact of injection parameters on mold in-cavity pressure and on the overall quality of green parts produced by low-pressure powder injection molding. The properties of two low-viscosity feedstocks (formulated from a water-atomized stainless-steel powder and wax-based binder system) were measured and implemented in an Autodesk Moldflow numerical model to quantify the molding pressures, which were finally validated using experimental real-scale injections. The results confirmed that an increase in mold temperature, an increase in feedstock temperature, and a decrease in solid loading decrease the mold in-cavity pressure, which was correlated with the feedstock viscosity. As a key result, real-scale injections confirmed that a minimum flow rate was required to avoid atypical high in-cavity pressure leading in several visual defects such as weld lines, flow marks, cracks, sinks, and incomplete filling. Due to differences in its thermal transfer properties, this flow rate threshold value decreases as the feedstock solid loading increases. For injection speeds higher than this value, the injection pressure measured experimentally was linearly correlated with the injection flow rate.

. The impedance change of an induction coil surrounding a metal tube adapter is investigated by the truncated region eigenfunction expansion (TREE) method. Conventional TREE method is inapplicable to this problem as a consequence of the numerical overflow of eigenfunctions of the air-metal multi-subdomain regions. The difficulty is surmounted by a normalization procedure for the numerical eigenfunctions obtained from the finite element method (FEM). Efficient algorithm is devised by the Clenshaw-Curtis quadrature rule for integrals involving the numerical eigenfunctions. Numerical results of the TREE and FEM simulation coincide very well in all cases, and the efficiency of the proposed method is also confirmed.

This work focuses on the research area of energy generation and harvesting for its transformation into electrical energy with the primary use of energy independent sensor systems in the power range P=10-10000W. These systems are applied in sensing quantities in the transport sector (bridge structures, transport infrastructure, and the others ). The proposed theoretical harvester models describing the transformation of motion energy into electrical energy, also provides theoretical and experimental models. Results obtained in the design and construction of a robust motion generator with primarily linear geometry system technology are presented. The expected output electrical power of the harvester is variable and is designed based on the linear motion generated by the engine- compressed air, small fuel system, etc. The fundamental design of the generator core has been continuously numerically modeled and an experimental setup has been developed to analyze specific parts and variations in order to validate the concept and achieve the most suitable parameters with the selected construction materials. A scaled down version of the model principle was tested in experiments and then the parameters and results were compared with the predicted theoretical analyses.

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.

Understanding the infiltration and solidification processes of liquid 5083Al alloy into Al2O3 three-dimensional reticulated porous ceramic (Al2O3(3D) RPC) is essential for optimizing the microstructure and properties of Al2O3(3D)/5083Al interpenetrating phase composites (IPCs) prepared by low-pressure infiltration process (LPIP). This study employs ProCAST software to simulate the infiltration and solidification processes of liquid 5083Al with pouring velocities (PV) of 0.4 m/s infiltrating into Al2O3(3D) RPC preforms with varying porosities at different pouring temperatures (PT) to prepare Al2O3(3D)/5083Al IPCs using LPIP. The results demonstrate that pore diameter of Al2O3(3D) RPC preforms and PT of liquid 5083Al significantly influence the of the infiltration. Solidification process analysis reveals that the Al2O3(3D) RPC preform with smaller pore diameters allows the lower pouring velocity of 5083Al to solidify faster compared to the preform with larger pore diameters. Al2O3(3D)/5083Al IPCs were prepared successfully from Al2O3(3D) RPC porosity of 15 PPI with liquid 5083Al at PV 0.4 m/s and PT 800 °C using LPIP, resulting in nearly fully dense composites, where both Al2O3(3D) RPCs and 5083Al interpenetrate throughout the microstructure. The infiltration and solidification defects were reduced under air pressure of 0.3 MPa (corresponding to PV of 0.4 m/s) during LPIP. Finite volume method simulations are in good agreement with experimental data, validating the suitability of the simplified model for Al2O3(3D) RPCs in the infiltration simulation.

We measure the accelerated cosmic expansion rate of the late-time Universe and derive constraints on cosmological parameters from the final baryon acoustic oscillation (BAO) uncorrelated measurements of different surveys after two decades of dedicated spectroscopic observation combined with the Pantheon sample of Type Ia supernovae, the Hubble diagram of gamma-ray bursts & quasars data, the Cosmic Chronometers data, and the recent Hubble constant value measurement from the Hubble Space Telescope and the SH0ES Team as an additional prior. In Λ cold dark matter (ΛCDM) scenario, the model fit yields Ωm=0.313±0.034 and ΩΛ=0.672±0.025. Combining BAO with CC+Pantheon+QSR+GRB data sets we get H0=69.21±1.22 km s−1 Mpc−1, rd=133.46±2.49 Mpc. For the a flat wCDM model, we get w=−1.111±0.030. For a dynamical evolution of dark energy in a flat w0waCDM cosmology, we get wa=−0.292±0.167. We apply the Akaike information criteria probe to compare the three models, and see that all models are favored and cannot be discarded from the current data.

This paper is concerned with developing an efficient numerical algorithm for fast implementation of the sparse grid method for computing the $d$-dimensional integral of a given function. The new algorithm, called the MDI-SG ({\em multilevel dimension iteration sparse grid}) method, implements the sparse grid method based on a dimension iteration/reduction procedure, it does not need to store the integration points, nor does it compute the function values independently at each integration point, instead, it reuses the computation for function evaluations as much as possible by performing the function evaluations at all integration points in a cluster and iteratively along coordinate directions. It is shown numerically that the computational complexity (in terms of CPU time) of the proposed MDI-SG method is of polynomial order $O(Nd^3 )$ or better, compared to the exponential order $O(N(\log N)^{d-1})$ for the standard sparse grid method, where $N$ denotes the maximum number of integration points in each coordinate direction. As a result, the proposed MDI-SG method effectively circumvents the curse of dimensionality suffered by the standard sparse grid method for high-dimensional numerical integration.

The practical necessity of strengthening of reinforced concrete structure elements requires a choice of a method which would provide cost-effectiveness, simplicity of construction and durability of the chosen solution. One of the ways is strengthening using the material belonging to the group of micro-reinforced concretes called ferrocement. The paper presents the comparative analysis of some of the results of RC beams exposed mainly to flexure, strengthened with ferrocement elements (strips) applied by gluing on the on the tensed side of RC beams. The results are obtained using the experimental, analytical and numerical research on the adequate models. Strengthening of the beams was performed in the form of four types of ferrocement strips (four different strip widths and four different numbers of wire mesh layers). The numerical non-linear analysis using finite elements method (FEM) was employed, with the introduction of the corresponding characteristics of constitutive material obtained by the experiments on the test specimens. The analytical researches were performed using known analytical calculation methods in order to formulate the model for the calculation of the ultimate and serviceability limit states of the strengthened cross-section. The results presented in the paper, through their comparative compatibility, confirmed the applicability of this strengthening method.