Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “implicit methods”

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 973 records · Page 54

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.↗

Derivation and evaluation of an approximate analysis for three-dimensional viscous subsonic flow with large secondary velocities

An approximate analysis is presented for calculating three-dimensional, low Mach number, laminar viscous flows in curved passages with large secondary flows and corner boundary layers. The analysis is based on the decomposition of the overall velocity field into inviscid and viscous components with the overall velocity being determined from superposition. An incompressible vorticity transport equation is used to estimate inviscid secondary flow velocities to be used as corrections to the potential flow velocity field. A parabolized streamwise momentum equation coupled to an adiabatic energy equation and global continuity equation is used to obtain an approximate viscous correction to the pressure and longitudinal velocity fields. A collateral flow assumption is invoked to estimate the viscous correction to the transverse velocity fields. The approximate analysis is solved numerically using an implicit ADI solution for the viscous pressure and velocity fields. An iterative ADI procedure is used to solve for the inviscid secondary vorticity and velocity fields. This method was applied to computing the flow within a turbine vane passage with inlet flow conditions of M = 0.1 and M = 0.25, Re = 1000 and adiabatic walls, and for a constant radius curved rectangular duct with R/D = 12 and 14 and with inlet flow conditions of M = 0.1, Re = 1000, and adiabatic walls.

Anderson, O. L.↗

Improved Euler analysis of advanced turboprop propeller flows

An implicit approximate factorization scheme in conjunction with a new boundary treatment was used to compute the inviscid flow field about an advanced high-speed propeller. The method of characteristics was used to apply impermeable boundary conditions. The convergence history of numerical calculations shows substantial decrease in the residual error decay. Euler solutions were computed for SR-3 propeller geometry and the results were compared with the experimental data and previous numerical results.

Yamamoto, O.↗

Fatigue life prediction of an intermetallic matrix composite at elevated temperatures

A strain-based fatigue life prediction method is proposed for an intermetallic matrix composite (IMC) under tensile cyclic loadings at elevated temperatures. Styled after the 'Universal Slopes' method, the model utilizes the composite's tensile properties to estimate fatigue life. Factors such as fiber volume ratio (Vf), number of plys and temperature dependence are implicitly incorporated into the model through these properties. The model constants are determined by using unidirectional fatigue data at temperatures of 425 and 815 C. Fatigue lives from two independent sources are used to verify the model at temperatures of 650 and 760 C. Cross-ply lives at 760 C are also predicted. It is demonstrated that the correlation between experimental and predicted lives is within a factor of two.

Bartolotta, P. A.↗

Some considerations on velocity vector accuracy in dust trajectory analysis

The relative contributions of comets and asteroids to the reservoir of dust in the interplanetary medium is not known. There are direct observations of dust released from comets and there is evidence to associate the IRAS dust bands with possible collisions of asteroids in the main belt. A means towards sorting out the parent sources has been proposed in the establishment of a dust collector in orbit about the Earth. The purpose of such a facility would be to collect not only cosmic dust particles intact but also the state vectors, as they arrive at the detector, the idea being that one may combine analytical laboratory analysis of the physics and chemistry of the captured particles with orbital data in order to help distinguish between bodies and identify parent bodies. The theoretical study of dust particle orbits in the solar system takes on greatly more importance if we use collected trajectory data. The orbital motion of dust when radiation and forces alone are acting is well understood. When gravitational forces due to the planets are included, the motion can become quite complex. In order to characterize the orbits of particles as they crossed the Earth's orbits, a study of the long-time dust orbital evolution was undertaken. We have considered various parameters associated with these dust orbits to see if one may in a general way discriminate between particles evolved from comets and asteroids. We proceed in this study as we have done previously. That is, we considered the dust particles as ideal black bodies, of density 1 gm/cc, spherical, with radii 10-100 microns. Particles of this size are affected by radiation forces, photon pressure, and Poynting-Robertson drag. Account was also taken of solar wind drag, which amounts to about 30 percent of the Poynting-Robertson drag negligible. The gravitational forces due to the planets are included, unlike in our previous study; the planetary orbits are those of true n-body interaction so that the possibility of secular resonance is included. Our method was to calculate explicitly by a numerical procedure the orbits of dust particles after they left their parent bodies. The motion is determined numerically with the implicit Runge-Kutta integrator using Gauss-Radau spacings.

Jackson, A. A.↗

Ablative Thermal Response Analysis Using the Finite Element Method

A review of the classic techniques used to solve ablative thermal response problems is presented. The advantages and disadvantages of both the finite element and finite difference methods are described. As a first step in developing a three dimensional finite element based ablative thermal response capability, a one dimensional computer tool has been developed. The finite element method is used to discretize the governing differential equations and Galerkin's method of weighted residuals is used to derive the element equations. A code to code comparison between the current 1-D tool and the 1-D Fully Implicit Ablation and Thermal Response Program (FIAT) has been performed.

Dec John A.↗

A semi-implicit dynamic relaxation algorithm for static nonlinear structural analysis

A semiimplicit dynamic relaxation technique for solution of the nonlinear structural equlibrium equation is presented. A previously presented basic transient response analysis algorithm is employed, permitting use of one solution method and one software module for both static and dynamic analyses. A theoretical comparison of the method with explicit dynamic relaxation techniques shows that it offers a substantially improved convergence property without additional computational overhead.

Park, K. C.↗

System Identification and POD Method Applied to Unsteady Aerodynamics

The representation of unsteady aerodynamic flow fields in terms of global aerodynamic modes has proven to be a useful method for reducing the size of the aerodynamic model over those representations that use local variables at discrete grid points in the flow field. Eigenmodes and Proper Orthogonal Decomposition (POD) modes have been used for this purpose with good effect. This suggests that system identification models may also be used to represent the aerodynamic flow field. Implicit in the use of a systems identification technique is the notion that a relative small state space model can be useful in describing a dynamical system. The POD model is first used to show that indeed a reduced order model can be obtained from a much larger numerical aerodynamical model (the vortex lattice method is used for illustrative purposes) and the results from the POD and the system identification methods are then compared. For the example considered, the two methods are shown to give comparable results in terms of accuracy and reduced model size. The advantages and limitations of each approach are briefly discussed. Both appear promising and complementary in their characteristics.

Tang, Deman↗

Numerical study of the effects of icing on viscous flow over wings

An improved hybrid method for computing unsteady compressible viscous flows is presented. This method divides the computational domain into two zones. In the outer zone, the unsteady full-potential equation (FPE) is solved. In the inner zone, the Navier-Stokes equations are solved using a diagonal form of an alternating-direction implicit (ADI) approximate factorization procedure. The two zones are tightly coupled so that steady and unsteady flows may be efficiently solved. Characteristic-based viscous/inviscid interface boundary conditions are employed to avoid spurious reflections at that interface. The resulting CPU times are less than 60 percent of that required for a full-blown Navier-Stokes analysis for steady flow applications and about 60 percent of the Navier-Stokes CPU times for unsteady flows in non-vector processing machines. Applications of the method are presented for a rectangular NACA 0012 wing in low subsonic steady flow at moderate and high angles of attack, and for an F-5 wing in steady and unsteady subsonic and transonic flows. Steady surface pressures are in very good agreement with experimental data and are essentially identical to Navier-Stokes predictions. Density contours show that shocks cross the viscous/inviscid interface smoothly, so that the accuracy of full Navier-Stokes equations can be retained with a significant savings in computational time.

Sankar, L. N.↗

Distributed memory compiler methods for irregular problems: Data copy reuse and runtime partitioning

Outlined here are two methods which we believe will play an important role in any distributed memory compiler able to handle sparse and unstructured problems. We describe how to link runtime partitioners to distributed memory compilers. In our scheme, programmers can implicitly specify how data and loop iterations are to be distributed between processors. This insulates users from having to deal explicitly with potentially complex algorithms that carry out work and data partitioning. We also describe a viable mechanism for tracking and reusing copies of off-processor data. In many programs, several loops access the same off-processor memory locations. As long as it can be verified that the values assigned to off-processor memory locations remain unmodified, we show that we can effectively reuse stored off-processor data. We present experimental data from a 3-D unstructured Euler solver run on iPSC/860 to demonstrate the usefulness of our methods.

Das, Raja↗

An investigation of several numerical procedures for time-asymptotic compressible Navier-Stokes solutions

The status of an investigation of four numerical techniques for the time-dependent compressible Navier-Stokes equations is presented. Results for free shear layer calculations in the Reynolds number range from 1000 to 81000 indicate that a sequential alternating-direction implicit (ADI) finite-difference procedure requires longer computing times to reach steady state than a low-storage hopscotch finite-difference procedure. A finite-element method with cubic approximating functions was found to require excessive computer storage and computation times. A fourth method, an alternating-direction cubic spline technique which is still being tested, is also described.

Rudy, D. H.↗

Non-linear effects on sound propagation in a variable area duct carrying high subsonic Mach number flows

The general one-dimensional non-linear equation for the acoustic velocity potential in a variable area duct carrying high subsonic Mach number flows is presented and solved numerically using an implicit finite difference scheme. For the linearized equation with no flow, the present scheme is compared with the exact solution and the fourth-order Runge-Kutta method with excellent agreement for dimensionless time periods. Non-linear solutions are more sensitive to Mach numbers and exciting amplitudes and less sensitive to exciting frequencies. In general, non-linear effects can be safely neglected for low Mach number flows but must be accounted for when high Mach number flows are encountered.

Lumsdaine, E.↗

Application of viscous-inviscid interaction methods to transonic turbulent flows

Two different viscous-inviscid interaction schemes were developed for the analysis of steady, turbulent, transonic, separated flows over axisymmetric bodies. The viscous and inviscid solutions are coupled through the displacement concept using a transpiration velocity approach. In the semi-inverse interaction scheme, the viscous and inviscid equations are solved in an explicitly separate manner and the displacement thickness distribution is iteratively updated by a simple coupling algorithm. In the simultaneous interaction method, local solutions of viscous and inviscid equations are treated simultaneously, and the displacement thickness is treated as an unknown and is obtained as a part of the solution through a global iteration procedure. The inviscid flow region is described by a direct finite-difference solution of a velocity potential equation in conservative form. The potential equation is solved on a numerically generated mesh by an approximate factorization (AF2) scheme in the semi-inverse interaction method and by a successive line overrelaxation (SLOR) scheme in the simultaneous interaction method. The boundary-layer equations are used for the viscous flow region. The continuity and momentum equations are solved inversely in a coupled manner using a fully implicit finite-difference scheme.

Lee, D.↗

Numerical simulation of the hypersonic flow around lifting vehicles

A method for solving the viscous hypersonic flow field around realistic configurations is presented. The numerical procedure for generating the required finite difference grid and the two-factored implicit flow solver are described. Results are presented for the shuttle orbiter and a generic wing-body configuration at hypersonic Mach numbers.

Rizk, Yehia↗

Numerical simulation of the hypersonic flow around lifting vehicles

A method of solving the viscous hypersonic flow field around realistic configurations is presented. The numerical procedure for generating the required finite difference grid and the two-factored implicit flow solver are described. Results are presented for the shuttle orbiter and a generic wing-body configuration at hypersonic Mach numbers.

Rizk, Yehia↗

Elliptic grid generation with orthogonality and spacing control on an arbitrary number of boundaries

A procedure for the generation of two and quasi-three-dimensional grids with control of orthogonality and spacing with respect to any and/or all boundaries of the domain is described. The elliptic grid generation equations of Thompson are solved implicitly. Control of the grid behavior is achieved through the introduction of forcing functions terms in the manner of Steger and Sorenson or in a modification of the method of Hilgenstock. The forcing function terms are constructed on the boundaries and propagated into the domain using transfinite Lagrangian bivariate interpolation. An anisotropic transfinite stencil is introduced and is shown to produce excellent grid behavior particularly in the vicinity of corner singularities. Emphasis is placed on the generation of viscous grids and the method is shown to be suited for use in the generation of grids for internal as well as external flow geometries. A FORTRAN program named PISCES has been written to implement the algorithm. Examples of grids for internal and external flows are given that highlight the characteristics and behavior of the algorithm.

White, J. A.↗

Three-dimensional incompressible Navier-Stokes solver using lower-upper symmetric-Gauss-Seidel algorithm

A numerical method based on the pseudocompressibility concept is developed for solving the three-dimensional incompressible Navier-Stokes equations using the lower-upper symmetric-Gauss-Seidel implicit scheme. Very high efficiency is achieved in a new flow solver, INS3D-LU code, by accomplishing the complete vectorizability of the algorithm on oblique planes of sweep in three dimensions.

Yoon, Seokkwan↗

Simulation of tail buffet using delta wing-vertical tail configuration

Computational simulation of the vertical tail buffet problem is accomplished using a delta wing-vertical tail configuration. Flow conditions are selected such that the wing primary-vortex cores experience vortex breakdown and the resulting flow interacts with the vertical tail. This multidisciplinary problem is solved successively using three sets of equations for the fluid flow, aeroelastic deflections and grid displacements. For the fluid dynamics part, the unsteady, compressible, full Navier-Stokes equations are solved accurately in time using an implicit, upwind, flux-difference splitting, finite-volume scheme. For the aeroelastic part, the aeroelastic equation for bending vibrations is solved accurately in time using the Galerkin method and the four-stage Runge-Kutta scheme. The grid for the fluid dynamics computations is updated every few time steps using a third set of interpolation equations. The computational application includes a delta wing of aspect ratio 1 and a rectangular vertical tail of aspect ratio 2, which is placed at 0.5 root-chord length downstream of the wing trailing edge. The wing angle of attack is 35 deg and the flow Mach number and Reynolds number are 0.4 and 10,000, respectively.

Kandil, Osama A.↗