Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational fluid dynamics simulation”

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 145 records · Page 8

CFD Simulations of the IHF 13-Inch Nozzle Flow: 55-deg Sphere-Cone Model, Manufactured Fences and Gaps

This paper reports computational analyses of testsin a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using 55-deg sphere-cone models placed in a free jet downstream of the 13-inch diameter conical nozzle in the Ames 60-MW Interaction Heating Facility. Some of the sphere-cone models include surface features such as manufactured fences and gaps intended to simulate effects of differential recession, all of which disturb the flow, producing augmented heating locally and downstream.Test calibration data were obtained using slug and Gardon gage stagnation calorimeters, and a sphere-cone calorimeter model with six Gardon gages and five pressure tabs.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 simulations take into account nonuniform total enthalpy and mass flux profiles at the nozzle inlet as well as the expansion waves emanating from the nozzle exit and their effects on the model flowfields.

Arc-jets↗

A machine learning-based interaction force model for non-spherical and irregular particles in low Reynolds number incompressible flows

In this study, interaction force of non-spherical particles in low Reynolds number gas-solid flow is investigated by neural network approaches. An artificial neural network (ANN) model is developed to correlate the non-spherical particle shape and the flow conditions with the interaction force. To define the particle shape, spherical harmonic expansion is applied. Furthermore, variational autoencoder model is then used to extract latent geometric features. The latent vector is utilized as an input with the Reynolds number for the ANN. The interaction force data, which is used as output data of the ANN, is obtained by particle resolved direct numerical simulation for 5200 non-spherical particles. The proposed model enables unsupervised extraction for non-spherical particle shapes and accurate predictions on the interaction force without heavy computation. This study provides the model that can explain complicated shapes of particles and be applied to a large scale, computational fluid dynamics simulation.

01 COAL, LIGNITE, AND PEAT↗

Pressure Response Optimization of an Eddy Current-Driven Flyer Plate Valve for the ITER Shattered Pellet Injection System

One technique for mitigating disruptions in a tokamak is shattered pellet injection (SPI). SPI is a process in which a large solid pellet consisting of deuterium, neon, or argon is desublimated in a pipe gun barsrel and launched downstream. Pellets are shattered just before entering the plasma by an impact with an angled tube. Injection of these materials into the plasma radiates stored thermal energy, limits current decay rates, suppresses the generation of runaway electrons, and dissipates runaway electrons if necessary. A critical element of the SPI system is a fast-acting valve that releases high-pressure gas to dislodge and accelerate pellets directly, or indirectly via a mechanical punch. A prototype valve sized for the ITER SPI system has been designed and fabricated. A pulsed high-voltage power supply energizes the valve’s internal magnetic coil, which induces eddy currents in the adjacent flyer plate resulting in a repulsive force between the flyer plate and the coil. The flyer plate action lifts a valve seat, allowing high-pressure gas to flow from the valve plenum to the downstream (breech) location of the pellet or mechanical punch. All of the valve’s internal components are designed to operate in ITER-level static background magnetic fields. Here, a study was conducted to optimize the downstream pressure response for a range of valve sizes and operating pressures. In particular, the study analyzes the breech pressure response associated with varying plenum pressures as well as varying breech volumes. A computational fluid dynamics simulation was built in STAR-CCM+ and validated against data from laboratory experiments. The resulting simulation outputs, in the form of downstream responses for a variety of initial plenum pressures and breech volumes, will be used as a complement to experimental data to ensure the pressure pulse is suitable for pellet survivability. These data, combined with studies on pellet shear strength and shock response, will be applied to optimization of overall operating parameters of the ITER SPI system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CFD Simulations of the IHF 13-Inch Nozzle Flow: 55-deg Sphere-Cone Model, Manufactured Fences and Gaps

This paper reports computational analyses of tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using 55-deg sphere-cone models placed in a free jet downstream of the 13-inch diameter conical nozzle in the Ames 60-MW Interaction Heating Facility. Some of the sphere-cone models include surface features such as manufactured fences and gaps intended to simulate effects of differential recession, all of which disturb the flow, producing augmented heating locally and downstream. Test calibration data were obtained using slug and Gardon gage stagnation calorimeters, and a sphere-cone calorimeter model with six Gardon gages and five pressure tabs. 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 simulations take into account nonuniform total enthalpy and mass flux profiles at the nozzle inlet as well as the expansion waves emanating from the nozzle exit and their effects on the model flowfields.

Arc-jets↗

CFD Simulations of the IHF 13-Inch Nozzle Flow: 55-deg Sphere-Cone Model, Manufactured Fences and Gaps

This paper reports computational analyses of tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using 55-deg sphere-cone models placed in a free jet downstream of the 13-inch diameter conical nozzle in the Ames 60-MW Interaction Heating Facility. Some of the sphere-cone models include surface features such as manufactured fences and gaps intended to simulate effects of differential recession, all of which disturb the flow, producing augmented heating locally and downstream. Test calibration data were obtained using slug and Gardon gage stagnation calorimeters, and a sphere-cone calorimeter model with six Gardon gages and five pressure tabs. 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 simulations take into account nonuniform total enthalpy and mass flux profiles at the nozzle inlet as well as the expansion waves emanating from the nozzle exit and their effects on the model flowfields.

Arc-jets↗

On Laminar to Turbulent Transition of Arc-Jet Flow in the NASA Ames Panel Test Facility

This paper provides experimental evidence and supporting computational analysis to characterize the laminar to turbulent flow transition in a high enthalpy arc-jet facility at NASA Ames Research Center. The arc-jet test data obtained in the 20 MW Panel Test Facility include measurements of surface pressure and heat flux on a water-cooled calibration plate, and measurements of surface temperature on a reaction-cured glass coated tile plate. Computational fluid dynamics simulations are performed to characterize the arc-jet test environment and estimate its parameters consistent with the facility and calibration measurements. The present analysis comprises simulations of the nonequilibrium flowfield in the facility nozzle, test box, and flowfield over test articles. Both laminar and turbulent simulations are performed, and the computed results are compared with the experimental measurements, including Stanton number dependence on Reynolds number. Comparisons of computed and measured surface heat fluxes (and temperatures), along with the accompanying analysis, confirm that that the boundary layer in the Panel Test Facility flow is transitional at certain archeater conditions.

Gokcen, Tahir↗

CFD Simulations of the IHF 13-Inch Nozzle Flow: 55° Sphere-Cone Model, Manufactured Fences and Gaps

This paper reports computational analyses of tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using 23.6-cm diameter, 55° sphere-cone models placed in a free jet downstream of the 33.0-cm diameter conical nozzle in the 60-MW Interaction Heating Facility. Some of the sphere-cone models include surface features such as manufactured fences and gaps intended to simulate effects of differential recession, all of which disturb the flow, producing augmented heating locally and downstream. Test calibration data were obtained using slug and Gardon gage stagnation calorimeters, and a sphere-cone calorimeter model with six Gardon heat flux gages and five pressure gages. 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. Various issues related to testing are considered: diffuser flow capture, flow characterization based on the calorimeter data, effects of the model surface recession, and prediction of surface quantities for the models with surface features. These simulations take into account nonuniform total enthalpy profiles at the nozzle inlet as well as the expansion waves emanating from the nozzle exit and their effects on the model flowfields.

Arc-jets↗

Detailed Simulations of a Three-Stage Supercritical Carbon Dioxide Axial Compressor with a Focus on the Shrouded Stator Cavities Flow

This paper describes the findings of detailed simulations performed to investigate the impact of seal teeth cavity leakage flow on the aerodynamic and thermal performance of a three-stage supercritical CO 2 axial compressor. The study compares a shrouded stator configuration (with cavities) to a cantilevered stator configuration (without cavities) to highlight their differences. High-fidelity computational fluid dynamics simulations were performed using non-linear harmonic (NLH) and mixing plane assumptions, considering various possible rotor/stator interface configurations for mixing plane calculations. The key performance parameters for each case were compared, and the best-performing configuration selected for further analysis. The individual stage performance parameters are also examined and compared between the cantilevered and shrouded configurations. It was observed that in the shrouded case, the leakage flow enters the cavity downstream of the stator trailing edge and gets entrained into the primary flow upstream of the stator, leading to boundary layer changes at the hub and degradation of stator and downstream rotor performance. Vortical flow structures were also observed in the stator wells, which tended to change the flow angles around the region, thereby affecting mixing and velocity distribution, resulting in a slight deterioration of compressor performance. Additionally, the study examines windage heating due to shear work from rotating walls, including the seal teeth surface. The amount of shear work done on the leakage flow and the corresponding rise in fluid temperature were quantified, tabulated, and further compared with a simple analytical model, showing good agreement between them and, hence, validating the numerical approach used.

42 ENGINEERING↗

Numerical simulation of frost formation and heat transfer on fin-and-tube heat exchangers in turbulent cross-flow

Frost formation in fin-and-tube heat exchangers in turbulent cross-flow presents significant challenges in industrial refrigeration applications, affecting heat transfer efficiency and operational reliability. The purpose of this work is to investigate frost deposition and growth on a staggered bank of a fin-and-tube freezer coil under turbulent forced convection conditions. The focus here is on investigating conditions that closely replicate real-world scenarios in large walk-in industrial freezers. Using a direct numerical simulation approach, we examine the flow dynamics and thermal behaviour in the presence of frost, considering turbulent regimes characterized by a Reynolds number in the range 1050 ≤ R e D , avg ≤ 4800 , with the characteristic length being the outer diameter of the tube and the velocity being the bulk fluid velocity between the plates (fins). Computational fluid dynamics simulations are employed to resolve the interactions between turbulent airflow and the frost layer. Our approach incorporates a modified immersed boundary method and a slow-time acceleration technique to address the complex dynamic interface between the continuously evolving frost layer and the flowing air stream. Our findings indicate that frost forms more on the sides of the finned surfaces (plates) and less on the tubes themselves. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗

High-Fidelity Simulations of a Tiltwing Vehicle for Urban Air Mobility

NASA is playing an important role in the development of the urban air mobility ecosystem by identifying key research areas and establishing a fleet of conceptual reference vehicles for urban air mobility applications. This paper delves deeper into the study of NASA’s six-passenger tiltwing air taxi aerodynamics using high-fidelity computational fluid dynamics simulations. NASA’s OVERFLOW flow solver has been loosely-coupled with the rotorcraft comprehensive analysis code CAMRAD II. Dual-time stepping, high-order accurate schemes, and the delayed detached-eddy simulation are employed. The tiltwing configuration is analyzed in airplane and helicopter mode, with flight conditions representative of an air taxi mission. Airloads, wake geometry, and performance are studied to identify the unique flow features of this configuration. NASA’s tiltwing vehicle is one of the conceptual designs intended to focus and guide NASA’s research activities in support of aircraft development for vertical take-off and landing air taxi operations.

RVLT↗

Investigation of flow regimes in arc plasma–gas interactions using a two-temperature arc in crossflow model

The perpendicular impingement of a gas stream on an electric arc, a configuration known as the arc in crossflow, is of primary relevance in the study of plasma-gas interactions, as well as in industrial applications such as circuit breakers and wire-arc spraying. The flow dynamics in the arc in crossflow are the result of coupled fluid-thermal-electromagnetic phenomena accompanied by large property gradients, which can produce significant deviations from Local Thermodynamic Equilibrium (LTE) among electrons and gas species. These characteristics can lead to the establishment of distinct flow regimes depending on the relative values of the controlling parameters of the system, such as inflow velocity, arc current, and inter-electrode spacing. A two-temperature non-LTE (NLTE) model is used to investigate the arc dynamics and the establishment of flow regimes in the arc in crossflow. The plasma flow model is implemented within a nonlinear Variational Multiscale (VMS) numerical discretization approach that is less dissipative, and hence better suited to capture unstable behavior, than traditional VMS methods commonly used in computational fluid dynamics simulations. The Reynolds and the Enthalpy dimensionless numbers, characterizing the relative flow strength and arc strength, respectively, are chosen as the controlling parameters of the system. Simulation results reveal the onset of dynamic behavior and the establishment of steady, periodic, quasi-periodic, and chaotic or potentially turbulent, regimes, as identified by distinct spatiotemporal fluctuations. The computational findings reveal the role of increasing the relative arc strength on enhancing flow stability by delaying the growth of fluctuating and unstable flow behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Computational Study of Plume Modeling For Space Launch System Abort Scenarios

This extended abstract details work on Space Launch System (SLS) core stage drag during abort scenarios that is currently being conducted. In the final paper, viscous computational fluid dynamic simulations will be used to study the axial force coefficient on the SLS core stage as a function of its distance to the upstream crew module and launch abort system, the freestream conditions, and the chemical model used to represent the launch abort system abort motor and attitude control motor plumes. In addition, the computational expense of each method will be studied in order to help determine the appropriate trade-off between fidelity and simulation turn around time.

CFD↗

On Characterization of Laminar to Turbulent Transitional Flow in the Truncated IHF Semi-Elliptical Nozzle

This paper reports computational analysis and experimental data to characterize the laminar to turbulent flow transition in 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 the flowfields over panel test articles, and comparisons with the measured test data. Both laminar and turbulent simulations are performed. Constant current studies are performed: the measured data are grouped at the same arc current levels with varying mass flow rates, and they are compared with corresponding computations. For all currents studied, these studies indicate that the semi-elliptical nozzle flow is transitional at relatively higher mass flow rates within the facility operating envelope.

Arc-jets↗

On Characterization of Laminar to Turbulent Transitional Flow in the Truncated IHF Semi-Elliptical Nozzle

This paper reports computational analysis and experimental data to characterize the laminar to turbulent flow transition in 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 the flowfields over panel test articles, and comparisons with the measured test data. Both laminar and turbulent simulations are performed. Constant current studies are conducted: the measured data are grouped at the same arc current levels with varying mass flow rates, and they are compared with corresponding computations. For all currents considered, these studies indicate that the semi-elliptical nozzle flow is transitional at relatively higher mass flow rates within the facility operating envelope. An approximate facility map for the transitional flow characterization is developed, showing laminar and transitional to turbulent flow regions.

Arc-jets↗

A phenomenological rate of injection model for predicting fuel injection with application to mixture formation in light-duty diesel engines

Fuel injection rate laws are one of the most important pieces of information needed when modeling engine combustion with computational fluid dynamics. In this study, a simple phenomenological model of a common-rail injector was developed and calibrated for the Bosch CRI2.2 platform. The model requires three tunable parameter fits, making it relatively easy to calibrate and suitable for injector modeling when high-fidelity information about the internal injector’s geometry and electrical circuit details are not available. Each injection pulse is modeled as a sequence of up to four stages: an injection needle mechanical opening transient, a full-lift viscous flow inertial transient, a Bernoulli steady-state stage, and a needle descent transient. Parameters for each stage are obtained as polynomial fits from measured injection rate properties. The model enforces total injected mass, and the intermediate stages are only introduced if the injection pulse duration is long enough. Experimental rates of injection from two separate campaigns on the same injector were used to calibrate the model. The model was first validated against measured injection rate laws featuring pilot injections, short partially premixed combustion pulses, and conventional diesel combustion injection strategies. Then, it was employed as an input to engine computational fluid dynamics simulations, which were run to simulate experiments of mixture formation in an optically accessible light-duty diesel engine. It was found that, though simple, this model is capable of predicting both pilot and main injection pulse mass flow rates well: the simulations yielded accurate predictions of in-cylinder equivalence ratio distributions from injection strategies for both partially premixed combustion and pilot injections. Also, once calibrated, the model produced appropriate results for a wide range of injected mass and rail pressure values. Finally, it was observed that usage of such a relatively simple model can be a good choice when high-fidelity injection rate input and highly detailed information of the injector’s geometry and operation are not available, particularly as noticeable discrepancies can be present also among different experimental campaigns on similar hardware.

42 ENGINEERING↗

FLOW Estimation and Rose Superposition (FLOWERS): an integral approach to engineering wake models

Abstract. Annual energy production (AEP) is often the objective function in wind plant layout optimization studies. The conventional method to compute AEP for a wind farm is to first evaluate power production for each discrete wind direction and speed using either computational fluid dynamics simulations or engineering wake models. The AEP is then calculated by weighted-averaging (based on the wind rose at the wind farm site) the power produced across all wind directions and speeds. We propose a novel formulation for time-averaged wake velocity that incorporates an analytical integral of a wake deficit model across every wind direction. This approach computes the average flow field more efficiently, and layout optimization is an obvious application to exploit this benefit. The clear advantage of this new approach is that the layout optimization produces solutions with comparable AEP performance yet is completed 2 orders of magnitude faster. The analytical integral and the use of a Fourier expansion to express the wind speed and wind direction frequency create a relatively smooth solution space for the gradient-based optimizer to excel in comparison to the existing weighted-averaging power calculation.

17 WIND ENERGY↗

Generic Wing-Body Aerodynamics Data Base

The wing-body aerodynamics data base consists of a series of CFD (Computational Fluid Dynamics) simulations about a generic wing body configuration consisting of a ogive-circular-cylinder fuselage and a simple symmetric wing mid-mounted on the fuselage. Solutions have been obtained for Nonlinear Potential (P), Euler (E) and Navier-Stokes (N) solvers over a range of subsonic and transonic Mach numbers and angles of attack. In addition, each solution has been computed on a series of grids, coarse, medium and fine to permit an assessment of grid refinement errors.

Holst, Terry L.↗

Automated and efficient local adaptive regression for principal component-based reduced-order modeling of turbulent reacting flows

Principal Component Analysis can be used to reduce the cost of Computational Fluid Dynamics simulations of turbulent reacting flows by reducing the dimensionality of the transported variables through projection of the thermochemical state onto a lower-dimensional manifold. However, because of the nonlinearity of the principal component source terms, nonlinear regression techniques must be utilized for the source terms in terms of the principal components. Unfortunately, widely available and utilized nonlinear regression techniques can have prohibitive computational requirements and/or accuracy that is highly dependent on user experience in ad hoc tuning of model architecture and hyperparameters. Here, in this work, a new nonlinear regression approach is proposed that is both computationally efficient and automated so does not require any user input. The approach is evaluated through a priori prediction of principal component source terms using data from a Direct Numerical Simulation of a turbulent nonpremixed n-heptane/air jet flame. In particular, the proposed framework consists of local regressions whose complexity is adapted according to the local nonlinearity of the data: local linear regression when accurate enough and local Artificial Neural Networks when nonlinear regression is required. The number of local clusters for local regression is determined automatically using the Davies-Bouldin index. In addition, Bayesian optimization is utilized for model training (i.e., to select the best architectures and hyperparameters of the nonlinear regressions in an unsupervised fashion), eliminating ad hoc hand-tuning and/or expensive grid searches. Overall, compared to a single, global neural network, the new local adaptive regression approach is shown to have comparable accuracy but 69% less training time due to the utilization of local linear regression and faster training of local neural networks.

42 ENGINEERING↗