Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “viscous effects”

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 397 records · Page 22

Development of Computational Aeroacoustics Code for Jet Noise and Flow Prediction

Accurate prediction of jet fan and exhaust plume flow and noise generation and propagation is very important in developing advanced aircraft engines that will pass current and future noise regulations. In jet fan flows as well as exhaust plumes, two major sources of noise are present: large-scale, coherent instabilities and small-scale turbulent eddies. In previous work for the NASA Glenn Research Center, three strategies have been explored in an effort to computationally predict the noise radiation from supersonic jet exhaust plumes. In order from the least expensive computationally to the most expensive computationally, these are: 1) Linearized Euler equations (LEE). 2) Very Large Eddy Simulations (VLES). 3) Large Eddy Simulations (LES). The first method solves the linearized Euler equations (LEE). These equations are obtained by linearizing about a given mean flow and the neglecting viscous effects. In this way, the noise from large-scale instabilities can be found for a given mean flow. The linearized Euler equations are computationally inexpensive, and have produced good noise results for supersonic jets where the large-scale instability noise dominates, as well as for the tone noise from a jet engine blade row. However, these linear equations do not predict the absolute magnitude of the noise; instead, only the relative magnitude is predicted. Also, the predicted disturbances do not modify the mean flow, removing a physical mechanism by which the amplitude of the disturbance may be controlled. Recent research for isolated airfoils' indicates that this may not affect the solution greatly at low frequencies. The second method addresses some of the concerns raised by the LEE method. In this approach, called Very Large Eddy Simulation (VLES), the unsteady Reynolds averaged Navier-Stokes equations are solved directly using a high-accuracy computational aeroacoustics numerical scheme. With the addition of a two-equation turbulence model and the use of a relatively coarse grid, the numerical solution is effectively filtered into a directly calculated mean flow with the small-scale turbulence being modeled, and an unsteady large-scale component that is also being directly calculated. In this way, the unsteady disturbances are calculated in a nonlinear way, with a direct effect on the mean flow. This method is not as fast as the LEE approach, but does have many advantages to recommend it; however, like the LEE approach, only the effect of the largest unsteady structures will be captured. An initial calculation was performed on a supersonic jet exhaust plume, with promising results, but the calculation was hampered by the explicit time marching scheme that was employed. This explicit scheme required a very small time step to resolve the nozzle boundary layer, which caused a long run time. Current work is focused on testing a lower-order implicit time marching method to combat this problem.

Keith, Theo G., Jr.↗

An Aerodynamic Analysis of a Spinning Missile with Dithering Canards

A generic spinning missile with dithering canards is used to demonstrate the utility of an overset structured grid approach for simulating the aerodynamics of rolling airframe missile systems. The approach is used to generate a modest aerodynamic database for the generic missile. The database is populated with solutions to the Euler and Navier-Stokes equations. It is used to evaluate grid resolution requirements for accurate prediction of instantaneous missile loads and the relative aerodynamic significance of angle-of-attack, canard pitching sequence, viscous effects, and roll-rate effects. A novel analytical method for inter- and extrapolation of database results is also given.

Meakin, Robert L.↗

Shock-Wave/Boundary-Layer Interactions in Hypersonic Low Density Flows

Results of numerical simulations of Mach 10 air flow over a hollow cylinder-flare and a double-cone are presented where viscous effects are significant. The flow phenomena include shock-shock and shock- boundary-layer interactions with accompanying flow separation, recirculation, and reattachment. The purpose of this study is to promote an understanding of the fundamental gas dynamics resulting from such complex interactions and to clarify the requirements for meaningful simulations of such flows when using the direct simulation Monte Carlo (DSMC) method. Particular emphasis is placed on the sensitivity of computed results to grid resolution. Comparisons of the DSMC results for the hollow cylinder-flare (30 deg.) configuration are made with the results of experimental measurements conducted in the ONERA RSCh wind tunnel for heating, pressure, and the extent of separation. Agreement between computations and measurements for various quantities is good except that for pressure. For the same flow conditions, the double- cone geometry (25 deg.- 65 deg.) produces much stronger interactions, and these interactions are investigated numerically using both DSMC and Navier-Stokes codes. For the double-cone computations, a two orders of magnitude variation in free-stream density (with Reynolds numbers from 247 to 24,7 19) is investigated using both computational methods. For this range of flow conditions, the computational results are in qualitative agreement for the extent of separation with the DSMC method always predicting a smaller separation region. Results from the Navier-Stokes calculations suggest that the flow for the highest density double-cone case may be unsteady; however, the DSMC solution does not show evidence of unsteadiness.

Moss, James N.↗

Introduction of the ASP3D Computer Program for Unsteady Aerodynamic and Aeroelastic Analyses

A new computer program has been developed called ASP3D (Advanced Small Perturbation 3D), which solves the small perturbation potential flow equation in an advanced form including mass-consistent surface and trailing wake boundary conditions, and entropy, vorticity, and viscous effects. The purpose of the program is for unsteady aerodynamic and aeroelastic analyses, especially in the nonlinear transonic flight regime. The program exploits the simplicity of stationary Cartesian meshes with the movement or deformation of the configuration under consideration incorporated into the solution algorithm through a planar surface boundary condition. The new ASP3D code is the result of a decade of developmental work on improvements to the small perturbation formulation, performed while the author was employed as a Senior Research Scientist in the Configuration Aerodynamics Branch at the NASA Langley Research Center. The ASP3D code is a significant improvement to the state-of-the-art for transonic aeroelastic analyses over the CAP-TSD code (Computational Aeroelasticity Program Transonic Small Disturbance), which was developed principally by the author in the mid-1980s. The author is in a unique position as the developer of both computer programs to compare, contrast, and ultimately make conclusions regarding the underlying formulations and utility of each code. The paper describes the salient features of the ASP3D code including the rationale for improvements in comparison with CAP-TSD. Numerous results are presented to demonstrate the ASP3D capability. The general conclusion is that the new ASP3D capability is superior to the older CAP-TSD code because of the myriad improvements developed and incorporated.

Batina, John T.↗

Application of the ASP3D Computer Program to Unsteady Aerodynamic and Aeroelastic Analyses

A new computer program has been developed called ASP3D (Advanced Small Perturbation - 3D), which solves the small perturbation potential flow equation in an advanced form including mass-consistent surface and trailing wake boundary conditions, and entropy, vorticity, and viscous effects. The purpose of the program is for unsteady aerodynamic and aeroelastic analyses, especially in the nonlinear transonic flight regime. The program exploits the simplicity of stationary Cartesian meshes with the movement or deformation of the configuration under consideration incorporated into the solution algorithm through a planar surface boundary condition. The paper presents unsteady aerodynamic and aeroelastic applications of ASP3D to assess the time dependent capability and demonstrate various features of the code.

Batina, John T.↗

A Non-Cut Cell Immersed Boundary Method for Use in Icing Simulations

This paper describes a computational fluid dynamic method used for modelling changes in aircraft geometry due to icing. While an aircraft undergoes icing, the accumulated ice results in a geometric alteration of the aerodynamic surfaces. In computational simulations for icing, it is necessary that the corresponding geometric change is taken into consideration. The method used, herein, for the representation of the geometric change due to icing is a non-cut cell Immersed Boundary Method (IBM). Computational cells that are in a body fitted grid of a clean aerodynamic geometry that are inside a predicted ice formation are identified. An IBM is then used to change these cells from being active computational cells to having properties of viscous solid bodies. This method has been implemented in the NASA developed node centered, finite volume computational fluid dynamics code, FUN3D. The presented capability is tested for two-dimensional airfoils including a clean airfoil, an iced airfoil, and an airfoil in harmonic pitching motion about its quarter chord. For these simulations velocity contours, pressure distributions, coefficients of lift, coefficients of drag, and coefficients of pitching moment about the airfoil's quarter chord are computed and used for comparison against experimental results, a higher order panel method code with viscous effects, XFOIL, and the results from FUN3D's original solution process. The results of the IBM simulations show that the accuracy of the IBM compares satisfactorily with the experimental results, XFOIL results, and the results from FUN3D's original solution process.

Sarofeen, Christian M.↗

Comparison of Numerically Simulated and Experimentally Measured Performance of a Rotating Detonation Engine

A quasi-two-dimensional, computational fluid dynamic (CFD) simulation of a rotating detonation engine (RDE) is described. The simulation operates in the detonation frame of reference and utilizes a relatively coarse grid such that only the essential primary flow field structure is captured. This construction and other simplifications yield rapidly converging, steady solutions. Viscous effects, and heat transfer effects are modeled using source terms. The effects of potential inlet flow reversals are modeled using boundary conditions. Results from the simulation are compared to measured data from an experimental RDE rig with a converging-diverging nozzle added. The comparison is favorable for the two operating points examined. The utility of the code as a performance optimization tool and a diagnostic tool are discussed.

combustion↗

Total Temperature Measurements Using a Rearward Facing Probe in Supercooled Liquid Droplet and Ice Crystal Clouds

Engine Icing Performance loss: rollback, surge, flameout, and even internal engine damagePartial melting and refreeze of ice inside engine core (Mason et al., 2006). Ingestion of ice crystals and aggregates, mixed-phase droplets, or supercooled liquid dropletsNeed to better understand the conditions and properties that lead to engine icing.Simulation and analysis (physical and computational, and modeling)Test facilities (PSL, NRC, ...). Thermal and computational models and analysisProbesMultiple probes (aerothermal probes and ice cloud characterization probes and techniques). Total temperatureTraditional total temperature probes (vented forward facing)Heated total temperature probes (Goodyear). Rearward facing (developmental). Total temperature relevance. Thermal interaction between the icing cloud and air flow impinging particles contribute to kinetic heating effect (Gent et al., 2000). Measurement considerations Temperature sensor accuracy. Incomplete recovery of total temperature. Thermal surfaces (sources and sinks). Flow effects (viscous losses). Debris contamination, including icing and ice ingestion.

Agui, Juan H.↗

Packed Bed Reactor Experiment (PBRE)-2: Pressure Drop Measurements in Microgravity

Single and two-phase flow through porous media are encountered on ground and in space in numerous applications related to life support systems, fuel cells, chemical/materials processing and transporting nutrients to plants. These systems operate differently in the microgravity environment encountered in space travel because the density differences no longer cause the phases to separate or “drain” under the body force of gravity. In the absence of gravity, the interfacial or capillary forces play a more significant role in determining the operational variables such as phase distribution, liquid holdup, and pressure drop especially when the liquid inertia and viscous forces are at minimum. Lower liquid flow rates in air-water two-phase flows represent the case for most of the space processes that depend on two-phase flows in porous media. The Packed Bed Reactor Experiment (PBRE) was developed and flown on the International Space Station (ISS) by NASA to extend a series of fundamental studies of gas-liquid flows through porous media. The goal of the this series of flight experiments (PBRE-1 and PBRE-2) is to better understand the hydrodynamics of two phase flows in porous media in flow regimes dominated by inertia and viscous effects and understanding the role of the interfacial forces and their role in determining the flow dynamics. The flight experiments were preceded by a series of low gravity experiments performed on the low gravity aircraft. These experiments led to the development of a semi-empirical two-phase pressure drop and flow pattern transition in support of a number of reactor beds planned for water reclamation processes onboard of the international space station. The PBRE flight experiment was designed to deliver a wide range of tightly controlled gas (nitrogen) and liquid (water) flows to one of two interchangeable test sections differing only in the type of internal packing material. In this paper, we describe the PBRE-2 experiment and present preliminary low gravity data on pressure drop and compare with semi-empirical model pressure drop predictions. This work is supported by the National Aeronautics and Space Administration (NASA) under grant no. 80NSSC20K0830

Suratman Number↗

Pressure Gradient and Gas Hold-Up in Two-Phase Flows through Porous Media in Microgravity

Single and two-phase flow through porous media are encountered on ground and in space in numerous applications related to life support systems, fuel cells, chemical/materials processing and transporting nutrients to plants. These systems operate differently in the microgravity environment encountered in space travel because the density differences no longer cause the phases to separate or “drain” under the body force of gravity. In the absence of gravity, the interfacial or capillary forces play a more significant role in determining the operational variables such as phase distribution, liquid holdup, and pressure drop especially when the liquid inertia and viscous forces are at minimum. Lower liquid flow rates in air-water two-phase flows represent the case for most of the space processes that depend on two-phase flows in porous media. The Packed Bed Reactor Experiment (PBRE) was developed and flown on the International Space Station (ISS) by NASA to extend a series of fundamental studies of gas-liquid flows through porous media. The goal of the this series of flight experiments (PBRE-1 and PBRE-2) is to better understand the hydrodynamics of two phase flows in porous media in flow regimes dominated by inertia and viscous effects and understanding the role of the interfacial forces and their role in determining the flow dynamics. The flight experiments were preceded by a series of low gravity experiments performed on the low gravity aircraft. These experiments led to the development of a semi-empirical two-phase pressure drop and flow pattern transition in support of a number of reactor beds planned for water reclamation processes onboard of the international space station. The PBRE flight experiment was designed to deliver a wide range of tightly controlled gas (nitrogen) and liquid (water) flows to one of two interchangeable test sections differing only in the type of internal packing material. In this presentation, we describe the PBRE-2 experiment and present preliminary low gravity data on pressure drop and gas hold-up and compare with semi-empirical model pressure drop predictions.

Brian Motil↗

The Influence of Adaptive Mesh Refinement on the Prediction of Vortex Interactions about a Generic Missile Airframe

The complex interaction of forebody and wing vortices significantly impacts missile aerodynamics. The formation of these vortices involves smooth regions of the geometry or geometric discontinuities like leading edges, trailing edges, tips, and corners. Regions of supersonic flow and complex shock topologies interact with boundary layers and vortices. Smooth-body separation and 3D viscous effects strain current Reynolds-averaged Navier-Stokes (RANS) techniques. The quantification and control of discretization error is critical to obtaining reliable simulation results and often turbulence model assessments are made in the presence of unquantified (and potentially large) discretization errors. Two mesh adaptation schemes are applied to steady RANS simulations. Multiscale unstructured mesh adaptation is applied to control interpolation error estimates of the Mach field, which resolves boundary layers, vortices, and shocks. A dual-mesh approach with overset communication is applied between an expert-crafted near-body unstructured mesh and an adaptive off-body Cartesian mesh refined with Q-criterion scaled by the strain tensor magnitude. A generic missile configuration is examined in a supersonic flow field to show the interaction of mesh adaptation and turbulence model. Turbulence model modifications for rotational correction and a quadratic constitutive relationship show a strong influence on adaptive mesh refinement and predicted rolling moment.

slender body aerodynamics↗

Robust Design Under Uncertainty of Hypersonic Inlets

The objective of this work is to optimize scramjet inlet designs under uncertainty at multiple fidelity levels. Optimizations are performed for both deterministic and stochastic conditions. For stochastic conditions, robust design is used to minimize the variance in system performance in addition to optimizing mean system performance. Uncertainty quantification is performed using non-intrusive polynomial chaos. Comparisons are made between low-fidelity analysis and high-fidelity analysis, as well as deterministic optimization and robust optimization. Results indicate that the general design trend of previous inlet optimization techniques is recovered for the external compression portion of an inlet. For the internal compression portion of an inlet, when viscous effects become more significant and boundary layer separation becomes more likely, the optimal trend departs from previous results. Multifidelity uncertainty quantification methods are found to substantially reduce the computation time of robust design while obtaining results that nearly match those of high-fidelity analysis. Robust optimization is found to decrease the standard deviation of throat Mach number by up to 40% compared to deterministic optimization. Changes in the design variables and total turning that lead to robust scramjet inlets are identified and linked back to fundamental physical principles through the theta-beta-Mach function.

Hypersonic↗

Midfidelity Model Verification for a Point-Absorbing Wave Energy Converter with Linear Power Takeoff: Preprint

In the preliminary design stage of a wave energy converter (WEC), fast and reliable simulation tools are required. High-fidelity numerical models are usually employed to study the wave-structure interaction, yet the computational cost is demanding. Instead, mid-fidelity models provide simulations in the order of real time. In this study, Uppsala University WEC operates in a relatively mild sea state and is modeled using WEC-Sim. The model is verified based on OpenFOAM simulations. To analyze the ability of mid-fidelity model to capture the WEC dynamics, the system is investigated separately in 1, 2 and 3 DoF. The contribution of viscous phenomena is examined. Moreover, linear and weakly non-linear solutions provided by WEC-Sim are studied. The results obtained indicate that the viscous effects in heave and surge motion can be neglected but not for the pitch. The weakly non-linear WEC-Sim solution successfully agrees with the CFD, whereas the linear solution could drive to misleading results.

CFD↗

Midfidelity Model Verification for a Point-Absorbing Wave Energy Converter with Linear Power Take-Off

In the preliminary design stage of a wave energy converter (WEC), fast and reliable simulation tools are required. High-fidelity numerical models are usually employed to study the wave-structure interaction, yet the computational cost is demanding. Instead, mid-fidelity models provide simulations in the order of real time. In this study, Uppsala University WEC operates in a relatively mild sea state and is modeled using WEC-Sim. The model is verified based on OpenFOAM simulations. To analyze the ability of mid-fidelity model to capture the WEC dynamics, the system is investigated separately in 1, 2 and 3 DoF. The contribution of viscous phenomena is examined. Moreover, linear and weakly non-linear solutions provided by WEC-Sim are studied. The results obtained indicate that the viscous effects in heave and surge motion can be neglected but not for the pitch. The weakly non-linear WEC-Sim solution successfully agrees with the CFD, whereas the linear solution could drive to misleading results.

CFD↗

An Efficient Three-Dimensional CFD-Based Numerical Wave Tank for a Wave Energy Converter in Extreme Irregular Waves: Preprint

A numerical wave tank approach for computational fluid dynamics (CFD) modelling of an extreme irregular seastate is presented. The technique couples a potential flow solution with a CFD solver for more efficient numerical predictions. This method has recently become attractive both for the research community and the industry working with offshore structures. The model is used to determine the response of a submerged pressure differential wave energy converter (WEC) in a fully nonlinear irregular wave condition using the high fidelity CFD code, STAR-CCM+. Potential flow based numerical models are commonly used to predict motions and performance of wave energy converters. Wave kinematics can deviate from potential flow predictions for extreme wave conditions; the excitation loads on an absorber can also be increasingly influenced by viscous effects, not captured in engineering level models. In these extreme conditions, a Reynolds-averaged Navier-Stokes CFD model can better predict motions and loads for a WEC. Long time series with varying random seed numbers can be used to identify singular extreme wave events from a stochastic irregular sea state. This approach simulates a more realistic wave series for a given sea state than a regular wave or a focused wave. However, it is computationally infeasible to run these long time series for three-dimensional (3D) CFD simulations. In this work, two-dimensional (2D) CFD simulations with a long domain allow the full development of an extreme nonlinear wave condition. The results are used to identify extreme events from a 50-year storm condition for the PacWave site off the coast of Oregon. A relatively short time window including this extreme event is then mapped to a 3D simulation using a user defined wave methodology. Convergence studies for domain length, wave forcing lengths, and time before the extreme event were conducted.

cfd↗