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,063 records · Page 59

All sources NASA NTRS DOE OSTI

NASA NTRS1989

The implementation of the graphics of program EAGLE: A numerical grid generation code on NASA Langley SNS computer system

Program EAGLE (Eglin Arbitrary Geometry Implicit Euler) Numerical Grid Generation System is a composite (multi-block) algebraic or elliptic grid generation system designed to discretize the domain in and/or around any arbitrarily shaped three dimensional regions. This system combines a boundary conforming surface generation scheme and includes plotting routines designed to take full advantage of the DISSPLA Graphics Package (Version 9.0). Program EAGLE is written to compile and execute efficiently on any Cray machine with or without solid state disk (SSD) devices. Also, the code uses namelist inputs which are supported by all Cray machines using the FORTRAN compiler CFT77. The namelist inputs makes it easier for the user to understand the inputs and operation of Program EAGLE. EAGLE's numerical grid generator is constructed in the following form: main program, EGG (executive routine); subroutine SURFAC (surface generation routine); subroutine GRID (grid generation routine); and subroutine GRDPLOT (grid plotting routines). The EAGLE code was modified to use on the NASA-LaRC SNS computer (Cray 2S) system. During the modification a conversion program was developed for the output data of EAGLE's subroutine GRID to permit the data to be graphically displayed by IRIS workstations, using Plot3D. The code of program EAGLE was modified to make operational subroutine GRDPLOT (using DI-3000 Graphics Software Packages) on the NASA-LaRC SNS Computer System. How to implement graphically, the output data of subroutine GRID was determined on any NASA-LaRC graphics terminal that has access to the SNS Computer System DI-300 Graphics Software Packages. A Quick Reference User Guide was developed for the use of program EAGLE on the NASA-LaRC SNS Computer System. One or more application program(s) was illustrated using program EAGLE on the NASA LaRC SNS Computer System, with emphasis on graphics illustrations.

Houston, Johnny L.↗

NASA NTRS1990

Numerical investigation of laminar-turbulent transition in a flat plate wake

Lamina-turbulent transition of high-deficit flat plate wakes is investigated by direct numerical simulations using the complete Navier Stokes equations. The simulations are based on a spatial model so that both the base flow and the disturbance flow can develop in the downstream direction. The Navier Stokes equations are used in a vorticity-velocity form and are solved using a combination of finite difference and spectral approximations. Fourier series are used in the spanwise direction. Second-order finite-differences are used to approximate the spatial derivatives in the streamwise and transverse directions. For the temporal discretion, a combination of ADI, Crank-Nicolson, and Adams-Bashforth methods is employed. The discretized velocity equations are solved using fast Helmholtz solvers. Code validation is accomplished by comparison of the numerical results to both linear stability and to experiments. Calculations of two- and/or three-dimensional sinuous and mode disturbances in the wake of flat plate are undertaken. For calculations of two-dimensional disturbances, the wake is forced at an amplitude level so that nonlinear disturbance development may be observed. In addition, the forcing amplitude is varied in order to determine its effect on the disturbance behavior. To investigate the onset of three-dimensionality, the wake is forced with a small-amplitude three-dimensional disturbance and a larger amplitude two-dimensional disturbance. The two-dimensional forcing amplitude is varied in order to determine its influence on the three-dimensional flow field.

Fasel, Hermann F.↗

NASA NTRS1990

Development of comprehensive numerical schemes for predicting evaporating gas-droplets flow processes of a liquid-fueled combustor

An existing Computational Fluid Dynamics code for simulating complex turbulent flows inside a liquid rocket combustion chamber was validated and further developed. The Advanced Rocket Injector/Combustor Code (ARICC) is simplified and validated against benchmark flow situations for laminar and turbulent flows. The numerical method used in ARICC Code is re-examined for incompressible flow calculations. For turbulent flows, both the subgrid and the two equation k-epsilon turbulence models are studied. Cases tested include idealized Burger's equation in complex geometries and boundaries, a laminar pipe flow, a high Reynolds number turbulent flow, and a confined coaxial jet with recirculations. The accuracy of the algorithm is examined by comparing the numerical results with the analytical solutions as well as experimented data with different grid sizes.

Chen, C. P.↗

NASA NTRS1990

Consistency and convergence for numerical radiation conditions

The problem of imposing radiation conditions at artificial boundaries for the numerical simulation of wave propagation is considered. Emphasis is on the behavior and analysis of the error which results from the restriction of the domain. The theory of error estimation is briefly outlined for boundary conditions. Use is made of the asymptotic analysis of propagating wave groups to derive and analyze boundary operators. For dissipative problems this leads to local, accurate conditions, but falls short in the hyperbolic case. A numerical experiment on the solution of the wave equation with cylindrical symmetry is described. A unified presentation of a number of conditions which have been proposed in the literature is given and the time dependence of the error which results from their use is displayed. The results are in qualitative agreement with theoretical considerations. It was found, however, that for this model problem it is particularly difficult to force the error to decay rapidly in time.

Hagstrom, Thomas↗

NASA NTRS1990

Numerical methods for systems of conservation laws of mixed type using flux splitting

The essentially non-oscillatory (ENO) finite difference scheme is applied to systems of conservation laws of mixed hyperbolic-elliptic type. A flux splitting, with the corresponding Jacobi matrices having real and positive/negative eigenvalues, is used. The hyperbolic ENO operator is applied separately. The scheme is numerically tested on the van der Waals equation in fluid dynamics. Convergence was observed with good resolution to weak solutions for various Riemann problems, which are then numerically checked to be admissible as the viscosity-capillarity limits. The interesting phenomena of the shrinking of elliptic regions if they are present in the initial conditions were also observed.

Shu, Chi-Wang↗

NASA NTRS1989

Theoretical and numerical study of a three-dimensional turbulent boundary layer

A boundary layer created on an infinite flat plate by a time-dependent freestream velocity vector whose magnitude is independent of time but whose direction changes at a constant angular velocity is theoretically studied using Reynolds-number scaling laws and numerical simulations performed over a range of Reynolds numbers. Results obtained with a higher-order version of existing theories of the Ekman layer are shown to agree well with the numerical results at three Reynolds numbers. The present results can be extrapolated to the case of high Reynolds numbers. The Reynolds-averaged equations reduce to a one-dimensional steady problem, making possible the easy and accurate testing of turbulence models.

Spalart, Philippe R.↗

NASA NTRS1989

Explosive east coast cyclogenesis - Numerical experimentation and model-based diagnostics

Numerical experimentation of explosive east-coast cyclogenesis is performed using the Florida State University Global Spectral Model (FSUGSM). The three cases examined here are the Presidents' Day storm of February 18-19, 1979 and the North Atlantic and Pacific bombs of January 18-20, 1979 which formed off the east coasts of the United States and Japan, respectively. The use of a global model provides a framework for studying the phenomena on the 3-5 day time scale. The forecast verifications of the numerical experiments indicate that the FSUGSM was able to adequately predict the phase, intensity, and synoptic-scale structure. These results justify the use of model data for diagnostic studies of the bomb. The model data are used to quantify the role of the adiabatic and diabatic forcing in the explosive cyclogenetic process, using surface pressure tendency to gage development.

Manobianco, John↗

NASA NTRS1990

A study of numerical methods for hyperbolic conservation laws with stiff source terms

In the present study of the behavior of typical numerical methods in the case of a model advection equation having a parameter-dependent source term, two approaches to the incorporation of the source terms are used: MacCormack-type predictor-corrector methods with flux limiters, and splitting methods in which the fluid dynamics and chemistry are handled in separate steps. The latter are found to perform slightly better. The model scalar equation is used to show that the incorrectness of the propagation speeds of discontinuities observed in the stiff case is due to the introduction of nonequilibrium values through numerical dissipation in the advection step.

Leveque, R. J.↗

NASA NTRS1990

Vortical flows over delta wings and numerical prediction of vortex breakdown

Navier-Stokes solutions of subsonic vortical flow over a 75 deg sweep delta wing with a sharp leading edge are presented. The sensitivity of the solution to the numerical scheme is examined using both a partially upwind scheme and a scheme with central differencing in all directions. At moderate angles of attack, no vortex breakdown is observed, whereas the higher angle-of-attack cases exhibit breakdown. The effect of numerical grid density is investigated, and solutions that are obtained with various grid densities are compared with experimental data. An embedded grid approach is implemented to enable higher resolution in selected isolated flow regions, such as the leeward-side surface, the leading-edge vortical flow, and the vortex breakdown region.

Ekaterinaris, J. A.↗

NASA NTRS1990

Numerical simulation of the incompressible internal flow through a tilting disk valve

A numerical simulation of the incompressible viscous flow through a prosthetic tilting disk heart valve is presented in order to demonstrate the current capability to model unsteady flows with moving boundaries. Both steady and unsteady flow calculations are performed by solving the incompressible Navier-Stokes equations in three-dimensional generalized curvilinear coordinates. In order to handle the moving boundary problems, the chimera grid embedding scheme which decomposes a complex computational domain into several simple subdomains is used. An algebraic turbulence model for internal flows is incorporated to reach the physiological values of Reynolds number. Good agreement is obtained between the numerical results and experimental measurements. It is found that the tilting disk valve causes large regions of separated flow, and regions of high shear.

Chang, I-Dee↗

NASA NTRS1989

Numerical analysis of natural convection in liquid droplets by phase change

A numerical analysis is performed on thermocapillary buoyancy convection induced by phase change in a liquid droplet. A finite-difference code is developed using an alternating-direction implicit (ADI) scheme. The intercoupling relation between thermocapillary force, buoyancy force, fluid property, heat transfer, and phase change, along with their effects on the induced flow patterns, are disclosed. The flow is classified into three types: thermocapillary, buoyancy, and combined convection. Among the three mechanisms, the combined convection simulates the experimental observations quite well, and the basic mechanism of the observed convection inside evaporating sessile drops is thus identified. It is disclosed that evaporation initiates unstable convection, while condensation always brings about a stable density distribution which eventually damps out all fluid disturbances. Another numerical model is presented to study the effect of boundary recession due to evaporation, and the 'peeling-off' effect (the removal of the surface layer of fluid by evaporation) is shown to be relevant.

Duh, J. C.↗

NASA NTRS1990

Numerical solutions of the linearized Euler equations for unsteady vortical flows around lifting airfoils

A linearized unsteady aerodynamic analysis is presented for unsteady, subsonic vortical flows around lifting airfoils. The analysis fully accounts for the distortion effects of the nonuniform mean flow on the imposed vortical disturbances. A frequency domain numerical scheme which implements this linearized approach is described, and numerical results are presented for a large variety of flow configurations. The results demonstrate the effects of airfoil thickness, angle of attack, camber, and Mach number on the unsteady lift and moment of airfoils subjected to periodic vortical gusts. The results show that mean flow distortion can have a very strong effect on the airfoil unsteady response, and that the effect depends strongly upon the reduced frequency, Mach number, and gust wave numbers.

Scott, James R.↗

NASA NTRS1990

Numerical simulation of an F-16A at angle of attack

The transonic flowfield around an F-16A fighter configuration at a moderate incidence angle is simulated by solving the Navier-Stokes equations on a single-block grid. The numerical solution matches experimental freestream conditions with a mach number of 0.85, 16 degrees angle of attack, and a characteristic Reynolds number of 12.75 million. MacCormack's explicit algorithm is used in conjuction with a local time step and consecutive mesh refinement procedure to accelerate numerical convergence. The Baldwin-Lomax algebraic model provides turbulent closure. Computed surface pressure distributions and the aircraft lift coefficient compare favorably with wind tunnel data. The drag coefficient in the simulation overpredicts the experimental value by 8 percent.

Huband, G. W.↗

NASA NTRS1989

Basic numerical methods

Some of the basic finite-difference schemes that can be used to solve the nonlinear equations that describe unsteady inviscid and viscous transonic flow are reviewed. Numerical schemes for solving the unsteady Euler and Navier-Stokes, boundary-layer, and nonlinear potential equations are described. Emphasis is given to the elementary ideas used in constructing various numerical procedures, not specific details of any one procedure.

Steger, Joseph L.↗

NASA NTRS1990

Numerical solutions of the linearized Euler equations for unsteady vortical flows around lifting airfoils

Scott, James R.↗

NASA NTRS1990

A numerical study of the separation wavenumber in the two-scale scattering approximation

The choice of the separation wavenumber Kd for the two-scale scattering approximation is examined by the comparison of two-scale calculations with exact numerical solutions using the moment method. It was found that the optimal choice of Kd is near k/2 when attempting to optimize for all incidence angles. When optimizing for specific incidence angles, Kd varied from k/5 at low incidence angles to just below k at higher incidence angles. These numerical results are in rough agreement with previous theoretical derivations of the separation wavenumber.

Durden, Stephen L.↗

NASA NTRS1990

A numerical study of the turbulent Ekman layer

The three-dimensional time-dependent turbulent flow in a neutrally stratified Ekman layer over a smooth flat surface was numerically simulated by directly solving the Navier-Stokes equations. Issues addressed using the direct numerical simulation (DNS) fields include the presence or absence of large-scale coherent structures ('longitudinal' or 'roll' vortices) in neutrally stratified Ekman-layer turbulence, the effects of the horizontal component of the angular velocity vector (i.e., latitude), and implications for models of the PBL. Experimental and DNS profiles are compared.

Coleman, G. N.↗

NASA NTRS1990

Two-dimensional atmospheric transport and chemistry model - Numerical experiments with a new advection algorithm

Extensive testing of the advective scheme proposed by Prather (1986) has been carried out in support of the California Institute of Technology-Jet Propulsion Laboratory two-dimensional model of the middle atmosphere. The original scheme is generalized to include higher-order moments. In addition, it is shown how well the scheme works in the presence of chemistry as well as eddy diffusion. Six types of numerical experiments including simple clock motion and pure advection in two dimensions have been investigated in detail. By comparison with analytic solutions, it is shown that the new algorithm can faithfully preserve concentration profiles, has essentially no numerical diffusion, and is superior to a typical fourth-order finite difference scheme.

Shia, Run-Lie↗