SEARCH · Engineering Papers
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.
Computer program determines the alternative geometries and associated design point performance of axial-flow turbines capable of satisfying specified design requirements. The program solves the flow field within the turbine without making the simplifying assumptions that result in restrictive designs.
Not Available
Computer program calculates leading- and trailing-edge circle radii, tangency angles on leading- and trailing-edge circles, and stagger angle of turbomachinery blade sections, using only spline points defining blade surfaces.
Explore the source record for details and available documents.
The operation of the MAST code, which computes transient solutions to the multiphase flow equations applicable to all-speed flows, is described. Two-phase flows are formulated based on the Eulerian-Lagrange scheme in which the continuous phase is described by the Navier-Stokes equation (or Reynolds equations for turbulent flows). Dispersed phase is formulated by a Lagrangian tracking scheme. The numerical solution algorithms utilized for fluid flows is a newly developed pressure-implicit algorithm based on the operator-splitting technique in generalized nonorthogonal coordinates. This operator split allows separate operation on each of the variable fields to handle pressure-velocity coupling. The obtained pressure correction equation has the hyperbolic nature and is effective for Mach numbers ranging from the incompressible limit to supersonic flow regimes. The present code adopts a nonstaggered grid arrangement; thus, the velocity components and other dependent variables are collocated at the same grid. A sequence of benchmark-quality problems, including incompressible, subsonic, transonic, supersonic, gas-droplet two-phase flows, as well as spray-combustion problems, were performed to demonstrate the robustness and accuracy of the present code.
Cartesian methods for computational fluid dynamics are considered to offer an accurate and robust approach for the simulation of aerodynamic flows around geometrically complex bodies. A part of the ongoing research in this domain is reviewed with the aim of providing insight into the fundamental challenges faced by the practitioners of this approach, and a guide to further research. The integration schemes used in Cartesian solvers are similar to those used in other approaches. The emphasis is on the geometric algorithms, surface modeling and boundary conditions required to design a successful Cartesian mesh scheme.
The objective of the effort was to numerically simulate viscous subsonic flow in a proposed elliptical two-duct version of the fuel side Hot Gas Manifold (HGM) for the Space Shuttle Main Engine (SSME). The numerical results were to complement both water flow and air flow experiments in the two-duct geometry performed at NASA-MSFC and Rocketdyne. The three-dimensional character of the HGM consists of two essentially different geometries. The first part of the construction is a concentric shell duct structure which channels the gases from a turbine exit into the second part comprised of two cylindrically shaped transfer ducts. The initial concentric shell portion can be further subdivided into a turnaround section and a bowl section. The turnaround duct (TAD) changes the direction of the mean flow by 180 degress from a smaller radius to a larger radius duct which discharges into the bowl. The cylindrical transfer ducts are attached to the bowl on one side thus providing a plane of symmetry midway between the two. Centerline flow distance from the TAD inlet to the transfer duct exit is approximately two feet. Details of the approach used to numerically simulate laminar or turbulent flow in the HGM geometry are presented. Computational results are presented and discussed.
Abstract Heat exchangers play a critical role in supercritical CO2 Brayton cycles by providing necessary waste heat recovery. Supercritical CO2 thermal cycles potentially achieve higher energy density and thermal efficiency operating at elevated temperatures and pressures. Accurate and computationally efficient estimation of heat exchanger performance metrics at these conditions is important for the design and optimization of sCO2 systems and thermal cycles. In this paper (Part II), a computationally efficient and accurate numerical model is developed to predict the performance of shell-and-tube heat exchangers (STHXs). Highly accurate correlations reported in Part I of this study are utilized to improve the accuracy of performance predictions, and the concept of volume averaging is used to abstract the geometry and reduce computation time. The numerical model is validated by comparison with computational fluid dynamics (CFD) simulations and provides high accuracy and significantly lower computation time compared to existing numerical models. A preliminary optimization study is conducted and the advantage of using supercritical CO2 as a working fluid for energy systems is demonstrated.
A joint workshop between the Air Force Office of Scientific Research, the Office of Naval Research, and the NASA Hypersonic Technology Project was held in June 2023 on the topic of high-speed boundary-layer transition. The purpose of this workshop was to gather subject matter experts in the areas of supersonic and hypersonic boundary-layer transition in the areas of geometries of interest, computational methods, facilities and diagnostics, flight experiments, flow control, and new approaches and to discuss the current state of the art in the field. As gaps and accomplishments were identified, roadmaps were also created to guide future investments in these areas. The top recommendations are (1) to define relevant common geometries of interest that can be worked by multiple research groups; (2) improve modeling of high-enthalpy flows involving chemistry and ablation, studies of receptivity of freestream disturbances and atmospheric particulates into the boundary layer, surface feature characterization and modeling of changes to geometry; (3) design and determine viable flow control methods; and (4) move towards incorporating boundary layer transition prediction into a standard practice for design of high-speed vehicles. To execute this last recommendation successfully, investments beyond fundamental (basic and early applied) research are needed.
Herein, we present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a shock-capturing Finite Volume (FV) scheme coupled via an hp adaptive mesh refinement (hp-AMR) strategy that offers high-order accurate resolution of the embedded geometries. The hp-AMR strategy is based on a multi-level block-structured domain partition in which each level is represented by block-structured Cartesian grids and the embedded geometry is represented implicitly by a level set function. The intersection of the embedded geometry with the grids produces the implicitly-defined mesh that consists of a collection of regular rectangular cells plus a relatively small number of irregular curved elements in the vicinity of the embedded boundaries. High-order quadrature rules for implicitly-defined domains enable high-order accuracy resolution of the curved elements with a cell-merging strategy to address the small-cell problem. The hp-AMR algorithm treats the system with a second-order finite volume scheme at the finest level to dynamically track the evolution of solution discontinuities while using dG schemes at coarser levels to provide high-order accuracy in smooth regions of the flow. On the dG levels, the methodology supports different orders of basis functions on different levels. The space-discretized governing equations are then advanced explicitly in time using high-order Runge-Kutta algorithms. Numerical tests are presented for two-dimensional and three-dimensional problems involving an ideal gas. The results are compared with both analytical solutions and experimental observations and demonstrate that the framework provides high-order accuracy for smooth flows and accurately captures solution discontinuities.
Porous electrodes are used as the core reactive component across electrochemical technologies. In flowing systems, controlling the fluid distribution, species transport, and reactive environment is critical to attaining high performance. However, conventional electrode materials like felts and papers provide few opportunities for precise engineering of the electrode and its microstructure. To address these limitations, architected electrodes composed of unit cells with spatially varying geometry determined via computational optimization are proposed. Resolved simulation is employed to develop a homogenized description of the constituent unit cells. These effective properties serve as inputs to a continuum model for the electrode when used in the negative half-cell of a vanadium redox flow battery. Porosity distributions minimizing power loss are then determined via computational design optimization to generate architected porosity electrodes. The architected electrodes are compared to bulk, uniform porosity electrodes and found to lead to increased power efficiency across operating flow rates and currents. The design methodology is further used to generate a scaled-up electrode with comparable power efficiency to the bench-scale systems. Finally, the variable porosity architecture and computational design methodology presented here thus offers a novel pathway for automatically generating spatially engineered electrode structures with improved power performance.
A solid freeform fabrication process and apparatus for making a three-dimensional reinforcement shape. The process comprises the steps of (1) operating a multiple-channel material deposition device for dispensing a liquid adhesive composition and selected reinforcement materials at predetermined proportions onto a work surface; (2) during the material deposition process, moving the deposition device and the work surface relative to each other in an X-Y plane defined by first and second directions and in a Z direction orthogonal to the X-Y plane so that the materials are deposited to form a first layer of the shape; (3) repeating these steps to deposit multiple layers for forming a three-dimensional preform shape; and (4) periodically hardening the adhesive to rigidize individual layers of the preform. These steps are preferably executed under the control of a computer system by taking additional steps of (5) creating a geometry of the shape on the computer with the geometry including a plurality of segments defining the preform shape and each segment being preferably coded with a reinforcement composition defining a specific proportion of different reinforcement materials; (6) generating programmed signals corresponding to each of the segments in a predetermined sequence; and (7) moving the deposition device and the work surface relative to each other in response to these programmed signals. Preferably, the system is also operated to generate a support structure for any un-supported feature of the 3-D preform shape.
Calculating the space-charge-limited-current density (SCLCD) for a complicated diode geometry often requires computationally expensive particle-in-cell (PIC) codes. Here, this paper addresses this issue by using the charge-free electric field $E_0$ calculated using COMSOL Multiphysics to determine local and global SCLCD. The SCLCD obtained by using the surface average of $|E_0|^2$ on the cathode recovers theoretical results for one-dimensional (1D) planar, cylindrical, and tip-to-tip geometries in appropriate limits. We next compared tip-to-tip calculations with the SCLCD obtained using the PIC code Empire. The SCLCD calculated using COMSOL agreed well with Empire for flatter 1D tip-to-tip geometries and diverged with increasing sharpness. Physically, Empire predicts lower SCLCD than COMSOL because the electrons spread due to concentrated space-charge at the tip, whereas theory assumes that the electrons follow the charge-free electric field lines. We further assess the behavior of the SCLCD for tips protruding from the centers of flat, circular plates of various areas. Larger plate areas with constant tip size recover the 1D planar SCLCD globally and 1D tip-to-tip SCLCD locally, while reducing the difference between Empire and COMSOL calculations since larger plates capture more of the emitted electrons, reducing SCLCD suppression due to beam spreading. These results show that charge-free electric field simulations can be used to determine the SCLCD without needing to simulate particle dynamics in PIC.
Paper will present recent developments in advanced material binder systems for solid freeform fabrication (SFF) technologies. The advantage of SFF is the capability to custom fabricate complex geometries directly from computer aided design data in layer- by-layer fashion, eliminated the need for traditional fixturing and tooling. Binders allow for the low temperature processing of 'green' structural materials, either metal, ceramic or composite, in traditional rapid prototyping machines. The greatest obstacle comes when green parts must then go through a sintering or burnout process to remove the binders and fully densify the parent material, without damaging or distorting the original part geometry. Critical issues and up-to-date assessments will be delivered on various material systems.
Preliminary verification and validation of an efficient Euler solver for adaptively refined Cartesian meshes with embedded boundaries is presented. The parallel, multilevel method makes use of a new on-the-fly parallel domain decomposition strategy based upon the use of space-filling curves, and automatically generates a sequence of coarse meshes for processing by the multigrid smoother. The coarse mesh generation algorithm produces grids which completely cover the computational domain at every level in the mesh hierarchy. A series of examples on realistically complex three-dimensional configurations demonstrate that this new coarsening algorithm reliably achieves mesh coarsening ratios in excess of 7 on adaptively refined meshes. Numerical investigations of the scheme's local truncation error demonstrate an achieved order of accuracy between 1.82 and 1.88. Convergence results for the multigrid scheme are presented for both subsonic and transonic test cases and demonstrate W-cycle multigrid convergence rates between 0.84 and 0.94. Preliminary parallel scalability tests on both simple wing and complex complete aircraft geometries shows a computational speedup of 52 on 64 processors using the run-time mesh partitioner.
The near-field behavior of a wingtip vortex flow is studied computationally using an incompressible flow solver for the Navier-Stokes equations based on the artificial compressibility method. Inaccuracies in current computational studies are addressed, especially, the role of numerical errors and transition/turbulence modeling. A subset problem is devised in order to make the study of vortex preservation more tractable. As part of the numerical checks, the flow is first run laminar while performing a systematic grid refinement study for the subset problem. Further studies on the numerical errors are conducted with the measured Reynolds stresses introduced into the momentum equations as source terms. As a preliminary study of turbulent flows, the one-equation Baldwin-Barth turbulence model is implemented as well as the approximation of the production term. The full-geometry case is computed using 1.1 million grid points. The results are compared with experiment.
Computational Fluid Dynamics (CFD) has considerably evolved in the last decade. There are many computer programs that can perform computations on viscous internal or external flows with chemical reactions. CFD has become a commonly used tool in the design and analysis of gas turbines, ramjet combustors, turbo-machinery, inlet ducts, rocket engines, jet interaction, missile, and ramjet nozzles. One of the problems of interest to NASA has always been the performance prediction for rocket and air-breathing engines. Due to the complexity of flow in these engines it is necessary to resolve the flowfield into a fine mesh to capture quantities like turbulence and heat transfer. However, calculation on a high-resolution grid is associated with a prohibitively increasing computational time that can downgrade the value of the CFD for practical engineering calculations. The Liquid Thrust Chamber Performance (LTCP) code was developed for NASA/MSFC (Marshall Space Flight Center) to perform liquid rocket engine performance calculations. This code is a 2D/axisymmetric full Navier-Stokes (NS) solver with fully coupled finite rate chemistry and Eulerian treatment of liquid fuel and/or oxidizer droplets. One of the advantages of this code has been the resemblance of its input file to the JANNAF (Joint Army Navy NASA Air Force Interagency Propulsion Committee) standard TDK code, and its automatic grid generation for JANNAF defined combustion chamber wall geometry. These options minimize the learning effort for TDK users, and make the code a good candidate for performing engineering calculations. Although the LTCP code was developed for liquid rocket engines, it is a general-purpose code and has been used for solving many engineering problems. However, the single zone formulation of the LTCP has limited the code to be applicable to problems with complex geometry. Furthermore, the computational time becomes prohibitively large for high-resolution problems with chemistry, two-equation turbulence model, and two-phase flow. To overcome these limitations, the LTCP code is rewritten to include the multi-zone capability with domain decomposition that makes it suitable for parallel processing, i.e., enabling the code to run every zone or sub-domain on a separate processor. This can reduce the run time by a factor of 6 to 8, depending on the problem.
The effect of forebody tangential slot blowing on the flowfield about an F/A-18 aircraft is investigated numerically using solutions of the Navier-Stokes equations. Computed solutions are obtained for an aircraft geometry which includes the fuselage, a wing with deflected leading-edge flap, empennage, and a faired-over engine inlet. The computational slot geometry corresponds to that used in full-scale wind-tunnel tests. Solutions are computed using flight test conditions and jet mass flow ratios equivalent to wind-tunnel test conditions. The effect of slot location is analyzed by computing two nontime-accurate solutions with a 16-in. slot located 3 in. and 11 in. aft of the nose of the aircraft. These computations resolve the trends observed in the full-scale wind-tunnel test data. The flow aft of the leading-edge extension vortex burst is unsteady. A time-accurate solution is obtained to investigate the flow characteristics aft of the vortex burst, including the effect of blowing on tail buffet.