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

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

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

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

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

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

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

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

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.

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

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

Wildland fires are responsible for large socio-economic impacts. Fires affect the environment, damage structures, threaten lives, cause health issues, and involve large suppression costs. These impacts can be mitigated via accurate fire spread forecast to inform the incident management team. We show that a fire forecast system based on a numerical weather prediction (NWP) model coupled with a wildland fire behavior model can provide this forecast. This is illustrated with the Chimney Tops II wildland fire responsible for large socio-economic impacts. The system is run at high horizontal resolution (111 m) over the region affected by the fire to provide a fine representation of the terrain and fuel heterogeneities and explicitly resolve atmospheric turbulence. Our findings suggest that one can use the high spatial resolution winds, fire spread and smoke forecast to minimize the adverse impacts of wildland fires.

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

A three-dimensional, unstructured grid, hydrodynamic and suspended-sediment transport model (i.e., SELFE-SED) was developed to simulate temporal and spatial variations of suspended sediment and was applied to the subtropical subalpine Tsuei-Feng Lake (TFL) of Taiwan. The model was validated with measured water level and suspended‑sediment concentration in 2009, 2010, and 2011. The overall model simulation results are in quantitative agreement with the observational data. The validated model was then applied to explore the most important parameter that affects the suspended-sediment concentration and to investigate the effect of wind stress on the mean current and suspended‑sediment distribution in this shallow lake. Modeling results of sensitivity analysis reveal that the settling velocity is a crucial parameter and erosion rate is less important in the suspended-sediment transport model. Remarkable lake circulation was found based on the strength of wind speed and wind direction. Strong wind would result in higher mean current in the top layer and suspended-sediment distribution in the top and bottom layers. This study demonstrated that the wind stress played a significant influence on mean circulation and suspended-sediment transport in a shallow lake.

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

Based on the goals of "high reliability, high frequency, rapid launch, and low cost" for space launch site, and drawing inspiration from the design principles of SpaceX's Starship launch test facility, an integrated dual-sided deflector system for convective cooling and thermal protection is presented. The interaction process between the gas jet and liquid water jet and its effect on the flow field environment are thoroughly studied using numerical calculation methods. Furthermore, Regarding the phase-change heat transfer issue in a compressible gas-liquid two-phase flow and the varying distribution of different bubble shapes and sizes at the gas-liquid interface, a modified Lee model is derived. The rationality and accuracy of the phase-change mass-transfer coefficient values can be improved by adjusting the equivalent diameter and shape factor of the bubbles. The research results demonstrate that compared to the classical Lee model, the modified Lee model can achieve a higher numerical accuracy in predicting the heat and mass transfer processes in gas-liquid two-phase flows. Through comparative analysis with traditional dual-sided deflector and conventional cooling system, the integrated dual-sided deflector system for convective cooling and thermal protection exhibits significant performance advantages in gas flow regulation, flow field environment improvement and erosion protection in the near-ground region of the launch site. Particularly, regarding the thermal environment of the dual-sided deflector, it not only achieves effective flow deflection, but also mitigates the degree of erosion caused by the gas jet on the deflector. This conclusion can provide theoretical references for the thermal protection design of commercial launch vehicle system at space launch site.

The ground consists of many layers of soil with different properties. The propagation speed and path of seismic waves are affected by different soil layers. It is necessary to understand that layered soil exhibits different dynamic behaviors and responses under the action of seismic waves. This study utilized weathered soil and silica sand as materials to create multi-layered soil conditions with varying degrees of compaction. By conducting 1 g shaking table test on multi-layered soil, the interaction and influence between different soil layers under different earthquake conditions were observed. The approach of numerical analysis aims to validate the experimental results and provide an in-depth understanding of the dynamic behaviors of multi-layered soil grounds during seismic events. The results of this study show that the amplification effect of seismic waves is related to the number of soil layers and the degree of compaction of the soil layers. This indicates that the multi-layered soil ground and the behavior of the soil layers play an important role in the propagation and impact of seismic waves, and this amplification effect is of great significance in actual seismic disaster risk assessment and design.

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

This paper introduces an innovative approach to numerically model Structure-Soil-Structure Interaction (SSSI) by integrating the Boundary Element Method (BEM) and the Finite Element Method (FEM) in a coupled manner. To assess the accuracy of the proposed method, a comparative study is undertaken, comparing its outcomes with those generated by the conventional FEM technique. Alongside accuracy, the computational efficiency aspect is crucial for the analysis of large-scale SSSI problems. Hence, the computational performance of the coupled BEM-FEM method undergoes a thorough examination and is compared with that of the standalone FEM method. The results from these comparisons illustrate the superior capabilities of the proposed method in comparison to the FEM method. The novel approach provides more reliable results compared to traditional FEM methods, serving as a valuable tool for engineers and researchers involved in structural analysis and design.

To identify the exchange coefficients of a nonlinear polynomial tricompartmental general system of α+β , we follow the procedure; firstly, the recommended solution is to introduce an adequate time t*>0 in a determinable manner. That is, after injecting quantity a into the main compartment, wait a moment for the exchange in the polynomial α+β order nonlinear general system to settle, and then compare this compartment to compartment 2 at t*. Secondly, applying the Taylor formula will linearize the system and identify the exchange coefficients. Finally, we will prove that the linearization method is stable.

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.

Employing numerical methods to simulate manufacturing processes and study their influential factors is increasing. The development of computers and engineering software has increased the ability and accuracy of numerical methods in simulating various aspects of different processes. One of the manufacturing processes that has been considered by many industries in recent years to join metals in solid-state with unique properties is called friction stir welding. It is challenging and, in some cases, impossible to experimentally study the various aspects of this process, such as temperature distribution, stress distribution, and material flow, due to severe plastic deformation in the weld zone. In this case, numerical methods are used to investigate these parameters and better understand the process. This study first investigates various numerical methods researchers use to simulate the process. The ability, pros, and cons of these methods to simulate the process are also examined. Then, the applications of numerical models in simulating the temperature distribution during the process as well as the effects of input parameters on the temperature history of the process, are deeply considered. Next, the application of numerical methods in material flow modeling during the process is investigated. Finally, the modeling of the microstructure of the welding zone is investigated using numerical methods that help significantly in predicting the weld microstructure.

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.

In this paper, we consider the fourth-order parabolic equations with integer and fractional order time derivatives with Neumann boundary conditions. The integer order time derivatives are approximated by backward Euler difference quotients, and the fractional order time derivatives are approximated by L1 interpolation. We propose the block-centered finite difference scheme for fourth-order parabolic equations of integer and fractional order time derivatives. We prove the stability of the block-centered finite difference scheme and the second-order convergence of the discrete L2 norms of the approximate solution and its derivatives of every order. Numerical examples are given to verify the effectiveness of the block-centered finite difference scheme.

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.

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.

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

This paper focuses on the combination of the method of moments and the wire-grid approximation as an effective computational technique for modeling symmetrical antennas with low computational cost and quite accurate results. The criteria and conditions for the use of wire-grid surface approximation from various sources are presented together with new recommendations for modeling symmetrical antenna structures using the wire-grid approximation. These recommendations are used to calculate the characteristics of biconical and horn antennas at different frequencies. The results obtained using different grid and mesh settings are compared to those obtained analytically. Moreover, the results are compared to those obtained using the finite difference time domain numerical method, as well as the measured ones. All results are shown to be in a good agreement. The used recommendations for building a symmetrical wire-grid of those symmetrical antenna elements provided the most advantageous parameters of the grid and mesh settings and the wire radius, which are able to give a quite accurate results with low computational cost. Additionally, the known equal area rule was modified for a rectangular grid form. The obtained radiation patterns of a conductive plate using both the original rule and the modified one are compared with the electrodynamic analysis results. It is shown that the use of the modified rule is more accurate when using a rectangle grid form.

The process of revitalising quiescent cells in order for them to proliferate plays a pivotal role in the repair of warn out tissues as well as for tissue homeostasis. This process is also crucial in the growth, development and well-being of higher multi-cellular organisms such as mammals. Deregulation of quiescent-proliferation transition is related to many diseases, in particular cancer. Recent studies have revealed that this process is regulated tightly by the Rb−E2F bistable switch mechanism. Based on experimental observations, in this study, we formulate a mathematical model to examine the effect of the growth factor concentration on the proliferation-quiescence transition in human cells. Working with a non-dimensionalised model, we prove positivity, boundedness and uniqueness of solutions. To understand model solution behaviour close to bifurcation points, we carry out bifurcation analysis, which is further verified by use of numerical bifurcation analysis, sensitivity analysis and numerical simulations. Indeed, bifurcation and numerical analysis of the model predicted a transition between stable, bistable and stable states, which are dependent on the growth factor concentration parameter (GF). The derived predictions confirm experimental observations.

Many authors proposed analytical methods for solving fully fuzzy Sylvester matrix equation (FFSME) based on Vec-operator and Kronecker product. However, these methods are restricted to nonnegative fuzzy numbers and cannot be extended to FFSME with near-zero fuzzy numbers. The main intention of this paper is to develop a new numerical method for solving FFSME with near-zero trapezoidal fuzzy numbers that provides a wider scope of trapezoidal fully fuzzy Sylvester matrix equation (TrFFSME) in scientific applications. This numerical method can solve the trapezoidal fully fuzzy Sylvester matrix equation with arbitrary coefficients and find all possible finite arbitrary solutions for the system. In order to obtain all possible fuzzy solutions, the TrFFSME is transferred to a system of non-linear equations based on newly developed arithmetic fuzzy multiplication between trapezoidal fuzzy numbers. The fuzzy solutions to the TrFFSME are obtained by developing a new two-stage algorithm. To illustrate the proposed method numerical example is solved.

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

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.

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.

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

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.

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.

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

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

This study focuses on the formulation and analysis of a numerical method to compute singularly perturbed functional differential equations (SP-FDEs) involving large mixed shifts. SP-FDEs are a class of differential equations with both shifts and small perturbation parameter. Such equations arise in various scientific fields to model phenomena with memory effects and perturbed dynamics. In this particular study, the considered SP-FDEs involve large negative shift (delay) and large positive shift (advance). Dealing with large mixed shifts in SP-FDEs gives rise to challenges in terms of stability analysis, numerical methods, and convergence of a solution. In this paper, the influence of the large mixed shifts is treated by choosing a uniform mesh discretized in such a way that the term with the shifts lie on the nodal points. Then, using the Numerov’s finite difference method, the continuous problem is transformed into a finite difference scheme and the influence of the perturbation parameter is handled by introducing a fitting factor in the difference scheme. Stability estimate and convergence analysis of the developed scheme are investigated and proved. Numerical experiments are carried out to demonstrate the validity and applicability of the developed scheme. It is shown that the numerical scheme obtained in this work is parameter-robust of second order convergence rate.

This paper presents an original method of additive manufacturing of cylindrical parts with variable circumference thickness, which allows the control of the deposition of molten material using an algorithm for decomposing the part geometry into volumetric elements with known dimensional configuration. In the absence of a post-processor capable of controlling additive manufacturing on a 5-axis numerical control machine, control of the deposition of molten material is done using parameterized programs, which can control both the feed speeds of the machine tool axes and the specific functions of the printing equipment. Additive manufacturing can make a positive contribution to sustainable development compared to traditional manufacturing technologies, thus making a positive contribution for a sustainable future. The aim of the work is to make it possible to 3D print parts with variable wall thickness using a CNC machining centre. To obtain the variable thickness layer we have implemented an original method of deposition with molten material (FDM), the coordination of the system composed of physical elements respectively programmable elements, is realized through control functions materialized in parameterized part programs, the generated outputs being the variable speed of the machine axes on a circular trajectory, the angular positioning, the filament advance.

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.

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.

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

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

Time evolution operators of a strongly ionizing medium are calculated by a time-dependent unitary transformation (TDUT) method. The TDUT method has been employed in quantum mechanical system composed of discrete states. This method is especially helpful for solving molecular rotational dynamics in quasi-adiabatic regimes because the strict unitary nature of the propagation operator allows us to set the temporal step size large; a tight limitation on the temporal step size ($\delta t <<1$) can be circumvented by the strict unitary nature. On the other hand, in a strongly ionizing system where the Hamiltonian is not Hermitian, the same approach cannot be directly applied because it is demanding to define a set of field-dressed eigenstates. In this study, the TDUT method was applied to the ionizing regime using the Kramers-Henneberger frame, in which the strong-field-dressed discrete eigenstates are given by the field-free discrete eigenstates in a moving frame. Although the present work verifies the method for a one-dimensional atom as a prototype, the method can be applied to three-dimensional atoms, and molecules exposed to strong laser fields.

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.

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.

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.

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.

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.

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.

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

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

In this paper, low circumferential reciprocating load foot-scale tests were performed on two truncated PHC B 600 130 tubular piles with bearing nodes to characterize the damage process and morphology of the specimens and to investigate the load-carrying performance of the members. The test results reveal that under the action of tensile-bending-shear loading, the bearing concrete in the node area buckles and is damaged, the anchored reinforcement in the node area yields, the constraint is weakened, an articulation point is formed, and the node rotational capacity increases. When the embedment depth increases from 200 mm to 300 mm, the ultimate bearing capacities of the positive and negative nodes increase by 57.60% and 54.60%, respectively. A numerical simulation is used to verify the test results. Considering the three types of piles with truncated nodes, the numerical simulation is used to analyze the node-bearing capacity at different embedment depths. Finally, two preferred node types are proposed as follows: pipe pile core-filled longitudinal reinforcement anchored into the bearing node and pipe pile body longitudinal reinforcement anchored into the bearing node + pipe pile core-filled longitudinal reinforcement anchored into the bearing node, with preferred embedment depths of 350 mm and 200 mm, respectively.

This paper proposes a lightweight joint with bamboo culm cross connection, which is constructed by screwing one side of the screw into the other side of the wooden plug, its simple construction makes it easy to manufacture and cost-effective. This form connects adjacent bamboo culms to form corner joints, and bolt connects adjacent bamboo culm frames to form edge joints, to explore the mechanical properties of the combination joints formed by corner joints and edge joints, by combining experiments and numerical simulations, the method discovered and explained the relationship and mechanism between the number and position of bolts and the stiffness of edge joints. A calculation model for frame-unit bamboo culm members considering semi-stiffness of joints is proposed, and the conversion equations between linear stiffness of corner and edge joint and their rotational stiffnesses are derived. By establishing the series and parallel relationship between the corner and edge joints in the assembled members, a generalized relational equation between the total stiffness of the assembled joint and the stiffness of the corner and edge joints is derived. A numerical simulation method for frame-unit bamboo culm members with semi-rigid joints is proposed, which can obtain the stiffness calculation formula of edge joints and the average rotational stiffness value of semi-rigid corner joints. A finite element analysis method based on semi-rigid joints is proposed for frame-unit bamboo culm members, and the internal force sketches and calculation methods for bamboo culms and joints are presented, which can be used for the internal force analysis and design of frame-unit bamboo culm lattice structures.

When the formation flying consists of satellites that are independent from each other and forms a virtual synthetic aperture telescope that replaces a single large telescope with all of its optical elements it is required mutual coordination between the satellites with high accuracy in order to achieve good optical observations. In this paper we consider the problem of developing a control system for a tetrahedral satellite formation flying representing a synthetic aperture telescope in the geostationary orbit for high resolution monitoring of fire in the forest area in infrared spectrum. Synthesizing the image with good quality requires to keep the configuration of formation with µm-class accuracy. Study of dynamics of passive uncontrolled motion of formation flying showed a significant deviation of the configuration from the required shape with a high frequency due to the action of gravitational forces. To keep the configuration in the required form with µm accuracy the analysis of efficiency of various controllers was carried out in the process of numerical simulation. The simulation results made it possible to highlight the features of using various approaches to the development of control system for the satellite formation flying representing a synthetic aperture telescope.

Bistatic scattering coefficient have important application of ocean active microwave remote sensing. In this paper, a general expression of the composite scattering coefficient for multiple dielectric targets of one-dimensional fractal sea surface is derived by MOM method. Secondly, the bistatic scattering coefficient is computed, there are two class missile targets and one class ship target of one-dimensional fractal sea surface, and the influence of the water dielectric constant, fractal dimension, wind speed, the target dielectric constant and spatial position is discussed, and the corresponding conclusions are given. The proposed algorithm is able to solve some composite scattering problems, such as “dielectric or conductive object + rough surface”, “dielectric or conductive floating object +rough ”and it has certain generalization and reference.

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.

This paper introduces a method to achieve optimal Gaussian estimation even when operating with low precision numerical representations of the full covariance matrix, such as 16-bit integer or 32-bit single precision formats, instead of more expensive double precision floating point representations (64 bits). The approximation's numerical error is accurately estimated in a conservative manner, resulting in minimal covariance inflation. When combined with a compressed estimator like the GCKF (Global Compressed Kalman Filter), which performs global updates at low frequency, this approach produces estimates that are nearly identical to the optimal ones. The deviation between the results obtained using this new approach and those obtained with the standard GCKF (which operates in double precision) is negligible. This new method is particularly relevant in high-dimensional estimation problems that require high-frequency operation. By leveraging the compressed operation of the GCKF, the estimation process applies the necessary approximation only during low-frequency global updates, using the single precision format instead of the more expensive double precision representations. The introduced numerical approximation is accurately estimated while maintaining a conservative approach, resulting in minimal covariance inflation. When combined with the GCKF, which performs global updates at low frequency, the estimates produced are almost identical to the optimal ones. This approach is particularly useful in high-dimensional estimation problems that require high-frequency operation. By representing the full covariance matrix in a low precision numerical format and applying appropriate scaling, memory usage is significantly reduced, offering substantial benefits for high-dimensional estimation problems and enabling real-time smoothing processes in such cases. By reducing the memory requirements for storing the full covariance matrix, the CPU avoids encountering a high number of page faults, resulting in improved processing performance. Lastly, the experimental section verifies the performance of the proposed approach in a single precision (32-bit floating point) GCKF SLAM process, assessing expected values and standard deviations of marginal PDFs.

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.

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.

This work is concerned with the existence and uniqueness of boundary value problems defined on semi-infinite intervals. These kinds of problems seldom admit exactly known solutions and, therefore, the theoretical information on their well-posedness is essential before attempting to derive an approximate solution by analytical or numerical means. Our utmost contribution in this context is the definition of a numerical test for investigating the existence and uniqueness of solutions of boundary problems defined on semi-infinite intervals. The main result is given by a theorem relating the existence and uniqueness question to the number of real zeros of a function implicitly defined within the formulation of the iterative transformation method. As a consequence, we can investigate the existence and uniqueness of solutions by studying the behaviour of that function. Within such a context the numerical test is illustrated by two examples where we find meaningful numerical results.

The work presents an approach to instrument the load sensing bearings for automotive applications for estimation of the loads acting on the wheels. The system comprises fiber-optic sensors based on addressed fiber Bragg structures (AFBS) with two symmetrical phase shifts. A mathematical model for load-deformation relation is presented, and the AFBS interrogation principle is described. The simulation includes (i) modeling of vehicle dynamics in a split-mu braking test, during which the longitudinal wheel loads are obtained, (ii) the subsequent estimation of bearing outer ring deformation using a beam model with simply supported boundary conditions, (iii) the conversion of strain into central wavelength shift of AFBS, and (iv) modeling of the beating signal at the photodetector. The simulation results show that the estimation error of the longitudinal wheel force from the strain data acquired from a single measurement point was 5.44% with root-mean-square error of 113.64 N. A prototype load sensing bearing was instrumented with a single AFBS sensor and mounted in a front right wheel hub of an experimental vehicle. The experimental setup demonstrated comparable results with the simulation during the braking test. The proposed system with load-sensing bearings is aimed at estimation of the loads acting on the wheels, which serve as input parameters for active safety systems, such as automatic braking, adaptive cruise control, or fully automated driving, in order to enhance their effectiveness and safety of the vehicle.

The boundary layer height (BLH) distinguishes the interface between the lower and free atmosphere, which is a key variable in numerical simulation, aerosol, and environmental pollution studies. Thus, the goal of this work is to propose a novel method conjuncting with numerical regularization to analyze climate characteristics of the marine boundary layer height (MBLH) using 2007-2011 GPS-RO data from the COSMIC mission. While traditionally, the difference method has been used to achieve this aim, herein, we propose an innovative method, in which the bending angle profile gradient was calculated using the numerical regulation method where the regulation parameters are determined by the double-parameter model function method. Then, the MBLH was determined by employing the maximum gradient method to ascertain the height corresponding to the smallest gradient. The results show a correlation between currents and the MBLH—a relationship that has not been previously demonstrated. A low MBLH is associated with seasons and regions where cold ocean currents are prevailing; whereas a high MBLH is observed in the seasons and places where warm currents are prevailing. This correlation was validated by comparing the obtained results with different occultation data—i.e., atmprf and echprf—which also established that atmprf is more sensitive to convective cloud top capture. In seas with active convection, the armprf calculated bending angle is higher than that from echprf. Subsequently, the standard deviation was used to express the MBLH confidence level. The results show that the MBLH standard deviation is highest in low latitudes and lowest in the middle and high latitudes. Furthermore, we analyzed the interannual MBLH variation trend, which displayed a seasonal variation and spatial distribution corresponding with the current and subsolar point. Finally, we conducted a case study in the South China Sea, and identified a distinctive seasonal change and downward trend.

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.

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.

A cosmological model was developed using the equation of state of photon gas, as well as cosmic time. The primary objective of this model is to see if determining the observed rotation speed of galactic matter is possible, without using dark matter (halo) as a parameter. To do so, a numerical application of the evolution of variables in accordance with cosmic time and a new state equation was developed to determine precise, realistic values for a number of cosmological parameters, such as energy of the universe , cosmological constant Λ, curvature of space k, energy density , age of the universe etc. Some assumptions were put forth in order to solve these equations. The current version of the model partially explains several of the observed phenomena that raise questions. Numerical application of the model has yielded the following results, among others: Initially, during the Planck era, at the very beginning of Planck time, , the universe contained a single photon at Planck temperature Planck energy in the Planck volume. During the photon inflation phase (before characteristic time ~10^-9 ), the number of original photons (alphatons) increased at each unit of Planck time and geometrical progression ~n³, where n is the quotient of cosmic time over Planck time . Then, the primordial number of photons reached a maximum of N~10⁸⁹, where it remained constant. It's primordial photons (alphatons) are still present today and represents the essential of the energy contained in the universe via the cosmological constante . Such geometric growth in the number of photons can bring a solution to the horizon problem through exchange and a photon energy volume that is in phase with that of the the volume energy of the universe. The age of the universe in cosmic time that is in line with positive energy conservation (in terms of conventional thermodynamics) and the creation of proton, neutron, electron, and neutrino masses, is ~76 [Gy] (observed . The predicted total mass (p, n, e and n), based on the Maxwell-Juttner relativistic statistical distribution, is ~7x10⁵⁰ . The predicted cosmic neutrino mass is ≤ 8.69x10^-32 (≤ 48.7 ) if based on observations of SN1987A. The temperature variation of the cosmic microwave background (CMB), as measured by Planck, can be said to be partially due to energy variations in the universe (DU/U) during the primordial baryon synthesis (energy jump from the creation of protons and neutrons). In this model, what is usually referred to as dark energy actually corresponds to the energy of the universe that has not been converted to mass, and which acts on the mass created by the energy-mass equivalence principle and the cosmological gravity field, , associated with the cosmological constant, which is high during primordial formation of the galaxies (<1 [Gy]). A look at the Casimir effect makes it possible to estimate a minimum Casimir pressure and thus determine our possible relative position in the universe at cosmic time 0,1813 (t₀/t_W=13,8 /76,1 ). Therefore, from the observed age of 13,8 [Gy], we can derive a possible cosmic age of 76,1 [Gy]. That energy of the universe, when taken into consideration during the formation of the first galaxies (< 1 [Gy]), provides a relatively adequate explanation of the non-Keplerian rotation of galactic masses. Indeed, such residual, non-baryonic energy, when considered in Newton’s gravity equation, adds the term , which can partially explain, without recourse to dark matter, the rotations of some galaxies, such as M33, UGC12591, UGC2885, NGC3198, NGC253, DDO161, UDG44, the MW and the Coma cluster. Today, in the MW, that cosmological gravity force is in the order of 10²⁶ times smaller than the conventional gravity force. The model predicts an acceleration of the mass in the universe (q~ 0,986); the energy associated with curvature is the driving force behind the expansion of the universe, rather than the energy associated with the cosmological constant . An equation to determine expansion is obtained using the energy form of the Friedmann equation relative to Planck power and cosmic time or Planck force acting at the frontier of the universe moving at c. This constant Planck force, from unknown sources, acts everywhere to the expansion of the universe as a stretching effect on the volume. Finally, the model partly explains the value a₀ of the MOND theory. Indeed, a₀ is not a true constant, but depends on the cosmological constant at the time the great structures were formed (~1 [Gy]), as well as an adjustment of the typical mass and dimension of those great structures, such as galaxies. The constant a₀ is a different expression of the cosmological gravity force as expressed by the cosmological constant, Λ, acting through the energy-mass equivalent during the formation of the structures. It does not put in question the value of G.

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.

In the process of use and manufacture, carbon fiber foam sandwich structures were often damaged by low-energy impact, resulting in performance degradation. Therefore, it was necessary to study the damage caused by low-speed impact of composite sandwich structures. In this paper, based on Hashin failure criterion, an equivalent finite element model of carbon fiber foam sandwich panel under low velocity impact was established. The model was used to simulate the damage of foam sandwich panel with [±45°/±45°/ (core) /±45°/±45°] ply structure under the impact energy of 10.58J, 21.17J, 31.75J and 42.34J. The simulation results of impact damage depth were compared with the experimental results. The error was less than 10%, which proved the rationality of the impact equivalent model. The model was used to predict and analyze the damage of foam sandwich panels with [±45°/ (core) /±45°], [±45°/ (0°, 90°) / (core) /±45°] and [±45°/ (0°, 90°) (core) / (0°, 90°) /±45°] ply structures under 21.17J impact energy. By comparing and analyzing the damage situation, impact force response time and impact velocity response time, the low energy impact resistance was analyzed. The results showed that increasing the number of ply structure of [±45°] can reduce the impact damage degree and improve the bearing capacity of sandwich panels.

This paper introduces a novel multi-strategy enhanced dung beetle optimization (MSDBO) algorithm that is designed to address several issues identified in the standard dung beetle optimization algorithm. Specifically, the MSDBO aims to enhance convergence speed, reduce susceptibility to local optima, and increase search accuracy. By incorporating three strategies: tent chaotic mapping for population initialization, the golden sinusoidal strategy for position updating, and the Lévy flight strategy for balancing exploration and exploitation, the standard dung beetle optimization algorithm is enhanced. The MSDBO algorithm is evaluated using twelve benchmark test functions and compared against five state-of-the-art algorithms. The results consistently show that MSDBO exhibits faster convergence speeds and more accurate solutions than the other algorithms across most of the test functions. In addition, MSDBO is also applied to optimize the parameters of a valve plate, including the close angle, cross angle, triangle groove sizes, and wrap angle. The optimization outcomes reveal that MSDBO effectively minimizes pressure ripples in the piston chamber, resulting in reduced flow rate fluctuations and noise emission compared to the initial design. This study highlights the potential of the MSDBO algorithm in tackling complex nonlinear engineering optimization problems.

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.

We study the tight focusing properties of generic bright and dark HyOPS beams in the presence of primary coma. The role of the polarization singularity index and handedness of the polarization of the HyOPS beams on the focused structure has been discussed. Results have been presented for the total intensity, component intensities, and component phase distributions for left-, and right-handed bright and dark star and lemon types singularities. The presence of primary coma distorted the focal plane intensity distributions for both positive and negative index generic C-points. Coma is known to disturb the circular symmetry of the focal plane intensity distribution. Similarly in tight focusing polarization is known to disturb the symmetry. Therefore, a beam with structured and inhomogeneous polarization distribution tightly focused under the influence of coma is a fit case to study. It is found that the presence of primary coma aberration in the focusing system produces a positional shift of the high-intensity peaks and a reduction of the intensity on one side of the center. As the strength of the primary coma increases, the focal plane intensity distributions shift more and more toward the right from the initial position. Unlike the scalar vortex case, in the case of HyOPS beams, the focal plane intensity distribution undergoes rotation, as the helicity of the HyOPS beams is inverted, in addition to shifting. All the component phase distributions are found to be embedded with phase vortices of charge ±1.

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.

In the southeastern coast of Brazil, the bays of Ilha Grande and Sepetiba are linked by the Ilha Grande Channel, where remarkably vigorous currents have been consistently observed. The sheer intensity of these currents defies a simple explanation through tidal forcing alone, prompting researchers to delve into potential influences such as basin resonance and the baroclinic effect. This study aims to elucidate the role of remote meteorological effects within this complex hydrodynamic system. A numerical approach with a coastal model nested within an ocean model was employed, enabling an in-depth examination of the intricate interplay between meteorological and tidal forcings. The study revealed a significant finding: the lag in signal propagation plays a pivotal role in determining how these signals impact the dynamics of the bays. The astronomical signal exhibits a minimal lag along the coast (1 minute), and leads to water level differences between the sea and the coastline, resulting in the generation of tidal currents at the bay entrances. On the other hand, the remote meteorological signal, with a stronger signal lag along the coast (4.92 hours), leads to the creation of water level difference between the bay’s entrances, inducing significant fluxes along the narrow Ilha Grande Channel.

This paper is concerned with developing an efficient numerical algorithm for fast implementation of the sparse grid method for computing the $d$-dimensional integral of a given function. The new algorithm, called the MDI-SG ({\em multilevel dimension iteration sparse grid}) method, implements the sparse grid method based on a dimension iteration/reduction procedure, it does not need to store the integration points, nor does it compute the function values independently at each integration point, instead, it reuses the computation for function evaluations as much as possible by performing the function evaluations at all integration points in a cluster and iteratively along coordinate directions. It is shown numerically that the computational complexity (in terms of CPU time) of the proposed MDI-SG method is of polynomial order $O(Nd^3 )$ or better, compared to the exponential order $O(N(\log N)^{d-1})$ for the standard sparse grid method, where $N$ denotes the maximum number of integration points in each coordinate direction. As a result, the proposed MDI-SG method effectively circumvents the curse of dimensionality suffered by the standard sparse grid method for high-dimensional numerical integration.

In crank presses, an important role is played by the presence of actuators with dwell, in which, with continuous movement of the input links, the output working link remains stationary at some part of the work. The use of a high-class structural group (fourth class) expands the functionality of the crank press, including it can ensure the link dwell of the working link on a given section of movement during the required period of time. As a result of the conducted research, four variants of actuators of the crank knee press were designed. Due to the rational choice of the positions of racks and crosshead guide, the best layout of the dimensions of the press mechanism has been achieved, changing the coordinates of the racks accordingly allows to design a mechanism with upper or lower dwells of the working links. When synthesizing (selecting) parameters such as: l_2 and h, it is possible to obtain a height of different duration. All this is shown on the basis of a numerical experiment based on kinematic analysis. In addition, the complex structure of the mechanism makes it possible to ensure the exact link dwell, which cannot be achieved due to the well-known knee presses based on the second-class mechanism [21].

Strong convective storms can be a serious threat to modern society, influencing both the economy and human life. Thunderstorms can bring much-needed rain after a dry spell, but the nature of thunderstorm mostly causes damage due to high wind associated with it and heavy intensity rain that mostly causes quick overlandflow and flashfloods. The paper focused on the 2011 Pukkelpop event in Belgium where due to heavy thunderstorm a lot of financial damage occurred, and five casualties. In this paper we look into the basic elements that enforces the creation of a thunderstorm like CAPE and windshear. A representative grid has been taken from the Alaro-0 model of 12.5 km resolution for the analysis of the thunderstorm over Belgium which is detailed out in the methodology. The model runs of Alaro-0 from Royal Meteorological Institute of Belgium (RMI) is used for historical and future years for CAPE and Shear. We look into STEnv which is a composite number that can provide some insight on occurrence of thunderstorm, and STEnv is when CAPE*Shear exceeds certain threshold. In this study the threshold has been chosen as 97%. Only summer months have been taken into account into this study as summer(s) has the most thunderstorm in Belgium. The analysis of CAPE and STEnv over historical period (1977-2005), near future (2041-2069) and far future (2071-2099) shows significant increase in occurrence of thunderstorm in the 2071-2099s, specially in extreme cases such as 97% and 99% quantile. CAPE (J/Kg) increases in all future scenarios compared to the historical period. In 2071-2099s, CAPE increases from 2100 J/Kg to 2300 J/Kg (90% quantile) and from 3500 J/Kg to 3900 J/Kg (99% quantile) compared to the historical period. The number of STEnv increases from 300 to 378 in top 10% storm case (90% quantile), 158 to 218 in top 5% storm (95% quantile), 97 to 153 in top 3% storm (97% quantile) and 33 to 66 in top 1% storm (99% quantile) with a 90% confidence interval overall. Shear mostly remains unchanged. The study also shows spatial variation of CAPE, shear and STEnv over Belgium for present and future scenarios in different extremities (90%, 95%, 97% and 99% quantiles). There is a good indication that increase in CAPE translates to increase in STEnv, which is an accepted measure to predict occurrence of thunderstorm. Future climate change also amplifies the occurrence of thunderstorm, which is coherent with relevant studies. The uncertainty related to the study can be reduced with more detailed dataset and involvement of numerical weather prediction models beyond data statistics.

Promoting the use of solar energy resources has always faced the challenges of instability and supply-demand mismatch. The key to solving these issues is to efficiently store and utilize solar energy resources using high-performance heat storage devices. This study designed a high-performance shell-and-tube phase change thermal storage device and established a numerical model using ANSYS software to summarize the device's dynamic melting law. To verify the accuracy of the numerical simulation, a performance testing platform for the phase change thermal storage device was built to investigate the impact of factors such as inlet water temperature, inlet water flow rate, type of heat storage, and initial temperature of the device, and to reveal the change law of the device's performance. The results show that the inlet water temperature has the most significant impact on the device's heat storage and release performance. When the device's heat storage or release is used for heating, changing the inlet water flow rate has a weak and limited effect on the device's performance. However, when the device's heat release is used to provide domestic hot water, increasing the makeup water temperature and reducing the inlet water flow rate can significantly improve the device's effective heat release. Furthermore, based on the experimental validation of the model's correctness, this study further simulated and studied the impact of different factors on the device's heat storage process to optimize its structural design and provide technical references for the device's actual operation and installation. The results show that the placement of fins has a negligible effect on the performance of the heat storage device while reducing the fin spacing and increasing the fin thickness can significantly improve the melting efficiency of the phase change material (PCM). Additionally, the heat storage characteristics of the device are significantly better in the vertical installation mode than in the horizontal installation mode. This study provides theoretical guidance and technical references for the design and use of phase change thermal storage devices.

(1) The article is devoted to the development of a theoretical and algorithmic basis for numerical modeling of the spontaneous potential method (SP) as applied to the study of sandy-argillaceous reservoirs. (2) In terms of coupled flows, we consider a physical-mathematical model of SP signals from an electrochemical source, with regard to the case of fluid-saturated shaly sandstone. (3) An algorithm for 2D finite-element modeling of SP signals was developed and implemented in software, along with its internal and external testing with analytical solutions. The numerical SP modeling was carried out, with determining the dependences on the reservoir thickness and porosity, the amount of argillaceous material and the type of minerals. We performed a comparative analysis of the simulated and field SP data, using the results of laboratory core examinations taken from wells in a number of fields in the Latitudinal Ob Region of Western Siberia. (4) The results of the study may be used either for the development of the existing SP techniques, by providing them with a consistent computational model, or for the design of new experimental approaches.

This study identifies new Alkali Activated Composite Mud Mortar to alternate Cement Mortar used for laying masonry works. This study handles both experimental and numerical analysis of strength and durability for mortar, especially bond strength along mortar and brick joints have studied detail and compared with ANSYS. The bond strength has estimated through the triplet brick prism of shear behavior. The strength and deformation characteristics of masonry constituents obtained from these tests depend on the actual composite conduct of masonry. There is always a good reception for environmentally friendly, low energy, new construction materials, and used in preparing mortar/brick unit combinations all around the world. Also, this study fulfills those requirements. Especially, masonry laid with the mud mortar meets all. Fly Ash, Ground Granulated Blast furnace Slag, and Quarry Dust have utilized the new composites' preparation with 8Molority, 10Molority developed excellent strength. The composite specimens are having compression and bond strength, ranging from 4.84-5.14 and 0.12-0.18N/mm2. The mixed ratio of materials has taken (0.5 parts of Fly Ash, 0.5 parts of Steel Slag, 1.75 parts of Soil, and 0.25 parts of Quarry Dust). The samples tested resistance against acidic, sulfate, and alkalinity indicated PH values less than 7 and 9 simultaneously. Compression strength varied from 10-20% and gained more and fit composite used for mortar in masonry works.

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.

5G heterogeneous network overlaid by millimeter-wave (mmWave) access employs mmWave meshed backhauling as a promising cost-efficient backhaul architecture. Due to the nature of mobile traffic distribution in practice which is both time-variant and spatially non-uniform, dynamic construction of mmWave meshed backhaul is prerequisite to support the varying traffic distribution. Focusing on such scenario of outdoor dynamic crowd (ODC), this paper proposes a novel method to control mmWave meshed backhaul for efficient operation of mmWave overlay 5G HetNet through Software-Defined Network (SDN) technology. Our algorithm is featured by two functionalities, i.e., backhauling route multiplexing for overloaded mmWave small cell base stations (SC-BSs) and mmWave SC-BSs’ ON/OFF status switching for underloaded spot. In this paper, the effectiveness of the proposed meshed network is confirmed by both numerical analyses and experimental results. Simulations are conducted over a practical user distribution modeled from measured data in realistic environments. Numerical results show that the proposed algorithm can cope with the locally intensive traffic and reduce energy consumption. Furthermore, a WiGig (Wireless Gigabit Alliance certified) device based testbed is developed for Proof-of-Concept (PoC) and preliminary measurement results confirm the proposed dynamic formation of the meshed network’s efficiency.

The development of computational acoustics allows simulation of sound generation and propagation in complex environment. In particular, meshfree methods are widely used to solve acoustics problems through arbitrarily distributed field points and approximation smoothness flexibility. As a Lagrangian meshfree method, smoothed particle hydrodynamics (SPH) method reduce the difficulty in solving problems with deformable boundaries, complex topologies, or multiphase medium. The traditional SPH method has been applied in acoustic simulation. This study presents the corrective smoothed particle method (CSPM), which is a combination of SPH kernel estimate and Taylor series expansion. The CSPM is introduced as a Lagrangian approach to improve accuracy in solving acoustic wave equations in the time domain. Moreover, a boundary treatment technique based on the hybrid meshfree and finite difference time domain (FDTD) method is proposed to represent different acoustic boundaries with particles. To model sound propagation in pipes with different boundaries, soft, rigid, and absorbing boundary conditions are built with this technique. Numerical results show that the CSPM algorithm is consistent and demonstrates convergence with exact solutions. Main computational parameters are discussed, and different boundary conditions are validated to be effective for benchmark problems in computational acoustics.

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

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

The paper is dedicated to the topical problem of examining permafrost state and the processes of its geocryological changes by means of geophysical methods. To monitor the cryolithozone, we propose and scientifically substantiate a new technique of pulsed electromagnetic cross-well sounding. Based on the vector finite-element method, we created a mathematical model of the cross-well sounding process with a pulsed source in a three-dimensional spatially heterogeneous medium. A high-performance parallel computing algorithm was developed and verified. Through realistic geoelectric models of permafrost with a talik under a highway, constructed following the results of electrotomography field data interpretation, we numerically simulated the pulsed sounding on the computing resources of the Siberian Supercomputer Center of SB RAS. The simulation results suggest the proposed system of pulsed electromagnetic cross-well monitoring to be characterized by a high sensitivity to the presence and dimensions of the talik. The devised approach can be oriented to addressing a wide range of issues related to monitoring permafrost rocks under civil and industrial facilities, buildings and constructions.

Fire is a significant threat to human life and civil infrastructures. Builders and architects are hankering for safer and sustainable alternatives of concrete that do not compromise with their design intent or fire safety requirements. The aim of the present work is to improve the residual compressive performance of concrete in the post-fire exposure, by incorporating by-products from urban residues. Based on sustainability and circular economy motivations, the attention is focused on rubber tire fly ash, aged brick powder, and plastic (PET) bottle residuals used as partial sand replacement. The selected by-products from urban residues are used for the preparation of Cement-Based Composites (CBCs) in two different proportions (10 % and 15 %). The thermal CBC behaviour is thus investigated under realistic fire scenarios (i.e., Direct Flame (DF) for 1 hour (1h)), by following the ISO 834 standard provisions, but necessarily resulting in non-uniform thermal exposure for the cubic specimens. The actual thermal exposure is further explored with a Finite Element (FE) model, giving evidence of thermal boundaries effects. The post-fire residual compressive strength of heated concrete and CBC samples is hence experimentally derived, and compared to unheated specimens in ambient conditions. From the experimental study, the enhanced post-fire performance of CBCs with PET bottle residual is generally found superior to other CBCs or concrete. The structure-property relation is also established, with the support of Scanning Electron Microscopy (SEM) micrographs. Based on existing empirical models of literature for the prediction of the compressive or residual compressive strength of standard concrete, newly developed empirical relations for both concrete and CBCs are assessed.

A novel concept of ion separation Anisotropic Diffusion Mobility Analyser (ADMA) with Individual Ion Counter (IIC) was shown by numerical modelling to be able to separate ions with high resolution and ultra-high sensitivity. An ADMA-IIC device has been developed and tested at sheath flow rates below 40 l/min with a number of chemicals relevant to airport security and military applications. The resolving power of the ADMA was observed to be in the range of 40 to 80. These data were close to the resolution obtained in numerical modelling predictions. It is greater than or comparable with the RP of currently available on-the-market devices for airport security and other applications. It was also confirmed that the ADMA technology enables very fast detection of analytes, at least 100 times faster than ion mobility spectrometers. It was shown that Brownian diffusion can be controlled and even overpowered by generating highly anisotropic diffusion in the ADMA device. The ultimate sensitivity of single-ion detection below the ppq range has been demonstrated.

Modular steel buildings (MSBs) interlocking inter-module connections (IMCs) have piqued researchers’ interest. However, the existing literature has not yet studied the optimisation of its plate thicknesses. This paper, therefore, focuses on optimising interlocking (IMCs) plate thickness in (MSBs), leveraging previous experimental and numerical simulation methodologies. Various numerical models for four MSBs interlocking (IMCs) with plate thicknesses (4 mm, 6 mm, 10 mm, and 12 mm) were developed under compression loading conditions. The study’s novelty lay in its comprehensive parametric analysis, which evaluated the slip phenomenon in (MSBs) connections and introduced a slip prediction model validated against empirical data. Further innovative machine learning approach was utilised to predict slip values based on applied force through a random forest regression model trained using the 'Treebagger' function. Sensitivity analysis and comparisons with alternative methods were utilised to underscore the reliability and applicability of the findings. The results showed that 11.03 mm was the optimal plate thickness for interlocking (IMCs) in modular steel buildings, with up to 8.08% reduced material costs. Increasing plate thickness boosts its deformation resistance by up to 50.75%. ‘TreeBagger’ random forest for anomaly detection and machine learning improves slip prediction up to 7% at higher force levels.

Due to the ongoing global warming on the Earth, permafrost degradation has been extensively taking place, which poses a substantial threat to civil and industrial facilities and infrastructure elements, as well as to the utilization of natural resources in the Arctic and high-latitude regions. In order to prevent the negative consequences of permafrost thawing under the foundations of constructions, various geophysical techniques for monitoring permafrost have been proposed and applied so far: temperature, electrical, seismic and many others. At the same time, non-stationary electromagnetic methods seem to have significant potential and a number of important practical advantages. We propose a cross-borehole exploration system for the highest localization of target objects in the cryolithozone. A novel mathematical apparatus for three-dimensional modeling of transient electromagnetic (TEM) signals by the vector finite element method has been developed. The original combination of the latter, the Sumudu integral transform and artificial neural networks makes it possible to examine spatially heterogeneous objects of the cryolithozone with a high contrast of geoelectric parameters, significantly reducing computational costs. We have studied the sensitivity of the TEM signals to the boundary between thawed and frozen rocks depending on the inter-borehole distance and spatial orientations of the transmitters and receivers. Numerical simulation results of the TEM monitoring of industrial facilities located on permafrost are considered. The formation of a talik has been shown to significantly manifest itself in the measured electromagnetic responses, which enables timely prevention of industrial disasters and environmental catastrophes.

1960s Europe saw a large number of residential neighborhoods built to house those migrating from the countryside. Today, more than 50 years later, these neighborhoods suffer high levels of functional, social, and technical obsolescence. In response to this, University of Seville developed the Aura Strategy as an intervention methodology to find global solutions to issues in outdated neighborhoods. To provide visibility to this aspect of the Aura strategy, the retrofit proposal presented in this article provides a solution to improve the roofing of buildings in a case study neighborhood (Polígono de San Pablo, Seville) and an analysis of the results in terms of energy and financial savings for local residents. The results show that for a population of roughly 18,000 (in 2018), net savings, including energy and retrofitting costs, ranging from nearly 6.5 to over 8.6 million can be made over the 20-year life cycle span. Likewise, the results obtained on the reduction of thermal loads indicate a 72% decrease in energy consumption, equivalent to a saving of close to 4,500 tons of greenhouse gas emissions for the district and the entire life cycle time period, with the consequent benefits on the impact on air quality and and the fight against climate change. climate change.

Understanding the dynamic interaction between piles and surrounding soil under vehicular impacts is essential for effectively designing and optimizing soil-embedded vehicle barrier systems. The complex behavior of pile-soil systems under impact loading, attributed to the soil’s nonlinear behavior and large deformation experienced by both components, presents significant simulation challenges. Popular computation techniques, such as the Updated Lagrangian Finite Element Method (UL-FEM), encounter difficulties in scenarios marked by large soil deformation, e.g., impacts involving rigid piles. While mesh-free particle and discrete element methods offer another option, their computational demands for field-scale pile-soil impact simulations are considerable. We introduce the erosion method to bridge this gap by integrating UL-FEM with an erosion algorithm designed for simulating large soil deformations during vehicular impacts. Validation against established physical impact tests confirmed the method’s effectiveness for flexible and rigid pile failure mechanisms. Additionally, this method was used to investigate the effects of soil mesh density, soil domain sizes, and boundary conditions on the dynamic impact response of pile-soil systems. Our findings provide guidelines for optimal soil domain size, mesh density, and boundary conditions. This investigation sets the stage for improved, computationally efficient techniques for the pile-soil impact problem, leading to better pile designs for vehicular impacts.

The effects of wastewater step-feeding on the performance of Horizontal Subsurface Flow Constructed Wetlands (HSF CWs) are numerically investigated. The purpose is to check if this alternative feeding technique increases the ability of HSF CWs to remove pollutants. Two methodologies are used: Initially, the Tanks-In-series (TIS) methodology, based on Finite Volume Method (FVM), is analyzed by using the volumetric degradation coefficient λ. In this case the operation of CW is similar with a series of Continuous Stirred-Tank Reactors (CSTR), under steady conditions. Then, the Step-Feeding (SF) procedure is presented, when the CW is operated like a Plug Flow Reactor (PFR). The results show that the SF procedure does not improve the performance of HSF CWs on removing Biochemical Oxygen Demand (BOD), operating under Mediterranean conditions.

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

The World Health Organization estimates that every year air pollution kills seven million people worldwide. As it is expected that climate change will affect future air quality patterns, the full understanding of the links between air pollution and climate change, and how they affect human health, are challenges of future research. In this scope, a methodology to assess the air quality impacts on health was developed. The WRF-CAMx modelling framework was applied for the medium-term future climate (considering the SSP24.5 scenario) and for the recent past (considered as baseline). Following the WHO recommendations, mortality health indicators were used to estimate health impacts of long-term exposures. For that, the Aveiro Region, in Portugal, was considered as a case study. Future climate results indicate the occurrence of higher temperatures, and lower total precipitation. Despite that, improvements in the main pollutants’ concentrations, and consequently in the reduction of the related premature deaths are foreseen, mainly due to the reduction of pollutants emissions imposed by the European legislation for the upcoming years. The applied approach constitutes an added value in this research field, being crucial to anticipate the effects of climate change on air quality and evaluate their impacts on human health.

This paper reviews how dynamically adaptive wavelet methods can be designed to simulate atmosphere and ocean dynamics in both flat and spherical geometries. We highlight the special features that these models must have in order to be valid for climate modelling applications. These include exact mass conservation and various mimetic properties that ensure the solutions remain physically realistic, even in the under-resolved conditions typical of climate models. Particular attention is paid to the implementation of complex topography in adaptive models. Using \textsc{wavetrisk} as an example, we explain in detail how to build a semi-realistic global atmosphere or ocean model of interest to the geophysical community. We end with a discussion of the challenges that remain to developing a realistic dynamically adaptive atmosphere or ocean climate models. These include scale-aware subgrid scale parameterizations of physical processes, such as clouds. Although we focus on adaptive wavelet methods, many of the topics we discuss are relevant for adaptive mesh refinement (AMR).

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