SEARCH · Engineering Papers

Results for “Statistical simulation”

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 1,207 records · Page 67

All sources NASA NTRS DOE OSTI

NASA NTRS1984

The effect of rotation on turbulence structure

Rotating turbulent channel flow has been simulated numerically using the large-eddy-simulation technique. With a three-dimensional, time-dependent, numerical simulation of turbulence, it was possible to obtain detailed flow-field information which had not been available previously from experiments. This information was used to study the structure and statistical properties of the flow with externally imposed rotation. It was found that in agreement with experimental observations, the Coriolis force introduced by the system rotation could stabilize/destabilize turbulence. The large-scale motion observed in the experiment was also found to exist in the computed flow field.

Kim, J.↗

NASA NTRS1988

Global environmental effects of impact-generated aerosols: Results from a general circulation model

Cooling and darkening at Earth's surface are expected to result from the interception of sunlight by the high altitude worldwide dust cloud generated by impact of a large asteroid or comet, according to the one-dimensional radioactive-convective atmospheric model (RCM) of Pollack et al. An analogous three-dimensional general circulation model (GCM) simulation obtains the same basic result as the RCM but there are important differences in detail. In the GCM simulation the heat capacity of the oceans, not included in the RCM, substantially mitigates land surface cooling. On the other hand, the GCM's low heat capacity surface allows surface temperatures to drop much more rapidly than reported by Pollack et al. These two differences between RCM and GCM simulations were noted previously in studies of nuclear winter; GCM results for comet/asteroid winter, however, are much more severe than for nuclear winter because the assumed aerosol amount is large enough to intercept all sunlight falling on Earth. In the simulation the global average of land surface temperature drops to the freezing point in just 4.5 days, one-tenth the time required in the Pollack et al. simulation. In addition to the standard case of Pollack et al., which represents the collision of a 10-km diameter asteroid with Earth, additional scenarios are considered ranging from the statistically more frequent impacts of smaller asteroids to the collision of Halley's comet with Earth. In the latter case the kinetic energy of impact is extremely large due to the head-on collision resulting from Halley's retrograde orbit.

Covey, C.↗

NASA NTRS1988

Analysis of homogeneous turbulent reacting flows

Full turbulence simulations at low Reynolds numbers were made for the single-step, irreversible, bimolecular reaction between non-premixed reactants in isochoric, decaying homogeneous turbulence. Various initial conditions for the scalar field were used in the simulations to control the initial scalar dissipation length scale, and simulations were also made for temperature-dependent reaction rates and for non-stoichiometric and unequal diffusivity conditions. Joint probability density functions (pdf's), conditional pdf's, and various statistical quantities appearing in the moment equations were computed. Preliminary analysis of the results indicates that compressive strain-rate correlates better than other dynamical quantities with local reaction rate, and the locations of peak reaction rates seem to be insensitive to the scalar field initial conditions.

Leonard, A. D.↗

NASA NTRS1989

Evolution of cosmic string networks

Results on cosmic strings are summarized including: (1) the application of non-equilibrium statistical mechanics to cosmic string evolution; (2) a simple one scale model for the long strings which has a great deal of predictive power; (3) results from large scale numerical simulations; and (4) a discussion of the observational consequences of our results. An upper bound on G mu of approximately 10(-7) emerges from the millisecond pulsar gravity wave bound. How numerical uncertainties affect this are discussed. Any changes which weaken the bound would probably also give the long strings the dominant role in producing observational consequences.

Albrecht, Andreas↗

NASA NTRS2011

Results of the Simulation and Assimilation of Doppler Wind Lidar Observations in Preparation for European Space Agency's Aeolus Mission

With the launch of the European Space Agency's Aeolus Mission in 2013, direct spaceborne measurements of vertical wind profiles are imminent via Doppler wind lidar technology. Part of the preparedness for such missions is the development of the proper data assimilation methodology for handling such observations. Since no heritage measurements exist in space, the Joint Observing System Simulation Experiment (Joint OSSE) framework has been utilized to generate a realistic proxy dataset as a precursor to flight. These data are being used for the development of the Gridpoint Statistical Interpolation (GSI) data assimilation system utilized at a number of centers through the United States including the Global Modeling and Assimilation Office (GMAO) at NASA/Goddard Space Flight Center and at the National Centers for Environmental Prediction (NOAA/NWS/NCEP) as an activity through the Joint Center for Satellite Data Assimilation. An update of this ongoing effort will be presented, including the methodology of proxy data generation, the limitations of the proxy data, the handling of line-of-sight wind measurements within the GSI, and the impact on both analyses and forecasts with the addition of the new data type.

McCarty, Will↗

NASA NTRS2011

The Simulation and Assimilation of Doppler Wind Lidar Observations in Support of Future Instruments

With the launch of the European Space Agency's Atmospheric Dynamics Mission (ADM-Aeolus) in 2011 and the call for the 3D-Winds mission in National Research Council's decadal survey, direct spaceborne measurements of vertical wind profiles are imminent via Doppler wind lidar technology. Part of the preparedness for such missions is the development of the proper data assimilation methodology for handling such observations. Since no heritage measurements exist in space, the Joint Observing System Simulation Experiment (Joint OSSE) framework is being utilized to generate a realistic proxy dataset as a precursor to flight. These data are being used for the development of the Gridpoint Statistical Interpolation (GSI) data assimilation system utilized at a number of centers through the United States including the Global Modeling and Assimilation Office (GMAO) at NASA/Goddard Space Flight Center and at the National Centers for Environmental Prediction (NOAA/NWS/NCEP). This effort will be presented, including the methodology of proxy data generation, the handling of line-of-sight wind measurements within the GSI, and the impact on both analyses and forecasts with the addition of the new data type.

McCarty, Will↗

NASA NTRS2011

Climate Change Impact Uncertainties for Maize in Panama: Farm Information, Climate Projections, and Yield Sensitivities

We present results from a pilot project to characterize and bound multi-disciplinary uncertainties around the assessment of maize (Zea mays) production impacts using the CERES-Maize crop model in a climate-sensitive region with a variety of farming systems (Panama). Segunda coa (autumn) maize yield in Panama currently suffers occasionally from high water stress at the end of the growing season, however under future climate conditions warmer temperatures accelerate crop maturation and elevated CO (sub 2) concentrations improve water retention. This combination reduces end-of-season water stresses and eventually leads to small mean yield gains according to median projections, although accelerated maturation reduces yields in seasons with low water stresses. Calibrations of cultivar traits, soil profile, and fertilizer amounts are most important for representing baseline yields, however sensitivity to all management factors is reduced in an assessment of future yield changes (most dramatically for fertilizers), suggesting that yield changes may be more generalizable than absolute yields. Uncertainty around General Circulation Model (GCM)s' projected changes in rainfall gain in importance throughout the century, with yield changes strongly correlated with growing season rainfall totals. Climate changes are expected to be obscured by the large inter-annual variations in Panamanian climate that will continue to be the dominant influence on seasonal maize yield into the coming decades. The relatively high (A2) and low (B1) emissions scenarios show little difference in their impact on future maize yields until the end of the century. Uncertainties related to the sensitivity of CERES-Maize to carbon dioxide concentrations have a substantial influence on projected changes, and remain a significant obstacle to climate change impacts assessment. Finally, an investigation into the potential of simple statistical yield emulators based upon key climate variables characterizes the important uncertainties behind the selection of climate change metrics and their performance against more complex process-based crop model simulations, revealing a danger in relying only on long-term mean quantities for crop impact assessment.

Maize↗

NASA NTRSDate not supplied

TPSAS-NF1676L-33626-DND

NASA is conducting a series of computational experiments to quantify atmospheric effects on low noise sonic booms. In the current study, simulated cruise nearfield pressure data from NASA?s X-59 Quiet Supersonic Technology aircraft was propagated from the aircraft to the ground at four cardinal headings through five years of realistic atmospheric profiles at 30 locations across the eastern USA. Statistical design of experiments was used to select the locations where primary sonic thump carpet widths (CW) and metric levels at the ground were computed. Atmospheric profiles were taken from the Climate Forecast System Reanalysis database, which contains reanalyzed atmospheric profiles four times daily. Decision tree analyses were performed to determine relative importance of predictors for CW and metric levels. Predictors included were latitude, longitude, date, time of day, season, climate, aircraft heading, and ground elevation. Results of this study indicate the propagation resolution needed to adequately characterize the distributions of CW and ground metric data, i.e., necessary separation distance between propagation locations and total number of atmospheres through which to propagate. This resolution will be used for a follow-on study of simulated X-59 carpets across the entire US mainland.

Will Doebler↗

NASA NTRS2007

A Simple Stochastic Model for Generating Broken Cloud Optical Depth and Top Height Fields

A simple and fast algorithm for generating two correlated stochastic twodimensional (2D) cloud fields is described. The algorithm is illustrated with two broken cumulus cloud fields: cloud optical depth and cloud top height retrieved from Moderate Resolution Imaging Spectrometer (MODIS). Only two 2D fields are required as an input. The algorithm output is statistical realizations of these two fields with approximately the same correlation and joint distribution functions as the original ones. The major assumption of the algorithm is statistical isotropy of the fields. In contrast to fractals and the Fourier filtering methods frequently used for stochastic cloud modeling, the proposed method is based on spectral models of homogeneous random fields. For keeping the same probability density function as the (first) original field, the method of inverse distribution function is used. When the spatial distribution of the first field has been generated, a realization of the correlated second field is simulated using a conditional distribution matrix. This paper is served as a theoretical justification to the publicly available software that has been recently released by the authors and can be freely downloaded from http://i3rc.gsfc.nasa.gov/Public codes clouds.htm. Though 2D rather than full 3D, stochastic realizations of two correlated cloud fields that mimic statistics of given fields have proved to be very useful to study 3D radiative transfer features of broken cumulus clouds for better understanding of shortwave radiation and interpretation of the remote sensing retrievals.

Prigarin, Sergei M.↗

NASA NTRS1982

The effects of horizontal body casting on blood volume, drug responsiveness, and +Gz tolerance in the rhesus monkey

To simulate the weightless condition, eight rhesus monkeys, instrumented with solid-state pressure transducers, were horizontally restrained in body casts for 28 days. Blood volume decreased an average of 13% after 14 days of restraint, due mainly to a drop in plasma volume. Aortic pressure and heart rate responses to norepinephrine and phenylephrine decreased after 14 days of restraint. The monkeys did not show a statistically significant decreased tolerance to a 90 deg sudden upright tilt after horizontal restraint. During the fifth week of casting, four animals were subjected to +Gz acceleration tests on a centrifuge. The acceleration tolerance of the casted monkeys was significantly reduced compared to four similarly instrumented control animals. These findings indicate that the cardiovascular deconditioning associated with simulated weightlessness results from an inability to maintain central blood volume during orthostatic stress.

Dickey, D. T.↗

NASA NTRS2015

Free-Energy Calculations. A Mathematical Perspective

Ion channels are pore-forming assemblies of transmembrane proteins that mediate and regulate ion transport through cell walls. They are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses. They are also essential to cardiac processes, muscle contraction and epithelial transport. Ion channels from lower organisms can act as toxins or antimicrobial agents, and in a number of cases are involved in infectious diseases. Because of their important and diverse biological functions they are frequent targets of drug action. Also, simple natural or synthetic channels find numerous applications in biotechnology. For these reasons, studies of ion channels are at the forefront of biophysics, structural biology and cellular biology. In the last decade, the increased availability of X-ray structures has greatly advanced our understanding of ion channels. However, their mechanism of action remains elusive. This is because, in order to assist controlled ion transport, ion channels are dynamic by nature, but X-ray crystallography captures the channel in a single, sometimes non-native state. To explain how ion channels work, X-ray structures have to be supplemented with dynamic information. In principle, molecular dynamics (MD) simulations can aid in providing this information, as this is precisely what MD has been designed to do. However, MD simulations suffer from their own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. To assess the reliability of MD simulations it is only natural to turn to the main function of channels - conducting ions - and compare calculated ionic conductance with electrophysiological data, mainly single channel recordings, obtained under similar conditions. If this comparison is satisfactory it would greatly increase our confidence that both the structures and our computational methodologies are sufficiently accurate. Channel conductance, defined as the ratio of ionic current through the channel to applied voltage, can be calculated in MD simulations by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. If the current is small, a voltage significantly higher than the experimental one needs to be applied to collect sufficient statistics of ion crossing events. Then, the calculated conductance has to be extrapolated to the experimental voltage using procedures of unknown accuracy. Instead, we propose an alternative approach that applies if ion transport through channels can be described with sufficient accuracy by the one-dimensional diffusion equation in the potential given by the free energy profile and applied voltage. Then, it is possible to test the assumptions of the equation, recover the full voltage/current dependence, determine the reliability of the calculated conductance and reconstruct the underlying (equilibrium) free energy profile, all from MD simulations at a single voltage. We will present the underlying theory, model calculations that test this theory and simulations on ion conductance through a channel that has been extensively studied experimentally. To our knowledge this is the first case in which the complete, experimentally measured dependence of the current on applied voltage has been reconstructed from MD simulations.

free energy↗

NASA NTRS1991

A statistical approach to the brittle fracture of a multi-phase solid

A stochastic damage model is proposed to quantify the inherent statistical distribution of the fracture toughness of a brittle, multi-phase solid. The model, based on the macrocrack-microcrack interaction, incorporates uncertainties in locations and orientations of microcracks. Due to the high concentration of microcracks near the macro-tip, a higher order analysis based on traction boundary integral equations is formulated first for an arbitrary array of cracks. The effects of uncertainties in locations and orientations of microcracks at a macro-tip are analyzed quantitatively by using the boundary integral equations method in conjunction with the computer simulation of the random microcrack array. The short range interactions resulting from surrounding microcracks closet to the main crack tip are investigated. The effects of microcrack density parameter are also explored in the present study. The validity of the present model is demonstrated by comparing its statistical output with the Neville distribution function, which gives correct fits to sets of experimental data from multi-phase solids.

Liu, W. K.↗

NASA NTRS2000

Operational Assimilation of GOES Data into a Mesoscale Model

A technique has been developed for assimilating GOES-derived skin temperature tendencies and insolation into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite- observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. The technique has been employed on a semi-operational basis at the Global Hydrology and Climate Center (GHCC) within the Penn State/National Center for Atmospheric Research (PSU/NCAR) Mesoscale Model (MM5) since 1 November 1998. We performed the assimilation on a model grid centered over the Southeastern US. In addition, a control run without assimilation was performed to provide insight into the performance of the assimilation technique. Bulk verification statistics (BIAS and RMSE) of surface air temperature and relative humidity of more than 250 case days has been performed to date. Results show that assimilation of the satellite data results reduces both the bias and RMSE for simulations of surface air temperature and relative humidity. We are working with forecasters at the National Weather Service Forecast Office located in Birmingham, AL to evaluate the impact of the assimilation on precipitation forecasts. In addition, work is currently underway to determine the ability of both the assimilation technique and Early Eta to sense the significant changes that occur in both leaf expansion and soil moisture conditions during the Spring-Summer of 1999. This will be accomplished via examination of the simulated surface energy budgets from both modeling systems.

Lapenta, William↗

NASA NTRS1994

Ozone formation during an episode over Europe: A 3-D chemical/transport model simulation

A 3-D regional photochemical tracer/transport model for Europe and the Eastern Atlantic has been developed based on the NASA/GISS CTM. The model resolution is 4x5 degrees latitude and longitude with 9 layers in the vertical (7 in the troposphere). Advective winds, convection statistics and other meteorological data from the NASA/GISS GCM are used. An extensive gas-phase chemical scheme based on the scheme used in our global 2D model has been incorporated in the 3D model. In this work ozone formation in the troposphere is studied with the 3D model during a 5 day period starting June 30. Extensive local ozone production is found and the relationship between the source regions and the downwind areas are discussed. Variations in local ozone formation as a function of total emission rate, as well as the composition of the emissions (HC/NO(x)) ratio and isoprene emissions) are elucidated. An important vertical transport process in the troposphere is by convective clouds. The 3D model includes an explicit parameterization of this process. It is shown that this process has significant influence on the calculated surface ozone concentrations.

Berntsen, Terje↗

NASA NTRS2019

Assimilation of Satellite-Based Snow Cover and Freeze/Thaw Observations Over the High Mountain Asia

Toward qualifying hydrologic changes in the High Mountain Asia (HMA) region, this study explores the use of a hyper-resolution (1 km) land data assimilation (DA)framework developed within the NASA Land Information System using the Noah Multi-parameterization Land Surface Model (Noah-MP) forced by the meteorological boundary conditions from Modern-Era Retrospective analysis for Research and Applications, Version 2 data. Two different sets of DA experiments are conducted:(1) the assimilation of a satellite-derived snow cover map (MOD10A1) and (2) the assimilation of the NASA MEaSUREs landscape freeze/thaw product from 2007 to 2008. The performance of the snow cover assimilation is evaluated via comparisons with available remote sensing-based snow water equivalent product and ground-based snow depth measurements. For example, in the comparison against ground-based snow depth measurements, the majority of the stations (13 of 14) show slightly improved goodness-of-fit statistics as a result of the snow DA, but only four are statistically significant. In addition, comparisons to the satellite-based land surface temperature products (MOD11A1 and MYD11A1) show that freeze/thaw DA yields improvements (at certain grid cells) of up to 0.58 K in the root-mean-square error (RMSE) and 0.77K in the absolute bias (relative to model-only simulations). In the comparison against three ground-based soil temperature measurements along the Himalayas, the bias and the RMSE in the 0-10 cm soil temperature are reduced (on average) by 10 and 7%,respectively. The improvements in the top layer of soil estimates also propagate through the deeper soil layers, where the bias and the RMSE in the 10-40 cm soil temperature are reduced (on average) by 9 and 6%, respectively. However, no statistically significant skill differences are observed for the freeze/thaw DA system in the comparisons against ground-based surface temperature measurements at mid-to-low altitude. Therefore, the two proposed DA schemes show the potential of improving the predictability of snow mass, surface temperature, and soil temperature states across HMA, but more ground-based measurements are still required, especially at high-altitudes, in order to document a more statistically significant improvement as a result of the two DA schemes.

High Mountain Asia↗

NASA NTRS1979

Advanced Meteorological Temperature Sounder (AMTS) simulations

Simulation studies are reported on temperature retrievals from AMTS and their effect on atmospheric analysis. Observations are simulated from radiosonde reports and observed cloud cover. Temperature retrievals are performed and RMS temperature and thickness errors are calculated relative to the radiosonde profiles and compared to similarly generated HIRS statistics. Significant improvement over HIRS is found throughout the atmosphere but especially in the stratosphere and lower troposphere.

Susskin, J.↗

NASA NTRS1982

Radiation effects studies for the high-resolution spectrograph

The generation and collection of charge carriers created during the passage of energetic protons through a silicon photodiode array are modeled. Pulse height distributions of noise charge collected during exposure of a digicon type diode array to 21 and 75 MeV protons were obtained. The magnitude of charge collected by a diode from each proton event is determined not only by diffusion, but by statistical considerations involving the ionization process itself. Utilizing analytical solutions to the diffusion equation for transport of minority carriers, together with the Vavilov theory of energy loss fluctuations in thin absorbers, simulations of the pulse height spectra which follow the experimental distributions fairly well are presented and an estimate for the minority carrier diffusion length L sub d is provided.

Smith, L. C.↗

NASA NTRS1984

Direct numerical simulation of curved turbulent channel flow

Low Reynolds number, mildly curved, turbulent channel flow has been simulated numerically without subgrid scale models. A new spectral numerical method developed for this problem was used, and the computations were performed with 2 million degrees of freedom. A variety of statistical and structural information has been extracted from the computed flow fields. These include mean velocity, turbulence stresses, velocity skewness, and flatness factors, space time correlations and spectra, all the terms in the Reynolds stress balance equations, and contour and vector plots of instantaneous velocity fields. The effects of curvature on this flow were determined by comparing the concave and convex sides of the channel. The observed effects are consistent with experimental observations for mild curvature. The most significant difference in the turbulence statistics between the concave and convex sides was in the Reynolds shear stress. This was accompanied by significant differences in the terms of the Reynolds shear stress balance equations. In addition, it was found that stationary Taylor-Gortler vortices were present and that they had a significant effect on the flow by contributing to the mean Reynolds shear stress, and by affecting the underlying turbulence.

Moser, R. D.↗