Results for “TIME SCALE” — Page 34

NASA NTRS2016

Investigation of Damping Physics and CFD Tool Validation for Simulation of Baffled Tanks at Variable Slosh Amplitude

Determination of slosh damping is a very challenging task as there is no analytical solution. The damping physics involves the vorticity dissipation which requires the full solution of the nonlinear Navier-Stokes equations. As a result, previous investigations were mainly carried out by extensive experiments. A systematical study is needed to understand the damping physics of baffled tanks, to identify the difference between the empirical Miles equation and experimental measurements, and to develop new semi-empirical relations to better represent the real damping physics. The approach of this study is to use Computational Fluid Dynamics (CFD) technology to shed light on the damping mechanisms of a baffled tank. First, a 1-D Navier-Stokes equation representing different length scales and time scales in the baffle damping physics is developed and analyzed. Loci-STREAM-VOF, a well validated CFD solver developed at NASA MSFC, is applied to study the vorticity field around a baffle and around the fluid-gas interface to highlight the dissipation mechanisms at different slosh amplitudes. Previous measurement data is then used to validate the CFD damping results. The study found several critical parameters controlling fluid damping from a baffle: local slosh amplitude to baffle thickness (A/t), surface liquid depth to tank radius (d/R), local slosh amplitude to baffle width (A/W); and non-dimensional slosh frequency. The simulation highlights three significant damping regimes where different mechanisms dominate. The study proves that the previously found discrepancies between Miles equation and experimental measurement are not due to the measurement scatter, but rather due to different damping mechanisms at various slosh amplitudes. The limitations on the use of Miles equation are discussed based on the flow regime.

Yang, H. Q.↗

NASA NTRS2018

Validity of Miles Equation in Predicting Propellant Slosh Damping in Baffled Tanks at Variable Slosh Amplitude

Determination of slosh damping is a very challenging task as there is no analytical solution. The damping physics involves the vorticity dissipation which requires the full solution of the nonlinear Navier-Stokes equations. As a result, previous investigations were mainly carried out by extensive experiments. A systematical study is needed to understand the damping physics of baffled tanks, to identify the difference between the empirical Miles equation and experimental measurements, and to develop new semi-empirical relations to better represent the real damping physics. The approach of this study is to use Computational Fluid Dynamics (CFD) technology to shed light on the damping mechanisms of a baffled tank. First, a 1-D Navier-Stokes equation representing different length scales and time scales in the baffle damping physics is developed and analyzed. Loci-STREAM-VOF, a well validated CFD solver developed at NASA MSFC, is applied to study the vorticity field around a baffle and around the fluid-gas interface to highlight the dissipation mechanisms at different slosh amplitudes. Previous measurement data is then used to validate the CFD damping results. The study found several critical parameters controlling fluid damping from a baffle: local slosh amplitude to baffle thickness (A/t), surface liquid depth to tank radius (d/R), local slosh amplitude to baffle width (A/W); and non-dimensional slosh frequency. The simulation highlights three significant damping regimes where different mechanisms dominate. The study proves that the previously found discrepancies between Miles equation and experimental measurement are not due to the measurement scatter, but rather due to different damping mechanisms at various slosh amplitudes. The limitations on the use of Miles equation are discussed based on the flow regime.

Yang, H. Q.↗

NASA NTRS2018

Validity of Miles Equation in Predicting Propellant Slosh Damping in Baffled Tanks at Variable Slosh Amplitude

Determination of slosh damping is a very challenging task as there is no analytical solution. The damping physics involves the vorticity dissipation which requires the full solution of the nonlinear Navier-Stokes equations. As a result, previous investigations were mainly carried out by extensive experiments. A systematical study is needed to understand the damping physics of baffled tanks, to identify the difference between the empirical Miles equation and experimental measurements, and to develop new semi-empirical relations to better represent the real damping physics. The approach of this study is to use Computational Fluid Dynamics (CFD) technology to shed light on the damping mechanisms of a baffled tank. First, a 1-D Navier-Stokes equation representing different length scales and time scales in the baffle damping physics is developed and analyzed. Loci-STREAM-VOF, a well validated CFD solver developed at NASA MSFC, is applied to study the vorticity field around a baffle and around the fluid-gas interface to highlight the dissipation mechanisms at different slosh amplitudes. Previous measurement data is then used to validate the CFD damping results. The study found several critical parameters controlling fluid damping from a baffle: local slosh amplitude to baffle thickness (A/t), surface liquid depth to tank radius (d/R), local slosh amplitude to baffle width (A/W); and non-dimensional slosh frequency. The simulation highlights three significant damping regimes where different mechanisms dominate. The study proves that the previously found discrepancies between Miles equation and experimental measurement are not due to the measurement scatter, but rather due to different damping mechanisms at various slosh amplitudes. The limitations on the use of Miles equation are discussed based on the flow regime.

Yang, H. Q.↗

DOE OSTI2023

Automated Symbolic Upscaling: 1. Model Generation for Extended Applicability Regimes

Abstract In porous media theory, upscaling techniques are fundamental to deriving rigorous Darcy‐scale models for flow and reactive transport in subsurface systems. Due to limitations in classical techniques, a number of ad hoc approaches have been proposed to address physical regimes in which reactive time scales are similar to, or faster than, diffusive time scales. In Part 1 of this two part series, we present a strategy for expanding the applicability of classical homogenization theory by generalizing the assumed closure form. We detail the implementation of this strategy on two reactive mass transport problems with moderately reactive physics. The strategy produces nontrivial homogenized models with emergent terms and effective parameters that couple reactive, diffusive, and advective transport. The differences in equation forms between the macroscopic and pore‐scale descriptions advise caution to further studies where the forms of macroscopic equations are assumed, as opposed to rigorously derived. Numerical validation is provided for each problem to show that the error estimates of homogenization theory are satisfied, and to justify the implemented strategy. In Part 2, the presented strategy is automated using symbolic computing to expedite its implementation.

Pietrzyk, Kyle↗

DOE OSTI2019

Auto-ignitive deflagration speed of methane (CH 4 ) blended dimethyl-ether (DME)/air mixtures at stratified conditions

Front propagation speeds from fully resolved unsteady one dimensional simulations with dimethyl-ether (DME)/methane (CH 4 )/air mixtures under engine relevant conditions are presented using complex kinetics and transport. Different time-scales of monochromatic inhomogeneities in DME concentration with varying DME/CH 4 blending ratios are simulated to unravel the fundamental aspects of auto-ignition and flame propagation under the influence of reactivity stratification. To understand the influence of different stratification time-scales on the flame-ignition interaction, two sets of conditions are simulated such that low temperature chemistry is present in only one of them. For a given amplitude of stratification, it is found that the instantaneous propagation speed is significantly affected by the level of CH 4 concentration in the binary fuel blend. Specifically, for cases with low temperature chemistry, at relatively smaller time-scales, the overall fluctuation in the instantaneous propagation speed is found to subside as the level of CH 4 concentration in the mixture is increased. However, for both sets of conditions, at comparatively larger time-scales, a rapid change in the instantaneous propagation speed is observed with an increase in the level of CH 4 concentration in the mixture. The intrinsic effects of stratification time-scales on the low temperature chemistry and the high temperature chemistry are further examined to assess the flame-ignition interaction. A displacement speed analysis is also carried out to elucidate the underlying combustion modes that are responsible for such a variation in flame response.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Modeling time-coincident ultrafast electron transfer and solvation processes at molecule-semiconductor interfaces

Kinetic models based on Fermi's Golden Rule are commonly employed to understand photoinduced electron transfer dynamics at molecule-semiconductor interfaces. Implicit in such second-order perturbative descriptions is the assumption that nuclear relaxation of the photoexcited electron donor is fast compared to electron injection into the semiconductor. This approximation breaks down in systems where electron transfer transitions occur on 100-fs time scale. Here, we present a fourth-order perturbative model that captures the interplay between time-coincident electron transfer and nuclear relaxation processes initiated by light absorption. The model consists of a fairly small number of parameters, which can be derived from standard spectroscopic measurements (e.g., linear absorbance, fluorescence) and/or first-principles electronic structure calculations. Insights provided by the model are illustrated for a two-level donor molecule coupled to both (i) a single acceptor level and (ii) a density of states (DOS) calculated for TiO{sub 2} using a first-principles electronic structure theory. These numerical calculations show that second-order kinetic theories fail to capture basic physical effects when the DOS exhibits narrow maxima near the energy of the molecular excited state. Overall, we conclude that the present fourth-order rate formula constitutes a rigorous and intuitive framework for understanding photoinduced electron transfer dynamics that occur on the 100-fs time scale.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Entanglement and correlations in fast collective neutrino flavor oscillations

Collective neutrino oscillations play a crucial role in transporting lepton flavor in astrophysical settings like supernovae and neutron star binary merger remnants, which are characterized by large neutrino densities. In these settings, simulations in the mean-field approximation show that neutrino-neutrino interactions can overtake vacuum oscillations and give rise to fast collective flavor evolution on time-scales t ∝ μ –1 , with μ proportional to the local neutrino density. In this work, we study the full out-of-equilibrium flavor dynamics in simple multi-angle geometries displaying fast oscillations in the mean field linear stability analysis. Focusing on simple initial conditions, we analyze the production of pair correlations and entanglement in the complete many-body-dynamics as a function of the number N of neutrinos in the system, for up to thousands of neutrinos. Similarly to simpler geometries with only two neutrino beams, we identify three regimes: stable configurations with vanishing flavor oscillations, marginally unstable configurations with evolution occurring on long time scales τ ≈ µ –1 $\sqrt{N}$, and unstable configurations showing flavor evolution on short time scales τ ≈ μ –1 log(N). Here, we present evidence that these fast collective modes are generated by the same dynamical phase transition which leads to the slow bipolar oscillations, establishing a connection between these two phenomena and explaining the difference in their time scales. We conclude by discussing a semi-classical approximation which reproduces the entanglement entropy at short to medium time scales and can be potentially useful in situations with more complicated geometries where classical simulation methods starts to become inefficient.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2014

Diffusion and spectroscopy of water and lipids in fully hydrated dimyristoylphosphatidylcholine bilayer membranes

Microscopic structure and dynamics of water and lipids in a fully hydrated dimyristoylphosphatidylcholine phospholipid lipid bilayer membrane in the liquid-crystalline phase have been analyzed with all-atom molecular dynamics simulations based on the recently parameterized CHARMM36 force field. The diffusive dynamics of the membrane lipids and of its hydration water, their reorientational motions as well as their corresponding spectral densities, related to the absorption of radiation, have been considered for the first time using the present force field. In addition, structural properties such as density and pressure profiles, a deuterium-order parameter, surface tension, and the extent of water penetration in the membrane have been analyzed. Molecular self-diffusion, reorientational motions, and spectral densities of atomic species reveal a variety of time scales playing a role in membrane dynamics. The mechanisms of lipid motion strongly depend on the time scale considered, from fast ballistic translation at the scale of picoseconds (effective diffusion coefficients of the order of 10{sup −5} cm{sup 2}/s) to diffusive flow of a few lipids forming nanodomains at the scale of hundreds of nanoseconds (diffusion coefficients of the order of 10{sup −8} cm{sup 2}/s). In the intermediate regime of sub-diffusion, collisions with nearest neighbors prevent the lipids to achieve full diffusion. Lipid reorientations along selected directions agree well with reported nuclear magnetic resonance data and indicate two different time scales, one about 1 ns and a second one in the range of 2–8 ns. We associated the two time scales of reorientational motions with angular distributions of selected vectors. Calculated spectral densities corresponding to lipid and water reveal an overall good qualitative agreement with Fourier transform infrared spectroscopy experiments. Our simulations indicate a blue-shift of the low frequency spectral bands of hydration water as a result of its interaction with lipids. We have thoroughly analyzed the physical meaning of all spectral features from lipid atomic sites and correlated them with experimental data. Our findings include a “wagging of the tails” frequency around 30 cm{sup −1}, which essentially corresponds to motions of the tail-group along the instantaneous plane formed by the two lipid tails, i.e., in-plane oscillations are clearly of bigger importance than those along the normal-to-the plane direction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

The Effect of Internal Variability on Ocean Temperature Adjustment in a Low-Resolution CESM Initial Condition Ensemble

Due to its large heat capacity and circulation, the ocean contributes significantly to global heat uptake, global heat transport, spatial temperature patterns, and variability. Quantifying ocean heat uptake across different temporal and spatial scales is important to quantify Earth's climate response to anthropogenic warming. Here we evaluate ocean adjustment time scales from two different fully coupled climate model ensembles using the Community Earth System Model. Both ensembles use the same model version, anthropogenic and natural forcings, and coupling configurations, but we initialize the ensembles in two different ways: (1) sampling joint internal variability of the ocean–atmosphere system (unique atmosphere and ocean conditions) and (2) sampling the internal variability of the atmosphere only (unique atmosphere, identical ocean conditions). Uncertainty due to internal variability is used as a proxy to quantify the time scales of ocean temperature adjustment at different depths and basins in Community Earth System Model. Time scales of equilibration are longer in the deep ocean than the upper ocean, highlighting the vertical structure of dynamic adjustment. The Atlantic equilibrates on shorter time scales (82 years above 1,000 m, 140 years below 1,000 m) relative to the Pacific (106 years above 1,000 m, 444 years below 1,000 m) in Community Earth System Model due to the large North Atlantic Deep Water formation and strong overturning circulation in the Atlantic. These results have broad implications for analyzing internal climate variability, ocean adjustment, and drift in global coupled model experiments and intercomparisons.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2003

Large Scale Water Vapor Sources Relative to the October 2000 Piedmont Flood

Very intense mesoscale or synoptic-scale rainfall events can occasionally be observed in the Mediterranean region without any deep cyclone developing over the areas affected by precipitation. In these perplexing cases the synoptic situation can superficially look similar to cases in which very little precipitation occurs. These situations could possibly baffle the operational weather forecasters. In this article, the major precipitation event that affected Piedmont (Italy) between 13 and 16 October 2000 is investigated. This is one of the cases in which no intense cyclone was observed within the Mediterranean region at any time, only a moderate system was present, and yet exceptional rainfall and flooding occurred. The emphasis of this study is on the moisture origin and transport. Moisture and energy balances are computed on different space- and time-scales, revealing that precipitation exceeds evaporation over an area inclusive of Piedmont and the northwestern Mediterranean region, on a time-scale encompassing the event and about two weeks preceding it. This is suggestive of an important moisture contribution originating from outside the region. A synoptic and dynamic analysis is then performed to outline the potential mechanisms that could have contributed to the large-scale moisture transport. The central part of the work uses a quasi-isentropic water-vapor back trajectory technique. The moisture sources obtained by this technique are compared with the results of the balances and with the synoptic situation, to unveil possible dynamic mechanisms and physical processes involved. It is found that moisture sources on a variety of atmospheric scales contribute to this event. First, an important contribution is caused by the extratropical remnants of former tropical storm Leslie. The large-scale environment related to this system allows a significant amount of moisture to be carried towards Europe. This happens on a time- scale of about 5-15 days preceding the Piedmont event. Second, water-vapor intrusions from the African Inter-Tropical Convergence Zone and evaporation from the eastern Atlantic contribute on the 2-5 day time-scale. The large-scale moist dynamics appears therefore to be one important factor enabling a moderate Mediterranean cyclone to produce heavy precipitation. Finally, local evaporation from the Mediterranean, water-vapor recycling, and orographically-induced low-level convergence enhance and concentrate the moisture over the area where heavy precipitation occurs. This happens on a 12-72 hour time-scale.

Turato, Barbara↗

NASA NTRS2010

Methods and apparatus for determining cardiac output

The present invention provides methods and apparatus for determining a dynamical property of the systemic or pulmonary arterial tree using long time scale information, i.e., information obtained from measurements over time scales greater than a single cardiac cycle. In one aspect, the invention provides a method and apparatus for monitoring cardiac output (CO) from a single blood pressure signal measurement obtained at any site in the systemic or pulmonary arterial tree or from any related measurement including, for example, fingertip photoplethysmography.According to the method the time constant of the arterial tree, defined to be the product of the total peripheral resistance (TPR) and the nearly constant arterial compliance, is determined by analyzing the long time scale variations (greater than a single cardiac cycle) in any of these blood pressure signals. Then, according to Ohm's law, a value proportional to CO may be determined from the ratio of the blood pressure signal to the estimated time constant. The proportional CO values derived from this method may be calibrated to absolute CO, if desired, with a single, absolute measure of CO (e.g., thermodilution). The present invention may be applied to invasive radial arterial blood pressure or pulmonary arterial blood pressure signals which are routinely measured in intensive care units and surgical suites or to noninvasively measured peripheral arterial blood pressure signals or related noninvasively measured signals in order to facilitate the clinical monitoring of CO as well as TPR.

Cohen, Richard J.↗

NASA NTRS1998

Boreal Winter Predictions with the GEOS-2 GCM: The Role of Boundary Forcing and Initial Conditions

Ensembles of atmospheric General Circulation Model (GCM) seasonal forecasts and long-term simulations (1980-94) are analyzed to assess the controlling influences of boundary forcing and memory of the initial conditions. Both the forecasts and simulations are carried out with version 2 of the Goddard Earth Observing System (GEOS-2) GCM forced with observed sea surface temperatures (SSTs). While much of the focus is on the seasonal time scale (January- March) and the Pacific North American (PNA) region, we also present results for other regions, shorter time scales, and other known modes of variability in the northern hemisphere extratropics. Forecasts of indices of some of the key large-scale modes of variability show that there is considerable variability in skill between different regions of the Northern Hemisphere. The eastern North Atlantic region has the poorest long lead forecast skill showing no skill beyond about 10 days. Skillful seasonal forecasts are primarily confined to the wave-like ENSO response emanating from the tropical Pacific. In the Northern Hemisphere, this is associated with the well-known Pacific/North American (PNA) pattern. Memory of the initial conditions is the major factor leading to skillful extratropical forecasts of lead time less than one month, while SST forcing is the only factor at the seasonal time scale. SST forcing contributes to skillful forecasts at sub- seasonal time scales only over the PNA region. The GEOS-2 GCM produces average (1980-94) signal to noise ratios which are less than one everywhere in the extratropics, except for the subtropical Pacific where they approach 1.5. When confined to the ENSO years, the maximum signal to noise ratios occur in the PNA region where they exceed three. An assessment of the sampling distribution of the forecasts suggests the model's ENSO response is very likely too weak. These results show some sensitivity to the uncertainties in the estimates of the SST forcing fields. In the North Pacific region, the sensitivity to SST forcing manifests itself primarily as changes in the variability of the PNA response, underscoring the need for an ensemble approach to the seasonal prediction problem.

Chang, Yehui↗

NASA NTRS1987

The formation of the exponential disk in spiral galaxies

A mechanism for the formation of the exponential stellar disk distribution in spiral galaxies is proposed. It is assumed that the gas which falls in to form the disk behaves like an accretion disk which evolves on a viscous time scale t(v) and forms stars on a time scale t(star). It is demonstrated that if t(v) roughly equals t(star) the resulting stellar disk has an exponential distribution independent of the disk rotation law and of the assumed viscosity prescription. Why the viscous and star formation time scales might be related is briefly commented on.

Lin, D. N. C.↗

NASA NTRS2010

The Causal Connection Between Disc and Power-Law Variability in Hard State Black Hole X-Ray Binaries

We use the XMM-Newton EPIC-pn instrument in timing mode to extend spectral time-lag studies of hard state black hole X-ray binaries into the soft X-ray band. \Ve show that variations of the disc blackbody emission substantially lead variations in the power-law emission, by tenths of a second on variability time-scales of seconds or longer. The large lags cannot be explained by Compton scattering but are consistent with time-delays due to viscous propagation of mass accretion fluctuations in the disc. However, on time-scales less than a second the disc lags the power-law variations by a few ms, consistent with the disc variations being dominated by X-ray heating by the power-law, with the short lag corresponding to the light-travel time between the power-law emitting region and the disc. Our results indicate that instabilities in the accretion disc are responsible for continuum variability on time-scales of seconds or longer and probably also on shorter time-scales.

Uttley, P.↗

DOE OSTI2022

Optical variability of quasars with 20-yr photometric light curves

ABSTRACT We study the optical gri photometric variability of a sample of 190 quasars within the SDSS Stripe 82 region that have long-term photometric coverage during ∼1998−2020 with SDSS, PanSTARRS-1, the Dark Energy Survey, and dedicated follow-up monitoring with Blanco 4m/DECam. With on average ∼200 nightly epochs per quasar per filter band, we improve the parameter constraints from a Damped Random Walk (DRW) model fit to the light curves over previous studies with 10–15 yr baselines and ≲ 100 epochs. We find that the average damping time-scale τDRW continues to rise with increased baseline, reaching a median value of ∼750 d (g band) in the rest frame of these quasars using the 20-yr light curves. Some quasars may have gradual, long-term trends in their light curves, suggesting that either the DRW fit requires very long baselines to converge, or that the underlying variability is more complex than a single DRW process for these quasars. Using a subset of quasars with better-constrained τDRW (less than 20 per cent of the baseline), we confirm a weak wavelength dependence of τDRW∝λ0.51 ± 0.20. We further quantify optical variability of these quasars over days to decades time-scales using structure function (SF) and power spectrum density (PSD) analyses. The SF and PSD measurements qualitatively confirm the measured (hundreds of days) damping time-scales from the DRW fits. However, the ensemble PSD is steeper than that of a DRW on time-scales less than ∼ a month for these luminous quasars, and this second break point correlates with the longer DRW damping time-scale.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Notch Fracture Toughness of Glasses: Dependence on Rate, Age, and Geometry

Understanding the fracture toughness (resistance) of glasses is a fundamental problem of prime theoretical and practical importance. Here we theoretically study its dependence on the loading rate, the age (history) of the glass, and the notch radius ρ. Reduced-dimensionality analysis suggests that the notch fracture toughness results from a competition between the initial, age- and history-dependent, plastic relaxation time scale τ$ ^{pl}_0 $ and an effective loading time scale τ ext ($ \dot{K} $ I , ρ ), where KI is the tensile stress-intensity-factor rate. The toughness is predicted to scale with ρ independently of ξ≡τ ext /τ$ ^{pl}_0 $ for ξ «1, to scale as T$ \sqrt{ρ} $log(ξ) for ξ«1 (related to thermal activation, where T is the temperature), and to feature a nonmonotonic behavior in the crossover region ξ~ O (1) (related to plastic yielding dynamics). These predictions are verified using 2D computations, providing a unified picture of the notch fracture toughness of glasses. Here, the theory highlights the importance of time-scale competition and far-from-steady-state elasto-viscoplastic dynamics for understanding the toughness and shows that the latter varies quite significantly with the glass age (history) and applied loading rate. Experimental support for bulk metallic glasses is presented, and possible implications for applications are discussed.

36 MATERIALS SCIENCE↗

NASA NTRS2008

Relationship between Mean Winds and Large 10-day Precipitation Anomalies over the Continental U.S.

To better understand the role of surface boundary forcing upon interannual and decadal variability of droughts and precipitation excesses over the continental U. S., an observational study has been carried out to better investigate the dynamical linkages between surface boundary forcing and droughts/floods. This study examines the relationship during the warm season between observed mean wind patterns on the 10-15 day time scale and anomalies in precipitation on that same time scale over the U. S. Each of the 7 cases of significant 10-day, warm-season drought over the U, S. since 1948 that has been investigated was accompanied by significant continental-scale ridging of the mean 10-day 250 mb wind pattern over the continental U.S. and a mid-Pacific relative maximum of 250 mb wind speed of at least 30 m/sec within 20 degrees of the international Date Line. This perhaps suggests the presence of a wave-guide phenomenon over the Pacific. The drought signal for each of these cases was still quite evident when averaged over the monthly time scale. While peak magnitudes of the precipitation anomalies were much larger for the 7 chosen cases of significant of precipitation excess were than for the drought cases none of the patterns of precipitation excess were as unambiguous or wide-spread as the drought cases. In only 2 to 3 cases was there evidence of a continental trough in the 250 mb winds. In fact, one case showed slight ridging over Canada and the extreme northern U.S. Strong southerly low-level jets, however, were evident in the 925 mb winds over south-central Texas or the nearby Gulf of Mexico in virtually all of the excessive precipitation cases, while the Great Plains LLJ was weak or displaced to the north of Texas in most of the drought cases. We will further investigate the asymmetry between droughts and precipitation excesses in their horizontal homogeneity and extent, in their temporal evolution, and in how the are forced by the combined effects of the mean 250 mb westerlies and the mean 925 mb southerly winds.

Helfand, H. Mark↗

NASA NTRS2020

Radiography in High Mass X-Ray Binaries: Micro-Structure of the Stellar Wind Through Variability of the Column Density

Context. In high mass X-ray binaries, an accreting compact object orbits a high mass star, which loses mass through a dense and in homogeneous wind. Aims. Using the compact object as an X-ray backlight, the time variability of the absorbing column density in the wind can be exploited in order to shed light on the micro-structure of the wind and obtain unbiased stellar mass-loss rates for high mass stars. Methods. We developed a simplified representation of the stellar wind where all the matter is gathered in spherical “clumps” that are radially advected away from the star. This model enables us to explore the connections between the stochastic properties of the wind and the variability of the column density for a comprehensive set of parameters related to the orbit and to the wind micro-structure, such as the size of the clumps and their individual mass. In particular, we focus on the evolution with the orbital phase of the standard deviation of the column density and of the characteristic duration of enhanced absorption episodes. Using the porosity length, we derive analytical predictions and compare them to the standard deviations and coherence time scales that were obtained. Results. We identified the favorable systems and orbital phases to determine the wind micro-structure. The coherence time scale of the column density is shown to be the self-crossing time of a single clump in front of the compact object. We thus provide a procedure to get accurate measurements of the size and of the mass of the clumps, purely based on the observable time variability of the column density. Conclusions. The coherence time scale grants direct access to the size of the clumps, while their mass can be deduced separately from the amplitude of the variability. We further show how monitoring the variability at superior conjunctions can probe the onset of the clump-forming region above the stellar photosphere. If the high column density variations in some high mass X-ray binaries are due to unaccreted clumps which are passing by the line-of-sight, this would require high mass clumps to reproduce the observed peak-to-peak amplitude and coherence time scales. These clump properties are marginally compatible with the ones derived from radiative-hydrodynamics simulations. Alternatively, the following components could contribute to the variability of the column density: larger orbital scale structures produced by a mechanism that has yet to be identified or a dense environment in the immediate vicinity of the accretor, such as an accretion disk, an outflow, or a spherical shell surrounding the magnetosphere of the accreting neutron star.

stars: mass loss↗