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 487 records · Page 27

The Isolated Synthetic Jet in Crossflow: A Benchmark for Flow Control Simulation

An overview of the data acquisition, reduction, and uncertainty of experimental measurements made of the flowfield created by the interaction of an isolated synthetic jet and a turbulent boundary layer is presented. The experimental measurements were undertaken to serve as the second of three computational fluid dynamics validation databases for Active Flow Control. The validation databases were presented at the NASA Langley Research Center Workshop on CFD Validation of Synthetic Jets and Turbulent Separation Control in March, 2004. Detailed measurements were made to document the boundary conditions for the flow and also for the phase-averaged flowfield itself. Three component Laser-Doppler Velocimetry, 2-D Particle Image Velocimetry, and Stereo Particle Image Velocimetry were utilized to document the phase-averaged velocity field and the turbulent stresses.

Schaeffler, Norman W.↗

Numerical simulation of viscous transonic airfoil flows

Numerical simulations of transonic airfoil flows using the Reynolds-averaged Navier-Stokes equations and various turbulence models are presented and compared with experimental data. Three different airfoils were investigated under varying flow conditions ranging from subcritical unseparated flows to supercritical separated flows. The turbulence models investigated consisted of three zero-equation models and one two-equation model. For unseparated flows involving weak viscous-inviscid interactions, the four models were comparable in their agreement with experiment. For separated flows involving strong viscous-inviscid interactions, the nonequilibrium zero-equation model of Johnson and King gave the best overall agreement with experiment.

Coakley, Thomas J.↗

Highly Variable Cycle Nozzle Concept: Validation of Flow and Noise Predictions

Results from experimental and numerical studies of highly Variable Cycle (HVC) exhaust model were presented. The model was designed and fabricated under a Supersonics NRA awarded to Rolls-Royce. The model had a lobed mixer for the core stream nozzle, and elliptic fan stream nozzle, and an ejector. Experiments included far-field acoustic array, phased array, and Particle Image Velocimetry (PIV) measurements. Numerical studies included flow simulations using the WIND-US code and far-field acoustic solutions using an acoustic analogy developed by Goldstein (2003) and Leib and Goldstein (2011). Far-field acoustic measurements showed increased noise levels over the round baseline nozzle when using non-static forward flight conditions. Phased array measurements showed noise sources near the ejector doors when tones were produced for small ejector door positions. Ejector door separation identified in the experiments was reproduced in the numerical flow simulations. Acoustic solutions were unable to match levels measured in the peak jet noise direction indicating additional development work is needed to predict noise from highly three-dimensional flows.

Halbig, Michael C.↗

An Initial Multi-Domain Modeling of an Actively Cooled Structure

A methodology for the simulation of turbine cooling flows is being developed. The methodology seeks to combine numerical techniques that optimize both accuracy and computational efficiency. Key components of the methodology include the use of multiblock grid systems for modeling complex geometries, and multigrid convergence acceleration for enhancing computational efficiency in highly resolved fluid flow simulations. The use of the methodology has been demonstrated in several turbo machinery flow and heat transfer studies. Ongoing and future work involves implementing additional turbulence models, improving computational efficiency, adding AMR.

Steinthorsson, Erlendur↗

Interacting scales and energy transfer in isotropic turbulence

The dependence of the energy transfer process on the disparity of the interacting scales is investigated in the inertial and far-dissipation ranges of isotropic turbulence. The strategy for generating the simulated flow fields and the choice of a disparity parameter to characterize the scaling of the interactions is discussed. The inertial range is found to be dominated by relatively local interactions, in agreement with the Kolmogorov assumption. The far-dissipation is found to be dominated by relatively non-local interactions, supporting the classical notion that the far-dissipation range is slaved to the Kolmogorov scales. The measured energy transfer is compared with the classical models of Heisenberg, Obukhov, and the more detailed analysis of Tennekes and Lumley. The energy transfer statistics measured in the numerically simulated flows are found to be nearly self-similar for wave numbers in the inertial range. Using the self-similar form measured within the limited scale range of the simulation, an 'ideal' energy transfer function and the corresponding energy flux rate for an inertial range of infinite extent are constructed. From this flux rate, the Kolmogorov constant is calculated to be 1.5, in excellent agreement with experiments.

Zhou, YE↗

Implicitly solving phase appearance and disappearance problems using two-fluid six-equation model

Phase appearance and disappearance issue presents serious numerical challenges in two-phase flow simulations using the two-fluid six-equation model. Numerical challenges arise from the singular equation system when one phase is absent, as well as from the discontinuity in the solution space when one phase appears or disappears. In this work, a high-resolution spatial discretization scheme on staggered grids and fully implicit methods were applied for the simulation of two-phase flow problems using the two-fluid six-equation model. A Jacobian-free Newton-Krylov (JFNK) method was used to solve the discretized nonlinear problem. An improved numerical treatment was proposed and proved to be effective to handle the numerical challenges. The treatment scheme is conceptually simple, easy to implement, and does not require explicit truncations on solutions, which is essential to conserve mass and energy. Various types of phase appearance and disappearance problems relevant to thermal-hydraulics analysis have been investigated, including a sedimentation problem, an oscillating manometer problem, a non-condensable gas injection problem, a single-phase flow with heat addition problem and a subcooled flow boiling problem. Successful simulations of these problems demonstrate the capability and robustness of the proposed numerical methods and numerical treatments. As a result, volume fraction of the absent phase can be calculated effectively as zero.

42 ENGINEERING↗

LOX/hydrogen coaxial injector atomization test program

Quantitative information about the atomization of injector sprays is needed to improve the accuracy of computational models that predict the performance and stability margin of liquid propellant rocket engines. To obtain this data, a facility for the study of spray atomization is being established at NASA-Lewis to determine the drop size and velocity distributions occurring in vaporizing liquid sprays at supercritical pressures. Hardware configuration and test conditions are selected to make the cold flow simulant testing correspond as closely as possible to conditions in liquid oxygen (LOX)/gaseous H2 rocket engines. Drop size correlations from the literature, developed for liquid/gas coaxial injector geometries, are used to make drop size predictions for LOX/H2 coaxial injectors. The mean drop size predictions for a single element coaxial injector range from 0.1 to 2000 microns, emphasizing the need for additional studies of the atomization process in LOX/H2 engines. Selection of cold flow simulants, measurement techniques, and hardware for LOX/H2 atomization simulations are discussed.

Zaller, M.↗

LOX/Hydrogen Coaxial Injector Atomization Test Program

Quantitative information about the atomization of injector sprays is needed to improve the accuracy of computational models that predict the performance and stability margin of liquid propellant rocket engines. To obtain this data, a facility for the study of spray atomization is being established at NASA-Lewis to determine the drop size and velocity distributions occurring in vaporizing liquid sprays at supercritical pressures. Hardware configuration and test conditions are selected to make the cold flow simulant testing correspond as closely as possible to conditions in liquid oxygen (LOX)/gaseous H2 rocket engines. Drop size correlations from the literature, developed for liquid/gas coaxial injector geometries, are used to make drop size predictions for LOX/H2 coaxial injectors. The mean drop size predictions for a single element coaxial injector range from 0.1 to 2000 microns, emphasizing the need for additional studies of the atomization process in LOX/H2 engines. Selection of cold flow simulants, measured techniques, and hardware for LOX/H2 atomization simulations are discussed.

Zaller, M.↗

Time-Lapse Gravity Monitoring of CO2 Migration Based on Numerical Modeling of a Faulted Storage Complex

In this study, the performance of both surface and borehole time-lapse gravity monitoring to detect CO2 leakage from a carbon storage site is evaluated. Several hypothetical scenarios of CO2 migration in a leaky fault, and thief zones at different depths at the Kimberlina site (California, USA) constitute the basis of the approach. The CO2 displacement is simulated using the TOUGH2 simulator applied to a detailed geological model of the site. The gravity responses to these CO2 plumes are simulated using forward modeling with sensors at ground surface and in vertical boreholes. Results of inversion on one scenario are also presented. The surface-based gravity responses obtained for the different leakage scenarios demonstrate that leakage can be detected at the surface in all the scenarios but the time to detection is highly variable (10 to 40 years) and dependent on the detection threshold considered. Borehole measurements of the vertical component of gravity provide excellent constrains in depth when they are located in proximity of the density anomaly associated with the presence of CO2, thus discriminating multiple leaks in different thief zones. Joint inversion of surface and borehole data can bring valuable information of the occurrence of leakages and their importance by providing a reasonable estimate of mass of displaced fluids. This study demonstrates the importance of combining multiphase flow simulations with gravity modeling in order to define if and when gravity monitoring would be applicable at a given storage site.

Time-lapse gravity monitoring, leak detection, CCS↗

A Moving Embedded Boundary Approach for the Compressible Navier-Stokes Equations in a Block-Structured Adaptive Refinement Framework

A computational technique has been developed to perform compressible flow simulations involving moving boundaries using an embedded boundary approach within the block-structured adaptive mesh refinement (SAMR) framework of AMReX [1], [91], [92]. We leverage the SAMR capability to obtain quantitatively accurate results whilst using robust, second-order finite volume schemes. A conservative, unsplit, cut-cell approach is utilized and a ghost-cell approach is developed for computing the flux on the moving, embedded boundary faces. A third-order least-squares formulation has been developed to compute the wall velocity gradients, and was found to significantly improve the performance of the solver in terms of the quantitative comparison of surface quantities such as the skin friction coefficient. Various test cases are performed to validate the method, and compared with analytical, experimental, and other numerical results in literature. Inviscid and viscous test cases are performed that span a wide regime of flow speeds - acoustic (harmonically pulsating sphere), smooth flows (expansion fan created by a receding piston) and flows with shocks (shock-cylinder interaction, shock-wedge interaction, pitching NACA 0012 airfoil and shock-cone interaction). A closed system with moving boundaries - an oscillating piston in a cylinder, showed that the percentage error in mass within the system decreases with refinement, demonstrating that the numerical scheme is conservative with grid refinement, but is not discretely conservative. Viscous test cases involve that of a horizontally moving cylinder at Re = 40, an inline oscillating cylinder at Re = 100, and a transversely oscillating cylinder at Re = 185. The judicious use of adaptive mesh refinement with appropriate refinement criteria to capture the regions of interest leads to well-resolved flow features, and good quantitative comparison is observed with the results available in literature.

adaptive refinement↗

Solute mixing regulates heterogeneity of mineral precipitation in porous media

Here, synchrotron X-ray microtomography was used to track the spatiotemporal evolution of mineral precipitation and the consequent alteration of the pore structure. Column experiments were conducted by injecting CaCl 2 and NaHCO 3 solutions into granular porous media either as a premixed supersaturated solution (external mixing) or as separate solutions that mixed within the specimen (internal mixing). The two mixing modes produced distinct mineral growth patterns. While internal mixing promoted transverse heterogeneity with precipitation at the mixing zone, external mixing favored relatively homogeneous precipitation along the flow direction. The impact of precipitation on pore water flow and permeability was assessed via 3-D flow simulations, which indicated anisotropic permeability evolution for both mixing modes. Under both mixing modes, precipitation decreased the median pore size and increased the skewness of the pore size distribution. Such similar pore-scale evolution patterns suggest that the clogging of individual pores depends primarily on local supersaturation state and pore geometry.

42 ENGINEERING↗

JACC: Leveraging HPC Meta-Programming and Performance Portability with the Just-in-Time and LLVM-based Julia Language

We present JACC (Julia for Accelerators), the first high-level, and performance-portable model for the just-in-time and LLVM-based Julia language. JACC provides a unified and lightweight front end across different back ends available in Julia, enabling the same Julia code to run efficiently on many HPC CPU and GPU targets. We evaluated the performance of JACC for common HPC kernels as well as for the most computationally demanding kernels used in applications, HPCCG, a supercomputing benchmark test for sparse domains, and HARVEY, a blood flow simulator to assist in the diagnosis and treatment of patients suffering from vascular diseases. We carried out the performance analysis on the most advanced US DOE supercomputers: Aurora, Frontier, and Perlmutter. Overall, we show that JACC has a negligible overhead versus vendor-specific solutions, reporting GPU speedups with no extra cost to programmability.

Valero-Lara, Pedro↗

A new approach to flow problems past a screen

A new method to simulate flows past a two-dimensional porous plate is proposed. The method is a simple application of a discrete vortex method which has been well established for flows past an inclined flat plate. In this method, a flow past a porous plate is described as superposition of the velocity potential of a uniform flow on that of a flow past a flat plate. The effect of porosity is replaced by that of a mass flow rate of the approaching freestream passing through the plate. Calculated flow features are compared with experiments. It is shown that the present method is effective in the simulation of the flow past a porous plate.

Inoue, O.↗

A Lattice Boltzmann Method for Turbomachinery Simulations

Lattice Boltzmann (LB) Method is a relatively new method for flow simulations. The start point of LB method is statistic mechanics and Boltzmann equation. The LB method tries to set up its model at molecular scale and simulate the flow at macroscopic scale. LBM has been applied to mostly incompressible flows and simple geometry.

Hsu, A. T.↗

Integration of Electromagnetic Geophysics Forward Simulation in Coupled Flow and Geomechanics for Monitoring a Gas Hydrate Deposit Located in the Ulleung Basin, East Sea, Korea

We investigate the feasibility of electromagnetic (EM) geophysics methods to detect the dissociation of gas hydrate specifically from a gas hydrate deposit located in the Ulleung Basin, East Sea, Korea via an integrated flow-geomechanics-EM geophysics simulation. To this end, coupled flow and geomechanics simulation is first performed with the multiple porosity model employed, where a mixed formulation with the finite volume (FV) and finite element (FE) methods are taken for the flow and geomechanics, respectively. From the saturation and porosity fields obtained from the coupled flow and geomechanics, the electrical conductivity model is established for the EM simulation. Solving the partial differential equation of electrical diffusion which is linearized using the 3D finite element method (FEM), the EM fields are then computed. For numerical experiments, particularly two approaches in the configuration for the EM methods are compared in this contribution: the surface-to-surface and the surface-to-borehole methods. When the surface-to-surface EM method is employed, the EM is found to be less sensitive, implying low detectability. Especially for the short term of production, the low detectability is attributed to the similarity of electrical resistivity between the dissociated gas (CH4) and hydrate as well as the specific dissociation pattern within the intercalated composites of the field. On the other hand, when the surface-to-borehole EM method is employed, its sensitivity to capture the produced gas flow is improved, confirming its detectability in monitoring gas flow. Hence, the EM geophysics simulation integrated with coupled flow and geomechanics can be a potential tool for monitoring gas hydrate deposits.

depressurization↗

Simulation of Fluid Flow and Collection Efficiency for a SEA Inc. Multi-Element Probe and Ice Crystal Detector Using GlennICE

Numerical simulation results of fluid flow and collection efficiency of the Science Engineering Associates Inc. Multi-Element Probe (Multiwire) and Ice Crystal Detector (ICD) are presented. Fluid flow simulations were conducted using NASA's FUN3D while collection efficiency simulations were conducted using NASA's LEWICE3D software and GlennICE software. For both probes, 3D unsteady flow results were time averaged. Simulations were computed for freestream velocities ranging from 85 to 185 m/s and freestream total pressures of 44.8 and 93.1 kPa. Collection efficiency results were computed for four spherical particle diameter sizes of 5, 20, 50, and 100 µm. GlennICE collection efficiency results for the multiwire were compared with previously published collection efficiency values calculated using LEWICE3D. Numerical collection efficiency results for the Ice Crystal Detector are presented for the first time.

aircraft, icing, Computational Fluid Dynamics↗

iMatTOUGH: An open-source Matlab-based graphical user interface for pre- and post-processing of TOUGH2 and iTOUGH2 models

TOUGH2 and iTOUGH2 are powerful models that simulate the heat and fluid flows in porous and fracture media, and perform parameter estimation, sensitivity analysis and uncertainty propagation analysis. However, setting up the input files is not only tedious, but error prone, and processing output files is time consuming. Here, we present an open source Matlab-based tool (iMatTOUGH) that supports the generation of all necessary inputs for both TOUGH2 and iTOUGH2 and visualize their outputs. The tool links the inputs of TOUGH2 and iTOUGH2, making sure the two input files are consistent. It supports the generation of rectangular computational mesh, i.e., it automatically generates the elements and connections as well as their properties as required by TOUGH2. The tool also allows the specification of initial and time-dependent boundary conditions for better subsurface heat and water flow simulations. The effectiveness of the tool is illustrated by an example that uses TOUGH2 and iTOUGH2 to estimate soil hydrological and thermal properties from soil temperature data and simulate the heat and water flows at the Rifle site in Colorado.

58 GEOSCIENCES↗