Search NASA⌕ Search

SEARCH · Search NASA

Results for “particle tracing”

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 91 records · Page 5

Navier-Stokes analysis of an oxidizer turbine blade with tip clearance with and without a mini-shroud

The Gas Generator Oxidizer Turbine (GGOT) Blade is being analyzed by various investigators under the NASA MSFC-sponsored Turbine Stage Technology Team design effort. The present work concentrates on the tip clearance region flow and associated losses; however, flow details for the passage region are also obtained in the simulations. The present calculations simulate the rotor blade row in a rotating reference frame with the appropriate coriolis and centrifugal acceleration term included in the momentum equations. The upstream computational boundary is located about one axial chord from the blade leading edge. The boundary conditions at this location have been determined by Pratt & Whitney using an Euler analysis without the vanes to obtain approximately the same flow profiles at the rotor as were obtained with the Euler stage analysis including the vanes. Inflow boundary layer profiles are then constructed assuming the skin friction coefficient at both the hub and the casing. The downstream computational boundary is located about one axial chord from the blade trailing edge, and the circumferentially averaged static pressure at this location was also obtained from the P&W Euler analysis. Results obtained for the 3-D baseline GGOT geometry at the full scale design Reynolds number show a region of high loss in the region near the casing. Particle traces in the near tip region show vortical flow behavior of the fluid which passes through the clearance region and exits at the downstream edge of the gap. In an effort to reduce clearance flow losses, the mini-shroud concept was proposed by the Pratt & Whitney design team. Calculations were performed on the GGO geometry with the mini-shroud. Results of these calculations indicate that the mini-shroud does not significantly affect the flow in the passage region, and although the tip clearance flow is different, the mini-shroud does not seem to prevent the above-mentioned vortical flow behavior. Since both flow distortion and total pressure losses are similar for both geometries, the addition of the mini-shroud does not seem to reduce the tip clearance flow effects.

Chan, Tony↗

An Interactive Tool for Discrete Phase Analysis in Two-Phase Flows

Under a NASA MSFC SBIR Phase 1 effort an interactive software package has been developed for the analysis of discrete (particulate) phase dynamics in two-phase flows in which the discrete phase does not significantly affect the continuous phase. This package contains a Graphical User Interface (based on the X Window system and the Motif tool kit) coupled to a particle tracing program, which allows the user to interactively set up and run a case for which a continuous phase grid and flow field are available. The software has been applied to a solid rocket motor problem, to demonstrate its ease of use and its suitability for problems of engineering interest, and has been delivered to NASA Marshall Space Flight Center.

Dejong, Frederik J.↗

Four-Nozzle Benchmark Wind Tunnel Model USA Code Solutions for Simulation of Multiple Rocket Base Flow Recirculation at 145,000 Feet Altitude

Multiple rocket exhaust plume interactions at high altitudes can produce base flow recirculation with attendant alteration of the base pressure coefficient and increased base heating. A search for a good wind tunnel benchmark problem to check grid clustering technique and turbulence modeling turned up the experiment done at AEDC in 1961 by Goethert and Matz on a 4.25-in. diameter domed missile base model with four rocket nozzles. This wind tunnel model with varied external bleed air flow for the base flow wake produced measured p/p(sub ref) at the center of the base as high as 3.3 due to plume flow recirculation back onto the base. At that time in 1961, relatively inexpensive experimentation with air at gamma = 1.4 and nozzle A(sub e)/A of 10.6 and theta(sub n) = 7.55 deg with P(sub c) = 155 psia simulated a LO2/LH2 rocket exhaust plume with gamma = 1.20, A(sub e)/A of 78 and P(sub c) about 1,000 psia. An array of base pressure taps on the aft dome gave a clear measurement of the plume recirculation effects at p(infinity) = 4.76 psfa corresponding to 145,000 ft altitude. Our CFD computations of the flow field with direct comparison of computed-versus-measured base pressure distribution (across the dome) provide detailed information on velocities and particle traces as well eddy viscosity in the base and nozzle region. The solution was obtained using a six-zone mesh with 284,000 grid points for one quadrant taking advantage of symmetry. Results are compared using a zero-equation algebraic and a one-equation pointwise R(sub t) turbulence model (work in progress). Good agreement with the experimental pressure data was obtained with both; and this benchmark showed the importance of: (1) proper grid clustering and (2) proper choice of turbulence modeling for rocket plume problems/recirculation at high altitude.

Dougherty, N. S.↗

Distributed interactive graphics applications in computational fluid dynamics

Implementation of two interactive, distributed graphics programs used in Computational Fluid Dynamics is discussed. Both programs run on a Cray 2 supercomputer and use a Silicon Graphics Iris workstation as the graphics front-end machine. The hardware and supporting software is from the Numerical Aerodynamic Simulation project. Using this configuration, the supercomputer does all of the numerically intensive work and the workstation allows the user to perform real-time interactive transformations on the displayed data. The first program was written originally as a distributed program which computes particle traces for fluid flow solutions existing on the supercomputer. The second is an older post-processing and plotting program which was modified to run in a distributed mode. Both programs have realized a large increase in capability as a distributed process. Some graphical results are presented.

Rogers, Stuart E.↗

Effect of B(sub y) on neutral line ridges and dynamical source ordering

In this paper we study the effect of a uniform magnetic field B(sub y) in the cross-tail y direction on the dynamics and distribution of energetic ions in a current sheet model including a neutral line. Martin and Speiser (1988) have previously shown that a 'ridge' in the velocity space distribution function is a remote neutral line signature in a two-dimensional field without B(sub y). Our results show that as B(sub y) is increased, using nominal tail parameters, there is very little change in the ridge signature for small B(sub y) (up to about 4 times B(sub z)). For intermediate values (up to about 10 times times B(sub Y)) the ridge becomes observable further from the neutral line, while close to the X line the ridge is strongly modified. For large B(sub y) (of the order of B(sub x)) the ridge evolves into the new neutral line signatures which depend strongly on whether the observations point is above or below the current sheet. These results are used to estimate B(sub y) less than 2 nT in the Active Magentospheric Particle Traces Explorers/Ion Release Module (AMPTE/IRM) event modeled by Speiser and Martin (1994), which is consistent with onboard magnetometer measurements. We further show that the cause of the new structures appearing at large B(sub Y) is the same source seperation effect that produced the original ridge and that all these signatures are relatively insensitive to a small parallel electric field. Finally, we point our a dynamical ordering which occurs at large B(sub y): Using high-resolution simulations, we find a common boundary separating regions of positive and negative initial x, y, and z position, as well as initial pitch angle. This boundary is purely dynamical in nature and is independent of modeled source distributions.

Martin, R. F., Jr.↗

Unsteady Simulation of the Viscous Flow About a V-22 Rotor and Wing in Hover

Results of an unsteady thin-layer Navier-Stokes simulation of a 0.658-scale V-22 rotor and wing configuration in hover are presented. All geometric components of the flapped-wing and rotor test rig, including rotor blades, are accurately modeled. Rotor motion and rotor/airframe interference effects are simulated directly using moving body overset grid methods. Tiltrotor hover aerodynamics are visualized via unsteady particle trace images. Wing download predictive ability is demonstrated. Simulation results are compared with experimental data.

Meakin, Robert L.↗

Subsonic Analysis of 0.04-Scale F-16XL Models Using an Unstructured Euler Code

The subsonic flow field about an F-16XL airplane model configuration was investigated with an inviscid unstructured grid technique. The computed surface pressures were compared to wind-tunnel test results at Mach 0.148 for a range of angles of attack from 0 deg to 20 deg. To evaluate the effect of grid dependency on the solution, a grid study was performed in which fine, medium, and coarse grid meshes were generated. The off-surface vortical flow field was locally adapted and showed improved correlation to the wind-tunnel data when compared to the nonadapted flow field. Computational results are also compared to experimental five-hole pressure probe data. A detailed analysis of the off-body computed pressure contours, velocity vectors, and particle traces are presented and discussed.

Lessard, Wendy B.↗

Aeropropulsion Technology (APT). Task 23 - Stator Seal Cavity Flow Investigation

The focus of NASA Contract NAS3-25950 Task 23 was to numerically investigate the flow through an axial compressor inner-banded stator seal cavity. The Allison/NASA developed ADPAC code was used to obtain all flow predictions. Flow through a labyrinth stator seal cavity of a high-speed compressor was modeled by coupling the cavity flow path and the main flow path of the compressor. A grid resolution study was performed to guarantee adequate grid spacing was used. Both unsteady rotor-stator-rotor interactions and steady-state isolated blade calculations were performed with and without the seal cavity present. A parameterized seal cavity study of the high-speed stator seal cavity collected a series of solutions for geometric variations. The parameter list included seal tooth gap, cavity depth, wheel speed, radial mismatch of hub flowpath, axial trench gap, hub corner treatments, and land edge treatments. Solution data presented includes radial and pitchwise distributions of flow variables and particle traces describing the flow character.

Heidegger, N. J.↗

Visualization of Unsteady Computational Fluid Dynamics

The current compute environment that most researchers are using for the calculation of 3D unsteady Computational Fluid Dynamic (CFD) results is a super-computer class machine. The Massively Parallel Processors (MPP's) such as the 160 node IBM SP2 at NAS and clusters of workstations acting as a single MPP (like NAS's SGI Power-Challenge array and the J90 cluster) provide the required computation bandwidth for CFD calculations of transient problems. If we follow the traditional computational analysis steps for CFD (and we wish to construct an interactive visualizer) we need to be aware of the following: (1) Disk space requirements. A single snap-shot must contain at least the values (primitive variables) stored at the appropriate locations within the mesh. For most simple 3D Euler solvers that means 5 floating point words. Navier-Stokes solutions with turbulence models may contain 7 state-variables. (2) Disk speed vs. Computational speeds. The time required to read the complete solution of a saved time frame from disk is now longer than the compute time for a set number of iterations from an explicit solver. Depending, on the hardware and solver an iteration of an implicit code may also take less time than reading the solution from disk. If one examines the performance improvements in the last decade or two, it is easy to see that depending on disk performance (vs. CPU improvement) may not be the best method for enhancing interactivity. (3) Cluster and Parallel Machine I/O problems. Disk access time is much worse within current parallel machines and cluster of workstations that are acting in concert to solve a single problem. In this case we are not trying to read the volume of data, but are running the solver and the solver outputs the solution. These traditional network interfaces must be used for the file system. (4) Numerics of particle traces. Most visualization tools can work upon a single snap shot of the data but some visualization tools for transient problems require dealing with time.

Haimes, Robert↗

Effects of low frequency waves and spiky electric fields in the magnetotail

Intense low frequency waves and large amplitude spiky electric fields are commonly observed in the active magnetotail. These fields are of dynamic significance for substorms because the amplitude of the waves is large enough to provide the dissipation necessary for reconnection. The waves modify the trajectories of ions in the magnetotail, resulting in enhanced energization and pitch angle scattering in comparison with the trajectories obtained without waves. The spikes may represent the nonlinear evolution of the waves. Observations are presented of low frequency waves and spiky electric fields measured in the plasma sheet. The relationship between the waves and the spikes is compared to wave evolution theories. The Lundquist number is evaluated for the magnetotail based on the wave observations and previously published studies of the current sheet scale sizes. The effects of the waves on the ion trajectories are discussed using the results from a particle tracing code.

Cattell, C.↗

Non-Linear Cosmological Power Spectra in Real and Redshift Space

We present an expression for the non-linear evolution of the cosmological power spectrum based on Lagrangian trajectories. This is simplified using the Zel'dovich approximation to trace particle displacements, assuming Gaussian initial conditions. The model is found to exhibit the transfer of power from large to small scales expected in self-gravitating fields. Some exact solutions are found for power-law initial spectra. We have extended this analysis into red-shift space and found a solution for the non-linear, anisotropic redshift-space power spectrum in the limit of plane-parallel redshift distortions. The quadrupole-to-monopole ratio is calculated for the case of power-law initial spectra. We find that the shape of this ratio depends on the shape of the initial spectrum, but when scaled to linear theory depends only weakly on the redshift-space distortion parameter, beta. The point of zero-crossing of the quadrupole, kappa(sub o), is found to obey a simple scaling relation and we calculate this scale in the Zel'dovich approximation. This model is found to be in good agreement with a series of N-body simulations on scales down to the zero-crossing of the quadrupole, although the wavenumber at zero-crossing is underestimated. These results are applied to the quadrupole-to-monopole ratio found in the merged QDOT plus 1.2-Jy-IRAS redshift survey. Using a likelihood technique we have estimated that the distortion parameter is constrained to be beta greater than 0.5 at the 95 percent level. Our results are fairly insensitive to the local primordial spectral slope, but the likelihood analysis suggests n = -2 un the translinear regime. The zero-crossing scale of the quadrupole is k(sub 0) = 0.5 +/- 0.1 h Mpc(exp -1) and from this we infer that the amplitude of clustering is sigma(sub 8) = 0.7 +/- 0.05. We suggest that the success of this model is due to non-linear redshift-space effects arising from infall on to caustic and is not dominated by virialized cluster cores. The latter should start to dominate on scales below the zero-crossing of the quadrupole, where our model breaks down.

Taylor, A. N.↗

Out-of-Core Streamline Visualization on Large Unstructured Meshes

It's advantageous for computational scientists to have the capability to perform interactive visualization on their desktop workstations. For data on large unstructured meshes, this capability is not generally available. In particular, particle tracing on unstructured grids can result in a high percentage of non-contiguous memory accesses and therefore may perform very poorly with virtual memory paging schemes. The alternative of visualizing a lower resolution of the data degrades the original high-resolution calculations. This paper presents an out-of-core approach for interactive streamline construction on large unstructured tetrahedral meshes containing millions of elements. The out-of-core algorithm uses an octree to partition and restructure the raw data into subsets stored into disk files for fast data retrieval. A memory management policy tailored to the streamline calculations is used such that during the streamline construction only a very small amount of data are brought into the main memory on demand. By carefully scheduling computation and data fetching, the overhead of reading data from the disk is significantly reduced and good memory performance results. This out-of-core algorithm makes possible interactive streamline visualization of large unstructured-grid data sets on a single mid-range workstation with relatively low main-memory capacity: 5-20 megabytes. Our test results also show that this approach is much more efficient than relying on virtual memory and operating system's paging algorithms.

Ueng, Shyh-Kuang↗

Parallel Visualization Co-Processing of Overnight CFD Propulsion Applications

An interactive visualization system pV3 is being developed for the investigation of advanced computational methodologies employing visualization and parallel processing for the extraction of information contained in large-scale transient engineering simulations. Visual techniques for extracting information from the data in terms of cutting planes, iso-surfaces, particle tracing and vector fields are included in this system. This paper discusses improvements to the pV3 system developed under NASA's Affordable High Performance Computing project.

Edwards, David E.↗

CIIs During the May 98 CMEs

The Thermal Ion Dynamics Experiment (TIDE) instrument on the ISTP/Polar spacecraft observed intense Cusp Ion Injections (CIIs) during the May 1998 series of Coronal Mass Ejections (CMEs). The energy signatures of the CIIs can be used to determine the magnetic field topology of the geomagnetic field and the Interplanetary Magnetic Field (IMF) at the reconnection site. Particle tracing using model magnetic fields and the TIDE on observations are used to locate the point of magnetopause crossing.

Chandler, Michael O.↗

Computational Analysis for Rocket-Based Combined-Cycle Systems During Rocket-Only Operation

A series of Reynolds-averaged Navier-Stokes calculations were employed to study the performance of rocket-based combined-cycle systems operating in an all-rocket mode. This parametric series of calculations were executed within a statistical framework, commonly known as design of experiments. The parametric design space included four geometric and two flowfield variables set at three levels each, for a total of 729 possible combinations. A D-optimal design strategy was selected. It required that only 36 separate computational fluid dynamics (CFD) solutions be performed to develop a full response surface model, which quantified the linear, bilinear, and curvilinear effects of the six experimental variables. The axisymmetric, Reynolds-averaged Navier-Stokes simulations were executed with the NPARC v3.0 code. The response used in the statistical analysis was created from Isp efficiency data integrated from the 36 CFD simulations. The influence of turbulence modeling was analyzed by using both one- and two-equation models. Careful attention was also given to quantify the influence of mesh dependence, iterative convergence, and artificial viscosity upon the resulting statistical model. Thirteen statistically significant effects were observed to have an influence on rocket-based combined-cycle nozzle performance. It was apparent that the free-expansion process, directly downstream of the rocket nozzle, can influence the Isp efficiency. Numerical schlieren images and particle traces have been used to further understand the physical phenomena behind several of the statistically significant results.

Steffen, C. J., Jr.↗

Transport of Energetic Ions in the Ring Current During Geomagnetic Storms

In the final year (plus no-cost extentions) of this grant, we have: Used the particle tracing code to perform a systematic study of the expected energy spectra over the full range of local times in the ring current using a variety of electric and magnetic field models. Shown that the Weimer electric field is superior to the Volland-Stern electric field in reproducing the observed energy spectra on the AMPTE CCE spacecraft. Redone our analysis of the pitch angle spectra of energetic ions during storms in the magnetosphere, using a larger data set, and a more reliable classification technique.

Kistler, Lynn M.↗

Scientific Visualization of Large Scale Unsteady Fluid Flow

In a numerical flow simulation. it is common to generate several thousand time steps of unsteady (time-dependent) flow data. Each time step may require tens to hundreds of megabytes for disk storage, and the total disk requirement for storing the unsteady flow data may be hundreds of gigabytes. Interactive visualization of unsteady flow data of this magnitude is presently impossible with the current hardware technology. This chapter describes the current approaches for unsteady flow visualization. An effective particle tracing technique for unsteady flow is also described. First, the life cycle of a typical numerical flow simulation is outlined. Several unsteady flow data sets from real-world problems are then given. The current approaches for visualizing unsteady flow are then described. There are many existing systems for flow visualization, and some of them are discussed. Streaklines depict time-varying phenomena that are sometimes difficult or impossible to see with other visualization techniques. The algorithms for computing streaklines are described. Several unsteady flow data sets have been visualized using streaklines, and the results are presented. Finally, some current issues in unsteady flow visualization are discussed.

Lane, David A.↗

Computational Modeling of Vortex Generators for Turbomachinery

In this work computational models were developed and used to investigate applications of vortex generators (VGs) to turbomachinery. The work was aimed at increasing the efficiency of compressor components designed for the NASA Ultra Efficient Engine Technology (UEET) program. Initial calculations were used to investigate the physical behavior of VGs. A parametric study of the effects of VG height was done using 3-D calculations of isolated VGs. A body force model was developed to simulate the effects of VGs without requiring complicated grids. The model was calibrated using 2-D calculations of the VG vanes and was validated using the 3-D results. Then three applications of VGs to a compressor rotor and stator were investigated: 1) The results of the 3-D calculations were used to simulate the use of small casing VGs used to generate rotor preswirl or counterswirl. Computed performance maps were used to evaluate the effects of VGs. 2) The body force model was used to simulate large part-span splitters on the casing ahead of the stator. Computed loss buckets showed the effects of the VGs. 3) The body force model was also used to investigate the use of tiny VGs on the stator suction surface for controlling secondary flows. Near-surface particle traces and exit loss profiles were used to evaluate the effects of the VGs.

Chima, R. V.↗