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 937 records · Page 52

Validation of A Multiphase Computational Fluid Dynamics Model for Vapor Pull-Through in Normal and Low Gravity

On-orbit fluid transfer such as refueling of propellant tanks and life-support systems can enable long-duration space missions. For safe and efficient liquid transfer operations, prior knowledge of liquid positioning and liquid-vapor interface behavior while draining in a reduced-gravity environment is required. Numerical models capable of predicting vapor ingestion (or vapor pull-through) can be used to design liquid transfer operations while reducing liquid residuals, mission risk and settling thrust required to prevent vapor ingestion. An experimental program conducted in the 2.2 Second Drop Tower facility at NASA Lewis Research Center in 1969 investigated the vapor ingestion phenomenon for a range of outflow rates and tank sizes providing a database for validation. This study presents a Computational Fluid Dynamics model capable of accurately predicting the vapor ingestion using the Volume-of-Fluid multiphase solver in commercial code STAR-CCM+. A description of the experimental setup and general trends from similar studies are presented. Comparisons of the numerical prediction and test data in normal and low gravity show good agreement and give confidence in pursuing design of full-scale propellant transfer systems.

Propellant Transfer↗

Reduced-Order Aerodynamic Modeling Based on CFD Frequency Responses from Multisine Inputs

A system identification analysis was performed to determine reduced-order models of a computational fluid dynamics (CFD) solver for linear aeroelastic analysis and control design. The application was to the FUN3D code and the flexible half-span wind tunnel test article, in transonic flow conditions, used in the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project. Multiple inputs (structural mode displacements and control surface deflections) were simultaneously excited with orthogonal phase-optimized multisines and multiple outputs (generalized aerodynamic forces) were recorded, from which the matrix of frequency responses were computed using a single CFD run. A state-space model was then fit to the frequency response data using a maximum-likelihood estimator.

System identification↗

Energy Exascale Computational Fluid Dynamics Simulations With the Spectral Element Method

Development and application of the open-source GPU-based fluid-thermal simulation code, NekRS, are described. Time advancement is based on an efficient kth-order accurate timesplit formulation coupled with scalable iterative solvers. Spatial discretization is based on the high-order spectral element method (SEM), which affords the use of fast, low-memory, matrix-free operator evaluation. Further, recent developments include support for nonconforming meshes using overset grids and for GPU-based Lagrangian particle tracking. Results of large-eddy simulations of atmospheric boundary layers for wind-energy applications as well as extensive nuclear energy applications are presented.

42 ENGINEERING↗

Computational fluid dynamics - The coming revolution

The development of aerodynamic theory is traced from the days of Aristotle to the present, with the next stage in computational fluid dynamics dependent on superspeed computers for flow calculations. Additional attention is given to the history of numerical methods inherent in writing computer codes applicable to viscous and inviscid analyses for complex configurations. The advent of the superconducting Josephson junction is noted to place configurational demands on computer design to avoid limitations imposed by the speed of light, and a Japanese projection of a computer capable of several hundred billion operations/sec is mentioned. The NASA Numerical Aerodynamic Simulator is described, showing capabilities of a billion operations/sec with a memory of 240 million words using existing technology. Near-term advances in fluid dynamics are discussed.

Graves, R. A., Jr.↗

Summary of Multiphysics Modeling for Sublimation Mass Transfer

This report summarizes the effort to develop a multiphysics modeling framework for sublimation mass transfer processes. The aim of the project was to develop a multiphysics code capable of predicting phase change between solid and gas in a closed container, and the movement of material within such a container when exposed to various exterior environmental conditions. Two modeling frameworks were developed towards this aim: one with a high-fidelity computational fluid dynamics (CFD) structure, and the other as a fast reduced-order model. Ultimately, both attempts were unsuccessful due to instabilities in the code and physical processes for which there is no adequate numerical representation. Both modeling attempts are briefly detailed before providing a brief survey of recent updates in the literature, and a recommendation for future work on this subject.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Center for Extended Magnetohydrodynamic Modeling (CEMM), Final Technical Report for UW-Madison Engineering 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 effort by the UW-Madison Engineering Physics component of CEMM improved the numerical methods in the NIMROD code with respect to modeling MHD interchange and two-fluid effects. It also improved the algorithm for solving large algebraic systems to enhance NIMROD's computational performance when running the code in parallel on increasing numbers of processor cores. The modeling efforts included a study of the nonlinear evolution of ballooning, and the findings were supported by analytical theory and computation. Numerical MHD modeling also examined the influences of two-fluid effects on ELMs and RMPs, where "two-fluid" refers to separation of electron dynamics from ion dynamics. Another related aspect considered the effects of a localized layer of electrical current density near the edge of the plasma confinement region. Two-fluid effects on resonant reconnecting instabilities were verified and applied to computations of sawtooth and drift-tearing dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

risingMicrobubbleLattice

This code is a simulation case to be run with OpenFOAM, an open-source computational fluid dynamics software. A gmsh mesh file is also included. Specifically, this simulation demonstrates the transport of a single microbubble rising through an ordered lattice due to an applied flow. The bubble deforms as it squeezes through the pores of the lattice. Is

Guo, Jack [Lawrence Livermore National Laboratory ↗

A Wave Energy Converter Design Load Case Study

This article presents an example by which design loads for a wave energy converter (WEC) might be estimated through the various stages of the WEC design process. Unlike previous studies, this study considers structural loads, for which, an accurate assessment is crucial to the optimization and survival of a WEC. Three levels of computational fidelity are considered. The first set of design load approximations are made using a potential flow frequency-domain boundary-element method with generalized body modes. The second set of design load approximations are made using a modified version of the linear-based time-domain code WEC-Sim. The final set of design load simulations are realized using computational fluid dynamics coupled with finite element analysis to evaluate the WEC’s loads in response to both regular and focused waves. This study demonstrates an efficient framework for evaluating loads through each of the design stages. In comparison with experimental and high-fidelity simulation results, the linear-based methods can roughly approximate the design loads and the sea states at which they occur. The high-fidelity simulations for regular wave responses correspond well with experimental data and appear to provide reliable design load data. The high-fidelity simulations of focused waves, however, result in highly nonlinear interactions that are not predicted by the linear-based most-likely extreme response design load method.

16 TIDAL AND WAVE POWER↗

NASA's CFD Validation Program

With computational fluid dynamics (CFD) becoming a productive research and design tool, the requirement to validate CFD codes has grown significantly. NASA had emphasized CFD validation activities since 1986 when a separate work element was formed to fund experimental activities related to validation. NASA's CFD and CFD validation programs are closely coordinated to ensure that experimental data bases are available as soon as possible for validating codes. In response to industry and academic requirements, four levels of experimental research have been defined as part of CFD validation with NASA's Aeronautics Advisory Committee (AAC) support although only the fourth level actually has the detailed information necessary for validating codes. Critical flow physics especially turbulence modeling are key to improved CFD codes. NASA has focused additional resources on transition and turbulence physics to meet these requirements. With improved turbulence models, CFD codes will be more accurate, robust, and efficient. However, with the level of detailed information available from CFD codes, highly accurate and detailed experiments are required to capture the critical information for validating codes. Advanced instrumentation especially non-intrusive instrumentation is required to acquire this information in validation experiments. The CFD validation program is being coordinated and managed to address these critical activities. A list of experiments which are currently being supported at least partially are included.

Satran, Dale R.↗

Determination of Rotordynamic Coefficients for Labyrinth Seals and Application to Rotordynamic Design Calculations

In today's rotordynamic calculations, the input parameters for a finite element analysis (FEA) determine very much the reliability of eigenvalue and eigenmode predictions. While modeling of an elastic structure by means of beam elements etc. is relatively straightforward to perform and the input data for journal bearings are usually known exactly enough, the determination of stiffness and damping for labyrinth seals is still the subject of many investigations. Therefore, the rotordynamic influence of labyrinths is often not included in FEA for rotating machinery because of a lack of computer programs to calculate these parameters. This circumstance can give rise to severe vibration problems especially for high performance turbines or compressors, resulting in remarkable economic losses. The forces generated in labyrinths can be described for small motions around the seal center with a linearized force-motion relationship. Several years ago, we started with the development of computer codes for the determination of rotordynamic seal coefficients. Our different approaches to evaluate the dynamic fluid forces generated by turbulent, compressible seal flow are introduced.

Weiser, P.↗

Turbine Technology Team - An overview of current and planned activities relevant to the National Launch System (NLS)

The current status of the activities and future plans of the Turbine Technology Team of the Consortium for Computational Fluid Dynamics is reviewed. The activities of the Turbine Team focus on developing and enhancing codes and models, obtaining data for code validation and general understanding of flows through turbines, and developing and analyzing the aerodynamic designs of turbines suitable for use in the Space Transportation Main Engine fuel and oxidizer turbopumps. Future work will include the experimental evaluation of the oxidizer turbine configuration, the development, analysis, and experimental verification of concepts to control secondary and tip losses, and the aerodynamic design, analysis, and experimental evaluation of turbine volutes.

Griffin, Lisa W.↗

Application of a distributed network in computational fluid dynamic simulations

A general-purpose 3-D, incompressible Navier-Stokes algorithm is implemented on a network of concurrently operating workstations using parallel virtual machine (PVM) and compared with its performance on a CRAY Y-MP and on an Intel iPSC/860. The problem is relatively computationally intensive, and has a communication structure based primarily on nearest-neighbor communication, making it ideally suited to message passing. Such problems are frequently encountered in computational fluid dynamics (CDF), and their solution is increasingly in demand. The communication structure is explicitly coded in the implementation to fully exploit the regularity in message passing in order to produce a near-optimal solution. Results are presented for various grid sizes using up to eight processors.

Deshpande, Manish↗

Helicopter Blade-Vortex Interaction Noise with Comparisons to CFD Calculations

A comparison of experimental acoustics data and computational predictions was performed for a helicopter rotor blade interacting with a parallel vortex. The experiment was designed to examine the aerodynamics and acoustics of parallel Blade-Vortex Interaction (BVI) and was performed in the Ames Research Center (ARC) 80- by 120-Foot Subsonic Wind Tunnel. An independently generated vortex interacted with a small-scale, nonlifting helicopter rotor at the 180 deg azimuth angle to create the interaction in a controlled environment. Computational Fluid Dynamics (CFD) was used to calculate near-field pressure time histories. The CFD code, called Transonic Unsteady Rotor Navier-Stokes (TURNS), was used to make comparisons with the acoustic pressure measurement at two microphone locations and several test conditions. The test conditions examined included hover tip Mach numbers of 0.6 and 0.7, advance ratio of 0.2, positive and negative vortex rotation, and the vortex passing above and below the rotor blade by 0.25 rotor chords. The results show that the CFD qualitatively predicts the acoustic characteristics very well, but quantitatively overpredicts the peak-to-peak sound pressure level by 15 percent in most cases. There also exists a discrepancy in the phasing (about 4 deg) of the BVI event in some cases. Additional calculations were performed to examine the effects of vortex strength, thickness, time accuracy, and directionality. This study validates the TURNS code for prediction of near-field acoustic pressures of controlled parallel BVI.

McCluer, Megan S.↗

Computational Fluid Dynamics Technology for Hypersonic Applications

Several current challenges in computational fluid dynamics and aerothermodynamics for hypersonic vehicle applications are discussed. Example simulations are presented from code validation and code benchmarking efforts to illustrate capabilities and limitations. Opportunities to advance the state-of-art in algorithms, grid generation and adaptation, and code validation are identified. Highlights of diverse efforts to address these challenges are then discussed. One such effort to re-engineer and synthesize the existing analysis capability in LAURA, VULCAN, and FUN3D will provide context for these discussions. The critical (and evolving) role of agile software engineering practice in the capability enhancement process is also noted.

Gnoffo, Peter A.↗

Thermal Stratification Analysis for Sodium Fast Reactors

The sodium fast reactor (SFR) is the most mature reactor concept of all the generation-IV nuclear systems and is a promising reactor design that is currently under development by several organizations. The majority of sodium fast reactor designs utilize a pool type arrangement which incorporates the primary coolant pumps and intermediate heat exchangers within the sodium pool. These components typically protrude into the pool thus reducing the risk and severity of a loss of coolant accidents. To further ensure safe operation under even the most severe transients a more comprehensive understanding of key thermal hydraulic phenomena in this pool is desired. One of the key technology gaps identified for SFR safety is determining the extent and the effects of thermal stratification developing in the pool during postulated accident scenarios such as a protected or unprotected loss of flow incident. In an effort to address these issues, detailed flow models of transient stratification in the pool during an accident can be developed. However, to develop the calculation models, and ensure they can reproduce the underlying physics, highly spatially resolved data is needed. This data can be used in conjunction with advanced computational fluid dynamic calculations to aid in the development of simple reduced dimensional models for systems codes such as SAM and SAS4A/SASSYS-1.

Schneider, James↗

Thermal Stratification Analysis for Sodium Fast Reactors

The sodium fast reactor (SFR) is the most mature reactor concept of all the generation-IV nuclear systems and is a promising reactor design that is currently under development by several organizations. The majority of sodium fast reactor designs utilize a pool type arrangement which incorporates the primary coolant pumps and intermediate heat exchangers within the sodium pool. These components typically protrude into the pool thus reducing the risk and severity of a loss of coolant accidents. To further ensure safe operation under even the most severe transients a more comprehensive understanding of key thermal hydraulic phenomena in this pool is desired. One of the key technology gaps identified for SFR safety is determining the extent and the effects of thermal stratification developing in the pool during postulated accident scenarios such as a protected or unprotected loss of flow incident. In an effort to address these issues, detailed flow models of transient stratification in the pool during an accident can be developed. However, to develop the calculation models, and ensure they can reproduce the underlying physics, highly spatially resolved data is needed. This data can be used in conjunction with advanced computational fluid dynamic calculations to aid in the development of simple reduced dimensional models for systems codes such as SAM and SAS4A/SASSYS-1.

Schneider, James↗

Development and Validation of SAM Multi-dimensional Flow Model for Thermal Mixing and Stratification Modeling

Thermal mixing and stratification in large pools or enclosures are very important phenomena that are critical to nuclear reactor safety. Because of the wide ranges of time and length scales associated with such phenomena, accurate modeling and simulation of thermal mixing and stratification remain as the key unresolved, challenging problems for reactor transient analyses. In traditional system analysis codes, simplified zero-dimensional (0-D) models are widely used for their high numerical efficiency, but they generally suffer from very limited prediction accuracies or range of applicability. Like traditional system analysis codes, the current version of the SAM code has implemented such simplified 0-D and one-dimensional mixing models. On the other hand, high-resolution Computational Fluid Dynamics (CFD) tools are often used to model complex thermal mixing and stratification phenomena. They are, however, generally numerically expensive, and they require large amounts of computational resources. It is therefore desirable to implement advanced and efficient thermal mixing and stratification modeling capabilities embedded in a system analysis code. This approach will improve the accuracy of reactor safety analyses when thermal mixing and stratification are involved, and also avoid using the large computational resources needed for high-resolution CFD analysis. Currently, with the support of the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation program, an effort has been launched to develop and implement a multi-dimensional flow model in the system analysis code SAM, and demonstrate its applications to model thermal mixing and stratification phenomena in large enclosures. The main outcomes of this research and development activity are summarized in this report, which presents an attempt to include a built-in advanced multi-dimensional flow model in a system analysis code with the focus on overcoming the simulation challenges of thermal mixing and stratification phenomena. In this report, we start with the introduction of existing SAM code capabilities to simulate thermal mixing and stratification phenomena, which is followed by a short summary of the multi-dimensional model implemented in the SAM code, including both the physical model and the Finite Element Method code implementation. In this study, two options were implemented in the SAM code to model turbulent flows: a relatively simple built-in turbulence model and an interface to accept externally computed turbulent viscosities (e.g., from a high- resolution CFD simulation). Code validation studies on this newly added capability were then carried out to compare SAM simulation results with experimental data from the SUPERCAVNA facility, which was designed to study the complex flow recirculation and thermal stratification phenomena relevant to sodium fast reactor designs. In this study, one transient and two steady-state test cases were used for code validation. Different approaches have been used to model the complex turbulence flow fields in the SUPERCAVNA facility. A highly simplified zero-equation turbulence model was first used, but it was determined that it is too simple to capture the complex turbulence flow fields in these test cases. Subsequently, the code validation continued with the use of turbulent viscosity data from high-resolution STAR-CCM+ CFD simulations to improve the accuracy of the results. Using this approach, the SAM simulation results showed very good agreement with both the SUPERCAVNA experimental data and STAR-CCM+ simulation results. In this report, we demonstrate the development, implementation, and successful validation of a multi-dimensional flow model in the SAM code, which aims to improve the simulation accuracy for complex thermal mixing and stratification phenomena. Lessons have also been learned, including that in cases where the flow fields are not well predicted by the zero-equation model, the iinput of turbulent viscosities from an external source can enhance the overall predictive capabilities needed to accurately capture complex thermal-hydraulic phenomena. Therefore, future research will be needed to further improve the code’s capabilities, such as by developing a more efficient and robust approach to capture the turbulence effects in the SAM code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

PSE-Based Aerodynamic Flow Transition Prediction Using Automated Unstructured CFD Integration

Accurate, robust, and efficient prediction of transition in viscous flows is a significant challenge in computational fluid dynamics. We present a coupled, high-fidelity, iterative approach that leverages the FUN3D flow solver and the LASTRAC stability code to predict transition in low-disturbance environments, initiated by the linear growth of boundary-layer instability modes. Our method integrates the ability of FUN3D to compute mixed laminar-transitional-turbulent mean flows via transition-sensitized Reynolds-Averaged Navier-Stokes equations with the ability of LASTRAC to perform linear stability analysis, all within an automated framework that requires no intermediate user involvement. Unlike conventional frameworks that rely on classical stability theory or reduced-order metamodels, our approach employs the parabolized stability equations to provide more accurate and reliable estimates of disturbance growth for multiple instability mechanisms, including Tollmien-Schlichting, Kelvin-Helmholtz, and crossflow modes. By accounting for the effects of mean-flow nonparallelism as well as the surface curvature, this approach lays the foundation for improved N-factor correlations for transition onset prediction in a broad class of flows. We apply this method to three distinct flow configurations: 1) flow over a zero-pressure-gradient flat plate, 2) the NLF-0416 airfoil with both natural and separation-induced transitions, and 3) a 6:1 prolate spheroid, where transition is primarily driven by crossflow instability. For the two-dimensional cases, a formulated intermittency distribution is used to model the transition zone in between the laminar and fully turbulent flows. The results include comparisons with experimental measurements, similar numerical approaches, and transport-equations-based models, demonstrating good agreement in surface pressure coefficients, transition onset locations, and skin-friction coefficients for all three configurations. Besides contributing a couple of new insights into boundary-layer transition in these canonical cases, this study provides a powerful tool for transition modeling in both research and design applications in aerodynamics.

Transition Prediction↗