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.

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

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.

ABSTRACT. In this paper, we suggest a biomathematical numerical method for analysing mRNA nucleotides sequences based on UA/CG Fibonacci numbers proportions. This method is used to evaluate then compare the spike genes related to the main SARS-CoV2 VARIANTS currently circulating within the world population. The 10 main results proposed to be reproduced by peers are: 1/ SARS-CoV2 genome and spike evolution in one year 2020-2021. 2/ SARS-CoV2 Origins. 3/ Comparing 11 reference variants spikes. 4/ analysing 32 CAL.20C California variant patients spikes. 5/ Toward a meta mRNA Fibonacci gene end message code. 6/ Analysing S501 UK, S484 South Africa and « 2 mutations » INDIA variants. 7/ Suggesting a possible variants spike mRNA palindrome symmetry metastructure improving mRNA stability then infectiousness. 8/ Analysing Fibonacci Metastructures in the mRNA coding for the vaccines PFIZER and MODERNA. 9/ Does the CG-rich modification of the synonymous codons of the spikes of the 2 mRNA vaccines affect the expression and quantity of SARS-CoV2 antibodies? 10/ The exceptional case of the Brazilian variant P.1. Particularly, we suggest the following conjecture at mRNA folding level : CONJECTURE of SARS-CoV2 VARIANTS: The growth of long Fibonacci structures in the shape of "podiums" for almost all of the variants studied (UK, California, South Africa, India, etc.) suggests the probable folding of the Spike mRNA in the form of a "hairpin", which can strengthen the cohesion and the lifespan of this mRNA. Finally, we show that these kinds of Fibonacci matastructures disapear TOTALLY by analysing the published mRNA sequences of PFIZER and MODERNA vaccines. One fact is certain, the two mRNAs of the Moderna and Pfizer vaccines will result in a low functionality of the spike vaccine. This is because their designers by seeking greater stability, have doped to build CG rich sequences which, as soon as they are inserted into the human host, will, paradoxically, seek to mutate, like SARS-CoV2 variants, towards CG ==> UA forms in order to improve their STABILITY and LIFETIME. We conclude using new biomathematics theoretical methods (Master code and numerical standing waves), and comparing the Spikes of the two vaccines Moderna and Pfizer, that there will be very probable differences in stability and shelf life of the two respective mRNAs vaccines. However, “State of the Art” analyzes will disclose that their two protein sequences are strictly identical. By modified their synonymous codons using different strategies, no one can guarantee that the quantity of antibodies generated will be identical in the two cases. We wish to draw attention to the great ADAPTATION power - at the global scale of their genomes - of the most infectious VARIANTS, such as the BRAZIL 20J / 501Y.V3 variant (P.1). This is very worrying for the VACCINES <==> VARIANTS run: We demonstrate how the Brazilian variant P.1 which becomes uncontrollable in Brazil in April 2021 has a level of organization of long metastructures of 17,711 bases covering the genome which is 3.6 more important than that of the 2 reference genomes SARS-CoV2 and worldwide D614G. We suggest that this high level of overall structure of this variant contributes to the stability of this genome and, might explain its greater contagiousness.

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.

Air filtration mechanisms in the composite filter media used in practical applications are important and challenging to understand because the component fibers could have various size scales and morphologies. In this work, a three-dimensional digital model of nanofiber-based filter media was reconstructed for the first time based on the X-ray tomography data for the cellulose substrate and the Focused Ion Beam-Scanning Electron Microscope (FIB-SEM) image analysis for the several microns thick (3.82-7.90 μm) electrospun polyvinylidene fluoride (PVDF) nanofiber membrane. Besides the high-resolution model where the details of the fibrous structures were fully resolved, another low-resolution model with approximated unresolved structures was also established. Filtration simulations utilizing these models were conducted considering the drag force, Brownian diffusion and aerodynamic slip. The simulated filtration efficiencies agreed well with the experiments for particles of 70-400 nm, including the most penetrating particle size (MPPS, 100-200 nm). Moreover, the structure-resolved models had higher accuracy but higher computational costs, while the unresolved simulations saved much running time but over-predicted the filtration efficiency, especially for smaller particles (<100 nm). Our study presents a comprehensive strategy for investigating the composite filter media with multiscale complex structures using a combination of advanced characterization technologies and modular simulation models.

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.

In this paper an integrated methodology for the coupling between 1D- and 3D-CFD simulation codes is presented, which has been developed to support the design and calibration of new diesel engines. The aim of the proposed methodology is to couple 1D engine models, which may be available in the early-stage engine development phases, with 3D predictive combustion simulations, in order to obtain reliable estimates of engine performance and emissions for newly designed automotive diesel engines. The coupling procedure features simulations performed in 1D-CFD by means of GT-SUITE and in 3D-CFD by means of Converge, executed within a specifically designed calculation methodology. An assessment of the coupling procedure has been performed by comparing its results with experimental data acquired on an automotive Diesel engine, considering different working points including both part load and full load conditions. Different multiple injection schedules have been evaluated for part-load operation, including pre and post injections. The proposed methodology, featuring detailed 3D chemistry modeling, was proven to be capable to properly assess pollutant formation, specifically to estimate NOx concentrations. Soot formation trend was also well-matched for most of the explored working points. The proposed procedure can therefore be considered as a suitable methodology to support the design and calibration of new Diesel engines, thanks to its ability to provide reliable engine performance and emissions estimations from the early-stage of a new engine development.

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.

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.

The steady streaming (SS) phenomenon is gaining increased attention in the microfluidics community, because it can generate net mass flow from the zero-mean vibration. We developed numerical simulation and experimental measurement tools to analyze this vibration induced flow, which has been challenging due to its unsteady nature. Validity of these analysis methods is confirmed by comparing the three-dimensional (3D) flow field induced around a cylindrical micropillar under circular vibration. In the numerical modeling, we directly solved the flow in the Lagrangian frame so that the substrate with a micropillar becomes stationary, and the result was converted to the Eulerian frame to compare them with the experimental results. The present approach enables to avoid the introduction of moving boundary or small perturbation approximation. The flow field obtained by the micro particle image velocimetry (PIV) measurement supported the three-dimensionality observed in the numerical results, which could be important for controlling the mass transport and manipulating particulate objects in the microfluidic systems.

Energy dissipation is one of the most important factors in choosing stepped spillways. However, few studies are about energy dissipation with different horizontal face angles. In this paper, the realizable k-ε turbulence model was used to study the flow field and energy dissipation rates of different stepped spillways with five horizontal face angles in the skimming flow regions. It shows that the free water surfaces and energy dissipation rates change with horizontal face angles and the energy dissipation rates initially decrease and then increase as the horizontal face angle increases. And the horizontal face angles cause the great changes of flow field.

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.

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

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

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

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

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.

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

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

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

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.

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

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

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

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

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.

Recently, NOAA has included GNSS (Global Navigation Satellite System) Radio Occultation (RO) data as one of the crucial long-term observables for weather and climate applications. To include more GNSS RO data in the numerical weather prediction system, the NOAA Commercial Weather Data Pilot program (CWDP) started to explore the commercial RO data available on the market. After two rounds of pilot studies, the CDWP decided to award the first Indefinite Delivery Indefinite Quantity (IDIQ) contract to GeoOptics and Spire Incs. in 2020. This study examines the quality of Spire data products for weather and climate applications. Spire RO data are collected from commercial CubeSats through careful comparison with the data from Formosa Satellite Mission 7–Constellation Observing System for Meteorology, Ionosphere, and Climate-2 (COSMIC-2), ERA-5, and high-quality radiosonde data. The results demonstrated that although with lower Signal-Noise-Ratio (SNR) in general, the pattern of the lowest penetration height for Spire is similar to those for COSMIC-2. The Spire and COSMIC-2 penetrate heights are between 0.6 and 0.8 km altitude at the tropical oceans. Although using different GNSS RO receivers, the precision of Spire STRASP receivers is of the same quality as those of COSMIC-2 Global Positioning System - GPS, GALILEO, and GLObal NAvigation Satellite System – GLONASS (TGRS) receivers. The retrieval accuracy from Spire is very compatible with those from COSMIC-2. We validated Spire temperature and water vapor profiles by comparing them with collocated radiosonde data. Generally, over the height region between 8 km and 16.5 km, the Spire temperature profiles match those from RS41 RAOB very well with temperature biases &lt; 0.02 K. Over the height range from 17.8 to 26.4 km, the temperature differences are ~-0.034 K with RS41 RAOB being warmer. We also estimated the error covariance matrix for Spire, COSMIC-2, and KOMPSAT-5. Results showed that the COSMIC-2 estimated error covariance values are slightly more significant over the oceans at the mid-latitudes (45oN-30oN and 30oS-45oS), which may also be owing to COSMIC-2 SNR being lower at those latitudinal zones.

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.

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.

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.

Polymethyl methacrylate (PMMA) polymer is widely used in various fields today. In order to thoroughly reveal the structural impact performance of PMMA materials, this paper is based on the comprehensive Johnson-Cook constitutive model and damage failure model, and accurately confirm the J-C constitutive and damage failure model parameters of PMMA through material test data, The dynamic process of steel bullet impacting PMMA plate structure is analyzed by finite element software ABAQUS. The calculation results show that the numerical simulation results in this paper are in good agreement with the experimental test data. Therefore, the feasibility and accuracy of the impact analysis of PMMA structures based on J-C constitutive and damage failure models in this paper are studied. Finally, based on the J-C constitutive model and damage failure model, the variation of the residual velocity of the bullet with the PMMA plate thickness is analyzed in depth, that is, the results show that the residual velocity of the bullet has a certain linear relationship with the plate thickness.

Taking a certain type of test centrifuge as the research object, the finite element model of the test centrifuge is established, the vibration characteristics and aerodynamic performance of the test centrifuge are analyzed, and the structural optimization design of the test centrifuge is carried out. In this paper, the load is applied according to the actual working condition of a certain type of test centrifuge. The vibration of the mounting seat of the test centrifuge is analyzed and the structure of the mounting seat is optimized. After optimization, the vibration of the mounting seat is 77.38 % lower than that of the original mounting seat. Then, the aerodynamic analysis of the test centrifuge is carried out. The analysis results show that the test centrifuge moves more smoothly under the whole package shell and the steady flow shell, and the resistance decreases and the shaft load decreases. Finally, the steady flow shell of the test centrifuge is optimized. The analysis shows that the increase of the width of the steady flow shell can reduce the resistance coefficient, which is helpful to the stability of the test centrifuge during operation and reduces the unbalanced response of the system caused by air resistance.

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 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 energy transfer process of laser selective melting is very complex. To study the effect of la-ser selective melting on the microstructure and properties of 18Ni-300 martensitic steel, ABAQUS was used to simulate the temperature of laser cladding 18Ni-300 martensitic steel at different time points and different laser power. The results show that the cross-section shape of the molten pool changes from round to oval With the increase of laser power, the higher the peak value of temperature time curve, the greater the temperature gradient; and the laser clad-ding experiment of 18Ni-300 martensitic steel was carried out, and the microstructure and me-chanical properties of the samples under different laser power were analyzed. The results show that with the increase of laser power, the grain size of the cladding layer becomes smaller and the microstructure becomes more compact; the hardness of the side surface of the sample is higher than that of the upper surface, and the tensile strength and elongation show a trend of first increasing and then decreasing.

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.

Melt mixing is a convenient method to prepare polymer nanocomposites, but the extent of the dispersion of the solid filler reached is often limited, and may compromise the anticipated performance of these materials during service. Since the efficiency of extensional flows on dispersion is now well recognized, several mixers were designed with the aim of inducing both shear and extensional flow components. This work combines experimental and numerical data to better understand the kinetics of the dispersion of graphite nanoplates in a polypropylene melt, using a mixing device that consists of a series of stacked rings with equal outer diameter and alternating larger and smaller inner diameters, thus creating a series of converging/diverging flows. Numerical simulation of the flow using the opensource OpenFOAM software assuming both inelastic and viscoelastic responses predicted the velocity, streamlines, flow type and shear and normal stress fields for the mixer. Experimental and computed data were combined to determine the trade-off between the local degree of dispersion of the PP/GnP nanocomposite, measured as Area ratio, and the absolute average value of the hydrodynamic stresses multiplied by the local cumulative residence time. From considerations based on a theoretical approach to dispersion, the cohesive strength of the GnP agglomerates studied was estimated to be in the range 3 - 10kPa.

While the knowledge about the development of understanding positive integers is rapidly growing, the development of understanding zero is not well-known. Here we tested several components of preschoolers’ understanding zero: whether they can use empty sets in numerical tasks, whether they can use empty sets soon after they understand the cardinality principle, whether they know what the word “zero” refers to, and whether they categorize zero as a number. The results show that preschoolers can handle empty sets in numerical tasks as soon as they understand the cardinality principle or even earlier, and some of them know that these sets are labeled as “zero.” However, they are unsure whether zero is a number. These results identify three components of knowledge about zero: operational knowledge, linguistic knowledge, and meta knowledge. To account for these results we propose that preschoolers might understand numbers as the properties of items or objects in a set. In this view, zero cannot be a number, because an empty set does not include any items, and the missing items cannot have any property, excluding also the number property. This model might explain why zero is handled correctly in numerical tasks, while it is not regarded to be a number.

Due to the accuracy of numerical calculation of fluid flow inside a hydrocyclone can be obtained using Computational Fluid Dynamics (CFD), highly modified super computers are used to simulate the fluid flow and track particle motion inside a hydrocyclone. This paper deals with the numerical study using three multiphase models viz. Volume of fluid, Mixture and Eulerian model. The dimensions of the hydrocyclone taken into consideration for numerical analysis is same as considered by Rajamani. Validation of axial and tangential velocities at different strategically decided axial stations, RMS axial and tangential velocity profiles of the hydrocyclone is done using Reynolds Stress Model (RSM). The hydrocyclone model has been designed in Creo 3.0 using the same dimensions which later was imported to CFD for meshing. Fine hexagonal mesh numbering up to 5 lacs were constructed to obtain optimum results. Fluid flow was allowed to be developed in ANSYS FLUENT 16.2. Entire simulation took 96 hours to generate results and track particle movements inside the hydrocyclone. The particle tracking has been done using three multiphase model. The first being the volume of fluid was used for validation purposes and the comparison of the Mixture and Eulerian model are the basic focus of this research work. Conclusive results indicate that usage of different multiphase model does not result in variation in particle motion. The slight variation in grade efficiency values are hardly noticeable. The Mixture model and Eulerian model predict lower separation efficiency as compared with Volume of fluid multiphase model.

In this manuscript, an automated framework dedicated to design space exploration and design optimization studies is presented. The framework integrates a set of numerical simulation, computer-aided design, numerical optimization, and data analytics tools using scripting capabilities. The tools used are open-source and freeware, and can be deployed on any platform. The main feature of the proposed methodology is the use of a cloud-based parametrical computer-aided design application, which allows the user to change any parametric variable defined in the solid model. We demonstrate the capabilities and flexibility of the framework using computational fluid dynamics applications; however, the same workflow can be used with any numerical simulation tool (e.g., a structural solver or a spread-sheet) that is able to interact via a command line interface or using scripting languages. We conduct design space exploration and design optimization studies using quantitative and qualitative metrics, and to reduce the high computing times and computational resources intrinsic to these kinds of studies, concurrent simulations and surrogate-based optimization are used.

Studies on the detection of layers with elevated black carbon aerosol (BC) concentrations and the formation conditions of these layers help understand the vertical distribution of BC concentrations, which will provide a basis for the assessment of climate effects and early BC pollution warnings. By using the Weather Research and Forecasting with Chemistry (WRF-Chem) numerical model, we performed a numerical simulation analysis on the authenticity of strong elevated BC concentration layers that were detected by an aircraft in the mixing layer over Harbin, China, which is a high-emission area, on a clear sunny afternoon in the early heating period of 2016. We then discuss possible problems and solutions when non-vertical paths are used to detect the vertical distribution of BC concentrations. Finally, we discuss the favorable conditions for the formation of elevated BC concentration layers by weak vertical flow. The results show that the horizontal variability of BC concentration in the mixing layer in the observation area in Harbin was sufficiently large during the measurement. This produced a false elevated layer, as detected by the aircraft during one round of spiral flight in the mixing layer. The root mean square of the horizontal distribution of BC concentration did not change with height in the mixing layer during the daytime, but it decreased with the thickness of the mixing layer and was higher in the mixing layer than in the free atmosphere. Therefore, the thinner the mixing layer, in which the vertical distribution of the BC concentration is detected in an inclined path, the stronger interference of the horizontal variability on the detected results. When a spiral flight detection path is used, the aircraft should fly at least two rounds in the mixing layer. In the daytime, due to strong turbulence in the mixing layer, weak vertical uplift is not favorable for the occurrence of elevated BC concentration layers in the mixing layer. In the nighttime, if weak vertical uplift is well matched with the BC concentration or its vertical gradient, elevated BC concentration layers can be formed in the atmosphere. Compared with upper layers far from the ground, nighttime elevated layers are easier to form in lower layers near the ground because high BC concentrations or large vertical gradients are more likely to occur in the lower layers. Both cases facilitate the occurrence of large vertical upward transport rates of BC.

This article is devoted to Feller's diffusion equation which arises naturally in probabilities and physics (e.g. wave turbulence theory). If discretized naively, this equation may represent serious numerical difficulties since the diffusion coefficient is practically unbounded and most of its solutions are weakly divergent at the origin. In order to overcome these difficulties we reformulate this equation using some ideas from the Lagrangian fluid mechanics. This allows us to obtain a numerical scheme with a rather generous stability condition. Finally, the algorithm admits an elegant implementation and the corresponding Matlab code is provided with this article under an open source license.

The assimilation of Soil Moisture and Ocean Salinity (SMOS) data into the ECMWF (European Centre for Medium Range Weather Forecasts) H-TESSEL (Hydrology revised - Tiled ECMWF Scheme for Surface Exchanges over Land) model is presented. SMOS soil moisture (SM) estimates have been produced specifically by training a neural network with SMOS brightness temperatures as input and H-TESSEL model SM simulations as reference. This can help the assimilation of SMOS information in several ways: (1) the neural network soil moisture (NNSM) data have a similar climatology to the model, (2) no global bias is present with respect to the model even if regional differences can exist. Experiments performing joint data assimilation (DA) of NNSM, 2 metre air temperature and relative humidity or NNSM-only DA are discussed. The resulting SM was evaluated against a large number of in situ measurements of SM obtaining similar results to those of the model with no assimilation, even if significant differences were found from site to site. In addition, atmospheric forecasts initialized with H-TESSEL runs (without DA) or with the analysed SM were compared to measure of the impact of the satellite information. Although, NNSM DA has an overall neutral impact in the forecast in the Tropics, a significant positive impact was found in other areas and periods, especially in regions with limited in situ information. The joint NNSM, T2m and RH2m DA improves the forecast for all the seasons in the Southern Hemisphere. The impact is mostly due to T2m and RH2m, but SMOS NN DA alone also improves the forecast in July- September. In the Northern Hemisphere, the joint NNSM, T2m and RH2m DA improves the forecast in April-September, while NNSM alone has a significant positive effect in July-September. Furthermore, forecasting skill maps show that SMOS NNSM improves the forecast in North America and in Northern Asia for up to 72 hours lead time.

The Lower Rio Grande Valley, South Texas is considered one of the more vulnerable coastal areas to flooding related with abrupt climate changes. From 1980-2017, there were 7 flooding events, 57 severe storm events, and 8 tropical cyclone events with losses exceeding $1 billion in the State of Texas, according to NOAA NCEI. Coastal flooding is typically a result of storm surge and heavy rainfall produced by hurricanes and tropical storms. In this study, the two-dimensional hydrodynamic flow circulation model is developed to predict the Lower Rio Grande Valley coastal area inundation due to the hurricane storm surge, especially in the case of Hurricane Beulah, 1967. The tropical cyclone properties and tidal constituents were assigned to the updated watershed geographic information with the bottom bathymetric and roughness data. For model validation, the Hurricane Dolly 2008 storm surge due to Hurricane Beulah at the coast and the storm surge reaches up to approximately 40 kilometers west from the coast through a natural river channel. This model can be used for a reliable engineering tool for the coastal hazard emergency management and disaster mitigation.

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

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

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

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.

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

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

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

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.

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.

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

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.

Large aluminum alloy ring forgings are the core components of heavy-duty rocket fuel storage tanks, and the large residual stress inside the rings leads to poor shape accuracy of large thin-walled parts. The initial stress of the 2219 aluminum alloy ring blank was tested using the drilling method, and the creep constitutive coefficient of the 2219 aluminum alloy was determined through stress relaxation tests. The numerical simulation processes of thermal stress relief (TSR), vibration stress relief (VSR), and thermal-vibration stress relief (TVSR) were compared and established. Through the correlation analysis between the actual measurement results of residual stress and the simulation results, it can be seen that the strong correlation in three directions at each measurement point accounts for over 37.5%, and the moderate correlation accounts for over 62.5%. This indicates that the numerical simulation model of 2219 aluminum alloy ring containing initial residual stress can accurately reflect the size and distribution of residual stress inside the actual ring. The simulation results show that the derived constitutive model can describe the stress relaxation process of TVSR by combining a single thermal time effect stress relaxation constitutive theory with a VSR plastic deformation material model. The simulation models established above were used to calculate the residual stress homogenization ability of three types of aging. The results showed that VSR, TSR, and TVSR can homogenize and reduce the residual stress field inside the ring, improve the distribution of residual stress inside the ring, and have a better overall homogenization ability of TVSR. The VSR control has a certain effect on reducing and homogenizing residual stress, but compared with TSR and TVSR, the reduction and homogenization ability of residual stress control is limited. The homogenization control effect TVSR>TSR>VSR, and the maximum equivalent stress homogenization rates of VSR, TSR, and TVSR are 52.8%, 80.6%, and 82.2%, respectively. Then, numerical simulation technology is used to study how the initial residual stress in the blank will cause the deformation of the ring during the thin-walled machining process. The roundness error theory of the minimum containment area method is applied to evaluate the deformation degree during the thin-walled numerical machining process, and the TVSR method is used for stress regulation. The deformation law of the thin-walled machining of the ring under different aging parameters is studied.

The main objective of this study is to find out the influence of cooperation and intra-specific competition in escaping predation through refuge of the prey population and the effect of the two intra-specific interactions on the dynamics of prey-predator systems. For this purposes, two mathematical models with Holling type-II functional response function have been proposed and analyzed. The first model includes cooperation among prey populations whereas the second one incorporates intra-specific competition. Existence conditions and stability of different equilibrium points of both models have been carried out for the qualitative behaviour of the systems. Refuge through intra-specific competition has a stabilizing role, whereas cooperation has a destabilizing role on the system dynamics. Periodic oscillations are observed in both the systems through Hopf bifurcation. From the analytical and numerical findings, we conclude that intra-specific competition affects the prey population and continuously checks the refuge class under a critical value, and thus it never becomes too large to cause predator extinction due to food scarcity. Conversely, cooperation leads maximum individuals to escape predation through the refuge so that predators will suffer from low predation success.

On the basis of the force and safety factor analysis of traditional retaining wall, a new analytical method of force and displacement of slope is suggested, the numerical theoretical solution of the stress distribution of the sliding body can be obtained under the conditions, which the stresses distribution are satisfied with the differential equilibrium equations, the boundary conditions, the compatibility equation and the macroscopic equilibrium equations. The interface stresses between the sliding body and the retaining wall are continuous, the theoretical solution of retaining wall stress distributions can be gotten, while the stress differential equilibrium equation, the compatibility equation, macroscopic force and moment balance equations are satisfied. The strain and displacement solutions can be obtained by using Duncan Chang and Hooke constitutive equations for the slope and the retaining wall respectively. The transfer station of landfill in the Guandukou Town of Badong County is taken as an example, the results of the sliding body and retaining wall analysis show: The stress and strain solutions of a slope and a retaining wall can be obtained by the proposed method, the forces between the retaining wall and the sliding body include shear stress and normal stress, and the stress distributions change with the deformation of the sliding body, the shape and material of retaining wall can be optimized according to the calculation results，it is feasible for the proposed analysis method of slope with retaining wall design to be run many years.

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

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.

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.

{\bf Abstract.} \\ {\bf Background:} Bangladesh hosts more than 800,000 Rohingya refugees from Myanmar. The low health immunity, lifestyle, access to good healthcare services, social-security cause this population to be at risk of far more direct effects of COVID-19 than the host population. Therefore, evidence-based forecasting of the COVID-19 burden is vital in this regard. In this study, we aimed to forecast the COVID-19 burden among the Rohingya refugees of Bangladesh to keep up with the disease outbreak's pace, health needs, and disaster preparedness. \noindent{\bf Methodology and Findings:} To estimate the possible consequences of COVID-19 in the Rohingya camps of Bangladesh, we used a modified Susceptible-Exposed-Infectious Recovered (SEIR) transmission model. All the values of different parameters used in this model were from the Bangladesh Government’s database and the relevant emerging literature. We addressed two different scenarios, i.e., the best-fitting model and good fitting model with unique consequences of COVID-19. Our best fitting model suggests that there will be good control over the transmission of the COVID-19 disease. At the end of December 2020, there will be only 169 confirmed COVID-19 cases in the Rohingya refugee camps. The average basic reproduction number ($ \mathcal{R}_{0} $) has been estimated to be 0.7563. \noindent{\bf Conclusion:} Our analysis suggests that, due to the extensive precautions from the Bangladesh government as well as other humanitarian organizations, the coronavirus disease will be under control if the maintenance continues like this. Although detailed and pragmatic preparedness should be adopted for the worst scenario.

A solution for the Prandtl-Blasius equation is essential to all kinds of boundary layer problems. This paper revisits this classic problem and presents a general Maple code as its numerical solution. The solutions were obtained from the Maple code, using the Runge-Kutta method. The study also considers convergence radius expanding and an approximate analytic solution is proposed by curve fitting. Similarly, the study resolves some boundary layer related problems and provide relevant Maple codes for these.

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.

Advancement in physics often results from analyzing numerical data and then creating a theoretical model that can explain and predict those data. In the field of lens design, the reverse is true: longstanding theoretical understanding is being overtaken by more powerful numerical methods.

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

A microstructure-based model is developed to study the effective anisotropic properties (magnetic, dielectric or thermal) of two phase particle-filled composites. The Green’s function technique and the effective field method are used to derive theoretically the homogenized (averaged) properties for a representative volume element containing isolated inclusion and an infinite, chain-structured particles. Those results are compared with the FE approximations conducted for the assumed representative volume element. In addition, as a special case, the Maxwell–Garnett model is retrieved when particle interactions are not taken into consideration. We shall also give some information on the optimal design of the effective anisotropic properties taking into account the shape of magnetic particles.

Marine cold air outbreaks (MCAOs) are large-scale phenomena in which cold air masses are advected over open ocean. It is well-known that these events are linked to the formation of polar lows and other mesoscale phenomena associated with high wind speeds, and that they therefore in some cases represent a hazard to maritime activities. However, it is still unknown whether MCAOs are generally conducive to higher wind speeds than normal. Here this is investigated by comparing the behaviour of ocean surface wind speeds during MCAOs in three atmospheric reanalysis products with different horizontal grid spacings, along with case studies using a convection-permitting numerical weather prediction model. The study regions are the Labrador Sea and the Greenland–Iceland–Norwegian (GIN) Seas, where MCAOs have been shown to be important for air–sea interaction and deep water formation. The main findings are: 1) Wind speeds during the most extreme MCAO events are stronger than normal and higher than wind speeds during less severe events; 2) The peak times of MCAO usually occur when baroclinic waves pass over the regions; and 3) Reanalyses with grid spacings of more than 50 km appear to underestimate winds driven by the large ocean–atmosphere energy fluxes during MCAOs. It is also shown that while the strong wind episodes during MCAOs generally last for just a few days, MCAOs can persist for up to 50 days. These findings demonstrate that it would be worthwhile to forecast MCAOs, and that it might be possible to do this beyond the standard weather forecasting range of up to 10 days.

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

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

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

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

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

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

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.

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

The objective of this study is to investigate the use of excavation behavior prediction and analysis software in the field of civil en-gineering. PLAXIS and SoilWorks, widely used in the engineering and academic communities, were employed as research tools for analyzing gravel layers and deriving appropriate ranges of soil elastic modulus for different software programs. By collecting excavation cases in the Xindian area of New Taipei City, Taiwan, and conducting deep excavation analysis using the displacement of retaining walls, this study provides recommended ranges for localized selection of soil elastic modulus. The suggested range for PLAXIS is 7,840 N/m² to 9,800 N/m², while Soil-Works suggests a range of 2,450 N/m² to 3,430 N/m². These ranges allow for reasonable estimation of maximum deformations during the final excavation stage. Based on the research findings, it is recommended that engineers refer to the provided ranges when selecting soil elastic modulus for excavation analysis in gravel layers using different soft-ware programs in the Xindian area of Taiwan to improve the accuracy of deformation predictions during the final excavation stage. These results serve as important references for engineers and contribute to the effective prediction and management of excavation behavior in civil engineering projects.

Creep is a fundamental property naturally exists in some types of rock, which is significant for the long-term stability of roadways during the mining process. In this paper, the long-term strength of coal and rock were obtained via laboratory experiments, and a Cvisc elasto-viscoplastic model was adopted and introduced in FLAC3D, based on the 3101 transport roadway in Hongqinghe Coal Mine, to investigate the influence of creep on the stability of deep-high stress roadway. The test results show that the long-term strength of 3-1 coal and sandy mudstone, respectively, is 18.65 MPa, 39.95 MPa. The plastic zone, the deformation, and the damage of roadway surrounding rock behave an obvious increase after being excavated for 720d as the creep model was chosen. The plastic zone failure with models shear-p (1090.7 m3), shear-n (381.7 m3), tension-n (98.4 m3), and tension-p (30.8 m3), the damage value has an obvious increment of increased 21.2% (0.053), the deformation increases in the order is two-sidewalls (1978 mm), roof (907 mm), floor (101 mm). The creep of the roadway can be divided into three stages: the accelerating stage, the decaying stage, and the stable stage. The creep speed of each stage is greatly affected by the presence or absence of anchor spray support: the creep speed of the bare roadway roof, roadway side, and floor stability, respectively, is 1.01, 1.02, and 0.12 mm/d; After anchor spray support, the creep velocity, correspondingly, decreased to 0.69, 0.37, and 0.12 mm/d, and the amount of surrounding rock damage decreased from 0.302 to 0.243. This indicates that the anchor spray support can significantly reduce the creep effect of the roadway. The Cvisc creep model was verified reliable and can provide guidance for deep high-stress coal roadway support.

The article presents a simplified method for determining the strength of a corrugated board packaging subjected to dynamic transport loads. The proposed algorithm consists of several calculation steps: (1) a static analysis of the compressive strength of the package, (2) an analysis of random vibrations in the frequency domain, used to determine the resonance frequencies, and (3) a dynamic analysis of the package loaded with computed resonant frequencies. For this purpose, numerical models of the static compression test of the packaging before and after the dynamic analysis of the package subjected to general transport loads were developed. In order to validate the model, the laboratory packaging compression tests were also performed for samples of boxes with 3-layer cardboard. Thanks to this, it was possible to verify the results of numerical simulations of compression tests for several box geometries. Which, in turn, allowed for the development of a method based on dynamic and post-dynamic (static) numerical analysis, allowing to determine with high accuracy the resistance of selected packaging to vibrations and dynamic loads. The results of the (validated experimentally) numerical analysis prove the usefulness of the simplified method presented here for precise estimation of the load capacity of various packages dynamically loaded during transport.

In engineering practice, one can often encounter issues related to optimization where the goal is to minimize material consumption, minimize stresses or deflections of the structure. In most cases, these issues are addressed with finite element analysis software and simple optimization algorithms. However, in the case of optimization of certain structures, it is not co straightforward. An example of such constructions are bubble deck ceilings, where in order to reduce the dead weight, air cavities are used, which are regularly arranged over the entire surface of the ceiling. In the case of these slabs, the flexural stiffness is not constant in all its cross-sections, which means that the use of structural finite elements (plate or shell) for static calculations is not possible, and therefore the optimization process becomes more difficult. The paper presents a minimization procedure of the weight of bubble deck slabs using numerical homogenization and sequential quadratic programming with constraints. Homogenization allows to determine the effective stiffnesses of the floor, which in the next step are sequentially corrected by changing the geometrical parameters of the floor and voids in order to achieve the assumed deflection. The presented procedure allows to minimize the use of material in a quick and effective way by automatically determining the optimal parameters describing the geometry of the bubble deck floor cross-section.

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.

New analytical solutions of the heat conduction equation are presented in cylindrical and spherical coordinates. Then these solutions are reproduced with high accuracy by recent explicit and unconditionally stable finite difference methods. After this, real experimental data from the literature regarding a heated cylinder are reproduced by the explicit numerical methods as well as by Finite Element Methods (FEM) ANSYS workbench. Convection and nonlinear radiation are also considered on the boundary of the cylinder, and it is shown that the explicit methods are much more accurate than the commercial software.

In order to study the effect of TiB2 particles on the mechanical properties of TiB2/6061Al composites, a series of TiB2/6061Al 3D representative volume element (RVE) were established based on SEM photos. This model took into account the ductile damage of the matrix and the traction separation behavior of the interface, and the linear damage evolution law was introduced to characterize the stiffness degradation behavior of the matrix elements. The mixed boundary conditions were used for RVE tensile experiments, and the accuracy of the predicted result was verified by the agreement of the experimental stress-strain curve. Results show that the addition of TiB2 particles can effectively promote the load-bearing capacity of the composite, but elongation is reduced. With the weight fraction of TiB2 increasing from 2.5% to 12.5%, the elastic modulus, yield strength, and tensile strength are increased by 8%, 10.37%, and 11.55% respectively, while the elongation decreased by 10%. The clustering rate of the TiB2 particles is also an important factor affecting the toughness of the composites. With the increase in the clustering rate of TiB2 particles from 20% to 80%, the load-bearing capacity of the composites is not improved, and the elongation of the composites is reduced by 8%. Moreover, the high strain region provides a fast expansion path for crack propagation and the particle spacing is a crucial factor that affects the stress field.

In this work, in spite of Milne’s recommendation using the three-point Newton–Cotes open formula (Milne’s rule) as a predictor rule and three-point Newton–Cotes closed formula (Simpson’s rule) as a corrector rule for 4-th differentiable functions with bounded derivatives. There is still a great need to introduce such formulas in other Lp spaces. Often, we need to approximate real integrals under the assumptions of the function involved. Because of that, the aim of this work is to introduce several Lp error estimates for the proposed perturbed Milne’s quadrature rule. Numerical experiments showing that our proposed quadrature rule is better than the classical Milne rule for certain types of functions are provided as well.

Infectious diseases continue to be a significant threat to human health and civilization, and finding effective methods to combat them is crucial. In this research, we investigate the impact of awareness campaigns and optimal control techniques on infectious diseases without proper vaccines. Specifically, we develop an SIRS type mathematical model that incorporates awareness campaigns through media and treatment for disease transmission dynamics and control. The model displays two equilibria, a disease-free equilibrium and an endemic equilibrium, and exhibits Hopf bifurcation when the bifurcation parameter exceeds its critical value, causing a switch in the stability of the system. We also propose an optimal control problem that minimizes the cost of control measures while achieving a desired level of disease control. By applying the minimum principle to the optimal control problem, we obtain analytical and numerical results that show how the infection rate of the disease affects the stability of the system and how awareness campaigns and treatment can maintain system stability. This study highlights the importance of awareness campaigns in controlling infectious diseases and demonstrates the effectiveness of optimal control theory in achieving disease control with minimal cost.

The shield tunnels inevitably endure various damage with the service time increasing. Meanwhile, steel corrugated plate has been used extensively in multiple conditions and proved effective to strengthen the segmental joints by full-scale tests. A numerical model is proposed to probe the feasibility of using new Stainless Steel Corrugated Plate (SSCP) to reinforce the shield tunnel segments. A new method, called virtual tracking element technology, is employed to achieve the simulation of realistic stress state of the segmental joint. Moreover, the segmental joint component analysis and parametric study are conducted based on the numerical model. The results demonstrate that: (1) the virtual tracking element technology proves valid and efficient to simulate the secondary stress state of segmental joint; (2) SSCP reinforcement is not fully utilized when the grade of segmental concrete is C50 and it has abundant safety margin for potential overload; (3) SSCP reinforcement is well-performed regardless of the burial depth and reinforcement in advance is recommended.

With the aim to improve cladding coating quality and prevent cracking, this paper established an extreme high-speed laser cladding thermo-mechanical coupling simulation model to study the evolution of the temperature field and residual stress distribution. Process parameters that impacted the macroscopic morphology of single-pass coatings were investigated. Numerical calculations and temperature field simulations were performed based on the process parameter data to validate the effect of the temperature gradient and cooling rate on the coating structure and the residual stress distribution. The results showed that a good coating quality could be achieved using a laser power of 2400 W, a cladding rate of 20 m/min, and a powder feeding rate of 20.32 g/min. The coatings’ cross-sectional morphology corresponded well with the temperature distribution predicted by numerical modeling of the melt pool. The microstructure of the molten coatings is affected by the temperature gradient and cooling rate, which vary greatly from the bottom to the middle to the top. Maximum residual stress appears between the bonding region of the coatings and the substrate, and the coatings themselves have significant residual stress in the form of tensile strains that are mostly distributed in the direction of the laser cladding speed.

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.

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

Flood routing in flood forecasting issue, calculation the height of flood bands, determining the river boundaries, and estimation of protective facilities for flood –exposed building is applicable. In many cases, due to the lack of measuring stations, the status of the upstream flood generating hydrograph is not known. The purpose of this study is to present an integrated method comprising of an optimization model and a hydrodynamic numerical model for flood modeling to determine the upstream hydrograph using the provided hydrograph at the downstream measuring station of a river. The routing procedure consists of three steps: (1) generating a hypothetical upstream hydrograph using genetic algorithm method; (2) hydrodynamic modeling using a numerical simulation model for flood routing according to the hypothetical hydrograph which is generated in the first step; (3) compare the calculated and observed hydrograph in downstream by using a fitness function. This recommended procedure was named as Reverse Flood Routing Method (RFRM) and was then applied to Karun River, the largest river in Iran. Comparing the generated upstream hydrograph by the RFRM model with the corresponding measured hydrograph at Ahvaz hydrometric station, as an ungauged river location, shows the high accuracy of the recommended model in this study.

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.