Digital twins are the foundation of autonomous off-road vehicles. Dynamic reduction methods are one of several ways to develop digital twins for off-road vehicles. The article commences with a comprehensive overview of the most widely used dynamic reduction methods and then introduces performance metrics for assessing their efficacies in the context of digital twins. The paper additionally includes a detailed mathematical derivation of the state-space representation for reduced-order finite element models. The state-space representation of the reduced finite element models facilitates their export to problem-solving environments for dynamic analysis. The state-space models are eventually solved utilizing the built-in libraries of numerical solvers in textual and graphical programming platforms. In addition, the article identifies the set of solvers that best suit the simulation of virtual models for off-road vehicles. This article also includes an evaluation of the simulation results of digital models with modes ranging from 0 to 30 Hz. In addition, the article demonstrates the lower bound of the frequency range necessary and sufficient to be retained in off-road vehicle virtual models. Finally, the paper presents the simulation outcomes for digital models of commercial off-road vehicles with custom-built control, electrical, and control systems.

Keyhole laser beam welding (LBW) of stainless steel sheets with a gap in between is numerically simulated with a three-dimensional, transient multi-physical model for laser material processing. At first, the model’s ability to reproduce experimental results on a relatively coarse computational mesh within reasonable computing time, so as to serve as process optimization tool, is presented. An example of process optimization, where a given set of weld seam quality criteria is fulfilled by iteratively optimizing a secondary laser beam, is shown. The relatively coarse mesh, in combination with a good model calibration to the experimental conditions, allows for sufficiently fast simulations to use this approach for optimization tasks. Finally, using a finer spatial and temporal discretization, the dynamic processes in the vicinity of the keyhole leading to the formation of pores are investigated. The physical phenomena predicted by the simulation are coherent with experimental observations found in literature.

Whether the universe is a computer simulation, or whether we wish to efficiently model our universe in a computer simulation, there would be benefits to modeling it in a fashion analogous to computer spreadsheet, each lattice cell can be conceived as containing all the mathematical formula necessary to continuously compute its state relative to changes in all its neighboring cells, and by progression, in relation to all the cells of entire space-time lattice. Alternatively, the “real” universe may itself be built on a space cell lattice, an irregular foam of space cells, in which each cell may be conceived as a multidimensional cell of distortable space, the shape of which fully describes (a) the four basic forces (gravity, electromagnetic, strong, weak) observed at that cell of space, and (b) the probability (or weight distribution) of any quantum states overlapping the cell and its neighbors. At an appropriate scale, it would appear that this conceptual model would resolve apparent conflicts between general relativity and quantum physics. It would also provide a new interpretation of Planck’s constant as description of the number of space cell events associated with any set of observable events. If formulae operating at a lattice cell level can be improve our ability to understand and model larger scale phenomena, this would be strong evidence in favor of the theory that mathematics is not just a human invention but rather an inherent feature of space-time itself.

Risk-based redispatch optimization is proposed as a methodology to support the Transmission System Operator (TSO) with preventive remedial actions obtained by extending the security constrained unit commitment/economic dispatch with constraints resulting from the risk assessed for the power system. Although being heuristic, the methodology is based on comprehensive dynamic security assessment as time-domain simulations are used, allowing to express the degree of all types of instabilities, e.g., caused by contingencies, in monetary terms. Therefore, the risk is assessed as the expected value of the cost incurred by the TSO. Such approach forms a new pathway to including risk in planning procedures already used by TSOs. Results obtained for the IEEE39 dynamic power system, with costs assigned to load shedding and generator tripping due to single transmission lines short-circuits, are shown as a reference case.

X-ray fluorescence is largely employed in the measurement of the thickness of coatings. Despite of its diffusion, the task is not straightforward because of the complex physics involved that results in high dependence on matrix effects. Thickness quantification is in practice accomplished using the Fundamental Parameters approach, adjusted with empirical measurements of standards with known composition and thickness. This approach has two major drawbacks: i) there are no standards for any possible coating and coating architecture and ii) even relying on standards, the quantification of unknown samples requires the precise knowledge of the matrix nature (e.g., in case of multilayer coatings the thickness and the composition of each underlayer). In this work, we describe a semiquantitative approach to coatings thickness measurement based on the construction of calibration curves through simulated XRF spectra built with Monte Carlo simulations. Simulations have been performed with the freeware software XMI-MSIM. We have assessed the accuracy of the methods by comparing the results with those obtained by i) XRF thickness determination with standards and ii) FIB-SEM cross-sectioning. Then we evaluated which parameters are critical in this kind of indirect thickness measurements.

In this paper, we present a comprehensive 4D study of the early stage of plastic deformation in a polycrystalline binary AlLi alloy. The entire microstructure is mapped with X-ray diffraction contrast tomography and a set of bulk grains is further studied via X-ray topotomography during mechanical loading. The observed contrast is analyzed with respect to the slip system activation and the evolution of the orientation spread is measured as a function of applied strain. The experimental observations are augmented by the mechanical response predicted by crystal plasticity finite element simulations to analyze the onset of plasticity in details. Simulation results show a general agreement of the individual slip system activation during loading and that comparison with experiments at the length scale of the grains may be used to fine tune the constitutive model parameters.

The heavy demand of water resources from the Zarqa River Basin (ZRB) has resulted in a base-flow reduction of the River from 5m³/s to less than 1m³/s. This paper aims at predicting Curve Numbers (CN) as a baseline scenario and proposing restoration scenarios for ZRB. The method includes classifying the soil type and land use, predicting CNs, and proposing CN restoration scenarios. Prediction of existing CNs will be in parallel with the runoff prediction using the US Army Corps of Engineers HEC-1 Model, and Rainfall-Runoff Model (RRM). The models have been set up at the land use distribution of 0.3% water body, 9.3% forest and orchard, 71% mixture of grass, weeds, and desert shrubs, 7.0% crops, 4.0% urban areas, and 8.4% bare soil. The results show that CN under dry condition are 59, 78 under a normal condition and 89 under a wet condition. During vegetation period, CN are 52, 72 and 86 for the dry, normal and wet condition respectively. The restoration scenarios; CN decreases runoff, and increases soil moisture when using the contours, terraces and crop residues. Analyzing results of CN scenarios will be a fundamental tool to support end-users related to their targets to achieve watershed restoration.

Previous research has shown that the effectiveness of selecting filter set from a large set of commercial broadband filters by vector analyzing method based on maximum linear independence (MLI). However, the traditional MLI is suboptimal due to predefining the first filter of the selected filter set being the maximum ℓ2 norm among all those of the filters. An exhaustive imaging simulation is conducted to investigate the features of the most competent filter set. In the simulation, every filter in a subset of all the filters is selected serving as the first filter in turn. From the results of the simulation, we found there are remarkable characteristics for the most competent filter set. Besides smaller condition number, the geometric features of the best-performed filter set comprise the distinct peak of the transmittance of the first filter, the generally uniform distributing of the peaks of the transmittance curve of the filters, the substantial overlapping of the transmittance curves with those of the adjacent filer sets. Therefore, the best-performed filter sets can be decided intuitively by simple vector analyzing and just a few experimental verifications. This work would be useful for optimizing spectral sensitivity of broadband multispectral imaging sensors or SFA sensors.

The water cushion depth of stilling basin with shallow-water cushion is a key factor that affects the flow regime of hydraulic jump in the basin. However, the specific depth at which the water cushion is considered as “shallow” has not be stated clearly for now, and only conceptual description is provided. This paper attempts to specify the best water cushion depth based on the flow regime of hydraulic jump and underflow speed; namely, in case of critical hydraulic jump in the basin, the best water cushion depth is located where the minimum distance to the bottom plate of the stilling basin is 1/5~1/4 of the water cushion depth. The theoretical analysis indicates, at different inclinations of discharge chute (θ) and depth ratios of inlet (m), instead of monotonic change, the Froude number (Fr) at inlet of the stilling basin with shallow-water cushion firstly reduces and then increases as the flow velocity at discharge chute inlet (V) increases; the parameters of inflection point (critical flow velocity and critical Fr) increase as the inclinations of discharge chute (θ) and depth ratios of inlet (m) increase. Such regularity is the theoretical basis for selecting representative study cases. The reliability of the large eddy simulation calculation results are verified by a model test; in the paper, 30 cases including five different Froude numbers and six shallow-water cushion depths are selected, for calculating the hydraulic factors such as flow profile, flow regime and flow velocity in the stilling basin with shallow-water cushion; and the varying pattern between the best depth of stilling basin with shallow-water cushion (depth-to-length ratio) and the inflow Froude number is obtained which indicates that the best depth of stilling basin with shallow-water cushion varies little as the change of the Froude number before reaching the critical Froude number; however, the best depth-to-length ratio of stilling basin with shallow-water cushion increases as the Froude number increases after the critical Froude number is reached. The study results in this paper are of reference significance to design and calculation of the stilling basin with shallow-water cushion.

In this paper we introduce a computational model for the simulation of hypervelocity impact (HVI) phenomena which is based on the Discrete Element Method (DEM). Our paper constitutes the first application of DEM to the modeling and simulating of impact events for velocities beyond 5 kms⁻¹. We present here the results of a systematic numerical study on HVI of solids. For modeling the solids, we use discrete spherical particles that interact with each other via potentials. In our numerical investigations we are particularly interested in the dynamics of material fragmentation upon impact. We model a typical HVI experiment configuration where a sphere strikes a thin plate and investigate the properties of the resulting debris cloud. We provide a quantitative computational analysis of the resulting debris cloud caused by impact and a comprehensive parameter study by varying key parameters of our model. We compare our findings from the simulations with recent HVI experiments performed at our institute. Our findings are that the DEM method leads to very stable, energy–conserving simulations of HVI scenarios that map the experimental setup where a sphere strikes a thin plate at hypervelocity speed. Our chosen interaction model works particularly well in the velocity range where the local stresses caused by impact shock waves markedly exceed the ultimate material strength.

The results of numerical experiments with the chemistry-climate model of the Institute of Numerical Mathematics of the Russian Academy of Sciences - Russian State Hydro-Meteorological University (INM-RSHU CCM) in comparison with the MERRA2 reanalysis data are presented. To assess the impact of sea surface temperature (SST) on stratospheric processes and the ozone layer, numerical experiments were carried out for scenarios corresponding to the El Niño and La Niña Southern Oscillation phases. It is shown that the El Niño phase contributes to the enhancement of heat and mass transfer to the polar stratosphere, which strengthen the Brewer-Dobson circulation and the flux of wave activity that affect the zonal wind. As a result, during El Niño, the probability of sudden stratospheric warming increases and the stratospheric polar vortex (SPV) weakens, re-sulting in an increase in total polar stratospheric ozone. The La Niña phase, on the contrary, is accompanied by a weakening of heat and mass transfer from the troposphere to the stratosphere, as a result of which the stratospheric polar vortex strengthens and the total column ozone decreases. It is shown that the effect associated with the Southern Oscillation is more pronounced in the Northern Hemisphere than in the Southern Hemisphere. Numerical experiments were also carried out for scenarios in which the sea surface temperature was set corresponding to the early 1980s, and the results were compared with calculations in which the SST was set corresponding to the late 2010s. This allowed us to assess the influence of the multiyear trend in the SST increase on stra-tospheric processes and the ozone layer. It was shown that the SST increase is most pronounced in the Arctic, which corresponds to the phenomenon of "Arctic amplification", and therefore, the heat and mass transfer to the stratosphere increases in the Northern Hemisphere, which contributes to some weakening of the SPV and increase of the total column ozone. But, compared to the Southern Oscillation, the influence of the SST trend on the SPV and total column ozone is much weaker and is observed mainly in the Northern Hemisphere.

This work is devoted to consideration and analysis of the application of molecular dynamics simulation (MDS) methods to the study of nanosized polymer polyvinylidene fluoride (PVDF) thin ferroelectric films (two-dimensional ferroelectrics) and their composites with graphene layers: 1) to study and calculations of the polarization switching time depending on the electric field and PVDF film thickness; 2) to study and calculations of the polarization switching time depending changes of the PVDF on PVDF-TrFE film; 3) to study the polarization switching time in PVDF under influence of graphene layers. All calculations at each MDS step were carried out using quantum semi-empirical methods PM3. A comparison and analysis of the results of these calculations and the kinetics of polarization switching within the framework of the Landau-Ginzburg-Devonshire theory for homogeneous switching in ferroelectric polymer films is carried out. The study of the composite heterostructures of the “graphene-PVDF” type and calculations of their polarization switching times were presented too. It is shown that replacing PVDF with PVDF-TrFE significantly changes the polarization switching times in these thin polymer films, and that introducing various graphene layers into the PVDF layered structure leads to both an increase and a decrease in the polarization switching time. Here everything also depends on the position and displacement of the ferroelectric coercive field depending on the system damping parameters.

In this paper, we study the effect of parasitic channel on electrical characteristics and power consumptions of vertically stacked gate-all-around (GAA) silicon nanosheet (NS) MOSFETs of the bulk and silicon-on-insulator (SOI) substrates. Through deliberate control of the bottom parasitic channel coverage, we systematically assess its influence on DC, AC, and RF characteristics including power consumption. The observed findings highlight a significant differentiation among the bottom parasitic channels: GAA (i.e., control sample with 100% coverage), omega-shaped gate NS (70% and 80% coverages), and FinFET (60% coverage). Compared with the device with a 60% coverage channel, the device with a 100% coverage channel possesses 3520% and 8030% normalized leakage current improvement for both NFET and PFET under the bulk substrate, respectively. Compared with bulk-substrate devices, the SOI-substrate devices further significantly improve the leakage current. Due to enhanced off-state current leakages in the 100% covered bulk-substrate device, the static power demonstrates a notable improvement (about 98.1%), compared to 60% covered devices. Hence, the utilization of a 60% covered channel will result in a degradation of both device and circuit performance. The observation of GAA NS MOSFET indicates that an increase in the bottom parasitic channel coverage ratio leads to a positive impact on power consumption and leakage currents. The examination presented here is useful in the fabrication of stacked GAA NS MOSFETs.

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.

A quasi-two-dimensional (quasi-2D) modelling approach is introduced to mimic transverse mixing of an inflow into a river from one of its banks, either an industrial outfall or a tributary. The concentrations of determinands in the inflow vary greatly from those in the river, leading to very long mixing lengths in the river downstream of the inflow location. Ideally, a two-dimensional (2D) model would be used on a small scale to capture the mixing of the two flow streams. However, for large-scale applications of several hundreds of kilometres of river length, such an approach demands too many computational resources and too much computational time, especially if the application will at some point require ensemble input from climate-change scenario data. However, a one-dimensional (1D) model with variables varying in the longitudinal flow direction but averaged across the cross-sections is too simple of an approach to capture the lateral mixing between different flow streams within the river. Hence, a quasi-2D method is proposed in which a simplified 1D solver is still applied but the discretisation of the model setup can be carried out in such a way as to enable a 2D representation of the model domain. The quasi-2D model setup also allows secondary channels and side lakes in floodplains to be incorporated into the discretisation. To show proof-of-concept, the approach has been tested on a stretch of the lower Athabasca River in Canada flowing through the oil sands region between Fort McMurray and Fort MacKay. A dye tracer and suspended sediments are the constituents modelled in this test case.

Globally, groundwater is the largest distributed storage of freshwater that plays an important role in an ecosystem’s sustainability in addition to aiding human adaptation to both climatic change and variability. However, this resource is not unlimited and its sustainability is highly dependent on its prudent use. Thus, efficient management of groundwater resources to prevent overexploitation, scarcity and drought has become a major challenge for researchers as well as water managers. To solve these challenges, many solutions such as simulation and optimisation models have been proffered through the use of historical data. Therefore, this has made efficient data gathering essential to maintain data-driven groundwater level resource management models from the observation site. The global evolution of the Internet of Things (IoTs), has increased the nature of data gathering for the management of groundwater resources. Recently, a number of research challenges such as the lack of computational efficiency and scalability due to uncertainties from input parameters to the groundwater level resource model have been revealed in the management of groundwater level resources. In addition, efficient data-driven groundwater level resource management relies hugely on information relating to changes in groundwater resource levels as well as its availability. At the moment, the groundwater managers lack an efficient and scalable groundwater management system to gather the required data. The literature revealed that the existing methods of collecting data lack efficiency to meet computational model requirements and meet management objectives. Although the IoTs enabled automated data processing systems are in existence by transmitting the generated data from IoT enabled devices into the cloud through the Internet. However, traditional IoT Internet is not scalable and efficient enough to process the generated vast IoT data. Thus, it is necessary to have an efficient and scalable IoT system to extract valuable information in real-time for groundwater level resource management. Unlike previous surveys which solely focussed on particular groundwater issues related to simulation and optimisation management models, nonetheless, this paper seeks to highlight the current groundwater level resources management models as well as the IoT contributions.

Current industrial trends bring new challenges in energy absorbing systems. Polymer materials as the traditional packaging material seem to be promising due to their low weight, structure and production price. Based on the review, the linear low-density polyethylene material was identified as the most promising material for absorbing impact energy. The current paper addresses the identification of the material parameters and the development of a Constitutive material model to be used in future design by virtual prototyping. The paper deals with the experimental measurement of the stress-strain relations of the linear low-density polyethylene under static and dynamic loading. The quasi-static measurement is realized in two perpendicular principal directions and is supplemented by a test measurement in the 45 degrees direction, i.e. exactly between the principal directions. The quasi-static stress-strain curves are analyzed as an initial step for dynamic strain rate dependent material behavior. The dynamic response is tested in the drop tower using a spherical impactor hitting the flat material multi-layered specimen at two different energy levels. The strain rate dependent material model is identified by optimizing the static material response obtained in the dynamic experiments. The material model is validated by the virtual reconstruction of the experiments and by comparing the numerical results to the experimental ones.

In order to realize the tram's low-floor structure, most of the trams that have been recently introduced adopt an independently rotating wheelset. In the case of trains driving in two regions with different gauges, an independently rotating wheelset may be applied to utilize the variable track technology. Since the independent rotation type wheelset has no rotational restraint of the left and right wheels, the difference in rotational speed between the outer and inner wheels occurs naturally during curved driving, and it is applied to railroad vehicles traveling in steep curve sections because it smoothly drives curved driving. However, the longitudinal creep force and the lateral restoring force are weakened as the left and right rotational constraints disappear. Lack of transverse direction restoring force weakens stability while causing continuous flange contact driving or zigzag phenomenon against disturbance. Under the conditions of driving in a steep curve, these railway vehicles generate excessive wear, noise, and lateral pressure, as well as deterioration of ride comfort and derailment. In order to overcome these drawbacks, a method has been proposed in which the torque of a motor mounted on each wheel is individually controlled to generate lateral restoring force or to improve driving performance through lateral displacement control using a yaw moment. In this paper, development using HILs was performed to check the performance and stability of the individual motor torque control technology before verifying by applying the individual motor torque control to the actual vehicle. HILs were constructed by combining a real-time dynamic analysis model of a railway vehicle with a drive motor to which real individual motor control was applied. Under the conditions of driving the test track where the actual test vehicle was tested, the analysis of the driving characteristics and the control characteristics of the disturbance was performed to confirm the proposed individual motor torque control performance.

The application of thermoelectric generators (TEG) on the nozzle of an aviation jet engine was studied by finite element TEG-simulations. Against the background of system-level requirements of the reference aircraft this work reports the resulting requirements on the TEG design with respect to applied thermoelectric (TE) element lengths and fill factors (F) within the TE modules in order to maintain a positive effect on the specific fuel consumption. Assuming a virtual optimized TE material and varying the convective heat transfer coefficients at the nozzle surfaces this work reports the achievable power output. System-level requirement on the gravimetric power density (> 100 Wkg-1) can only be met for F ≤ 21%. Extrapolating TEG coverage to the full nozzle surface, the power output reaches 1.65 kW per engine. Assessment of further potential is demonstrated by a parametric study on the fill factor, convective heat transfer coefficients, and materials performance. This study confirms a feasible design range for TEG installation on the aircraft nozzle with a positive impact on the fuel consumption. This application translates into a reduction of operational costs, allowing for an economically efficient installation of TEG in consideration of the cost-specific power output of modern thermoelectric materials.

Controlling functional properties of matter and combine them for engineering a functional device is nowadays a common direction of scientific community. For instance, heterogeneous magnetic nanostructures can make use of different types of geometrical and compositional modulations to achieve the control of the magnetization reversal along with the nano-entities and thus enabling the fabrication of spintronic, magnetic data storage and sensing devices, among others. In this work, diameter modulated FeNi nanowires are fabricated paying special effort to obtain sharp transition regions between two segments of different diameters (from about 450 nm to 120 nm), enabling precise control over the magnetic behavior of the sample. Micromagnetic simulations performed on single bi-segmented nanowires predict a double step magnetization reversal where the wide segment magnetization switches near 200 Oe through a vortex domain wall, while at 500 Oe the magnetization of the narrow one is reversed through a corkscrew like mechanism. Finally, these results are confirmed with magneto-optic Kerr effect measurements at the transition of isolated bi-segmented nanowires. Furthermore, macroscopic vibrating sample magnetometry is used to demonstrate that the magnetic decoupling of nanowire segments is the main phenomenon occurring over the entire fabricated nanowires.

We present observations from evaluation of internal environmental quality of industrial halls with priority on daylighting in combination with the integral lighting. The physical parameters related to indoor lighting in large industrial halls in winter and summer periods were analyzed using in situ measurements and computational methods. These are part of a comprehensive research on indoor environmental quality of industrial halls with the aims of saving energy and providing a comfortable environment for the workers while improving the productivity. The results showed that the procedures used for evaluation of residential or office buildings may not be used for industrial buildings. We also observed that the criteria of occupants’ comforts for indoor industrial buildings may differ from those of other kinds of buildings. Based on these results, an adequate attention is required for designing the industrial buildings. For this reason, appropriate evaluation methods and criteria should be created. We found the measured values of daylight factor very close to the skylight component of the total illumination. The skylight component was observed on average 30% that of the measured daylight factor values. Although the daylight is not emphasized when designing the industrial buildings and its contribution is small, but it is very important for the workers psychology and physiology. The workers must feel a connection with the exterior environment; otherwise, their productivity decreases.

This paper offers an integrated approach for correlated storage assignment strategy (CSAS) to optimize travel distance by considering item correlation, picking frequency, zoning, and slot constraints to improve operational efficiency in distribution centers. We simplify the mixed integer non-linear programming (MINLP) model and incorporate integrated heuristic procedures for faster convergence. We introduce positive and negative centroid deviations as techniques to guide the model convergence and explore different scenarios. In the second stage of item assignment, we prioritize items within zones using a ranking formula that is optimal for both single-item and multiple-item orders. By analyzing the distribution of correlated item across zones and the impact of their picking frequency on optimizing travel distance, we propose an improved CSAS that minimize travel distance by strategically placing items based on their proximity to the depot. Our model reduces trips in distant zones by prioritizing larger average order sizes, while maximizing trips for frequently ordered smaller lists in closer zones. This method significantly reduces travel distance by minimizing the need for pickers in distant zones to frequently travel to the depot. Simulation results show reductions of up to 36.75% for fitted data and 33.59% for recent data.

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

Based on the engineering background of Daliuta mining area, the fracture distribution and water movement within the weakly cemented overlay strata are studied through liquid-solid coupling physical simulation and COMSOL numerical simulation test. 3D fracture network of overburden strata in mining is plotted based on the result of the physical simulation test. The results showed that the ‘fractured zone’ of the weakly cemented overburden strata in the western area develops more evidently than that in the eastern area ---- wide-opened fractures form at both ends of the working area and this leads to a wide distribution of fractures in the overlying strata.After the initial caving of the roof, height of the penetrating fractures increases rapidly and the separated strata space becomes the main channel for water flow. Under the action of water flow and rock strata movement, mining-induced fractures in overlying rock display cyclic change of expansion, penetration and closure. When the working face is fully mined, the penetrating fractures in the overlying strata and mining-induced fractures up in the working face become the main passages for water flow. The experimental result provides a reference for underground reservoirs layout and maintenance and water-conserved mining in western China.

Statistical modeling of lifetime data plays an essential role in a wide range of practical fields, such as health and engineering. There have been a lot of studies done to develop statistical models that can better describe health data than traditional models. For the first time, we pioneer a novel family of continuous probability distributions called the generalized odd beta prime generalized (GOBP-G) family of distributions. The cumulative distribution and probability density functions of the new family are presented. A new generalization of the Weibull distribution called "generalized odd beta prime-Weibull" (GOBPW) is proposed using the pioneered GOBP-G family. The mixture representations of the new distribution are defined and derived. Some formal statistical properties of the GOBPW distribution, such as the moments, moment generating function, incomplete moments, information generating function, entropies, stress-strength function, quantile function, and order statistics, are derived. The estimation of the parameters of the proposed distribution is evaluated using the maximum likelihood estimation approach. Different cancer disease data sets, such as the bladder, head and neck, acute bone, and blood cancers, are used to illustrate the applicability and usefulness of the new model and were compared using several statistical accuracy measures with that of well-established extended Weibull distributions, which are the beta modified Weibull distribution, Kumaraswamy modified Weibull distribution, gamma generalized modified Weibull distribution, gamma log-logistic Weibull distribution, and beta log-logistic Weibull distribution. The results show that the proposed model gives better results than the competitive models. This study could guide the relevant stakeholders in choosing a suitable statistical model for the health data instead of relying on traditional models to enhance decision-making.

It is not easy to dynamically analyze thermal bridges that require multidimensional analysis in building energy simulations, which are mostly one-dimensional platforms. To solve this problem, many studies have been conducted and, recently, a study was conducted to model the thermal bridge based on the data by approaching this in a similar way to steady-state analysis, showing high accuracy. This was an early-stage study, which is only applicable when the indoor temperature is constant. By extending this study, a thermal bridge model that can be applied even when the indoor temperature changes over time is proposed and validated. Since the governing equation, the heat diffusion equation, is linear, the key idea is to create and apply two thermal bridge transfer function models by expressing the heat flow entering the room as a linear combination of the transfer function for indoor temperature and the transfer function for outdoor temperature. For the proposed thermal bridge model, the NRMSE of the model itself showed a high accuracy of 99.9%, and in the verification through annual simulation using the model, the NRMSE showed an accuracy of 88.8%.

Microturbines can be used not only in models and education but also to propel UAVs. However, their wider adoption is limited by their relatively low efficiency and durability. Validated simulation models are required to monitor their performance, improve their lifetime, and design engine control systems. This study aims at developing a numerical model of a micro gas turbine for prediction and prognostics of engine performance. To build a reliable zero-dimensional model, the available compressor and turbine maps were scaled to the available test bench data with the least squares method, to meet the performance of the engine achieved during bench and flight tests. A steady-state aeroengine model was implemented in GSP and compared with experimental operating points. The selected flight data was then used as input for the transient engine model. The exhaust gas temperature (EGT) and fuel flow were chosen as the two key parameters to validate the model, comparing the numerical predicted values with the experimental ones. The observed difference between the model and the flight data was lower than 3% for both EGT and fuel flow.

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

The present work deals with the development of finite element methodology for obtaining the stress distributions in thick cylindrical HK40 stainless steel pipe that carry high temperature fluids. The material properties and loading are assumed to be random variables. Thermal stresses that are generated along radial, axial and tangential directions are computed generally using analytical expressions which are very complex. To circumvent such an issue, the probability theory and mathematical statistics have been applied to many engineering problems which allows to determine the safety both quantitatively and objectively based on the concepts of reliability. Monte Carlo simulation methodology is used to study the probabilistic characteristics of thermal stresses which is used for estimating the probabilistic distributions of stresses against the variations arising due to material properties and load. A 2-D Probabilistic finite element code is developed in MATLAB and the deterministic solution is compared with ABAQUS solutions. The values of stresses that are obtained from the variation of elastic modulus are found to be low as compared to the case where the load alone is varying. The probability of failure of the pipe structure is predicted against the variations in internal pressure and thermal gradient. These finite element framework developments are useful for the life estimation of piping structures in high temperature applications and subsequently quantifying the uncertainties in loading and material properties.

Combining wind and solar photovoltaic (PV) generation can provide complementary renewable power production, but depends on correlated resources. This study analyzed 10 years of wind data from Naama, Algeria to evaluate the potential for evening wind generation to offset the loss of solar at sunset. Average wind speeds showed a distinct increase during evening hours, coinciding with the decrease in solar irradiance. Wind turbine simulations using a 1.5 MW turbine and the wind data showed sufficient resources for profitable power production after sunset. Statistical analyses confirmed significantly higher wind speeds and simulated power output in evening vs daylight periods (p<0.05). The Pearson correlation coefficient between evening wind speeds and decreasing solar irradiance was 0.63, supporting a strong positive relationship. These findings indicate Naama has adequate wind resources to deploy economically viable wind power capacity that can complement existing solar infrastructure and provide renewable electricity after dark , .

Neural or mental simulation of actions is a powerful tool for allowing cognitive agents to develop Prospection Capabilities that are crucial for learning and memorizing key aspects for challenging skills. In previous studies we developed an approach based on the animation of the redundant human Body Schema based on the Passive Motion Paradigm (PMP). In this paper we show that this approach can be easily extended to hyper-redundant serpentine robots as well as to hybrid configurations where the serpentine robot is functionally integrated with a traditional skeletal infrastructure. A simulation model is analyzed in detail showing that it incorporates spatio-temporal features discovered in the biomechanical studies of biological hydrostats like the elephant trunk or the octopus tentacles. It is proposed that such generative internal model can be the basis for a cognitive architecture appropriate for serpentine robots, independent of the underlying design and control technologies. Although robotic hydrostats have received a lot of attention in the last decades, the great majority of research activities have been focused on the actuation/sensorial/material technologies that can support the design of hyperredundant soft robots as well as the related control methodologies. The cognitive level of analysis has been limited to motion planning, without addressing synergy formation and metal time travel. This is where the contribution of the paper is focused on.

Space-based gravitational wave detection is extremely sensitive to disturbances. The Keplerian configuration cannot accurately reflect the variations in spacecraft configuration. Planetary gravitational disturbances are one of the main sources. Numerical simulation is an effective method to investigate the impact of perturbation on spacecraft orbits. This study shows that, in the context of the Taiji project, Earth’s gravity is an essential factor in the change in heliocentric formation configuration, contributing to the relative acceleration between spacecrafts in the order of O(10−6)m·s−2. Considering 00:00:00 on October 27, 2032 as the initial orbiting moment, under the influence of Earth’s gravitational perturbation, the maximum relative change in armlengths and variation rates of armlengths for Taiji is 1.6×105km, 32m·s−1 respectively, compared with the unperturbed Keplerian orbit. Additionally, by considering the gravitational perturbations of Venus and Jupiter, the armlength and relative velocity for Taiji are reduced by 16.01% and 17.45%, respectively, compared with when only considering that of Earth. The maximum amplitude of the formation motion indicator changes with the orbit entry time. Results show that the relative velocity increase between the spacecrafts is minimal when the initial orbital moment occurs in July. Moreover, the numerical simulation results are inconsistent when using different ephemerises. The differences between ephemeris DE440 and DE430 are smaller than those between DE440 and DE421.

Non-Linear Simulation by Harmonic Balance Techniques of Load Modulated Power Amplifier Driven by Random Modulated Signals. Comparison with measured results.

This paper investigates urban microclimate in relation to the typical major streets and mobility layout in Turin, i.e. boulevards with lateral minor roads, called “controviali”, an almost unique urban mobility layout. The aim is to examine the possibility of intervention in relation to the type and scope of mobility for microclimate mitigation and climate adaptation purposes. To this, this paper evaluates to which extent urban morphology and mobility features contribute to determining microclimatic conditions, to finally identify site-specific prioritised strategies to cope with its mitigation. This research is structured in two phases. The first phase is dedicated to urban-scale analysis of the streets with “controviali”, resulting in the creation of a database containing their properties and typology, to finally derive six clusters of boulevards. The second phase involves the selection of one representative case study per cluster and simulation of their microclimate conditions, performed by the software ENVI-met. The results are discussed in relation to the parameters proven in the prior literature review to have a major impact on the microclimate. The assessment of the results indicated which elements affect the climatic conditions the most, allowing for the proposal of different prioritised mitigation strategies in the conclusion.

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 difficult current environmental situation, caused by construction industry residues containing ceramic materials could be improved by using these materials as recycled aggregates in mortars, with their processing causing a reduction in their use in landfill, contributing to recycling and also minimizing the consumption of virgin materials. Although some research is currently being carried out into recycled mortars, little is known about their stress-strain (σ-ε); therefore this work will provide the experimental results obtained from recycled mortars with recycled ceramic aggregates (with contents of 0, 10, 20, 30, 50 and 100%), such as: the density, the compression strength, as well as the σ-ε curves representative of their behavior. The values obtained from the analysis process of the results are those of: σ (elastic ranges and failure maximum), ε (elastic ranges and failure maximum), and Resilience and Toughness; in order to finally obtain, through numerical analysis, the equations to predict their behavior (related to their recycled content). At the end of the investigation it is established that mortars with recycled ceramic aggregate contents of up to 20% could be assimilated just like mortars with the usual aggregates, and the prediction equations produced could be used in cases of similar applications.

We propose a non-invasive approach to control the quality of spheroids and their fusion into a complex bioconstruct. The proposed method is based on the union of the nutrient concentration change calculated using Fick's law and the "reaction-diffusion" equations, taking into account absorption coefficient with specifying the exact fusion geometry using implicit Function Repre-sentation (FRep) functions. The proposed approach allows us to analyze the viability of cells within the spheroid, predict the fusion of spheroids, and accurately model complex heteroge-neous biostructures as a future task for our research. These results will significantly accelerate the development of such a promising field of additive biotechnologies as 3D bioprinting.

Summary: Studying biological systems generally relies on computational modeling and simulation, e.g., for model-driven discovery and hypothesis testing. Progress in standardization efforts led to the development of interrelated file formats to exchange and reuse models in systems biology, such as SBML, the Simulation Experiment Description Markup Language (SED-ML), or the Open Modeling EXchange format (OMEX). Conducting simulation experiments based on these formats requires efficient and reusable implementations to make them accessible to the broader scientific community and to ensure the reproducibility of the results. The Systems Biology Simulation Core Library (SBSCL) provides interpreters and solvers for these standards as a versatile open-source API in Java™. The library simulates even complex bio-models and supports deterministic Ordinary Differential Equations (ODEs); Stochastic Differential Equations (SDEs); constraint-based analyses; recent SBML and SED-ML versions; exchange of results, and visualization of in silico experiments; open modeling exchange formats (COMBINE archives); hierarchically structured models; and compatibility with standard testing systems, including the Systems Biology Test Suite and published models from the BioModels and BiGG databases. Availability and implementation: SBSCL is freely available at https://draeger-lab.github.io/SBSCL/.

The famous singlet correlations of a composite quantum system consisting of two spatially separated components exhibit notable features of two kinds. The first kind consists of striking certainty relations: perfect correlation and perfect anti-correlation in certain settings. The second kind consists of a number of symmetries, in particular, invariance under rotation, as well as invariance under exchange of components, parity, or chirality. In this note, I investigate the class of correlation functions that can be generated by classical composite physical systems when we restrict attention to systems which reproduce the certainty relations exactly, and for which the rotational invariance of the correlation function is the manifestation of rotational invariance of the underlying classical physics. I call such correlation functions classical EPR-B correlations. It turns out that the other three (binary) symmetries can then be obtained "for free": they are exhibited by the correlation function, and can be imposed on the underlying physics by adding an underlying randomisation level. We end up with a simple probabilistic description of all possible classical EPR-B correlations in terms of a "spinning coloured disk" model, and a research programme: describe these functions in a concise analytic way.

Nuclear technologies have a strong potential and a unique role to play in delivering reliable low carbon energy for a future net-zero society. However, to assure the sustainability required for the long-term success, nuclear will need to deliver innovative solutions as proposed in iMAGINE. One of the most attractive features, but also a key challenge for the envisaged highly integrated nuclear energy system is the need for a demand driven salt clean-up system. The work described provides an insight into the interplay between a potential salt clean-up system and the reactor operation in a plutonium started core in a dynamic approach. The results presented will help to optimize the parameters for the salt clean-up process as well as to understand the differences which appear between a core started with enriched Uranium and Plutonium as the fissile material. The integrated model is used to investigate the effects of the initial fissile material on core size, achievable burnup, and long term operation. Different approaches are tested to achieve a higher burnup in the significantly smaller Pu driven core. The effects of different clean-up system throughputs on the concentration of fission products in the reactor salt and its consequences are discussed for general molten salt reactor design. Finally, an investigation of how a plutonium loaded core could be used to provide fuel for future reactors through fuel salt splitting is presented with the outcome that one Pu started reactor of the same size as a uranium started core could deliver fuel for 1.5 new cores due to enhanced breeding. The results provide an essential understanding for the progress of iMAGINE as well as the basis for inter-disciplinary work required for optimizing iMAGINE.

Aortic disease has a significant impact on quality of life. Involvement of the aortic arch requires preservation of blood supply to the brain during surgery. Deep hypothermic circulatory arrest is an established technique for this purpose although neurological injury remains high. Additional techniques have been used to reduce the risk although controversy still remains. A three-way cannulation approach, including both carotid arteries and the femoral artery or the ascending aorta, has been used successfully for aortic arch replacement and redo procedures. We have de-veloped circuits of the circulation to simulate blood flow during this type of cannulation set up. The aim is to analyse using CARDIOSIM© cardiovascular simulation platform, how the haemodynamic and energetic parameters are affected and the benefit derived with particular reference to the cerebral circulation.

Electrifying off-grid and isolated islands in the Philippines remains one of the challenges that hinders community development. One of the solutions seen to ensure energy security, expand energy access and promote a low carbon future in this isolated islands is the use of renewable energy sources. This study wishes to determine the nearshore wave energy resource during monsoon seasons in Cuyo Island using a 40-year wave hindcast and 9-year on-site wind speed data to develop high resolution wave energy model using SWAN wave model, and assessed its annual energy production through matching with wave energy devices. Results shows that average significant wave height (Hs), peak period (Tp) and wave power density (Pd) during northeast monsoon are Hs = 1.35 m, Tp = 4.79 s and Pd = 4.05 kW/m respectively, while southwest monsoon which is sheltered by the mainland resulted to a lesser outcome, Hs = 0.52 m, Tp = 3.37 s and Pd = 0.34 kW/m. While the simulated model was observed to overestimate the wave energy resource (Bias = 0.398, RMSE = 0.54 and SI = 1.34), it has a strong relationship with the observed values (average r = 0.9). Its annual energy production is highest at Station 5, with AEPWaveBouy = 43.761 MWh, AEP-Pelamis = 216.786 MWh and AEPWave Dragon = 2462.66 MWh. At present, the minimum requirement for a wave energy development to be feasible is 5 kW/m, which in this case, Cuyo Island falls short, but with the continuous evolution of wave energy converters, applications on milder re-sources will soon materialized.

In this review, we trace the evolution of the quantum spin-wave theory treating non-collinear spin configurations. Non-collinear spin configurations are consequences of the frustration created by competing interactions. They include simple chiral magnets due to competing nearest-neighbor (NN) and next-NN interactions and systems with geometry frustration such as the triangular antiferromagnet and the Kagomé lattice. We review here spin-wave results of such systems and also systems with the Dzyaloshinskii-Moriya interaction. Accent is put on these non-collinear ground states which have to be calculated before applying any spin-wave theory to determine the spectrum of the elementary excitations from the ground states. We mostly show results from a self-consistent Green’function theory to calculate the spin-wave spectrum and the layer magnetizations at finite T in two and three dimensions as well as in thin films with surface effects. Some new unpublished results are also included. Analytical details and the validity of the method are shown and discussed.

There is an increased interest in the district-scale energy transition within interdisciplinary research community. Agent-based modelling presents a suitable approach to address variety of questions related to policies, technologies, processes, and the different stakeholder roles that can foster such transition. This state-of-the-art review focuses on the application of agent-based modelling for exploring policy interventions that facilitate the decarbonisation (i.e., energy transition) of districts and neighbourhoods while considering stakeholders’ social characteristics and interactions. We systematically select and analyse peer-reviewed literature and discuss the key modelling aspects, such as model purpose, agents and decision-making logic, spatial and temporal aspects, and empirical grounding. The analysis reveals that the most established agent-based models’ focus on innovation diffusion (e.g., adoption of solar panels) and dissemination of energy-saving behaviour among a group of buildings in urban areas. We see a considerable gap in exploring the decisions and interactions of agents other than residential households, such as commercial and even industrial energy consumers (and prosumers). Moreover, measures such as building retrofits and conversion to district energy systems involve many stakeholders and complex interactions between them that up to now have hardly been represented in the agent-based modelling environment.

Computation of Common Middle Point seismic sections and their subsequent time migration and diffraction imaging provides very important knowledge about the internal structure of 3D heterogeneous geological media and are key elements for successive geological interpretation. Full-scale numerical simulation, that computes all single shot seismograms, provides a full understanding of how the features of the image reflect the properties of the subsurface prototype. Unfortunately, this kind of simulations of 3D seismic surveys for realistic geological media needs huge computer resources, especially for simulation of seismic waves’ propagation through multiscale media like cavernous fractured reservoirs. In order to significantly reduce the query of computer resources we propose to model these 3D seismic cubes directly rather than the shot-by-shot simulation with subsequent CMP stacking. To do that we modify the well known "exploding reflectors principle" for 3D heterogeneous multiscale media by use of the finite-difference technique on the base of grids locally refined in time and space. To be able to simulate realistic models and acquisition we develop scalable parallel software, which needs reasonable computational costs. Numerical results for simulation of Common Middle Points sections and their time migration are presented and discussed.

Full-order state-space models represent the starting point for the development of advanced control methods for wind turbine systems (WTSs). Regarding existing control-oriented WTS models, two research gaps must be noted: (i) There exists no full-order WTS model in form of one overall ordinary differential equation that considers all dynamical effects which significantly influence the electrical power output; (ii) all existing reduced-order WTS models are subject to rather arbitrary simplifications and are not validated against a full-order model. Therefore, in this paper, two full-order nonlinear state-space models (of 11th and 9th-order in the (a, b, c)- and (d, q)-reference frame, resp.) for variable-speed variable-pitch permanent magnet synchronous generator WTSs are derived. The full-order models cover all relevant dynamical effects with significant impact on the system’s power output, including the switching behavior of the power electronic devices. Based on the full-order models, by a step-by-step model reduction procedure, two reduced-order WTS models are deduced: A non-switching (averaging) 7th-order WTS model and a non-switching 3rd-order WTS model. Comparative simulation results reveal that all models capture the dominant system dynamics properly. The full-order models allow for a detailed analysis covering the high frequency oscillations in the instantaneous power output due to the switching in the power converters. The reduced-order models provide a time-averaged instantaneous power output (which still correctly reflects the energy produced by the WTS) and come with a drastically reduced complexity making those models appropriate for large-scale power grid controller design.

Multi-junction solar cells comprised of stacked III-V semiconductor junctions represent the highest-efficiency photovoltaic technology, with recent demonstrations exceeding 47% efficiency . Optimizing the design and thickness of each junction is critical for maximizing performance . This work utilizes Silvaco TCAD tools to systematically optimize a 5-junction cell based on AlInP, AlGaInP, AlGaInAs, GaInP, GaAs, InGaAs, and Ge similar to recent record cells . The junction thicknesses are varied using a predictive profiler to sample the parameter space . For each combination, the spectral absorption and I-V characteristics are simulated to determine the efficiency . Statistical analysis identifies the optimal thickness set that maximizes performance to 26% under 1 sun illumination .

Left ventricular (LV) longitudinal function is mechanically coupled to the elasticity of the ascending aorta (AA). The pathophysiologic link between stiff AA and reduced longitudinal strain and subsequent deterioration in longitudinal LV systolic function is likely relevant in heart failure with preserved ejection fraction (HFpEF). A proposed therapeutic effect of freeing the LV apex and allowing for LV inverse longitudinal shortening was studied in-silico utilizing Living Left Heart Human Model (Dassault Systémes Simulia Corporation). LV function was evaluated in a model with A) Elastic AA; B) Stiff AA; and C) Stiff AA with free LV apex. The cardiac model simulation demonstrated that freeing the apex caused inverse LV longitudinal shortening that can abolish the deleterious mechanical effect of a stiff AA on LV function. A stiff AA and impairment of LV longitudinal strain are common in patients HFpEF. The hypothesis-generating model strongly suggests that freeing the apex and inverse longitudinal shortening may improve LV function in HFpEF patients with a stiff AA.

Molecular dynamics simulations were performed to examine the impact of temperature (ranging from 250K to 350K) on the structural behavior of cellulose chains when dissolved in an aqueous cuprammonium hydroxide Cuam solvent system. By monitoring the root mean square deviation (RMSD) and radius of gyration (Rg) values at each temperature, we found that the geometric deformation of the cellulose chains changed significantly at T=300K. In contrast, Rg and RMSD values remained stable at temperatures in the range of 250K, 270K and the range of 330K, 350K. Calculating the interaction energy between the solvent and cellulose chains revealed that temperature significantly influences the cellulose dissolution process in the Cuam solvent system and that T=300K is an efficient temperature for its dissolution. The number of intra- and interchain hydrogen bonds was also calculated as a function of temperature, and the analysis confirmed that there is no breakage of these bonds at temperatures below 300 K (i,e: 250K, 270K) and above 310 K (I,e, : 330K and 350K) either between two native chains or between a native chain and another derivative.

Several factors impact the durability of concrete bridge decks, including traffic loads, fatigue, temperature changes, environmental stress, and maintenance activities. Detecting problems such as corrosion, delamination, or concrete degradation early on can lower maintenance costs. Nondestructive evaluation (NDE) techniques can detect these issues at early stages. Each NDE method, meanwhile, has limitations that reduce the accuracy of the assessment. In this study, multiple NDE technologies were combined with machine learning algorithms to improve the interpretation of half-cell potential (HCP) and electrical resistivity (ER) measurements. A parametric study was performed to analyze the influence of five parameters on HCP and ER measurements, such as degree of saturation, corrosion length, delamination depth, concrete cover, and moisture condition of delamination. The results were obtained through finite element simulations and used to build two machine learning algorithms, a classification algorithm and a regression algorithm, based on the Random Forest methodology. The algorithms were tested using data collected from a bridge deck in the BEAST® facility. Both machine learning algorithms were effective in improving the interpretation of ER and HCP measurements using data from multiple NDE technologies.

To enable the deployment of battery-electric vehicles (BEV) as passenger cars in the private transport sector, a suitable charging infrastructure is crucial. In this paper a methodology for efficient spatial distribution of charging infrastructure is evaluated by investigating a scenario with a market penetration of BEVs of 100 percent (around 1.3 million vehicles). It aims towards the development of various charging infrastructure scenarios - including public and private charging - which are suitable to cover the charging demand. Therefore, these scenarios are investigated in detail with focus on number of public charging points, their spatial distribution, the available charging power and the necessary capital costs. For the creation of those charging infrastructure scenarios a placement model is developed. It uses the data of a MATSim (Multi-Agent Transport Simulation) traffic simulation of the metropolitan area of Berlin to evaluate and optimize different distributions of charging infrastructure. The model uses a genetic algorithm and the principle of multi-objective optimization. The capital cost of the charging points and the mean detour car drivers must cover additionally are used as optimization criteria. Using these criteria should lead to cost efficient infrastructure solutions which provide high usability at the same time. The main advantage of the method selected is that multiple optimal solutions with different characteristics can be found and suitable solutions can be selected by using other criteria subsequently. The optimized charging infrastructure solutions show capital costs between 624 and 2950 million euro. Users must cover an additionally mean detour of 254m to 590m per charging process to reach an available charging point. According to the results a suitable ratio between charging points and vehicles is between 11:1 and 5:1. A share of fast charging infrastructure (>50kW) of less than ten percent seems to be sufficient, if it is situated at main traffic routes and highly frequented places.

A new attack against the Kirchhoff-Law-Johnson-Noise (KLJN) key distribution system is introduced. The attack is based on 1) Utilizing the dc-voltage-source - which we put at Alice’s end in our mathematical modeling of the attack-that could exist due to errors, imbalances, Electromagnetic interference and etc.2) On studying the number of samples per bit in the security key that the measured Alice/Bob voltages exceeds or falls below a threshold voltage, respectively. The threshold voltage is the average between dc voltages across low and high resistors- generated by a dc-voltage source. We count the number of samples the voltage at Bob’s end (containing both the noise and dc components) exceeds the threshold voltage and how many times it falls below the threshold to judge whether the resistor can be guessed as low or high for every cycle. Also for a pre-specified key-length we count the number of high resistance estimations per bit –at Bob’s end-according to the previous criterion and the non-successful estimations per bit to judge the final guessed resistor value at Alice’s end and Bob’s end, where if we have more bits with most of its measured samples are above the threshold voltage then we will predict Bob’s resistance as high resistance, otherwise we predict Bob’s resistance to be low resistance. The Simulation was conducted and the attack proved that it is successful unless the temperature increased dramatically to ranges more than a threshold temperature ~ Kelvin that increases when the number of samples per bit increases.

The local and global order in dense packings of linear, semi-flexible polymers of tangent hard spheres are studied by employing extensive Monte Carlo simulations at increasing volume fractions. Chain stiffness is controlled by a tunable harmonic potential for the bending angle whose intensity dictates the rigidity of the polymer backbone as a function of the bending constant and equilibrium angle. The studied angles range from acute to obtuse ones, reaching the limit of rod-like polymers. We analyze how packing density and chain stiffness affect the ability of chains to self-organize at the local and global levels. The former corresponds to crystallinity as quantified by the Characteristic Crystallographic Element (CCE) norm descriptor, while the latter is computed through the scalar orientational order parameter. In all cases, we identify the critical volume fraction for the phase transition and gauge the established crystal morphologies, developing a complete phase diagram as a function of packing density and equilibrium bending angle. A plethora of structures is obtained, ranging from random hexagonal closed packed morphologies of mixed character and almost perfect face centered cubic (FCC) and hexagonal close-packed (HCP) crystals at the level of monomers, to nematic mesophases, with prolate and oblate mesogens at the level of chains. For rod-like chains, hysteresis is observed between the establishment of long-range nematic order and crystallization, while for right-angle chains both transitions are synchronized. A comparison is also provided against the analogous packings of monomeric and fully flexible chains of hard spheres.

Computationally modelling a nuclear reactor startup core for a benchmark against the existing models is highly desirable for an independent assessment informing nuclear engineers and energy policymakers. This work presents a startup core model of the UK’s first Evolutionary Pressurised Water Reactor (EPR) based on Monte Carlo simulations of particle collisions using Serpent 2, a continuous-energy Monte Carlo reactor physics burnup code. Coupling between neutronics and thermal-hydraulic conditions with the fuel depletion is incorporated into the multi-dimensional branches, obtaining the thermal flux and fission rate (power) distributions radially and axially from the three dimensional (3D) single assembly level to a 3D full core. Shannon entropy is employed to characterise the convergence of the fission source distribution, with 3 billion neutron histories tracked by parallel computing. Source biasing is applied for the variance reduction. Benchmarking the proposed Monte Carlo 3D full-core model against the traditional deterministic transport computation suite used by the UK Office for Nuclear Regulation (ONR), a reasonably good agreement within statistics is demonstrated for the safety-related reactivity coefficients, which creates trust in the EPR safety report.

Information is the currency of life, but the origin of prebiotic information remains a mystery. We propose transitional pathways from the cosmic building blocks of life to the complex prebiotic organic chemistry that led to the origin of information systems. The prebiotic information system, specifically the genetic code, is segregated, linear, and digital and probably appeared during biogenesis four billion years ago. In the peptide/RNA world, lipid membranes randomly encapsulated amino acids, RNA, and protein molecules, drawn from the prebiotic soup, to initiate a molecular symbiosis inside the protocells. This endosymbiosis led to the hierarchical emergence of several requisite components of the translation machine: tRNAs, aaRS, mRNAs, and ribosomes. When assembled in the right order, the translation machine created biosynthetic polypeptides, a process that transferred information from mRNAs to proteins. This was the beginning of the prebiotic information age. The molecular attraction between tRNA and amino acids led to different stages of the translation machines and the genetic code. tRNA is an ancient molecule that designed and built mRNA for storing the information of its cognate amino acid. Each mRNA strand became the storage device for the genetic information that encoded the amino acid sequences in triplet nucleotides. As information appeared in the digital languages of the codon within mRNA, and the genetic code for protein synthesis evolved, the prebiotic chemistry then became more organized and directional. The origin of the genetic code is enigmatic; herein we propose an evolutionary explanation: the demand for a wide range of specific enzymes in the peptide/RNA world was the main selective pressure for the origin of information-directed protein synthesis. We review three main concepts on the origin and evolution of the genetic code: the stereochemical theory, the coevolution theory, and the adaptive theory. These three theories are compatible with our coevolution model of the translation machines and the genetic code. We suggest biosynthetic pathways as the origin of the specific translation machines which provided the framework for the origin of the genetic code. During translation, the genetic code developed in three stages coincident with the refinement of the translation machines: GNC code developed by the pre-tRNA/pre-aaRS /pre-mRNA machine, SNS code by the tRNA/aaRS/mRNA machine, and finally the universal genetic code by the tRNA/aaRS/mRNA/ribosome machine. Our hypothesis provides the logical and incremental steps for the origin of the programmed protein synthesis. In order to understand the prebiotic information system better, we converted letter codons into numerical codons in the Universal Genetic Code Table. We have developed a software called CATI (Codon-Amino Acid-Translator-Imitator) to translate randomly chosen numerical codons into corresponding amino acids and vice versa. This conversion has granted us insight into how the translation might have worked in the peptide/RNA world. There is great potential in the application of numerical codons to bioinformatics such as barcoding, DNA mining, or DNA fingerprinting. We constructed the likely biochemical pathways for the origin of translation and the genetic code using the Model-View-Controller (MVC) software framework, and the translation machinery step-by-step. Using AnyLogic software we were able to simulate and visualize the entire evolution of the translation machines and the genetic code. The results indicate that the emergence of the information age from the peptide/RNA world was a watershed event in the origin of life about four billion years ago.

Surgical navigation has evolved as a vital tool in maxillofacial surgery, offering precise and patient-specific data. This study explores the clinical applications and accuracy of intraoperative tool tracking in maxillofacial surgery. The research includes 37 patients with various pathologies who underwent surgery assisted by a surgical navigation system using tracked instruments. The study showcases two representative cases: one involving coronoid hyperplasia with mouth opening deficit and another addressing nose-orbital-ethmoidal-frontal ossifying fibroma resection. The results indicate that surgical navigation with tracked instruments provides high precision (<1.5 mm error), reduced surgical time, and less invasive approaches. The study highlights the potential for reproducible outcomes and increased safety, especially in complex cases. Despite some limitations, the synergy between surgical navigation and tracked instruments offers a promising approach in maxillofacial surgery, expanding its applications beyond current practices.

The Late Heavy Bombardment Period (4.1 to 3.8 billion years ago) of heightened impact cratering activity on young Earth is likely the driving force for the origin of life. During the Eoarchean, asteroids such as carbonaceous chondrites delivered the building blocks of life and water to early Earth. Asteroid collisions created innumerable hydrothermal crater lakes in the Eoarchean crust which inadvertently became the perfect cradle for prebiotic chemistry. These hydrothermal crater lakes were filled with cosmic water and the building blocks of life. forming a thick prebiotic soup. The unique combination of exogenous delivery of extraterrestrial building blocks of life, and the endogenous biosynthesis in hydrothermal impact crater lakes very likely gave rise to life. A new symbiotic model for the origin of life within the hydrothermal crater lakes is here proposed. In this scenario, life arose around four billion years ago through five hierarchical stages of increasing molecular complexity: cosmic, geologic, chemical, information, and biological. During the prebiotic synthesis, membranes first appeared in the hydrothermal crater lakes, followed by the simultaneous origin of RNA and protein molecules, creating the RNA/protein world. These proteins were noncoded protein enzymes that facilitated chemical reactions. RNA molecules formed in the hydrothermal crater basin by polymerization of the nucleotides on the montmorillonite mineral substrate. Similarly, the initial synthesis of abiotic protein enzymes was mediated by the condensation of amino acids on pyrite surfaces. The regular wet-dry cycles within the crater lakes assisted further concentration, condensation, and polymerization of the RNAs and proteins. Lipid membranes randomly encapsulated amino acids, RNA, and protein molecules from the prebiotic soup to initiate a molecular symbiosis inside the protocells, this led to the hierarchical emergence of several cell components. As the role of protein enzymes became essential for catalytic process in the RNA/protein world, Darwinian selection from noncoded to coded protein synthesis led to translation systems and the genetic code, heralding the information stage. In this stage, the biochemical pathways suggest the successive emergence of translation machineries such as tRNAs, aaRS, mRNAs, and of ribosomes for protein synthesis. The molecular attraction between tRNA and amino acid led to the emergence of translation machinery and the genetic code. tRNA is an ancient molecule that created mRNA for the purpose of storing amino acid information like a digital strip. Each mRNA strand became the storage device for genetic information that encoded the amino acid sequences in triplet nucleotides. As information became available in the digital languages of the codon within mRNA, biosynthesis became less random and more organized and directional. The original translation machinery was simpler before the emergence of the ribosome than that of today. We review three main concepts on the origin and evolution of the genetic code: the stereochemical theory, the coevolution theory, and adaptive theory. We believe that these three theories are not mutually exclusive, but are compatible with our coevolution model of translations machines and the genetic code. We suggest biosynthetic pathways as the origin of the translation machine that provided the framework for the origin of the genetic code. During translation, the genetic code developed in three stages coincident with the refinement of the translation machinery: GNC code with four codons and four amino acids during interactions of pre-tRNA/pre-aaRS /pre-mRNA, SNS code consisting of 16 codons and 10 amino acids appeared during the tRNA/aaRS/mRNA interaction, and finally the universal genetic code evolved with the emergence of the tRNA/aaRS/mRNA/ribosome machine. The universal code consists of 64 codons and 20 amino acids, with a redundancy that minimizes errors in translation. To address the question of the origin of the biological information system in the RNA/protein world, we converted letter codons into numerical codons in the Universal Genetic Code Table. We developed a software called CATI (Codon-Amino Acid-Translator-Imitator) to translate randomly chosen numerical codons into corresponding amino acids and vice versa, gaining insight into how translation might have worked in the RNA/protein world. We simulated the likely biochemical pathways for the origin of translation and the genetic code using the Model-View-Controller (MVC) software framework, and the translation machinery step-by-step. We used AnyLogic software to simulate and visualize the evolution of the translation machines and the genetic code. We conclude that the emergence of the information age from the RNA/protein world was a watershed event in the origin of life about four billion years ago.

Rapid developments in the field of underwater photogrammetry have given scientists1the ability to produce accurate 3-dimensional (3D) models which are now increasingly used in the representation and study of local areas of interest. This paper addresses the lack of systematic analysis of 3D reconstruction and navigation fusion strategies, as well as associated error evaluation of models produced at larger scales in GPS-denied environments using a monocular camera (often in deep-sea scenarios). Based on our prior work on automatic scale estimation of Structure from Motion (SfM)-based 3D models using laser scalers, an automatic scale accuracy framework is presented. The confidence level for each of the scale error estimates is independently assessed through the propagation of the uncertainties associated with image features and laser spot detections using a Monte Carlo simulation. The number of iterations used in the simulation was validated through the analysis of the final estimate behaviour. To facilitate the detection and uncertainty estimation of even greatly attenuated laser beams, an automatic laser spot detection method, mitigating the effects of scene texture, was developed, with the main novelty of estimating the uncertainties based on the recovered characteristic shapes of laser spots with radially decreasing1 intensities. The effects of four different reconstruction strategies resulting from the combinations of Incremental/GlobalSfM, and thea priori/a posterioriuse of navigation data were analyzed using two distinct survey scenarios captured during the SUBSAINTES 2017 cruise (doi: 10.17600/17001000). The study demonstrates that surveys with multiple overlaps of non-sequential images result in a nearly identical solution regardless of the strategy (SfM or navigation fusion), while surveys with weakly connected sequentially acquired images are prone to produce broad-scale deformation (doming effect) when navigation is not included in the optimization. Thus the scenarios with complex survey patterns substantially benefit from using multi-objective BA navigation fusion. In all cases, the errors in the models are inferior to 5%, with errors often being around 1%. The effects of combining data from multiple surveys were also evaluated. The introduction of additional vectors in the optimization of multi-survey problems successfully accounted for offset changes present in the underwater USBL-based navigation data and thus minimize the effect of contradicting navigation priors. Our results also illustrate the importance of collecting a multitude of evaluation data at different locations and moments during the survey.

Experiments to measure a single molecule/particle, i.e. an individual molecule/particle, at room temperature or under physiological conditions without immobilization, for example on a surface, or without significant hydrodynamic flow have so far failed. This failure has given impetus to the underlying theory of Brownian molecular motion towards its stochastic due to diffusion. Quantifying the thermodynamic jitter of molecules/particles inspires many and forms the theoretical basis of single-molecule/single-particle biophysics and biochemistry. For the first time, our simulation results for a live cell (cytoplasm) show that the tracks of individual single molecules are localized in Brownian motion, while there is fanning out in fractal diffusion (anomalous diffusion).

This article presents a comprehensive study on enhancing grid resilience through advanced forecasting and optimization techniques in the context of power outages. Power outages pose significant challenges to modern societies, affecting various sectors such as industries, households, and critical infrastructures. The research combines statistical analysis, machine learning algorithms, and optimization methods to address this issue to develop a holistic approach for predicting and mitigating power outage events. The proposed methodology involves the use of Monte Carlo simulations in MATLAB for future outage prediction, Long Short-Term Memory (LSTM) networks for forecasting solar irradiance and load profiles, and a hybrid LSTM-Particle Swarm Optimization (PSO) model to improve accuracy. Furthermore, the role of Battery State of Charge (SoC) in enhancing system resilience is explored. The study also assesses the techno-economic advantages of a grid-tied microgrid integrated with solar panels and batteries over conventional grid systems. The results highlight the potential of the proposed approach in strengthening grid resilience, reducing downtime, and fostering sustainable energy utilization.

The demand for improving the nuclear waste management has since long been identified as one of the major hurdles for widespread use of nuclear energy. Nuclear waste management, through partitioning and transmutation (P&T) has been researched since the 1990ies with partitioning being a pre-requisite for the process. Recently, an innovative approach of reactors directly operating on spent nuclear fuel, iMAGINE, has been proposed which could deliver on the aims of P&T as a side effect to more efficient nuclear energy production in the future. The results shown here confirm that it is possible to achieve the aims of P&T without prior partitioning, but it is certainly less efficient since the high concentration of MAs, required for efficient burning, are impossible to obtain in a short operational time. The proposed nuclear waste management approach will be a long-term effort when it is accomplished without partitioning/separation technologies. However, none of the analyses contradicts this effort – the key points are a) when the technology for treating the waste is possible and reliable, the time horizon will not be a major concern b) the waste management is now intrinsically linked with energy production instead of requiring dedicated costly facilities, delivering a promising economic basis c) the waste management is now associated with long-term energy production and massively improved resource utilization. In summary, the objectives of P&T are achievable without prior partitioning – this is a ground-breaking result. However, the approach will require continued nuclear energy production for a very long time. In case of a nuclear phase-out decision, a much more efficient transmutation scenario would have to be envisaged as already discussed several years ago for Germany.

The aim of this study is to obtain the Bayes estimators and the maximum likelihood estimators (MLEs) for the unknown parameters of the Rayleigh Weibull (RW) distribution based on progres-sive type-II censored samples. The approximate Bayes estimators are calculated using the idea of Lindley and Tierney-Kadane's approximation method under the squared-error loss function when the Bayes estimators are not handed in explicit forms. In this study, the approximate Bayes esti-mates are compared with the maximum likelihood estimates in the aspect of the estimated risks (ERs) using Monte Carlo simulation. In addition, the coverage probabilities of the parametric bootstrap estimates are calculated. Real lifetime data sets belonging to the cancer types as bladder cancer, head and neck cancer, and leukemia are used to illustrate the emprical results belonging to the approximate Bayes estimates, the maximum likelihood estimates, and the parametric bootstrap intervals.

This paper presents a new methodology to derive and analyze strategies for a fully decarbonized urban transport system which combines conceptual vehicle design, a large-scale agent-based transport simulation, operational cost analysis, and life cycle assessment for a complete urban region. The holistic approach evaluates technical feasibility, system cost, energy demand, transportation time and sustainability-related impacts of various decarbonization strategies. In contrast to previous work, the consequences of a transformation to fully decarbonized transport system scenarios are quantified across all traffic segments, considering procurement, operation and disposal. The methodology can be applied to arbitrary regions and transport systems. Here, the metropolitan region of Berlin is chosen as a demonstration case. First results are shown for a complete conversion of all traffic segments from conventional propulsion technology to battery electric vehicles. The transition of private individual traffic is analyzed regarding technical feasibility, energy demand and environmental impact. Commercial goods, municipal traffic and public transport are analyzed with respect to system cost and environmental impacts. We can show a feasible transition path for all cases with substantially lower greenhouse gas emissions. Based on current technologies and today’s cost structures our simulation shows a moderate increase in total systems cost of 13-18%.

Carbon dioxide capture and storage (CCS) is a necessary requirement for high-emitting CO2 industries to significantly reduce volumes of greenhouse gases released to the atmosphere and mitigate the climate change. Geological CO2 storage into depleted oil and gas fields is the fastest and most accessible option for CCS deployment allowing for re-purposing of existing infrastructure and utilizing the significant knowledge about the subsurface acquired during the field production operations. Location of such depleted fields in neighbourhood of high-emitting CO2 industries is additional advantage available for matured on-shore Europe fields. Considering these advantages, the oil and gas operators are now evaluating different possibilities for CO2 sequestration projects for the fields approaching end of production. This article describes an integrated approach to reservoir simulations focused on evaluating CO2 injection pilot at one of such matured fields operated by MND and located in the Czech Republic. The potential CO2 injection site in focus is a naturally fractured carbonate reservoir. This oil-bearing formation has a gas cap at the top and a limited aquifer at the bottom of the reservoir and was produced mainly by pressure depletion with quite limited pressure support from water injection. The article summarizes results of efforts taken by the multi-disciplinary team to develop and apply an integrated approach to reservoir simulation started from geological modelling of the naturally fractured reservoir and integrating results of laboratory studies and their interpretation (geomechanical and geochemistry laboratory experiments) and dynamic field data analysis (pressure transient analysis including time-lapse) into the history matching of the reservoir model and simulations of pilot CO2 injection. Attention is given to evaluating and modelling stress-sensitive fracture properties and safe injection envelope preventing induced fracturing as well as impact of potential salt precipitation in the near wellbore area. These effects are considered in the context of pilot CO2 injection and addressed in reservoir simulations of the injection scenarios. Single porosity and permeability reservoir simulations with dominating fracture flow and black-oil formulation with CO2 simulated as a solvent were used in this study and arguments for such choice of the simulation approach for the site in focus are shortly discussed. The reservoir simulations have indicated larger site injection capacity than required for the pilot injection and allowed to evaluate CO2 migration pathway within the reservoir. The application of the approach to the site in focus has also revealed large uncertainties, related to fracture description and geomechanical evaluations, resulted in an uncertain safe injection envelope. These uncertainties should be addressed in further studies in preparation to the pilot. The article concludes with an overview of the outcomes of the integrated approach and its application to the field in focus including discussion of issues and uncertainties revealed.

Long-range interactions are relevant for a large variety of quantum systems in quantum optics and condensed matter physics. In particular, the control of quantum-optical platforms promises to gain deep insights in quantum-critical properties induced by the long-range nature of interactions. From a theoretical perspective, long-range interactions are notoriously complicated to treat. Here, we give an overview of recent advancements to investigate quantum magnets with long-range interactions focusing on two techniques based on Monte Carlo integration. First, the method of perturbative continuous unitary transformations where classical Monte Carlo integration is applied within the embedding scheme of white graphs. This linked-cluster expansion allows to extract high-order series expansions of energies and observables in the thermodynamic limit. Second, stochastic series expansion quantum Monte Carlo which enables calculations on large finite systems. Finite-size scaling can then be used to determine physical properties of the infinite system. In recent years, both techniques have been applied successfully to one- and two-dimensional quantum magnets involving long-range Ising, XY, and Heisenberg interactions on various bipartite and non-bipartite lattices. Here, we summarise the obtained quantum-critical properties including critical exponents for all these systems in a coherent way. Further, we review how long-range interactions are used to study quantum phase transitions above the upper critical dimension and the scaling techniques to extract these quantum critical properties from the numerical calculations.

The article is devoted to the preliminary concept of the Future Planetary Defense System (FPDS), emphasizing Astroballistics. This paper is intended to support international efforts to improve the planetary security of the Earth. The work covers three areas of knowledge: Astronautics, Astrodynamics, and Astroballistics. The most important part of the presented article is dynamic, contact combat modeling against small, deformable celestial bodies. For these purposes, the original, proprietary hydrocode of the Free Particle Method (HEFPM-G) with gravity was used. The main aim of combat is to redirect the Potentially Hazardous Objects (PHOs) to orbits safe for the Earth or destroy them. This concept's first task is to find, prepare and use dynamic three-dimensional models of celestial bodies' motion and spacecraft or human-crewed spaceships in the Solar System's relativistic frame. The second task is to prepare the FPDS' architecture and computer simulation space missions' initial concepts in the internal part of the solar system. The third and main task covers simulating, using hydrocodes, selected methods of fighting 100 m and 140 m diameter asteroids. The order of the article is as follows. The first part of the article presents an architecture and FPDS' modus operandi. Preliminary design and development of FPDS' space missions, including navigation, mission dynamics simulation, is prepared using an open-source space mission analysis and design tool. E.g., Asterank and Trajectory Browser or GMAT are presented in the second part. The third part of the article is devoted to computer hydrocodes (HEFPM-G) and the modeling and simulation of asteroid–asteroid collision, laser radiation effects on an asteroid, and FPDS spacecraft's warhead contact interaction on the small celestial body. The authors formulated the Main Planetary Defense Problem (MPDP) in this paper. The proposition of this problem solving has been realized by preparing its concepts, architecture, and modus operandi of the FPDS mission. Finally, a series of realistic simulations were made using hydrocode to deflect or destroy dangerous asteroids. The summary and conclusions can be found in the fourth part of the article.

For evaluating probabilities of arbitrary random events with respect to a given multivariate probability distribution, specific techniques are of great interest. An important two-dimensional high risk limit law is the Gauss-exponential distribution whose probabilities can be dealt with based upon the Gauss-Laplace law. The latter will be considered here as an element of the newly introduced family of (p,q)-spherical distributions. Based upon a suitably defined non-Euclidean arc-length measure on (p,q)-circles we prove geometric and stochastic representations of these distributions and correspondingly distributed random vectors, respectively. These representations allow dealing with the new probability measures similarly like with elliptically contoured distributions and more general homogeneous star-shaped ones. This is demonstrated at hand of a generalization of the Box-Muller simulation method. En passant, we prove an extension of the sector and circle number functions.