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

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

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

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.

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.

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

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

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.

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.

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

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

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.

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

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.

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

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

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

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

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

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

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.

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

Present study offers mathematical calculations of the roll-coating procedure lubricated with upper Convected Maxwell Fluid. An incompressible isothermal viscoelastic fluid is considered with roll and the porous web having uniform velocities. By employing Lubrication Approximation Theory, the desired equations of motion for the fluid concerned over porous web are modelled and analyzed. The suction rate on the web and injection rate at the roll surface are anticipated proportionate. Results for velocity profile and pressure gradient are obtained analytically. Fluid parameters of industrial significance i.e. detachment point, pressure, sheet /roll separating force, power contribution and coating thickness are also calculated numerically. Substantial and monotonic increase is witnessed in these quantities with the increase of flow parameters.

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

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.

In the present trend of constructing taller and longer structures, the application of lightweight aggregate concrete is becoming an increasing important advanced solution in the modern construction industry. In engineering practice, the analysis of lightweight concrete elements is performed using the same algorithms used for normal concrete elements. As an alternative to traditional engineering methods, nonlinear numerical algorithms based on constitutive material models may be used. The paper presents a comparative analysis of curvature calculations for flexural lightweight concrete elements, incorporating analytical code methods EN 1992-1 and ACI 318-14, as well as a numerical analysis using the constitutive model of cracked tensile lightweight concrete recently proposed by the authors. To evaluate the adequacy of the theoretical predictions, experimental data of 51 lightweight concrete beams tested during five different programmes were collected. A comparison of theoretical and experimental results showed that the most accurate predictions are obtained using numerical analysis and the constitutive model proposed by the authors. In the future, the latter algorithm can be used as a reliable tool for improving the design standard methods or numerical modelling of lightweight concrete elements subjected to short-term loading.

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

In the modern construction industry, lightweight aggregate concrete (LWAC) is often used in the production of load-bearing structural members. LWAC can be up to 40% lighter by volume in comparison to normal strength concrete. On the other hand, the lack of adequate numerical models often limits the practical application of innovative building materials, such as lightweight concrete, in real projects. This trend is due to the uncertainties in design standard methods and calculation errors, the level of which is generally unacceptable to civil engineers in terms of safety and reliability. In the present paper, a comparative numerical deformation analysis of a full-scale bridge deck slab and girder has been carried out. Using the physical model proposed by the authors and the finite element software ATENA, the deformations of full–scale lightweight and traditional reinforced concrete elements under short-term effects of permanent and variable loads was compared. Depending on the safety and serviceability limit requirements, it was found that the amount of longitudinal reinforcement in lightweight reinforced concrete elements can be reduced compared to normal reinforced concrete elements with the same parameters. The results of the numerical analysis show that the deformation analysis model proposed by the authors can be a reliable tool for the design of lightweight concrete flexural members by selecting the optimum geometrical and reinforcement parameters limited by the stiffness condition.

Non steady applied magnetic field impact on a liquid metals has good prospects for industry. For a better understanding of heat and mass transfer processes under these circumstances, numerical simulations are needed. A combination of finite elements and volumes methods was used to calculate the flow and solidification of liquid metal under electromagnetic influence. Validation of numerical results was carried out by means of measuring with ultrasound Doppler velocimetry technique, as well as with neutron radiography snapshots of the position and shape of the solid/liquid interface. As a result of the first part of the work, a numerical model of electromagnetic stirring and solidification was developed and validated. This model could be an effective tool for analyzing the electromagnetic stirring during the solidification process. In the second part, the dependences of the velocity pulsation amplitude and the melt velocity maximum value on the magnetic field pulsation frequency are obtained. It was found numerically the ability of the pulsating force action to develop higher values of the liquid metal velocity at a frequency close to the MHD resonance. Obtained characteristics give a more detailed description of the electrically conductive liquid behaviour under action of pulsating traveling magnetic field.

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

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

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

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

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

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

From June 18 to 19, 2013, the Ésera river in the Pyrenees, Northern Spain, caused widespread damage due to flooding as a result of torrential rains and sustained snowmelt. We estimate the contribution of snow melt to total discharge applying a snow energy balance to the catchment. Precipitation is derived from sparse local measurements and the WRF-ARW model over three nested domains, down to a grid cell size of 2 km. Temperature profiles, precipitation and precipitation gradient are well simulated, although with a possible displacement regarding the observations. Snowpack melting was correctly reproduced and verified in three instrumented sites, and according to satellite images. We found that the hydrological simulations agree well with measured discharge. Snowmelt represented 33% of total runoff during the main flood event and 23% at peak flow. The snow energy balance model indicates that most of the energy for snow melt during the day of maximum precipitation came from turbulent fluxes. This approach forecast correctly peak flow and discharge during normal conditions at least 24h in advance and could give an early warning of the extreme event 2.5 days before.

In this work the existence and nonexistence of stationary radial solutions to the elliptic partial differential equation arising in the molecular beam epitaxy are studied. Since we are interested in radial solutions we arrive at the following fourth-order differential equation

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

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

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

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

Alzheimer's disease (AD) is the first cause of dementia worldwide and represents a public health challenge. Current diagnostic methods still rely on extended interviews and paper tests. We aim to create a novel, quick cognitive-screening tool on a numerical tablet. This program, built and edited with Unity®, runs on Android® for the Samsung Galaxy Tab S7 FE®. Composed of seven tasks inspired by the Mini-Mental Status Examination and the Montréal Cognitive Assessment, it browses several cognitive functions. The architectural design of this tablet application is distinguished by its multifaceted capabilities, encompassing not only seamless offline functionality but also a mechanism to ensure the singularity of data amalgamated from diverse sites. Additionally, a paramount emphasis is placed on safeguarding the confidentiality of patient information in the healthcare domain. Furthermore, the application empowers individual site managers by allowing them to access and peruse specific datasets, thereby enhancing their operational efficacy and decision-making processes. We performed a usability assessment among 24 healthy patients with a final F-SUS score of "excellent". Participants perceived the tool as simple to use and achieved the test in a mean time of 142 seconds, confirming that a short assessment on a numerical tablet is possible.

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

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

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

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

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

Abstract—In the present paper the static and fatigue analysis of plates with holes is conducted. Structures are made of 2D woven roving composites. The parametrization of convex holes is proposed. Experimental results of specimens without holes and with different shapes of notches under static and fatigue tensile loads are presented. The fatigue behaviour of structures is analysed with the use of the low cycle fatigue method. The attention is focused on the nonlinear deformations of 2D woven roving structures leading to the scatter of results. The analysis is supplemented by the numerical finite element modeling. The agreement between experimental and finite element modeling is very good. Since we do not intend to repeat in the abstract the review of the existing literature I would like to emphasize only that the present work is a continuation and extension of the results discussed in the literature – see the next sections of the paper. The finite element modeling of plates with holes is presented in details by Muc [1], Muc, Ulatowska [2] and Muc, Romanowicz [3] including also the problems of the accuracy. The mesh correctness and accuracy is discussed in Ref [2] where the results are compared to the analytical results illustrated in the Soviet Union literature. The experimental results and the form of specimens are shown by Pandita et al. [4] and they are fundaments for our present experimental investigations. Two parametrical description of static and fatigue behaviour of plates with holes characterizes the goal and the novelty of the paper. The analysis demonstrates how those parameters affect on the variations of static and fatigue failure loads. In addition, it should be emphasized that the fuzzy approach is employed herein to reduce the number of experimental data and to avoid statistical analysis of experimental results – see e.g. Muc [5]. In the mesoscale modeling and using the fuzzy approach it is possible to achieve a very good description of static and fatigue behavior of plates with convex holes subjected to tension.

In this work, new refinements of some numerical radius inequalities are proved. Namely, new improvements and refinements purify the recent inequalities of some famous inequalities concerning the numerical radius of Hilbert space operators. The proven inequalities in this work are not only an improvement over old inequalities, but rather they are stronger than them. Several examples that support the validity of our results are established as well.

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

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

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.

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

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

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

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

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

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.

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

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.

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

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

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.

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

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

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

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

We investigated preschool-aged children’s understanding of early fractional tasks and how that performance correlates with fine motor skills and use of gestures while counting. Participants were 33 preschoolers aged 4 to 5 in two Southeastern public elementary schools. Children were tested individually in an interview-like setting. Mathematics tasks were presented in a paper and pencil format and the Grooved Pegboard test assessed fine motor skills. Finally, utilization of gestures was evaluated by taking a behavioral rating of the child’s hand morphology, accuracy of gestures, and synchrony of gestures and spoken word while performing a counting task. Results indicate performance on fractional reasoning tasks significantly predicts both fine motor ability and accuracy of gestures.

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.

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

In this work, some numerical radius inequalities based on the recent Dragomir extension of Furuta's inequality are obtained. Some particular cases are also provided.

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.

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

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

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

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.

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.

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

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

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

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

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

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

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.

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

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

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

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

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

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

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