Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Simulation”

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

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

At least 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Mature thunderstorm cloud-top structure and dynamics - A three-dimensional numerical simulation study

The present investigation is concerned with results from an initial set of comparative experiments in a project which utilize a three-dimensional convective storm model. The modeling results presented are related to four comparative experiments, designated Cases A through D. One of two scientific questions considered involves the dynamical processes, either near the cloud top or well within the cloud interior, which contribute to organize cloud thermal patterns such as those revealed by IR satellite imagery for some storms having strong internal cloud-scale rotation. The second question is concerned with differences, in cloud-top height and temperature field characteristics, between thunderstorms with and without significant internal cloud-scale rotation. The four experiments A-D are compared with regard to both interior and cloud-top configurations in the context of the second question. A particular strong-shear experiment, Case B, is analyzed to address question one.

Schlesinger, R. E.↗
NASA NTRS1984

A preliminary study of numerical simulation of thermosolutal convection of interest to crystal growth

Calculations were performed with computer models using three types of finite difference methods of thermosolutal convection: horizontal heating of a container filled with a stably stratified solution, finger convection in a container, and finger convection in a horizontally infinite channel. The importance of including thermosolutal convection in models of crystal growth is emphasized, and the difficulties in doing so are demonstrated. It is pointed out that these difficulties, due primarily to the fine structure of the convection, may be partly overcome by the use of fine grids and implicit time stepping methods.

Miller, T. L.↗
NASA NTRS1984

Numerical simulations of fission

Hydrodynamic simulation techniques to the fission problem of the lunar origin were applied. It is indicated that, for fluids with the compressibility of stars, dynamic fission instabilities lead to spiral arm ejection of mass and angular momentum in the form of a ring or disk of debris, not as a single body. Some quantitative aspects of these results seem favorable to lunar origin by fission. Fission instabilities in fluid objects with a smaller degree of compressibility, more closely approximating terrestrial material are considered. Although the qualitative features are similar, there are significant quantitative differences for the stiffer equation of state. The implications of the results for the fission hypothesis of lunar origin are discussed. Evolutions illustrating possible approaches to instability are presented.

Durisen, R. H.↗
NASA NTRS1984

Numerical Simulation of North Atlantic Sea Ice Variability, 1951 - 1980

A two-level dynamic-thermodynamic sea ice model is used to simulate the growth, drift and decay of sea ice in the Northern Hemisphere during a 30-year period, 1951 to 1980. The model is run with a daily timestep on a 222 km grid and is forced by interanually varying fields of geostrophic wind and temperature-derived thermodynamic fluxes. The objective is a quantitative description of large-scale sea ice variability in terms of the dynamic and thermodynamic processes responsible for the fluctuations, especially in the North Atlantic where sea ice represents a substantial input of fresh water. The fields of ice velocity and thickness contain strong seasonal as well as interannual variability. The mean drift pattern results in thicknesses of 4 to 5 m offshore of northern Canada and Greenland, while winter thicknesses of approximately 2 m are typical of Alaskan. Eurasian and East Greenland coastal waters. The 30-year mean fields are characterized by ecessive ice in the North Atlantic during winter and by a summer retreat that is more rapid than observed.

Walsh, J. E.↗
NASA NTRS1985

A numerical simulation of the NFAC (National Full-scale Aerodynamics Complex) open-return wind tunnel inlet flow

The flow into an open return wind tunnel inlet was simulated using Euler equations. An explicit predictor-corrector method was employed to solve the system. The calculation is time-accurate and was performed to achieve a steady-state solution. The predictions are in reasonable agreement with the experimental data. Wall pressures are accurately predicted except in a region of recirculating flow. Flow-field surveys agree qualitatively with laser velocimeter measurements. The method can be used in the design process for open return wind tunnels.

Kaul, U. K.↗
NASA NTRS1984

Numerical simulation of the transonic flowfield for wing/nacelle configurations

An efficient grid-interfacing zonal algorithm has been developed for computing the three-dimensional transonic flow field about wing/nacelle multicomponent configurations. The algorithm uses the full-potential formulation and the AF2 fully-implicit approximate factorization scheme. The flow field position is computed using a component-adaptive grid approach in which separate grids are employed for the individual components in the multicomponent configuration, where each component grid is optimized for a particular geometry such as the wing or nacelle. The wing and nacelle component grids are allowed to overlap, and flow field information is transmitted from one grid to another through the overlap region using trivariate interpolation. This paper presents a discussion of the computational methods used to generate both the wing and nacelle component grids, the technique used to interface the component grids, and the method used to obtain the inviscid multicomponent flow field solution. Computed results and correlations with experiment are presented to illustrate application of the analysis.

Atta, E. H.↗
NASA NTRS1983

Numerical simulations of plasma double layers

The results of analytical studies of quasi-static electric fields along geomagnetic field lines are discussed. The calculations were targeted at the structure, generation mechanisms and stability parameters. The field consists of two oppositely charged layers, either weakly or strongly charged, with an electric field between. Existence conditions are defined for the double layer field and balancing requirements are explored. Details of the simulation techniques, i.e., particle in cell and Vlasov simulations, for studying the double layer are outlined, noting that both periodic and quasi-periodic simulations are used. Solutions to Poisson's equation for fixed and floating point boundary conditions are generated. Finally, attention is also given to oblique and two-dimensional magnetic double layers.

Goertz, C. K.↗
NASA NTRS1985

A numerical simulation of auroral ionospheric electrodynamics

A computer simulation of auroral ionospheric electrodynamics in the altitude range 80 to 250 km has been developed. The routine will either simulate typical electron precipitation profiles or accept observed data. Using a model background ionosphere, ion production rates are calculated from which equilibrium electron densities and the Hall and Pedersen conductivities may be determined. With the specification of suitable boundary conditions, the entire three-dimensional current system and electric field may be calculated within the simulation region. The results of the application of the routine to a typical inverted-V precipitation profile are demonstrated. The routine is used to explore the observed anticorrelation between electric field magnitude and peak energy in the precipitating electron spectrum of an auroral arc.

Mallinckrodt, A. J.↗
NASA NTRS1985

The numerical simulation of steady transonic rotational flow using a dual potential formulation

A finite-difference method is presented that simulates steady transonic rotational flow of an inviscid fluid by representing the velocity field as the sum of scalar and vector potentials. This dual potential velocity decomposition extends the validity of the scalar (full) velocity potential to include vorticity. The inclusion of a vector potential also permits an alternate treatment of lift that does not require a circulation wake cut. This is accomplished by specifying the vector potential as a constant on the airfoil surface in order to satisfy a Kutta condition. The governing equations are solved as iteratively decoupled scalar equations using approximate factorization techniques, and the overall efficiency approaches that of the full potential equation. The governing equations are able to convect entropy and vorticity throughout the flow field and are equivalent to the Euler equations in continuous flow domains, however at shocks the Rankine-Hugoniot entropy jump must be supplied. An entropy correction method is presented and verified with transonic airfoil solutions of the Euler equations.

Chaderjian, N. M.↗
NASA NTRS1985

A numerical simulation of the NFAC (National Full-scale Aerodynamics Complex) open-return wind tunnel inlet flow

The flow into an open return wind tunnel inlet was simulated using Euler equations. An explicit predictor-corrector method was employed to solve the system. The calculation is time-accurate and was performed to achieve a steady-state solution. The predictions are in reasonable agreement with the experimental data. Wall pressures are accurately predicted except in a region of recirculating flow. Flow-field surveys agree qualitatively with laser velocimeter measurements. The method can be used in the design process for open-return wind tunnels.

Kaul, U. K.↗
NASA NTRS1984

Some features of auroral electric fields as seen in 2D numerical simulations

Results of 2D plasma simulations are presented and related to auroral observations. The formation of V-shaped potentials is studied with a 2 1/2 dimensional electrostatic particle-in-cell code for a magnetized plasma. It is shown that amplitudes for perpendicular electric fields are larger than for parallel electric fields, and for Te less than 100 eV, the amplitudes are comparable to the electric fields associated with the electrostatic shocks observed from the S3-3 satellite. The excitation of electrostatic ion-cyclotron EIC waves which occurs in the region below the parallel potential drop is discussed. In auroral plasmas EIC waves are observed above the V-shaped double layers in association with ion beams and field-aligned currents. The results also show that oppositely directed electric fields in the center and at the edges of the simulation region produce oppositely directed currents. Precipitating auroral ions in association with electron inverted-V events are seen by the DMSP-F6 satellite.

Thiemann, H.↗
NASA NTRS1985

Numerical simulation of positive-potential conductors in the presence of a plasma and a secondary-emitting dielectric

A 2 1/2-dimensional, cylindrically symmetric particle in cell program is used to simulate the interaction of a plasma with a positively biased disk covered by a dielectric material (Kapton). Several runs using different combinations of bias voltages and secondary electron yields are discussed. Special attention is paid to the role of the secondaries as the dielectric surface changes from a negative potential to become increasingly positive. This is termed the 'snapover' process. The contribution due to secondaries is examined by plotting the distribution (number versus position) of the particles emitted by the impacting on the dielectric. These plots show that secondaries are drawn back into the dielectric as well as the conductor in the case of snapover. Mechanisms and experimental implications for these particular distributions are discussed.

Kessel, R. L.↗
NASA NTRS1985

Numerical simulation of transient inviscid gas flows in a shock tube

Time-dependent upwind high resolution schemes for solving the Euler equations were developed and applied to simulate 1-D and 2-D transient inviscid gas flows in a shock tube. Using obstacles of different geometries, a series of calculations were carried out to investigate the transient complex shock-wave diffraction phenomena. Shock-wave and objects interactions with shock Mach-number Ms ranging from 2 to 20 were considered. Comparison with analytical and available experimental results indicate good agreement. Display of detailed flow structures including multiple Mach shocks, slip surfaces, and vortex are also given.

Yang, J. Y.↗
NASA NTRS1985

Numerical simulation of viscous supersonic flow over a generic fighter configuration

A procedure is presented, as well as some results, to calculate the flow over a generic fighter configuration. A parabolized marching Navier-Stokes code is used to obtain the solution over a wing-canopy body. The flow conditions simulate supersonic cruise with a freestream Mach number of 2.169 and angles of attack of 4 and 10 deg. The body surface is considered to be an adiabatic wall and the flow is assumed to be turbulent for the given Reynolds number.

Chaussee, D. S.↗
NASA NTRS1985

Numerical simulation of El Chichon volcanic cloud

The stratospheric volcanic cloud from the eruption of El Chichon, Mexico, on April 4, 1982 was observed routinely by a Nd:YAG lidar system from April 18, 1982 at Kyushu University. There were two phenomena detected by lidar observation of the El Chichon volcanic cloud: (1) The main enhanced layer was formed in the summer easterly wind region over Fukuoka (33.5 deg N, 130.4 deg E). The layer had two high gradient top and bottom boundaries. The layer width slowly broadened until September 1982, when the easterly wind changed to the westerly wind. A similar phenomenon was observed after the eruption of Mt. St. Helens in May 1980. (2) The stratospheric integrated backscattering coefficient (IBC) reached a peak value on May 3, 1982. It then gradually decreased until August, but reincreased significantly from September to December 1982. After January 1983, IBC decreased slowly. This remarkable reincrease had not been seen in past increases which had been observed by lidar at Kyushu University. These phenomena are briefly discussed.

Shibata, T.↗
NASA NTRS1985

Results from a limited area mesoscale numerical simulation for 10 April 1979

Results are presented from a nine-hour limited area fine mesh (35-km) mesoscale model simulation initialized with SESAME-AVE I radiosonde data for Apr. 10, 1979 at 2100 GMT. Emphasis is on the diagnosis of mesoscale structure in the mass and precipitation fields. Along the Texas/Oklahoma border, independent of the short wave, convective precipitation formed several hours into the simulation and was organized into a narrow band suggestive of the observed April 10 squall line.

Kalb, M. W.↗
NASA NTRS1986

Numerical simulation of backward step and jet exhaust flows

Preparatory to a continuing study of efficient, high resolution flow structure capture on patched mesh systems, turbulent flowfield solutions for a tactical missile base flow problem are computed on two first level composite grids of differing topologies. One mesh is of the wraparound character and the other is a step mesh on the Cartesian topology. A simplified Baldwin-Lomax algebraic eddy viscosity model is applied in the wake region. Results for the two meshes are compared with each other and with other computations reported in the literature. The totality of results suggest both mesh topology and resolution are important factors in the accuracy of predictive results as well as turbulence modeling. Related to interests in high altitude laminar flows bluff bodies such as the AOTV, the paper also presents results from a solution from laminar flow over a backward facing step and makes comparison with the available classical experiment of Chapman, et al.

Lombard, C. K.↗