Search NASA⌕ Search

SEARCH · Search NASA

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 235 records · Page 13

Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer

The results from direct numerical simulations of a spatially evolving, supersonic, flat-plate turbulent boundary-layer flow, with free-stream Mach number of 2.25 are presented. The simulated flow field extends from a transition region, initiated by wall suction and blowing near the inflow boundary, into the fully turbulent regime. Distributions of mean and turbulent flow quantities are obtained and an analysis of these quantities is performed at a downstream station corresponding to Re(sub x)= 5.548 x10(exp 6) based on distance from the leading edge.

Gatski, T. B.↗

Numerical Simulation of the Water Cycle Change Over the 20th Century

We have used numerical models to test the impact of the change in Sea Surface Temperatures (SSTs) and carbon dioxide (CO2) concentration on the global circulation, particularly focusing on the hydrologic cycle, namely the global cycling of water and continental recycling of water. We have run four numerical simulations using mean annual SST from the early part of the 20th century (1900-1920) and the later part (1980-2000). In addition, we vary the CO2 concentrations for these periods as well. The duration of the simulations is 15 years, and the spatial resolution is 2 degrees. We use passive tracers to study the geographical sources of water. Surface evaporation from predetermined continental and oceanic regions provides the source of water for each passive tracer. In this way, we compute the percent of precipitation of each region over the globe. This can also be used to estimate precipitation recycling. In addition, we are using the passive tracers to independently compute the global cycling of water (compared to the traditional, Q/P calculation).

Bosilovich, Michael G.↗

Direct numerical simulation of incompressible axisymmetric flows

In the present work, we propose to conduct direct numerical simulations (DNS) of incompressible turbulent axisymmetric jets and wakes. The objectives of the study are to understand the fundamental behavior of axisymmetric jets and wakes, which are perhaps the most technologically relevant free shear flows (e.g. combuster injectors, propulsion jet). Among the data to be generated are various statistical quantities of importance in turbulence modeling, like the mean velocity, turbulent stresses, and all the terms in the Reynolds-stress balance equations. In addition, we will be interested in the evolution of large-scale structures that are common in free shear flow. The axisymmetric jet or wake is also a good problem in which to try the newly developed b-spline numerical method. Using b-splines as interpolating functions in the non-periodic direction offers many advantages. B-splines have local support, which leads to sparse matrices that can be efficiently stored and solved. Also, they offer spectral-like accuracy that are C(exp O-1) continuous, where O is the order of the spline used; this means that derivatives of the velocity such as the vorticity are smoothly and accurately represented. For purposes of validation against existing results, the present code will also be able to simulate internal flows (ones that require a no-slip boundary condition). Implementation of no-slip boundary condition is trivial in the context of the b-splines.

Loulou, Patrick↗

Numerical simulations of the nonlinear kink modes in linearly stable supersonic slip surfaces

High-resolution numerical simulations of supersonic slip surfaces were performed using the Cray 2 at the Minnesota Supercomputer Center to confirm earlier analytic nonlinear stability calculations of such structures. The present study presents a numerical solution for the problem solved analytically by Artola and Majda (1987) and follows growing modes well into the large-amplitude regime. The response of an equal-density Mach-4 slip surface to a variety of impinging nonlinear sound wave trains is computed using the piecewise-parabolic method (PPM). The nature of the nonlinear kink modes observed in the simulations was found to be similar to that discussed by Artola and Majda. Most of the speeds moved in either direction, with speeds near the predicted value. It is suggested that the stationary modes will eventually dominate the flow at much later times, and it is argued that the stationary modes are more disruptive than the propagating kink modes.

Pedelty, Jeffrey A.↗

Efficient numerical simulation of electron states in quantum wires

A new algorithm is presented for the numerical simulation of electrons in a quantum wire as described by a two-dimensional eigenvalue problem for Schroedinger's equation coupled with Poisson's equation. Initially, the algorithm employs an underrelaxed fixed point iteration to generate an approximation which is reasonably close to the solution. Subsequently, this approximate solution is employed as an initial guess for a Jacobian-free implementation of an approximate Newton method. In this manner the nonlinearity in the model is dealt with effectively. The effectiveness of this approach is demonstrated in a set of numerical experiments which study the electron states on the cross section of a quantum wire structure based on III-V semiconductors at 4.2 and 77 K.

Kerkhoven, Thomas↗

Determination of thermospheric quantities from simple ionospheric observations using numerical simulation

Measured ionospheric electron content and peak electron concentration data are introduced into a numerical simulation of the ionosphere to yield values of induced plasma drifts and exospheric neutral temperatures consistent with the observations. Data collected on 23-24 March 1970 on the East Coast of the U.S.A. are analyzed and the results are in agreement with incoherent radar measurements at Millstone Hill, Massachusetts. Neutral winds and meridional exospheric temperature gradients that give rise to the computed plasma drifts are calculated through the use of a dynamic model of the thermosphere.

Antoniadis, D. A.↗

Analysis of Numerical Simulation Database for Pressure Fluctuations Induced by High-Speed Turbulent Boundary Layers

Direct numerical simulations (DNS) of Mach 6 turbulent boundary layer with nominal freestream Mach number of 6 and Reynolds number of Re(sub T) approximately 460 are conducted at two wall temperatures (Tw/Tr = 0.25, 0.76) to investigate the generated pressure fluctuations and their dependence on wall temperature. Simulations indicate that the influence of wall temperature on pressure fluctuations is largely limited to the near-wall region, with the characteristics of wall-pressure fluctuations showing a strong temperature dependence. Wall temperature has little influence on the propagation speed of the freestream pressure signal. The freestream radiation intensity compares well between wall-temperature cases when normalized by the local wall shear; the propagation speed of the freestream pressure signal and the orientation of the radiation wave front show little dependence on the wall temperature.

Duan, Lian↗

Numerical simulation of a disk-shaped electron accelerating electrostatic probe

This paper presents a method utilizing the numerical plasma simulation particle-in-cell technique for studying electron accelerating probe characteristics. A probe configuration other than the commonly used plane, cylindrical, or spherical case will be studied; the probe configuration used is that of a disk. Cylindrical geometry is used in the simulation model. One of the advantages of a numerical simulation is the ease with which the details of the system can be examined. Therefore, in addition to the current-voltage characteristics of the static probe, other results are presented such as the time evolution of the plasma properties and studies of the sheath region surrounding the probe.

Nonnast, J. H.↗

A One Dimensional, Time Dependent Inlet/Engine Numerical Simulation for Aircraft Propulsion Systems

The NASA Lewis Research Center (LeRC) and the Arnold Engineering Development Center (AEDC) have developed a closely coupled computer simulation system that provides a one dimensional, high frequency inlet/engine numerical simulation for aircraft propulsion systems. The simulation system, operating under the LeRC-developed Application Portable Parallel Library (APPL), closely coupled a supersonic inlet with a gas turbine engine. The supersonic inlet was modeled using the Large Perturbation Inlet (LAPIN) computer code, and the gas turbine engine was modeled using the Aerodynamic Turbine Engine Code (ATEC). Both LAPIN and ATEC provide a one dimensional, compressible, time dependent flow solution by solving the one dimensional Euler equations for the conservation of mass, momentum, and energy. Source terms are used to model features such as bleed flows, turbomachinery component characteristics, and inlet subsonic spillage while unstarted. High frequency events, such as compressor surge and inlet unstart, can be simulated with a high degree of fidelity. The simulation system was exercised using a supersonic inlet with sixty percent of the supersonic area contraction occurring internally, and a GE J85-13 turbojet engine.

Doug Garrard↗

Direct Numerical Simulation of Smooth-Body Flow Separation Around a Ramp

Spanwise-periodic computation of a turbulent flow past a two-dimensional smooth ramp geometry is performed in the form of a direct numerical simulation. The Reynolds number based on the ramp height is about 147,000. A straight section that precedes the smooth ramp allows the incoming turbulent boundary layer to grow under a weak favorable pressure gradient. The boundary layer introduced at the domain inlet has a momentum-thickness based Reynolds number of 2000. The turbulent boundary layer nearing the ramp first interacts with a relatively stronger favorable pressure gradient, followed by a strong adverse pressure gradient. Consequently, the boundary layer experiences a modest acceleration before decelerating and separating. Analysis of the data over this region hints at the formation of an internal layer beneath the accelerated boundary layer. The analysis also reveals that this internal layer forms the origin of the free shear layer that emerges in the deceleration region and separates. The streamwise extent of the separated region is comparable to the ramp length, while the viscous layer thickness near reattachment is about the same as the ramp height; hence, the boundary layer undergoing separation and subsequent reattachment in the present configuration experiences its thickness being amplified by about tenfold. The reattached flow continues to develop further under a diminishing pressure gradient in the recovery region in a similar fashion to a zero pressure gradient turbulent boundary layer.

Turbulence↗

Comparison of Nonlinear Random Response Using Equivalent Linearization and Numerical Simulation

A recently developed finite-element-based equivalent linearization approach for the analysis of random vibrations of geometrically nonlinear multiple degree-of-freedom structures is validated. The validation is based on comparisons with results from a finite element based numerical simulation analysis using a numerical integration technique in physical coordinates. In particular, results for the case of a clamped-clamped beam are considered for an extensive load range to establish the limits of validity of the equivalent linearization approach.

Rizzi, Stephen A.↗

Direct numerical simulation of transition: The spatial approach

The objective of this paper is to review the efforts in spatial direct numerical simulations for transition modeling. Much recent success has been realized in the development of more efficient numerical algorithms as well as a robust downstream boundary condition. Efforts at explaining complex physical phenomena through the use of simulations are reviewed.

Reed, Helen L.↗

Numerical simulation of the flow field around a complete aircraft

The present effort represents a first attempt of numerical simulation of the flow field around a complete aircraft-like, lifting configuration utilizing the Reynolds averaged Navier-Stokes equations. The numerical solution generated for the experimental aircraft concept X24C-10D at a Mach number of 5.95 not only exhibited accurate prediction of detailed flow properties but also of the integrated aerodynamic coefficients. In addition, the present analysis demonstrated that a page structure of data collected into cyclic blocks is an efficient and viable means for processing the Navier-Stokes equations on the CRAY XMP-22 computer with external memory device.

Shang, J. S.↗

Numerical simulation of a conducting disk partially covered by an insulator interacting with a plasma

A numerical simulation has been made of the 'pinhole effect' which produces the characteristic 'S-shaped' current-voltage curve. A disk-shaped conducting probe immersed in a plasma is modeled using a particle-in-cell (PIC) code. A probe partially covered by a very thin insulating layer is considered, as well as a probe mounted on an insulating disk. The simulation uses a cylindrical particle mover and allows for a variable number of particles in the system. The simulation space grid uses three different mesh sizes, the coarsest being away from the probe and the finest near the probe, in order to accurately calculate the trajectories of the simulation particles contributing the current to the probe and to the surface charge density on the dielectric. The calculation of the electrostatic potential is done self-consistently using successive over-relaxation (SOR). Backscattering and secondary electron emission are included for the case of positive probe voltage.

Chaky, R. C.↗

Numerical simulations of adiabatic axisymmetric accretion flow. I - A new mechanism for the formation of jets

Numerical simulations of the uniform axisymmetric flow past a gravitating sphere have been studied. It is found that the structure of the flow is extremely sensitive to the boundary condition at the surface of the gravitating object. For the case in which the boundary is totally absorbing, a steady state flow is reached. However, for a boundary which is not totally absorbing, steady state flows are not obtained. The morphology of the flow is also sensitive to the Mach number at infinity and to the ratio of the free-fall velocity at the surface of the gravitating object to the flow velocity at inifinity. A new mechanism for the formation of jets is identified in which a fraction of the accretion energy is tapped to drive an anisotropic supersonic outflow with collimation provided by a combination of the inertia of matter which surrounds the beam and the development of multiple shock structures.

Fryxell, B. A.↗

Rocket engine numerical simulator

The topics are presented in viewgraph form and include the following: a rocket engine numerical simulator (RENS) definition; objectives; justification; approach; potential applications; potential users; RENS work flowchart; RENS prototype; and conclusion.

Davidian, Ken↗

Towards the Identification of the Keeper Erosion Cause(s): Numerical Simulations of the Plasma and Neutral Gas Using the Global Cathode Model OrCa2D-II

Numerical simulations with the time-dependent Orificed Cathode (OrCa2D-II) computer code show that classical enhancements of the plasma resistivity can not account for the elevated electron temperatures and steep plasma potential gradients measured in the plume of a 25-27.5 A discharge hollow cathode. The cathode, which employs a 0.11-in diameter orifice, was operated at 5.5 sccm without an applied magnetic field using two different anode geometries. It is found that anomalous resistivity based on electron-driven instabilities improves the comparison between theory and experiment. It is also estimated that other effects such as the Hall-effect from the self-induced magnetic field, not presently included in OrCa2D-II, may contribute to the constriction of the current density streamlines thus explaining the higher plasma densities observed along the centerline.

plasma modeling↗