Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model initialization”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

The WRF Lightning Forecasting Algorithm: Sensitivities to Microphysics and Boundary Layer Physics Schemes

As with WRF model simulated convection in general, LFA output is sensitive to cloud and PBL physics. WRF convection is also sensitive to model initialization uncertainties and unresolved errors.  Large sensitivity exists in HRRR-like Thompson microphysics and PBL, for which original LFA calibration constant needs to be multiplied by 1.13 to give proper FRD amplitudes. Thompson microphysics scheme also shows poorest correlations relative to reference data, suggesting low predictability of HRRR LFA output. Need to validate HRRR LFA against GLM observations.

weather prediction↗
NASA NTRS2020

Lower-Stratospheric Tracer Trends and Variability from Reanalyses

Previous studies have robustly identified a decrease since 1998 in lower stratospheric ozone in the Northern Hemisphere (NH). While this ozone decrease is qualitatively explained as resulting from changes in the large-scale circulation, there is not yet a quantitative mechanistic explanation of these changes. Here we explore the drivers of recent ozone changes using two different configurations of the Goddard Earth Observing System (GEOS) general circulation model. The first configuration of GEOS includes a full chemistry module and is constrained with meteorological fields from the MERRA-2 reanalysis. This configuration (M2GMI) is used to analyze an idealized tracer that covaries closely with ozone on interannual and decadal timescales, revealing that recent ozone decreases in the NH subtropics are associated with a poleward expansion of upwelling in the NH LS, with reduced (enhanced) downwelling over northern subtropics (midlatitudes). The second configuration of GEOS is a free-running version of the GEOS Chemistry-Climate Model (CCM) that is used to perform a ten-member ensemble of free-running simulations. Comparisons of the two configurations reveal that, while the free-running model can produce negative ozone changes in the NH LS, the magnitude of these changes is significantly weaker, relative to both M2GMI and MERRA-2; moreover, these weaker ozone decreases are consistent with weaker simulated changes in the residual circulation. Finally, we examine the GEOS model results in the broader context of the hindcast simulations performed as part of Phase 1 of the Chemistry Climate Modeling Initiative. We show that the majority of the free-running simulations considered here also exhibit weaker long-term residual circulation changes, compared to reanalyses.

Wargan, Kris↗
NASA NTRS2020

Mechanisms Linked to Recent Ozone Decreases in the Northern Hemisphere Lower Stratosphere

Previous studies have robustly identified a decrease since 1998 in lower stratospheric (LS) ozone in the Northern Hemisphere (NH). While this ozone decrease is qualitatively explained as resulting from changes in the large‐scale circulation, there is not yet a quantitative mechanistic explanation of these changes. Here we explore the drivers of recent ozone changes using two different configurations of the Goddard Earth Observing System (GEOS) general circulation model. The first configuration of GEOS includes a full chemistry module and is constrained with meteorological fields from the Modern‐Era Retrospective analysis for Research and Applications, Version 2 (MERRA‐2). This configuration (M2GMI) is used to analyze an idealized tracer that covaries closely with ozone on interannual and decadal timescales, revealing that recent ozone decreases in the NH subtropics are associated with a poleward expansion of upwelling in the NH LS, with reduced (enhanced) downwelling over northern subtropics (midlatitudes). The second configuration of GEOS is a free‐running version of the GEOS Chemistry‐Climate Model (CCM) that is used to perform a ten‐member ensemble of free‐running simulations. Comparisons of the two configurations reveal that, while the free‐running model can produce negative ozone changes in the NH LS, the magnitude of these changes is significantly weaker, relative to both M2GMI and MERRA‐2; moreover, these weaker ozone decreases are consistent with weaker simulated changes in the residual circulation. Finally, we examine the GEOS model results in the broader context of the hindcast simulations performed as part of Phase 1 of the Chemistry Climate Modeling Initiative. We show that the majority of the free‐running simulations considered here also exhibit weaker long‐term residual circulation changes, compared to reanalyses.

Ozone↗
NASA NTRSDate not supplied

Assessing the Viability of Using GEOS-Forecast Product for Landslides Forecasting: A Step Toward Early Warning System

Landslides across the globe are mostly triggered by extreme rainfall events affecting infrastructure, transportation and livelihoods. The risks are rarely quantified due to lack of data, analytical skills and limited modeling techniques. Knowledge of local to global scale landslide risks provides communities and national agencies the ability to adapt disaster management practices to mitigate and recover from these hazards. In order to minimize the risks and improve characterization of community resilience to landslides, it is vital to have reliable information about the factors triggering landslides such as rainfall, well ahead in time. Forecasting potential landslide activity and impacts can be achieved through reliable precipitation forecast models. However, it is challenging because of the temporal and spatial variability of precipitation, an important factor in triggering landslides. Evaluation of the precipitation field, associated errors, and sampling uncertainties is integral for development of efficient and reliable landslide forecasting and early warning system. This study develops a methodology to assess the viability of using a precipitation field provided by a global model and its potential integration in the landslide forecasting system. The study focuses on the comparison between the IMERG (Integrated Multi-satellitE Retrievals for Global Precipitation Mission) and GEOS (NASA Goddard Earth Observing System)-Forecast product over contiguous United States (CONUS). GEOS model assimilates new observations every 6 hours, at 00, 06, 12, and 18 UTC. The framework is tested on the GEOS-Forecast Model initialized at 00 UTC using daily IMERG early product as reference using both categorical and continuous statistics. The categorical statistics includes the probability of detection (POD), success ratio (SR), critical success index (CSI), and the hit bias. Continuous statistics such as correlation, normalized standard deviation, and root-mean-square error are also evaluated. Overall, GEOS-Forecast precipitation field over the analysis period (~1 year) show underestimation with respect to IMERG early for the daily accumulated rainfall. However, the probability distribution function and cumulative distribution function of both show similar patterns. In terms of correlations, POD, SR, CSI, hit bias, the performance varies with respect to the rainfall threshold used.

Sana Khan↗
NASA NTRS2020

Tidal Disruptions of Main Sequence Stars - I. Observable Quantities and their Dependence on Stellar and Black Hole Mass

This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main-sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses(MM0.1510)and six black hole masses( MM105 105BH7). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of Må for Må/Me3. Independent of Må, relativistic effects increase the critical pericenter distance for full disruption events by up to a factor of∼3 relative to the Newtonian prediction. The probability of a direct capture is also independent of Må; at MBH/Me;5×106 this probability is equal to the probability of a complete disruption. The breadth of the debris energy distribution ΔE can differ from the standard estimate by factors of0.35−2, depending on Må and MBH, implying a corresponding change(∝(ΔE)−3/2)in the characteristic mass-return timescale. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly universal relation between the star’s angular momentum and the fraction of Må remaining. Within the “empty loss-cone” regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions.

Taeho Ryu↗
NASA NTRS2020

Tidal Disruptions of Main Sequence Stars - II. Simulation Methodology and Stellar Mass Dependence of the Character of Full Tidal Disruptions

This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main-sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses (0.15 < M_/M<10)and six black hole masses( ́10<MM105 105BH7). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of Må for Må/Me3. Independent of Må, relativistic effects increase the critical pericenter distance for full disruption events by up to a factor of ∼3 relative to the Newtonian prediction. The probability of a direct capture is also independent of Må; at MBH/Me5×106this probability is equal to the probability of a complete disruption. The breadth of the debris energy distribution ΔE can differ from the standard estimate by factors of 0.35−2, depending on Må and MBH, implying a corresponding change(∝(ΔE)−3/2)in the characteristic mass-return timescale. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly universal relation between the star’s angular momentum and the fraction of Må remaining. Within the “empty loss-cone” regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions

Taeho Ryu↗
NASA NTRSDate not supplied

Characterizing Martian X-ray Amorphous Materials through Terrestrial Analogs

X-ray amorphous materials (i.e., lacking long-range crystallographic order) have been identified on the martian surface from orbit and in-situ. Initial models of orbital IR spectral data identified volcanic glasses [e.g., 1], but subsequent interpretations suggested that amorphous silicates are dominantly secondary in nature and formed from water-rock interactions [e.g., 2-4]. ThermalIR spectra from the Mini-Thermal Emission Spectrometer on the Spirit rover show evidence for the amorphous secondary product opaline silica[e.g., 5].X-ray diffraction (XRD) measurements by the CheMin instrument on the Mars Science Laboratory rover have identified 20-70 wt.% amorphous materials in every sample analyzed in Gale crater to date [e.g., 6]. Mass balance calculations using CheMin results and bulk elemental measurements by the Alpha Particle X-ray Spectrometer suggest the amorphous component in Gale crater is variably enriched in Si, Fe, and S[e.g., 7,8].The compositions are not consistent with pure volcanic glass, indicating water must have been involved in the formation of amorphous materials in Gale crater. We seek to better constrain the conditions under which amorphous materials on Mars formed by studying the composition and short-range atomic order of amorphous materials in terrestrial analog environments via XRD and transmission electron microscopy (TEM).

E. B. Rampe↗
NASA NTRSDate not supplied

Characterizing Martian X-Ray Amorphous Materials Through Terrestrial Analogs

X-ray amorphous materials (i.e., lacking long-range crystallographic order) have been identified on the martian surface from orbit and in-situ. Initial models of orbital IR spectral data identified volcanic glasses [e.g., 1], but subsequent interpretations suggested that amorphous silicates are dominantly secondary in nature and formed from water-rock interactions [e.g., 2-4]. Thermal IR spectra from the Mini-Thermal Emission Spectrometer on the Spirit rover show evidence for the amorphous secondary product opaline silica [e.g., 5]. X-ray diffraction (XRD) measurements by the CheMin instrument on the Mars Science Laboratory rover have identified 20-70 wt.% amorphous materials in every sample analyzed in Gale crater to date [e.g., 6]. Mass balance calculations using CheMin results and bulk elemental measurements by the Alpha Particle X-ray Spectrometer suggest the amorphous component in Gale crater is variably enriched in Si, Fe, and S [e.g., 7,8]. The compositions are not consistent with pure volcanic glass, indicating water must have been involved in the formation of amorphous materials in Gale crater. We seek to better constrain the conditions under which amorphous materials on Mars formed by studying the composition and short-range atomic order of amorphous materials in terrestrial analog environments via XRD and transmission electron microscopy (TEM).

E B Rampe↗
NASA NTRS2022

Assimilation of GPM-retrieved Ocean Surface Meteorology Data for Two Snowstorm Events during ICE-POP 2018

As a component of the National Aeronautics and Space Administration (NASA) Weather Focus Area and Global Precipitation Measurement (GPM) Ground Validation participation in the International Collaborative Experiments for PyeongChang 2018 Olympic and Paralympic Winter Games (ICE-POP 2018) field research and forecast demonstration programs, hourly ocean surface meteorology properties were retrieved from the GPM microwave observations for January – March 2018. In this study, the retrieved ocean surface meteorological products – 2-m temperature, 2-m specific humidity, and 10-m wind speed were assimilated into a regional numerical weather prediction (NWP) framework to explore the application of these observations for two heavy snowfall events during the ICE-POP 2018: 27-28 February, and 7-8 March 2018. The Weather Research and Forecasting (WRF) model and the community Gridpoint Statistical Interpolation (GSI) were used to conduct high resolution simulations and data assimilation experiments. The results indicate that the data assimilation has a large influence on surface thermodynamic and wind fields in the model initial condition for both events. With cycled data assimilation, significantly positive influence of the retrieved surface observation was found for the March case with improved quantitative precipitation forecast and reduced error in temperature forecast. A slightly smaller yet positive impact was also found in the forecast of the February case.

assimilation↗
NASA NTRSDate not supplied

TROPOMI Geometry-dependent Lambertian-Equivalent surface Reflectivity (GLER) Product for Improved Trace-Gas Retrieval

Accurate information about the reflectivity of the Earth's surface is required for most satellite retrievals of atmospheric composition, and this information is generally taken from monthly surface reflectivity climatology that neglects angular dependence. Previously we introduced Geometry-dependent Lambertian-equivalent surface reflectivity (GLER) which captures solar and satellite viewing angle dependence as well as daily and seasonal changes. GLER is calculated from simulations of Rayleigh-only top-of-atmosphere (TOA) radiances over non-Lambertian surfaces. We use NASA's Moderate Resolution Imaging Spectroradiometer (MODIS) bidirectional reflectance distribution function (BRDF) product over land and the wind-dependent Cox–Munk wave-facet slope distribution including water-leaving radiance over water to accounts for surface BRDF. We have developed global GLER product, previously for the Ozone Monitoring Instrument (OMI) and recently for Sentinel-5 Precursor (S5P) TROPOspheric Monitoring Instrument (TROPOMI) with several new improvements and updates. We have implemented the near real time daily V006 MODIS MCD43C1 BRDF data and gap-filled with a daily BRDF coefficient climatology created from 2002-2017 V006 MCD43GF data. The NASA’s Global Modelling Initiative hourly 0.25 x 0.25 deg Replay simulations are used for more accurate determination of pixel specific terrain pressure. To improve detection of seasonal snow/ice scenes, we use the 4-km snow cover product from the Interactive Multi-sensor Snow and Ice Mapping System (IMS). Finally, we use an improved version of the vector linearized discrete ordinate radiative transfer (VLIDORT) for update of the top-of-atmosphere (TOA) radiance look-up-tables (LUTs). We demonstrate how the use of GLER is beneficial to TROPOMI’s high spatial resolution (up to 3.5 km x 3.5 km) measurements to monitor atmospheric trace gas pollutants down to the sub-city scale.

TROPOMI↗
NASA NTRSDate not supplied

Verification of TROPOMI NO2 Product Using OMI NO2 algorithm

We evaluated S5P TROPOspheric Monitoring Instrument (TROPOMI) operational nitrogen dioxide (NO2) product by comparing with the NO2 retrievals from Ozone Monitoring Instrument (OMI) onboard NASA’s Aura satellite. We compared spatially matched NO2 vertical column density (VCD) data from OMI and TROPOMI to identify any discrepancies between the two operational products. We also applied OMI Air-Mass Factor (AMF) algorithm to TROPOMI Slant Column Density (SCD) data (v02.03.01) aiming to create long-term merged NO2 Earth Science Data Record (https://disc.gsfc.nasa.gov/datasets/TROPOMI_MINDS_NO2_1.1/summary ). The algorithm (version 4.0) uses new Geometry-dependent Lambertian Equivalent surface Reflectivity (GLER) product available for each TROPOMI pixel. GLER is pre-calculated using the vector linearized discrete ordinate radiative transfer (VLIDORT) model, which uses as input high-resolution bidirectional reflectance distribution function (BRDF) information from NASA's Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) over land and the wind-dependent Cox–Munk wave-facet slope distribution over water, the latter with a contribution from the water-leaving radiance based on MODIS gap-filled in-water chlorophyll-a data. The GLER data, combined with consistently retrieved cloud parameters, provide improved information for the calculation of the scattering weight profiles. Additional AMF improvements result from using high resolution (0.25o latitude x 0.25o longitude) a priori NO2 profile shapes and other auxiliary information from the Global Modeling Initiative (GMI) Replay simulation sampled using S5P orbital simulator. Our TROPOMI_MINDS_NO2 algorithm employs the stratosphere-troposphere separation scheme, de-striping, and surface snow/ice treatment consistently with the latest OMI NO2 re-processing (version 4). We evaluate the TROPOMI_MINDS_NO2 product by comparing with the S5P operational NO2 product as well as with independent NO2 observations from ground-based Pandora and aircraft observations. Our results indicate that the new TROPOMI_MINDS_NO2 retrievals are generally higher than the S5P operational NO2 product over polluted regions and show improved agreement with independent validation data.

TROPOMI↗
NASA NTRS2022

Global Total Ozone Recovery Trends Attributed to Ozone-Depleting Substance (ODS) Changes Derived From Five Merged Ozone Datasets

We report on updated trends using different merged zonal mean total ozone datasets from satellite and ground-based observations for the period from 1979 to 2020. This work is an update of the trends reported in Weber et al. (2018) using the same datasets up to 2016. Merged datasets used in this study include NASA MOD v8.7 and NOAA Cohesive Data (COH) v8.6, both based on data from the series of Solar Backscatter Ultraviolet (SBUV), SBUV-2, and Ozone Mapping and Profiler Suite (OMPS) satellite instruments (1978–present), as well as the Global Ozone Monitoring Experiment (GOME)-type Total Ozone – Essential Climate Variable (GTO-ECV) and GOME-SCIAMACHY-GOME-2 (GSG) merged datasets (both 1995–present), mainly comprising satellite data from GOME, SCIAMACHY, OMI, GOME-2A, GOME-2B, and TROPOMI. The fifth dataset consists of the annual mean zonal mean data from ground-based measurements collected at the World Ozone and Ultraviolet Radiation Data Centre (WOUDC). Trends were determined by applying a multiple linear regression (MLR) to annual mean zonal mean data. The addition of 4 more years consolidated the fact that total ozone is indeed slowly recovering in both hemispheres as a result of phasing out ozone-depleting substances (ODSs) as mandated by the Montreal Protocol. The near-global (60° S–60° N) ODS-related ozone trend of the median of all datasets after 1995 was 0.4 ± 0.2 (2σ) %/decade, which is roughly a third of the decreasing rate of 1.5 ± 0.6 %/decade from 1978 until 1995. The ratio of decline and increase is nearly identical to that of the EESC (equivalent effective stratospheric chlorine or stratospheric halogen) change rates before and after 1995, confirming the success of the Montreal Protocol. The observed total ozone time series are also in very good agreement with the median of 17 chemistry climate models from CCMI-1 (Chemistry-Climate Model Initiative Phase 1) with current ODS and GHG (greenhouse gas) scenarios (REF-C2 scenario). The positive ODS-related trends in the Northern Hemisphere (NH) after 1995 are only obtained with a sufficient number of terms in the MLR accounting properly for dynamical ozone changes (Brewer–Dobson circulation, Arctic Oscillation (AO), and Antarctic Oscillation (AAO)). A standard MLR (limited to solar, Quasi-Biennial Oscillation (QBO), volcanic, and El Niño–Southern Oscillation (ENSO)) leads to zero trends, showing that the small positive ODS-related trends have been balanced by negative trend contributions from atmospheric dynamics, resulting in nearly constant total ozone levels since 2000.

Total Column Ozone Trends↗
NASA NTRSDate not supplied

Coupled Aerosol-Chemistry Simulations of the January 2022 Eruptions of Hunga Tonga-Hunga Ha’apai in the NASA GEOS Earth System Model

The January 2022 eruptions of the underwater Hunga Tonga-Hunga Ha’apai volcano injected more than 100 Tg of water vapor and approximately 0.5 Tg of sulfur dioxide into the stratosphere. Injected materials reached as high as ~55 km in altitude, while the main plume from the eruption travelled west over Australia and the Indian Ocean between about 20 – 30 km altitude. We investigate the transport and the impact of the erupted materials in coupled aerosol-chemistry simulations performed with the NASA Goddard Earth Observing System (GEOS) Earth system model. A sulfur mechanism introduced into the Global Modeling Initiative (GMI) stratospheric-tropospheric chemistry package allows for an interactive simulation of the water vapor-chemistry impacts, and the large water vapor perturbation results in rapid conversion of sulfur dioxide to sulfate aerosol in our aerosol mechanism. We report on the results of a multi-year ensemble of simulations performed with this system that include a control ensemble (no eruption), a water vapor-only injection ensemble, and a water vapor and sulfur dioxide injection ensemble. Simulations of the near-field aerosol and chemistry transport are compared to available measurements from OMPS-LP, OMPS-NM, CALIOP, MLS, and other sensors. The extended impact of the volcanic materials on the stratospheric composition and chemistry over the next several years is further investigated in forecast simulations.

Peter Colarco↗
NASA NTRSDate not supplied

Selecting an Encapsulant for an Aerospace Superconducting Machine

Achieving larger benefits from electrified propulsion in single aisle aircraft (the largest commercial aircraft market segment) necessitates that a significant fraction of their 20 MW or greater propulsion system be electrified. 5 to possibly 20 of the MW-scale electric machines investigated by NASA would be required to meet this need. Approximately 10-20 of NASA’s MW-scale electric machines would be required to reach the necessary power levels. Therefore, electric machines of 5 MW or greater are of interest. However, conventionally cooled high power density MW-scale electric machines still produce a significant amount of waste heat, which makes obtaining these power levels problematic. Highly efficient fully superconducting electric machines (or those with a cryogenically cooled, high purity copper or aluminum stator) are a possible solution but have multiple barriers to be overcome. One such barrier is the durability of superconducting windings. Superconductors are typically brittle, and encapsulants used to hold windings in place are a significant source of stress because of the mismatches in coefficient of thermal expansion at temperatures between room temperature and operational temperatures (20 to 77 K). These encapsulants also form a thermal barrier that can be detrimental to cooling the superconductors. This paper describes an initial study of commercially available encapsulants as it relates to superconducting machine applications along with initial modeling of thermal stress in superconducting windings for a superconducting or cryogenically cooled stator in a 5 MW electric machine.

encapsulant↗
NASA NTRSDate not supplied

Selecting an Encapsulant for an Aerospace Superconducting Machine

Achieving larger benefits from electrified propulsion in single aisle aircraft (the largest commercial aircraft market segment) necessitates that a significant fraction of their 20 MW or greater propulsion system be electrified. These aircraft would require up to 10-20 of today’s state-of-the art machines, which have been demonstrated up to 2.5 MW for aerospace use, and could benefit further from 5 MW+ designs. However, when conventionally cooled, high power density machines still produce a significant amount of waste heat, translating into heavy thermal management and detracting from system benefits. Highly efficient fully superconducting electric machines (or those with a cryogenically cooled, high purity copper or aluminum stator) are a possible solution but have multiple barriers to overcome. One such barrier is the durability of superconducting windings. Superconductors are typically brittle, and encapsulants used to hold windings in place are a significant source of stress due to mismatches in coefficient of thermally induced contraction between ambient and operational temperatures (20 to 77 K). Lack of compliance between the materials at operational temperatures also adds significant stress. These encapsulants also form a thermal barrier that can be detrimental to cooling the superconductors. This work describes an initial study of commercially available encapsulants as it relates to superconducting machine applications, along with initial modeling of thermal stress in superconducting windings for a superconducting or cryogenically cooled stator in a 5 MW electric machine.

encapsulant↗
NASA NTRSDate not supplied

Numerical Predictions of the Flow and Heat Transfer Characteristics in the Film Boiling Regime During Tube Quenching

Cryogenic fluid management plays a major role in refueling of spacecrafts while in space for NASA’s future human space exploration missions. Due to the low boiling points of cryogens, storage, transport and handling of these fluids becomes difficult and may result in inefficient operation of the space propulsion systems. For refueling applications in space, the cryogenic fluids have to be transported across different locations and hence, the transfer of cryogenic fluids through pipes become critical. The cryogenic chill-down process is characterized by different regimes of flow boiling, viz., film boiling, transition boiling and nucleate boiling. The prediction of these regimes in a single CFD framework available in the literature is challenging and the present work attempts to address this challenge by initially modeling the film boiling regime accurately and to incorporate an user-defined function for transition and nucleate boiling at a later stage. Hence, the aim of the present work is to numerically model and validate the film boiling regime of the chilldown curve for liquid nitrogen experiments available in the literature. The validations are carried out at different inlet mass fluxes to have a robust simulation methodology. A dispersed mixture model is used to predict the vapor-liquid interface dynamics with the phase change phenomena modeled using the Lee model.

line chilldown↗
DOE OSTI2019

Development of an Underground Thermal Battery for Enabling Ground Source Heat Pump Applications and Shaping Electric Demand of Buildings

An innovative underground thermal battery (UTB) can provide not only a stable heat sink or heat source for efficient heat pump operation, but also a thermal energy storage to overcome the mismatch between intermittent renewable power supply and fluctuating thermal demands of a building. This article introduces the recent developments of the UTB, including concept design, lab tests, and modeling. Initial results indicate that the UTB has a potential to enable both efficient air-conditioning and active demand side management without sacrificing the desired comfort condition and convenience in the building.

25 ENERGY STORAGE↗
DOE OSTI2025

MSU IETC ML for Modbus (AN EDGE)

This study explores machine learning for decoding Modbus RTU data using K-Nearest Neighbors (KNN) models. An initial KNN model trained on 8,000 packets achieved 95.15% accuracy. Although ML improves generalization, accuracy still falls short of deterministic methods. These findings have implications for Modbus traffic analysis, intrusion detection in industrial networks, and adaptive error correction in real-time monitoring systems. By refining ML-based decoding, future work could enable more efficient anomaly detection and predictive maintenance in industrial automation and cybersecurity applications.

Communication Protocol↗