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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Models for multiple heating mechanisms

Current meteoritic and astrophysical evidence (especially FU Orionis events) suggests that large-scale gas-dynamics shocks penetrating to the nebular midplane represent a promising class of multiple heating mechanisms for the formation of chondrules. Initial numerical models have shown that shock heating of chondrule precursor grains in optically thick clouds of grains is especially efficient. Shocks generated by infalling molecular cloud clumps and shocks separating different mass elements in a time dependently and axially asymmetrically evolving disk are possibilities deserving further attention. Meteoritic constraints, astrophysical considerations and multiple heating mechanisms are addressed.

Hood, L. L.↗
NASA NTRS1987

Simulation of Airplane and Rocket Trajectories

Simulation and Optimization of Rocket Trajectories program (SORT) contains comprehensive mathematical models for simulating aircraft dynamics, freely falling objects, and many types of ballistic trajectories. Provides high-fidelity, three-degrees-of-freedom simulation for atmospheric and exoatmospheric flight. It numerically models vehicle subsystems and vehicle environment. Used for wide range of simulations. Written in machine-independent FORTRAN 77.

Wahbah, Magdy M.↗
NASA NTRS1991

Accretion and evolution of solar system bodies

We use a combination of analytical and numerical methods to study dynamical processes involved in the formation of planets and smaller bodies in the solar system. Our goal was to identify and understand critical processes and to link them in a numerical model of planetesimal accretion. We study effects of these processes by applying them in the context of the standard model of solar system formation, which involves accretion of the terrestrial planets and cores of the giant planet from small planetesimals. The principal focus of our research effort is the numerical simulation of accretion of a swarm of planetesimals into bodies of planetary size. Our computer code uses a Monte Carlo method to determine collisional interactions within the swarm. These interactions are not determined simply by a relative velocity, but rather by explicit distributions of keplerian orbital elements. The planetesimal swarm is divided into a number of zones in semimajor axis, which are allowed to interact. The present version of our code has the capability of following detailed distributions of size, eccentricity, and inclination in each zone.

Weidenschilling, S. J.↗
NASA NTRS1993

A consistent model reduction of measured modal parameters for reduced-order active control

The problem of synthesizing reduced-order linear models of vibrating structures for the design of fixed-order dynamic feedback control is investigated. The present technique builds on a recently developed procedure for constructing an objective set of mass and stiffness matrices from measured modal parameters that are akin to the Craig-Bampton synthesized ones obtained from finite element models. The constructed mass and stiffness matrices are determined directly from the identification of experimental data, however, rather than through correlation or reconciliation of a finite element model. A model truncation criterion is then applied to the identified minimum-order mass and stiffness model to satisfy certain observability/controllability requirements for the reduced model. Numerical examples illustrate the effectiveness of the proposed technique for synthesizing reduced-order controllers from system realizations of experimental data. The dynamic performance of the resulting closed-loop models is assessed using the known full-order structural dynamics and compared with existing model reduction techniques.

Alvin, K. F.↗
NASA NTRS1997

Turbulent Premixed Flame Propagation in Microgravity

A facility in which turbulent Couette flow could be generated in a microgravity environment was designed and built. To fit into the NASA Lewis drop tower the device had to be very compact. This means that edge effects and flow re-circulation were expected to affect the flow. The flow was thoroughly investigated using LDV and was found to be largely two dimensional away from the edges with constant turbulence intensities in the core. Slight flow asymmetries are introduced by the non symmetric re-circulation of the fluid outside the test region. Belt flutter problems were remedied by adding a pair of guide plates to the belt. In general, the flow field was found to be quite similar to previously investigated Couette flows. However, turbulence levels and associated shear stresses were higher. This is probably due to the confined re-circulation zone reintroducing turbulence into the test section. An estimate of the length scales in the flow showed that the measurements were able to resolve nearly all the length scales of interest. Using a new LES method for subgrid combustion it has been demonstrated that the new procedure is computational feasible even on workstation type environment. It is found that this model is capable of capturing the propagation of the premixed names by resolving the flame in the LES grid within 2-3 grid points. In contrast, conventional LES results in numerical smearing of the flame and completely inaccurate estimate of the turbulent propagation speed. Preliminary study suggests that there is observable effect of buoyancy in the 1g environment suggesting the need for microgravity experiments of the upcoming experimental combustion studies. With the cold flow properties characterized, an identical hot flow facility is under construction. It is assumed that the turbulence properties ahead of the flame in this new device will closely match the results obtained here. This is required since the hot facility will not enable LDV measurements. The reacting flow facility is also being constructed with planned drop tower experiments in mind. Therefore, issues related to safety and structural integrity are being take into account. Further development of the numerical model will also be carried out to include finite-rate kinetics for representative premixed cases. More detail analysis of the flame structure and propagation nature will be investigated. Simulations will also be compared to the flame properties observed in the experiments.

Menon, S.↗
NASA NTRS2005

Propagation of Interplanetary Disturbances in the Outer Heliosphere

Contents include the following: 1. We have developed a one-dimensional, spherically symmetric, multi-fluid MHD model that includes solar wind protons and electrons, pickup ions, and interstellar neutral hydrogen. This model advances the existing solar wind models for the outer heliosphere in two important ways: one is that it distinguishes solar wind protons from pickup ions, and the other is that it allows for energy transfer from pickup ions to the solar wind protons. Model results compare favorably with the Voyager 2 observations. 2. 2. Solar wind slowdown and interstellar neutral density. The solar wind in the outer heliosphere is fundamentally different from that in the inner heliosphere since the effects of interstellar neutrals become significant. 3. ICME propagation from the inner to outer heliosphere. Large coronal mass ejections (CMEs) have major effects on the structure of the solar wind and the heliosphere. The plasma and magnetic field can be compressed ahead of interplanetary CMEs. 4. During the current solar cycle (Cycle 23), several major CMEs associated with solar flares produced large transient shocks which were observed by widely-separated spacecraft such as Wind at Earth and Voyager 2 beyond 60 AU. Using data from these spacecraft, we use the multi-fluid model to investigate shock propagation and interaction in the heliosphere. Specifically, we studied the Bastille Day 2000, April 2001 and Halloween 2003 events. 5. Statistical properties of the solar wind in the outer heliosphere. In a collaboration with L.F. Burlaga of GSFC, it is shown that the basic statistical properties of the solar wind in the outer heliosphere can be well produced by our model. We studied the large-scale heliospheric magnetic field strength fluctuations as a function of distance from the Sun during the declining phase of a solar cycle, using our numerical model with observations made at 1 AU during 1995 as input. 6. Radial heliospheric magnetic field events. The heliospheric magnetic field (HMF) direction, on average, conforms well to the Parker spiral.

Wang, Chi↗
NASA NTRS2007

Optimized Heat Interception for Cryogen Tank Support

We consider means for using the cooling available in boil-off gas to intercept heat conducted through the support structure of a cryogen tank. A one-dimensional model of the structure coupled to a gas stream gives an analytical expression for heat leak in terms of flow rate for temperature independent-properties and laminar flow. A numerical model has been developed for heat transfer on a thin cylindrical tube with an attached vent line. The model is used to determine the vent path layout that will minimize heat flow into the cryogen tank. The results are useful for a number of applications, but the one of interest in this study is the minimization of the boil-off in large cryopropellant tanks in low Earth and low lunar orbit.

Canavan, Edgar R.↗
NASA NTRS2001

Inorganic Chlorine Partitioning in the Summer Lower Stratosphere: Modeled and Measured [ClONO2/HCl] During POLARIS

We examine inorganic chlorine (Cly) partitioning in the summer lower stratosphere using in situ ER-2 aircraft observations made during the Photochemistry of Ozone Loss in the Arctic Region in Summer (POLARIS) campaign. New steady state and numerical models estimate [ClONO2]/[HCl] using currently accepted photochemistry. These models are tightly constrained by observations with OH (parameterized as a function of solar zenith angle) substituting for modeled HO2 chemistry. We find that inorganic chlorine photochemistry alone overestimates observed [ClONO2]/[HCl] by approximately 55-60% at mid and high latitudes. On the basis of POLARIS studies of the inorganic chlorine budget, [ClO]/[ClONO2], and an intercomparison with balloon observations, the most direct explanation for the model-measurement discrepancy in Cly partitioning is an error in the reactions, rate constants, and measured species concentrations linking HCl and ClO (simulated [ClO]/[HCl] too high) in combination with a possible systematic error in the ER-2 ClONO2 measurement (too low). The high precision of our simulation (+/-15% 1-sigma for [ClONO2]/[HCl], which is compared with observations) increases confidence in the observations, photolysis calculations, and laboratory rate constants. These results, along with other findings, should lead to improvements in both the accuracy and precision of stratospheric photochemical models.

Voss, P. B.↗
NASA NTRS1992

Inertial oscillation of a vertical rotating draft with application to a supercell storm

An analytic model (vertical rotating draft) which includes the gross features of a supercell storm on an f-plane, undergoes an inertial oscillation that appears to have been overlooked in previous analytic and numerical models. The oscillation is nonlinear and consists of a long quiescent phase and a short intense phase. During the intense phase, the rotating draft has the following features of a supercell: the diameter of the core contracts as it spins up and expands as it spins down; if vertical wind shear is included, the track of the rotating draft turns to the right (an anticyclonic rotating draft turns to the left); this turning point is followed by a predominantly upward flow; and the horizontal pressure gradient is very small (a property of most tornadoless supercells). The rapid spin-up during the intense phase and the high Rossby numbers obtainable establish the ability of the Coriolis force to spin up single cyclonic or anticyclonic supercells by means of this inertial oscillation. This surprising result has implications for numerical supercell simulations, which generally do not rely on the Coriolis force as a source of rotation. The physics and mathematics of the inertial oscillation are given, and the solution is applied to a documented supercell.

Supercell thunderstorm↗
NASA NTRS1991

Basic studies of baroclinic flows

Computations were completed of transition curves in the conventional annulus, including hysteresis effect. The model GEOSIM was used to compute the transition between axisymmetric flow and baroclinic wave flow in the conventional annulus experiments. Thorough testing and documentation of the GEOSIM code were also completed. The Spacelab 3 results from the Geophysical Fluid Flow Cell (GFFC) were reviewed and numerical modeling was performed of many of the cases with horizontal temperature gradients as well as heating from below, with different rates of rotation. A numerical study of the lower transition to axisymmetric flow in the baroclinic annulus was performed using GEOSIM.

Miller, Tim L.↗
NASA NTRS1996

Bubble Dynamics on a Heated Surface

In this work, we study the combined thermocapillary and natural convective flow generated by a bubble on a heated solid surface. The interaction between gas and vapor bubbles with the surrounding fluid is of interest for both space and ground-based processing. On earth, the volumetric forces are dominant, especially, in apparatuses with large volume to surface ratio. But in the reduced gravity environment of orbiting spacecraft, surface forces become more important and the effects of Marangoni convection are easily unmasked. In order to delineate the roles of the various interacting phenomena, a combined numerical-experimental approach is adopted. The temperature field is visualized using Mach-Zehnder interferometry and the flow field is observed by a laser sheet flow visualization technique. A finite element numerical model is developed which solves the two-dimensional momentum and energy equations and includes the effects of bubble surface deformation. Steady state temperature and velocity fields predicted by the finite element model are in excellent qualitative agreement with the experimental results. A parametric study of the interaction between Marangoni and natural convective flows including conditions pertinent to microgravity space experiments is presented. Numerical simulations clearly indicate that there is a considerable difference between 1-g and low-g temperature and flow fields induced by the bubble.

Kassemi, Mohammad↗
NASA NTRS2015

NASA Galactic Cosmic Radiation Environment Model: Badhwar-O'Neill (2014)

The Badhwar‐O'Neill (BON) Galactic Cosmic Ray (GCR) flux model is used by NASA to certify microelectronic systems and in the analysis of radiation health risks for human space flight missions. Of special interest to NASA is the kinetic energy region below 4.0 GeV/n due to the fact that exposure from GCR behind shielding (e.g., inside a space vehicle) is heavily influenced by the GCR particles from this energy domain. The BON model numerically solves the Fokker‐Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration under the assumption of a spherically symmetric heliosphere. The model utilizes a GCR measurements database from various particle detectors to determine the boundary conditions. By using an updated GCR database and improved model fit parameters, the new BON model (BON14) is significantly improved over the previous BON models for describing the GCR radiation environment of interest to human space flight.

O'Neill, P. M.↗
NASA NTRS1999

Numerical Analysis of Temperature Gradients and Interface Shape During Directional Solidification of Al and Al-Cu Alloy Under Microgravity Conditions

Numerical modeling was undertaken to analyze the influence of radial thermal gradient on solid/liquid (s/1) interface shape and convection patterns during solidification of pure Al and Al-4 wt% Cu alloy. The objective of the numerical task was to predict the influence of convective velocity on an insoluble particle near a s/l interface. These predictions would then be used to define the minimum gravity level (g) required to investigate the fundamental physics of interaction between a particle and a s/I interface. To satisfy this objective, steady state calculations were performed for different gravity levels and orientations with the gravity vector. ne furnace configuration used in this analysis is the proposed International Space Station Furnace, Quench Module Insert (QMI) 1. Results from a thermal model of the furnace core were used as initial boundary conditions for solidification modeling. General model of binary alloy solidification was based on the finite element code FIDAP. It was found that for the worst case orientation of 90 degrees with the gravity vector and a g level of 10(exp -4)g(sub o) (g(sub o) = 9.8 m/s(exp 2)) the dominant forces acting on the particle would be the fundamental drag and interfacial forces.

Bune, Andris V.↗
NASA NTRS1999

Noninvasive assessment of left atrial maximum dP/dt by a combination of transmitral and pulmonary venous flow

OBJECTIVES: The study assessed whether hemodynamic parameters of left atrial (LA) systolic function could be estimated noninvasively using Doppler echocardiography. BACKGROUND: Left atrial systolic function is an important aspect of cardiac function. Doppler echocardiography can measure changes in LA volume, but has not been shown to relate to hemodynamic parameters such as the maximal value of the first derivative of the pressure (LA dP/dt(max)). METHODS: Eighteen patients in sinus rhythm were studied immediately before and after open heart surgery using simultaneous LA pressure measurements and intraoperative transesophageal echocardiography. Left atrial pressure was measured with a micromanometer catheter, and LA dP/dt(max) during atrial contraction was obtained. Transmitral and pulmonary venous flow were recorded by pulsed Doppler echocardiography. Peak velocity, and mean acceleration and deceleration, and the time-velocity integral of each flow during atrial contraction was measured. The initial eight patients served as the study group to derive a multilinear regression equation to estimate LA dP/dt(max) from Doppler parameters, and the latter 10 patients served as the test group to validate the equation. A previously validated numeric model was used to confirm these results. RESULTS: In the study group, LA dP/dt(max) showed a linear relation with LA pressure before atrial contraction (r = 0.80, p < 0.005), confirming the presence of the Frank-Starling mechanism in the LA. Among transmitral flow parameters, mean acceleration showed the strongest correlation with LA dP/dt(max) (r = 0.78, p < 0.001). Among pulmonary venous flow parameters, no single parameter was sufficient to estimate LA dP/dt(max) with an r2 > 0.30. By stepwise and multiple linear regression analysis, LA dP/dt(max) was best described as follows: LA dP/dt(max) = 0.1 M-AC +/- 1.8 P-V - 4.1; r = 0.88, p < 0.0001, where M-AC is the mean acceleration of transmitral flow and P-V is the peak velocity of pulmonary venous flow during atrial contraction. This equation was tested in the latter 10 patients of the test group. Predicted and measured LA dP/dt(max) correlated well (r = 0.90, p < 0.0001). Numerical simulation verified that this relationship held across a wide range of atrial elastance, ventricular relaxation and systolic function, with LA dP/dt(max) predicted by the above equation with r = 0.94. CONCLUSIONS: A combination of transmitral and pulmonary venous flow parameters can provide a hemodynamic assessment of LA systolic function.

NASA Discipline Cardiopulmonary↗
NASA NTRS2001

Flame-Vortex Interactions Imaged in Microgravity - To Assess the Theory Flame Stretch

The goals of this research are to: 1) Assess the Theory of Flame Stretch by operating a unique flame-vortex experiment under microgravity conditions in the NASA Glenn 2.2 Second Drop Tower (drops to identify operating conditions have been completed); 2) Obtain high speed shadowgraph images (500-1000 frames/s) using the drop rig (images were obtained at one-g, and the NASA Kodak RO camera is being mounted on the drop rig); 3) Obtain shadowgraph and PIV images at 1-g while varying the effects of buoyancy by controlling the Froude number (completed); 4) Numerically model the inwardly-propagating spherical flame that is observed in the experiment using full chemistry and the RUN 1DL code (completed); 5) Send images of the flame shape to Dr. G. Patniak at NRL who is numerically simulating the entire flame-vortex interaction of the present experiment (data transfer completed); and 6) Assess the feasibility of obtaining PIV velocity field images in the drop rig, which would be useful (but not required) for our assessment of the Theory of Flame Stretch (PIV images were obtained at one-g using same low laser power that is available from fiber optic cable in drop tower). The motivation for the work is to obtain novel measurement needed to develop a physically accurate model of turbulent combustion that can help in the control of engine pollutants. The unique experiment allows, for the first time, the detailed study of a negatively-curved (negatively stretched) flame, which is one of the five fundamental types of premixed flames. While there have been studies of flat flames, positively-curved (outwardly-propagating) cases and positively-strained (counterflow) cases, this is the first detailed study of a negatively-curved (inwardly-propagating) flame. The first set of drops in the 2.2 Second Drop Tower showed that microgravity provides more favorable conditions for achieving inwardly-propagating flames (IPFs) than 1-g. A vortex interacts with a flame and creates a spherical pocket, which burns inwardly. Shadowgraphs at 1000 frames/sec quantify the Markstein number and flame speed. A Low-Laser Power PIV System was developed and is being added to the drop package. Numerical computations were required to explain why the Markstein numbers measured for the inwardly-propagating flames differ from those of outward propagating flames; this is an important research issue in the assessment of the Theory of Flame Stretch. The RUN-1DL code (developed by Prof. B. Rogg) was run for IPF and OPFs with complex methane and propane chemistry. Results confirmed that Ma for the IPFs are larger than for OPFs as was observed experimentally. Physical reasons for these new findings about the Theory of Flame Stretch are being determined from the experiments and the computations. Several journal papers have been published; the drop package is described in the AIAA Journal, while the one-g results appear in three other journal papers.

Driscoll, James F.↗
NASA NTRS1973

A field investigation and numerical simulation of coastal fog

A field investigation of the microphysical and micrometeorological features of fogs occurring near Los Angeles and Vandenberg, California was conducted. Observations of wind speed and direction, temperature, dew point, vertical wind velocity, dew deposition, drop-size distribution, liquid water content, and haze and cloud nucleus concentration were obtained. These observations were initiated in late evening prior to fog formation and continued until the time of dissipation in both advection and radiation fogs. Data were also acquired in one valley fog and several dense haze situations. The behavior of these parameters prior to and during fog are discussed in detail. A two-dimensional numerical model was developed to investigate the formation and dissipation of advection fogs under the influence of horizontal variations in surface temperature. The model predicts the evolution of potential temperature, water vapor content, and liquid water content in a vertical plane as determined by vertical turbulent transfer and horizontal advection. Results are discussed from preliminary numerical experiments on the formation of warm-air advection fog and dissipation by natural and artificial heating from the surface.

Mack, E. J.↗
NASA NTRS1973

Behavior of nitric oxide formed by the space shuttle in the mesosphere

A numerical model of eddy diffusive transport of emitted gases from the space shuttle wake, including chemical reactions between the emitted constituents and the ambient atmosphere, has been constructed for 75 km altitude. The numerical methods involve explicit solution of the diffusion equation and Runge-Kutta method for the chemical reactions. The time required to reach background levels of nitric oxide concentration of 7 x 10 to the 7th power molecules/cc has been calculated. This relaxation time depends strongly on atmospheric conditions.

Eberstein, I. J.↗