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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Implicit finite-difference simulations of steady and unsteady transonic flows

Implicit methods for several fluid dynamic formulations have been developed and applied to steady-state and low-frequency transonic flows. The basic steps involved in the construction of implicit schemes include: selection of linearly stable accurate implicit difference operators, time-linearization of nonlinear terms, and approximate factorization of the implicit operators into easily solved systems of equations. The proposed schemes are found very efficient for the simpler formulations.

Ballhaus, W. F.↗
NASA NTRS1982

Transient phenomena of shock-induced turbulent separation for a spikebody and stalling airfoil at transonic and supersonic speeds

The time-dependent, compressible, Reynolds-averaged, full Navier-Stokes equations are applied to solve an axisymmetric flow around a forward-facing stepbody (spikebody) at supersonic speeds and a stalling airfoil at transonic speeds. Important transient and unsteady phenomena, not yet well understood, are examined, and significant new findings of the present solution to the phenomena are discussed. The phenomena described in detail are as follows: The evolution of the shock wave pressure built up by the impact of the pressure waves, one from the trailing edge; the separation of the flow as influenced by the shock wave; the location of the reversed flow, the separation point, and the reattachment point; and the transient (or unsteady) phenomena of the flow pulsation, oscillation, and stalling of the body and airfoil wake flow. The numerical results show that the transient flow instability is caused by a supersonic jet induced in the separation bubble by the shock-bifurcation (lambda shock) mechanism between the separation shock and the reflected shock. Pulsation and stall phenomena are caused by a sudden increase in the leading-edge pressure due to the jet and the separation bubble interacting along the stagnation point flow.

Yoshikawa, K. K.↗
NASA NTRS1985

The science of computing - The evolution of parallel processing

The present paper is concerned with the approaches to be employed to overcome the set of limitations in software technology which impedes currently an effective use of parallel hardware technology. The process required to solve the arising problems is found to involve four different stages. At the present time, Stage One is nearly finished, while Stage Two is under way. Tentative explorations are beginning on Stage Three, and Stage Four is more distant. In Stage One, parallelism is introduced into the hardware of a single computer, which consists of one or more processors, a main storage system, a secondary storage system, and various peripheral devices. In Stage Two, parallel execution of cooperating programs on different machines becomes explicit, while in Stage Three, new languages will make parallelism implicit. In Stage Four, there will be very high level user interfaces capable of interacting with scientists at the same level of abstraction as scientists do with each other.

Denning, P. J.↗
NASA NTRS1986

On the instability of a three-dimensional attachment-line boundary layer - Weakly nonlinear theory and a numerical approach

The instability of a three-dimensional attachment-line boundary layer is considered in the nonlinear regime. Using weakly nonlinear theory, it is found that, apart from a small interval near the (linear) critical Reynolds number, finite-amplitude solutions bifurcate subcritically from the upper branch of the neutral curve. The time-dependent Navier-Stokes equations for the attachment-line flow have been solved using a Fourier-Chebyshev spectral method and the subcritical instability is found at wavenumbers that correspond to the upper branch. Both the theory and the numerical calculations show the existence of supercritical finite-amplitude (equilibrium) states near the lower branch which explains why the observed flow exhibits a preference for the lower branch modes. The effect of blowing and suction on nonlinear stability of the attachment-line boundary layer is also investigated.

Hall, P.↗
NASA NTRS1987

Finite-volume Euler and Navier-Stokes solvers for three-dimensional and conical vortex flows over delta wings

A unified central-difference finite-volume Euler and Navier-Stokes solver with four-stage Runge-Kutta time stepping is presented. The computer code developed for this purpose is capable of solving the standard set and nonstandard sets (zero-total-pressure loss) of Euler equations and the thin-layer and full Navier-Stokes equations. Applications are presented for conical supersonic flows with weak shocks using the standard and nonstandard sets of Euler equations, and the thin-layer and full Navier-Stokes equations for sharp and round-edged delta wings. Applications are also presented for three-dimensional transonic and subsonic flows using the standard set of Euler equations for sharp-edged delta wings. The computational results of the different sets of equations are compared with each other and with the experimental results and conclusions on the validity of these sets to these applications, are presented.

Kandil, Osama A.↗
NASA NTRS1988

A diagonally inverted LU implicit multigrid scheme

A new Diagonally Inverted LU Implicit scheme is developed within the framework of the multigrid method for the 3-D unsteady Euler equations. The matrix systems that are to be inverted in the LU scheme are treated by local diagonalizing transformations that decouple them into systems of scalar equations. Unlike the Diagonalized ADI method, the time accuracy of the LU scheme is not reduced since the diagonalization procedure does not destroy time conservation. Even more importantly, this diagonalization significantly reduces the computational effort required to solve the LU approximation and therefore transforms it into a more efficient method of numerically solving the 3-D Euler equations.

Yokota, Jeffrey W.↗
NASA NTRS1988

Thermo-viscoelastic analysis of composite materials, volume 1

Advanced composite materials, especially graphite/epoxy, are being applied to aircraft structures in order to improve performance and save weight. An important consideration in composite design is the residual strength of a structure containing holes, delaminations, or interlaminar damage when subjected to compressive loads. Recent studies have revealed the importance of viscoelastic effects in polymer-based composites. The viscoelastic effect is particularly significant at elevated temperature/moisture conditions since the matrix material is strongly affected by the environment. The solution of viscoelastic problems in composites was limited to special cases which can be solved by classical lamination theory. A finite element procedure is presented for calculating time-dependent stresses and strains in composite structures with general configurations and complicated boundary conditions. Using this procedure the in-plane and interlaminar stress distributions and histories in notched and unnotched composites were obtained for mechanical and thermal loads. Both two-dimensional and three-dimensional viscoelastic problems are analyzed. The effects of layup orientation and load spectrum on creep response and stress relaxation were also studied.

Lin, K. Y.↗
NASA NTRS1988

A diagonally inverted LU implicit multigrid scheme

A new Diagonally Inverted LU Implicit scheme is developed within the framework of the multigrid method for the three-dimensional unsteady Euler equations. The matrix systems that are to be inverted in the LU scheme are treated by local diagonalizing transformations that decouple them into systems of scalar equations. Unlike the Diagonalized ADI method, the time accuracy of the LU scheme is not reduced since the diagonalizing procedure does not destroy time conservation. Even more importantly, this diagonalization significantly reduces the computational effort required to solve the LU approximation and therefore transforms it into a more efficient method of numerically solving the three-dimensional Euler equations.

Yokota, J. W.↗
NASA NTRS1990

Numerical Methods For Chemically Reacting Flows

Issues related to numerical stability, accuracy, and resolution discussed. Technical memorandum presents issues in numerical solution of hyperbolic conservation laws containing "stiff" (relatively large and rapidly changing) source terms. Such equations often used to represent chemically reacting flows. Usually solved by finite-difference numerical methods. Source terms generally necessitate use of small time and/or space steps to obtain sufficient resolution, especially at discontinuities, where incorrect mathematical modeling results in unphysical solutions.

Leveque, R. J.↗
NASA NTRS1989

A current disruption mechanism in the neutral sheet for triggering substorm expansions

Two main areas were addressed in support of an effort to understand mechanism responsible for the broadband electrostatic noise (BEN) observed in the magnetotail. The first area concerns the generation of BEN in the boundary layer region of the magnetotail whereas the second area concerns the occassional presence of BEN in the neutral sheet region. For the generation of BEN in the boundary layer region, a hybrid simulation code was developed to perform reliable longtime, quiet, highly resolved simulations of field aligned electron and ion beam flow. The result of the simulation shows that broadband emissions cannot be generated by beam-plasma instability if realistic values of the ion beam parameters are used. The waves generated from beam-plasma instability are highly discrete and are of high frequencies. For the plasma sheet boundary layer condition, the wave frequencies are in the kHz range, which is incompatible with the observation that the peak power in BEN occur in the 10's of Hz range. It was found that the BEN characteristics are more consistent with lower hybrid drift instability. For the occasional presence of BEN in the neutral sheet region, a linear analysis of the kinetic cross-field streaming instability appropriate to the neutral sheet condition just prior to onset of substorm expansion was performed. By solving numerically the dispersion relation, it was found that the instability has a growth time comparable to the onset time scale of substorm onset. The excited waves have a mixed polarization in the lower hybrid frequency range. The imposed drift driving the instability corresponds to unmagnetized ions undergoing current sheet acceleration in the presence of a cross-tail electric field. The required electric field strength is in the 10 mV/m range which is well within the observed electric field values detected in the neutral sheet during substorms. This finding can potentially account for the disruption of cross-tail current and its diversion to the ionosphere to form the substorm current wedge. Furthermore, a number of features associated with substorm expansion onset can be understood based on this substorm onset scenario.

Lui, A. T. Y.↗
NASA NTRS1989

Nonlinear interactions in mixing layers and compressible heated round jets

The nonlinear interactions between a fundamental instability mode and both its harmonics and the changing mean flow are studied using the weakly nonlinear stability theory of Stuart and Watson, and numerical solutions of coupled nonlinear partial differential equations. The first part focuses on incompressible cold (or isothermal; constant temperature throughout) mixing layers, and for these, the first and second Landau constants are calculated as functions of wavenumber and Reynolds number. It is found that the dominant contribution to the Landau constants arises from the mean flow changes and not from the higher harmonics. In order to establish the range of validity of the weakly nonlinear theory, the weakly nonlinear and numerical solutions are compared and the limitation of each is discussed. At small amplitudes and at low-to-moderate Reynolds numbers, the two results compare well in describing the saturation of the fundamental, the distortion of the mean flow, and the initial stages of vorticity roll-up. At larger amplitudes, the interaction between the fundamental, second harmonic, and the mean flow is strongly nonlinear and the numerical solution predicts flow oscillations, whereas the weakly nonlinear theory yields saturation. In the second part, the weakly nonlinear theory is extended to heated (or nonisothermal; mean temperature distribution) subsonic round jets where quadratic and cubic nonlinear interactions are present, and the Landau constants also depend on jet temperature ratio, Mach number and azimuthal mode number. Under exponential growth and nonlinear saturation, it is found that heating and compressibility suppress the growth of instability waves, that the first azimuthal mode is the dominant instability mode, and that the weakly nonlinear solution describes the early stages of the roll-up of an axisymmetric shear layer. The receptivity of a typical jet flow to pulse type input disturbance is also studied by solving the initial value problem and then examining the behavior of the long-time solution.

Jarrah, Yousef Mohd↗
NASA NTRS1990

Coupled BE/FE/BE approach for scattering from fluid-filled structures

NASHUA is a coupled finite element/boundary element capability built around NASTRAN for calculating the low frequency far-field acoustic pressure field radiated or scattered by an arbitrary, submerged, three-dimensional, elastic structure subjected to either internal time-harmonic mechanical loads or external time-harmonic incident loadings. Described here are the formulation and use of NASHUA for solving such structural acoustics problems when the structure is fluid-filled. NASTRAN is used to generate the structural finite element model and to perform most of the required matrix operations. Both fluid domains are modeled using the boundary element capability in NASHUA, whose matrix formulation (and the associated NASTRAN DMAP) for evacuated structures can be used with suitable interpretation of the matrix definitions. After computing surface pressures and normal velocities, far-field pressures are evaluated using an asymptotic form of the Helmholtz exterior integral equation. The proposed numerical approach is validated by comparing the acoustic field scattered from a submerged fluid-filled spherical thin shell to that obtained with a series solution, which is also derived here.

Everstine, Gordon C.↗
NASA NTRS1990

Numerical simulation of 3-D Benard convection with gravitational modulation

In this numerical study, randomly and sinusoidally modulated gravitational fields imposed on three-dimensional Rayleigh-Benard convection are investigated in an effort to understand the effects of vibration (G-Jitter) on fluid systems. The time-dependent, Navier-Stokes equations and the energy equation with Boussinesq approximations are solved by a semi-implicit, pseudospectral procedure. An analysis of energy balances indicates that with increasing modulation amplitude, transition from synchronous to relaxation oscillation goes through the subharmonic response. Random modulations are found to be less stabilizing than sinusoidal and are shown to impose three-dimensionality on the flow for some parameter ranges both at terrestrial and zero base gravity conditions.

Biringen, S.↗
NASA NTRS1991

Alternative smoothing algorithms for on-line estimation problems

In recent years, much emphasis has been placed on the Kalman filter to solve on-line estimation problems. Smoothers, which determine state estimates at a certain time t using measurements between 0 and T, have been used primarily in off-line estimation schemes. This paper presents alternative forms of smoothers that can be used for on-line estimation problems. The alternative smoothing algorithms provide state estimates which are more refined than Kalman filtered state estimates. The alternative smoothing algorithms are illustrated in two examples. The first example considers a linear discrete-time system while the second example considers a nonlinear continuous-time system with discrete-time measurements.

Baumgartner, Eric T.↗
NASA NTRS1991

Non-linear dynamics of a geared rotor-bearing system with multiple clearances

The nonlinear frequency response characteristics of a geared rotor-bearing system are examined. A three-DOF dynamic model is developed which includes nonlinearities associated with radial clearances in the radial rolling element bearings and backlash between a spur gear pair; linear time-invariant gear meshing stiffness is assumed. The corresponding linear system problem is also solved, and predicted natural frequencies and modes match with finite element method results. The applicability of both analytical and numerical solution techniques to the multi-DOF nonlinear problem is investigated. Satisfactory agreement is found between the theory presented here and the available experimental data. Nonlinear modal interactions and differences between internal static transmission error excitation and external torque excitation are discussed. In addition, parametric studies are performed to understand the effect of system parameters such as bearing stiffness to gear mesh stiffness ratio, alternating to mean force ratio and radial bearing preload to mean force ratio on the nonlinear dynamic behavior.

Kahraman, A.↗
NASA NTRS1991

Response of a confined gas to volumetric heating in the absence of gravity. I - Slow transients

A one-dimensional model for bulk motion induced by a transient volumetric heat source in a confined gas at zero gravity is considered. Rational approximation methods are used to derive a quantitative theory for the gas response to a spatiallly distributed, time-dependent internal power deposition. The resulting low Mach number compressible flow equations are solved by using perturbation methods. Solutions are given for a conduction-free core and thin conductive boundary layers adjacent to the end walls. It is found that may any spatially nonuniform power deposition will cause fluid motion. Net mass transport in the closed container will occur for certain spatially distributed heating. The model mimics the thermal effects of an exothermic gas phase reaction in vapor transport experiments conducted in space. The solutions demonstrate that thermally induced mass transport can be as large as diffusive mass transport in a typical experiment.

Herczynski, A.↗
NASA NTRS1990

Analysis of unsteady compressible viscous layers

The development of an analysis to predict the unsteady compressible flows in blade boundary layers and wakes is presented. The equations that govern the flows in these regions are transformed using an unsteady turbulent generalization of the Levy-Lees transformation. The transformed equations are solved using a finite difference technique in which the solution proceeds by marching in time and in the streamwise direction. Both laminar and turbulent flows are studied, the latter using algebraic turbulence and transition models. Laminar solutions for a flat plate are shown to approach classical asymptotic results for both high and low frequency unsteady motions. Turbulent flat-plate results are in qualitative agreement with previous predictions and measurements. Finally, the numerical technique is also applied to the stator and rotor of a low-speed turbine stage to determine unsteady effects on surface heating. The results compare reasonably well with measured heat transfer data and indicate that nonlinear effects have minimal impact on the mean and unsteady components of the flow.

Power, G. D.↗