3D Numerical Model Research Articles

How Switchbacks Can Maintain a Longer Time in the Interplanetary Space

Parker Solar Probe, the closest spacecraft to the Sun, has renewed our understanding of the solar corona and the interplanetary space. One of its important findings is the prevalence of switchbacks, which display localized magnetic reversals along the otherwise Parker spirals. While some switchbacks might disappear quickly, others can maintain for a long period of time, and there are indications that many switchbacks strengthen from the solar corona to the interplanetary space despite their magnetic tension force, which tends to straighten the magnetic field lines. Therefore, how these switchbacks could be maintained for a long period of time remains a mystery. In this paper, we employed a 3D data-driven global full magnetohydrodynamics numerical model to explore the evolution of switchbacks formed in the dynamic corona. Our simulations indicate that two factors can affect the lifetime of a switchback. One factor is the combination of angle and leg length ensures that the switchback with greater curvature after reconnection can last longer, and the greater the angle, the more magnetic field lines that can be reconnected, and thus the longer the duration. We call this influencing factor flux tube shape factor. The other factor is the velocity shear, i.e., when the solar wind at the convex-outward turning of a switchback is faster than that at the concave-outward turning, the switchback becomes enhanced during propagation, and it weakens when the velocity difference is opposite.

Lab-scale validation and parametric numerical investigation on hydrogen production in a porous media as a byproduct of in-situ combustion

A numerical 1D model was developed and validated using the measured data on a porous tube lab experiment for hydrogen generation as an in-situ combustion (ISC) byproduct, and a detailed parametric numerical study was performed. This work explores the impacts of various water injection parameters on the generation of hydrogen, the oil recovery factor (ORF), and carbon monoxide using a thorough parametric analysis. The methodology included simulating the injection of mutually enriched air and water into the combustion tube, focusing on parameters such as water temperature, quality, and flow rates. Key findings reveal that the oxygen ratio in the oxidizer substantially influences hydrogen production and ORF, with hydrogen generation increasing from 1.65 × 103 cm3 to 3.0 × 103 cm3, which is around 82 % when the oxygen percentage is raised from 50 % to 95 %. The results showed that employing wet combustion instead of dry combustion increased the hydrogen production rate by roughly four times. Variable water temperature has insignificant impacts on the hydrogen production rate and ORF. Increasing the steam quality has an opposed effect on the hydrogen generation rate. Additionally, escalating the injected water flow rate from 3000 cm³/day to 15000 cm³/day boosts hydrogen production from 10000 cm³ to 53400 cm³, respectively.

In silico assessment of the reliability and performance of artificial sphincter for urinary incontinence.

The standard artificial urinary sphincter (AUS) is an implantable device for the treatment of urinary incontinence by applying a pressure loading around the urethra through an inflatable cuff, often inducing no-physiological stimulation up to tissue degenerative phenomena. A novel in silico approach is proposed to fill the gap of the traditional procedures by providing tools to quantitatively assess AUS reliability and performance based on AUS-urethra interaction. The approach requires the development of 3D numerical models of AUS and urethra, and experimental investigations to define their mechanical behaviors. Computational analyses are performed to simulate the urethral lumen occlusion by AUS inflation under different pressures, and the lumen opening by applying an intraluminal pressure progressively increased under the AUS action (Abaqus Explicit solver). The AUS reliability is evaluated in terms of tissue stimulation by the mechanical fields potentially responsible for vasoconstriction and tissue damage, while the performance by the intraluminal pressure that causes the lumen opening for a specific occlusive pressure, showing the maximum urethral pressure for which continence is guaranteed. The present study implemented the procedure considering the gold standard AMS 800 and a novel patented AUS. Results provided the comparison between two sphincteric devices and the evaluation of the influence of different building materials and geometrical features on the AUS functionality. The approach was developed for the AUS, but it could be adapted also to artificial sphincters for the treatment of other anatomical dysfunctions, widening the analyzable device configurations and reducing experimental and ethical efforts.

Influence of bottom protection structures group on the hydraulic characteristics in sharp bends

Bottom protection (BP) is a commonly used engineering measure in river training. A 3D hydrodynamic numerical model is established to study the influence of the overall arrangements of BP structures on the hydraulic characteristics of the sharp bends. Under the action of the BP structures group, the maximum longitudinal flow velocity is closer to the water surface vertically and closer to the concave bank laterally compared with no BP. The high shear stress zones on the riverbed at the exit straight section of the bend do not occur. The so-called primary circulation after the apex of the bend gradually weakens and disappears. If the coverage rate of BP is increased to 1.5 times that of the sequential BP, the average longitudinal flow velocity near the bottom of BP decreases by 10% and the average bed shear stress at the outlet section of the bend decreases by 35%. In the staggered BP condition, the average bed shear stress at the outlet section of the bend is reduced by 23% compared with the sequential BP. Moreover, the results show that the protection effects of the dense BP and the staggered BP are better than that of sequential BP.

Numerical Simulation of Flow Field Structure of Vanadium Redox Flow Battery and its Optimization on Mass Transfer Performance

The structural design of the flow channel of a redox flow battery directly affects ion transport efficiency, electrode overpotential, and stack performance during charge-discharge cycles. A tapered hierarchical interdigitated flow field design that has independent flow channel structures for different levels of flow was developed in this work. Especially, the secondary branch channels are the tapered type and the corresponding cross-sections are gradually reduced along the flow direction, which is beneficial for improving the flow rate at the end of channels and enhancing mass transfer. The performances of a vanadium redox flow battery with interdigitated flow field, hierarchical interdigitated flow field, and tapered hierarchical interdigitated flow field were evaluated through 3D numerical model. The results showed that at 240 mA cm−2 and 6 ml s−1, the pump-based efficiency of the hierarchical interdigitated flow field increased by 4%-7% compared with the interdigitated flow field. Furthermore, the pump-based efficiency with tapered hierarchical interdigitated flow field increased by 1.6%-3% compared with the hierarchical interdigitated flow field. This indicates that the tapered hierarchical interdigitated flow field shows further advantages in redox flow battery applications.

Application of artificial countermeasures to enhance desilting efficiency in a reservoir under normal and extreme events

ABSTRACT The influx of flood-induced turbidity flows with high sediment concentrations into the reservoir can create sediment deposition issues, potentially jeopardizing water supply reliability. Efficient management of reservoir water resources is imperative for sustainable utilization. Evaluating the need for reservoir desilting involves comprehending the inflow and outflow dynamics of water-sediment discharges, the capacity of existing outlet structures, reservoir operation objectives, and available desilting technologies. However, given the altering hydrological conditions due to the impacts of climate change, there is a necessity for adaptive methods to ensure reservoir storage maintenance. Consequently, this study introduces the utilization of a flushing channel with artificial guiding structures to enhance desilting efficiency. A 2D numerical model offers convenient computational tools for assessing water-sediment transport behaviors across various operational scenarios based on reservoir management strategies. Data obtained from field observations is collected and analyzed to calibrate and verify the numerical model. By employing 2D numerical modeling, sediment concentration and desilting efficiency are calculated to support effective desilting operations within a reservoir. Simulation results indicate that the integration of additional artificial guiding structures at the bottom of the reservoir can enhance sedimentation mitigation outcomes by approximately 17.63% under normal conditions and by 5.27% under extreme hydrological conditions.

Toward Building a Virtual Laboratory to Investigate Rainfall Microphysics at Process Scales

Abstract A 3D numerical model was built to serve as a virtual microphysics laboratory (VML) to investigate rainfall microphysical processes. One key goal for the VML is to elucidate the physical basis of warm precipitation processes toward improving existing parameterizations beyond the constraints of past physical experiments. This manuscript presents results from VML simulations of classical tower experiments of raindrop collisional collection and breakup. The simulations capture large raindrop oscillations in shape and velocity in both horizontal and vertical planes and reveal that drop instability increases with diameter due to the weakening of the surface tension compared with the body force. A detailed evaluation against reference experimental datasets of binary collisions over a wide range of drop sizes shows that the VML reproduces collision outcomes well including coalescence, and disk, sheet, and filament breakups. Furthermore, the VML simulations captured spontaneous breakup, and secondary coalescence and breakup. The breakup type, fragment number, and size distribution are analyzed in the context of collision kinetic energy, diameter ratio, and relative position, with a view to capture the dynamic evolution of the vertical microstructure of rainfall in models and to interpret remote sensing measurements. Significance Statement Presently, uncertainty in precipitation estimation and prediction remains one of the grand challenges in water cycle studies. This work presents a detailed 3D simulator to characterize the evolution of drop size distributions (DSDs), without the space and functional constraints of laboratory experiments. The virtual microphysics laboratory (VML) is applied to replicate classical tower experiments from which parameterizations of precipitation processes used presently in weather and climate models and remote sensing algorithms were derived. The results presented demonstrate that the VML is a robust tool to capture DSD dynamics at the scale of individual raindrops (precipitation microphysics). VML will be used to characterize DSD dynamics across scales for environmental conditions and weather regimes for which no measurements are available.

Using a three-dimensional computational fluid dynamics model and an agglomerate model to investigate the effect of varying agglomerate parameters and output voltages on proton exchange membrane fuel cell performance

Modeling of proton exchange membrane (PEM) fuel cells is attracting more attention as fuel cell technology continues to develop. In this study, we considered a hybrid model that combines an agglomerate model based on the agglomeration of catalyst particles and the coverage-dependent kinetic equation of platinum oxide for ORR, and another 3D numerical model of a PEM fuel cell based on computational fluid dynamics (CFD). The obtained results from our developed models were validated with experimental results from literature. In fact, we investigated the effects of changing the agglomerate radius (Ragg), the ionomer volume fraction within the agglomerate (Yi,agg), the effective agglomerate surface area (Ai,agg), the distribution of the gases and the temperature on the cell performances.The results revealed that the cell performances are strongly influenced by changing Ragg and Yi,agg for medium and high current densities: The activation loss increases with increasing Ragg and decreasing Yi,agg. Also, Ai,agg increases with decreasing Ragg and increasing Yi,agg. In addition, the PEM fuel cell's power output is significantly enhanced when Ragg is decreased and Yi,agg is increased, the optimal power being obtained for values of Ragg=100nm and Yi,agg = 0.6. The numerical results also showed that decreasing the output voltage from 0.95V to 0.35V can accelerate the electrochemical reaction process.

Shaking table test and numerical analyses of a multi‐story traditional tower‐style building

AbstractThe traditional tower‐style building (TTSB) is an innovative structural form constructed in modern cities imitating the overall appearance of ancient timber pagodas, and it is an extraordinarily cultured high‐rise construction. The limited stipulations for high‐rise TTSBs in the seismic design code pose challenges in assessing the reliable performance of the unique structure subjected to severe earthquakes. This paper presents the shaking table test and numerical analyses of a 1/15‐scale 13‐story TTSB specimen with an integral tower height of 5.36 m, which is a steel frame‐braced core‐tube structure consisting of seven bright floors (built‐out stories evident from the outside, BF) and six dim floors (built‐in stories not apparent from the outside, DF), with transition columns and stiffening trusses around the exterior perimeter. The experimental results showed that the tested tower‐style building had excellent seismic performance and reliable structural integrity. It only experienced minor damage when subjected to extremely high‐intensity motions. The interstory drift, dynamic strain, and floor acceleration response of the core‐tube region with eccentric steel braces were more significant under severe excitation than those with cross‐centered symmetric steel braces. The vertical reaction increased at the upper floor of the tower under vertical acceleration, and differences in the dynamic response of the middle and upper floors were much more apparent after the test. Moreover, 3D numerical simulation models of the tested tower were established and validated against the test responses. Successively, the validated numerical model was used to investigate the influence of the transition column at different floors on the peak interstory drift response and the relevant strain distribution, and the proposal for a proper position of the transition column was recommended at the end.

Predicting weld pool metrics in laser welding of aluminum alloys using data-driven surrogate modeling: A FEA-DoE-GPRN hybrid approach

Multi-physics computational models based on finite element analysis, offer detailed insights into the dynamics and metrics in the weld pool formed by laser welding. Conversely, data-driven surrogate models provide a cost-effective means to predict desired responses. These models establish statistical or mathematical correlations with input–output data, eliminating the need for additional simulations during design optimization. This study proposes a data-driven surrogate model, employing the Gaussian process regression network (GPRN), to predict weld pool metrics, such as weld width and depth of penetration in laser welding of aluminum alloy. A 3D computational fluid dynamics-based numerical model is initially constructed and experimentally validated to predict weld pool metrics. Subsequent experimental runs, guided by the design of experiments, include various configurations of process parameter settings. The developed numerical model computes weld pool metrics for each experimental run, forming a dataset for training and testing the GPRN model. The GPRN model is evaluated against simulated data, showing adequacy with a mean square error of 1.7 µm and mean absolute percentage error of 10−7, with experimental validation further confirming its accuracy, revealing a minimum error of 1.7%, a maximum error of 8%, and an average error of 3%. The key contribution and novelty of this study lie in the development of the hybrid data-driven model, which accurately predicts weld pool metrics while minimizing experimental and computational efforts.

3D Numerical Modeling for Investigating Structural Controls on Orogenic Gold Mineralization, Sanshandao Gold Belt, Eastern China

3D Numerical Modeling for Investigating Structural Controls on Orogenic Gold Mineralization, Sanshandao Gold Belt, Eastern China

Behavior of rock bolt plates made of engineering fiber composites: Experimental investigations, 3D surface characterization and numerical models

Bearing plates made from plastic composites can be used as an alternative to their steel counterparts in rock bolt or soil nail applications. To achieve this goal, an existing recycled high-density polyethylene bearing plate was investigated and later modified to improve its engineering properties. Laboratory studies were conducted to determine the failure load of the existing and modified plates, and a numerical model was developed for complementary analysis. The results of both efforts clearly showed that the existing bearing plate was not adequate in terms of strength and creep properties, as it quickly yielded with large displacements at relatively low loads. In order to enhance the strength of the plate, both geometric and material modifications are made by our research group to obtain a more efficient plate. Numerical models were used to determine the frame layout, and a series of analyses were performed to evaluate the effects of frame thickness, number and arrangement. Once the design was optimized and finalized, a mold was created to match the new geometry for manufacturing new plates through injection molding. A test setup was also established in the laboratory and numerous compression tests were performed on the manufactured new plates. The measured load-displacement behavior of plates made of polyethylene and polyamide with a variety of additives were discussed separately. It was determined that the new plastic plates reinforced with polyamide through various additives have the potential to reach a strength up to 200–240 kN, which is at least two times higher than the existing one, with distinct economic advantages.

Harbor Sedimentation Management Using Numerical Modeling and Exploratory Data Analysis

Sedimentation in the harbors’ basins is an environmental phenomenon that frequently disrupts safe shipping and necessitates costly dredging operations. The layout of harbors and the permeability of protective structures such as breakwaters influence sediment transport within harbor basins. Thus, through a multistep framework, this study investigates the sedimentation management issues for the Egyptian proposed Ezbet Elborg fishing harbor based on field measurements and a numerical morphodynamic coastal modeling system (CMS). First, field measurements were analyzed and evaluated for acquiring a full grasp of the research area’s bathymetry and hydrodynamics. Second, a two-dimensional (2D) numerical simulation CMS model was set up and calibrated against field measurements wherein the developed CMS model highly correlated with actual measurements by 97%. CMS results demonstrate that the predominant NNW wave with the formed longshore current on both the harbor’s sides affects sediment accumulation within the harbor’s basin. Third, 100 simulations for the proposed harbor including different structural modulation scenarios affecting the sedimentation issue were investigated via the calibrated CMS model. Finally, an exploratory data analysis (EDA) is performed via correlation matrix and ANOVA test for the CMS’s scenarios’ results to gain an in-depth view of the relation between the harbors’ layout and the structural characteristics with the sedimentation volumes. Results showed that breakwaters’ orientation affects sediment accumulation more than its length. Also, breakwater permeability and basin width are significantly affecting sediment accumulation. Ultimately, the current study makes a substantial contribution to integrated coastal structure management (ICSM) by helping coastal stakeholders to mitigate the negative impacts of the harbors’ sediment deposition aiming at sustaining both environmental and economic aspects.

Detailed modeling and numerical analysis of thermo mechanical stresses in the crystalline silicon and thin film PV modules under varying climatic conditions

The present manuscript provides an in-depth analysis of the modeling and simulation of thermo-mechanical stresses in photovoltaic panels using finite element (FEM) analysis. It covers important aspects such as the choice of the model, including governing equations, boundary conditions, and simplifying assumptions. The study mainly focuses on evaluating the impact of environmental conditions on the thermo-mechanical stresses of first-generation c-Si and second-generation thin-film CdTe and a-Si PV panels in five different climatic zones in Algeria. The proposed transient 3D numerical model is validated against experimental measurements by other researchers and the tabulated values of the nominal operating temperature of the cell (NOCT). The obtained results indicate that the first principle stress for the c-Si cell is lower compared to CdTe and a-Si cells; this can be explained by looking at the thickness, as well as the Young modulus of the solar cell. Consequently, in regions or countries with hot climatic conditions or large temperature gradients, it is recommended to use PV panels based on first-generation c-Si cells with glass-glass as the front and back covers.

GOLD PLASMONIC ARRAY STRUCTURES FOR SENSING APPLICATIONS

This article is devoted to the theoretical study of the plasmonic properties of periodically arranged arrays of gold nanoparticles. The Comsol Multiphysics software, which is based on the finite element method, was used to build 3D numerical models for the simulation and conduct research. In this work the electric field distribution and optical characteristics of the spherical gold nanoparticles array were studied. Individual localized surface plasmon resonance modes are influenced when metallic nanoparticles are in the close proximity and as a result the electric near- fields can couple, resulting in a new hybrid mode. We mainly focused here on the investigation of two crucial questions, particularly, influences of the gap between the nanoparticles and the refractive index of the surrounding medium on the resulting optical response of the gold nanoparticles arrays. The array of periodically arragement gold nanoparticles is characterized by an enhanced local electric field between the nanoparticles, which is inversely proportional to the gap between the particles. The field strength and optical properties (reflection, transmission, and absorption) can be conveniently manipulated by changing the gap between particles. In additional, their potential applications as sensetive plasmonic sensors element have been considered. The studied structure has a significant potential for practical applications due to its wide range of the operating wavelengths and ease of the high-throughput fabrication. In the course of the study, it was established that the change in the distance between the surface of nanoparticles by 1 nm leads to a significant shift in the spectral transmission and reflection curves on the spectral range. In addition, these studies showed that an increase in the distance between the surfaces of nanoparticles leads to the decrease in the near-field interaction between gold nanoparticles in the array. Therefore, the obtained results can be successfully used in the manufacture of highly sensitive plasmon sensors with the possibility of controlling the sensitivity and the working spectral range.

Quantitative analysis of asphalt concrete’s tension-compression asymmetry effects on pavement response through 3D numerical modeling with dual viscoelastic model

Asphalt concrete (AC) exhibits noticeable tension-compression (TC) asymmetry, but it is typically considered isotropic in pavement design. This study aims to quantitatively evaluate the effect of AC’s TC asymmetry on pavement response under loading through numerical modeling. To achieve this objective, a temperature-dependent dual viscoelastic constitutive model was applied to incorporate AC’s TC asymmetry into pavement modeling. Besides, three pavement structure models, including one with thick AC layers, one with thin AC layers, and one with a Portland cement concrete (PCC) base, were developed. The responses of the three pavement structures under traffic and environmental loading conditions were simulated. Modeling results showed that AC’s TC asymmetry can significantly increase the vertical strain in AC, leading to higher stress concentration and larger deformation in AC layers. Unlike the conventional understanding that tensile strain concentrates at the bottom of AC layers, high tensile strains were observed in the top AC layers, especially for the pavement with thick AC layers or a PCC base, which may lead to top-down cracking. High tensile strains were observed on the granular subbase for the pavement with thin AC layers, which may induce bottom-up cracking. Besides, AC’s TC asymmetry also significantly increases the shear strains, especially the horizontal shear strain in AC layers, which may result in debonding and shoving distresses in AC layers. It was also noticed that the pavement response highly depends on its temperature and vehicular speed. A higher temperature or lower vehicular speed leads to more significant AC’s TC asymmetry as well as its effects on asphalt pavement’s response. The outcomes of this study are expected to help enhance the design and maintenance of more durable asphalt pavement.

Unraveling dispersion and buoyancy dynamics around radial A + B → C reaction fronts: microgravity experiments and numerical simulations

Radial Reaction–Diffusion–Advection (RDA) fronts for A + B → C reactions find wide applications in many natural and technological processes. In liquid solutions, their dynamics can be perturbed by buoyancy-driven convection due to concentration gradients across the front. In this context, we conducted microgravity experiments aboard a sounding rocket, in order to disentangle dispersion and buoyancy effects in such fronts. We studied experimentally the dynamics due to the radial injection of A in B at a constant flow rate, in absence of gravity. We compared the obtained results with numerical simulations using either radial one– (1D) or two–dimensional (2D) models. We showed that gravitational acceleration significantly distorts the RDA dynamics on ground, even if the vertical dimension of the reactor and density gradients are small. We further quantified the importance of such buoyant phenomena. Finally, we showed that 1D numerical models with radial symmetry fail to predict the dynamics of RDA fronts in thicker geometries, while 2D radial models are necessary to accurately describe RDA dynamics where Taylor–Aris dispersion is significant.

3D Numerical Modeling of Negative Skin Friction Effects on Tapered Piles in Soft Soils

3D Numerical Modeling of Negative Skin Friction Effects on Tapered Piles in Soft Soils

The influence of viscous slab rheology on numerical models of subduction

Abstract. Numerical models of subduction commonly use diffusion and dislocation creep laws from laboratory deformation experiments to determine the rheology of the lithosphere. The specific implementation of these laws varies from study to study, and the impacts of this variation on model behavior have not been thoroughly explored. We run simplified 2D numerical models of free subduction in SULEC, with viscoplastic slabs following (1) a diffusion creep law, (2) a dislocation creep law, and (3) both simultaneously, as well as several variations of model 3 with reduced resistance to bending. We compare the results of these models to a model with a constant-viscosity slab to determine the impact of the implementation of different lithospheric flow laws on subduction dynamics. In creep-governed models, higher subduction velocity causes a longer effective slab length, increasing slab pull and asthenospheric drag, which, in turn, affect subduction velocity. Numerical and analogue models implementing constant-viscosity slabs lack this feedback but still capture morphological patterns observed in more complex models. Dislocation creep is the primary deformation mechanism throughout the subducting lithosphere in our models. However, both diffusion creep and dislocation creep predict very high viscosities in the cold core of the slab. At the trench, the effective viscosity is lowered by plastic failure, rendering effective slab thickness the primary control on bending resistance and subduction velocity. However, at depth, plastic failure is not active, and the viscosity cap is reached in significant portions of the slab. The resulting high slab stiffness causes the subducting plate to curl under itself at the mantle transition zone, affecting patterns in subduction velocity, slab dip, and trench migration over time. Peierls creep and localized grain size reduction likely limit the stress and viscosity in the cores of real slabs. Numerical models implementing only power-law creep and neglecting Peierls creep are likely to overestimate the stiffness of subducting lithosphere, which may impact model results in a variety of respects.

Numerical study of thermal performance of a PCM in a modifying tube-bundle latent heat thermal storage

The main aim of this study is to numerically investigate the thermal performance of a newly designed modifying tube-bundle Latent Heat Thermal Storage (LHTS). The new LHTS design proposed a half-cylindrical tube attached to the Heat Transfer Fluid tube, which implies making the PCM indirect contact with the HTF tube. This reduces the heat transfer resistance by eliminating the metal construction of the tubes. For this purpose, a 3D numerical model was developed based on mass, momentum and energy conservations using ANSYS Fluent software. After validating the model developed by comparing numerical results with those obtained experimentally by a published study, the PCM liquid fraction and amount of power stored during the melting and solidification were studied. Water and lauric acid as a working fluid and PCM were implemented in the model. In addition, three different Heat Transfer Fluid (HTF) inlet temperatures 15°C,20°Cand25°C, and HTF’s flow rates in terms of Reynolds number 1000,1500and2000 were tested. The thermal behaviour of the present LHSU design was compared with that of the same volume double pipe LHTS unit. Results showed that the melting and solidification processes of the proposed storage were very uniform, with convection and conduction heat dominating PCM melting and solidification, respectively. The increase in the heat transfer contact area between the HTF and the tube bundle contributed significantly to the PCM melting and solidification acceleration, making each tube bundle behave as a single storage unit. Therefore, the time required for complete melting and solidifying of the PCM was half (50%) of the time needed for the same amount in the double pipe LHTS. Similarly, the melting and solidification time in the preset storage was shortened by about 22% and 41.6%, respectively, compared with the regular tube LHTS unit.