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,225 records · Page 68

All sources NASA NTRS DOE OSTI

NASA NTRS1980

The role of microturbulence on collisionless reconnection

The linear, non-linear and anomalous transport properties associated with various microinstabilities driven by cross field currents in reconnecting geometries are reviewed. An assessment of their role in collisionless tearing based on analytic theory, computer simulations and experimental evidence, supports the dominant role of lower hybrid waves. The relevance of microturbulence on macroscopic stationary and time dependent models of merging is presented. It is concluded that a fluid-numerical simulation approach that includes (at each space and time step) the effects of anomalous transport in a self consistent manner, similar to the one used for laboratory collisionless shocks, represents the best method for studying and modeling the details of the reconnection process.

Papadopoulos, K.↗

NASA NTRS1979

The solution of the three-dimensional viscous-compressible Navier-Stokes equations on a vector computer

The development of a vectorized computer code for the solution of the three-dimensional viscous-compressible Navier-Stokes equations is described. The code is applied on the CDC STAR-100 vector computer which is capable of achieving high result rates when a high degree of parallelism is present in the computations. The computational technique is an explicit time-split MacCormack predictor-corrector algorithm. Since a large volume of data is processed and virtual memory utilized, a data management scheme based on interleaving is used. The program has been applied to obtain the solution of the laminar supersonic flow about a family of three-dimensional corners. The equations of motion are expressed in a generalized form relative to a uniform rectangular computational domain. The metric coefficient and boundary conditions must be supplied for the corresponding physical domain. For calculations with 30,000 grid points, a computational rate of 0.00015 seconds per grid point per time step is observed.

Smith, R. E.↗

NASA NTRS1979

A fast semi-implicit algorithm for problems of mixed type

Certain physical processes are modeled by partial differential equations which are parabolic over part of the domain and elliptic over the remainder. A family of semi-implicit algorithms which are well suited to initial-boundary value problems of this mixed type is discussed. One important feature of these algorithms is the use of an approximate inverse for the solution of the implicit linear system. A strong error analysis results in an estimate of the total error as a function of approximate inverse error e and time step h.

Frederickson, P. O.↗

NASA NTRS1980

A model of CO-CH4 global transport/chemistry. II - Preliminary calculation by 1 dimensional model

Calculations of the one-dimensional CO-CH4 transport/chemistry model were performed to investigate the sensitivity of the OH vertical profile to various factors. The boundary conditions at the earth's surface and the tropopause are discussed. The results show that the liquid-phase or soil-phase resistance dominates the mass transfer of both CO and CH4 between the atmosphere and the earth's surface. At the tropopause, the CO and CH4 fluxes are set approximately equal to the value at the previous time step in the numerical calculation. Photolysis rate coefficients that increase with altitude modified the OH vertical profile so that the maximum value was at midtroposphere altitudes. Increasing the CO concentration caused a decrease of the OH concentration, and the OH profile was very sensitive to the NO(x) profiles. The calculated OH profiles are compared with those from other investigations and also with some available measurements. Although the available measurements of the OH profiles are limited, the calculated OH profiles agree fairly well with them.

Kitada, T.↗

NASA NTRS1981

Sample data effects of high-pass filters

Four mathematical models of linear first and second order high pass washout filters were analyzed. These models were Euler's Integration, Zero-Order Hold, Bilinear Transformation, and Second-Order Adams-Bashforth Integration. Bode responses for each model at various sample rates were compared to the continuous filter response. Higher sample rates produced Bode responses approaching the continuous response and the Bilinear Transformation model produced the best responses over the frequency spectrum and sample rates. Pole location analysis of each model in the z-plane showed the Bilinear Transformation and Zero-Order Hold models gave stable poles regardless of time step size, whereas the other models did not display stable poles. A near constant gain error over the entire frequency spectrum was discovered in the Zero-Order Hold cases and a correction gain was calculated for the first-order high-pass filter case.

Chin, D. O.↗

NASA NTRS1980

Time dependent deformation and stress in the lithosphere

Efficient computer programs incorporating frontal solution and time stepping procedure were developed for the modelling of geodynamic problems. This scheme allows for investigating the quasi static phenomena including the effects of the rheological structure of a tectonically active region. From three dimensional models of strike slip earthquakes, it was found that lateral variation of viscosity affects the characteristics of surface deformations. The vertical deformation is especially informative about the viscosity structure in a strike slip fault zone. A three dimensional viscoelastic model of a thrust earthquake indicated that the transient disturbance on plate velocity due to a great plate boundary earthquake is significant at intermediate distances, but becomes barely measurable 1000 km away from the source.

Yang, M.↗

NASA NTRS1981

A Nonlinear Propulsion System Simulation Technique for Piloted Simulators

In the past, propulsion system simulations used in flight simulators have been extremely simple. This resulted in a loss of simulation realism since significant engine and aircraft interactions were neglected and important internal engine parameters were not computed. More detailed propulsion system simulators are needed to permit evaluations of modern aircraft propulsion systems in a simulated flight environment. A real time digital simulation technique has been developed which provides the capabilities needed to evaluate propulsion system performance and aircraft system interaction on manned flight simulators. A parameter correlation technique is used with real and pseudo dynamics in a stable integration convergence loop. The technique has been applied to a multivariable propulsion system for use in a piloted NASA flight simulator program. Cycle time is 2.0 ms on a Univac 1110 computer and 5.7 ms on the simulator computer, a Xerox Sigma 8. The model is stable and accurate with time steps up to 50 ms. The program evaluated the simulation technique and the propulsion system digital control. The simulation technique and model used in that program are described and results from the simulation are presented.

Mihaloew, J. R.↗

NASA NTRS1981

Approximate factorization for incompressible flow

For computational solution of the incompressible Navier-Stokes equations, the approximate factorization (AF) algorithm is used to solve the vectorized momentum equation in delta form based on the pressure calculated in the previous time step. The newly calculated velocities are substituted into the pressure equation (obtained from a linear combination of the continuity and momentum equation), which is then solved by means of line SOR. Computational results are presented for the NACA 66 sub 3 018 airfoil at Reynolds numbers of 1000 and 40,000 and attack angles of 0 and 6 degrees. Comparison with wind tunnel data for Re = 40,000 indicates good qualitative agreement between measured and calculated pressure distributions. Quantitative agreement is only fair, however, with the calculations somewhat displaced from the measurements. Furthermore, the computed velocity profiles are unrealistically thick around the airfoil, due to the excessive amount of artificial viscosity needed for stability. Based on the performance of the algorithm with regard to stability, it is concluded that AF/SOR is suitable for calculations at Reynolds numbers less than 10,000. Speedwise, the method is faster than point SOR by at least a factor of two.

Bernard, R. S.↗

NASA NTRS1981

Some aspects of algorithm performance and modeling in transient 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 algorithms with variable time steps, known as the GEAR package, is described. Four test problems, used for evaluating and comparing various algorithms, were 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 wing of the space shuttle orbiter. 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 NTRS1980

A clear-sky spectral solar radiation model for snow-covered mountainous terrain

A clear-sky spectral solar radiation model for direct and diffuse fluxes, combined with topographic calculations from digital terrain data, computes either incident, net, or reflected solar radiation at any point on a snow surface in mountainous terrain. The radiation may be integrated over any wavelength range from 250 to 5000 nm, or over any time step. Atmospheric attenuation parameters are ozone, water vapor, the Angstrom turbidity coefficient and exponent, and the absorptance to reflectance ratio of the atmospheric aerosols. The model derives these, from measurements which may contain both systematic and random errors, by finding the least squares solution to an overdetermined set of nonlinear equations. For calculations over a specified area, it employs table look-up procedures, so that computation speed for the spectral model approaches that for a lumped model. Thus it may be useful as part of a snow surface energy budget calculation over a drainage basin.

Dozier, J.↗

NASA NTRS1980

The application of Seasat-1 radar altimetry to continental shelf circulation modeling

Seasat-1 radar altimetry data have been used to verify a continental shelf circulation model subjected to tidal forcing on the open boundary. The model makes use of the semi-implicit mode of time integration, removing the surface gravity wave time step restriction. Both the altimetry and the model predictions are compared with empirically determined tidal fluctuations and generally good agreement is obtained.

Cornillon, P.↗

NASA NTRS1980

On time discretizations for spectral methods

New methods are introduced for the time integration of the Fourier and Chebyshev methods of solution for dynamic differential equations. These methods are unconditionally stable, even though no matrix inversions are required. Time steps are chosen by accuracy requirements alone. For the Fourier method both leapfrog and Runge-Kutta methods are considered. For the Chebyshev method only Runge-Kutta schemes are tested. Numerical calculations are presented to verify the analytic results. Applications to the shallow water equations are presented.

Gottlieb, D.↗

NASA NTRS1981

Optimal time splitting for two- and three-dimensional Navier-Stokes equations with mixed derivatives

A new explicit, time splitting algorithm has been developed for finite difference modelling of the full two and three-dimensional time-dependent, compressible, viscous Navier-Stokes equations of fluid mechanics. The scheme is optimal in the sense that the split operators achieve their maximum allowable time step, i.e., the corresponding Courant number. The algorithm allows a conservation-form formulation. Stability is proven analytically and verified numerically. In proving stability it was shown that all nine matrix coefficients of the Navier-Stokes equations are simultaneously symmetrizable by a similarity transformation. Two such transformations and their resulting symmetric matrix coefficients are presented explicitly.

Abarbanel, S.↗

NASA NTRS1981

Numerical method for gas dynamics combining characteristic and conservation concepts

An efficient implicit numerical method that solves the compressible Navier-Stokes equations in arbitrary curvilinear coordinates by the finite-volume technique is presented. An intrinsically dissipative difference scheme and a fully implicit treatment of boundary conditions, based on characteristic and conservation concepts, are used to improve stability and accuracy. Efficiency is achieved by using a diagonal form of the implicit algorithm and spatially varying time-steps. Comparisons of various schemes and methods are presented for one- and two-dimensional flows, including transonic separated flow past a thick circular-arc airfoil in a channel. The new method is equal to or better than a version of MacCormack's hybrid method in accuracy and it converges to a steady state up to an order of magnitude faster.

Coakley, T. J.↗

NASA NTRS1980

Influence of rotor blade tip shape on tip vortex shedding - An unsteady, inviscid analysis

An unsteady potential flow panel method is described based on a time-stepping procedure and using planar quadrilateral panels to represent the surface of thick blades. Each panel has a constant source and doublet distribution and a central control point where an internal Dirichlet boundary condition is applied. The method includes detail paneling around the tip edge and offers the facility of prescribing a tip-edge separation. Preliminary results are shown for four blade tip shapes which are treated as semi-span wings oscillating in pitch about a mean angle of attack of 6 deg. Further work is planned for the wake model before proceeding to higher angles of attack.

Maskew, B.↗

NASA NTRS1981

Numerical stability of an explicit finite difference scheme for the solution of transient conduction in composite media

A theoretical evaluation of the stability of an explicit finite difference solution of the transient temperature field in a composite medium is presented. The grid points of the field are assumed uniformly spaced, and media interfaces are either vertical or horizontal and pass through grid points. In addition, perfect contact between different media (infinite interfacial conductance) is assumed. A finite difference form of the conduction equation is not valid at media interfaces; therefore, heat balance forms are derived. These equations were subjected to stability analysis, and a computer graphics code was developed that permitted determination of a maximum time step for a given grid spacing.

Campbell, W.↗

NASA NTRS1981

Development of an engineering analysis of progressive damage in composites during low velocity impact

A computerized, analytical methodology was developed to study damage accumulation during low velocity lateral impact of layered composite plates. The impact event was modeled as perfectly plastic with complete momentum transfer to the plate structure. A transient dynamic finite element approach was selected to predict the displacement time response of the plate structure. Composite ply and interlaminar stresses were computed at selected time intervals and subsequently evaluated to predict layer and interlaminar damage. The effects of damage on elemental stiffness were then incorporated back into the analysis for subsequent time steps. Damage predicted included fiber failure, matrix ply failure and interlaminar delamination.

Humphreys, E. A.↗

NASA NTRS1982

A numerical method for integrating the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes

An extension of the method of moments is developed for the numerical integration of the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes. The number density function n sub k (x,t) in each separate droplet packet between droplet mass grid points (x sub k, x sub k+1) is represented by an expansion in orthogonal polynomials with a given weighting function. In this way droplet number concentrations, liquid water contents and other moments in each droplet packet are conserved and the problem of solving the kinetic equations is replaced by one of solving a set of coupled differential equations for the number density function moments. The method is tested against analytic solutions of the corresponding kinetic equations. Numerical results are obtained for different coalescence/breakup and condensation/evaporation kernels and for different initial droplet spectra. Also droplet mass grid intervals, weighting functions, and time steps are varied.

Emukashvily, I. M.↗