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 631 records · Page 35

Mechanical Response of Castlegate Sandstone under Hydrostatic Cyclic Loading

The stress history of rocks in the subsurface affects their mechanical and petrophysical properties. Rocks can often experience repeated cycles of loading and unloading due to fluid pressure fluctuations, which will lead to different mechanical behavior from static conditions. This is of importance for several geophysical and industrial applications, for example, wastewater injection and reservoir storage wells, which generate repeated stress perturbations. Laboratory experiments were conducted with Castlegate sandstone to observe the effects of different cyclic pressure loading conditions on a common reservoir analogue. Each sample was hydrostatically loaded in a triaxial cell to a low effective confining pressure, and either pore pressure or confining pressure was cycled at different rates over the course of a few weeks. Fluid permeability was measured during initial loading and periodically between stress cycles. Samples that undergo cyclic loading experience significantly more inelastic (nonrecoverable) strain compared to samples tested without cyclic hydrostatic loading. Permeability decreases rapidly for all tests during the first few days of testing, but the decrease and variability of permeability after this depend upon the loading conditions of each test. Cycling conditions do affect the mechanical behavior; the elastic moduli decrease with the increasing loading rate and stress cycling. The degree of volumetric strain induced by stress cycles is the major control on permeability change in the sandstones, with less compaction leading to more variation from measurement to measurement. The data indicate that cyclic loading degrades permeability and porosity more than static conditions over a similar period, but the petrophysical properties are dictated more by the hydrostatic loading rate rather than the total length of time stress cycling is imposed.

58 GEOSCIENCES↗

Constraining the destruction rate of K 40 in stellar nucleosynthesis through the study of the Ar 40 ( p , n ) K 40 reaction

Background: We present that K 40 plays a significant role in the radiogenic heating of Earth-like exoplanets, which can affect the development of a habitable environment on their surfaces. The initial amount of K 40 in the interior of these planets depends on the composition of the interstellar clouds from which they formed. Within this context, nuclear reactions that regulate the production of K 40 during stellar evolution can play a critical role. Purpose: In this study, we constrain for the first time the astrophysical reaction rate of K 40 ( n , p ) Ar 40 , which is responsible for the destruction of K 40 during stellar nucleosynthesis. We provide to the nuclear physics community high-resolution data on the cross section and angular distribution of the Ar 40 ( p , n ) K 40 reaction. These are important to various applications involving Ar 40 . The associated reaction rate of the Ar 40 ( p , n ) K 40 process addresses a reaction rate gap in the Joint Institute for Nuclear Astrophysics REACLIB database in the region of intermediate-mass isotopes. Methods: We performed differential cross-section measurements on the Ar 40 ( p , n ) K 40 reaction, for six energies in the center-of-mass system between 3.2 and 4.0 MeV and various angles between 0 ° and 135 ° . The experiment took place at the Edwards Accelerator Laboratory at Ohio University using the beam swinger target location and a standard neutron time-of-flight technique. We extracted total and partial cross sections by integrating the double differential cross sections we measured. Results: The total and partial cross sections varied with energy due to the contribution from isobaric analog states and Ericson type fluctuations. The energy-averaged neutron angular distributions were symmetrical relative to 90 ° . Based on the experimental data, local transmission coefficients were extracted and were used to calculate the astrophysical reaction rates of Ar 40 ( p , n ) K 40 and K 40 ( n , p ) Ar 40 reactions. The new rates were found to vary significantly from the theoretical rates in the REACLIB library. We implemented the new rates in network calculations to study nucleosynthesis via the slow neutron capture process, and we found that the produced abundance of K 40 is reduced by up to 10% compared to calculations with the library rates. At the same time, the above result removes a significant portion of the previous theoretical uncertainty on the K 40 yields from stellar evolution calculations. Conclusions: Our results support a destruction rate of K 40 in massive stars via the K 40 ( n , p ) Ar 40 reaction that is larger compared to previous estimates. The rate of K 40 destruction via the K 40 ( n , p ) Ar 40 reaction now has a dramatically reduced uncertainty based on our measurement. Lastly, this result directly affects the predicted stellar yields of K 40 from nucleosynthesis, which is a critical input parameter for the galactic chemical evolution models that are currently employed for the study of significant properties of exoplanets.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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

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

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

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↗

Differential equations for cosmological correlators

Cosmological fluctuations retain a memory of the physics that generated them in their spatial correlations. The strength of correlations varies smoothly as a function of external kinematics, which is encoded in differential equations satisfied by cosmological correlation functions. In this work, we provide a broader perspective on the origin and structure of these differential equations. As a concrete example, we study conformally coupled scalar fields in a power-law cosmology. The wavefunction coefficients in this model have integral representations, with the integrands being the product of the corresponding flat-space results and “twist factors” that depend on the cosmological evolution. Similar twisted integrals arise for loop amplitudes in dimensional regularization, and their recent study has led to the discovery of rich mathematical structures and powerful new tools for computing multi-loop Feynman integrals in quantum field theory. The integrals of interest in cosmology are also part of a finite-dimensional basis of master integrals, which satisfy a system of first-order differential equations. We develop a formalism to derive these differential equations for arbitrary tree graphs. The results can be represented in graphical form by associating the singularities of the differential equations with a set of graph tubings. Upon differentiation, these tubings grow in a local and predictive fashion. In fact, a few remarkably simple rules allow us to predict — by hand — the equations for all tree graphs. While the rules of this “kinematic flow” are defined purely in terms of data on the boundary of the spacetime, they reflect the physics of bulk time evolution. We also study the analogous structures in tr ϕ 3 theory, and see some glimpses of hidden structure in the sum over planar graphs. This suggests that there is an autonomous combinatorial or geometric construction from which cosmological correlations, and the associated spacetime, emerge.

Cosmological models↗

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

Understanding Life History Diversity of a Wild Steelhead Population and Managing for Resiliency

Abstract Diversity in steelhead Oncorhynchus mykiss (anadromous Rainbow Trout) is often characterized by varying time spent in freshwater and saltwater environments. However, ontogenetic shifts are not often highlighted as an important diversity characteristic in salmonids. Larger water bodies can play a key role in population stability by buffering limited resources in smaller tributary systems. Based on freshwater and saltwater ages of adult steelhead in Fish Creek, Idaho, there were 12 different life history trajectories observed. Juvenile abundance varied from 12,083 (SE = 1,225.7) to 62,463 (SE = 1,753.8), with ages ranging from young of the year (age 0) to 4 years. Adult abundance fluctuated from 17 (SE < 0.1) to 499 (SE = 67.1), and total age of adults varied from 3 to 7 years. Diversity was most evident in juvenile movement and rearing strategies. We observed four types of movement and rearing pathways in which juveniles either directly migrated to the ocean from Fish Creek or resided for up to three winters in main-stem river habitat in the Clearwater River basin prior to ocean migration. Most juveniles emigrated in the fall at age 2, and most ocean-emigrating smolts were age 3. Subsequently, most juveniles resided for at least one additional winter period before ocean migration. Juvenile diversity was also reflected in adult ages; the mean proportion that spent 3 years in freshwater was 0.55 (SE = 0.03), and the mean proportion that spent 2 years in freshwater was 0.41 (SE = 0.03). Adult sex ratio was female biased (mean ± SE = 0.66 ± 0.02), with a higher proportion of females (0.89 ± 0.03) than males (0.65 ± 0.05) residing in the ocean for at least 2 years. Density-dependent mechanisms were evident in juvenile rearing habitat in Fish Creek for cohorts of high female abundance. Population productivity potentials can be met when tributary and main-stem rearing habitat can be utilized, and this will be essential in attaining recovery goals.

Dobos, Marika E.↗

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

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

Improving coastal water level estimation by merging nadir-only satellite altimetry data into a hydrodynamic model

Providing robust real time flood warnings is of paramount importance to coastal communities. Although state-of-the-art hydrodynamic models are capable of robustly predicting Coastal Water Levels (CWL), unresolved drivers affecting level fluctuations are often not represented by the model governing equations. This work evaluates a novel method to improve the performance of the ADvanced CIRCulation (ADCIRC) hydrodynamic model by assimilating observations from four nadir-only satellite altimetry missions against a set of National Oceanic and Atmospheric Administration (NOAA) gauge stations located across the entire U.S. East Coast. Two different types of simulations were performed – Open Loop (OL) and Data Assimilation (DA). Five different simulations were performed where four different satellite altimetry observations were assimilated individually and combined with two different scenarios – with and without considering the data quality flags. Results indicate that, despite their limited spatial coverage, merging nadir-only observations into ADCIRC from the newly launched Surface Water and Ocean Topography (SWOT)’s nadir altimeter can improve the model performance at 76% of the gauge locations, whereas Sentinel-6 improves 73% of the total locations, Jason-3 74%, and SARAL 21%. Furthermore, combining observations from SWOT-nadir, Jason-3, and Sentinel-6 can improve the ADCIRC performance at more than 80% of the gauge locations for 107-day simulation. Nadir-only satellite altimetry observations can be useful for improving the model performance even if flagged as “poor quality” near the coast. When the flagged data are disregarded, SWOT can improve ADCIRC at 78%, Sentinel-6 at 73%, Jason-3 at 53%, and SARAL at 21% of the gauge locations. The ability to improve the model simulations largely depends on the availability of a satellite overpass nearby. Therefore, model performance can be further enhanced if satellite observations are available during a storm surge event, stressing the importance of frequent satellite overpasses.

Aafnan Bhuiyan, Soelem↗

Scalability and Effectiveness of Smart Charge Management

The rise in electric vehicle (EV) adoption presents growing challenges for power grids, particularly from simultaneous residential charging, which can cause voltage fluctuations and increase feeder peak loads. Baltimore Gas and Electric (BGE), with support from the U.S. Department of Energy, initiated a pilot program to evaluate managed residential EV charging through Smart Charge Management (SCM). This study analyzes real-world charging behavior data from the pilot and feeder-level base loads from BGE to simulate residential charging scenarios through 2035 across the Washington, DC–Baltimore region. Grid impacts under unmanaged charging are compared to three SCM strategies: TOU-immediate, TOU-distributed, and Load Balancing. Results show that the magnitude of peak reduction is highly feeder-dependent. Some feeders achieve reductions of more than 40% at high enrollment levels, while others show improvements closer to 10–15%. This heterogeneity reflects differences in baseline feeder load shapes, EV penetration, and plug-in behavior across customers. Results also highlight trade-offs between shifting load away from peak periods and minimizing secondary demand peaks, offering practical insights for future utility program design.

Electric vehicle↗

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↗

Image Correlation Pattern Optimization for Micro-Scale In-Situ Strain Measurements

The accuracy and precision of digital image correlation (DIC) is a function of three primary ingredients: image acquisition, image analysis, and the subject of the image. Development of the first two (i.e. image acquisition techniques and image correlation algorithms) has led to widespread use of DIC; however, fewer developments have been focused on the third ingredient. Typically, subjects of DIC images are mechanical specimens with either a natural surface pattern or a pattern applied to the surface. Research in the area of DIC patterns has primarily been aimed at identifying which surface patterns are best suited for DIC, by comparing patterns to each other. Because the easiest and most widespread methods of applying patterns have a high degree of randomness associated with them (e.g., airbrush, spray paint, particle decoration, etc.), less effort has been spent on exact construction of ideal patterns. With the development of patterning techniques such as microstamping and lithography, patterns can be applied to a specimen pixel by pixel from a patterned image. In these cases, especially because the patterns are reused many times, an optimal pattern is sought such that error introduced into DIC from the pattern is minimized. DIC consists of tracking the motion of an array of nodes from a reference image to a deformed image. Every pixel in the images has an associated intensity (grayscale) value, with discretization depending on the bit depth of the image. Because individual pixel matching by intensity value yields a non-unique scale-dependent problem, subsets around each node are used for identification. A correlation criteria is used to find the best match of a particular subset of a reference image within a deformed image. The reader is referred to references for enumerations of typical correlation criteria. As illustrated by Schreier and Sutton and Lu and Cary systematic errors can be introduced by representing the underlying deformation with under-matched shape functions. An important implication, as discussed by Sutton et al., is that in the presence of highly localized deformations (e.g., crack fronts), error can be reduced by minimizing the subset size. In other words, smaller subsets allow the more accurate resolution of localized deformations. Contrarily, the choice of optimal subset size has been widely studied and a general consensus is that larger subsets with more information content are less prone to random error. Thus, an optimal subset size balances the systematic error from under matched deformations with random error from measurement noise. The alternative approach pursued in the current work is to choose a small subset size and optimize the information content within (i.e., optimizing an applied DIC pattern), rather than finding an optimal subset size. In the literature, many pattern quality metrics have been proposed, e.g., sum of square intensity gradient (SSSIG), mean subset fluctuation, gray level co-occurrence, autocorrelation-based metrics, and speckle-based metrics. The majority of these metrics were developed to quantify the quality of common pseudo-random patterns after they have been applied, and were not created with the intent of pattern generation. As such, it is found that none of the metrics examined in this study are fit to be the objective function of a pattern generation optimization. In some cases, such as with speckle-based metrics, application to pixel by pixel patterns is ill-conditioned and requires somewhat arbitrary extensions. In other cases, such as with the SSSIG, it is shown that trivial solutions exist for the optimum of the metric which are ill-suited for DIC (such as a checkerboard pattern). In the current work, a multi-metric optimization method is proposed whereby quality is viewed as a combination of individual quality metrics. Specifically, SSSIG and two auto-correlation metrics are used which have generally competitive objectives. Thus, each metric could be viewed as a constraint imposed upon the others, thereby precluding the achievement of their trivial solutions. In this way, optimization produces a pattern which balances the benefits of multiple quality metrics. The resulting pattern, along with randomly generated patterns, is subjected to numerical deformations and analyzed with DIC software. The optimal pattern is shown to outperform randomly generated patterns.

Bomarito, G. F.↗

Improving the Optical Quality Factor of the WGM Resonator

Resonators usually are characterized with two partially dependent values: finesse (F) and quality factor (Q). The finesse of an empty Fabry-Perot (FP) resonator is defined solely by the quality of its mirrors and is calculated as F=piR(exp 1/2)/(1-R). The maximum up-to-date value of reflectivity R approximately equal to 1 - 1.6 x 10(exp -6) is achieved with dielectric mirrors. An FP resonator made with the mirrors has finesse F=1.9 x 10(exp 6). Further practical increase of the finesse of FP resonators is problematic because of the absorption and the scattering of light in the mirror material through fundamental limit on the reflection losses given by the internal material losses and by thermodynamic density fluctuations on the order of parts in 109. The quality factor of a resonator depends on both its finesse and its geometrical size. A one-dimensional FP resonator has Q=2 F L/lambda, where L is the distance between the mirrors and lambda is the wavelength. It is easy to see that the quality factor of the resonator is unlimited because L is unlimited. F and Q are equally important. In some cases, finesse is technically more valuable than the quality factor. For instance, buildup of the optical power inside the resonator, as well as the Purcell factor, is proportional to finesse. Sometimes, however, the quality factor is more valuable. For example, inverse threshold power of intracavity hyperparametric oscillation is proportional to Q(exp 2) and efficiency of parametric frequency mixing is proportional to Q(exp 3). Therefore, it is important to know both the maximally achievable finesse and quality factor values of a resonator. Whispering gallery mode (WGM) resonators are capable of achieving larger finesse compared to FP resonators. For instance, fused silica resonators with finesse 2.3 x 10(exp 6) and 2.8 x 10(exp 6) have been demonstrated. Crystalline WGM resonators reveal even larger finesse values, F=6.3 x 10(exp 6), because of low attenuation of light in the transparent optical crystals. The larger values of F and Q result in the enhancement of various nonlinear processes. Low-threshold Raman lasing, optomechanical oscillations, frequency doubling, and hyperparametric oscillations based on these resonators have been recently demonstrated. Theory predicts a possibility of nearly 10(exp 14) room-temperature optical Q-factors of optical crystalline WGM resonators, which correspond to finesse levels higher than 10(exp 9). Experiments have shown numbers a thousand times lower than that. The difference occurs due to media imperfections. To substantially reduce the optical losses caused by the imperfections, a specific, multi-step, asymptotic processing of the resonator is implemented. The technique has been initially developed to reduce microwave absorption in dielectric resonators. One step of the process consists of mechanical polishing performed after high temperature annealing. Several steps repeat one after another to lead to significant reduction in optical attenuation and, as a result, to the increase of Q-factor as well as finesse of the resonator which demonstrates a CaF2 WGM resonator with F greater than 10(exp 7) and Q greater than 10(exp 11).

Savchenkov, Anatoliy↗

Experimental and Theoretical Studies of Pulsating Turbulent Flow

The objective of this investigation was to study the effects of small amplitude sinusoidal pulsations on fully developed turbulent flow in a tube from both experimental and theoretical viewpoints. Theoretical models for the macroscopic behavior of pulsating turbulent tube flow were developed for the two cases of very low and very high pulsation frequencies. The models are based on assumptions of quasi-steady and frozen eddy viscosity flow behavior, respectively. The models successfully predict unsteady velocity profiles, thereby supporting the currently proposed definitions of frequency regimes in pulsating turbulent flow. Experimental measurements were made of the time-dependent pressure drop and velocity profiles over the range of frequency-to-Reynolds number ratios from 0.0095 to 0.24. The two macroscopic models developed in this study predict unsteady velocity profiles which are in moderately good agreement with the experiments in their respective frequency regimes, and a previously developed quasi-steady model is found to predict experimental velocity profiles well in both the quasisteady and the frozen eddy viscosity frequency regimes. The effect of flow pulsations on the dissipation of turbulence energy in the vicinity of the wall was measured in the lower transition frequency regime. The long-time averaged dissipation was observed to be unchanged from the steady flow dissipation, within the accuracy of the experiment. A theoretical model of the periodic viscous sublayer was also developed and applied to pulsating flow in a tube, in order to investigate the effects of flow pulsations on the rate of production of turbulence in the region of the wall. The periodic viscous sublayer model predicts sublayer growth periods in steady flow which agree with the published experimental data. When the model is applied to pulsating flow, the response of the sublayer growth period falls into three frequency regimes, the parameters of which are in approximate agreement with the frequency regimes which are defined on the basis of macroscopic flow behavior. The sublayer renewal cycle exhibits quasi-steady flow behavior when the sublayer growth period is much less than the pulsation period, transition behavior when these two periods are approximately equal, and frozen eddy viscosity behavior when the sublayer period is much longer than the pulsation period. The effect of the sublayer growth and renewal cycle on the level of turbulence was investigated by two methods. The velocity fluctuations seen by a point velocity probe located close to the wall were predicted from the model in one method and the rate of turbulence production was estimated from the frequency of sublayer renewal events in the other.

Kingston, G. C.↗

Electronic 1/f noise as a probe of dimensional effects on spin-glass dynamics

Over the past decade, spin-glass simulations have improved to the point that they now access time- and length-scales comparable to experiments at the mesoscale. A recent series of thin-film field-cooled/zero-field-cooled magnetization (FC/ZFC) experiments demonstrated activated spin dynamics, with a temperature-independent activation energy proportional to the logarithm of the film thickness and with coefficients in remarkable agreement with the simulation. These measurements require the application of small magnetic fields, which has been shown to affect the spin-glass energy landscape. Measurements of the 1/f noise in metallic spin-glasses have been previously shown to be a sensitive probe of the spin dynamics, and the measurements can be made without applying a magnetic field. In this mini-review, we review these techniques and discuss how transport measurements can fit into the current landscape of spin-glass measurements. We compare previous measurements to more recent measurements on similar films, made with ostensibly different cooling protocols, and compare both the previous and recent measurements to the magnetometry. The transport measurements—taken over a wider range of temperature than magnetometry—suggest that the maximum spin-glass energy barrier height is temperature-dependent, not fixed, possibly due to two-dimensional dynamics. We discuss this possibility, along with future measurements, which may be able to resolve this mystery.

1/f noise↗