SEARCH · Engineering Papers

Results for “Numerical”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,657 records · Page 92

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Time domain numerical calculations of unsteady vortical flows about a flat plate airfoil

A time domain numerical scheme is developed to solve for the unsteady flow about a flat plate airfoil due to imposed upstream, small amplitude, transverse velocity perturbations. The governing equation for the resulting unsteady potential is a homogeneous, constant coefficient, convective wave equation. Accurate solution of the problem requires the development of approximate boundary conditions which correctly model the physics of the unsteady flow in the far field. A uniformly valid far field boundary condition is developed, and numerical results are presented using this condition. The stability of the scheme is discussed, and the stability restriction for the scheme is established as a function of the Mach number. Finally, comparisons are made with the frequency domain calculation by Scott and Atassi, and the relative strengths and weaknesses of each approach are assessed.

Hariharan, S. I.↗

NASA NTRS1992

A numerical assessment of rough surface scattering theories. I - Horizontal polarization. II - Vertical polarization

A numerical evaluation is presented of the regime of validity for various rough surface scattering theories against numerical results obtained by employing the method of moments. The contribution of each theory is considered up to second order in the perturbation expansion for the surface current. Considering both vertical and horizontal polarizations, the unified perturbation method provides best results among all theories weighed.

Rodriguez, Ernesto↗

NASA NTRS1992

A finite-volume numerical method to calculate fluid forces and rotordynamic coefficients in seals

A numerical method to calculate rotordynamic coefficients of seals is presented. The flow in a seal is solved by using a finite-volume formulation of the full Navier-Stokes equations with appropriate turbulence models. The seal rotor is perturbed along a diameter such that the position of the rotor is a sinusoidal function of time. The resulting flow domain changes with time, and the time-dependent flow in the seal is solved using a space conserving moving grid formulation. The time-varying fluid pressure reaction forces are then linked with the rotor center displacement, velocity and acceleration to yield the rotordynamic coefficients. Results for an annular seal are presented, and compared with experimental data and other more simplified numerical methods.

Athavale, M. M.↗

NASA NTRS1992

Analytic and numerical models of wave-CISK with conditional heating

Wave-CISK with conditional heating is investigated in the equatorial zonal-height plane by analytic and numerical techniques. For two- and three-level models, previous results are extended to give additional evidence that the most unstable mode has a single wet region of infinitesimal width. A three-level model has qualitatively similar behavior as the two-level model except that propagating solutions are possible due to coalescence of internal vertical modes. Phase speeds with conditional heating are found to be slightly greater than those for unconditional heating. The structure has one circulation cell in the vertical and is asymmetric in longitude with stronger motion on the leading edge. Growth rate is inversely proportional to the width of the single wet region. That width can be limited by second-order diffusion. A general integral relationship between growth rate, viscosity, phase speed, and heating is derived. The main conclusion is that the linear wave-CISK catastrophe is modified by conditional heating but not eliminated. The preferred mode of instability has one wet region, but it occurs on the smallest possible scale. It is likely that numerical models that use conditional heating are sensitive to resolution, especially for the commonly used spectral truncations, unless there is sufficiently strong damping at the smallest scales.

Crum, Francis X.↗

NASA NTRS1992

Thermal stability during rotation in space - A scaling and numerical analysis

An analysis of thermal stability in a fluid during rotation is presented. The conditions resulting in thermal stability are examined and numerically verified. It is shown that at least one thermal field exists for any rotating fluid, with or without a gravitational background field, where convection does not occur. The numerical model used is fully nonlinear and 3D. Results are presented on the basis of both steady-state and fully transient models. The thermally stable condition is further investigated in the context of the average net acceleration and the inhomogeneous part of the acceleration field (gradient of acceleration). It is shown that in the thermally stable condition the flow that would result from the average net acceleration is equal and opposite to the flow that would result from the inhomogeneous part of the acceleration field. It is concluded that the thermally stable condition is actually a flow mode transition at Ad = 1.

Arnold, William A.↗

NASA NTRS1992

Numerical study of the small scale structures in Boussinesq convection

Two-dimensional Boussinesq convection is studied numerically using two different methods: a filtered pseudospectral method and a high order accurate Essentially Nonoscillatory (ENO) scheme. The issue whether finite time singularity occurs for initially smooth flows is investigated. The numerical results suggest that the collapse of the bubble cap is unlikely to occur in resolved calculations. The strain rate corresponding to the intensification of the density gradient across the front saturates at the bubble cap. We also found that the cascade of energy to small scales is dominated by the formulation of thin and sharp fronts across which density jumps.

Weinan, E.↗

NASA NTRS1992

On an origin of numerical diffusion: Violation of invariance under space-time inversion

The invariant properties of the convection equation du/dt + adu/dx = 0 (where d is the partial differential operator) with respect to spatial reflection, time reversal, and space-time inversion are studied. Generally, a finite-difference analog of this equation may possess some or none of these properties. It is shown that, under certain conditions, the von Neumann amplification factor of an analog satisfies a special relation for each invariant property this analog possesses. Particularly, an analog is neutrally stable and thus free of numerical diffusion if it possesses the invariant property related to space-time inversion. It is also explained why generally (1) an upwind scheme possesses neither the invariant property related to spatial reflection nor that related to space-time inversion, and (2) an explicit scheme possesses neither the invariant property related to time reversal nor that related to space-time inversion. Extension to the viscous case and a remarkable connection between the current work and a new numerical framework for solving conservation laws are also discussed.

Chang, Sin-Chung↗

NASA NTRS1992

Numerical investigation of the Earth's rotation during a complete precession cycle

A theory for the long-term rotational motion of the quasi-rigid Earth was constructed by numerical integration. The theory spans 72,000 years centered about 1968 A.D., and provides accurate rotational and positional data for the Earth in the recent past and the near future. The physical model is termed dynamically consistent because developments for the active forces and torques are truncated based solely on their magnitudes regardless of their origin. The model includes all appropriate forces and torques due to the geopotential and tidal effects as well as lunisolar and planetary contributions. The elastic and inelastic deformations due to tidal action were too small to affect the mass properties of the Earth at the truncation level of the model. However, long-term dissipative effects of the tidal forces and torques were not negligible. These considerations gave the model its quasi-rigid characterization. The numerical output provided both rotational and orbital-element data. The data were fitted throughout the 72,000-year range using Chebyshev polynomial series.

Richardson, David L.↗

NASA NTRS1991

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗

NASA NTRS1992

A brief description of a new numerical framework for solving conservation laws: The method of space-time conservation element and solution element

A new numerical method for solving conservation laws is being developed. It differs substantially from the well established methods, i.e., finite difference, finite volume, finite element, and spectral methods, in both concept and methodology. It is much simpler than a typical high resolution method. No flux limiter or any technique related to characteristics is involved. No artificial viscosity or smoothing is introduced, and no moving mesh is used. Yet this method is capable of generating highly accurate shock tube solutions. The slight numerical overshoot and/or oscillations generated can be removed if a simple averaging formula initially used is replaced by a weighted formula. This modification has little effect on other parts of the solution. Because of its simplicity, generalization of this new method for multi-dimensional problems is straightforward.

Chang, Sin-Chung↗

NASA NTRS1993

Numerical simulations of galactic wakes

We have performed a series of numerical simulations which examine the simple but realistic case of a galaxy with a smooth potential corresponding to a King profile, and with no interstellar medium. This corresponds to the case of an early-type cluster galaxy that has previously been stripped of its gas. The simulations provide a numerical basis for future work involving galaxies with interstellar media. We use a time-explicit finite difference code (Lufkin & Hawley 1992, ApJ submitted) to obtain solutions to the nonlinear equations for adiabatic flow in axisymmetry. The computational grid covers a region extending 500 kpc from a stationary galaxy radially and along the symmetry axis. We use a graded mesh, allowing for full resolution near the galaxy. Because the assumed potential has a finite depth, no artificial inner boundary condition is necessary; reflection symmetry at the axis is assumed. With a total grid size of 64 radial x 128 vertical, each run requires approximately 5 minutes of cpu time on the NCSA Cray-2.

Lufkin, Eric A.↗

NASA NTRS1993

A new observational and numerical study of tidal interactions in M81-M82-NGC3077 system

A nearby system of interacting galaxies M81-M82-NGC3077 triplet (D = 3.3 Mpc; Freeman & Madore 1988) has been studied using multi-wavelength observations and numerical simulations to obtain a comprehensive understanding on the dynamics and the consequences of tidal interactions in a group environment. The VLA 12-field Mosaic H I observations of 2 x 1.5 deg. region have revealed a vast array of H I filaments which suggests that the severity and extent of tidal disruptions far exceed the previous estimates. A tidal remnant of the former H I disk of M82 extending up to 30 kpc (in projection) is identified for the first time, and the pervasive effects of the tidal disruption are traced into the inner disk by optical and CO observations, including a kinematic trace of a large scale bar potential (Yun, Ho, & Lo 1992). The H I disk of M81 is traced out to 40 kpc in radius, and a large scale (l approx. 20 kpc) velocity anomaly ('High Velocity Trough'), which may be a remnant of a gaseous collision, is found within the disk of M81. The large H I bridge between M81 and NGC 3077 (van der Hulst 1979) is also found to extend approx. greater than 50 kpc further, bending around NGC 3077, toward M82. The total H I detected in this experiment, 5.6 x 10(exp 9) solar mass, represents the majority of the single-dish flux (Appleton, Davies, & Stephenson 1981) and suggests that the bulk of H I found in the region belongs to the three galaxies and the tidal filaments. The impact and details of the tidal interactions have been further examined through the use of numerical techniques. The 'restricted 3-body' approach was used to simulate the observed distribution of tidal H I streamers connecting the three galaxies, and the success of the simulation is further strengthened by the accurate predictions on the gas kinematics.

Yun, M. S.↗

NASA NTRS1993

Numerical simulation of free shear flows: Towards a predictive computational aeroacoustics capability

Implicit and explicit spatial differencing techniques with fourth order accuracy have been developed. The implicit technique is based on the Pade compact scheme. A Dispersion Relation Preserving concept has been incorporated into both of the numerical schemes. Two dimensional Euler computation of a spatially-developing free shear flow, with and without external excitation, has been performed to demonstrate the capability of numerical schemes developed. Results are in good agreement with theory and experimental observation regarding the growth rate of fluctuating velocity, the convective velocity, and the vortex-pairing process.

Soh, W. Y.↗

NASA NTRS1993

Meso-beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

The first section is on 3-D numerical modeling of terrain-induced circulations and covers the following: (1) additional insights into gravity wave generation mechanisms based on the control simulation; (2) ongoing nested-grid numerical simulations; (3) work to be completed during the remainder of FY-93; and (4) work objectives for FY-94. The second section is on linear theory and theoretical modeling and covers the following: (1) the free response of a uniform barotropic flow to an initially stationary unbalanced (ageostrophic) zonal wind anomaly; and (2) the free response of a uniform barotropic flow to an initially stationary balanced zonal wind anomaly.

Lin, Yuh-Lang↗

NASA NTRS1992

Numerical modeling of runback water on ice protected aircraft surfaces

A numerical simulation for 'running wet' aircraft anti-icing systems is developed. The model includes breakup of the water film, which exists in regions of direct impingement, into individual rivulets. The wetness factor distribution resulting from the film breakup and the rivulet configuration on the surface are predicted in the numerical solution procedure. The solid wall is modeled as a multilayer structure and the anti-icing system used is of the thermal type utilizing hot air and/or electrical heating elements embedded with the layers. Details of the calculation procedure and the methods used are presented.

Al-Khalil, Kamel M.↗

NASA NTRS1993

Numerical simulations of three-dimensional laminar flow over a backward facing step; flow near side walls

This paper reports the results of numerical simulations of steady, laminar flow over a backward-facing step. The governing equations used in the simulations are the full 'compressible' Navier-Stokes equations, solutions to which were computed by using a cell-centered, finite volume discretization. The convection terms of the governing equations were discretized by using the Advection Upwind Splitting Method (AUSM), whereas the diffusion terms were discretized using central differencing formulas. The validity and accuracy of the numerical solutions were verified by comparing the results to existing experimental data for flow at identical Reynolds numbers in the same back step geometry. The paper focuses attention on the details of the flow field near the side wall of the geometry.

Steinthorsson, Erlendur↗

NASA NTRS1992

Testing and numerical modeling of hypervelocity impact damaged Space Station multilayer insulation

Results are presented of experiments measuring the degradation of the insulating capabilities of the multilayer insulation (MLI) of the Space Station Freedom, when subjected to hypervelocity impact damage. A simple numerical model was developed for use in an engineering design environment for quick assessment of thermal effect of the impact. The model was validated using results from thermal vacuum tests on MLI with simulated damage. The numerical model results agreed with experimental data.

Rule, William K.↗

NASA NTRS1992

Numerical simulation of a low-emission gas turbine combustor using KIVA-II

A modified version of the KIVA-II code was used to obtain a multidimensional numerical solution for the turbulent two-phase chemically reacting flows inside a staged turbine combustor (STC). The STC under consideration is equipped with an advanced airblast fuel nozzle and encompasses a fuel nozzle (FN), a rich-burn (RB) zone, a converging connecting pipe, a quick-quench (QQ) zone, a diverging connecting pipe, and a lean-combustion (LC) zone. The STC was divided into two subsystems, namely, FN/RB zone and QQ/LC zones, and the numerical solutions were obtained separately for each subsystem. Preliminary data characterize the major features of the flow and temperature fields inside the STC. Information on velocity, temperature, and some critical species in the FN/RB zone is presented. In the QQ/LC zones, formation of the co- and counter-rotating bulk flow and the sandwiched-ring-shaped temperature field can be clearly seen. The calculations of the mass-weighted standard deviation and the pattern factor of temperature indicated that the mixing performance of the STC is very promising.

Yang, S. L.↗