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.

A galvalume corrugated sheet was utilized as formwork for a reinforced concrete beam in flexure. A numerical model was validated to the experimentally obtained data, and further adopted to simulate the behavior of this composite structure under elevated temperatures. The properties and constitutive stress-strain data of the basic materials were obtained from experiments, and superimposed into the finite element model. The study concluded that the load carrying capacity of the member decreased as a direct function on temperature increase, and the cracking moment was very sensitive to the temperature fluctuation. The elevated temperatures also altered the failure mode.

This paper presents the results of aerodynamical performance research focused on maintaining the thermal comfort and increasing the energy efficiency of a typical social housing unit located in a high-density urban area. Bioclimatic design strategies are used to develop a sustainable and economic technology in existing housing clusters in Mexico City. A full-scale prototype, built on campus facilities, was used to study the flow conditions around the design. All scaled prototypes implement similar criterion. Furthermore, a scaled prototype is evaluated within a low speed wind tunnel installation. Additionally, numerical simulations were performed at transient state based on previous physical measurements and historical local climatic conditions to determine preferable modifications.

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.

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.

It is difficult to obtain necessary information accurately from Social Networking Service (SNS) while raising children, and it is thought that there is a certain demand for the development of a system that presents appropriate information to users according to the child's developmental stage. There are still few examples of research on knowledge extraction that focuses on childcare. This research aims to develop a system that extracts and presents useful knowledge for people who are actually raising children, using texts about childcare posted on Twitter. In many systems numbers in text data are just strings like words and are normalized to zero or simply ignored. In this paper, we created a set of tweet texts and a set of profiles created according to the developmental stages of infants from "0-year-old child" to "6-year-old child". For each set, we used ML algorithms such as NB (Naive Bayes), LR (Logistic Regression), ANN (Approximate Nearest Neighbor algorithms search), XGboost, RF (random forest), decision trees, and SVM (Support Vector Machine) to compare with BERT, a neural language model, to construct a classification model that predicts numbers from "0" to "6" from sentences. The accuracy rate predicted by the BERT (Bidirectional Encoder Representations from Transformers) classifier was slightly higher than that of the NB, LR, ANN, XGboost, RF, decision trees, and SVM classifiers, indicating that the BERT classification method was better.

Since the discovery of graphene, due to its excellent optical, thermal, mechanical and electrical properties, it has a broad application prospect in energy, materials, biomedicine, electromagnetism and other fields. A great quantity of researches on the physical mechanism of graphene has been applied to engineering in electromagnetism and optics. To study the properties of graphene, different kinds of numerical methods have been developed for simulating the electromagnetic field effects of graphene, and equivalent boundary conditions have been employed to replace graphene in the computational domain. In this work, the numerical methods with equivalent boundary conditions are reviewed, and some examples are provided to illustrate their applicability.

It is difficult to obtain necessary information accurately from Social Networking Service (SNS) while raising children, and it is thought that there is a certain demand for the development of a system that presents appropriate information to users according to the child's developmental stage. There are still few examples of research on knowledge extraction that focuses on childcare. This research aims to develop a system that extracts and presents useful knowledge for people who are actually raising children, using texts about childcare posted on Twitter. In many systems numbers in text data are just strings like words and are normalized to zero or simply ignored. In this paper, we created a set of tweet texts and a set of profiles created according to the developmental stages of infants from "0-year-old child" to "6-year-old child". For each set, we used Support Vector Machine (SVM), and Bidirectional Encoder Representations from Transformers (BERT), a neural language model, to construct a classification model that predicts numbers from "0" to "6" from sentences. The accuracy rate predicted by the BERT classifier was slightly higher than that of the SVM classifier, indicating that the BERT classification method was better.

Herein we foremost detailed the numerical modeling of the double absorber layer- methyl ammonium lead iodide– carbon nitride layer solar cell and subsequently provided in-depth insight on the active layer associated with dominant radiative and non-radiative recombination losses limiting the efficiency ( ) of the solar cell. Under recombination kinetics phenomena, we explored the influence of Radiative recombination, Auger recombination, Shockley Read Hall recombination, the energy distribution of defects; Band Tail recombination (Hoping Model), Gaussian distribution, metastable defect states including single donor (0/+), single acceptor (-/0), Double Donor (0/+/2+), double acceptor (2/-/0-), and the interface layer defects on the output characteristics of the solar cell. Setting defect (or trap) density to with uniform energy distribution of defects for all the layers, we achieved the of 24. 16 %. A considerable enhancement in power conversion efficiency was perceived as we reduced the trap density to for the absorber layers. Further, it was observed that for the absorber layer with double donor defect states, the active layer should be carefully synthesized to reduce crystal order defects to keep the total defect density as low as to achieve efficient device characteristics

In this study, we have considered for numerical solution of a Poisson equation in a domain. we have presented a first and second order finite difference method for solving the system of the intial boundary value. We have presented the numerical solution of the Poisson equations in a two-dimensional finite region. Numerical Solution and exact solution for the presented method has been compared and the figures presented. The present error analysis tables demonstrate the efficiency of the method.

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

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

This code provides computational facilities to simulate current versus time during the charging of Li-ion cells at desire constant voltage by considering multiscale physical phenomena. This code only considers a powder of active materials (at microscale or nanoscale) and a small part of electrolyte around it as a half cell. Then it is extended to a complete cell by applying correct boundary conditions. This code is very useful by modifying code parameters to understand the effect of the complex shape of active materials powder (surface area and powder size), kind of electrolyte, and the applied voltages on the charging response of Li-ion cell. As a summary, a microscale approach to the design of Li-ion cells has been provided via this code.

There are many criterion to generalize the concept of numerical radius; one of the most recent interesting generalization is what so called the generalized Euclidean operator radius. Simply, it is the numerical radius of multivariable operators. In this work, several new inequalities, refinements and generalizations are established for this kind of numerical radius.

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

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

Nonlinear thermal transport of non-Newtonian polymer flows is an increasingly important area in materials engineering. Motivated by new developments in this area which entail more re-fined and mathematical frameworks, the present analysis investigates the boundary layer ap-proximation and heat transfer persuaded by a symmetrical cylindrical surface positioned hori-zontally. To simulate thermal relaxation impacts the bioconvection cross nanofluid flow Buongiorno model is deployed. The impact of the magnetic field applied on the nanofluid using the heat generation and melting phenomena are inspected. The nonlinear effect of thermosolutal buoyant forces is incorporated into the proposed model. The novel mathematical equations comprise thermophoresis and Brownian diffusion effects. Via robust transformation techniques, the primitive resulting partial equations for momentum, energy, concentration, and motile living microorganism are rendered into nonlinear ordinary equations with convective boundary postulates. An explicit and efficient numerical solver procedure in Mathematica 11.0 programming platform is developed to engage the nonlinear equations. The consequence of multiple governing parameters on dimensionless fluid profiles is inspected through plotted visuals and tables. Finally, outcomes declared on the surface drag force, heat, and mass transfer coefficients are exhibited through 3D visuals via different influential parameters.

Taking Shanxi Province, a typical energy region in China, as the research object, data on economic development quality and environmental pollution intensity from 2001 to 2021 were selected. Linear regression, numerical simulation, and Pearson correlation coefficient were used to analyze the evolution characteristics of economic development quality and environmental pollution intensity in Shanxi Province in the past 20 years, and to explore the relationship between economic development and environmental pollution. The results indicate that since 2001, Shanxi Province has made long-term progress in economic development, with GDP increasing by nearly 10 times and an average annual growth rate of around 7%. In the past 20 years, the main pollutants have shown a trend of increasing first and then decreasing steadily, with a turning point occurring during the "12th Five Year Plan" period. This shows that the environmental policies and investments of the Chinese and Shanxi provincial governments in the past 10 years have been effective in the new era. The numerical simulation curve results show that the per capita GDP exhibits a classic inverted "U" curve relationship with wastewater and SO2 emissions, with a turning point occurring at around 20000 yuan per capita GDP; However, there is a monotonic decreasing trend with Chemical oxygen demand and ammonia nitrogen emissions, and a monotonic increasing trend with solid waste generation. There is no turning point yet. The correlation analysis results further support the fitting curve conclusion.

The stability of the filling roof as an important bearing unit in the stope of the access stope in the filling mining method is of great significance to guarantee the safe and efficient production of the mine. Arrangement of mining approach in downward cemented filling stope is the key factor af-fecting stability of filling body roof. Based on combination of laboratory test, theoretical analysis and numerical simulation, the influence of different mining approach arrangement on stability of filling body roof is analyzed. The weak filling surface is formed between adjacent mining paths. The mechanical strength of the weak filling surface is significantly reduced by laboratory experiments. The relationship between distribution of the weak filling surface, azimuth angle and stability of filling roof is further studied by numerical simulation. The results show that the stability of the filler roof is the best when the upper and lower stratified mining approach is arranged vertically.

In this paper, the research progress on dynamic behaviors of RC structural members is reviewed from experimental, theoretical and numerical perspectives. Firstly, the basic overview, measurement methods, main conclusions and current limitations of available dynamic loading tests are presented. Then the theoretical studies on the dynamic constitutive models of RC materials, the dynamic increase factor (DIF) model for concrete and reinforcing steel and the proposed modified models of dynamic behavior parameter at the structural member level are summarized. Finally, the available modelling approach and the method for incorporating the dynamic effect in numerical simulation of RC structures are reviewed. Moreover, a brief introduction is made to the dynamic hysteretic model established on the experimental data, which provides an alternative approach to consider the dynamic effect beside the commonly used DIF method. This paper provides valuable reference for experimental studies and numerical simulations on the dynamic behaviors of RC structures, and also puts forward some issues that need to be solved in the future research works..

In the construction of shield crossing existing mined tunnel without load, it is imperative to develop corresponding design standards that reflect actual engineering force characteristics to ensure the successful completion of tunnel construction. This research uses The MIDAS-GTS finite element software to facilitate the creation of a numerical model of the shield structure for an air-push-over mine tunnel project in Changsha, China, investigating the stress field’s evolution during shield construction, and calculating the maximum positive and negative bending moments and maximum axial forces for different structures and other force states under various construction conditions. This study's findings informed the design and construction optimisation of the shield tunnelling empty-push method. The outcomes of this numerical simulation led to several key findings: (1) The soil density exerted a significantly greater impact on the internal forces of the initial support structure than both the tunnel depth and soil Poisson's ratio. Additionally, a sudden shift in internal forces occurred within the 300-350 mm range when the lining thickness was altered. (2) Factors such as the tunnel depth, soil density, soil Poisson's ratio μ, and lining thickness have a similarly influenced internal forces of the segment and the initial support. Notably, the backfill layer thickness significantly affected the segment’s maximum axial force causing an abrupt change of approximately 300 mm. (3) Under the shield machine equipment’s weight constraint, it is essential to control the guide rail’s thickness, to prevent it from becoming overly large.

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 dynamic development of electromobility poses challenges to designers regarding not only the efficiency of energy transformation but also the battery life, which is influenced by the stability of its operating temperature. Designing cooling systems is connected not only with the optimization of energy management but also with other environmental parameters, such as noise emission. The paper presents the numerical optimization of an innovative radiator for use in electric buses in terms of energy consumption and noise emission. The results of the numerical studies were verified in laboratory and field conditions, showing a very good convergence of the model with the results of the experiments.

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.

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.

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.

Pore fluid pressure is very important to the generation, migration and accumulation of petroleum. The previous studies indicate that the Chezhen sag region is typically characterized by pore fluid overpressure (the difference between pore fluid pressure and hydrostatic pressure). To date, the formation mechanism of pore overpressure and the accumulation regularity of “upper source-lower reservoir” type in this region remains unknown. In order to investigate these problems, a 2-D finite element model is established with regards to the fluid-solid coupling in Chegu 25 block of Chezhen depression, fully considering the existing seismic, logging data and regional tectonic stress environment. Meanwhile, we calculate the abnormal overpressure generated at the source rock during hydrocarbon generation and the process of hydrocarbon migration and accumulation along the faults, analyzing the dynamic conditions of hydrocarbon downward accumulation. The results show that overpressure can accelerate migration of hydrocarbon and improve the efficiency of hydrocarbon accumulation. However, when the overpressure is too large, tensile fractures and shear fractures may occur, resulting in hydrocarbon dissipation, and changing the results of oil and gas accumulation. The overpressure at the source rock is mainly caused by hydrocarbon generation, while the overpressure at the reservoir is primarily created by unbalanced compaction. As the dominant channel of hydrocarbon migration, overpressure will change the direction and path of hydrocarbon migration in the fault. Therefore, the “upper source-lower reservoir” hydrocarbon accumulation model is strongly explained by the high permeability of faults and the presence overpressure. In summary, the simulated overpressure results are in good agreement with the mud weight equivalent overpressure and the drill stem tests (DSTs).

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

During the recent decades, global warming by greenhouse gas evolution has attracted worldwide attention and evermore strict regulations thereon have become institutionalized as international policies. In the process, more environment-friendly power generation technologies have been developed utilizing fossil fuels with a view to timely commercialization. As one of such “clean coal” technology, Integrated Gasification Fuel Cell system is a promising power generation means and COS hydrolysis reactor is installed downstream of coal syngas to remove acidic gas constituents such as H2S and COS. The most significant design parameters affecting performance of the COS hydrolysis reactor were selected to be GHSV, (catalytic) reaction temperature, and length ratio and numerical modeling was performed considering heat and fluid flow transfer as well as chemical reaction kinetics. Effect of the selected design parameters on the variation of conversion rate and reactant gas mixture concentration were comprehensively investigated to predict performance of COS hydrolysis reactor. Stochastic modeling of reactor performance was finally performed using Monte Carlo simulation and linear regression fitting.

This study investigates the impact of thin flame combustor design on hydrogen flame characteristics and combustion performance through numerical simulations. Variations in flame shape and combustibility between 100% methane and 100% hydrogen combustion are analyzed. Three combustor header shapes (flat, concave, and convex) are modeled to assess their influence on flame behavior. Results show distinct flow patterns, with concave headers promoting strong central flows and convex headers dispersing the flow outward. Temperature field analysis indicates that hydrogen flames have higher temperatures and shorter quenching distances compared to methane flames. A comparative analysis of combustion products is conducted to evaluate combustion performance and NOx emissions. Consequently, the concave header has high combustibility and increases temperatures and NOx fraction in hydrogen combustion.

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.

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

To solve the problem of deformation and stress in segments when a full-face hard rock tunnel boring machine (TBM) tunnel for an urban subway passes through existing buildings, a three-dimensional mechanical calculation model of TBM tunneling under existing buildings was established based on numerical analysis to comprehensively study the spatial attributes of tunnel segments during TBM tunneling. The Qingdao Metro Line 1 TBM tunnel in the hard rock stratum of Xiaobei District passes through existing buildings. The results indicated that (1) the stress and deformation of the segments were different at different positions. The horizontal and vertical stress and displacement of the segments in front of the building were significantly affected by the building. When the segment was not in the area below the horizontal plane of the building, the vertical stress of the segment was significantly greater than the two stresses in the horizontal direction. (2) The vertical displacements of the vault and the arch bottom were the largest, and the influence of vertical accessories was stronger when the segment was closer to the building.(3) When the tunnel was excavated to approximately one tunnel diameter in front of the building, the building was affected. If the surrounding rock is unstable, the excavation face may collapse and the ground surface may sink, possibly leading to inclination, cracking, and uneven settlement of the building.

The mountains in northern Xinjiang of China were studied during the snowmelt season. Multi-source fusion live data Chinese Land Data Assimilation System (CLDAS,0.05°× 0.05°,hourly data) were used as real data, and Central Meteorological Observatory guidance forecast (SCMOC,0.05°× 0.05 °,forecasting the following 10 days in 3 h intervals) was used as forecast data, both issued by the China Meteorological Administration. The dynamic linear regression and the average filter correction algorithms were selected to revise the original forecast products SCMOC.Based on conventional temperature forecast information,we designed four temperature rise prediction algorithms for essential factors affecting snowmelt.Temperature rise prediction algorithms included the daily maximum temperature algorithm, daily temperature rise range algorithm, snowmelt temperature algorithm, and daily snowmelt duration algorithm.Four temperature rise prediction algorithms values were calculated of each prediction product.Root-mean-square error algorithm and temperature prediction accuracy algorithm were used to compare and test each prediction algorithm value from the time sequence and spatial distribution.Comprehensive tests show that the forecast product revised by the average filter algorithm was superior to revised by the dynamic linear regression algorithm as well as the original forecast product.Through these algorithms, the more suitable temperature rise forecast value for each grid point in the study area can be obtained at different prediction times. The comprehensive and accurate temperature forecast value in mountainous snowmelt season can provide an accurate theoretical basis for effective prediction of the runoff in snowmelt areas and prevention of snowmelt flooding.

On the aero-engine combustor liner, a phenomenon of overheating resembling a hot spot exists, accompanied by a preferential peeling of the YSZ thermal barrier coating that will negatively affect the service life of the component. The temperature gradient of the hot spot may affect the YSZ fracture pattern.In this paper, the morphological and property changes of YSZ sintering is investigated experimentally, and a coupled analysis of the YSZ Crack Propagation Mode under the hot spot is performed using ABAQUS. The results show that the pore fractal size D of YSZ increases by 15%. Inside the hot spot region, the coating has a shear stress of 75MPa. By inserting Cohesive units globally in the model to simulate the random crack initiation and propagation, unlike the thinning of YSZ layered peeling caused by uniform superheating, hot spot will cause the crack initiation at the tip of the pores inside the coating and the oblique propagation, eventually forming an oblique crack connection through the pores. When the temperature gradient reaches 30K/mm, the crack propagation is 40% greater than in a uniform temperature field; consequently, the YSZ in the hot spot boundary region with a greater temperature gradient results in more severe bulk peeling.

Hydrogel microneedles are a promising technology for the delivery of different types of medicines locally and painlessly, as well as ISF extraction. As the hydrogel microneedles are inserted into the tissue, they swell and release drugs. To improve the effectiveness of this technology in delivering medicine at controlled and desirable doses and intervals, a deep understanding of the mechanism of drug delivery inside the microneedles is required. In this work, drug diffusion inside a tapered microneedle is investigated using numerical simulation. The microneedle is divided into many small elements, and the mass transfer equation of meloxicam is solved in each element over time. The skin is simulated as the sink in the microneedle surface for drug absorption. Simulations are performed for different sizes of microneedles. For a microneedle with a height of 500 µm and a base diameter of 250 µm, the drug completely penetrates the skin within 3.2 seconds. The rate of drug diffusion from the tip of the microneedle is higher than diffusion from the side area near the microneedle base. The obtained data demonstrate that in addition to the height and the base diameter, the microneedle’s aspect ratio, h/d, also affects the time of drug diffusion. We present a nonlinear equation to predict the time of complete drug diffusion as a function of the microneedle geometrical parameter, including the height and base diameter. The proposed equation calculates the total drug diffusion time with an error of less than 7% for all studied cases. Predicting drug diffusion patterns inside microneedles can be helpful in the biomedical field, especially in the drug-controlled release system for the optimization of drug delivery.

Projected future sea level rise and marine storminess is a serious threat for beaches as they induce beach flooding and erosion. Among other factors, wave runup play an important role in beach evolution and must be robustly assessed. However, little attention has been paid to the uncertainties associated to its characterization and how do they compare to other sources of uncertainty. As a first step towards this goal, we have quantified the impact of several sources of error in the estimation of wave runup on sandy beaches. Understanding what factors are more influential in the quality of the results will help to determine the main sources of uncertainty in beach flooding projections. To reach that goal, a calibrated state-of-the-art numerical modelling system has been setup for two beaches in the Mallorca islands (NW Mediterranean). The system has been forced with the best available information of nearshore incoming waves and has been validated against observations to define the benchmark quality. To determine the key factors affecting the quality of the system’s results, different systems configurations have been tested with different degrees of complexity.Our results show that using the most sophisticated modelling system with the best information on boundary conditions, bottom bathymetry, and submerged vegetation leads to a swash RMSE comparable to the standard deviation of the observed swash. We have also found that the choice of lateral boundary conditions (i.e., source of information for the incoming waves) can double the RMSE and induce large biases. Our results also show that using a simple empirical approach usually underestimates the wave runup. However, in locations with vegetated seabed there is a compensation error, and the empirical approach can lead to acceptable results if forced by nearshore waves. In addition, we have compared the error estimates obtained with uncertainties associated to projected sea level rise. Our findings suggest that the uncertainty associated with wave runup modelling should be considered in the assessment of the total uncertainty of future beach flooding, since the analysis performed indicates that this uncertainty can account for between 16% to 60% of the uncertainties linked to mean sea level projections.

Dynamic Realistic Failure Process Analysis (DRFPA) is a numerical software based on FEM it is firstly validated by comparing with the experimental test. Such validation is carried out under two consistent concrete specimens under uniaxial compression. After validation, DRFPA is used to study the influence of two parallel pre-existing crack inclination angles on the fracture of concrete under dynamic loads. The evolution of wing cracks and secondary cracks is investigated by observing the initiation, coalescence, and propagation of them. Meanwhile, the relationship between failure mode and inclination angles is discussed via analyzing the acoustic emission events. Finally, it is concluded the relationship between inclination angles of pre-cracks and performance/stability of concrete structure.

For the numerical simulation of components made of short fiber-reinforced composites the correct prediction of the deformation including the elastic and plastic behavior and its spatial distribution is essential. When using purely deterministic modeling approaches the information of the probabilistic microstructure is not included in the simulation process. One possible approach for the integration of stochastic information is the use of random fields. In this study numerical simulations of tensile test specimens are conducted utilizing a finite deformation elastic-ideal plastic material model. A selection of the material parameters covering the elastic and plastic domain are represented by cross-correlated second-order Gaussian random fields to incorporate the probabilistic nature of the material parameters. To validate the modeling approach tensile tests until failure are carried out experimentally, that confirm the assumption of spatially distributed material behavior in both the elastic and plastic domain. Since the correlation lengths of the random fields cannot be determined by pure analytic treatments, additionally numerical simulations are performed for different values of the correlation length. The numerical simulations endorse the influence of the correlation length on the overall behavior. For a correlation length of 5mm a good conformity with the experimental results is obtained. Therefore, it is concluded, that the presented modeling approach is suitable to predict the elastic and plastic deformation of a set of tensile test specimens made of short fiber-reinforced composite sufficiently.

Relativistic atomic and molecular structure calculations start from very different premises to quantum electrodynamics (QED). QED predicts radiative corrections that can nly be approximated crudely in today's atomic structure programs. This paper describes rigorous simple algorithms for including accurate vacuum polarization and self-energy corrections in the GRASP programs.

Our aim is to implement and test some less known numerical methods that might be used to solve an advection equation. The description is restricted to information necessary to implement and test the presented methods and it is given only for one-dimensional case.

This paper examines the influence of physical parameters on the collapse dynamics of a spherical bubble filled with diatomic gas ($\kappa=7/5$). The present numerical investigation shows that each physical parameter affects the bubble collapse dynamics differently. After comparing the contribution of each physical parameter, it appears that, of all the parameters, the surrounding liquid environment affects the bubble collapse dynamics the most. Meanwhile, surface tension has the weakest influence and can be ignored in the bubble collapse dynamics. However, surface tension must be retained in the initial analysis since this, as well as the pressure difference jointly control initial bubble formation. As an essential part of this study, a general Maple code is provided.

This paper revisits a classic problem in physics - Hertz elastic dynamics of two colliding elastic spheres. This study obtains impact period in terms of hypergeometric function and successfully combines Deresiewicz's three segmental solutions into one single solution. Our numerical investigation confirms that Deresiewicz's inversion is a good approximation. As an essential part of this study, a general Maple code is provided.

The effects of the surface waves generated by the wind have a significant effect on the currents. A wave current coupled model plays an important role in the design of offshore structures. The interaction between fluids such as incompressible ocean waves and current and offshore structures is significant with many real-time applications in offshore engineering. These coupled models can be applied to Offshore Floating Production Operating and offloading (FPSO), Wind or current turbines and offshore pipelines. The complex issues related to the design are analyzed by using Computational Fluid Dynamics, which requires an investigation of the multiphase flow between wave and current and the structure which is considered restrictive due to the computational cost. If viscous effects are neglected then the single-phase flow models have been recommended, where wave-current interaction have been modelled successfully. Models have been developed where velocities and pressure are computed and the results can be verified with the experimental results available in the literature. In this study the existing numerical methods, mesh types are discussed along with their coupling methods. Here single-phase and multiphase models with small and medium movement are reviewed and their applications are highlighted.

Computational Fluid Dynamics (CFD) simulations, based on Reynolds Averaged Navier Stokes (RANS) models, are a useful tool for a wide range of coastal and offshore applications, providing a high fidelity representation of the underlying hydrodynamic processes. Generating input waves in the CFD simulation is performed by a numerical wavemaker (NWM), with a variety of different NWM methods existing for this task. While NWMs, based on impulse source methods, have been widely applied for wave generation in depth averaged, shallow water models, they have not seen the same level of adoption in the more general RANS based CFD simulations, due to difficulties in relating the required impulse source function to the resulting free surface elevation for non-shallow water cases. This paper presents an implementation of an impulse source wavemaker, which is able to self-calibrate the impulse source function to produce a desired wave series in deep or shallow water at a specific point in time and space. Example applications are presented, for a numerical wave tank (NWT), based on the opensource CFD software OpenFOAM, for wave packets in deep and shallow water, highlighting the correct calibration of phase and amplitude. Also, the suitability for cases requiring very low reflection from NWT boundaries is demonstrated. Possible issues in the use of the method are discussed and guidance for good application is given.

The application of instrumented indentation to assess material properties like Young’s modulus and micro-hardness has become a standard method. In recent developments, indentation experiments and simulations have been combined to inverse methods, from which further material parameters as yield strength, work hardening rate, and tensile strength can be determined. In this work, an inverse method is introduced by which material parameters for cyclic plasticity, i.e. kinematic hardening parameters, can be determined. To accomplish this, cyclic Vickers indentation experiments are combined with finite element simulations of the indentation with unknown material properties, which are then determined by inverse analysis. To validate the proposed method, these parameters are subsequently applied to predict the uniaxial stress-strain response of a material with success. The method has been validated successfully for a quenched and tempered martensitic steel and for technically pure copper, where an excellent agreement between measured and predicted cyclic stress-strain-curves has been achieved. Hence, the proposed inverse method based on cyclic nanoindentation, as a quasi-non-destructive method, could complement or even substitute the resource-intensive conventional fatigue testing in the future for some applications.

Cloud-based CNC is an emerging paradigm of Industry 4.0 where computer numerical control (CNC) functionalities are moved to the cloud and provided to manufacturing machines as a service. Among many benefits, C-CNC allows manufacturing machines to leverage advanced control algorithms running on cloud computers to boost their performance at low cost, without need for major hardware upgrades. However, a fundamental challenge of C-CNC is how to guarantee safety and reliability of machine control given variable Internet quality of service, especially on public Internet networks. We propose a three-tier redundant architecture to address this challenge. We then prototype tier one of the architecture on a 3D printer successfully controlled via C-CNC over public Internet connections, and discuss follow-on research opportunities.

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.

Different inlet structures have a significant impact on the internal flow characteristics of a solid-magnesium powder water ramjet engine. Based on the magnesium-water reaction model, a computational fluid dynamics (CFD) method is applied to establish a numerical simulation method for the internal flow field of the engine, and the internal flow characteristics of the engine under different inlet structure conditions are studied. The simulation results show that high-temperature gas can effectively promote the ignition of magnesium powder at the top of the combustion chamber, while accelerating the evaporation of the first inlet water and increasing the combustion rate of magnesium powder. The secondary inlet has the most significant effect on the temperature inside the combustion chamber. When the secondary inlet flow rate increases towards the top of the combustion chamber, it increases the amount of heat absorbed by the evaporating water at the top of the chamber, thereby reducing the temperature at the top of the combustion chamber. However, when the flow rate is low, it results in insufficient oxidizer at the top of the combustion chamber, which is unfavorable for the combustion of magnesium powder.

A time integration scheme based on semi-Lagrangian advection and Laplace transform adjustment has been implemented in a baroclinic primitive equation model. The semi-Lagrangian scheme makes it possible to use large time steps. However, errors arising from the semi-implicit scheme increase with the time step size. In contrast, the errors using the Laplace transform adjustment remain relatively small for typical time steps used with semi-Lagrangian advection. Numerical experiments confirm the superior performance of the Laplace transform scheme relative to the semi-implicit reference model. The algorithmic complexity of the scheme is comparable to the reference model, making it computationally competitive, and indicating its potential for integrating weather and climate prediction models.

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.

Transient flows result in unbalanced forces and high pressure in pipelines. Under these conditions, the combined effects of flow-induced forces along with sudden pipe displacements can create cracks in the pipeline, especially at the junctions. This situation consequently results in water leakage and reduced operational efficiency of the pipeline. In this study, displacements and stresses in a buried pressurized water transmission pipe installed by pipe jacking method are investigated using numerical modeling and considering interactions between fluid, pipe, and soil. The analyses were performed consecutively under no-flow, steady flow, and transient flow conditions, in order to investigate the effects of flow conditions on displacements and stresses in the system. Analyses of the results show that displacements and stresses in the jointed concrete pipes are significant under transient flow conditions. Moreover, because of pressure transient effects, maximum tensile stresses exceed the tensile strength of concrete at the junctions, leading to cracks and consequent water leakage.

This work exclusively focuses on the assessment of the thermal performance of a multifunctional facade panel incorporating PCM in foam layers, recurring to a hot box heat flux meter method to determine the thermal transmittance (U-value). The experimental setup is based on the steady-state approach using climatic chambers, assuring a stable thermal environment. Even small fractions of PCM achieved a small reduction in thermal amplitude. Numerical simulations using Ansys Fluent were developed to evaluate the performance of PCM use over a wide range of temperature boundary conditions and operating modes. These numerical models were calibrated and validated using the results of experimental tests, achieving a correlation factor of 0.9674, thus accurately representing a real-world scenario. The decrement factor (f) was used to analyse the data. It was identified that the efficiency of the panel and size of the optimum region increased with the PCM fraction growth. The simulated behaviour was optimum when the input mean temperature is of 20 °C for a room temperature of between 18-20 °C. The results show the significant potential of the multi-layered panel and the thermal regulator effect of the PCM incorporated on indoor space temperature to reach good thermal comfort levels.

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.

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

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.

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.

The flavor, aroma and color of coffee may change due to mechanical damage, reducing its qual-ity. To measure the mechanical response of materials, compression tests can be performed, de-termining the effect of force application on the mechanical behavior of the fruit at different stages of ripeness. In this context, this study aimed to analyze the deformation, strain energy and equivalent von Mises stress of coffee fruits at mature, semi-mature and immature maturation stages and compressed by collapse forces. Compression in three directions (x, y and z) was sim-ulated on coffee fruit models using the finite element method, with a compression support in the direction opposite to the force application axis. The numerical simulation of the compression process allowed to verify that the more mature the fruit, the greater the associated mean defor-mation (2.20 mm mm-1, 0.78 mm mm-1 and 0.88 mm mm-1), the lower the mean strain energy (0.07 mJ, 0.21 mJ and 0. 34 mJ) and the lower the mean equivalent von Mises stress (0.25 MPa, 1.03 MPa and 1.25 MPa), corresponding to ripe, semi-ripe and immature fruits; and these proper-ties may be used to analyze how processes in the production chain affect the amount of mechan-ical damage.

The article presents the results of testing the conditions of closing foundry voids during the hot forging operation of an ingot made of zirconium with 1% Nb alloy and use of the physical and numerical modelling, which is a continuation of research presented in a previous thematically related article published in the journal Materials. The study of the impact of forging operation parameters on the rheology zirconium with 1% Nb alloy was carried out on the Gleeble 3800 device. Using the commercial FORGE®NxT 2.1 program, a numerical analysis of the influence of thermo-mechanical parameters of the hot elongation operation in trapezoidal flat and rhombic trapezoidal anvils on the closure offoundry voids was performed. The analysis of the obtained test results was used to formulate recommendations on the technology of hot forging and the geometry of anvils, which ensure the closure foundry voids during the operation of hot forging

This paper presents some aspects of command and control of the new generation of machines-tools with numerical control. A problem faced by numerical control programmers is related to the length (relatively long) of the programs, with implications in the costs and duration of the manufacturing preparation. A method of reducing the main program and other redundant expenses is the introduction and use of processing cycles and subroutines for automatic calculation of the cutting tool trajectory, a fact that also allows the optimization and simplification of numerical control programming on flexible order manufacturing systems numeric. The method involves the use of predefined cycles and appropriate subroutines in order to improve the efficiency of the programming activity in a balanced way. It is a new and rapid method for optimize the structure of program using the numerical subroutine for command the movement of tool, in our case, with application for turning and milling.

In the present study, a simple approximation expression is given for the relationship between the period and amplitude of a simple pendulum under magnetic action. The analytical solution presented for the given problem. Two numerical quadrature methods Simpson's and Boole's method were utilized to demonstrate a new approximation of the problem. The results of the numerical quadrature have been compared to the exact solution. Absolute and relative mistakes of the problem have been presented. The Matlab program 2013R has created a numerical method that is used to analyze the outcome, It has been determined that the comparison's outcomes attest to the method's suitability and correctness. Moreover, the results show that numerical solution is suitable for the problem.

The twisted Wang transform distribution family, defined as the composition of parameter shifted inverse CDF function with an original CDF function, is found to be most suitable for matching low shape factor distributions, characterizing hard to fit or to simulate reinsurance portfolio losses for some perils from our previous study. Among them, the best form for matching a hard-to-fit empirical loss distribution for a specific peril, is the Exponential Fractional Extra Power 0 Distribution in (0,1) with CDF:.The simplest yet still a good form of this family is the Transformed Hyperbolic Tangent Distribution with CDF:,which has analytical formulas for the moments. The twisted Wang transform distribution family is compared and confirmed to be superior to all other well-known distribution families through extensive numerical optimization practice, distribution forms guesses, and computer-aided exploration experiments.

For embankment slope stability, soil cohesivity is one of the most important shear strength parameters. The effect of cohesiveness on the slope stability of a homogeneous embankment dam under rapid drawdown loading conditions was examined in this study. With the use of numerical modeling in GeoStudio, different situations were explored based on cohesiveness (0 kN/m², 5 kN/m², 10 kN/m^2, and 15 kN/m²) under a 1 m per day drawdown rate. The factor of safety value obtained from the long-term steady-state condition under 15 kN/m² cohesion was equivalent to a 116.8% increase from the one obtained under 0 kN/m² cohesion. The factor of safety values obtained after subjecting the embankment to different soil cohesion levels yielded a p-value of 1.91 x 10^-41, according to the Analysis of Variance. The calculated p-value (alpha value) is less than 0.05, suggesting that the differences between the examined cohesiveness values based on the list of the factor of safety values are statistically significant. The findings derived in this study show that it is significant to capture the effect of material characteristics during the design phase of an embankment dam.

Wave is a common environmental load that often causes serious damages to offshore structures. In addition, the stability for the submarine artificial slope is also affected by the wave loading. Although the landslide of submarine slopes induced by the waves has received wide attention, the research on the influence of solitary-wave is rare. In this study, a 2-D integrated numerical model is developed to investigate the stability of the foundation trench under the solitary wave loading. The Reynolds-Averaged-Stokes (RANS) equations are used to simulate the propagation of a solitary wave, while the current is realized by setting boundary inlet/outlet velocity. The pore pressure induced by the solitary wave is calculated by Darcy’s law and the seabed is characterized by Mohr-Coulomb constitutive model. Firstly, the wave model is validated through the comparison between analytical solution and experimental data. The initial consolidation state of slope under hydrostatic pressure is achieved as the initial state. Then, the factor of stability (FOS) for the slope corresponding to different distance between wave crest and slope top is calculated with the strength reduction method. The minimum of FOS is defined as the stability index for the slope with specific slope ratio during the process of dynamic wave loading. The parametric study is conducted to examine the effects of soil strength parameters, slope ratio and current direction. At last, the influence of upper slope ratio in a two-stage slope is also discussed.

Automotive aerodynamics comprises of the study of aerodynamics of road vehicles. Its main goals are reducing drag, minimizing noise emission, improving fuel economy, preventing undesired lift forces and minimizing other causes of aerodynamic instability at high speeds. The Ahmed body has the form of a highly simplified car, consisting of a blunt nose with rounded edges fixed onto a box-like middle section and a rear end that has an upper slanted surface, the angle of which can be varied. It retains vital features of real vehicles in order to study the flow fields around it and the related turbulence models which characterizes the actual flow at elevated Reynolds number. In the present study, the aerodynamic behavior of this body is investigated numerically by the aid of commercial CFD tool: Ansys Fluent. The results of the simulation are validated with available experimental data and results of the simulations from other literatures. The numerical data were obtained for a fixed free stream velocity of 25 m/s at the inlet. The simulations were performed at a fixed slant angle of 25 degree and zero yaw angle. The present study focuses on how local refinement of mesh inside the concerned body and the outside, helps affect the results and for which grid dependency test is the primary objective of this paper. The present study also helps demonstrate how the drag of the body behaves, which is mainly the effect of pressure drag force generated at the rear portion of the body. The study also focuses on important properties like the velocity magnitude at different locations for different meshing cases, and to capture the flow pattern in the front or near the wake region. The study can be further helpful to future researchers in determining resistance, fuel efficiency etc. helping designers to optimize in specialized areas for better efficiency.

A conventional 2D numerical model is improved by incorporating three submodels to consider different effects of secondary flow and a module for cohesive sediment transport. The model is applied to a meandering reach of Yangtze River to investigate secondary flow effects on cohesive sediment deposition, and a preferable submodel is selected based on the flow simulation results. Sediment simulation results indicate that the improved model predictions are in better agreement with the measurements in planar distribution of deposition as the increased sediment deposits caused by secondary current on the convex bank have been well predicted. Secondary flow effects on predicted amount of deposition become more obvious during the period when the sediment load is low and velocity redistribution induced by the bed topography is evident. Such effects vary with the settling velocity and critical shear stress for deposition of cohesive sediment. The bed topography effects can be reflected by the secondary flow submodels.

Numerical modelling is becoming a major tool for supporting environmental studies at different scales, thanks to the capability of up-to-date codes in reproducing the natural behaviour in a quite reliable manner. In evaluating the habitat diversity of anthropized rivers, however, many issues are rising because of the intrinsic complexity of the processes involved. Using a reach of the Po River in Italy as a case study, the present works aims to provide an estimate of the changes of the Eco-Environmental Diversity as a response to different constant flow discharges. The goals are achieved by means of two solvers of the iRIC suite, applied in sequence to firstly simulate the fluvial hydrodynamics and subsequently provide an estimate of the habitat conditions. Despite the simplifications intrinsically present in the models and the ones introduced for practical purposes, the results pointed out that the reduction of the flow discharge recently observed can threat the overall biological status of the river. Because of the modelling uncertainties, on the other side, these preliminary outcomes show the need for more research, both in terms of data acquisition and numerical schematization, for adequately evaluate the effects of transient hydrology on the river ecosystems. Moreover, additional field surveys are necessary to calibrate and validate the used model for having sufficiently reliable estimates.

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

Structural vibration induced by dynamic load or natural vibration is a nonnegligible factor in failure analysis. Based on vibrating boundary condition, impact resistance of shape memory alloy reinforced composites is investigated. In this investigation, modified Hashin’s failure criterion, Brinson’s model and visco-hyperelastic model are implemented into the numerical model to charactering the mechanical behavior of glass fiber/epoxy resin laminates, SMAs and interphase, respectively. First, fixed boundary condition is maintained in simulation to verify the accuracy of material parameters and procedures by comparing with experimental data. Then, a series of vibrating boundaries with different frequencies and amplitudes are applied during the simulation process to reveals the effect on impact resistances. The statistics of absorbed energy and contact force indicate that impact resistance of the composite under high frequency and large amplitude is lower than that under low frequency and small amplitude, and summarized by a mathematical expression.

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

Drafting formations have been long recognized as highly effective for reducing drag and enhancing athletic performance, particularly in race walking events. The precise spacing and positioning of the race walkers are critical to optimizing the effectiveness of drafting. In this study, the drag reduction of 15 drafting formations is investigated using wind-tunnel experiments and CFD numerical simulations. The results show excellent consistency in drag reduction rate between the two methods, with differences being within 10%. This can be attributed to spacing replacing body shape differences as the primary factor influencing drag reduction. Optimal double, triple, and quadruple drafting formations produce the same results in both wind-tunnel experiments and CFD simulation, resulting in drag reductions of 67%, 66%, and 81% (wind-tunnel) and 65%, 72%, and 85% (CFD). The sources of drag differences in the two methods are discussed from various aspects. The flow field obtained through CFD analysis is used to examine the mechanism of drag reduction, revealing that drafting formations have a significant shielding effect on incoming air, which reduces the number and speed of airflow impacting the core race walker. This shielding effect is identified as the primary cause of drag reduction. Using an empirical model for mechanical power output, optimal double, triple, and quadruple drafting formations enhance sport economy (4.4-5.7%), speed (3.61-4.67%), and performance (173.8-223.3s) compared to race walking alone. The findings can serve as a reference for race walkers' positioning strategy and provide insights for considering drafting formations in various running events.

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.

Since the first case of 2019 novel coronavirus disease (COVID-19) detected on Jan 30, 2020, in India, the number of cases rapidly increased to 3819 cases including 106 deaths as of 5 April 2020. Taking this into account, in the present work, we are studying a Bats-Hosts-Reservoir-People transmission fractional-order COVID-19 model for simulating the potential transmission with the thought of individual social response and control measures by the government. The real data available about infectious cases from $14^{th}$ March to $26^{th}$ March 2020 is analysed and accordingly various parameters of the model are estimated or fitted. The Picard successive approximation technique and Banach's fixed point theory have been used for verification of the existence and stability criteria of the model. Numerical computations are done utilizing the iterative Laplace transform method. In the end, we illustrate the obtained results graphically. The purpose of this study is to estimate the effectiveness of preventive measures, predicting future outbreaks and potential control strategies using the mathematical model.

Façade elements are a building component that satisfies multiple performance parameters. Among other things, “advanced façades” take advantage of hybrid solutions, such as assembling laminated materials. In addition to the enhanced mechanical properties that are typical of optimally composed hybrid structural components, these systems are energy-efficient, durable, and offer lighting comfort and optimal thermal performance, an example of which is the structural solution developed in collaboration with the University of Zagreb and the University of Ljubljana within the Croatian Science Foundation VETROLIGNUM project. The design concept involves the mechanical interaction of timber and glass load-bearing members without sealing or bonding the glass-to-timber surfaces. Following earlier research efforts devoted to the structural analysis and optimization of thus-assembled hybrid Cross-Laminated Timber (CLT)-glass façade elements, in this paper, special focus is given to a thermal and energy performance investigation under ordinary operational conditions. A simplified numerical model representative of a full-size building is first presented by taking advantage of continuous ambient records from a Live-Lab mock-up facility in Zagreb. Afterwards, a more detailed Finite Element (FE) numerical analysis is carried out at the component level to further explore the potential of CLT–glass façade elements. The collected numerical results show that CLT–glass composite panels can offer stable and promising thermal performance for façades similar to national and European standard requirements.

Corrugated cardboard is an ecological material, mainly because, in addition to virgin cellulose fibers also the fibers recovered during recycling process are used in its production. However, the use of recycled fibers causes slight deterioration of the mechanical properties of the corrugated board. In addition, converting processes such as printing, die-cutting, lamination, etc. cause micro-damage in the corrugated cardboard layers. In this work, the focus is precisely on the crushing of corrugated cardboard. A series of laboratory experiments were conducted, in which the different types of single-walled corrugated cardboards were pressed in a fully controlled manner to check the impact of the crush on the basic material parameters. The amount of crushing (with a precision of 10 micrometers) was controlled by a precise FEMat device, for crushing the corrugated board in the range from 10 to 70 % of its original thickness. In this study, the influence of crushing on bending, twisting and shear stiffness as well as a residual thickness and edge crush resistance of corrugated board was investigated. Then, a procedure based on a numerical homogenization, taking into account a partial delamination in the corrugated layers to determine the degraded material stiffness was proposed. Finally, using the empirical-numerical method, a simplified calculation model of corrugated cardboard was derived, which satisfactorily reflects the experimental results.

Facade elements are known to represent a building component with multiple performance parameters to satisfy. Among others, “advanced facades” take advantage of hybrid solutions, like the assemblage of laminated materials. In addition to enhanced mechanical properties that are typical of optimally composed hybrid structural components, these systems are energy-efficient, durable, and offer lightening comfort and optimal thermal performance. This is the case of the structural solution developed in joint research efforts of the University of Zagreb and the University of Ljubljana, within the Croatian Science Foundation VETROLIGNUM project. The design concept involves the mechanical interaction of timber and glass load-bearing members, without sealing or bonded glass-to-timber surfaces. Laminated glass infilled timber frames are recognized as a new generation of structural members with relevant load-carrying capacity (and especially the enhancement of earthquake resistance of framed systems), but also energy-efficient and cost-effective solutions.

Tornadoes remain an active subject of observational and numerical research due to the damage and fatalities they cause worldwide as well as poor understanding of their behavior, such as what processes result in their genesis and what determines their longevity and morphology. A numerical model executed on a supercomputer run at high resolution can serve as a powerful tool for exploring the rapidly evolving tornado-scale features within a simulated storm, but saving large amounts data for meaningful analysis can result in unacceptably slow model performance, an unwieldy amount of saved data, and saved data spread across millions of files. In this paper, a system for efficiently saving and managing hundreds of terabytes of compressed model output is described in order to support a supercomputer simulation of a tornadic supercell thunderstorm. The challenges of managing a simulation spanning over a quarter trillion grid volumes across the Blue Waters supercomputer are also described. The simulated supercell produces a long-track EF5 tornado, and the near-tornado environment is described during tornadogenesis, where single upward-growing vortex undergoes several vortex mergers before transitioning into a multiple vortex tornado.

The goal of this lab is to look at how the characteristics of cementitious materials are affected by bending loads combined with dry and wet cycles of sulfate. In this paper, three sets of specimens with different water–cement ratios were designed. This study applied 20% and 40% of the ultimate bending load as a continuous bending load to the cementitious material specimens for a sulfate erosion test. It was discovered that the specimens' porosity, flexural strength, and compressive strength all tended to improve initially before declining as the cementitious material was degraded in the sulfate solution aged. The flexural strength of the cementitious materials declined as the bending load increased, and their internal sulfate ion concentration rose. A transport model of sulfate ions in cementitious materials under the coupling effect of bending load and sulfate erosion was created by combining the improved version of Fick’s second law and chemical reaction kinetics. The simulation results show that the ion transport properties of sulfate ions in cementitious materials are in accordance with the experimental laws. These research results formed the basis of the unified standards for the coupling of load dry and wet cycle systems, as well as evaluating the durability of cementitious materials under corrosive environments.

First of all the paper presents a solid solution of piezoelectric ceramic that was synthesized according to the general formula Pb(1-x)Srx(Ti0.48Zr0.52)(1-y)NbyO3 with x = 0.05 and y = 0.02, using wet ceramic processing technology and using an oxide mix as prime material. The effects of dopants (Sr2+ and Nb5+) on phase constitution, on microstructure and on the dielectric and piezoelectric properties were determinate. The Zr/Ti ratio was chosen near the morphotropic phase boundary of the PZT system in studied composition. The XRD data revealed that the PZT doped composition had tetragonal perovskite structure. Secondly the paper presents the design of the novel piezo-motor based on a surface wave which translates the linear extension of different piezoelectric segments of a piezoelectric cylinder into a rotational bending movement. This rotational bending of the piezoelectric cylinder is then transformed into a continuous rotation of the rotor through a calculated contact. The design of the motor takes advantage of the high piezoelectric constants of the developed material in an optimal way in order to increase the energetic efficiency. A brief mathematical model of electromechanical answer of piezoelectric materials is presented as it was used in the modeling of the material during the numerical simulations.

The formation of hydrogen sulfide (H₂S) in petroleum reservoirs by anaerobic microbial activity is called biogenic souring of reservoirs and is a risk in the oil and gas industry as the compound is extremely toxic, flammable and corrosive. In this paper, we present a workflow and the tools to assess biogenic souring from a pragmatic engineering perspective. The retention of H₂S in the reservoir due to the reactions with iron-bearing rock minerals (e.g. siderite) is shown in a theoretical approach here and supported with literature data. Cases are given for two fields under secondary (waterflooding) and tertiary flooding with microbial enhanced oil recovery (MEOR). The use of the Monte Carlo method as a numerical modelling tool to incorporate uncertainties in the measured physical/chemical/ biochemical data is demonstrated as well. A list of studies conducted with different chemicals alone or in combination with various biocides to mitigate biogenic souring provides an overview of potential inhibitors as well as possible applications. Furthermore, the results of static and dynamic inhibition tests using molybdate are presented in more detail due to its promising mitigation ability. Finally, a three-step workflow for the risk assessment of biogenic souring and its possible mitigation is presented and discussed.

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.

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.

Polydimethylsiloxane (PDMS) is one of the most popular elastomers and has been used in different fields, especially in biomechanics research. Among the many interesting features of this material, its hyperelastic behaviour stands out, which allows the use of PDMS in various applications, like the ones that mimic soft tissues. However, the hyperelastic behaviour is not linear and needs detailed analysis, especially the characterization of shear strain. In this work, two approaches, numerical and experimental, were proposed to characterize the effect of shear strain on PDMS. The experimental method was implemented as a simple shear testing associated with 3D digital image correlation and was made using two specimens with two thicknesses of PDMS (2 and 4 mm). A finite element software was used to implement the numerical simulations, in which four different simulations using the Mooney-Rivlin, Yeoh, Gent, and Polynomial hyperelastic constitutive models were performed. These approaches showed that the maximum value of shear strain occurred in the central region of the PDMS, and higher values emerged for the 2 mm PDMS thickness. Qualitatively, in the central area of the specimen, the numerical and experimental results have similar behaviours and the values of shear strain are close. For higher values of displacement and thicknesses, the numerical simulation results move further away from experimental values.

Fractional differential mathematical model unfolding the dynamics of the COVID-19 pandemic in India is presented and explored in this paper. The purpose of this study is to estimate the future outbreak of disease and potential control strategies using mathematical models in India as a whole country as well as in some of the states of the country. This model is calibrated based on reported cases of infections over the month of April 2020 in India. We have used iterative fractional complex transform method to find approximate solutions of the model having modified Riemann Liouville fractional differential operator. We have also carried out a comparative analysis between actual and estimated cumulative cases graphically, moreover, most sensitive parameters for basic reproduction number$(R_0)$ are computed and their effect on transmission dynamics of COVID-19 pandemic is investigated in detail.

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

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

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.

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.

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.

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.

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.

This paper presents a numerical method for solving a two-dimensional subdiffusion equation with a Caputo fractional derivative. The problem considered assumes symmetry in both the equation’s solution domain and the boundary conditions, allowing for a reduction of the two-dimensional equation to a one-dimensional one. The proposed method is an extension of the fractional Crank-Nicolson method, based on the discretization of the equivalent integral-differential equation. To validate the method, the obtained results were compared with a solution obtained through Laplace transform. The analytical solution in the image of the Laplace transform was inverted using the Gaver-Wynn-Rho algorithm implemented in the specialized mathematical computing environment, Wolfram Mathematica. The results clearly show the mutual convergence of the solutions obtained by the two methods.

In recent years discrete element method (DEM) has gradually been applied to the traditional fluidization simulation of fine particles in micro fluidized bed (MFB). The application of DEM in simulating fast fluidization of fine particles in MFB has not yet received attention. This article presents a drag model that relies on the surrounding environment of particles, namely particle circumstance-dependent drag modle or PCDD model. Fast fluidization in MFB of fine particles is simulated using DEM based on the PCDD model. It is found that the gas-solid two-phase flow in MFB meets the common law of fast fluidization in circulating fluidized bed (CFB) to a certain extent, and also presents a special law that is significantly different from that in CFB. Simulations indicate that the local structure in MFB exhibits particle aggregation which is a natural property of fast fluidized, forming a structure where continuous dilute phase and dispersed concentrated phase coexist. The formation and fragmentation of particle clusters in different local regions of the bed have time synchronization. There exists strong effect of solid back-mixing in MFB, leading to relatively low outlet solid flux. The gas back-mixing effect is, however, not so distinct. The axial porosity shows a monotonically increasing distribution with the bed height, but does not strictly follow the single exponential distribution. The solid volume fraction at the bottom of the bed is significantly lower than the correlated value in CFB. The axial heterogeneous distribution of the cross-sectional average porosity in the lower half of the bed is also weakened. The radial porosity shows a distribution pattern of higher in the central region and lower in the sidewall region. Compared with the correlation results in CFB, the porosity near the central region in MFB is relatively low, while the porosity near the sidewall region is relatively high. The geometric size of the container, especially the radial size, is much smaller than that of a CFB. The relative area of contact between the MFB wall and particles is larger, and the wall friction factor becomes more significant. This may be the main reason for the time synchronization of particle agglomeration, severe particle back-mixing, low outlet solid flux, monotonically increasing axial porosity and weak core-annular structure of radial porosity.

Nowadays, there are implemented devices whose purpose is to perform massive computations by saving resources at the time they reduce the latency of arithmetic operations. These devices are usually GPUs, FPGAs and other specialised devices such as "Coral". Neural networks, digital filters and numerical simulators take advantage of the massively parallel operations of such devices. One way to reduce the amount of resources used is to limit the size of the registers that store data. This has led to the proliferation of numeric formats with a length of less than 32 bits, known as short floating point or SFP. We have developed several SFP’s for use in our neural network accelerator design, allowing for different levels of accuracy. We use a 16-bit format for data transfer and different formats can be used simultaneously for internal operations. The internal operations can be performed in 16-bit, 20-bit and 24-bit. The use of registers larger than 16-bit allows the preservation of fractional information while increasing precision. By leveraging some of the FPGA’s arithmetic resources, our design outperforms designs implemented from scratch and is competitive with specialized arithmetic circuits already implemented in the FPGA.

The amplified Arctic warming is one of several factors influencing atmospheric dynamics. In this work, we consider a series of numerical experiments to identify the direct role of the Arctic sea ice reduction process in forming climatic trends in the northern hemisphere. Aimed at this, we used two more or less independent mechanisms of ice reduction. The first is traditionally associated with increasing the concentration of carbon dioxide in the atmosphere from the historic level of 360 ppm to 450 ppm and 600 ppm. This growth increases air temperature and decreases the ice volume. The second mechanism is associated with a reduction in the reflectivity of ice and snow. We assume that comparing the results of these two experiments allows us to judge the direct role of ice reduction. The most prominent consequences of ice reduction, as a result, were the weakening of temperature gradient at the tropopause level in mid-latitudes, the slower zonal wind at 50-60∘N, intensification of wave activity in Europe, Western America, and Chukotka, and its weakening in the south of Siberia and Kazakhstan. We also consider how climate change may alter regimes such as blocking and stationary Rossby waves. The study used the INM-CM48 climate system model .

The need for quick and easy deflection calculations of various prefabricated slabs causes that simplified procedures and numerical tools are used more and more often. Modelling of full 3D finite element (FE) geometry of such plates is not only uneconomical but often requires the use of complex software and advanced numerical knowledge. Therefore, numerical homogenization is an excellent tool, which can be easily employed to simplify a model, especially when accurate modelling is not necessary. Homogenization allows for simplifying a computational model and replacing a complicated composite structure with a homogeneous plate. Here, a numerical homogenization method based on strain energy equivalence is derived. Using the method proposed, the structure of the prefabricated concrete slabs reinforced with steel spatial trusses is homogenized to a single plate element with an effective stiffness. There is a complete equivalence between the full 3D FE model built with solid elements combined with truss structural elements and the simplified homogenized plate FE model. The method allows for the correct homogenization of any complex composite structures made of both solid and structural elements, without the need to perform advanced numerical analyses. The only requirement is a correctly formulated stiffness matrix of a representative volume element (RVE) and appropriate formulation of the transformation between kinematic constrains on RVE boundary and generalized strains.

The adoption of Pump as Turbines (PATs) both in small scale hydroelectric plants and water supply systems has brought different advantages, the most recognized being cost-effectiveness as compared to other hydroturbines. However, due to their lack of flow control ability, their intolerance to off-design operations constitutes one tough shortfall. Moreover, this technology’s newness leads to PAT flow dynamics still being ununderstood even to date. Therefore, this study intends to numerically investigate the mixed flow PAT’s pump mode flow dynamics for five operating conditions expanding from optimum (1QBEP) to deep part-load (0.41 QBEP) conditions. Moreover the effect of runner blade angle on the same has been investigated where three angles namely -2°, 0°, and 2° have been considered. PAT flow stability was found to deteriorate as the flow decreased, where associated pressure pulsation level worsened at different flow zones. In addition, the blade angle increase led correspondingly increasing flow unsteadiness and pressure pulsation levels, where the pulsation frequencies from rotor-stator interactions were dominant for most of flow zones. This study’s findings are of a crucial importance to both scientific and engineering communities as they contribute to thorough understanding of PAT flow dynamics.