Search NASA⌕ Search

SEARCH · Search NASA

Results for “wake function”

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 253 records · Page 14

Calculation of wall and free turbulent-shear flows at supersonic speeds

Supersonic turbulent flows are simulated numerically by solving the Reynolds-averaged full Navier-Stokes equations by an implicit finite-volume method. This flux-vector-split upwind scheme uses approximate factorization followed by line-Gauss-Seidel relaxations. The effects of turbulence are modeled by two eddy viscosity models. High-Reynolds-number form of the k-epsilon model is coupled with a wall-function to avoid excessive refinement of the grid in the low-Reynolds number regions. The k-epsilon equations are solved by the explicit-implicit MacCormack method. The algebraic Baldwin-Lomax model is also being used as an inexpensive alternative for the cases which do not experience massive separation. Several examples of two-dimensional solutions are given to illustrate both wall and free turbulent shear flows which include fluid dynamic phenomena, such as shocks, boundary layers, shear layers, wakes, separations and recirculations. The results compared with experimental data show good agreeent.

Baysal, O.↗

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗

Meteoroid Environment Modeling: the Meteoroid Engineering Model and Shower Forecasting

INTRODUCTION: The meteoroid environment is often divided conceptually into meteor showers and the sporadic meteor background. It is commonly but incorrectly assumed that meteoroid impacts primarily occur during meteor showers; instead, the vast majority of hazardous meteoroids belong to the sporadic complex. Unlike meteor showers, which persist for a few hours to a few weeks, sporadic meteoroids impact the Earth's atmosphere and spacecraft throughout the year. The Meteoroid Environment Office (MEO) has produced two environment models to handle these cases: the Meteoroid Engineering Model (MEM) and an annual meteor shower forecast. The sporadic complex, despite its year-round activity, is not isotropic in its directionality. Instead, their apparent points of origin, or radiants, are organized into groups called "sources". The speed, directionality, and size distribution of these sporadic sources are modeled by the Meteoroid Engineering Model (MEM), which is currently in its second major release version (MEMR2) [Moorhead et al., 2015]. MEM provides the meteoroid flux relative to a user-provided spacecraft trajectory; it provides the total flux as well as the flux per angular bin, speed interval, and on specific surfaces (ram, wake, etc.). Because the sporadic complex dominates the meteoroid flux, MEM is the most appropriate model to use in spacecraft design. Although showers make up a small fraction of the meteoroid environment, they can produce significant short-term enhancements of the meteoroid flux. Thus, it can be valuable to consider showers when assessing risks associated with vehicle operations that are brief in duration. To assist with such assessments, the MEO issues an annual forecast that reports meteor shower fluxes as a function of time and compares showers with the time-averaged total meteoroid flux. This permits missions to do quick assessments of the increase in risk posed by meteor showers.

Moorhead, Althea V.↗

Meteoroid Environment Modeling: The Meteoroid Engineering Model and Shower Forecasting

The meteoroid environment is often divided conceptually into meteor showers and the sporadic meteor background. It is commonly but incorrectly assumed that meteoroid impacts primarily occur during meteor showers; instead, the vast majority of hazardous meteoroids belong to the sporadic complex. Unlike meteor showers, which persist for a few hours to a few weeks, sporadic meteoroids impact the Earth's atmosphere and spacecraft throughout the year. The Meteoroid Environment Office (MEO) has produced two environment models to handle these cases: the Meteoroid Engineering Model (MEM) and an annual meteor shower forecast. The sporadic complex, despite its year-round activity, is not isotropic in its directionality. Instead, their apparent points of origin, or radiants, are organized into groups called "sources". The speed, directionality, and size distribution of these sporadic sources are modeled by the Meteoroid Engineering Model (MEM), which is currently in its second major release version (MEMR2) [Moorhead et al., 2015]. MEM provides the meteoroid flux relative to a user-provided spacecraft trajectory; it provides the total flux as well as the flux per angular bin, speed interval, and on specific surfaces (ram, wake, etc.). Because the sporadic complex dominates the meteoroid flux, MEM is the most appropriate model to use in spacecraft design. Although showers make up a small fraction of the meteoroid environment, they can produce significant short-term enhancements of the meteoroid flux. Thus, it can be valuable to consider showers when assessing risks associated with vehicle operations that are brief in duration. To assist with such assessments, the MEO issues an annual forecast that reports meteor shower fluxes as a function of time and compares showers with the time-averaged total meteoroid flux. This permits missions to do quick assessments of the increase in risk posed by meteor showers. Section II describes MEM in more detail and describes our current efforts to improve its characteristics for a future release. Section III describes the annual shower forecast and highlights recent improvements made to its algorithm and inputs.

Moorhead, Althea V.↗

Upgraded viscous flow analysis of multi-element airfoils

A description of an improved version of the NASA/Lockheed multi-element airfoil analysis computer program is presented. The improvements include several major modifications of the aerodynamic model as well as substantial changes of the computer code. The modifications of the aerodynamic model comprise the representation of the boundary layer and wake displacement effects with an equivalent source distribution, the prediction of wake parameters with Green's lag-entrainment method, the calculation of turbulent boundary layer separation with the method of Nash and Hicks, the estimation of the onset of confluent boundary layer separation with a modified form of Goradia's method, and the prediction of profile drag with the formula of Squire and Young. The modifications of the computer program for which the structured approach to computer software development was employed are also described. Important aspects of the structured program development such as the functional decomposition of the aerodynamic theory and its numerical implementation, the analysis of the data flow within the code, and the application of a pseudo code are discussed.

Brune, G. W.↗

Laser Doppler velocimeter system simulation for sensing aircraft wake vortices

A hydrodynamic model of aircraft vortex wakes in an irregular wind shear field near the ground is developed and used as a basis for modeling the characteristics of a laser Doppler detection and vortex location system. The trailing vortex sheet and the wind shear are represented by discrete free vortices distributed over a two-dimensional grid. The time dependent hydrodynamic equations are solved by direct numerical integration in the Boussinesq approximation. The ground boundary is simulated by images, and fast Fourier Transform techniques are used to evaluate the vorticity stream function. The atmospheric turbulence was simulated by constructing specific realizations at time equal to zero, assuming that Kolmogoroff's law applies, and that the dissipation rate is constant throughout the flow field. The response of a simulated laser Doppler velocimeter is analyzed by simulating the signal return from the flow field as sensed by a simulation of the optical/electronic system.

Thomson, J. A. L.↗

Upgraded viscous flow analysis of multielement airfoils

A description of an improved version of the NASA/Lockheed multielement airfoil analysis computer program is presented. The improvements include several major modifications of the aerodynamic model as well as substantial changes of the computer code. The modifications of the aerodynamic model comprise the representation of the boundary layer and wake displacement effects with an equivalent source distribution, the prediction of wake parameters with Green's lag-entrainment method, the calculation of turbulent boundary layer separation with the method of Nash and Hicks, the estimation of the onset of confluent boundary layer separation with a modified form of Goradia's method, and the prediction of profile drag with the formula of Squire and Young. The paper further describes the modifications of the computer program for which the structured approach to computer software development was employed. Important aspects of the structured program development such as the functional decomposition of the aerodynamic theory and its numerical implementation, the analysis of the data flow within the code, and the application of a pseudo code are discussed.

Brune, G. W.↗

Nitric Oxide Formation by Meteoroids in the Upper Atmosphere

The process of nitric oxide formation during atmospheric entry of meteoroids is analyzed theoretically. An ablating meteoroid is assumed to be a point source in a uniform flow with a continuum regime evolving in its wake. The amount of nitric oxide produced by high-temperature reactions of air in the continuum regime is calculated by numerical integration of chemical-rate equations. This is accomplished by assuming that flow properties are constant across the reacting region, the radius of the region being determined from considerations of shock-wave formation and molecular diffusion. The results, when summed over the observed mass, velocity, and entry-angle distributions of meteoroids, provide annual global production rates of nitric oxide as a function of altitude. The peak production of nitric oxide is found to occur at altitudes between 9 x 10(exp 4) and 10(exp 5) m, the total annual rate being about 4 x 10(exp 7) kg. The present results suggest that the large concentration of nitric oxide observed below 9.5 x 10(exp 4) m could be attributed to meteoroids instead of photodissociation of nitrogen into metastable, 2D-state atoms, as has been previously hypothesized.

Menees, Gene P.↗

Wake Cycle Robustness of the Mars Science Laboratory Flight Software

The Mars Science Laboratory (MSL) is a spacecraft being developed by the Jet Propulsion Laboratory (JPL) for the purpose of in-situ exploration on the surface of Mars. The objective of MSL is to explore and quantitatively assess a local region on the Martian surface as a habitat for microbial life, past or present. This objective will be accomplished through the assessment of the biological potential of at least one target environment, the characterization of the geology and geochemistry of the landing region, an investigation of the planetary process relevant to past habitability, and a characterization of surface radiation. For this purpose, MSL incorporates a total of ten scientific instruments for which functions are to include, among others, atmospheric and descent imaging, chemical composition analysis, and radiation measurement. The Flight Software (FSW) system is responsible for all mission phases, including launch, cruise, entry-descent-landing, and surface operation of the rover. Because of the essential nature of flight software to project success, each of the software modules is undergoing extensive testing to identify and correct errors.

Whitehill, Robert↗

An algebraic turbulence model for three-dimensional viscous flows

An algebraic turbulence model is proposed for use with three-dimensional Navier-Stokes analyses. It incorporates features of both the Baldwin-Lomax and Cebeci-Smith models. The Baldwin-Lomax model uses the maximum of a function f(y) to determine length and velocity scales. An analysis of the Baldwin-Lomax model shows that f(y) can have a spurious maximum close to the wall, causing numerical problems and non-physical results. The proposed model uses integral relations to determine delta(*) u(sub e) and delta used in the Cebeci-Smith mode. It eliminates a constant in the Baldwin-Lomax model and determines the two remaining constants by comparison to the Cebeci-Smith formulation. Pressure gradient effects, a new wake model, and the implementation of these features in a three-dimensional Navier-Stokes code are also described. Results are shown for a flat plate boundary layer, an annular turbine cascade, and endwall heat transfer in a linear turbine cascade. The heat transfer results agree well with experimental data which shows large variations in endwall Stanton number contours with Reynolds number.

Chima, R. V.↗

An algebraic turbulence model for three-dimensional viscous flows

An algebraic turbulence model is proposed for use with three-dimensional Navier-Stokes analyses. It incorporates features of both the Baldwin-Lomax and Cebeci-Smith models. The Baldwin-Lomax model uses the maximum of a function f(v) to determine length and velocity scales. An analysis of the Baldwin-Lomax model shows that f(v) can have a spurious maximum close to the wall, causing numerical problems and non-physical results. The proposed model uses integral relations to determine delta(*) u(sub e) and delta used in the Cebeci-Smith mode. It eliminates a constant in the Baldwin-Lomax model and determines the two remaining constants by comparison to the Cebeci-Smith formulation. Pressure gradient effects, a new wake model, and the implementation of these features in a three-dimensional Navier-Stokes code are also described. Results are shown for a flat plate boundary layer, an annular turbine cascade, and endwall heat transfer in a linear turbine cascade. The heat transfer results agree well with experimental data which shows large variations in endwall Stanton number contours with Reynolds number.

Chima, R. V.↗

Extragalactic Radio Jets and Intergalactic Medium 2

During the mid 1990s, the ROSAT satellite provided the highest resolution X-ray imaging, with useful X-ray spectroscopy, and excellent sensitivity. ROSAT"s PSPC and HRI instruments gave us our first high-quality data on the AGN phenomenon in nearby galaxies - the central theme of this research project. Initially this project concentrated on separating the AGN-related component of the X-ray emission from thermal radiation from the surrounding atmospheres in radio galaxies. There was ample reason to believe that this separation would be possible, although earlier work had taken the view that the X-ray emission from radio galaxies is either wholly AGN-related (and hence of sub-arcsec scale), or wholly from the galaxy and cluster atmospheres (and hence of scales - 10 arcsec or more in the low-redshift radio galaxies that we chose to study). First with the PSPC, and then with the HRI, we proved that a wide range of AGN core X ray emissivities could be found, and that the generic radio galaxy produces both AGN-related and atmosphere-related X-radiation. We demonstrated that there is a close relationship between the core X-ray power and the core radio power, and found instances in which the hot atmospheres of the galaxies must be participating in cooling flows, and other cases where those atmospheres are relatively stable. In some cases, there is a clear active relationship between the extended radio emission and the X-ray atmospheres: thus in NGC 326, we interpret the distorted radio structure as evidence of buoyancy as the radio plasma rises through the cluster atmosphere. In 3C 449, we can see that the radio plasma is displacing the X-ray emitting gas. In a further case (NGC 1265), we hoped that the motion of a radio galaxy through a cluster atmosphere would be apparent though the wake that it might establish - unfortunately, other structure in the atmosphere of the Perseus cluster tends to dominate our X-ray image, and no useful limits could be placed on the strength of the wake. In the light of so-called "unified" models of AGN, it is expected that the low-power radio galaxies studied under these programmes are the unbeamed counterparts of BL Lac objects. Accordingly, we predicted that BL Lac objects, as a class, should show X-ray halos from atmospheres similar to those associated with radio galaxies. The difficulty with studying this is that BL Lac objects tend to be X-ray bright, and the halos then vanish under the wings of the point response function. We found a BL Lac with less of a dynamic range problem, mapped it with the ROSAT HRI, and duly found an atmosphere (Hardcastle et al. 1999). Unfortunately for unification models, this atmosphere is exceptionally dense, so that it must be participating in an intense cooling flow. The implications of this for unification models are far from clear, and further work is under way.

Birkinshaw, Mark↗

Cassini Thermal Observations of Saturn's Main Rings: Implications for Particle Rotation and Vertical Mixing

In late 2004 and 2005 the Cassini composite infrared spectrometer (CIRS) obtained spatially resolved thermal infrared radial scans of Saturn's main rings (A, B and C, and Cassini Division) that show ring temperatures decreasing with increasing solar phase angle, (alpha), on both the lit and unlit faces of the ring plane. These temperature differences suggest that Saturn's main rings include a population of ring particles that spin slowly, with a spin period greater than 3.6 h, given their low thermal inertia. The A ring shows the smallest temperature variation with (alpha), and this variation decreases with distance from the planet. This suggests an increasing number of smaller, and/or more rapidly rotating ring particles with more uniform temperatures, resulting perhaps from stirring by the density waves in the outer A ring and/or self-gravity wakes. The temperatures of the A and B rings are correlated with their optical depth, (tau), when viewed from the lit face, and anti-correlated when viewed from the unlit face. On the unlit face of the B ring, not only do the lowest temperatures correlate with the largest (tau), these temperatures are also the same at both low and high a, suggesting that little sunlight is penetrating these regions. The temperature differential from the lit to the unlit side of the rings is a strong, nearly linear, function of optical depth. This is consistent with the expectation that little sunlight penetrates to the dark side of the densest rings, but also suggests that little vertical mixing of ring particles is taking place in the A and B rings.

Cassini composite infrared spectrometer (CIRS)↗

Crack Growth Behavior in the Threshold Region for High Cycle Fatigue Loading

The research results described in this paper presents a new understanding of the behavior of fatigue crack growth in the threshold region. It is believed by some crack growth experts that the ASTM load shedding test method does not produce true or valid threshold properties. The concern involves the observed fanning of threshold region da/dN data plots for some materials in which the low R-ratio data fans out or away from the high R-ratio data. This data fanning or elevation of threshold values is obviously caused by an increase in crack closure in the low R-ratio tested specimens. This increase in crack closure is assumed by some investigators to be caused by a plastic wake on the crack surfaces that was created during the load shedding test phase. This study shows that the increase in crack closure is the result of an extensive occurrence of crack bifurcation behavior in some materials, particularly in aluminum alloys, when the crack tip cyclic yield zone size becomes less than the grain size of the alloy. This behavior is related to the high stacking fault energy (SFE) property of aluminum alloys which results in easier slip characteristics. Therefore, the particular fanning behavior in aluminum alloys is a function of intrinsic dislocation property of the materials and that the fanned data represents valid material properties. However, for corrosion sensitive steel alloys used in this study the fanning was caused by a build-up of iron oxide at the crack tip from fretting corrosion.

Forman, Royce G.↗

Numerical simulation of spacecraft charging by impact-induced plasmas during a cometary flyby

A numerical model is developed for the interaction of a cometary probe, such as Giotto, with its environment, i.e., dust particles and gas. The spacecraft was set on a course to pass the comet at a velocity of 69 km/sec, so a chance existed that a potential field would form around the spacecraft and block lower energy particles from reaching the spacecraft instruments. The motion of electrons and ions is traced as a function of time to examine the evolution of the electric field, electric potential and the total space charge distributions on the surface of the spacecraft and its environment. Account is taken of the density of the particles and gas molecules at various distances from the comet, the collision energies involved, and the Giotto geometry. A solution is defined for the Poisson equation to describe the evolution of the plasma around Giotto, including the effects of ion collisions with the Al bumper protecting the spacecraft. The simulation predicts formation of an ion wake behind Giotto and the evolution of a positive potential on the order of 10 V around the spacecraft, i.e., sufficient for a positive potential barrier near the surface of the spacecraft.

Thiemann, H.↗

The Effect of Trailing Edge Grid Resolution in Large-Eddy Simulations

Ceramic Matrix Composite (CMC) with environmental barrier coating can tolerate significantly higher temperatures compared with the conventional metal alloy blades and thus can help achieve higher thermal efficiency in gas turbine engines. Due to the complex fabrication process of CMC blades, larger blade trailing edge thicknesses and larger leading-edge radii of curvature are expected. These features significantly alter aerodynamic performance especially near the trailing edge and thus the loss profile. In previous studies, the authors used the large eddy simulation (LES) as well as the Reynolds Averaged Navier-Stokes (RANS) coupled with an intermittency function-based transition model to simulate an aerodynamics performance of three different CMC blades and compared the predictions against the data acquired at NASA Glenn Transonic Turbine Blade Cascade Rig. Although we were able to successfully predict different characteristics of aerodynamic performance among three different blades (e.g., pressure loading and integrated losses), LES consistently predicted a small recirculation bubble at the end of the suction side. This was not inconsistent with the physics of the flow, but the pressure loading, locally, near the trailing edge, and as a result, the measured wake profiles did not match the experimental data. To delve further into these issues, a new LES grid with a much finer mesh near the trailing edge was generated where the spanwise extent of the grid was increased by 50%. An LES simulation for the CMC blade with the largest trailing edge thicknesses at the high-freestream turbulence (Tu=13%) condition showed that using such a mesh is necessary for LES to accurately capture the transition phenomena on the suction side. The resulting total pressure coefficient profile and the shift of the wake profile were greatly improved.

CMC↗

The Effect of Trailing Edge Grid Resolution in Large-Eddy Simulations

Ceramic Matrix Composite (CMC) with environmental barrier coating can tolerate significantly higher temperatures compared with the conventional metal alloy blades and thus can help achieve higher thermal efficiency in gas turbine engines. Due to the complex fabrication process of CMC blades, larger blade trailing edge thicknesses and larger leading-edge radii of curvature are expected. These features significantly alter aerodynamic performance especially near the trailing edge and thus the loss profile. In previous studies, the authors used the large eddy simulation (LES) as well as the Reynolds Averaged Navier-Stokes (RANS) coupled with an intermittency function-based transition model to simulate an aerodynamics performance of three different CMC blades and compared the predictions against the data acquired at NASA Glenn Transonic Turbine Blade Cascade Rig. Although we were able to successfully predict different characteristics of aerodynamic performance among three different blades (e.g., pressure loading and integrated losses), LES consistently predicted a small recirculation bubble at the end of the suction side. This was not inconsistent with the physics of the flow, but the pressure loading, locally, near the trailing edge, and as a result, the measured wake profiles did not match the experimental data. To delve further into these issues, a new LES grid with a much finer mesh near the trailing edge was generated where the spanwise extent of the grid was increased by 50%. An LES simulation for the CMC blade with the largest trailing edge thicknesses at the high-freestream turbulence (Tu=13%) condition showed that using such a mesh is necessary for LES to accurately capture the transition phenomena on the suction side. The resulting total pressure coefficient profile and the shift of the wake profile were greatly improved.

LES↗

The Design of Pumpjets for Hydrodynamic Propulsion

A procedure for use in the design of a wake adapted pumpjet mounted on the aft end of a body of revolution is presented. To this end, a pumpjet is designed for the Akron airship. The propulsor mass flow is selected to minimize kinetic energy losses through the duct and in the discharge jet. The shaft speed and disk size are selected to satisfy specified limits of cavitation performance and to provide acceptable blade loading. The streamtubes which pass through a propulsor mounted on a tapered afterbody follow essentially conical surfaces. A method is provided for defining these surfaces as a function of shroud geometry, rotor head distribution, and the energy distribution of the ingested mass flow. The three-dimensional effects to which the conical flow subjects the cylindrical blade design sections are described and a technique is presented which permits incorporation of these effects in the blade design procedure.

Bruce, E. P.↗