Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Positivity-preserving numerical schemes for multidimensional advection

This report describes the construction of an explicit, single time-step, conservative, finite-volume method for multidimensional advective flow, based on a uniformly third-order polynomial interpolation algorithm (UTOPIA). Particular attention is paid to the problem of flow-to-grid angle-dependent, anisotropic distortion typical of one-dimensional schemes used component-wise. The third-order multidimensional scheme automatically includes certain cross-difference terms that guarantee good isotropy (and stability). However, above first-order, polynomial-based advection schemes do not preserve positivity (the multidimensional analogue of monotonicity). For this reason, a multidimensional generalization of the first author's universal flux-limiter is sought. This is a very challenging problem. A simple flux-limiter can be found; but this introduces strong anisotropic distortion. A more sophisticated technique, limiting part of the flux and then restoring the isotropy-maintaining cross-terms afterwards, gives more satisfactory results. Test cases are confined to two dimensions; three-dimensional extensions are briefly discussed.

Leonard, B. P.↗
NASA NTRS1992

Investigation of forced unsteady separated flows using velocity-vorticity form of Navier-Stokes equations

The phenomenon of forced unsteady separation and eruption of boundary-layer vorticity is a highly-complex, high-Reynolds number flow phenomenon, which abruptly leads to the formation of a dynamic stall vortex as demonstrated earlier by the authors for a NACA 0015 airfoil undergoing constant rate pitch-up motion. This, as well as the results of other researchers, have convincingly demonstrated a complex vortical structure within the state of unsteady separation prior to the evolution of dynamic stall. This phenomenon of vortex eruption, although observed in studying dynamic stall phenomena, is also associated with transition from laminar to turbulence flow and its generic nature has been stressed by many researchers including the present investigators. An unsteady Navier-Stokes (NS) analysis is developed for arbitrarily maneuvering bodies using velocity-vorticity variables; this formulation is nearly form-invariant under a generalized non-inertial coordinate transformation. A fully-implicit uniformly second-order accurate method is used, with the nonlinear convective terms approximated using a biased third-order upwind differencing scheme to be able to simulate higher-Re flows. No explicit artificial dissipation is added. The numerical method is fully vectorized and currently achieves a computational index of 7 micro-seconds per time step per mesh point, using a single processor on a CRAY Y-MP. The simulation results show that the energetic free shear from the leading edge is responsible for the wall viscous layer to abruptly erupt near the center of the counterclockwise rotating eddy in the unsteady boundary layer. Primary, secondary, tertiary and quaternary vortices have been observed before the dynamic stall vortex evolves and gathers its maximum strength. This study will discuss the simulation results of Reynolds number up to Re = 45,000 and will also discuss the efforts of initial acceleration in a specific maneuver, on the evolution of the stall vortex.

Ghia, K. N.↗
NASA NTRS1993

Development of an unstructured solution adaptive method for the quasi-three-dimensional Euler and Navier-Stokes equations

A general solution adaptive scheme based on a remeshing technique is developed for solving the two-dimensional and quasi-three-dimensional Euler and Favre-averaged Navier-Stokes equations. The numerical scheme is formulated on an unstructured triangular mesh utilizing an edge-based pointer system which defines the edge connectivity of the mesh structure. Jameson's four-stage hybrid Runge-Kutta scheme is used to march the solution in time. The convergence rate is enhanced through the use of local time stepping and implicit residual averaging. As the solution evolves, the mesh is regenerated adaptively using flow field information. Mesh adaptation parameters are evaluated such that an estimated local numerical error is equally distributed over the whole domain. For inviscid flows, the present approach generates a complete unstructured triangular mesh using the advancing front method. For turbulent flows, the approach combines a local highly stretched structured triangular mesh in the boundary layer region with an unstructured mesh in the remaining regions to efficiently resolve the important flow features. One-equation and two-equation turbulence models are incorporated into the present unstructured approach. Results are presented for a wide range of flow problems including two-dimensional multi-element airfoils, two-dimensional cascades, and quasi-three-dimensional cascades. This approach is shown to gain flow resolution in the refined regions while achieving a great reduction in the computational effort and storage requirements since solution points are not wasted in regions where they are not required.

Jiang, Yi-Tsann↗
NASA NTRS1993

Navier-Stokes analysis of three-dimensional flow and heat transfer inside turbine blade rows

A numerical method for solving the three-dimensional, Navier-Stokes equations for unsteady, viscous flow and heat transfer through multiple turbomachinery blade rows is presented. The method solves the fully three-dimensional Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed. A third-order accurate upwinding scheme is used to approximate convection terms while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow and heat transfer field of the High Pressure Fuel side Turbo-Pump (HPFTP) of the Space Shuttle Main Engine (SSME). The stage calculation is performed by coupling the stator and the rotor flow fields at each time step through an over-laid grid. Numerical results for the complete geometry with the vane trailing edge cutback are presented and compared with the available experimental data.

Hah, C.↗
NASA NTRS1992

An effective multigrid method for high-speed flows

The use is considered of a multigrid method with central differencing to solve the Navier-Stokes equations for high speed flows. The time dependent form of the equations is integrated with a Runge-Kutta scheme accelerated by local time stepping and variable coefficient implicit residual smoothing. Of particular importance are the details of the numerical dissipation formulation, especially the switch between the second and fourth difference terms. Solutions are given for 2-D laminar flow over a circular cylinder and a 15 deg compression ramp.

Swanson, R. C.↗
NASA NTRS1993

Discreteness noise versus force errors in N-body simulations

A low accuracy in the force calculation per time step of a few percent for each particle pair is sufficient for collisionless N-body simulations. Higher accuracy is made meaningless by the dominant discreteness noise in the form of two-body relaxation, which can be reduced only by increasing the number of particles. Since an N-body simulation is a Monte Carlo procedure in which each particle-particle force is essentially random, i.e., carries an error of about 1000 percent, the only requirement is a systematic averaging-out of these intrinsic errors. We illustrate these assertions with two specific examples in which individual pairwise forces are deliberately allowed to carry significant errors: tree-codes on supercomputers and algorithms on special-purpose machines with low-precision hardware.

Hernquist, Lars↗
NASA NTRS1993

The effect of Reynolds number and turbulence on airfoil aerodynamics at -90 degrees incidence

A method has been developed for calculating the viscous flow about airfoils in with and without deflected flaps at -90 deg incidence. This method provides for the solution of the unsteady incompressible Navier-Stokes equations by means of an implicit technique. The solution is calculated on a body-fitted computational mesh using a staggered grid method. The vorticity is defined at the node points, and the velocity components are defined at the mesh-cell sides. The staggered-grid orientation provides for accurate representation of vorticity at the node points and the continuity equation at the mesh-cell centers. The method provides for the direct solution of the flow field and satisfies the continuity equation to machine zero at each time-step. The method is evaluated in terms of its ability to predict two-dimensional flow about an airfoil at -90 degrees incidence for varying Reynolds number and different boundary layer models. A laminar and a turbulent boundary layer models. A laminar and a turbulent boundary layer model are considered in the evaluation of the method. The variation of the average loading and surface pressure distribution due to flap deflection, Reynolds number, and laminar or turbulent flow are presented and compared with experimental results. The comparisons indicate that the calculated drag and drag reduction caused by flap deflection and the calculated average surface pressure are in excellent agreement with the measured results at a similar Reynolds number.

Stremel, Paul M.↗
NASA NTRS1993

An installed nacelle design method using multiblock Euler solver

An efficient multiblock Euler design method was developed for designing a nacelle installed on geometrically complex airplane configurations. This approach employed a design module based on a direct iterative surface curvature method developed at NASA-Langley. A multiblock Euler flow solver was used for computing flow around complex geometries. The flow solver used a finite-volume formulation with explicit time-stepping to solve the Euler equations. It used a multiblock version of the multigrid method to accelerate the convergence of the calculations. The design module successively updated the surface geometry to reduce the difference between the computed and target pressure distributions. In the flow solver, the change in surface geometry was simulated by applying surface transpiration boundary conditions to avoid repeated grid generation during design iterations. Smoothness of the designed surface was ensured by alternate application of streamwise and circumferential smoothings.

Chen, H. C.↗
NASA NTRS1993

Acceleration of convergence and spectrum transformation of implicit finite difference operators associated with Navier-Stokes equations

Implicit difference operator spectra are presently computed by applying eigensystem analysis techniques to finite-difference formulations of 2D Euler and Navier-Stokes equations, and attention is given to these iterative methods' convergence and stability characteristics by taking into account the effects of grid geometry, time-step, numerical viscosity, and boundary conditions. On the basis of the eigenvalue distributions for various flow configurations, the feasibility of applying such convergence-acceleration techniques as eigenvalue annihilation and relaxation is discussed. Spectrum-shifting is applied to NASA-Ames' ARC2D flow code, achieving a 20-33 percent efficiency.

Saleem, M.↗
NASA NTRS1993

Direct calculations of waves in fluid flows using a high-order compact difference scheme

The solution of the unsteady Euler equations by a sixth-order compact difference scheme combined with a fourth-order Runge-Kutta method is investigated. Closed-form expressions for the amplification factors and their corresponding dispersion correlations are obtained by Fourier analysis of the fully discretized, two-dimensional Euler equations, and the numerical dissipation, dispersion, and anisotropic effects are assessed. It is found that the CFL limit for stable calculations is about 0.8. For a CFL number equal to 0.6, the smallest wavelength which is resolved without numerical damping is about 6 to 8 grid nodes. For phase speeds corresponding to acoustic waves, the corresponding time period is resolved by about 200 to 300 time steps. Three numerical examples of waves in compressible flow are included.

Yu, Sheng-Tao↗
NASA NTRS1993

Three-dimensional unstructured grid Euler method applied to turbine blades

Flow through a turbine annular cascade is calculated using a three-dimensional Euler method based on unstructured tetrahedral meshes. The equations are integrated in time using an explicit Runge-Kutta time-stepping scheme. The inviscid flux terms are discretized using a cell-centered finite-volume formulation with upwind flux-difference splitting. The tetrahedral meshes around the turbine blade are generated using an advancing-front technique with forced geometric periodicity between the blades. Good agreement is obtained between the present calculation and the experiment for both surface pressure distribution and flow behavior in the passage between the blades, demonstrating the capability of the present methodology for turbomachinery flow applications.

Kwon, Oh J.↗
NASA NTRS1993

An algorithm for solving the perturbed gas dynamic equations

The present application of a compact, higher-order central-difference approximation to the linearized Euler equations illustrates the multimodal character of these equations by means of computations for acoustic, vortical, and entropy waves. Such dissipationless central-difference methods are shown to propagate waves exhibiting excellent phase and amplitude resolution on the basis of relatively large time-steps; they can be applied to wave problems governed by systems of first-order partial differential equations.

Davis, Sanford↗
NASA NTRS1993

Application of an unstructured grid flow solver to planes, trains and automobiles

Rampant, an unstructured flow solver developed at Fluent Inc., is used to compute three-dimensional, viscous, turbulent, compressible flow fields within complex solution domains. Rampant is an explicit, finite-volume flow solver capable of computing flow fields using either triangular (2d) or tetrahedral (3d) unstructured grids. Local time stepping, implicit residual smoothing, and multigrid techniques are used to accelerate the convergence of the explicit scheme. The paper describes the Rampant flow solver and presents flow field solutions about a plane, train, and automobile.

Spragle, Gregory S.↗
NASA NTRS1993

A global multilevel atmospheric model using a vector semi-Lagrangian finite-difference scheme. I - Adiabatic formulation

An adiabatic global multilevel primitive equation model using a two time-level, semi-Lagrangian semi-implicit finite-difference integration scheme is presented. A Lorenz grid is used for vertical discretization and a C grid for the horizontal discretization. The momentum equation is discretized in vector form, thus avoiding problems near the poles. The 3D model equations are reduced by a linear transformation to a set of 2D elliptic equations, whose solution is found by means of an efficient direct solver. The model (with minimal physics) is integrated for 10 days starting from an initialized state derived from real data. A resolution of 16 levels in the vertical is used, with various horizontal resolutions. The model is found to be stable and efficient, and to give realistic output fields. Integrations with time steps of 10 min, 30 min, and 1 h are compared, and the differences are found to be acceptable.

Bates, J. R.↗
NASA NTRS1993

Fast dynamos with finite resistivity in steady flows with stagnation points

Results are presented of a kinematic fast dynamo problem for two classes of steady incompressible flows: the ABC flow and the spatially aperiodic flow of Lau and Finn (1992). The numerical method used to find the solutions is described, together with convergence studies with respect to the time step and the number of points N of the spatial grid. It is shown that the growth rate and frequency can be extrapolated to N = infinity. Results are presented indicating that fast kinematic dynamos can exist in both these flows and that chaotic flow is a necessary condition. It was found that, for the ABC flow with A = B = C, there are two dynamo modes: an oscillating mode and a purely growing mode.

Lau, Yun-Tung↗
NASA NTRS1993

The application of preconditioning in viscous flows

The present time-derivative preconditioning algorithm is effective in flow conditions ranging from inviscid to very diffusive flows, as well as low subsonic to supersonic flow velocities. By means of a preconditioning matrix that (1) introduces well-conditioned eigenvalues and (2) avoids nonphysical time reversals for viscous flows, a mechanism is obtained which controls the inviscid and viscous time-step parameters at very diffusive flows. These capabilities are demonstrated for a variety of sample problems; convergence rates of solutions that are indistinguishable from those obtained without preconditioning are shown to be accelerated by as much as two orders of magnitude.

Choi, Y.-H.↗
NASA NTRS1993

Global/local methods for probabilistic structural analysis

A probabilistic global/local method is proposed to reduce the computational requirements of probabilistic structural analysis. A coarser global model is used for most of the computations with a local more refined model used only at key probabilistic conditions. The global model is used to establish the cumulative distribution function (cdf) and the Most Probable Point (MPP). The local model then uses the predicted MPP to adjust the cdf value. The global/local method is used within the advanced mean value probabilistic algorithm. The local model can be more refined with respect to the g1obal model in terms of finer mesh, smaller time step, tighter tolerances, etc. and can be used with linear or nonlinear models. The basis for this approach is described in terms of the correlation between the global and local models which can be estimated from the global and local MPPs. A numerical example is presented using the NESSUS probabilistic structural analysis program with the finite element method used for the structural modeling. The results clearly indicate a significant computer savings with minimal loss in accuracy.

Millwater, H. R.↗
NASA NTRS1993

High order finite difference and multigrid methods for spatially evolving instability in a planar channel

The fourth-order finite-difference scheme with fully implicit time-marching presently used to computationally study the spatial instability of planar Poiseuille flow incorporates a novel treatment for outflow boundary conditions that renders the buffer area as short as one wavelength. A semicoarsening multigrid method accelerates convergence for the implicit scheme at each time step; a line-distributive relaxation is developed as a robust fast solver that is efficient for anisotropic grids. Computational cost is no greater than that of explicit schemes, and excellent agreement with linear theory is obtained.

Liu, C.↗