Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Observing System 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 379 records · Page 21

Assessing the Atmospheric Response to Subgrid Surface Heterogeneity in the Single-Column Community Earth System Model, Version 2 (CESM2)

Land-atmosphere interactions are central to the evolution of the atmospheric boundary layer and the subsequent formation of clouds and precipitation. Existing global climate models represent these connections with bulk approximations on coarse spatial scales, but observations suggest that small-scale variations in surface characteristics and co-located turbulent and momentum fluxes can significantly impact the atmosphere. Recent model development efforts have attempted to capture this phenomenon by coupling existing representations of subgrid-scale (SGS) heterogeneity between land and atmosphere models. Such approaches are in their infancy and it is not yet clear if they can produce a realistic atmospheric response to surface heterogeneity. Here, we implement a parameterization to capture the effects of SGS heterogeneity in the Community Earth System Model (CESM2), and compare single-column simulations against high-resolution Weather Research and Forecasting (WRF) large-eddy simulations (LESs), which we use as a proxy for observations. The CESM2 experiments increase the temperature and humidity variances in the lowest atmospheric levels, but the response is weaker than in WRF-LES. In part, this is attributed to an underestimate of surface heterogeneity in the land model due to a lack of SGS meteorology, a separation between deep and shallow convection schemes in the atmosphere, and a lack of explicitly represented mesoscale secondary circulations. These results highlight the complex processes involved in capturing the effects of SGS heterogeneity and suggest the need for parameterizations that communicate their influence not only at the surface but also vertically.

54 ENVIRONMENTAL SCIENCES↗

Bridging adsorption behavior of confined CH 4 -CO 2 binary mixtures across scales

An accurate understanding of the competitive adsorption of CH 4 -CO 2 binary mixtures in nano-confined systems is critical for engineering CO 2 storage in shale gas reservoirs. Due to difficulties in making reliable experimental observations in nano-scale, atomistic simulations (ASs), such as the Grand Canonical Monte Carlo (GCMC) method, provide a viable approach to studying the adsorption behavior of confined fluids. ASs are, however, limited in the size of the compositional domain due to the high computational cost. This work proposes a framework that combines AS and the lattice Boltzmann (LB) method to bridge the physics of confined fluids across scales. The Peng–Robinson equation of state (PR-EoS) produces fugacity coefficients, which serve as input for conducting multi-component GCMC simulations. These GCMC simulations explore the competitive adsorption behavior of CH 4 -CO 2 in nano-slits at various composition, pressure, and channel-width conditions. Both components generate adsorption layers with high densities near the walls with CO 2 preferentially adsorbing compared to CH 4 on the organic walls of carbon sheets. At the mesoscale, a pseudopotential model represents the intermolecular forces in multi-component, multiple-relaxation-time LB simulations. The LB simulations are in good agreement with the GCMC results, allowing us to obtain values for tunable LB parameters. We then extend the use of LB to simulate adsorption behavior in complex networks with nano-sized channels. The phase behavior and fluid properties in the complex geometries of nano-channels differ from nano-slits and bulk systems. Furthermore, the bridging of physics from GCMC (microscale) to LB (mesoscale) via the macroscale PR-EoS connects the adsorption behavior of binary systems across scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phase-field modeling of thermally-grown oxide and damage evolution in environmental barrier coatings

Silicon carbide-based ceramic matrix composites protected by environmental barrier coatings (EBCs) present a promising materials solution for next-generation gas turbines. Developing more robust and efficient EBCs is therefore of significant technological importance. During the service in high-temperature oxidative environments, there is a thermally grown oxide (TGO) layer, spontaneously formed in the EBC system. TGO is recognized as a critical factor for the degradation and failure of EBCs, yet the detailed mechanisms of TGO growth and its effect on EBC failure remain unclear. In this study we develop a comprehensive chemo-mechano-phase-field model to simulate growth of the TGO in EBCs, factoring in creep and deformation, and especially the cracking behaviors. The volume expansion due to TGO growth and the resulting large inelastic deformation are addressed by using our recently developed, so-called incremental realization of inelastic deformation (IRID) algorithm, in combination with an adapted Hu-Chen spectral solver for elasticity. Simulations of TGO growth are performed considering different growth modes of TGOs determined mainly by the ratio of oxidant permeability in the topcoat to that in the TGO itself. Large-scale three-dimensional (3D) simulations are performed to model the formation of interconnecting vertical/channel cracks (often called ‘mud cracks’). The simulated crack morphology are in excellent agreement with the experimental observations from the literature. The simulations also provide insights into the cracking of EBCs and its dependence on the structure and constituent properties of the coating system. Furthermore, these results demonstrate the developed damage model can be a useful tool for design of more durable EBCs.

Cracking↗

As a Matter of State: The Role of Thermodynamics in Magnetohydrodynamic Turbulence

Turbulence simulations play a key role in advancing the general understanding of the physical properties of turbulence and in interpreting astrophysical observations of turbulent plasmas. For the sake of simplicity, however, turbulence simulations are often conducted in the isothermal limit. Given that the majority of astrophysical systems are not governed by isothermal dynamics, we aim to quantify the impact of thermodynamics on the physics of turbulence, through varying adiabatic index, γ, combined with a range of optically thin cooling functions. Here, we present a suite of ideal magnetohydrodynamics simulations of thermally balanced stationary turbulence in the subsonic, super-Alfvénic, high ${\beta }_{{\rm{p}}}$ (ratio of thermal to magnetic pressure) regime, where turbulent dissipation is balanced by two idealized cooling functions (approximating linear cooling and free–free emission) and examine the impact of the equation of state by considering cases that correspond to isothermal, monatomic, and diatomic gases. We find a strong anticorrelation between thermal and magnetic pressure independent of thermodynamics, whereas the strong anticorrelation between density and magnetic field found in the isothermal case weakens with increasing γ. Similarly, the linear relation between variations in density and thermal pressure with sonic Mach number becomes steeper with increasing γ. This suggests that there exists a degeneracy in these relations with respect to thermodynamics and Mach number in this regime, which is dominated by slow magnetosonic modes. These results have implications for attempts to infer (e.g.,) Mach numbers from (e.g.,) Faraday rotation measurements, without additional information regarding the thermodynamics of the plasma. However, our results suggest that this degeneracy can be broken by utilizing higher-order moments of observable distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗

Computing diffraction patterns of microstructures from phase-field simulations

The diffraction pattern of a material contains information not only on the crystal structures of its constituting phases, but also on its mesoscale spatial distributions of phases, grains, and ferroelastic, ferroelectric, and ferromagnetic domains. While diffraction patterns from experiments such as X-ray diffraction are presented in the reciprocal or Fourier space, mesoscale microstructure models such as the phase-field method naturally produce real-space images of spatial distribution of chemical composition, structural, and ferroic domains. Although one could rather readily compute the Fourier amplitudes of chemical and structural domain distributions generated by mesoscale simulations, they only contain information about the length scale and alignment of the real-space chemical and structure domains. Therefore, a direct comparison between diffraction experiments and mesoscale microstructure simulations is not possible. Here, we develop a theoretical approach to directly compute the crystal diffraction patterns of microstructures predicted by phase-field simulations. In particular, we consider five representative examples of microstructure patterns involving purely compositional domains, a single pair of tetragonal twin structures, multiple twin variants in a hexagonal system, ferroelectric polar vortices, and polycrystalline grains. The results are compared with previous experimental observations as well as X-ray diffraction experiments performed in the present study. Furthermore, the theoretical framework allows one to directly connect material microstructures and diffraction patterns predicted from phase-field simulations and the corresponding diffraction patterns from experiments, and thus providing guidance to experimental diffraction characterization and interpretation of microstructures.

36 MATERIALS SCIENCE↗

Extended kinetic lattice grand canonical Monte Carlo simulation method for transport of multicomponent ion mixtures through a model nanopore system

Abstract An extended version of the original kinetic lattice grand canonical Monte Carlo simulation method combined with mean field theory (KLGCMC/MF) ( J. Chem. Phys. 2007 , 127, 024706) is presented for the study of transport of multicomponent ion mixtures through a model nanopore. Comparison of the extended KLGCMC/MF (eKLGCMC/MF) simulation results with Poisson–Nernst–Planck (PNP) calculations is also made to confirm the validity of the extended simulation approach. Unlike the original version of KLGCMC/MF simulation method that treats only a binary ionic solution with one cation and one anion species, this extended version can deal with a system that includes ternary ion mixtures. A diffusion probability algorithm is also added to the extended version of the simulation method to describe the inhomogeneous diffusivity of ions that is often observed in the ion permeation through nanopores. Both Legendre and Chebyshev polynomials of the second kind were tested as a basis set for the basis set expansion (BSE) method with which to calculate the reaction field energy in the eKLGCMC/MF simulation. It turned out that the Legendre polynomials perform better than the Chebyshev polynomials, and as a result, the Legendre polynomials were implemented in the current version of eKLGCMC/MF simulation algorithm. The presented eKLGCMC/MF simulation method with new features finds its potential applications in nanopore systems where the correlation between ion species with the same sign of charges plays a key role such as oscillating ion currents or anomalous mole fraction effects.

Choi, Inhyeok↗

The Interplay of Solvation and Polarization Effects on Ion Pairing in Nanoconfined Electrolytes

The nature of ion–ion interactions in electrolytes confined to nanoscale pores has important implications for energy storage and separation technologies. However, the physical effects dictating the structure of nanoconfined electrolytes remain debated. Here we employ machine-learning-based molecular dynamics simulations to investigate ion–ion interactions with density functional theory level accuracy in a prototypical confined electrolyte, aqueous NaCl within graphene slit pores. We find that the free energy of ion pairing in highly confined electrolytes deviates substantially from that in bulk solutions, observing a decrease in contact ion pairing but an increase in solvent-separated ion pairing. These changes arise from an interplay of ion solvation effects and graphene’s electronic structure. Notably, the behavior observed from our first-principles-level simulations is not reproduced even qualitatively with the classical force fields conventionally used to model these systems. The insight provided in this work opens new avenues for predicting and controlling the structure of nanoconfined electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characterizing Ice Nucleating Particles Over the Southern Ocean Using Simultaneous Aircraft and Ship Observations

Abstract Supercooled liquid clouds are ubiquitous over the Southern Ocean (SO), even to temperatures below −20°C, and comprise a large fraction of the marine boundary layer (MBL) clouds. Earth system models and reanalysis products have struggled to reproduce the observed cloud phase distribution and occurrence of cloud ice in the region. Recent simulations found the microphysical representation of ice nucleation and growth has a large impact on these properties, however, measurements of SO ice nucleating particles (INPs) to validate simulations are sparse. This study presents measurements of INPs from simultaneous aircraft and ship campaigns conducted over the SO in austral summer 2018, which include the first in situ observations in and above cloud in the region. Our results confirm recent observations that INP concentrations are uniformly lower than measurements made in the late 1960s. While INP concentrations below and above cloud are similar, higher ice nucleation efficiency above cloud supports model simulations that the dominant INP composition varies with height. Model parameterizations based solely on aerosol properties capture the mean relationship between INP concentration and temperature but not the observed variability, which is likely related to the only modest correlations observed between INPs and environmental or aerosol metrics. Including wind speed in addition to activation temperature in a marine INP parameterization reduces bias but does not explain the large range of observed INP concentrations. Direct and indirect inference of marine INP size suggests MBL INPs, at least during Austral summer, are dominated by particles with diameters smaller than 500 nm.

Meteorology & Atmospheric Sciences↗

An analytical model of “Electron-Only” magnetic reconnection rates

Abstract “Electron-only” reconnection, which is both uncoupled from the surrounding ions and much faster than standard reconnection, is arguably ubiquitous in turbulence. One critical step to understanding the rate in this novel regime is to model the outflow speed that limits the transport of the magnetic flux, which is super ion Alfvénic but significantly lower than the electron Alfvén speed based on the asymptotic reconnecting field. Here we develop a simple model to determine this limiting speed by taking into account the multiscale nature of reconnection, the Hall-mediated electron outflow speed, and the pressure buildup within the small system. The predicted scalings of rates and various key quantities compare well with fully kinetic simulations and can be useful for interpreting the observations of NASA’s Magnetospheric-Multiscale (MMS) mission and other ongoing missions.

Physics↗

Accurately simulating core-collapse self-interacting dark matter halos

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations have motivated current efforts to explain surprisingly compact DM halos. If DM is not collisionless, but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we sought to better understand the numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially, our aim was to better understand the evolution of satellite halos. To do so, we ran SIDM N-body simulations of a low-mass halo in isolation and within an external gravitational potential. The simulation set-up was motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a SIDM kernel size that is too large can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our most highly resolved simulation (N = 5 × 10 7 ) we have made the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, for example that the energy conservation error is better kept well below 1%. This allows simulations to be run of halos that become concentrated enough to explain observations of GD-1-like stellar streams or strong gravitational lensing systems.

dark matter↗

Characterization of the DUNE photodetectors and study of the event burst phenomenon

The Deep Underground Neutrino Experiment (DUNE) is an upcoming neutrino physics experiment that will answer some of the most compelling questions in particle physics and cosmology. The DUNE far detectors employ silicon photomultipliers (SiPMs) to detect light produced by charged particles interacting in a large liquid argon time projection chamber (LArTPC). The SiPMs are photosensors consisting of an array of single-photon avalanche diodes (SPAD) operating in Geiger mode. Their high sensitivity and dynamic range, as well as the possibility to fill large surfaces with high-granularity sensors, makes them an ideal choice for DUNE. An international consortium of research groups is currently engaged in a systematic comparison of the performances of the SiPM models that have been custom developed for DUNE by two manufacturers. Such detailed studies, which include gain measurements and a structure study of the dark count rate at 77 K, are meant to determine the best choice of the photodetection system for DUNE, as well as characterize the response of the chosen detectors for the DUNE simulation. Moreover, an investigation of a newly observed phenomenon, in which quick bursts of tens of events at close range are collected in individual SiPMs, is being carried out, which potentially impacts the design of future models and their implementation in particle physics experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Random Quantum Circuits

Quantum circuits—built from local unitary gates and local measurements—are a new playground for quantum many-body physics and a tractable setting to explore universal collective phenomena far from equilibrium. These models have shed light on longstanding questions about thermalization and chaos, and on the underlying universal dynamics of quantum information and entanglement. In addition, such models generate new sets of questions and give rise to phenomena with no traditional analog, such as dynamical phase transitions in quantum systems that are monitored by an external observer. Quantum circuit dynamics is also topical in view of experimental progress in building digital quantum simulators that allow control of precisely these ingredients. Randomness in the circuit elements allows a high level of theoretical control, with a key theme being mappings between real-time quantum dynamics and effective classical lattice models or dynamical processes. Many of the universal phenomena that can be identified in this tractable setting apply to much wider classes of more structured many-body dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Lattice Effective Field Theory Simulations of Nuclei

Lattice effective field theory applies the principles of effective field theory in a lattice framework where space and time are discretized. Nucleons are placed on the lattice sites, and the interactions are tuned to replicate the observed features of the nuclear force. Monte Carlo simulations are then employed to predict the properties of nuclear few- and many-body systems. Here, we review the basic methods and several theoretical and algorithmic advances that have been used to further our understanding of atomic nuclei.

nuclear lattice effective field theory↗

Understanding the Roles of Convective Trigger Functions in the Diurnal Cycle of Precipitation in the NCAR CAM5

The wrong diurnal cycle of precipitation is a common weakness of current global climate models (GCMs). To improve the simulation of the diurnal cycle of precipitation and understand what physical processes control it, we test a convective trigger function described in Xie et al. (2019) with additional optimizations in the NCAR Community Atmosphere Model version 5 (CAM5). The revised trigger function consists of three modifications: 1) replacing the Convective Available Potential Energy (CAPE) trigger with a dynamic CAPE (dCAPE) trigger, 2) allowing convection to originate above the top of planetary boundary layer (i.e., the unrestricted air parcel launch level - ULL), and 3) optimizing the entrainment rate and threshold value of the dynamic CAPE generation rate for convection onset based on observations. Results from 1°-resolution simulations show that the revised trigger can alleviate the long-standing GCM problem of too early maximum precipitation during the day and missing the nocturnal precipitation peak that is observed in many regions, including the U.S. southern Great Plains (SGP). The revised trigger also improves the simulation of the propagation of precipitation systems downstream of the Rockies and the Amazon region. A further composite analysis over the SGP unravels the mechanisms through which the revised trigger affects convection. Additional sensitivity tests show that both the peak time and the amplitude of the diurnal cycle of precipitation are sensitive to the entrainment rate and dCAPE threshold values.

54 ENVIRONMENTAL SCIENCES↗

Role of Tropical Cyclones in Determining ENSO Characteristics

El Niño-Southern Oscillation (ENSO) can effectively modulate global tropical cyclone (TC) activity, but the role TCs may play in determining ENSO characteristics remains unclear. Here we investigate the impact of TC winds on ENSO using a suite of Earth system model experiments where we insert TC winds, extracted from a TC-permitting high-resolution simulation, into a low-resolution model configuration with nearly no intrinsic TCs. The presence of TC winds in the model increases ENSO power and shifts ENSO frequency closer to what we observe. TCs lead to an increase of strong to extreme El Niño events seen in observations and not simulated in the low-resolution model without intrinsic TCs, mainly through enhanced zonal advection feedback and thermocline feedback. Our results indicate that TCs play a fundamental role in producing the ENSO characteristics we experience today in the climate system and point to a two-way climatological interaction between TCs and ENSO.

54 ENVIRONMENTAL SCIENCES↗

The GFDL Variable-Resolution Global Chemistry-Climate Model for Research at the Nexus of US Climate and Air Quality Extremes

We present a variable-resolution global chemistry-climate model (AM4VR) developed at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL) for research at the nexus of US climate and air quality extremes. AM4VR has a horizontal resolution of 13 km over the US, allowing it to resolve urban-to-rural chemical regimes, mesoscale convective systems, and land-surface heterogeneity. With the resolution gradually reducing to 100 km over the Indian Ocean, we achieve multi-decadal simulations driven by observed sea surface temperatures at 50% of the computational cost for a 25-km uniform-resolution grid. In contrast with GFDL's AM4.1 contributing to the sixth Coupled Model Intercomparison Project at 100 km resolution, AM4VR features much improved US climate mean patterns and variability. In particular, AM4VR shows improved representation of: precipitation seasonal-to-diurnal cycles and extremes, notably reducing the central US dry-and-warm bias; western US snowpack and summer drought, with implications for wildfires; and the North American monsoon, affecting dust storms. AM4VR exhibits excellent representation of winter precipitation, summer drought, and air pollution meteorology in California with complex terrain, enabling skillful prediction of both extreme summer ozone pollution and winter haze events in the Central Valley. AM4VR also provides vast improvements in the process-level representations of biogenic volatile organic compound emissions, interactive dust emissions from land, and removal of air pollutants by terrestrial ecosystems. We highlight the value of increased model resolution in representing climate–air quality interactions through land-biosphere feedbacks. AM4VR offers a novel opportunity to study global dimensions to US air quality, especially the role of Earth system feedbacks in a changing climate.

54 ENVIRONMENTAL SCIENCES↗

The interplay of excitonic delocalization and vibrational localization in optical lineshapes: A variational polaron approach

The dynamics of molecular excitonic systems are complicated by a competition between electronic coupling (which drives delocalization) and vibrational-electronic (vibronic) interactions (which tend to encourage electronic localization). A particular challenge of molecular systems is that they typically possess a large number of independent vibrations, with frequencies often spanning the entire spectrum of relevant electronic energy gaps. Recent spectroscopic observations and numerical simulations on a water-soluble chlorophyll-binding protein (WSCP) reveal a transition between two regimes of vibronic behavior, a Redfield-like regime in which low-frequency vibrations respond to a delocalized excitonic state, and a Förster-like regime where high-frequency vibrations act as incoherent excitations on individual pigments. Although numerical simulations can reproduce these effects, there is a need for a simple, systematic theory that accurately describes the smooth transition between these two regimes in experimental spectra. Here we address this challenge by generalizing the variational polaron transform approach of [Bloemsma et al., Chem. Phys. 481, 250 (2016)] to include arbitrary bath densities for systems with or without symmetry. We benchmark this theory against both numerical matrix-diagonalization methods and experimental 77 K fluorescence spectra for two WSCP variants, obtaining quite satisfactory agreement in both cases. Here, we apply this theory to offer an explanation for the large loss in apparent electronic coupling in the WSCP Q57K mutant and to examine the likely impact of the interplay between excitonic delocalization and vibrational localization on vibrational sideband shapes and apparent coupling strengths in high-resolution optical spectra for chlorophyll-protein complexes such as WSCP.

14 SOLAR ENERGY↗

Impacts of sea ice leads on sea salt aerosols and atmospheric chemistry in the Arctic

Abstract. The processes contributing to Arctic cold-season (November–April) sea salt aerosols (SSAs) remain uncertain. Observations from coastal Alaska suggest that emissions from open leads in sea ice, which are not included in climate models, may play a dominant role. Their Arctic-wide significance has not yet been quantified. Here, we create an emission parameterization of SSAs from leads by combining satellite data of lead area (the Advanced Microwave Scanning Radiometer–Earth Observation System (AMSR-E) product) and a chemical transport model (GEOS-Chem) to quantify pan-Arctic SSA emissions from leads during the cold season from 2002 to 2008 and to predict their impacts on atmospheric chemistry, evaluating the results of our simulated SSAs against in situ observations. The AMSR-E product detects large leads with certainty (> 3 km in size), and, hence, our study is limited to quantifying emissions from large leads. Lead emissions vary seasonally and interannually. Simulated total monthly SSA emissions increase by 1.1 %–1.8 % (≥60° N latitude) and 5.6 %–7.5 % (≥75° N) for the 2002–2008 cold seasons. SSA concentrations primarily increase at the location of leads, where standard model concentrations are low. GEOS-Chem overestimates SSA concentrations at Arctic sites compared to ground observations, even when lead emissions are not included, suggesting underestimation of SSA sinks and/or uncertainties in SSA emissions from blowing snow and the open ocean. Multi-year monthly mean surface bromine atom (Br) concentrations increase by 2.8 %–8.8 % due to SSAs from leads for the 2002–2008 cold seasons. Changes in ozone concentrations are negligible. While leads contribute < 10 % to Arctic-wide SSA emissions in the years 2002–2008, these emissions occur in regions of low background aerosol concentrations. Leads may increase in frequency under future climate change, which could increase SSA emissions from leads.

Emme, Erin J. (ORCID:0000000268406352)↗