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.

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.

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

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.

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.

The properties of iterative splitting methods with serial versions have been analyzed since recent years, see [1] and [3]. We extend the iterative splitting methods to a class of parallel versions, which allow to reduce the computational time and keep the benet of the higher accuracy with each iterative step. Parallel splitting methods are nowadays important to solve large problems, which can be splitted in subproblems and computed independently with the dierent processors. We present a novel parallel iterative splitting method, which is based on the multi-splitting methods, see [2], [10] and [15]. Such a exibilisation with multisplitting methods allow to decompose large iterative splitting methods and recover the benet of their underlying waveform-relaxation (WR) methods. We discuss the convergence results of the parallel iterative splitting methods, while we could reformulate such an error to a summation of the individual WR methods. We discuss the numerical convergence of the serial and parallel iterative splitting methods and present dierent numerical applications to validate the benet of the parallel versions.

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

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

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

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.

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.

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.

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

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

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

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.

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

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

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

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.

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.

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

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

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.

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.

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.

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

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

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

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

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

The object of this work is to determine the correlation on the Nusselt number on the individual rows of a four-row tubular finned heat exchanger with continuous fins with a staggered tube arrangement using CFD modelling. Correlations for calculating Darcy-Weisbach friction factors on individual tube rows were also determined. Relationships for the Nusselt number and friction factor derived for the entire exchanger based on CFD modelling were compared with those available in the literature determined using experimental data. The maximum relative differences between the Nusselt number for a four-row exchanger determined experimentally and by CFD modelling are in the range from 22% for a Reynolds number based on a tube's outside diameter of 1,000 to 30% for a Reynolds number of 13,000. The maximum relative differences between the friction factor for a four-row exchanger determined experimentally and by CFD modelling are in the range of 50% for a Reynolds number based on a tube outer diameter of 1,000 to 10% for a Reynolds number of 13,000. The CFD modeling performed shows that in the range of Reynolds numbers based on hydraulic diameters from 150 to 1,400, the Nusselt number for the first row in a four-row finned heat exchanger is about 22% to 15% higher than the average Nusselt number for the entire exchanger. In the range of Reynolds number changes based on hydraulic diameter from 2,800 to 6,000, the Nusselt numbers on the first and second rows of tubes are close to each other. Correlations on Nusselt numbers and friction factors derived for individual tube rows can be used in the design of plate-fin and tube heat exchangers used in equipment such as air-source heat pumps, automotive radiators, air-conditioning systems and in air hot-liquid coolers. In particular, the correlations can be used to select the optimum number of tube rows in the exchanger.

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

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.

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

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.

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.

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.

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.

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

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.

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.

Data obtained by Direct Numerical Simulations (DNS) of pressure-driven turbulent channel flow are studied in the range 180 ≤ Reτ≤ 10,000. Reynolds number effects on the Mean Velocity Profile (MVP) and second order statistics are analyzed with a view of finding logarithmic behavior in the overlap region or even further from the wall, well in the boundary layer’s outer region. Values of Kármán constant for the MVPs and the Townsend-Perry constants for the streamwise and spanwise fluctuation variances are determined for the Reynolds numbers considered. A data-driven model of the MVP, proposed and validated for zero pressure gradient flow over a flat plate, is employed for pressure-driven channel flow by appropriately adjusting Coles’ strength of the wake function parameter, Π. There is excellent agreement between the analytic model predictions of MVP and the DNS-computed MVP as well as of the Reynolds shear stress profile. The skin friction coefficient Cf is calculated analytically. The agreement between the analytical model predictions and the DNS-based computed discrete values of Cf is excellent.

In the present study, a simple approximation expression is given for the relationship between the period and amplitude of a simple pendulum under magnetic action. The analytical solution presented for the given problem. Two numerical quadrature methods Simpson's and Boole's method were utilized to demonstrate a new approximation of the problem. The results of the numerical quadrature have been compared to the exact solution. Absolute and relative mistakes of the problem have been presented. The Matlab program 2013R has created a numerical method that is used to analyze the outcome, It has been determined that the comparison's outcomes attest to the method's suitability and correctness. Moreover, the results show that numerical solution is suitable for the problem.

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.

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

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.

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.

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.

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.

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

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

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

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.

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.

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.

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.

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.

This study concerns the transport and retention of polydisperse micron-sized (16 ± 6 µm) of microplastics (MPls) in porous media under varying flow rate conditions. Sorption kinetics were modeled using first-order reversible and irreversible kinetic sorption models, with sensitivity analyses providing insights into each sorption parameter's effect. Both numerical modeling and experimental measurements were employed to assess sand filter retention rates. The impact of flow rate on sorption reveals variations in distribution coefficient (Kd), mass transfer coefficient (β), and irreversible sorption rate (K1). Lower flow rates are correlated with higher Kd and β values, indicating an increase in sorption and diminished mass transfer rates. The findings revealed that an increase in Kd resulted in a more gradual sorption process with a decrease in peak concentration, whereas changes in β influenced the rate of sorption and peak concentration to a lesser extent compared to Kd. Lower K1 values are associated with higher peak concentrations and decreased retention efficiency. Retention rates were evaluated by a numerical model and found as 28 ± 1% at a flow rate of 31 ml min⁻¹ and 17 ± 1% at 65 ml min⁻¹. The introduction of MPls into soil environments has been noted to modify transport dynamics into soil. As a result, these alterations effects hydrological characteristics of soil, thereby impacting quality of groundwater and agricultural output. The mean absolute error (MAE) of 6% between modeled and observed retention rates suggests minor discrepancies. This study highlights the importance of examining retention efficiency and the accuracy of numerical models in porous media during MPl transport.

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.

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

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.

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.

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

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

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

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.

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

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

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.

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.

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.

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.

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.

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.

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

The steel grid structure with bolted ball joints is gaining popularity in offshore platforms and aquaculture cages, but the corrosion resistance of these joints is a critical safety concern. This article introduces an innovative method to enhance the sealing performance of bolted ball joints, preventing seawater ingress and corrosion. The approach involves creating a sealed surface within the contact gap between the sleeve and connecting components. Circular grooves and sealing washers are added to both ends of the sleeve. A two-dimensional finite element analysis model of the bolted ball joint with the sealing structure was created using SOLIDWORKS and ANSYS. The study analyzes the sealing gasket's contact pressure at various compression levels and evaluates its performance with bubble tests for air tightness. Research results show a linear relationship between contact pressure and compression rate, with a sealing pressure of 5.95MPa achieved at a 14% compression rate. This simple, cost-effective sealing method effectively prevents corrosion from seawater, making it ideal for marine environments.

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

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.

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

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

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.

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

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.

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

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

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

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

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

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

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

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

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

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.

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