SEARCH · Search NASA

Results for “Multiscale simulations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

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

Oblique propagation and temperature effects on the resonant right-hand ion beam instability

The resonant right-hand instability (RHI) is often the dominant mode driven by reflected ions upstream of Earth’s quasi-parallel bow shock. In the tradition of Peter Gary, this paper further explores the right-hand instability using numerical solutions of the plasma dispersion relation and non-linear kinetic simulations, with parameters inspired by observations from NASA’s Magnetospheric Multiscale (MMS) mission. Agreement is found between the ion distributions in the particle-in-cell simulations and Magnetospheric Multiscale spacecraft data, which show the gyrophase bunching characteristic of the instability. The non-linear structures created by right-hand instability tend to be stronger when the plasma beta is lower. These structures have sizes of around 100 to 200 ion inertial lengths perpendicular to the magnetic field, presenting planet-sized disturbances to the magnetosphere. 2d and 3D hybrid particle-in-cell simulations show that modes with a range of propagation angles oblique to the magnetic field are excited, providing a ground to understand previous statistical studies of observed foreshock waves.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2023

Oblique Propagation and Temperature Effects on the Resonant Right-Hand Ion Beam Instability

Ari Le↗

DOE OSTI2017

Optimizing Discharge Capacity of Li-O 2 Batteries by Design of Air-Electrode Porous Structure: Multifidelity Modeling and Optimization

We develop a new mathematical framework to study the optimal design of air electrode microstructures for lithium-oxygen (Li-O2) batteries. It can eectively reduce the number of expensive experiments for testing dierent air-electrodes, thereby minimizing the cost in the design of Li-O2 batteries. The design parameters to characterize an air-electrode microstructure include the porosity, surface-to-volume ratio, and parameters associated with the pore-size distribution. A surrogate model (also known as response surface) for discharge capacity is rst constructed as a function of these design parameters. The surrogate model is accurate and easy to evaluate such that an optimization can be performed based on it. In particular, a Gaussian process regression method, co-kriging, is employed due to its accuracy and eciency in predicting high-dimensional responses from a combination of multidelity data. Specically, a small amount of data from high-delity simulations are combined with a large number of data obtained from computationally ecient low-delity simulations. The high-delity simulation is based on a multiscale modeling approach that couples the microscale (pore-scale) and macroscale (device-scale) models. Whereas, the low-delity simulation is based on an empirical macroscale model. The constructed response surface provides quantitative understanding and prediction about how air electrode microstructures aect the discharge performance of Li-O2 batteries. The succeeding sensitivity analysis via Sobol indices and optimization via genetic algorithm ultimately oer a reliable guidance on the optimal design of air electrode microstructures. The proposed mathematical framework can be generalized to investigate other new energy storage techniques and materials.

25 ENERGY STORAGE↗

DOE OSTI2023

Ab Initio Quantum Information Processor Design with Single-Molecule Magnets: A Multiscale Modeling Approach (Final Report)

This final report summarizes the team's efforts to develop a multiscale modeling approach that ranges from different levels of ab-initio quantum chemistry simulations to effective models and time-dependent external control, and to use this approach to systematically design quantum information processors with TbPc 2 single-molecule magnets. The impact of the work is two-fold: (i) New quantum chemistry simulation techniques capable of treating complex, multiscale problems such as the TbPc 2 molecule were developed, and (ii) the prospects for building quantum processors based on single-molecule magnets coupled by superconducting transmission line resonators were analyzed. The outcomes of this project revealed that current technology is at the cusp of being able to realize the main components of such a processor, and they highlighted the need to achieve stronger molecule-resonator interactions to enhance the viability of this approach. The multiscale modeling techniques developed during this project are general and transferable to other molecules and will thus have a broad impact on the field of quantum chemistry. Methods for controlling and simulating many coupled qubits developed here will also impact other quantum information technologies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2017

Mesoscopic-microscopic spatial stochastic simulation with automatic system partitioning

The reaction-diffusion master equation (RDME) is a model that allows for efficient on-lattice simulation of spatially resolved stochastic chemical kinetics. Compared to off-lattice hard-sphere simulations with Brownian dynamics or Green’s function reaction dynamics, the RDME can be orders of magnitude faster if the lattice spacing can be chosen coarse enough. However, strongly diffusion-controlled reactions mandate a very fine mesh resolution for acceptable accuracy. It is common that reactions in the same model differ in their degree of diffusion control and therefore require different degrees of mesh resolution. This renders mesoscopic simulation inefficient for systems with multiscale properties. Mesoscopic-microscopic hybrid methods address this problem by resolving the most challenging reactions with a microscale, off-lattice simulation. However, all methods to date require manual partitioning of a system, effectively limiting their usefulness as “black-box” simulation codes. In this paper, we propose a hybrid simulation algorithm with automatic system partitioning based on indirect a priori error estimates. We demonstrate the accuracy and efficiency of the method on models of diffusion-controlled networks in 3D.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Molecular dynamics-based multiscale damage initiation model for CNT/epoxy nanopolymers

This study presents a methodology that accurately simulates the brittle behavior of epoxy polymers initiating at the molecular level due to bond elongation and subsequent bond dissociation. The system investigated in this study comprises a combination of crystalline carbon nanotubes (CNTs) dispersed in epoxy polymer molecules. Molecular dynamics (MD) simulations are performed with an appropriate bond order-based force field to capture deformation-induced bond dissociation between atoms within the simulation volume. During deformation, the thermal vibration of molecules causes the elongated bonds to re-equilibrate; thus, the effect of mechanical deformation on bond elongation and scission cannot be captured effectively. This issue is overcome by deforming the simulation volume at zero temperature—a technique adopted from the concept of quasi-continuum and demonstrated successfully in the authors’ previous work. Results showed that a combination of MD deformation tests with ultra-high strain rates at near-zero temperatures provides a computationally efficient alternative for the study of bond dissociation phenomenon in amorphous epoxy polymer. In this paper, the ultra-high strain rate deformation approach is extended to the CNT-epoxy system at various CNT weight fractions and the corresponding bond disassociation energy extracted from the simulation volume is used as input to a low-fidelity continuum damage mechanics (CDM) model to demonstrate the bridging of length scales and to study matrix failure at the microscale. The material parameters for the classical CDM model are directly obtained from physics-based atomistic simulations, thus improving the accuracy of the multiscale approach.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Multiscale investigation of the mechanism of biomass deconstruction in the dimethyl isosorbide/water Co-solvent pretreatment system

In the context of promoting a circular bioeconomy, the development of green and efficient lignocellulosic biomass pretreatment technologies so as to realize high value-added biomass utilization is of intense interest. We demonstrated the potential of the bio-based green solvent dimethyl isosorbide (DMI) for the fractionation of Eucalyptus biomass with excellent performance. Here, to investigate the mechanisms involved in biomass fractionation, microimaging and microspectroscopic techniques were employed together with molecular dynamics (MD) simulation and COSMO-RS quantum chemical calculations to derive multiscale information. Both the microstructure and regional chemistry of the cell wall vary significantly with the volume ratio of DMI/H 2 O. The strongest effects were found at DMI/H 2 O = 9:1 and showed visible cell wall tearing cracks and cell wall deformation and collapse as well as the lowest values of cell wall thickness and circularity. From the MD simulations, lignin exhibits collapsed-like structure in pure H 2 O with low solvent accessibility surface area (SASA) and radius of gyration (R g ). In contrast, lignin in DMI/H 2 O shows extended structure with high SASA and solvent interactions dominated by van der Waals forces, with maximal contact in the 9:1 (v/v) system. Further, the COSMO-RS calculated sigma (σ-) potential suggests the intermolecular interactions in DMI and DMI/H 2 O co-solvent are weak, leading to stronger interaction with lignin and correspondingly higher lignin dissolution. The radial distribution functions and σ-potential all show that again DMI/H 2 O at 9:1 is an optimal volume ratio for high lignin dissolution. Finally, this study provides a solvent-ratio dependent mechanism for the action of polar aprotic solvents in the deconstruction of biomass.

09 BIOMASS FUELS↗

DOE OSTI2019

Qualitative Analysis, Structure-Preserving Discretizations and Scalable Solvers for Multiphysics Problems. Final Technical Report

Many problems in physics, engineering, meteorology, computer science, and other fields involve a certain level of coupling between different physical quantities. With the development of advanced computational modeling, simulation, and hardware, simulations for many important multiphysics/multiscale coupling effects are now more tractable. In this project, we have developed robust, efficient, and practical numerical algorithms for such simulation software that are capable of tackling these complex mathematical models in many fundamental ways such as stability, wellposedness, and consistency. We expect that such an integrated theoretical and algorithmic framework will have a substantial impact on numerous multiscale and multiphysics modeling and simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

HL-2A's ELM cycle simulations by integrating BOUT++'s drift MHD and transport code

A new integrating model has been developed to couple tokamak edge multiscale magnetohydrodynamic (MHD) events and transport simulations, such as edge-localized mode (ELM) cycles. As a proof of principle, we first start from a set of three-field two-fluid model equations, which includes the pressure, current, and vorticity. Here, the equations are separated into the slowly evolving part of the axisymmetric component by taking a time average of the axisymmetric component. The time-averaged fluxes, which are quadratic in fluctuating quantities, act as driven terms for the time-averaged axisymmetric quantities that determine the plasma transport, and therefore the large-scale evolution of the plasma profiles. Then the HL-2A's ELM cycles are simulated using the model. Good agreements of ELM size and pedestal recovery time have been achieved for the solutions obtained from the coupled simulation compared with experiment. For one ELM cycle simulation, the coupled code can achieve a speedup of a factor of up to 30 over standalone code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

A review of predictive nonlinear theories for multiscale modeling of heterogeneous materials

Since the beginning of the industrial age, material performance and design have been in the midst of innovation of many disruptive technologies. Today's electronics, space, medical, transportation, and other industries are enriched by development, design and deployment of composite, heterogeneous and multifunctional materials. As a result, materials innovation is now considerably outpaced by other aspects from component design to product cycle. In this article, we review predictive nonlinear theories for multiscale modeling of heterogeneous materials. Deeper attention is given to multiscale modeling in space and to computational homogenization in addressing challenging materials science questions. Moreover, we discuss a state-of-the-art platform in predictive image-based, multiscale modeling with co-designed simulations and experiments that executes on the world's largest supercomputers. Such a modeling framework consists of experimental tools, computational methods, and digital data strategies. Once fully completed, this collaborative and interdisciplinary framework can be the basis of Virtual Materials Testing standards and aids in the development of new material formulations. Moreover, it will decrease the time to market of innovative products.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Multiscale mechanical design of the lightweight, stiff, and damage-tolerant cuttlebone: A computational study

Cuttlebone, the endoskeleton of cuttlefish, offers an intriguing biological structural model for designing low-density cellular ceramics with high stiffness and damage tolerance. Cuttlebone is highly porous (porosity ~93%) and lightweight (density less than 20% of seawater), constructed mainly by brittle aragonite (95 wt%), but capable of sustaining hydrostatic water pressures over 20 atmospheres and exhibits energy absorption capability under compression comparable to many metallic foams (~4.4 kJ/kg). Here, in this work, we computationally investigate how such remarkable mechanical efficiency is enabled by the multiscale structure of cuttlebone. Using the common cuttlefish, Sepia Officinalis, as a model system, we first conducted high-resolution synchrotron micro-computed tomography (µ-CT) and quantified the cuttlebone's multiscale geometry, including the 3D asymmetric shape of individual walls, the wall assembly patterns, and the long-range structural gradient of walls across the entire cuttlebone (ca. 38 chambers). The acquired 3D structural information enables systematic finite-element simulations, which further reveal the multiscale mechanical design of cuttlebone: at the wall level, wall asymmetry provides optimized energy absorption while maintaining high structural stiffness; at the chamber level, variation of walls (number, pattern, and waviness amplitude) contributes to progressive damage; at the entire skeletal level, the gradient of chamber heights tailors the local mechanical anisotropy of the cuttlebone for reduced stress concentration. Our results provide integrated insights into understanding the cuttlebone's multiscale mechanical design and provide useful knowledge for the designs of lightweight cellular ceramics.

36 MATERIALS SCIENCE↗

NASA NTRS2012

The Geometric Factor of Electrostatic Plasma Analyzers: A Case Study from the Fast Plasma Investigation for the Magnetospheric Multiscale mission

We report our findings comparing the geometric factor (GF) as determined from simulations and laboratory measurements of the new Dual Electron Spectrometer (DES) being developed at NASA Goddard Space Flight Center as part of the Fast Plasma Investigation on NASA's Magnetospheric Multiscale mission. Particle simulations are increasingly playing an essential role in the design and calibration of electrostatic analyzers, facilitating the identification and mitigation of the many sources of systematic error present in laboratory calibration. While equations for laboratory measurement of the Geometric Factpr (GF) have been described in the literature, these are not directly applicable to simulation since the two are carried out under substantially different assumptions and conditions, making direct comparison very challenging. Starting from first principles, we derive generalized expressions for the determination of the GF in simulation and laboratory, and discuss how we have estimated errors in both cases. Finally, we apply these equations to the new DES instrument and show that the results agree within errors. Thus we show that the techniques presented here will produce consistent results between laboratory and simulation, and present the first description of the performance of the new DES instrument in the literature.

Magnetospheric Multiscale Mission↗

DOE OSTI2020

Multiscale modeling of microwave-heated multiphase systems

Modeling of microwave heating of multiphase systems requires resolving the electromagnetic field from the single particle to the entire cavity. We introduce a multiscale methodology for computationally affordable simulations of microwave heating of multiphase systems consisting of one phase dispersed in a continuum phase. The methodology homogenizes the original multiphase system, by taking advantage of the large separation of length scales, to calculate the effective permittivity, effective thermal conductivity, and volumetric power absorbed in each phase. The methodology is rigorously assessed against particle-resolved detailed numerical simulations of model systems. We demonstrate a significant reduction in computational cost while retaining the accuracy of detailed simulations. Lastly, the present methodology enables high throughput exploration of novel multiphase system designs utilizing microwave heating.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Theory of x-ray photon correlation spectroscopy for multiscale flows

Complex multiscale flows associated with instabilities and turbulence are commonly induced under high-energy density (HED) conditions, but accurate measurement of their transport properties has been challenging. X-ray photon correlation spectroscopy (XPCS) with coherent x-ray sources can, in principle, probe material dynamics to infer transport properties using time autocorrelation of density fluctuations. Here we develop a theoretical framework for utilizing XPCS to study material diffusivity in multiscale flows. We extend single-scale shear flow theories to broadband flows using a multiscale analysis that captures shear and diffusion dynamics. Our theory is validated with simulated XPCS for Brownian particles advected in multiscale flows. We demonstrate the versatility of the method over several orders of magnitude in timescale using sequential-pulse XPCS, single-pulse x-ray speckle visibility spectroscopy (XSVS), and double-pulse XSVS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2019

Multiscale Feature-feature Interactions Control Patterns of Hyporheic Exchange in a Simulated Headwater Mountain Stream

Recent predictions of hyporheic exchange at the basin-scale assume individual features control exchange independently of each other, which has been demonstrated in relatively uniform, low-gradient rivers. However, this assumption may not hold in steep catchments where both the type and size of individual features may vary over short distances, leading to irregular patterns of feature dominance on hyporheic exchange flows. Also, steep longitudinal gradients support substantial downvalley flows in the subsurface, which may create feedbacks between adjacent features. In this study, we test the extent to which features interact with one another and whether they can be aggregated to make reach-scale predictions in a headwater mountain stream. Using systematic manipulations of a 2-D stream centerline model and spectral analyses, we test for the presence of both feature-feature and multiscale interactions. Our results show that changing the height of individual step-pool features can alter hyporheic flow fields in neighboring, and sometimes distant, features. Spectral analyses revealed two scales of streambed topography—a local scale of single features and an intermediate scale that spanned multiple local-scale features. All features produced hyporheic exchange, but turnover of deeper hyporheic water only occurred at a few key locations where local- and intermediate-scale features amplified each other. Further, shallow bedrock increases the ratio of local- and intermediate-scale flowpaths to regional-scale flowpaths. Conceptual models portraying hyporheic exchange as a series of nested flowpaths should include the interactions among streambed topographic features in mountain streams. These results have implications for upscaling, field experiments, and stream restoration in steep catchments.

54 ENVIRONMENTAL SCIENCES↗

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↗

NASA NTRSDate not supplied

3D Phase-Field Simulations of Pattern Formation During Freeze Casting

We present the results of a combined experimental and phase-field modeling study of pattern formation during freeze casting. Those studies use unidirectional freezing of simple binary liquid mixtures of water and sugars (sucrose and trehalose) in a temperature gradient, which suffice to produce hierarchical templated structures similar to those observed in more complex multi-component freeze-cast systems including lamellae, undulated ridges, and more exotic “jellyfish-like” substructures. Multiscale 3D phase-field simulations reproduce remarkably well those structures quantitatively and identify key properties of the ice-water interface that control their formation. They further reveal that lamellae form as a result of a novel symmetry-breaking secondary instability of partially faceted cellular structures and pinpoint additional secondary instability mechanisms giving rise to smaller-scale substructures.

Kaihua Ji↗