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.

468 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

CFD Based Design of a Filming Injector for N+3 Combustors

An effort was undertaken to perform CFD analysis of fluid flow in Lean-Direct Injection (LDI) combustors with axial swirl-venturi elements for next-generation LDI-3 combustor design. The National Combustion Code (NCC) was used to perform non-reacting and two-phase reacting flow computations for a newly-designed pre-filming type fuel injector LDI-3 injector, in a single-injector and a five-injector array configuration. All computations were performed with a consistent approach of mesh-optimization, spray-modeling, ignition and kinetics-modeling. Computational predictions of the aerodynamics of the single-injector were used to arrive at an optimized main-injector design that meets effective area and fuel-air mixing criteria. Emissions (EINOx) characteristics were predicted for a medium-power engine cycle condition, and will be compared with data when it is made available from experimental measurements. The use of a PDF-like turbulence-chemistry interaction model with NCC's Time-Filtered Navier-Stokes (TFNS) solver is shown to produce a significant impact on the CFD results, when compared with a laminar-chemistry TFNS approach for the five-injector computations.

combustion chambers↗
NASA NTRSDate not supplied

Aeroacoustic Computations of a Transonic Truss-Braced Wing Aircraft: Part 1 – Aero-dynamic and Airframe Noise Simulations

Computational results are presented for an 8%-scale, full-span, transonic truss-braced wing (TTBW) model simulated as installed in the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel and in free-air conditions. The simulations were conducted with the lat-tice Boltzmann solver PowerFLOW® to capture the time-accurate characteristics of the flow. The aerodynamic behavior of the aircraft was investigated in the landing configuration, with high-lift devices and landing gear deployed, as well as in the clean (cruise) configuration, with these components stowed. Analyses were performed on local flow quantities, global forces, and time-averaged surface pressures. Aerodynamic quantities were shown to be sensitive to mesh resolution levels, driven by small geometric features inherent to the TTBW model. Flow features of the TTBW model were examined, with the wing/strut configuration of this model presenting unique behaviors generally not found in conventional transport aircraft. Near-field, time-dependent flow quantities obtained from the scale-resolving simulations were used in conjunction with a Ffowcs-Williams and Hawkings integral approach to predict the far-field airframe noise signature of this advanced concept. The effects of permeable data sur-face end caps on the far-field noise spectrum in the flyover direction were determined to be negligible.

Transonic Truss-Braced Wing↗
NASA NTRS1990

Multiprocessing on supercomputers for computational aerodynamics

Very little use is made of multiple processors available on current supercomputers (computers with a theoretical peak performance capability equal to 100 MFLOPs or more) in computational aerodynamics to significantly improve turnaround time. The productivity of a computer user is directly related to this turnaround time. In a time-sharing environment, the improvement in this speed is achieved when multiple processors are used efficiently to execute an algorithm. The concept of multiple instructions and multiple data (MIMD) through multi-tasking is applied via a strategy which requires relatively minor modifications to an existing code for a single processor. Essentially, this approach maps the available memory to multiple processors, exploiting the C-FORTRAN-Unix interface. The existing single processor code is mapped without the need for developing a new algorithm. The procedure for building a code utilizing this approach is automated with the Unix stream editor. As a demonstration of this approach, a Multiple Processor Multiple Grid (MPMG) code is developed. It is capable of using nine processors, and can be easily extended to a larger number of processors. This code solves the three-dimensional, Reynolds averaged, thin-layer and slender-layer Navier-Stokes equations with an implicit, approximately factored and diagonalized method. The solver is applied to generic oblique-wing aircraft problem on a four processor Cray-2 computer. A tricubic interpolation scheme is developed to increase the accuracy of coupling of overlapped grids. For the oblique-wing aircraft problem, a speedup of two in elapsed (turnaround) time is observed in a saturated time-sharing environment.

Yarrow, Maurice↗
NASA NTRS2017

Jet Noise Prediction Using Hybrid RANS/LES with Structured Overset Grids

Hybrid RANS/LES simulations using the structured overset grid approach and a low dissipation finite-difference method within the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are presented for jet noise prediction. The simulations are part of a validation effort to demonstrate jet noise prediction capability, to assess noise characteristics of next generation quiet supersonic commercial jets. Results are compared with experimental data acquired in the Small Hot Jet Acoustic Rig in the Aeroacoustic Propulsion Laboratory at NASA Glenn Research Center. Details of the structured overset grids, numerical discretization, and turbulence model are provided. Near-field comparisons to PIV data and far-field comparisons to microphone data are discussed. Excellent agreement of time-averaged mean quantities along the jet center-line and lip-line are obtained. Good agreement with root-mean-squared (RMS) quantities along the jet center-line are also obtained. An over prediction of lip-line RMS at the nozzle exit is observed leading to a small over-prediction of side-line far-field noise (less than 2 dB). Good agreement with far-field noise in the Mach wave radiation directions is achieved.

Housman, Jeffrey A.↗
NASA NTRS1993

Development of computational methods for heavy lift launch vehicles

The research effort has been focused on the development of an advanced flow solver for complex viscous turbulent flows with shock waves. The three-dimensional Euler and full/thin-layer Reynolds-averaged Navier-Stokes equations for compressible flows are solved on structured hexahedral grids. The Baldwin-Lomax algebraic turbulence model is used for closure. The space discretization is based on a cell-centered finite-volume method augmented by a variety of numerical dissipation models with optional total variation diminishing limiters. The governing equations are integrated in time by an implicit method based on lower-upper factorization and symmetric Gauss-Seidel relaxation. The algorithm is vectorized on diagonal planes of sweep using two-dimensional indices in three dimensions. A new computer program named CENS3D has been developed for viscous turbulent flows with discontinuities. Details of the code are described in Appendix A and Appendix B. With the developments of the numerical algorithm and dissipation model, the simulation of three-dimensional viscous compressible flows has become more efficient and accurate. The results of the research are expected to yield a direct impact on the design process of future liquid fueled launch systems.

Yoon, Seokkwan↗
NASA NTRS2007

Numerical Investigation of Plasma Flows in Magnetic Nozzles

Magnetic nozzles are used in many laboratory experiments in which plasma flows are to be confined, cooled. accelerated, or directed. At present, however, there is no generally accepted theoretical description that explains the phenomena of plasma detachment from an externally-imposed magnetic field. This is an important problem in the field of plasma propulsion, where the ionized gas must detach from the applied, solenoidal magnetic field to realize thrust production. In this paper we simulate a plasma flowing in the presence of an applied magnetic field using a multidimensional numerical simulation tool that includes theoretical models of the various dispersive and dissipative processes present in the plasma. This is an extension of the simulation tool employed in previous work by Sankaran et al. The new tool employs the same formulation of the governing equation set. but retains the axial and radial components of magnetic field and the azimuthal component of velocity that were neglected in other works. We aim to compare the computational results with the various proposed magnetic nozzle detachment theories to develop an understanding of the physical mechanisms that cause detachment. An applied magnetic field topology is obtained using a magnetostatic field solver and this field is superimposed on the time-dependent magnetic field induced in the plasma to provide a self-consistent field description. The applied magnetic field and model geometry match those found in experiments by Kurtki and Okada. We model this geometry because there ts a substantial amount of experimental data that can be compared to our computations, allowing for validation of the model. In addition, comparison of the simulation results with the experimentally obtained plasma parameters will provide insight into the mechanisms that lead to plasma detachment, revealing how the 3 scale with different input parameters.

Polzin, K. A.↗
NASA NTRS1995

Dual-Code Solution Strategy for Chemically-Reacting Hypersonic Flows

A new procedure seeks to combine the thin-layer Navier-Stokes solver LAURA with the parabolized Navier-Stokes solver UPS for the aerothermodynamic solution of chemically-reacting air flow fields. The interface protocol is presented and the method is applied to two slender, blunted shapes. Both axisymmetric and three-dimensional solutions are included with surface pressure and heat transfer comparisons between the present method and previously published results. The case of Mach 25 flow over an axisymmetric six degree sphere-cone with a non-catalytic wall is considered to 100 nose radii. A stability bound on the marching step size was observed with this case and is attributed to chemistry effects resulting from the non-catalytic wall boundary condition. A second case with Mach 28 flow over a sphere-cone-cylinder-flare configuration is computed at both two and five degree angles of attack with a fully-catalytic wall. Surface pressures are seen to be within five percent with the present method compared to the baseline LAURA solution and heat transfers are within 10 percent. The effect of grid resolution is investigated in both the radial and streamwise directions. The procedure demonstrates significant, order of magnitude reductions in solution time and required memory for the three-dimensional case in comparison to an all thin-layer Navier-Stokes solution.

Wood, William A.↗
NASA NTRS2009

Numerical Investigation of Plasma Flows in Magnetic Nozzles

Magnetic nozzles are used in many laboratory experiments in which plasma flows are to be confined, cooled, accelerated, or directed. At present, however, there is no generally accepted theoretical description that explains the phenomena of plasma expansion in and detachment from an externally-imposed magnetic field. The latter is an especially important problem in the field of plasma propulsion, where the ionized gas must detach from the applied, solenoidal magnetic field to realize thrust production. In this paper we simulate a plasma flowing in the presence of an applied magnetic field using a multidimensional numerical simulation tool that includes theoretical models of the various dispersive and dissipative processes present in the plasma. This is an extension of the simulation tool employed in previous work by Sankaran et al. The new tool employs the same formulation of the governing equation set, but retains the axial and radial components of magnetic field and the azimuthal component of velocity that were neglected. We aim to compare the computational results with the various proposed magnetic nozzle detachment theories to develop an understanding of the physical mechanisms that cause detachment. An applied magnetic field topology is obtained using a magnetostatic field solver, and this field is superimposed on the time-dependent magnetic field induced in the plasma to provide a self-consistent field description. The applied magnetic field and model geometry match those found in experiments by Kuriki and Okada. 4 A schematic showing the setup used in those experiments is shown. We model this geometry because there is a substantial amount of experimental data that can be compared to our computations, allowing for validation of the model. In addition, comparison of the simulation results with the experimentally obtained plasma parameters will provide insight into the mechanisms that lead to plasma detachment, revealing how they scale with different input parameters.

Sankaran, Kamesh↗
NASA NTRS1999

Ames Optimized TCA Configuration

Configuration design at Ames was carried out with the SYN87-SB (single block) Euler code using a 193 x 49 x 65 C-H grid. The Euler solver is coupled to the constrained (NPSOL) and the unconstrained (QNMDIF) optimization packages. Since the single block grid is able to model only wing-body configurations, the nacelle/diverter effects were included in the optimization process by SYN87's option to superimpose the nacelle/diverter interference pressures on the wing. These interference pressures were calculated using the AIRPLANE code. AIRPLANE is an Euler solver that uses a unstructured tetrahedral mesh and is capable of computations about arbitrary complete configurations. In addition, the buoyancy effects of the nacelle/diverters were also included in the design process by imposing the pressure field obtained during the design process onto the triangulated surfaces of the nacelle/diverter mesh generated by AIRPLANE. The interference pressures and nacelle buoyancy effects are added to the final forces after each flow field calculation. Full details of the (recently enhanced) ghost nacelle capability are given in a related talk. The pseudo nacelle corrections were greatly improved during this design cycle. During the Ref H and Cycle 1 design activities, the nacelles were only translated and pitched. In the cycle 2 design effort the nacelles can translate vertically, and pitch to accommodate the changes in the lower surface geometry. The diverter heights (between their leading and trailing edges) were modified during design as the shape of the lower wing changed, with the drag of the diverter changing accordingly. Both adjoint and finite difference gradients were used during optimization. The adjoint-based gradients were found to give good direction in the design space for configurations near the starting point, but as the design approached a minimum, the finite difference gradients were found to be more accurate. Use of finite difference gradients was limited by the CPU time limit available on the Cray machines. A typical optimization run using finite difference gradients can use only 30 to 40 design variables and one optimization iteration within the 8 hour queue limit for the chosen grid size and convergence level. The efficiency afforded by the adjoint method allowed for 50-120 design variables and 5-10 optimization iterations in the 8 hour queue. Geometric perturbations to the wing and fuselage were made using the Hicks/Henne (HH) shape functions. The HH functions were distributed uniformly along the chords of the wing defining sections and lofted linearly. During single-surface design, constraints on thickness and volume at selected wing stations were imposed. Both fuselage camber and cross-sectional area distributions were permitted to change during design. The major disadvantage to the use of these functions is the inherent surface waviness produced by repeated use of such functions. Many smoothing operations were required following optimization runs to produce a configuration with reasonable smoothness. Wagner functions were also used on the wing sections but were never used on the fuselage. The Wagner functions are a family of increasingly oscillatory functions that have also been used extensively in airfoil design. The leading and trailing edge regions of the wing were designed by use of polynomial and monomial functions respectively. Twist was attempted but was abandoned because of little performance improvement available from changing the baseline twist.

Cliff, Susan E.↗
NASA NTRS2016

Applied Time Domain Stability Margin Assessment for Nonlinear Time-Varying Systems

The baseline stability margins for NASA's Space Launch System (SLS) launch vehicle were generated via the classical approach of linearizing the system equations of motion and determining the gain and phase margins from the resulting frequency domain model. To improve the fidelity of the classical methods, the linear frequency domain approach can be extended by replacing static, memoryless nonlinearities with describing functions. This technique, however, does not address the time varying nature of the dynamics of a launch vehicle in flight. An alternative technique for the evaluation of the stability of the nonlinear launch vehicle dynamics along its trajectory is to incrementally adjust the gain and/or time delay in the time domain simulation until the system exhibits unstable behavior. This technique has the added benefit of providing a direct comparison between the time domain and frequency domain tools in support of simulation validation. This technique was implemented by using the Stability Aerospace Vehicle Analysis Tool (SAVANT) computer simulation to evaluate the stability of the SLS system with the Adaptive Augmenting Control (AAC) active and inactive along its ascent trajectory. The gains for which the vehicle maintains apparent time-domain stability defines the gain margins, and the time delay similarly defines the phase margin. This method of extracting the control stability margins from the time-domain simulation is relatively straightforward and the resultant margins can be compared to the linearized system results. The sections herein describe the techniques employed to extract the time-domain margins, compare the results between these nonlinear and the linear methods, and provide explanations for observed discrepancies. The SLS ascent trajectory was simulated with SAVANT and the classical linear stability margins were evaluated at one second intervals. The linear analysis was performed with the AAC algorithm disabled to attain baseline stability margins. At each time point, the system was linearized about the current operating point using Simulink's built-in solver. Each linearized system in time was evaluated for its rigid-body gain margin (high frequency gain margin), rigid-body phase margin, and aero gain margin (low frequency gain margin) for each control axis. Using the stability margins derived from the baseline linearization approach, the time domain derived stability margins were determined by executing time domain simulations in which axis-specific incremental gain and phase adjustments were made to the nominal system about the expected neutral stability point at specific flight times. The baseline stability margin time histories were used to shift the system gain to various values around the zero margin point such that a precise amount of expected gain margin was maintained throughout flight. When assessing the gain margins, the gain was applied starting at the time point under consideration, thereafter following the variation in the margin found in the linear analysis. When assessing the rigid-body phase margin, a constant time delay was applied to the system starting at the time point under consideration. If the baseline stability margins were correctly determined via the linear analysis, the time domain simulation results should contain unstable behavior at certain gain and phase values. Examples will be shown from repeated simulations with variable added gain and phase lag. Faithfulness of margins calculated from the linear analysis to the nonlinear system will be demonstrated.

Kiefer, J. M.↗
NASA NTRS2006

Reentry-Vehicle Shape Optimization Using a Cartesian Adjoint Method and CAD Geometry

A DJOINT 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. 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 (e.g., geometric parameters that control the shape). Classic aerodynamic applications of gradient-based optimization include the design of cruise configurations for transonic and supersonic flow, as well as the design of high-lift systems. are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric computer-aided design (CAD). In previous work on Cartesian adjoint solvers, Melvin et al. 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 two-dimensional Euler equations using a ghost-cell method to enforce the wall boundary conditions. In Refs. 18 and 19, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm were the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The accuracy of the gradient computation was verified using several three-dimensional test cases, which included design variables such as the free stream parameters and the planform shape of an isolated wing. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Factors under consideration include the computation of mesh sensitivities that provide a reliable approximation of the objective function gradient, as well as the computation of surface shape sensitivities based on a direct-CAD interface. We present detailed gradient verification studies and then focus on a shape optimization problem for an Apollo-like reentry vehicle. The goal of the optimization is to enhance the lift-to-drag ratio of the capsule by modifying the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design.

Nemec, Marian↗
NASA NTRS1989

Numerical simulation of rotorcraft

The objective of the research is to develop and validate accurate, user-oriented viscous CFD codes (with inviscid options) for three-dimensional, unsteady aerodynamic flows about arbitrary rotorcraft configurations. Unsteady, three-dimensional Euler and Navier-Stokes codes are developed, adapted, and extended to rotor-body combinations. Flow solvers are coupled with zonal grid topologies, including rotating and nonrotating blocks. Special grid clustering and wave-fitting techniques were developed to capture low-level radiating acoustic waves. Significant progress was made in computing the propagation of acoustic waves due to the interaction of a concentrated vortex and a helicopter airfoil. The need for higher-order schemes was firmly established in relatively inexpensive two-dimensional calculations. In three dimensions, the number of grid points required to capture the low-level acoustic waves becomes very large, so that large supercomputer memory becomes essential. Good agreement was obtained between the numerical results obtained with a thin-layer Navier-Stokes code and experimental data from a model rotor. In addition, several nonrotating configurations that are sometimes proposed to simulate rotor blade tips in conventional wind tunnels were examined, and the complex flow around the radical tip shape of the world's fastest helicopter is under investigation. These studies demonstrate the flexibility and power of CFD to gain physical insight, study novel ideas, and examine various possibilities that might be difficult or impossible to set up in physical experiments. As a prelude to studies of rotor-body aerodynamic interactions, a preliminary grid topology and moving-interface strategy were developed. A new Euler/Navier-Stokes code using these techniques computes the vortical wake directly, rather than modeling it, as in most previous rotorcraft studies. Several hover cases were run for conventional and advanced-geometry blades. Numerical schemes using multi-zones and/or adaptive grids appear to be necessary to simulate the complex vortical flows in rotor wakes.

Mccroskey, William J.↗
NASA NTRSDate not supplied

Wall-modeled LES of the NASA Juncture Flow Experiment

We perform wall-modeled large eddy simulation (WMLES) of the NASA Juncture Flow experiment (Kegerise & Neuhart, NASA/TM–2019–20286) using an equilibrium wall model and unstructured finite volume solver, to assess its predictive ability for complex geometries. The flow condition simulated corresponds to 5 degrees angle of attack with a freestream Mach number of 0.189, and a Reynolds number based on the mean crank chord of 2.4 million. All simulations neglect the effects of wind tunnel walls, sting and mast present in the experiment. To assess the performance of WMLES on simpler flows with very coarse grids, we simulate flow in a turbulent channel at a friction Reynolds number, Reτ≈2000, and flow through a square duct at anReτ≈1000. The duct flow results using the equilibrium wall model indicate that the skin friction is inaccurate near the corner for grid topologies and resolutions typically used in WMLES, having implications for more complicated flows such as the juncture flow. For the juncture flow simulation, two different approaches are investigated. The first approach consists of a truncated-domain simulation wherein the inflow plane of the computational domain is placed at≈0.2chord length of the wing, and we prescribe the mean flow from a separate Reynolds-averaged Navier-Stokes (RANS) solution along with synthetic turbulence to initiate realistic unsteadiness in the domain. The second approach involves simulating the entire geometry with trips to trigger transition to turbulence. The truncated- and full-domain simulations contain about 62 and 90 million cells, respectively; with 8-10 points per boundary layer thickness in the wall-parallel directions, 14-20 points in the wall-normal direction and a near-wall viscous spacing (∆n+1) of≈100. Preliminary results are encouraging overall in terms of the prediction of wall pressure, wall skin friction, velocity and stresses; but indicate that further work is required to improve the predictions in the separation bubble and wing-fuselage corner regions. The grid resolution used in this study is still fairly coarse, and the results should be interpreted as work in progress.

Prahladh S Iyer↗
NASA NTRS2003

Computational Evaluation of the Steady and Pulsed Jet Effects on the Performance of a Circulation Control Wing Section

Circulation Control Wing (CCW) technology is a very effective way of achieving very high lift coefficients needed by aircraft during take-off and landing. This technology can also be used to directly control the flow field over the wing. Compared to a conventional high-lift system, a Circulation Control Wing (CCW) can generate the required values of lift coefficient C(sub L,max) during take-off/landing with fewer or no moving parts and much less complexity. Earlier designs of CCW configurations used airfoils with a large radius rounded trailing edge to maximize the lift benefit. However, these designs also produced very high drag. These high drag levels associated with the blunt, large radius trailing edge can be prohibitive under cruise conditions when Circulation Control is no longer necessary. To overcome this difficulty, an advanced CCW section, i.e., a circulation hinged flap was developed to replace the original rounded trailing edge CC airfoil. This concept developed by Englar is shown. The upper surface of the CCW flap is a large-radius arc surface, but the lower surface of the flap is flat. The flap could be deflected from 0 degrees to 90 degrees. When an aircraft takes-off or lands, the flap is deflected as in a conventional high lift system. Then this large radius on the upper surface produces a large jet turning angle, leading to high lift. When the aircraft is in cruise, the flap is retracted and a conventional sharp trailing edge shape results, greatly reducing the drag. This kind of flap does have some moving elements that increase the weight and complexity over an earlier CCW design. But overall, the hinged flap design still maintains most of the Circulation Control high lift advantages, while greatly reducing the drag in cruising condition associated with the rounded trailing edge CCW design. In the present work, an unsteady three-dimensional Navier-Stokes analysis procedure has been developed and applied to this advanced CCW configuration. The solver can be used in both a 2-D and a 3-D mode, and can thus model airfoils as well as finite wings. The jet slot location, slot height, and the flap angle can all be varied easily and individually in the grid generator and the flow solver. Steady jets, pulsed jets, the leading edge and trailing edge blowing can all be studied with this solver.

Liu, Yi↗
NASA NTRS1991

An Explicit Upwind Algorithm for Solving the Parabolized Navier-Stokes Equations

An explicit, upwind algorithm was developed for the direct (noniterative) integration of the 3-D Parabolized Navier-Stokes (PNS) equations in a generalized coordinate system. The new algorithm uses upwind approximations of the numerical fluxes for the pressure and convection terms obtained by combining flux difference splittings (FDS) formed from the solution of an approximate Riemann (RP). The approximate RP is solved using an extension of the method developed by Roe for steady supersonic flow of an ideal gas. Roe's method is extended for use with the 3-D PNS equations expressed in generalized coordinates and to include Vigneron's technique of splitting the streamwise pressure gradient. The difficulty associated with applying Roe's scheme in the subsonic region is overcome. The second-order upwind differencing of the flux derivatives are obtained by adding FDS to either an original forward or backward differencing of the flux derivative. This approach is used to modify an explicit MacCormack differencing scheme into an upwind differencing scheme. The second order upwind flux approximations, applied with flux limiters, provide a method for numerically capturing shocks without the need for additional artificial damping terms which require adjustment by the user. In addition, a cubic equation is derived for determining Vegneron's pressure splitting coefficient using the updated streamwise flux vector. Decoding the streamwise flux vector with the updated value of Vigneron's pressure splitting improves the stability of the scheme. The new algorithm is applied to 2-D and 3-D supersonic and hypersonic laminar flow test cases. Results are presented for the experimental studies of Holden and of Tracy. In addition, a flow field solution is presented for a generic hypersonic aircraft at a Mach number of 24.5 and angle of attack of 1 degree. The computed results compare well to both experimental data and numerical results from other algorithms. Computational times required for the upwind PNS code are approximately equal to an explicit PNS MacCormack's code and existing implicit PNS solvers.

Korte, John J.↗
NASA NTRS1999

Performance and Application of Parallel OVERFLOW Codes on Distributed and Shared Memory Platforms

The presentation discusses recent studies on the performance of the two parallel versions of the aerodynamics CFD code, OVERFLOW_MPI and _MLP. Developed at NASA Ames, the serial version, OVERFLOW, is a multidimensional Navier-Stokes flow solver based on overset (Chimera) grid technology. The code has recently been parallelized in two ways. One is based on the explicit message-passing interface (MPI) across processors and uses the _MPI communication package. This approach is primarily suited for distributed memory systems and workstation clusters. The second, termed the multi-level parallel (MLP) method, is simple and uses shared memory for all communications. The _MLP code is suitable on distributed-shared memory systems. For both methods, the message passing takes place across the processors or processes at the advancement of each time step. This procedure is, in effect, the Chimera boundary conditions update, which is done in an explicit "Jacobi" style. In contrast, the update in the serial code is done in more of the "Gauss-Sidel" fashion. The programming efforts for the _MPI code is more complicated than for the _MLP code; the former requires modification of the outer and some inner shells of the serial code, whereas the latter focuses only on the outer shell of the code. The _MPI version offers a great deal of flexibility in distributing grid zones across a specified number of processors in order to achieve load balancing. The approach is capable of partitioning zones across multiple processors or sending each zone and/or cluster of several zones into a single processor. The message passing across the processors consists of Chimera boundary and/or an overlap of "halo" boundary points for each partitioned zone. The MLP version is a new coarse-grain parallel concept at the zonal and intra-zonal levels. A grouping strategy is used to distribute zones into several groups forming sub-processes which will run in parallel. The total volume of grid points in each group are approximately balanced. A proper number of threads are initially allocated to each group, and in subsequent iterations during the run-time, the number of threads are adjusted to achieve load balancing across the processes. Each process exploits the multitasking directives already established in Overflow.

Djomehri, M. Jahed↗
NASA NTRSDate not supplied

Magnetohydrodynamics (MHD) Aerocapture System for Enabling Faster-Larger Planetary Science & Human Exploration Missions

Since our completing the NIAC Phase I NIAC on this Advanced Aerocapture System, NASA Langley Research Center has funded or supported a number of studies and code enhancements through its Center Innovation Fund (CIF) and NASA’s NSTGRO and Internship Programs to mature the analysis capabilities and quantify the merits of the MHD Aerocapture System technology. These efforts have resulted in a plug and play analysis capability for assessing MHD aerocapture system performance for arrival at many planetary bodies of interest. Our efforts have especially focused on the potential mass savings for improving the capacity for science observations at Neptune and Triton. A re-cent Forbes article published “‘Orbital mechanics is probably going to decide for us whether we go to Uranus or Neptune because we need to flyby Jupiter,’ said Kunio Sayanagi at Hampton University, Virginia, who also worked on the Neptune Odyssey proposal…. Exactly when a mission can be sent to Uranus, or Neptune, depends on the relative position of Jupiter, which can help give a spacecraft a gravitational slingshot. That drastically shortens the cruise phase.” [1] Since shortening the cruise phase is important for these science missions, any mass savings enabled by the MHD Aerocapture System could be reallocated to increasing Thermal Protection System mass to allow faster arrival speeds and/or for onboarding additional payloads such as science instruments, batteries, or propellant for conducting more science for longer durations in the desired orbits. The analysis steps and codes for conducting trades and sizing vehicles for aerocapture are as follows: Step 1 is to conduct aeroheating analysis using LAURA of the selected entry vehicle shape to identify locations on the forebody where ionization and flow velocity are sufficient for producing Lo-rentz forces. LAURA is a multiblock structured grid finite-volume CFD solver developed at the NASA Langley Research Center. [2] LAURA has been used for aerothermal analysis support of the entry, de-scent and landing (EDL) phase of interplanetary missions over the last three decades [3-7]. Step 2 is to port the LAURA results into CFDWARP to calcu-late electrical and thermal conductivities of ionized flow for sizing MHD patch system and calculating Lorentz forces needed for controls analysis. CFDWARP is a CFD code that uses advanced nu-merical methods that enable the simulation of the full coupling between the aerodynamics, the magne-tohydrodynamics, and the non-neutral plasma sheaths. CFDWARP has the unique capability to simulate efficiently the non-neutral sheaths (near the electrodes) in coupled form with the quasi-neutral bulk MHD flow [8-11]. Step 3 is to link re-sults from LAURA and CFDWARP into POST2 for calculating entry trajectories and comparing MHD control results with other aerodynamic control strategies. The Program to Optimize Simulated Tra-jectories II (POST2) is a generalized point mass, discrete parameter targeting and optimization pro-gram. POST2 provides the capability to target and optimize point mass trajectories for multiple pow-ered or un-powered vehicles near an arbitrary rotat-ing, oblate planet [12]. Step 4: TPS sizing was per-formed using the Fully Implicit Ablation and Ther-mal-response code (FIAT) tool which computes the transient one-dimensional thermal response and surface thermochemistry of a multilayer stackup of thermal protection, bonding, and structural materi-als subject to aeroheating on one surface [13]. The sizing and margining methodology used was based on the approach documented by Mahzari and Milos [14] for the dual-layered heatshield for extreme entry environment technology (DL-HEEET) TPS concept. TPS analysis utilizes trajectory information from POST2. Using this step-wise plug and play MHD Aerocapture performance assessment process, our analysis targets a Neptune aerocapture trajectory that will place the spacecraft in an observation orbit for Triton. [15]. Magnetohydrodynamic (MHD) control of a 4.5-meter diameter MSL-style capsule resulted in TPS mass savings of nearly 2000 kg when using an MHD system mass of under 200 kg. The flight path for a vehicle using the MHD control strategy has a much lower heat rate and heat load compared to the conventional aerodynamic aerocapture strategies known as bank angle con-trolled (BAC) and direct force controlled (DFC). Both BAC and DFC have heat rates significantly greater than 1500 W/cm2 typically used as an upper limit for PICA. Thus, DL-HEEET TPS concept was required for the BAC and DFC control strategies. However, considering the more benign environ-ments for the MHD case, additional TPS concepts with improved mass efficiency were also assessed. PICA was considered for the MHD controlled strat-egy since the maximum heat rate was well within the limits (<1500 W/cm2) of PICA. TPS sizing re-sulted in a significant mass reduction. The PICA layer for this sizing case was about 7.8 cm. As a point of reference, the Mars 2020 mission, which used this same PICA concept, had a PICA thickness of 3.18 cm [16]. The trajectories used for the TPS sizing originat-ed from the POST2 simulations. The current, I, to an electromagnet configuration can be manipulated to allow for active control of the vehicle. Manipula-tion of the current, I, changes the magnetic field, B, which affects the Lorentz force and therefore the MHD drag force on the vehicle. Our analysis in-cluded both open-loop and close-loop control. Closed-loop control will enable improved overall performance when taking into account mission level uncertainties, such as interplanetary delivery errors and atmospheric modeling uncertainties. The open-loop and closed-loop MHD control cases do not dip as deep into the atmosphere as the aerodynamic cases. Three types of aerodynamic-only approach-es are investigated: bank angle modulation (BAM), director force control (DFC), and Drag Modulated. BAM and DFC make use of vehicle aerodynamic angles to steer the vehicle. Thus, changing the aer-odynamic forces acting on the vehicle for control, aerodynamic drag modulated case requires a vary-ing drag area to modulate the drag force. The MHD drag modulated case modulates MHD generated drag force that adds to the aerodynamic drag. This higher atmospheric activation of drag forces by the MHD patch results in significantly less heat flux on the vehicle. The MHD technology will enable shorter cruise times and deceleration of larger payloads for increasing the capacity for science at the Ice Giants or for returning astronauts to Earth from cislunar space or from Mars. The purpose of this presentation is to provide more details about this work and to highlight plans for further research and development including a flight demonstration.

R. W. Moses↗
NASA NTRS2019

Supporting Crew Autonomy in Deep Space Exploration: Preliminary Onboard Capability Requirements and Proposed Research Questions. Technical Report of the Autonomous Crew Operations Technical Interchange Meeting

Communication delays are a critical challenge posed by long duration deep space exploration. Space missions historically have relied on an ever-present Mission Control Center (MCC) to direct operations in near real-time. As unanticipated anomalies that defeat fault detection and resolution systems do arise, the lack of real-time communication will significantly weaken what the MCC support represents: a reliable safety net for the flight crew through its deep and diverse areas of expertise and investigative resources. As a consequence, future space vehicles and habitats need to be equipped with capabilities to support the flight crew to operate with little or no ground support. Considerations must be given to vehicle and mission designs that will fortify the traditionally ground-centered safety net and forge new support systems, when communication delays exist. In August 2018, NASA’s Human Research Program, through its Human Factors and Behavioral Performance Element, convened a Technical Interchange Meeting (TIM) on Autonomous Crew Operations at NASA Ames Research Center. The goal of the meeting was to gather input from NASA centers, industry, academia, and branches of the Department of Defense (DoD) to address how intelligent technologies can be applied to augment onboard capabilities to support crew anomaly response. The TIM featured 24 presentations by 29 speakers and hosted a total of 59 attendees, including 43 from 5 NASA centers (Ames, Johnson, Langley, Marshall, and Jet Propulsion Lab) and 4 from the DoD (3 from Army Research Lab and 1 from Naval Postgraduate School), with remaining attendees from academia (e.g., UC Davis, CMU) and industry (e.g., IBM, Siemens). Discussions were centered around three themes: standards and guidelines, lessons learned in analog environments, and technologies. To help provide a framework for discussion, a concept matrix describing anomaly response processes was created prior to the TIM (Figure 1, page 6). The matrix captures the steps involved (monitoring and detection, diagnosis, solution development and evaluation, solution implementation and verification, resolution documentation) as well as the resources and capabilities required to support these steps (data, knowledge, analysis, synthesis, resource management). A wallpaper size printout of the matrix was utilized at the TIM to solicit attendee inputs along the three themes; the activity garnered 108 submissions of ideas. Overall, what emerged from TIM discussions was a picture of mismatch between crew anomaly response needs and support that can be provided by existing intelligent technologies. The needs are broad, spanning multiple steps and processes/resources, with many of which lacking support from existing technologies, such as knowledge management throughout the steps of problem solving (especially in resolution documentation) and manpower management. The solutions provided by existing intelligent technologies are specific to the steps/processes that they are designed to support and constrained to solving only problems similar to those that have occurred before. What is lacking from technologies is typically made up by humans, specifically their complex critical thinking, creative problem solving, and domain expertise. In the end, the TIM highlighted the pressing need to support responses to onboard anomalies during autonomous crew operations, particularly those that have eluded the system tests, inspection, and other assurance processes. Such anomalies can potentially threaten crew and vehicle safety, as well as significantly impact overall operations with additional workload. These fairly rare events are difficult to anticipate and prepare for, given the state-of-the-art in intelligent technologies. This is true even for anomalies that stem from “unknown knowns”—cases in which there is sufficient external information to characterize the problem but the overall pattern fails to be recognized by the problem solver, or in which the internal knowledge needed to solve a problem is held tacitly and potentially accessible by the problem solver but not articulated. It follows that the ability to tackle anomalies lies not only with the availability of relevant information and knowledge but also their accessibility in times of need. To that end, we propose research questions along the following three broad themes: • How intelligent technologies can help make relevant knowledge and information available? • How intelligent technologies can help make relevant knowledge and information accessible? • How intelligent technologies can help support the crew operating as a team in anomaly response processes?

autonomous crew operations↗