Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “direct solver”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Computational strategies for three-dimensional flow simulations on distributed computer systems

This research effort is directed towards an examination of issues involved in porting large computational fluid dynamics codes in use within the industry to a distributed computing environment. This effort addresses strategies for implementing the distributed computing in a device independent fashion and load balancing. A flow solver called TEAM presently in use at Lockheed Aeronautical Systems Company was acquired to start this effort. The following tasks were completed: (1) The TEAM code was ported to a number of distributed computing platforms including a cluster of HP workstations located in the School of Aerospace Engineering at Georgia Tech; a cluster of DEC Alpha Workstations in the Graphics visualization lab located at Georgia Tech; a cluster of SGI workstations located at NASA Ames Research Center; and an IBM SP-2 system located at NASA ARC. (2) A number of communication strategies were implemented. Specifically, the manager-worker strategy and the worker-worker strategy were tested. (3) A variety of load balancing strategies were investigated. Specifically, the static load balancing, task queue balancing and the Crutchfield algorithm were coded and evaluated. (4) The classical explicit Runge-Kutta scheme in the TEAM solver was replaced with an LU implicit scheme. And (5) the implicit TEAM-PVM solver was extensively validated through studies of unsteady transonic flow over an F-5 wing, undergoing combined bending and torsional motion. These investigations are documented in extensive detail in the dissertation, 'Computational Strategies for Three-Dimensional Flow Simulations on Distributed Computing Systems', enclosed as an appendix.

Sankar, Lakshmi N.↗
NASA NTRS2001

Multigrid approaches to non-linear diffusion problems on unstructured meshes

The efficiency of three multigrid methods for solving highly non-linear diffusion problems on two-dimensional unstructured meshes is examined. The three multigrid methods differ mainly in the manner in which the nonlinearities of the governing equations are handled. These comprise a non-linear full approximation storage (FAS) multigrid method which is used to solve the non-linear equations directly, a linear multigrid method which is used to solve the linear system arising from a Newton linearization of the non-linear system, and a hybrid scheme which is based on a non-linear FAS multigrid scheme, but employs a linear solver on each level as a smoother. Results indicate that all methods are equally effective at converging the non-linear residual in a given number of grid sweeps, but that the linear solver is more efficient in cpu time due to the lower cost of linear versus non-linear grid sweeps.

Mavriplis, Dimitri J.↗
NASA NTRS2004

Test and Analyses of a Composite Multi-Bay Fuselage Panel Under Uni-Axial Compression

A composite panel containing three stringers and two frames cut from a vacuum-assisted resin transfer molded (VaRTM) stitched fuselage article was tested under uni-axial compression loading. The stringers and frames divided the panel into six bays with two columns of three bays each along the compressive loading direction. The two frames were supported at the ends with pins to restrict the out-of-plane translation. The free edges of the panel were constrained by knife-edges. The panel was modeled with shell finite elements and analyzed with ABAQUS nonlinear solver. The nonlinear predictions were compared with the test results in out-of-plane displacements, back-to-back surface strains on stringer flanges and back-to-back surface strains at the centers of the skin-bays. The analysis predictions were in good agreement with the test data up to post-buckling.

Li, Jian↗
NASA NTRS2004

Adjoint Algorithm for CAD-Based Shape Optimization Using a Cartesian Method

Adjoint solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape optimization. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (geometric parameters that control the shape). More recently, emerging adjoint applications focus on the analysis problem, where the adjoint solution is used to drive mesh adaptation, as well as to provide estimates of functional error bounds and corrections. The attractive feature of this approach is that the mesh-adaptation procedure targets a specific functional, thereby localizing the mesh refinement and reducing computational cost. Our focus is on the development of adjoint-based optimization techniques for a Cartesian method with embedded boundaries.12 In contrast t o implementations on structured and unstructured grids, Cartesian methods decouple the surface discretization from the volume mesh. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin et developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation. Central to automated shape optimization algorithms is the issue of geometry modeling and control. The need to optimize complex, "real-life" geometry provides a strong incentive for the use of parametric-CAD systems within the optimization procedure. In previous work, we presented an effective optimization framework that incorporates a direct-CAD interface. In this work, we enhance the capabilities of this framework with efficient gradient computations using the discrete adjoint method. We present details of the adjoint numerical implementation, which reuses the domain decomposition, multigrid, and time-marching schemes of the flow solver. Furthermore, we explain and demonstrate the use of CAD in conjunction with the Cartesian adjoint approach. The final paper will contain a number of complex geometry, industrially relevant examples with many design variables to demonstrate the effectiveness of the adjoint method on Cartesian meshes.

Nemec, Marian↗
NASA NTRS2019

Measured and Simulated Acoustic Signature of a Full-Scale Aircraft with Airframe Noise Reduction Technology Installed

Microphone phased-array and pole-mounted microphone data gathered during the NASA Acoustics Research Measurements flight tests were used to benchmark results from companion full-scale aeroacoustics simulations. Conducted with the lattice Boltzmann solver PowerFLOW®, the simulations predicted the acoustic behavior of various tested aircraft configurations. Emphasis was placed on those flown during the third flight test - a Fowler flap-equipped Gulfstream G-III with and without noise abatement technology on the main landing gear. Direct comparisons between experimental and synthetic microphone phasedarray data were achieved by applying the same processing and deconvolution technique to both sets of data. To extend the validation of the computations to the metric used for noise certification, the Effective Perceived Noise Level, a high-fidelity digital model of the nose landing gear, which was excluded from earlier computations, was developed and integrated into the G-III aircraft geometry. The acoustic study presented here demonstrates that the simulated beamform maps and corresponding integrated farfield spectra accurately predict the locations and strengths of the prominent airframe noise sources present on the G-III aircraft.

Khorrami, Mehdi R.↗
NASA NTRSDate not supplied

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗
NASA NTRS1995

Large-eddy simulation of flow past a circular cylinder

Some of the most challenging applications of large-eddy simulation are those in complex geometries where spectral methods are of limited use. For such applications more conventional methods such as finite difference or finite element have to be used. However, it has become clear in recent years that dissipative numerical schemes which are routinely used in viscous flow simulations are not good candidates for use in LES of turbulent flows. Except in cases where the flow is extremely well resolved, it has been found that upwind schemes tend to damp out a significant portion of the small scales that can be resolved on the grid. Furthermore, it has been found that even specially designed higher-order upwind schemes that have been used successfully in the direct numerical simulation of turbulent flows produce too much dissipation when used in conjunction with large-eddy simulation. The objective of the current study is to perform a LES of incompressible flow past a circular cylinder at a Reynolds number of 3900 using a solver which employs an energy-conservative second-order central difference scheme for spatial discretization and compare the results obtained with those of Beaudan & Moin (1994) and with the experiments in order to assess the performance of the central scheme for this relatively complex geometry.

Mittal, R.↗
NASA NTRS1992

Flutter analysis of a supersonic cascade in time domain using an ADI Euler solver

The aeroelastic stability of a two-dimensional cascade oscillating in supersonic axial flow is analyzed in the time domain. The aeroelastic model consists of a single degree of freedom typical section structural model for each blade of the cascade and an unsteady two-dimensional cascade aerodynamic model based on the Euler equations. The Euler equations are solved using a time accurate Alternating Direction Implicit (ADI) solution scheme. The aeroelastic equations are integrated in time. The effect of interblade phase angle is included in the aeroelastic analysis by an appropriate choice of initial and boundary conditions. Flutter predictions are obtained from the time response of a flat plate cascade in single degree of freedom pitching motion. The results correlate well with those obtained from a separate frequency domain flutter analysis for all values of interblade phase angles considered. Flutter results are then presented for cascades having airfoil sections representative of a supersonic throughflow fan. The validity of the time integration method for a cascade of airfoils at various interblade phase angles is demonstrated.

Reddy, T. S. R.↗
NASA NTRS2013

Computational Studies of Magnetic Nozzle Performance

An extensive literature review of magnetic nozzle research has been performed, examining previous work, as well as a review of fundamental principles. This has allow us to catalog all basic physical mechanisms which we believe underlie the thrust generation process. Energy conversion mechanisms include the approximate conservation of the magnetic moment adiabatic invariant, generalized hall and thermoelectric acceleration, swirl acceleration, thermal energy transformation into directed kinetic energy, and Joule heating. Momentum transfer results from the interaction of the applied magnetic field with currents induced in the plasma plume., while plasma detachment mechanisms include resistive diffusion, recombination and charge exchange collisions, magnetic reconnection, loss of adiabaticity, inertial forces, current closure, and self-field detachment. We have performed a preliminary study of Hall effects on magnetic nozzle jets with weak guiding magnetic fields and weak expansions (p(sub jet) approx. = P(sub background)). The conclusion from this study is that the Hall effect creates an azimuthal rotation of the plasma jet and, more generally, creates helical structures in the induced current, velocity field, and magnetic fields. We have studied plasma jet expansion to near vacuum without a guiding magnetic field, and are presently including a guiding magnetic field using a resistive MHD solver. This research is progressing toward the implementation of a full generalized Ohm's law solver. In our paper, we will summarize the basic principle, as well as the literature survey and briefly review our previous results. Our most recent results at the time of submittal will also be included. Efforts are currently underway to construct an experiment at the University of Michigan Plasmadynamics and Electric Propulsion Laboratory (PEPL) to study magnetic nozzle physics for a RF-thruster. Our computational study will work directly with this experiment to validate the numerical model, in order to study magnetic nozzle physics and optimize magnetic nozzle design. Preliminary results from the PEPL experiment will also be presented.

Ebersohn, Frans H.↗
NASA NTRS2015

Tetrahedral-Mesh Simulation of Turbulent Flows with the Space-Time Conservative Schemes

Direct numerical simulations of turbulent flows are predominantly carried out using structured, hexahedral meshes despite decades of development in unstructured mesh methods. Tetrahedral meshes offer ease of mesh generation around complex geometries and the potential of an orientation free grid that would provide un-biased small-scale dissipation and more accurate intermediate scale solutions. However, due to the lack of consistent multi-dimensional numerical formulations in conventional schemes for triangular and tetrahedral meshes at the cell interfaces, numerical issues exist when flow discontinuities or stagnation regions are present. The space-time conservative conservation element solution element (CESE) method - due to its Riemann-solver-free shock capturing capabilities, non-dissipative baseline schemes, and flux conservation in time as well as space - has the potential to more accurately simulate turbulent flows using unstructured tetrahedral meshes. To pave the way towards accurate simulation of shock/turbulent boundary-layer interaction, a series of wave and shock interaction benchmark problems that increase in complexity, are computed in this paper with triangular/tetrahedral meshes. Preliminary computations for the normal shock/turbulence interactions are carried out with a relatively coarse mesh, by direct numerical simulations standards, in order to assess other effects such as boundary conditions and the necessity of a buffer domain. The results indicate that qualitative agreement with previous studies can be obtained for flows where, strong shocks co-exist along with unsteady waves that display a broad range of scales, with a relatively compact computational domain and less stringent requirements for grid clustering near the shock. With the space-time conservation properties, stable solutions without any spurious wave reflections can be obtained without a need for buffer domains near the outflow/farfield boundaries. Computational results for the isotropic turbulent flow decay, at a relatively high turbulent Mach number, show a nicely behaved spectral decay rate for medium to high wave numbers. The high-order CESE schemes offer very robust solutions even with the presence of strong shocks or widespread shocklets. The explicit formulation in conjunction with a close to unity theoretical upper Courant number bound has the potential to offer an efficient numerical framework for general compressible turbulent flow simulations with unstructured meshes.

Chang, Chau-Lyan↗
NASA NTRS2021

A Study into Validating a Coupled Method of Characteristics and Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

Direct Simulation Monte Carlo Analysis Code↗
NASA NTRS1994

Canonical-variables multigrid method for steady-state Euler equations

A novel approach for the solution of inviscid flow problems for subsonic compressible flows is described. The approach is based on canonical forms of the equations, in which subsystems governed by hyperbolic operators are separated from those governed by elliptic ones. The discretizations used as well as the iterative techniques for the different subsystems are inherently different. Hyperbolic parts, which describe, in general, propagation phenomena, are discretized using upwind schemes and are solved by marching techniques. Elliptic parts, which are directionally unbiased, are discretized using h-elliptic central discretizations, and are solved by pointwise relaxations together with coarse grid acceleration. The resulting discretization schemes introduce artificial viscosity only for the hyperbolic parts of the system; thus a smaller total artificial viscosity is used, while the multigrid solvers used are much more efficient. Solutions of the subsonic compressible Euler equations are achieved at the same efficiency as the full potential equation.

Taasan, Shlomo↗
NASA NTRSDate not supplied

Computational Analysis of X-57 Rate Derivatives and Comparison to Wind Tunnel Data

Rate derivates of the wind tunnel model of X-57 were evaluated using hysteresis from computational fluid dynamics simulation and compared with that computed using the wind tunnel experimental data. To accomplish the task, hysteresis of the pitch, roll, and yaw oscillations were produced numerically using two different CFD solvers: STARCCM+ and Launch Ascent and Vehicle Analysis. The rate derivatives were computed using same parameter identification technique on CFD-based and wind tunnel-based hysteresis data. Results show CFD-based pitch rate derivatives have large error, up to 73%, while the lateral and direction derivates compared well with errors within 5% with respect to the that computed using the wind tunnel data. This may be due to lack of wind tunnel wall and sting model in the CFD simulation as well as lack of wall correction in the wind tunnel data which has larger impact on pitch than lateral and directional oscillations. Geometrical features such as the shape, slope of extremities, and enclosed area were compared between CFD-based and wind tunnel-based hysteresis. Results show that area enclosed by the hysteresis had largest influence on the accuracy of the computed rate derivatives given that shape of the hysteresis were roughly equivalent.

Seung Y. Yoo↗
NASA NTRS1986

On the structure of parallelism in a highly concurrent PDE solver

A parallel multigrid algorithm for solving elliptic partial differential equations is developed and evaluated. A V-cycle multigrid method is altered to increase the degree of parallelism. A numerical analysis of the resulting concurrent-iteration multigrid algorithm is performed; its architectural implications are considered; highly parallel systems without shared memory are examined (including mesh-connected arrays, mesh-shuffle-connected systems, permutation networks, and direct VLSI embeddings); and the results of numerical experiments are presented in tables and graphs.

Gannon, D.↗
NASA NTRS1987

A study of multi-body aerodynamic interference at transonic Mach numbers

A wind tunnel experiment involving single, double, and triple combinations of mutually interfering generic, unfinned aircraft stores has been conducted. Each combination of stores was tested at Mach numbers from 0.60 to 1.20 and at angles of attack from 0 to 25 deg for the single store and from 0 to 6 deg for the double and triple store configurations. Extensive axial and circumferential pressure and flow visualization data at each store location were obtained. Euler solutions for each configuration at 0 deg incidence have been generated and compared with experimental data. This comparison indicates an Euler flow solver can yield accurate predictions of the location and magnitude of multibody interference provided an appropriate grid is used and the viscous effects associated with these configurations remain small. The data indicate multibody interference in the transonic region increases as the freestream Mach number approaches 1 from either direction, and subsides as the Mach number moves away from sonic conditions. This interference is characterized by a large, localized reduction in pressure on the inboard surfaces of the bodies which results in forces that draw the configuration closer together.

Cottrell, Charles J.↗
NASA NTRS1996

OVERFLOW-Interaction with Industry

A Navier-Stokes flow solver, OVERFLOW, has been developed by researchers at NASA Ames Research Center to use overset (Chimera) grids to simulate the flow about complex aerodynamic shapes. Primary customers of the OVERFLOW flow solver and related software include McDonnell Douglas and Boeing, as well as the NASA Focused Programs for Advanced Subsonic Technology (AST) and High Speed Research (HSR). Code development has focused on customer issues, including improving code performance, ability to run on workstation clusters and the NAS SP2, and direct interaction with industry on accuracy assessment and validation. Significant interaction with NAS has produced a capability tailored to the Ames computing environment, and code contributions have come from a wide range of sources, both within and outside Ames.

Buning, Pieter G.↗
NASA NTRS2003

Performance Analysis of a Hybrid Overset Multi-Block Application on Multiple Architectures

This paper presents a detailed performance analysis of a multi-block overset grid compu- tational fluid dynamics app!ication on multiple state-of-the-art computer architectures. The application is implemented using a hybrid MPI+OpenMP programming paradigm that exploits both coarse and fine-grain parallelism; the former via MPI message passing and the latter via OpenMP directives. The hybrid model also extends the applicability of multi-block programs to large clusters of SNIP nodes by overcoming the restriction that the number of processors be less than the number of grid blocks. A key kernel of the application, namely the LU-SGS linear solver, had to be modified to enhance the performance of the hybrid approach on the target machines. Investigations were conducted on cacheless Cray SX6 vector processors, cache-based IBM Power3 and Power4 architectures, and single system image SGI Origin3000 platforms. Overall results for complex vortex dynamics simulations demonstrate that the SX6 achieves the highest performance and outperforms the RISC-based architectures; however, the best scaling performance was achieved on the Power3.

Djomehri, M. Jahed↗
NASA NTRS2008

Enhancements on the Convex Programming Based Powered Descent Guidance Algorithm for Mars Landing

In this paper, we present enhancements on the powered descent guidance algorithm developed for Mars pinpoint landing. The guidance algorithm solves the powered descent minimum fuel trajectory optimization problem via a direct numerical method. Our main contribution is to formulate the trajectory optimization problem, which has nonconvex control constraints, as a finite dimensional convex optimization problem, specifically as a finite dimensional second order cone programming (SOCP) problem. SOCP is a subclass of convex programming, and there are efficient SOCP solvers with deterministic convergence properties. Hence, the resulting guidance algorithm can potentially be implemented onboard a spacecraft for real-time applications. Particularly, this paper discusses the algorithmic improvements obtained by: (i) Using an efficient approach to choose the optimal time-of-flight; (ii) Using a computationally inexpensive way to detect the feasibility/ infeasibility of the problem due to the thrust-to-weight constraint; (iii) Incorporating the rotation rate of the planet into the problem formulation; (iv) Developing additional constraints on the position and velocity to guarantee no-subsurface flight between the time samples of the temporal discretization; (v) Developing a fuel-limited targeting algorithm; (vi) Initial result on developing an onboard table lookup method to obtain almost fuel optimal solutions in real-time.

Guidance↗