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 631 records · Page 35

All sources NASA NTRS DOE OSTI

NASA NTRS1982

An implicit solution of the three-dimensional Navier-Stokes equations for an airfoil spanning a wind tunnel

An implicit finite-difference algorithm is developed for the numerical solution of the incompressible three dimensional Navier-Stokes equations in the non-conservative primitive-variable formulation. The flow field about an airfoil spanning a wind-tunnel is computed. The coordinate system is generated by an extension of the two dimensional body-fitted coordinate generation techniques of Thompson, as well as that of Sorenson, into three dimensions. Two dimensional grids are stacked along a spanwise coordinate defined by a simple analytical function. A Poisson pressure equation for advancing the pressure in time is arrived at by performing a divergence operation on the momentum equations. The pressure at each time-step is calculated on the assumption that continuity be unconditionally satisfied. An eddy viscosity coefficient, computed according to the algebraic turbulence formulation of Baldwin and Lomax, simulates the effects of turbulence.

Moitra, A.↗

NASA NTRS1982

Operator split methods in the numerical solution of the finite deformation elastoplastic dynamic problem

The spatial formulation of the elastoplastic dynamic problem for finite deformations is considered. A thermodynamic argument leads to an additive decomposition of the spatial rate of deformation tensor and allows an operator split of the evolutionary equations of the problem into elastic and plastic parts. This operator split is taken as the basis for the definition of a global product algorithm. In the context of finite element discretization the product algorithm entails, for every time step, the solution of a nonlinear elastodynamic problem followed by the application of plastic algorithms that operate on the stresses and internal variables at the integration points and bring in the plastic constitutive equations. Suitable plastic algorithms are discussed for the cases of perfect and hardening plasticity and viscoplasticity. The proposed formalism does not depend on any notion of smoothness of the yield surface and is applicable to arbitrary convex elastic regions, with or without corners. The stabiity properties of the global product algorithm are shown to be identical to those of the algorithm used for the integration of the nonlinear elastodynamic problem. Numerical examples illustrate the accuracy of the method.

Pinsky, P. M.↗

NASA NTRS1982

Aerodynamic characteristics of airfoils with ice accretions

Results of a wind tunnel test to evaluate the performance of an airfoil with simulated rime ice are presented with theoretical comparisons. A NACA 65A413 airfoil was tested in the OSU 6 x 22 inch Transonic Airfoil Wind Tunnel at a Reynolds number near three million and Mach numbers from 0.20 to 0.80. The model was tested in four configurations to determine the aero-dynamic effects of the roughness and shape of a rime ice accretion. The simulated rime ice shape was obtained analytically using a time-stepping dry ice accretion computer code. Lift, drag, moment coefficients, and pressure distributions for the clean and simulated rime ice cases are reported. The measured degradation in airfoil performance is compared to an analytical method which uses existing airfoil analysis computer codes with empirical corrections for the surface roughness. A discussion of the empirical surface roughness correction and uses of other airfoil computer methods is included.

Bragg, M. B.↗

NASA NTRS1981

Some aspects of algorithm performance and modeling in transient thermal analysis of structures

The status of an effort to increase the efficiency of calculating transient temperature fields in complex aerospace vehicle structures is described. The advantages and disadvantages of explicit and implicit algorithms are discussed. A promising set of implicit algorithms with variable time steps, known as the GEAR package is described. Four test problems, used for evaluating and comparing various algorithms, have been selected and finite-element models of the configurations are described. These problems include a Space Shuttle frame component, an insulated cylinder, a metallic panel for a thermal protection system, and a model of the Space Shuttle Orbiter wing. Results generally indicate a preference for implicit over explicit algorithms for solution of transient structural heat transfer problems when the governing equations are 'stiff' (typical of many practical problems such as insulated metal structures).

Adelman, H. M.↗

NASA NTRS1981

Multigrid simulation of asymptotic curved-duct flows using a semi-implicit numerical technique

Asymptotic flows inside curved ducts of rectangular as well as polar cross section are analyzed using the Navier-Stokes equations in terms of the axial velocity and vorticity and the cross-flow stream function. Numerical solutions of the three second-order coupled elliptic partial differential equations governing this flow are obtained efficiently using the coupled alternating-direction implicit (ADI) method as well as the multigrid strongly-implicit (SI) scheme. For the flow configuration studied, the ADI method is found to be more sensitive to the time steps used than is the SI scheme. Use of the multigrid-coupled-strongly-implicit (MG-SI) scheme makes it possible to efficiently obtain fine-grid solutions for configurations having strong secondary flow. It is shown that, for this asymptotic curved-duct flow, the similarity parameter of significance is the Dean's number K rather than the Reynolds number Re. Results are obtained for curved ducts with square cross sections for K up to 900, which here corresponds to Re = 9,000 for this internal flow configuration.

Ghia, K. N.↗

NASA NTRS1982

A vectorized, finite-volume, adaptive grid algorithm applied to planetary entry problems

An adaptive grid, finite-volume method has been applied to problems in planetary entry for computing complete flowfields. The adaption algorithm is implicit in nature and is keyed to resolve user specified gradients. The finite-volume algorithm is explicit, utilizing a maximum time step advancement at each grid point to accelerate convergence to the steady state. The present version of the code is for the laminar flow of a perfect gas. The role of the adaption algorithm in resolving various features of blunt body/wake flow for planetary entry conditions is emphasized.

Gnoffo, P. A.↗

NASA NTRS1982

Utilization of normal mode initial conditions for detecting errors in the dynamics part of primitive equation global models

When a global atmospheric basic state has constant angular velocity and its temperature varies with altitude only, there exist normal mode solutions to the linearized global primitive equations. The use of these normal modes, which have known behavior in time, is superior to the use of the Rossby-Haurwitz wave as initial conditions for detecting errors in the dynamics part of primitive equation global models. With these initial conditions, integration through only one time step is sufficient to detect many formulation and coding errors. Other tests are still required for detecting problems of nonlinear instability and conservation of integral properties, however.

Chao, W. C.↗

NASA NTRS1982

A vectorized, finite-volume, adaptive-grid algorithm for Navier-Stokes calculations

An adaptive grid, finite-volume method has been used to solve the Navier-Stokes equations for complete (forebody and afterbody) flowfields around blunt bodies. The code, which is applicable for axisymmetric or two-dimensional flows, allows the mesh to adjust during the computation to provide a closer spacing of mesh points in regions of high gradients, thus minimizing the number of required computational points. The solution technique is explicit, utilizing a maximum time-step advancement at each grid point to accelerate convergence to the steady state. The code has been fully vectorized for efficient solution on the CYBER 203 computer. A very flexible rezoning routine is used to concentrate mesh points anywhere in the field, either by a user-defined weighting function or by allowing high gradient regions to adjust the grid. The grid adjustment routine is implicit in nature and represents a very small portion of the total computational cost. Currently, the code runs in approximately 0.000016 seconds per grid point per iteration.

Gnoffo, P. A.↗

NASA NTRS1982

On the solution of creep induced buckling in general structure

This paper considers the pre and post buckling behavior of general structures exposed to high temperature fields for long durations wherein creep effects become significant. The solution to this problem is made possible through the use of closed upper bounding constraint surfaces which enable the development of a new time stepping algorithm. This permits the stable and efficient solution of structural problems which exhibit indefinite tangent properties. Due to the manner of constraining/bounding successive iterates, the algorithm developed herein is largely self adaptive, inherently stable, sufficiently flexible to handle geometric material and boundary induced nonlinearity, and can be incorporated into either finite element or difference simulations. To illustrate the capability of the procedure, as well as, the physics of creep induced pre and post buckling behavior, the results of several numerical experiments are included.

Padovan, J.↗

NASA NTRS1982

Application of two-point difference schemes to the conservative Euler equations for one-dimensional flows

An implicit finite-difference method is presented for obtaining steady-state solutions to the time-dependent, conservative Euler equations for flows containing shocks. The method uses a two-point central difference scheme with dissipation added at supersonic points via the retarded density concept. Application of the method to the one-dimensional nozzle flow equations for various combinations of subsonic and supersonic boundary conditions show the method to be very efficient. Residuals are typically reduced to machine zero in approximately 35 time steps for 50 mesh points. It is shown that the scheme offers certain advantages over the more widely-used three-point schemes, especially in regard to application of boundary conditions.

Wornom, S. F.↗

NASA NTRS1983

Finite element formulation for transient heat treat problems

The macrothermomechanical behavior of materials subjected to rapid thermal or mechanical loading such as occurs in most heat treatments is described. The equations are developed for Lagrangian, Eulerian, and intermediary kinematic descriptions and are independent of the constitutive laws and the equation of state; they can be solved numerically for a specified material and boundary conditions. The coupled transport effects between dissipation and energy are included. The conventional linearized stability approach indicates the numerical procedure to be stable, with certain restriction on the time step size.

Mullen, R. L.↗

NASA NTRS1983

A system study for the application of microcomputers to research flight test techniques

The onboard simulator is a three degree of freedom aircraft behavior simulator which provides parameters used by the interception procedure. These parameters can be used for verifying closed loop performance before flight. The air to air intercept mode is a software package integrated in the simulation process that generates a target motion and performs a tracking procedure that predicts the most likely next target position, for a defined time step. This procedure also updates relative position parameters and gives adequate fire commands. A microcomputer based on an aircraft spin warning system periodically samples the assymetric thrust and yaw rate of an airplane and then issues voice synthesized warnings and /or suggests to the ilot how to respond to the situation.

Smyth, R. K.↗

NASA NTRS1983

Efficient computerized model for dynamic analysis of energy conversion systems

In searching for the optimum parameters that minimize the total life cycle cost of an energy conversion system, various combinations of components are examined and the resulting system performance and associated economics are studied. The systems performance and economics simulation computer program (SPECS) was developed to fill this need. The program simulates the fluid flow, thermal, and electrical characteristics of a system of components on a quasi-steady state basis for a variety of energy conversion systems. A unique approach is used in which the set of characteristic equations is solved by the Newton-Raphson technique. This approach eliminates the tedious iterative loops which are found in comparable programs such as TRNSYS or SOLTES-1. Several efficient features were also incorporated such as the centralized control and energy management scheme, and analogous treatment of energy flow in electrical and mechanical components, and the modeling of components of similar fundamental characteristics using generic subroutines. Initial tests indicate that this model can be used effectively with a relatively small number of time steps and low computer cost.

Hughes, R. D.↗

NASA NTRS1982

Iterative spectral methods and spectral solutions to compressible flows

A spectral multigrid scheme is described which can solve pseudospectral discretizations of self-adjoint elliptic problems in O(N log N) operations. An iterative technique for efficiently implementing semi-implicit time-stepping for pseudospectral discretizations of Navier-Stokes equations is discussed. This approach can handle variable coefficient terms in an effective manner. Pseudospectral solutions of compressible flow problems are presented. These include one dimensional problems and two dimensional Euler solutions. Results are given both for shock-capturing approaches and for shock-fitting ones.

Hussaini, M. Y.↗

NASA NTRS1982

An implicit formulation for the one-dimensional two-phase multi-interface Stefan Problem

In recent years there has been an intensified effort in the analysis of the Stefan Problem. Particular interest is directed towards determining the location of the interface between the solid and liquid phases of the substance being analyzed. A method has been developed using an implicit finite difference scheme which provides an accurate calculation of the phase change interface locations for a problem with simultaneous multiple phases. This is achieved through an iterative technique which requires the phase change energy to be correctly distributed within the nodes which are undergoing a phase transition. The new method predicts the interface locations without requiring time step size control or node size modification at the interfaces.

Williams, S. D.↗

NASA NTRS1983

An enhanced version of an implicit code for the Euler equations

A two-dimensional implicit finite-difference code is applied to the inviscid Euler equations to compute transonic flow past airfoils in order to provide well-documented standard test cases for the general user community. The code is an improved version of Steger's 1976 implicit code. Enhancements include the use of up-wind differencing in supersonic regions before shocks and a variable time step to accelerate convergence. An airfoil grid generation routine based on algebraic techniques is employed. The grids are clustered near shocks to improve resolution. Computed results are compared with other numerical results from the literature.

Pulliam, T. H.↗

NASA NTRS1983

Numerical simulation of the Weddell Sea pack ice

In order to determine the degree to which a coupled dynamic thermodynamic model can reproduce the seasonal cycle of the Antarctic Sea ice in the Weddell Sea, a series of two-year simulations of the Weddell Sea ice pack were carried out. The simulations employed Hibler's (1979) model and used a one-day time step on an 18 x 15 grid with a resolution of 222 km. Daily atmospheric data from 1979 were used to drive the simulations. The simulations yielded a seasonal cycle of ice with maximum and minimum extents close to that observed. Except for portions of the western Weddell, the advance of the ice is found to be primarily thermodynamic in nature, while the rapid decay depends critically on the presence of both leads and lateral ice advection. In early summer the ice motion causes a residual tongue of ice to extend eastward from the Antarctic Peninsula in agreement with observations. Mean ice thicknesses are consistent with observations, as are the mean drift rates of about 5 km/day.

Hibler, W. D., III↗

NASA NTRS1982

Computer modeling of thermoelectric generator performance

Features of the DEGRA 2 computer code for simulating the operations of a spacecraft thermoelectric generator are described. The code models the physical processes occurring during operation. Input variables include the thermoelectric couple geometry and composition, the thermoelectric materials' properties, interfaces and insulation in the thermopile, the heat source characteristics, mission trajectory, and generator electrical requirements. Time steps can be specified and sublimation of the leg and hot shoe is accounted for, as are shorts between legs. Calculations are performed for conduction, Peltier, Thomson, and Joule heating, the cold junction can be adjusted for solar radition, and the legs of the thermoelectric couple are segmented to enhance the approximation accuracy. A trial run covering 18 couple modules yielded data with 0.3% accuracy with regard to test data. The model has been successful with selenide materials, SiGe, and SiN4, with output of all critical operational variables.

Chmielewski, A. B.↗