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

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.

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.

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

To improve the effect of multi-point stretch forming of sheet metal, it is proposed in this paper to replace fixed ball head with swinging ball head. According to the multi-point dies with different arrangements, this research establishes the finite element models of the following stretch forming, i.e. fixed ball heads with conventional arrangement, swinging ball heads with conventional arrangement, swinging ball heads with declining staggered arrangement, and swinging ball heads with parallel staggered arrangement and then numerical simulation is performed. The simulation results show that by replacing fixed ball head with swinging ball head, the surface indentation of the formed part was effectively suppressed, the stress and tension strain distribution of the formed part was improved and the forming quality was improved; that the thickness of the elastic pad was reduced, the springback was reduced and the forming accuracy was improved; and that when the ball head was applied to multi-point die with staggered arrangement, better forming result was achieved, where the best forming result was achieved in combining the swinging ball heads with the multi-point die with parallel staggered arrangement. The forming experiments were carried out, and the experimental results were consistent with the trend of numerical simulation results, which verified the correctness of the numerical simulation.

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

Several problems of radiative transfer are yet unsolved because of the difficulties of the calculations involved in them, especially if the intervening shapes are geometrically complex. The main goal of our investigation in this domain is to convert the formulas that were previously derived, into a graphical interface based on the projected solid-angle principle. Such procedure is now feasible by virtue of several widely diffused programs for Algorithms Aided Design (AAD). Accuracy and reliability of the process is controlled by means of the analytical software DianaX developed at an earlier stage by the authors. With this new approach the often cumbersome procedure of lighting and thermal exchange calculations can be simplified and made available for the neophyte, with the undeniable advantage of reduced computer time.

Structural vibration induced by dynamic load or natural vibration is a nonnegligible factor in failure analysis. Based on vibrating boundary condition, impact resistance of shape memory alloy reinforced composites is investigated. In this investigation, modified Hashin’s failure criterion, Brinson’s model and visco-hyperelastic model are implemented into the numerical model to charactering the mechanical behavior of glass fiber/epoxy resin laminates, SMAs and interphase, respectively. First, fixed boundary condition is maintained in simulation to verify the accuracy of material parameters and procedures by comparing with experimental data. Then, a series of vibrating boundaries with different frequencies and amplitudes are applied during the simulation process to reveals the effect on impact resistances. The statistics of absorbed energy and contact force indicate that impact resistance of the composite under high frequency and large amplitude is lower than that under low frequency and small amplitude, and summarized by a mathematical expression.

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

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

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

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

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

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

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.

Wave is a common environmental load that often causes serious damages to offshore structures. In addition, the stability for the submarine artificial slope is also affected by the wave loading. Although the landslide of submarine slopes induced by the waves has received wide attention, the research on the influence of solitary-wave is rare. In this study, a 2-D integrated numerical model is developed to investigate the stability of the foundation trench under the solitary wave loading. The Reynolds-Averaged-Stokes (RANS) equations are used to simulate the propagation of a solitary wave, while the current is realized by setting boundary inlet/outlet velocity. The pore pressure induced by the solitary wave is calculated by Darcy’s law and the seabed is characterized by Mohr-Coulomb constitutive model. Firstly, the wave model is validated through the comparison between analytical solution and experimental data. The initial consolidation state of slope under hydrostatic pressure is achieved as the initial state. Then, the factor of stability (FOS) for the slope corresponding to different distance between wave crest and slope top is calculated with the strength reduction method. The minimum of FOS is defined as the stability index for the slope with specific slope ratio during the process of dynamic wave loading. The parametric study is conducted to examine the effects of soil strength parameters, slope ratio and current direction. At last, the influence of upper slope ratio in a two-stage slope is also discussed.

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

This review is devoted to present the history of CARDIOSIM© software simulator platform, which was developed in Italy to simulate the human cardiovascular and respiratory system. The first version of CARDIOSIM© was developed at the Institute of Biomedical Technologies of the National Research Council in Rome. The first platform version published in 1991 ran on PC with disk operating system (MS-DOS) and was developed using the Turbo Basic language. The last version runs on PC with Microsoft Windows 10 operating system; it is implemented in Visual Basic and C++ languages. The platform has a modular structure consisting of seven different general sections, which can be assembled to reproduce different pathophysiological conditions. The software simulator can reproduce the most important circulatory phenomena in terms of pressure and volume relationships. It represents the whole circulation using a lumped-parameter model and enables the simulation of different cardiovascular conditions according to Starling’s law of the heart and a modified time-varying elastance model. Different mechanical ventilatory and circulatory devices have been implemented in the platform including thoracic artificial lung, ECMO, IABP, pulsatile and continuous right and left ventricular assist devices, biventricular pacemaker and biventricular assist devices. CARDIOSIM© is used in clinical and educational environment.

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

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.

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

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

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

One of the most rapidly advancing areas of deep learning research aims at creating models that learn to disentangle the latent factors of variation from a data distribution. However, modeling joint probability mass functions is usually prohibitive, which motivates the use of conditional models assuming that some information is given as input. In the domain of numerical cognition, deep learning architectures have successfully demonstrated that approximate numerosity representations can emerge in multi-layer networks that build latent representations of a set of images with a varying number of items. However, existing models have focused on tasks requiring to conditionally estimate numerosity information from a given image. Here we focus on a set of much more challenging tasks, which require to conditionally generate synthetic images containing a given number of items. We show that attention-based architectures operating at the pixel level can learn to produce well-formed images approximately containing a specific number of items, even when the target numerosity was not present in the training distribution.

The emerging industry of offshore wind turbines mounted on floating bases has garnered significant attention from both academia and industry. The desire to understand the complex physics of these floating structures has led to the development of numerical and physical modelling techniques. While physical testing has traditionally been employed, there is a growing focus on cost-effective and accurate high-fidelity numerical modelling as a potential alternative or supplement. However, commonly used numerical engineering tools in the offshore industry are considered mid- to low-fidelity and may lack the desired precision for floating offshore wind turbines (FOWTs). Given the complexity of these simulation codes, it is crucial to validate their accuracy. To address this, the International Energy Agency (IEA) Wind Technology Collaboration Programme initiated various research endeavors, including the Offshore Code Comparison Collaboration (OC3), Offshore Code Comparison Collaboration Continuation (OC4), Offshore Code Comparison Collaboration Continuation with Correlation (OC5), and the recent Offshore Code Comparison Collaboration Continued with Correlation and Uncertainty (OC6) projects. This study offers a comprehensive survey of the simulation tools available for FOWTs which were part of OC projects, focusing particularly on horizontal axis wind turbines (HAWTs) and highlighting their capabilities and fundamental theories.

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.

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

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.

For embankment slope stability, soil cohesivity is one of the most important shear strength parameters. The effect of cohesiveness on the slope stability of a homogeneous embankment dam under rapid drawdown loading conditions was examined in this study. With the use of numerical modeling in GeoStudio, different situations were explored based on cohesiveness (0 kN/m², 5 kN/m², 10 kN/m^2, and 15 kN/m²) under a 1 m per day drawdown rate. The factor of safety value obtained from the long-term steady-state condition under 15 kN/m² cohesion was equivalent to a 116.8% increase from the one obtained under 0 kN/m² cohesion. The factor of safety values obtained after subjecting the embankment to different soil cohesion levels yielded a p-value of 1.91 x 10^-41, according to the Analysis of Variance. The calculated p-value (alpha value) is less than 0.05, suggesting that the differences between the examined cohesiveness values based on the list of the factor of safety values are statistically significant. The findings derived in this study show that it is significant to capture the effect of material characteristics during the design phase of an embankment dam.

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

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.

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.

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

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.

Cladding is typically used to protect components from wear and corrosion while also improving the aesthetic value and reliability of the substrate. The cladding process induces significant residual stresses due to the temperature difference between the substrate and the clad layer. However, these residual stresses could be effectively utilized by modifying processes and geometrical parameters. This paper introduces a novel methodology for using the weld-cladding process as a cost-effective alternative to various existing reinforcement techniques. The numerical analyses are performed to maximize the reinforcement of a cylindrical tool. The investigation of how the weld cladding develops compressive stresses on the specimen in response to a change in the weld beads and the welding sequence is presented. For the benchmark shape, experimental verification of the numerical model is performed. The impact of the distance between the weld beads and the effect of the tool diameter is numerically investigated. Furthermore, the variation in compressive stresses due to temperature fluctuations during the extrusion process has been evaluated. The results showed that adequate compressive stresses are generated on the welded parts through the cladding process after cooling. Hence, the targeted reinforcement of the substrate can be achieved by optimizing the welding sequence and process parameters.

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

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

Recently, more wind turbine systems are being installed in deep waters far from the coast. Several concepts of floating wind turbine systems (FWTSs) are developed, among which the semisubmersible platform, due to its applicability in different water depths, good hydrodynamic performance, and facility in the installation process, constitutes the most explored technology compared to the others. However, a significant obstacle to the industrialization of this technology is a design of a cost-effective FWTS, which can be achieved by optimizing the geometry, size and weight of the floating platform, along with the mooring system. This is only possible by selecting a method capable of accurately analyzing the FWTS coupled hydro-aero-structural dynamics at each design stage. Accordingly, this paper aims at providing a detailed overview of the most common coupled numerical and physical methods, including their basic assumptions, formulations, limitations, and costs, used for analyzing the dynamics of FWTSs, mainly those supported by a semisubmersible, to assist the choice of the most suitable method at each design phase of FWTSs. Finally, the article discusses possible future research to address challenges in modeling FWTSs dynamics that persist to date.

Pericytes are increasingly recognized as being important in the control of blood-brain barrier permeability and vascular flow. Research on this important cell type has been hindered by widespread confusion regarding the phenotypic identity and nomenclature of pericytes and other perivascular cell types. In addition, pericyte heterogeneity and mouse-human species differences have contributed to confusion. Herein we summarize our present knowledge on the identification of pericytes and pericyte subsets in humans, primarily focusing on recent findings in humans and nonhuman primates. Precise identification and definition of pericytes and pericyte subsets in humans may help us to better understand pericyte biology and develop new therapeutic approaches specifically targeting disease-associated pericyte subsets.

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

A key feature of pulmonary arterial hypertension (PAH) is the hyperplastic proliferation exhibited by the vascular smooth muscle cells from patients (HPASMC). The growth inducers FOXM1 and PLK1 are highly upregulated in these cells. The mechanism by which these two proteins direct aberrant growth in these cells is not clear. Herein we identify cyclin dependent kinase 1 (CDK1) also termed cell division cycle protein 2 (CDC2), as having a primary role in promoting progress of the cell cycle leading to proliferation in HPASMC. HPASMC obtained from PAH patients and pulmonary arteries from Sugen/hypoxia rats were investigated for their expression of CDC2. Protein levels of CDC2 were much higher in PAH than in cells from normal donors. Knocking down FOXM1 or PLK1 protein expression with siRNA or pharmacological inhibitors lowered the cellular expression of CDC2 considerably. However, knockdown of CDC2 with siRNA or inhibiting its activity with RO-3306 did not reduce the protein expression of FOXM1 or PLK1. Expression of CDC2 and FOXM1 reached its maximum at G1/S, while PLK1 reached its maximum at G2/M phase of the cell cycle. The expression of other CDKs such as CDK2, CDK4, CDK6, CDK7 and CDK9 did not change in PAH HPASMC. Moreover, inhibition via Wee1 inhibitor adavosertib or siRNAs targeting Wee1, Myt1, CDC25A, CDC25B, or CDC25C led to dramatic decreases in CDC2 protein expression. Lastly, we found CDC2 expression at the RNA and protein level to be upregulated in pulmonary arteries during disease progression Sugen/hypoxia rats. In sum, our present results illustrate that the increased expression of FOXM1 and PLK1 in PAH leads directly to increased expression of CDC2 resulting in a potentiated growth hyperactivity of PASMC from patients with pulmonary hypertension. Our results further suggest that the regulation of CDC2, or associated regulatory proteins, will prove beneficial in the treatment of this disease.

This study employs a high-fidelity numerical framework to determine the plastic material flow patterns and temperature distributions that lead to void formation during friction stir welding (FSW), and to relate the void morphologies to the underlying alloy material properties and process conditions. Three aluminum alloys, viz., 6061-T6, 7075-T6, and 5053-H18 were investigated under varying traverse speeds. The choice of aluminum alloys enables investigation of a wide range of thermal and mechanical properties. The numerical simulations were validated using experimental observations of void morphologies in these three alloys. Temperatures, plastic strain rates, and material flow patterns are considered. The key results from this study are: (1) The predicted stir zone and void morphology are in good agreement with the experimental observations, (2) The temperature and plastic strain-rate maps in the steady-state process conditions show a strong dependency on the alloy type and traverse speeds, (3) The material velocity contours provide a good insight into the material flow in the stir zone for the FSW process conditions that result in voids as well as those that do not result in voids. The numerical model and the ensuing parametric studies presented in this work provide a framework for understanding material flow under different process conditions in aluminum alloys, and potentially in other alloys. Furthermore, the utility of the numerical model for making quantitative predictions and investigating different process parameters to reduce void formation is demonstrated.

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.

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.

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

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

In this paper, we investigate the accuracy and robustness of three classes of methods for solving two-phase incompressible flows on a staggered grid. Here, the unsteady two-phase flow equations are simulated by finite volumes and penalty methods using implicit and monolithic approaches (such as the augmented Lagrangian and the fully coupled methods), where all velocity components and pressure variables are solved simultaneously (as opposed to segregated methods). The interface tracking is performed with a Volume-of-Fluid (VOF) method, using the Piecewise Linear Interface Construction (PLIC) technique. The home code Fugu is used for implementing the various methods. Our target application is the simulation of two-phase flows at high density and viscosity ratios, which are known to be challenging to simulate. The resulting strategies of monolithic approaches will be proven to be considerably better suited for these two-phase cases, they also allow to use larger time step than segregated methods.

Chronic pulmonary diseases such as asthma, COPD, and Idiopathic pulmonary fibrosis are significant causes of mortality and morbidity worldwide. Currently, there is no radical treatment for many chronic pulmonary diseases, and the treatment options focus on relieving the symptoms and improving lung function. Therefore, efficient therapeutic agents are highly needed. Bronchial epithelial cells and airway smooth muscle cells and their crosstalk play a significant role in the pathogenesis of these diseases. Thus, targeting the interactions of these two cell types could open the door to a new generation of effective therapeutic options. However, the studies on how these two cell types interact and how their crosstalk adds up to respiratory diseases are not well established. With the rise of modern research tools and technology, such as lab-on-chip, organoids, co-culture techniques, and advanced immunofluorescence imaging, a substantial degree of evidence about these cell interactions emerged. Hence, this contribution aims to summarize the growing evidence of bronchial epithelial cells and airway smooth muscle cells crosstalk under normal and pathophysiological conditions. The review first deliberates the effects of both healthy and stressed epithelial cells on airway smooth muscle cells, taking into account three themes; contraction, migration, and proliferation. Then, it discusses the impact of airway smooth muscle cells on the epithelium in inflammatory settings. Later, it examines the role of airway smooth muscle cells in the early development of bronchial epithelial cells and their recovery after injury.

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

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

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

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

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

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

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

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

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

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.

Gene electrotransfer is one of the main non-viral methods for intracellular delivery of plasmid DNA, wherein pulsed electric fields are used to transiently permeabilize the cell membrane allowing enhanced transmembrane transport. By localizing the electric field over small portions of the cell membrane using nanostructured substrates, it is possible to increase considerably the gene electrotransfer efficiency while preserving cell viability. In this study, we design an electrotransfer approach based on commercially available cell culture inserts with polyethylene-terephthalate (PET) porous substrate. We first use multiscale numerical modelling to determine the pulse parameters, substrate pore size, and other factors that are expected to result in successful gene electrotransfer. Based on numerical results we design a simple device combining an insert with substrates containing pores with 0.4 um and 1.0 um diameter, a multiwell plate, and a pair of wire electrodes. We test the device in three mammalian cell lines and obtain transfection efficiencies similar to those achieved with bulk electroporation, but with low-voltage pulses that do not require the use of expensive electroporators. Our combined theoretical and experimental analysis calls for further systematic studies that will investigate the influence of substrate pore size and porosity on gene electrotransfer efficiency and cell viability.

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

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

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.

Cellular senescence is characterized by permanent proliferation and migration arrest, senescence-associated secretory phenotype (SASP), and oxidative stress. Senescent vascular smooth muscle cells (VSMCs) contribute to cardiovascular diseases and atherosclerotic plaque instability. To establish and characterize a model of replicative senescence (RS), human aortic VSMCs were serially passaged to represent different stages of RS such as young proliferating cells and old/senescent non-proliferating cells. More than 50% of old cells stained positive for the senescence-associated β-galactosidase compared to 20% of young cells. Old cells have a slower proliferation rate, a migratory activity reduced by 50%, but increased levels of TP53 and of cell cycle inhibitors p21/p16 expression, and accumulate in the G1 phase. Old cells showed a flattened appearance and enlarged and regular nuclei, and downregulation of the expression of LMNB1 and HMGB1. Old cells showed also an increased expression of SASP molecules (IL1β, IL6, IL8, and MMP3). Moreover, among a set of 12 manually selected long non-coding RNAs (lncRNAs), we detected significant upregulation of PURPL and NEAT1. We observed also increased levels of RRAD mRNA. The detection of novel molecular markers of senescence, such as RRAD, PURPL, and NEAT1, could be helpful for future studies on potential anti-aging factors.

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.

Mutations in the LMNA gene (encoding lamin A/C proteins) cause several human cardiac diseases, including dilated cardiomyopathies (LMNA-DCM). The main clinical risks in LMNA-DCM patients are sudden cardiac death and progressive left ventricular ejection fraction deterioration, and therefore most human and animal studies have sought to define the mechanisms through which LMNA mutations provoke cardiac alterations, with particular focus on cardiomyocytes. To investigate if LMNA mutations also cause vascular alterations that might contribute to the etiopathogenesis of LMNA-DCM, we generated and characterized Lmnaflox/floxSM22αCre mice, which constitutively lack lamin A/C in vascular smooth muscle cells (VSMCs), cardiac fibroblasts, and cardiomyocytes. Like mice with whole body or cardiomyocyte-specific lamin A/C ablation, Lmnaflox/floxSM22αCre mice recapitulated the main hallmarks of human LMNA-DCM, including ventricular systolic dysfunction, cardiac conduction defects, cardiac fibrosis, and premature death. These alterations were associated with hyperactivation of Smad3 and elevated expression of the proapoptotic protein caspase 3 in the heart. Lmnaflox/floxSM22αCre mice also exhibited perivascular fibrosis in the coronary arteries, and a switch of aortic VSMCs from the ‘contractile’ to the ‘synthetic’ phenotype. Ex vivo wire myography in isolated aortic rings revealed impaired maximum contraction capacity and an altered response to vasoconstrictor and vasodilator agents in Lmnaflox/floxSM22αCre mice. To our knowledge, our results provide the first evidence of phenotypic alterations in VSMCs that might contribute significantly to the pathophysiology of some forms of LMNA-DCM. Future work addressing the mechanisms underlying vascular defects in LMNA-DCM may open new therapeutic avenues for these diseases.

One of the most important subjects of hydraulic engineering is the reliable estimation of the transverse distribution in rectangular channel of bed and wall shear stresses. This study makes use of the Tsallis entropy, Genetic Programming (GP) and (ANFIS) methods to assess the shear stress distribution (SSD) in rectangular channel. To evaluate the results of the Tsallis entropy, GP and ANFIS models, laboratory observations were used in which shear stress was measured using an optimized Preston tube. This is then used to measure the SSD in various aspect ratios in the rectangular channel. To investigate the shear stress percentage, 10 data series with a total of 112 different data for were used. The results of the sensitivity analysis show that the most influential parameter for the SSD in smooth rectangular channel is the dimensionless parameter B/H, Where the transverse co-ordinate is B, and the flow depth is H. With the parameters (b/B), (B/H) for the bed and (z/H), (B/H) for the wall as inputs, the modeling of the GP was better than the other one. Based on the analysis, it can be concluded that the use of GP and ANFIS algorithms is more effective in estimating shear stress in smooth rectangular channels than the Tsallis entropy-based equations.

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.

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

This paper proposes and demonstrates, in full scale, a novel type of energy geostructure (“the Climate Road”) that combines a ground source heat pump (GSHP) with a sustainable urban drainage system (SUDS) by utilizing the gravel roadbed simultaneously as energy source and rainwater retarding basin. The Climate Road measures 50m x 8m x 1m (length, width, depth) and has 800 m of geothermal piping embedded in the roadbed, serving as the heat collector for a GSHP that supplies a nearby kindergarten with domestic hot water and space heating. Model analysis of operational data from 2018-2021 indicates sustainable annual heat production levels around 0.6 MWh per meter road, with a COP of 2.9-3.1. The continued infiltration of rainwater to the roadbed increases the amount of extractable heat by an estimated 17% compared to the case of zero infiltration. Using the developed model for scenario analysis we find that draining rainwater from three single family houses and storing 30% of the annual heating consumption in the roadbed, increases the predicted extractable energy by 56% compared to zero infiltration with no seasonal energy storage. The Climate Road is capable of supplying three single family houses with heating, cooling and rainwater management year-round.

This paper offers a comprehensive examination of the process involved in developing and automating supervised end-to-end machine learning workflows for forecasting and classification purposes. It offers a complete overview of the components (i.e. feature engineering, model selection, etc), principles (i.e. bias-variance decomposition, model complexity, overfitting, model sensitivity to feature assumptions and scaling, output interpretability, etc), models (i.e. neural networks, regression models, etc), methods (i.e. Cross-Validation, data augmentation, etc), metrics (i.e. Mean Squared Error, F1-score, etc) and tools that rule most supervised learning applications with numerical and categorical data, as well as their integration, automation, and deployment. The end goal and contribution of this paper is the education and guidance of the non-AI expert academic community over complete and rigorous machine learning pipelines and data science practices, from problem scoping to design and state-of-the-art automation tools, including basic principles and reasoning in the choice of methods. The paper delves into the critical stages of supervised machine learning workflow development, many of which are often omitted by researchers due to brevity, and covers foundational concepts essential for understanding and optimizing a functional machine learning workflow, thereby offering a holistic view of task-specific application development for applied researchers who are not AI experts.

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

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.

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.

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.

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

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

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

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.

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.

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

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.

The concepts of fractional calculus and the related techniques have been gaining momentum in circuit system fields, and an increasing body of evidence suggests that fractional-order characteristics are widely distributed in electronic components. In this paper, we consider the fractional characteristics of non-solid aluminum electrolytic capacitors and focus on the nonlinear dynamic characteristics of peak-current-mode controlled (PCMC) dc-dc boost converters with such components. In the work, we establish a piecewise smooth non-commensurate fractional-order model to describe the converter, in which the fractional-order equivalent impedance model of capacitors is adopted and the detailed control loop is considered. In order to analyze the converter, we use the fractional Adams-Bashforth-Moulton typed method (F-ABM) and stroboscopic map technique to obtain the time-domain solutions and the bifurcation diagram. The results of bifurcation diagrams indicate that, with the change of different parameters, the converter enters chaotic state through period doubling bifurcation. In order to provide empirical evidence for the nonlinear dynamic analysis, we consider both circuit-level simulation and experiments in the work, the results of which further confirm the results of theoretical analysis.

The potentially active A subunit of coagulation factor XIII (FXIII-A) is an intracellular transglutaminase expressed in various cell types including platelets and monocytes/macrophages. It is involved in stabilizing protein structures by cross-linking through Nε-(?-L-glutamyl)-L-lysyl iso-peptide bonds. Macrophages are major cellular constituents of the atherosclerotic plaque and are important in determining its structural/functional features. Two of their important functions are the accumulation of oxidized LDL in the lipid core, and by cross-linking structural proteins they may stabilize the plaque and protect the thrombi of atherogenic origin against fibrinolytic degradation. It is important to know whether these functions operate in parallel utilizing the same cellular compartments. First, we showed that monocyte-derived human macrophages significantly increase their FXIII-A content when up-taking oxidized LDL. This phenomenon is very likely independent of the process of transformation into foam cells, as the transformation of vascular smooth muscle cells into foam cells fails to result in the expression of FXIII-A. FXIII containing macrophage-like cells are abundant in the plaque and FXIII-A is also present in the extracellular core. Several cells co-stained for FXIII-A and for Oil Red O suggest that expression of FXIII-A and lipid up-take are common features of macrophages present in the atherosclerotic plaque.

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

H-adaptivity is an effective tool to introduce local mesh refinement in FEM-based numerical simulation of crack propagation. The implementation of h-adaptivity could benefit the numerical simulation of fatigue or accidental load scenarios involving large structures such as ship hulls. In engineering applications, the element deletion method is frequently used to represent cracks. However, the element deletion method has some drawbacks such as strong mesh dependency and loss of mass or energy. In order to mitigate this problem, the element splitting method could be applied. In this study, a numerical method called ‘h-adaptive element splitting’ (h-AES) is introduced. The h-AES method is applied in FEM programs by combining h-adaptivity with the element splitting method. Two examples using the h-AES method to simulate cracks in large structures under linear-elastic fracture mechanics scenario are presented. The numerical results are verified against analytical solutions. Based on the examples, the h-AES method is proven to be able to introduce mesh refinement in large-scale numerical models that consist of structured coarse meshes. By employing the mesh refinement introduced in this paper, very small cracks are well represented in large structures.

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.

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

The role of fractional calculus in circuit systems has received increased attention in recent years. In order to evaluate the effectiveness of time-domain calculation methods in the analysis of fractional-order piecewise smooth circuit systems, an experimental prototype is developed and the effects of three typical calculation methods in different test scenarios are compared and studied in this paper. It is proved that Oustaloup’s rational approximation method usually overestimates the peak-to-peak current and brings in pulse-voltage phenomenon in piecewise smooth test scenarios, while the results of two iterative recurrence-form numerical methods are in good agreement with the experimental results. The study results are dedicated to provide a reference for efficiently deploying calculation methods in fractional-order piecewise smooth circuit systems. Some quantitative analysis results are concluded in this paper.

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

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

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

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

Asthma is a common chronic respiratory disease that affects millions of people worldwide, and its prevalence continues to increase. Vitamin D has been proposed as a potential environmental factor in asthma pathogenesis, due to its immunomodulatory effects. This systematic review aimed to evaluate the effect of vitamin D supplementation on airway remodeling in asthma. We conducted a comprehensive search of electronic databases and identified experimental studies and randomized controlled trials that investigated the impact of vitamin D supplementation on airway remodeling in asthma. The studies included in this review suggest that vitamin D inhibits airway smooth muscle cell contraction and remodeling, reduces inflammation, and regulates collagen synthesis in the airways, and modulates the action of bronchial fibroblasts. However, one study suggests that TGF-β1 can impair vitamin D-induced and constitutive airway epithelial host defense mechanisms. Overall, vitamin D appears to influence many of the pathways central to asthma progression, indicating a promising potential in asthma management and prevenation of asthma progression.

This paper studies the evaporative heat transfer characteristics of R744 at low temperatures in a horizontal smooth tube as a cascade refrigeration system (CRS) among hybrid cascade refrigeration systems (HCRSs). There is a lack of research on the low-temperature evaporative heat transfer characteristics of R744 under the operating conditions of evaporators used in actual CRSs used in supermarkets. Therefore, this study aims to provide basic data on the evaporative heat transfer characteristics of R744 in the evaporators of refrigerators used in supermarkets. The tube used in the evaporation experiment conducted herein was a smooth horizontal copper tube with an inner diameter and length of 11.46 mm and 8000 mm, respectively. The experimental pa-rameters were as follows: heat fluxes of 12–21.5 kW/m2, mass fluxes of 75–225 kg/(m2·s), and saturation temperatures of −50–−30 °C. The main results are summarized as follows. (1) When designing the R744 evaporator, the mass and heat fluxes must be maximized within the operating conditions, and the saturation temperature must be designed to be as low as possible. (2) The evaporative heat transfer coefficient of R744 can be predicted well by using the correlation formula of Chen at the evaporation temperature of −40 °C in the CRS.

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

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.

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

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