Search NASA⌕ Search

SEARCH · Search NASA

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 343 records · Page 19

Numerical modeling of rotor flows with a conservative form of the full-potential equations

A computer program has been developed to solve a three-dimensional conservative formulation of the full-potential equation. Its ability to solve transonic, unsteady rotor flows is demonstrated by comparison to forward flight non-lifting pressure data at low to moderate advance ratios. A 'split potential' formulation has been added to the code which incorporates known vorticity fields into the full-potential calculation. Using this methodology, rotor wake contributions have been incorporated into the computer code. Pressure results are presented for lifting rotors in hover. These results are compared to experimental data as well as to other predictions.

Strawn, R. C.↗

A numerical model of magnetosphere-ionosphere coupling Preliminary results

A three-dimensional simulation model was developed to study magnetosphere-ionosphere coupling in the auroral region. One-fluid MHD equations are adopted to model the magnetosphere, and current density continuity equations were solved consistently to model the ionosphere. In the preliminary simulation runs described here, the electrodynamics of region 1 field-aligned currents were modeled. Initially, the electric field is taken to be in the magnetospheric equatorial plane. Alfven waves then propagate down to the ionosphere, accompanying a field-aligned current and exciting the electrostatic potential in the ionosphere by electrodynamic coupling. The results of the preliminary runs directly correspond with the fundamental characteristics of global magnetosphere-ionosphere coupling. These characteristics include an ionospheric electrostatic potential which varies in its development in time, depending upon the ratio of ionospheric resistance to magnetospheric impedance, the flowing of field-aligned currents into the ionosphere on the dawnside and out on the duskside, and the distribution along the geomagnetic field of the amplitude of the field-aligned current density in proportion to the field intensity.

Watanabe, K.↗

A detailed numerical model of a supersonic reacting mixing layer

A current research effort is underway at the NASA Langley Research Center to achieve a detailed understanding of important phenomena present when a supersonic flow undergoes chemical reaction. A computer program has been developed to study the details of such flows. The program has been constructed to consider the multicomponent diffusion and convection of important species, the finite-rate reaction of these species, and the resulting interaction between the fluid mechanics and chemistry. Code results from the analysis of a spatially developing and reacting mixing layer are presented, and conclusions are drawn regarding the structure of the evolving layer and its associated flame.

Drummond, J. P.↗

A feasibility study on the numerical modeling of interior noise fields

New developments in electronics and computer control have made active methods of noise suppression feasible. To date, however, active methods of noise control have been applied principally in geometrically uniform spaces with localized sound sources. The purpose of this paper is to demonstrate how finite element techniques may be used to model sound fields and active control inside geometrically complex spaces with non-uniform acoustic boundaries and multiple spatially separated sources. A closed-form solution is derived for the case of a one-dimensional tube with a vibrating piston at each end. Results are compared with finite element solutions and the problem is then generalized to two-dimensions with non-uniform geometry and non-uniform acoustic boundaries.

Abrahamson, A. L.↗

Comparison of vertical velocities analyzed by a numerical model and measured by a VHF wind profiler

The use of wind profilers for measuring vertical velocities in the troposphere and lower stratosphere is potentially of great interest for verification of forecasts, diagnosis of mesoscale circulations, and studies of wave motions. The studies of profiler vertical velocities to date have shown that the observed patterns of ascent and subsidence are reasonable when compared to the synoptic conditions. However, difficulties arise when a direct verification of the profiler vertical winds is sought. Since no other technique can measure the vertical velocities over the same height range and with the same claimed accuracy as the profilers, direct comparisons are impossible. The only alternative is to compare the measurements to analyzed vertical velocity fields. Here, researchers compare vertical measurements made with the SOUSY VHF radar over a period of 11 days at the beginning of November 1981 to the analyzed vertical velocities produced by the European Center for Medium-range Weather Forecasting (ECMWF) model for grid points near the radar site.

Larsen, M. F.↗

Intercomparison of numerical models of flaring coronal loops

The proposed Benchmark Problem consists of an infinitesimal magnetic flux tube containing a low-beta plasma. The field strength is assumed to be so large that the plasma can move only along the flux tube, whose shape remains invariant with time (i.e., the fluid motion is essentially one-dimensional). The flux tube cross section is taken to be constant over its entire length. In planar view the flux tube has a semi-circular shape, symmetric about its midpoint s = s sub max and intersecting the chromosphere-corona interface (CCI) perpendicularly at each foot point. The arc length from the loop apex to the CCI is 10,000 km. The flux tube extends an additional 2000 km below the CCI to include the chromosphere, which initially has a uniform temperature of 8000 K. The temperature at the top of the loop was fixed initially at 2 X 1 million K. The plasma is assumed to be a perfect gas (gamma = 5/3), consisting of pure hydrogen which is considered to be fully ionized at all temperatures. For simplicity, moreover, the electron and ion temperatures are taken to be everywhere equal at all times (corresponding to an artificially enhanced electron-ion collisional coupling). While there was more-or-less unanimous agreement as to certain global properties of the system behavior (peak temperature reached, thermal-wave time scales, etc.), no two groups could claim satisfactory accord when a more detailed comparison of solutions was attempted.

Kopp, R. A.↗

Numerical modelling of cryogenic propellant behavior in low-G

A partial survey is presented of recent research, sponsored by the NASA Lewis Research Center, into the computational modelling of cryogenic propellant behavior in a low gravity environment. This presentation is intended to provide insight into some of the specific problems being studied and into how these studies are part of an integrated plan to develop predictive capabilities. A brief description of the computational models developed to analyze jet induced mixing in cryogenic propellant tankage is presented along with representative results. Similar information is presented for a recent examination of on-orbit self-pressurization. A study of propellant reorientation has recently been initiated and preliminary results are included. The presentation concludes with a list of ongoing efforts and projected goals.

Hochstein, John I.↗

A surface temperature and moisture parameterization for use in mesoscale numerical models

A modified multi-level soil moisture and surface temperature model is presented for use as in defining lower boundary conditions in mesoscale weather models. Account is taken of the hydraulic and thermal diffusion properties of the soil, their variations with soil type, and the mixing ratio at the surface. Techniques are defined for integrating the surface input into the multi-level scheme. Sample simulation runs were performed with the modified model and the original model defined by Pielke, et al. (1977, 1981). The models were applied to regional weather forecasting over soils composed of sand and clay loam. The new form of the model avoided iterations necessary in the earlier version of the model and achieved convergence at reasonable profiles for surface temperature and moisture in regions where the earlier version of the model failed.

Tremback, C. J.↗

Time-dependent numerical modeling of dust halo formation at comets

The evolution of gas and dust distributions following a spatially and temporally localized comet outburst was calculated using a hybrid kinetic-hydrodynamic method. In the inner coma the time-dependent continuity, momentum, and energy equations of the dusty gas flow were solved simultaneously using 12 dust sizes. Beyond 300 km a three-dimensional kinetic model was used to calculate the trajectory of each individual dust grain. It was found that following the onset of the comet outburst a gas-dust blast wave propagates outward in the inner coma. About 15 minutes after the increased gas and dust production was initiated at the nucleus, a new equilibrium was reached in the inner coma. The most important feature of this new steady state was the significant increase of the dust terminal velocities. These higher terminal velocity values resulted in larger apex distances for dust particles emitted during the outburst. The dust particles spend a relatively long time near their apex points; therefore, the outburst generates long-lasting distinct dust envelopes in front of the regular dust coma.

Gombosi, T. I.↗

Numerical modeling of the vortex/airfoil interaction

A modeling of the vortex-airfoil interaction is presented in which the finite-area of the real vortices is taken into consideration. Two vortex models are used. In the first, a disturbed piece of vorticity layer is simulated by four rows of discrete vortices of small strength. In the second, a number of discrete vortices is arranged within a circle. The first model may simulate a shear layer or a wake, while the second, a well-formed vortex. The method was applied to the calculation of the pressure induced on the surface of the airfoil by the interacting vortex. Both models give similar results. It was found that for large distances of the vortex from the surface of the airfoil, the consideration or not of the finite-area of the vortex is not a significant factor in determining the induced pressure field. However, when the distance of the vortex from the surface is reduced, its shape is distorted and the induced pressure pulses have lower amplitude than the ones induced by an equivalent point vortex. In the limit, where the vortex impinges on the leading edge of the airfoil, it is split into two and the time dependent pressure coefficient takes even negative values at some time intervals.

Panaras, Argyris G.↗

Numerical modeling of proton and electron bursts in a realistic magnetotail

Particle trajectories of the magnetotail ion and electron bursts were simulated using the Beard-Hirschi-Propp (1982) magnetic field model and data generated by IMP 8. Particle trajectories were computed both in the absence of the electric field, and with electric field introduced locally perpendicular to the magnetic field. The results indicate that the motions of energetic particles near the nominal neutral sheet can be divided into two types. Proton and electron parallel trajectories have most of their trajectory motion parallel to the magnetic field lines and travel long distances along field lines. Perpendicular particle trajectories travel much shorter distances with much of the trajectory motion perpendicular to the magnetic field and with proton trajectories interacting with the neutral sheet region. Only the premidnight trajectories connect to reasonable source locations, which explains the preferential occurrence of earthward energetic proton burst anisotropies in the premidnight sector.

Kutchko, Frank J.↗

A numerical modeling study of a Montana thunderstorm. I - Model results versus observations involving nonelectrical aspects. II - Model results versus observations involving electrical aspects

Model results and the observed cloud behavior are examined in terms of nonelectrical and electrical aspects of a thunderstorm. The characteristics of the two-dimensional, time-dependent atmospheric electricity model used to simulate the cloud observations of July 19, 1981 in Miles City, Montana are described. The interactions of the dynamics and microphysics of the cloud with the charge separation mechanisms are analyzed. It is observed that the model accurately represents many of the observed characteristics of the cloud; however, the cloud base height, maximum liquid water content, and the time from first formation of precipitation until it reaches the ground are not accurately modeled. It is found that the model adequately represents the electrical field structure of the cloud and the electrical field strengths.

Helsdon, John H., Jr.↗

Comments on the challenge of using mesoscale data in mesoscale numerical models

The dangers of addressing the initialization issues for limited-area mesoscale models by extending the lessons learned during the development of global analysis and prediction systems are discussed. Lack of impact with data inserts at one time suggests that the lateral boundary conditions imposed on the limited-area models might force the model simulation toward a preferred solution, work against the new data being inserted into the model and, therefore, limit the potential impact that this data can have on the model system. The second potential pitfall involves the imposition of balance constraints on the data that are being inserted into the model to compute winds from temperature data and/or temperature from wind data.

Uccellini, Louis W.↗

Numerical modelling of the classical nova outburst

A mechanism is described that promises to explain how nova outbursts take place on white dwarf of 1 solar mass or less and for accretion rates of 4 x 10 to the -10 solar mass/yr or greater.

Kutter, G. S.↗

A Numerical Model of Unsteady, Subsonic Aeroelastic Behavior

A method for predicting unsteady, subsonic aeroelastic responses was developed. The technique accounts for aerodynamic nonlinearities associated with angles of attack, vortex-dominated flow, static deformations, and unsteady behavior. The fluid and the wing together are treated as a single dynamical system, and the equations of motion for the structure and flow field are integrated simultaneously and interactively in the time domain. The method employs an iterative scheme based on a predictor-corrector technique. The aerodynamic loads are computed by the general unsteady vortex-lattice method and are determined simultaneously with the motion of the wing. Because the unsteady vortex-lattice method predicts the wake as part of the solution, the history of the motion is taken into account; hysteresis is predicted. Two models are used to demonstrate the technique: a rigid wing on an elastic support experiencing plunge and pitch about the elastic axis, and an elastic wing rigidly supported at the root chord experiencing spanwise bending and twisting. The method can be readily extended to account for structural nonlinearities and/or substitute aerodynamic load models. The time domain solution coupled with the unsteady vortex-lattice method provides the capability of graphically depicting wing and wake motion.

Strganac, Thomas W.↗

The use of available potential energy to evaluate the impact of satellite data on numerical model analysis during FGGE

To evaluate the effect of the FGGE satellite observing system, the following two data sets were compared by examining the available potential energy (APE) and extratropical cyclone activity within the entire global domain during the first Special Observing Period: (1) the complete FGGE IIIb set, which incorporates satellite soundings, and (2) a NOSAT set which incorporates only conventional data. The time series of the daily total APEs indicate that NOSAT values are larger than the FGGE values, although in the Northern Hemisphere the differences are negligible. Analyses of cyclone scale features revealed only minor differences between the Northern Hemisphere FGGE and NOSAT analyses. On the other hand, substantial differences were revealed in the two Southern Hemisphere analyses, where the satellite soundings apparently add detail to the FGGE set.

Horn, Lyle H.↗

A numerical model for supersonic reacting mixing layers

A current research effort is underway at the NASA Langley Research Center to achieve a detailed understanding of important phenomena present when a supersonic flow undergoes a chemical reaction. A computer program has been developed to study the details of such flows. The program has been constructed to consider the multicomponent diffusion and convection of important species, the finite-rate reaction of these species, and the resulting interaction between the fluid mechanics and chemistry. Code results from the analysis of a spatially developing and reacting mixing layer are presented, and conclusions are drawn regarding the structure of the evolving layer and its associated flame.

Drummond, J. Philip↗

Review and assessment of the database and numerical modeling for turbine heat transfer

The objectives of the HOST Turbine Heat Transfer subproject were to obtain a better understanding of the physics of the aerothermodynamic phenomena and to assess and improve the analytical methods used to predict the flow and heat transfer in high-temperature gas turbines. At the time the HOST project was initiated, an across-the-board improvement in turbine design technology was needed. A building-block approach was utilized and the research ranged from the study of fundamental phenomena and modeling to experiments in simulated real engine environments. Experimental research accounted for approximately 75 percent of the funding while the analytical efforts were approximately 25 percent. A healthy government/industry/university partnership, with industry providing almost half of the research, was created to advance the turbine heat transfer design technology base.

Gladden, H. J.↗