SEARCH · Engineering Papers

Results for “Coupled models”

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.

SourceResearch connectionFrom yearThrough year

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

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

NASA NTRS2000

TAILSIM Users Guide

The TAILSIM program uses a 4th order Runge-Kutta method to integrate the standard aircraft equations-of-motion (EOM). The EOM determine three translational and three rotational accelerations about the aircraft's body axis reference system. The forces and moments that drive the EOM are determined from aerodynamic coefficients, dynamic derivatives, and control inputs. Values for these terms are determined from linear interpolation of tables that are a function of parameters such as angle-of-attack and surface deflections. Buildup equations combine these terms and dimensionalize them to generate the driving total forces and moments. Features that make TAILSIM applicable to studies of tailplane stall include modeling of the reversible control System, modeling of the pilot performing a load factor and/or airspeed command task, and modeling of vertical gusts. The reversible control system dynamics can be described as two hinged masses connected by a spring. resulting in a fifth order system. The pilot model is a standard form of lead-lag with a time delay applied to an integrated pitch rate and/or airspeed error feedback. The time delay is implemented by a Pade approximation, while the commanded pitch rate is determined by a commanded load factor. Vertical gust inputs include a single 1-cosine gust and a continuous NASA Dryden gust model. These dynamic models. coupled with the use of a nonlinear database, allow the tailplane stall characteristics, elevator response, and resulting aircraft response, to be modeled. A useful output capability of the TAILSIM program is the ability to display multiple post-run plot pages to allow a quick assessment of the time history response. There are 16 plot pages currently available to the user. Each plot page displays 9 parameters. Each parameter can also be displayed individually. on a one plot-per-page format. For a more refined display of the results the program can also create files of tabulated data. which can then be used by other plotting programs. The TAILSIM program was written straightforwardly assuming the user would want to change the database tables, the buildup equations, the output parameters. and the pilot model parameters. A separate database file and input file are automatically read in by the program. The use of an include file to set up all common blocks facilitates easy changing of parameter names and array sizes.

Hiltner, Dale W.↗

NASA NTRS2009

Validation of Community Models: Identifying Events in Space Weather Model Timelines

I develop and document a set of procedures which test the quality of predictions of solar wind speed and polarity of the interplanetary magnetic field (IMF) made by coupled models of the ambient solar corona and heliosphere. The Wang-Sheeley-Arge (WSA) model is used to illustrate the application of these validation procedures. I present an algorithm which detects transitions of the solar wind from slow to high speed. I also present an algorithm which processes the measured polarity of the outward directed component of the IMF. This removes high-frequency variations to expose the longer-scale changes that reflect IMF sector changes. I apply these algorithms to WSA model predictions made using a small set of photospheric synoptic magnetograms obtained by the Global Oscillation Network Group as input to the model. The results of this preliminary validation of the WSA model (version 1.6) are summarized.

MacNeice, Peter↗

NASA NTRS2012

Aerosol Direct, Indirect, Semidirect, and Surface Albedo Effects from Sector Contributions Based on the IPCC AR5 Emissions for Preindustrial and Present-day Conditions

The anthropogenic increase in aerosol concentrations since preindustrial times and its net cooling effect on the atmosphere is thought to mask some of the greenhouse gas-induced warming. Although the overall effect of aerosols on solar radiation and clouds is most certainly negative, some individual forcing agents and feedbacks have positive forcing effects. Recent studies have tried to identify some of those positive forcing agents and their individual emission sectors, with the hope that mitigation policies could be developed to target those emitters. Understanding the net effect of multisource emitting sectors and the involved cloud feedbacks is very challenging, and this paper will clarify forcing and feedback effects by separating direct, indirect, semidirect and surface albedo effects due to aerosols. To this end, we apply the Goddard Institute for Space Studies climate model including detailed aerosol microphysics to examine aerosol impacts on climate by isolating single emission sector contributions as given by the Coupled Model Intercomparison Project Phase 5 (CMIP5) emission data sets developed for Intergovernmental Panel on Climate Change (IPCC) AR5. For the modeled past 150 years, using the climate model and emissions from preindustrial times to present-day, the total global annual mean aerosol radiative forcing is -0.6 W/m(exp 2), with the largest contribution from the direct effect (-0.5 W/m(exp 2)). Aerosol-induced changes on cloud cover often depends on cloud type and geographical region. The indirect (includes only the cloud albedo effect with -0.17 W/m(exp 2)) and semidirect effects (-0.10 W/m(exp 2)) can be isolated on a regional scale, and they often have opposing forcing effects, leading to overall small forcing effects on a global scale. Although the surface albedo effects from aerosols are small (0.016 W/m(exp 2)), triggered feedbacks on top of the atmosphere (TOA) radiative forcing can be 10 times larger. Our results point out that each emission sector has varying impacts by geographical region. For example, the single sector most responsible for a net positive radiative forcing is the transportation sector in the United States, agricultural burning and transportation in Europe, and the domestic emission sector in Asia. These sectors are attractive mitigation targets.

Bauer, Susanne E.↗

NASA NTRS2012

A Canonical Response in Rainfall Characteristics to Global Warming: Projections by IPCC CMIP5 Models

Changes in rainfall characteristics induced by global warming are examined based on probability distribution function (PDF) analysis, from outputs of 14 IPCC (Intergovernmental Panel on Climate Change), CMIP (5th Coupled Model Intercomparison Project) models under various scenarios of increased CO2 emissions. Results show that collectively CMIP5 models project a robust and consistent global and regional rainfall response to CO2 warming. Globally, the models show a 1-3% increase in rainfall per degree rise in temperature, with a canonical response featuring large increase (100-250 %) in frequency of occurrence of very heavy rain, a reduction (5-10%) of moderate rain, and an increase (10-15%) of light rain events. Regionally, even though details vary among models, a majority of the models (>10 out of 14) project a consistent large scale response with more heavy rain events in climatologically wet regions, most pronounced in the Pacific ITCZ and the Asian monsoon. Moderate rain events are found to decrease over extensive regions of the subtropical and extratropical oceans, but increases over the extratropical land regions, and the Southern Oceans. The spatial distribution of light rain resembles that of moderate rain, but mostly with opposite polarity. The majority of the models also show increase in the number of dry events (absence or only trace amount of rain) over subtropical and tropical land regions in both hemispheres. These results suggest that rainfall characteristics are changing and that increased extreme rainfall events and droughts occurrences are connected, as a consequent of a global adjustment of the large scale circulation to global warming.

Lau, William K. M.↗

NASA NTRS2015

A Priori Analysis of a Compressible Flamelet Model using RANS Data for a Dual-Mode Scramjet Combustor

In an effort to make large eddy simulation of hydrocarbon-fueled scramjet combustors more computationally accessible using realistic chemical reaction mechanisms, a compressible flamelet/progress variable (FPV) model was proposed that extends current FPV model formulations to high-speed, compressible flows. Development of this model relied on observations garnered from an a priori analysis of the Reynolds-Averaged Navier-Stokes (RANS) data obtained for the Hypersonic International Flight Research and Experimentation (HI-FiRE) dual-mode scramjet combustor. The RANS data were obtained using a reduced chemical mechanism for the combustion of a JP-7 surrogate and were validated using avail- able experimental data. These RANS data were then post-processed to obtain, in an a priori fashion, the scalar fields corresponding to an FPV-based modeling approach. In the current work, in addition to the proposed compressible flamelet model, a standard incompressible FPV model was also considered. Several candidate progress variables were investigated for their ability to recover static temperature and major and minor product species. The effects of pressure and temperature on the tabulated progress variable source term were characterized, and model coupling terms embedded in the Reynolds- averaged Navier-Stokes equations were studied. Finally, results for the novel compressible flamelet/progress variable model were presented to demonstrate the improvement attained by modeling the effects of pressure and flamelet boundary conditions on the combustion.

Quinlan, Jesse R.↗

NASA NTRS2021

Improved Representation of Atmospheric Dynamics in CMIP6 Models Removes Climate Sensitivity Dependence on Hadley Cell Climatological Extent

The persistent inter-model spread in the response of global mean surface temperature to increased CO. (known as the “Equilibrium Climate Sensitivity”, or “ECS”) is a crucial problem across model generations. This work examines the influence of the models’ present-day atmospheric circulation climatologies, and the accompanying climatological cloud radiative effects, in explaining that spread. We analyze the Coupled Model Intercomparison Project Phase 6 (CMIP6) models and _nd that they simulate a more poleward, and thus more realistic, edge of the Hadley Circulation (HC) in the SouthernHemisphere than the CMIP5 models, although the climatological shortwave cloud radiative effects are similar in the two generations of models. A few CMIP5 models with extreme equatorward biases in the HC edge exhibited high ECS due to strong Southern midlatitude shortwave cloud radiative warming in response to climate change, suggesting an ECS dependence on HC position. We find that such constraint no longer holds for the CMIP6 models, due to absence of models with extreme HC climatologies. In spite of this, however, the CMIP6 models show an increased spread in ECS, with more models in the high ECS range. Additionally, an improved representation of the climatological jet dynamics does not lead to a new emergent constraint in the CMIP6 models either.

Equilibrium Climate Sensitivity↗

NASA NTRS2013

Reducing Spread in Climate Model Projections of a September Ice-Free Arctic

This paper addresses the specter of a September ice-free Arctic in the 21st century using newly available simulations from the Coupled Model Intercomparison Project Phase 5 (CMIP5). We find that large spread in the projected timing of the September ice-free Arctic in 30 CMIP5 models is associated at least as much with different atmospheric model components as with initial conditions. Here we reduce the spread in the timing of an ice-free state using two different approaches for the 30 CMIP5 models: (i) model selection based on the ability to reproduce the observed sea ice climatology and variability since 1979 and (ii) constrained estimation based on the strong and persistent relationship between present and future sea ice conditions. Results from the two approaches show good agreement. Under a high-emission scenario both approaches project that September ice extent will drop to approx. 1.7 million sq km in the mid 2040s and reach the ice-free state (defined as 1 million sq km) in 2054-2058. Under a medium-mitigation scenario, both approaches project a decrease to approx.1.7 million sq km in the early 2060s, followed by a leveling off in the ice extent.

simulation↗

DOE OSTI2021

Plasma-catalyst modeling for materials selection: challenges and opportunities in nitrogen oxidation

Reliable kinetic models are essential to rationalize observations and to guide the selection of appropriate plasma-catalyst combinations. Kinetic models coupling plasma- and heterogeneous-catalytic-chemistries, however, are primitive. Here, we combine reduced plasma chemistries and density functional theory parameterized surface reactions to model NO formation from N2 oxidation, an alternative nitrogen fixation process. We first show plasma excited species can enhance turnover frequencies of both Pt and Au. The enhancing potential is a function of plasma species and their densities. We then compare NO production between plasma/Pt and plasma/Au with an integral reactor model and explore the dependence of NO production as a function of plasma conditions, reactor configurations and chemical compositions. Plasma/Pt generates more NO than plasma/Au under most thermal and plasma parametric conditions. Plasma/Au outperforms plasma/Pt under high plasma excitations at low temperatures. Furthermore, these results highlight that selection of plasma-catalyst combinations should consider the properties of plasmas, materials and the coupling of the two.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Mean-Field Ring Polymer Rates Using a Population Dividing Surface

Mean-field Ring Polymer Molecular Dynamics (MF-RPMD) offers a computationally efficient method for the simulation of reaction rates in multi-level systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, the dividing surface must be chosen to explicitly sample kinked ring polymer configurations where at least one bead is in a different electronic state than the others. Building on this, we introduce a population difference coordinate and a kink-constrained dividing surface, and we test the accuracy of the resulting mean-field rate theory on a series of bilinear coupling model systems. Here, we demonstrate that this new MF-RPMD rate approach is efficient to implement and quantitatively accurate for models over a wide range of driving forces, coupling strengths, and temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Simulation of hardened cement degradation and estimation of uncertainty in predicted failure times with peridynamics

Modeling the degradation of cement-based infrastructure due to aqueous environmental conditions continues to be a challenge. In order to develop a capability to predict concrete infrastructure failure due to chemical degradation, here we created a chemomechanical model of the effects of long-term water exposure on cement paste. The model couples the mechanical static equilibrium balance with reactive–diffusive transport and incorporates fracture and failure via peridynamics (a meshless simulation method). The model includes fundamental aspects of degradation of ordinary Portland cement (OPC) paste, including the observed softening, reduced toughness, and shrinkage of the cement paste, and increased reactivity and transport with water induced degradation. This version of the model focuses on the first stage of cement paste decalcification, the dissolution of portlandite. Given unknowns in the cement paste degradation process and the cost of uncertainty quantification (UQ), we adopt a minimally complex model in two dimensions (2D) in order to perform sensitivity analysis and UQ. We calibrate the model to existing experimental data using simulations of common tests such as flexure, compression and diffusion. Then we calculate the global sensitivity and uncertainty of predicted failure times based on variation of eleven unique and fundamental material properties. We observed particularly strong sensitivities to the diffusion coefficient, the reaction rate, and the shrinkage with degradation. Also, the predicted time of first fracture is highly correlated with the time to total failure in compression, which implies fracture can indicate impending degradation induced failure; however, the distributions of the two events overlap so the lead time may be minimal. Extension of the model to include the multiple reactions that describe complete degradation, viscous relaxation, post-peak load mechanisms, and to three dimensions to explore the interactions of complex fracture patterns evoked by more realistic geometry is straightforward and ongoing.

36 MATERIALS SCIENCE↗

DOE OSTI2023

The Roles of Tidal Marshes in the Estuarine Biochemical Processes: A Numerical Modeling Study

Observations suggest that the existence of tidal marsh can alter the oxygen and nutrient dynamics in adjacent water bodies, but assessing the impacts of large tidal marshes on an estuary is challenging. In this study, we use a modeling approach to investigate the roles of tidal marshes on the estuarine biochemical processes. The marsh model, which simulates the ecological functions of marshes at seasonal and annual time-scales, is embedded inside an unstructured-grid three-dimensional hydrodynamic and eutrophication model (SCHISM-ICM). This modeling system simulates the growth and metabolism of the tidal marshes and links biological processes to nutrient dynamics in the water column and sediment. This model dynamically simulates nutrient recycling and physical transport of the materials between marshes and open water through wetting-drying processes. This coupled model system is validated and successfully applied to the York River Estuary. Model results suggest that tidal marshes influence the local diurnal dissolved oxygen (DO) cycle by exporting dissolved organic carbon and high sediment oxygen demand in the marsh system through the tidal exchange. The high deposition rates of organics and diurnal DO cycle enhance the sediment release of phosphorus. On the other hand, marshes tend to decrease dissolved inorganic nitrogen in the water column by settling particulate nutrients and enhancing the denitrification process. The study demonstrates that tidal marshes exert substantial impacts on the estuarine biochemical processes. The developed tidal marsh model enhances eutrophication modeling and advances the understanding of the feedback effects between marsh biogeochemistry and estuarine eutrophication processes on a systemic scale.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Impacts of climate change on subannual hydropower generation: a multi-model assessment of the United States federal hydropower plant

Abstract Hydropower is a low-carbon emission renewable energy source that provides competitive and flexible electricity generation and is essential to the evolving power grid in the context of decarbonization. Assessing hydropower availability in a changing climate is technically challenging because there is a lack of consensus in the modeling representation of key dynamics across scales and processes. Focusing on 132 US federal hydropower plants, in this study we evaluate the compounded impact of climate and reservoir-hydropower models’ structural uncertainties on monthly hydropower projections. In particular, instead of relying on one single regression-based hydropower model, we introduce another conceptual reservoir operations-hydropower model in the assessment framework. This multi-model assessment approach allows us to partition uncertainties associated with both climate and hydropower models for better clarity. Results suggest that while at least 70% of the uncertainties at the annual scale and 50% at the seasonal scale can be attributed to the choice of climate models, up to 50% of seasonal variability can be attributed to the choice of hydropower models, particularly in regions over the western US where the reservoir storage is substantial. The analysis identifies regions where multi-model assessments are needed and presents a novel approach to partition uncertainties in hydropower projections. Another outcome includes an updated evaluation of Coupled Model Intercomparison Project Phase 5 (CMIP5)-based federal hydropower projection, at the monthly scale and with a larger ensemble, which can provide a baseline for understanding future assessments based on CMIP6 and beyond.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1989

Effects of implementing the Simple Biosphere Model in a general circulation model

The Simple Biosphere Model (SiB) of Sellers et al., was designed to simulate the interactions between the earth's land surface and the atmosphere by treating the vegetation explicitly and realistically, thereby incorporating the biophysical controls on the exchanges of radiation, momentum, sensible and latent heat between the two systems. This paper describes the steps taken to implement SiB in a modified version of the National Meteorological Center's global spectral general circulation model (GCM) and explores the impact of the implementation on the simulated land surface fluxes and near-surface meteorological conditions. The coupled model (SiB-GCM) was used to produce summer and winter simulations. The same GCM was used with a conventional hydrological model (Ctl-GCM) to produce comparable 'control' summer and winter simulations for comparison. It was found that SiB-GCM produced a more realistic partitioning of energy at the land surface than Ctl-GCM. Generally, SiB-GCM, produced more sensible heat flux and less latent heat flux over vegetated land than did Ctl-GCM and this resulted in a much deeper daytime planetary boundary layer and reduced precipitation rates over the continents in SiB-GCM. In the summer simulation, the 200 mb jet stream was slightly weakened in the SiB-GCM relative to the Ctl-GCM results and analyses made from observations.

Sato, N.↗

NASA NTRS2000

Computational Fluid Dynamics Modeling of Nickel Hydrogen Batteries

An electrochemical Ni-H2 battery model has been expanded to include thermal effects. A thermal energy conservation equation was derived from first principles. An electrochemical and thermal coupled model was created by the addition of this equation to an existing multiphase, electrochemical model. Charging at various rates was investigated and the results validated against experimental data. Reaction currents, pressure changes, temperature profiles, and concentration variations within the cell are predicted numerically and compared with available data and theory.

Cullion, R.↗

NASA NTRS2022

Southern Ocean Solar Reflection Biases in Cmip6 Models Linked to Cloud Phase and Vertical Structure Representations

Over the Southern Ocean (SO, 40˚S-70˚S), climate models have consistently underestimated solar reflection. Here we evaluate the relationship between cloud profiles, cloud phase and radiation over the SO in Coupled Model Intercomparison Project Phase 6 (CMIP6) models against Clouds and the Earth's Radiant Energy System (CERES) and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO). We find that the lack of solar reflection is slightly improved in CMIP6 models compared to CMIP5’s, attributable to a better representation of cloud fraction and phase. We show that clouds have a different vertical structure and radiative effect south and north of where the 0˚C isotherm meets the surface (~55˚S). Although the models capture the greater vertical extent of clouds south of 55˚S, they fail to reproduce the observed increase in solar reflection, which we pinpoint to cloud phase biases. Increasing CMIP6 supercooled liquid cloud opacity should help reduce their persistent SW biases.

Southern Ocean↗

DOE OSTI2025

Influence of Lake Ice Biases in Reanalysis Data on Downscaled Climate Simulations over the Great Lakes Region

This data package contains observation-based and model-simulated datasets (all provided in NetCDF format) for evaluating how wintertime lake-ice representation affects regional weather and climate over the Laurentian Great Lakes (freshwater lake ecosystem) during the high–ice-cover winter of 2009. The observational component includes: (1) Stage IV gridded precipitation at 4 km, hourly resolution for January–February 2009 over the Great Lakes region (radar–gauge multisensor precipitation analyses); (2) Great Lakes Surface Environmental Analysis (GLSEA) satellite-derived lake-ice coverage at 1.3 km, daily resolution for the 2009 winter months, providing ice coverage over Lakes Superior, Michigan, Huron, Erie, and Ontario; and (3) in situ measurements at the Standard Rock site on Lake Superior from the Great Lakes Evaporation Network (GLEN) at hourly resolution, including near-surface atmospheric variables and sensible and latent heat fluxes (air–lake exchange) at a fixed point location. The modeling component provides corresponding fields from two simulations, both archived at 4 km, hourly resolution: a standalone Weather Research Forecasting model (WRF) run driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5 (ERA5), and a two-way coupled model using WRF and the Finite Volume Community Ocean Model (WRF-FVCOM, a 3-D hydrodynamic lake model). These outputs include variables relevant to air–lake interaction and lake-effect processes (e.g., near-surface temperature, humidity, wind, precipitation, and surface turbulent fluxes), enabling direct comparison with the observational datasets. Users can analyze and visualize these NetCDF files with common tools such as Python (e.g., xarray, netCDF4, numpy, pandas), NCO/CDO, Panoply, or ncview; NetCDF variables can also be converted to other formats (e.g., CSV, GeoTIFF) using these utilities.

EARTH SCIENCE > ATMOSPHERE > ATMOSPHERIC TEMPERATU↗

NASA NTRS2001

Iron Catalyst Chemistry in High Pressure Carbon Monoxide Nanotube Reactor

The high-pressure carbon monoxide (HiPco) technique for producing single wall carbon nanotubes (SWNT) is analyzed using a chemical reaction model coupled with properties calculated along streamlines. Streamline properties for mixing jets are calculated by the FLUENT code using the k-e turbulent model for pure carbon monixide. The HiPco process introduces cold iron pentacarbonyl diluted in CO, or alternatively nitrogen, at high pressure, ca. 30 atmospheres into a conical mixing zone. Hot CO is also introduced via three jets at angles with respect to the axis of the reactor. Hot CO decomposes the Fe(CO)5 to release atomic Fe. Cluster reaction rates are from Krestinin, et aI., based on shock tube measurements. Another model is from classical cluster theory given by Girshick's team. The calculations are performed on streamlines that assume that a cold mixture of Fe(CO)5 in CO is introduced along the reactor axis. Then iron forms clusters that catalyze the formation of SWNTs from the Boudouard reaction on Fe-containing clusters by reaction with CO. To simulate the chemical process along streamlines that were calculated by the fluid dynamics code FLUENT, a time history of temperature and dilution are determined along streamlines. Alternative catalyst injection schemes are also evaluated.

Scott, Carl D.↗

NASA NTRS2022

The Dynamical Role of the Chesapeake Bay on the Local Ozone Pollution Using Mesoscale Modeling—A Case Study

This study investigated the dynamic influence of the Chesapeake Bay (CB) on local ozone (O3) concentration and distribution using a weather forecasting model. The Weather Research and Forecasting model coupled with Chemistry (WRF–Chem) was employed to simulate O3 production and transportation near the CB. Baseline (water) as well as sensitivity (nowater) model experiments of bay circulation were carried out with and without bay water by changing the water surface from water to land (loam). First, the model performance simulating O3 was evaluated using the baseline experiment results and AirNow surface wind and O3 observations. The results showed that the model overestimates surface O3 by up to 20–30%. Further, the comparisons of the baseline and sensitivity experiments revealed higher O3 mixing ratios, primarily due to the resulting bay breeze circulation. These increases, after considering model overestimation, represent a mean bay dynamics circulation-induced contribution of up to 10% at night and 5% during the day. Furthermore, the boundary layer over northern CB, where it is at its narrowest width, was higher (by 1.2 km on average) during daytime due to higher surface temperatures observed. The boundary layer depth difference between the northern, central, and southern regions of the bay leads to a differential in the role of bay circulation dynamics in the observed O3 increase. The relatively wider swath of water surface over southern CB resulted in a lower boundary layer depth and stronger breeze circulation and this circulation contributed to O3 concentrations. Moreover, since the case selected has a minimal bay breeze circulation, the associated surface ozone enhancements represent what is expected at least at a minimum.

Chesapeak Bay↗