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 451 records · Page 25

Probabilistic Methods for Structural Reliability and Risk

A probabilistic method is used to evaluate the structural reliability and risk of select metallic and composite structures. The method is a multiscale, multifunctional and it is based on the most elemental level. A multi-factor interaction model is used to describe the material properties which are subsequently evaluated probabilistically. The metallic structure is a two rotor aircraft engine, while the composite structures consist of laminated plies (multiscale) and the properties of each ply are the multifunctional representation. The structural component is modeled by finite element. The solution method for structural responses is obtained by an updated simulation scheme. The results show that the risk for the two rotor engine is about 0.0001 and the composite built-up structure is also 0.0001.

Chamis, Christos C.↗

Probabilistic Methods for Structural Reliability and Risk

A probabilistic method is used to evaluate the structural reliability and risk of select metallic and composite structures. The method is a multiscale, multifunctional and it is based on the most elemental level. A multifactor interaction model is used to describe the material properties which are subsequently evaluated probabilistically. The metallic structure is a two rotor aircraft engine, while the composite structures consist of laminated plies (multiscale) and the properties of each ply are the multifunctional representation. The structural component is modeled by finite element. The solution method for structural responses is obtained by an updated simulation scheme. The results show that the risk for the two rotor engine is about 0.0001 and the composite built-up structure is also 0.0001.

Chamis, Christos C.↗

Peridynamic Model for Single-Layer Graphene Obtained from Coarse Grained Bond Forces

An ordinary state-based peridynamic material model is proposed for single sheet graphene. The model is calibrated using coarse grained molecular dynamics simulations. The coarse graining method allows the dependence of bond force on bond length to be determined, including the horizon. The peridynamic model allows the horizon to be rescaled, providing a multiscale capability and allowing for substantial reductions in computational cost compared with molecular dynamics. The calibrated peridynamic model is compared to experimental data on the deflection and perforation of a graphene monolayer by an atomic force microscope probe.

36 MATERIALS SCIENCE↗

Modeling dust mineralogical composition: sensitivity to soil mineralogy atlases and their expected climate impacts

Soil dust aerosols are a key component of the climate system, as they interact with short- and long-wave radiation, alter cloud formation processes, affect atmospheric chemistry and play a role in biogeochemical cycles by providing nutrient inputs such as iron and phosphorus. The influence of dust on these processes depends on its physicochemical properties, which, far from being homogeneous, are shaped by its regionally varying mineral composition. The relative amount of minerals in dust depends on the source region and shows a large geographical variability. However, many state-of-the-art Earth system models (ESMs), upon which climate analyses and projections rely, still consider dust mineralogy to be invariant. The explicit representation of minerals in ESMs is more hindered by our limited knowledge of the global soil composition along with the resulting size-resolved airborne mineralogy than by computational constraints. In this work we introduce an explicit mineralogy representation within the state-of-the-art Multiscale Online Nonhydrostatic AtmospheRe CHemistry (MONARCH) model. We review and compare two existing soil mineralogy datasets, which remain a source of uncertainty for dust mineralogy modeling and provide an evaluation of multiannual simulations against available mineralogy observations. Soil mineralogy datasets are based on measurements performed after wet sieving, which breaks the aggregates found in the parent soil. Our model predicts the emitted particle size distribution (PSD) in terms of its constituent minerals based on brittle fragmentation theory (BFT), which reconstructs the emitted mineral aggregates destroyed by wet sieving. Our simulations broadly reproduce the most abundant mineral fractions independently of the soil composition data used. Feldspars and calcite are highly sensitive to the soil mineralogy map, mainly due to the different assumptions made in each soil dataset to extrapolate a handful of soil measurements to arid and semi-arid regions worldwide. For the least abundant or more difficult-to-determine minerals, such as iron oxides, uncertainties in soil mineralogy yield differences in annual mean aerosol mass fractions of up to ~ 100 %. Although BFT restores coarse aggregates including phyllosilicates that usually break during soil analysis, we still identify an overestimation of coarse quartz mass fractions (above 2 µm in diameter). In a dedicated experiment, we estimate the fraction of dust with undetermined composition as given by a soil map, which makes up ~ 10 % of the emitted dust mass at the global scale and can be regionally larger. Changes in the underlying soil mineralogy impact our estimates of climate-relevant variables, particularly affecting the regional variability of the single-scattering albedo at solar wavelengths or the total iron deposited over oceans. All in all, this assessment represents a baseline for future model experiments including new mineralogical maps constrained by high-quality spaceborne hyperspectral measurements, such as those arising from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) mission.

54 ENVIRONMENTAL SCIENCES↗

Modeling Dust Mineralogical Composition: Sensitivity to Soil Mineralogy Atlases and Their Expected Climate Impacts

Soil dust aerosols are a key component of the climate system, as they interact with short- and long-wave radiation, alter cloud formation processes, affect atmospheric chemistry and play a role in biogeochemical cycles by providing nutrient inputs such as iron and phosphorus. The influence of dust on these processes depends on its physicochemical properties, which, far from being homogeneous, are shaped by its regionally varying mineral composition. The relative amount of minerals in dust depends on the source region and shows a large geographical variability. However, many state-of-the-art Earth system models (ESMs), upon which climate analyses and projections rely, still consider dust mineralogy to be invariant. The explicit representation of minerals in ESMs is more hindered by our limited knowledge of the global soil composition along with the resulting size-resolved airborne mineralogy than by computational constraints. In this work we introduce an explicit mineralogy representation within the state-of-the-art Multiscale Online Nonhydrostatic AtmospheRe CHemistry (MONARCH) model. We review and compare two existing soil mineralogy datasets, which remain a source of uncertainty for dust mineralogy modeling and provide an evaluation of multiannual simulations against available mineralogy observations. Soil mineralogy datasets are based on measurements performed after wet sieving, which breaks the aggregates found in the parent soil. Our model predicts the emitted particle size distribution (PSD) in terms of its constituent minerals based on brittle fragmentation theory (BFT), which reconstructs the emitted mineral aggregates destroyed by wet sieving. Our simulations broadly reproduce the most abundant mineral fractions independently of the soil composition data used. Feldspars and calcite are highly sensitive to the soil mineralogy map, mainly due to the different assumptions made in each soil dataset to extrapolate a handful of soil measurements to arid and semi-arid regions worldwide. For the least abundant or more difficult-to-determine minerals, such as iron oxides, uncertainties in soil mineralogy yield differences in annual mean aerosol mass fractions of up to ∼ 100 %. Although BFT restores coarse aggregates including phyllosilicates that usually break during soil analysis, we still identify an overestimation of coarse quartz mass fractions (above 2 µm in diameter). In a dedicated experiment, we estimate the fraction of dust with undetermined composition as given by a soil map, which makes up ∼ 10 % of the emitted dust mass at the global scale and can be regionally larger. Changes in the underlying soil mineralogy impact our estimates of climate-relevant variables, particularly affecting the regional variability of the single-scattering albedo at solar wavelengths or the total iron deposited over oceans. All in all, this assessment represents a baseline for future model experiments including new mineralogical maps constrained by high-quality spaceborne hyperspectral measurements, such as those arising from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) mission.

Soil dust↗

Recent advances in machine learning towards multiscale soft materials design

The multiscale design of soft materials requires an ensemble of computational techniques spanning quantum-chemistry to molecular dynamics to continuum modeling. The recent emergence of machine-learning (ML) and modern optimization algorithms has accelerated material property prediction, as well as stimulated the development of hybrid ML/molecular modeling methodologies capable of providing physical insights unobtainable from purely physics-based modeling and intuition. Such hybrid techniques also have important ramifications for the ML-enhanced interpretation of results from simulations and experiments alike. Leveraging ML techniques for the design of chemical or morphological structures based on a target property or functionality represents an exciting goal for the general area of soft materials, including polymers, liquid crystals, colloids, or biomolecules, to name a few representative classes of systems. Furthermore, we provide a perspective on recent work using ML techniques of relevance for the multiscale design of soft materials and outline potential future directions of interest to the soft materials community.

36 MATERIALS SCIENCE↗

Fluid modeling of low-temperature plasmas

Fluid models are essential for understanding and predicting low-temperature plasma (LTP) behavior in various scientific and industrial settings. This paper provides an introductory tutorial on fluid modeling of LTPs, covering model formulation, implementation, and computational simulations. The tutorial focuses on five main components of the formulation of LTP fluid models: fluid flow, energy, chemistry, electromagnetism, and material properties, as well as in essential aspects of model implementations, including multiscale phenomena, multiphysics coupling, and numerical convergence. Designed for students and early-career researchers, this work offers a practical foundation for developing and using fluid models, from in-house computational codes to commercial software, bridging fundamental theory with real-world applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Exploring far-from-equilibrium ultrafast polarization control in ferroelectric oxides with excited-state neural network quantum molecular dynamics

Ferroelectric materials exhibit a rich range of complex polar topologies, but their study under far-from-equilibrium optical excitation has been largely unexplored because of the difficulty in modeling the multiple spatiotemporal scales involved quantum-mechanically. To study optical excitation at spatiotemporal scales where these topologies emerge, we have performed multiscale excited-state neural network quantum molecular dynamics simulations that integrate quantum-mechanical description of electronic excitation and billion-atom machine learning molecular dynamics to describe ultrafast polarization control in an archetypal ferroelectric oxide, lead titanate. Far-from-equilibrium quantum simulations reveal a marked photo-induced change in the electronic energy landscape and resulting cross-over from ferroelectric to octahedral tilting topological dynamics within picoseconds. The coupling and frustration of these dynamics, in turn, create topological defects in the form of polar strings. The demonstrated nexus of multiscale quantum simulation and machine learning will boost not only the emerging field of ferroelectric topotronics but also broader optoelectronic applications.

36 MATERIALS SCIENCE↗

Description and Evaluation of the Multiscale Online Nonhydrostatic AtmospheRe CHemistry Model (NMMB-MONARCH) Version 1.0: Gas-Phase Chemistry at Global Scale

This paper presents a comprehensive description and benchmark evaluation of the tropospheric gas-phase chemistry component of the Multiscale Online Nonhydrostatic AtmospheRe CHemistry model (NMMBMONARCH), formerly known as NMMB/BSC-CTM, that can be run on both regional and global domains. Here, we provide an extensive evaluation of a global annual cycle simulation using a variety of background surface stations (EMEP, WDCGG and CASTNET), ozonesondes (WOUDC, CMD and SHADOZ), aircraft data (MOZAIC and several campaigns), and satellite observations (SCIAMACHY and MOPITT).We also include an extensive discussion of our results in comparison to other state-of-the-art models. We note that in this study, we omitted aerosol processes and some natural emissions (lightning and volcano emissions). The model shows a realistic oxidative capacity across the globe. The seasonal cycle for CO is fairly well represented at different locations (correlations around 0.3-0.7 in surface concentrations), although concentrations are underestimated in spring and winter in the Northern Hemisphere, and are overestimated throughout the year at 800 and 500 hPa in the Southern Hemisphere. Nitrogen species are well represented in almost all locations, particularly NO2 in Europe (root mean square error - RMSE - below 5 ppb). The modeled vertical distributions of NOx and HNO3 are in excellent agreement with the observed values and the spatial and seasonal trends of tropospheric NO2 columns correspond well to observations from SCIAMACHY, capturing the highly polluted areas and the biomass burning cycle throughout the year. Over Asia, the model underestimates NOx from March to August, probably due to an underestimation of NOx emissions in the region. Overall, the comparison of the modeled CO and NO2 with MOPITT and SCIAMACHY observations emphasizes the need for more accurate emission rates from anthropogenic and biomass burning sources (i.e., specification of temporal variability).

atmospheric chemistry↗

Particle Acceleration by Dense Impulsive Structures Moving in Ambient Magnetospheric Plasma. 3-D Hybrid Kinetic Modeling and MMS Observations

High resolution observations of dense plasma impulsive structures moving through an ambient background magnetospheric flows were captured by the Magnetospheric Multiscale mission. The observations show particle heating and acceleration, shock-like wave formation, and whistler wave excitation inside the interface between the dense impulsive plasma structures and the ambient plasma. A multiscale hybrid kinetic simulation provides an explanation of the observed wave-particle interactions with the assumption that the dense plasma structures may be represented by plasma clouds which are formed at the magnetopause layer due to reconnection processes. Plain Language Summary: Dense, impulsive plasma structures moving through a background plasma were captured by the NASA Magnetospheric Multiscale mission. The observations show that the dense structures can generate strong perturbations in the electromagnetic field and shock-like waves. Interactions between these electromagnetic waves and the particles results in particle acceleration. 3-D hybrid kinetic modeling (particle description for ions and fluid description for electrons) was used to investigate the plasma physics of the observed structures. It was assumed that the plasma clouds were formed by magnetic field reconnection inside the magnetopause, which is the interface between the solar wind particles and the cold low-density magnetospheric plasma. The work helps us understand the plasma environment at the interface between the Earth and solar wind, near planetary moons, within astrophysical explosions, and possibly at the interface between the solar wind and local interstellar medium.

A. S. Lipatov↗

Improvement of the $\mathrm{BISON U_3Si_2}$ modeling capabilities based on multiscale developments to modeling fission gas behavior

Uranium silicide (U 3 Si 2 ) is a concept explored as a potential alternative to UO 2 fuel used in light water reactors (LWRs) since it may improve accident tolerance and economics due to its higher thermal conductivity and increased uranium density. U 3 Si 2 has been previously used in research reactors in the form of dispersion fuel, but operated at lower temperatures than commercial LWRs. The research reactor data illustrated that significant gaseous swelling occurs as the fuel burnup increases. Therefore, it is imperative to understand the fission gas behavior of U 3 Si 2 under higher temperature LWR operating conditions. In this work, molecular dynamics and phase-field modeling techniques are used to reduce the uncertainty in select modeling assumptions made in developing the fission gas behavior model for U 3 Si 2 in the BISON fuel performance code. These lower length scale informed models are then utilized in the validation of BISON U 3 Si 2 modeling capabilities to simulate the ATF-1 experiments irradiated in the Advanced Test Reactor (ATR). Sensitivity analysis (SA) and uncertainty quantification (UQ) are included as part of the validation process to identify where further experiments and lower length scale modeling would be beneficial. Here, the multiscale modeling approach utilized in this work can be applied to new fuel concepts being explored for both LWRs and advanced reactors (e.g., uranium nitride, uranium carbide).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Thermoviscoelastic modelling of high strain thin-ply composites by means of multiscale plate and beam model

High strain thin-ply (HS-TPC) technology is being increasingly adopted for high-performanceaerospace applications. Albeit many of these structures such as deployable booms can bemodeled as one-dimensional beam problems, there is a lack of themoviscoelastic beam mod-els to efficiently and accurately simulate HS-TPC. This work will use mechanics of struc-ture genome (MSG) to construct linear thermoviscoelastic beam models that can homogenizethree-dimensional heterogeneous materials made of constituents with time- and temperature-dependent behavior. The formulation derives the transient strain energy based on integralformulation for thermorheologically simple materials subject to finite temperature changeswith the restriction that the strain is small. A lenticular boom is used as a numerical exampleto verify the MSG-based linear thermoviscoelastic beam model against MSG-based shell/platemodel, which has already been validated against experimental data provided by NASA, anddirect numerical simulations performed in a finite element commercial package.

Orzuri Rique Garaizar↗

Coupled SAM/Griffin Model of a Reference Fluoride-Salt-Cooled High-Temperature Reactor for Multi-Physics Simulations

A multi-physics coupled simulation model of a reference pebble bed fluoride-salt-cooled high-temperature reactor (PB-FHR) has been developed with SAM and Griffin computer codes for transient safety analysis. The reference problem of a prototypical reactor design serves as the foundation for the U.S. NRC (Nuclear Regulatory Commission) to verify the adequacy of computer codes and evaluation models for specific reactor types. In this work, the previously developed SAM model for PB-FHR primary loop has been updated for the coupled simulation. The updated SAM PB-FHR model includes a 2-D axial symmetric core region and external core components in 0-D/1-D. In addition to the primary loop, a detailed model of the RCCS (reactor cavity cooling system) is added. The 2-D and 1-D domains are tightly coupled using the single-solve approach developed in SAM. In the pebble bed region, the SAM multiscale explicit pebble model is applied to calculate the pebble and TRISO fuel kernel temperatures. The Griffin model used in this work is based on a model developed at Idaho National Laboratory in collaboration with the U.S. NRC. The Griffin neutronics model and SAM thermal hydraulics model is coupled with the Comprehensive Reactor Analysis Bundle (CRAB or alternately BlueCRAB) application. Both steady-state and transient scenarios are simulated to demonstrate the model's suitability for multi-physics simulations of PB-FHR transients.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The impact of capillary heterogeneity on CO 2 flow and trapping across scales

Capillary heterogeneity has been identified over the last decade as a key control on subsurface CO 2 flow behavior during geological CO 2 sequestration. These heterogeneities can be formed in all sedimentary rocks, ranging from slight variations in the sand grain sizes to extensive sequences of interbedded sands, shales, and limestones. Capillary heterogeneity has been largely, although not entirely, overlooked in subsurface flow modeling because it is assumed to only directly influence fluid redistribution over scales of centimeters to meters. However, even small-scale fluid movements can result in dramatic impacts on the mobility and trapping of the CO 2 over kilometers. Therefore, neglecting capillary heterogeneity at multiple scales could potentially lead to errors in modeling and predicting field-scale plume migration. In this review paper, we aim to provide a consistent overview to (1) establish that capillary heterogeneity can have a major impact on CO 2 plume migration, (2) establish the respective length scales at which capillary heterogeneity matters, and (3) provide guidance for numerical modeling. This review covers pertinent literature and extracts key observations from the core to the field scales. Experimental studies have shown that millimeter-decimeter scale capillary heterogeneity can cause the so-called capillary heterogeneity trapping in addition to pore-scale residual trapping. Even at such a small scale, capillary heterogeneity can already lead to complex upscaled constitutive relationships, such as flow-rate dependent and anisotropic relative permeability, which affects field-scale CO 2 migration even when field-scale heterogeneities are present. Under gravity-dominated flow regimes, centimeter-meter scale capillary heterogeneity can entrap a significant amount of CO 2 at field scale, not just after imbibition but also during drainage. In certain cases, the presence of capillary heterogeneity can even completely stop the vertical movement of the CO 2 plume, hence greatly reducing leakage risks. At meter-kilometer scale, the influence of capillary heterogeneity is more pronounced and can hinder or redirect CO 2 migration in both lateral and vertical directions. The impact of capillary heterogeneity across multiple spatial scales poses a great challenge in modeling CO 2 migration at field scale, because it is practically impossible to build a field-scale earth model with grid blocks at millimeter scale. We recommend a hierarchical modeling approach to address this challenge. At field scale, earth models are built to capture geological features and heterogeneities in high but still practical grid resolutions. For each facies or rock type of the field-scale model, high- resolution meter-scale “conceptual” models are built with millimeter-scale grid blocks to capture representative fine-scale bedding geometries and heterogeneities in various environments of deposition, bridging the gap from subcore scale to the size of a field-scale simulation grid block. Upscaling is then used to preserve the smaller-scale flow dynamics of various rock types in field-scale simulations. Here, future work is needed to (1) refine, improve, and validate the hierarchical modeling approach; (2) build libraries of fine-scale bedding models for facies in various environments of deposition; (3) quantify multiscale capillary heterogeneity effects under subsurface uncertainties; (4) gain learning from different storage formations; and (5) establish best practices that balance accuracy and computational speed.

Capillary heterogeneity↗

BattPhase—A Convergent, Non-Oscillatory, Efficient Algorithm and Code for Predicting Shape Changes in Lithium Metal Batteries Using Phase-Field Models: Part I. Secondary Current Distribution

Electrochemical models at different scales and varying levels of complexity have been used in the literature to study the evolution of the anode surface in lithium metal batteries. This includes continuum, mesoscale (phase-field approaches), and multiscale models. Thermodynamics-based equations have been used to study phase changes in lithium batteries using phase-field approaches. However, grid convergence studies and the effect of additional parameters needed to simulate these models are not well-documented in the literature. In this paper, using a motivating example of a moving boundary model in one- and two-dimensions, we show how one can formulate phase-field models, implement algorithms for the same and analyze the results. An open-access code with no restrictions is provided as well. This article concludes with some thoughts on the computational efficiency of phase-field models for simulating dendritic growth.

25 ENERGY STORAGE↗

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↗

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↗

Comparing Convective Self‐Aggregation in Idealized Models to Observed Moist Static Energy Variability Near the Equator

Abstract Idealized convection‐permitting simulations of radiative‐convective equilibrium have become a popular tool for understanding the physical processes leading to horizontal variability of tropical water vapor and rainfall. However, the applicability of idealized simulations to nature is still unclear given that important processes are typically neglected, such as lateral water vapor advection by extratropical intrusions, or interactive ocean coupling. Here, we exploit spectral analysis to compactly summarize the multiscale processes supporting convective aggregation. By applying this framework to high‐resolution reanalysis data and satellite observations in addition to idealized simulations, we compare convective‐aggregation processes across horizontal scales and data sets. The results affirm the validity of the radiative‐convective equilibrium simulations as an analogy to the real world. Column moist static energy tendencies share similar signs and scale selectivity in convection‐permitting models and observations: Radiation increases variance at wavelengths above 1,000 km, while advection damps variance across wavelengths, and surface fluxes mostly reduce variance between 1,000 and 10,000 km.

Beucler, Tom↗