Search NASA⌕ Search

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.

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

At least 1,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Electrostatic shock dynamics in superthermal plasmas

The propagation of ion acoustic shocks in nonthermal plasmas is investigated, both analytically and numerically. An unmagnetized collisionless electron-ion plasma is considered, featuring a superthermal (non-Maxwellian) electron distribution, which is modeled by a {kappa}-(kappa) distribution function. Adopting a multiscale approach, it is shown that the dynamics of low-amplitude shocks is modeled by a hybrid Korteweg-de Vries-Burgers (KdVB) equation, in which the nonlinear and dispersion coefficients are functions of the {kappa} parameter, while the dissipative coefficient is a linear function of the ion viscosity. All relevant shock parameters are shown to depend on {kappa}: higher deviations from a pure Maxwellian behavior induce shocks which are narrower, faster, and of larger amplitude. The stability profile of the kink-shaped solutions of the KdVB equation against external perturbations is investigated. The spatial profile of the shocks is found to depend upon the dispersion and the dissipation term, and the role of the interplay between dispersion and dissipation is elucidated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

A multiscale restriction-smoothed basis method for high contrast porous media represented on unstructured grids

A wide variety of multiscale methods have been proposed in the literature to reduce runtime and provide better scaling for the solution of Poisson-type equations modeling flow in porous media. We present a new multiscale restricted-smoothed basis (MsRSB) method that is designed to be applicable to both rectilinear grids and unstructured grids. Like many other multiscale methods, MsRSB relies on a coarse partition of the underlying fine grid and a set of local prolongation operators (multiscale basis functions) that map unknowns associated with the fine grid cells to unknowns associated with blocks in the coarse partition. These mappings are constructed by restricted smoothing: Starting from a constant, a localized iterative scheme is applied directly to the fine-scale discretization to compute prolongation operators that are consistent with the local properties of the differential operators. The resulting method has three main advantages: First of all, both the coarse and the fine grid can have general polyhedral geometry and unstructured topology. This means that partitions and good prolongation operators can easily be constructed for complex models involving high media contrasts and unstructured cell connections introduced by faults, pinch-outs, erosion, local grid refinement, etc. In particular, the coarse partition can be adapted to geological or flow-field properties represented on cells or faces to improve accuracy. Secondly, the method is accurate and robust when compared to existing multiscale methods and does not need expensive recomputation of local basis functions to account for transient behavior: Dynamic mobility changes are incorporated by continuing to iterate a few extra steps on existing basis functions. This way, the cost of updating the prolongation operators becomes proportional to the amount of change in fluid mobility and one reduces the need for expensive, tolerance-based updates. Finally, since the MsRSB method is formulated on top of a cell-centered, conservative, finite-volume method, it is applicable to any flow model in which one can isolate a pressure equation. Herein, we only discuss single and two-phase incompressible models. Compressible flow, e.g., as modeled by the black-oil equations, is discussed in a separate paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Influence of Surface and Bulk Water Ice on the Reactivity of a Water-forming Reaction

On the surface of icy dust grains in the dense regions of the interstellar medium, a rich chemistry can take place. Due to the low temperature, reactions that proceed via a barrier can only take place through tunneling. The reaction H+H{sub 2}O{sub 2}⟶H{sub 2}O+OH is such a case with a gas-phase barrier of ∼26.5 kJ mol{sup −1}. Still, the reaction is known to be involved in water formation on interstellar grains. Here, we investigate the influence of a water ice surface and of bulk ice on the reaction rate constant. Rate constants are calculated using instanton theory down to 74 K. The ice is taken into account via multiscale modeling, describing the reactants and the direct surrounding at the quantum mechanical level with density functional theory (DFT), while the rest of the ice is modeled on the molecular mechanical level with a force field. We find that H{sub 2}O{sub 2} binding energies cannot be captured by a single value, but rather they depend on the number of hydrogen bonds with surface molecules. In highly amorphous surroundings, the binding site can block the routes of attack and impede the reaction. Furthermore, the activation energies do not correlate with the binding energies of the same sites. The unimolecular rate constants related to the Langmuir–Hinshelwood mechanism increase as the activation energy decreases. Thus, we provide a lower limit for the rate constant and argue that rate constants can have values up to two order of magnitude larger than this limit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Weakly Nonlinear Magnetorotational Instability in a Global, Cylindrical Taylor–Couette Flow

We conduct a global, weakly nonlinear analysis of the magnetorotational instability (MRI) in a Taylor–Couette flow. This is a multiscale, perturbative treatment of the nonideal, axisymmetric MRI near threshold, subject to realistic radial boundary conditions and cylindrical geometry. We analyze both the standard MRI, initialized by a constant vertical background magnetic field, and the helical MRI, with an azimuthal background field component. This is the first weakly nonlinear analysis of the MRI in a global Taylor–Couette geometry, as well as the first weakly nonlinear analysis of the helical MRI. We find that the evolution of the amplitude of the standard MRI is described by a real Ginzburg–Landau equation (GLE), whereas the amplitude of the helical MRI takes the form of a complex GLE. This suggests that the saturated state of the helical MRI may itself be unstable on long spatial and temporal scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

WE-AB-202-11: Radiobiological Modeling of Tumor Response During Radiotherapy Based On Pre-Treatment Dynamic PET Imaging Data

Purpose: To evaluate the ability of a multiscale radiobiological model of tumor response to predict mid-treatment hypoxia images, based on pretreatment imaging of perfusion and hypoxia with [18-F]FMISO dynamic PET and glucose metabolism with [18-F]FDG PET. Methods: A mechanistic tumor control probability (TCP) radiobiological model describing the interplay between tumor cell proliferation and hypoxia (Jeong et al., PMB 2013) was extended to account for intra-tumor nutrient heterogeneity, dynamic cell migration due to nutrient gradients, and stromal cells. This extended model was tested on 10 head and neck cancer patients treated with chemoradiotherapy, randomly drawn from a larger MSKCC protocol involving baseline and mid-therapy dynamic PET scans. For each voxel, initial fractions of proliferative and hypoxic tumor cells were obtained by finding an approximate solution to a system of linear equations relating cell fractions to voxel-level FDG uptake, perfusion (FMISO K{sub 1}) and hypoxia (FMISO k{sub 3}). The TCP model then predicted their evolution over time up until the mid treatment scan. Finally, the linear model was reapplied to predict each lesion’s median hypoxia level (k{sub 3}[med,sim]) which in turn was compared to the FMISO k{sub 3}[med] measured at mid-therapy. Results: The average k3[med] of the tumors in pre-treatment scans was 0.0035 min{sup −1}, with an inter-tumor standard deviation of σ[pre]=0.0034 min{sup −1}. The initial simulated k{sub 3}[med,sim] of each tumor agreed with the corresponding measurements within 0.1σ[pre]. In 7 out of 10 lesions, the mid-treatment k{sub 3}[med,sim] prediction agreed with the data within 0.3σ[pre]. The remaining cases corresponded to the most extreme relative changes in k{sub 3}[med]. Conclusion: This work presents a method to personalize the prediction of a TCP model using pre-treatment kinetic imaging data, and validates the modeling of radiotherapy response by predicting changes in median hypoxia values during treatment. Variations from predicted response may be a useful biomarker, which should be further explored. Partially supported by NIH grant #1 R01 CA157770-01A1 and a grant from Varian Corporation.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2023

Comparison of Dust Optical Depth From Multi-Sensor Products and MONARCH (Multiscale Online Non-hydrostatic AtmospheRe CHemistry) Dust Reanalysis Over North Africa, the Middle East, and Europe

Aerosol reanalysis datasets are model-based, observationally constrained, continuous 3D aerosol fields with a relatively high temporal frequency that can be used to assess aerosol variations and trends, climate effects, and impacts on socioeconomic sectors, such as health. Here we compare and assess the recently published MONARCH (Multiscale Online Non-hydrostatic AtmospheRe CHemistry) high-resolution regional desert dust reanalysis over northern Africa, the Middle East, and Europe (NAMEE) with a combination of ground-based observations and space-based dust retrievals and products. In particular, we compare the total and coarse dust optical depth (DOD) from the new reanalysis with DOD products derived from MODIS (MODerate resolution Imaging Spectroradiometer), MISR (Multi-angle Imaging SpectroRadiometer), and IASI (Infrared Atmospheric Sounding Interferometer) spaceborne instruments. Despite the larger uncertainties, satellite-based datasets provide a better geographical coverage than ground-based observations, and the use of different retrievals and products allows at least partially overcoming some single-product weaknesses in the comparison. Nevertheless, limitations and uncertainties due to the type of sensor, its operating principle, its sensitivity, its temporal and spatial resolution, and the methodology for retrieving or further deriving dust products are factors that bias the reanalysis assessment. We, therefore, also use ground-based DOD observations provided by 238 stations of the AERONET (AErosol RObotic NETwork) located within the NAMEE region as a reference evaluation dataset. In particular, prior to the reanalysis assessment, the satellite datasets were evaluated against AERONET, showing moderate underestimations in the vicinities of dust sources and downwind regions, whereas small or significant overestimations, depending on the dataset, can be found in the remote regions. Taking these results into consideration, the MONARCH reanalysis assessment shows that total and coarse-DOD simulations are consistent with satellite- and ground-based data, qualitatively capturing the major dust sources in the area in addition to the dust transport patterns. Moreover, the MONARCH reanalysis reproduces the seasonal dust cycle, identifying the increased dust activity that occurred in the NAMEE region during spring and summer. The quantitative comparison between the MONARCH reanalysis DOD and satellite multi-sensor products shows that the reanalysis tends to slightly overestimate the desert dust that is emitted from the source regions and underestimate the transported dust over the outflow regions, implying that the model's removal of dust particles from the atmosphere, through deposition processes, is too effective. More specifically, small positive biases are found over the Sahara desert (0.04) and negative biases over the Atlantic Ocean and the Arabian Sea (−0.04), which constitute the main pathways of the long-range dust transport. Considering the DOD values recorded on average there, such discrepancies can be considered low, as the low relative bias in the Sahara desert (< 50 %) and over the adjacent maritime regions (< 100 %) certifies. Similarly, over areas with intense dust activity, the linear correlation coefficient between the MONARCH reanalysis simulations and the ensemble of the satellite products is significantly high for both total and coarse DOD, reaching 0.8 over the Middle East, the Atlantic Ocean, and the Arabian Sea and exceeding it over the African continent. Moreover, the low relative biases and high correlations are associated with regions for which large numbers of observations are available, thus allowing for robust reanalysis assessment.

Michail Mytilinaios↗
DOE OSTI2020

On Simulating the Formation of Structured, Crystalline Systems via Non-classical Pathways

Observations in crystal growth and assembly from recent in situ methods suggest alternative, non-classical crystallization pathways play an important role in the determination of the micro- and meso- structures in crystalline systems. These processes display parallels that cross-cut multiple disciplines investigating crystallization across four orders of magnitude in size scales and widely differing environments, hinting that alternative crystal growth pathways may be a fundamental scheme in natural crystal formation. Using a system of short-range attractive microbeads, we demonstrate that the addition of a small concentration of sub-species incommensurate with the lattice spacing of the dominant species results in a stark change in crystal size and morphology. These changes are attributed to the presence of fleeting, amorphous-like configurations of beads that ultimately change the melting and growth dynamics in preferred directions. From these real-time observations, we hypothesize the amorphous mineral precursors present in biological mineralized tissues undergo similar non-classical crystallization processes resulting in the complex structures found in biomineralization.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A Langevin approach to multi-scale modeling

In plasmas, distribution functions often demonstrate long anisotropic tails or otherwise significant deviations from local Maxwellians. The tails, especially if they are pulled out from the bulk, pose a serious challenge for numerical simulations as resolving both the bulk and the tail on the same mesh is often challenging. A multi-scale approach, providing evolution equations for the bulk and the tail individually, could offer a resolution in the sense that both populations could be treated on separate meshes or different reduction techniques applied to the bulk and the tail population. In this paper, we propose a multi-scale method which allows us to split a distribution function into a bulk and a tail so that both populations remain genuine, non-negative distribution functions and may carry density, momentum, and energy. The proposed method is based on the observation that the motion of an individual test particle in a plasma obeys a stochastic differential equation, also referred to as a Langevin equation. Finally, this allows us to define transition probabilities between the bulk and the tail and to provide evolution equations for both populations separately.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2019

Revealing the Multiscale Nature of Turbulence in Space Plasmas with an Innovative Swarm of Spacecraft

Abstract: Turbulence is fundamentally a multiscale phenomena, with energy, mass, and momentum being transported across both spatial and temporal scales. In a magnetized collisionless plasma, there are unresolved questions about the structure of the turbulent transfer of energy, as well as how the energy is extracted from the cascade and dissipated as heat in the constituent charged particles. The plasma in the solar wind acts as an accessible natural laboratory to study these processes, and much progress has indeed been made since the dawn of the space age in understanding turbulence via in situ observations of turbulent plasmas near Earth. However, to date, these observations have been limited to a single, or at best, a tight cluster of points, leading to ambiguities in at what scales energy is contained, how it is transported, and by what mechanism it is dissipated. In this presentation, we describe a heliophysics mission concept aimed at understanding turbulence that is enabled by a swarm of small satellites. The proposed "HelioSwarm" mission will measure turbulent fields and flows and charged particles simultaneously at many points spanning size and time scales from the fluid to sub-ion regime. In doing so, we will be able to disentangle how the turbulence depends on time and space, directly observe the change in internal energy in the plasma, and definitively capture the dynamic relation between turbulence and structures. While the processes under examination are universal, arising throughout our solar system and universe, they are difficult to reproduce in either terrestrial laboratories or numerical simulations, meaning that the only a multipoint observatory in the near-Earth heliosphere will be able to study them in sufficient detail to discern the underlying physics. In this talk, we highlight the enabling role that small satellites play in providing closure on these long-standing but critically important science questions.

Spence, H. E.↗
DOE OSTI2021

Arrested Phase Separation in Double-Exchange Models: Large-Scale Simulation Enabled by Machine Learning

In this work, we present large-scale dynamical simulations of electronic phase separation in the single-band double-exchange model based on deep-learning neural-network potentials trained from small-size exact diagonalization solutions. We uncover an intriguing correlation-induced freezing behavior as doped holes are segregated from half filled insulating background during equilibration. While the aggregation of holes is stabilized by the formation of ferromagnetic clusters through Hund’s coupling between charge carriers and local magnetic moments, this stabilization also creates confining potentials for holes when antiferromagnetic spin-spin correlation is well developed in the background. The dramatically reduced mobility of the self-trapped holes prematurely disrupts further growth of the ferromagnetic clusters, leading to an arrested phase separation. Implications of our findings for phase separation dynamics in materials that exhibit colossal magnetoresistance effect are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Predicted Fracture Tendency of Naturally Occurring Aluminum Surface Coatings under Tensile Loading

Naturally occurring coatings on aluminum metal, such as its oxide or hydroxide, serve to protect the material from corrosion. Understanding the conditions under which these coatings mechanically fail is therefore expected to be an important aspect of predictive models for aluminum component lifetimes. To this end, we develop and apply a molecular dynamics (MD) modeling framework for conducting tension tests that is capable of isolating factors governing the mechanical strength as a function of coating chemistry, defect morphology, and variables associated with the loading path. We consider two representative materials, including γ-Al 2 O 3 and γ-Al(OH) 3 (i.e., oxide and hydroxide), both of which form readily as aluminum surface coatings. Our results indicate that defects have a significant bearing on the strength of aluminum oxide, with grain boundaries serving to reduce the strain at failure from ε zz = 0.300 to 0.219, relative to perfect single crystal. Our simulations also predict that porosity lowers the elastic stiffness and yield strength of the oxide. Relative to perfect crystal, we find porosity factors of 5%, 10% and 20% decrease the yield stress by 26%, 36% and 53%, respectively. MD predicts that perfect hydroxide and oxide single crystal have respective strains at failure of 0.08 and 0.31 under tensile uniaxial strain loading, and that the corresponding yield stresses are respectively 1.6 and 11.1 GPa. These data indicate that the hydroxide is substantially more susceptible to mechanical failure than the oxide. Our results, coupled with literature findings that indicate hot and humid conditions favor formation of hydroxide and defective oxide coatings, indicate the potential for a complicated dependence of aluminum corrosion susceptibility and stress corrosion cracking on aging history.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Structural complexity of snapshots of two-dimensional Fermi-Hubbard systems

The development of quantum gas microscopy for two-dimensional optical lattices has provided an unparalleled tool to study the Fermi-Hubbard model (FHM) with ultracold atoms. Spin-resolved projective measurements, or snapshots, have played a significant role in quantifying correlation functions, theory verification, and thus the uncovering of underlying physical phenomena such as antiferromagnetism at commensurate filling on bipartite lattices and other charge and spin correlations, as well as dynamical properties at various densities. Here we employ a recent concept, the multiscale structural complexity, and show that when computed for the snapshots (of single spin species, local moments, or total density) it can provide a theory-free property, immediately accessible to experiments. Specifically, after benchmarking results for Ising and $XY$ models, we study the structural complexity of snapshots of the repulsive FHM in the two-dimensional square lattice as a function of doping and temperature. We generate projective measurements using determinant quantum Monte Carlo and compare their complexities against those from the experiment. We demonstrate that these complexities are linked to relevant physical observables such as the entropy and double occupancy. Their behaviors capture the development of correlations and relevant length scales in the system. Furthermore, we provide an open-source code in python which can be implemented into data analysis routines in experimental settings for the square lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Long-range hydrogen-binding effects of carbide interfaces in iron

A micromechanics model was developed to evaluate the elastic binding energy between carbide precipitates and hydrogen interstitials using Eshelby's equivalent inclusion method. In this work, density functional theory (DFT) simulations were performed to obtain the material-specific quantities, e.g., lattice constants and the elastic constants, for the continuum model. Using this model, we find that for coherent carbide precipitates, hydrogen atoms are more likely to bind on the broad surfaces of the disk-like precipitates, which is consistent with experimental observations. For semicoherent and incoherent precipitates, our model suggests that it is possible for semicoherent precipitates to have significant hydrogen binding capability while there is no hydrogen-binding capability of incoherent precipitates, which also agrees with experimental findings. In addition, several factors that influence the binding energies between hydrogen atoms and carbide precipitates were quantitatively analyzed, including the precipitate size, morphology, orientation, and interface. These collective results include both the position and the value of the strongest hydrogen-binding interaction for a wide range of carbide stoichiometries, which contributes to our understanding of hydrogen trapping in steel-based materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Design and Manufacturing of High Performance, Reduced Charge Heat Exchangers

Air-to-fluid heat exchangers (HXs) play a critical role as the main heat transfer component in Heating, Ventilation, Air-Conditioning, and Refrigeration (HVAC&R) systems. However, their airside thermal resistance significantly inhibits their overall performance. Furthermore, these HXs must be continually more compact to meet the latest refrigerant charge limits to reduce emissions. Recent literature suggests that traditional HX geometries (e.g., round or flat tubes with fins), have reached their limits, and more sophisticated shape- and topology-optimized designs are required to achieve the next jump in performance. This research sheds light on the next generation of air-to-refrigerant HXs and aims to address several practical issues to commercialization such as novelty challenges (improved performance for significant charge reduction; modeling expertise & time investment), manufacturing challenges (non-round tube manufacturing; tube-header integration; product qualification, e.g., burst pressure testing, extreme operational environment, etc.), and operational challenges (flow maldistribution, fouling & wetting, noise & vibration). For example, a >20% improvement on one (or more) HX-level performance metrics (e.g., envelope volume, airside pressure drop, face area, capacity, refrigerant charge, weight, cost, etc.) must be achieved before a HX design is considered for commercialization. We present a new, comprehensive and experimentally validated air-to-refrigerant HX optimization framework with simultaneous thermal-hydraulic performance and mechanical strength considerations for novel, non-round, shape- and topology-optimized tubes capable of optimizing single and two-phase HX designs for any refrigerant choice and performance requirement with significant engineering time savings compared to conventional design practices. The framework was exercised for a wide range of applications and refrigerants, resulting in HXs which achieved greater than 20% improved performance, 20% reductions in size, and 25% reductions in refrigerant charge. To enable non-round tube bundle use in next generation HVAC&R equipment, novel manufacturing techniques were investigated, including the development of conventional manufacturing methods for small diameter, non-round tubes and novel tube-header integration strategies. In total, ten HX prototypes were manufactured, nine using conventional methods directly attributed to this project and one using advanced additive manufacturing methods. The five-year manufacturing feasibility of the proposed HXs was found to have a good outlook. The non-round tube HX simulated performance was validated through comprehensive experimental testing, including nine in-house component-level tests, one independent component-level test at an industry partner laboratory, and in-house system-level tests of using a commercially-available, residential packaged A/C unit which was retrofitted with a non-round tube prototype HX. It was found that HX designs proposed by the new framework can successfully predict experimental thermal-hydraulic performance within ±10-20% the first time with no manual design changes, eliminating the need for time-consuming and expensive prototyping efforts. This work will accelerate design and time to market for next generation HXs while simultaneously facilitating industry transition to new refrigerants at lower charge.

42 ENGINEERING↗
DOE OSTI2014

Aerosol indirect effect on the grid-scale clouds in the two-way coupled WRF–CMAQ: model description, development, evaluation and regional analysis

This study implemented first, second and glaciation aerosol indirect effects (AIE) on resolved clouds in the two-way coupled Weather Research and Forecasting Community Multiscale Air Quality (WRF–CMAQ) modeling system by including parameterizations for both cloud drop and ice number concentrations on the basis of CMAQ-predicted aerosol distributions and WRF meteorological conditions. The performance of the newly developed WRF–CMAQ model, with alternate Community Atmospheric Model (CAM) and Rapid Radiative Transfer Model for GCMs (RRTMG) radiation schemes, was evaluated with observations from the Clouds and the See http://ceres.larc.nasa.gov/ . Earth's Radiant Energy System (CERES) satellite and surface monitoring networks (AQS, IMPROVE, CASTNET, STN, and PRISM) over the continental US (CONUS) (12 km resolution) and eastern Texas (4 km resolution) during August and September of 2006. The results at the Air Quality System (AQS) surface sites show that in August, the normalized mean bias (NMB) values for PM 2.5 over the eastern US (EUS) and the western US (WUS) are 5.3% (-0.1%) and 0.4% (-5.2%) for WRF–CMAQ/CAM (WRF–CMAQ/RRTMG), respectively. The evaluation of PM 2.5 chemical composition reveals that in August, WRF–CMAQ/CAM (WRF–CMAQ/RRTMG) consistently underestimated the observed SO 4 2- by -23.0% (-27.7%), -12.5% (-18.9%) and -7.9% (-14.8%) over the EUS at the Clean Air Status Trends Network (CASTNET), Interagency Monitoring of Protected Visual Environments (IMPROVE) and Speciated Trends Network (STN) sites, respectively. Both configurations (WRF–CMAQ/CAM, WRF–CMAQ/RRTMG) overestimated the observed mean organic carbon (OC), elemental carbon (EC) and and total carbon (TC) concentrations over the EUS in August at the IMPROVE sites. Both configurations generally underestimated the cloud field (shortwave cloud forcing, SWCF) over the CONUS in August due to the fact that the AIE on the subgrid convective clouds was not considered when the model simulations were run at the 12 km resolution. This is in agreement with the fact that both configurations captured SWCF and longwave cloud forcing (LWCF) very well for the 4 km simulation over eastern Texas, when all clouds were resolved by the finer resolution domain. The simulations of WRF–CMAQ/CAM and WRF–CMAQ/RRTMG show dramatic improvements for SWCF, LWCF, cloud optical depth (COD), cloud fractions and precipitation over the ocean relative to those of WRF default cases in August. The model performance in September is similar to that in August, except for a greater overestimation of PM 2.5 due to the overestimations of SO 4 2- , NH 4 + , NO 3 - , and TC over the EUS, less underestimation of clouds (SWCF) over the land areas due to the lower SWCF values, and fewer convective clouds in September. Finally, this work shows that inclusion of indirect aerosol effect treatments in WRF–CMAQ represents a significant advancement and milestone in air quality modeling and the development of integrated emissions control strategies for air quality management and climate change mitigation.

indirect aerosol forcing↗
DOE OSTI2017

An intermediate-scale model for thermal hydrology in low-relief permafrost-affected landscapes

Integrated surface/subsurface models for simulating the thermal hydrology of permafrost-affected regions in a warming climate have recently become available, but computational demands of those new process-rich simu- lation tools have thus far limited their applications to one-dimensional or small two-dimensional simulations. We present a mixed-dimensional model structure for efficiently simulating surface/subsurface thermal hydrology in low-relief permafrost regions at watershed scales. The approach replaces a full three-dimensional system with a two-dimensional overland thermal hydrology system and a family of one-dimensional vertical columns, where each column represents a fully coupled surface/subsurface thermal hydrology system without lateral flow. The system is then operator split, sequentially updating the overland flow system without sources and the one-dimensional columns without lateral flows. We show that the app- roach is highly scalable, supports subcycling of different processes, and compares well with the corresponding fully three-dimensional representation at significantly less computational cost. Those advances enable recently developed representations of freezing soil physics to be coupled with thermal overland flow and surface energy balance at scales of 100s of meters. Furthermore developed and demonstrated for permafrost thermal hydrology, the mixed-dimensional model structure is applicable to integrated surface/subsurface thermal hydrology in general.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Predicted Melt Curve and Liquid Shear Viscosity of RDX up to 30 GPa

Recent grain-scale simulations of HMX and TATB have shown that predictions for hot spot formation in high explosives are particularly sensitive to accurate determinations of the pressure-dependent melt curve and the shear viscosity of the liquid phase. These physics terms are poorly constrained beyond ambient pressure for the explosive RDX. Here we adopt an all-atom modeling approach using molecular dynamics (MD) simulations to predict the melt curve of RDX near to detonation conditions (30 GPa) and determine the shear viscosity of the liquid as a function of temperature and pressure above the melt curve. Phase-coexistence simulations were used to determine the melt curve, which is predicted to vary by almost 1100 K as the pressure increases from 0 GPa to 30 GPa. Equilibrium MD simulations and the Green-Kubo formalism were used to obtain the pressure-temperature-dependent shear viscosity. The shear viscosity of RDX is predicted to be of similar magnitude to the viscosity of TATB at low GPa-range pressures, and to be roughly an order of magnitude lower than the viscosity of HMX. The temperature dependence of the shear viscosity is Arrhenius at a given pressure, and the exponential pre-factor and activation term exhibit a strong, yet complicated, pressure dependence. An empirical pressure-temperature dependent function for RDX shear viscosity is developed that simultaneously captures a wide range of MD predictions while taking an analytic form that extrapolates smoothly beyond the fitted regime. The relative strength of the pressure and temperature dependencies of these two physics terms is found to be of similar magnitude for RDX, HMX, and TATB, which motivates incorporating these results in future RDX grain-scale modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Influence of loading control on strain bursts and dislocation avalanches at the nanometer and micrometer scale

Through three-dimensional discrete dislocation dynamics simulations, we show that by tuning the mode of external loading, the collective dynamics of dislocations undergo a transition from driven avalanches under stress control to quasiperiodic oscillations under strain control. We directly correlate measured intermittent plastic events with internal dislocation activities and collective dynamics. Under different loading modes, the roles of the weakest dislocation source and the defect population trend are significantly different. This finding raises new possibilities of controlling correlated dislocation activities and obtaining a low defect density in nanostructured devices by tuning external constraints. In addition, the effect of machine stiffness comes to light. The statistical analysis of the burst magnitude is revisited and carefully discussed. Self-organized criticality and scale-free statistics of strain bursts are obeyed under stress control. However, this behavior is shown to break down under strain control. Rapid stress drops under pure strain control force truncation of dislocation avalanches, leading to a dynamical transition to quasiperiodic oscillations.

3-dimensional systems↗