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.

Motivated by the way animals react in nature, we introduce a fully-connected neural-network model, which incorporates the external pattern presented to the system as a fundamental tool of the recognition process. In this new scenario, in the absence of this external field, memories are not attractors inside basins of attraction, as in usual attractor neural networks, although basins may be created according to the external pattern, thus allowing storing a much larger number of memories. The key point consists in calibrating the influence of the external pattern such as to cancel the noise generated by those memories not correlated with the external pattern.We illustrate how this proposal works by including this new contribution in the standard Hopfield model, showing a significant increase in its recognition capacity (typically by a factor $10^{2}$). As an additional feature of this model, we show its ability to react promptly to changes in the external environment, a crucial attribute of living beings. This procedure can be applied to a wide variety of neural-network models to increase considerably their recognition and reaction capabilities.

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

FRAP technique have been used for decades to measure movements of molecules in 2D. Data obtained by FRAP experiments in cell plasma membranes are assumed to be described through means of two parameters, a diffusion coefficient D (as defined in a pure Brownian model) and a mobile fraction M. Nevertheless, it has also been shown that recoveries can be nicely fit using anomalous sub-diffusion. FRAP at variable radii has been developed using the Brownian diffusion model to access geometrical characteristics of the surrounding landscape of the molecule. Here we performed numerical simulations of continuous time random walk (CTRW) anomalous subdiffusion and interpreted them in the context of variable radii FRAP. These simulations were compared to experimental data obtained at variable radii on living cells using the PH domain of the membrane binding protein EFA6 (exchange factor for ARF6, a small G protein). This protein domain is an excellent candidate to explore the structure of the interface between cytosol and plasma membrane in cells. By direct comparison of our numerical simulations to the experiments, we show that this protein does not exhibit anomalous diffusion in BHK cells. The non Brownian PH-EFA6 dynamics observed here is more related to spatial heterogeneities such as cytoskeleton fences effects.

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

The theory says that number interacts with spaces based on the meaning of mathematics, in this way it is possible to obtain numerical knowledge. The theory shows the possibility of using spaces related to each other to efficiently discover whether results and origins are even or odd, but always applying rules to make this work. the study uses equivalent representations of odd and even numbers. Thus, it is possible to discover the possibilities of origins and their results.

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

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

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

Wireline Formation Testing is an important technique in the oil field of exploration and development. Not only it can make real fluid samples of the formation directly obtained to know exactly oil existed in formation or not; but it can also show flowing pressure change to see the production capacity of the formation. So it is important measurement for formation evaluation during the drilling process and supports for activities of the exploration and development in the oil field. A numerical simulation model in this article is researched and established based on the Finite Volume Method with considering the influence of sensitive parameters such as reservoir heterogeneity, probe suction area, and mud-filtrate invasion depth during the drilling. The model is capable to design and evaluate for formation fluid sampling operations with calculating hydrocarbon content and flowing pressure. Furthermore, through case application, the performance and effect in the process of Wireline Formation Testing was demonstrated. The results indicate that this technology can be served as an effective auxiliary tool for fluid sampling operations with the function of optimizing fluid sampling measures. It can improve the accuracy of predicting indicators such as hydrocarbon content and breakthrough time during the sampling process. This study provides important backups and technical guidance for professionals in geological exploration and oilfield development.

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

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

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

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.

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

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

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

In the case of solid slabs made of reinforced concrete that are usually subjected to bending, large areas of the structure are stressed well below their load-bearing capacity. 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.

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

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

Fractured media are seldom considered viable for storing excess thermal energy seasonally, compared to the more commonly used porous aquifers, due to their heterogeneity, lower storage capacity, and potential for unpredictable flow paths. This study employs COMSOL Multiphysics to create a numerical model of a fractured granite reservoir located below the Bedretto underground laboratory. The energy efficiency of both single-well and doublet systems was evaluated over 25 cycles, each comprising four 91-day phases: injection, storage, production, and resting. Results indicated that the single-well system was more efficient, with energy efficiency increasing with the flow rate. For an injection temperature of 60 °C, the system achieved energy efficiencies of up to 59.3% and 67.7%, after 10 and 25 cycles, respectively, aligning with porous Aquifer Thermal Energy Storage (ATES) system efficiencies at this temperature. Exploration of electricity production using water injected at 120 °C revealed it to be uneconomical and ineffective. However, the most promising use of the retrieved water was found to be direct ground heating for a greenhouse. This technology offers a potential alternative for thermal energy storage in regions lacking suitable porous aquifers, although identifying appropriate locations with the necessary geological characteristics may present significant challenges.

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

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

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

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

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

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

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

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

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.

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

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

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

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

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.

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

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

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

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

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.

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

Liquid propellants are fast becoming attractive for pulsed plasma thrusters due to their high efficiency and low contamination issues. However, the complete plasma interaction and acceleration processes are still not very clear. Present paper develops a multi-layer numerical model for liquid propellant PPTs. The model proposes a possible acceleration mechanism for liquid fed pulsed plasma thrusters and accurately predicts the propellant utilization capabilities and estimations for the fraction of propellant gas that is completely ionized and accelerated to high exit velocities. Validation of the numerical model and the assumptions on which the model is based on is achieved by comparing the experimental results from two different liquid-fed thrusters developed at the University of Tokyo. Results show 50% of liquid propellant injected is completely ionized and accelerated to high exit velocities (> 50Km/s), whereas, neutral gas contribute to only 7% of the total specific impulse and accelerated to low exit velocity (< 4Km/s). The model shows an accuracy up-to 92%. Optimization methods are briefly discussed to ensure efficient propellant utilization and performance. The model acts as a tool to understand the background physics and to optimize the performance for liquid-fed PPTs.

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

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

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

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

We are concerned with the use of some classical spectral collocation methods as well as with the new software system Chebfun in order to compute high order (index) eigenpairs of singular as well as regular Schrodinger eigenproblems. We want to highlight both the qualities as well as the shortcomings of these methods and evaluate them vis-a-vis the usual ones. In order to resolve a boundary singularity we use Chebfun with the simple domain truncation technique. Although this method is equally easy to apply with spectral collocation, things are more nuanced in the case of these methods. A special technique to introduce boundary conditions as well as a coordinate transform which maps an unbounded domain to a nite one are the ingredients. A challenging set of "hard" benchmark problems, for which usual numerical methods (f. d., f. e. m., shooting etc.) fail, are analysed. In order to separate "good"and "bad"eigenvalues we estimate the drift of the set of eigenvalues of interest with respect to the order of approximation and/or scaling of domain parameter. It automatically provides us with a measure of the error within which the eigenvalues are computed and a hint on numerical stability. We pay a particular attention to problems with almost multiple eigenvalues as well as for problems with a mixed (continuous) spectrum. In the latter case we try to numerically highlight its existence. Special attention will be paid to the higher eigenpairs (the pair of eigenvalue and the corresponding eigenfunction approximated by an eigenvector spanning its nodal values).

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

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

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

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

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

Distributions with bounded support show considerable sparsity over those with unbounded support, despite the fact that there are a number of real-world contexts where observations take values from a bounded range (proportions, percentages, and fractions are typical examples). For proportion modelling, a flexible family of two-parameter distribution functions associated with the exponential distribution is proposed here. Mathematical and statistical properties of the novel distribution are examined, including quantiles, mode, moments, hazard rate function, and its characterization. The parameter estimation procedure using the maximum likelihood method was carried out, and applications to environmental and engineering data were also considered. To this end, various statistical tests are used, along with some other information criterion indicators to determine how well the model fits the data. The proposed model is found to be the most efficient plan in most cases for the datasets considered.

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

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

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

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

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

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

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

Based on the geologic conditions of natural gas hydrate (NGH) accumulations in the Shenhu area, northern slope of the South China Sea, this paper used a method of combining reservoir stimulation technique (RST) with depressurization to investigate the conversion efficiency of marine NGH reservoirs in different intrinsic permeability and initial NGH saturation conditions, and analyze the influence of variably-stimulation effect on marine NGH conversion efficiency in different accumulation conditions, provided a reference scheme for improving the NGH conversion efficiency in the Shenhu area. In this work, we performed calculations for the variations in CH4 production rate and cumulative volume of CH4 in different initial NGH saturation, intrinsic permeability and stimulation effect conditions, variance analysis and range analysis methods were used to analyze the significance of these key factors and their interaction, and investigated the sensitivity of stimulation effect on NGH conversion efficiency, respectively. The simulation results showed that stimulation effect have a significant influence on NGH conversion efficiency, and the influence of interaction between these three factors were not obvious. Possibly most important, we clarified an optimum stimulation effect for higher NGH conversion efficiency under different accumulation conditions, especially in low-permeability and high-saturation, which corresponding stimulation effect were different.

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.

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.

Recently, there has been increasing attention toward inhaled nanoparticles (NPs) to develop in-halation therapies for diseases associated with the pulmonary system and investigate the toxic ef-fects of hazardous environmental particles on human lung health. Taking advantage of microfluidic technology for cell culture applications, lung-on-a-chip devices with great potential in replicating the lung air-blood barrier (ABB) have opened new research insights in preclinical pathology and therapeutic studies associated with aerosol NPs. Surprisingly, the air interface in such devices, which makes them tremendously unique among the other common liquid-based cellular mi-crosystems and enables exposure of cultivated cells to the drug/toxic aerosols, has been largely disregarded. Accordingly, there exists a significant research gap in the comprehensive under-standing of the NPs’ dynamics in the lung-on-a-chip devices, which is momentous to control the injection of aerosols to the target area and attain a desired physiochemical efficacy. To address this research gap, a numerical parametric study is presented to provide a deep insight into the dynamic behavior of the airborne NPs in a gas-liquid dual-channel lung-on-a-chip device with a porous membrane separating the channels. A finite element multi-physics model is developed to enable a particle tracing investigation simultaneously in both air and medium phases to replicate the in vivo conditions avoiding laborious and costly experimental trials. The results include the impact of fluid flow as well as geometrical properties on the distribution, deposition, and translocation of the NPs with diameters ranging from 10 nm to 900 nm in the lung-on-a-chip. The obtained results strongly suggest the aerosol injection of NPs instead of the aqueous solution for more efficient deposition on the substrate of the air channel and higher translocation to the media channel. The current study proposes to optimize the affecting parameters to control the injection and delivery of aerosol par-ticles into the lung-on-chip device depending on the objectives of biomedical investigations and provides optimized values for some specific cases. Therefore, this study can assist scientists and researchers in complementing their experimental investigation in future preclinical studies on pulmonary pathology associated with inhaled hazardous and toxic environmental particles, as well as therapeutic studies for developing inhalation drug delivery.

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.

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

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

Solid asperity interactions are common and inevitable under severe loading conditions for any lubricated contact. Heavy-duty machine components (gear, bearings etc.) generally operate under Mixed Lubrication, where uneven surface features contact each other when generated fluid pressure is not enough to support the external load. The Reynolds equation is commonly used to simulate smooth lubricated contacts numerically. In rough lubricated interfaces where opposite surface asperities make contact, the Reynolds equation alone cannot accurately predict pressure using the traditional numerical simulation method. In this paper, Lubrication-Contact interface conditions (LCICs) have been implemented and extended to solve the multiple asperity contact problem using the full-multigrid approach. The developed novel algorithm has significantly accelerated the solution process and improved the accuracy and efficiency of pressure calculation for fluid-solid sub-interactions that can occur under ML regions. The Finite Difference Method (FDM) results have been compared with Computational Fluid Dynamics (CFD) simulation to validate the newly developed model. Hence, the proposed optimized solution method will provide valuable insight to researchers and industry engineers interested in simulating the ML problem where the effect of the fluid-solid interface can be captured effectively to improve reliability in the calculation of the life expectancy of the lubricated parts.

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

Keywords: underground engineering; numerical simulation; excavation length effect, major principal stress; displacement; damage initiation; CPU time

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

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

In this paper we introduce the hypo-q-norms on a Cartesian product of algebras of bounded linear operators on Hilbert spaces. A representation of these norms in terms of inner products, the equivalence with the q-norms on a Cartesian product and some reverse inequalities obtained via the scalar reverses of Cauchy-Buniakowski-Schwarz inequality are also given. Several bounds for the norms δ_p, ϑ_p and the real norms η_r,p and θ_r,p are provided as well.

As long as the non-contact digital printing is not a common standard in the corrugated packaging industry, corrugated board crushing is a real issue that affects the load capacity of the boxes. Crushing mainly occurs during the converting of corrugated board (e.g. analog flexographic printing or laminating) and is a process that cannot be avoided. However, as show in this study, it can be controlled. In this work, extended laboratory tests were carried out on the crushing of double-walled corrugated board. The influence of fully controlled crushing (with a precision: ±10 μm) in the range from 10 to 70 % on different laboratory measurements was checked. Most of the typical mechanical tests were performed e.g. edge crush test, four-point bending test, shear stiffness test, torsional stiffness test, etc. on reference and crushed specimens. The residual thickness reduction of the crushed samples was also controlled. All empirical observations and performed measurements were the basis for building an analytical model of crushed corrugated board. The proven and verified model was then used to study the crushing effect of the selected corrugated board on the efficiency of simple packages with various dimensions.

Predictive simulation of concurrent debris flow using only pre-disaster information has proven to be difficult as a result of problems in predicting the location of debris-flow initiation (i.e., slope failure). However, because catchment topography has concave characteristics, with all channels in a catchment joining each other as they flow downstream, it is possible to predict damage to downstream area using relatively inaccurate initiation points. Based on this, this paper presents methodologies employing debris-flow initiation points generated randomly using statistical slope failure prediction. A many-case simulation across numerous initiation points was performed to quantify the effect of slope-failure location in terms of deviations in the predicted water level and terrain deformation. It was found that the relative standard deviation diminished as the points approached the downstream area, indicating a location-based predictability effect.

Biological diversity (BD) explored by the biological systematics is a complexly organized natural phenomenon and can be partitioned in several aspects defined with references to various causal factors structuring biota. These BD aspects are studied by particular research programs based on specific taxonomic theories (TT). They provide in total a framework for comprehending the structure of the biological systematics and its multi-aspect relations to other fields of biology. General principles of individualizing BD aspects and construing TT as quasi-axiomatics are briefly considered. It stressed that each TT is characterized by a specific combination of interrelated ontological and epistemological premises most adequate to the BD aspect a TT deals with. The following contemporary research programs in systematics are recognized and characterized in brief: phenetic, rational (with several subprograms), numerical, typological (with several subprograms), biosystematic, biomorphic, phylogenetic (with several subprograms), evo-devo. From a scientific pluralism perspective, all these research programs related to particular naturally defined BD aspects are of the same biological and scientific significance and no one of them can pretend to take a privileged position. They elaborate “locally” natural classifications that can be united by a kind of generalized faceted classification.

Rigid-hulled inflatable boats are extremely practical and popular nowadays. They offer a effective conciliation among usability and costs. Their stable and seaworthy behaviour is guaranteed by performing hydroplaning hulls coupled with unsinkable inflated tubes. At the same time, their design is often based on tradition and preconceptions. Rarely, the design assumptions are validated by the reality or, even, by deeper investigations. In this article, both numerical methods and experimental mechanics techniques are proposed as an essential way for supporting the designers in their decisive tasks. Three different situations are detailed where a numerical or an experimental approach shows its benefit inside the engineering design process: firstly permitting to investigate the behaviour of materials driving the fiberglass selection; then measuring the levels of stress and strain in the hull during sailing; finally, using all available information as a base for developing numerical models of the hull slamming in waves. Even if the discussion is focused on a rigid inflatable boat, large part of its considerations is relevant beyond this special case.

This study reports the development of a numerical method to simulate two-phase flows of Newtonian fluids that are incompressible, immiscible, and isothermal. The interface in the simulation was located and reconstructed using the geometric volume of fluid (VOF) method. The implementation of the Piecewise-Linear Interface Calculation (PLIC) scheme of the VOF method was performed to solve the three-dimensional (3D) interface transport during the dynamics of two-phase flows. In this method, the interface is approximated by a line segment in each interfacial cell. The balance of forces at the interface was accounted for using the continuum interfacial force (CSF) model. To solve the Navier-Stokes equations, meshless finite difference schemes from the HiG-Flow computational fluid dynamics software were employed. The 3D PLIC-VOF HiG-Flow algorithm was used to simulate several benchmark two-phase flows for the purpose of validating the numerical implementation. First, the performance of the PLIC implementation was evaluated by conducting two standard advection numerical tests: the 3D shearing flow test and the 3D deforming field test. Good agreement is obtained for the 3D interface shape using both the 3D PLIC-VOF HiG-Flow algorithm and those found on the scientific literature. In addition, the absolute error of the volume tracking advection calculation obtained by our 3D PLIC-VOF HiG-Flow algorithm was found to be smaller than the one found in the scientific literature for both the 3D shearing and 3D deforming flow tests. Lastly, the 3D bubble rising problem was simulated for different fluid densities and viscosity’s ratios. Again, good agreement is obtained for the 3D interface shape using both the newly implemented algorithm and the OpenFOAM, DROPS and NaSt3D software. The bubble’s rise velocity and evolution of the bubble’s center of mass is also computed with the 3D PLIC-VOF HiG-Flow algorithm and found to be in agreement with these software.

Lamb waves occur in thin-walled structures in two wave modes, the symmetric and antisymmetric mode. Their oscillation on the structures‘ surfaces is either in phase (symmetric) or shifted by a phase angle of π (antisymmetric). In this work, a method is developed to compare the surfaces‘ oscillation phase relation. It is based on the evaluation of time signals regarding the instantaneous phase angle using the continuous wavelet transformation and as a comparative method the short-time Fourier transformation. For this purpose, numerical simulations utilizing the finite element method provide time signals from the top and bottom surface of different thin-walled structures. They differ with respect to their material settings and laminate configurations. The numerically obtained time signals are evaluated by the developed methods. The occurring oscillation phase differences on the top and bottom surface are studied and both methods are compared. Subsequently, the oscillation phase is evaluated experimentally for the wave propagation in a fiber metal laminate. It is shown that the method based on the continuous wavelet transformation is suitable for the evaluation of oscillation phase relations in time signals. Additionally, it is proven that fiber metal laminates show only two phase relations which indicates the occurrence of Lamb waves.

This work investigates how integrated polyimide inlays with applied sensor bodies influence the guided ultrasonic wave propagation in narrow glass fiber-reinforced polymer specimens. Preliminary numerical simulations indicate that in a damping-free specimen, the inlays show reflections for the S0-mode propagation. Hence, an air-coupled ultrasonic technique and a 3D laser Doppler vibrometer measurement are used to measure different parts of the propagating waves’ displacement field after burst excitation at different frequencies. No significant reflections on the inlay can be seen in the experiments. However, it is shown that the reflections from the strip specimen’s narrow width cause periodical reflections that superimpose with the excited wave fronts. A continuous wavelet transformation in the time-frequency domain filters discontinuities from the measurement signal and is used for reconstruction of the time signals. The reconstructed signals are used in a spatial continuous wavelet transformation to identify the occurring wavelengths and hence to prove the assumption of reflections from the narrow edges. Since the amplitude of the reflections identified in the numerical data at the polyimide inlays are an order of magnitude smaller than the excited wave packages, it is concluded that material damping of the epoxy resin matrix extincts possible reflections from the inlays.

: In this paper we present an optimization of the planar manufacturing scheme for stretch-free, shape induced metal interconnects to simplify the fabrication with the aim of maximizing the flexibility in the structure regarding to stress and strain. The Formation of trenches between silicon islands are actively used in the lithographic process to create arc shape structures by spin coating resist into the trenches. The resulting resist form is used as a template for the metal lines, which are structured on top. Because this arc shape is beneficial for the flexibility of these bridges. The trench depth as a key parameter for the stress distribution is investigated by applying numerical simulations. The simulated results show that the increase in penetration depth of the metal bridge into the trench increases the tensile load which is converted into a shear Force Q(x), that usually leads to increased strains the structure can generate. For the fabrication the filling of the trenches with resist is optimized by varying the spin speed. Compared to the theoretical resistance, the current-voltage measurements of the metal bridges show a similar behavior and almost every structural variation is capable of functioning as a flexible electrical interconnect in a complete island-bridge array.

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

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

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

We are witnessing a progressive divestment of some institutions with strong traditions and skills in physical modelling and their consequent impoverishment, to the detriment of numerical modelling. For many reasons, the economic imperatives and the exponential growth of computational means and numerical methods should certainly not be excluded. In this work, we aimed to highlight the new requirements of the recent sophisticated developments in physical modelling, precisely due to the new needs imposed on them by mathematical and numerical modelling and the growing risks in civil construction works. In this context, reflections are reported, justified by scientific and real-world examples, on the need for maintenance and reinforcement of investments in physical modelling, both to support the scientific community and to design buildings of significant economic, social and environmental impact.

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

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

This study develops a predictor-corrector algorithm for the numerical simulation of IVPs involving singular generalized fractional derivatives with Mittag-Leffler kernels. The proposed algorithm converts the considered IVP into a Volterra-type integral equation and then uses Trapezoidal rule to obtain approximate solutions. Numerical approximate solutions of some singular generalized fractional derivative with Mittag-Leffler kernels models have been presented to demonstrate the efficiency and accuracy of the proposed algorithm. The algorithm describes the influence of the fractional derivative parameters on the dynamics of the studied models. The suggested method is expected to be effectively employed in the field of simulating generalized fractional derivative models

Lithium-metal batteries offer high energy density, but are challenged by dendrite growth, which can lead to short circuits and battery failure. Various models with varying degrees of accuracy and computational cost have been developed to understand and predict dendrite growth. In this study, a simple simulation model to predict macroscale dendrite growth on lithium metal electrodes is presented. The model uses a 3D single-particle diffusion-limited aggregation (DLA) algorithm, biased by the electric field, to simulate dendrite growth. The model was validated using 3D MRI of dendrite growth in a lithium symmetric cell. It was used to simulate deposition between consecutive measured dendrite structures. The model also introduces an important parameter, namely the biasing factor c, which determines the balance between diffusion and electric field effects. This parameter was calibrated using the overlap and displacement of deposition sites metrics, achieving optimal results at c = 0.7. The model predictions were in good agreement with experimental data confirmed by fractal dimension analysis. This model provides a useful tool for predicting and studying dendrite formation in lithium and other metal batteries.

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

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.

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.

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

In the steeply inclined and extremely thick coal seams (SIETCS) mined using the horizontal sublevel top-coal caving (HSTCC) method, the uncertainty of gas emission is a safety threat to the mining operations. In order to reduce the occurrence of accidents, the determination of gas emission is crucial. In this paper, we first proposed a prediction model for workers at the floor coalbed to calculate gas emission on site. We then put forward a finite element numerical simulation for researchers to predict gas emission from the floor coalbed. At last, we measured gas emitted from the floor coalbed of SIETCS in Wudong Coal Mine in a specific mining period and used the data to verify the applicability and accuracy of these two gas-emission prediction methods. The results showed that the gas emission from Wudong Coal Mine was 1.08 m³/min calculated based on the prediction model and 1.07 m³/min obtained using the user-defined integration method. Both methods have their own advantages, disadvantages and applicable objects, and are important in predicting gas emission from SIETCS mined using HSTCC method.

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

Coastal wave modeling and forecasting are essential in oceanography, sustainable marine energy, and ocean engineering. Precise forecasting of wave speed and direction are crucial for offshore operations, marine energy, risk management, environmental management, coastal and sustainable maritime management. The coastline of Brunei Darussalam can generate between 15 and 126 Giga Watt of wave energy. In this experimental research, we used two numerical approaches, the finite difference, and spectral element methods, to model and simulate wave speed and direction. We calculated the mean error between numerical and analytical solutions. We proposed a novel, promising, univariate time series forecasting model by combining the Long Short-Term (LSTM) with KerasTuner hyperparameters tuning and optimization techniques. This method helps us to improve the accuracy and efficiency of time-series forecasting.. The experimental data was computed from high-precision Acoustic Doppler Current Profiler (ADCP) sensor data. This research is part of the preliminary feasibility analysis of wave energy production in Brunei Darussalam and net zero commitment for a sustainable environment. Seven independent forecast experiments were performed for wave speed and direction in degree and radian units for 1, 3, 6, 8, 10, 12, and 24 hours. Mean squared error (MSE) was adopted as a metric for both training and testing. The experimental results reveal that the wave speed forecast has the lowest MSE compared to direction, regardless of the unit of measure, but has a longer duration. In addition, the direction forecast in the degree unit has the lowest errors compared to the unit of radians; the latter has a longer running time than the former. The model has delivered optimal results throughout the experiments with minor training and test errors. We conducted a thorough evaluations on two benchmark time series datasets, which include the study dataset and air quality index dataset, to validate the performance of the proposed model with other models. The proposed model outperforms cutting-edge forecasting models, such as the conventional LSTM, ARIMA, and Prophet. The model has the slightest forecast error compared to the existing literature’s result.

Recently, Indium Tin Oxide, a highly transparent, well conductive, and CMOS-compatible material, has been paying strong attention to the thermo-optic controlled silicon photonics industry because it allows a miniature of the gap between the core silicon and the heater, thus enabling reducing the electric power consumption and enhancing the switching speed. In this article, we propose an ultralow loss and small-size ITO microheater for the phase shift tuning. The designated microheater is manipulated in realizing a numerical co-design of compact and high bandwidth three-mode converter and three-mode selective router. Simulation results illustrate the 3-dB bandwidth for the three-mode converter and three-mode selective router as much as 100-nm and 40-nm during crosstalk under -25 dB, respectively. Besides, co-designed devices attain relatively large fabrication tolerances corresponding to width and height tolerances of ±50 nm and ±5 nm. In addition, the proposed devices consumed less than 90 mW total power consumption and took a fast switching time below 8 μs. Moreover, both two co-designs can be integrated into an estimated compact footprint of 8 μm2160 μm. Such excellent performances demonstrate the attractive potential of ITO as low-loss thermo-optic phase shifters and open an alternative way for enabling ultrafast and high-speed mode division multiplexing systems and very large-scale photonic integrated circuits.

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