SEARCH · Engineering Papers

Results for “time step”

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,207 records · Page 67

All sources NASA NTRS DOE OSTI

NASA NTRS1976

Implicit finite-difference procedures for the computation of vortex wakes

Implicit finite-difference procedures for the primitive form of the incompressible Navier-Stokes and the compressible Euler equations are used to compute vortex wake flows. The partial differential equations in strong conservation-law form are transformed to cluster grid points in regions with large changes in vorticity. In addition to clustering, fourth-order accurate, spatial difference operators are used to help resolve the flow-field gradients. The use of implicit time-differencing permits large time steps to be taken since temporal variations are typically small. Computational efficiency is achieved by approximate factorization. Both two-dimensional and preliminary three-dimensional calculations are described and qualitatively compared with existing experimental data.

Steger, J. L.↗

NASA NTRS1976

A NASTRAN implementation of the doubly asymptotic approximation for underwater shock response

A detailed description is given of how the decoupling approximation known as the doubly asymptotic approximation is implemented with NASTRAN to solve shock problems for submerged structures. The general approach involves locating the nonsymmetric terms (which couple structural and fluid variables) on the right hand side of the equations. This approach results in coefficient matrices of acceptable bandwidth but degrades numerical stability, requiring a smaller time step size than would otherwise be used. It is also shown how the structure's added (virtual) mass matrix, is calculated with NASTRAN.

Everstine, G. C.↗

NASA NTRS1977

Evaluation of subgrid-scale turbulence models using a fully simulated turbulent flow

An exact turbulent flow field was calculated on a three-dimensional grid with 64 points on a side. The flow simulates grid-generated turbulence from wind tunnel experiments. In this simulation, the grid spacing is small enough to include essentially all of the viscous energy dissipation, and the box is large enough to contain the largest eddy in the flow. The method is limited to low-turbulence Reynolds numbers, in our case R sub lambda = 36.6. To complete the calculation using a reasonable amount of computer time with reasonable accuracy, a third-order time-integration scheme was developed which runs at about the same speed as a simple first-order scheme. It obtains this accuracy by saving the velocity field and its first-time derivative at each time step. Fourth-order accurate space-differencing is used.

Clark, R. A.↗

NASA NTRS1976

Dissipative, forced turbulence in two-dimensional magnetohydrodynamics

The equations of motion for turbulent two-dimensional magnetohydrodynamic flows are solved in the presence of finite viscosity and resistivity, for the case in which external forces (mechanical and/or magnetic) act on the fluid. The goal is to verify the existence of a magnetohydrodynamic dynamo effect which is represented mathematically by a substantial back-transfer of mean square vector potential to the longest allowed Fourier wavelengths. External forces consisting of a random part plus a fraction of the value at the previous time step are employed, after the manner of Lilly for the Navier-Stokes case. The regime explored is that for which the mechanical and magnetic Reynolds numbers are in the region of 100 to 1000. The conclusions are that mechanical forcing terms alone cannot lead to dynamo action, but that dynamo action can result from either magnetic forcing terms or from both mechanical and magnetic forcing terms simultaneously.

Fyfe, D.↗

NASA NTRS1976

Simulation of three-dimensional separated flows with vortex filaments

A Lagrangian vorticity method for numerical simulation of three-dimensional separated flows with vortex filaments, about solid bodies, is discussed. The method is an extension of an earlier one for three-dimensional rotational flows away from solid boundaries e.g., jets, vortex rings, and aircraft trailing vortices after initial rollup. In this method an harmonic contribution to the velocity field is computed at each time step to ensure tangency of the velocity field at the boundary. The mechanics of the boundary layer are approximated with sufficient accuracy so that the separation lines are located correctly on the surface of the body.

Leonard, A.↗

NASA NTRS1977

An implicit factored scheme for the compressible Navier-Stokes equations

An implicit finite difference scheme is developed for the numerical solution of the compressible Navier-Stokes equations in conservation-law form. The algorithm is second-order-time accurate, noniterative, and spatially factored. In order to obtain an efficient factored algorithm, the spatial cross-derivatives are evaluated explicitly. However, the algorithm is unconditionally stable and, although a three-time-level scheme, requires only two-time-levels of data storage. The algorithm is constructed in a 'delta' form (i.e., increments of the conserved variables and fluxes) that provides a direct derivation of the scheme and leads to an efficient computational algorithm. In addition, the delta form has the advantageous property of a steady-state (if one exists) independent of the size of the time step. Numerical results are presented for a two-dimensional shock boundary-layer interaction problem.

Beam, R. M.↗

NASA NTRS1977

Implicit finite difference simulation of flow about arbitrary geometries with application to airfoils

Finite difference procedures are used to solve either the Euler equations or the 'thin layer' Navier-Stokes equations subject to arbitrary boundary conditions. An automatic grid generation program is employed, and because an implicit finite difference algorithm for the flow equations is used, time steps are not severely limited when grid points are finely distributed. Computational efficiency and compatibility to vectorized computer processors is maintained by use of approximate factorization techniques. Computed results for both inviscid and viscous flow about airfoils are described and compared to various known solutions.

Steger, J. L.↗

NASA NTRS1978

Predicting surface heat flux

Report presents technique involving single embedded thermocouple used to predict flux and temperature for high-or-low conductivity materials that have temperature-and-pressure dependent properties. Technique solves for heat rate and temperature at given surface at each time step, rather than for entire history. Newton-Raphson technique is used for temperature solutions; then quadratic fit is employed. Report contains representative graphs and tables.

Curry, D. M.↗

NASA NTRS1977

CAVE: A computer code for two-dimensional transient heating analysis of conceptual thermal protection systems for hypersonic vehicles

A digital computer code CAVE (Conduction Analysis Via Eigenvalues), which finds application in the analysis of two dimensional transient heating of hypersonic vehicles is described. The CAVE is written in FORTRAN 4 and is operational on both IBM 360-67 and CDC 6600 computers. The method of solution is a hybrid analytical numerical technique that is inherently stable permitting large time steps even with the best of conductors having the finest of mesh size. The aerodynamic heating boundary conditions are calculated by the code based on the input flight trajectory or can optionally be calculated external to the code and then entered as input data. The code computes the network conduction and convection links, as well as capacitance values, given basic geometrical and mesh sizes, for four generations (leading edges, cooled panels, X-24C structure and slabs). Input and output formats are presented and explained. Sample problems are included. A brief summary of the hybrid analytical-numerical technique, which utilizes eigenvalues (thermal frequencies) and eigenvectors (thermal mode vectors) is given along with aerodynamic heating equations that have been incorporated in the code and flow charts.

Rathjen, K. A.↗

NASA NTRS1978

Numerical solution of the Navier-Stokes equations for super-sonic flows with strong shocks

The numerical solution of the full Navier-Stokes Equations for viscous flows with high Mach numbers and a strong detached bow shock was obtained. Two dimensional flows around a circular cylinder, and a circular cylinder with an aft-body in the form of a fairing, were considered. The solution of the compressible N.S. equations was accomplished by the method of finite differences. An implicit scheme of solution, the S.O.R., was used with the optimum acceleration parameters determined by trial and error. The tensor notation was used in writing the N-S Equations transformed into general curvilinear coordinates. The equations for the generation of the coordinate system were solved, followed by the solution of the N.S. equations, at the end of a set of given number of time steps. "Wiggles", constituted the one major problem that needed to be overcome. These oscillations give rise to quantities such as negative temperatures, which ultimately caused the computational program to break down. Certain dissipative finite-difference schemes damped these oscillations.

Devarayalu, K.↗

NASA NTRS1978

Supersonic flow over ablated nosetips using an unsteady, implicit numerical procedure

The axisymmetric supersonic flow over passive, that is, nonablating, indented nosetips of reentry vehicles is determined using an unsteady implicit numerical algorithm which solves either the inviscid Euler equations or the 'thin-layer' Navier-Stokes equations. A nonorthogonal independent variable transformation is used to map the distorted physical region, containing multiple zones of embedded subsonic flow into a rectangular computational domain at whose boundaries the required permeable or impermeable boundary conditions are simulated. Use of the implicit algorithm results in faster convergence to the steady state because of a larger allowable time step over conventional explicit schemes. The numerical results obtained compare favorably with existing experimental data for very mildly and severely indented blunt nosetips.

Kutler, P.↗

NASA NTRS1978

On acceleration of MacCormack's scheme

An acceleration of MacCormack's scheme due to Desideri and Tannehill is analyzed. It is found that for hyperbolic problems one cannot improve upon the efficiency of MacCormack's method. For parabolic problems the time step can be chosen arbitrarily large without loss of stability by an appropriate choice of the acceleration parameters. When applied to the heat equation this method is equivalent to both the Dufort-Frankel scheme and to MSOR.

Gottlieb, D.↗

NASA NTRS1977

Comparisons among second-order, fourth-order, and pseudospectral advection techniques in a mesoscale numerical model

The reported investigation is concerned with mesoscale squall-line simulation utilizing numerical techniques. Particular attention is given to fourth-order results and their implications for the dynamical evolution of mesoscale squall-line systems which contain severe local storms. An adiabatic inviscid set of prognostic equations is utilized in a z coordinate system for the fundamental experiment. The fourth-order advection is employed for all time-dependent equations. The Euler-backward time marching scheme is utilized with a 60 second time step. The horizontal mesh length is 42 km. Horizontal diffusion is accomplished by utilizing a smoother-desmoother. Lateral boundary conditions are designed to reduce the development of strong gradients in dependent variables near the boundaries by bringing interior values from the grid to the boundaries.

Kaplan, M.↗

NASA NTRS1977

The multiple nesting of mesoscale, submesoscale, and nonhydrostatic microscale numerical models for a case of mesocyclogenesis and severe storm development

A fundamental approach involving the integration of two different numerical models over several different spatial mesh lengths is considered. Model I is a hydrostatic, incompressible, and adiabatic primitive (Newtonian) equation set with no boundary layer forcing, i.e., surface heat or moisture flux or external frictional stresses. Typically this model is initialized from conventional radiosonde data and run for a 12-15 hour forecast period over a 40 km horizontal mesh. If phase I of the tornado developmental processes develops, the data is 'frozen' in space and time, interpolated to a 20 km mesh, and integrated for a shorter time period (3 hours) with a shorter time step. Model II is a nonhydrostatic compressible 18-level complement of model I. If model II indicates the development of the intensification of phase V dynamical processes, the data is frozen and interpolated to finer mesh lengths until the tornadic circulation system is approximated

Kaplan, M. L.↗

NASA NTRS1977

Severe weather prediction studies

A numerical forecast model is described which uses three-dimensional primitive equations to predict the evolution of the wind, temperature and moisture fields along with the subsynoptic-scale convective activity. There are 15 vertical levels set on a longitude-latitude horizontal grid. In this study, the horizontal grid distance is 70 km (0.625 deg latitude and 0.816 deg longitude at 40 deg N) and the time step is 90 sec. Orography is included and terrain coordinates are adopted. The forecast model is used to test two tornado-producing storms with very different synoptic structures. A 12 hr forecast is carried out for both cases. The rms errors of the forecast indicate much smaller errors in the middle troposphere than in the lower troposphere. The results demonstrate that a rather fine-mesh limited-area numerical model is capable of providing a valuable short-range forecast of the circulation and the mesoscale precipitation features.

Chang, C.-B.↗

NASA NTRS1973

User's instructions for the cardiovascular Walters model

The model is a combined, steady-state cardiovascular and thermal model. It was originally developed for interactive use, but was converted to batch mode simulation for the Sigma 3 computer. The model has the purpose to compute steady-state circulatory and thermal variables in response to exercise work loads and environmental factors. During a computer simulation run, several selected variables are printed at each time step. End conditions are also printed at the completion of the run.

Croston, R. C.↗

NASA NTRS1979

The performance of a seasonal global climatic model

A seasonal global climatic model developed by Sellers is studied. Several changes are made to correct errors in the original model and to allow its use in a time dependent simulation mode. The major changes are in the infrared radiation formulation and in the size of the time steps. The seasonal cycles of surface temperature, wind, ice, snow, albedo, horizontal heat transport, and vertical flux of solar and infrared radiation are compared to recent observations. The seasonal cycle of temperature is fairly well reproduced, but has too small an amplitude in the Northern Hemisphere high latitudes. An overestimate of albedo at the poles is found, which is related to errors in some of the other variables. The sensitivity of the model to CO2 and solar constant (O) changes is summarized. Doubling the CO2 amount causes a 2.45 K increase in global average surface temperature (T). Increasing Q by 1% increases T by 3.14K. Decreasing 0 by 1% decreases T by 4.65K. An annual average version of the model is more sensitive to changes in Q than the model with the seasonal cycle.

Robock, A.↗

NASA NTRS1980

Supersonic flow over three-dimensional ablated nosetips using an unsteady implicit numerical procedure

The three-dimensional supersonic flow over passive, that is, nonablating, indented nosetips of reentry vehicles is determined using an unsteady implicit numerical algorithm which solves either the inviscid Euler equations or the 'thin-layer' Navier-Stokes equations. A nonorthogonal independent variable transformation is used to map the distorted physical domain, containing multiple zones of embedded subsonic flow and separated flow regions into a rectangular computational volume at whose boundaries the required permeable or impermeable boundary conditions are simulated. Use of the implicit algorithm results in faster convergence to the steady state because of a larger allowable time step over conventional explicit schemes. The numerical results obtained compare favorably with existing numerical solutions and experimental data for simple spheres which validates the program. Results are also presented for analytically defined indented bodies for both laminar and turbulent flow conditions that demonstrate the program's capability for computing such flows.

Kutler, P.↗