Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Fluid Dynamics Conduction”

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 325 records · Page 18

Experimental Investigation of Subsonic Turbulent Boundary Layer Flow Over a Wall-Mounted Axisymmetric Hill

An important goal for modern fluid mechanics experiments is to provide datasets which present a challenge for Computational Fluid Dynamics simulations to reproduce. Such "CFD validation experiments" should be well-characterized and well-documented, and should investigate flows which are difficult for CFD to calculate. It is also often convenient for the experiment to be challenging for CFD in some aspects while simple in others. This report is part of the continuing documentation of a series of experiments conducted to characterize the flow around an axisymmetric, modified-cosine-shaped, wall-mounted hill named "FAITH" (Fundamental Aero Investigates The Hill). Computation of this flow is easy in some ways - subsonic flow over a simple shape - while being complex in others - separated flow and boundary layer interactions. The primary set of experiments were performed on a 15.2 cm high, 45.7 cm base diameter machined aluminum model that was tested at mean speeds of 50 m/s (Reynolds Number based on height = 500,000). The ratio of model height to boundary later height was approximately 3. The flow was characterized using surface oil flow visualization, Cobra probe to determine point-wise steady and unsteady 3D velocities, Particle Image Velocimetry (PIV) to determine 3D velocities and turbulence statistics along specified planes, Pressure Sensitive Paint (PSP) to determine mean surface pressures, and Fringe Imaging Skin Friction (FISF) to determine surface skin friction magnitude and direction. A set of pathfinder experiments were also performed in a water channel on a smaller scale (5.1 cm high, 15.2 cm base diameter) sintered nylon model. The water channel test was conducted at a mean test section speed of 3 cm/s (Reynolds Number of 1500), but at the same ratio of model height to boundary layer thickness. Dye injection from both the model and an upstream rake was used to visualize the flow. This report summarizes the experimental set-up, techniques used, and data acquired. It also describes some details of the dataset that is being constructed for use by other researchers, especially the CFD community.

fluid machanics cfd validation wind tunnel↗

Investigation of Physical Mechanisms for Jet Noise Reduction by Plug Nozzle Porosity

This paper presents the results of a computational study on porous-plug nozzles for their eventual application in reducing takeoff noise levels for supersonic civil transport. The study explores a wide range of porous-plug geometries, encompassing variations in hole size, plug length and porosity. Computational fluid dynamics (CFD) analyses are employed in concert with an advanced automated mesh refinement (AMR) scheme, Sketch-to-Solution (S2S), to investigate porous plug nozzle design strategies that may diminish shock strength, broad-band shock noise (BBSN), flow separation, and other undesirable characteristics associated with high takeoff noise levels. A selected subset of the porous-plug nozzles that were analyzed computationally were manufactured at subscale using stereolithography 3D printing. These nozzles were tested experimentally, obtaining far-field noise spectral data and schlieren flow visualization pictures. Detailed comparisons between CFD predictions and experimental results were conducted. These comparisons focus on gaining a better understanding of the underlying physical mechanisms responsible for noise reduction with the porous plugs relative to corresponding solid plugs.

computational fluid dynamics↗

Prototype energy models for data centers

Data centers in the United States consume about two percent of the nation’s electricity. Because heat gains from IT equipment drive cooling demand, data centers offer unique opportunities for energy savings. However, no prototype energy model for data centers is available in the suite of existing U.S. Department of Energy’s Commercial Prototype Building Models. Here we present the development of two new data center prototype models and their implementation in OpenStudio and EnergyPlus. The small-size data center model represents a computer room in a building served by computer room air conditioners (CRACs); while the large-sized model represents stand-alone data centers served by computer room air handlers (CRAHs) with a central chiller plant. For each data center model, two levels of IT equipment (ITE) load density were considered, to cover the wide range of IT power density of data centers: 40 and 100 W/ft 2 (430 and 1076 W/m 2 ) for the computer room, and 100 and 500 W/ft 2 (1076 and 5382 W/m 2 ) for the stand-alone data center. All other assumptions, such as building envelope, lighting, HVAC efficiencies and schedules, were based on the minimal requirements of ASHRAE Standard 90.1 at various vintages. We introduced a novel concept of supply and return air approach temperatures to capture the essential effects of non-uniform airflow and temperature distribution in data centers. The approach temperatures were pre-computed by computational fluid dynamics (CFD) simulations for various configurations of ITE loads and airflow containment management in data centers. A new feature was developed in EnergyPlus to implement the approach temperature method. A case study was conducted to demonstrate the use of the data center models. The two data center models cover all U.S. climate zones and can be used to evaluate energy saving measures for data centers, as well as to support development of data center energy efficiency codes and standards.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

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)↗

Effect of including variable gas properties and entrained particles in the flow analysis of the ASRM nozzle

The objective of this study is to determine the effect of using constant gas properties for the analysis of the Advanced Solid Rocket Motor (ASRM) nozzle and to gain more understanding concerning those types of analysis which might require this additional complexity. Kinetics data for viscosity, thermal conductivity, specific heat, and specific heat ratio are extracted from the Solid Propellant Rocket Motor Performance Prediction Computer Program (SPP) and tabulated as a function of temperature. These tables are added to the Aerovisc computational fluid dynamics (CFD) code in place of the constant gas property values. The results of a CFD analysis of the ASRM 48 inch motor with constant gas properties will be compared with an analysis which uses variable gas properties. Mach number, surface pressure, and torque plots are presented. A full scale ASRM analysis using SPP with and without particle flow is presented.

Clayton, Curtis D.↗

Prediction of forces and moments for hypersonic flight vehicle control effectors

This research project includes three distinct phases. For completeness, all three phases of the work are briefly described in this report. The goal was to develop methods of predicting flight control forces and moments for hypersonic vehicles which could be used in a preliminary design environment. The first phase included a preliminary assessment of subsonic/supersonic panel methods and hypersonic local flow inclination methods for such predictions. While these findings clearly indicated the usefulness of such methods for conceptual design activities, deficiencies exist in some areas. Thus, a second phase of research was conducted in which a better understanding was sought for the reasons behind the successes and failures of the methods considered, particularly for the cases at hypersonic Mach numbers. This second phase involved using computational fluid dynamics methods to examine the flow fields in detail. Through these detailed predictions, the deficiencies in the simple surface inclination methods were determined. In the third phase of this work, an improvement to the surface inclination methods was developed. This used a novel method for including viscous effects by modifying the geometry to include the viscous/shock layer.

Maughmer, Mark D.↗

Improved ceramic slip casting technique

A primary concern in modern fluid dynamics research is the experimental verification of computational aerothermodynamic codes. This research requires high precision and detail in the test model employed. Ceramic materials are used for these models because of their low heat conductivity and their survivability at high temperatures. To fabricate such models, slip casting techniques were developed to provide net-form, precision casting capability for high-purity ceramic materials in aqueous solutions. In previous slip casting techniques, block, or flask molds made of plaster-of-paris were used to draw liquid from the slip material. Upon setting, parts were removed from the flask mold and cured in a kiln at high temperatures. Casting detail was usually limited with this technique -- detailed parts were frequently damaged upon separation from the flask mold, as the molded parts are extremely delicate in the uncured state, and the flask mold is inflexible. Ceramic surfaces were also marred by 'parting lines' caused by mold separation. This adversely affected the aerodynamic surface quality of the model as well. (Parting lines are invariably necessary on or near the leading edges of wings, nosetips, and fins for mold separation. These areas are also critical for flow boundary layer control.) Parting agents used in the casting process also affected surface quality. These agents eventually soaked into the mold, the model, or flaked off when releasing the case model. Different materials were tried, such as oils, paraffin, and even an algae. The algae released best, but some of it remained on the model and imparted an uneven texture and discoloration on the model surface when cured. According to the present invention, a wax pattern for a shell mold is provided, and an aqueous mixture of a calcium sulfate-bonded investment material is applied as a coating to the wax pattern. The coated wax pattern is then dried, followed by curing to vaporize the wax pattern and leave a shell mold of the calcium sulfate-bonded investment material. The shell mold is cooled to room temperature, and a ceramic slip is poured therein. After a ceramic shell of desired thickness has set up in the shell mold, excess ceramic slip is poured out. While still wet, the shell mold is peeled from the ceramic shell to expose any delicate or detailed parts, after which the ceramic shell is cured to provide a complete, detailed, precision ceramic article without parting lines.

Buck, Gregory M.↗

Transonic Investigation of Two-Dimensional Nozzles Designed for Supersonic Cruise

An experimental and computational investigation has been conducted to determine the off-design uninstalled drag characteristics a of a two-dimensional convergent-divergent nozzle designed for a supersonic cruise civil transport. The main objective of this investigation was to determine the effects of varying nozzle external flap curvature and sidewall boattail angle and curvature on nozzle drag The experimental investigation was conducted in the Langley 16-Foot Transonic Tunnel at Mach numbers from 0.80 to 1.20 at nozzle pressure ratios up to nine. Three-dimensional simulations of nozzle performance were obtained with the computational fluid dynamics code PAB using turbulence closure and nonlinear Reynolds stress modeling. The results of this investigation indicate that excellent correlation between experimental and predicted results was obtained for the nozzle with a moderate amount of boattail curvature. The nozzle with an external flap having a sharp shoulder (no curvature) had the lowest nozzle pressure drag. At a Mach number of 1.2, sidewall pressure drag doubled as sidewall boattail angle was increased from 4 to 8 deg. Reducing the height of the sidewall caused large decreases in both the sidewall and flap pressure drags.

Capone, Francis J.↗

A Computational Study of the Abrupt Wing Stall (AWS) Characteristics for Various Fighter Jets: F/A-18E and F-16C - Part 1

Steady state Computational Fluid Dynamic (CFD) simulations are used to gain an understanding of the physics behind the Abrupt Wing Stall (AWS) phenomenon and to arrive at static Figures of Merit (FOMs). Navier-Stokes simulations are conducted using the NASA Langley developed TetrUSS simulation suite which is based on tetrahedral, unstructured grids. The physics of the AWS phenomenon is understood by comparing CFD simulation results on two aircraft; a pre-production F/A-18E configuration which exhibits AWS phenomenon under certain geometric and flow conditions, and an F-16C aircraft configuration that does not. The CFD code is first validated against two sets of experimental data to build confidence in its use for the problem of AWS. An attempt is then made to understand the possible causes of AWS by analyzing and comparing the detailed flow fields between the two configurations under a variety of flow conditions. Based on this approach, a number of static Figures of Merit are developed to predict the potential existence of AWS. The FOMs include the break in the lift and wing root bending moment versus angle-of-attack (AoA) curves and the rate of change of sectional lift with respect to AoA. A companion paper, in Part II, describes a similar CFD study on two other aircraft, the AV-8B Harrier and F/A-18C. Results from both these studies, as well as other CFD studies conducted as part of the AWS program are used to recommend a CFD procedure for predicting the existence of AWS in future aerospace designs.

Parikh, Paresh↗

Slosh-Induced Ullage Collapse Validation by Computational Fluid Dynamics and Application to In-space Propellant Dynamics

Ullage collapse predictions made using Computational Fluid Dynamics (CFD) were validated then made for an in-space cryogenic propellant unsettling application. Ullage collapse is the decrease in propellant tank pressure due to cooling of the ullage gas, which can be significant for cryogenic propulsion systems. Vehicles may be developed that are robust to temperature and pressure changes caused by heat leaks, e.g., propellant tanks with high maximum operation pressures, technologies yielding active thermal conditioning, and systems with ample pressurant for conditioning. These methods are robust because additional mass, cost, and complexity are used to address cryogenic fluid management (CFM). As vehicles are optimized for increased payload capacity, tracking of events such as ullage collapse and self-pressurization in tanks will become increasingly important for propellant inventory calculations. Some ullage collapse events can be studied through ground experimentation, but propellant unsettling in and around in-space coast mission phases must be studied in the absence of an acceleration field. Computational Fluid Dynamics is a powerful tool to address this need, enabling prediction of complex physics prior to flight. Validation of the chosen CFD tool and modeling methodology was done using data from a previously conducted subscale slosh-induced ullage collapse ground test performed in part by the Japan Aerospace Exploration Agency (JAXA). A small tank was partially filled with liquid nitrogen and pressurized with gaseous nitrogen at various temperatures. The tank was oscillated laterally to induce high amplitude slosh dynamics resulting in wave breakup. Successful qualitative validation of slosh wave dynamics and quantitative validation of ullage collapse promoted confidence in extending ullage collapse predictions to a flight application. Mechanisms of ullage collapse caused by main engine cutoff while in-space were identified. Both slosh dynamics leading to engine shutdown and body dynamics following shutdown were shown to contribute to propellant unsettling and ullage collapse.

J M Brodnick↗

Investigation of Multiphase Flow in a Packed Bed Reactor Under Microgravity Conditions

In this paper we study the two-phase flow phenomena in a packed bed reactor using an integrated experimental and numerical method. The cylindrical bed is filled with uniformly sized spheres. In the experiment water and air are injected into the bed simultaneously. The pressure distribution along the bed will be measured. The numerical simulation is based on a two-phase flow solver which solves the Navier-Stokes equations on Cartesian grids. A novel coupled level set and moment of fluid method is used to construct the interface. A sequential method is used to position spheres in the cylinder. Preliminary experimental results showed that the tested flow rates resulted in pulse flow. The numerical simulation revealed that air bubbles could merge into larger bubbles and also could break up into smaller bubbles to pass through the pores in the bed. Preliminary results showed that flow passed through regions where the porosity is high. Comparison between the experimental and numerical results in terms of pressure distributions at different flow injection rates will be conducted. Comparison of flow phenomena under terrestrial gravity and microgravity will be made.

computational fluid dynamics↗

Simulation of Channel Flow with Square Ribs for Blanket First-Wall Cooling: Geometry-Specific Tuning of k-ω Model Using Adjoint Method

Cooling of the plasma-facing first wall is challenging in the design of blanket components because of the high heat flux (on the order of 𝑀𝑊/𝑚2) from the plasma, especially when a low thermal mass medium like helium is chosen as the coolant. Therefore, heat transfer enhancement in which the convective heat transfer rate is augmented by the addition of turbulence-promoting structures becomes a key initiative for providing sufficient cooling capability with helium. Previously, computational fluid dynamics simulations had been performed on pipe flows with different transverse and longitudinal ribbed geometries at Oak Ridge National Laboratory to compare the enhancement performance among different ribbed geometries. Rib shape morphing had been conducted to obtain an optimized rib profile. In the work presented here, the adjoint method is adopted in the ANSYS Fluent solver for turbulence model augmentation, and the Generalized k-ω (GEKO) turbulence model is employed because of its ability of tuning the turbulence model. The Nusselt number and pressure drop obtained from the channel flow with bottom ribbed wall experiments are used as the targets. Sensitivity analysis provides information as guidance to improve the turbulence model accuracy. The augmented GEKO model is tuned for the studied ribbed channel geometry and flow conditions, providing improved predictive accuracy within this context. Extension to other configurations offers potential but may require additional tuning and validation.

Xu, Tracy [ORNL] (ORCID:0009000193700887)↗

Computational Fluid Dynamics Modeling of Solar Thermal Dry Reforming of Methane in a Parabolic Trough

Computational fluid dynamics simulations of solar-thermal dry reforming of methane using a parabolic trough configuration were performed. Parametric simulations of different combinations of gas flow rate, receiver tube emissivity, and geometric concentration ratio were conducted to determine configurations that could achieve the required catalyst temperatures of at least 700 °C to achieve high conversion of CH4 and CO2 to H2 and CO. Results showed that the concentration ratio of the parabolic trough collector had to be increased from ~70 to ~120 and the receiver-tube emissivity had to be reduced to ~0.2 to achieve bulk average catalyst temperatures of greater than 700 °C. Lower gas flow rates also reduced enthalpic heat losses and increased catalyst temperatures.

Ho, Clifford↗

Finite element algorithms for compressible flow computation on a supercomputer

Methods of applying computational fluid dynamics techniques to engineering problems are discussed. The linking factors between theoretical and applied research are the shape and weight functions which produce the Galerkin finite element schemes. The flow is unsteady, viscous, heat conducting and compressible, and steady-state flows are the asymptotic limit of unsteady flows. All flows are described by the time-averaged Navier-Stokes equations (NSE) with mass and energy conservation. Local curvilinear intrinsic coordinates are applied to discretize the NSE in arbitrary geometric domains. A modified weighted residuals approach defines the discrete analogs of the physical systems modeled and a general interpolants method is used to derive families of numerical models, both implicit and explicit, finite difference and finite elemnt, from a single point of departure. Time-dependent element approximations are achieved with a progressive assembly of generalized elements method. Hyperbolic steady-state Euler algorithms and quasi-parabolic are spatial marching algorithms for solving the discretized equations. Application of the techniques is illustrated through calculation of the compressible flow around a hypersonic flight vehicle at 100,000 ft at zero angle of attack at speeds of Mach 4-10.

Spradley, L. W.↗

High-Fidelity Multidisciplinary Design Optimization Methodology with Application to Rotor Blades

A multidisciplinary design optimization procedure has been developed and applied to rotorcraft simulations involving tightly-coupled, high-fidelity computational fluid dynamics and comprehensive analysis. A discretely-consistent, adjoint-based sensitivity analysis available in the fluid dynamics solver provides sensitivities arising from unsteady turbulent flows on unstructured, dynamic, overset meshes, while a complex-variable approach is used to compute structural sensitivities with respect to aerodynamic loads. The multidisciplinary sensitivity analysis is conducted through integrating the sensitivity components from each discipline of the coupled system. Accuracy of the coupled system for high-fidelity rotorcraft analysis is verified; simulation results exhibit good agreement with established solutions. A constrained gradient-based design optimization for a HART-II rotorcraft configuration is demonstrated. The computational cost for individual components of the multidisciplinary sensitivity analysis is assessed and improved.

Wang, Li↗

High-fidelity Multidisciplinary Sensitivity Analysis and Design Optimization for Rotorcraft Applications

A multidisciplinary sensitivity analysis of rotorcraft simulations involving tightly coupled high-fidelity computational fluid dynamics and comprehensive analysis solvers is presented and evaluated. A sensitivity-enabled fluid dynamics solver and a nonlinear flexible multibody dynamics solver are coupled to predict aerodynamic loads and structural responses of helicopter rotor blades. A discretely consistent adjoint-based sensitivity analysis available in the fluid dynamics solver provides sensitivities arising from unsteady turbulent flows and unstructured dynamic overset meshes, while a complex-variable approach is used to compute structural sensitivities with respect to aerodynamic loads. The multidisciplinary sensitivity analysis is conducted through integrating the sensitivity components from each discipline of the coupled system. Accuracy of the coupled system is validated by conducting simulations for a benchmark rotorcraft model and comparing solutions with established analyses and experimental data. Sensitivities of lift computed by the multidisciplinary sensitivity analysis are verified by comparison with the sensitivities obtained by complex-variable simulations. Finally the multidisciplinary sensitivity analysis is applied to a constrained gradient-based design optimization for a HART-II rotorcraft configuration.

Wang, Li↗

The use of computer models to predict temperature and smoke movement in high bay spaces

The Building and Fire Research Laboratory (BFRL) was given the opportunity to make measurements during fire calibration tests of the heat detection system in an aircraft hangar with a nominal 30.4 (100 ft) ceiling height near Dallas, TX. Fire gas temperatures resulting from an approximately 8250 kW isopropyl alcohol pool fire were measured above the fire and along the ceiling. The results of the experiments were then compared to predictions from the computer fire models DETACT-QS, FPETOOL, and LAVENT. In section A of the analysis conducted, DETACT-QS AND FPETOOL significantly underpredicted the gas temperature. LAVENT at the position below the ceiling corresponding to maximum temperature and velocity provided better agreement with the data. For large spaces, hot gas transport time and an improved fire plume dynamics model should be incorporated into the computer fire model activation routines. A computational fluid dynamics (CFD) model, HARWELL FLOW3D, was then used to model the hot gas movement in the space. Reasonable agreement was found between the temperatures predicted from the CFD calculations and the temperatures measured in the aircraft hangar. In section B, an existing NASA high bay space was modeled using the CFD model. The NASA space was a clean room, 27.4 m (90 ft) high with forced horizontal laminar flow. The purpose of this analysis is to determine how the existing fire detection devices would respond to various size fires in the space. The analysis was conducted for 32 MW, 400 kW, and 40 kW fires.

Notarianni, Kathy A.↗

Scalability of Localized Arc Filament Plasma Actuators

Temporal flow control of a jet has been widely studied in the past to enhance jet mixing or reduce jet noise. Most of this research, however, has been done using small diameter low Reynolds number jets that often have little resemblance to the much larger jets common in real world applications because the flow actuators available lacked either the power or bandwidth to sufficiently impact these larger higher energy jets. The Localized Arc Filament Plasma Actuators (LAFPA), developed at the Ohio State University (OSU), have demonstrated the ability to impact a small high speed jet in experiments conducted at OSU and the power to perturb a larger high Reynolds number jet in experiments conducted at the NASA Glenn Research Center. However, the response measured in the large-scale experiments was significantly reduced for the same number of actuators compared to the jet response found in the small-scale experiments. A computational study has been initiated to simulate the LAFPA system with additional actuators on a large-scale jet to determine the number of actuators required to achieve the same desired response for a given jet diameter. Central to this computational study is a model for the LAFPA that both accurately represents the physics of the actuator and can be implemented into a computational fluid dynamics solver. One possible model, based on pressure waves created by the rapid localized heating that occurs at the actuator, is investigated using simplified axisymmetric simulations. The results of these simulations will be used to determine the validity of the model before more realistic and time consuming three-dimensional simulations are conducted to ultimately determine the scalability of the LAFPA system.

Brown, Clifford A.↗