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 595 records · Page 33

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗

A High-Resolution Capability for Large-Eddy Simulation of Jet Flows

A large-eddy simulation (LES) code that utilizes high-resolution numerical schemes is described and applied to a compressible jet flow. The code is written in a general manner such that the accuracy/resolution of the simulation can be selected by the user. Time discretization is performed using a family of low-dispersion Runge-Kutta schemes, selectable from first- to fourth-order. Spatial discretization is performed using central differencing schemes. Both standard schemes, second- to twelfth-order (3 to 13 point stencils) and Dispersion Relation Preserving schemes from 7 to 13 point stencils are available. The code is written in Fortran 90 and uses hybrid MPI/OpenMP parallelization. The code is applied to the simulation of a Mach 0.9 jet flow. Four-stage third-order Runge-Kutta time stepping and the 13 point DRP spatial discretization scheme of Bogey and Bailly are used. The high resolution numerics used allows for the use of relatively sparse grids. Three levels of grid resolution are examined, 3.5, 6.5, and 9.2 million points. Mean flow, first-order turbulent statistics and turbulent spectra are reported. Good agreement with experimental data for mean flow and first-order turbulent statistics is shown.

DeBonis, James R.↗

Liquid-phase mixing of bipropellant doublets

Experimental results of unlike doublet mixing are correlated with an analytically derived equation predicting fluid cavitation. The correlation relates the minimum orifice pressure drop required to initiate cavitation, with the system back pressure, cold flow simulant vapor pressure, and the orifice flow discharge and contraction coefficients. Stream flow instabilities are also visually correlated with the onset of cavitation and orifice discharge coefficient measurements. The influence of cavitation on the characteristic phenomenon of hydraulic flip is observed for both circular and noncircular shaped orifices. For certain intermediate orifice lengths, some noncircular shapes are shown to produce more fully developed flows (shorter recovery lengths) and therefore a more cohesive jet, which in turn yields slightly higher cold flow mixing uniformities than circular shaped orifices of equal absolute length. The particular noncircular shaped elements evaluated are shown to be more sensitive to liquid stream misimpingement than the corresponding circular orifices.

Hoehn, F. W.↗

Vector potential methods

Vector potential and related methods, for the simulation of both inviscid and viscous flows over aerodynamic configurations, are briefly reviewed. The advantages and disadvantages of several formulations are discussed and alternate strategies are recommended. Scalar potential, modified potential, alternate formulations of Euler equations, least-squares formulation, variational principles, iterative techniques and related methods, and viscous flow simulation are discussed.

Hafez, M.↗

Microgravity liquid propellant management

The requirement to settle or to position liquid fluid over the outlet end of a spacecraft propellant tank prior to main engine restart, poses a microgravity fluid behavior problem. Resettlement or reorientation of liquid propellant can be accomplished by providing optimal acceleration to the spacecraft such that the propellant is reoriented over the tank outlet without any vapor entrainment, any excessive geysering, or any other undersirable fluid motion for the space fluid management under microgravity environment. The most efficient technique is studied for propellant resettling through the minimization of propellant usage and weight penalties. Both full scale and subscale liquid propellant tank of Space Transfer Vehicle were used to simulate flow profiles for liquid hydrogen reorientation over the tank outlet. In subscale simulation, both constant and impulsive resettling acceleration were used to simulate the liquid flow reorientation. Comparisons between the constant reverse gravity acceleration and impulsive reverse gravity acceleration to be used for activation of propellant resettlement shows that impulsive reverse gravity thrust is superior to constant reverse gravity thrust.

Hung, R. J.↗

CFD-based Design of LOX Pump Inlet Duct for Reduced Dynamic Loads

Numerical simulations have been completed for a variety of designs for a 90 deg elbow duct. The objective is to identify a design that minimizes the dynamic load entering a LOX turbopump located at the elbow exit. Designs simulated to date indicate that simpler duct geometries result in lower losses. Benchmark simulations have verified that the compressible flow code used in this study is applicable to these incompressible flow simulations.

Rothermel, Jeffry↗

CFD-Based Design of Turbopump Inlet Duct for Reduced Dynamic Loads

Numerical simulations have been completed for a variety of designs for a 90 deg elbow duct. The objective is to identify a design that minimizes the dynamic load entering a LOX turbopump located at the elbow exit. Designs simulated to date indicate that simpler duct geometries result in lower losses. Benchmark simulations have verified that the compressible flow codes used in this study are applicable to these incompressible flow simulations.

Rothermel, Jeffry↗

Computational Fluid Dynamics Analysis Method Developed for Rocket-Based Combined Cycle Engine Inlet

Renewed interest in hypersonic propulsion systems has led to research programs investigating combined cycle engines that are designed to operate efficiently across the flight regime. The Rocket-Based Combined Cycle Engine is a propulsion system under development at the NASA Lewis Research Center. This engine integrates a high specific impulse, low thrust-to-weight, airbreathing engine with a low-impulse, high thrust-to-weight rocket. From takeoff to Mach 2.5, the engine operates as an air-augmented rocket. At Mach 2.5, the engine becomes a dual-mode ramjet; and beyond Mach 8, the rocket is turned back on. One Rocket-Based Combined Cycle Engine variation known as the "Strut-Jet" concept is being investigated jointly by NASA Lewis, the U.S. Air Force, Gencorp Aerojet, General Applied Science Labs (GASL), and Lockheed Martin Corporation. Work thus far has included wind tunnel experiments and computational fluid dynamics (CFD) investigations with the NPARC code. The CFD method was initiated by modeling the geometry of the Strut-Jet with the GRIDGEN structured grid generator. Grids representing a subscale inlet model and the full-scale demonstrator geometry were constructed. These grids modeled one-half of the symmetric inlet flow path, including the precompression plate, diverter, center duct, side duct, and combustor. After the grid generation, full Navier-Stokes flow simulations were conducted with the NPARC Navier-Stokes code. The Chien low-Reynolds-number k-e turbulence model was employed to simulate the high-speed turbulent flow. Finally, the CFD solutions were postprocessed with a Fortran code. This code provided wall static pressure distributions, pitot pressure distributions, mass flow rates, and internal drag. These results were compared with experimental data from a subscale inlet test for code validation; then they were used to help evaluate the demonstrator engine net thrust.

Source record↗

DNS/LES Simulations of Separated Flows at High Reynolds Numbers

Direct numerical simulations (DNS) and large-eddy simulations (LES) simulations of flow through a periodic channel with a constriction are performed using the dynamic Smagorinsky model at two Reynolds numbers of 2800 and 10595. The LES equations are solved using higher order compact schemes. DNS are performed for the lower Reynolds number case using a fine grid and the data are used to validate the LES results obtained with a coarse and a medium size grid. LES simulations are also performed for the higher Reynolds number case using a coarse and a medium size grid. The results are compared with an existing reference data set. The DNS and LES results agreed well with the reference data. Reynolds stresses, sub-grid eddy viscosity, and the budgets for the turbulent kinetic energy are also presented. It is found that the turbulent fluctuations in the normal and spanwise directions have the same magnitude. The turbulent kinetic energy budget shows that the production peaks near the separation point region and the production to dissipation ratio is very high on the order of five in this region. It is also observed that the production is balanced by the advection, diffusion, and dissipation in the shear layer region. The dominant term is the turbulent diffusion that is about two times the molecular dissipation.

Balakumar, P.↗

Accurate simulation of vortical flows

Even with recent supercomputers having a large memory, Navier-Stokes simulations for vortical flows do not provide satisfactory results because of the lack of grid resolution to accurately simulate the strength of separation vortices. To overcome this problem, a zonal method is proposed to increase the number of grid points locally. Interface scheme which is critical for an efficient and stable zonal method is based on the Fortified Navier-Stokes concept. Application to both two-dimensional conical and three-dimensional delta wing problems indicates this simple zonal method can improve the accuracy of vortical flow simulations.

Fujii, Kozo↗

dfnWorks: A discrete fracture network framework for modeling subsurface flow and transport

DFNWORKS is a parallelized computational suite to generate three-dimensional discrete fracture networks (DFN) and simulate flow and transport. Developed at Los Alamos National Laboratory over the past five years, it has been used to study flow and transport in fractured media at scales ranging from millimeters to kilometers. The networks are created and meshed using DFNGEN, which combines FRAM (the feature rejection algorithm for meshing) methodology to stochastically generate three-dimensional DFNs with the LaGriT meshing toolbox to create a high-quality computational mesh representation. The representation produces a conforming Delaunay triangulation suitable for high performance computing finite volume solvers in an intrinsically parallel fashion. Flow through the network is simulated in dfnFlow, which utilizes the massively parallel subsurface flow and reactive transport finite volume code PFLOTRAN. A Lagrangian approach to simulating transport through the DFN is adopted within DFNTRANS to determine pathlines and solute transport through the DFN. Example applications of this suite in the areas of nuclear waste repository science, hydraulic fracturing and CO 2 sequestration are also included.

42 ENGINEERING↗

Direct simulations of compressible wall-bounded turbulence

Several direct numerical simulations of high-speed turbulent Couette flow were performed with a new spectral code. Mach numbers up to three and a Reynolds number of 3000 were used. A new time-integration scheme was developed to handle Mach numbers above 1.5, which require greater accuracy and stability than lower Mach numbers. At low Mach number, the large streamwise eddies found by M. J. Lee in high incompressible Couette flow simulations were reproduced. At higher Mach numbers these structures still exist, but they become considerably less organized (although the disorganization may be a function of the spanwise box size). While the same types of vortical structures seen in the incompressible flow are observed at higher Mach numbers, a new structure involving the divergence of the velocity is also observed. This structure is generally associated with low shear areas next to the walls, but it has not been determined whether it is a cause or an effect of the low shear. A 'nonphysical' simulation was performed to determine by what mechanism the Mach number affects the flow. It appears that pressure gradient (acoustic) effects are more important than variable viscosity effects in determining the wall shear, but the size of vortical structures is determined more by the local kinematic viscosity. Low-order mean statistics are provided to help quantify these effects.

Buell, Jeffrey C.↗

The transport of suspended sediment in the wake of a marine hydrokinetic turbine: Simulations via a validated Discrete Random Walk (DRW) model

The interaction of sediment with Marine Hydrokinetic (MHK) Turbines is studied via turbulent flow simulation coupled with a model for particle dynamics. The Discrete Random Walk (DRW) model is commonly used for simulations of particle dispersion in environmental flows. A method to improve the accuracy of the DRW model in conjunction with RANS simulations of the flow around MHK turbines is presented. A key issues for the use of the DRW model in marine applications is identified, finding the right characteristic eddy lifetime scale, and a simple methodology based on G.I. Taylor's classical dispersion theory is develop to calibrate this physical model parameter. The model is validated against experiments in the literature and the physics of suspended sediment transport in the wake of an MHK turbine is studied with this improved method. The qualitative changes observed in the particle trajectories agree well with theoretical and empirical observations and the quantitive trends allow for comparison of multiple particle cases under different flow conditions, confirming that this implementation of CFD can become a powerful tool to select the most challenging cases, among the thousands of potential studies on the environmental effects of MHK energy, for further study with experiments and pilot projects.

54 ENVIRONMENTAL SCIENCES↗

Time reversal symmetry breaking and odd viscosity in active fluids: Green–Kubo and NEMD results

Active fluids, which are driven at the microscale by non-conservative forces, are known to exhibit novel transport phenomena due to the breaking of time reversal symmetry. Recently, Epstein and Mandadapu [arXiv:1907.10041 (2019)] obtained Green-Kubo relations for the full set of viscous coefficients governing isotropic chiral active fluids, including the so-called odd viscosity, invoking Onsager's regression hypothesis for the decay of fluctuations in active non-equilibrium steady states. In this Communication, we test these Green-Kubo relations using molecular dynamics simulations of a canonical model system consisting of actively torqued dumbbells. We find the resulting odd and shear viscosity values from the Green-Kubo relations to be in good agreement with values measured independently through non-equilibrium molecular dynamics flow simulations. This provides a test of the Green-Kubo relations and lends support to the application of the Onsager regression hypothesis in relation to viscous behaviors of active matter systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two inviscid computational simulations of separated flow about airfoils

Two inviscid computational simulations of separated flow about airfoils are described. The basic computational method is the line relaxation finite-difference method. Viscous separation is approximated with inviscid free-streamline separation. The point of separation is specified, and the pressure in the separation region is calculated. In the first simulation, the empiricism of constant pressure in the separation region is employed. This empiricism is easier to implement with the present method than with singularity methods. In the second simulation, acoustic theory is used to determine the pressure in the separation region. The results of both simulations are compared with experiment.

Barnwell, R. W.↗

Experiments On Flow In A Coronary Artery

Report describes experiments on simulated flow of blood in atherosclerotic human coronary artery. Experiments performed on polyurethane cast made from S-shaped coronary artery of cadaver. Sucrose solution with viscosity of blood pumped through cast at physiologically realistic rates, and flow made pulsatile by mechanism alternately compressing and releasing elastic tube just upstream of cast.

Back, Lloyd H.↗

More-Accurate Model of Flows in Rocket Injectors

An improved computational model for simulating flows in liquid-propellant injectors in rocket engines has been developed. Models like this one are needed for predicting fluxes of heat in, and performances of, the engines. An important part of predicting performance is predicting fluctuations of temperature, fluctuations of concentrations of chemical species, and effects of turbulence on diffusion of heat and chemical species. Customarily, diffusion effects are represented by parameters known in the art as the Prandtl and Schmidt numbers. Prior formulations include ad hoc assumptions of constant values of these parameters, but these assumptions and, hence, the formulations, are inaccurate for complex flows. In the improved model, these parameters are neither constant nor specified in advance: instead, they are variables obtained as part of the solution. Consequently, this model represents the effects of turbulence on diffusion of heat and chemical species more accurately than prior formulations do, and may enable more-accurate prediction of mixing and flows of heat in rocket-engine combustion chambers. The model has been implemented within CRUNCH CFD, a proprietary computational fluid dynamics (CFD) computer program, and has been tested within that program. The model could also be implemented within other CFD programs.

Hosangadi, Ashvin↗