Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiscale Modeling and Simulation”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

PAB3D: Its History in the Use of Turbulence Models in the Simulation of Jet and Nozzle Flows

This is a review paper for PAB3D s history in the implementation of turbulence models for simulating jet and nozzle flows. We describe different turbulence models used in the simulation of subsonic and supersonic jet and nozzle flows. The time-averaged simulations use modified linear or nonlinear two-equation models to account for supersonic flow as well as high temperature mixing. Two multiscale-type turbulence models are used for unsteady flow simulations. These models require modifications to the Reynolds Averaged Navier-Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two-equation (k-epsilon) model with a RANS/LES transition function, dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier-Stokes (PANS) formulation. All of these models are implemented in the three-dimensional Navier-Stokes code PAB3D. This paper discusses computational methods, code implementation, computed results for a wide range of nozzle configurations at various operating conditions, and comparisons with available experimental data. Very good agreement is shown between the numerical solutions and available experimental data over a wide range of operating conditions.

Abdol-Hamid, Khaled S.↗
DOE OSTI2023

Multiscale Reactive Model for 1,3,5-Triamino-2,4,6-trinitrobenzene Inferred by Reactive MD Simulations and Unsupervised Learning

When high-energy-density materials are subjected to thermal or mechanical insults at extreme conditions (shock loading), a coupled response between the thermo-mechanical and chemical behaviors is systematically induced. Herein we develop a reaction model for the fast chemistry of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) at the mesoscopic scale, where the chemical behavior is determined by underlying microscopic reactive simulations. The slow carbon cluster formation is not discussed in the present work. All-atom reactive molecular dynamics (MD) simulations are performed with the ReaxFF potential, and a reduced-order chemical kinetics model for TATB is fitted to isothermal and adiabatic simulations of single crystal chemical decomposition. Unsupervised machine learning techniques based on non-negative matrix factorization are applied to MD trajectories to model the decomposition kinetics of TATB in terms of a four-component model. The associated heats of reaction are fit to the temperature evolution from adiabatic decomposition trajectories. Using a chemical species analysis, we show that non-negative matrix factorization captures the main chemical decomposition steps of TATB and provides an accurate estimation of their evolution with temperature. The final analytical formulation, coupled to a diffusion term, is incorporated into a continuum formalism, and simulation results are compared one-to-one against MD simulations of 1D reaction propagation along different crystallographic directions and with different initial temperatures. A good agreement is found for both the temporal and spatial evolution of the temperature field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Numerical Simulation of Container Breach and Airborne Release of Solids Due to Mechanical Insults

Throughout U.S. Department of Energy (DOE) complexes, safety engineers employ the five-factor formula to calculate the source term (ST) that includes parameters of airborne release fraction (ARF), respirable fraction (RF) and damage ratio (DR). Limited experimental data on fragmentation of solids, such as ceramic pellets (i.e., PuO2), and container breach due to mechanical insults (i.e., drop and forklift impact), can be supplemented by modeling and simulation using high fidelity computational tools to estimate these parameters. This paper presents the use of Sandia National Laboratories' SIERRA solid mechanics (SM) finite element code to investigate the behavior of the widely utilized waste container (such as 7A Drum) subject to a range of free fall impact and puncture scenarios. The resulting behavior of the container is assessed, and the estimates are presented for bounding DRs from calculated breach areas for the various accident conditions considered. This paper also describes a novel multiscale constitutive model recently implemented in SIERRA/SM that simulates the fracture of brittle materials such as PuO 2 and determines ARF during the fracture process. Furthermore, comparisons are made between model predictions and simple bench-top experiments.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2024

Leveraging MARVEL and SPHERE to Demonstrate NEAMS Thermal Hydraulics Codes

The Single Primary Heat Extraction and Removal Emulator (SPHERE) and Microreactor Applications Research Validation and Evaluation (MARVEL) were utilized to demonstrate thermal hydraulics codes of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. A recent long-duration test of a sodium heat pipe conducted at the SPHERE facility was utilized for the validation of the heat pipe code Sockeye. The first 10 hours of operation featured a startup of the heat pipe from room temperature; this portion was modeled with Sockeye and compared to the experimental data. Measurements included several thermocouple measurements across the length of the heat pipe surface. Two of Sockeye’s heat pipe models were applied to experiment, as well as several models of the cooling at the condenser end of the heat pipe. Sockeye results gave a decent match to the experimental data, with notable discrepancy in the condenser section, where experimental data indicated a significant inactive length of the heat pipe. The MARVEL microreactor was modeled through a coupled simulation of System Analysis Module (SAM) and SubChannel Module (SCM). This multiscale coupled simulation was performed as an exercise to demonstrate the compatibility and applicability of these codes for thermal-hydraulic analysis of microreactors. The coupling between the SCM and SAM codes for multiscale modeling of MARVEL, was achieved through a domain overlapping approach. In this example a transient SAM simulation provided boundary conditions to SCM, and SCM provided pressure drop to SAM.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Development and Validation of MALAMUTE model for Electric Field Assisted Sintering of Structural Materials

Fusion power plant designs feature extreme material performance requirements for structural material candidates. In addition to conventional alloys, more advanced composites and oxide dispersion strengthened (ODS) alloys are being explored, however, achieving the desired microstructures to maximize performance using traditional manufacturing methods can be challenging. The advanced manufacturing (AM) electric field-assisted sintering (EFAS) technique offers improved control over the final microstructure through higher heating and cooling rates and moderate pressures. Modeling and simulation tools show promise in elucidating the process-structure-property-performance (PSPP) correlation for AM-produced parts, including the EFAS process. An inherently multiscale process, the EFAS technique aligns well with the multiscale modeling capability of the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE)[cite]. We present here an electro-thermo-mechanical approach to modeling the EFAS process using the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) code. Prediction of the field and gradient distributions across the EFAS tooling is required to accurately describe the conditions for the lower-scale microstructural evolution models. In this work we present the MALAMUTE model developed to predict the electrical potential, temperature, and mechanical stress distribution across the EFAS graphite tooling and part at the larger engineering-scale. Validation of the MALAMUTE engineering-scale model is completed using data from experimental densification and pre-densified runs of iron powder via EFAS at 1000oC. These runs were conducted using a Thermal Technology DCS-5 EFAS system. Data collected during the experiment runs include the direct current (DC) supplied to the graphite tooling, the temperature of the graphite tooling as measured with a pyrometer, and the force applied to the top of the graphite tooling stack, and the data were recorded every 10 seconds. Our validation approach used the current and force data from the EFAS run as boundary condition inputs to the MAMALUTE simulation; the temperature data were used to evaluate the MALAMUTE EFAS model prediction. Results of the MALAMUTE simulations are employed to connect the external pyrometer temperature measurement to the temperature profile across the part undergoing consolidation. We investigate the impact of material property variation and mesh deformation on the temperature profile as predicted by MALAMUTE. We conclude by highlighting projects where the MALAMUTE EFAS modeling and simulation capabilities will be used to assist experimental design.

36 - MATERIALS SCIENCE↗
DOE OSTI2021

Emergent Recursive Multiscale Interaction in Complex Systems

This project studied the potential for multiscale group dynamics in complex social systems, including emergent recursive interaction. Current social theory on group formation and interaction focuses on a single scale (individuals forming groups) and is largely qualitative in its explanation of mechanisms. We combined theory, modeling, and data analysis to find evidence that these multiscale phenomena exist, and to investigate their potential consequences and develop predictive capabilities. In this report, we discuss the results of data analysis showing that some group dynamics theory holds at multiple scales. We introduce a new theory on communicative vibration that uses social network dynamics to predict group life cycle events. We discuss a model of behavioral responses to the COVID-19 pandemic that incorporates influence and social pressures. Finally, we discuss a set of modeling techniques that can be used to simulate multiscale group phenomena.

42 ENGINEERING↗
DOE OSTI2022

Assessing the Impact of Lightning NOx Emissions in CMAQ Using Lightning Flash Data from WWLLN over the Contiguous United States

Comparison of lightning flash data from the National Lightning Detection Network (NLDN) and from the World Wide Lightning Location Network (WWLLN) over the contiguous United States (CONUS) for the 2016–2018 period reveals temporally and spatially varying flash rates that would influence lightning NO x (LNO x ) production due to known detection efficiency differences especially during summer months over land (versus over ocean). However, the lightning flash density differences between the two networks show persistent seasonal patterns over geographical regions (e.g., land versus ocean). Since the NLDN data are considered to have higher accuracy (lightning detection with >95% efficiency), we developed scaling factors for the WWLLN flash data based on the ratios of WWLLN to NLDN flash data over time (months of year) and space. In this study, sensitivity simulations using the Community Multiscale Air Quality (CMAQ) model are performed utilizing the original data sets (both NLDN and WWLLN) and the scaled WWLLN flash data for LNO x production over the CONUS. The model performance of using the different lightning flash datasets for ambient O 3 and NO x mixing ratios that are directly impacted by LNO x emissions and the wet and dry deposition of oxidized nitrogen species that are indirectly impacted by LNO x emissions is assessed based on comparisons with ground-based observations, vertical profile measurements, and satellite products. During summer months, the original WWLLN data produced less LNO x emissions (due to its lower lightning detection efficiency) compared to the NLDN data, which resulted in less improvement in model performance than the simulation using NLDN data as compared to the simulation without any LNO x emissions. However, the scaled WWLLN data produced LNO x estimates and model performance comparable with the NLDN data, suggesting that scaled WWLLN may be used as a substitute for the NLDN data to provide LNO x estimates in air quality models when the NLDN data are not available (e.g., due to prohibitive cost or lack of spatial coverage).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

SO2 over Central China: Measurements, Numerical Simulations and the Tropospheric Sulfur Budget

SO2 in central China was measured in situ from an aircraft and remotely using the Ozone Monitoring Instrument (OMI) from the Aura satellite; results were used to develop a numerical tool for evaluating the tropospheric sulfur budget - sources, sinks, transformation and transport. In April 2008, measured ambient SO2 concentrations decreased from approx.7 ppbv near the surface to approx. 1 ppbv at 1800 m altitude (an effective scale height of approx.800 m), but distinct SO2 plumes were observed between 1800 and 4500 m, the aircraft's ceiling. These free tropospheric plumes play a major role in the export of SO2 and in the accuracy of OMI retrievals. The mean SO2 column contents from aircraft measurements (0.73 DU, Dobson Units) and operational OMI SO2 products (0.63+/-0.26 DU) were close. The OMI retrievals were well correlated with in situ measurements (r = 0.84), but showed low bias (slope = 0.54). A new OMI retrieval algorithm was tested and showed improved agreement and bias (r = 0.87, slope = 0.86). The Community Multiscale Air Quality (CMAQ) model was used to simulate sulfur chemistry, exhibiting reasonable agreement (r = 0.62, slope = 1.33) with in situ SO2 columns. The mean CMAQ SO2 loading over central and eastern China was 54 kT, approx.30% more than the estimate from OMI SO2 products, 42 kT. These numerical simulations, constrained by observations, indicate that ",50% (35 to 61 %) of the anthropogenic sulfur emissions were transported downwind, and the overall lifetime of tropospheric SO2 was 38+/-7 h.

He, Hao↗
DOE OSTI2021

Unraveling the Ion Adsorption Kinetics in Microporous Carbon Electrodes: A Multiscale Quantum-Continuum Simulation and Experimental Approach

With the understanding that sorption in porous carbon electrodes is crucial to many environmental and energy technologies, such as capacitive deionization (CDI), supercapacitor energy storage, and activated carbon filters. In each of these examples, a practical model that can describe ion electrosorption kinetics is highly desirable for accelerating material design. Here, we proposed a multiscale model to study the ion electrosorption kinetics in porous carbon electrodes by combining quantum mechanical simulations with continuum approaches. Our model integrates the Butler-Volmer (BV) equation for sorption kinetics and a continuously stirred tank reactor (CSTR) formulation with atomistic calculations of ion hydration and ion-pore interactions based on density functional theory (DFT). We validated our model experimentally by using ion mixtures in a flow-through electrode CDI device and developed an in-line UV absorption system to provide unprecedented resolution of individual ions in the separation process. We showed that the multiscale model captures unexpected experimental phenomena that cannot be explained by the traditional ion electrosorption theory. The proposed multiscale framework provides a viable approach for modeling separation processes in systems where pore sizes and ion hydration effects strongly influence the sorption kinetics, which can be leveraged to explore possible strategies for improving carbon-based and, more broadly, pore-based technologies

36 MATERIALS SCIENCE↗
DOE OSTI2023

Variational quantum simulation of the critical Ising model with symmetry averaging

Here we investigate the use of deep multiscale entanglement renormalization ansatz (DMERA) circuits as a variational ansatz. We use the exactly solvable one-dimensional critical transverse-field Ising model as a test bed. Numerically exact simulation of the quantum circuit ansatz can in this case be carried out to hundreds of qubits by exploiting efficient classical algorithms for simulating matchgate circuits. We find that, for this system, the DMERA strongly outperforms a standard quantum approximate optimization algorithm (QAOA)–style ansatz, and that a major source of systematic error in correlation functions approximated using the DMERA is the breaking of the translational and Kramers-Wannier symmetries of the transverse-field Ising model. We are able to reduce this error by up to four orders of magnitude by symmetry averaging, without incurring additional cost in qubits or circuit depth. Here, we propose that this technique for mitigating systematic error could be applied to noisy intermediate-scale quantum (NISQ) simulations of physical systems with other symmetries.

1-dimensional spin chains↗
DOE OSTI2019

Effect of Diffusion Length in Modeling of Equiaxed Dendritic Solidification under Buoyancy Flow in a Configuration of Hebditch–Hunt Experiment

Modeling of equiaxed solidification is vital for understanding the solidification process of metallic alloys. In this work, an extended literature review is given for the models currently used for equiaxed solidification simulations. Based on this analysis, we present a three-phase multiscale equiaxed solidification model in which some approximations regarding solute transport at microscopic scale are put together in a new way and incorporated into macroscopic transport equations. For the latter, a term relating to the momentum exchange between the two phases is revised, and a modification for the grain packing algorithm is proposed. A modernized model for equiaxed dendrite growth is tested using a case of solidification of Sn-5 wt pct Pb alloy in a parallelepiped cavity that mimics the Hebditch–Hunt experiment. The results obtained using two approaches to calculate diffusion length are presented and compared both with each other and with numerical results from elsewhere. It is demonstrated that diffusion length has a crucial effect on the final segregation pattern.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A review of microwave-assisted process intensified multiphase reactors

Microwaves provide alternative heating and allow process intensification due to their rapid, volumetric, and selective nature. Recognizing the central role of multiphase reactors in chemical industry, a recent surge in employing microwaves is observed. We review the recent experimental and modeling investigations of microwave heating of multiphase reactors with emphasis on chemical engineering applications. Here we demonstrate that there is accumulated evidence for improved performance via microwave heating and a clear opportunity for further process intensification. In most of the cases, this improved performance stems from a temperature gradient between two phases. We discuss the ongoing debate on the mechanism by which microwaves affect chemical processes exacerbated by the inability of measuring the temperature distribution in a microwave cavity. We outline recent progress in this direction in monolith reactors and needs for future work. We underscore the lack of detailed modeling and simulation tools even for single-phase systems and emphasize the imperative for multiscale predictive modeling to bridge the experimental-modeling gap. Promising results are shown by a few recently published modeling studies that can predict the experimental measurements in complex multiphase reactors. A combination of experimental and modeling tools can provide a comprehensive picture of the microwave multiphase reactors as well as a means toward scale-up and optimization.

42 ENGINEERING↗
DOE OSTI2025

Binding and Translocation of Substrate Allosterically Promotes Functional Interactions Within the AlkB–AlkG Electron Transfer Complex

The alkane monooxygenase AlkB and rubredoxin AlkG form an electron transfer complex that hydroxylates terminal alkanes to produce alcohols. The recent cryoEM study of Fontimonas thermophila AlkB-AlkG complex revealed its architecture, including a dodecane (D12) substrate at the active site. However, FtAlkBG molecular mechanism of action of remains unknown. Here, in this study, we examined its dynamics and interactions by multiscale computations, including molecular dynamics simulations, elastic network models, and QM/MM of the oxygen activation mechanism at the AlkB catalytic site. D12 maintained stable interactions within the catalytic site during two MD runs, coordinated by hydrophobic residues L263-L264, I267, I133. A third extended run revealed that D12 could translocate to a membrane-exposed site near S49/F46 along a hydrophobic channel gated by I54. During this translocation, D12 was temporarily stabilized at intermediate sites IS1 (lined by I27/L30-G31/G50/L53-I54/P59/S124/A127-V128) and IS2 (I33-G34/L37/L45-F46/S49) before nearly exiting the protein, and diffused back to the active site, assisted by L30. Substrate binding and translocation across those intermediate sites affects the coupling between the iron centers in AlkBG, and interfacial interactions between AlkB-AlkG. The channel was further connected to the cytosol, near two surface-exposed arginines, potentially allowing for O 2 passage. The allosteric effects between D12 putative entry site, catalytic site and AlkB-AlkG interface were analyzed by ENM-based methods which confirmed the cooperative perturbation-responses and strongly correlated movements of residues belonging to those distal regions. Our study provides new mechanistic insights into key sites and their interactions that could be targeted for developing AlkB-variants with desirable alkane conversion functions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Multiscale Machine-Learned Modeling Infrastructure RAS

MuMMI RAS is the application component of the Multiscale Machine-Learned Modeling Infrastructure (MuMMI). It simulates RAS protein interactions at three scales of resolution coupled with ML-based selection and in-situ feedback. MuMMI RAS is particularly configured for massive scale, running thousands of simultaneous jobs with several terabytes of data on Lassen and Summit.

Moon, JosephY.↗
DOE OSTI2020

Imaging transient solidification behavior

Solidification processing offers the first opportunity to control microstructure, properties, and performance in metallic alloy components. Until recently, microstructural evaluations were limited to post-solidification characterization by destructive methods. We review the development of time-resolved, in situimaging techniques capable of capturing solid-liquid interfacial evolution in metallic alloys with high spatial and temporal resolution underdiverse solidification conditions relevant for applications ranging from conventional directional solidification, crystal growth, and casting,to welding and additive manufacturing. These experiments enable direct visualization of transient behaviors that would otherwise remain unknown, uniquely providinginsights into the physics that impact microstructure and defect development and strategies for microstructuralcontrol and defect mitigation. Understanding microstructural evolution and the characteristics that form under various solidification conditions is essential for the development of multiscale, experimentallyinformed predictive modeling. This is highlightedby solidification simulations that utilize in situmeasurements of solidification dynamics from state-of-the-art experimental techniques.

36 MATERIALS SCIENCE↗
NASA NTRS2023

Highly Resolved Mapping of NO 2 Vertical Column Densities From GeoTASO Measurements Over a Megacity and Industrial Area During the KORUS-AQ Campaign

The Korea–United States Air Quality (KORUS-AQ) campaign is a joint study between the United States National Aeronautics and Space Administration (NASA) and the South Korea National Institute of Environmental Research (NIER) to monitor megacity and transboundary air pollution around the Korean Peninsula using airborne and ground-based measurements. Here, tropospheric nitrogen dioxide (NO 2 ) slant column density (SCD) measurements were retrieved from Geostationary Trace and Aerosol Sensor Optimization (GeoTASO) L1B data during the KORUS-AQ campaign (2 May to 10 June 2016). The retrieved SCDs were converted to tropospheric vertical column densities using the air mass factor (AMF) obtained from a radiative transfer calculation with trace gas profiles and aerosol property inputs simulated with the Community Multiscale Air Quality (CMAQ) model and surface reflectance data obtained from the Moderate Resolution Imaging Spectroradiometer (MODIS). For the first time, we examine highly resolved (250 m × 250 m resolution) tropospheric NO 2 over the Seoul and Busan metropolitan regions and the industrial region of Anmyeon. We reveal that the maximum NO 2 vertical column densities (VCDs) were 4.94×10 16 and 1.46×10 17 molec. cm −2 at 09:00 and 15:00 LT over Seoul, respectively, 6.86×10 16 and 4.89×10 16 molec. cm −2 in the morning and afternoon over Busan, respectively, and 1.64×10 16 molec. cm −2 over Anmyeon. The VCDs retrieved from the GeoTASO airborne instrument were correlated with those obtained from the Ozone Monitoring Instrument (OMI) (r=0.48), NASA's Pandora Spectrometer System (r=0.91), and NO 2 mixing ratios obtained from in situ measurements (r=0.07 in the morning, r=0.26 in the afternoon over the Seoul, and r>0.56 over Busan). Based on our results, GeoTASO is useful for identifying NO 2 hotspots and their spatial distribution in highly populated cities and industrial areas.

Gyo-Hwang Choo↗
DOE OSTI2023

Modeling the 4D discharge of lithium-ion batteries with a multiscale time-dependent deep learning framework

The lithium-ion battery (LIB) field is moving towards the direction of investigating spatially resolved physical phenomena in the 3D porous microstructure of electrodes. These pore-scale simulations give new insights into the local dynamics of lithiation/de-lithiation and charge transport, Nevertheless, the computational time of these simulations limits the integration of these models in optimization workflows of cycling conditions or electrode manufacturing processes. Machine learning models present a way of assessing in real-time the performance of materials. While several successful techniques for replicating simulations with machine learning have been proposed, this case study presents a more demanding problem, due to the necessity of understanding the behavior of heterogeneous 3D local data, as it evolves in time: this poses both a scientific and a technical challenge. To this end, we propose an autoregressive multiscale convolutional neural network model to predict relevant quantities at the pore-scale in the solid phase: the lithium concentration (in the active material) and potential (in the active material and carbon binder). Here, these are ultimately used to reconstruct the battery discharge curve. 3D images of the electrode microstructures are the input to the network, trained with a dataset of finite element method simulations to predict the discharge behavior of the cathode side in lithium ion batteries. We propose this machine learning model as a proof-of-concept of the applicability of multiscale networks for time-dependent physics problems. The trained model exhibits very high accuracy (with errors lower than 2 %) in forecasting the discharge behavior of new unseen cathodes.

25 ENERGY STORAGE↗