Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale models”

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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Coupling of CTF and RELAP5-3D Within an Enhanced Fidelity Nuclear Power Plant Simulator

A robust and accurate multiphysics engineering simulator is being developed to model the core behavior and system response of pressurized water reactors. This simulator relies on the NESTLE and CTF computer codes to model the neutronics and thermal hydraulics (TH), respectively, inside the core on a nodal scale and on the Reactor Excursion and Leak Analysis Program—Three Dimensional (RELAP5-3D) to model the entire nuclear steam supply system. The RELAP5-3D model includes highly detailed nodalization and multidimensional flow modeling throughout the vessel. Previously, pin-resolved data generated via the Virtual Environment for Reactor Analysis core simulator were used to improve the accuracy of the NESTLE core predictions. The engineering simulator being developed as part of this work uses the 3KEYMASTER platform to couple the enhanced NESTLE model to a nodal-fidelity CTF model to balance run time with accuracy; NESTLE provides node-dependent powers to CTF, and CTF provides node-dependent coolant densities and fuel temperatures to NESTLE.An overlapping domain approach is used for the core TH in which RELAP5-3D provides core boundary conditions based on the system response and CTF provides a node-dependent coolant heating rate to the RELAP5-3D core solution. In the preliminary TH demonstration discussed in this paper, CTF and RELAP5-3D provided similar steady-state core predictions, indicating the hydraulic compatibility between the codes, as well as reasonable and expected behavior under hypothetical transient conditions. This provides an initial step in ongoing efforts toward a robust, multiscale TH/neutronics engineering simulator capability.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2006

Multiscale Analysis of Delamination of Carbon Fiber-Epoxy Laminates with Carbon Nanotubes

A multi-scale analysis is presented to parametrically describe the Mode I delamination of a carbon fiber/epoxy laminate. In the midplane of the laminate, carbon nanotubes are included for the purposes of selectively enhancing the fracture toughness of the laminate. To analyze carbon fiber epoxy carbon nanotube laminate, the multi-scale methodology presented here links a series of parameterizations taken at various length scales ranging from the atomistic through the micromechanical to the structural level. At the atomistic scale molecular dynamics simulations are performed in conjunction with an equivalent continuum approach to develop constitutive properties for representative volume elements of the molecular structure of components of the laminate. The molecular-level constitutive results are then used in the Mori-Tanaka micromechanics to develop bulk properties for the epoxy-carbon nanotube matrix system. In order to demonstrate a possible application of this multi-scale methodology, a double cantilever beam specimen is modeled. An existing analysis is employed which uses discrete springs to model the fiber bridging affect during delamination propagation. In the absence of empirical data or a damage mechanics model describing the effect of CNTs on fracture toughness, several tractions laws are postulated, linking CNT volume fraction to fiber bridging in a DCB specimen. Results from this demonstration are presented in terms of DCB specimen load-displacement responses.

Riddick, Jaret C.↗
DOE OSTI2020

Multi-scale model based design of membrane reactor/separator processes for intensified hydrogen production through the water gas shift reaction

This paper aims to first quantify the impact of various diffusion models (Maxwell-Stefan, Wilke, Dusty-Gas) on the predictions of a multi-scale membrane reactor/separator mathematical model, and to then demonstrate this model's use for the design and process intensification of membrane reactor/separator systems for hydrogen production. This multi-scale model captures velocity, temperature and species' concentration profiles along the catalyst pellet's radial direction, and along the reactor's axial direction, by solving the momentum, energy, and species transport equations, accounting for convection, conduction, reaction, and diffusion mechanisms. In the first part of work, the effect of pellet-scale design parameters (mean pore diameter, volumetric porosity, tortuosity factor, etc.) and various species' flux models on the model predictions is studied. In the second part, the study focuses on the comparison, in terms of their process intensification characteristics, of various hydrogen production processes. These include a conventional high-temperature shift reactor (HTSR)/low-temperature shift reactor (LTSR) sequence, a novel HTSR/membrane separator (MS)/LTSR/MS sequence, and a process that involves low-temperature shift membrane reactors-LTSMR in a series.

08 HYDROGEN↗
DOE OSTI2021

Persistent, “Mysterious” Seismoacoustic Signals Reported in Oklahoma State during 2019

Here, we report on the source of seismoacoustic pulses that were observed across the state of Oklahoma (OK) during summer of 2019, and the subject of national media coverage and speculation. Seismic network data collected across four U.S. states and interviews with witnesses to the pulse’s effect on residential structures demonstrate that they were triggered by routine ammunition disposal operations conducted by McAlester Army Ammunition Plant (McAAP). During these operations, conventional explosives destroy obsolete munitions stored in pits through a controlled sequence of electronically timed shots that occur over tens of minutes. Despite noise-abatement efforts that reduce coupling of acoustic energy with air, some lower frequency, subaudible (infrasonic) sound radiates from these shots as discrete pulses. We use nine months of blast log documents, seismic network records, analyst picks, and physical modeling to demonstrate that seismic stations as far as 640 km from McAAP sample these pulses, which record seasonal patterns in stratospheric and tropospheric winds, as well as the dynamic formation of waveguides and shadow zones. Digital short-term average to long-term average detectors that we augment with dynamic thresholds and time-binning operations identify these pulses with a fair probability, when compared with visual observations. Our analyses thereby provide estimates of observation rates for both partial and full sequences of these pulses, as well as single shots. We suggest that disposal operations can exploit existing, composite seismic networks to predict where residents are likely to witness blasting. Crucially, our data also show that dense seismic networks can record multiscale atmospheric processes in the absence of infrasound arrays.

58 GEOSCIENCES↗
DOE OSTI2022

Clustering of inertial particles in turbulent flow through a porous unit cell

We report direct numerical simulation is used to investigate effects of turbulent flow in the confined geometry of a face-centred cubic porous unit cell on the transport, clustering and deposition of fine particles at different Stokes numbers (St = 0.01, 0.1, 0.5, 1, 2) and at a pore Reynolds number of 500. Particles are advanced using one-way coupling and the collision of particles with pore walls is modelled as perfectly elastic with specular reflection. Tools for studying inertial particle dynamics and clustering developed for homogeneous flows are adapted to take into account the embedded, curved geometry of the pore walls. The pattern and dynamics of clustering are investigated using the volume change of Voronoi tesselation in time to analyse the divergence and convergence of the particles. Similar to the case of homogeneous, isotropic turbulence, the cluster formation is present at large volumes, while cluster destruction is prominent at small volumes and these effects are amplified with the Stokes number. However, unlike homogeneous, isotropic turbulence, the formation of a large number of very small volumes was observed at all Stokes numbers and attributed to the collision of particles with the pore wall. Multiscale wavelet analysis of the particle number density indicates that the peak of the energy density spectrum, representative of enhanced particle clustering, shifts towards larger scales with an increase in the Stokes number. Scale-dependent skewness and flatness quantify the intermittent void and cluster distribution, with cluster formation observed at small scales for all Stokes numbers, and void regions at large scales for large Stokes numbers.

42 ENGINEERING↗
NASA NTRS2016

Energy Limits of Electron Acceleration in the Plasma Sheet During Substorms: A Case Study with the Magnetospheric Multiscale (MMS) Mission

We present multipoint observations of earthward moving dipolarization fronts and energetic particle injections from NASAs Magnetospheric Multiscale mission with a focus on electron acceleration. From a case study during a substorm on 02 August 2015, we find that electrons are only accelerated over a finite energy range, from a lower energy threshold at approx. 7-9 keV up to an upper energy cutoff in the hundreds of keV range. At energies lower than the threshold energy, electron fluxes decrease, potentially due to precipitation by strong parallel electrostatic wavefields or initial sources in the lobes. Electrons at energies higher than the threshold are accelerated cumulatively by a series of impulsive magnetic dipolarization events. This case demonstrates how the upper energy cutoff increases, in this case from approx. 130 keV to >500 keV, with each depolarization/injection during sustained activity. We also present a simple model accounting for these energy limits that reveals that electron energization is dominated by betatron acceleration.

Turner, D. L.↗
DOE OSTI2020

Computer Simulation of Proton Transport in Fuel Cell Membranes (Final Report)

This DOE-supported research grant focused on understanding the nature of proton transport in complex systems such as proton exchange membranes (PEMs). The most unique aspect of the research was the development and implementation of a novel multiscale reactive molecular dynamics (MS-RMD) methodology. In this approach, covalent bonds can dynamically break and form, allowing one to accurately treat the proton hopping process essential to capturing the physics of proton transport. The Voth group applied this method to proton exchange membrane systems, providing insight into their proton transport mechanism. They found that protons can diffuse most rapidly in the water-rich regions, but that protons actually spend so little time in such regions that transport along the hydrophobic – hydrophilic interface controls the membrane performance. In addition, the group worked to increase understanding of acidic solutions, developing methods for simulating and interpreting experimental infrared vibrational spectroscopy for excess protons (acidic solutions). The group also implemented novel tools for developing proton transport reactive MD models using a relative entropy minimization scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2021

The Dynamics of a High Mach Number Quasi-Perpendicular Shock: MMS Observations

Shock parameters at Earth’s bow shock, in rare instances, can approach the Mach numbers predicted at astrophysical shocks and supernova remnants. We present our analysis of a high Alfv ́en Mach number (MA= 27) shock, by utilizing multipoint measurements from the Magnetospheric Multiscale (MMS) spacecraft during a crossing of Earth’s quasi-perpendicular bow shock. We find that the shock dynamics are mostly driven by reflected ions, perturbations that they generate, and nonlinear amplification of the perturbations. Our analyses indicate that reflected ions create modest magnetic enhancements upstream of the shock front which evolve in a nonlinear manner as they traverse the shock foot. They can transform into proto-shocks that propagate at small angles to the magnetic field and towards the bow shock. The nonstationary bow shock shows signatures of both reformation and surface ripples. Our observations indicate that although shock reformation occurs, the main shock layer never disappears. These observations are at high plasmaβ, a parameter regime which has not been well explored by numerical models.

H Madanian↗
DOE OSTI2022

Machine learning approaches to streamline and enhance the analysis of multiscale imaging data for bioaerosol and soil particles

Bioaerosol and soil particles are ubiquitous in the environment. They are multicomponent and complex in nature displaying mixed inorganic and organic components. The way components are mixed in a bioaerosol sample is referred to as its mixing state. Soil particles are also a mixture of inorganic (mineral) and organic (soil organic matter) components. Bioaerosol particles contribute to a major fraction of coarse mode atmospheric particles, especially in the tropical areas, contributing up to 80 % of the particle mass concentration. The mixing state of particles is crucial to evaluate because it impacts several important environmental processes such as warm and cold cloud formation and radiation budget. Mixing states in aerosols are accompanied by chemical reactions across solid-liquid-gas interfaces. In this study, we utilized elemental compositions and microcopy images of thousands of atmospheric particles acquired by computer-controlled scanning electron microscope equipped with an energy-dispersive x-ray spectrometer to compute the mixing state of atmospheric particles. A 2D convolutional neural network (CNN), also known as convnet, was used to model the relationship between low resolution imaging data and higher resolution spectroscopy data, with the former as training input and the latter as target output. Two types of CNNs were implemented and tested; a basic CNN and an Inception-v3 network. For binary classification, the basic CNN achieved an accuracy of 84.29 % across all atom types, and the Inception-v3-like network achieved an accuracy 85.51 %. This study demonstrates the applicability of deep learning to handle large amounts of imaging/chemical spectroscopy data efficiently and evaluate particle mixing state from a range of environmental samples.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

The Lunar GNSS Receiver Experiment (LuGRE)

The Lunar GNSS Receiver Experiment (LuGRE) is a joint NASA-Italian Space Agency (ASI) payload on the Firefly Blue Ghost Mission 1 (BGM1) with the goal to demonstrate GNSS-based positioning, navigation, and timing at the Moon. LuGRE was chosen by the NASA Commercial Lunar Payload Services (CLPS) program as one of ten payloads on its “19D” task order for delivery to the lunar surface in 2023. The LuGRE payload consists of a weak-signal GNSS receiver, a high-gain L-band patch antenna, a low-noise amplifier, and an RF filter. The receiver will track GPS L1 C/A and L5, and Galileo E1 and E5a signals and will return pseudorange, carrier phase, and Doppler measurements to the ground. It will also calculate least-squares point solutions and Kalman-filter based navigation solutions onboard. In addition, the receiver features the capability to record raw I/Q baseband samples for downlink and ground processing. LuGRE will build on the legacy of prior missions in the Space Service Volume (SSV) including the initial experiments by AMSAT-OSCAR 40 and others, the GOES-R series of geostationary weather satellites, and the NASA Magnetospheric Multiscale (MMS) mission currently operating on GPS-based navigation at nearly 50% of lunar distance. Further, LuGRE will be one of the very first demonstrations of GNSS signal reception and navigation in the lunar environment and on the lunar surface, paving the way for operational use by future lunar missions such as Orion, Gateway, robotic and human landers, and surface rovers. Ultimately, all LuGRE science data will be released to a public data archive for the benefit of the GNSS and space communities. This paper provides a detailed overview of the LuGRE payload, including its design, concept of operations, and its predicted ability to meet its core science objectives. The baseline science investigations and priorities are outlined. Simulated performance results are shown based on the latest calibrated models including signal strength, signal availability, onboard navigation performance and convergence properties, and ground-based post-processed navigation performance.

LuGRE↗
DOE OSTI2026

Multiphysics simulation of TRISO fuel compacts and the effects of homogenization on silver release predictions

This work studies the impact of explicit and homogenized modeling approaches on the multiphysics simulation of TRistructural ISOtropic (TRISO) fuel compacts in prismatic High Temperature Gas Reactors (HTGRs) and the silver release predictions. TRISO fuel compacts exhibit complex double heterogeneity that significantly affects heat conduction, neutron transport, and fission products release. In this work, we use Cardinal, a multiphysics tool based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, to couple neutron transport and heat conduction. OpenMC is used for neutron transport and the MOOSE heat transfer module is used for the heat conduction. Then, we use BISON for the silver release predictions. Two different modeling approaches—explicit modeling of individual TRISO particles and homogenized representation using effective thermal properties—are compared at high TRISO packing fractions (20% and 40%) across varying power densities. Results demonstrate that homogenization significantly under-predicts peak temperatures, for the case of high power/TRISO, there is a difference of 83.73 K in the maximum temperature. Additionally, homogenization underestimates the silver release predictions compared to explicit modeling. Specifically, the temperature and power density differences lead to significant differences in silver release predictions. This work demonstrates the importance of accurately modeling heterogeneity of the TRISO particles to reliably predict fission product release and assess reactor safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Multi-GPU immersed boundary method hemodynamics simulations

Large-scale simulations of blood flow that resolve the 3D deformation of each comprising cell are increasingly popular owing to algorithmic developments in conjunction with advances in compute capability. Among different approaches for modeling cell-resolved hemodynamics, fluid structure interaction (FSI) algorithms based on the immersed boundary method are frequently employed for coupling separate solvers for the background fluid and the cells within one framework. GPUs can accelerate these simulations; however, both current pre-exascale and future exascale CPU-GPU heterogeneous systems face communication challenges critical to performance and scalability. In this paper, we describe, to our knowledge, the largest distributed GPU-accelerated FSI simulations of high hematocrit cell-resolved flows with over 17 million red blood cells. We compare scaling on a fat node system with six GPUs per node and on a system with a single GPU per node. Through comparison between the CPU- and GPU-based implementations, we identify the costs of data movement in multiscale multi-grid FSI simulations on heterogeneous systems and show it to be the greatest performance bottleneck on the GPU.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Toward Addressing the Challenge to Predict the Heat Capacities of RDX and HMX Energetic Materials

Availability of heat capacity as function of pressure and temperature is an essential prerequisite for development of a computational multiscale strategy capable to address the evolution of microstructure and energy release in advanced high energy density materials. In the case of 1,3,5-trinitro-1,3,5-triazinane (RDX) and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) systems as two of the most studied energetic materials, there are substantial gaps in experimental data, with available heat capacities values distributed only in a region close to standard ambient conditions. In this study we demonstrate how these major experimental limitations can be addressed in the case of the RDX and HMX systems based on the combined use of classical and quantum mechanical calculations. We show that by considering ideal gas properties evaluated using quantum mechanical methods, and residual properties obtained from molecular simulations using fully flexible atomistic force field models, excellent agreement can be obtained for the predicted heat capacities to the most recent experimental values. An important advantage of the current computational methodology is that it allows evaluation of both constant-volume and constant-pressure heat capacities for a broad interval of temperatures and pressures, which encompasses solid and liquid phases conditions. In the case of the solid α and γ phases of RDX and the β phase of HMX, the predicted results follow closely both the available experimental data at standard ambient conditions and the results obtained using density functional theory calculations at high pressures, a regime where experimental data are not available. A perspective to expand the current methodology is also discussed.

36 MATERIALS SCIENCE↗
NASA NTRS2015

High Frequency Design Considerations for the Large Detector Number and Small Form Factor Dual Electron Spectrometer of the Fast Plasma Investigation on NASA's Magnetospheric Multiscale Mission

Each half of the Dual Electron Spectrometer (DES) of the Fast Plasma Investigation (FPI) on NASA's Magnetospheric MultiScale (MMS) mission utilizes a microchannel plate Chevron stack feeding 16 separate detection channels each with a dedicated anode and amplifier/discriminator chip. The desire to detect events on a single channel with a temporal spacing of 100 ns and a fixed dead-time drove our decision to use an amplifier/discriminator with a very fast (GHz class) front end. Since the inherent frequency response of each pulse in the output of the DES microchannel plate system also has frequency components above a GHz, this produced a number of design constraints not normally expected in electronic systems operating at peak speeds of 10 MHz. Additional constraints are imposed by the geometry of the instrument requiring all 16 channels along with each anode and amplifier/discriminator to be packaged in a relatively small space. We developed an electrical model for board level interactions between the detector channels to allow us to design a board topology which gave us the best detection sensitivity and lowest channel to channel crosstalk. The amplifier/discriminator output was designed to prevent the outputs from one channel from producing triggers on the inputs of other channels. A number of Radio Frequency design techniques were then applied to prevent signals from other subsystems (e.g. the high voltage power supply, command and data handling board, and Ultraviolet stimulation for the MCP) from generating false events. These techniques enabled us to operate the board at its highest sensitivity when operated in isolation and at very high sensitivity when placed into the overall system.

Microchannel plates high frequency performance↗
DOE OSTI2024

Scalable Multiphysics Block Preconditioning for Low Mach Number Compressible Resistive MHD with Application to Magnetic Confinement Fusion

This study investigates multiphysics block preconditioners that are critical in devising scalable Newton–Krylov iterative solvers for longer time-scale fully implicit fluid plasma models. The specific model of interest is the visco-resistive, low Mach number, compressible magnetohydrodynamics (MHD) model. This model describes the dynamics of conducting fluids in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics. The specific application of interest that motivates this study is the macroscopic simulation of longer time-scale stability and disruptions of magnetic confinement fusion devices, specifically the ITER Tokamak. The computational solution of the governing balance equations for mass, momentum, heat transfer, and magnetic induction for resistive MHD systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena as well as the significant range of time and length scales that the interactions of these physical mechanisms produce. To handle the range of time and spatial scales of interest, a fully implicit unstructured variational multiscale finite element formulation is employed. For the scalable solution of the Newton linearized systems, fully coupled block preconditioners are designed to leverage algebraic multigrid subsolves. In conclusion, results are presented for the strong and weak scaling of the method as well as the robustness of these techniques for a large range of Lundquist numbers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Digital Twins for Materials

Digital twins are emerging as powerful tools for supporting innovation as well as optimizing the in-service performance of a broad range of complex physical machines, devices, and components. A digital twin is generally designed to provide accurate in-silico representation of the form (i.e., appearance) and the functional response of a specified (unique) physical twin. This paper offers a new perspective on how the emerging concept of digital twins could be applied to accelerate materials innovation efforts. Specifically, it is argued that the material itself can be considered as a highly complex multiscale physical system whose form (i.e., details of the material structure over a hierarchy of material length) and function (i.e., response to external stimuli typically characterized through suitably defined material properties) can be captured suitably in a digital twin. Accordingly, the digital twin can represent the evolution of structure, process, and performance of the material over time, with regard to both process history and in-service environment. This paper establishes the foundational concepts and frameworks needed to formulate and continuously update both the form and function of the digital twin of a selected material physical twin. The form of the proposed material digital twin can be captured effectively using the broadly applicable framework of n-point spatial correlations, while its function at the different length scales can be captured using homogenization and localization process-structure-property surrogate models calibrated to collections of available experimental and physics-based simulation data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Tracking ion intercalation into layered Ti 3 C 2 MXene films across length scales

Enhancing the energy stored and power delivered by layered materials relies strongly on improved understanding of the intricate interplay of electrolyte ions, solvents, and electrode interactions as well as the role of confinement. Here we report a highly integrated study with multiscale theory/modelling and experiments to track the intercalation of aqueous Li + , Na + , K + , Cs + , and Mg 2+ ions into Ti 3 C 2 MXene. The integrated analysis of experiments assisted by theory/modelling allows for a deep understanding of energy storage processes highlighting the importance of the dynamics of cations, their positionings between MXene sheets, their effects on mechanical properties and capacitive energy storage. Computational simulations and operando calorimetry measurements prove the processes involving cation dehydration and H+ rehydration, showing a good correlation for heat variations between experiments and theory. Operando liquid AFM mapped energy dissipation of ions appears non-uniformly across the MXene surface, indicating heterogeneities of ions inside the MXene and confirming partially the ion behaviour obtained in theory. We directly demonstrate that the average distance between the cation and MXene surface follows a modified two-sided Helmholtz model when plotted versus the open circuit potential capacitance, revealing a different electrical double layer mechanism in confinement. This new fundamental understanding lays the foundation for improved functional devices utilizing electrodes and membranes made of two-dimensional materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Improving High Resolution Offshore Wind Resource Assessments and Forecasts using Observations in the MA/RI Lease Areas

The third Wind Forecast Improvement Project (WFIP3) sought to improve understanding of the physical phenomena in the atmosphere and ocean that dictate the structure and variability of wind and thermodynamic fields within the Marine Atmospheric Boundary Layer (MABL). WFIP3 conducted a comprehensive 18-month observational study over the Northeast U.S. outer continental shelf, a high use coastal zone, using a 3-D multiscale sensor array to highly resolve the temporal, vertical, and horizontal structure of the coupled atmospheric and oceanic boundary layers. Multiple land-based study sites adjacent to the coastal ocean observed surface meteorology and vertical profiles of atmospheric properties via passive infrared and microwave radiometers, active lidars and radars, and radiosondes. At sea, an array of surface flux buoys and two vertical profiling lidar buoys observed both atmospheric and oceanic properties, augmented by land-based oceanographic radar systems and routine ship-based surveys. Intensive observations of the MABL over the ocean was done from an air-sea interaction flux tower and extended deployments of a large barge platform. WFIP3 focused on mesoscale and sub-mesoscale flows -- including sea breezes, low-level jets, low-level clouds, and coastal storms -- and the ability of advanced numerical model parameterizations to represent them within fully coupled oceanic and atmospheric modeling systems and foundational weather forecast models. Numerous critical forecasting phenomena were observed, however the project was terminated prior to the completion of the field observational period and the analysis period.

54 ENVIRONMENTAL SCIENCES↗