Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational fluid dynamics”

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 235 records · Page 13

Cryogenic Design and Thermal Analysis of EIC Central Detector (ePIC) Solenoid Magnet (MARCO)

The Electron Ion Collider (EIC) physics program utilizes a 2.0 T superconducting magnet at the heart of its ePIC detector system. This approximately 3.5 m long and 2.84 m diameter warm bore magnet has a 20 tons cold mass which is conduction-cooled using liquid helium at 4.5 K. A closed loop active thermosiphon system is chosen to facilitate the cooling and to maintain a minimum of 2 K temperature margin on the peak operating temperature (4.7 K) of the superconductor. Here, this paper presents the cryogenic design of the cooling system and the thermal analysis of the solenoid. A Computational fluid dynamics (CFD) model was developed to quantify the performance of the two–phase closed thermosiphon system and predict the temperature gradient on the cold mass.

Gopinath, Sandesh [Thomas Jefferson National Accel↗

Background-Oriented Schlieren Velocimetry of Helium Coolant Flow in Additively Manufactured Channels

High-pressure helium gas cooling is an attractive solution for thermal management of the fusion blanket first wall, as this coolant is chemically and neutronically inert and separable from hydrogenic species. However, due to the low thermal mass of helium, geometric optimization of these channels is required to provide sufficient cooling at manageable flow rates and pumping burdens. Increasingly, analysis and optimization of these coolant channels rely on computational fluid dynamics (CFD) simulations, and these require relevant experimental data for turbulence model validation. Toward this end, a high-pressure helium gas flow visualization system has been employed to image the flow of helium in flow channels with one-sided heating, mimicking the blanket first wall environment. Flow of helium at 4 MPa pressure and flow rates up to 68 g/s (Reynolds number 57 000) is supplied to rectangular channel test sections, with uniform heating applied to the bottom wall of the channel at heat fluxes varied between roughly 50 and 130 kW/m2. A high-speed camera is used to image index of refraction (IOR) gradients in the fluid via background-oriented schlieren (BOS), and temperature and pressure instrumentation are used to characterize thermal-hydraulic performance of each channel. Cross correlation of time-resolved BOS images is then used to calculate time-averaged 2-D helium velocity fields. Flow in additively manufactured (AM) channels is examined in this manner, including both featureless channels and those containing baffling as a heat transfer enhancement. The flow distribution seen in the featureless case differs significantly from that seen in prior simulations, whereas the flow in the baffled case shows the predicted behavior of flow forced along the heated wall. This augmented flow distribution is seen to increase the heat transfer coefficient in the baffled test section. Here, strategies are discussed for ongoing and future validation of these simulations, with the aim of model deployment for blanket cooling design and optimization.

Additive manufacturing↗

Flow and Heat Transfer Experimental Study for 3D-Printed Solar Receiving Tubes With Helical Fins at Internal Surface

3D-printing technology was applied to fabricate novel solar thermal collection tubes that have internal heat transfer enhancement fins and external surfaces with high solar absorptivity and low emissivity due to the ability to use different materials in one tube. Helical fins were selected to introduce circumferential flow and thus minimize the circumferential temperature difference of the tube that receives sunlight on one side. The structures of the helical fins were previously optimized from computational fluid dynamics (CFD) analysis with the objective of low entropy production rate by looking for high heat transfer coefficient and relatively lower pressure loss. High-temperature alloy, Inconel-718, was used to 3D print the tubes, which can resist corrosion for the potential application of molten chloride salts as heat transfer fluid. Experimental tests were carried out using water as the heat transfer fluid with the high heat flux provided by a tubular furnace heater. The tested Reynolds number ranges from 3.9 × 10 3 to 6.1 × 10 4 . Heat transfer coefficients of up to 2.8 times that of the smooth tube could be obtained with the expense of increased pressure loss compared to that of the smooth tube. The total system entropy generation can be significantly reduced due to the benefit of heat transfer enhancement that is greater than the expenses of the increased pressure loss. The experimental results of the 3D-printed heat transfer tubes confirmed the CFD-based results of fin optimization. Furthermore, the novel heat transfer tube is recommended for application in concentrating solar power systems.

14 SOLAR ENERGY↗

Chemical Reactor Network Modeling of Ammonia Rich-Quench-Lean Combustion Using a Partially Stirred Reactor Approach

Ammonia is a promising alternative to hydrogen with high energy density and favorable storage and transport characteristics. However, low flammability and a propensity for high nitrogen oxide (NO x ) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NO x emissions with ammonia. In this approach, the rich stage serves to oxidize a portion of the fuel while thermally decomposing as much of the remaining ammonia as possible, generating hydrogen. In the second (lean) stage, air is rapidly introduced, burning out the hydrogen and residual ammonia. Two-stage RQL combustion of ammonia has been investigated in the open literature both experimentally and numerically. In general, idealized chemical reactor network (CRN) models predict NO x concentrations below those of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing nonadiabatic operation. The typical CRN model is comprised of a perfectly-stirred-reactor (PSR), followed by a plug-flow-reactor (PFR), meant to represent the flame, and postflame zones, respectively. In the two-stage RQL approach two PSR-PFR networks are arranged sequentially, corresponding to the rich and lean stages, with secondary air injection in between. In the authors' past work, this arrangement has demonstrated the significant sensitivity of exit NO x to the rich stage equivalence ratio, while the amount of secondary air injection was shown to be less critical. In this paper, the CRN model is extended to (1) include the impacts of heat loss and (2) utilize a partially-stirred-reactor (PaSR) approach to study the impacts of mixing on emissions performance. Varying amounts of heat loss are applied to the rich relaxation zone to understand emissions performance and changes to optimization of equivalence ratio and residence time. Premixed and nonpremixed configurations are considered in the rich stage PaSR, with varying degrees of mixing intensity to study the interaction between mixing, transport, and kinetic timescales. Critically, the impact of mixing between hot products and secondary air injection is studied to understand practical injector needs. Results show unburnt ammonia leaving the rich stage as a primary contributor to NO x emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have been shown to create conditions that are conducive to increased N 2 O formation in the lean stage. In conclusion, the results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia.

Combustion↗

Robust Solution Verification Experiments on Nonuniform Meshes

The activities of verification, validation, and uncertainty quantification (VVUQ) provide a comprehensive means to assess the credibility of computational models. Within VVUQ, solution verification assesses numerical errors and evaluates whether the simulation is sufficiently accurate for its intended applications. As computational modeling gains traction in the development of complex, high-consequence systems, the need for robust solution verification intensifies, particularly because experimental data for these systems are often limited. This work examines improvements in the robustness of Richardson extrapolation (RE), a method commonly used in solution verification to study the discretization error of computational models using a power law. Nonuniform mesh refinement is discussed alongside other pollutants that affect the robustness of the power law model. Maximum likelihood estimation (MLE) is proposed as a robust strategy to address the uncertainty generated by nonuniform mesh refinement. An exploratory computational fluid dynamics (CFD) study of a 2D planar Poiseuille flow is conducted to determine if nonuniform mesh noise can be modeled with this MLE approach for more robust RE.

Weinmeister, Justin [ORNL] (ORCID:0000000160090237↗

Design, Analysis, and Experimental Testing of Hydrogen Lean Direct Injection Nozzles at Elevated Pressure

Abstract There are many challenges of commissioning a hydrogen combustor into future gas turbine engines; especially regarding achieving emissions goals. Previously, Escudero et al. and Tran et al. conducted a study to adapt the liquid fuel Lean Direct Injection (LDI) concept from Jet-A to gaseous natural gas-hydrogen blends and pure hydrogen [1], [2]. Experimental data was collected at atmospheric conditions using a Box Behnken design of experiments. The design of experiments suggested that biasing the air split in favor of the inner air circuit and increasing the swirl strength of this inner air passage resulted in improved NOx emissions, while the inverse was true for stability, which was quantified by studying the lean blowoff point (LBO) [1], [2]. The trends revealed by the original experiment [1], [2] provided a design direction for further iterations of the experimental hardware. The study presented herein describes the further investigation of such LDI injectors through experimental methods and computational fluid dynamic (CFD) simulations at atmospheric conditions, which were used to identify potential flow behaviors driving enhanced emissions performance. Further evaluation of select injectors from both studies was then conducted at elevated pressures up to 6 atmospheres. The results from both experiments are presented in this study, which include flame observations, emissions measurements, and operational challenges. NOx emissions results are reported on a volume basis in ppmvd corrected to 15% O2 and corrected for fuel. A predictive model for relating NOx emissions to test conditions at atmospheric conditions show high significance to adiabatic flame temperature while little to no significance to fuel composition for the best performing configurations. The results illustrate the connection between atmospheric testing and testing elevated pressures. The design direction indicated by the initial tests and CFD results in promising configurations for implementation into a Multi-point LDI array.

08 HYDROGEN↗

A Simulated Evaluation of Powder Flowability Through a Partially Obstructed Consumable in Blown Powder Directed Energy Deposition Systems

Abstract In the interest of continued industrialization of metal additive manufacturing in modern production environments, cost is often referenced as a primary deterrent to new adopters. Conventional economic models for additive systems, processes, and supply chains often focus on specific process applications with little generalizability, or they neglect significant costs associated with production such as machine maintenance and consumable part replacement. Compounding the latter issue are substantial knowledge gaps in consumable part wear characterization for additive and other convergent manufacturing systems. In coaxial blown powder directed energy deposition systems, gas atomized metal powder is wasted during material deposition at a rate that is partly dependent on present wear phenomena in a consumable nozzle housed in the cladding head assembly. The price and lead time required to replace the nozzle incentivizes its reuse even when visibly worn. Often this initiates a process quality decline in the form of underbuilt geometry and internal defects due to losses in powder catchment efficiency. While depositing H13 steel using a hybrid manufacturing machine tool equipped with such a deposition system, a unique partial clog with a bridge-like structure formed at the consumable nozzle exit when supporting argon gas flows failed mid-process. To further understand coaxial multi-phase powder flow in the event of support gas failure, a computational fluid dynamics simulation is tailored to relevant process parameters, H13 powder material profile, and machine operator observations collected after the incident. The resulting differences in powder flow compared to control gas flow parameters is presented and discussed. The powder flowability and performance of the clogged nozzle is then assessed by using an optical profilometer to extract the profile of the clog and recreate the clog geometry within the simulation environment. In past work this simulation has been experimentally validated for a 316L steel powder material profile and used specifically for analyzing powder stream geometry and catchment efficiency. After the initial powder flow characterization, the clog is removed, and the nozzle is reprofiled. After removing the obstructing clog, the newly unobstructed nozzle geometry, the original off the shelf nozzle geometry, and additional nozzle profiles exploring different consumable refurbishment strategies are reevaluated in the simulation. Powder catchment efficiency for all variant nozzle geometries and relevant flow variables are compared and discussed, along with potential mitigation strategies for optimizing powder flowability with worn consumables. This work expands on the known morphology of blown powder obstructions and wear defects present in consumable coaxial nozzles while discussing pragmatic simulation driven responses to unanticipated subsystem failure in hybrid manufacturing machining platforms.

DeWitte, Lisa↗

High- and Mid-Fidelity Modeling Comparison for a Floating Marine Turbine System

There is a lack of suitable numerical tools, particularly open-source tools, that can be used for designing and optimizing marine turbine systems. The National Renewable Energy Laboratory has added features to their widely used mid-fidelity wind turbine modeling code, OpenFAST, to enable modeling of axial-flow marine turbines. This necessitated the addition of several physical effects relevant to marine turbines that are neglected for wind turbines. These include buoyancy, added mass and inertial loads, wave-current superposition, and changes to the coordinate systems. This updated version of OpenFAST allows for the modeling of both fixed and floating marine turbines at a speed comparable to real time. While efficient for large sets of load cases and design studies, mid-fidelity codes make simplifying assumptions that may impact their accuracy. High-fidelity computational fluid dynamics (CFD) simulations can capture more flow effects with fewer assumptions and provide detailed body pressure mapping and flow-field information. It is important to compare predictions between mid-fidelity and high-fidelity codes, both to verify the models and to understand the limitations. A floating marine turbine system was modeled both with OpenFAST and with the commercial CFD code STAR-CCM+. The CFD model used a three-dimensional unsteady Reynolds-averaged Navier-Stokes solver for a volume-of-fluid numerical wave and current tank. The blade-resolved simulations used the sliding-interface technique for the spinning rotor and an overset grid to accommodate the rigid-body motion of the floating system. The mooring system was modeled with a custom coupling of the CFD solver with the open-source code MoorDyn. This improves upon the existing quasi-static catenary solver in STAR-CCM+, which lacks seabed contact or line-to-line connections. Simulation results for a floating marine turbine are compared between OpenFAST and CFD, highlighting the capabilities of the mid-fidelity code and identifying the areas where a high-fidelity approach is needed.

CFD↗

HELIUM LEAK TEST MODELING OF A SPENT NUCLEAR FUEL CANISTER

The U.S. Department of Energy (DOE) is considering the development of one or more federal consolidated interim storage facilities (CISFs) to be used to store commercial spent nuclear fuel (SNF) at locations in the U.S. One of the first technical challenges of a CISF is performing an inspection of SNF canisters upon their receipt to confirm they can be placed into the CISF’s licensed storage configuration. The canister receipt inspection is critical to CISF site operations. The test is conceived as being a helium (He) leak check, intended to confirm that the confinement boundary of a SNF canister is intact. SNF canisters are filled with He when they are sealed, so detection of a He leak indicates that a through-wall flaw has occurred in the canister confinement boundary. Other measurements are planned to occur upon canister receipt in addition to the He leak check such as krypton-85 measurements, which would indicate confinement breaches of one or more fuel rods in addition to a breach of the SNF canister. However, the He leak check has been identified as one of such high importance and has such significant technical challenges that a full-scale demonstration is needed to confirm the He leak test’s viability and to assist in planning relative to its operational requirements. A modeling methodology for simulating the He detection test was developed to help inform the test plan and the design of the test vessels. To develop the modeling methodology a detailed computational fluid dynamics (CFD) benchmark model was constructed to compare against leak rate test data from a transportation package for radioactive material. This report is focused on modeling efforts to simulate the benchmark leak test.

Suffield, Sarah R.↗

Using 3D Measurements from MRI to Improve VVUQ

Magnetic resonance velocimetry (MRV) is a diagnostic that can measure 3D, three-component turbulent velocity fields in arbitrarily complex flow configurations. The approach uses magnetic resonance imaging (MRI), which is a system commonly available in radiology departments or medical research centers. MRV acquires signal from hydrogen protons in water channel flows. Despite the obvious utility of the measurements for investigating a variety of flows, comparisons with computational fluid dynamic (CFD) simulations largely retain traditional metrics, comparing velocity profiles and planes at discrete regions within the flow rather than using the 3D nature of the measured field. In this effort, MRV was conducted in a simple water channel with six spanwise-centered periodic obstacles roughly shaped as cubes. At a fully turbulent Reynolds number of 15,000, the channel includes two flows: a main streamwise flow and, in the wake of the second obstacle, an injected flow oriented perpendicular to the streamwise flow. The flow geometry includes partial obstacles on the side walls and is inherently 3D because of the interaction of the two streams and the wake features, making the flow challenging to completely measure with traditional optical techniques. Aside from line profiles and planar comparisons, two 3D metrics are used to compare the data with a steady Reynolds-averaged Navier–Stokes (RANS) simulation result. A brief discussion about the comparison is provided, including comments about uncertainty.

Benson, Mike [ORNL] (ORCID:000000023210116X)↗

CFDverify

Estimating the discretization error of computational fluid dynamics (CFD) or other scientific codes as part of solution verification is often a non-trivial part of the analysis process. Methods can be complicated and may include assumptions/qualifications that need to be checked during analysis. CFD analysts, therefore, are likely to make errors when trying to conduct this necessary analysis on their own in not knowing about the best method for their problem, not correctly implementing a method, or in not have diagnostic tools to determine if the method was correctly applied.

Weinmeister, Justin [Oak Ridge National Laboratory↗

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 ↗

quinoa v1.1

Adaptive computational fluid dynamics

Pandare, Aditya [Los Alamos National Laboratory]↗

CFD-Enhanced Calibration of a Multihole Probe for Small Uncrewed Aircraft Systems

A robust calibration method with a 9-hole probe (9HP) for inertial wind vector measurements from a small uncrewed aircraft system (sUAS) is presented. Calibration accuracy is improved by using computational fluid dynamics (CFD) to estimate corrections, such as those for test section blockage, that are generally applied in wind tunnel testing. A method for estimating experimental bias is presented to account for flow variations over the pressure taps as the 9HP is repositioned in the test section. Installation effects from flow over the airframe and upwash produced by the wings are estimated from CFD simulations of the entire airframe at the expected cruise speed. An initial CFD analysis has demonstrated a quantifiable linear relationship between the measured and actual angle of attack as measured from the relative wind frame. In conclusion, the objective of these additional steps is to increase the accuracy of 9HP calibration.

54 ENVIRONMENTAL SCIENCES↗

Development of prechamber enabled mixing-controlled combustion strategy for ultra-low methane emissions from lean burn natural gas engines

This numerical study explores the optimization of Prechamber Enabled Mixing-Controlled Combustion (PC-MCC) using natural gas in heavy-duty engines, aiming to enhance combustion efficiency to minimize methane slip and NOx emissions. The approach involves a prechamber ignition system, distinct from conventional spark ignition (SI) systems, to initiate combustion of direct injected natural gas. By leveraging the robust ignition characteristics of the prechamber, the PC-MCC method demonstrates significant potential in achieving efficient combustion akin to diesel engines but with lower greenhouse gas emissions. The research evaluates the effects of various geometric and operational parameters on the combustion process and emissions, including prechamber volume, nozzle diameter, direct injector (DI) geometry, and engine operating strategies. Computational Fluid Dynamics (CFD) simulations are utilized, focusing on a heavy-duty, single-cylinder engine modeled after the Caterpillar C9.3B engine. Key findings indicate that a prechamber volume of 3 cc, coupled with a nozzle diameter of 2.75 mm for two prechamber holes, strikes an optimal balance between combustion efficiency and emissions reduction. This configuration ensures robust combustion across a range of operating conditions while maintaining methane slip within targeted limits. Further investigation into DI geometry shows the significance of the injector umbrella angle and nozzle diameter in shaping the fuel-air mixing and combustion dynamics. An umbrella angle of 130° and a nozzle diameter of 300 microns are identified as optimal, promoting rapid and efficient combustion with minimized methane and NOx emissions. The study also investigates the impact of injection timing and pressure, highlighting their roles in controlling combustion timing and influencing emissions levels. Advanced injection timing is found to be crucial in achieving the desired low methane slip, whereas retarded injection timing assists to reduce NOx emissions while having a slight increase in methane emissions. Operating strategies incorporating various levels of Exhaust Gas Recirculation (EGR) are assessed for their effectiveness in further reducing emissions. The research demonstrates that a judicious combination of internal hot EGR and careful calibration of DI pressure and SOI timing can achieve significant reductions in NOx emissions while keeping methane slip under control. Specifically, an internal EGR level of 15%–25%, combined with DI pressures of 200–300 bar and injection timings at or after top dead center, is recommended. These findings contribute valuable insights into the development of advanced combustion techniques for natural gas engines, offering a viable pathway to reduce methane slip without compromising engine efficiency or performance. The PC-MCC system presents a promising solution for the future of heavy-duty natural gas engine technology to reduce methane emissions.

Nsaif, Osama↗

Development of a two-zone zero-dimensional mixing model for actively fueled prechamber engines

In this numerical work, a zero-dimensional, two-zone mixing model was developed and validated to capture the mixture preparation process of an actively fueled prechamber engine. The model was developed for ethanol-gasoline fuel blends but is applicable to various fuels of interest relevant to prechamber applications. The model was exercised to demonstrate the effects of boundary conditions, fuel composition, and residual content on the mixture stoichiometry predictions. The 0D model showed reasonable agreement in prechamber equivalence ratio predictions with higher fidelity computational fluid dynamics modeling, but the accuracy of model was constrained to the governing assumptions of the injection process. The modeling assumptions surrounding the non-dimensional injection of mass into the prechamber were scrutinized to identify key areas of model improvement. The results indicate heat transfer effects play a crucial role in the evaporation of the fuel within the confined auxiliary chamber, especially when considering fuels with high latent heat of vaporization. Additionally, the presence of fuel films that collect on the prechamber surfaces complicated the modeled evaporative cooling loses as the modeling suggests the films evaporate primarily by means of conduction from the boundary. Furthermore, the developed model serves as a simplified framework to generate calibration maps of prechamber fueling or run real time as an open loop control strategy such that desired prechamber mixture proportions are effectively targeted.

Zero-dimensional↗

Forest structural complexity and ignition pattern influence simulated prescribed fire effects

Background: Forest structural characteristics, the burning environment, and the choice of ignition pattern each influence prescribed fire behaviors and resulting fire effects; however, few studies examine the influences and interactions of these factors. Understanding how interactions among these drivers can influence prescribed fire behavior and effects is crucial for executing prescribed fires that can safely and effectively meet management objectives. To analyze the interactions between the fuels complex and ignition patterns, we used FIRETEC, a three-dimensional computational fluid dynamics fire behavior model, to simulate fire behavior and effects across a range of horizontal and vertical forest structural complexities. For each forest structure, we then simulated three different prescribed fires each with a unique ignition pattern: strip-head, dot, and alternating dot. Results: Forest structural complexity and ignition pattern affected the proportions of simulated crown scorch, consumption, and damage for prescribed fires in a dry, fire-prone ecosystem. Prescribed fires in forests with complex canopy structures resulted in increased crown consumption, scorch, and damage compared to less spatially complex forests. The choice of using a strip-head ignition pattern over either a dot or alternating-dot pattern increased the degree of crown foliage scorched and damaged, though did not affect the proportion of crown consumed. We found no evidence of an interaction between forest structural complexity and ignition pattern on canopy fuel consumption, scorch, or damage. Conclusions: We found that forest structure and ignition pattern, two powerful drivers of fire behavior that forest managers can readily account for or even manipulate, can be leveraged to influence fire behavior and the resultant fire effects of prescribed fire. These simulation findings have critical implications for how managers can plan and perform forest thinning and prescribed burn treatments to meet risk management or ecological objectives.

54 ENVIRONMENTAL SCIENCES↗

High-Fidelity Numerical Wave Tank Verification & Validation Study: Wave Generation Through Paddle Motion: Preprint

This paper presents a numerical benchmark study of wave propagation due to a paddle motion using different high-fidelity numerical models, which are capable of replicating the nearly actual physical wave tank testing. A full time series of the measured wave generation paddle motion which was used to generate wave propagation in the physical wave tank will be utilized in each of the models contributed by IEA OES Task 10's participants, which includes both computational fluid dynamics (CFD) and smooth hydrodynamic particle (SPH). The high-fidelity simulations of the physical wave testcase will allow for the evaluation of the initial transient effects from wave ramp-up and its evolution in the wave tank over time for two representative regular waves with varying levels of nonlinearity. A couple of interesting metrics like the predicted wave surface elevation at select wave probes, wave period, and phase-shift in time will be assessed to evaluate the relative accuracy of numerical models versus experimental data within specified time intervals. These models will serve as a guide for modelers in the wave energy community and provide a base case to allow further and more detailed numerical modeling of the fixed Kramer Sphere Cases under wave excitation force wave tank testing.

HYDRO ENERGY,TIDAL AND WAVE POWER↗