Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,099 records · Page 61

Fine-Scale Comparison of TOMS Total Ozone Data with Model Analysis of an Intense Midwestern Cyclone

High-resolution (approx. 40 km) along-track total column ozone data from the Total Ozone Mapping Spectrometer (TOMS) instrument are compared with a high-resolution mesoscale numerical model analysis of an intense cyclone in the Midwestern United States. Total ozone increased by 100 DU (nearly 38%) as the TOMS instrument passed over the associated tropopause fold region. Complex structure is seen in the meteorological fields and compares well with the total ozone observations. Ozone data support the meteorological analysis showing that stratospheric descent was confined to levels above approx. 600 hPa; significant positive potential vorticity at lower levels is attributable to diabetic processes. Likewise, meteorological fields show that two pronounced ozone streamers extending north and northeastward into Canada at high levels are not bands of stratospheric air feeding into the cyclone; one is a channel of exhaust downstream from the system, and the other apparently previously connected the main cyclonic circulation to a southward intrusion of polar stratospheric air and advected eastward as the cut-off cyclone evolved. Good agreement between small-scale features in the model output and total ozone data underscores the latter's potential usefulness in diagnosing upper tropospheric/lower stratospheric dynamics and kinematics.

Olsen, Mark A.↗

Modeling and Simulation of Trapping Mechanisms of Granular Media in Space

This paper describes the modeling and simulation of trapped granular media, within the context of the Granular Imager project. We describe the physics of trapped granular media in space, and the methodologies used to stably confine and shape such a medium using electromagnetic fields. The numerical models have also been validated with results in the literature, obtaining excellent agreement. The results of the numerical tests indicate that it is possible, with structural arrangements of rings and plates at different levels of electrostatic potential, to stably confine one or more charged particles, when driven by voltages that can be modulated in time and space.

Quadrelli, Marco B.↗

Use of satellite data and modeling to assess the influence of stratospheric processes on the troposphere

Current research focused on the examination of the effects of ozone heating and Newtonian cooling on the linear stability of transient planetary waves. A finite difference numerical model was developed in which the basic state wind, temperature, and ozone distributions were constructed using ground-based and satellite data. Preliminary results indicate that ozone dynamics interaction may play a more important role in the stability and maintenance of planetary waves than previously thought, particularly in summer when the mean solar zenith angle is smallest and thus ozone heating largest. Examination of the finite amplitude interactions among radiation, ozone, and dynamics was carried out. Self-consistent, coupled equations governing the weakly nonlinear interactions are currently being analyzed to provide a better understanding of wave-mean flow interactions and ozone transport in a continuously stratified model of the troposphere-stratosphere coupled system.

Nathan, Terrence↗

Microphysics, Meteorology, Microwave and Modeling of Mediterranean Storms: The M(sup 5) Problem

Comprehensive understanding of the microphysical nature of Mediterranean storms requires a combination of in situ meteorological data analysis and radar-passive microwave data analysis, effectively integrated with numerical modeling studies at various scales, particularly from synoptic scale down to mesoscale. The microphysical properties of and their controls on severe storms are intrinsically related to meteorological processes under which storms have evolved, processes which eventually select and control the dominant microphysical properties themselves. Insofar as hazardous Mediterranean storms, highlighted by the September 25-28/1992 Genova flood event, the October 5-7/1998 Friuli flood event, and the October 13-15/2000 Piemonte flood event (all taking place in northern Italy), developing a comprehensive microphysical interpretation requires an understanding of the multiple phases of storm evolution and the heterogeneous nature of precipitation fields within the storm domains. This involves convective development, stratiform transition and decay, orographic lifting, and sloped frontal lifting proc esses. This also involves vertical motions and thermodynamical instabilities governing physical processes that determine details of the liquid/ice water contents, size distributions, and fall rates of the various modes of hydrometeors found within the storm environments. This paper presents detailed 4-dimensional analyses of the microphysical elements of the three severe Mediterranean storms identified above, investigated with the aid of SSM/I and TRMM satellite measurements (and other remote sensing measurements). The analyses are guided by nonhydrostatic mesoscale model simulations at high resolution of the intense rain producing portions of the storm environments. The results emphasize how meteorological controls taking place at the large scale, coupled with localized terrain controls, ultimately determine the most salient features of the bulk microphysical properties of the storms. These results have bearing on precipitation remote sensing from space, and the role of modeling in designing precipitation retrieval algorithms.

Smith, Eric A.↗

The Madden-Julian Oscillation and its Impact on Northern Hemisphere Weather Predictability during Wintertime

The Madden-Julian Oscillation (MJO) is known as the dominant mode of tropical intraseasonal variability and has an important role in the coupled-atmosphere system. This study used twin numerical model experiments to investigate the influence of the MJO activity on weather predictability in the midlatitudes of the Northern Hemisphere during boreal winter. The National Aeronautics and Space Administration (NASA) Goddard laboratory for the Atmospheres (GLA) general circulation model was first used in a 10-yr simulation with fixed climatological SSTs to generate a validation data set as well as to select initial conditions for active MJO periods and Null cases. Two perturbation numerical experiments were performed for the 75 cases selected [(4 MJO phases + Null phase) _ 15 initial conditions in each]. For each alternative initial condition, the model was integrated for 90 days. Mean anomaly correlations in the midlatitudes of the Northern Hemisphere (2O deg N_60 deg.N) and standardized root-mean-square errors were computed to validate forecasts and control run. The analyses of 500-hPa geopotential height, 200-hPa Streamfunction and 850-hPa zonal wind component systematically show larger predictability during periods of active MJO as opposed to quiescent episodes of the oscillation.

Jones, Charles↗

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program was to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, with the funding from this NASA program, we concentrated on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we developed unique global models that allowed us to study the coupling between the different regions. These results are highlighted in the next section. Another important aspect of our NASA Theory Program concerned the effect that localized 'structure' had on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkland current patterns) or time variations in these input due to storms and substorms. Also, some of the plasma flows that we predicted with our macroscopic models could be unstable, and another one of our goals was to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another goal of our NASA Theory Program was to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This could involve a detailed comparison of kinetic, semi-kinetic, and hydrodynamic predictions for a given polar wind scenario or it could involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations provides insight into when the various models can be used with confidence.

Schunk, Robert W.↗

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable, and another one of our goals is to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulation). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, R. W.↗

Gravity and suspension effects on the dynamics of controlled structures

The effects of gravity on the actuators, sensors and structural plant of a controlled flexible structure are investigated. These influences include gravity stiffening, gravity induced deformations, gravity influence on inertial sensors and actuators, and the interaction with the suspension system. For each of the four influences, a simplified model is derived, from which the appropriate non-dimensional parameter can be identified, and the magnitude of the influence estimated. The construction of a detailed numerical model, which includes the four gravity influences, is outlined. This modeling procedure is then applied to a simplified model of the Middeck Active Control Eperiment (MACE). The influence on the poles and input-output transfer function is examined for several variations in the plant. The modeling procedure is also applied to the MACE experimental hardware. Experimental and model-based transfer functions are compared to demonstrate the significant improvement by incorporating gravity effects.

Rey, Daniel A.↗

Microphysical Timescales in Clouds and their Application in Cloud-Resolving Modeling

Computational phenomena (i.e., spurious supersaturation and negative mixing ratio of cloud water) usually exist in cloud-resolving models when the time step for explicit integration is larger than a microphysical timescale in clouds. In this paper, the microphysical timescales in clouds are studied, showing that the timescale of water vapor condensation (or cloud water evaporation) is smaller than 10 s - the order of a typical time step for cloud-resolving models. To avoid spurious computational phenomena in cloud-resolving modeling, it is suggested that moist entropy be used as a prognostic thermodynamic variable, and temperature be diagnosed from that and other prognostic variables. A simple numerical model with moist entropy as a prognostic variable, for example, is presented to show that spurious computational phenomena are removed when moist entropy is used as a prognostic variable.

Zeng, Xi-Ping↗

Modeling of an Atmospheric-Boundary-Layer Profile in Support of Experiments in the NASA Langley Transonic Dynamics Tunnel

This paper presents the results of the Computational Fluid Dynamics (CFD) analyses to model two atmospheric-boundary-layer (ABL) profiles inside the NASA Langley Research Center Transonic Dynamics Tunnel (TDT). The CFD models include tunnel walls and all additional hardware needed to simulate the ABL profiles. This hardware includes Irwin spires and floor-mounted roughness elements. The numerical simulations show that the application of higher-fidelity numerical methods is necessary to compute boundary-layer and turbulenceintensity profiles that match experimental data. The ABL profiles are computed both inside the numerical model of the TDT and in a classical free-air model. Both Unsteady Reynolds Averaged Navier-Stokes (URANS) with Spalart-Allmaras (SA) turbulence model and Modified Delayed Detached Eddy (MDDES) simulation methods are used. The results show that the MDDES-simulation results match the experimental data very well while URANS-SA does not. The results also show that the wind-tunnel walls have a significant effect on ABL prediction.

Chwalowski, Pawel↗

Comment on 'Convergence characteristics of a vortex-lattice method for nonlinear configuration aerodynamics'

It is proposed that the study of Rusak et al. (1985), which reports numerical modeling sensitivities on longitudinal force/moment properties for a vortex-lattice method incorporating free vortex filaments to represent the leading-edge vortex separation, employs a formula that is strongly affected by the particular points of analysis chosen. This results in a narrowly applicable curve fit, where numerical sensitivities of the theory are inappropriately traded off against physical effects that are not modeled in that theory. Attention is also given to questionable drag estimate computations.

Luckring, J. M.↗

Using Coupled Models to Study the Effects of River Discharge on Biogeochemical Cycling and Hypoxia in the Northern Gulf of Mexico

We describe emerging capabilities to understand physical processes and biogeoehemical cycles in coastal waters through the use of satellites, numerical models, and ship observations. Emerging capabilities provide significantly improved ability to model ecological systems and the impact of environmental management actions on them. The complex interaction of physical and biogeoehemical processes responsible for hypoxic events requires an integrated approach to research, monitoring, and modeling in order to fully define the processes leading to hypoxia. Our efforts characterizes the carbon cycle associated with river plumes and the export of organic matter and nutrients form coastal Louisiana wetlands and embayments in a spatially and temporally intensive manner previously not possible. Riverine nutrients clearly affect ecosystems in the northern Gulf of Mexico as evidenced in the occurrence of regional hypoxia events. Less known and largely unqualified is the export of organic matter and nutrients from the large areas of disappearing coastal wetlands and large embayments adjacent to the Louisiana Continental Shelf. This project provides new methods to track the river plume along the shelf and to estimate the rate of export of suspended inorganic and organic paniculate matter and dissolved organic matter form coastal habitats of south Louisiana.

COUPLED MODEL↗

Verification of thermal analysis codes for modeling solid rocket nozzles

One of the objectives of the Solid Propulsion Integrity Program (SPIP) at Marshall Space Flight Center (MSFC) is development of thermal analysis codes capable of accurately predicting the temperature field, pore pressure field and the surface recession experienced by decomposing polymers which are used as thermal barriers in solid rocket nozzles. The objective of this study is to provide means for verifications of thermal analysis codes developed for modeling of flow and heat transfer in solid rocket nozzles. In order to meet the stated objective, a test facility was designed and constructed for measurement of the transient temperature field in a sample composite subjected to a constant heat flux boundary condition. The heating was provided via a steel thin-foil with a thickness of 0.025 mm. The designed electrical circuit can provide a heating rate of 1800 W. The heater was sandwiched between two identical samples, and thus ensure equal power distribution between them. The samples were fitted with Type K thermocouples, and the exact location of the thermocouples were determined via X-rays. The experiments were modeled via a one-dimensional code (UT1D) as a conduction and phase change heat transfer process. Since the pyrolysis gas flow was in the direction normal to the heat flow, the numerical model could not account for the convection cooling effect of the pyrolysis gas flow. Therefore, the predicted values in the decomposition zone are considered to be an upper estimate of the temperature. From the analysis of the experimental and the numerical results the following are concluded: (1) The virgin and char specific heat data for FM 5055 as reported by SoRI can not be used to obtain any reasonable agreement between the measured temperatures and the predictions. However, use of virgin and char specific heat data given in Acurex report produced good agreement for most of the measured temperatures. (2) Constant heat flux heating process can produce a much higher heating rate than the radiative heating process. The results show that heating rates of 125 C/s (225 F/s) can be achieved. (3) The electrical resistance of the FM 5055 samples were about 150 Omega in the virgin state, and decreased to about 35 Omega when charred. A reliable scheme must be developed to electrically insulate the composite from the heater foil in order to prevent any current leakage into the sample which can result in volumetric heating of the char zone of the composite.

Keyhani, M.↗

Convection equation modeling: A non-iterative direct matrix solution algorithm for use with SINDA

The determination of the boundary conditions for a component-level analysis, applying discrete finite element and finite difference modeling techniques often requires an analysis of complex coupled phenomenon that cannot be described algebraically. For example, an analysis of the temperature field of a coldplate surface with an integral fluid loop requires a solution to the parabolic heat equation and also requires the boundary conditions that describe the local fluid temperature. However, the local fluid temperature is described by a convection equation that can only be solved with the knowledge of the locally-coupled coldplate temperatures. Generally speaking, it is not computationally efficient, and sometimes, not even possible to perform a direct, coupled phenomenon analysis of the component-level and boundary condition models within a single analysis code. An alternative is to perform a disjoint analysis, but transmit the necessary information between models during the simulation to provide an indirect coupling. For this approach to be effective, the component-level model retains full detail while the boundary condition model is simplified to provide a fast, first-order prediction of the phenomenon in question. Specifically for the present study, the coldplate structure is analyzed with a discrete, numerical model (SINDA) while the fluid loop convection equation is analyzed with a discrete, analytical model (direct matrix solution). This indirect coupling allows a satisfactory prediction of the boundary condition, while not subjugating the overall computational efficiency of the component-level analysis. In the present study a discussion of the complete analysis of the derivation and direct matrix solution algorithm of the convection equation is presented. Discretization is analyzed and discussed to extend of solution accuracy, stability and computation speed. Case studies considering a pulsed and harmonic inlet disturbance to the fluid loop are analyzed to assist in the discussion of numerical dissipation and accuracy. In addition, the issues of code melding or integration with standard class solvers such as SINDA are discussed to advise the user of the potential problems to be encountered.

Schrage, Dean S.↗

Mechanical energy flow models of rods and beams

It has been proposed that the flow of mechanical energy through a structural/acoustic system may be modeled in a manner similar to that of flow of thermal energy/in a heat conduction problem. If this hypothesis is true, it would result in relatively efficient numerical models of structure-borne energy in large built-up structures. Fewer parameters are required to approximate the energy solution than are required to model the characteristic wave behavior of structural vibration by using traditional displacement formulations. The energy flow hypothesis is tested in this investigation for both longitudinal vibration in rods and transverse flexural vibrations of beams. The rod is shown to behave approximately according to the thermal energy flow analogy. However, the beam solutions behave significantly differently than predicted by the thermal analogy unless locally-space-averaged energy and power are considered. Several techniques for coupling dissimilar rods and beams are also discussed. Illustrations of the solution accuracy of the methods are included.

Wohlever, J. C.↗

The r-process and neutrino-heated supernova ejecta

As a neutron star is formed by the collapse of the iron core of a massive star, its Kelvin-Helmholtz evolution is characterized by the release of gravitational binding energy as neutrinos. The interaction of these neutrinos with heated material above the neutron star generates a hot bubble in an atmosphere that is nearly in hydrostatic equilibrium and heated, after approximately 10 s, to an entropy of S/N(sub AS)k greater than or approximately = 400. The neutron-to-proton ratio for material moving outward through this bubble is set by the balance between neutrino and antineutrino capture on nucleons. Because the electron antineutrino spectrum at this time is hotter than the electron neutrino spectrum, the bubble is neutron-rich (0.38 less than or approximately = Y(sub e) less than or approximately = 0.47). Previous work using a schematic model has shown that these conditions are well suited to the production of heavy elements by the r-process. In this paper we have advanced the numerical modeling of a 20 solar mass 'delayed' supernova explosion to the point that we can follow the detailed evolution of material moving through the bubble at the late times appropiate to r-process nucleosynthesis. The supernova model predicts a final kinetic energy for the ejecta of 1.5 x 10(exp 51) ergs and leaves behind a remnant with a baryon mass of 1.50 solar mass (and a gravitational mass of 1.445 solar mass). We follow the thermodynamic and compositional evolution of 40 trajectories in rho(t), T(t), Y(sub e)(t) for a logarithmic grid of mass elements for the last approximately = 0.03 solar mass to be ejected by the proto-neutron star down to the last less than 10(exp -6) solar mass of material expelled at up to approximately = 18 s after core collapse. We find that an excellent fit to the solar r-process abundance distribution is obtained with no adjustable parameters in the nucleosynthesis calculations. Moreover, the abundances are produced in the quantities required to account for the present Galactic abundances. However, at earlier times, this one-dimensional model ejects too much material with entropies S/N(sub A)k approximately 50 and Y(sub e) approximately 0.46. This leads to an acceptable over production of N = 50 nuclei, particularly Sr-88, Y-89, and Zr-90, relative to their solar abundances. We speculate on various means to avoid the early overproduction and/or ejection of N = 50 isotonic nuclei while still producing and ejecting the correct amount of r-process material.

Woosley, S. E.↗

Concerning the Charging of an Exploration Craft on and near a Small Asteroid

Introduction: An object immersed in an airless plasma environment will experience a natural process of surface charging in order to acheieve current balance, or zero net electric current to the object. It has been shown in recent computer simulations that the small-body plasma environment is very complex [1], considering effects of photoemission, topography, and formation of a plasma wake. For this work we consider an exploration craft (or astronaut) immersed within a plasma environment near an asteroid, which exhibits widely varying solar wind and photoelectric particle fluxes and continuously evolving illumination conditions. Objective: We aim to determine how an explo-ration craft or astronaut suit accumulates charge while located in the "nightside" asteroid wake where the particle fluxes are reduced, and in the dayside near-surface photoelectron sheath, by combining an object charging model [2] with kinetic simulations of a near-asteroid plasma environment [1]. We consider an astronaut floating near the asteroid while not in contact with the surface, as well as an astronaut moving along the surface using their hands/gloves to crawl along. Results: The modeling results suggest that remediation of triboelectric charge via accumulation of plasma currents is an important factor to consider when designing future NEA mission infrastructure, especially if repeated and frequent contact with the surface is planned. In shadowed regions such as the location shown in Fig. 1a, the plasma currents are so low (and the effective charge-remediation timescale so long, e.g. minutes to hours) that repeated contact with the surface tribocharges the glove in an uncontrollable fashion, as shown for two representative electron temperatures in Fig. 2a. The resulting buildup of significant negative charge would eventually initiate some other "current of last resort" [4] such as transport of positively-charged dust, field-emission from the glove, or significant alteration of environmental ion currents within the wake. In contrast, the few-meters-thick dayside photoelectron sheath in which the astronaut of Fig. 1b is immersed in is so rich in electrons (and hence so electrically conductive) that accumulated tribocharge dissipates almost instantaneously (e.g. in less than a ms) as shown in Fig. 2b. As our model astronaut orbits the NEA they would experience plasma currents and associated charge re-mediation times spanning many orders of magnitude, and the fusion between our numerical models provides a detailed understanding of the charging hazards possibly associated with contact-based NEA exploration.

Asteroid↗

Meteorology, Macrophysics, Microphysics, Microwaves, and Mesoscale Modeling of Mediterranean Mountain Storms: The M8 Laboratory

Comprehensive understanding of the microphysical nature of Mediterranean storms can be accomplished by a combination of in situ meteorological data analysis and radar-passive microwave data analysis, effectively integrated with numerical modeling studies at various scales, from synoptic scale down through the mesoscale, the cloud macrophysical scale, and ultimately the cloud microphysical scale. The microphysical properties of and their controls on severe storms are intrinsically related to meteorological processes under which storms have evolved, processes which eventually select and control the dominant microphysical properties themselves. This involves intense convective development, stratiform decay, orographic lifting, and sloped frontal lifting processes, as well as the associated vertical motions and thermodynamical instabilities governing physical processes that affect details of the size distributions and fall rates of the various types of hydrometeors found within the storm environment. Insofar as hazardous Mediterranean storms, highlighted in this study by three mountain storms producing damaging floods in northern Italy between 1992 and 2000, developing a comprehensive microphysical interpretation requires an understanding of the multiple phases of storm evolution and the heterogeneous nature of precipitation fields within a storm domain. This involves convective development, stratiform transition and decay, orographic lifting, and sloped frontal lifting processes. This also involves vertical motions and thermodynamical instabilities governing physical processes that determine details of the liquid/ice water contents, size disi:ributions, and fall rates of the various modes of hydrometeors found within hazardous storm environments.

Starr, David O.↗