Search NASA⌕ Search

SEARCH · Search NASA

Results for “Modeling & Simulation (M&S)”

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 109 records · Page 6

Modeling of Highly Instrumented Honeywell Turbofan Engine Tested with Ice Crystal Ingestion in the NASA Propulsion System Laboratory

The Propulsion Systems Laboratory (PSL), an altitude test facility at NASA Glenn Research Center, has been used to test a highly instrumented turbine engine at simulated altitude operating conditions. This is a continuation of the PSL testing that successfully duplicated the icing events that were experienced in a previous engine (serial LF01) during flight through ice crystal clouds, which was the first turbofan engine tested in PSL. This second model of the ALF502R-5A serial number LF11 is a highly instrumented version of the previous engine. The PSL facility provides a continuous cloud of ice crystals with controlled characteristics of size and concentration, which are ingested by the engine during operation at simulated altitudes. Several of the previous operating points tested in the LF01 engine were duplicated to confirm repeatability in LF11. The instrumentation included video cameras to visually illustrate the accretion of ice in the low pressure compressor (LPC) exit guide vane region in order to confirm the ice accretion, which was suspected during the testing of the LF01. Traditional instrumentation included static pressure taps in the low pressure compressor inner and outer flow path walls, as well as total pressure and temperature rakes in the low pressure compressor region. The test data was utilized to determine the losses and blockages due to accretion in the exit guide vane region of the LPC. Multiple data points were analyzed with the Honeywell Customer Deck. A full engine roll back point was modeled with the Numerical Propulsion System Simulation (NPSS) code. The mean line compressor flow analysis code with ice crystal modeling was utilized to estimate the parameters that indicate the risk of accretion, as well as to estimate the degree of blockage and losses caused by accretion during a full engine roll back point. The analysis provided additional validation of the icing risk parameters within the LPC, as well as the creation of models for estimating the rates of blockage growth and losses.

Simulation↗

Microburst vertical wind estimation from horizontal wind measurements

The vertical wind or downdraft component of a microburst-generated wind shear can significantly degrade airplane performance. Doppler radar and lidar are two sensor technologies being tested to provide flight crews with early warning of the presence of hazardous wind shear. An inherent limitation of Doppler-based sensors is the inability to measure velocities perpendicular to the line of sight, which results in an underestimate of the total wind shear hazard. One solution to the line-of-sight limitation is to use a vertical wind model to estimate the vertical component from the horizontal wind measurement. The objective of this study was to assess the ability of simple vertical wind models to improve the hazard prediction capability of an airborne Doppler sensor in a realistic microburst environment. Both simulation and flight test measurements were used to test the vertical wind models. The results indicate that in the altitude region of interest (at or below 300 m), the simple vertical wind models improved the hazard estimate. The radar simulation study showed that the magnitude of the performance improvement was altitude dependent. The altitude of maximum performance improvement occurred at about 300 m.

Vicroy, Dan D.↗

External Sources of Water for Mercury's Putative Ice Deposits

Radar images have revealed the possible presence of ice deposits in Mercury's polar regions. Although thermal models indicate that water ice can be stable in permanently shaded regions near Mercury's poles, the ultimate source of the water remains unclear. We use stochastic models and other theoretical methods to investigate the role of external sources in supplying Mercury with the requisite amount of water. By extrapolating the current terrestrial influx of interplanetary dust particles to that at Mercury, we find that continual micrometeoritic bombardment of Mercury over the last 3.5 byr could have resulted in the delivery of (3-60) x 10(exp 16) grams of water ice to the permanently shaded regions at Mercury's poles (equivalent to an average ice thickness of 0.8-20 m). Erosion by micrometeoritic impact on exposed ice deposits could reduce the above value by about a half. For comparison, the current ice deposits on Mercury are believed to be somewhere between approx. 2 and 20 m thick. Using a Monte Carlo model to simulate the impact history of Mercury, we find that asteroids and comets can also deliver an amount of water consistent with the observations. Impacts from Jupiter-family comets over the last 3.5 billion years can supply (0.1-200) x 10(exp 16) grams of water to Mercury's polar regions (corresponding to ice deposits 0.05-60 m thick), Halley-type comets can supply (0.2-20) x 10(exp 16) grams of water to the poles (0.07-7 m of ice), and asteroids can provide (0.4-20) x 10(exp 16) grams of water to the poles (0.1-8 m of ice). Although all these external sources are nominally sufficient to explain the estimated amount of ice currently at Mercury's poles, impacts by a few large comets and/or asteroids seem to provide the best explanation for both the amount and cleanliness of the ice deposits on Mercury. Despite their low population estimates in the inner solar system, Jupiter-family comets are particularly promising candidates for delivering water to Mercury because they have a larger volatile content than asteroids and more favorable orbital and impact characteristics than Halley-type comets.

Moses, Julianne I.↗

Simulated interannual variability of the Greenland Sea deep water formation and its connection to surface forcing

A fully prognostic Arctic ice-ocean model is used to study the interannual variability of deepwater formation in the Greenland Sea Gyre based on the simulations for the Arctic ice-ocean system for the period 1955 and 1960 - 1985. The model uses monthly climatology for thermodynamic forcing components (such as air temperature and cloudiness), together with constant annual net precipitation and river runoff. The daily wind forcing is derived from analyzed sea level air pressures from the National Center for Atmospheric Research (NCAR). In summary, the model shows that the occurence of deep convection in the Greenland Sea Gyre is controlled by the extensive Fram Strait ice export and/or local wind conditions in the Greenland Sea. In the latter case the weakening of the local wind curl allows the Polar Front to move eastward. The movement of the Polar Front causes adverse ice conditions, often together with much larger than normal ice export from the Arctic, such as in 1968, which can block convection in the gyre. The density difference between upper and lower layers is investigated as an indication of water mass formation through convection, occurring as strong diffusion in the model. The model-simulated density difference between the average top 100 m and deep levels reveals that the period 1960 - 1985 had only a few distinct years with weak stratification, and, especially, the model predicts no deep convection since the nid-1970s. The common factor for the years of the weakest decrease of the model-predicted heat content of the upper 2000 m which can, to a high degree, be explained by local heat loss.

Haekkinen, Sirpa↗

Effect of high-latitude ionospheric convection on Sun-aligned polar caps

A coupled magnetospheric-ionospheric (M-I) magnetohydrodynamic (MHD) model has been used to simulate the formation of Sun-aligned polar cap arcs for a variety of interplanetary magnetic field (IMF) dependent polar cap convection fields. The formation process involves launching an Alfven shear wave from the magnetosphere to the ionosphere where the ionospheric conductance can react self-consistently to changes in the upward currents. We assume that the initial Alfven shear wave is the result of solar wind-magnetosphere interactions. The simulations show how the E region density is affected by the changes in the electron precipitation that are associated with the upward currents. These changes in conductance lead to both a modified Alfven wave reflection at the ionosphere and the generation of secondary Alfven waves in the ionosphere. The ensuing bouncing of the Alfven waves between the ionosphere and magnetosphere is followed until an asymptotic solution is obtained. At the magnetosphere the Alfven waves reflect at a fixed boundary. The coupled M-I Sun-aligned polar cap arc model of Zhu et al.(1993a) is used to carry out the simulations. This study focuses on the dependence of the polar cap arc formation on the background (global) convection pattern. Since the polar cap arcs occur for northward and strong B(sub y) IMF conditions, a variety of background convection patterns can exist when the arcs are present. The study shows that polar cap arcs can be formed for all these convection patterns; however, the arc features are dramatically different for the different patterns. For weak sunward convection a relatively confined single pair of current sheets is associated with the imposed Alfven shear wave structure. However, when the electric field exceeds a threshold, the arc structure intensifies, and the conductance increases as does the local Joule heating rate. These increases are faster than a linear dependence on the background electric field strength. Furthermore, above the threshold, the single current sheet pair splits into multiple current sheet pairs. For the fixed initial ionospheric and magnetospheric conditions used in this study, the separation distance between the current pairs was found to be almost independent of the background electric field strength. For either three-cell or distorted two-cell background convection patterns the arc formation favored the positive B(sub y) case in the northern hemisphere.

Sojka, J. J.↗

Improved Simulations of Biomass Burning Aerosol Optical Properties and Lifetimes in the NASA GEOS Model during the ORACLES-I Campaign

In order to improve aerosol representation in the NASA Goddard Earth Observing System (GEOS) model, we evaluated simulations of the aerosol properties and transport over the southern African biomass burning source and outflow region using observations made during the first deployment of the ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) field campaign in September 2016. An example case study of September 24 was analyzed in detail, during which aircraft-based in-situ and remote sensing observations showed the presence of a multi-layered smoke plume structure with significant vertical variation in single scattering albedo (SSA). Our baseline GEOS simulations were not able to represent the observed SSA variation, nor the observed organic aerosol-to-black carbon ratio (OA : BC). Analyzing the simulated smoke age suggests that the higher altitude, less absorbing smoke plume was younger (~4 days), while the lower altitude and more absorbing smoke plume was older (~7 days). We hypothesize a chemical or microphysical loss process exists to explain the change in aerosol absorption as the smoke plume ages, and we apply a simple 6-day e-folding loss rate to the model hydrophilic biomass burning OA to mimic this process. Adding this loss process required some adjustment to the model assumed scaling factors of aerosol emissions to conserve the regional aerosol loading and further improve the simulated OA : BC ratio. Accordingly, we have increased our biomass burning emissions of OA by 60 % and biomass burning BC by 15 %. We also utilized the ORACLES airborne observations to better constrain the simulation of aerosol optical properties, adjusting the assumed particle size, hygroscopic growth, and absorption. Our final GEOS model simulation with additional OA loss and updated optics showed a better performance in simulating AOD and SSA compared to independent ground and space-based retrievals for the entire month of September 2016, including the OMI Aerosol Index. In terms of radiative implications of our model adjustments, the final GEOS simulation suggested a decreased atmospheric warming of about 10 % (~2 W m -2 ) over the south-east Atlantic region and above the stratocumulus cloud decks compared to the model baseline simulations. These results improve the representation of the smoke age, transport, and optical properties in Earth system models.

aerosols↗

Artificial stereo presentation of meteorological data fields

The innate capability to perceive three-dimensional stereo imagery has been exploited to present multidimensional meteorological data fields. Variations on an artificial stereo technique first discussed by Pichel et al. (1973) are used to display single and multispectral images in a vivid and easily assimilated manner. Examples of visible/infrared artificial stereo are given for Hurricane Allen and for severe thunderstorms on 10 April 1979. Three-dimensional output from a mesoscale model also is presented. The images may be viewed through the glasses inserted in the February 1981 issue of the Bulletin of the American Meteorological Society, with the red lens over the right eye. The images have been produced on the interactive Atmospheric and Oceanographic Information Processing System (AOIPS) at Goddard Space Flight Center. Stereo presentation is an important aid in understanding meteorological phenomena for operational weather forecasting, research case studies, and model simulations.

Hasler, A. F.↗

Premonsoon Aerosol Characterization and Radiative Effects Over the Indo-Gangetic Plains: Implications for Regional Climate Warming

The Himalayas have a profound effect on the South Asian climate and the regional hydrological cycle, as it forms a barrier for the strong monsoon winds and serves as an elevated heat source, thus controlling the onset and distribution of precipitation during the Indian summer monsoon. Recent studies have suggested that radiative heating by absorbing aerosols, such as dust and black carbon over the Indo-Gangetic Plains (IGP) and slopes of the Himalayas, may significantly accelerate the seasonal warming of the Hindu Kush-Himalayas-Tibetan Plateau (HKHT) and influence the subsequent evolution of the summer monsoon. This paper presents a detailed characterization of aerosols over the IGP and their radiative effects during the premonsoon season (April-May-June) when dust transport constitutes the bulk of the regional aerosol loading, using ground radiometric and spaceborne observations. During the dust-laden period, there is a strong response of surface shortwave flux to aerosol absorption indicated by the diurnally averaged forcing efficiency of -70 W/sq m per unit optical depth. The simulated aerosol single-scattering albedo, constrained by surface flux and aerosol measurements, is estimated to be 0.89+/- 0.01 (at approx.550 nm) with diurnal mean surface and top-of-atmosphere forcing values ranging from -11 to -79.8 W/sq m and +1.4 to +12 W/sq m, respectively, for the premonsoon period. The model-simulated solar heating rate profile peaks in the lower troposphere with enhanced heating penetrating into the middle troposphere (5-6 km), caused by vertically extended aerosols over the IGP with peak altitude of approx.5 km as indicated by spaceborne Cloud-Aerosol Lidar with Orthogonal Polarization observations. On a long-term climate scale, our analysis, on the basis of microwave satellite measurements of tropospheric temperatures from 1979 to 2007, indicates accelerated annual mean warming rates found over the Himalayan-Hindu Kush region (0.21 C/decade+/-0.08 C/decade) and underscores the potential role of enhanced aerosol solar absorption in the maximum warming localized over the western Himalayas (0.26 C/decade f 0.09 C/decade) that significantly exceed the entire HKHT and global warming rates. We believe the accelerated warming rates reported here are critical to both the South Asian summer monsoon and hydro-glaciological resource variability in the Himalayan-Hindu Kush snowpack and therefore to the densely populated downstream regions.

Gautam, Ritesh↗

Science-return modeling and simulation

In this paper, an on-going effort to integrate science-return modeling and simulation with JPL's institutional mission system design process will be presented.

modeling simulation↗

Flux Comparison of Master-8 and Ordem 3.1 Modelled Space Debris Population

With ESA’s Meteoroid And Space debris Terrestrial Environment Reference (MASTER-8) model and NASA’s Orbital Debris Engineering Model(ORDEM) 3.1, the two premier orbital debris engineering models have been officially released. The two models come with significant enhancements and now represent state-of-the-art orbital debris modelling for their respective agencies. Both models provide the community with estimates of the space debris environment from low Earth orbit (LEO) up to at least geostationary altitude. The MASTER population is an event-based simulation of all known events that generate debris and objects that are part of the U.S. Space Surveillance Network (SSN) catalog, which provides coverage of objects with diameters down to approximately 10cm in LEO and 1 m in geosynchronous Earth orbit (GEO). Different models are used to simulate the artificial objects and their orbit evolution over time. These models are called “sources” since they assign an origin to each individual object and consist of fragments, solid rocket motor (SRM) remainders, sodium-potassium(NaK) droplets, paintflakes, ejecta, and multi-layer insulation (MLI) fragments. The objects from each source are characterized by having individual release mechanisms, as well as orbital distributions, material composition, size, and mass distributions. Dedicated radar and telescope observation data is used to calibrate the model for objects larger than 1 cm in LEO and larger than 10 cm in GEO. For calibrating the small-sized objects, below 1 cm, impact data from returned surfaces are analyzed. Because the >1cm object population is dominated by fragments, the fragmentation event database was updated to include new events, as well as re-evaluate past events. Special attention was drawn to re-evaluating theFengyun-1C anti-satellite test from 2007 and Cosmos-Iridium collision event from 2009 since these events shape the fragment population because of their severity. After 2009, the two largest fragmentations in terms of number of tracked debris are the Briz-M explosion in 2012 and the NOAA-16explosion in 2015. In total, there are 261 confirmed fragmentations in the database up to November 2016.The baseline population for ORDEM 3.1 is based on the U.S. SSN catalog, and observational datasets from radar, in situ, and optical sources provide a foundation from which the model populations are statistically extrapolated to smaller size regions. These regions are not well-covered by the SSN catalog yet may pose the greatest threat to operational spacecraft. The NASA Standard Satellite Breakup Model is used to generate fragments greater than 1mm from collisions and explosions, and these fragment populations are scaled using ground-based radar data. Specific major debris-producing events, including the Fengyun-1C, Iridium 33, and Cosmos2251 debris clouds, and unique populations, such as NaK droplets, were re-examined, modelled, and added to the ORDEM environment separately. Optical measurement data is used to model the GEO population down to 10 cm. The debris environment is propagated using NASA’s LEO-to-GEO Environment Debris model, and future explosions of intact objects and collisions involving objects greater than 10 cm are assessed statistically. The environment from a few millimetres down to 10 𝜇m is modelled using a special degradation model where small particles are generated from intact spacecraft and rocket bodies, then the populations are scaled to fit in situ cratering data from Space Shuttle returned surfaces. Fragments smaller than 10 cm are differentiated based on material density categories, i.e., high-, medium-, and low-density, to better characterize the potential debris risk posed to upper stages and spacecraft. This paper will discuss the MASTER and ORDEM approaches for modelling populations and compare fluxes for specific orbits, including sun-synchronous, ISS-altitude, geosynchronous transfer, and GEO. In the end, a conclusion is drawn towards the importance of having multiple fundamentally different, yet validated, models to estimate the space debris population.

Andre Horstmann↗

Computational Investigation of the Aerodynamic Effects on Fluidic Thrust Vectoring

A computational investigation of the aerodynamic effects on fluidic thrust vectoring has been conducted. Three-dimensional simulations of a two-dimensional, convergent-divergent (2DCD) nozzle with fluidic injection for pitch vector control were run with the computational fluid dynamics code PAB using turbulence closure and linear Reynolds stress modeling. Simulations were computed with static freestream conditions (M=0.05) and at Mach numbers from M=0.3 to 1.2, with scheduled nozzle pressure ratios (from 3.6 to 7.2) and secondary to primary total pressure ratios of p(sub t,s)/p(sub t,p)=0.6 and 1.0. Results indicate that the freestream flow decreases vectoring performance and thrust efficiency compared with static (wind-off) conditions. The aerodynamic penalty to thrust vector angle ranged from 1.5 degrees at a nozzle pressure ratio of 6 with M=0.9 freestream conditions to 2.9 degrees at a nozzle pressure ratio of 5.2 with M=0.7 freestream conditions, compared to the same nozzle pressure ratios with static freestream conditions. The aerodynamic penalty to thrust ratio decreased from 4 percent to 0.8 percent as nozzle pressure ratio increased from 3.6 to 7.2. As expected, the freestream flow had little influence on discharge coefficient.

Deere, K. A.↗

Status of Solar Sail Technology Within NASA

In the early 2000s, NASA made substantial progress in the development of solar sail propulsion systems for use in robotic science and exploration of the solar system. Two different 20-m solar sail systems were produced and they successfully completed functional vacuum testing in NASA Glenn Research Center's (GRC's) Space Power Facility at Plum Brook Station, Ohio. The sails were designed and developed by ATK Space Systems and L Garde, respectively. The sail systems consist of a central structure with four deployable booms that support the sails. These sail designs are robust enough for deployment in a one-atmosphere, one-gravity environment and were scalable to much larger solar sails perhaps as large as 150 m on a side. Computation modeling and analytical simulations were also performed to assess the scalability of the technology to the large sizes required to implement the first generation of missions using solar sails. Life and space environmental effects testing of sail and component materials were also conducted. NASA terminated funding for solar sails and other advanced space propulsion technologies shortly after these ground demonstrations were completed. In order to capitalize on the $30M investment made in solar sail technology to that point, NASA Marshall Space Flight Center (MSFC) funded the NanoSail-D, a subscale solar sail system designed for possible small spacecraft applications. The NanoSail-D mission flew on board the ill-fated Falcon-1 Rocket launched August 2, 2008, and due to the failure of that rocket, never achieved orbit. The NanoSail-D flight spare will be flown in the Fall of 2010. This paper will summarize NASA's investment in solar sail technology to-date and discuss future opportunities

Johnson, Les↗

A High-Resolution Merged Wind Dataset for DYNAMO: Progress and Future Plans

In order to support research on optimal data assimilation methods for the Cyclone Global Navigation Satellite System (CYGNSS), launching in 2016, work has been ongoing to produce a high‐resolution merged wind dataset for the Dynamics of the Madden Julian Oscillation (DYNAMO) field campaign, which took place during late 2011/early 2012. The winds are produced by assimilating DYNAMO observations into the Weather Research and Forecasting (WRF) three‐dimensional variational (3DVAR) system. Data sources from the DYNAMO campaign include the upper‐air sounding network, radial velocities from the radar network, vector winds from the Advanced Scatterometer (ASCAT) and Oceansat‐2 Scatterometer (OSCAT) satellite instruments, the NOAA High Resolution Doppler Lidar (HRDL), and several others. In order the prep them for 3DVAR, significant additional quality control work is being done for the currently available TOGA and SMART‐R radar datasets, including automatically dealiasing radial velocities and correcting for intermittent TOGA antenna azimuth angle errors. The assimilated winds are being made available as model output fields from WRF on two separate grids with different horizontal resolutions ‐ a 3‐km grid focusing on the main DYNAMO quadrilateral (i.e., Gan Island, the R/V Revelle, the R/V Mirai, and Diego Garcia), and a 1‐km grid focusing on the Revelle. The wind dataset is focused on three separate approximately 2‐week periods during the Madden Julian Oscillation (MJO) onsets that occurred in October, November, and December 2011. Work is ongoing to convert the 10‐m surface winds from these model fields to simulated CYGNSS observations using the CYGNSS End‐To‐End Simulator (E2ES), and these simulated satellite observations are being compared to radar observations of DYNAMO precipitation systems to document the anticipated ability of CYGNSS to provide information on the relationships between surface winds and oceanic precipitation at the mesoscale level. This research will improve our understanding of the future utility of CYGNSS for documenting key MJO processes.

Lang, Timothy J.↗