SEARCH · Search NASA

Results for “time-dependent numerical methods”

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 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1977

Boundary-fitted curvilinear coordinate systems for solution of partial differential equations on fields containing any number of arbitrary two-dimensional bodies

A method is presented for automatic numerical generation of a general curvilinear coordinate system with coordinate lines coincident with all boundaries of a general multi-connected two-dimensional region containing any number of arbitrarily shaped bodies. No restrictions are placed on the shape of the boundaries, which may even be time-dependent, and the approach is not restricted in principle to two dimensions. With this procedure the numerical solution of a partial differential system may be done on a fixed rectangular field with a square mesh with no interpolation required regardless of the shape of the physical boundaries, regardless of the spacing of the curvilinear coordinate lines in the physical field, and regardless of the movement of the coordinate system in the physical plane. A number of examples of coordinate systems and application thereof to the solution of partial differential equations are given. The FORTRAN computer program and instructions for use are included.

Thompson, J. F.↗

NASA NTRS1972

Effects of circular geometry in simulation of convection in rotating spacecraft tanks

For computational simulation of the convection and mixing of stratified fluid in a rotating tank (such as used in Apollo flights) with time-dependent rotation, the Navier-Stokes convection problem was formulated for a circular tank configuration. The final equations results from a general approximate theory for combined forced and contained natural convection in a time-dependent rotating system. The equations are cast in terms of vorticity and stream function in a form convenient for computation, with a transformed coordinate system, and appropriate boundary conditions are derived. Accurate representations for the cryogenic supercritical oxygen thermodynamic properties are used in the computations, and an efficient numerical finite difference scheme and computational method are employed.

Martin, E. D.↗

NASA NTRS1975

Interaction of strong shock wave with turbulent boundary layer

The reported investigation represents an extension of the time-dependent solution of separated laminar flows based on the complete Navier-Stokes equations reported by MacCormack (1971) and Carter (1973). The current study includes turbulence models in conjunction with the compressible flow equations. The calculations start with a uniform flow except for values imposed along the upstream and outer boundaries. The basic numerical method is discussed along with questions concerning the exponential accuracy and the resolution of the viscous sublayer in a compressed region.

Baldwin, B. S.↗

NASA NTRS1977

Non-linear analysis of vibrations of irregular plates

A numerical perturbation method is used to investigate the forced vibrations of irregular plates. Nonlinear terms associated with the midplane stretching are retained in the analysis. The numerical part of the method involves the use of linear, finite element techniques to determine the free oscillation mode shapes and frequencies and to obtain the linear midplane stress resultants caused by the midplane stretching. Representing the solution as an expansion in terms of these linear mode shapes, these modes and the resultants are used to determine the equations governing the time-dependent coefficients of this expansion. These equations are solved by using the method of multiple scales. Specific solutions are given for the main-resonant vibrations of an elliptical plate in the presence of internal resonances. The results indicate that modes other than the driven mode can be drawn into the steady state response. Though the excitation is composed of a single harmonic, the response may not be periodic. Moreover, the particular types of responses that can occur are highly dependent on the mode being excited and are sensitive to small geometrical changes.

Lobitz, D. W.↗

NASA NTRS1977

TOMCAT - A code for numerical generation of boundary-fitted curvilinear coordinate systems on fields containing any number of arbitrary two-dimensional bodies

A method for automatic generation of boundary-fitted curvilinear coordinate systems, where the transformed coordinates are solutions of an elliptic differential system in the physical plane, and where the coordinate lines are coincident with all boundaries of a general multiply-connected, two-dimensional region containing any number of arbitrarily shaped bodies, and is described along with a suitable computer code for implementing the method. Any partial differential system can be solved on the boundary-fitted coordinate system by appropriate transformations. The transformed equations are approximated by finite differences and solved numerically in the transformed plane. All computations, whether for generating coordinate system or then solving the transformed equations, can be done on a rectangular field with square mesh with no interpolation required on the boundaries. The physical boundaries may even be time-dependent.

Thompson, J. F.↗

NASA NTRS1978

Numerical simulation of MHD shock waves in the solar wind

The effects of the interplanetary magnetic field on the propagation speed of shock waves through an ambient solar wind are examined by numerical solutions of the time-dependent nonlinear equations of motion. The magnetic field always increases the velocity of strong shocks. Although the field may temporarily slow down weak shocks inside 1 AU, it eventually also causes weak shocks to travel faster than they would without the magnetic field at larger distances. Consistent with the increase in the shock velocity, the gas pressure ratio across a shock is reduced considerably in the presence of the magnetic field. The numerical method is used to simulate (starting at 0.3 AU) the large deceleration of a shock observed in the lower corona by ground-based radio instrumentation and the more gradual deceleration of the shock in the solar wind observed by the Pioneer 9 and Pioneer 10 spacecraft.

Steinolfson, R. S.↗

NASA NTRS1978

Finite elements and finite differences for transonic flow calculations

The paper reviews the chief finite difference and finite element techniques used for numerical solution of nonlinear mixed elliptic-hyperbolic equations governing transonic flow. The forms of the governing equations for unsteady two-dimensional transonic flow considered are the Euler equation, the full potential equation in both conservative and nonconservative form, the transonic small-disturbance equation in both conservative and nonconservative form, and the hodograph equations for the small-disturbance case and the full-potential case. Finite difference methods considered include time-dependent methods, relaxation methods, semidirect methods, and hybrid methods. Finite element methods include finite element Lax-Wendroff schemes, implicit Galerkin method, mixed variational principles, dual iterative procedures, optimal control methods and least squares.

Hafez, M. M.↗

NASA NTRS1955

Generalized indical forces on deforming rectangular wings in supersonic flight

A method is presented for determining the time-dependent flow over a rectangular wing moving with a supersonic forward speed and undergoing small vertical distortions expressible as polynomials involving spanwise and chordwise distances. The solution for the velocity potential is presented in a form analogous to that for steady supersonic flow having the familiar "reflected area" concept discovered by Evvard. Particular attention is paid to indicial-type motions and results are expressed in terms of generalized indicial forces. Numerical results for Mach numbers equal to 1.1 and 1.2 are given for polynomials of the first and fifth degree in the chordwise and spanwise directions, respectively, on a wing having an aspect ratio of 4.

Lomax, Harvard↗

NASA NTRS1979

Multi-level adaptive finite element methods. 1: Variation problems

A general numerical strategy for solving partial differential equations and other functional problems by cycling between coarser and finer levels of discretization is described. Optimal discretization schemes are provided together with very fast general solvers. It is described in terms of finite element discretizations of general nonlinear minimization problems. The basic processes (relaxation sweeps, fine-grid-to-coarse-grid transfers of residuals, coarse-to-fine interpolations of corrections) are directly and naturally determined by the objective functional and the sequence of approximation spaces. The natural processes, however, are not always optimal. Concrete examples are given and some new techniques are reviewed. Including the local truncation extrapolation and a multilevel procedure for inexpensively solving chains of many boundary value problems, such as those arising in the solution of time-dependent problems.

Brandt, A.↗

NASA NTRS1974

Analysis of dynamic stall using unsteady boundary-layer theory

The unsteady turbulent boundary layer and potential flow about a pitching airfoil are analyzed using numerical methods to determine the effect of pitch rate on the delay in forward movement of the rear flow reversal point. An explicit finite difference scheme is used to integrate the unsteady boundary layer equations, which are coupled at each instant of time to a fully unsteady and nonlinear potential flow analysis. A substantial delay in forward movement of the reversal point is demonstrated with increasing pitch rate, and it is shown that the delay results partly from the alleviation of the gradients in the potential flow, and partly from the effects of unsteadiness in the boundary layer itself. The predicted delay in flow-reversal onset, and its variation with pitch rate, are shown to be in reasonable agreement with experimental data relating to the delay in dynamic stall. From the comparisons it can be concluded (a) that the effects of time-dependence are sufficient to explain the failure of the boundary layer to separate during the dynamic overshoot, and (b) that there may be some link between forward movement of the reversal point and dynamic stall.

Scruggs, R. M.↗

NASA NTRS1974

Analysis of three-dimensional unsteady flow around oscillating wings

A method based on the Navier-Stokes equations was developed for determining analytically the three-dimensional unsteady flow patterns around oscillating wings. The Helmholz vorticity transport equations were discretized in three-dimensional finite element form from a variational formulation and integrated numerically. At each time step of the numerical integration the velocity field was calculated from the representation of the three-dimensional wing by a system of optimized distribution of vortices in space. During the numerical integration of the vorticity transport equations the time-dependent boundary conditions on the wing were specified as external constraint conditions. Examples of obtained results describing the three-dimensional unsteady flow around a wing were presented.

Bratanow, T.↗

NASA NTRS1974

Numerical simulation of high Reynolds number transonic flows

A code has been developed for simulating high Reynolds number transonic flow fields of arbitrary configuration. An explicit finite-difference method with time splitting is used to solve the time-dependent equations for compressible turbulent flow. A nonorthogonal computational mesh of arbitrary configuration facilitates the description of the flow field. The code is applied to simulate the flow over an 18 percent thick circular-arc biconvex airfoil at zero angle of attack and free-stream Mach number of 0.775. A simple mixing-length model is used to describe the turbulence and chord Reynolds numbers of 1, 2, 4, and 10 million are considered. The solution describes in sufficient detail both the shock-induced and trailing-edge separation regions, and provides the profile and friction drag.

Deiwert, G. S.↗

NASA NTRS1975

Finite difference methods for transient signal propagation in stratified dispersive media

Explicit difference equations are presented for the solution of a signal of arbitrary waveform propagating in an ohmic dielectric, a cold plasma, a Debye model dielectric, and a Lorentz model dielectric. These difference equations are derived from the governing time-dependent integro-differential equations for the electric fields by a finite difference method. A special difference equation is derived for the grid point at the boundary of two different media. Employing this difference equation, transient signal propagation in an inhomogeneous media can be solved provided that the medium is approximated in a step-wise fashion. The solutions are generated simply by marching on in time. It is concluded that while the classical transform methods will remain useful in certain cases, with the development of the finite difference methods described, an extensive class of problems of transient signal propagating in stratified dispersive media can be effectively solved by numerical methods.

Lam, D. H.↗

NASA NTRS1975

An investigation of several numerical procedures for time-asymptotic compressible Navier-Stokes solutions

The status of an investigation of four numerical techniques for the time-dependent compressible Navier-Stokes equations is presented. Results for free shear layer calculations in the Reynolds number range from 1000 to 81000 indicate that a sequential alternating-direction implicit (ADI) finite-difference procedure requires longer computing times to reach steady state than a low-storage hopscotch finite-difference procedure. A finite-element method with cubic approximating functions was found to require excessive computer storage and computation times. A fourth method, an alternating-direction cubic spline technique which is still being tested, is also described.

Rudy, D. H.↗

NASA NTRS1976

Correlation between multispectral photography and near-surface turbidities

Four-band multispectral photography obtained from an aerial platform at an altitude of about 10,000 feet has been utilized to measure near-surface turbidity at numerous sampling sites in the Ross Barnett Reservoir, Mississippi. Correlation of the photographs with turbidity measurements has been accomplished via an empirical mathematical model which depends upon visual color recognition when the composited photographs are examined on either an I squared S model 600 or a Spectral Data model 65 color-additive viewer. The mathematical model was developed utilizing least-squares, iterative, and standard statistical methods and includes a time-dependent term related to sun angle. This model is consistent with information obtained from two overflights of the target area - July 30, 1973 and October 30, 1973 - and now is being evaluated with regard to information obtained from a third overflight on November 8, 1974.

Wertz, D. L.↗

NASA NTRS1977

Numerical solution of the viscous hypersonic flow past blunted cones at angle of attack

Hypersonic viscous flow over spherically blunted cones of large half angle is computed at small angles of attack in the plane of symmetry of the flow field. Time-dependent viscous shock-layer equations in body-oriented coordinates are used to describe the flow field. The shock wave is treated as a discontinuity, across which the Rankine-Hugoniot relations are used to compute the flow conditions behind the shock. A time-marching second-order finite-difference method is used to solve the equations for a perfect gas. The local CFL (Courant-Friedrich-Lewy) time increment is used to advance the solution in time at each grid point. A fourth-order damping is used to damp the oscillations in the flow quantities. The numerical results of the present analysis for quantities such as shock standoff distance, surface-pressure distribution, and heating rates compare well with the existing theoretical and experimental results.

Kumar, A.↗

NASA NTRS1979

Laminar and turbulent flows over spherically blunted cone and hyperboloid with massive surface blowing

Numerical solutions are presented for the flow over a spherically blunted cone and hyperboloid with massive surface blowing. Time-dependent viscous shock-layer equations are used to describe the flow field. The boundary conditions on the body surface include a prescribed blowing-rate distribution. The governing equations are solved by a time-asymptotic finite-difference method. Results presented here are only for a perfect gas-type flow at zero angle of attack. Both laminar and turbulent flow solutions are obtained. It is found that the effect of the surface blowing on the laminar flow field is to smooth out the curvature discontinuity at the sphere-cone juncture point, which results in a positive pressure gradient over the body. The shock slope increases on the downstream portion of the body as the surface blowing rate is increased. The turbulent flow with surface blowing is found to redevelop a boundary-layer-like region near the surface. The effects of this boundary-layer region on the flow field and heating rates are discussed.

Kumar, A.↗

NASA NTRS1979

Laminar and turbulent flows over a spherically blunted cone with massive surface blowing

Numerical solutions are presented for the flow over a spherically blunted cone with massive surface blowing. Time-dependent viscous shock layer equations are used to describe the flowfield. The boundary conditions on the body surface include a prescribed blowing rate distribution. The governing equations are solved by a time-asymptotic finite-difference method. Results presented here are only for a perfect gas-type flow at zero angle of attack. Both laminar and turbulent flow solutions are obtained. It is found that the surface blowing smooths out the effect of the curvature discontinuity at the sphere-cone juncture point on the laminar flowfield and results in a negative pressure gradient over the body. The shock slope increases on the downstream portion of the body as the surface blowing rate is increased. The turbulent flow with surface blowing is found to redevelop a boundary-layer-like region near the surface. The effects of this boundary-layer-like region on the flowfield and heating rates are discussed.

Kumar, A.↗