Engineering Papers⌕ Search

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.

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

At least 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Numerical modeling of chemical vapor deposition (CVD) in a horizontal reactor

In the present numerical prediction of the deposition rate of silicon from silane in a CVD process, the conservation equations for mass, momentum, energy, and chemical species are solved on a staggered grid using the SIMPLE algorithm, while the rate of chemical reactions in the gas phase and on the susceptor surface is obtained from an Arrhenius rate equation. Predicted deposition rates as a function of position along the susceptor with and without the gas phase chemical reaction are compared with the available experimental and numerical data; agreement is excellent except at the leading edge of the susceptor, where the deposition rate is overpredicted.

Sheikholeslami, M. Z.↗
NASA NTRS1988

SIRTF controller simulation - Instability masked by numerical integration

In the course of a simulation study of a candidate design for the Space Infrared Telescope Facility (SIRTF), an unusual phenomenon was observed. The uncompensated control system was unstable, but a numerical simulation with the fixed-step-size classical fourth-order Runge-Kutta method gave a stable response. This phenomenon is described in the setting in which it occurred. The Runge-Kutta simulation model is analyzed as a discrete linear system and shown to be stable, thus corroborating the numerical results.

Rajan, N.↗
NASA NTRS1988

Numerical analysis of laminar and turbulent incompressible flows using the finite element Fluid Dynamics Analysis Package (FIDAP)

The purpose of the study is the evaluation of the numerical accuracy of FIDAP (Fluid Dynamics Analysis Package). Accordingly, four test problems in laminar and turbulent incompressible flows are selected and the computational results of these problems compared with other numerical solutions and/or experimental data. These problems include: (1) 2-D laminar flow inside a wall-driven cavity; (2) 2-D laminar flow over a backward-facing step; (3) 2-D turbulent flow over a backward-facing step; and (4) 2-D turbulent flow through a turn-around duct.

Sohn, Jeong L.↗
NASA NTRS1988

On the coupling of hyperbolic and parabolic systems: Analytical and numerical approach

The coupling of hyperbolic and parabolic systems is discussed in a domain Omega divided into two distinct subdomains omega(+) and omega(-). The main concern is to find the proper interface conditions to be fulfilled at the surface separating the two domains. Next, they are used in the numerical approximation of the problem. The justification of the interface conditions is based on a singular perturbation analysis, i.e., the hyperbolic system is rendered parabolic by adding a small artifical viscosity. As this goes to zero, the coupled parabolic-parabolic problem degenerates into the original one, yielding some conditions at the interface. These are taken as interface conditions for the hyperbolic-parabolic problem. Actually, two alternative sets of interface conditions are discussed according to whether the regularization procedure is variational or nonvariational. It is shown how these conditions can be used in the frame of a numerical approximation to the given problem. Furthermore, a method of resolution is discussed which alternates the resolution of the hyperbolic problem within omega(-) and of the parabolic one within omega(+). The spectral collocation method is proposed, as an example of space discretization (different methods could be used as well); both explicit and implicit time-advancing schemes are considered. The present study is a preliminary step toward the analysis of the coupling between Euler and Navier-Stokes equations for compressible flows.

Gastaldi, Fabio↗
NASA NTRS1988

Some observations on boundary conditions for numerical conservation laws

Four choices of outflow boundary conditions are considered for numerical conservation laws. All four methods are stable for linear problems, for which examples are presented where either a boundary layer forms or the numerical scheme, together with the boundary condition, is unstable due to the formation of a reflected shock. A simple heuristic argument is presented for determining the suitability of the boundary condition.

Kamowitz, David↗
NASA NTRS1989

Numerical simulation and comparison of symmetrical/supercritical airfoils for the near tip region of a helicopter in forward flight

Aerodynamic loads on a multi-bladed helicopter rotor in forward flight at transonic tip conditions are calculated. The unsteady, three-dimensional, time-accurate compressible Reynolds-averaged thin layer Navier-Stokes equations are solved in a rotating coordinate system on a body-conformed, curvilinear grid of C-H topology. Detailed boundary layer and global numerical comparisons of NACA-0012 symmetrical and CAST7-158 supercritical airfoils are made under identical forward flight conditions. The rotor wake effects are modeled by applying a correction to the geometric angle of attack of the blade. This correction is obtained by computing the local induced downwash velocity with a free wake analysis program. The calculations are performed on the Numerical Aerodynamic Simulation Cray 2 and the VPS32 (a derivative of a Cyber 205 at the Langley Research Center) for a model helicopter rotor in forward flight.

Badavi, F. F.↗
NASA NTRS1989

Numerical model of solar dynamic radiator for parametric analysis

Growth power requirements for Space Station Freedom will be met through addition of 25 kW solar dynamic (SD) power modules. The SD module rejects waste heat from the power conversion cycle to space through a pumped-loop, multi-panel, deployable radiator. The baseline radiator configuration was defined during the Space Station conceptual design phase and is a function of the state point and heat rejection requirements of the power conversion unit. Requirements determined by the overall station design such as mass, system redundancy, micrometeoroid and space debris impact survivability, launch packaging, costs, and thermal and structural interaction with other station components have also been design drivers for the radiator configuration. Extensive thermal and power cycle modeling capabilities have been developed which are powerful tools in Station design and analysis, but which prove cumbersome and costly for simple component preliminary design studies. In order to aid in refining the SD radiator to the mature design stage, a simple and flexible numerical model was developed. The model simulates heat transfer and fluid flow performance of the radiator and calculates area mass and impact survivability for many combinations of flow tube and panel configurations, fluid and material properties, and environmental and cycle variations. A brief description and discussion of the numerical model, it's capabilities and limitations, and results of the parametric studies performed is presented.

Rhatigan, Jennifer L.↗
NASA NTRS1989

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 NTRS1988

Spur-type instability observed on numerically simulated vortex filaments

An instability observed on vortex filaments during numerical simulations of the three-dimensional, time-dependent dynamics of vortex wakes is studied to determine when and why it occurs. It is concluded that the observed instability is a consequence of the use of straight-line vortex segments of finite length to model continuously curving vortex filaments. The instability appears to occur only when the link length is a sizable fraction of the vortex span and, therefore, is not expected in an experiment. Guidelines are then given that help avoid numerical instabilities when vortex filaments are used in flow simulations.

Rossow, Vernon J.↗
NASA NTRS1988

Providing operational guidance for the development of sea breeze thunderstorms at the Kennedy Space Center - An experiment using a mesoscale numerical model

The effectiveness of a mesoscale numerical model to provide improved local forecast guidance is evaluated with respect to sea breeze convection storms at the Kennedy Space Center. The model and operational forecast guidance production are described. A case study is presented for sea breeze convection storms and lightning events on July 1, 1986. It is found that the mesoscale numerical model outperforms purely subjective predictions of sea breeze convection. The range of applications for the model are considered.

Lyons, Walter A.↗
NASA NTRS1988

Plume formation and lithosphere erosion - A comparison of laboratory and numerical experiments

The mechanics of thermal plume formation and intrusion into the lithosphere are investigated using a combination of laboratory and numerical simulations. The sequence of events leading to lithospheric thinning and uplift by thermal plumes is established, and some numerical estimates of the time scales for each stage in this process are derived that are applicable to the mantle. It is demonstrated that the two-dimensional finite element computations successfully reproduce the qualitative features seen in the experiments, with a quantitative discrepancy of typically 30 percent or less. The results of some calculations on plume formation and intrusion into model lithospheres with a variety of rheologies are presented.

Olson, Peter↗
NASA NTRS1989

Low Reynolds number numerical solutions of chaotic flow

Numerical computations of two-dimensional flow past an airfoil at low Mach number, large angle of attack, and low Reynolds number are reported which show a sequence of flow states leading from single-period vortex shedding to chaos via the period-doubling mechanism. Analysis of the flow in terms of phase diagrams, Poincare sections, and flowfield variables are used to substantiate these results. The critical Reynolds number for the period-doubling bifurcations is shown to be sensitive to mesh refinement and the influence of large amounts of numerical dissipation. In extreme cases, large amounts of added dissipation can delay or completely eliminate the chaotic response. The effect of artificial dissipation at these low Reynolds numbers is to produce a new effective Reynolds number for the computations.

Pulliam, Thomas H.↗
NASA NTRS1989

Numerical models of two complex hypersonic flowfields

Numerical solutions of the compressible Navier-Stokes equations assuming a perfect gas and laminar flow are applied to two hypersonic experiments conducted at NASA Langley Research Center. The first is cooling by injection of a gas jet through the nose of an ogive-cone, and the second is the aerothermal environment in the gap formed by the wing and elevon section of a test model of the Space Shuttle. The simulations demonstrate that the model predicts accurate pressures and shock stand-off for the gas jet and accurate pressures for the wing-elevon cove. Experimental heating rates for the gas jet indicate turbulent mixing to be occurring, an effect not modeled in the numerical solution.

Macaraeg, M. G.↗
NASA NTRS1989

A numerical method for predicting hypersonic flowfields

The flow about a body traveling at hypersonic speed is energetic enough to cause the atmospheric gases to chemically react and reach states in thermal nonequilibrium. The prediction of hypersonic flowfields requires a numerical method capable of solving the conservation equations of fluid flow, the chemical rate equations for specie formation and dissociation, and the transfer of energy relations between translational and vibrational temperature states. Because the number of equations to be solved is large, the numerical method should also be as efficient as possible. The proposed paper presents a fully implicit method that fully couples the solution of the fluid flow equations with the gas physics and chemistry relations. The method flux splits the inviscid flow terms, central differences of the viscous terms, preserves element conservation in the strong chemistry source terms, and solves the resulting block matrix equation by Gauss Seidel line relaxation.

Maccormack, Robert W.↗
NASA NTRS1988

Application of similitude principle to the numerical simulation of solar atmospheric dynamics

Numerical simulation has become an essential part of solar physics because the complex nonlinear characteristics of solar phenomena makes analytical solutions difficult to obtain. Realistic simulation of the birth and decay of an active region still is not possible because of the wide range of spatial and time scales that must be considered. Therefore, proper scaling rules must be recognized for the development of appropriate models. In this paper, the similitude principle is applied to develop scaling rules. It is found that these rules are highly dependent on the physical nature of the specific problem under consideration. A set of 'similitude critiques' is presented for some specific physical conditions. Numerical examples of coronal dynamic response and active region dynamics are used to demonstrate these ideas.

Wu, S. T.↗
NASA NTRS1989

Numerical experiments on stability and transition in plane channel flow

Various secondary and tertiary instabilities in plane channael flow are explored via time-dependent numerical simulations using the incompressible Navier-Stokes equations. Comparisons are made between transitional flows at Reynolds numbers 1500, 5000, and 8000. The lambda vortex, detached shear layer, and inverted vortex regions are identified, and the origin of the latter is explained. The laminar breakdown of the Re = 1500 flow is computed with high resolution, and the nature of its ensuing hairpin eddies is clarified by numerical particle paths. The potential of center-mode, rather than wall-mode transitions is proposed, and the resulting flow structure is described.

Zang, Thomas A.↗
NASA NTRS1989

Numerical simulation of flow through a two-strut scramjet inlet

A three-dimensional, Reynolds-averaged Navier-Stokes code has been used to numerically analyze flow through a two-strut, supersonic combustion ramjet (scramjet) inlet configuration. It solves the governing equations in full conservation form using either a fully explicit or explicit-implicit method. An algebraic, two-layer eddy-viscosity model is used for turbulent flow calculations. The analysis allows inclusion of end effects that are caused by the aft placement of the cowl on the underside of the inlet. A special grid has been developed to accommodate the struts embedded in the inlet flowfield. Detailed numerical results are presented here for the two-strut configuration, and a comparison is made with the available experimental results.

Kumar, Ajay↗
NASA NTRS1989

On some numerical schemes for transonic flow problems

New second-order-accurate finite-difference approximations for a class of nonlinear PDEs of mixed type, which includes the two-dimensional low-frequency transonic small-disturbance equation (TSD) and full-potential equation (FP), are presented. For the TSD equation, the scheme is implemented via a time-splitting algorithm; the inclusion of flux limiters keeps the total variation nonincreasing and eliminates spurious oscillations near shocks. Global-linear-stability, TVD, and entropy-stability results are proved. Numerical results for the flow over a thin airfoil are presented. Current techniques used to solve the TSD equation may easily be extended to second-order accuracy by this method. For the FP equation, the new scheme requires no subsonic/supersonic switching and no numerical flux biasing. Global linear stability for all values of the Mach number is proved.

Mostrel, Marco Mosche↗