Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse power flow”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Dependence of wind-farm-induced gravity waves and wind farm performance on non-dimensional atmospheric parameters and simulation configuration

This large-eddy simulation (LES) study examines how wind-farm-induced atmospheric gravity waves (AGWs) and wind farm performance depend on non-dimensional atmospheric parameters and simulation configuration. A hypothetical aligned wind farm of actuator disks is simulated under neutral surface conditions, with a stable capping inversion and a mildly stable free atmosphere, to assess the effects of stratification beyond the atmospheric boundary layer (ABL) on ABL flow. Simulation set-ups fully resolving AGWs are validated to minimize spurious wave generation and reflection from the domain boundaries. The validated set-up is then used to analyze AGW types and characteristics, as well as stratification impacts under conventionally neutral boundary layer (CNBL) conditions. These conditions are governed by four non-dimensional parameters: the Froude numbers of the free atmosphere and capping inversion (Fr, Fr i ), and the aspect ratios of the ABL and wind farm (H̃ i , S h ). Simulation configurations that fully resolve AGWs – capturing at least one wavelength both horizontally and vertically – yield the most realistic stratification effects on ABL flow, whereas partial or unresolved configurations produce non-physical, channel-like behavior. A coherent description of the AGW phenomena is provided, highlighting the central role of capping inversion displacement in linking ABL fluctuations with AGWs. Trapped waves are confined within the capping inversion, while interfacial and internal waves aloft are identified as the AGW types most relevant to wind farm performance. The wavy inversion, analogous to an interfacial wave, forms converging and diverging zones that drive power fluctuations across the farm. The interfacial wavelength, measured over the wind farm, corresponds to one diverging, one converging, and one mildly diverging zone. As the interfacial wavelength decreases with Fr i , multiple convergence–divergence zones develop under sub-critical conditions (Fr i <1.0), while for super-critical conditions (Fr i > 1.0), the wavelength approaches the farm length. Wave amplitude increases with decreasing H̃ i (i.e., shallower capping inversions). Wind farm performance is most sensitive to H̃ i : shallow boundary layers increase blockage and reduce efficiency, while deeper layers enhance efficiency. Increasing Fr and Fr i mitigates blockage, and increasing S h mainly improves wake recovery. Although local power fluctuations arise from AGWs, overall wind farm efficiency remains nearly constant with Fr and Fr i , improving primarily with larger H̃ i and S h .

17 WIND ENERGY↗
DOE OSTI2026

Alignment switching in 3D-printed smectic liquid crystal elastomers

Extrusion-based additive manufacturing has emerged as a powerful platform for designing shape-morphing materials through controlled orientation. However, existing approaches primarily rely on a single mode of flow-induced alignment, limiting orientation programmability. Herein, we present a direct-ink-writing approach for smectic liquid crystal elastics that exploits two distinct alignment modes within a single ink. The smectic ink exhibits shear- and temperature-dependent orientation switching, enabling molecular alignment either perpendicular or parallel to the print direction. Combined rheological, X-ray, and molecular dynamics analyses reveal that this alignment inversion arises from the preservation or collapse of smectic layers under flow. This reversible switching encodes both contractile and elongational actuation within individual filaments, greatly expanding the design freedom of printed liquid crystal elastomers. We demonstrate 2D and 3D structures with diverse programmed shape transformations, highlighting the potential of this platform for adaptive soft actuators and architected functional materials.

Lee, Jin Hyeong [Pusan National University, Busan,↗
DOE OSTI2021

A unified theory of zonal flow shears and density corrugations in drift wave turbulence

A unified theory of zonal flow shears and density corrugations in drift wave turbulence is presented. Polarization and density advection beat excitation are studied in combination with modulational response. Noise is driven by two-time flux correlation. While the effective zonal flow eddy viscosity can go negative, the zonal diffusivity is positive definite. There is no inverse cascade of density corrugation. The connection between avalanches and corrugations is discussed. Here, the zonal cross-correlation is identified and calculated. Conditions for alignment of zonal shears and corrugation gradients are determined, and the implications for staircase structure are discussed. We show that the synergy of beat noise and modulational effects is stronger than either alone. Strong zonal flows can be excited well below the modulational instability threshold. In the context of L–H transition, zonal noise quenches turbulence overshoot by eliminating the threshold for zonal flow excitation. The power threshold for L–H transition is lowered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Mechanisms of Dynamic Near-Wake Modulation of a Utility-Scale Wind Turbine

The current study uses large eddy simulations to investigate the transient response of a utility-scale wind turbine wake to dynamic changes in atmospheric and operational conditions, as observed in previous field-scale measurements. Most wind turbine wake investigations assume quasi-steady conditions, but real wind turbines operate in a highly stochastic atmosphere, and their operation (e.g. blade pitch, yaw angle) changes constantly in response. Furthermore, dynamic control strategies have been recently proposed to optimize wind farm power generation and longevity. Therefore, improved understanding of dynamic wake behaviours is essential. First, changes in blade pitch are investigated and the wake expansion response is found to display hysteresis as a result of flow inertia. The time scales of the wake response to different pitch rates are quantified. Next, changes in wind direction with different time scales are explored. Under short time scales, the wake deflection is in the opposite direction of that observed under quasi-steady conditions. Finally, yaw changes are implemented at different rates, and the maximum inverse wake deflection and time scale are quantified, showing a clear dependence on yaw rate. To gain further physical understanding of the mechanism behind the inverse wake deflection, the streamwise vorticity in different parts of the wake is quantified. The results of this study provide guidance for the design of advanced wake flow control algorithms. The lag in wake response observed for both blade pitch and yaw changes shows that proposed dynamic control strategies must implement turbine operational changes with a time scale of the order of the rotor time scale or slower.

flow-structure interactions↗
DOE OSTI2021

Multi-Fidelity Active Subspaces for Wind Farm Uncertainty Quantification

Wind plants operate in stochastic environments characterized by complex turbulent flow dynamics and high-dimensional random variables. A key step in uncertainty quantification studies is sensitivity analysis and dimension reduction that can facilitate the development of surrogate models to be used for forward and inverse propagation or optimization under uncertainty. Prior work has shown active subspaces are an effective tool for identifying important directions in the space of stochastic inputs; however, they have only been applied to single-fidelity wind plant models. In this study, we investigate the efficacy of a multi-fidelity active subspace method for analyzing the uncertainty in wind plant power output. The multi-fidelity active subspace estimator offers the promise of increased accuracy in identifying active subspaces as compared to a single-fidelity estimator for the same computational cost, or a reduction in cost for the same accuracy. This makes the study of uncertainty in larger wind plants and with higher fidelity physics tractable. The multi-fidelity active subspace method is applied to gridded and existing wind plant layouts with single and multiple inflow conditions and its performance for surrogate modeling and uncertainty propagation is compared against a single-fidelity active subspace method. This multi-fidelity approach yields substantial computational speedups of 2x - 3.4x across the test cases along with acceptable accuracy in surrogate modeling and computing statistical moments.

active subspace↗
DOE OSTI2023

Dark matter halo mass functions and density profiles from mass and energy cascade

Abstract Halo abundance and structure play a central role for modeling structure formation and evolution. Without relying on a spherical or ellipsoidal collapse model, we analytically derive the halo mass function and cuspy halo density (inner slope of −4/3) based on the mass and energy cascade theory in dark matter flow. The hierarchical halo structure formation leads to halo or particle random walk with a position-dependent waiting time $$\tau _g$$ τ g . First, the inverse mass cascade from small to large scales leads to the halo random walk in mass space with $$\tau _g\propto m_h^{-\lambda }$$ τ g ∝ m h - λ , where $$m_h$$ m h is the halo mass and $$\lambda$$ λ is a halo geometry parameter with predicted value of 2/3. The corresponding Fokker-Planck solution for halo random walk in mass space gives rise to the halo mass function with a power-law behavior on small scale and exponential decay on large scale. This can be further improved by considering two different $$\lambda$$ λ for haloes below and above a critical mass scale $$m_h^*$$ m h ∗ , i.e. a double- $$\lambda$$ λ halo mass function. Second, a double- $$\gamma$$ γ density profile can be derived based on the particle random walk in 3D space with a position-dependent waiting time $$\tau _g \propto \Phi (r)^{-1} \propto r^{-\gamma }$$ τ g ∝ Φ ( r ) - 1 ∝ r - γ , where $$\Phi$$ Φ is the gravitational potential and r is the particle distance to halo center. Theory predicts $$\gamma =2/3$$ γ = 2 / 3 that leads to a cuspy density profile with an inner slope of −4/3, consistent with the predicted scaling laws from energy cascade. The Press-Schechter mass function and Einasto density profile are just special cases of proposed models. The small scale permanence can be identified due to the scale-independent rate of mass and energy cascade, where density profiles of different halo masses and redshifts converge to the $$-4/3$$ - 4 / 3 scaling law ( $$\rho _h \propto r^{-4/3}$$ ρ h ∝ r - 4 / 3 ) on small scales. Theory predicts the halo number density scales with halo mass as $$\propto m_h^{-1.9}$$ ∝ m h - 1.9 , while the halo mass density scales as $$\propto m_h^{4/9}$$ ∝ m h 4 / 9 . Results were compared against the Illustris simulations. This new perspective provides a theory for nearly universal halo mass functions and density profiles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Experimental demonstration of transient CHI start-up using a floating biased electrode configuration

Transient coaxial helicity injection (CHI) current start-up in a spherical tokamak (ST) has been demonstrated for the first time using a reactor-relevant floating single biased (FSB) electrode configuration on QUEST. Unlike previous applications of CHI on HIT-II and on NSTX in which the entire ST vessel was split in two using two toroidal ceramic rings as vacuum breaks, in the FSB electrode configuration on QUEST, only the lower outer divertor plate is electrically insulated from the rest of the vessel. This is the CHI injector region. A CHI discharge is generated by driving current, using an external power supply, along magnetic field lines that connect the inner and outer divertor plates in the CHI injector region. The CHI plasma evolution and closed flux formation on QUEST are very similar to those on NSTX and have generated toroidal currents of 150 kA with 100 kA of current flowing on closed flux surfaces. An alternate, force-free model approach has been used to derive the CHI scaling relations. The resulting CHI scaling relations from the model are consistent with the conventional CHI scaling relations, which state that the CHI generated toroidal current is proportional to the CHI injector flux and the CHI injector current is inversely proportional to the toroidal field. These results from QUEST bode well for the application of CHI in a new generation of compact high field STs and low aspect ratio tokamaks in which the space for the central solenoid is very restricted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Inverse energy transfer in decaying, three-dimensional, non-helical magnetic turbulence due to magnetic reconnection

ABSTRACT It has been recently shown numerically that there exists an inverse transfer of magnetic energy in decaying, non-helical, magnetically dominated, magnetohydrodynamic turbulence in three dimensions (3D). We suggest that magnetic reconnection is the underlying physical mechanism responsible for this inverse transfer. In the two-dimensional (2D) case, the inverse transfer is easily inferred to be due to smaller magnetic islands merging to form larger ones via reconnection. We find that the scaling behaviour is similar between the 2D and 3D cases, i.e. the magnetic energy evolves as t−1, and the magnetic power spectrum follows a slope of k−2. We show that on normalizing time by the magnetic reconnection time-scale, the evolution curves of the magnetic field in systems with different Lundquist numbers collapse on to one another. Furthermore, transfer function plots show signatures of magnetic reconnection driving the inverse transfer. We also discuss the conserved quantities in the system and show that the behaviour of these quantities is similar between the 2D and 3D simulations, thus making the case that the dynamics in 3D could be approximately explained by what we understand in 2D. Lastly, we also conduct simulations where the magnetic field is subdominant to the flow. Here, too, we find an inverse transfer of magnetic energy in 3D. In these simulations, the magnetic energy evolves as t−1.4 and, interestingly, a dynamo effect is observed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Acceleration mechanisms of energetic ion debris in laser-driven tin plasma EUV sources

Laser-driven tin plasmas are driving new-generation nanolithography as sources of extreme ultraviolet (EUV) radiation centered at 13.5 nm. A major challenge facing industrial EUV source development is predicting energetic ion debris produced during the plasma expansion that may damage the sensitive EUV channeling multilayer optics. Gaining a detailed understanding of the plasma dynamics and ion acceleration mechanisms in these sources could provide critical insights for designing debris mitigation strategies in future high-power EUV sources. Here, we develop a fully kinetic model of tin-EUV sources using one-dimensional particle-in-cell simulations to study ion debris acceleration, which will be valuable for cross-validation of radiation-hydrodynamic simulations. An inverse-bremsstrahlung heating operator is used to model the interaction of a tin target with an Nd:YAG laser, and thermal conduction is included through a Monte Carlo Coulomb collision operator. While the large-scale evolution is in reasonable agreement with analogous hydrodynamic simulations, the significant timescale for collisional equilibration between electrons and ions allows for the development of prominent two-temperature features. A collimated flow of energetic ions is produced with a spectrum that is significantly enhanced at high energies compared to fluid simulations. The dominant acceleration mechanism is found to be a large-scale electric field supported mainly by the electron pressure gradient, which is enhanced in the kinetic simulations due to the increased electron temperature. We discuss the implications of these results for future modeling of tin-EUV sources and the development of debris mitigation schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Control Design for a Marine Hydrokinetic Cycloturbine Vehicle

Abstract Marine Hydrokinetic (MHK) cycloturbines generate sustainable power by exploiting tidal currents. By powering the turbines and using pitching foils for control, a vehicle comprised of MHK cycloturbines also has the ability to station keep and maneuver. The vehicle consists of four counter-rotating cycloturbines, with hydrofoils oriented perpendicular to the flow in a paddlewheel configuration. Lift and drag generated from these foils sum together to produce thrust. An experimentally tuned simulation model that solves the six-degrees-of-freedom rigid body equations of motion for the MHK vehicle subject to hydrodynamic, hydrostatic, and propulsive forces is used to aid the design of vehicle controllers. Global feedback controllers are initially designed by applying classical control methods to an approximate linear model of the system dynamics. A higher performing nonlinear controller is designed using the nonlinear dynamic inversion (NDI) method. NDI accounts for the nonlinearities of the MHK system and therefore is suitable for a wide range of operating conditions. The response of the classical and NDI controllers to speed, depth, roll, pitch, and yaw commands are evaluated and compared in simulation. The classical controller outperforms the NDI controller for small amplitude maneuvers, although the degradation with NDI is minor. However, in the nonlinear operating regime the NDI controller outperforms the classical controller and the classical controller exhibits instability.

Automation & Control Systems↗
DOE OSTI2023

Evolution and feedback of AGN jets of different cosmic ray composition

Jet feedback from active galactic nuclei (AGN) is one of the most promising mechanisms for suppressing cooling flows in cool-core clusters. However, the composition of AGN jets and bubbles remains uncertain; they could be thermally dominated, or dominated by cosmic ray proton (CRp), cosmic ray electron (CRe), or magnetic energy. In this work, we investigate the evolution and feedback effects of CRp and CRe dominated jets by conducting 3D magnetohydrodynamic simulations of AGN jet-inflated bubbles in the intracluster medium using the FLASH code. We present the evolution of their energies, dynamics, and heating, and model their expected cavity-power versus radio-luminosity relation (P cav –L R ). We find that bubbles inflated by CRe dominated jets follow a very similar dynamical evolution to CRp dominated bubbles even though CRe within bubbles suffer significantly stronger synchrotron and inverse-Compton cooling. This is because, as CRe lose their energy, the jet-inflated bubbles quickly become thermally dominated within ∼30 Myr. Their total energy stops decreasing with CR energy and evolves similarly to CRp dominated bubbles. The ability of CRe and CRp dominated bubbles to heat the intracluster medium is also comparable; the cold gas formed via local thermal instabilities is well suppressed in both cases. The CRp and CRe bubbles follow different evolutionary trajectories on the P cav –L R plane, but the values are broadly consistent with observed ranges for FR-I sources. We also discuss observational techniques that have potential for constraining the composition of AGN jets and bubbles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Laser-Plasma Spatiotemporal Cyanide Spectroscopy and Applications

This article reports new measurements of laser-induced plasma hypersonic expansion measurements of diatomic molecular cyanide (CN). Focused, high-peak-power 1064 nm Q-switched radiation of the order of 1 TW/cm 2 generated optical breakdown plasma in a cell containing a 1:1 molar gas mixture of N 2 and CO 2 at a fixed pressure of 1.1 × 10 5 Pascal and in a 100 mL/min flow of the mixture. Line-of-sight (LOS) analysis of recorded molecular spectra indicated the outgoing shockwave at expansion speeds well in excess of Mach 5. Spectra of atomic carbon confirmed increased electron density near the shockwave, and, equally, molecular CN spectra revealed higher excitation temperature near the shockwave. Results were consistent with corresponding high-speed shadowgraphs obtained by visualization with an effective shutter speed of 5 nanoseconds. In addition, LOS analysis and the application of integral inversion techniques allow inferences about the spatiotemporal plasma distribution.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electrochemically driven carbon dioxide separation

This project explored the viability of a Ni(OH) 2 based hydroxide exchange membrane carbon capture (HEMCC) device for direct air capture. It is built off an H 2 fuel cell based HEMCC previously developed in the Yushan Yan group at University of Delaware. Taking the same membrane-based, electrochemically driven pH gradient concept, similar flux performance was able to be seen in the Ni(OH) 2 system as the H 2 system when using comparable current densities. The project produced two types of membrane electrode assemblies (MEA). The first was a traditional MEA with two electrodes and a membrane, the second was a flow-through membrane MEA. Consistent performance was achieved with the traditional MEA with an energy cost of 1.1 MWh∙ton -1 and flux of 82 kg∙m -2 ∙yr -1 . This was the most stable of the two designs. This project investigated strategies to improve performance with a flow-through membrane design. Two designs were made, one with a cast phase inversion membrane, and one with powdered membrane layer. The template phase inversion membrane achieved low pressure drop but had performance limitations due to a skin layer of membrane limiting CO 2 gas transport. The powdered membrane layer had better performance but higher pressure drop. Finally, using the powdered membrane, commercial battery materials were able to be used in order to achieve higher flux. This showed that the flow through membrane design does have the capability to overcome flux inefficiencies in the traditional MEA. Moreover, the process design of the system is proposed and given in this project. The process mass and energy balances were calculated for a reference plant of 1000 t/yr CO 2 capture, which contains several subsystems, e.g. air processing subsystem, electrical subsystem, and CO 2 purification and compression subsystem is designed and evaluated. Based on our calculation, When the power consumption for stack is 1 MWh t -1 CO 2 , Additional 382 kWh t -1 CO 2 will be consumed by other subsystem of the plant, i.e. The total power consumption of 1.38 MWh t -1 CO 2 , lower than the final milestone of 1.5 MWh t -1 CO 2 in this project.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

General relativistic MHD simulations of non-thermal flaring in Sagittarius A*

Sgr A* exhibits regular variability in its multiwavelength emission, including daily X-ray flares and roughly continuous near-infrared (NIR) flickering. The origin of this variability is still ambiguous since both inverse Compton and synchrotron emission are possible radiative mechanisms. The underlying particle distributions are also not well constrained, particularly the non-thermal contribution. In this work, we employ the GPU-accelerated general relativistic magnetohydrodynamics code H-AMR to perform a study of flare flux distributions, including the effect of particle acceleration for the first time in high-resolution 3D simulations of Sgr A*. For the particle acceleration, we use the general relativistic ray-tracing code bhoss to perform the radiative transfer, assuming a hybrid thermal+non-thermal electron energy distribution. We extract ~60 h light curves in the sub-millimetre, NIR and X-ray wavebands, and compare the power spectra and the cumulative flux distributions of the light curves to statistical descriptions for Sgr A* flares. Our results indicate that non-thermal populations of electrons arising from turbulence-driven reconnection in weakly magnetized accretion flows lead to moderate NIR and X-ray flares and reasonably describe the X-ray flux distribution while fulfilling multiwavelength flux constraints. These models exhibit high rms per cent amplitudes, $\gtrsim 150{{\ \rm per\ cent}}$ both in the NIR and the X-rays, with changes in the accretion rate driving the 230 GHz flux variability, in agreement with Sgr A* observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Translation and rotation equivariant normalizing flow (TRENF) for optimal cosmological analysis

ABSTRACT Our Universe is homogeneous and isotropic, and its perturbations obey translation and rotation symmetry. In this work, we develop translation and rotation equivariant normalizing flow (TRENF), a generative normalizing flow (NF) model which explicitly incorporates these symmetries, defining the data likelihood via a sequence of Fourier space-based convolutions and pixel-wise non-linear transforms. TRENF gives direct access to the high dimensional data likelihood p(x|y) as a function of the labels y, such as cosmological parameters. In contrast to traditional analyses based on summary statistics, the NF approach has no loss of information since it preserves the full dimensionality of the data. On Gaussian random fields, the TRENF likelihood agrees well with the analytical expression and saturates the Fisher information content in the labels y. On non-linear cosmological overdensity fields from N-body simulations, TRENF leads to significant improvements in constraining power over the standard power spectrum summary statistic. TRENF is also a generative model of the data, and we show that TRENF samples agree well with the N-body simulations it trained on, and that the inverse mapping of the data agrees well with a Gaussian white noise both visually and on various summary statistics: when this is perfectly achieved the resulting p(x|y) likelihood analysis becomes optimal. Finally, we develop a generalization of this model that can handle effects that break the symmetry of the data, such as the survey mask, which enables likelihood analysis on data without periodic boundaries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Photoredox Catalysis with Spin Magnetic Field Effects: A Scheme for Chiral Resolution

Quantities that break both mirror symmetry and time-reversal symmetry, such as the orbital angular momentum, are known to connect molecular chirality with an applied magnetic field. This concept has led to observations such as magneto-chiral dichroism and chirality-induced spin selectivity (CISS). However, being small, these effects often require additional amplification procedures such as flow chemistry to achieve bulk enantioseparation. In this work, we demonstrate how the magnetic field effect on photogenerated radical pairs, which also breaks time-reversal symmetry, can be harnessed for enantiopurification. Fundamental to this process is the collective decay of the singlet and triplet radical-pair states made possible by an applied magnetic field. Because opposite enantiomers exhibit spin-orbit coupling matrix elements of opposite signs, the singlet and triplet decay channels interfere constructively in one enantiomer. Meanwhile, molecules of the other enantiomer are funnelled into the first enantiomer through excited-state chirality inversion, achieving (dynamic kinetic) chiral resolution. Using an axially chiral binaphthyl derivative and a borane photosensitizer as prototype, we predict an appreciable enantiomeric excess (e.e.) of 90% to be possible at steady state, attained within hundreds of milliseconds when irradiated by a laser. Importantly, our analytical results showcase regimes of perfect enantioselectivity (100% e.e.), accessible by further chemical optimization of the photosensitizer for which general strategies are discussed. Altogether, this work illustrates a so-far untapped but powerful control knob for photoredox catalysis based on spin chemistry principles.

Magnetic properties↗
DOE OSTI2023

Using Dark Fiber and Distributed Acoustic Sensing to Characterize a Geothermal System in the Imperial Valley, Southern California

The Imperial Valley, CA, is a tectonically active transtensional basin located south of the Salton Sea; the area hosts numerous geothermal fields, including significant hidden hydrothermal resources without surface manifestations. Development of inexpensive, rugged, and highly sensitive exploration techniques for undiscovered geothermal systems is critical for accelerating geothermal power deployment as well as unlocking a low-carbon energy future. We present a case study utilizing distributed acoustic sensing (DAS) and ambient noise interferometry for geothermal reservoir imaging, utilizing unlit fiber-optic telecommunication infrastructure (dark fiber). The study exploits two days of passive DAS data acquired in early November 2020 over a ~28-km section of fiber from Calipatria, CA to Imperial, CA. We apply ambient noise interferometry to retrieve coherent signals from DAS records and develop a bin stacking technique to attenuate the effects from persistent localized noise sources and to enhance retrieval of coherent surface waves. As a result, we are able to obtain high-resolution two-dimensional (2D) S wave velocity ($V_s$) structure to 3 km depth, based on joint inversion of both the fundamental and higher overtones. We observe a previously unmapped high $V_s$ and low $V_p$/$V_s$ ratio feature beneath the Brawley geothermal system, which we interpret to be a zone of hydrothermal mineralization and lower porosity. This interpretation is consistent with a host of other measurements including surface heat flow, gravity anomalies, and available borehole wireline data. These results demonstrate the potential utility of DAS deployed on dark fiber for geothermal system exploration and characterization in the appropriate geological settings.

Vp/Vs imaging↗
DOE OSTI2020

Inverse Compton Cooling in the Coronae of Simulated Black Hole Accretion Flows

Here we present a formulation for a local cooling function to be employed in the diffuse, hot corona region of 3D GRMHD simulations of accreting black holes. This new cooling function calculates the cooling rate due to inverse Compton scattering by considering the relevant microphysics in each cell in the corona and approximating the radiation energy density and Compton temperature thereby integrating over the thermal seed photon flux from the disk surface. The method either assumes the ion and electron temperatures are equal (1T) or calculates them separately (2T) using an instantaneous equilibrium approach predicated on the actual relevant rate equations (Coulomb and Compton). The method is shown to be consistent with a more detailed ray-tracing calculation where the bulk of the cooling occurs, but is substantially less costly to perform. As an example, we apply these methods to a harm3d simulation of a 10M⊙, non-spinning black hole, accreting at nominally 1% the Eddington value. This new approach leads to radiative efficiency values sime65% above Novikov–Thorne, with a larger fraction of total cooling in the corona as compared to simulations performed using the original target-temperature cooling function. Time-averaged post-processing reveals that the continuum spectral observations predicted from these simulations are qualitatively similar to actual X-ray binary data, especially so for the 1T approach, which yields a harder power-law component (Γ = 2.25) compared to the 2T version (Γ = 2.53).

79 ASTRONOMY AND ASTROPHYSICS↗