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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Phase-space methods for neutrino oscillations: Extension to multibeams

The phase-space approach (PSA), which was originally introduced in Lacroix [] to describe neutrino flavor oscillations for interacting neutrinos emitted from stellar objects is extended to describe arbitrary numbers of neutrino beams. The PSA is based on mapping the quantum fluctuations into a statistical treatment by sampling initial conditions followed by independent mean-field evolution. A new method is proposed to perform this sampling that allows treating an arbitrary number of neutrinos in each neutrino beams. We validate the technique successfully and confirm its predictive power on several examples where a reference exact calculation is possible. We show that it can describe many-body effects, such as entanglement and dissipation induced by the interaction between neutrinos. Due to the complexity of the problem, exact solutions can only be calculated for rather limited cases, with a limited number of beams and/or neutrinos in each beam. The PSA approach considerably reduces the numerical cost and provides an efficient technique to accurately simulate arbitrary numbers of beams. Examples of PSA results are given here, including up to 200 beams with time-independent or time-dependent Hamiltonians. We anticipate that this approach will be useful to bridge exact microscopic techniques with more traditional transport theories used in neutrino oscillations. It will also provide important reference calculations for future quantum computer applications where other techniques are not applicable to classical computers. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Dephasing in Fluxonium Qubits from Coherent Quantum Phase Slips

Phase slips occur across all Josephson junctions (JJs) at a rate that increases with the impedance of the junction. In superconducting qubits composed of JJ-array superinductors—such as fluxonium—phase slips in the array can lead to decoherence. In particular, phase-slip processes at the individual array junctions can coherently interfere, each with an Aharonov-Casher phase that depends on the offset charges of the array islands. These coherent quantum phase slips (CQPS) perturbatively modify the qubit frequency, and therefore charge noise on the array islands will lead to dephasing. By varying the impedance of the array junctions, we design a set of fluxonium qubits in which the expected phase-slip rate within the JJ array changes by several orders of magnitude. We characterize the coherence times of these qubits and demonstrate that the scaling of CQPS-induced dephasing rates agrees with our theoretical model. Furthermore, we perform noise spectroscopy of two qubits in regimes dominated by either CQPS or flux noise. We find that the noise power spectrum associated with CQPS dephasing appears to be featureless at low frequencies and not 1 / f . Numerical simulations indicate that this behavior is consistent with charge noise generated by charge-parity fluctuations within the array. Our findings broadly inform JJ-array-design trade-offs, relevant for the numerous superconducting-qubit designs employing JJ-array superinductors. Published by the American Physical Society 2024

Physics↗
NASA NTRS1997

Improved Estimates of Clear Sky Longwave Flux and Application to the Tropical Greenhouse Effect

The first objective of this investigation is to eliminate the clear-sky offset introduced by the scene-identification procedures developed for the Earth Radiation Budget Experiment (ERBE). Estimates of this systematic bias range from 10 to as high as 30 W/sq m. The initial version of the ScaRaB data is being processed with the original ERBE algorithm. Since the ERBE procedure for scene identification is based upon zonal flux averages, clear scenes with longwave emission well below the zonal mean value are mistakenly classified as cloudy. The erroneous classification is more frequent in regions with deep convection and enhanced mid- and upper-tropospheric humidity. We will develop scene identification parameters with zonal and/or time dependence to reduce or eliminate the bias in the clear- sky data. The modified scene identification procedure could be used for the ScaRaB-specific version of the Earth-radiation products. The second objective is to investigate changes in the clear-sky Outgoing Longwave Radiation (OLR) associated with decadal variations in the tropical and subtropical climate. There is considerable evidence for a shift in the climate state starting in approximately 1977. The shift is accompanied by higher SSTs in the equatorial Pacific, increased tropical convection, and higher values of atmospheric humidity. Other evidence indicates that the humidity in the tropical troposphere has been steadily increasing over the last 30 years. It is not known whether the atmospheric greenhouse effect has increased during this period in response to these changes in SST and precipitable water. We will investigate the decadal-scale fluctuations in the greenhouse effect using Nimbus-7, ERBE, and ScaRaB measurements spaning 1979 to the present. The data from the different satellites will be intercalibrated by comparison with model calculations based upon ship radiosonde observations. The fluxes calculated from the radiation model will also be used for validation of the ScaRaB fluxes.

Collins, W. D.↗
NASA NTRS1978

Large eddy simulation of incompressible turbulent channel flow

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

Moin, P.↗
NASA NTRS2006

Solar Forced Dansgaard/Oeschger Events?

Climate records for the last ice age (which ended 11,500 years ago) show enormous climate fluctuations in the North Atlantic region - the so-called Dansgaard/Oeschger events. During these events air temperatures in Greenland changed on the order of 10 degrees Celsius within a few decades. These changes were attributed to shifts in ocean circulation which influences the warm water supply from lower latitudes to the North Atlantic region. Interestingly, the rapid warmings tend to recur approximately every 1500 years or multiples thereof. This has led researchers to speculate about an external cause for these changes with the variable Sun being one possible candidate. Support for this hypothesis came from climate reconstructions, which suggested that the Sun influenced the climate in the North Atlantic region on these time scales during the last approximately 12,000 years of relatively stable Holocene climate. However, Be-10 measurements in ice cores do not indicate that the Sun caused or triggered the Dansgaard/Oeschger events. Depending on the solar magnetic shielding more or less Be-10 is produced in the Earth's atmosphere. Therefore, 10Be can be used as a proxy for solar activity changes. Since Be-10 can be measured in ice cores, it is possible to compare the variable solar forcing directly with the climate record from the same ice core. This removes any uncertainties in the relative dating, and the solar-climate link can be reliably studied. Notwithstanding that some Dansgaard/Oeschger warmings could be related to increased solar activity, there is no indication that this is the case for all of the Dansgaard/Oeschger events. Therefore, during the last ice age the Be-10 and ice core climate data do not indicate a persistent solar influence on North Atlantic climate.

Muscheler, R.↗
NASA NTRS2008

An Analytical Model for the Prediction of a Micro-Dosimeter Response Function

A rapid analytical procedure for the prediction of a micro-dosimeter response function in low Earth orbit (LEO), correlated with the Space Transportation System (STS, shuttle) Tissue Equivalent Proportional Counter (TEPC) measurements is presented. The analytical model takes into consideration the energy loss straggling and chord length distribution of the detector, and is capable of predicting energy deposition fluctuations in a cylindrical micro-volume of arbitrary aspect ratio (height/diameter) by incoming ions through both direct and indirect (ray) events. At any designated (ray traced) target point within the vehicle, the model accepts the differential flux spectrum of Galactic Cosmic Rays (GCR) and/or trapped protons at LEO as input. On a desktop PC, the response function of TEPC for each ion in the GCR/trapped field is computed at the average rate of 30 seconds/ion. The ionizing radiation environment at LEO is represented by O'Neill fs GCR model (2004), covering charged particles in the 1 less than or equal to Z less than or equal to 28. O'Neill's free space GCR model is coupled with the Langley Research Center (LaRC) angular dependent geomagnetic cutoff model to compute the transmission coefficient in LEO. The trapped proton environment is represented by a LaRC developed time dependent procedure which couples the AP8MIN/AP8MAX, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment is represented by the extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. The charged particle transport calculations correlated with STS 51 and 114 flights are accomplished by using the most recent version (2005) of the LaRC deterministic High charge (Z) and Energy TRaNsport (HZETRN) code. We present the correlations between the TEPC model predictions (response function) and TEPC measured differential/integral spectra in the lineal energy (y) domain for both GCR and trapped protons, with the conclusion that the model correctly accounts for the increase in flux at low y values where energetic ions are the primary contributor. We further discuss that, even with the incorporation of angular dependency in the cutoffs, comparison of the GCR differential/integral flux between STS 51 and 114 TEPC measured data and current calculations indicates that there still exists an underestimation by the simulations at low to mid range y values. This underestimation is partly related the exclusion of the secondary pion particle production from the current version of HZETRN.

Badavi, Francis F.↗
DOE OSTI2022

The economic value of photovoltaic performance loss mitigation in electricity spot markets

Photovoltaic (PV) performance is affected by reversible and irreversible losses. These can typically be mitigated through responsive and proactive operations and maintenance (O&M) activities. However, to generate profit, the cost of O&M must be lower than the value of the recovered electricity. This value depends both on the amount of recovered energy and on the electricity prices, which can vary significantly over time in spot markets. The present work investigates the impact of the electricity price variability on the PV profitability and on the related O&M activities in Italy, Portugal, and Spain. Here, it is found that the PV revenues varied by 1.6 × to 1.8 × within the investigated countries in the last 5 years. Moreover, forecasts predict higher average prices in the current decade compared to the previous one. These will increase the future PV revenues by up to 60% by 2030 compared to their 2015–2020 mean values. These higher revenues will make more funds available for better maintenance and for higher quality components, potentially leading to even higher energy yield and profits. Linearly growing or constant price assumptions cannot fully reproduce these expected price trends. Furthermore, significant price fluctuations can lead to unexpected scenarios and alter the predictions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRS1996

Nonlinear Response of Layer Growth Dynamics in the Mixed Kinetics-Bulk-Transport Regime

In situ high-resolution interferometry on horizontal facets of the protein lysozyme reveal that the local growth rate R, vicinal slope p, and tangential (step) velocity v fluctuate by up to 80% of their average values. The time scale of these fluctuations, which occur under steady bulk transport conditions through the formation and decay of step bunches (macrosteps), is of the order of 10 min. The fluctuation amplitude of R increases with growth rate (supersaturation) and crystal size, while the amplitude of the v and p fluctuations changes relatively little. Based on a stability analysis for equidistant step trains in the mixed transport-interface-kinetics regime, we argue that the fluctuations originate from the coupling of bulk transport with nonlinear interface kinetics. Furthermore, step bunches moving across the interface in the direction of or opposite to the buoyancy-driven convective flow increase or decrease in height, respectively. This is in agreement with analytical treatments of the interaction of moving steps with solution flow. Major excursions in growth rate are associated with the formation of lattice defects (striations). We show that, in general, the system-dependent kinetic Peclet number, Pe(sub k) , i.e., the relative weight of bulk transport and interface kinetics in the control of the growth process, governs the step bunching dynamics. Since Pe(sub k) can be modified by either forced solution flow or suppression of buoyancy-driven convection under reduced gravity, this model provides a rationale for the choice of specific transport conditions to minimize the formation of compositional inhomogeneities under steady bulk nutrient crystallization conditions.

Vekilov, Peter G.↗
NASA NTRS2022

Periodic Solar Wind Structures Observed in Measurements of Elemental and Ionic Composition in situ at L1

Mesoscale periodic structures observed in solar wind plasma serve as an important diagnostic tool for constraining the processes that govern the formation of the solar wind. These structures have been observed in situ and in remote data as fluctuations in proton and electron density. However, only two events of this type have been reported regarding the elemental and ionic composition. Composition measurements are especially important in gaining an understanding of the origin of the solar wind as the composition is frozen into the plasma at the Sun and does not evolve as it advects through the heliosphere. Here, we present the analysis of four events containing mesoscale periodic solar wind structure during which the Iron and Magnesium number density data, measured by the Solar Wind Ion Composition Spectrometer (SWICS)on board the Advanced Composition Explorer spacecraft, are validated at statistically significant count levels. We use a spectral analysis method specifically designed to extract periodic signals from astrophysical time series and apply it to the SWICS 12-minute native resolution dataset. We find variations in the relative abundance of elements with low first ionization potential, mass dependencies, and charge state during time intervals in which mesoscale periodic structures are observed. These variations are linked to temporal or spatial variations in solar source regions and put constraints on the solar wind formation mechanisms that produce them. Techniques presented here are relevant for future, higher-resolution studies of data from new instruments such as Solar Orbiter’s Heavy Ion Sensor.

Irena Gershkovich↗
DOE OSTI2021

Singularities in nearly uniform one-dimensional condensates due to quantum diffusion

Dissipative systems often exhibit wavelength-dependent loss rates. One prominent example is Rydberg polaritons formed by electromagnetically induced transparency, which have long been a leading candidate for studying the physics of interacting photons and also hold promise as a platform for quantum information. In this system, dissipation is in the form of quantum diffusion, i.e., proportional to k 2 (k being the wavevector) and vanishing at long wavelengths as k→0. Here, we show that one-dimensional condensates subject to this type of loss are unstable to long-wavelength density fluctuations in an unusual manner: after a prolonged period in which the condensate appears to relax to a uniform state, local depleted regions quickly form and spread ballistically throughout the system. We connect this behavior to the leading-order equation for the nearly uniform condensate—a dispersive analog to the Kardar-Parisi-Zhang equation—which develops singularities in finite time. Furthermore, we show that the wavefronts of the depleted regions are described by purely dissipative solitons within a pair of hydrodynamic equations, with no counterpart in lossless condensates. We close by discussing conditions under which such singularities and the resulting solitons can be physically realized.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1998

TOPEX/Poseidon orbit maintenance for the first five years

The TOPEX/Poseidon orbit maintenance strategy was changed following launch to include the effects of observed unmodeled, and hence anomalous, along-track accelerations. The anomalous force causes the semi-major axis, a, to either increase (called "boost") or decrease ("deboost" or "decay") depending on the satellite attitude and solar array pitch angle offset. Although this force is the most uncertain parameter in ground track prediction, it has been used as a passive technique for orbit maintenance, thereby reducing the number of propulsive maneuvers, enhancing maneuver spacing, and to place maneuvers at convenient times. This passive technique was first demonstrated in May 1993. The TOPEX/Poseidon orbit has been uniquely maintained using both passive (non-propulsive) and active (propulsive) maneuvers. Furthermore, the orbit has been maintained using only the passive technique since the ninth orbit maintenance maneuver on January 15, 1996. Only nine orbit maintenance maneuvers have been required to maintain the ground track, including verification site over flights, since achieving the operational orbit on September 21, 1992 (mission requirement: 95% within +/- l km). During this period, a has varied within 7,714,429 +/- 7 m, while the inclination i periodically fluctuated in the range 66.0408 deg. +/- 0.0040 deg. The frozen orbit (required e < 0.001 and omega approximately equals to 90 deg.) has been maintained without any dedicated eccentricity maneuvers. The frozen eccentricity vector has completed two periodic cycles and it is currently tracing its third cycle (period approximately equals 26.7 months).

Bhat, R. S.↗
DOE OSTI2024

Rephasing spectral diffusion in time-bin spin-spin entanglement protocols

Generating high-fidelity spin-spin entanglement is an essential task of quantum repeater networks for the distribution of quantum information across long distances. Solid-state-based spin-photon interfaces are promising candidates to realize nodes of a quantum network, but are often limited by spectral diffusion of the optical transition, which results in phase errors on the entangled states. Here, we introduce a method to correct phase errors from quasi-static frequency fluctuations after the entangled state is generated, by shelving the emitters in the excited state to refocus the unknown phase. For quasi-static frequency fluctuations, the fidelity is determined only by the lifetime of the excited state used for shelving, making it particularly suitable for systems with a long-lived shelving state with correlated spectral diffusion. Such a shelving state may be found in Kramers doublet systems such as rare-earth emitters and color centers in Si or SiC interfaced with nanophotonic cavities with a strongly frequency-dependent Purcell enhancement. Furthermore, the protocol can be used to generate high-fidelity entangled spin pairs without reducing the rate of entanglement generation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Particle-number distribution in large fluctuations at the tip of branching random walks

Here, we investigate properties of the particle distribution near the tip of one-dimensional branching random walks at large times t , focusing on unusual realizations in which the rightmost lead particle is very far ahead of its expected position, but still within a distance smaller than the diffusion radius ~$\sqrt{t}$. Our approach consists in a study of the generating function $G_{Δx}(λ) = Σ_n$ ${λ^n}p_n(Δx)$ for the probabilities $p_n(Δx)$ of observing $\textit{n}$ particles in an interval of given size $Δ\textit{x}$ from the lead particle to its left, fixing the position of the latter. This generating function can be expressed with the help of functions solving the Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation with suitable initial conditions. In the infinite-time and large-$Δ\textit{x}$ limits, we find that the mean number of particles in the interval grows exponentially with $Δ\textit{x}$, and that the generating function obeys a nontrivial scaling law, depending on $Δ\textit{x}$ and λ through the combined variable $[Δx — f(λ)]^3 / Δx^2$, where $\textit{f}$(λ) ≡ – ln(1 – λ) – ln [– ln(1 – λ)]. From this property, one may conjecture that the growth of the typical particle number with the size of the interval is slower than exponential, but, surprisingly enough, only by a subleading factor at large Δ$\textit{x}$. The scaling we argue is consistent with results from a numerical integration of the FKPP equation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1999

Dynamics of Single Chains of Suspended Ferrofluid Particles

We present an experimental study of the dynamics of isolated chains made of super-paramagnetic particles under the influence of a magnetic field. The motivation of this work is to understand if the chain fluctuations exist and, if it does, how does the fluctuation affect chain aggregation. We find that single chains strongly fluctuate and that the characteristic frequency of their fluctuations is inversely proportional to the magnetic field strength. The higher the field the lower the characteristic frequency of the chain fluctuations. In the high magnetic field limit, chains behave like rigid rods without any internal motions. In this work, we used ferrofluid particles suspended in water. These particles do not have any intrinsic magnetization. Once a magnetic field is applied, a dipole moment is induced in each particle, proportional to the magnetic field. A dipolar magnetic interaction then occurs between particles. If dipole-dipole magnetic energy is higher than the thermal energy, the result is a structure change inside the dipolar fluid. The ratio of these two energies is expressed by a coupling constant lambda as: lambda = (pi(a(exp 3))(chi(exp 2))(mu(sub 0))(H(sub 0))(exp 2))/18kT Where a is the particle radius, mu(sub 0) is the vacuum magnetic permeability, H(sub 0) the applied magnetic field, k the Boltzmann constant and T the absolute temperature. If lambda > 1, magnetic particles form chains along the field direction. The lateral coalescence of several chains may form bigger aggregates especially if the particle volume fraction is high. While many studies and applications deal with the rheological properties and the structural changes of these dipolar fluids, this work focuses on the understanding of the chain dynamics. In order to probe the chain dynamics, we used dynamic light scattering (DLS) in self-beating mode as our experimental technique. The experimental geometry is such that the scattering plane is perpendicular to the magnetic field. Therefore, only motions in this plane are probed. A very dilute sample of a ferrofluid emulsion with a particle volume fraction of 10(exp -5) is used in this experiment. We chose such a low volume fraction to avoid multiple light scattering as well as lateral chain-chain aggregation. DLS measures the dynamic structure factor S(q,t) of the sample (q is the scattering wave vector, t is the time). In the absence of the magnetic field, identical particles of ferrofluid droplets are randomly distributed and S(q,t) reduces to exp(-q(exp 2)2D(sub 0)t). D(sub 0)=(kT/(6(pi)(eta)(a)) is the diffusion coefficient of Brownian particles (where Xi = (6(pi)(eta)(a)) is the Stokes frictional coefficient of a spherical particle in a fluid of viscosity eta). If interactions or polydispersity can not be ignored, an effective diffusion coefficient is introduced. Formally, D(sub eff) is defined as: D(sub eff) = - q(exp -2) partial derivative of (ln(S(q,t)) with respect to time, as t goes to 0. D(sub eff) reduces to D(sub 0) if no interactions and only a few particles size are present. Therefore, we can use DLS to measure particle size. The particle radius was found to be a=0.23 mu m with 7% of polydispersity. In this case, if we vary the scattering angle theta (and so q) we do not have any change in the measured diffusion coefficient: it is q-independent. When a magnetic field is applied, particles aggregate into chains if lambda > 1. We first studied the kinetics of the chain formation when lambda = 406. At a fixed scattering angle, we measured diffusion coefficient D(sub eff) as a function of time. Experimentally, we find that D(sub eff) decreases monotonously with time. Physically, this means that chains are becoming longer and longer. Since we are only sensitive to motions in the scattering plane and since chains have their main axis perpendicular to this plane, the measured diffusion coefficient is the trans-verse diffusion coefficient. We can relate D(sub eff) to the mean number of particles per chain N(t) at a given time and to the diffusion coefficient of an isolated particle D(sub 0) as D(sub eff)=f(N(t))D(sub 0). Since f(N) is known from other recent work, N can be expressed as a function of the time. We found a square root dependency: N(t) proportional to the square root of t. As expected for very low volume fraction, this behavior is characteristic of a diffusion-limited aggregation as suggested by several authors and by our previous work. In this study, we focus on the dependence of the effective diffusion coefficient on the scattering angle and the magnetic field strength. After the magnetic field is applied (lambda = 406) for a long time, typically 6 hours, kinetics of chain formation becomes very slow. Chain size does not vary much over the next hour period. Thus, we can perform different interesting experiments. First, at a fixed magnetic field, we measure the effective diffusion coefficient as a function of the scattering angle (from 5 to 130 deg). Our results show that the measured diffusion coefficient increases linearly with the scattering angle: D(sub eff) proportional to q. If we do the same experiment for different lambda values, D(sub eff) depends on lambda as D(sub eff) proportional to lambda(exp -1/2). We also find for different lambda values that the same asymptotic D(sub eff) value is obtained when q approaches zero. The angle dependency of D(sub eff) suggests that an additional motion exists besides chain drifting. Chain size is constant during experiment, which was verified by measuring the same diffusion coefficient at the beginning and at the end of the angle switching. If chains are rigid, D(sub eff) is independent of q. Therefore, we found that D(sub eff) not only measures the motion of the entire chain but also its internal fluctuations. These internal motions are the fluctuations of the particles in the chain. To understand the q dependency of D(sub eff), let us look at the probing length used. In our study, the characteristic length scale probed is l=2pi/q which is in the range of 0.9<l/a<20. When l is much larger than the particle radius (l/a going to infinity) we are mainly sensitive to the center of mass diffusion of the chain. D(sub eff) is then the diffusion of the entire chain and depends only on N but not on q and lambda. The value D(sub eff) allows us to obtain the number of particles per chain N. In the opposite limit (i.e. l/a < 1), we are sensitive to motions on the size of individual particles. The main contributions to the measured diffusion coefficient come from internal motions of the chains (i.e. particles' fluctuations inside the chain). We investigated also the effect of the chain size on D(sub eff). For the same value of the magnetic interaction lambda and different chain sizes, we found that the slope of D(sub eff) versus q decreases when N increases as a power law D(sub eff) proportional to N(exp -0.7). This is an important behavior because this means that the longer the chain the slower its internal fluctuations besides the whole chain motion.

Cutillas, S.↗
DOE OSTI2021

Relaxation timescales and electron-phonon coupling in optically pumped $\mathrm{YBa_2Cu_3O_{6+\mathcal{x}}}$ revealed by time-resolved Raman scattering

Time-resolved measurements provide a new way to disentangle complex interactions in quantum materials due to their different timescales. We used pump-probe Raman scattering to investigate the apical oxygen vibration in $\mathrm{YBa_2Cu_3O_{6+\mathcal{x}}}$ under nonequilibrium conditions. Time dependence of the phonon population demonstrated strong electron-phonon coupling. Most importantly, the phonon shifts to a higher energy due to transient smearing of the Fermi surface in a remarkable agreement with diagrammatic theory. Here we also discuss insights into photoinduced superconductivity reported at lower doping that follow from these results.

36 MATERIALS SCIENCE↗
NASA NTRS1983

Ozone transport by diabatic and planetary wave circulations on a beta plane

A beta plane model is used to study the transport of ozone and a tracer distributed like ozone. In the first experiment, ozone is transported with a time dependent diabatic circulation in a beta channel. Increase in the ozone column density occurs principally through downward advection from the high altitude source region into the lower stratosphere. Compared with observations, the maximum in the total ozone field occurs too late and is too large. Transient planetary waves produce large temporary increases in total ozone near the pole during sudden warming events, but permanent increases are much smaller. Warmings which penetrate below the chemical transition region (major warmings) are most effective in changing the zonal mean ozone distribution. Combining the sudden warming and diabatic circulations simulates many of the observed features of the seasonal changes in total ozone at high latitudes. The diabatic circulation provides the major enhancement of the ozone column density while the transient planetary waves move ozone to high latitudes. This result suggests sudden warmings can play an important role in our understanding of the fluctuations in total ozone, and the final warming may be the key event which determines the date and magnitude of the spring maximum.

Rood, R. B.↗
DOE OSTI2022

Probing eddy size and its effective mixing length in stably stratified roughness sublayer flows

Stably stratified roughness sublayer flows are ubiquitous yet remain difficult to represent in models and to interpret using field experiments. Here, continuous high-frequency potential temperature profiles from the forest floor up to 6.5 times the canopy height observed with distributed temperature sensing (DTS) are used to link eddy topology to roughness sublayer stability correction functions and coupling between air layers within and above the canopy. The experiments are conducted at two forest stands classified as hydrodynamically sparse and dense. Near-continuous profiles of eddy sizes (length scales) and effective mixing lengths for heat are derived from the observed profiles using a novel conditional sampling approach. The approach utilizes potential temperature isoline fluctuations from a statically stable background state. The transport of potential temperature by an observed eddy is assumed to be conserved (adiabatic movement) and we assume that irreversible heat exchange between the eddy and the surrounding background occurs along the (vertical) periphery of the eddy. This assumption is analogous to Prandtl's mixing-length concept, where momentum is transported rapidly vertically and then equilibrated with the local mean velocity gradient. A distinct dependence of the derived length scales on background stratification, height above ground, and canopy characteristics emerges from the observed profiles. Implications of these findings for (1) the failure of Monin–Obukhov similarity in the roughness sublayer and (2) above-canopy flow coupling to the forest floor are examined. The findings have practical applications in terms of analysing similar DTS data sets with the proposed approach, modelling roughness sublayer flows, and interpreting nocturnal eddy covariance measurements above tall forested canopies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Omnigenous umbilic stellarators

To better understand the dependence of the magnetic field structure in the plasma edge on the plasma boundary shape, in the context of X-point and island divertor designs, we define and develop a class of stellarators called umbilic stellarators. These equilibria are characterised by a single continuous high-curvature edge on the plasma boundary that goes around multiple times toroidally before meeting itself. We develop a technique that allows us to simultaneously optimise the plasma boundary along with a curve lying on the boundary on which we impose a high curvature while imposing omnigenity – a property of the magnetic field that ensures trapped particle confinement throughout the plasma volume. We find that umbilic stellarators naturally tend to favour piecewise omnigenity instead of omnigenity with a specific helicity. After generating omnigenous umbilic stellarators, we design coil sets for some of them and explore the fieldline structure in the edge and its sensitivity to small fluctuations in the plasma. Finally, using single-stage optimisation, we simultaneously modify the plasma and coil shape and propose an experiment to modify an existing tokamak to a finite-𝛽 stellarator using this technique and explore a potentially simpler way to convert a limited tokamak into a diverted stellarator.

fusion plasma↗