SEARCH · Engineering Papers

Results for “solvers”

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,495 records · Page 83

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Numerical study of the effects of icing on fixed and rotary wing performance

The sectional and total aerodynamic load characteristics and performance degradation of swept wings and helicopter rotors have been studied using a three-dimensional, compressible Navier-Stokes solver. Correlations of predictions with experimental data for swept wings with and without leading-edge ice formation show the ability of the present computational technique to accurately predict both the distributed surface pressures and integrated sectional loads. The leading-edge flow separation and reattachment on the wing surface associated with the leading-edge ice are also captured well showing a vortex formation and the spanwise migration of the flow inside the separated flow region. In the case of the helicopter rotors in hover, the rotor thrust loss and the torque penalties due to the leading-edge ice formation are numerically demonstrated.

Kwon, Oh J.↗

NASA NTRS1991

A two-equation turbulence model for compressible reacting flows

A two-equation turbulence model (k-epsilon) has been modified to model compressible reacting flows. A compressibility correction model based on modeling the dilatational terms in the Reynolds stress equations has been used. Computations are carried out using the Navier-Stokes solver SPARK which includes a finite rate chemistry model for hydrogen-air combustion. Initially, planar compressible mixing layers are simulated. Comparisons of the predictions with available experimental data and the predictions of a compressible Reynolds stress closure indicate that the model is well suited for the study of such flows. The decrease in the growth rate of the mixing layer with increasing convective Mach number is well predicted. The model is then tested on a three-dimensional supersonic flow problem where solid boundaries are present. An axisymmetric nonreacting mixing layer flow involving dissimilar gases is evaluated next using the model and the predictions are compared with available experimental data. Finally, two reacting turbulent flow configurations have been studied and comparison of predictions and experimental data has been done. The comparisons indicate that the model predicts the overall features of the test cases very well.

Narayan, J. R.↗

NASA NTRS1991

Inviscid steady/unsteady flow calculations

The solution of the Euler equations using a flux splitting procedure is considered for low subsonic to high supersonic flows. Steady and unsteady, internal and external flow fields, are computed. For transient flows, a direct sparse matrix solver is applied to compute the flow field at each instant of time. Oscillation free normal and oblique shocks are captured. Unstart and restart of a simplified two-dimensional inlet is investigated.

Pordal, H. S.↗

NASA NTRS1990

Development of unstructured grid methods for steady and unsteady aerodynamic analysis

The current status of the development of unstructured grid methods in the Unsteady Aerodynamic Branch at NASA-Langley is described. These methods are being developed for steady and unsteady aerodynamic applications. The flow solvers that were developed for the solution of the unsteady Euler and Navier-Stokes equations are highlighted and selected results are given which demonstrate various features of the capability. The results demonstrate 2-D and 3-D applications for both steady and unsteady flows. Comparisons are also made with solutions obtained using a structured grid code and with experimental data to determine the accuracy of the unstructured grid methodology. These comparisons show good agreement which thus verifies the accuracy.

Batina, John T.↗

NASA NTRS1991

Standard test case for a low-speed, turbulent junction vortex flow

The mean flow structure upstream, around within, and in the near wake of a turbulent junction or horseshoe vortex is reported for an incompressible subsonic flow. Measurements of the primitive variables of velocity and pressure are reported on all surfaces bounding the region of the vortex flow and on three transverse and one streamwise plane within the flowfield itself for comparisons between the measured and any calculated flow variables. Detailed surface flow visualizations and some direct force measurements of surface shear stress are also available. The data is a highly detailed, coherent, self-consistent set offered to the computational fluid mechanics community as a standard test case for the evaluation of the capability of numerical solvers intended for predicting the flowfield in such a complex, separated, three-dimensional turbulent flow. This data base is available for copying to user supplied tapes or for transmission via BITNET, as well as in two National Technical Information Service (NTIS) reports.

Pierce, F. J.↗

NASA NTRS1991

The aerodynamic characteristics of vortex ingestion for the F/A-18 inlet duct

A Reduced Navier-Stokes (RNS) solution technique was successfully combined with the concept of partitioned geometry and mesh generation to form a very efficient 3-D RNS code aimed at the analysis-design engineering environment. Partitioned geometry and mesh generation is a pre-processor to augment existing geometry and grid generation programs which allows the solver to (1) recluster an existing gridlife mesh lattice, and (2) perturb an existing gridfile definition to alter the cross-sectional shape and inlet duct centerline distribution without returning to the external geometry and grid generator. The present results provide a quantitative validation of the initial value space marching 3-D RNS procedure and demonstrates accurate predictions of the engine face flow field, with a separation present in the inlet duct as well as when vortex generators are installed to suppress flow separation. The present results also demonstrate the ability of the 3-D RNS procedure to analyze the flow physics associated with vortex ingestion in general geometry ducts such as the F/A-18 inlet. At the conditions investigated, these interactions are basically inviscid like, i.e., the dominant aerodynamic characteristics have their origin in inviscid flow theory.

Anderson, Bernhard H.↗

NASA NTRS1991

Shedding patterns of the near-wake vortices behind a circular cylinder

The unsteady incompressible Navier-Stokes equations have been accurately solved for the laminar flow past a circular cylinder in the Reynolds number range 50-200. A direct elliptic solver called the SEVP is used to rapidly advance the streamfunction in time, facilitating the overall convergence to the fully periodic or quasi-steady state. A new integral-series method is developed for the far-field streamfunction condition on a finite two-dimensional computational domain. The use of fourth-order Hermitain relations for the convection terms in the conservation-form vorticity transport equation has also contributed to the good comparison of the present results with the earlier experimental data. The vortex-shedding patterns visualized by the experimentalist are numerically reproduced here in the given Reynolds number range. Discussions that may be helpful in interpreting the behavior of the shedding frequency are presented in the main text.

Sa, Jong-Youb↗

NASA NTRS1990

Trajectory optimization for real-time guidance. I - Time-varying LQR on a parallel processor

A key algorithmic element of a real-time trajectory optimization hardware/software implementation, the quadratic program (QP) solver element, is presented. The purpose of the effort is to make nonlinear trajectory optimization fast enough to provide real-time commands during guidance of a vehicle such as an aeromaneuvering orbiter. Many methods of nonlinear programming require the solution of a QP at each iteration. In the trajectory optimization case the QP has a special dynamic programming structure, a LQR-like structure. QP algorithm speed is increased by taking advantage of this special structure and by parallel implementation.

Psiaki, Mark L.↗

NASA NTRS1991

Parallel computation in a three-dimensional elastic-plastic finite-element analysis

CRAY 2 and CRAY Y-MP evaluation runs are undertaken for a three-dimensional, elastoplastic FEM code implementation of the 'autotasking' parallel-processing technique, which is used in all code components except the matrix-equation solver. For a typical example problem in the case of a four-processor CRAY 2, a speedup factor of 2.1 was found achievable in a dedicated environment, and 1.7 in a multiuser environment; for a three-processor CRAY Y-MP, speedup in a multiuser environment was about 2.4. Only minimal effort was required to implement autotasking in this implementation code.

Shivakumar, K. N.↗

NASA NTRS1991

Transonic adaptive flutter suppression using approximate unsteady time domain aerodynamics

A digital adaptive controller is applied to the active flutter suppression problem of a wing under time varying flight conditions in subsonic and transonic flow. Linear quadratic controller gain at each time step is obtained using an iterative Riccati solver. The digital adaptive optimal controller is robust with respect to the unknown external loads. Flutter and divergence instabilities are simultaneously suppressed using a trailing-edge control surface and displacement sensing. A new transonic unsteady aerodynamic approximation methodology is developed which enables one to carry out the rapid calculation required for transonic aeroservoelastic applications. This approximation is based on a combination of unsteady subsonic aerodynamics combined with a transonic correction procedure. Aeroservoelastic transient time response is obtained using Roger's approximation, state transition matrices and an iterative time marching algorithm. The aeroservoelastic system in the time domain is modeled using a deterministic ARMA model together with a parameter estimator. Transonic flutter boundaries of a wing structure are computed, in the time domain, using an estimated aeroelastic system matrix and are in good agreement with experimental data for the low transonic Mach number range.

Pak, Chan-Gi↗

NASA NTRS1990

SIAM Conference on Parallel Processing for Scientific Computing, 4th, Chicago, IL, Dec. 11-13, 1989, Proceedings

Attention is given to such topics as an evaluation of block algorithm variants in LAPACK and presents a large-grain parallel sparse system solver, a multiprocessor method for the solution of the generalized Eigenvalue problem on an interval, and a parallel QR algorithm for iterative subspace methods on the CM2. A discussion of numerical methods includes the topics of asynchronous numerical solutions of PDEs on parallel computers, parallel homotopy curve tracking on a hypercube, and solving Navier-Stokes equations on the Cedar Multi-Cluster system. A section on differential equations includes a discussion of a six-color procedure for the parallel solution of elliptic systems using the finite quadtree structure, data parallel algorithms for the finite element method, and domain decomposition methods in aerodynamics. Topics dealing with massively parallel computing include hypercube vs. 2-dimensional meshes and massively parallel computation of conservation laws. Performance and tools are also discussed.

Dongarra, Jack↗

NASA NTRS1990

Domain decomposition methods in aerodynamics

Compressible Euler equations are solved for two-dimensional problems by a preconditioned conjugate gradient-like technique. An approximate Riemann solver is used to compute the numerical fluxes to second order accuracy in space. Two ways to achieve parallelism are tested, one which makes use of parallelism inherent in triangular solves and the other which employs domain decomposition techniques. The vectorization/parallelism in triangular solves is realized by the use of a recording technique called wavefront ordering. This process involves the interpretation of the triangular matrix as a directed graph and the analysis of the data dependencies. It is noted that the factorization can also be done in parallel with the wave front ordering. The performances of two ways of partitioning the domain, strips and slabs, are compared. Results on Cray YMP are reported for an inviscid transonic test case. The performances of linear algebra kernels are also reported.

Venkatakrishnan, V.↗

NASA NTRS1990

A flux-split solution procedure for unsteady flow calculations

The solution of reduced Navier Stokes (RNS) equations is considered using a flux-split procedure. Unsteady flow in a two dimensional engine inlet is computed. The problems of unstart and restart are investigated. A sparse matrix direct solver combined with domain decomposition strategy is used to compute the unsteady flow field at each instant of time. Strong shock-boundary layer interaction, time varying shocks and time varying recirculation regions are efficiently captured.

Pordal, H. S.↗

NASA NTRS1991

Computational aspects of chemically reacting flows

An implicit, time-dependent approach is taken to solve the inviscid conservative equations in a vertex-based finite-volume formulation. The convective fluxes are estimated by a modified Riemann solver, and higher-order accuracy is obtained by means of the MUSCL. Implicit techniques are adopted to deal with the equation stiffness as a result of the clustered grid and chemical source term. The real-gas Roe matrix and the ADI and LU schemes have been tested on the sample problems: hypersonic flow over blunt configurations, such as the forebody of a double ellipse, and the near wake of a conic aerobrake. Temperature distribution and standoff distance are validated against the shock-fitting, central-differencing results. Equilibrium temperature may be recovered by artificially increasing the chemical rates at the wall. The present method satisfactorily resolves the flow structure, including the strong expansion behind the flat base. The relative merits and performance of the real-gas TVD and implicit schemes are discussed briefly.

Li, C. P.↗

NASA NTRS1991

Multitasking the INS3D-LU code on the Cray Y-MP

This paper presents the results of multitasking the INS3D-LU code on eight processors. The code is a full Navier-Stokes solver for incompressible fluid in three dimensional generalized coordinates using a lower-upper symmetric-Gauss-Seidel implicit scheme. This code has been fully vectorized on oblique planes of sweep and parallelized using autotasking with some directives and minor modifications. The timing results for five grid sizes are presented and analyzed. The code has achieved a processing rate of over one Gflops.

Fatoohi, Rod↗

NASA NTRS1991

Grid convergence for adaptive methods

The ability of adaptive methods to obtain accurate results is examined using two different Euler solvers for a near-sonic flow containing several important flow features. It is shown that the accuracy obtained can be greatly affected by the lack of resolution of smooth portions of the flow caused by adapting only to the more prevalent flow features such as discontinuities. In particular, common methods of adaptation can lead to results in which shocks are well resolved but whose locations are highly inaccurate due to the lack of resolution of the smoother regions. An explanation for this behavior is given and a correction is proposed.

Warren, Gary P.↗

NASA NTRS1991

Conservative multizonal interface algorithm for the 3-D Navier-Stokes equations

A conservative zonal interface algorithm using features of both structured and unstructured mesh CFD technology is presented. The flow solver within each of the zones is based on structured mesh CFD technology. The interface algorithm was implemented into two three-dimensional Navier-Stokes finite volume codes and was found to yield good results.

Klopfer, G. H.↗

NASA NTRS1991

A numerical investigation of turbulent baseflow

The turbulent base flow of a cone with a sting in hypersonic flow is investigated using a time-dependent Navier-Stokes code called UWIN. A PNS solver provides the upstream flow conditions which drive the development of the base flow. The computational method is described along with its application to the base flow problem. The grid generation and boundary conditions are described. The results and the summary conclude this note.

Chaussee, Denny S.↗