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

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.

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

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

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

Ultra-large cold-water pipes (CWP) may present structural instability phenomenon under complex and variable operating conditions, and the current theoretical analysis rarely takes into account the influence of multiple key parameters at the same time, resulting in some errors between the prediction results and the actual deformation behavior of the pipe. Considering this issue, this paper establishes a computational analysis model of the mechanical performance of the CWP based on the Euler-Bernoulli beam theory, considering the effect of waves, the clump weight at the bottom, internal flow and sea current, and gives a time-dependent solution of the mechanical response of the pipe with the action of multiple key parameters using the differential quadrature method (DQM). Additionally, the correctness and validity of the theoretical model are well verified by comparing the numerical solution with the theoretical results. Finally, to study the effects of operating conditions on the mechanical response of the pipe, the effects of waves, the clump weight, internal flow and sea current on the lateral deflection and maximum bending moment of the pipeline are discussed based on numerical simulation and orthogonal experimental method. The results show that changing the wave strength and internal flow velocity has little effect on the lateral deflection of the CWP, and the clump weight at the bottom can effectively suppress the lateral deflection of the pipeline, but increasing the current velocity can significantly increase the lateral deflection of the pipeline, which can lead to the instability of the pipe. Nevertheless, the effects of waves, the clump weight, internal flow and sea current on the maximum bending moment of the pipeline are basically the same, all of them increase sharply at first, and then decrease gradually until they level off, and the differences in the effects are mainly reflected in the different locations of the pipe sections. This paper suggests the design guidance of CWP in aspect of dynamic response depending on the operating condition.

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

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

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

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

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

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

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

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

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

High density of buildings and the large traffic volume of cities give rise to narrow spaces and high passenger flows in most subway stations. When evacuation is required during an emergency, these problems may trigger issues in operational safety. Therefore, it is imperative to comprehensively assess subway station emergency evacuation times and capacities during the operating process to ensure the design is safe, science-based, and rational. To further enhance safe evacuation capacity, this paper proposes an overall evacuation time model that considers the multielement charac-teristics of subway stations in multiple segments. The evacuation route is decomposed into five stages according to the critical nodes in the evacuation process. An overall emergency evacuation time model is established based on the diversity of bottlenecks in the five stages, integrating elements such as the horizontal movement velocity of passengers, subway equipment parameters, and human density. Taking Xi’an Wulukou Subway Station as an example, this paper verifies the outcomes of the theoretical model against the Pathfinder software and conducts additional anal-yses of the evacuation conditions of stairs and exits. The results show that the error between the theoretical emergency evacuation time and simulation evacuation time was 5.4%, The emergency evacuation model established in this study boasts strong robustness and stability.

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

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

The conventional explanation of delayed-choice experiments seems to violate our causal intuition. This apparent violation is caused by a misinterpretation of the conventional formulation of quantum mechanics. I reanalyze these experiments using advanced and time-symmetric formulations of quantum mechanics. All three formulations give the same experimental predictions, but the advanced and time-symmetric formulations violate our causal intuition that effects only happen after causes. I explore reasons why our causal intuition may be wrong at the quantum level. I also suggest how conventional causation might be recovered in the classical limit, and speculate on cosmological boundary conditions.

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

In the present paper four passive micromixers designs (G1, G2, G3 and G4) inspired on distillation columns trays were proposed. The devices performance was assessed by numerical simulations. The mixing performance was investigated for a Reynolds number range from 0.01 to 100 and channel height of 200 µm, 500 µm and 1000 µm for oil/ethanol flow. G1 and G4 designs provided high mixing index (> 0.975). The G1 device achieved superior mixing performance with a moderate pressure drop (about 0.5 MPa) due to the induced flow recirculation pattern for a relative high flow rate of 0.21 L/min, highlighting the potential use of such microdevice for scale-up and numbering-up of microdevices in modular chemical plant processing.

The structure of turbulent flow over non-flat surfaces is a topic of major interest in practical applications in both engineering and geophysical settings. A lot of work has been done in the fully rough regime at high Reynolds numbers where the effect on the outer layer turbulence structure and the resulting friction drag is well documented. It turns out that surface topology plays a significant role on the flow drag especially in the transitional roughness regime and therefore, hard to characterize. Survey of literature shows that roughness function depends on the interaction of roughness height, flow Reynolds number and topology shape. In addition, if the surface topology contains large enough scales then it can impact the outer layer dynamics and in turn modulate the total frictional force. Therefore, it is important to understand the mechanisms underlying drag increase from systematically varied surface undulations in order to better interpret quantifications based on mean statistics such as roughness function. In this study, we explore the mechanisms that modulate the turbulence structure over a two-dimensional (2D) sinusoidal wavy surface with a fixed amplitude, but varying slope. To accomplish this, we model the turbulent flow between two infinitely wide 2D wavy plates at a friction Reynolds number, $Re_{\tau}=180$. We pursue two different but related flavors of analysis. The first one adopts a roughness characterization flavor of such wavy surfaces. The second one focuses on understanding the non-equilibrium near surface turbulence structure and their impact on roughness characterization. Analysis of the different statistical quantifications show strong dependence on wave slope for the roughness function indicating drag increase due to enhanced turbulent stresses resulting from increased production of vertical velocity variance from the surface undulations.

MCrAlX (M: Ni or Co or both, X: minor elements) coatings have been widely used to protect hot components in gas turbines against oxidation and hot corrosion at high temperatures. Understanding the influence of the X-elements on oxidation behaviour is important in the design of durable MCrAlX coatings. In this study, NiCoCrAlX coatings doped with Y+Ru and Ce, respectively, were deposited on Inconel-792 substrate by HVOF. The samples were subjected to isothermal oxidation test in laboratory air at 9000, 1000, 1100 ºC and cyclic oxidation test between 100 ºC and 1100 ºC with 1 hour dwell time at 1100 ºC. It was observed that the coating with Ce shows a much higher oxidation rate than the coating with Y+Ru under both isothermal and cyclic oxidation tests. Besides, β-depletion due to interdiffusion between coating and substrate was significantly lower with Ru addition. Simulations results using both a moving phase boundary model and an established oxidation-diffusion model show that Ru stabilize β grains which reduces β-depletion of coating due to substrate interdiffusion. This paper presents a comprehensive study of the influence of Ce and Ru on oxidation behaviour including investigation of the microstructure evolution influenced by oxidation time at coating surface and coating-substrate interface combining experiment and simulation results.

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

Paper is a material whose mechanical properties are highly dependent on humidity and temperature. Which naturally also builds the relationship between the stiffness and strength of corrugated board and changing weather conditions. In this paper, attention is focused on the dependence of the physical properties of the cardboard on changes in humidity and temperature, which undergo dynamic fluctuations both during the production of corrugated board and during its storage. Two techniques were used to test this effect, namely numerical homogenization and global sensitivity analysis. Both methods were implemented to determine the theoretical relationships between the change in humidity and/or temperature in each layer of corrugated board and its global bending, compression and shear stiffness. The procedure was used to analyze different types of 5-ply and 3-ply cardboard. The obtained results allowed to build a complete map of the relationship between the change in humidity of selected layers and the strength characteristics of the full assembly.

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

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

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

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

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

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

With the popularization of natural gas and the requirements for environmental protection, the development and utilization of natural gas is particularly important. The status of natural gas in China's oil and gas exploration and development is constantly improving, and the country is paying more and more attention to the exploitation and utilization of natural gas. The Upper Paleozoic tight sandstone in the Ordos Basin is characterized by low porosity, low permeability and large area of concealed gas reservoirs. By injecting CO2 into the formation, the recovery of natural gas can be improved, and at the same time, the stable storage of CO2 can be achieved to achieve a win-win situation of CO2 emission reduction and utilization. Injecting greenhouse gas CO2 into gas reservoirs for storage and improving recovery has also become a hot research issue. In order to improve the recovery efficiency of tight sandstone gas reservoir, this paper takes the complex tight sandstone of Upper Paleozoic in Ordos Basin as the research object, through indoor physical simulation experiments, carried out the influence of displacement rate, fracture dip angle, core permeability, core dryness and wetness on CO2 gas displacement efficiency and storage efficiency, and analyzed the influence of different factors on CO2 gas displacement efficiency and storage efficiency to improve the recovery and storage efficiency. The research results show that under different conditions, when the CO2 injection pore volume is less than 1PV, the relationship between the CH4 recovery rate and the CO2 injection pore volume is linear, and the tilt angle is 45 °. When the CO2 injection pore volume exceeds 1PV, the CH4 recovery rate increases slightly with the increase of displacement speed, the recovery rate of CO2 displacement CH4 is between 87% - 97%, and the CO2 breakthrough time is 0.7PV-0.9PV. In low-permeability and low-speed displacement cores, the diffusion of CO2 molecules is more significant. The lower the displacement speed is, the earlier the breakthrough time of CO2 is, and the final recovery of CH4 slightly decreases. Gravity has a great impact on CO2 storage and enhanced recovery. The breakthrough of high injection and low recovery of CO2 is earlier, and the recovery of CH4 is about 3.3% lower than that of low injection and high recovery. The bound water makes the displacement phase CO2 partially dissolved in the formation water, and the CO2 breakthrough lags about 0.1PV. Ultimately, CH4 recovery factor and CO2 storage rate are higher than those of dry core displacement. The research results provide theoretical data support for CO2 injection to improve recovery and storage efficiency in complex tight sandstone gas reservoirs.

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

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

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

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

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

In order to reveal the mechanism and influencing factors of high pressure jet cavitation of nozzle in submerged environment, this study focused on the evolutionary process of cavitation bubbles and combined finite volume method and mixed multi-phase flow model to analyze the cavitation, velocity distribution and experimental cavitation intensity of fishing net cleaning equipment. Results show that the cavitation inception, growth and collapse mainly occurred in the peripheral region of the flow field. Ring-shaped cavitation erosion zone appeared on the test sample target. A lot of small dense erosion pits were densely distributed in the ring-shaped erosion zone, erosion marks were observed in the center. The cavitation erosion intensity was greatly affected by the nozzle structure. As the diameter of nozzle increased from 0.6 mm to 1 mm, the maximum gas volume fraction increased by 8.53%. The nozzle outlet enlargement angle greatly increased the cavitation intensity. The nozzle with an outlet angle of 30° exhibited the optimal cavitation erosion performance . The cavitation volume fraction of the nozzle with short necking structure was slightly larger than that of the nozzle with long necking structure at the same level in the necking length rang of 3mm to 7mm. In terms of the influence of nozzle structure on the cavitation erosion effect, the nozzle diameter D ranked the first, followed by the outlet angle α, and the necking length L was at the last.

The Gleeble-1500D thermal simulation test machine was used to conduct the isothermal compression test on 21-4N at the strain rate ( ) of 0.01-10s-1, the deformation temperature (T) of 1273-1453K and the maximum deformation is 0.916. The data of the stress-strain ( - )were obtained. Based on the - data, the Johnson-Cook (J-C), modified J-C, Arrhenius and Back-Propagation Artificial Neural Network (BP-ANN) models were established. The accuracy of four models were verified, analyzed and compared. The results show that J-C model has a higher accuracy only under reference deformation conditions. When the deformation condition changes greatly, the accuracy of J-C model is significantly reduced. The coupling effect of T and of modified J-C model is considered, and the prediction accuracy is greatly improved The Arrhenius model introduces Zener-Hollomon (Z) to represent the coupling effect of T and , it has a fairly high prediction accuracy. And it can predict flow stress ( ) accurately at different conditions. The accuracy of BP-ANN model is the highest, but its learning rate is low, the learning and memory are unstable. It has no memory for the weights and thresholds of the completed training. So, there are certain limitations of it in use. Finally, a FEM of the isothermal compression experiment for four models were established, and the distribution of the equivalent stress field, equivalent strain field and temperature field with the deformation degree of 60% were obtained.

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

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

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

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

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

The damage evolution rules of coal under true triaxial combined static-dynamic loads are important to understand and reveal the mechanism of rockburst in coal mines. The numerical simulation was carried out using FLAC3D to explore the coupling effect of static and dynamic stress, especially the influence of dynamic load parameters. The results show that the increase of dynamic load amplitude not only makes the plastic zone grow exponentially and shorten the damage development time, but also decreases the residual strength of the coal sample, which leads to a much more violent of the dynamic failure. The stress wave can propagate more stable in the coal sample with the increase of dynamic load frequency and the development time of the plastic zone also show synchronous fluctuating features, and when the frequency is like the intrinsic frequency of the coal sample, the dynamic load has a significant stimulate effect. The static and dynamic loads have different damage effect to the coal sample under true triaxial stress path, the weight decreases in order of the minimum principal stress, the maximum principal, the amplitude of the dynamic load and at last the dynamic stress frequency by the entropy weight method analysis.

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.

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.

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.

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

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

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.

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.

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.

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.

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.

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.

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.

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.

African indigenous herbs and medicinal plants are extensively reported for their protective and therapeutic potential against Diabetes. However, most of the constituent phytochemicals that may account for the therapeutic activity are largely uninvestigated. Exploring novel phytochemical–target interactions in silico could help to provide useful insights into the mechanism underpinning the biological activities of these plants and their constituent phytochemicals. The aim of this study was to explore the interactions of terpene struc-tures previously reported from selected African antidiabetic with two emerging drug targets in diabetes. Structure-based virtual screening was used to screen 107 terpene structures against DPP-4 and PTP1B enzymes. The MD simulation, MM-GBSA free energy calculation and DFT were used to clarify the inter-actions of the in silico hits with the target enzymes. Molecular docking, post-docking Prime MM-GBSA analysis revealed the top terpenes with high binding affinity with the active site regions of DPP-4 and PTP1B. Ensemble docking revealed three triterpenes viz: cucurbitacin B (T1), 6-Oxoisoiguesterin (T4) and 20-Epi-isoiguesterinol (T2) as in silico hits which exhibit strong interaction potential with critical residues that define the catalytic triad (Ser630 and His740); oxyanion cavity (Ser631); hydrophobic S1 pocket (Tyr662) and the charged S2 pocket (Arg125) in the active site region of DPP-4. The interactions of these compounds with DPP-4 exhibited structural stability and conformational flexibility during 100 ns full atomistic MD simulation as indicated by the structure dynamics parameters including RMSD, RoG, SASA and hydrogen bond number. The post-MD MM/GBSA calculations further revealed the stability of the triterpene-DPP-4 complexes. Furthermore, Frontier molecular orbital calculations showed that, the triterpenes possess high interaction potential with the enzyme. The triterpenes showed desirable ADMET properties and drug-likeness. Therefore, cucurbitacin B, 6-Oxoisoiguesterin and 20-Epi-isoiguesterinol are recommended for experimental validation.

Acoustic emission (AE) is caused by the sudden release of energy by a material as a result of material degradation related to deformation, crack, or fault within a solid. The same situation also occurs in leaks caused by turbulence in the fluid around the leak. In this study, analytical modeling for AE due to leakage through a circular pinhole in a gas storage cylinder was performed. The displacement fields responsible for AE, excited by the concentrated force (CF) associated with the turbulent flow though the pinhole, were derived by solving the Navier-Lamé equation. The CF as an excitation source was formulated in terms of a fluctuating Reynolds stress (FRS) and spatial Green’s function. In particular, a series of experiments were conducted under different operating conditions to explore the characteristics of the AE signals due to leak in a gas cylinder. Finally, the simulation and experimental results were compared to verify the accuracy of the simulation results.

For both conventional and renewable energy conversion processes, computational fluid dynamics (CFD) has been used to address more energy-related challenges in recent decades. Using CFD to investigate vertical-axis wind turbines has become more common in recent years. The main goals of this application have been to more accurately predict the turbine's performance and to comprehend the complicated nature of the complex turbulent flow. The Vertical Axis Wind Turbine (VAWT) simulation for energy-generating applications has several intricate components. One of them is being aware of the chaotic flow that occurs during the first stages of the starting process, which influences overall effectiveness. In this article, the performance of the wind turbine was increased by using a passive flow control approach. The numerical research was carried out using Large Eddy Simulation for four alternative tip speed ratios in both cases, the classic and the optimized, equipped with a vortex trap on the exterior of the blades. The power and torque coefficient variations, as well as the velocity magnitude contours, show that the starting process may begin with a significant improvement in efficiency when flow control is used.

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

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

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.

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.

There are ~200.000 households living in the detached houses and ger (yurt) with the small-coal stove which burns a raw coal in the Ulaanbaatar city. The proper heating system and improvement on the energy efficiency of residential dwellings is vital important for Ulaanbaatar city to reduce an air pollution as well as for the operation of current central energy system. This study shows the experimental results of two Gers with two different heating system and different thermal insulation. Technical availability of the system consisting of ETS heater with daytime charging schedule and areal PV system is also examined by using a simulation of MATLAB programming. As a result of the experiment, the indoor comfort level and energy efficiency of the Ger with added insulation and ETS heater with nighttime charging was enhanced to compare with the reference Ger. The Ger with added insulation and ETS heater consumed 3169 kWh for electric appliances and 5989 kWh for the heating season. The simulation was resulted that the ETS heater to be installed in the Ger 2 had prospected to increase the PV self-consumption rate by 76% because of the daytime charging schedule of ETS heater. The PV system supplied 31% of the total energy consumption with remained 69% of the main grid.

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.

Mathematical modeling is the most important tool for constructing the theory of different kinds of two-phase flows. The review is devoted to the analysis of the advent of the approach of mathematical modeling of two-phase flows, where solid particles act mainly as the dispersed phase. The main problems and features of the study of gas-solid two-phase flows are excluded. The main characteristics of gas flows with solid particles are given and the classification of two-phase flows is developed on their basis. The Lagrangian and Euler approaches to modeling the motion of a dispersed phase (particles) are described. Much attention is paid to the consideration of numerical simulation methods that describe descriptions of turbulent gas flow at different hierarchical levels (RANS, LES, DNS), different levels of description of interphase interactions (one-way coupling (OWC), two-way coupling (TWC) and four-way coupling (FWC)), as well as at different levels of interface resolution (partial-point (PP) and particle-resolved (PR)). Excluded are examples of some studies carried out on the basis of the identified approaches, as well as their exclusion for mathematical modeling of various classes of gas-solid two-phase flows.

Agricultural drainage canals that connect upstream fish spawning areas to downstream rivers and lakes serve as crucial habitats for migrating fish. However, disconnections, such as drops and chutes, have been constructed in these canals due to agricultural modernization and flood control measures, hindering the movement of fish with limited swimming abilities. Portable fishways offer a promising solution to reconnect water bodies in agricultural canals, as they can be easily removed during high water discharges to avoid impeding the canals' drainage function. In addition to experimental assessments of fishway functionality, employing a hydrodynamic model to explore effective placement strategies for portable fishways is essential to maximize their effectiveness. This study presents a method to determine the best horizontal location of a portable fishway in an agricultural drainage canal using two-dimensional hydrodynamic simulations within the specified cases. The applicability of this method is demonstrated by addressing the positioning challenge of a portable fishway on a chute in an agricultural drainage canal in Japan. Results indicate that the proposed method allows for the selection of a suitable location, considering preferable hydraulic conditions both within the portable fishway and around its entrance.

The monkeypox virus (MPXV) has become a major threat due to the increasing global caseload and the ongoing multi-country outbreak in non-endemic territories. Due to limited research in this avenue and the lack of intervention strategies, the present study was aimed to virtually screen bioactive phytochemicals against envelope proteins of MPXV via rigorous computational approaches. Molecular docking and molecular dynamic (MD) simulations were used to investigate the binding affinity of 12 phytochemicals against three envelope proteins of MPXV, viz., D13, A26, and H3. Silibinin, oleanolic acid, and ursolic acid were computationally identified as potential phytochemicals that showed strong binding affinity towards all the tested structural proteins of MPXV through molecular docking. The stability of the docked complexes was also confirmed by MD simulations. ADME analysis also computationally confirmed the drug-like properties of the phytochemicals, thereby asserting their suitability for consumption. Hence, this study envisions the candidature of bioactive phytochemicals as promising inhibitors against the MPXV, serving as template molecules that could further be experimentally evaluated for their efficacy against monkeypox.

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.

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

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

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

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

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.

Subsurface temperature is the key parameter in geothermal exploration. An accurate estimation of the reservoir temperature is of high importance and usually done either by interpolation of borehole temperature measurement data or numerical modeling. However, temperature measurements at depths which are of interest for deep geothermal applications (usually deeper than 2 km) are generally sparse. A pure interpolation of such sparse data always involves large uncertainties and usually neglects knowledge of the 3D reservoir geometry or the rock and reservoir properties governing the heat transport. Classical numerical modeling approaches at regional scale usually only include conductive heat transport and do not reflect thermal anomalies along faults created by convective transport. These thermal anomalies however are usually the target of geothermal exploitation. Kriging with trend does allow including secondary data to improve the interpolation of the primary one. Using this approach temperature measurements of depths larger than 1,000 m of the federal state of Hessen/Germany have been interpolated in 3D. A 3D numerical conductive temperature model was used as secondary information. This way the interpolation result reflects thermal anomalies detected by direct temperature measurements as well as the geological structure. This results in a considerable quality increase of the subsurface temperature estimation.

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

The paper presents the prediction of hydrodynamic characteristics of different rudder -plane configurations behind a submarine in deep water through the Reynolds-Averaged Navier-Stokes method in Ansys Fluent Solver. The hydrodynamic force and moment acting on the Cross-plane, X-plane, and Y-plane rudder configurations are analyzed for estimating hydrodynamic coefficientsfor predicting the advantages of each configuration to control the ship. The obtained results show that X-force of the X-plane rudder configuration is larger than the force acting on the Cross-plane rudder and Y-plane rudder configurations. However, Y-force, and Z-force of the X-plane rudder configuration are significantly greater than the Y-force of the left configurations. The same tendency can be seen in the moment of the X-plane rudder about the y-axis and z-axis. The roll moment induced by the Y-plane and Cross-plane rudder configurations is significantly larger than the case of the X-plane rudder configuration.

The research work explores the impact of temperature on Silicon photovoltaic (PV) panels considering Nigeria as case study. It is found that high solar radiation in Nigeria increases surface temperature of PV panel above 25oC of optimal operating temperature of PV panels. Redundant energy gain from incident solar energy is dissipated at the rear of PV panels as heat, which reduces the efficiency of PV panel. Cooling mechanism is needed to cool the PV panels to increase the efficiency. In this study we demonstrated a unique design of a bifunctional photovoltaic-thermal (PVT) system by employing heat exchanger at the rear of PV panel and water is used as a working fluid that is circulated through the heat exchanger to remove the excess heat. It is observed from the simulation results that Maiduguri peak annual electrical power output gives 1907kWh/kWp which is the highest, because of high solar radiation average of 727W/m2 in a year. And for Makurdi the peak annual electrical power output of 1542kWh/kWp while for Port Harcourt the peak power output of 1355kWh/kWp respectively. Interestingly, Polycrystalline Si-PV surface temperature decrease from 49.25oC to 38.38oC. Electrical power increased from 1526.83W to 1566.82W in a day and efficiency increase from 13.99% to 15.01%.

This research delves into the impact of complex low-speed wake flow structures on the aerodynamic characteristics of objects downstream. It employs the proper orthogonal decomposition (POD) method and the domain precursor simulation method to compare traditional methods and validate this approach. The study generates several flow structures of parallel dual-cylinder wakes with different scales and spacing. The variations in the aerodynamic coefficient of three downstream objects at various times passing through wakes of varying scales are appropriately compared and analyzed. The study establishes that the wake with a cylinder spacing of G=1.5 has a more compact and concentrated modal structure than that with a cylinder spacing of G=0.35. Smaller objects are more responsive to the wake flow structure with a spacing of G=1.5, whereas larger objects respond more to the flow structure with a spacing of G=0.35. The achieved results also reveal that the aerodynamic force coefficients of objects passing through the wakefield at different times are closely related to the temporal characteristics of the wake flow structure with different scales.

The current work investigates the auxetic tensile deformation behavior of the inversehoneycomb structure with 5 × 5 cells made of biodegradable poly(butylene adipate-coterephthalate) (PBAT). Fused deposition modeling, an additive manufacturing method, produced such specimens. Residual stress (RS) and warpage, more or less, always exist in such specimens due to layer-by-layer fabrication, i.e., repeated heating and cooling. The RS influences the auxetic deformation behavior, but its measurement is challenging due to the very fine structure. Instead, the finite-element(FE)-based process simulation realized by an ABAQUS plug-in numerically predicts the RS and warpage. The predicted warpage shows a negligible slight deviation compared to the design topology. This process simulation also delivers the temperature evolution of a small volume material, revealing the local cyclic heating and cooling. The achieved RS serves as the initial condition for the FE model used to investigate the auxetic tensile behavior. With the outcomes from FE calculation without considering RS at hand, the effect of the RS on the deformation behavior is discussed for the global force-displacement curve, the structural Poisson’s ratio evolution, the deformed structural status, the stress distribution, and evolution, where the first three and the warpage are also compared with experimental results. Furthermore, the FE simulation can easily provide the global stress-strain flow curve with the total stress calculated from the elemental ones.

Keyhole laser beam welding (LBW) of stainless steel sheets with a gap in between is numerically simulated with a three-dimensional, transient multi-physical model for laser material processing. At first, the model’s ability to reproduce experimental results on a relatively coarse computational mesh within reasonable computing time, so as to serve as process optimization tool, is presented. An example of process optimization, where a given set of weld seam quality criteria is fulfilled by iteratively optimizing a secondary laser beam, is shown. The relatively coarse mesh, in combination with a good model calibration to the experimental conditions, allows for sufficiently fast simulations to use this approach for optimization tasks. Finally, using a finer spatial and temporal discretization, the dynamic processes in the vicinity of the keyhole leading to the formation of pores are investigated. The physical phenomena predicted by the simulation are coherent with experimental observations found in literature.

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

In this paper, possibilities for the bandwidth and wavelength utilization in future next-generation passive optical networks are presented and reasons for realizing and utilizing of extended dynamic wavelength and bandwidth algorithms for NG-PON2 networks are analyzed. Next, principles of the effective dynamic bandwidth allocation are introduced in details, focused on the importance of the decision criterion optimization. For achieving better bandwidth utilization of dedicated wavelengths, this paper is focused on the novel effective dynamic bandwidth allocation algorithm with adaptive allocation of wavelengths to ONU units and on the optimization of the decision criterion. For analyzing the novel extended DWBA algorithm utilizing for various wavelength allocation cases possible in NG-PON2 networks, a simulation program with the enhancement of the EDBA algorithm is realized. Finally, an optimization of the decision criterion defining a minimum bandwidth utilization of the actual wavelength is executed for NG-PON2 networks based on the hybrid TWDM technique.