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,153 records · Page 64

Testing New Programming Paradigms with NAS Parallel Benchmarks

Over the past decade, high performance computing has evolved rapidly, not only in hardware architectures but also with increasing complexity of real applications. Technologies have been developing to aim at scaling up to thousands of processors on both distributed and shared memory systems. Development of parallel programs on these computers is always a challenging task. Today, writing parallel programs with message passing (e.g. MPI) is the most popular way of achieving scalability and high performance. However, writing message passing programs is difficult and error prone. Recent years new effort has been made in defining new parallel programming paradigms. The best examples are: HPF (based on data parallelism) and OpenMP (based on shared memory parallelism). Both provide simple and clear extensions to sequential programs, thus greatly simplify the tedious tasks encountered in writing message passing programs. HPF is independent of memory hierarchy, however, due to the immaturity of compiler technology its performance is still questionable. Although use of parallel compiler directives is not new, OpenMP offers a portable solution in the shared-memory domain. Another important development involves the tremendous progress in the internet and its associated technology. Although still in its infancy, Java promisses portability in a heterogeneous environment and offers possibility to "compile once and run anywhere." In light of testing these new technologies, we implemented new parallel versions of the NAS Parallel Benchmarks (NPBs) with HPF and OpenMP directives, and extended the work with Java and Java-threads. The purpose of this study is to examine the effectiveness of alternative programming paradigms. NPBs consist of five kernels and three simulated applications that mimic the computation and data movement of large scale computational fluid dynamics (CFD) applications. We started with the serial version included in NPB2.3. Optimization of memory and cache usage was applied to several benchmarks, noticeably BT and SP, resulting in better sequential performance. In order to overcome the lack of an HPF performance model and guide the development of the HPF codes, we employed an empirical performance model for several primitives found in the benchmarks. We encountered a few limitations of HPF, such as lack of supporting the "REDISTRIBUTION" directive and no easy way to handle irregular computation. The parallelization with OpenMP directives was done at the outer-most loop level to achieve the largest granularity. The performance of six HPF and OpenMP benchmarks is compared with their MPI counterparts for the Class-A problem size in the figure in next page. These results were obtained on an SGI Origin2000 (195MHz) with MIPSpro-f77 compiler 7.2.1 for OpenMP and MPI codes and PGI pghpf-2.4.3 compiler with MPI interface for HPF programs.

Jin, H.↗

Coupled Monte Carlo and thermal-hydraulics modeling of a prismatic gas reactor fuel assembly using Cardinal

Cardinal is a MOOSE application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics to the MOOSE framework, closing the neutronics and thermal-fluid gaps in conducting tightly-coupled, high-resolution multiscale and multiphysics analyses. By leveraging MOOSE's interfaces for wrapping external codes, Cardinal overcomes many challenges encountered in earlier multiphysics coupling works, such as file-based I/O or overly-restrictive geometry mapping requirements. In this work, we leverage a subset of the multiphysics interfaces in Cardinal to perform coupling of OpenMC neutron transport, MOOSE heat conduction, and THM thermal-fluids for steady-state modeling of a prismatic gas reactor fuel assembly.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Thermohydrodynamic Analysis of Cryogenic Liquid Turbulent Flow Fluid Film Bearings

This report describes a thermohydrodynamic analysis and computer programs for the prediction of the static and dynamic force response of fluid film bearings for cryogenic applications. The research performed addressed effectively the most important theoretical and practical issues related to the operation and performance of cryogenic fluid film bearings. Five computer codes have been licensed by the Texas A&M University to NASA centers and contractors and a total of 14 technical papers have been published.

San Andres, Luis↗

Benchmark Wall Heat Flux Data for a GO2/GH2 Single Element Injector

Computational Fluid Dynamics (CFD) is becoming an important component of injector design in the rocket industry. Injector designers who use CFD in the design process need to understand the accuracy level of the particular code being used for certain aspects of the design. This paper presents a recent effort to acquire benchmark quality data to be used for CFD code validation. Detailed chamber wall temperature and heat flux data was acquired for a gaseous oxygen, gaseous hydrogen single element shear coaxial injector in a 1.5 inch diameter copper heat sink chamber at Penn State University. The data was taken using both coaxial and water cooled heat flux gauges. Tests were run using hot gases generated from both fuel and oxidizer preburners. Tests were conducted over a chamber pressure range of 300 to 750 psia. Data analysis and uncertainty information will also be presented.

Pal, Sibtosh↗

Equivalent Source Method Applied to Launch Acoustic Simulations

Aeroacoustic simulations of the launch environment are described. A hybrid computational fluid dynamics (CFD)/computational aeroacoustic (CAA) approach is developed in order to accurately and efficiently predict the sound pressure level spectrum on the launch vehicle and surrounding structures. The high-fidelity CFD code LAVA (Launch Ascent and Vehicle Analysis), is used to generate pressure time history at select locations in the flow field. A 3D exterior Helmholtz solver is then used to iteratively determine a set of monopole sources which mimic the noise generating mechanisms identified by the CFD solver. The acoustic pressure field generated from the Helmholtz solver is then used to evaluate the sound pressure levels.

Housman, Jeffrey A.↗

Parallel High-Order Anisotropic Meshing Using Discrete Metric Tensors

This paper presents a metric-aligned meshing algorithm that relies on the Lp-Centroidal Voronoi Tesselation approach. A prototype of this algorithm was first presented at the Scitech conference of 2018 and this work is an extension to that paper. At the end of the previously presented work, a set of problems were mentioned which we are trying to address in this paper. First, we show a significant improvement in code performance since we were limited to present relatively benign (analytical) test cases. Second, we demonstrate here that we are able to rely on discrete metric data that is delivered by a Computational Fluid Dynamics (CFD) solver. Third, we demonstrate how to generate high-order curved elements that are aligned with the underlying discrete metric field.

Ekelschot, Dirk↗

Fluid modeling of low-temperature plasmas

Fluid models are essential for understanding and predicting low-temperature plasma (LTP) behavior in various scientific and industrial settings. This paper provides an introductory tutorial on fluid modeling of LTPs, covering model formulation, implementation, and computational simulations. The tutorial focuses on five main components of the formulation of LTP fluid models: fluid flow, energy, chemistry, electromagnetism, and material properties, as well as in essential aspects of model implementations, including multiscale phenomena, multiphysics coupling, and numerical convergence. Designed for students and early-career researchers, this work offers a practical foundation for developing and using fluid models, from in-house computational codes to commercial software, bridging fundamental theory with real-world applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermal Stratification Modeling for Sodium-Cooled Fast Reactors: A Status Update

Thermal hydraulic behavior in the upper plenum of pool- type sodium-cooled fast reactors (SFRs) is a major concern, as many design challenges are concentrated in this region. As SFR designs aim for licensing and commercialization, it is important to accurately analyze and predict the thermal-hydraulic behavior in this region during accident scenarios, specifically thermal stratification. Thermal stratification models are currently a major source of uncertainty in most system codes for all types of power plants. Most system codes, including SAS4A/SASSYS-1, a system level code developed by Argonne National Laboratory (Argonne), use very coarse meshes that cannot capture the complexities of the stratification phenomena. While the commonly employed lumped-volume based models for thermal stratification are able to run in a matter of seconds, they result in approximate results and can only handle simple cases. Other 2-D and 3-D methods, such as computational fluid dynamics (CFD) models, can analyze simple configurations with higher fidelity, but come with a relatively large computational expense. Finding a modeling solution that is both accurate and computationally efficient has proven difficult. This paper provides details of a review and gap analysis of the various modeling approaches proposed to date and explores a path forward for future thermal stratification modeling efforts, with a focus on developing new models for the SAS4A/SASSYS-1 system code.

Morgan, Sarah↗

Progress in the validation of rotor aerodynamic codes using field data

Within the framework of the fourth phase of the International Energy Agency (IEA) Wind Task 29, a large comparison exercise between measurements and aeroelastic simulations has been carried out featuring three simulation cases in axial, sheared and yawed inflow conditions. Results were obtained from more than 19 simulation tools originating from 12 institutes, ranging in fidelity from blade element momentum (BEM) to computational fluid dynamics (CFDs) and compared to state-of-the-art field measurements from the 2 MW DanAero turbine. More than 15 different variable types ranging from lifting-line variables to blade surface pressures, loads and velocities have been compared for the different conditions, resulting in over 250 comparison plots. The result is a unique insight into the current status and accuracy of rotor aerodynamic modeling. For axial flow conditions, a good agreement was found between the various code types, where a dedicated grid sensitivity study was necessary for the CFD simulations. However, compared to wind tunnel experiments on rotors featuring controlled conditions, it remains a challenge to achieve good agreement of absolute levels between simulations and measurements in the field. For sheared inflow conditions, uncertainties due to rotational and unsteady effects on airfoil data result in the CFD predictions standing out above the codes that need input of sectional airfoil data. However, it was demonstrated that using CFD-synthesized airfoil data is an effective means to bypass this shortcoming. For yawed flow conditions, it was observed that modeling of the skewed wake effect is still problematic for BEM codes where CFD and free vortex wake codes inherently model the underlying physics correctly. The next step is a comparison in turbulent inflow conditions, which is featured in IEA Wind Task 47. Doing this analysis in cooperation under the auspices of the IEA Wind Technology Collaboration Program (TCP) has led to many mutual benefits for the participants. The large size of the consortium brought ample manpower for the analysis where the learning process by combining several complementary experiences and modeling techniques gave valuable insights that could not be found when the analysis is carried out individually.

17 WIND ENERGY↗

A Study of Cavitation-Ignition Bubble Combustion

We present the results of an experimental and computational study of the physics and chemistry of cavitation-ignition bubble combustion (CIBC), a process that occurs when combustible gaseous mixtures are ignited by the high temperatures found inside a rapidly collapsing bubble. The CIBC process was modeled using a time-dependent compressible fluid-dynamics code that includes finite-rate chemistry. The model predicts that gas-phase reactions within the bubble produce CO and other gaseous by-products of combustion. In addition, heat and mechanical energy release through a bubble volume-expansion phase are also predicted by the model. We experimentally demonstrate the CIBC process using an ultrasonically excited cavitation flow reactor with various hydrocarbon-air mixtures in liquid water. Low concentrations (< 160 ppm) of carbon monoxide (CO) emissions from the ultrasonic reactor were measured, and found to be proportional to the acoustic excitation power. The results of the model were consistent with the measured experimental results. Based on the experimental findings, the computational model, and previous reports of the "micro-diesel effect" in industrial hydraulic systems, we conclude that CIBC is indeed possible and exists in ultrasonically- and hydrodynamically-induced cavitation. Finally, estimates of the utility of CIBC process as a means of powering an idealized heat engine are also presented.

Nguyen, Quang-Viet↗

Bichromatic Wave Selection for Validation of the Difference-Frequency Transfer Function for the OC6 Validation Campaign: Preprint

The focus of the Offshore Code Comparison Collaboration, Continuation, with Correlation and unCertainity (OC6) project, which operates under the International Energy Agency (IEA) Wind Task 30, is to refine the accuracy of engineering tools used to design offshore wind turbines. In support of this work, a new validation campaign is being developed that seeks to better understand the nonlinear wave loading that excites floating wind systems at their low-frequency, rigid-body modes in surge and pitch. The validation data will be used in a three-way validation between simplified engineering tools and higher-fidelity tools, such as computational fluid dynamics (CFD). Irregular wave spectrums, which are traditionally used to examine the nonlinear wave interaction with offshore structures, are too computationally expensive to be simulated in CFD tools, and so bichromatic wave cases will be used instead. This paper reviews the process used to choose the bichromatic wave pairs to be used in the campaign to validate the second-order difference-frequency quadratic and potential loads at the surge and pitch natural frequencies of a floating semisubmersible.

17 WIND ENERGY↗

Improvements to the Blade Element Momentum Formulation of OpenFAST for Skewed Inflows

In this work, we modify the blade element momentum algorithm of OpenFAST to improve its predictions under large skewed inflow conditions. We use the well-known Glauert's skew correction and introduce continuous extension of the model for high-thrust conditions. We present the rationale behind Glauert's empirical model and discuss the different conventions possible for the axial induction factor. We verify the model against the higher-fidelity lifting-line vortex method and blade-resolved computational fluid dynamics, and we observe that the new implementation enhances the accuracy and reliability of OpenFAST's aerodynamic modeling capabilities in conditions involving large skew angles. For the parametric studies run using the different codes, we find that the power changes with the skew angle as cos 1.7 (θ skew ) and the thrust as cos 0.65 (θ skew ). An analysis of the azimuthal variation of the induced velocities in the rotor plane reveals that current redistribution models used in blade element momentum codes may need to be refined.

17 WIND ENERGY↗

Computational Fluid Dynamics Study of a Cross-Flow Marine Hydrokinetic Turbine and the Combined Influence of Struts and Helical Blades

A computational fluid dynamics study was performed for a cross-flow marine hydrokinetic turbine. The analysis was done in three dimensions and used the unsteady Reynolds-averaged Navier-Stokes solver in the commercial code STAR-CCM+. The base turbine configuration is the RivGen(R) Turbine, designed by the Ocean Renewable Power Company. A convergence and uncertainty analysis was performed for both the spatial and temporal discretization; this was done using the base configuration, which features support struts and helical foils. Both struts and helical blades introduce three-dimensional flow effects, influencing the complex flow phenomenon of dynamic stall. The study compares the relative impact of struts on power performance and blade loading for both helical and straight blades, and found that for this turbine the relative loss in power from struts was lower with helical blades.

CFD↗

Evolution and status of CFD techniques for scramjet applications

Computational Fluid Dynamics (CFD) methods are presently applied to the design and analysis of scramjet propulsion systems for missiles and aircraft operating at speeds in excess of Mach 5, in forms ranging from simple inviscid codes for parametric studies of inlet geometries to Navier-Stokes codes incorporating chemistry models which address complex scramjet combustion phenomena. Attention is given to the development history of CFD for the hypersonic regime, as well as to scramjet flowfield modeling in inlets, combustors and nozzles, and prospective development trends.

White, M. E.↗

Investigation of the transient fuel preburner manifold and combustor

A computational fluid dynamics (CFD) model with finite rate reactions, FDNS, was developed to study the start transient of the Space Shuttle Main Engine (SSME) fuel preburner (FPB). FDNS is a time accurate, pressure based CFD code. An upwind scheme was employed for spatial discretization. The upwind scheme was based on second and fourth order central differencing with adaptive artificial dissipation. A state of the art two-equation k-epsilon (T) turbulence model was employed for the turbulence calculation. A Pade' Rational Solution (PARASOL) chemistry algorithm was coupled with the point implicit procedure. FDNS was benchmarked with three well documented experiments: a confined swirling coaxial jet, a non-reactive ramjet dump combustor, and a reactive ramjet dump combustor. Excellent comparisons were obtained for the benchmark cases. The code was then used to study the start transient of an axisymmetric SSME fuel preburner. Predicted transient operation of the preburner agrees well with experiment. Furthermore, it was also found that an appreciable amount of unburned oxygen entered the turbine stages.

Wang, Ten-See↗

SmaggIce User Guide

SmaggIce (Surface Modeling And Grid Generation for Iced Airfoils) is one of NASNs aircraft icing research codes developed at the Glenn Research Center. It is a software toolkit used in the process of aerodynamic performance prediction of iced airfoils. It includes tools which complement the 2D grid-based Computational Fluid Dynamics (CFD) process: geometry probing; surface preparation for gridding: smoothing and re-discretization of geometry. Future releases will also include support for all aspects of gridding: domain decomposition; perimeter discretization; grid generation and modification.

Baez, Marivell↗

A Perspective on Computational Aerothermodynamics at NASA

The evolving role of computational aerothermodynamics (CA) within NASA over the past 20 years is reviewed. The presentation highlights contributions to understanding the Space Shuttle pitching moment anomaly observed in the first shuttle flight, prediction of a static instability for Mars Pathfinder, and the use of CA for damage assessment in post-Columbia mission support. In the view forward, several current challenges in computational fluid dynamics and aerothermodynamics for hypersonic vehicle applications are discussed. Example simulations are presented to illustrate capabilities and limitations. Opportunities to advance the state-of-art in algorithms, grid generation and adaptation, and code validation are identified.

Gnoffo, Peter A.↗

right-sized dimple evaluator

The code (aka computer program written as a Matlab script) uses a unique set of n independent equations to solve for n turbulence variables. The code requires the input of a characteristic dimension, a characteristic fluid velocity, the fluid dynamic viscosity, and the fluid density. Most importantly, the code estimates the size of three key turbulent eddies: Kolmogorov, Taylor, and integral. Based on the eddy sizes, dimples dimensions are prescribed such that the key eddies (principally Taylor, and sometimes Kolmogorov), can be generated by the dimple rim and flow unimpeded through the dimple’s concave cavity. It is hypothesized that turbulent eddies are generated by the dimple rim at the dimple-surface interface. The newly-generated eddies in turn entrain the movement of surrounding regions of fluid, creating more mixing. The eddies also generate lift near the wall surrounding the dimple, as they accelerate and reduce pressure in the regions near and at the dimple cavity, thereby minimizing the fluid drag.

Rodriguez, Sal↗