Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time dependent fluctuations”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Experimental Investigation of Separated and Transitional Boundary Layers Under Low-Pressure Turbine Airfoil Conditions

Modern low-pressure turbine airfoils are subject to increasingly stronger pressure gradients as designers impose higher loading in an effort to improve efficiency and to reduce part count. The adverse pressure gradients on the suction side of these airfoils can lead to boundary-layer separation, particularly under cruise conditions. Separation bubbles, notably those which fail to reattach, can result in a significant degradation of engine efficiency. Accurate prediction of separation and reattachment is hence crucial to improved turbine design. This requires an improved understanding of the transition flow physics. Transition may begin before or after separation, depending on the Reynolds number and other flow conditions, has a strong influence on subsequent reattachment, and may even eliminate separation. Further complicating the problem are the high free-stream turbulence levels in a real engine environment, the strong pressure gradients along the airfoils, the curvature of the airfoils, and the unsteadiness associated with wake passing from upstream stages. Because of the complicated flow situation, transition in these devices can take many paths that can coexist, vary in importance, and possibly also interact, at different locations and instances in time. The present work was carried out in an attempt to systematically sort out some of these issues. Detailed velocity measurements were made along a flat plate subject to the same nominal dimensionless pressure gradient as the suction side of a modern low-pressure turbine airfoil ('Pak-B'). The Reynolds number based on wetted plate length and nominal exit velocity, Re, was varied from 50;000 to 300; 000, covering cruise to takeoff conditions. Low, 0.2%, and high, 7%, inlet free-stream turbulence intensities were set using passive grids. These turbulence levels correspond to about 0.2% and 2.5% turbulence intensity in the test section when normalized with the exit velocity. The Reynolds number and free-stream turbulence level do not have a significant effect on the location of boundary-layer separation unless they are high enough to induce transition upstream of separation. The location and extent of the transition zone, in contrast, depend strongly on Re and TI. The beginning of reattachment closely follows the onset of transition. Under low free-stream turbulence conditions the boundary layer is laminar at separation and then begins to exhibit fluctuations in a finite frequency band in the shear layer over the separation bubble. These fluctuations are due to instability waves. The fluctuations grow in magnitude, higher harmonics are generated, and finally lead to a breakdown to turbulence. Transition begins in the shear layer, but quickly spreads to the near wall region and causes the boundary layer to reattach. The transition is rapid and the resulting turbulence contains a full range of high and low frequencies. Under high free-stream turbulence conditions, slowly growing low-frequency fluctuations are induced in the pretransitional boundary layer by the free-stream. The separation bubbles are considerably thinner than in the low TI cases, resulting in thinner boundary layers at the end of the test wall. At Re=50,000 and 100,000, the pre-transitional boundary layer separates at about the same location as in the low TI cases. Transition occurs through a bypass mode, begins upstream of the corresponding low-TI location, and proceeds in a manner similar to that of an attached boundary layer. Under high TI at Re=200,000 and 300,000, transition begins before separation. The boundary layer may separate, but if it does the separation bubble is very short and does not significantly affect the downstream development of the boundary layer. A comparison is made to previous work in a simulated cascade.

Hultgren, Lennart S.↗
NASA NTRS1991

Incoherent vs. coherent behavior in the normal state of copper oxide superconductors

The self-consistent quantum fluctuations around the mean-field Hartree-Fock state of the Hubbard model provide a very good description of the ground state and low temperature properties of a 2-D itinerant antiferromagnet. Very good agreement with numerical calculations and experimental data is obtained by including the one- and two-loop spin wave corrections to various physical quantities. In particular, the destruction of the long-range order above the Neel temperature can be understood as a spontaneous generation of a length-scale epsilon(T), which should be identified as the spin correlation length. For finite doping, the question of the Hartree-Fock starting point becomes a more complex one since an extra hole tends to self-trap in antiferromagnetic background. Such quantum defects in an underlying antiferromagnetic state can be spin-bags or vortex-like structures and tend to suppress the long-range order. If motion of the holes occurs on a time-scale shorter than the one associated with the motion of these quantum defects of a spin background, one obtains several important empirical features of the normal state of CuO superconductors like linear T-dependence of resistivity, the cusp in the tunneling density of states, etc. As opposed to a familiar Fermi-liquid behavior, the phenomenology of the above system is dominated by a large incoherent piece of a single hole propagator, resulting in many unusual normal state properties.

Tesanovic, Zlatko↗
NASA NTRS2017

Improvement of Automated POST Case Success Rate Using Support Vector Machines

During early conceptual design of complex systems, concept down selection can have a large impact upon program life-cycle cost. Therefore, any concepts selected during early design will inherently commit program costs and affect the overall probability of program success. For this reason it is important to consider as large a design space as possible in order to better inform the down selection process. For conceptual design of launch vehicles, trajectory analysis and optimization often presents the largest obstacle to evaluating large trade spaces. This is due to the sensitivity of the trajectory discipline to changes in all other aspects of the vehicle design. Small deltas in the performance of other subsystems can result in relatively large fluctuations in the ascent trajectory because the solution space is non-linear and multi-modal [1]. In order to help capture large design spaces for new launch vehicles, the authors have performed previous work seeking to automate the execution of the industry standard tool, Program to Optimize Simulated Trajectories (POST). This work initially focused on implementation of analyst heuristics to enable closure of cases in an automated fashion, with the goal of applying the concepts of design of experiments (DOE) and surrogate modeling to enable near instantaneous throughput of vehicle cases [2]. Additional work was then completed to improve the DOE process by utilizing a graph theory based approach to connect similar design points [3]. The conclusion of the previous work illustrated the utility of the graph theory approach for completing a DOE through POST. However, this approach was still dependent upon the use of random repetitions to generate seed points for the graph. As noted in [3], only 8% of these random repetitions resulted in converged trajectories. This ultimately affects the ability of the random reps method to confidently approach the global optima for a given vehicle case in a reasonable amount of time. With only an 8% pass rate, tens or hundreds of thousands of reps may be needed to be confident that the best repetition is at least close to the global optima. However, typical design study time constraints require that fewer repetitions be attempted, sometimes resulting in seed points that have only a handful of successful completions. If a small number of successful repetitions are used to generate a seed point, the graph method may inherit some inaccuracies as it chains DOE cases from the non-global-optimal seed points. This creates inherent noise in the graph data, which can limit the accuracy of the resulting surrogate models. For this reason, the goal of this work is to improve the seed point generation method and ultimately the accuracy of the resulting POST surrogate model. The work focuses on increasing the case pass rate for seed point generation.

Zwack, Matthew R.↗
NASA NTRS2009

Exploring Connections between Global Climate Indices and African Vegetation Phenology

Variations in agricultural production due to rainfall and temperature fluctuations are a primary cause of food insecurity on the continent in Africa. Agriculturally destructive droughts and floods are monitored from space using satellite remote sensing by organizations seeking to provide quantitative and predictive information about food security crises. Better knowledge on the relation between climate indices and food production may increase the use of these indices in famine early warning systems and climate outlook forums on the continent. Here we explore the relationship between phenology metrics derived from the 26 year AVHRR NDVI record and the North Atlantic Oscillation index (NAO), the Indian Ocean Dipole (IOD), the Pacific Decadal Oscillation (PDO), the Multivariate ENSO Index (MEI) and the Southern Oscillation Index (SOI). We explore spatial relationships between growing conditions as measured by the NDVI and the five climate indices in Eastern, Western and Southern Africa to determine the regions and periods when they have a significant impact. The focus is to provide a clear indication as to which climate index has the most impact on the three regions during the past quarter century. We found that the start of season and cumulative NDVI were significantly affected by variations in the climate indices. The particular climate index and the timing showing highest correlation depended heavily on the region examined. The research shows that climate indices can contribute to understanding growing season variability in Eastern, Western and Southern Africa.

Brown, Molly E.↗
NASA NTRS1994

Ion transport and loss in the Earth's quiet ring current. 2: Diffusion and magnetosphere-ionosphere coupling

We have studied the transport and loss of H(+), He(+), and He(++) ions in the Earth's quiet time ring current (1 to 300 keV/e, 3 to 7 R(sub E), Kp less than 2+, absolute value of Dst less than 11, 70 to 110 degs pitchangles, all LT) comparing the standard radial diffusion model developed for the higher-energy radiation belt particles with measurements of the lower energy ring current ions in a previous paper. Large deviations of that model, which fit only 50% of the data to within a factor of 10, suggested that another transport mechanism is operating in the ring current. Here we derive a modified diffusion coefficient corrected for electric field effects on ring current energy ions that fit nearly 80% of the data to within a factor of 2. Thus we infer that electric field fluctuations from the low-latitude to midlatitude ionosphere (ionospheric dynamo) dominated the ring current transport, rather than high-latitude or solar wind fluctuations. Much of the remaining deviation may arise from convective electric field transport of the E less than 30 keV particles. Since convection effects cannot be correctly treated with this azimuthally symmetric model, we defer treatment of the lowest-energy ions to a another paper. We give chi(exp 2) contours for the best fit, showing the dependence of the fit upon the internal/external spectral power of the predicted electric and magnetic field fluctuations.

Sheldon, R. B.↗
NASA NTRS1994

Hydro-isostatic deflection and tectonic tilting in the central Andes: Initial results of a GPS survey of Lake Minchin shorelines

Sufficiently large lake loads provide a means of probing rheological stratification of the crust and upper mantle. Lake Minchin was the largest of the late Pleistocene pluvial lakes in the central Andes. Prominent shorelines, which formed during temporary still-stands in the climatically driven lake level history, preserve records of lateral variations in subsequent net vertical motions. At its maximum extent the lake was 140 m deep and spanned 400 km N-S and 200 km E-R. The load of surficial water contained in Lake Minchin was sufficient to depress the crust and underlying mantle by 20-40 m, depending on the subjacent rheology. Any other differential vertical motions will also be recorded as departures from horizontality of the shorelines. We recently conducted a survey of shoreline elevations of Lake Minchin with the express intent of monitoring the hydro-isostatic deflection and tectonic tilting. Using real-time differential Global Positioning System (GPS), we measured topographic profiles across suites of shorelines at 15 widely separated locations throughout the basin. Horizontal and vertical accuracies attained are roughly 30 and 70 cm, respectively. Geomorphic evidence suggests that the highest shoreline was occupied only briefly (probably less than 200 years) and radiocarbon dates on gastropod shells found in association with the shore deposits constrain the age to roughly 17 kyr. The basin-side pattern of elevations of the highest shoreline is composed of two distinct signals: (27 +/- 1) m of hydro-isostatic deflection due to the lake load, and a planar tilt with east and north components of (6.8 +/- 0.4) 10(exp -5) and 9-5.3 +/- 0.3) 10(exp -5). This rate of tilting is too high to be plausibly attributed to steady tectonism, and presumably reflects some unresolved combination of tectonism plus the effects of oceanic and lacustrine loads on a laterally heterogeneous substrate. The history of lake level fluctuations is still inadequately known to allow detailed inferences of crust and mantle rheology. However, it is already clear that the effective elastic plate thickness is closer to 40 km than the 60-70 km crustal thickness in the central Andes and the effective viscosity is less than 5 10(exp 20) Pa s.

Bills, Bruce G.↗
NASA NTRS1995

The effects of redshifts and focusing on the spectrum of an accretion disk in the galactic center black hole candidate Sagittarius A(sup *)

There are firm indications that Sgr A(sup *), a compact, nonthermal radio source at the Galactic center, may be powered by the dissipation of gravitational energy as gas trapped from an ambient wind descends down the potential well, first through a quasi-spherical inflow (extending out to approximately 3 x 10(exp 16) cm) and then through a small accretion disk at less than or approximately = 5-10 Schwarzschild radii. Earlier three-dimensional Bondi-Hoyle numerical simulations have indicated that fluctuations in the accreted specific angular momentum can lead to a variability in the disk flux on a timescale of years. With greatly improved flux measurements at K and H, and the hint of a approximately 10 minute modulation in the IR luminosity, it is crucial to model the disk emission much more precisely than has been attempted thus far. In this Letter we take into account the effects of Doppler and gravitational redshifts, the light-travel time factor, and the light bending near the black hole to determine the measurable spectrum of Sgr A(sup *) in the increasingly important 10(exp 13) Hz less than or approximately = v less than or approximately = 10(exp 16) Hz frequency range. We find that the relativistic disk spectrum is much softer than its Newtonian counterpart, with a predicted UV flux roughly an order of magnitude smaller than had previously been anticipated. In addition, we find that when the physical conditions in the disk are taken to be consistent with the properties of the quasi-spherical infall (specifically, in terms of the accretion rate and disk size), only a slowly spinning or Schwarzschild black hole appears to fit the observations. Our calculations also reveal that the disk flux is much more weakly dependent on the observer's inclination angle than had been suspected on the basis of earlier Newtonian estimates.

Hollywood, J. M.↗
NASA NTRS1998

Experimental Techniques Verified for Determining Yield and Flow Surfaces

Structural components in aircraft engines are subjected to multiaxial loads when in service. For such components, life prediction methodologies are dependent on the accuracy of the constitutive models that determine the elastic and inelastic portions of a loading cycle. A threshold surface (such as a yield surface) is customarily used to differentiate between reversible and irreversible flow. For elastoplastic materials, a yield surface can be used to delimit the elastic region in a given stress space. The concept of a yield surface is central to the mathematical formulation of a classical plasticity theory, but at elevated temperatures, material response can be highly time dependent. Thus, viscoplastic theories have been developed to account for this time dependency. Since the key to many of these theories is experimental validation, the objective of this work (refs. 1 and 2) at the NASA Lewis Research Center was to verify that current laboratory techniques and equipment are sufficient to determine flow surfaces at elevated temperatures. By probing many times in the axial-torsional stress space, we could define the yield and flow surfaces. A small offset definition of yield (10 me) was used to delineate the boundary between reversible and irreversible behavior so that the material state remained essentially unchanged and multiple probes could be done on the same specimen. The strain was measured with an off-the-shelf multiaxial extensometer that could measure the axial and torsional strains over a wide range of temperatures. The accuracy and resolution of this extensometer was verified by comparing its data with strain gauge data at room temperature. The extensometer was found to have sufficient resolution for these experiments. In addition, the amount of crosstalk (i.e., the accumulation of apparent strain in one direction when strain in the other direction is applied) was found to be negligible. Tubular specimens were induction heated to determine the flow surfaces at elevated temperatures. The heating system induced a large amount of noise in the data. By reducing thermal fluctuations and using appropriate data averaging schemes, we could render the noise inconsequential. Thus, accurate and reproducible flow surfaces (see the figure) could be obtained.

Lerch, Brad A.↗
NASA NTRS1994

Surface temperatures in the polar regions from Nimbus 7 temperature humidity infrared radiometer

Monthly surface temperatures in the Arctic and Antarctic regions have been derived from the 11.5 micrometer thermal infrared channel of the Nimbus 7 temperature humidity infrared radiometer (THIR) for a whole year in 1979 and for a winter and a summer month from 1980 through 1985. The data set shows interannual variability and provides spatial details that allow identification of temperature patterns over sea ice and ice sheet surfaces. For example, the coldest spot in the southern hemisphere is observed to be consistently in the Antarctic plateau in the southern hemisphere, while that in the northern hemisphere is usually located in Greenland, or one of three other general areas: Siberia, the central Arctic, or the Canadian Archipelago. Also, in the southern hemisphere, the amplitude of the seasonal fluctuation of ice sheet temperatures is about 3 times that of sea ice, while in the northern hemisphere, the corresponding fluctuations for the two surfaces are about the same. The main sources of error in the retrieval are cloud and other atmospheric effects. These were minimized by first choosing the highest radiance value from the set of measurements during the day taken within a 30 km by 30 km grid of each daily map. Then the difference of daily maps was taken and where the difference is greater than a certain threshold (which in this case is 12 C), the data element is deleted. Overall, the monthly maps derived from the resulting daily maps are spatially and temporally consistent, are coherent with the topograph y of the Antarctic continent and the location of the sea ice edge, and are in qualitative agreement with climatological data. Quantitatively, THIR data are in good agreement with Antarctic ice sheet surface air temperature station data with a correlation coefficient of 0.997 and a standard deviation of 2.0 C. The absolute values are not as good over the sea ice edges, but a comparison with Russian 2-m drift station temperatures shows very high correlation (with correlation coefficient at 0.998) and a standard deviation of 1.1 C. Overall, the rms error is estimated to be from 1 to 2 C, depending on the surface, while the average bias when compared with in situ data is less than 2 C.

Comiso, Josefino C.↗
NASA NTRS1999

Step Bunching: Influence of Impurities and Solution Flow

Step bunching results in striations even at relatively early stages of its development and in inclusions of mother liquor at the later stages. Therefore, eliminating step bunching is crucial for high crystal perfection. At least 5 major effects causing and influencing step bunching are known: (1) Basic morphological instability of stepped interfaces. It is caused by concentration gradient in the solution normal to the face and by the redistribution of solute tangentially to the interface which redistribution enhances occasional perturbations in step density due to various types of noise; (2) Aggravation of the above basic instability by solution flowing tangentially to the face in the same directions as the steps or stabilization of equidistant step train if these flows are antiparallel; (3) Enhanced bunching at supersaturation where step velocity v increases with relative supersaturation s much faster than linear. This v(s) dependence is believed to be associated with impurities. The impurities of which adsorption time is comparable with the time needed to deposit one lattice layer may also be responsible for bunching; (4) Very intensive solution flow stabilizes growing interface even at parallel solution and step flows; (5) Macrosteps were observed to nucleate at crystal corners and edges. Numerical simulation, assuming step-step interactions via surface diffusion also show that step bunching may be induced by random step nucleation at the facet edge and by discontinuity in the step density (a ridge) somewhere in the middle of a face. The corresponding bunching patterns produce the ones observed in experiment. The nature of step bunching generated at the corners and edges and by dislocation step sources, as well as the also relative importance and interrelations between mechanisms 1-5 is not clear, both from experimental and theoretical standpoints. Furthermore, several laws controlling the evolution of existing step bunches have been suggested, though unambiguous conclusions are still missing. Addressing these issues is the major goal of the present project. The theory addressing the above problem, experimental methods, several figures which include: (1) the spatial wave numbers at which the system is neutrally stable as a function of growth velocity for linear kinetics and supersaturation for nonlinear kinetics; (2) a schematic of the experiment of lysozyme crystal growing under conditions of natural convection; (3) fluctuations in time, t, of the normal growth rate, R(t), vicinal slope, p(t) and Fourier Spectra of R(t), discussions and conclusions are presented.

Chernov, A. A.↗
NASA NTRS2012

Time-Dependent Searches for Point Sources of Neutrinos with the 40-String and 22-String Configurations of IceCube

This paper presents four searches for flaring sources of neutrinos using the IceCube neutrino telescope. For the first time, a search is performed over the entire parameter space of energy, direction and time with sensitivity to neutrino flares lasting between 20 microseconds and a year duration from astrophysical sources. Searches which integrate over time are less sensitive to flares because they are affected by a larger background of atmospheric neutrinos and muons that can be reduced by the use of additional timing information. Flaring sources considered here, such as active galactic nuclei, soft gamma ray repeaters and gamma-ray bursts, are promising candidate neutrino emitters. Two searches are untriggered in the sense that they look for any possible flare in the entire sky and from a predefined catalog of sources from which photon flares have been recorded. The other two searches are triggered by multi-wavelength information on flares from blazars and from a soft gamma-ray repeater. One triggered search uses lightcurves from Fermi-LAT which provides continuous monitoring. A second triggered search uses information where the flux states have been measured only for short periods of time near the flares. The untriggered searches use data taken by 40 strings of IceCube between Apr 5, 2008 and May 20, 2009. The triggered searches also use data taken by the 22-string configuration of IceCube operating between May 31, 2007 and Apr 5, 2008. The results from all four searches are compatible with a fluctuation of the background.

Stamatikos, M.↗
NASA NTRS2012

Time-Dependent Searches for Point Sources of Neutrinos with the 4O-String and 22-String Configurations of IceCube

This paper presents four searches for flaring sources of neutrinos using the IceCube neutrino telescope. For the first time, a search is performed over the entire parameter space of energy, direction, and time with sensitivity to neutrino flares lasting between 20 micro-s and a year duration from astrophysical sources. Searches that integrate over time are less sensitive to flares because they are affected by a larger background of atmospheric neutrinos and muons that can be reduced by the use of additional timing information. Flaring sources considered here, such as active galactic nuclei, soft gamma-ray repeaters, and gamma-ray bursts, are promising candidate neutrino emitters. Two searches are "untriggered" in the sense that they look for any possible flare in the entire sky and from a predefined catalog of sources from which photon flares have been recorded. The other two searches are triggered by multi-wavelength information on flares from blazars and from a soft gamma-ray repeater. One triggered search uses lightcurves from Fermi-LAT which provides continuous monitoring. A second triggered search uses information where the flux states have been measured only for short periods of time near the flares. The untriggered searches use data taken by 40 strings of IceCube between 2008 April 5 and 2009 May 20. The triggered searches also use data taken by the 22-string configuration of IceCube operating between 2007 May 31 and 2008 April 5. The results from all four searches are compatible with a fluctuation of the background.

Abbasi, R.↗
NASA NTRS1996

An Atmospheric Tape Recorder: The Imprint of Tropical Tropopause Temperatures on Stratospheric Water Vapor

We describe observations of tropical stratospheric water vapor q that show clear evidence of large-scale upward advection of the signal from annual fluctuations in the effective 'entry mixing ratio' q(sub E) of air entering the tropical stratosphere. In other words, air is 'marked,' on emergence above the highest cloud tops, like a signal recorded on an upward moving magnetic tape. We define q(sub E) as the mean water vapor mixing ratio, at the tropical tropopause, of air that will subsequently rise and enter the stratospheric 'overworld' at about 400 K. The observations show a systematic phase lag, increasing with altitude, between the annual cycle in q(sub E) and the annual cycle in q at higher altitudes. The observed phase lag agrees with the phase lag calculated assuming advection by the transformed Eulerian-mean vertical velocity of a q(sub E) crudely estimated from 100-hPa temperatures, which we use as a convenient proxy for tropopause temperatures. The phase agreement confirms the overall robustness of the calculation and strongly supports the tape recorder hypothesis. Establishing a quantitative link between q(sub E) and observed tropopause temperatures, however, proves difficult because the process of marking the tape depends subtly on both small- and large-scale processes. The tape speed, or large-scale upward advection speed, has a substantial annual variation and a smaller variation due to the quasi-biennial oscillation, which delays or accelerates the arrival of the signal by a month or two in the middle stratosphere. As the tape moves upward, the signal is attenuated with an e-folding time of about 7 to 9 months between 100 and 50 hPa and about 15 to 18 months between 50 and 20 hPa, constraining possible orders of magnitude both of vertical diffusion K(sub z) and of rates of mixing in from the extratropics. For instance, if there were no mixing in, then K(sub z) would be in the range 0.03-0.09 m(exp 2)/s; this is an upper bound on K(sub z).

Mote, Philip W.↗
NASA NTRSDate not supplied

Exploring Radiation, Water and Energy Closure in Key Deep Tropical Regions Over the Last 3 Decades

The closure of global energy and water budgets at various time scales depends heavily upon the on the processes of water and energy exchange in the deep tropics. A key process observed is the change in convection regimes from disorganized convection to organized deep convection. During these periods an imbalance in key tropical regions is observed of water budget terms between observations and atmospheric reanalysis systems, implying that processes are not simulated well. For instance, periods of imbalance in the water budget may be at least partly attributed to precipitation estimates due to the transition between cloud and precipitation systems in which ice cloud affects microwave radar estimates (Lietmann-Niimi et al. 2024, submitted). A recent paper by Hsaio et al (2023, pre-published version) found that the transition process is characterized by wind shear and longwave cloud radiative feedback signatures as well. This is just one example how these processes also affect the radiative fluxes and thus the energy budget at the TOA. In this presentation, we assess the net radiative budget and components in the deep tropics by comparing multiple data products (GEWEX SRB Rel4IP, CERES EBAF/SYN1deg Ed4.2, ISSCP FH, etc.) and describing the temporal variability of energy balance terms across the deep tropics for the period from 1988 to near present in six key oceanic regions. These six key tropical oceanic regions are West Indian Ocean, the south Indian Ocean east, tropical western, central and eastern Pacific Ocean, and Tropical Atlantic Ocean. In these regions we assess the sensitivity of the fluxes as a function of fluctuations in cloud types in response to various larger scale atmospheric processes (e.g. El Niño, Indian Ocean Dipole). We assess this variability of the radiative flux anomalies (including the net TOA, surface and the atmospheric divergence fluxes) against water vapor divergence, cloud properties (include ISCCP “Weather States, Tselioudis et al. 2021) and larger scale wind shear to assess the uncertainties of energy closure in these regions. These are contrasted against the modeled radiative fluxes from ERA-5 and MERRA-2. Ocean buoy measurements of radiative fluxes are utilized to help assess data radiative flux uncertainty over the nearly 40 year period. In general, the variability for overlap periods of these various data products agrees well, but there are significant differences in the net fluxes that vary according to the rendering of surface, atmospheric, aerosol and cloud properties. These differences also affect closure terms and inferences on systematic uncertainties within the observational estimates. We conclude with recommendations for continuing work

Radiation budget↗
NASA NTRS2004

Numerical Study of Flow Past a Circular Cylinder Using RANS, Hybrid RANS/LES and PANS Formulations

Two multiscale type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds Averaged Navier-Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two-equation (k(sub epsilon)) model with a RANS/LES transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier-Stokes (PANS) model, where the unresolved kinetic energy parameter (f(sub k)) is allowed to vary as a function of grid spacing and the turbulence length scale. Solutions from these models are compared to RANS results and experimental data for a stationary and rotating cylinder. The parameter f(sub k) varies between zero and one and has the characteristic to be equal to one in the viscous sub layer, and when the RANS turbulent viscosity becomes smaller than the LES viscosity. The formulation, usage methodology, and validation example are presented to demonstrate the enhancement of PAB3D's time-accurate and turbulence modeling capabilities. The models are compared to RANS results and experimental data for turbulent separated flows (TS) and laminar separated flows (LS) around stationary and rotating cylinders. For a stationary cylinder, the TS case is accurately simulated using the general two-equation k(sub epsilon) turbulence model (eddy-viscosity model). PAB3D accurately predicts the drag coefficient (CD), lift coefficient (CL) and the Strouhal number (St). The LS case was a challenge for the RANS computation with an eddy-viscosity turbulence model. The RANS/LES and PANS performed well and showed marked improvements over the RANS solution. The modified PANS model was the most accurate. For the rotating cylinder, the spin ratio varied from zero to one, and the PANS results were in good agreement with published experimental data. RANS/LES and PANS capture both temporal and spatial fluctuations and produce large-scale structures that do not occur in the RANS simulation. The current results show promise for the capability of PANS in simulating unsteady and complex flow phenomena.

Elmiligui, Alaa A.↗
NASA NTRS1994

Reionization in a cold dark matter universe: The feedback of galaxy formation on the intergalactic medium

We study the coupled evolution of the intergalactic medium (IGM) and the emerging structure in the universe in the context of the cold dark matter (CDM) model, with a special focus on the consequences of imposing reionization and the Gunn-Peterson constraint as a boundary condition on the model. We have calculated the time-varying density of the IGM by coupling our detailed, numerical calculations of the thermal and ionization balance and radiative transfer in a uniform, spatially averaged IGM of H and He, including the mean opacity of an evolving distribution of gas clumps which correspond to quasar absorption line clouds, to the linearized equations for the growth of density fluctuations in both the gaseous and dark matter components in a CDM universe. We use the linear growth equations to identify the fraction of the gas which must have collapsed out at each epoch, an approach similar in spirit to the so-called Press-Schechter formalism. We identify the IGM density with the uncollapsed baryon fraction. The collapsed fraction is postulated to be a source of energy injection into the IGM, by radiation or bulk hydrodynamical heating (e.g., via shocks) or both, at a rate which is marginally enough to satisfy the Gunn-Peterson constraint at z less than 5. Our results include the following: (1) We find that the IGM in a CDM model must have contained a substantial fraction of the total baryon density of the universe both during and after its reionization epoch. (2) As a result, our previous conclusion that the observed Quasi-Stellar Objects (QSOs) at high redshift are not sufficient to ionize the IGM enough to satisfy the Gunn-Peterson constraint is confirmed. (3) We predict a detectable He II Gunn-Peterson effect at 304(1 + z) A in the spectra of quasars at a range of redshift z greater than or approx. 3, depending on the nature of the sources of IGM reionization. (4) We find, moreover, that a CDM model with high bias parameter b (i.e., b greater than or approx. 2) cannot account for the baryon content of the universe at z approximately 3 observed in quasar absorption line gas unless Omega (sub B) significantly exceeds the maximum value allowed by big bang nucleocynthesis. (5) For a CDM model with bias parameter within the allowed range of (lower) values, the lower limit to Omega(sub B) imposed by big bang nucleosynthesis (Omega(sub B) h(sup 2) greater than or equal to 0.01) combines with our results to yield the minimum IGM density for the CDM fodel. For CDM with b = 1 (Cosmic Background Explorer (COBE) normalization), we find Omega(sub IGM)(sup min) (z approximately 4) approx. equal 0.02-0.03, and Omega(sub IGM)(sup min)(z approximately 0) approx. equal 0.005-0.03, depending upon the nature of the sources of IGM reionization. (6) In general, we find that self-consistent reionization of the IGM by the collapsed baryon fraction has a strong effect on the rate of collapse. (7) As a further example, we show that the feedback effect on the IGM of energy release by the collapsed baryon fraction may explain the slow evolution of the observed comoving QSO number density between z = 5 and z = 2, followed by the sharp decline after z = 2.

Shapiro, Paul R.↗
NASA NTRS1973

On Lidar Sounding of the Atmosphere to Estimate Static and Dynamic Characteristics of Aerosol Inhomogeneities

1. A possible application of intensity fluctuations of a pulse light signal reflected by atmospheric aerosols is analyzed by the correlation method to evaluate static (medium sizes, shape) and dynamic (speed and direction of movement, lifetime) characteristics of aerosol inhomogeneities. The aerosol inhomogeneities are assumed to be expanded, pressed, disintegrated and originated constantly in accordance with random laws, the set of inhomogeneities as a whole traveling together with air masses and having predominant movement in wind direction. It is shown that the characteristics of aerosol inhomogeneities considered can be expressed by the coefficients of the correlation function expansion of the reflected signal fluctuation intensity in Tailor series. 2/ Correlation systems for evaluating static and dynamic characteristics of driving objects can be divided into two types according to the kind and quantity of used information: the systems with coordinates of the information removal "points" to be fixed in space, and the systems with a parallel simultaneous information removal at discrete moments of time. The systems for determination of wind direction considered in are the examples of the first type system. However, the operating information removal for two points is insufficient to estimate completely static and dynamic characteristics of inhomogenities, their quantity ought to be increased up to three of them for two-dimensional problem and up to four of them for three-dimensional problem as it is usually done in the ionospheric studies. The second type systems are used for the investigation of a medium shape and speed of the clouds according to photographs made from satellites. These systems are also used for solution of navigation problems. The use of optical quantum generators with a scanning beam is seen to increase greatly the working information removal in comparison with the first type systems. Nevertheless, scanning rate is not sufficient sometimes in order· to consider a general picture of aerosol inhomogeneities to be stationary. In this connection the use of the systems of second type treatment becomes a matter of essential difficulty. 3. Aerosol inhomogeneities simulation has been carried out on the basis of the digital computer experiments with the aim of estimating static and dynamic characteristics of inhomogeneities by an optical beam in the atmosphere at different scanning procedures. The dependence of determination accuracy of these characteristics on the type of chosen laws of aerosol particle distribution in the atmosphere, the parameters of inhomogeneities geometry, their speed and the law of scanning have been obtained.

Zuev, V. E.↗
NASA NTRS2002

Charge Generation and Propagation in Igneous Rocks

Various electrical phenomena have been reported prior to or concurrent with earthquakes such as resistivity changes, ground potentials, electromagnetic (EM), and luminous signals. Doubts have been raised as to whether some of these phenomena are real and indeed precursory. One of the reasons for uncertainty is that, despite decades of intense work, there is still no physically coherent model. Using low- to medium-velocity impacts to measure electrical signals with microsecond time resolution, it has now been observed that when dry gabbro and diorite cores are impacted at relatively low velocities, approximately 100 m/s, highly mobile charge carriers are generated in a small volume near the impact point. They spread through the rocks, causing electric potentials exceeding +400 mV, EM, and light emission. As the charge cloud spreads, the rock becomes momentarily conductive. When a dry granite block is impacted at higher velocity, approximately 1.5 km/s, the propagation of the P and S waves is registered through the transient piezoelectric response of quartz. After the sound waves have passed, the surface of the granite block becomes positively charged, suggesting the same charge carriers as observed during the low-velocity impact experiments, expanding from within the bulk. During the next 2-3 ms the surface potential oscillates, indicating pulses of electrons injected from ground and contact electrodes. The observations are consistent with positive holes, e.g., defect electrons in the O(2-) sublattice, traveling via the O 2p-dominated valence band of the silicate minerals. Before activation, the positive holes lay dormant in the form of electrically inactive positive hole pairs (PHP), chemically equivalent to peroxy links, O3X/OO\XO3, with X=Si(4+), Al(3+), etc. PHPs are introduced into the minerals by way of hydroxyl,O3X-OH, which all nominally anhydrous minerals incorporate when crystallizing in H2O-laden environments. The fact that positive holes can be activated by low-energy impacts, and their attendant sound waves, suggests that they can also be activated by microfracturing. Depending on where in the stressed rock volume the charge carriers are activated, they will form rapidly moving or fluctuating charge clouds that may account for earthquake-related electrical signals and EM emission. Wherever such charge clouds intersect the surface, high fields are expected, causing electric discharges and earthquake lights.

Freund, Friedemann↗