Results for “computational fluid dynamics simulation” — Page 5

DOE OSTI2022

Impact of large-scale effects on mass transfer and concentration polarization in Reverse Osmosis membrane systems

Here, we present well-resolved computational fluid dynamics simulations of a large-scale reverse osmosis membrane-spacer configuration (1 m). Our computational model solves the flow and transport equations with variable solute-dependent properties. We utilize a high resolution computational mesh to resolve all relevant length scales associated with spacer-induced mixing and thin concentration boundary layers. An important contribution of this work is the development of a modified mass-transfer correlation that accounts for the development of the concentration boundary layer along the channel. A set of 2D axisymmetric simulations were performed for a spiral wound module layer with varying cross-flow conditions and spacer diameters which indicate a significant entrance length effect for concentration profile development at lower flow rates while mixing effects dominate at higher flow rates. The mass-transfer correlations at higher flow rates compare well with published correlations while a surrogate model for Sherwood number was obtained that depends on an additional similarity variable that accounted for entrance length effects at lower flow rates. Finally, a large-scale membrane-spacer design relevant to high-pressure reverse osmosis is studied with a non-uniform arrangement of spacers, which indicate a substantial saving in pressure drop (~40%) compared to traditional uniformly spaced pattern with minor variations (~2%) in concentration polarization, product water quality (~1%) and water recovery (~7%) compared to a uniform spacer pattern.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Comparison of a Full-Scale and a 1:10 Scale Low-Speed Two-Stroke Marine Engine Using Computational Fluid Dynamics

International marine shipping is a growing component of international trade; a vast majority of all the world’s goods are being transported on large ocean-going vessels. The International Maritime Organization (IMO) introduced the Energy Efficiency Design Index in 2013, a regulatory framework of associated metrics for reducing emissions of CO 2 per tonne-mile from shipping by approximately 10% each decade. Therefore, decarbonizing the maritime sector requires the development of new fuel sources. Because of the extremely large physical size of the internal combustion engines present in shipping vessels, experimental iterative development of the engine and fuel system is cost-prohibitive. Thus, the ability to perform combustion system development in a scaled platform that can be more easily operated and modeled computationally is of interest. To that end, scaling relationships are needed to translate the results from a smaller engine to a larger counterpart. Scaling studies to date have been restricted to low scaling ratios, four-stroke light-duty engines, and under-resolved computational fluid dynamic simulations that likely do not accurately capture the physics of scaling. In this work, computational models of a 1:10 scale and a full-scale two-stroke crosshead low-speed marine engine were created and validated against experiments obtained in a real 1:10 scale engine installed at Oak Ridge National Laboratory. Further, due to the large size of the full-scale engine, the model required large high-performance computing resources to be evaluated. The availability of high-performance computing resources at the Department of Energy’s Leadership Computing Facilities is an enabler of the current work. The results of the small- and large-scale engine simulations were compared to analyze the effectiveness of the appropriate scaling laws under these extreme scaling ratio conditions.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRSDate not supplied

Observations of Keyhole Porosity and Comparisons to Analytical Models for Ti-6Al-4V Powder Bed Fusion

Keyhole porosity defects are a common concern in powder bed fusion (PBF) processing. Keyhole porosity prediction models have generally fallen into two categories – high fidelity computational fluid dynamics simulations and low fidelity analytical or empirical models based on keyhole vibration dynamics. The low computational cost of low fidelity models allows them to better approach part-scale predictions. However, studies on low fidelity modeling techniques are limited by the lack of experimental data to assess the validity of the calibration over a wide range of processing conditions. This work extracts and quantifies keyhole porosity across 14 laser velocities, two laser powers, and six+ repetitions for a total of 176 independent trials. The measured porosity data are compared to low fidelity keyhole models to assess their success in predictive porosity occurrence. This work impacts the field by providing independent validation of keyhole porosity models for PBF for use in part-scale defect prediction.

powder bed fusion↗

NASA NTRS2016

Extension of Miles Equation for Ring Baffle Damping Predictions to Small Slosh Amplitudes and Large Baffle Widths

The Miles equation has long been used to predict slosh damping in liquid propellant tanks due to ring baffles. The original work by Miles identifies defined limits to its range of application. Recent evaluations of the Space Launch System identified that the Core Stage baffle designs resulted in violating the limits of the application of the Miles equation. This paper describes the work conducted by NASA/MSFC to develop methods to predict slosh damping from ring baffles for conditions for which Miles equation is not applicable. For asymptotically small slosh amplitudes or conversely large baffle widths, an asymptotic expression for slosh damping was developed and calibrated using historical experimental sub-scale slosh damping data. For the parameter space that lies between region of applicability of the asymptotic expression and the Miles equation, Computational Fluid Dynamics simulations of slosh damping were used to develop an expression for slosh damping. The combined multi-regime slosh prediction methodology is shown to be smooth at regime boundaries and consistent with both sub-scale experimental slosh damping data and the results of validated Computational Fluid Dynamics predictions of slosh damping due to ring baffles.

West, Jeff↗

DOE OSTI2025

Thermal-Fluid and Thermal-Structural Response of the T-Tube Modular Divertor to Spatiotemporally Varying Heat Loads

Tungsten (W) is the leading candidate for divertor target plates because of its high melting point (>3000°C), thermal conductivity, and ultimate tensile stress. While W and its alloys are the only solid materials that can survive the high heat fluxes incident on the divertor, W’s low-ductility high ductile-to-brittle transition temperature of ~600°C and relatively low recrystallization temperature (RT) of ~1300°C pose structural (among other) challenges. The objective of this work is to estimate the thermal-fluid and thermal-structural performance of the helium (He)-cooled T-tube divertor, which was originally developed by the Advanced Reactor Innovation and Evaluation Study (ARIES) using numerical simulations. Here, predictions of temperature distributions across the plasma-facing structural component and surface pressures from computational fluid dynamics simulations are used to determine stress distributions using commercial structural finite element modeling software over a range of fusion-relevant conditions. The maximum allowable incident heat fluxes are determined based on the temperature limits imposed by the ITER elastic Structural Design Criteria for In-vessel Components (SDC-IC) and the maximum RT over a range of He mass flow rates and presented in the form of performance design charts. Our recent work found that thermal- structural criteria accounting for the low ductility of W in a finger-type modular divertor constrain the maximum incident heat fluxes to values well below the ITER specifications, and those based on considering only the RT demonstrate that integrated thermal-fluid and elastic structural performance evaluation are required for accurate assessment of divertor performance. This novel analysis of the T-tube considers how nonuniform and transient incident heat fluxes affect its thermal-fluid and thermal-structural performance, as well as the effect of volumetric heating, which can be as great as 27% of the power incident on the divertor surface. The W tile of the T-tube, with its relatively large plasma-facing area of ~15 cm 2 , will likely experience significant spatial variations in incident heat flux. This work therefore assesses whether steady-state incident heat flux profiles with a peak of 10 MW/m 2 and maximum heat flux gradients of 200 MW/m 2 per m exceed the structural limits imposed by the ITER elastic SDC-IC and the maximum RT over a range of fusion-relevant conditions. The effect of transient heat fluxes typical of plasma detachment and reattachment from the target plate due, for example, to gas injection are also evaluated

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRSDate not supplied

On Characterization of Flow Disturbances in Arc-Jet Testing

This paper reports computational simulations and analysis of flow characterization tests in high enthalpy arc-jet facilities at NASA Ames Research Center. Flow disturbances and their characterization are presented through case studies, and their implications for testing are discussed. The case studies include free-jet test configurations from three different arc heaters and nozzles: the 60-MW IHF 30-inch conical nozzle, the 10-MW TP3 15-inch conical nozzle, and the 20-MW AHF 12-inch conical nozzle. For all cases, test articles are placed in the jet exiting the conical nozzle, and the existence of a flow disturbance is confirmed through flow survey data using pitot pressure and heat flux probes and accompanying analysis. The paper focuses on cases where the effects of the disturbances on test article surface quantities are important. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzles and test box, including the models tested. Comparisons of computations with the experimental measurements are presented. The computations that reproduce the probe sweep data approximately are essential to interpret the arc-jet test data accurately, while providing insights into several observed calorimeter anomalies caused by flow disturbances.

Arc-jet testing↗

NASA NTRSDate not supplied

On Characterization of Flow Disturbances in Arc-Jet Testing

This paper reports computational simulations and analysis of flow characterization tests in high enthalpy arc-jet facilities at NASA Ames Research Center. Flow disturbances and their characterization are presented through case studies, and their implications for testing are discussed. The case studies include free-jet test configurations from three different arc heaters and nozzles: the 60-MW IHF 30-inch conical nozzle, the 10-MW TP3 15-inch conical nozzle, and the 20-MW AHF 12-inch conical nozzle. For all cases, test articles are placed in the jet exiting the conical nozzle, and the existence of a flow disturbance is confirmed through flow survey data using pitot pressure and heat flux probes and accompanying analysis. The paper focuses on cases where the effects of the disturbances on test article surface quantities are important. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzles and test box, including the models tested. Comparisons of computations with the experimental measurements are presented. The computations that reproduce the probe sweep data approximately are essential to interpret the arc-jet test data accurately, while providing insights into several observed calorimeter anomalies caused by flow disturbances.

Arc-jet testing↗

DOE OSTI2021

A Scoping Study to Develop a Computational Fluid Dynamics Based Model to Predict Radiological Materials Packaging Temperatures within a Generic Staging Building

The objective of this work is to perform computational fluid dynamics simulations of a ventilated radiological-material-package staging building to determine the effect of including or excluding a variety of physical effects and computational methods on both the predicted package temperatures and required computational resources. The generic building contains 640 drum-packages containing heat-generating, radiological material supported on four racks that are eight levels tall. Additionally, there is a forced ventilation system, lighting, and insulated walls. Computational models were constructed that included or excluded (a) shelving, (b) effects of unsteadiness, and (c) radiation heat transfer. Simulations with each drum modeled separately were compared to simpler simulations with sets of four drums represented by an equivalent box-package. These models employed between 106 to 107 elements. Steady state simulations predicted package temperatures that were within 0.1°C of simulations that included transient effects and required only one-eighth the computational resources. Calculations that excluded shelving predicted temperatures within 0.6°C of simulations that included shelving and required one-fourth the computational resources. Excluding radiation heat transfer systematically increased temperatures by around 1.5°C but reduced computational resources by a factor of four. The box-package model reduced the computational resources by a factor of 3, but systematically predicted higher temperatures by around 1.1°C. These results will be used to develop an efficient computational fluid dynamics model to assess the ability of different staging building designs to prevent the temperature of package components from exceeding specified limits.

42 ENGINEERING↗

NASA NTRS2019

Extraction of Dynamic Inflow Models for Coaxial and Tandem Rotors from CFD Simulations

The dynamic inflow coupling with rotor/body dynamics is crucial in the analysis of stability and control law design for helicopters. Over the past several decades, finite-state inflow models for single rotor configurations in hover, forward flight, and maneuver have developed (Ref.1-3). By capturing the interference effects between rotors, the extension of pressure potential finite state inflow model has promising result for coaxial rotor configuration (Ref.4-6). Recently, the focus of the dynamic inflow modeling has shifted to tandem rotor configurations (Ref.7, 8). The development of the dynamic inflow models for tandem rotor configuration still have some limitations due to the lack of knowledge of rotor-to-rotor interference, and rotor-wake interference. Experimental methods, and computational fluid dynamics methods are commonly used to understand the rotor performance and rotor airload variations, and measure or predict inflow velocity distributions at the rotor desk. The inflow distributions are subsequently used to improve the dynamic inflow models. Tandem rotor configurations have been studied experimentally and computationally for several decades (Ref.9-12). Sweet (Ref.10) observed that a tandem rotor with 76-percent-radius overlap required 14% more induced power at hovering condition, relative to an isolated rotor of equivalent disk area. Sweet also found that, above a shaft-to-shaft distance of 1.03 diameter, the performance of the tandem rotor was nearly the same as two isolated rotors. The objective of the present study is to apply computational fluid dynamics simulations of tandem rotors for the extraction of dynamic inflow models. The extended methodology is first validated by comparing the computed induced power against test data. Subsequently inflow distributions and wake structures are analyzed.

Extraction of Dynamic Inflow Models↗

NASA NTRS2016

Flow Characterization Studies of the 10-MW TP3 Arc-Jet Facility: Probe Sweeps

This paper reports computational simulations and analysis in support of calibration and flow characterization tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted in the NASA Ames 10-MW TP3 facility using flat-faced stagnation calorimeters at six conditions corresponding to the steps of a simulated flight heating profile. Data were obtained using a conical nozzle test configuration in which the models were placed in a free jet downstream of the nozzle. Experimental surveys of arc-jet test flow with pitot pressure and heat flux probes were also performed at these arc-heater conditions, providing assessment of the flow uniformity and valuable data for the flow characterization. Two different sets of pitot pressure and heat probes were used: 9.1-mm sphere-cone probes (nose radius of 4.57 mm or 0.18 in) with null-point heat flux gages, and 15.9-mm (0.625 in) diameter hemisphere probes with Gardon gages. The probe survey data clearly show that the test flow in the TP3 facility is not uniform at most conditions (not even axisymmetric at some conditions), and the extent of non-uniformity is highly dependent on various arc-jet parameters such as arc current, mass flow rate, and the amount of cold-gas injection at the arc-heater plenum. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested. Comparisons of computations with the experimental measurements show reasonably good agreement except at the extreme low pressure conditions of the facility envelope.

Goekcen, Tahir↗

NASA NTRS2016

Flow Characterization Studies of the 10-MW TP3 Arc-Jet Facility: Probe Sweeps

This paper reports computational simulations and analysis in support of calibration and flow characterization tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted in the NASA Ames 10-MW TP3 facility using flat-faced stagnation calorimeters at six conditions corresponding to the steps of a simulated flight heating profile. Data were obtained using a conical nozzle test configuration in which the models were placed in a free jet downstream of the nozzle. Experimental surveys of arc-jet test flow with pitot pressure and heat flux probes were also performed at these arc-heater conditions, providing assessment of the flow uniformity and valuable data for the flow characterization. Two different sets of pitot pressure and heat probes were used: 9.1-mm sphere-cone probes (nose radius of 4.57 mm or 0.18 in) with null-point heat flux gages, and 15.9-mm (0.625 in) diameter hemisphere probes with Gardon gages. The probe survey data clearly show that the test flow in the TP3 facility is not uniform at most conditions (not even axisymmetric at some conditions), and the extent of non-uniformity is highly dependent on various arc-jet parameters such as arc current, mass flow rate, and the amount of cold-gas injection at the arc-heater plenum. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested. Comparisons of computations with the experimental measurements show reasonably good agreement except at the extreme low pressure conditions of the facility envelope.

CFD↗

NASA NTRS2014

Parallel Adaptive High-Order CFD Simulations Characterizing Cavity Acoustics for the Complete SOFIA Aircraft

This paper presents one-of-a-kind MPI-parallel computational fluid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft of a Boeing 747SP. These simulations focus on how the unsteady flow field inside and over the cavity interferes with the optical path and mounting of the telescope. A temporally fourth-order Runge-Kutta, and spatially fifth-order WENO-5Z scheme was used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh refinement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32,000 cores and 4 billion cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregularities caused by the highly complex geometry. Limits to scaling beyond 32K cores are identified, and targeted code optimizations are discussed.

Acoustics↗

NASA NTRS2023

Anisotropic Mach Cone Aligned Mesh Adaptation for Low Boom Simulations

An automated off-body Mach cone aligned structured curvilinear grid generation procedure is presented for near-field computational fluid dynamics simulations. This procedure combines output-based indicators and mesh redistribution to perform anisotropic mesh adaptation while maintaining Mach cone alignment. Automation is achieved through a novel direction-based adaptation indicator formulation. The adaptation procedure is demonstrated on the JAXA Wing Body geometry and X-59 C608 demonstrator model from the Second and Third AIAA Sonic Boom Prediction Workshops, respectively. It is demonstrated that anisotropic mesh adaptation may result in a greater than fifty percent reduction in resource usage required to achieve the same level of accuracy as uniform and user constructed Mach cone aligned grids for near-field pressure signatures, ground-level overpressure signatures, and loudness metrics.

CST↗

NASA NTRS2023

Anisotropic Mach Cone Aligned Mesh Adaptation for Low Boom Simulations

An automated off-body Mach cone aligned structured curvilinear grid generation procedure is presented for near-field computational fluid dynamics simulations. This procedure combines output-based indicators and mesh redistribution to perform anisotropic mesh adaptation while maintaining Mach cone alignment. Automation is achieved through a novel direction-based adaptation indicator formulation. The adaptation procedure is demonstrated on the JAXA Wing Body geometry and X-59 C608 demonstrator model from the Second and Third AIAA Sonic Boom Prediction Workshops, respectively. It is demonstrated that anisotropic mesh adaptation may result in a greater than fifty percent reduction in resource usage required to achieve the same level of accuracy as uniform and user constructed Mach cone aligned grids for near-field pressure signatures, ground-level overpressure signatures, and loudness metrics.

CST↗

NASA NTRS2015

Visualization of a Sweeping Jet by Laser Speckle Retro-reflective Background Oriented Schlieren

Sweeping jet actuators are currently being studied as a method of active flow control to increase the performance of wings and tails. Computational Fluid Dynamics simulations show that the jet does not sweep at a constant rate, instead spending more time at the minimum and maximum of its extent. As part of a test done to provide dynamic flow measurements for Computation Fluid Dynamics validation in NASA Ames Fluid Mechanics Lab, Particle Image Velocimetry (PIV) was performed on a single sweeping jet actuator designed by NASA and Boeing that was inset in a small airfoil. However, even at low jet pressures, the PIV could not be acquired at high enough frequencies to both spatially and temporally resolve the stickiness of the jet motion. Therefore, a small portion of the PIV laser was picked off and used to acquire flow visualization at higher frequencies through a new technique, laser speckle Retroreflective Background Oriented Schlieren (RBOS).

sweeping jet actuator↗

NASA NTRS2015

Computational Simulations of the 10-MW TP3 Arc-Jet Facility Flow

This paper reports computational simulations and analysis in support of calibration tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using stagnation calorimeters and two different blunted wedge models with calibration plates at a wide range of conditions in the NASA Ames 10-megawatt TP3 (Test Position 3) facility. Data were obtained using four different conical nozzles with the same test configuration in which the models were placed in a free jet downstream of the nozzle. Experimental surveys of arc-jet test flow with pitot and null-point heat flux probes were also performed at several arc-heater conditions, providing assessment of the flow uniformity and valuable data for the flow characterization. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested, and comparisons with the experimental measurements. These computational simulations provide estimates of the arc-jet test environment parameters that are not measured but are needed to evaluate the performance of thermal protection system materials, along with further valuable insights into the arc-jet testing environment. Simulation results are used to estimate centerline total enthalpy, surface shear, boundary layer thickness, and boundary layer edge Mach number and to verify that specific test requirements from the Orion program are met.

CFD↗

DOE OSTI2026

Hydrodynamic Modelling for Tidal Energy Kites: Cooperative Research and Development Final Report, CRADA Number CRD-24-30123

The primary goal of the work is to improve accuracy of hydrodynamic loads acting on the lifting surface of the kite turbine due to complex air-water-lifting body interaction in a low-order numerical model using high-fidelity numerical input. The result of this study will help the Participant better design control algorithms for the kite turbine so that the turbine can harvest more energy while maintaining stability of the platform as well as its reliability. The use of hydrodynamic coefficients/load components from computational fluid dynamic simulation will inform better inputs for a medium-fidelity dynamic model (specifically an OrcaFlex model) which is intended for use in rapid design iterations in the design process.

16 TIDAL AND WAVE POWER↗

NASA NTRSDate not supplied

CFD Simulations of the Truncated IHF Semi-Elliptical Nozzle Flow: Convective Calibration Plate

This paper reports computational simulations of a high-enthalpy panel test facility at NASA Ames Research Center. The recently developed facility configuration uses a truncated version of an existing semi-elliptical nozzle in the 60-MW Interaction Heating Facility, and it provides testing capability for panel test articles as large as 43.2 cm x 43.2 cm (17 in x 17 in). A water-cooled calibration plate, which is attached to the semi-elliptical nozzle bottom surface and instrumented with twenty Gardon-gage heat flux calorimeters and twelve pressure gages, is used to characterize the convective heating environment over the panel test articles. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, and flowfield over panel test articles. Both laminar and turbulent simulations are performed, and the computational results are compared with the calibration plate measurements. Comparisons of computed and measured heat fluxes indicate that the boundary layer in the semi-elliptical nozzle is transitional at certain conditions (relatively higher arc-heater pressures or mass flow rates).

Arc-jets↗