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 1,531 records · Page 85

Prediction and Predictability of the Madden Julian Oscillation in the NASA GEOS-5 Seasonal-to-Subseasonal System

In this study, we examine the prediction skill and predictability of the Madden Julian Oscillation (MJO) in a recent version of the NASA GEOS-5 atmosphere-ocean coupled model run at at 1/2 degree horizontal resolution. The results are based on a suite of hindcasts produced as part of the NOAA SubX project, consisting of seven ensemble members initialized every 5 days for the period 1999-2015. The atmospheric initial conditions were taken from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), and the ocean and the sea ice were taken from a GMAO ocean analysis. The land states were initialized from the MERRA-2 land output, which is based on observation-corrected precipitation fields. We investigated the MJO prediction skill in terms of the bivariate correlation coefficient for the real-time multivariate MJO (RMM) indices. The correlation coefficient stays at or above 0.5 out to forecast lead times of 26-36 days, with a pronounced increase in skill for forecasts initialized from phase 3, when the MJO convective anomaly is located in the central tropical Indian Ocean. A corresponding estimate of the upper limit of the predictability is calculated by considering a single ensemble member as the truth and verifying the ensemble mean of the remaining members against that. The predictability estimates fall between 35-37 days (taken as forecast lead when the correlation reaches 0.5) and are rather insensitive to the initial MJO phase. The model shows slightly higher skill when the initial conditions contain strong MJO events compared to weak events, although the difference in skill is evident only from lead 1 to 20. Similar to other models, the RMM-index-based skill arises mostly from the circulation components of the index. The skill of the convective component of the index drops to 0.5 by day 20 as opposed to day 30 for circulation fields. The propagation of the MJO anomalies over the Maritime Continent does not appear problematic in the GEOS-5 hindcasts implying that the Maritime Continent predictability barrier may not be a major concern in this model. Finally, the MJO prediction skill in this version of GEOS-5 is superior to that of the current seasonal prediction system at the GMAO; this could be partly attributed to a slightly better representation of the MJO in the free running version of this model and partly to the improved atmospheric initialization from MERRA-2.

Achuthavarier, Deepthi↗

Simulations of the monthly mean atmosphere for February 1976 with the GISS model

Three monthly mean simulations of the global atmosphere were computed for February 1976 with the GISS model from observed initial conditions on the first day of the month. In a replication experiment, two of these computations generated slightly different monthly mean states, apparently due to the schedule of interruptions on the computer. The root-mean-square errors of replication over the Northern Hemisphere were found to be about 2 mb, 20 m and 1 K for sea level pressure, 500 mb height and 850 mb temperature, respectively. The monthly mean 500 mb forecast results for February 1976 over the Northern Hemisphere are consistent with those from earlier GISS model experiments and again indicate some predictive skill at that level. Use of the observed monthly mean sea surface temperature (SST) field for February 1976 in place of the climatological SST field for February resulted in slightly improved simulations over the globe and Northern Hemisphere,

Spar, J.↗

Multiscale Modeling of Solar Coronal Magnetic Reconnection

Magnetic reconnection is widely believed to be the primary process by which the magnetic field releases energy to plasma in the Sun's corona. For example, in the breakout model for the initiation of coronal mass ejections/eruptive flares, reconnection is responsible for the catastrophic destabilizing of magnetic force balance in the corona, leading to explosive energy release. A critical requirement for the reconnection is that it have a "switch-on' nature in that the reconnection stays off until a large store of magnetic free energy has built up, and then it turn on abruptly and stay on until most of this free energy has been released. We discuss the implications of this requirement for reconnection in the context of the breakout model for CMEs/flares. We argue that it imposes stringent constraints on the properties of the flux breaking mechanism, which is expected to operate in the corona on kinetic scales. We present numerical simulations demonstrating how the reconnection and the eruption depend on the effective resistivity, i.e., the effective Lundquist number, and propose a model for incorporating kinetic flux-breaking mechanisms into MHO calculation of CMEs/flares.

Antiochos, Spiro K.↗

Similarity Criteria for the Study of Removal of Spherical Non-metallic Inclusions in Physical Models of Continuous Casting Tundishes: A More Fundamental Approach

Physical modeling is widely used for studying steelmaking processes. The accuracy of this technique depends on the use of appropriate similarity criteria during experimental design. The present work deduced rigorously the similarity criteria for the study of spherical non-metallic inclusions in isothermal physical models of tundishes. Initially, macroscopic quantities were used to simplify the multiphase problem. Then, the similarity criteria were obtained from the dimensionless form of the exact equation of motion for an inclusion in Stokes regime. The validity range of these criteria was discussed and they were compared with common literature practice. Some practical considerations were also presented. This work clarifies various implicit assumptions usually adopted in literature. Additionally, a restrictive criterion to choose the density of the inclusion simulating particles was obtained, which can be satisfied, for example, by hollow glass particles in water models. This criterion is necessary to adequately describe the inclusion motion in regions where fluid elements suffer strong acceleration (e.g., inside nozzles, near flow controllers or in the entering fluid jet). Moreover, it was established an upper limit for the concentration of particles injected in the inlet nozzle of physical models, which do not simulate the phenomenon of inclusion aggregation.

36 MATERIALS SCIENCE↗

Weather Research and Forecasting Model Wind Sensitivity Study at Edwards Air Force Base, CA

This abstract describes work that will be done by the Applied Meteorology Unit (AMU) in assessing the success of different model configurations in predicting "wind cycling" cases at Edwards Air Force Base, CA (EAFB), in which the wind speeds and directions oscillate among towers near the EAFB runway. The Weather Research and Forecasting (WRF) model allows users to choose among two dynamical cores - the Advanced Research WRF (ARW) and the Non-hydrostatic Mesoscale Model (NMM). There are also data assimilation analysis packages available for the initialization of the WRF model - the Local Analysis and Prediction System (LAPS) and the Advanced Regional Prediction System (ARPS) Data Analysis System (ADAS). Having a series of initialization options and WRF cores, as well as many options within each core, creates challenges for local forecasters, such as determining which configuration options are best to address specific forecast concerns. The goal of this project is to assess the different configurations available and determine which configuration will best predict surface wind speed and direction at EAFB.

Watson, Leela R.↗

Observations for Model Intercomparison Project (Obs4MIPs): status for CMIP6

The Observations for Model Intercomparison Project (Obs4MIPs) was initiated in 2010 to facilitate the use of observations in climate model evaluation and research, with a particular target being the Coupled Model Intercomparison Project (CMIP), a major initiative of the World Climate Research Programme (WCRP). To this end, Obs4MIPs (1) targets observed variables that can be compared to CMIP model variables; (2) utilizes dataset formatting specifications and metadata requirements closely aligned with CMIP model output; (3) provides brief technical documentation for each dataset, designed for nonexperts and tailored towards relevance for model evaluation, including information on uncertainty, dataset merits, and limitations; and (4) disseminates the data through the Earth System Grid Federation (ESGF) platforms, making the observations searchable and accessible via the same portals as the model output. Taken together, these characteristics of the organization and structure of obs4MIPs should entice a more diverse community of researchers to engage in the comparison of model output with observations and to contribute to a more comprehensive evaluation of the climate models. At present, the number of obs4MIPs datasets has grown to about 80; many are undergoing updates, with another 20 or so in preparation, and more than 100 are proposed and under consideration. A partial list of current global satellite-based datasets includes humidity and temperature profiles; a wide range of cloud and aerosol observations; ocean surface wind, temperature, height, and sea ice fraction; surface and top-of-atmosphere longwave and shortwave radiation; and ozone (O 3 ), methane (CH 4 ), and carbon dioxide (CO 2 ) products. A partial list of proposed products expected to be useful in analyzing CMIP6 results includes the following: alternative products for the above quantities, additional products for ocean surface flux and chlorophyll products, a number of vegetation products (e.g., FAPAR, LAI, burned area fraction), ice sheet mass and height, carbon monoxide (CO), and nitrogen dioxide (NO 2 ). While most existing obs4MIPs datasets consist of monthly-mean gridded data over the global domain, products with higher time resolution (e.g., daily) and/or regional products are now receiving more attention. Along with an increasing number of datasets, obs4MIPs has implemented a number of capability upgrades including (1) an updated obs4MIPs data specifications document that provides additional search facets and generally improves congruence with CMIP6 specifications for model datasets, (2) a set of six easily understood indicators that help guide users as to a dataset's maturity and suitability for application, and (3) an option to supply supplemental information about a dataset beyond what can be found in the standard metadata. With the maturation of the obs4MIPs framework, the dataset inclusion process, and the dataset formatting guidelines and resources, the scope of the observations being considered is expected to grow to include gridded in situ datasets as well as datasets with a regional focus, and the ultimate intent is to judiciously expand this scope to any observation dataset that has applicability for evaluation of the types of Earth system models used in CMIP.

54 ENVIRONMENTAL SCIENCES↗

A Synoptic System for Capturing Ecosystem Control Points Across Terrestrial‐Aquatic Interfaces

Interconnected landscape features such as terrestrial‐aquatic interfaces play an outsized role in biogeochemical cycles as ecosystem control points, but it is notoriously challenging to characterize these. Here, we document a synoptic sensor network design that is (a) flexible to accommodate diverse ecosystem interfaces and gradients, (b) adaptable to monitoring and modeling needs of small and large projects alike, (c) standardized for intercomparability across sites and field experiments, and (d) adequately replicated to capture heterogeneity of each parameter monitored. This real‐time monitoring of surface water, groundwater, soil, and vegetation supports configuration and evaluation of models that span upland, wetland, open water strata, and transitions between them. We established the network at seven sites along the Chesapeake Bay and Lake Erie coastlines, including large‐scale flood manipulation experiments in both regions. A central design element is “one data logger program to rule them all”—a collection of sensor‐specific modules deployed on 40 loggers controlling ∼2,000 sensors, with the goal of streamlining maintenance, debugging, and reproducible data processing. The network generates ∼6 M observations per month, capturing system dynamics at the broad spatial and fine temporal scales needed to initialize and benchmark models; measurement frequency can be modified remotely to capture events. This network design has also revealed behaviors not represented in Earth system models, such as transient groundwater oxygen pulses. Completely documented and open source, this standardized, flexible, and efficient sensor network design can reduce barriers to understanding environmental changes and ecosystem responses across systems and scales.

Ward, Nicholas D. [Pacific Northwest National Labo↗

Assessing mechanical response of CO 2 storage into a depleted carbonate reef using a site-scale geomechanical model calibrated with field tests and InSAR monitoring data

Geomechanical risks of injection have raised concerns regarding secure CO 2 storage. In this work, a combined monitoring and modeling approach is used to assess the stress changes and surface uplift associated with CO 2 injection into a depleted carbonate reef of the Michigan basin. A site-scale geomechanical model is built by assigning mechanical properties of formations using well-log and experimental data. Gravity load is applied to the model to estimate the vertical component of stress as well as different lateral boundary displacement scenarios to estimate horizontal stresses. We used a poroelastic pressure-dependent model (instead of a linear elastic mechanical earth model) to calibrate initial stresses using hydraulic fracture test data measured at depleted reservoir status. Multi-phase fluid flow-geomechanical simulations are performed to estimate the poroelastic response during (1) primary depletion (2) field-scale CO 2 injection phase (3) a hypothetical forecast scenario in which well bottom hole pressure (BHP) reach 45000 KPa. The predicted surface uplift is less than 1 mm at the end of the field-scale CO 2 injection phase which is in good agreement with Interferometric Synthetic Aperture Radar (InSAR) uplift measurement. Although the InSAR data shows an insignificant uplift, hydromechanical modeling of injection shows that CO 2 injection still causes reservoir deformation emphasizing the role of carbonate overburden and reservoir formation mechanical properties and limited size of reef on diminishing the surface deformation. Modeling indicates poroelastic response of caprock matters to estimate uplift. The lower permeability of the top two layers provides additional barrier to large uplift. Also, history of subsidence due to production should be accounted to predict uplift due to a follow up injection correctly. This report shows the significance of combining a calibrated geomechanical model with field measured stresses and monitoring data to be used as a tool to ensure the safety of CO 2 storage.

42 ENGINEERING↗

A knowledge-informed large language model framework for U.S. nuclear power plant shutdown initiating event classification for probabilistic risk assessment

Identifying and classifying shutdown initiating events (SDIEs) is critical for developing shutdown probabilistic risk assessment for nuclear power plants. Existing computational approaches cannot achieve satisfactory performance due to the challenges of unavailable large, labeled datasets, imbalanced event types, and label noise. To address these challenges, we propose a hybrid pipeline that integrates a knowledge-informed machine learning model to prescreen non-SDIEs and a large language model (LLM) to classify SDIEs into four types. In the prescreening stage, we proposed a set of 44 SDIE text patterns that consist of the most salient keywords and phrases from six SDIE types. Text vectorization based on the SDIE patterns generates feature vectors that are highly separable by using a simple binary classifier. The second stage builds Bidirectional Encoder Representations from Transformers (BERT)-based LLM, which learns generic English language representations from self-supervised pretraining on a large dataset and adapts to SDIE classification by fine-tuning it on an SDIE dataset. The proposed approaches are evaluated on a dataset with 10,928 events using precision, recall ratio, F 1 score, and average accuracy. In conclusion, the results demonstrate that the prescreening stage can exclude more than 97% non-SDIEs, and the LLM achieves an average accuracy of 95.1% for SDIE classification.

99 - GENERAL AND MISCELLANEOUS↗

Saturn's icy satellites - Thermal and structural models

Thermal history models which assume formation as homogeneous ice-silicate mixtures are constructed for the small, icy Saturnian satellites Mimas, Tethys, Dione, Rhea, and Iapetus, including the effects of radiogenic and accretional heating, conductive and subsolidus convective heat transfer, and lithosphgeric growth. Accretional heating is not likely to have melted the water ice in the interiors of these bodies, and solid state creep of the ice-dominated material precludes melting by radiogenic heating. The four largest satellites are sufficiently large and rich in heat-producing silicates to possess a solid state convection system beneath a rigid lithosphere, irrespective of initial conditions. The model thermal histories are qualitatively consistent with the present appearances of these satellites.

Ellsworth, K.↗

Duke Energy Carbon-Free Resource Integration Study: Capacity Expansion Findings and Production Cost Modeling Plan [Slides]

This presentation highlights initial results from the NREL's Low-Carbon Integration study for Duke Energy. Specifically the presentation focuses on preliminary findings from the ReEDS capacity expansion model for the Carolinas. In addition, the presentation also introduces initial plans for the production cost modeling stage of the work.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Nucleosynthesis of an 11.8 M ⊙ Supernova with 3D Simulation of the Inner Ejecta: Overall Yields and Implications for Short-lived Radionuclides in the Early Solar System

Based on a 3D supernova simulation of an 11.8 M ⊙ progenitor model with an initial solar composition, we study the nucleosynthesis using tracers covering the innermost 0.1 M ⊙ of the ejecta. These ejecta are mostly proton-rich and contribute significant amounts of 45 Sc and 64 Zn. The production of heavier isotopes is sensitive to the electron fraction and hence the neutrino emission from the proto-neutron star. The yields of these isotopes are rather uncertain due to the approximate neutrino transport used in the simulation. In order to obtain the total yields for the whole supernova, we combine the results from the tracers with those for the outer layers from a suitable 1D model. Using the yields of short-lived radionuclides (SLRs), we explore the possibility that an 11.8 M ⊙ supernova might have triggered the formation of the solar system and provided some of the SLRs measured in meteorites. In particular, we discuss two new scenarios that can account for at least the data on 41 Ca, 53 Mn, and 60 Fe without exceeding those on the other SLRs.

Astronomy & Astrophysics↗

The Reanalysis for Stratospheric Trace Gas Studies

The "Reanalysis for Stratospheric Trace Gas Studies" (ReSTS) project intends to provide a state-of-the-art meteorological dataset for use in chemistry-transport models (CTMs). The initial focus is on the period from May 1991 until April 1995, which encompasses the initial few years of observations from NASA's Upper Atmosphere Research Satellite (UARS) as well as the abrupt increase and slow decay of volcanic aerosol loading from Mount Pinatubo. The reanalysis is being performed using NASA's operational "GEOS-4" data assimilation system. The results presented in this paper will begin with an overview of the meteorology of the ReSTS period, validating the ReSTS data against independent analyses, followed by an evaluation of the system statistics (the "Observation minus Forecast" for different input data types). Some sensitivities to the assumptions made in the assimilation process will be presented, including aspects of the forecast model and the statistical model used in the observation-forecast merging. Aspects of trace-gas transport will be discussed, focussing on the on-line water vapor distribution (which is a quasi-inert trace gas in the stratosphere), off-line transport studies of ozone and selected long-lived species, as well as assimilated ozone. The results will be discussed in the context of transport studies and our ability to use reanalyzed datasets to examine long-term variations in the meteorology and composition of the middle atmosphere.

Pawson, Steven↗

Multivariate Dynamical Modeling to Investigate Human Adaptation to Space Flight: Initial Concepts

The array of physiological changes that occur when humans venture into space for long periods presents a challenge to future exploration. The changes are conventionally investigated independently, but a complete understanding of adaptation requires a conceptual basis founded in intergrative physiology, aided by appropriate mathematical modeling. NASA is in the early stages of developing such an approach.

Shelhamer, Mark↗

Multivariate Dynamic Modeling to Investigate Human Adaptation to Space Flight: Initial Concepts

The array of physiological changes that occur when humans venture into space for long periods presents a challenge to future exploration. The changes are conventionally investigated independently, but a complete understanding of adaptation requires a conceptual basis founded in integrative physiology, aided by appropriate mathematical modeling. NASA is in the early stages of developing such an approach.

Shelhamer, Mark↗

A review of thermo-hydro-mechanical modeling of coupled processes in fractured rock: From continuum to discontinuum perspective

Coupled thermo-hydro-mechanical (THM) processes in fractured rock are playing a crucial role in geoscience and geoengineering applications. Diverse and conceptually distinct approaches have emerged over the past decades in both continuum and discontinuum perspectives leading to significant progress in their comprehending and modeling. This review paper offers an integrated perspective on existing modeling methodologies providing guidance for model selection based on the initial and boundary conditions. By comparing various models, one can better assess the uncertainties in predictions, particularly those related to the conceptual models. The review explores how these methodologies have significantly enhanced the fundamental understanding of how fractures respond to fluid injection and production, and improved predictive capabilities pertaining to coupled processes within fractured systems. It emphasizes the importance of utilizing advanced computational technologies and thoroughly considering fundamental theories and principles established through past experimental evidence and practical experience. The selection and calibration of model parameters should be based on typical ranges and applied to the specific conditions of applications. The challenges arising from inherent heterogeneity and uncertainties, nonlinear THM coupled processes, scale dependence, and computational limitations in representing field scale fractures are discussed. Realizing potential advances on computational capacity calls for methodical conceptualization, mathematical modeling, selection of numerical solution strategies, implementation, and calibration to foster simulation outcomes that intricately reflect the nuanced complexities of geological phenomena. Future research efforts should focus on innovative approaches to tackle the hurdles and advance the state-of-the-art in this critical field of study.

Coupling scheme↗

Chaotic behavior and fluid dynamics

It is shown that while fluids are very complicated systems having many degrees of freedom, the basic mechanism underlying the observed chaotic behavior of fluid turbulence may not require the consideration of many degrees of freedom. Simple dynamical systems display erratic and chaotic behavior reminiscent of turbulence, the key to which is the 'sensitivity to initial conditions' introduced by Ruelle (1978). Simple models which are sensitive to initial data and have permitted workers to study chaotic behavior are described, with attention to the Lorenz equations, Landau's idea of a continuous transition to turbulence via an infinite cascade of bifurcations, and one-dimensional maps.

Yorke, J. A.↗

Systematic estimation of forecast and observation error covariances in four-dimensional data assimilation

A two-part algorithm is presented for reliably computing weather forecast model and observational error covariances during data assimilation. Data errors arise from instrumental inaccuracies and sub-grid scale variability, whereas forecast errors occur because of modeling errors and the propagation of previous analysis errors. A Kalman filter is defined as the primary algorithm for estimating the forecast and analysis error convariance matrices. A second algorithm is described for quantifying the noise covariance matrices of any degree to obtain accurate values for the observational error covariances. Numerical results are provided from a linearized one-dimensional shallow-water model. The results cover observational noise covariances, initial instrumental errors and erroneous model values.

Dee, D. P.↗