SEARCH · Engineering Papers

Results for “numerics”

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,135 records · Page 63

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Effect of airfoil (trailing-edge) thickness on the numerical solution of panel methods based on the Dirichlet boundary condition

The practical limit of airfoil thickness ratio for which acceptable engineering results are obtainable with the Dirichlet boundary-condition-based numerical methods is investigated. This is done by studying the effect of thickness on the calculated pressure distribution near the trailing edge and by comparing the aerodynamic coefficients with available exact solutions. The first objective of this study, owing to the wide use of such computational methods, is to demonstrate the numerical symptoms that occur when the body or wing thickness approaches zero and to increase the awareness of potential users of these methods. Additionally, an effort is made to obtain the practical limits of the trailing-edge thickness where such problems will appear in the flow solution, and to propose some possible cures for very thin airfoils or those with cusped trailing edges.

Yon, Steven↗

NASA NTRS1992

Direct numerical simulation of laminar breakdown in high-speed, axisymmetric boundary layers

Temporal direct numerical simulation of laminar breakdown via subharmonic secondary instability in high-speed axisymmetric boundary layers has been accomplished using a highly accurate, fully explicit algorithm which combines spectral collocation and high-order compact-difference techniques. Numerical test cases confirm that subharmonic secondary instability is confirmed to be a viable path to transition in high-speed boundary-layer flow. Secondary instability is shown to account for peaks in the Reynolds stresses at or near the critical layer which are not possible from the second-mode primary instability alone. Reynolds stresses spatially reconstructed from the temporal model via the Gaster transformation show a 'spreading angle' of about 12 deg, in qualitative agreement with experimental findings. The rate of broadening of the Reynolds stress peak is a strongly nonlinear phenomenon which cannot be reproduced by secondary instability theory.

Pruett, C. D.↗

NASA NTRS1992

A comparative study between experimental results and numerical predictions of multi-wall structural response to hypervelocity impact

Over the last three decades, multiwall structures have been analyzed extensively, primarily through experiment, as a means of increasing the protection afforded to spacecraft structure. However, as structural configurations become more varied, the number of tests required to characterize their response increases dramatically. As an alternative, numerical modeling of high-speed impact phenomena is often being used to predict the response of a variety of structural systems under impact loading conditions. This paper presents the results of a preliminary numerical/experimental investigation of the hypervelocity impact response of multiwall structures. The results of experimental high-speed impact tests are compared against the predictions of the HULL hydrodynamic computer code. It is shown that the hypervelocity impact response characteristics of a specific system cannot be accurately predicted from a limited number of HULL code impact simulations. However, if a wide range of impact loadings conditions are considered, then the ballistic limit curve of the system based on the entire series of numerical simulations can be used as a relatively accurate indication of actual system response.

Schonberg, William P.↗

NASA NTRS1992

Numerical models of mantle convection

An overview of numerical methods describing the structure and dynamics of the mantle is presented with attention given to novel 3D modeling techniques. The paper reviews 3D spherical and Cartesian models for constant viscosity emphasizing the assumptions regarding style of convection, time dependence, and implications for the mantle. Similarly treated are 3D Cartesian models with temperature-dependent viscosities, and briefly examined are models that are based on compressibility, nonlinear viscosity, or plates. Extensive illustrations are presented detailing: (1) temperature variations from models of 3D thermal convection in spherical shells; (2) thermal anomalies in equatorial cross sections; and (3) temperature variations in a spherical shell heated from within. The discussion relates the numerical results of the models with real mantle-convection events, and the simulations are shown to yield increasingly realistic representations of material behavior.

Schubert, G.↗

NASA NTRS1992

Numerical simulations using a dynamic solution-adaptive grid algorithm, with applications to unsteady internal flows

An investigation into the numerical simulation of unsteady flows is undertaken using a two-stage Runge-Kutta scheme coupled with the dynamic solution-adaptive grid algorithm developed by the authors. The inviscid fluxes are described by a modified Advective Upwind Split Method to eliminate the need for artificial dissipation. A well-documented numerical example containing moving discontinuities is presented that demonstrates the ability of the coupled grid/solver scheme to accurately capture unsteady flowfield phenomena. Applications are to a typical inlet diffuser configuration at Mach 3.0 with excessive back pressure inducing inlet unstart.

Benson, Rusty A.↗

NASA NTRS1992

The void spectrum in two-dimensional numerical simulations of gravitational clustering

An algorithm for deriving a spectrum of void sizes from two-dimensional high-resolution numerical simulations of gravitational clustering is tested, and it is verified that it produces the correct results where those results can be anticipated. The method is used to study the growth of voids as clustering proceeds. It is found that the most stable indicator of the characteristic void 'size' in the simulations is the mean fractional area covered by voids of diameter d, in a density field smoothed at its correlation length. Very accurate scaling behavior is found in power-law numerical models as they evolve. Eventually, this scaling breaks down as the nonlinearity reaches larger scales. It is shown that this breakdown is a manifestation of the undesirable effect of boundary conditions on simulations, even with the very large dynamic range possible here. A simple criterion is suggested for deciding when simulations with modest large-scale power may systematically underestimate the frequency of larger voids.

Kauffmann, Guinevere↗

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.↗