Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solve time”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Time-dependent solutions of viscous incompressible flows in moving co-ordinates

A time-accurate solution method for the incompressible Navier-Stokes equations in generalized moving coordinates is presented. A finite-volume discretization method that satisfies the geometric conservation laws for time-varying computational cells is used. The discrete equations are solved by a fractional-step solution procedure. The solution is second-order-accurate in space and first-order-accurate in time. The pressure and the volume fluxes are chosen as the unknowns to facilitate the formulation of a consistent Poisson equation and thus to obtain a robust Poisson solver with favorable convergence properties. The method is validated by comparing the solutions with other numerical and experimental results. Good agreement is obtained in all cases.

Rosenfeld, Moshe↗
NASA NTRS1994

On the nonlinear stability of viscous modes within the Rayleigh problem on an infinite flat plate

The stability has been investigated of the unsteady flow past an infinite flat plate when it is moved impulsively from rest, in its own plane. For small times the instantaneous stability of the flow depends on the linearized equations of motion which reduce in this problem to the Orr-Sommerfeld equation. It is known that the flow for certain values of Reynolds number, frequency and wave number is unstable to Tollmien-Schlichting waves, as in the case of the Blasius boundary layer flow past a flat plate. With increase in time, the unstable waves only undergo growth for a finite time interval, and this growth rate is itself a function of time. The influence of finite amplitude effects is studied by solving the full Navier-Stokes equations. It is found that the stability characteristics are markedly changed both by the consideration of the time evolution of the flow, and by the introduction of finite amplitude effects.

Webb, J. C.↗
NASA NTRS1970

Dynamic characteristics of a variable-mass flexible missile

The general motion of a variable mass flexible missile with internal flow and aerodynamic forces is considered. The resulting formulation comprises six ordinary differential equations for rigid body motion and three partial differential equations for elastic motion. The simultaneous differential equations are nonlinear and possess time-dependent coefficients. The differential equations are solved by a semi-analytical method leading to a set of purely ordinary differential equations which are then solved numerically. A computer program was developed for the numerical solution and results are presented for a given set of initial conditions.

Meirovitch, L.↗
NASA NTRS1974

Numerical computation of two dimensional viscous blunt body flows with an impinging shock, part 2

Two-dimensional viscous blunt body flows with an impinging shock have been computed using a time-dependent finite-difference method which solves the complete set of Navier-Stokes equations for a compressible flow. For low Reynolds number flows, the entire flow field, including the bow shock and impinging shock, has been captured in the computation. For higher Reynolds number flows, the bow shock is treated as a discontinuity across which the Rankine-Hugoniot equations are applied, while the boundary layer and interaction regions are captured as before. Using this latter shock-fitting approach, a Type III shock interaction flow field has been computed with flow conditions corresponding to the space shuttle orbiter freestream conditions at 61 km (200,000 ft).

Holst, T. L.↗
NASA NTRS1975

Suppression of nonlinear oscillations in combustors with partial length acoustic liners

An analytical model is formulated for a three-dimensional nonlinear stability problem in a rocket motor combustion chamber. The chamber is modeled as a right circular cylinder with a short (multi-orifice) nozzle, and an acoustic linear covering an arbitrary portion of the cylindrical periphery. The combustion is concentrated at the injector and the gas flow field is characterized by a mean Mach number. The unsteady combustion processes are formulated using the Crocco time lag model. The resulting equations are solved using a Green's function method combined with numerical evaluation techniques. The influence of acoustic liners on the nonlinear waveforms is predicted. Nonlinear stability limits and regions where triggering is possible are also predicted for both lined and unlined combustors in terms of the combustion parameters.

Espander, W. R.↗
NASA NTRS1974

Numerical computation of viscous blunt body flows with a planar impinging shock

Two- and three-dimensional, viscous blunt body flows with planar impinging shocks are computed using an explicit, time-dependent, finite-difference method to solve the complete set of Navier-Stokes equations. The bow shock is treated as a discontinuity, while all interior shock layer detail such as shear layers, shock waves, jets and the wall boundary layer are automatically captured in the solution. Numerical results are presented for cases in which planar shock waves of different strengths and orientations are allowed to impinge on the flow field surroundings an infinite cylinder resulting in two- and three-dimensional shock interference patterns. The numerical results are compared with experiment.

Holst, T. L.↗
NASA NTRS1975

Numerical computation of two-dimensional viscous blunt body flows with an impinging shock

Two-dimensional, viscous, blunt body flows with an impinging shock wave are computed using a time-dependent, finite-difference method to solve the complete set of Navier-Stokes equations. The bow shock wave is treated as a discontinuity, while all interior shock layer detail such as shear layers, shock waves, jets, and the wall boundary layer are automatically captured in the solution. Numerical results are presented for cases in which shock waves of different strengths are allowed to impinge on the flow field surrounding a circular cylinder resulting in different shock interference patterns. The two-dimensional results are compared qualitatively with existing three-dimensional experiments.

Tannehil, J. C.↗
NASA NTRS1975

Transient airload computer analysis for simulating wind induced impulsive noise conditions of a hovering helicopter rotor

A numerical analysis was developed to determine the airloads on helicopter rotors operating under near-hovering flight conditions capable of producing impulsive noise. A computer program was written in which the solutions for the rotor tip vortex geometry, inflow, aeroelastic response, and airloads are solved in a coupled manner at sequential time steps, with or without the influence of an imposed steady ambient wind or transient gust. The program was developed for future applications in which predicted airloads would be incorporated in an acoustics analysis to attempt to predict and analyze impulsive noise (blade slap). The analysis was applied to a hovering full-scale rotor for which impulsive noise was recorded in the presence of ambient wind. The predicted tip vortex coordinates are in reasonable agreement with the test data, and the blade airload solutions converged to a periodic behavior for an imposed steady ambient wind conditions.

Hall, G. F.↗
NASA NTRS1976

An analysis of auroral E region neutral winds based on incoherent scatter radar observations at Chatanika

The physical origin of auroral E-region neutral winds during geomagnetic disturbances is investigated by comparing incoherent-backscatter radar observations made at Chatanika, Alaska, during such disturbances with detailed theoretical calculations of neutral velocities for these conditions. In the analysis, an ion-drag collisional-coupling model is adopted as the primary driving force for the neutral winds. Electric fields and electron number densities observed by the radar facility are used in the calculations, and the coupled ion and neutral momentum equations are solved numerically as functions of altitude and time in the ion-drag approximation. The altitude structure of the velocity field is examined, and the detailed response of the neutral gas to the ion-drag mechanism is delineated. Weighted altitude averages of the calculated velocities are obtained, the accuracy of the procedure for deriving neutral winds from incoherent-scatter radar observations is evaluated, and the comparison is carried out. It is concluded that both ion drag and pressure gradients appear to be required to explain the neutral winds observed in the auroral E-region during geomagnetic disturbances.

Comfort, R. H.↗
NASA NTRS1977

Numerical solution of axisymmetric boattail flow fields with plume simulators

Turbulent separating flows over axisymmetric afterbody-boattail configurations with solid sting plume simulators are computed with a time-dependent finite-difference method to solve the compressible Navier-Stokes equations. The Reynolds stress terms are replaced with a two-layer eddy viscosity model including a relaxation formula to model the nonequilibrium effects of the separated flow. The mesh is alined with the boattail body through an analytic transformation which accommodates a wide variety of boattail geometries. Numerical results for a series of boattail geometries over a wide range of Reynolds number (140,000 to 140 million) are presented and, when possible, compared with experimental data or independent numerical results.

Holst, T. L.↗
NASA NTRS1978

Numerical solution of two-dimensional turbulent blunt body flows with an impinging shock

An implicit finite-difference method has been developed to compute two-dimensional, turbulent, blunt body flows with an impinging shock wave. The full time-averaged Navier-Stokes equations are solved with algebraic eddy viscosity and turbulent Prandtl number models employed for shear stress and heat flux. The irregular-shaped bow shock is treated as a discontinuity across which the Rankine-Hugoniot equations are applied. A Type III turbulent shock interference flow field has been computed and the numerical results compare favorably with existing experimental data. In addition, comparisons are made between the present implicit code and a previous explicit code.

Tannehill, J. C.↗
NASA NTRS1979

An exploratory study of a finite difference method for calculating unsteady transonic potential flow

A method for calculating transonic flow over steady and oscillating airfoils was developed by Isogai. The full potential equation is solved with a semi-implicit, time-marching, finite difference technique. Steady flow solutions are obtained from time asymptotic solutions for a steady airfoil. Corresponding oscillatory solutions are obtained by initiating an oscillation and marching in time for several cycles until a converged periodic solution is achieved. The method is described in general terms and results for the case of an airfoil with an oscillating flap are presented for Mach numbers 0.500 and 0.875. Although satisfactory results are obtained for some reduced frequencies, it is found that the numerical technique generates spurious oscillations in the indicial response functions and in the variation of the aerodynamic coefficients with reduced frequency. These oscillations are examined with a dynamic data reduction method to evaluate their effects and trends with reduced frequency and Mach number. Further development of the numerical method is needed to eliminate these oscillations.

Bennett, R. M.↗
NASA NTRS1979

Amplitude effects on the dynamic performance of hydrostatic gas thrust bearings

A strip gas film bearing with inherently compensated inlets is analyzed to determine the effect of disturbance amplitude on its dynamic performance. The governing Reynolds' equation is solved using finite-difference techniques. The time dependent load capacity is represented by a Fourier series up to and including the third harmonics. For the range of amplitudes investigated the linear stiffness was independent of the amplitude, and the linear damping was inversely proportional to (1 - epsilon-squared) to the 1.5 power where epsilon is the amplitude relative to the film thickness.

Stiffler, A. K.↗
NASA NTRS1980

Numerical calculation of steady inviscid full potential compressible flow about wind turbine blades

The air flow through a propeller-type wind turbine rotor is characterized by three-dimensional rotating cascade effects about the inner portions of the rotor blades and compressibility effects about the tip regions of the blades. In the case of large rotor diameter and/or increased rotor angular speed, the existence of small supersonic zones terminated by weak shocks is possible. An exact nonlinear mathematical model (called a steady Full Potential Equation - FPE) that accounts for the above phenomena has been rederived. An artificially time dependent version of FPE was iteratively solved by a finite volume technique involving an artificial viscosity and a three-level consecutive mesh refinement. The exact boundary conditions were applied by generating a boundary conforming periodic computation mesh.

Dulikravich, D. S.↗
NASA NTRS1981

A numerical simulation of two-dimensional inlet flow fields

Inlet flow fields for airbreathing missiles are calculated by the adaptation of a two-dimensional computational method developed for the flow around airfoils. A supersonic free stream is assumed to allow the forebody calculation to be uncoupled from the inlet calculation. The inlet calculation employs an implicit, time-marching finite-difference procedure to solve the thin-layer Navier-Stokes equations formulated in body-fitted coordinates. Because the method can be used for a flow field with both subsonic and supersonic regions, it is applicable to subcritical as well as supercritical inlet operation.

Biringen, S.↗
NASA NTRS1983

Some recent applications of XTRAN3S

A time marching finite difference code, XTRAN3S that solves the three dimensional transonic small perturbation equation for flow over isolated wings was developed. During initial applications of the program, problems were encountered in the prediction of unsteady forces. The use of a revised grid and force calculation scheme improved those predictions. Comparisons are made between predicted and experimental pressure data for a rectangular supercritical wing. Comparisons of steady and unsteady data at M sub infinity = 0.700 show good agreement between calculated and experimental values. A comparison of steady data at M sub infinity 0.825 shows poor agreement between calculations and experiment. Program difficulties were encountered with swept and tapered configurations.

Seidel, D. A.↗
NASA NTRS1983

Some recent applications of XTRAN3S

A time-marching finite difference code, XTRAN3S, that solves the three-dimensional transonic small perturbation equation for flow over isolated wings has recently been developed. During initial applications of the program, problems were encountered in the prediction of unsteady forces. The use of a revised grid and force calculation scheme improved those predictions. Comparisons are made between predicted and experimental pressure data for a rectangular supercritical wing. Comparisons of steady and unsteady data at freestream Mach number = 0.700 show good agreement between calculated and experimental values. A comparison of steady data at freestream Mach number = 0.825 shows poor agreement between calculations and experiment. Program difficulties have been encountered with swept and tapered configurations.

Seidel, D. A.↗