Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow 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 163 records · Page 9

Application of an Upwind High Resolution Finite-Differencing Scheme and Multigrid Method in Steady-State Incompressible Flow Simulations

The analysis and design of a submarine propulsor requires the ability to predict the characteristics of both laminar and turbulent flows to a higher degree of accuracy. This report presents results of certain benchmark computations based on an upwind, high-resolution, finite-differencing Navier-Stokes solver. The purpose of the computations is to evaluate the ability, the accuracy and the performance of the solver in the simulation of detailed features of viscous flows. Features of interest include flow separation and reattachment, surface pressure and skin friction distributions. Those features are particularly relevant to the propulsor analysis. Test cases with a wide range of Reynolds numbers are selected; therefore, the effects of the convective and the diffusive terms of the solver can be evaluated separately. Test cases include flows over bluff bodies, such as circular cylinders and spheres, at various low Reynolds numbers, flows over a flat plate with and without turbulence effects, and turbulent flows over axisymmetric bodies with and without propulsor effects. Finally, to enhance the iterative solution procedure, a full approximation scheme V-cycle multigrid method is implemented. Preliminary results indicate that the method significantly reduces the computational effort.

NAVIER STOKES EQUATIONS↗

Unsteady Flow Simulations in Support of the SSME HEX Turning Vane Cracking Investigation with the ATD HPOTP

Unsteady flow computations are being performed with the P&W (ATD) and the Rocketdyne baseline configurations of the SSME LO2 turbine turnaround duct (TAD) and heat exchanger (HEX). The work is in support of the HEX inner turning vane cracking investigation. Fatigue cracking has occurred during hot firings with the P&W configuration on the HEX inner vane, and it appears the fix will involve changes to the TAD splitter vane position and to the TAD inner wall curvature to reduce the dynamic loading on the inner vane. Unsteady flow computations on the P&W baseline and fix and on the Rocketdyne baseline reference follow steady-flow screening computations done by MSFC/ED32 on several trial configurations arriving at the fix. The P&W TAD inlet velocity profile has a strong radial velocity component that directs the flow toward the inner wall and raises the local velocity a factor of two and the dynamic pressure a factor, of four. The fix is intended to redistribute the flow more evenly across the HEX inner and outer vanes like the Rocketdyne baseline reference. Vane buffeting at frequencies around 4,000 Hz is the leading suspected cause of the problem. Our simulations (work in progress) are being done with the USA 2D axisymmetric code approximating the flow as axisymmetric u+v 2D (axial, u, and radial, v, components only). The HEX coils are included in the model to make sure the fix does not adversely affect the HEX environment. Turbulent kinetic energy, k, levels where k = 1/2 v' rms2 are locally as high as 10,000 ft2/sec2 for the P&W baseline at the engine interface (between the TAD and HEX) at the HEX inner vane location. However, k is less than 8,000 on the HEX outer vane and only about 4,500 on the HEX inner vane for the Rocketdyne baseline. Unsteady turbulence intensity, v'rms/v, and pressure, p', are being computed in the present computations to compare with steady-flow Reynolds-averaged computations where p'rms = const (pk) for overall rms random turbulence from 0.1 to 12,000 Hz frequency. Random overall static, p'rms fluctuations as large as 1.7 psi are estimated from k on the HEX inner vane for the P&W baseline configuration but only about 0.7 psi for the Rocketdyne configuration.

Dougherty, N. S.↗

Groundwater flow simulation of the Savannah River Site general separations area

The most recent groundwater flow model of the General Separations Area, Savannah River Site, is referred to as the “GSA/PORFLOW” model. GSA/PORFLOW was developed in 2004 by porting an existing General Separations Area groundwater flow model from the FACT code to the PORFLOW code. The preceding “GSA/FACT” model was developed in 1997 using characterization and monitoring data through the mid-1990’s. Both models were manually calibrated to field data. Significantly more field data have been acquired since the 1990’s and model calibration using mathematical optimization software has become routine and recommended practice. The current task involved updating the GSA/PORFLOW model using selected field data current through at least 2015, and use of the PEST code to calibrate the model and quantify parameter uncertainty. This new GSA groundwater flow model is named “GSA2016” in reference to the year in which most development occurred. The GSA2016 model update is intended to address issues raised by the DOE Low-Level Waste (LLW) Disposal Facility Federal Review Group (LFRG) in a 2008 review of the E-Area Performance Assessment, and by the Nuclear Regulatory Commission in reviews of tank closure and Saltstone Disposal Facility Performance Assessments.

54 ENVIRONMENTAL SCIENCES↗

Internal Nozzle Flow Simulations of Gasoline-Like Fuels under Diesel Operating Conditions

Spray formation in internal combustion engines with direct injection is strictly correlated with internal nozzle flow characteristics, which are in turn influenced by fuel physical properties and injector needle motion. This paper pre-sents a series of 3D simulations that model the in-nozzle flow in a 5-hole mini-sac diesel injector. Two gasoline-like naphtha fuels, namely full-range and light naphtha, were tested under operating conditions typical of diesel applica-tions and were compared with n-dodecane, selected from a palette used as diesel surrogates. Validated methodolo-gies from our previous work were employed to account for realistic needle motion. The multi-phase nature of the problem was described by the mixture model assumption with the Volume of Fluid method. Cavitation effects were included by means of the Homogeneous Relaxation Model and turbulence closure was achieved with the Standard k-ε model in an Unsteady Reynolds-Averaged Navier-Stokes formulation. The results revealed that injector perfor-mance and propensity to cavitation are influenced by the fuel properties. Analyses of several physical quantities were carried out to highlight the fuel-to-fuel differences in terms of mass flow rate, discharge coefficients, and fuel vapor volume fraction inside the orifices. A series of parametric investigations was also performed to assess the fuel response to varied fuel injection temperature, injection pressure, and cross-sectional orifice area. For all cases, the strict correlation between cavitation magnitude and saturation pressure was confirmed. Owing to their higher volatil-ity, the two gasoline-like fuels were characterized by higher cavitation across all the simulated conditions. Occur-rence of cavitation was mostly found at the needle seat and at the orifice inlets during the injection event’s transient, when very small gaps exist between the needle and its seat. This behavior tended to disappear at maximum needle lift, where cavitation was absent for all fuels. Differences in mass flow rate between the naphtha fuels and n-dodecane were measured and ascribed to the different densities of the three fuels. Nevertheless, they were found to be smaller than expected, owing to the lower viscosity of the gasoline-like fuels. This beneficial influence of the lower viscosity was shown to be less effective at higher temperature, where the relative viscosity differences de-creased.

Torelli, R.↗

Flow Simulation of Supersonic Inlet with Bypass Annular Duct

A relaxed isentropic compression supersonic inlet is a new concept that produces smaller cowl drag than a conventional inlet, but incurs lower total pressure recovery and increased flow distortion in the (radially) outer flowpath. A supersonic inlet comprising a bypass annulus to the relaxed isentropic compression inlet dumps out airflow of low quality through the bypass duct. A reliable computational fluid dynamics solution can provide considerable useful information to ascertain quantitatively relative merits of the concept, and further provide a basis for optimizing the design. For a fast and reliable performance evaluation of the inlet performance, an equivalent axisymmetric model whose area changes accounts for geometric and physical (blockage) effects resulting from the original complex three-dimensional configuration is proposed. In addition, full three-dimensional calculations are conducted for studying flow phenomena and verifying the validity of the equivalent model. The inlet-engine coupling is carried out by embedding numerical propulsion system simulation engine data into the flow solver for interactive boundary conditions at the engine fan face and exhaust plane. It was found that the blockage resulting from complex three-dimensional geometries in the bypass duct causes significant degradation of inlet performance by pushing the terminal normal shock upstream.

Kim, HyoungJin↗

Implicit Thermochemical Nonequilibrium Compressible Flow Simulations on Unstructured Grids Using GPUs

As next-generation exascale-class systems arrive, existing software must be updated accordingly to effectively utilize these systems. For high concurrency and energy efficiency, many of these systems utilize GPU architectures. In this work, we present a CUDA C++ implementation of FUN3D's thermochemical nonequilibrium capability for turbulent flows. Efficiency is demonstrated at scale using the Summit system at the Oak Ridge Leadership Computing Facility which is representative of future exascale systems. This work enables faster, higher fidelity, and scale-resolving simulations of thermochemical nonequilibrium flows including reentry, hypersonics, and combustion.

CFD, GPU, HPC, Hypersonics, Chemistry↗

Three-dimensional flow simulation about the AFE vehicle in the transitional regime

The direct-simulation Monte Carlo technique is used to analyze the hypersonic rarefied flow about the three-dimensional NASA Aeroassist Flight Experiment vehicle. Results are given for typical transitional flows encountered during the vehicle's atmospheric entry from altitudes of 200-100 km with an entry velocity of 9.9 km/s. It is found that dissociation is important at altitudes of 110 km and below, and that transitional effects are significant even at an altitude of 200 km.

Celenligil, M. Cevdet↗

Demonstration of an automated CFD system for three-dimensional flow simulations

In this paper the capabilities of an automated CFD system which is currently available at NLR are demonstrated. Transonic flow around the AS28G wing/body configuration and hypersonic flow through a generic three-dimensional mixed-compression airbreathing inlet are simulated. An assessment of the level of automation of the current CFD-system is made. The problem-turnaround time lies within the order of a week for both applications.

Vanderburg, J. W.↗

Incompressible viscous flow simulations of the NFAC wind tunnel

The capabilities of an existing 3-D incompressible Navier-Stokes flow solver, INS3D, are extended and improved to solve turbulent flows through the incorporation of zero- and two-equation turbulence models. The two-equation model equations are solved in their high Reynolds number form and utilize wall functions in the treatment of solid wall boundary conditions. The implicit approximate factorization scheme is modified to improve the stability of the two-equation solver. Applications to the 3-D viscous flow inside the 80 by 120 feet open return wind tunnel of the National Full Scale Aerodynamics Complex (NFAC) are discussed and described.

Champney, Joelle Milene↗

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

Transonic flow simulation of prop-fan configurations

This paper presents the results of a study that used a three-dimensional transonic analysis method to analyze flow around a wing-mounted prop-fan configuration. The transonic analysis method employs a cell-oriented finite volume approach to solve the full potential equation in conservative form over a numerically generated surface-fitted grid. The propeller slipstream effects are simulated using a linearized transpiration boundary condition. The results of analyses of flow over clean wing/body, wing/body/nacelle, and wing/body nacelle with slipstream are compared to test data.

Samant, S. S.↗

Hydrocarbon Cold-Flow Simulation of LOX/LH2 Rocket Aerodynamics

In FY19, this project developed a test method which employed heated ethane in the physical simulation of rocket plume aerodynamics. Based on this work, EUS partnered with significant leveraged funding. A portable, self-contained test bed was also fabricated to support validation of the new method. In FY20, this project focused on the formal validation of ethane’s performance as an aerodynamic simulant as well as the development of the analytical and hardware infrastructure required to carry out that task. The test bed was activated, a series of heated-ethane tests on one or more rocket/diffuser configurations were conducted, and characterizations of the performance of the ethane simulation compared to hot-fire were performed. A series of 127 tests were conducted on 13 aerodynamic configurations to evaluate ethane’s performance as a simulant. It was found that the new methodology can reproduce rocket plume flow fields with +/- 0.5% pressure error at ~3% physical scale. The result is a cost-per-test reduction of ~99% compared to the hot-fire subscale testing typically used to obtain equivalent data. Some initial efforts have taken place, and future efforts are being planned with LaRC supersonic retro propulsion and MSFC SLS aerodynamics teams, which are interested in using the ethane system to wind-tunnel test human-rated Mars lander retrorockets and SLS Booster. Details of the prototypical experiments were published as NASA TM 2020-5009122,and have been presented in JANNAF Virtual Conference paper, "New Developments in Retro propulsion Testing for Mars Entry, Descent and Landing".

Daniel Jones↗

Multistage Turbomachinery Flows Simulated Numerically

At the NASA Lewis Research Center, a comprehensive assessment was made of the predictive capability of the average passage flow model as applied to multistage axial-flow compressors. This model, which describes the time-averaged flow field within a typical passage of a blade row embedded in a multistage configuration, is being widely used throughout U.S. aircraft industry as an integral part of their design systems. Rotor flow-angle deviation. In this work, detailed data taken within a four and one-half stage large low-speed compressor were used to assess the weaknesses and strengths of the predictive capabilities of the average passage flow model. The low-speed compressor blading is of modern design and employs stator end-bends. Measurements were made with slow- and high response instrumentation. The high-response measurements revealed the velocity components of both the rotor and stator wakes. From the measured wake profiles, we found that the flow exiting the rotors deviated from the rotor exit metal angle to a lesser degree than was predicted by the average passage flow model. This was found to be due to blade boundary layer transition, which recently has been shown to exist on multistage axial compressor rotor and stator blades, but was not accounted for in the average passage model. Consequently, a model that mimics the effects of blade boundary layer transition, Shih k-epsilon model, was incorporated into the average passage model. Simulations that incorporated this transition model showed a dramatic improvement in agreement with data. The altered model thus improved predictive capability for multistage axial-flow compressors, and this was verified by detailed experimental measurement.

Hathaway, Michael D.↗

Flow simulations for nacelle-propeller configurations using Euler equations

Euler codes for both axisymmetric and general three-dimensional nacelle-propeller flow analysis have been developed. Surface-fitted grids are generated either by a finite difference method or by an algebraic method. The propeller is represented by an actuator disk along a computational plane where proper boundary conditions are assigned to simulate the propeller power loading. Computed results for a NASA SR3 propeller as well as a NASA turboprop configuration are compared with test data. Good agreement has been achieved through the present simulation method.

Yu, N. J.↗

Multispecies Hypersonic Flow Simulator

MHYSA is a conservative finite-difference code to solve the 1D, 2D, and 3D multispecies Navier-Stokes equations for hypersonic flow applications. It supports Cartesian grids only and has the option of including immersed boundaries. - Solves 1D Euler and 2D and 3D Navier-Stokes equations. - Solves the PDEs over Cartesian grids. - Written entirely in C and uses the MPICH library. It also uses OpenMP threads but this is a work-inprogress. - Can be compiled with PETSc (http://www.mcs.anl.gov/petsc/), if available, where it can use PETSc's time integration module TS (http://www.mcs.anl.gov/petsc/petsc-current/src/ts/). - For 3-dimensional simulations, the immersed boundaries can be used to solve over non-Cartesian geometries. MHYSA has been developed to be scalable, and apart from the usual functionalities to solve a system of PDEs on distributed architectures, it provides scalable file I/O functions.

Ghosh, Debojyoti↗

A mass–momentum consistent coupling for mesh-adaptive two-phase flow simulations

Here, we present a novel mass-momentum consistent coupling between a geometric volume-of-fluid scheme and an incompressible flow solver with differing directional-splitting approaches. The advection of the volume fraction is performed using a direction-split algorithm, whereas the momentum advection algorithm uses a traditional unsplit, fractional-step approach. Both algorithms employ finite-volume discretizations based on Cartesian meshes. In solving the mass-momentum consistency problem, momentum fluxes at the cell faces are weighted by the density fluxes based on the already advected volume fraction. The success of our approach lies on introducing a Favre-averaged velocity interpolation at the liquid/gas interface along with a minmod slope limiter. Mesh-convergence studies show that when the minmod slope limiter is used, the two-phase solver retains an accuracy between first and second order, but when a purely upwind scheme is considered, its accuracy drops to first order. Finally, after considering several validation problems, the solver is shown to agree well with reference numerical and experimental data while retaining its robustness and efficiency.

97 MATHEMATICS AND COMPUTING↗