SEARCH · Engineering Papers

Results for “Finite-differences”

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 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Navier-Stokes simulation of transonic wing flow fields using a zonal grid approach

The transonic Navier-Stokes code was used to simulate flow fields about isolated wings for workshop wind-tunnel and free-air cases using the thin-layer Reynolds-averaged Navier-Stokes equations. An implicit finite-difference scheme based on a diagonal version of the Beam-Warming algorithm was used to integrate the governing equations. A zonal grid approach was used to allow efficient grid refinement near the wing surface. The flow field was sensitive to the turbulent transition model, and flow unsteadiness was observed for a wind-tunnel case but not for the corresponding free-air case. The specification of experimental pressure at the wind-tunnel exit plane is the primary reason for the difference of these two numerical solutions.

Chaderjian, Neal M.↗

NASA NTRS1988

A transonic interactive boundary-layer theory for laminar and turbulent flow over swept wings

A 3-D laminar and turbulent boundary-layer method is developed for compressible flow over swept wings. The governing equations and curvature terms are derived in detail for a nonorthogonal, curvilinear coordinate system. Reynolds shear-stress terms are modeled by the Cebeci-Smith eddy-viscosity formulation. The governing equations are descretized using the second-order accurate, predictor-corrector finite-difference technique of Matsuno, which has the advantage that the crossflow difference formulas are formed independent of the sign of the crossflow velocity component. The method is coupled with a full potential wing/body inviscid code (FLO-30) and the inviscid-viscous interaction is performed by updating the original wing surface with the viscous displacement surface calculated by the boundary-layer code. The number of these global iterations ranged from five to twelve depending on Mach number, sweep angle, and angle of attack. Several test cases are computed by this method and the results are compared with another inviscid-viscous interaction method (TAWFIVE) and with experimental data.

Woodson, Shawn H.↗

NASA NTRS1988

Transient heat transfer from shrinking LOX-drop

In order to achieve prescribed experimental conditions in wind tunnels, the nitrogen-oxygen mixture is enriched by injection of liquid oxygen (LOX) upstream of methane burners. The objective of the study is to determine the LOX drop evaporation rate, which is dominated by heat transfer from the air stream (mass transfer mechanism can be later coupled with the solution from the present work). Since the concentration of drops in the surrounding gas is quite high, the concept of infinite medium cannot be utilized. Drop evaporation, as part of system's mass balance, is an important source-term in the finite-difference 2/3-dimensional network terms of laboratory-system (Eulerian) coordinates, while individual drop behavior, including its thermal history, is analyzed in Lagrangian coordinates.

Kondic, Nenad N.↗

NASA NTRS1989

Two-dimensional numerical simulation of a Stirling engine heat exchanger

The first phase of an effort to develop multidimensional models of Stirling engine components is described; the ultimate goal is to model an entire engine working space. More specifically, parallel plate and tubular heat exchanger models with emphasis on the central part of the channel (i.e., ignoring hydrodynamic and thermal end effects) are described. The model assumes: laminar, incompressible flow with constant thermophysical properties. In addition, a constant axial temperature gradient is imposed. The governing equations, describing the model, were solved using Crank-Nicloson finite-difference scheme. Model predictions were compared with analytical solutions for oscillating/reversing flow and heat transfer in order to check numerical accuracy. Excellent agreement was obtained for the model predictions with analytical solutions available for both flow in circular tubes and between parallel plates. Also the heat transfer computational results are in good agreement with the heat transfer analytical results for parallel plates.

Ibrahim, Mounir B.↗

NASA NTRS1989

Low-frequency sound propagation modeling over a locally-reacting boundary using the parabolic approximation

There is substantial interest in the analytical and numerical modeling of low-frequency, long-range atmospheric acoustic propagation. Ray-based models, because of frequency limitations, do not always give an adequate prediction of quantities such as sound pressure or intensity levels. However, the parabolic approximation method, widely used in ocean acoustics, and often more accurate than ray models for lower frequencies of interest, can be applied to acoustic propagation in the atmosphere. Modifications of an existing implicit finite-difference implementation for computing solutions to the parabolic approximation are discussed. A locally-reacting boundary is used together with a one-parameter impedance model. Intensity calculations are performed for a number of flow resistivity values in both quiescent and windy atmospheres. Variations in the value of this parameter are shown to have substantial effects on the spatial variation of the acoustic signal.

Robertson, J. S.↗

NASA NTRS1988

The origin of pristine KREEP - Effects of mixing between UrKREEP and the magmas parental to the Mg-rich cumulates

The fact that pristine KREEP basalts generally have bulk-rock molar MgO/(MgO+FeO) ratios lower than average for the lunar crust despite their extraordinarily high incompatible element contents is discussed. Pristine KREEP basalts also have primitive Ni contents compared to mare basalts with comparable REE contents. It is suggested that these facts can be explained through mixing between ancient KREEP precursor materials and primitive Mg-rich melts. Finite-difference models of this mixing, followed by anorthosite assimilation and fractional crystallization provide satisfactory fits to the composition of these basalts. It is suggested that this is a further confirmation of the magmasphere hypothesis.

Warren, Paul H.↗

NASA NTRS1988

Numerical simulation of the vortical flow over a round-edged double-delta wing

Simulations of three-dimensional vortical flows over a thin double-delta wing with an aspect ratio of 2.05 have been performed. Steady-state solutions to the unsteady incompressible Navier-Stokes equations are obtained using a new implicit upwind-relaxation finite-difference scheme. The method is second-order accurate spatially and naturally dissipative. Numerical results indicate that key features of both vortical interaction and vortex breakdown are successfully simulated. The computed lift coefficients and lateral trajectories of the vortical cores are in good agreement with the experimental data.

Hsu, C.-H.↗

NASA NTRS1988

Unsteady transonic aerodynamics and aeroelastic calculations at low-supersonic freestreams

A computational procedure is presented to simulate transonic unsteady flows and corresponding aeroelasticity of wings at low-supersonic freestreams. The flow is modeled by using the transonic small-perturbation theory. The structural equations of motions are modeled using modal equations of motion directly coupled with aerodynamics. Supersonic freestreams are simulated by properly accounting for the boundary conditions based on pressure waves along the flow characteristics in streamwise planes. The flow equations are solved using the time-accurate, alternating-direction implicit finite-difference scheme. The coupled aeroelastic equations of motion are solved by an integration procedure based on the time-accurate, linear-acceleration method. The flow modeling is verified by comparing calculations with experiments for both steady and unsteady flows at supersonic freestreams. The unsteady computations are made for oscillating wings. Comparisons of computed results with experiments show good agreement. Aeroelastic responses are computed for a rectangular wing at Mach numbers ranging from subtransonic to upper-transonic (supersonic) freestreams. The extension of the transonic dip into the upper transonic regime is illustrated.

Guruswamy, Guru P.↗

NASA NTRS1988

Backward differentiation approximations of nonlinear differential/algebraic systems

Finite-difference approximations of dynamical systems modeled by nonlinear, semiexplicit, differential/algebraic equations are analyzed. Convergence for the backward differentiation method is proved for index two and index three problems when the numerical initial values obey certain constraints. The appropriate asymptotic convergence rates and the leading error terms are determined.

Brenan, Kathryn E.↗

NASA NTRS1988

Unsteady transonic airfoil computation using the integral solution of full-potential equation

The shock-capturing integral-equation scheme developed by Kandil and Hu (1987) for the analysis of steady transonic flow over airfoils is extended to the unsteady case. The full potential formulation of the governing equations is reviewed; the solution method is outlined; and results for a NACA 0012 airfoil in forced pitching oscillation at Mach 0.755 are presented in extensive graphs and briefly characterized. The present technique is shown to require significantly less computation time than finite-difference or finite-volume methods, and to give shock-motion predictions in good agreement with those of an implicit finite-volume Euler solver; the surface-pressure peaks are slightly underpredicted.

Kandil, O. A.↗

NASA NTRS1988

Nonoscillatory solution of the steady-state inviscid Burgers' equation by mathematical programming

In order to obtain the physically relevant discontinuous numerical solution, the steady-state inviscid Burgers' equation is singularly perturbed through the addition of a small amount of viscosity. A 'cell-centered' finite-difference scheme is proposed which employs two points for the inviscid part and four points for the viscous one. While difficulties are experienced in the capture of interior layers centered at node points, computational results for interior layers centered between node points, and for boundary layers, exhibit accurate nonoscillatory solutions whose discontinuities are captured in one cell on both coarse and fine grids.

Lavery, John E.↗

NASA NTRS1988

Navier-Stokes simulation of transonic wing flow fields using a zonal grid approach

Chaderjian, Neal M.↗

NASA NTRS1988

Viscous shock-layer solutions for the low-density hypersonic flow past long slender bodies

Results are obtained for the surface pressure, drag, heat-transfer, and skin-friction coefficients for hyperboloids and sphere cones. Body half angles from 5 to 22.5 degrees are considered for various low-density flow conditions. Recently obtained surface-slip and shock-slip equations are employed to account for the low-density effects. The method of solution employed for the viscous shock-layer (VSL) equations is a partially coupled spatial-marching implicit finite-difference technique. The flow cases analyzed include highly cooled long slender bodies in high Mach number flows. The present perfect-gas VSL calculations compare quite well with available experimental data. Results have also been obtained from the steady-state Navier-Stokes (NS) equations by successive approximations. Comparison between the NS and VSL results indicates that VSL equations even with body and shock-slip boundary conditions may not be adequate in the stagnation region at altitudes greater than about 75 km for the cases analyzed here.

Gupta, R. N.↗

NASA NTRS1988

Chemical nonequilibrium and viscous flow computation for conic aerobrake bodies

Three-dimensional analyses are presented for the viscous, reactive flow over a complete entry-body configuration with a wide-angle conic surface. The predictive method uses a split approach that solves iteratively the Navier-Stokes and the continuity equations of chemical species. The finite-difference formulation and the computational grid are adapted to the bow shock and the conformally mapped body such that the velocity components are in the computational spherical-polar space. Combinations of several conic forebody and afterbody configurations have been studied using wind-tunnel, Space Shuttle, and aerobraking orbital transfer vehicle (AOTV) entry conditions. The effects of the borebody bluntness and of finite-rate chemical reactions on the shock layer, the wall catalycity on the boundary layer, the shear-layer impingement on the afterbody, and the base-flow environment are discussed.

Li, C. P.↗

NASA NTRS1988

Navier-Stokes simulation of unsteady three-dimensional blade-vortex interactions

The unsteady, viscous, three-dimensional flow field of a helicopter rotor blade encountering a passing vortex is calculated by solving the thin layer Navier-Stokes equations by a finite-difference numerical procedure. A prescribed vortex method is adopted to preserve the structure of the interacting concentrated vortex. The test cases considered correspond to the experimental model rotor test conditions of Caradonna et al.

Srinivasan, G. R.↗

NASA NTRS1988

The role of eigensystems in the stability and convergence rate of numerical time-marching systems

A set of linear homogeneous nonautonomous ordinary differential equations for which a general solution can be written is developed and discussed. The stability and convergence properties of standard finite-difference methods used to integrate these equations are also considered.

Lomax, H.↗

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↗

NASA NTRS1989

Review of FD-TD numerical modeling of electromagnetic wave scattering and radar cross section

Applications of the finite-difference time-domain (FD-TD) method for numerical modeling of electromagnetic wave interactions with structures are reviewed, concentrating on scattering and radar cross section (RCS). A number of two- and three-dimensional examples of FD-TD modeling of scattering and penetration are provided. The objects modeled range in nature from simple geometric shapes to extremely complex aerospace and biological systems. Rigorous analytical or experimental validatons are provided for the canonical shapes, and it is shown that FD-TD predictive data for near fields and RCS are in excellent agreement with the benchmark data. It is concluded that with continuing advances in FD-TD modeling theory for target features relevant to the RCS problems and in vector and concurrent supercomputer technology, it is likely that FD-TD numerical modeling will occupy an important place in RCS technology in the 1990s and beyond.

Taflove, Allen↗