Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “supercell modeling”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Inertial Oscillation of a Vertical Rotating Draft with Application to a Supercell Storm: Video Supplement to NASA-TP-3230

In this video (8 min., color, sound, VHS), animation depicts the inertial oscillation of a new mathematical model ('vertical rotating draft') for spinning up a single supercell storm. The oscillation consists of a long quiescent phase when the draft is large in diameter and rotates anticyclonically and a short intense phase when the draft is small and cyclonic. During the intense phase, the rotating draft resembles a supercell. The physical basis for the oscillation is depicted by tracking air parcels in the draft as they move along inertial circles (projected on a horizontal plane), where the horizontal pressure gradient is zero and the Coriolis force balances the centrifugal force. A side view of the oscillation shows that contraction and expansion are linked, respectively, to buoyantly driven compressible downdraft and updraft. An aerial view tracks the draft as it moves above the surface of the Earth and turns to the right during the intense phase. Radar echoes from a supercell storm are superimposed for comparison. The data appear to support only the intense phase. A critical experiment would measure the predominantly downward flow that theoretically occurs before the right turn in a supercell track and causes contraction and spin-up.

Supercell thunderstorm↗
DOE OSTI2026

Machine learning approach for vibronically renormalized electronic band structures

Here, we present a machine learning (ML) method for efficient computation of vibrational thermal expectation values of physical properties from first principles. Our approach is based on the nonperturbative frozen phonon formulation in which stochastic Monte Carlo algorithm is employed to sample configurations of nuclei in a supercell at finite temperatures based on a first-principles phonon model. A deep-learning neural network is trained to accurately predict physical properties associated with sampled phonon configurations, thus bypassing the time-consuming ab initio calculations. To incorporate the point-group symmetry of the electronic system into the ML model, group-theoretical methods are used to develop a symmetry-invariant descriptor for phonon configurations in the supercell. We apply our ML approach to compute the temperature dependent electronic energy gap of silicon based on density functional theory (DFT). We show that, with less than a hundred DFT calculations for training the neural network model, an order of magnitude larger number of sampling can be achieved for the computation of the vibrational thermal expectation values. Our work highlights the promising potential of ML techniques for finite temperature first-principles electronic structure methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Quantitative surface studies of moiré intercalated Pb structures under Gr/SiC

Quantitative growth studies using surface diffraction and scanning tunneling microscopy are carried out to determine the intercalation conditions for Pb under Gr/SiC, with initial phase a mixture of buffer-layer and single-layer graphene. Intercalation is a kinetically complex, multistage process with intercalated Pb moving between different subsurface locations with annealing, eventually to the location of higher stability. The transfer is determined by monitoring the difference of a given spot intensity from its clean value, as a function of growth parameters. In addition, a broad, low-intensity (10×10) gr moiré phase is observed exclusively originating from intercalated Pb. Well-defined growth conditions are identified for the moiré spot to be observed reversibly. The (10×10) gr phase follows changes of the Pb(10) spot during growth and annealing experiments, instead of changes of the other Gr/SiC spots indicating preferred nucleation of Pb islands on the (10×10) gr domains. These structural results are discussed in connection to different models proposed for the (10×10) gr supercell and the subsurface Pb coverage. In conclusion, mass transport of subsurface Pb is unusually high, possibly because of the very high compression of the intercalated Pb, when compared to metallic Pb(111).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1987

The terminal area simulation system. Volume 2: Verification cases

The numerical simulation of five case studies are presented and are compared with available data in order to verify the three-dimensional version of the Terminal Area Simulation System (TASS). A spectrum of convective storm types are selected for the case studies. Included are: a High-Plains supercell hailstorm, a small and relatively short-lived High-Plains cumulonimbus, a convective storm which produced the 2 August 1985 DFW microburst, a South Florida convective complex, and a tornadic Oklahoma thunderstorm. For each of the cases the model results compared reasonably well with observed data. In the simulations of the supercell storms many of their characteristic features were modeled, such as the hook echo, BWER, mesocyclone, gust fronts, giant persistent updraft, wall cloud, flanking-line towers, anvil and radar reflectivity overhang, and rightward veering in the storm propagation. In the simulation of the tornadic storm a horseshoe-shaped updraft configuration and cyclic changes in storm intensity and structure were noted. The simulation of the DFW microburst agreed remarkably well with sparse observed data. The simulated outflow rapidly expanded in a nearly symmetrical pattern and was associated with a ringvortex. A South Florida convective complex was simulated and contained updrafts and downdrafts in the form of discrete bubbles. The numerical simulations, in all cases, always remained stable and bounded with no anomalous trends.

Proctor, F. H.↗
DOE OSTI2025

Classical signatures of quenched and thermal disorder in the dynamics of correlated spin systems

Neutron scattering is frequently used to look for evidence of features indicative of quantum-entangled phases of matter such as continua from fractionalisation or quantised excitations. However, the non-specificity of these features and difficulty of both fully quantum treatments and semiclassical models of disorder, make the diagnosis of such states problematic. Here, we demonstrate the feasibility of semiclassical treatments of disordered systems for supercells of ~10 000 spins. By examining a number of classically disordered models we show the presence of quantised excitations, broad continua and anomalous damping originating from quenched disorder or large classical degeneracies.

36 MATERIALS SCIENCE↗
NASA NTRS2014

Terminal Area Simulation System User's Guide - Version 10.0

The Terminal Area Simulation System (TASS) is a three-dimensional, time-dependent, large eddy simulation model that has been developed for studies of wake vortex and weather hazards to aviation, along with other atmospheric turbulence, and cloud-scale weather phenomenology. This document describes the source code for TASS version 10.0 and provides users with needed documentation to run the model. The source code is programed in Fortran language and is formulated to take advantage of vector and efficient multi-processor scaling for execution on massively-parallel supercomputer clusters. The code contains different initialization modules allowing the study of aircraft wake vortex interaction with the atmosphere and ground, atmospheric turbulence, atmospheric boundary layers, precipitating convective clouds, hail storms, gust fronts, microburst windshear, supercell and mesoscale convective systems, tornadic storms, and ring vortices. The model is able to operate in either two- or three-dimensions with equations numerically formulated on a Cartesian grid. The primary output from the TASS is time-dependent domain fields generated by the prognostic equations and diagnosed variables. This document will enable a user to understand the general logic of TASS, and will show how to configure and initialize the model domain. Also described are the formats of the input and output files, as well as the parameters that control the input and output.

Switzer, George F.↗
NASA NTRS2008

Extreme Helicity and Intense Convective Towers in Hurricane Bonnie

Helicity was calculated in Hurricane Bonnie (1998) using tropospheric-deep dropsonde soundings from the NASA Convection and Moisture Experiment. Large helicity existed downshear of the storm center with respect to the ambient vertical wind shear. It was associated with veering, semicircular hodographs created by strong, vortex-scale, radial-vertical flow induced by the shear. The most extreme values of helicity, among the largest ever reported in the literature, occurred in the vicinity of deep convective cells in the downshear-left quadrant. These cells reached as high as 17.5 km and displayed the temporal and spatial scales of supercells. Convective available potential energy (CAPE) averaged 861 J/kg downshear, but only about one-third as large upshear. The soundings nearest the deep cells were evaluated using two empirical supercell parameters that make use of CAPE, helicity, and/or shear. These parameters supported the possible existence of supercells as a consequence of the exceptional helicity combined with moderate but sufficient CAPE. Ambient vertical wind shear exceeded 12 m/s for 30 h, yet the hurricane maintained 50 m/s maximum winds. It is hypothesized that the long-lived convective cells enabled the storm to resist the negative impact of the shear. Supercells in large-helicity, curved-hodograph environments appear to provide a useful conceptual model for intense convection in the hurricane core. Helicity calculations might also give some insight into the behavior of vortical hot towers, which share some characteristics with supercells.

Molinari, John↗
DOE OSTI2025

Investigating the effects of local environment on nitrogen vacancies in high-entropy metal nitrides

High-entropy metal nitrides are an important material class in a variety of applications, and the role of nitrogen vacancies is of great importance for understanding their stability and mechanical properties. Here, we study six different high-entropy nitrides with eight different metal species to build a predictive model of the nitrogen-vacancy formation energy. We construct sets of supercells that maximize the number of unique nitrogen environments for a given chemistry, and then use density-functional theory to calculate the energy density for all nitrogen sites, and the vacancy formation energies for the highest, lowest, and a median subset based on the energy densities. The energy density of nitrogen sites correlates with the vacancy formation energies, for binary, ternary, and high-entropy nitrides. A linear regression model predicts the vacancy formation energies using only the nearest-neighbor composition; across our eight metals, we find the largest vacancy formation energies next to Hf, then Zr, Ti, V, Cr, Ta, Nb, and the lowest near Mo. Additionally, we see that binary nitride data show qualitatively similar vacancy formation energy trends for high-entropy nitrides; however, the binary data alone are insufficient to predict the complex nitride behavior. Our model is both predictive and easily interpretable, and correlates with experimental data.

DeSilva, Charith R. [Univ. of Illinois at Urbana-C↗
DOE OSTI2025

Multiscale Modeling Framework Using Element‐Based Galerkin Methods for Moist Atmospheric Limited‐Area Simulations

This paper presents a multiscale modeling framework (MMF) to model moist atmospheric limited-area weather. The MMF resolves large-scale convection using a coarse grid while simultaneously resolving local features through numerous fine local grids and coupling them seamlessly. Both large- and small-scale processes are modeled using the compressible Navier-Stokes equations within the Nonhydrostatic Unified Model of the Atmosphere (NUMA), and are discretized using a continuous element-based Galerkin method (spectral elements) with high-order basis functions. Consequently, the large-scale and small-scale models share the same dynamical core but have the flexibility to be adjusted individually. The proposed MMF method is tested in 2D and 3D idealized limited-area weather problems involving storm clouds produced by squall line and supercell simulations. Numerical results from the MMF showed enhanced representation of cloud processes compared to the coarse model.

Kang, Soonpil [Naval Postgraduate School, Monterey↗
NASA NTRS2016

WRF-Chem Simulations of Lightning-NOx Production and Transport in Oklahoma and Colorado Thunderstorms Observed During DC3

The focus of this analysis is on lightning-generated nitrogen oxides (LNOx) and their distribution for two thunderstorms observed during the Deep Convective Clouds and Chemistry (DC3) field campaign in May-June 2012. The Weather Research and Forecasting Chemistry (WRF-Chem) model is used to perform cloud-resolved simulations for the May 29-30 Oklahoma severe convection, which contained one supercell, and the June 6-7 Colorado squall line. Aircraft and ground-based observations (e.g., trace gases, lightning and radar) collected during DC3 are used in comparisons against the model-simulated lightning flashes generated by the flash rate parameterization schemes (FRPSs) incorporated into the model, as well as the model-simulated LNOx predicted in the anvil outflow. Newly generated FRPSs based on DC3 radar observations and Lightning Mapping Array data are implemented in the model, along with previously developed schemes from the literature. The results of these analyses will also be compared between storms to investigate which FRPSs were most appropriate for the two types of convection and to examine the variation in the LNOx production. The simulated LNOx results from WRF-Chem will also be compared against other previously studied mid-latitude thunderstorms.

Lightning↗
NASA NTRS2006

The Effects of Lightning NO(x) Production during the July 21 EULINOX Storm studied with a 3-D Cloud-scale Chemical Transport Model

The July 21,1998 thunderstonn observed during the European Lightning Nitrogen Oxides Project (EULINOX) project was simulated using the three-dimensional Goddard Cumulus Ensemble (GCE) model. The simulation successfully reproduced a number of observed storm features including the splitting of the original cell into a southern cell which developed supercell characteristics, and a northern cell which became multicellular. Output from the GCE simulation was used to drive an offline cloud-scale chemical transport model which calculates tracer transport and includes a parameterization of lightning NO(x) production which uses observed flash rates as input. Estimates of lightning NO(x) production were deduced by assuming various values of production per intracloud and production per cloud-to-ground flash and comparing the results with in-cloud aircraft observations. The assumption that both types of flashes produce 360 moles of NO per flash on average compared most favorably with column mass and probability distribution functions calculated from observations. This assumed production per flash corresponds to a global annual lightning NOx source of 7 Tg N per yr. Chemical reactions were included in the model to evaluate the impact of lightning NO(x), on ozone. During the storm, the inclusion of lightning NOx in the model results in a small loss of ozone (on average less than 4 ppbv) at all model levels. Simulations of the chemical environment in the 24 hours following the storm show on average a small increase in the net production of ozone at most levels resulting from lightning NO(x), maximizing at approximately 5 ppbv per day at 5.5 km. Between 8 and 10.5 km, lightning NO(x) causes decreased net ozone production.

Ott, Lesley E.↗
NASA NTRS1989

Simulated appearance of an evolving supercell hail storm at SSM/I microwave frequencies

The contribution of distinct classes of hydrometeors generated from a cloud model to top-of-atmosphere microwave brightness temperatures is studied at the four separate frequencies of the Special Sensor Microwave Imager (SSM/I). The relative radiative role of individual modes of hydrometeors during the evolution of a severe thunderstorm is examined using a microwave radiative transfer model coupled with a three-dimensional cloud model capable of simulating the essential dynamical and microphysical features of an intense hail store. Consideration is given to plans to use the simulation results with actual multichannel measurements from SSM/I to provide a framework for multispectral precipitation retrieval algorithms.

Smith, Eric A.↗
DOE OSTI2025

Ultrahigh-pressure crystallographic passage towards metallic hydrogen

The structural evolution of molecular hydrogen H 2 under multi-megabar compression and its relation to atomic metallic hydrogen is a key unsolved problem in condensed-matter physics. Although dozens of crystal structures have been proposed by theory, only one, the simple hexagonal-close-packed (hcp) structure of only spherical disordered H 2 , has been previously confirmed in experiments. Through advancing nano-focused synchrotron X-ray probes, here we report the observation of the transition from hcp H 2 to a post-hcp structure with a six-fold larger supercell at pressures above 212 GPa, indicating the change of spherical H 2 to various ordered configurations. Theoretical calculations based on our XRD results found a time-averaged structure model in the space group $P\bar{6}2c$ with alternating layers of spherically disordered H 2 and new graphene-like layers consisting of H 2 trimers (H 6 ) formed by the association of three H 2 molecules. Here, this supercell has not been reported by any previous theoretical study for the post-hcp phase, but is close to a number of theoretical models with mixed-layer structures. The evidence of a structural transition beyond hcp establishes the trend of H 2 molecular association towards polymerization at extreme pressures, giving clues about the nature of the molecular-to-atomic transition of metallic hydrogen. Considering the spectroscopic behaviours that show strong vibrational and bending peaks of H 2 up to 400 GPa, it would be prudent to speculate the continuation of hydrogen molecular polymerization up to its metallization.

Ji, Cheng [Center for High Pressure Science and Te↗
DOE OSTI2026

Transition Metal Taggants in UO 2 from First Principles

The incorporation of transition metals into nuclear fuel has gained attention both to improve fuel properties and as a possible nuclear forensics tool. Recent experimental studies by Ulrich et al. and Adorno Lopes et al. have investigated Ni and Fe as candidate transition metal dopants for potential nuclear forensics purposes and found that there is minimal alteration to key UO 2 fuel properties. In the present work, we performed density functional theory (DFT) investigations into possible defect structures for Ni and Fe incorporation. The dynamic stability of these defect structures was validated by calculating the phonon density of states. We found that transition metal incorporation likely occurs via substitution at U sites in the fluorite crystal structure with a nearby O vacancy for charge-balancing, which agrees with experimentally proposed structures. Additionally, this defect structure does not cause long-range alterations to the crystal parameters—an important consideration for use as a nuclear fuel taggant. Future needed work involves computational investigation of additional defect concentrations using larger supercells, additional transition metal charge states, and thermal effects. Such investigations can be carried forward into sintering models for a more complete understanding of the suitability of Ni and Fe as fuel taggants.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRSDate not supplied

Machine Learning Approaches for Rare-Earth Silicate Environmental Barrier Coating Thermochemical and Thermomechanical Property Predictions

Environmental barrier coatings (EBCs) are a necessary enabling technology for the transition from superalloys to silicon carbide (SiC) ceramic matrix composites (CMCs) in gas turbine engines for increased efficiency and decreased fuel costs. SiC-based CMCs are prone to oxidation-based degradation in the engine hot section, and rare-earth (RE) silicates are promising candidates for EBCs due to their close thermal expansion match to the composite substrate and oxidation resistance. However, the design of EBCs is hindered by the large chemical space of candidate materials and the difficulty in obtaining material properties for engineering optimization. This is especially difficult as research continues into mixed-cation or “high-entropy” RE silicates. First-principles computational methods such as density functional theory (DFT) are highly effective at calculating material properties to guide coating design but are limited by their computational cost. Atomistic simulations have the potential to both accelerate property calculations and expand the properties able to be calculated due to their lower computational compared to DFT. However, they require interatomic potentials (IAPs) specific to the material system of interest, and, to our knowledge, there are no suitable IAPs for RE silicates. Machine learning (ML) is a promising technique to accelerate material property predictions indirectly by generating IAPs for atomistic simulations or via direct prediction. In this work, we present two ML approaches to accelerate the calculation of RE silicate properties relevant to EBC design: 1) a ML-derived interatomic potential (IAP) for atomistic simulations of yttrium disilicate (Y2Si2O7) from DFT training data, and 2) a neural network (NN) model to directly predict thermochemical properties of RE silicates and oxides directly from easily obtainable unit cell parameters. Classical MD simulations using the IAP yield lattice properties and bond lengths in good agreement with both DFT and experimental results from x-ray diffraction. Thermodynamic properties calculated using the finite-displacement phonon method and quasi-harmonic approximation were orders of magnitude faster than DFT with good agreement to the DFT results. The IAP was also used to calculate properties such as coefficient of thermal expansion (CTE) that require large simulation supercells and are therefore difficult with DFT. The IAP correctly predicted the anisotropic nature of the CTE in three different phases of Y2Si2O7. The NN model predicts constant pressure heat capacity, Cp, orders of magnitude faster than DFT calculations, which can enable its use as a surrogate model for multiscale simulations. The two methods presented in this work demonstrate the utility of ML for accelerating the prediction of RE silicate properties, which can in turn accelerate EBC design and optimization.

machine learning↗
DOE OSTI2026

Machine learning guided prediction of solute segregation at coherent and semi-coherent metal/oxide interfaces

Investigation of semi-coherent metal/oxide interfaces with misfit dislocations using density functional theory (DFT) is computationally intensive to the point of being prohibitive, as it involves several hundreds to many thousands of atoms. In this study, we examined the solute segregation behavior at the Fe/Y 2 O 3 interface—a model interface for cladding applications in nuclear fission reactors—using a combination of DFT calculations and machine learning (ML) approaches. Both coherent and semi-coherent interfaces were considered. ML models were trained on DFT-calculated segregation energies to identify the key chemical, geometric and strain energy related features that govern solute segregation behavior at coherent Fe/Y 2 O 3 interfaces. Furthermore, it was found that ML models when trained on DFT calculated segregation energy of elements at a coherent interface, comprising of about a hundred-atom supercell, can predict the segregation energy of elements at a semi-coherent Fe/Y 2 O 3 interface (with multiple hundreds of atoms) at a fraction of computational cost (1/35th), with an accuracy comparable to DFT calculations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2017

High Ice Water Concentrations in the 19 August 2015 Coastal Mesoconvective System

During August 2015, NASA's DC-8 research aircraft was flown into High Ice Water Content (HIWC) events as part of a three-week campaign to collect airborne radar data and to obtain measurements from microphysical probes. Goals for this flight campaign included improved characterization of HIWC events, especially from an airborne radar perspective. This paper focuses on one of the flight days, in which a coastal mesoscale convective system (MCS) was investigated for HIWC conditions. The system appears to have been maintained by bands of convection flowing in from the Gulf of Mexico. These convective bands were capped by a large cloud canopy, which masks the underlying structure if viewed from an infrared sensing satellite. The DC-8 was equipped with an IsoKinetic Probe that measured ice concentrations of up to 2.3 g m(exp -3) within the cloud canopy of this system. Sustained measurements of ice crystals with concentrations exceeding 1 g m(exp -3) were encountered for up to ten minutes of flight time. Airborne Radar reflectivity factors were found to be weak within these regions of high ice water concentrations, suggesting that Radar detection of HIWC would be a challenging endeavor. This case is then investigated using a three-dimensional numerical cloud model. Profiles of ice water concentrations and radar reflectivity factor demonstrate similar magnitudes and scales between the flight measurements and model simulation. Also discussed are recent modifications to the numerical model's ice-microphysics that are based on measurements during the flight campaign. The numerical model and its updated ice-microphysics are further validated with a simulation of a well-known case of a supercell hailstorm measured during the Cooperative Convective Precipitation Experiment. Differences in HIWC between the continental supercell and the coastal MCS are discussed.

Proctor, Fred H.↗