Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence methods”

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 793 records · Page 44

Fast laser-induced aerosol formation for visualization of gas flows

A technique for aerosol seeding of gas flows by laser-induced particle formation is demonstrated using a pulsed Nd:YAG laser (1.06 microns) for optical breakdown of a mixture of SF6 and H2 in an inert carrier gas. It is noted that, contrary to the smoke-wire approach, the laser-induced particles form first in zones of high turbulence, since mixing enhances coagulation. The method also allows seeding to be performed in locations hardly accessible otherwise and is mechanically nonintrusive. Finally, a study of the mixture and the breakdown effects indicates that for H2:SF6 ratios between 3:1 and 15:1 the particle formation is only limited by the physics of the gas/particle conversion.

Hassa, C.↗

Numerical simulation of compressible, turbulent, two-phase flow

A computer program for numerically simulating compressible, turbulent, two-phase flows is described and applied. Special attention is given to flows in which dust is ingested into the turbulent boundary layer behind shock waves moving over the earth's surface. it is assumed that the two phases are interpenetrating continua which are coupled by drag forces and heat transfer. The particle phase is assumed to be dilute, and turbulent effects are modeled by zero- and two-equation eddy viscosity models. An important feature of the turbulence modeling is the treatment of surface boundary conditions which control the ingestion of particles into the boundary layer by turbulent friction and diffusion. The numerical method uses second-order implicit upwind differencing of the inviscid terms of the equations and second-order central differencing of the viscous terms. A diagonal form of the implicit algorithm is used to improve efficiency, and the transformation to a curvilinear coordinate system is accomplished by the finite volume techniques. Applications to a series of representative flows include a two-phase nozzle flow, the steady flow of air over a sand bed, and the air flow behind a normal shock wave in uniform motion over a sand bed. Results of the latter two applications are compared with experimental results.

Coakley, t. J.↗

Large-scale computations in fluid mechanics; Proceedings of the Fifteenth Summer Seminar on Applied Mathematics, University of California, La Jolla, CA, June 27-July 8, 1983. Parts 1 & 2

Papers are presented on such topics as the use of semi-Lagrangian advective schemes in meteorological modeling; computation with high-resolution upwind schemes for hyperbolic equations; dynamics of flame propagation in a turbulent field; a modified finite element method for solving the incompressible Navier-Stokes equations; computational fusion magnetohydrodynamics; and a nonoscillatory shock capturing scheme using flux-limited dissipation. Consideration is also given to the use of spectral techniques in numerical weather prediction; numerical methods for the incorporation of mountains in atmospheric models; techniques for the numerical simulation of large-scale eddies in geophysical fluid dynamics; high-resolution TVD schemes using flux limiters; upwind-difference methods for aerodynamic problems governed by the Euler equations; and an MHD model of the earth's magnetosphere.

Engquist, B. E.↗

An improved viscid/inviscid interaction procedure for transonic flow over airfoils

A new interacting boundary layer approach for computing the viscous transonic flow over airfoils is described. The theory includes a complete treatment of viscous interaction effects induced by the wake and accounts for normal pressure gradient effects across the boundary layer near trailing edges. The method is based on systematic expansions of the full Reynolds equation of turbulent flow in the limit of Reynolds numbers, Reynolds infinity. Procedures are developed for incorporating the local trailing edge solution into the numerical solution of the coupled full potential and integral boundary layer equations. Although the theory is strictly applicable to airfoils with cusped or nearly cusped trailing edges and to turbulent boundary layers that remain fully attached to the airfoil surface, the method was successfully applied to more general airfoils and to flows with small separation zones. Comparisons of theoretical solutions with wind tunnel data indicate the present method can accurately predict the section characteristics of airfoils including the absolute levels of drag.

Melnik, R. E.↗

Liquid crystals for high-altitude in-flight boundary layer flow visualization

A novel method for visualization of laminar to turbulent boundary layer transition in flight is presented, in which liquid crystal coatings are used to indicate the transition by a change of color (in response to the changes in the wall shear and the temperature). The liquid crystals applied to the wing surfaces were either cholesteric crystals formulated to be active between temperatures of 30 and 100 F or a sample of chiral nematic liquid crystals formulated to be active between -18 and -22 F. Liquid crystal coatings were shown to provide visualization data at flight altitudes and temperatures at which data could not be obtained previously. Flight evaluations were successfully conducted on a T-34C and a Learjet research aircraft at Mach numbers of up to 0.8 and altitudes of up to 50,000 ft. The time response during oscillating sideslip tests was very rapid. In wind tunnel tests, a hysteresislike behavior was observed in the liquid crystal coatings, possibly because of the lower dynamic pressures in wind tunnels as compared to flight.

Gall, P. D.↗

Fuel-air mixing and combustion in a two-dimensional Wankel engine

A two-equation turbulence model, an algebraic grid generalization method, and an approximate factorization time-linearized numerical technique are used to study the effects of mixture stratification at the intake port and gaseous fuel injection on the flow field and fuel-air mixing in a two-dimensional rotary engine model. The fuel distribution in the combustion chamber is found to be a function of the air-fuel mixture fluctuations at the intake port. It is shown that the fuel is advected by the flow field induced by the rotor and is concentrated near the leading apex during the intake stroke, while during compression, the fuel concentration is highest near the trailing apex and is lowest near the rotor. It is also found that the fuel concentration near the trailing apex and rotor is small except at high injection velocities.

Shih, T. I.-P.↗

Vortex/boundary-layer interactions: Data report, volume 1

This report summarizes the work done under NASA Grant NAGw-581, Vortex/Boundary Layer Interactions. The experimental methods are discussed in detail and numerical results are presented, but are not fully interpreted. This report should be useful to anyone who wishes to make further use of the data (available on floppy disc or magnetic tape) for the development of turbulence models or the validation of predictive methods. Journal papers are in course of preparation.

Cutler, A. D.↗

Assessment of higher order turbulence models for complex two- and three-dimensional flowfields

A numerical method is presented to solve the three-dimensional Navier-Stokes equations in combination with a full Reynolds-stress turbulence model. Computations will be shown for three complex flowfields. The results of the Reynolds-stress model will be compared with those predicted by two different versions of the k-omega model. It will be shown that an improved version of the k-omega model gives as accurate results as the Reynolds-stress model.

Menter, Florian R.↗

Efficiency and reliability enhancements in propulsion flowfield modeling

The implementation of traditional CFD algorithms in practical propulsion related flowfields often leads to dramatic reductions in efficiency and/or robustness. The present research is directed at understanding the reasons for this deterioration and finding methods to circumvent it. Work to date has focussed on low Mach number regions, viscous dominated regions, and high grid aspect ratios. Time derivative preconditioning, improved definition of the local time stepping, and appropriate application of boundary conditions are employed to decrease the required time to obtain a solution, while maintaining accuracy. A number of cases having features typical of rocket engine flowfields are computed to demonstrate the improvement over conventional methods. These cases include laminar and turbulent high Reynolds number flat plate boundary layers, flow over a backward-facing step, a diffusion flame, and wall heat-flux calculations in a turbulent converging-diverging nozzle. Results from these cases show convergence that is virtually independent of the local Mach number and the grid aspect ratio, which translates to a convergence speed-up of up to several orders of magnitude over conventional algorithms. Current emphasis is in extending these results to three-dimensional flows with highly stretched grids.

Buelow, Philip E. O.↗

Parameterization and scaling of Arctic ice conditions in the context of ice-atmosphere processes

This report summarizes achievements during year three of our project to investigate the use of ERS-1 SAR data to study Arctic ice and ice/atmosphere processes. The project was granted a one year extension, and goals for the final year are outlined. The specific objects of the project are to determine how the development and evolution of open water/thin ice areas within the interior ice pack vary under different atmospheric synoptic regimes; compare how open water/thin ice fractions estimated from large-area divergence measurements differ from fractions determined by summing localized openings in the pack; relate these questions of scale and process to methods of observation, modeling, and averaging over time and space; determine whether SAR data might be used to calibrate ice concentration estimates from medium and low-rate bit sensors (AVHRR and DMSP-OLS) and the special sensor microwave imager (SSM/I); and investigate methods to integrate SAR data for turbulent heat flux parametrization at the atmosphere interface with other satellite data.

Barry, R. G.↗

Advanced Multi-phase Flow CFD Model Development for Solid Rocket Motor Flowfield Analysis

A Navier-Stokes code, finite difference Navier-Stokes (FDNS), is used to analyze the complicated internal flowfield of the SRM (solid rocket motor) to explore the impacts due to the effects of chemical reaction, particle dynamics, and slag accumulation on the solid rocket motor (SRM). The particulate multi-phase flowfield with chemical reaction, particle evaporation, combustion, breakup, and agglomeration models are included in present study to obtain a better understanding of the SRM design. Finite rate chemistry model is applied to simulate the chemical reaction effects. Hermsen correlation model is used for the combustion simulation. The evaporation model introduced by Spalding is utilized to include the heat transfer from the particulate phase to the gase phase due to the evaporation of the particles. A correlation of the minimum particle size for breakup expressed in terms of the Al/Al2O3 surface tension and shear force was employed to simulate the breakup of particles. It is assumed that the breakup occurs when the Weber number exceeds 6. A simple L agglomeration model is used to investigate the particle agglomeration. However, due to the large computer memory requirements for the agglomeration model, only 2D cases are tested with the agglomeration model. The VOF (Volume of Fluid) method is employed to simulate the slag buildup in the aft-end cavity of the redesigned solid rocket motor (RSRM). Monte Carlo method is employed to calculate the turbulent dispersion effect of the particles. The flowfield analysis obtained using the FDNS code in the present research with finite rate chemical reaction, particle evaporation, combustion, breakup, agglomeration, and VOG models will provide a design guide for the potential improvement of the SRM including the use of materials and the shape of nozzle geometry such that a better performance of the SRM can be achieved. The simulation of the slag buildup in the aft-end cavity can assist the designer to improve the design of the RSRM geometry.

Liaw, Paul↗

Development of a two-equation turbulence model for hypersonic flows: Evaluation of a low Reynolds number correction to the Kappa - epsilon two equation compressible turbulence model - Volume 1

The objective of the current research is the development of an improved k-epsilon two-equation compressible turbulence model for turbulent boundary layer flows experiencing strong viscous-inviscid interactions. The development of an improved model is important in the design of hypersonic vehicles such as the National Aerospace Plane (NASP) and the High Speed Civil Transport (HSCT). Improvements have been made to the low Reynolds number functions in the eddy viscosity and dissipation of solenoidal dissipation of the k-epsilon turbulence mode. These corrections offer easily applicable modifications that may be utilized for more complex geometries. The low Reynolds number corrections are functions of the turbulent Reynolds number and are therefore independent of the coordinate system. The proposed model offers advantages over some current models which are based upon the physical distance from the wall, that modify the constants of the standard model, or that make more corrections than are necessary to the governing equations. The code has been developed to solve the Favre averaged, boundary layer equations for mass, momentum, energy, turbulence kinetic energy, and dissipation of solenoidal dissipation using Keller's box scheme and the Newton spatial marching method. The code has been validated by removing the turbulent terms and comparing the solution with the Blasius solution, and by comparing the turbulent solution with an existing k-epsilon model code using wall function boundary conditions. Excellent agreement is seen between the computed solution and the Blasius solution, and between the two codes. The model has been tested for both subsonic and supersonic flat-plate turbulent boundary layer flow by comparing the computed skin friction with the Van Driest II theory and the experimental data of Weighardt; by comparing the transformed velocity profile with the data of Weighardt, and the Law of the Wall and the Law of the Wake; and by comparing the computed results of an adverse pressure gradient with the experimental data of Fernando and Smits. Good agreement is obtained with the experimental correlations for all flow conditions.

Knight, Doyle D.↗

Evaluation of Grid Modification Methods for On- and Off-Track Sonic Boom Analysis

Grid modification methods have been under development at NASA to enable better predictions of low boom pressure signatures from supersonic aircraft. As part of this effort, two new codes, Stretched and Sheared Grid - Modified (SSG) and Boom Grid (BG), have been developed in the past year. The CFD results from these codes have been compared with ones from the earlier grid modification codes Stretched and Sheared Grid (SSGRID) and Mach Cone Aligned Prism (MCAP) and also with the available experimental results. NASA's unstructured grid suite of software TetrUSS and the automatic sourcing code AUTOSRC were used for base grid generation and flow solutions. The BG method has been evaluated on three wind tunnel models. Pressure signatures have been obtained up to two body lengths below a Gulfstream aircraft wind tunnel model. Good agreement with the wind tunnel results have been obtained for both on-track and off-track (up to 53 degrees) cases. On-track pressure signatures up to ten body lengths below a Straight Line Segmented Leading Edge (SLSLE) wind tunnel model have been extracted. Good agreement with the wind tunnel results have been obtained. Pressure signatures have been obtained at 1.5 body lengths below a Lockheed Martin aircraft wind tunnel model. Good agreement with the wind tunnel results have been obtained for both on-track and off-track (up to 40 degrees) cases. Grid sensitivity studies have been carried out to investigate any grid size related issues. Methods have been evaluated for fully turbulent, mixed laminar/turbulent and fully laminar flow conditions.

Nayani, Sudheer N.↗

Measurements of few-mode fiber photonic lanterns in emulated atmospheric conditions for a low earth orbit space to ground optical communication receiver application

Photonic lanterns are being evaluated as a component of a scalable photon counting real-time optical ground receiver for space-to-ground photon-starved communication applications. The function of the lantern as a component of a receiver is to efficiently couple and deliver light from the atmospherically distorted focal spot formed behind a telescope to multiple small-core fiber-coupled single-element super-conducting nanowire detectors. This architecture solution is being compared to a multimode fiber coupled to a multi-element detector array. This paper presents a set of measurements that begins this comparison. This first set of measurements are a comparison of the throughput coupling loss at emulated atmospheric conditions for the case of a 60 cm diameter telescope receiving light from a low earth orbit satellite. The atmospheric conditions are numerically simulated at a range of turbulence levels using a beam propagation method and are physically emulated with a spatial light modulator. The results show that for the same number of output legs as the single-mode fiber lantern, the few-mode fiber lantern increases the power throughput up to 3.92 dB at the worst emulated atmospheric conditions tested of D/r(sub 0)=8.6. Furthermore, the coupling loss of the few-mode fiber lantern approaches the capability of a 30 micron graded index multimode fiber chosen for coupling to a 16 element detector array.

Tedder, Sarah A.↗

Measurements of Few-Mode Fiber Photonic Lanterns in Emulated Atmospheric Conditions for a Low Earth Orbit Space to Ground Optical Communication Receiver Application

Photonic lanterns are being evaluated as a component of a scalable photon counting real-time optical ground receiver for space-to-ground photon-starved communication applications. The function of the lantern as a component of a receiver is to efficiently couple and deliver light from the atmospherically distorted focal spot formed behind a telescope to multiple small-core fiber-coupled single-element super-conducting nanowire detectors. This architecture solution is being compared to a multimode fiber coupled to a multi-element detector array. This paper presents a set of measurements that begins this comparison. This first set of measurements are a comparison of the throughput coupling loss at emulated atmospheric conditions for the case of a 60 cm diameter telescope receiving light from a low earth orbit satellite. The atmospheric conditions are numerically simulated at a range of turbulence levels using a beam propagation method and are physically emulated with a spatial light modulator. The results show that for the same number of output legs as the single-mode fiber lantern, the few mode fiber lantern increases the power throughput up to 3.92 dB at the worst emulated atmospheric conditions tested of D/r0=8.6. Furthermore, the coupling loss of the few mode fiber lantern approaches the capability of a 30 micron graded index multimode fiber chosen for coupling to a 16 element detector array.

Tedder, Sarah A.↗

Wall-Modeled Large Eddy Simulations of Transonic Buffet Over a Supercritical Airfoil

A series of scale-resolving simulations of flow over the ONERA OAT15A airfoil have been performed at an angle of attack of 3.5◦, just past the onset of buffet. The focus of this study is to document the sensitivity of the wall-modeled large eddy simulation (WMLES) methodology for curvilinear structured overset grids within the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework to mesh spacing, mesh distribution, and domain size. A secondary purpose of the study is to compare the results from WMLES to unsteady Reynolds-averaged Navier Stokes (URANS) simulations and hybrid RANS-LES (HRLES) within the same LAVA solver framework. The study provides a unique perspective regarding comparisons between different turbulence modeling approaches, time-integration methods, and computational performance since many of the same numerical routines are used for all three types of simulations. The results are compared with experiments and previous numerical studies of the same geometry and flow conditions.

TTT↗

High-resolution turbulent simulations using the Connection Machine-2

The spectral method provides an efficient algorithm for solving the 3D incompressible Navier-Stokes equations in periodic boundaries. Most people, so far, have used vectorized machines, such as the CRAY-2, to implement fast Fourier transformations and time integrations in the spectral calculations. In this paper, new results are presented using the spectral calculations on the Connection Machine-2 with a parallel algorithm. The large memory of the Connection Machine-2 and the parallel algorithm allows, of the first time, to implement a 512-cubed mesh resolution for high Reynolds number flows. The computational speed of the present code is about 30 percent faster than the fastest CRAY-2 simulations with four processors. Parallel machines, such as the Connection Machine-2, will possibly provide new computational power for understanding the intermittency and cascade mechanism in fluid turbulence.

Chen, Shiyi↗

Large eddy simulation of incompressible turbulent channel flow

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Moin, P.↗