Search NASA⌕ Search

SEARCH · Search NASA

Results for “data flow statistics”

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 217 records · Page 12

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.↗

Lagrangian and Eulerian statistics obtained from direct numerical simulations of homogeneous turbulence

Direct numerical simulation is used to study dispersion in decaying isotropic turbulence and homogeneous shear flow. Both Lagrangian and Eulerian data are presented allowing direct comparison, but at fairly low Reynolds number. The quantities presented include properties of the dispersion tensor, isoprobability contours of particle displacement, Lagrangian and Eulerian velocity autocorrelations and time scale ratios, and the eddy diffusivity tensor. The Lagrangian time microscale is found to be consistently larger than the Eulerian microscale, presumably due to the advection of the small scales by the large scales in the Eulerian reference frame.

Squires, Kyle D.↗

Roles of Engineering Correlations in Hypersonic Entry Boundary Layer Transition Prediction

Efforts to design and operate hypersonic entry vehicles are constrained by many considerations that involve all aspects of an entry vehicle system. One of the more significant physical phenomenon that affect entry trajectory and thermal protection system design is the occurrence of boundary layer transition from a laminar to turbulent state. During the Space Shuttle Return To Flight activity following the loss of Columbia and her crew of seven, NASA's entry aerothermodynamics community implemented an engineering correlation based framework for the prediction of boundary layer transition on the Orbiter. The methodology for this implementation relies upon similar correlation techniques that have been is use for several decades. What makes the Orbiter boundary layer transition correlation implementation unique is that a statistically significant data set was acquired in multiple ground test facilities, flight data exists to assist in establishing a better correlation and the framework was founded upon state of the art chemical nonequilibrium Navier Stokes flow field simulations. Recent entry flight testing performed with the Orbiter Discovery now provides a means to validate this engineering correlation approach to higher confidence. These results only serve to reinforce the essential role that engineering correlations currently exercise in the design and operation of entry vehicles. The framework of information related to the Orbiter empirical boundary layer transition prediction capability will be utilized to establish a fresh perspective on this role, and to discuss the characteristics which are desirable in a next generation advancement. The details of the paper will review the experimental facilities and techniques that were utilized to perform the implementation of the Orbiter RTF BLT Vsn 2 prediction capability. Statistically significant results for multiple engineering correlations from a ground testing campaign will be reviewed in order to describe why only certain correlations were selected for complete implementation to support the Shuttle Program. Historical Orbiter flight data on early boundary layer transition due to protruding gap fillers will be described in relation to the selected empirical correlations. In addition, Orbiter entry flight testing results from the BLT Flight Experiment will be discussed in relation to these correlations. Applicability of such correlations to the entry design problem will be reviewed, and finally a perspective on the desirable characteristics for a next generation capability based on high fidelity physical models will be provided.

Campbell, Charles H.↗

Ensemble Data Assimilation Without Ensembles: Methodology and Application to Ocean Data Assimilation

Two methods to estimate background error covariances for data assimilation are introduced. While both share properties with the ensemble Kalman filter (EnKF), they differ from it in that they do not require the integration of multiple model trajectories. Instead, all the necessary covariance information is obtained from a single model integration. The first method is referred-to as SAFE (Space Adaptive Forecast error Estimation) because it estimates error covariances from the spatial distribution of model variables within a single state vector. It can thus be thought of as sampling an ensemble in space. The second method, named FAST (Flow Adaptive error Statistics from a Time series), constructs an ensemble sampled from a moving window along a model trajectory. The underlying assumption in these methods is that forecast errors in data assimilation are primarily phase errors in space and/or time.

Data Assimilation↗

Statistical uncertainties of the v n { 2 k } harmonics from Q cumulants

Analytic formulas to calculate statistical uncertainties of v n { 2 k } harmonics extracted from the Q cumulants are presented. The Q cumulants are multivariate polynomial functions of the weighted means of 2m-particle azimuthal correlations, $\langle\langle 2 m \rangle\rangle$. Variances and covariances of the $\langle\langle 2 m \rangle\rangle$ are included in the analytic formulas of the uncertainties that can be calculated simultaneously with the calculations of the v n { 2 k } harmonics. The calculations are performed using a simple toy model, which roughly simulates elliptic flow azimuthal anisotropy with magnitudes around 0.05. The results are compared with the results obtained by the many data subsets, and by the bootstrapping method. The first one is a common way of estimation of the statistical uncertainties of the v n { 2 k } harmonics in a real experiment. In order to increase precision in the measurement of the v n { 2 k } harmonics, a large number of 15 000 subsets and the same number of the resampling in the bootstrap method is used. Unlike the other ways of the analytic calculation of the v n { 2 k } statistical uncertainties, our proposal that includes the use of squared weights in the calculation of both the variances and covariances, gives the best agreement with the results obtained from the subsets and bootstrap method. Additionally, a recurrence equation between Q cumulants of any order is also presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Role of Declining Snow Cover in the Desiccation of the Great Salt Lake, Utah, using MODIS Data

The Great Salt Lake (GSL) in Utah has been shrinking since the middle of the 19th Century, leading to decreased area and volume, and increased salinity. We use satellite data products from the Terra and Aqua MODerate-resolution Imaging Spectroradiometer (MODIS) and the Landsat-7 and -8 satellites, along with meteorological and streamflow data, and modeled data products to study the relationship between changing snow-cover conditions and the decline of the GSL since 2000 in the context of the historical record of lake levels. The GSL basin includes much of the snow-dominated Wasatch and Uinta mountain ranges to the east of the lake. Snowmelt feeds the Bear, Jordan, and Weber rivers which are the three main rivers that flow into the lake. Snowmelt-timing maps, derived from a new MODIS standard snow-cover product, MOD10A1F, show that snow melted ~9.5 days earlier in the GSL basin during the study period, extending from 2000 – 2018. Air temperatures derived from 26 meteorological stations and surface temperatures measured by the Aqua MODIS land-surface temperature (LST) products, MYD21A1D and MYD21A1N, show trends of increasing temperature of ~0.94°C (a=0.05), and ~2.18°C, respectively, with most of the LST trends in the GSL basin being statistically significant (a=0.05). Increasing air temperatures in the basin have led to less precipitation falling as snow, lower snow depth (by ~34.5 mm (=0.01)) and snow-water equivalent (0.02 mm (a=0.01)), and earlier snowmelt. Also during the study period, Global Land surface Evaporation Amsterdam Model data show evaporation increasing by ~3.2 mm/yr, with trends in much of the basin being statistically significant (a=0.05). Trends calculated from the various products are generally in agreement indicating higher temperatures, greater evaporation, less snowfall and snow-on-the ground, and earlier snowmelt. Earlier snowmelt contributes to increasing evaporative loss from water flowing toward the lake. Furthermore, a lower mountain snowpack and less precipitation falling as snow (versus rain) is associated with lower stream discharge even if overall precipitation stays the same. The surface-water temperature of the GSL also increased over the study period by ~ 0.69°C, according to the MODIS LST data products, and the surface-water elevation of the lake dropped by ~1.7 m between 2000 and 2018 based on United States Geological Survey measurements, and the areal extent of the lake decreased by ~901 km2 as measured using Landsat imagery. Desiccation of the lake is associated with deleterious effects on wildlife, recreational activities, and some local industries. And, importantly, an expanding lake bed can also fuel dust storms that promote dangerous air quality along the Wasatch Front. This work elucidates the key role that satellite remote sensing can play in documenting earlier snowmelt and other changes in the GSL basin that influence the ongoing decline of the Great Salt Lake.

Dorothy K Hall↗

Modeling flow and transport in fracture networks using graphs

Fractures form the main pathways for flow in the subsurface within low-permeability rock. For this reason, accurately predicting flow and transport in fractured systems is vital for improving the performance of subsurface applications. Fracture sizes in these systems can range from millimeters to kilometers. Although modeling flow and transport using the discrete fracture network (DFN) approach is known to be more accurate due to incorporation of the detailed fracture network structure over continuum-based methods, capturing the flow and transport in such a wide range of scales is still computationally intractable. Furthermore, if one has to quantify uncertainty, hundreds of realizations of these DFN models have to be run. To reduce the computational burden, we solve flow and transport on a graph representation of a DFN. We study the accuracy of the graph approach by comparing breakthrough times and tracer particle statistical data between the graph-based and the high-fidelity DFN approaches, for fracture networks with varying number of fractures and degree of heterogeneity. Due to our recent developments in capabilities to perform DFN high-fidelity simulations on fracture networks with large number of fractures, we are in a unique position to perform such a comparison. We show that the graph approach shows a consistent bias with up to an order of magnitude slower breakthrough when compared to the DFN approach. We show that this is due to graph algorithm's underprediction of the pressure gradients across intersections on a given fracture, leading to slower tracer particle speeds between intersections and longer travel times. We present a bias correction methodology to the graph algorithm that reduces the discrepancy between the DFN and graph predictions. We show that with this bias correction, the graph algorithm predictions significantly improve and the results are very accurate. In conclusion, the good accuracy and the low computational cost, with O(10 4 ) times lower times than the DFN, makes the graph algorithm an ideal technique to incorporate in uncertainty quantification methods.

58 GEOSCIENCES↗

GEOS-3 radar altimeter study for the South Atlantic Bight

Three years of radar altimeter data from GEOS-3 for the South Atlantic Bight were processed. Mean monthly topographic maps were produced which estimate geostrophic flow in the region. Statistical distribution of the surface wind speed and significant wave height as a function of both space and time are presented.

Leitao, C. D.↗

Thermal regime of the continental lithosphere

From studies of the global heat flow data set, it has been generalized, with respect to the continental lithosphere, that there is a negative correlation between heat flow and the lithosphere's tectonic edge, and that the lithosphere's thermal evolution is similar to that of the ocean basins, resulting in a 'stable geotherm' in both environments. It is presently noted that a regional study perspective for heat flow data leads to doubts concerning the general applicability of either statement. Rao et al. (1982) have demonstrated that the data are not normally distributed, and that it is not possible to establish a negative correlation between heat flow and age in a rigorous statistical fashion. While some sites of stable continental blocks may have a geotherm that is by chance similar to that for old ocean basins, this need not hold true generally, and many stable continental terranes will be characterized by geotherms very different from those for old ocean basins.

Morgan, P.↗

Transport of passive scalars in a turbulent channel flow

A direct numerical simulation of a turbulent channel flow with three passive scalars at different molecular Prandtl numbers is performed. Computed statistics including the turbulent Prandtl numbers are compared with existing experimental data. The computed fields are also examined to investigate the spatial structure of the scalar fields. The scalar fields are highly correlated with the streamwise velocity; the correlation coefficient between the temperature and the streamwise velocity is as high as 0.95 in the wall region. The joint probability distributions between the temperature and velocity fluctuations are also examined; they suggest that it might be possible to model the scalar fluxes in the wall region in a manner similar to the Reynolds stresses.

Kim, John↗

Transport of passive scalars in a turbulent channel flow

A direct numerical simulation of a turbulent channel flow with three passive scalars at different molecular Prandtl numbers is performed. Computed statistics including the turbulent Prandtl numbers are compared with existing experimental data. The computed fields are also examined to investigate the spatial structure of the scalar fields. The scalar fields are highly correlated with the streamwise velocity; the correlation coefficient between the temperature and the streamwise velocity is as high as 0.95 in the wall region. The joint probability distributions between the temperature and velocity fluctuations are also examined; they suggest that it might be possible to model the scalar fluxes in the wall region in a manner similar to the Reynolds stresses.

Kim, John↗

Interfacial and Transport Property Measurements of Liquid Sodium within a Porous Structure

The sodium thermal electrochemical converter (Na-TEC) is a heat engine that generates electricity via the isothermal expansion of sodium ions through a beta”-alumina solid-electrolyte. The transfer of liquid sodium from a low-pressure condenser to a high-pressure evaporator is necessary to complete the thermodynamic cycle for this conversion process. A unique sodium capillary pump for the Na-TEC is proposed, whereby low-pressure sodium vapor is condensed within a non-wetting (i.e. contact angle > π/2) porous sample. The curvature created at this non-wetting condensation interface supplies a directed force upon the liquid sodium that effectively “pushes” it towards a high pressure region, rather than a traditional capillary wick that “pulls” the liquid. Porous structures made of stainless steel 316 are commercially available, and this material is non-wetting to sodium at the relevant temperature range of the Na-TEC condenser (500K – 673K). An experimental set-up is used to measure several interfacial and transport properties of the condensed sodium within this porous structure. These properties include the temperature dependent Laplace pressure generated at the phase change interface and the permeability of the porous structure. For the former, liquid sodium is in contact with the porous structure and is progressively pressurized with argon gas in increments of ~ 5 kPa until the sodium “breaks through” and the liquid begins to flow. This breakthrough pressure is determined when electrical contact is established between the flowing liquid metal and a set of electrodes cemented into the flow path. The velocity of the flowing sodium can also be measured with these electrodes, and the permeability is then determined by fitting velocity data to the classical Darcy law. Statistically significant results of these measurements, taken at multiple sodium temperatures, are presented. Furthermore, contact angle measurements of sessile sodium drops on stainless steel 316 are demonstrated. This contact angle information is used to define an effective pore radius for the porous structure, which is then validated with microscopy. All the physical parameters described herein are used as inputs in a condensation model to predict the mass flowrate of condensing sodium under a given pressure gradient across the porous structure. This condensation model is used to design the capillary pump for the Na-TEC.

14 SOLAR ENERGY↗

Experiments with Turbulent Jets at Mach Number 0.9

A systematic investigation of the structure of turbulent jets before their interaction with shock or expansion waves was undertaken during the last year. In particular compressibility and density effects in circular jets issuing in still air were investigated experimentally. Jets with nitrogen, helium, and krypton gases at 0.3, 0.6, and 0.9 Mach numbers were investigated in detail. Particle Image Velocimetry technique was developed, tested, and used to obtain qualitative information of the two-dimensional velocity field on a plane inside the flow field, which was illuminated by a laser sheet. The motion of particles was recorded by a CCD camera, which was appropriately triggered to capture two images within a fraction of a microsecond. Statistical averaging of the data at each location reduced the large amount of acquired data. It was found that the spreading rate of the jets was reduced with increased Mach numbers or increased density ratio. It was also found that decay rates of centerline Mach numbers are higher in gases with reduced density ratio. Mach number fluctuations appear to decrease with increasing Mach number of the flow. It has been proposed that the reason for this behavior is the reduction of vortex stretching activities with increased Mach number.

Agui, Juan↗

Critical points of the cosmic velocity field and the uncertainties in the value of the Hubble constant

The existence of critical points for the peculiar velocity field is a natural feature of the correlated vector field. These points appear at the junctions of velocity domains with different orientations of their averaged velocity vectors. Since peculiar velocities are the important cause of the scatter in the Hubble expansion rate, we propose that a more precise determination of the Hubble constant can be made by restricting analysis to a subsample of observational data containing only the zones around the critical points of the peculiar velocity field, associated with voids and saddle points. On large-scales the critical points, where the first derivative of the gravitational potential vanishes, can easily be identified using the density field and classified by the behavior of the Hessian of the gravitational potential. We use high-resolution N-body simulations to show that these regions are stable in time and hence are excellent tracers of the initial conditions. Furthermore, we show that the variance of the Hubble flow can be substantially minimized by restricting observations to the subsample of such regions of vanishing velocity instead of aiming at increasing the statistics by averaging indiscriminately using the full data sets, as is the common approach.

79 ASTRONOMY AND ASTROPHYSICS↗

Biasing errors and corrections

The dependence of laser velocimeter measurement rate on flow velocity is discussed. Investigations outlining that any dependence is purely statistical, and is nonstationary both spatially and temporally, are described. Main conclusions drawn are that the times between successive particle arrivals should be routinely measured and the calculation of the velocity data rate correlation coefficient should be performed to determine if a dependency exists. If none is found, accept the data ensemble as an independent sample of the flow. If a dependency is found, the data should be modified to obtain an independent sample. Universal correcting procedures should never be applied because their underlying assumptions are not valid.

Meyers, James F.↗

Characteristics of inhomogeneous jets in confined swirling air flows

An experimental program to study the characteristics of inhomogeneous jets in confined swirling flows to obtain detailed and accurate data for the evaluation and improvement of turbulent transport modeling for combustor flows is discussed. The work was also motivated by the need to investigate and quantify the influence of confinement and swirl on the characteristics of inhomogeneous jets. The flow facility was constructed in a simple way which allows easy interchange of different swirlers and the freedom to vary the jet Reynolds number. The velocity measurements were taken with a one color, one component DISA Model 55L laser-Doppler anemometer employing the forward scatter mode. Standard statistical methods are used to evaluate the various moments of the signals to give the flow characteristics. The present work was directed at the understanding of the velocity field. Therefore, only velocity and turbulence data of the axial and circumferential components are reported for inhomogeneous jets in confined swirling air flows.

So, R. M. C.↗

Supersonic Traveling Crossflow Wave Characteristics in Ground and Flight Tests

This paper continues analysis of data acquired in a series of supersonic crossflow transition experiments on a 65° swept- wing model tested in both wind tunnel and flight. Flow visualization data are summarized to provide context for the swept-wing instability mechanisms of interest. Unsteady pressure measurements obtained near the wing model leading edge (LE) are studied using 2nd- and 3rd-order statistical methods to extract traveling crossflow wave characteristics from these data. Comparisons between tunnel and flight data with boundary-layer (BL) Linear Stability Theory (LST) predictions provide insights into BL transition phenomena similarities and differences observed in tests using the same swept-wing model configurations. Comparisons of unsteady pressure data results for wing LE configurations with and without distributed-roughness-element (DRE) patterns installed, suggest traveling crossflow disturbance growth is reduced when this type of flow control is applied. A brief stud y of different background surface roughness levels on the wing LE suggests the potential for another flow control approach to promote laminar flow on swept wings using streamwise-biased surface finishes.

Lewis R Owens↗

Search for the chiral magnetic effect in Au+Au collisions at $\sqrt{s_{NN}}$ = 27 GeV with the STAR forward event plane detectors

A decisive experimental test of the Chiral Magnetic Effect (CME) is considered one of the major scientific goals at the Relativistic Heavy-Ion Collider (RHIC) towards understanding the nontrivial topological fluctuations of the Quantum Chromodynamics vacuum. In heavy-ion collisions, the CME is expected to result in a charge separation phenomenon across the reaction plane, whose strength could be strongly energy dependent. The previous CME searches have been focused on top RHIC energy collisions. In this Letter, we present a low energy search for the CME in Au+Au collisions at $\sqrt{s_{NN}}$ = 27 GeV. We measure elliptic flow scaled charge-dependent correlators relative to the event planes that are defined at both mid-rapidity |η| < 1.0 and at forward rapidity 2.1 < |η| < 5.1. We compare the results based on the directed flow plane ( $Ψ$ 1 ) at forward rapidity and the elliptic flow plane ( $Ψ$ 2 ) at both central and forward rapidity. The CME scenario is expected to result in a larger correlation relative to $Ψ$ 1 than to $Ψ$ 2 , while a flow driven background scenario would lead to a consistent result for both event planes. In 10-50% centrality, results using three different event planes are found to be consistent within experimental uncertainties, suggesting a flow driven background scenario dominating the measurement. We obtain an upper limit on the deviation from a flow driven background scenario at the 95% confidence level. This work opens up a possible road map towards future CME search with the high statistics data from the RHIC Beam Energy Scan Phase-II.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗