Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD modeling”

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 271 records · Page 15

Computational Fluid Dynamics Uncertainty Analysis Applied to Heat Transfer over a Flat Plate

There have been few discussions on using Computational Fluid Dynamics (CFD) without experimental validation. Pairing experimental data, uncertainty analysis, and analytical predictions provides a comprehensive approach to verification and is the current state of the art. With pressed budgets, collecting experimental data is rare or non-existent. This paper investigates and proposes a method to perform CFD uncertainty analysis only from computational data. The method uses current CFD uncertainty techniques coupled with the Student-T distribution to predict the heat transfer coefficient over a at plate. The inputs to the CFD model are varied from a specified tolerance or bias error and the difference in the results are used to estimate the uncertainty. The variation in each input is ranked from least to greatest to determine the order of importance. The results are compared to heat transfer correlations and conclusions drawn about the feasibility of using CFD without experimental data. The results provide a tactic to analytically estimate the uncertainty in a CFD model when experimental data is unavailable

Computational↗

Data Analysis Techniques for Fan Performance in Highly-Distorted Flows from Boundary Layer Ingesting Inlets

The design of a unique distortion-tolerant fan for a high-bypass ratio boundary-layer ingesting propulsion system has been completed and a rig constructed and tested in the NASA Glenn 8’x6’ wind tunnel. Processing the data from the experiment presented some interesting challenges because of the complexity of the experimental setup and the flow through the test rig. The experiment was run in three phases, each of which employed a unique complement of inlet throat and fan face instrumentation to avoid the blockage that would have resulted from simultaneously installing all of the rakes. The measurement from the individual test points were subsequently combined to compute the overall stage performance. A CFD model of the experiment was used to gain understanding of the flow field and to test some of the techniques proposed for interpolating and extrapolating the measurements into regions where measurements were not made. This capability became extremely useful when it was discovered that there was an unexpected total temperature distortion in the tunnel. The CFD model was modified by inserting a total temperature profile at the upstream boundary that mimicked the measured distortion where measurements were available and that CFD solution was used to investigate methods to infer the complete total temperature field at the fan face.

Hardin, Larry W.↗

Advanced CO 2 Capture Solvent Systems for Dynamic Power Generation: Quarterly Research Performance Progress Report, QR4 (Q4FY24)

We developed an integrated Computational Fluid Dynamics (CFD) model to simulate the multi-physics coupled cooling process of mixed gas by cold water within a Direct Contact Cooler (DCC) equipped with a rotating packing bed (RPB). The model captures the interactions between fluid dynamics, heat transfer, mass transport, and phase transitions, while accounting for key operational variables such as RPB rotational speed and the mass flow rates of both liquid and gas. The CFD model has been validated using experimental data, specifically by comparing predicted outflow gas and liquid temperatures to measured results. Our findings demonstrate the significant effects of RPB rotational speed and mass flow rates on cooling performance, providing valuable insights for optimizing DCC efficiency in industrial applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Coupling coarse mesh CFD with fine mesh CFD for modeling for Molten Salt Reactors in the Virtual Test Bed VTB

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims at developing a simulation tool kit to accelerate the development and deployment of nuclear power technologies. NEAMS multiphysics tools have been designed to provide numerical simulation support for the design and licensing of GEN IV reactors. Pronghorn is NEAMS’s coarse-mesh computational fluid dynamics (CFD) tool, which is designed to run 3D core transients in GEN IV reactors at a reduced computational cost. To increase their accuracy, coarse-mesh CFD simulations require calibrated closure coefficients. One way of computing these coefficients is via the Nek5000, NEAMS’s high-fidelity CFD tool. This article discusses our current research lines in informing Pronghorn closure coefficients via Nek5000 to enable multiphysics simulations of the core cavity of the molten-salt fast reactor. We present an application in which Pronghorn mixing length turbulent viscosity has been calibrated from Nek5000 simulations. The resulting Pronghorn thermal-hydraulics model is then coupled to Griffin, the NEAMS neutron transport solver, to solve for the steady-state configuration of the molten-salt fast reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Prediction of the Thrust Performance and the Flowfield of Liquid Rocket Engines

In an effort to improve the current solutions in the design and analysis of liquid propulsive engines, a computational fluid dynamics (CFD) model capable of calculating the reacting flows from the combustion chamber, through the nozzle to the external plume, was developed. The Space Shuttle Main Engine (SSME) fired at sea level, was investigated as a sample case. The CFD model, FDNS, is a pressure based, non-staggered grid, viscous/inviscid, ideal gas/real gas, reactive code. An adaptive upwinding differencing scheme is employed for the spatial discretization. The upwind scheme is based on fourth order central differencing with fourth order damping for smooth regions, and second order central differencing with second order damping for shock capturing. It is equipped with a CHMQGM equilibrium chemistry algorithm and a PARASOL finite rate chemistry algorithm using the point implicit method. The computed flow results and performance compared well with those of other standard codes and engine hot fire test data. In addition, the transient nozzle flowfield calculation was also performed to demonstrate the ability of FDNS in capturing the flow separation during the startup process.

Wang, T.-S.↗

Experimentation and Modeling of Jet A Thermal Stability in a Heated Tube

High performance aircraft typically use hydrocarbon fuel to regeneratively cool the airframe and engine components. As the coolant temperatures increase, the fuel may react with dissolved oxygen forming deposits that limit the regenerative cooling system performance. This study investigates the deposition of Jet A using a thermal stability experiment and computational fluid dynamics (CFD) modeling. The experimental portion of this study is performed with a high Reynolds number thermal stability (HiRets) tester in which fuel passes though an electrically heated tube and the fuel outlet temperature is held constant. If the thermal stability temperature of the fuel is exceeded, deposits form and adhere to the inside of the tube creating an insulating layer between the tube and the fuel. The HiRets tester measures the tube outer wall temperatures near the fuel outlet to report the effect of deposition occurring inside the tube. Final deposits are also estimated with a carbon burn off analysis. The CFD model was developed and used to simulate the fluid dynamics, heat transfer, chemistry, and transport of the deposit precursors. The model is calibrated to the experiment temperature results and carbon burn-off deposition results. The model results show that the dominant factor in deposition is the heated wall temperature and that most of the deposits are formed in the laminar sublayer. The models predicted a 7.0E-6 kilograms per square meter-sec deposition rate, which compared well to the carbon burn-off analysis deposition rate of 1.0E-6 kilograms per square meter-sec.

Khodabandeh, Julia W.↗

Quench process modeling & simulation in the heat-treatment of critical aerospace components

To manufacture light-weight, advanced metal alloy components for gas turbine engines, quench heat-treatment processes are typically used. By quenching the component from elevated temperatures, the alloy sometimes undergoes a solid-state phase transformation which produces special microstructures with the required, enhanced mechanical properties. However, the quenching can also lead to cracks forming in the component. Addressing the quench cracking problems adds a significant burden to the cost, schedule, and energy demand of manufacture. Currently, optimizing the quench process to mitigate or avoid the cracking is performed largely by trial-and-error, relying heavily on costly experimental (thermocouple)trials to understand the local thermal gradients which cause the cracks to form. In this first part (Phase 1) of the work, high-performance computing is employed to establish the ability of modern CFD (computational fluid dynamics) to alleviate or wholly replace the experimental quenching trials by virtual testing. A Baseline CFD model is defined and its accuracy established to be comparable to(and which usually exceeds) the accuracy of existing HTC (heat-transfer coefficient) based simulation methods of quenching. As a first-principles based approach, “calibration” of the Baseline CFD model is independent of the quench process itself, but instead relies on the accuracy of the underlying (modeled),generic two-phase fluid processes which cannot be currently resolved by CFD for large, industrial-scale cases. A novel, high-fidelity DNS capability has been developed and verified to examine and further improve upon the mean-field closure submodels on which the Baseline CFD approach is based.1

99 GENERAL AND MISCELLANEOUS↗

Quench process modeling & simulation in the heat-treatment of critical aerospace components

To manufacture light-weight, advanced metal alloy components for gas turbine engines, quench heat-treatment processes are typically used. By quenching the component from elevated tempera-tures, the alloy sometimes undergoes a solid-state phase transformation which produces special microstructures with the required, enhanced mechanical properties. However, the quenching can also lead to cracks forming in the component. Addressing the quench cracking problems adds a significant burden to the cost, schedule, and energy demand of manufacture. Currently, optimizing the quench process to mitigate or avoid the cracking is performed largely by trial-and-error, relying heavily on costly experimental (thermocouple) trials to understand the local thermal gradients which cause the cracks to form. In this first part (Phase 1) of the work, high-performance computing is employed to establish the ability of modern CFD (computational fluid dynamics) to alleviate or wholly replace the experimental quenching trials by virtual testing. A baseline CFD model is defined and its accuracy established to be comparable to (and which usually exceeds) the accuracy of existing HTC (heat-transfer correlation) based simulation methods of quenching. As a first-principles based approach, “calibration” of the Baseline CFD model is independent of the quench process itself, but instead relies on the accuracy of the underlying (modeled), generic two-phase fluid processes which cannot be currently resolved by CFD for large, industrial-scale cases. A novel, high-fidelity DNS capability has been developed and verified to examine and further improve upon the mean-field closure submodels on which the Baseline CFD approach is based.

36 MATERIALS SCIENCE↗

Coupling coarse-mesh CFD with fine-mesh CFD for modeling molten-salt reactors in the Virtual Test Bed (VTB)

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims at developing a simulation tool kit to accelerate the development and deployment of nuclear power technologies. NEAMS multiphysics tools have been designed to provide numerical simulation support for the design and licensing of GEN IV reactors. Pronghorn is NEAMS's coarse-mesh computational fluid dynamics (CFD) tool, which is designed to run 3D core transients in GEN IV reactors at a reduced computational cost. To increase their accuracy, coarse-mesh CFD simulations require calibrated closure coefficients. One way of computing these coefficients is via the Nek5000, NEAMS's high-fidelity CFD tool. This article discusses our current research lines in informing Pronghorn closure coefficients via Nek5000 to enable multiphysics simulations of the core cavity of the molten-salt fast reactor. We present an application in which Pronghorn mixing length turbulent viscosity has been calibrated from Nek5000 simulations. The resulting Pronghorn thermal-hydraulics model is then coupled to Griffin, the NEAMS neutron transport solver, to solve for the steady-state configuration of the molten-salt fast reactor. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Computational study on the impact of gasoline-ethanol blending on autoignition and soot/NO x emissions under low-load gasoline compression ignition conditions

Here, in the present work, computational fluid dynamics (CFD) simulations of a single-cylinder gasoline compression ignition (GCI) engine are performed to investigate the impact of gasoline-ethanol blending on autoignition, nitrogen oxide (NO x ), and soot emissions under low-load conditions. In order to represent the test gasoline (RD5-87), a four-component toluene primary reference fuel (TPRF)+ethanol (ETPRF) surrogate (with 10% ethanol by volume; E10) is employed. A three-dimensional (3D) engine CFD model employing finite-rate chemistry with a skeletal kinetic mechanism (including NO x sub-mechanism), adaptive mesh refinement (AMR), and hybrid method of moments (HMOM) is adopted to capture the in-cylinder combustion phenomena and soot/NO x emissions. The engine CFD model is validated against experimental data for three gasoline-ethanol blends: E10, E30 and E100, with varying ethanol content by volume. Model validation is carried out for a broad range of start-of-injection (SOI) timings (−21, −27, −36, and −45 crank angle degrees (°CA) after top-dead-center (aTDC)) with respect to in-cylinder pressure, heat release rate, combustion phasing, NO x and soot emissions. For relatively later injection timings (−21 and −27 °CA aTDC), E30 yields higher amount of soot than E10; while the trend reverses for early injection cases (−36 and −45 °CA aTDC ). On the other hand, E100 yields the lowest amount of soot among all fuels irrespective of SOI timing. Further, E10 shows a non-monotonic trend in soot emissions with SOI timing: SOI-36>SOI-45>SOI-21>SOI-27, while soot emissions from E30 exhibit monotonic decrease with advancing SOI timing. NO x emissions from various fuels follow a trend of E10>E30>E100. On the other hand, NO x emissions increase as SOI timing is advanced for all fuels, with an anomaly for E10 and E100 where NO x decreases when SOI is advanced beyond −36 °CA aTDC. Detailed analysis of the numerical results is performed to investigate the soot/NO x emission trends and elucidate the impact of chemical composition and physical properties on autoignition and emissions characteristics.

Computational fluid dynamics↗

Particle Size and Density Measurement

This project has developed a miniature gravity settling classifier (flow cell) and flow system that separates slurry particles by settling rate into six different Cuts. The settling rates for particles in each Cut were determined from COMSOL computational fluid dynamics (CFD) modeling of the flow in the flow cell together with particle tracking for particles of differing settling rates. For particles collected in each Cut, the particle concentration and size distribution were measured. The effective particle densities from each Cut are calculated from their respective CFD modeled settling rates using a relationship between the unhindered settling rate, particle size distribution (PSD), and effective density. The particle size and density distribution (PSDD) of the slurry sample is then calculated by combining the size and density distributions determined for the particles collected in each of the settling rate Cuts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Computational Fluid Dynamics Simulations of Glass Vitrification Refractory Coupon Tests

The Waste Treatment and Immobilization Plant (WTP) at the Hanford site is nearing the start of the Direct-Feed Low-Activity Waste (DFLAW) operations. DFLAW is destined to convert a pretreated low activity waste portion of the 56 million gallons of tank waste into a stable solid glass. In the subsequent decade completion of the high-level waste (HLW) facility is anticipated. Sustained operational missions of both LAW and HLW melter facilities are expected over multiple decades. In high-temperature glass melters, the refractory lining corrodes over time, which could potentially be an issue for longer term operations, this refractory corrosion is higher at the level of the glass-air interface due to surface tension driven flow. The glass viscosity, melt pool temperature, and glass chemical composition can impact the rate at which the refractory corrodes. This rate is important to quantify for the various waste glasses to be produced at the WTP since the integrity of the refractory should not be a limiting factor affecting the lifetime of the melter. To this end, a series of glasses representative of the first batches of waste glass produced by the WTP will be melted in small-scale crucibles with Monofrax® K-3 coupons inserted. The corrosion of the K-3 will be measured in the melt and at the meltline (or neckline). A model for the corrosion rate will be constructed and implemented into a previously developed framework for a computational fluid dynamics (CFD) model of the full-scale WTP. To assist with experimental design and validate the implementation of the model in the full-scale melter, CFD simulations of the small-scale crucible tests were performed. The bubbling that occurs in the small-scale crucible is initially validated here with a model that uses silicone oil at room temperature. The viscosity of the oil ranges from 1 to 100 Pa•s, which corresponds to operating glass pool temperatures near 1150 °C down to idling temperatures near 950 °C. The simulation results show good agreement with the bubble sizes that form during experiments. CFD modeling of the crucible setup was used to determine bubbling characteristics to match the range of near-wall velocities expected in the full-scale WTP. This study presents the initial CFD modeling results, corrosion testing plan, and some preliminary corrosion samples with an outline for the next steps for the development of the corrosion model.

Abboud, Alexander W. [Idaho National Lab]↗

Grid generation and surface modeling for CFD

When computing the flow around complex three dimensional configurations, the generation of the mesh is the most time consuming part of any calculation. With some meshing technologies this can take of the order of a man month or more. The requirement for a number of design iterations coupled with ever decreasing time allocated for design leads to the need for a significant acceleration of this process. Of the two competing approaches, block-structured and unstructured, only the unstructured approach will allow fully automatic mesh generation directly from a CAD model. Using this approach coupled with the techniques described in this paper, it is possible to reduce the mesh generation time from man months to a few hours on a workstation. The desire to closely couple a CFD code with a design or optimization algorithm requires that the changes to the geometry be performed quickly and in a smooth manner. This need for smoothness necessitates the use of Bezier polynomials in place of the more usual NURBS or cubic splines. A two dimensional Bezier polynomial based design system is described.

Connell, Stuart D.↗

Heat transfer enhancement due to cold cap motion from bubbling in a waste glass melter

In this study, a computational fluid dynamics (CFD) model was developed to model the motion of a solid cold cap in a waste glass melter. Forced convection bubblers at the base of the melter release air into the molten glass, which forms large bubbles that travel upwards to the cold cap and augment heat transfer from the glass to the cold cap. The CFD model employs the Navier- Stokes equations to solve for the fluctuating flowfield using a rigid body motion dynamic fluid body interaction module. This allows for movement of the floating body in response to the bubbling forces calculated at each time step. Further, the heat flux delivered to the cold cap by the convective bubbling is studied as a function of the bubbling rate. Results for the moving cold cap are compared with the computed heat flux trends for a stationary cold cap. The heat flux delivered to the cold cap from the glass is 25% higher for the case with the moving cold cap. The heat flux was found to be proportional to v0.6 as opposed to v0.9 (where v is the normalized bubbling rate) for the stationary cold cap.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Challenge Problem 1: Preliminary Model Development and Assessment of Flexible Heat Transfer Modeling Approaches

This report presents the modeling progress of a university consortium studying flexible modeling for turbulent heat transfer. In the context of advanced nuclear reactor design, where challenges include non-unity Prandtl fluids, free and mixed convection, and transitional flow, computational fluid dynamics methods are underutilized due to low confidence in modeling approaches and high computational cost. This report evaluates the performance of engineering CFD models in diabatic, buoyant turbulent flow. It finds that all turbulence models including several variants of the k-ε and k-ω models struggle to predict accurate turbulent momentum and heat transfer in such flows. The Nusselt numbers have been compared between the models and the DNS data, where calculations have been performed for each case and trends have shown a good agreement between DNS estimated Nusselt numbers and available correlations as well as experimental data. A novel DNS correlation for high Pr cases as they are transitioning from mixed convection to natural convection has been proposed. While Nusselt number errors relative to DNS range from 20% to 50%, the models capture similar trends to DNS with respect to Nusselt suppression and amplification under varying levels of buoyancy effect. This report also highlights model form error as a significant contributor to CFD predictions and proposes a framework for quantifying model error and improving confidence in CFD calculations. To improve the predictive capability of engineering CFD models, data-driven approaches for turbulence models are investigated. Theoretical frameworks based on the invariant tensor / vector basis neural networks for prediction of Reynolds stresses and turbulent heat fluxes are employed. The models are developed using direct numerical simulations data for forced convection flows of different fluids in vertical planar channel domain. The framework is implemented in spectral element solvers Nek5000 / nekRS and has shown a potential for future development and consideration of mixed convection flows.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Burner Rig with an Unattached Duct for Evaluating the Erosion Resistance of Thermal Barrier Coatings

Extensive computational fluid dynamics (CFD) modeling backed by experimental observation has demonstrated the feasibility of using an unattached duct to increase the velocity and spatial spread of erodent particles exiting from a burner rig. It was shown that gas velocity and temperature are mostly retained if the inner diameter of the unattached duct equaled the exit diameter of the burner rig nozzle. For particles having a mean diameter of 550 millimeters, the modeled velocity attained at a distance 2.0 in. (50.8 millimeters) beyond the exit of a 12 in. (305 millimeters) long duct was approximately twice as large as the velocity the same distance from the nozzle when the duct was not present. For finer particles, the relative enhancement was somewhat less approximately 1.5 times greater. CFD modeling was also used to guide the construction of a device for slowing down the velocity of the particles being injected into the burner rig. This device used a simple 45 degree fitting to slow the particle velocity in the feed line from 20 meters per second, which is in the range needed to convey the particles, to about 3 meters per second just as they are injected into the burner. This lower injection velocity would lessen the severity of the collision of large particles with the wall of the burner liner opposite the injection port, thereby reducing potential damage to the burner liner by high-velocity particles.

Miller, Robert A.↗

Hypersonic Inflatable Aerodynamic Decelerator (HIAD) Technology Development Overview

The successful flight of the Inflatable Reentry Vehicle Experiment (IRVE)-3 has further demonstrated the potential value of Hypersonic Inflatable Aerodynamic Decelerator (HIAD) technology. This technology development effort is funded by NASA's Space Technology Mission Directorate (STMD) Game Changing Development Program (GCDP). This paper provides an overview of a multi-year HIAD technology development effort, detailing the projects completed to date and the additional testing planned for the future. The effort was divided into three areas: Flexible Systems Development (FSD), Mission Advanced Entry Concepts (AEC), and Flight Validation. FSD consists of a Flexible Thermal Protection Systems (FTPS) element, which is investigating high temperature materials, coatings, and additives for use in the bladder, insulator, and heat shield layers; and an Inflatable Structures (IS) element which includes manufacture and testing (laboratory and wind tunnel) of inflatable structures and their associated structural elements. AEC consists of the Mission Applications element developing concepts (including payload interfaces) for missions at multiple destinations for the purpose of demonstrating the benefits and need for the HIAD technology as well as the Next Generation Subsystems element. Ground test development has been pursued in parallel with the Flight Validation IRVE-3 flight test. A larger scale (6m diameter) HIAD inflatable structure was constructed and aerodynamically tested in the National Full-scale Aerodynamics Complex (NFAC) 40ft by 80ft test section along with a duplicate of the IRVE-3 3m article. Both the 6m and 3m articles were tested with instrumented aerodynamic covers which incorporated an array of pressure taps to capture surface pressure distribution to validate Computational Fluid Dynamics (CFD) model predictions of surface pressure distribution. The 3m article also had a duplicate IRVE-3 Thermal Protection System (TPS) to test in addition to testing with the Aerocover configuration. Both the Aerocovers and the TPS were populated with high contrast targets so that photogrammetric solutions of the loaded surface could be created. These solutions both refined the aerodynamic shape for CFD modeling and provided a deformed shape to validate structural Finite Element Analysis (FEA) models. Extensive aerothermal testing has been performed on the TPS candidates. This testing has been conducted in several facilities across the country. The majority of the testing has been conducted in the Boeing Large Core Arc Tunnel (LCAT). HIAD is continuing to mature testing methodology in this facility and is developing new test sample fixtures and control methodologies to improve understanding and quality of the environments to which the samples are subjected. Additional testing has been and continues to be performed in the NASA LaRC 8ft High Temperature Tunnel, where samples up to 2ft by 2ft are being tested over representative underlying structures incorporating construction features such as sewn seams and through-thickness quilting. With the successful completion to the IRVE-3 flight demonstration, mission planning efforts are ramping up on the development of the HIAD Earth Atmospheric Reenty Test (HEART) which will demonstrate a relevant scale vehicle in relevant environments via a large-scale aeroshell (approximately 8.5m) entering at orbital velocity (approximately 7km/sec) with an entry mass on the order of 4MT. Also, the Build to Print (BTP) hardware built as a risk mitigation for the IRVE-3 project to have a "spare" ready to go in the event of a launch vehicle delivery failure is now available for an additional sub-orbital flight experiment. Mission planning is underway to define a mission that can utilize this existing hardware and help the HIAD project further mature this technology.

Hughes, Stephen J.↗