Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

A Multiscale Modeling Approach to Cryocompressed Hydrogen Storage Pressure Vessels – Part II: Constitutive Modeling and Finite Element Analysis

Subjected to thermomechanical cycling, composite cryogenic hydrogen (H2) storage pressure vessels experience high stresses in the carbon fiber (CF)/epoxy overwrap which can lead to vessel failure due to a combination of degradation mechanisms such as matrix cracking, fiber/matrix debonding, delamination, and fiber rupture. The present paper is the second part of a set of two articles on this subject that addresses analysis and design of these pressure vessels. Predictive finite element (FE) modeling capabilities have been used to support a material acceptance process that aims at providing guidance to evaluate specialty resins, vessel liner options, and CF composites through thermomechanical testing. We have applied a multiscale modeling approach recently developed (B.N. Nguyen et al., International Journal of Hydrogen Energy, 2019, https://doi.org/10.1016/j.ijhydene.2019.09.200) to perform three-dimensional (3D) FE analysis and design of a cryo-compressed H2 storage pressure vessel for its operation in a large temperature range from room to cryogenic temperatures. The developed approach termed the micro-meso-macro approach determines the laminar stresses in different layers of the vessel in addition to constituent (i.e., fiber and matrix) stresses that are important to the design of the composite overwrap for its layup and material combination to reduce the risk of vessel failure.

Nguyen, Ba Nghiep↗

Assessing pore network heterogeneity across multiple scales to inform CO2 injection models

Geologic heterogeneity is a key feature that must be considered when translations of scaled data are performed. This paper presents the assessment of geologic heterogeneity using a multiscale workflow that includes image analysis-based methods coupled with well log analysis to provide data in which fractals and machine learning methods estimate the carbon dioxide (CO 2 ) storage resource potential of a reservoir. The heterogeneity of rock properties of the complex Bell Creek reservoir in Montana, USA, was explored at the pore scale (~nm to mm), core scale (~mm to m), and well scale (~cm to m). The data used in this study included advanced image analysis of micro-CT (computed tomography) images (pore scale), thin sections (pore scale), plugs and core images (core scale) and well logs (well scale). The micro-CT images were segmented using a U-net segmentation approach into objects of pores and grains. Further, the segmented images were reconstructed into subvolumes of different sizes. Physical properties (porosity and permeability) and fractal dimensions were calculated for the various subvolumes, and Lorenz coefficient (Lc) values, a single parameter to describe the degree of heterogeneity within a pay zone section, were calculated from thin-section images and well logs. Porosity and fractal dimension values were used to estimate the 188-µm threshold of representative elementary volume (REV) in this study. Both the Lc and fractal dimension values were found to be negatively correlated. When these two parameters are combined, it is possible to discern differences in the complex porous networks of the samples analyzed in this study.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Quantitative Measurement of Li-Ion Concentration and Diffusivity in Solid-State Electrolyte

Here, we present a quantitative method to measure the concentration and diffusivity of Li ion in a solid-state electrolyte with nanoscale depth resolution. We designed a standard sample with different depths of trenches and used time-of-flight secondary ion mass spectroscopy, inductively coupled plasma optical emission spectroscopy, three-dimensional (3D) optical microscopy, and electrochemical strain microscopy to obtain effective Li-ion concentration and the corresponding diffusivity as a function of depth. As a result, we verified the dependence of electromechanical strain induced by ionic oscillation on mobile-ion concentration. We also showed the diffusivity change as a function of effective concentration and discussed the difference in diffusivity near the surface and inside the bulk. Our method can be widely used for electromechanical strain in various ionic conductors, rendering our technique universal for quantitative analysis of ionic dynamics at a multiscale.

36 MATERIALS SCIENCE↗

Reduced models for ETG transport in the tokamak pedestal

This paper reports on the development of reduced models for electron temperature gradient (ETG) driven transport in the pedestal. Model development is enabled by a set of 61 nonlinear gyrokinetic simulations with input parameters taken from pedestals in a broad range of experimental scenarios. The simulation data have been consolidated in a new database for gyrokinetic simulation data, the multiscale gyrokinetic database (MGKDB), facilitating the analysis. The modeling approach may be considered a generalization of the standard quasilinear mixing length procedure. The parameter η, the ratio of the density to temperature gradient scale length, emerges as the key parameter for formulating an effective saturation rule. With a single order-unity fitting coefficient, the model achieves an error of 15%. A similar model for ETG particle flux is also described. We also present simple algebraic expressions for the transport informed by an algorithm for symbolic regression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Species Transport Framework Development in SAM for System-Level Tritium Source Term Analysis

The SAM code is under development as a modern system-level modeling and simulation tool for advanced non–light water reactor safety analyses, with recent efforts to add capabilities to evaluate radiological source term risks in these novel reactor concepts. By leveraging the established system-level multiphysics thermal-hydraulic models in SAM, a framework for tightly coupled species transport modeling has been integrated into the code for engineering-scale source term evaluation. This species transport framework was first applied to the simulation of tritium, which is a well-known source term in conventional light water reactors. Tritium poses a unique risk in salt-cooled reactors, especially those with lithium-bearing salts such as the fluoride salt–cooled high-temperature reactor (FHR) concept, as tritium is generated in the salt coolant in significant quantities due to neutron interactions. A compounding factor is the increased mobility of tritium at high temperatures, which is able to permeate through metals while also potentially being retained in graphite pebbles and structures. Engineering-scale models for the tritium transport pathways in a FHR have been developed using the new species transport framework in SAM. The capabilities are assessed through analytical verification problems and validated with data from a graphite retention experiment. In conclusion, the system-level model is demonstrated by performing an initial estimate of baseline tritium generation and flows in a generic reference SAM FHR model, setting a foundation for future studies of source term transient analysis with the potential for further multiscale and multiphysics integration.

SAM↗

Streaming Data in HPC Workflows Using ADIOS

The “IO Wall” problem, in which the gap between computation rate and data access rate grows continuously, poses significant problems to scientific workflows which have traditionally relied upon using the filesystem for intermediate storage between workflow stages. One way to avoid this problem in scientific workflows is to stream data directly from producers to consumers and avoiding storage entirely. However, the manner in which this is accomplished is key to both performance and usability. This paper presents the Sustainable Staging Transport, an approach which allows direct streaming between traditional file writers and readers with few application changes. SST is an ADIOS “engine”, accessible via standard ADIOS APIs, and because ADIOS allows engines to be chosen at run-time, many existing file-oriented ADIOS workflows can utilize SST for direct application-to-application communication without any source code changes. This paper describes the design of SST and presents performance results from various applications that use SST, for feeding model training with simulation data with substantially higher bandwidth than the theoretical limits of Frontier’s file system, for strong coupling of separately developed applications for multiphysics multiscale simulation, or for in situ analysis and visualization of data to complete all data processing shortly after the simulation finishes.

Podhorszki, Norbert [ORNL] (ORCID:000000019647542X↗

Analyzing the Multiscale Impacts of Implementing Energy-Efficient HVAC Improvements Through Energy Audits and Economic Input–Output Analysis

Abstract Heating, ventilation, and air-conditioning (HVAC) systems are usually an industry’s highest consumer of energy, most of which goes toward space cooling in buildings. Industrial energy-efficiency audits not only benefit manufacturers but also generate significant economic and environmental benefits to localities, states, and the nation. This article analyzes the micro- and macro scale impacts of implementing energy-efficient HVAC systems by integrating the industrial building energy data with the macroeconomic regional economic flow model. Micro-scale data include 10 years of historical energy, cost, and carbon dioxide savings achieved from energy-efficient HVAC implementation offered to manufacturers through industrial energy audits. The data were integrated into the macroeconomic modeling framework to illuminate the cascading regional economic impacts of implementing energy-efficient HVAC recommendations in manufacturing facilities. Results show that if recommendations had been implemented throughout all manufacturers in the region, $656 M energy costs would have been directly saved, 7.8 million metric tons of carbon dioxide emissions would have been avoided, and 4387 jobs could have been created, resulting in a total annual economic impact of $899 M stemming from direct, indirect, and induced impacts. The results offer insight into how industrial energy systems can be designed and provide models for how communities can accomplish a net-zero society.

Energy & Fuels↗

Spatial Optimization of Multiscale Biorefinery Deployment for a Diversified Bioeconomy in the United States

Strategic biorefinery siting is critical for a diversified bioeconomy, yet industry, policy, and research often focus on either large-scale biofuel plants or smaller-scale specialty bioproduct facilities, with limited coordination across scales. We address this gap by modeling biorefinery deployment spanning a 28-fold difference in capacity. We developed an open-source, spatially explicit framework integrating techno-economic analysis with logistics and refinery cost surrogate models to evaluate multiscale miscanthus-derived biorefineries across the rainfed U.S. for the production of ethanol, succinic acid, lactic acid, potassium sorbate, and acrylic acid. Overall costs change little as feedstock density increases, while transport distances decrease by ∼30 to 67% (∼100 km) and siting flexibility improves. Specifically, a 5-fold feedstock density increase (2% to 10% of suitable land) reduces minimum selling prices by <10% (e.g., 0.27 USD·gal –1 for ethanol). This limited economic sensitivity suggests dense planting is not required for competitive deployment, particularly for smaller-scale facilities. Representing collection areas as irregular rather than circular expands the feasible space under low-density scenarios. While large-scale refineries anchor regional supply chains, smaller facilities retain spatial flexibility even when large refineries are established. These findings highlight the importance of spatial representation and multiscale coordination for robust, regionally tailored biomanufacturing networks to advance renewable carbon integration without extensive land conversion.

biorefinery siting↗

Stochastic nonlinear analysis of unidirectional fiber composites using image-based microstructural uncertainty quantification

Here, we present a data-driven nonlinear uncertainty quantification and propagation framework to study the microstructure-induced stochastic performance of unidirectional (UD) carbon fiber reinforced polymer (CFRP) composites. The proposed approach integrates (1) microscopic image characterization, (2) stochastic microstructure reconstruction, and (3) efficient multiscale finite element simulations enabled by self-consistent clustering (SCA) analysis. To model the complex microstructural variability, the proposed UQ methods take the non-Gaussian uncertainty sources into account through a distribution-free sampling approach leveraging nonparametric and asymptotic statistical tools. A hierarchical conditional sampling strategy enables the simultaneous sampling of multiple sources of uncertainties. Our approach provides insights into the impact of microstructural variabilities, which are shown to have an increasing impact on the nonlinear responses of UD CFRP parts under progressive compression loading and ultimately on the failure rate over time. We discover that before CFRP parts start to fail, a characteristic time period emerges with distinctive uncertainty distributions specific to the microstructure variability and the probability of failure. Identifying the failure time period is crucial to the reliability prediction, which is an essential component of CFRP design.

36 MATERIALS SCIENCE↗

Slow Electron Holes in the Earth's Magnetosheath

Here, we present a statistical analysis of electrostatic solitary waves observed aboard Magnetospheric Multiscale spacecraft in the Earth's magnetosheath. Applying single-spacecraft interferometry to several hundred solitary waves collected in about 2-minute interval, we show that almost all of them have the electrostatic potential of positive polarity and propagate quasi-parallel to the local magnetic field with plasma frame velocities of the order of 100 km/s. The solitary waves have typical parallel half-widths from 10 to 100 m that is between 1 and 10 Debye lengths and typical amplitudes of the electrostatic potential from 10 to 200 mV that is between 0.01% and 1% of local electron temperature. The solitary waves are associated with quasi-Maxwellian ion velocity distribution functions, and their plasma frame velocities are comparable with ion thermal speed and well below electron thermal speed. We argue that the solitary waves of positive polarity are slow electron holes and estimate the time scale of their acceleration, which occurs due to interaction with ions, to be of the order of one second. The observation of slow electron holes indicates that their lifetime was shorter than the acceleration time scale. We argue that multi-spacecraft interferometry applied previously to these solitary waves is not applicable because of their too-short spatial scales. The source of the slow electron holes and the role in electron-ion energy exchange remain to be established.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiscale petrophysical characterization and flow unit classification of the Minnelusa eolian sandstones

Integration of petrophysical and geological information is critical to simulation of subsurface carbon storage (GCS). In this sense, two depositional facies were identified from the core description and well-log interpretation, namely massive (MS) and cross-bedded (CB) facies groups. Additionally, pore-scale characteristics were studied by a combination of techniques, e.g. Nuclear Magnetic Resonance (NMR) and mercury intrusion capillary pressure (MICP). Scanning electron microscope (SEM) and petrographic analyses show that the pore structure is dominantly controlled by the depositional environment and dolomite cementation. NMR-T 2 distributions of MS and CB facies show triple and quadruple modes, respectively. In addition, MICP of high- and low-permeability MS facies samples, and their CB facies group mixtures were collected. The MS sample pore-throat size distribution is uni-modal, while the triple-modal characteristic of the mixtures indicates heterogeneous pore structures at the sub-core scale for CB facies. The reliably estimates of porosity and permeability for both facies groups via NMR techniques and the MLR (Multiple Linear Regression) approach demonstrate the applicability of these techniques to eolian sandstone. Moreover, irreducible water saturation via the T 2 -cutoff method correlates strongly with T 2LM instead of porosity. Finally, the rock quality index and flow zone indicator were calculated based on Combinable Magnetic Resonance (CMR) log interpretations. This provides direct connection to properties measured in the well. Four flow units were classified for both facies groups. Finally, results show that better reservoir quality with significant heterogeneities is observed in the CB facies. This study highlights the importance integrating multiscale petrophysical properties including facies, pore architecture and diagenesis analysis with core- to log-scale property characterization. The results herein validate our reservoir characterization and flow unit classification in eolian reservoirs.

58 GEOSCIENCES↗

Multiscale and multiphysics FEA simulation and materials optimization for laser ultrasound transducers

In this study, the relationship between the nanocomposite design and the laser ultrasound transducer (LUT) characteristics was investigated through simulations in multiple scale levels for material behavior, device response, and acoustic wave propagation in media. First, the effects of the nanoparticle size and concentration on the effective properties of composites were quantitatively investigated with the finite element analysis (FEA) method. Second, the effective properties of the nanocomposite were assigned to the layer, which is modeled as a homogeneous material, in the FEA for the LUT simulating the energy conversion from the incident laser to the acoustic wave. Finally, the ultrasound propagation in the water was calculated by a theoretical wave propagation model. The FEA-based prediction was compared with the experimental data in the literature and a theoretical analysis for LUT based on Thermal-Acoustic coupling. As a result, the ultrasound waves on the transducer surface and at a distance in the water could be predicted. Based on the hierarchically integrated prediction procedure, the optimal conditions of the photoacoustic nanocomposites were investigated through the parametric study with the particle size and concentration as variables. The results guide the material designs optimized for different device characteristics, such as high pressure and broad bandwidth.

36 MATERIALS SCIENCE↗

Tritium Transport Modeling at the Pore Scale in Ceramic Breeder Materials Using TMAP8

Fusion reactors depend on the blanket material to breed and release tritium at the same rate or faster than it is consumed by the fusion reaction. Cellular ceramic breeders (CCBs) are dense materials that can maintain a high tritium breeding ratio while promoting tritium release because of highly connected pores. Assessing the tritium breeding capabilities of these materials requires a combination of extensive experimental and modeling efforts. In this work, we develop and calibrate a multiphysics model of tritium transport. This novel model accounts for ceramic and pore diffusion, trapping and detrapping, and several surface reactions at the pore surface. We perform a sensitivity analysis and calibrate the model by comparing its predictions against experimental measurements of deuterium absorption. The calibrated model is then used to model tritium absorption in samples with different pore microstructures to investigate the effect of pore interconnectivity on tritium absorption. The model is part of the development of the multiscale, multiphysics framework for tritium transport [i.e., the Tritium Migration Analysis Program (TMAP8)], which is itself built on top of the finite-element multiphysics framework multiphysics object-oriented simulation environment (MOOSE). This study demonstrates some of TMAP8’s capabilities and is the first step toward assessing the tritium breeding capabilities of ceramic breeder material designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TRITIUM TRANSPORT MODELING AT THE PORE SCALE IN CERAMIC BEEDER MATERIALS USING TMAP8

Fusion reactors depend on the blanket material to breed and release tritium at the same rate or faster than it is consumed by the fusion reaction. Cellular ceramic breeders are dense materials that can maintain a high tritium breeding ratio while promoting tritium release thanks to highly connected pores. Assessing the tritium breeding capabilities of these materials requires a combination of extensive experimental and modeling efforts. In this work, we develop and calibrate a multiphysics model of tritium transport. This novel model accounts for ceramic and pore diffusion, trapping and detrapping, and several surface reactions at the pore surface. We perform a sensitivity analysis and calibrate the model by comparing its predictions against experimental measurements of deuterium absorption. The calibrated model is then used to model tritium absorption in samples with different pore microstructures to investigate the effect of pore interconnectivity on tritium absorption. The model is part of the development of the multiscale, multiphysics framework for tritium transport (i.e., the Tritium Migration Analysis Program [TMAP8]), which is itself built on top of the finite element multiphysics framework MOOSE. This study demonstrates some of TMAP8’s capabilities, and is the first step toward assessing the tritium breeding capabilities of ceramic breeder material designs.

08 HYDROGEN↗

Kinetic characterization of 3D magnetic reconnection: A transformative step (Final Report)

We investigated magnetic reconnection in the Earth’s magnetopause, magnetotail, and the bow shock, by means of fully kinetic simulations and analysis of space data obtained by NASA’s Magnetospheric Multiscale mission. We elucidated signatures of particle energization in particle distribution functions in the vicinity of reconnection X-line, properties of reconnection in shock-driven turbulence, and wave excitations and nonlinear structures. We advanced the understanding of magnetic reconnection from the viewpoint of particle kinetics in the diffusion region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Dependence of Electron-only Magnetic Reconnection on Plasma Beta and Magnetic Shear: Theory and Observations

We investigate the onset conditions for electron-only magnetic reconnection through an analysis of 86 reconnection events observed by the Magnetospheric Multiscale mission in magnetosheath and magnetospheric current sheets. Our study presents the first observational evidence that the occurrence of electron-only reconnection depends on both the magnetic-field shear angle (θ) across the current sheet and the cross-sheet electron beta jump, quantified either by the inflow electron beta difference (Δβ e ) or by a modified electron beta jump (Δβ$^*_e$) defined by the cross-sheet electron pressure difference normalized by the average of the inflow magnetic pressure in the two inflow regions. We extend the diamagnetic-drift suppression mechanism from standard reconnection to electron-only reconnection by deriving relations applicable to electron-scale current sheets. Using scaling estimates, the electron-only reconnection suppression relation is derived as Δβ$^*_e$ ≳ 2(L/d e ) tan (θ/2), while a natural extension of the ion-coupled reconnection suppression criterion is Δβ e ≳ 2(L/d e ) tan (θ/2). At low Δβ e and Δβ$^*_e$ values, reconnection is detected across a large range of magnetic shear (from low to high), whereas high Δβ e and Δβ$^*_e$ conditions restrict reconnection to large shear configurations: shear angles as small as θ ≈ 5° occur when Δβ e (or Δβ$^*_e$) ≲0.1. However, only events exceeding θ ≳ 70° are found when Δβ e ≳ 1.0; this apparent restriction weakens when employing Δβ$^*_e$, as the corresponding data points shift toward smaller values of the modified electron beta jump. These findings quantitatively support theoretical predictions of reconnection suppression in high-β e plasmas at low shears, attributed to X-line drift induced by pressure gradients across electron-scale current sheets.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiscale modeling of metal-hydride interphases—quantification of decoupled chemo-mechanical energies

Abstract The quantification of interphase properties between metals and their corresponding hydrides is crucial for modeling the thermodynamics and kinetics of the hydrogenation processes in solid-state hydrogen storage materials. In particular, interphase boundary energies assume a pivotal role in determining the kinetics of nucleation, growth, and coarsening of hydrides, alongside accompanying morphological evolution during hydrogenation. The total interphase energy arises from both chemical bonding and mechanical strains in these solid-state systems. Since these contributions are usually coupled, it is challenging to distinguish via conventional computational approaches. Here, a comprehensive atomistic modeling methodology is developed to decouple chemical and mechanical energy contributions using first-principles calculations, of which feasibility is demonstrated by quantifying chemical and elastic strain energies of key interfaces within the FeTi metal-hydride system. Derived materials parameters are then employed for mesoscopic micromechanical analysis, predicting crystallographic orientations in line with experimental observations. The multiscale approach outlined verifies the importance of the chemo-mechanical interplay in the morphological evolution of growing hydride phases, and can be generalized to investigate other systems. In addition, it can streamline the design of atomistic models for the quantitative evaluation of interphase properties between dissimilar phases and allow for efficient predictions of their preferred phase boundary orientations.

36 MATERIALS SCIENCE↗

A predictive modeling tool for damage analysis and design of hydrogen storage composite pressure vessels

In this study, a predictive modeling tool is developed for damage analysis and design of hydrogen (H 2 ) storage composite pressure vessels. It integrates micromechanics of matrix cracking into a continuum damage mechanics (CDM) description for damage evolution, and three-dimensional (3D) finite element (FE) modeling of the vessel structural response. At the scale of the composite layer (mesoscale), the temperature-dependent stiffness reduction law in terms of the damage variable for transverse matrix cracking is computed using an Eshelby-Mori-Tanaka approach for the initial composite thermoelastic properties and a self-consistent model for the stiffness reduction as a function of the damage variable. While transverse matrix cracking obeying a damage evolution relation can progressively evolve from an initiation to a saturation state, fiber failure is predicted by a micromechanical fiber rupture criterion that accounts for the fiber strength and matrix stress. The implementation of this integrated multiscale modeling model into a 3D FE formulation enables damage analysis and design of H 2 storage composite pressure vessels. The developed tool is illustrated through 3D damage analyses of a cryogenically compressed H 2 storage vessel model subjected to thermomechanical loadings to investigate effects of the helical layer fiber orientation and loading scenario on damage development, vessel integrity and burst pressure.

08 HYDROGEN↗