Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational fluid dynamics code”

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 1,027 records · Page 57

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)↗

A summary of the computational fluid dynamics analysis of the hypersonic flights of Pegasus

A computational fluid dynamics (CFD) analysis of the hypersonic flights of the NASA's Pegasus launch vehicle was carried out using the PARC code. The flight data included heat transfer rates and sparse pressure coefficients primarily in the fillet region of the vehicle. Since the CFD simulations did not account for ablation at the vehicle surface, the test heat transfer results were obtained for the low Mach number simulations; as the Mach number was increased, the ablative effects at the vehicle surface became more apparent and contributed to the error.

Fricker, Darren↗

Numerical investigations in three-dimensional internal flows

In the present reporting period, the 3D STUFF code was used to solve the underbody flow for the waverider configuration vehicle. In order to start the space-marched version of the code, the time-marched version (TUFF) was used to solve the forward portion of the underside of the forebody flowfield. A grid was generated which went from the tip of the nose to the location of the cowl lip. This includes all of the inlet ramp system. Previous indications were that three-dimensional effects could be expected on the ramps of such an aircraft. For purposes of the present study, no sidewalls were assumed. The sidewalls were eliminated to simplify the calculations and to show the potential effects of three-dimensional flow in the absence of a full sidewall. Further, this flow was also analyzed using the newly released OVERFLOW code and comparisons between the two codes were made. In addition to these 3D calculations, 2D calculations using the OVERFLOW code were also obtained for the Mach 5 inlet model in this reporting period. Comparisons between the experimental data, previous computational fluid dynamics (CFD) results and those from OVERFLOW were made.

Rose, William C.↗

Numerical Analysis of 2-D and 3-D MHD Flows Relevant to Fusion Applications

Here, the analysis of many fusion applications such as liquid-metal blankets requires application of computational fluid dynamics (CFD) methods for electrically conductive liquids in geometrically complex regions and in the presence of a strong magnetic field. A current state of the art general purpose CFD code allows modeling of the flow in complex geometric regions, with simultaneous conjugated heat transfer analysis in liquid and surrounding solid parts. Together with a magnetohydrodynamics (MHD) capability, the general purpose CFD code will be a valuable tool for the design and optimization of fusion devices. This paper describes an introduction of MHD capability into the general purpose CFD code CFX, part of the ANSYS Workbench. The code was adapted for MHD problems using a magnetic induction approach. CFX allows introduction of user-defined variables using transport or Poisson equations. For MHD adaptation of the code three additional transport equations were introduced for the components of the magnetic field, in addition to the Poisson equation for electric potential. The Lorentz force is included in the momentum transport equation as a source term. Fusion applications usually involve very strong magnetic fields, with values of the Hartmann number of up to tens of thousands. In this situation a system of MHD equations become very rigid with very large source terms and very strong variable gradients. To increase system robustness, special measures were introduced during the iterative convergence process, such as linearization using source coefficient for momentum equations. The MHD implementation in general purpose CFD code was tested against benchmarks, specifically selected for liquid-metal blanket applications. Results of numerical simulations using present implementation closely match analytical solutions for a Hartmann number of up to 1500 for a 2-D laminar flow in the duct of square cross section, with conducting and nonconducting walls. Results for a 3-D test case are also included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Control Law Design in a Computational Aeroelasticity Environment

A methodology for designing active control laws in a computational aeroelasticity environment is given. The methodology involves employing a systems identification technique to develop an explicit state-space model for control law design from the output of a computational aeroelasticity code. The particular computational aeroelasticity code employed in this paper solves the transonic small disturbance aerodynamic equation using a time-accurate, finite-difference scheme. Linear structural dynamics equations are integrated simultaneously with the computational fluid dynamics equations to determine the time responses of the structure. These structural responses are employed as the input to a modern systems identification technique that determines the Markov parameters of an "equivalent linear system". The Eigensystem Realization Algorithm is then employed to develop an explicit state-space model of the equivalent linear system. The Linear Quadratic Guassian control law design technique is employed to design a control law. The computational aeroelasticity code is modified to accept control laws and perform closed-loop simulations. Flutter control of a rectangular wing model is chosen to demonstrate the methodology. Various cases are used to illustrate the usefulness of the methodology as the nonlinearity of the aeroelastic system is increased through increased angle-of-attack changes.

Newsom, Jerry R.↗

A Performance-Portable MultiGPU Implementation of 3D Euler Equations using ProtoX and IRIS

Computational scientists often face challenges when developing and optimizing code for high-performance computing (HPC), especially when trying to leverage GPUs. Given the heterogeneity of the nodes that comprise many modern HPC facilities, considerable demand exists for performance portable solutions for the core computational kernels used in many scientific computing libraries. In this work, we demonstrate a fourth-order finite volume method–based implementation of the Euler equations, which are an integral part of computational fluid dynamics. Our performance-portable multiGPU implementation for Euler equations uses ProtoX to generate kernels and IRIS for portability. ProtoX is a domain-specific language that uses a structured-grid partial differential equation library called Proto as its front end and the SPIRAL code generation system as its back end to generate optimized kernels for different architectures. Optimized kernels generated by ProtoX are orchestrated through the IRIS intelligent runtime system to provide portability. Two levels of optimizations within the IRIS runtime— directed acyclic graph fusion and task fusion—are explored to efficiently utilize computing resources in a multiGPU environment. Performance improvement through these optimizations is showcased by comparing the base ProtoX-IRIS implementation on AMD GPUs (Frontier node) and on NVIDIA GPUs (NVIDIA DGX-1).

Mankad, Het↗

Design of Optimum Ducts Using an Efficient 3-D Viscous Computational Flow Analysis

Design of fluid dynamically efficient ducts is addressed through the combination of an optimization analysis with a three-dimensional viscous fluid dynamic analysis code. For efficiency, a parabolic fluid dynamic analysis was used. Since each function evaluation in an optimization analysis is a full three-dimensional viscous flow analysis requiring 200,000 grid points, it is important to use both an efficient fluid dynamic analysis and an efficient optimization technique. Three optimization techniques are evaluated on a series of test functions. The Quasi-Newton (BFGS, eta = .9) technique was selected as the preferred technique. A series of basic duct design problems are performed. On a two-parameter optimization problem, the BFGS technique is demonstrated to require half as many function evaluations as a steepest descent technique.

Madabhushi, Ravi K.↗

A Wall-Distance Method for Turbulence Modeling

The distance from a grid point to the closest wall surface, wall distance, is a funda- mental quantity in turbulence modeling. Efficiency of wall-distance calculations has become more critical as the size of computational grids has significantly increased in recent years. This paper reports on an initial implementation of a new search-based wall-distance method that is suitable for general unstructured computational fluid dynamics (CFD) grids and tailored for requirements specific for turbulence modeling. The method represents a two-step approach to calculate the wall distance. In the first step, the wall distance is approximated for each grid point as the minimum distance from this point to a vertex of a triangular face at the wall. The point-to-vertex distance calculation is relatively inexpensive but may lead to a significant error in the wall-distance ap- proximation, especially for grid points near the wall. In the second step, for grid points located within a predefined distance ( threshold ) from the wall, the wall distance is computed as the minimum distance to wall faces. As a result, the wall distance is exact for all grid points within the threshold. This two-step approach reduces the computational cost yet achieves high and controllable accuracy in the evaluation of the wall distance. Algorithmic enhancements are presented to improve efficiency of wall-distance computations. Comprehensive assessment of the new method is reported for large-scale unstructured CFD grids generated for the Fifth AIAA CFD High-Lift Prediction Workshop. The performance of the new wall-distance method compares favorably with performance of two established methods implemented in high-performance CFD codes.

Wall Distance↗

Contributions of TetrUSS to Project Orion

The NASA Constellation program has relied heavily on Computational Fluid Dynamics simulations for generating aerodynamic databases and design loads. The Orion Project focuses on the Orion Crew Module and the Orion Launch Abort Vehicle. NASA TetrUSS codes (GridTool/VGRID/USM3D) have been applied in a supporting role to the Crew Exploration Vehicle Aerosciences Project for investigating various aerodynamic sensitivities and supplementing the aerodynamic database. This paper provides an overview of the contributions from the TetrUSS team to the Project Orion Crew Module and Launch Abort Vehicle aerodynamics, along with selected examples to highlight the challenges encountered along the way. A brief description of geometries and tasks will be discussed followed by a description of the flow solution process that produced production level computational solutions. Four tasks conducted by the USM3D team will be discussed to show how USM3D provided aerodynamic data for inclusion in the Orion aero-database, contributed data for the build-up of aerodynamic uncertainties for the aero-database, and provided insight into the flow features about the Crew Module and the Launch Abort Vehicle.

Mcmillin, Susan N.↗

Transonic flow theory of airfoils and wings

There are plans to use the supercritical wing on the next generation of commercial aircraft so as to economize on fuel consumption by reducing drag. Computer codes have served well in meeting the consequent demand for new wing sections. The possibility of replacing wind tunnel tests by computational fluid dynamics is discussed. Another approach to the supercritical wing is through shockless airfoils. A novel boundary value problem in the hodograph plane is studied that enables one to design a shockless airfoil so that its pressure distribution very nearly takes on data that are prescribed.

Garabedian, P. R.↗

Optimal dynamic remapping of parallel computations

A large class of computations are characterized by a sequence of phases, with phase changes occurring unpredictably. The decision problem was considered regarding the remapping of workload to processors in a parallel computation when the utility of remapping and the future behavior of the workload is uncertain, and phases exhibit stable execution requirements during a given phase, but requirements may change radically between phases. For these problems a workload assignment generated for one phase may hinder performance during the next phase. This problem is treated formally for a probabilistic model of computation with at most two phases. The fundamental problem of balancing the expected remapping performance gain against the delay cost was addressed. Stochastic dynamic programming is used to show that the remapping decision policy minimizing the expected running time of the computation has an extremely simple structure. Because the gain may not be predictable, the performance of a heuristic policy that does not require estimnation of the gain is examined. The heuristic method's feasibility is demonstrated by its use on an adaptive fluid dynamics code on a multiprocessor. The results suggest that except in extreme cases, the remapping decision problem is essentially that of dynamically determining whether gain can be achieved by remapping after a phase change. The results also suggest that this heuristic is applicable to computations with more than two phases.

Nicol, David M.↗

Analysis of a hypersonic waverider research vehicle with a hydrocarbon scramjet engine

The results of a feasibility study of a hypersonic waverider research vehicle with a hydrocarbon scramjet engine are presented. The integrated waverider/scramjet geometry is first optimized with a vehicle synthesis code to produce a maximum product of the lift-to-drag ratio and the cycle specific impulse, hence cruise range. Computational fluid dynamics (CFD) is then employed to provide a nose-to-tail analysis of the system at the on-design conditions. Some differences are noted between the results of the two analysis techniques. A comparison of experimental, engineering analysis and CFD results on a waverider forebody are also included for validation.

Molvik, Gregory A.↗

NASA Lewis Research Center Workshop on Forced Response in Turbomachinery

A summary of the NASA Lewis Research Center (LeRC) Workshop on Forced Response in Turbomachinery in August, 1993 is presented. It was sponsored by the following NASA organizations: Structures, Space Propulsion Technology, and Propulsion Systems Divisions of NASA LeRC and the Aeronautics and Advanced Concepts & Technology Offices of NASA Headquarters. In addition, the workshop was held in conjunction with the GUIde (Government/Industry/Universities) Consortium on Forced Response. The workshop was specifically designed to receive suggestions and comments from industry on current research at NASA LeRC in the area of forced vibratory response of turbomachinery blades which includes both computational and experimental approaches. There were eight presentations and a code demonstration. Major areas of research included aeroelastic response, steady and unsteady fluid dynamics, mistuning, and corresponding experimental work.

Stefko, George L.↗

Unsteady Flow in a Supersonic Turbine with Variable Specific Heats

Modern high-work turbines can be compact, transonic, supersonic, counter-rotating, or use a dense drive gas. The vast majority of modern rocket turbine designs fall into these Categories. These turbines usually have large temperature variations across a given stage, and are characterized by large amounts of flow unsteadiness. The flow unsteadiness can have a major impact on the turbine performance and durability. For example, the Space Transportation Main Engine (STME) fuel turbine, a high work, transonic design, was found to have an unsteady inter-row shock which reduced efficiency by 2 points and increased dynamic loading by 24 percent. The Revolutionary Reusable Technology Turbopump (RRTT), which uses full flow oxygen for its drive gas, was found to shed vortices with such energy as to raise serious blade durability concerns. In both cases, the sources of the problems were uncovered (before turbopump testing) with the application of validated, unsteady computational fluid dynamics (CFD) to the designs. In the case of the RRTT and the Alternate Turbopump Development (ATD) turbines, the unsteady CFD codes have been used not just to identify problems, but to guide designs which mitigate problems due to unsteadiness. Using unsteady flow analyses as a part of the design process has led to turbine designs with higher performance (which affects temperature and mass flow rate) and fewer dynamics problems. One of the many assumptions made during the design and analysis of supersonic turbine stages is that the values of the specific heats are constant. In some analyses the value is based on an average of the expected upstream and downstream temperatures. In stages where the temperature can vary by 300 to 500 K, however, the assumption of constant fluid properties may lead to erroneous performance and durability predictions. In this study the suitability of assuming constant specific heats has been investigated by performing three-dimensional unsteady Navier-Stokes simulations for a supersonic turbine stage.

Dorney, Daniel J.↗

Enhanced Control Effector Designs for Airbreathing Transatmospheric Vehicles

A study was conducted to evaluate the potential effectiveness of a moveable cowl-trailing-edge design for airbreathing hypersonic single-stage-to-orbit (SSTO) configurations, which can be extended or deflected from the nominal position in order to provide additional pitching moment capability. This additional pitching moment capability may reduce the necessary deflection angle from conventional control surfaces and the associated trim drag penalty. Calculations for a generic SSTO configuration with baseline and modified cowl trailing edge geometries were performed at Mach 6 and 10 initially with an engineering analysis code in order to examine several design parametrics. In order to more accurately model geometries and flow physics, 2-D viscous computational fluid dynamics (CFD) predictions were obtained. FInally, a limited set of 3D CFD predictions were obtained at Mach 6 in order to show the effects of modeling 3D flow fields as well as the full 3D vehicle geometry. Comparisons of aftbody surface pressures and force and moment predictions show differences between initial predictions and 2-D CFD solutions due to geometry modeling and calculation method differences. The 3-D CFD predictions confirm the trends observed in the 2-D solutions and provide additional information on 3D effects. These analyses shows that cowl-trailing-edge extensions were effective in providing additional (nose-down) pitching moment increments as well as increased thrust compared to the baseline geometry. These effects reduce the control surface deflection angle required to trim the vehicle and the associated trim drag. Cowl-trailingedge deflections were not as effective, generating more nose-up pitching moment and decreased thrust compared to non-deflected cases.

Cockrell, Charles E.↗

Center for Extended Magnetohydrodynamic Modeling (CEMM, UW-Madison Physics)

The Center for Extended Magnetohydrodynamic Modeling (CEMM) focused on the further development, verification, validation, and application of two scientific computation programs, the NIMROD and M3D-C1 codes, for modeling macroscopic dynamics in magnetized plasma. Both NIMROD and M3D-C1 solve the mathematical equations that model two-fluid 3D magnetohydrodynamic (MHD). They are designed and optimized to be efficient and accurate for toroidal confinement physics calculations. The effort improved the fidelity of the physics models, the realism of the treatment of the surrounding vessel, coils, and other structures, the accuracy and efficiency of the algorithms, and the ability to scale calculations to thousands of processors in order to take advantage of leading-edge computers. The primary application areas of Resonant Magnetic Perturbations (RMP), Edge Harmonic Oscillations (EHO), Edge Localized Modes (ELMs), Disruption Studies, Sawteeth and Stationary States, and Kinetic MHD and Neoclassical Tearing Modes (NTM) are of concern for the international ITER experiment, which is under construction in France. The component of CEMM that is described in this report applied recently developed aspects of NIMROD to sawtooth and EHO dynamics, which occur in tokamak plasma confinement systems. In particular, the effort researched energetic-particle, two-fluid, and plasma-flow effects. The work included an effort to validate energetic-particle modeling using laboratory measurements and previous analysis of giant sawtooth oscillations, where energetic particles temporarily stabilize the system. A side effort that grew out of this work was benchmarking of the slab-geometry ion temperature gradient (ITG) instability computed by NIMROD's two-fluid modeling against analytical results. These fluid-model results were also compared with results from kinetic modeling to help quantify the range of validity for the two-fluid model. The EHO simulations show that sheared plasma flow plays a crucial role in saturating the underlying MHD instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Progress in computational unsteady aerodynamics

After vigorous development for over twenty years, Computational Fluid Dynamics (CFD) in the field of aerospace engineering has arrived at a turning point toward maturity. This paper discusses issues related to algorithm development for the Euler/Navier Stokes equations, code validation and recent applications of CFD for unsteady aerodynamics. Algorithm development is a fundamental element for a good CFD program. Code validation tries to bridge the reliability gap between CFD and experiment. Many of the recent applications also take a multidisciplinary approach, which is a future trend for CFD applications. As computers become more affordable, CFD is expected to be a better scientific and engineering tool.

Obayashi, Shigeru↗

Modeling Tools for Propulsion Analysis and Computational Fluid Dynamics on the Internet

The existing RocketWeb(TradeMark) Internet Analysis System (httr)://www.iohnsonrockets.com/rocketweb) provides an integrated set of advanced analysis tools that can be securely accessed over the Internet. Since these tools consist of both batch and interactive analysis codes, the system includes convenient methods for creating input files and evaluating the resulting data. The RocketWeb(TradeMark) system also contains many features that permit data sharing which, when further developed, will facilitate real-time, geographically diverse, collaborative engineering within a designated work group. Adding work group management functionality while simultaneously extending and integrating the system's set of design and analysis tools will create a system providing rigorous, controlled design development, reducing design cycle time and cost.

Muss, J. A.↗