Lattice Spring Model & Simulation Based Parameter Inference
Presentation of mathematical content of SDRD project related to subsidence and shock physics
SEARCH · Engineering Papers
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.
Presentation of mathematical content of SDRD project related to subsidence and shock physics
Current assessment of aerosol radiative effect is hindered by our incomplete knowledge of aerosol optical properties, especially absorption, and our current inability to quantify physical and microphysical processes. In this research, we investigate direct aerosol radiative effect over heavy aerosol loading areas (e.g., Indo-Gangetic Plains, South/East Asia) and its feedbacks on the South Asian climate during the pre-monsoon season (March-June) using the Purdue Regional Climate Model (PRCM) with prescribed aerosol data derived by the NASA Goddard Earth Observing System Model (GEOS-5). Our modeling domain covers South and East Asia (60-140E and 0-50N) with spatial resolutions of 45 km in horizontal and 28 layers in vertical. The model is integrated from 15 February to 30 June 2008 continuously without nudging (i.e., only forced by initial/boundary conditions). Two numerical experiments are conducted with and without the aerosol-radiation effects. Both simulations are successful in reproducing the synoptic patterns on seasonal-to-interannual time scales and capturing a pre-monsoon feature of the northward rainfall propagation over Indian region in early June which shown in Tropical Rainfall Measuring Mission (TRMM) observation. Preliminary result suggests aerosol-radiation interactions mainly alter surface-atmosphere energetics and further result in an adjustment of the vertical temperature distribution in lower atmosphere (below 700 hPa). The modifications of temperature and associated rainfall and circulation feedbacks on the regional climate will be discussed in the presentation.
Explore the source record for details and available documents.
Biomass sorghum (Sorghum bicolor L. Moench) is an annual C 4 grass that has emerged as a candidate bioenergy crop but has not been widely grown in the United States. Corn (Zea mays L.), another annual C 4 grass, has been produced on a large scale in the United States for biofuels. Iowa leads the nation in both corn and ethanol production. The high productivity of corn in Iowa creates the research question: could biomass sorghum be as or more productive than corn in the state in terms of bioenergy? Efforts to use crop models to fill the gaps left by field experimentation on biomass sorghum have also been limited thus far. To address our research question, we collected biophysical data on biomass sorghum grown in Iowa for two growing seasons, used them to develop a biomass sorghum module in an agroecosystem model (Agro-IBIS), and predict the potential performance of biomass sorghum across Iowa compared to maize. Despite dry conditions in 2019 and 2020, average biomass sorghum yields were 17.20 Mg ha –1 . By comparison, average corn aboveground biomass was higher in 2020 (22.02 Mg ha –1 ). Soil cores indicated average belowground biomass of 1.46 Mg ha –1 , with roots concentrated near the surface (73% of biomass above 50 cm in a 1 m core). When biomass sorghum model parameters were calibrated with measured values, model output was in close agreement with measured biomass (slope = 0.932, R 2 = 0.91) and evapotranspiration (slope = 0.757, R 2 = 0.64). Subsequent regional simulations revealed a notable latitudinal gradient in biomass sorghum yield, with a strong linear relationship between yield and seasonal growing degree-days (R 2 = 0.89). When these yields were compared to simulated corn aboveground biomass yields, only 3.4% of the state had biomass sorghum yields that were significantly higher than corn. This number was reduced to 0.3% when comparing the crops in terms of energy ethanol yield. Thus, we conclude biomass sorghum must be improved to be competitive with corn as an annual biofuel crop in Iowa. This study provides a baseline against which to compare advancements in sorghum breeding for biomass and stress tolerance in Iowa.
Growing evidence indicates that the atmospheric and oceanic circulation experiences a systematic poleward shift in a warming climate. However, the complexity of the climate system, including the coupling between the ocean and the atmosphere, natural climate variability and land-sea distribution, tends to obfuscate the causal mechanism underlying the circulation shift. Here, using an idealised coupled aqua-planet model, we explore the mechanism of the shifting circulation, by isolating the contributing factors from the direct CO 2 forcing, the indirect ocean surface warming, and the wind-stress feedback from the ocean dynamics. We find that, in contrast to the direct CO 2 forcing, ocean surface warming, in particular an enhanced subtropical ocean warming, plays an important role in driving the circulation shift. This enhanced subtropical ocean warming emerges from the background Ekman convergence of surface anomalous heat in the absence of the ocean dynamical change. Additionally, it expands the tropical warm water zone, causes a poleward shift of the mid-latitude temperature gradient, hence forces a corresponding shift in the atmospheric circulation and the associated wind pattern. The shift in wind, in turn drives a shift in the ocean circulation. Our simulations, despite being idealised, capture the main features of the observed climate changes, for example, the enhanced subtropical ocean warming, poleward shift of the patterns of near-surface wind, sea level pressure, storm tracks, precipitation and large-scale ocean circulation, implying that increase in greenhouse gas concentrations not only raises the temperature, but can also systematically shift the climate zones poleward.
The NASA Crew Launch Vehicle is a two-stage orbital launcher designed to meet NASA's current as well as future needs for human space flight. In order to free the designers to explore more possibilities during the design phase, a need exists for the ability to quickly perform simulation on both the baseline vehicle as well as the vehicle after proposed changes due to mission planning, vehicle configuration and avionics changes, proposed new guidance and control algorithms, and any other contingencies the designers may wish to consider. Further, after the vehicle is designed and built, the need will remain for such analysis in the event of future mission planning. An easily reconfigurable, modular, nonlinear six-degree-of-freedom simulation matching NASA Marshall's in-house high-fidelity simulator is created with the ability to quickly perform simulation and analysis of the Crew Launch Vehicle throughout the entire launch profile. Simulation results are presented and discussed, and an example comparison fly-off between two candidate controllers is presented.
ABSTRACT Fuzzy dark matter (FDM), comprised of ultralight ($m \sim 10^{-22}\,{\rm eV}$) boson particles, has received significant attention as a viable alternative to cold dark matter (CDM), as it approximates CDM on large scales (${\gtrsim}1$ Mpc) while potentially resolving some of its small-scale problems via kiloparsec-scale quantum interference. However, the most basic FDM model, with one free parameter (the boson mass), is subject to a tension: small boson masses yield the desired cores of dwarf galaxies but underpredict structure in the Lyman-α forest, while large boson masses render FDM effectively identical to CDM. This Catch-22 problem may be alleviated by considering an axion-like particle with attractive particle self-interactions. We simulate an idealized FDM halo with self-interactions parametrized by an energy decay constant $f \sim 10^{15}~\rm {GeV}$ related to the axion symmetry-breaking conjectured to solve the strong-CP problem in particle physics. We observe solitons, a hallmark of FDM, condensing within a broader halo envelope, and find that the density profile and soliton mass depend on self-interaction strength. We propose generalized formulae to extend those from previous works to include self-interactions. We also investigate a critical mass threshold predicted for strong interactions at which the soliton collapses into a compact, unresolved state. We find that the collapse happens quickly, and its effects are initially contained to the central region of the halo.
The CALIPSO aerosol algorithms currently assign one lidar ratio (LR) value globally for each of the seven tropospheric aerosol types. In this study, MODIS total column aerosol optical depths (AODs) are used to constrain collocated CALIOP backscatter profiles in a Fernald inversion that infers aerosol LRs for CALIOP-classified marine and dusty marine aerosols. The GOCART aerosol model is leveraged to estimate the sea salt volume fractions (SSVFs) that are collocated with the CALIOP+MODIS LR retrievals. An inverse empirical relationship is found between the SSVFs and LRs (i.e., smaller SSVFs and larger LRs near coastlines, but the opposite in the remote oceans). This SSVF/LR relationship is applied to create regional and seasonal hybrid (i.e., retrieval & model-assisted) climatological LR maps so as to develop more robust LR selections for marine and dusty marine aerosols in the CALIPSO algorithms. These analyses also provide critical LR information for the next generation of spaceborne elastic backscatter lidars.
The CALIPSO aerosol algorithms currently assign one lidar ratio (LR) value globally for each of the seven tropospheric aerosol types. In this study, MODIS total column aerosol optical depths (AODs) are used to constrain collocated CALIOP backscatter profiles in a Fernald inversion that infers aerosol LRs for CALIOP-classified marine and dusty marine aerosols. The GOCART aerosol model is leveraged to estimate the sea salt volume fractions (SSVFs) that are collocated with the CALIOP+MODIS LR retrievals. An inverse empirical relationship is found between the SSVFs and LRs (i.e., smaller SSVFs and larger LRs near coastlines, but the opposite in the remote oceans). This SSVF/LR relationship is applied to create regional and seasonal hybrid (i.e., retrieval & model-assisted) climatological LR maps so as to develop more robust LR selections for marine and dusty marine aerosols in the CALIPSO algorithms. These analyses also provide critical LR information for the next generation of spaceborne elastic backscatter lidars.
The U.S. Department of Energy’s (DOE’s) Vehicle Technologies Office (VTO) and Hydrogen and Fuel Cell Technologies Office (HFTO) aim to develop sustainable, affordable and efficient technologies for transportation of goods and people. Translating investments in advanced transportation component technologies and powertrains to estimate vehicle-level fuel savings potential is critical for understanding DOE’s impact. In this work, we simulated technologies funded by VTO and HFTO for light duty vehicles. The simulations were performed across: Multiple powertrain configurations (i.e., conventional, power-split, extended-range electric vehicle, battery electric drive, and fuel-cell vehicles), Vehicle classes (i.e., compact car, midsize car, small sport utility vehicle [SUV], midsize SUV, and pickup trucks); and Fuels (i.e., gasoline, diesel, hydrogen, and battery electricity). These various technologies are assessed for six different timeframes: laboratory years 2015, 2020, 2025, 2030, and 2045. A delay of 5 years is assumed between laboratory year and model year (year technology is introduced into production). Finally, uncertainties are included for both technology performance and cost aspects by considering two cases: Low case, aligned with DOE technology manager estimates of expected original equipment manufacturer (OEM) improvements based on regulations, business as usual; and High case, aligned with aggressive technology advancements based on R&D targets developed through support by VTO & HFTO. These scenarios are not intended as predictions of future performances. The energy and cost impact of different technologies were estimated using Autonomie (www.autonomie.net), Argonne vehicle system simulation tool. Autonomie is a state-of-the-art vehicle system simulation tool used to assess the energy consumption, performance and cost of multiple advanced vehicle technologies across classes (from light to heavy duty), powertrains (from conventional to HEVs, FCEVs, PHEVs and BEVs), components and control strategies. Autonomie is packaged with a complete set of vehicle models for a wide range of vehicle classes, powertrain configurations and component technologies, including vehicle level and component level controls. These controls were developed and calibrated using dynamometer test data. Autonomie has been used to support a wide range of studies including analyzing various component technologies, sizing powertrains components for different vehicle requirements, comparing the benefits of powertrain configurations, optimizing both heuristic and route based vehicle energy control and predicting transportation energy use when paired with a traffic modeling tool such as POLARIS. This report documents the assumptions and estimates the vehicle-level energy consumption benefits and associated technology costs for the various types of light duty vehicles. All details of vehicle assumptions and simulation results are available in the spreadsheets accompanying this report.
Tropical cyclone Paka formed during the first week of December 1997 and underwent three periods of rapid intensification over the following two weeks. During one of these periods, which initiated early on December 10, Paka's Dvorak-measured windspeed increased from 23 to 60 m/s over a 48-hr period. On December 18, during the last rapid deepening episode, Paka became a supertyphoon with a maximum wind speed of about 80 m/s. In this study, the Penn State/NCAR Mesoscale Model (MM5) with improved physics (i.e., cloud microphysics, radiation, land-soil-vegetation-surface processes, and TOGA COARE flux scheme) and a multiple level nesting technique (135, 45 and 15 km horizontal resolution) will be used to simulate supertyphoon Paka. We performed two runs initialized with Goddard Earth Observing System (GEOS) data sets. The first GEOS data set does not incorporate either TRMM (tropical rainfall measuring mission satellite) or SSM/I (sensor microwave imager) observed rainfall fields into the GEOS's assimilation system while the second one does. Preliminary results show that the MM5 simulated surface pressure deepened by more than 25 mb (45 km resolution domain) in the run initialized with the GEOS data set incorporating TRMM and SSM/I derived rainfall, compared to the one initialized without. However, the track and precipitation patterns are quite similar between the runs. In our presentation, we will show the impact of TRMM rainfall upon the MM5 simulation of Paka at various horizontal resolutions. We will also examine the physical processes associated with initial explosive development by comparing MM5 simulated rainfall and latent heat release. In addition, budget (vorticity, PV, momentum and heat) calculations and sensitivity tests will be performed to examine the upper-tropospheric and SST mechanisms responsible for the explosive development of Paka.
Increasing the energy storage capability of ionic liquid supercapacitors will require better understanding of ion-electrode interactions. We have probed the influence of these interactions on the structure and differential capacitance of of an ionic liquid ([EMIM][BF4]) at an ideal graphite interface as a function of model accuracy. Of note, differential capacitance is determined through newly derived and validated fluctuation formulas. In terms of model accuracy, we test electrostatic techniques, electrode charging techniques, and electrolyte interatomic potentials. For electrostatic summations, we employ high cost, high fidelity techniques as well as less expensive, approximate techniques for summation in slab geometry. For electrode charging, uniform, constant-charge and environmentally responsive, constant-potential conditions are employed. For the ionic liquid, constant charge and atomically polarizable models are employed. We comment on the role of model accuracy on the structure and energetics of the electric double layer as well as on the magnitude and shape of differential capacitance.
The electrical coupling between the solar wind, magnetosphere, and ionosphere is studied. The coupling is analyzed using observations of high-latitude ion convection measured by the Sondre Stromfjord radar in Greenland and a computer simulation. The computer simulation calculates the ionospheric electric potential distribution for a given configuration of field-aligned currents and conductivity distribution. The technique for measuring F-region in velocities at high time resolution over a large range of latitudes is described. Variations in the currents on ionospheric plasma convection are examined using a model of field-aligned currents linking the solar wind with the dayside, high-latitude ionosphere. The data reveal that high-latitude ionospheric convection patterns, electric fields, and field-aligned currents are dependent on IMF orientation; it is observed that the electric field, which drives the F-region plasma curve, responds within about 14 minutes to IMF variations in the magnetopause. Comparisons of the simulated plasma convection with the ion velocity measurements reveal good correlation between the data.
A computer model of a LiBr-H2O single-stage absorption machine was developed. The model, utilizing a given set of design data such as water-flow rates and inlet or outlet temperatures of these flow rates but without knowing the interior characteristics of the machine (heat transfer rates and surface areas), can be used to predict or simulate off-design performance. Results from 130 off-design cases for a given commercial machine agree with the published data within 2 percent.
Sulfur dioxide-derived cloud condensation nuclei are expected to enhance the planetary albedo, thereby cooling the planet. This effect might counteract the global warming expected from enhanced greenhouse gases. A detailed treatment of the relationship between fossil fuel burning and the SO2 effect on cloud albedo is implemented in a two-dimensional model for assessing the climate impact. Using a conservative approach, results show that the cooling induced by the SO2 emission can presently counteract 50 percent of the CO2 greenhouse warming. Since 1980, a strong warming trend has been predicted by the model: 0.15 C during the 1980-1990 period alone. The model predicts that by the year 2060 the SO2 cooling reduces climate warming by 0.5 C or 25 percent for the Intergovernmental Panel on Climate Change (IPCC) business as usual (BAU) scenario and 0.2 C or 20 percent for scenario D (for a slow pace of fossil fuel burning). The hypothesis is examined that the different responses between the Northern Hemisphere and the Southern Hemisphere can be used to validate the presence of the SO2-induced cooling.
Cloud microphysics are inevitable affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds. Therefore, size distribution parameterized as spectral bin microphysics are needed to explicitly study the effect of atmospheric aerosol concentration on cloud development, rainfall production, and rainfall rates convective clouds. Recently, two detailed spectral-bin microphysical schemes were implemented into the Goddard Cumulus Ensembel (GCE) model. The formulation for the explicit spectral-bim microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles [i.e., pristine ice crystals (columnar and plate-like), snow (dendrites and aggregates), groupel and frozen drops/hall] Each type is described by a special size distribution function containing many categories (i.e., 33 bins). Atmospheric aerosols are also described using number density size-distribution functions.A spectral-bin microphysical model is very expensive from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep cloud systems in the west Pacific warm pool region and in the mid-latitude using identical thermodynamic conditions but with different concentrations of CCN: a low "clean" concentration and a high "dirty" concentration. Besides the initial differences in aerosol concentration, preliminary results indicate that the low CCN concentration case produces rainfall at the surface sooner than the high CCN case but has less cloud water mass aloft. Because the spectral-bim model explicitly calculates and allows for the examination of both the mass and number concentration of cpecies in each size category, a detailed analysis of the instantaneous size spectrum can be obtained for the two cases. It is shown that since the low CCN case produces fever droplets, larger size develop due to greater condencational and collectional growth, leading to a broader size spectrum in comparison to the high CCN case.
Cloud microphysics are inevitably affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds. Therefore, size distributions parameterized as spectral bin microphysics are needed to explicitly study the effects of atmospheric aerosol concentration on cloud development, rainfall production, and rainfall rates for convective clouds. Recently, two detailed spectral-bin microphysical schemes were implemented into the Goddard Cumulus Ensemble (GCE) model. The formulation for the explicit spectral-bin microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles [i.e.,pristine ice crystals (columnar and plate-like), snow (dendrites and aggregates), graupel and frozen drops/hail]. Each type is described by a special size distribution function containing many categories (i.e. 33 bins). Atmospheric aerosols are also described using number density size-distribution functions.A spectral-bin microphysical model is very expensive from a from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region using identical thermodynamic conditions but with different concentrations of CCN: a low "clean" concentration and a high "dirty" concentration. Besides the initial differences in aerosol concentration, preliminary results indicate that the low CCN concentration case produces rainfall at the surface sooner than the high CCN case but has less cloud water mass aloft. Because the spectral-bin model explicitly calculates and allows for the examination of both the mass and number concentration of species in each size categor, a detailed analysis of the instantaneous size spectrum can be obtained for the two cases. It is shown that since the low CCN case produces fewer droplets, larger sized develop due to the greater condensational and collectional growth, leading to a broader size spectrum in comparison to the high CCN case.
Cloud microphysics is inevitably affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds. Therefore, size distributions parameterized as spectral bin microphysics are needed to explicitly study the effects of atmospheric aerosol concentration on cloud development, rainfall production, and rainfall rates for convective clouds. Recently, two detailed spectral-bin microphysical schemes were implemented into the Goddard Cumulus Ensembel (GCE) model. The formulation for the explicit spectral-bin microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles [i.e. pristine ice crystals (columnar and plate-like), snow (dendrites and aggregates), graupel and frozen drops/hail]. Each type is described by a special size distribution function containing many categories (i.e. 33 bins). Atmospheric aerosols are also described using number density size distribution functions. A spectral-bin microphysical model is very expensive from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region and in the mid-latitude continent with different concentrations of CCN: a low "c1ean"concentration and a high "dirty" concentration. In addition, differences and similarities between bulk microphysics and spectral-bin microphysical schemes will be examined and discussed.