Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLARE”

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 145 records · Page 8

A Coupled Analysis of Atmospheric Mass Loss and Tidal Evolution in XUV Irradiated Exoplanets: The TRAPPIST-1 Case Study

Exoplanets residing close to their stars can experience evolution of both their physical structures and their orbits due to the influence of their host stars. In this work, we present a coupled analysis of dynamical tidal dissipation and atmospheric mass loss for exoplanets in X-ray and ultraviolet (XUV) irradiated environments. As our primary application, we use this model to study the TRAPPIST-1 system and place constraints on the interior structure and orbital evolution of the planets. We start by reporting on an ultraviolet continuum flux measurement (centered around ∼1900 Å) for the star TRAPPIST-1, based on 300 ks of Neil Gehrels Swift Observatory data, and which enables an estimate of the XUV-driven thermal escape arising from XUV photodissociation for each planet. We find that the X-ray flaring luminosity, measured from our X-ray detections, of TRAPPIST-1 is 5.6 × 10{sup −4} L {sub *}, while the full flux including non-flaring periods is 6.1 × 10{sup −5} L {sub *}, when L {sub *} is TRAPPIST-1's bolometric luminosity. We then construct a model that includes both atmospheric mass loss and tidal evolution and requires the planets to attain their present-day orbital elements during this coupled evolution. We use this model to constrain the ratio Q{sup ′}=3Q/2k{sub 2} for each planet. Finally, we use additional numerical models implemented with the Virtual Planet Simulator VPLanet to study ocean retention for these planets using our derived system parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

The Anatomy of an Unusual Edge-on Protoplanetary Disk. I. Dust Settling in a Cold Disk

As the earliest stage of planet formation, massive, optically thick, and gas-rich protoplanetary disks provide key insights into the physics of star and planet formation. When viewed edge-on, high-resolution images offer a unique opportunity to study both the radial and vertical structures of these disks and relate this to vertical settling, radial drift, grain growth, and changes in the midplane temperatures. In this work, we present multi-epoch Hubble Space Telescope and Keck scattered light images, and an Atacama Large Millimeter/submillimeter Array 1.3 mm continuum map for the remarkably flat edge-on protoplanetary disk SSTC2DJ163131.2–242627, a young solar-type star in ρ Ophiuchus. We model the 0.8 μm and 1.3 mm images in separate Markov Chain Monte Carlo (MCMC) runs to investigate the geometry and dust properties of the disk using the MCFOST radiative transfer code. In scattered light, we are sensitive to the smaller dust grains in the surface layers of the disk, while the submillimeter dust continuum observations probe larger grains closer to the disk midplane. An MCMC run combining both data sets using a covariance-based log-likelihood estimation was marginally successful, implying insufficient complexity in our disk model. The disk is well characterized by a flared disk model with an exponentially tapered outer edge viewed nearly edge-on, though some degree of dust settling is required to reproduce the vertically thin profile and lack of apparent flaring. A colder than expected disk midplane, evidence for dust settling, and residual radial substructures all point to a more complex radial density profile to be probed with future, higher-resolution observations.

47 OTHER INSTRUMENTATION↗

The High-energy Spectrum of the Nearby Planet-hosting Inactive Mid-M Dwarf LHS 3844

To fully characterize the atmospheres, or lack thereof, of terrestrial exoplanets, we must include the high-energy environments provided by their host stars. The nearby mid-M dwarf LHS 3844 hosts a terrestrial world that lacks a substantial atmosphere. We present a time-series UV spectrum of LHS 3844 from 1131 to 3215 Å captured by HST/COS. We detect one flare in the FUV that has an absolute energy of 8.96 ± 0.77 × 10{sup 28} erg and an equivalent duration of 355 ± 31 s. We extract the flare and quiescent UV spectra separately. For each spectrum, we estimate the Lyα flux using correlations between UV line strengths. We use Swift-XRT to place an upper limit on the soft X-ray flux and construct a differential emission model to estimate flux that is obscured by the interstellar medium. We compare the differential emission model flux estimates in the XUV to other methods that rely on scaling from the Lyα, Si iv, and N v lines in the UV. The XUV, FUV, and NUV flux of LHS 3844 relative to its bolometric luminosity is log{sub 10}(L {sub band}/L {sub Bol}) =−3.65, −4.16, and −4.48, respectively, for the quiescent state. These values agree with trends in high-energy flux as a function of stellar effective temperature found by the MUSCLES survey for a sample of early-M dwarfs. Many of the most spectroscopically accessible terrestrial exoplanets orbit inactive mid-to-late M dwarfs like LHS 3844. Measurements of M dwarf high-energy spectra are preferable for exoplanet characterization but are not always possible. The spectrum of LHS 3844 is a useful proxy for the current radiation environment for these worlds.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Seeking Echoes of Circumstellar Disks in Kepler Light Curves

Light echoes of flares on active stars offer the opportunity for direct detection of circumstellar dust. We revisit the problem of identifying faint echoes in postflare light curves, focusing on debris disks from ongoing planet formation. Starting with simulations, we develop an algorithm for estimating the radial extent and total mass from disk echo profiles. We apply this algorithm to light curves from over 2100 stars observed by NASA’s Kepler mission, selected for multiple short-lived flares in either the long-cadence or short-cadence data sets. While flux uncertainties in light curves from individual stars preclude useful mass limits on circumstellar disks, catalog-averaged light curves yield constraints on disk mass that are comparable to estimates from known debris disks. The average mass in micron- to millimeter-sized dust around the Kepler stars cannot exceed 10% of an Earth mass in exo-Kuiper belts or 10% of a lunar mass in the terrestrial zone. We group stars according to IR excess, based on WISE W1–W3 color, as an indicator for the presence of circumstellar dust. The mass limits are greater for stars with strong IR excess, a hint that echoes are lurking not far beneath the noise in postflare light curves. With increased sensitivity, echo detection will let time-domain astronomy complement spectroscopic and direct-imaging studies in mapping how, when, and where planets form.

79 ASTRONOMY AND ASTROPHYSICS↗

Turbulent Model of Crab Nebula Radiation

We construct a turbulent model of the Crab Nebula's nonthermal emission. The present model resolves a number of long-standing problems of the Kennel–Coroniti model: (i) the sigma problem, (ii) the hard spectrum of radio electrons, (iii) the high peak energy of gamma-ray flares, (iv) and the spatial evolution of the infrared (IR) emission. The Nebula contains two populations of injected particles: Component-I, accelerated at the wind termination shock via the Fermi-I mechanism; and Component-II, accelerated in reconnecting turbulence in highly magnetized (σ >> 1) plasma in the central part of the Crab Nebula. The reconnecting turbulence in Component-II extends from radio to gamma-rays: it accelerates radio electrons with a hard spectrum, destroys the large-scale magnetic flux (and thus resolves the sigma problem), and occasionally produces gamma-ray flares (from the largest-scale reconnection events). The model reproduces the broadband spectrum of the Crab Nebula, from low-frequency synchrotron emission in radio to inverse-Compton emission at TeV energies, as well as the spatially resolved evolution of the spectral indices in the IR and optical bands.

79 ASTRONOMY AND ASTROPHYSICS↗

Radiation and Polarization Signatures from Magnetic Reconnection in Relativistic Jets. I. A Systematic Study

Blazars are relativistic magnetized plasma outflows from supermassive black holes that point very close to our line of sight. Their emission is nonthermal-dominated and highly variable across the entire electromagnetic spectrum. Relativistic magnetic reconnection has been proposed as the driver of particle acceleration during blazar flares. While recent particle-in-cell (PIC) simulations have self-consistently studied the evolution of magnetic reconnection and particle acceleration therein, the resulting radiation signatures have not been systematically explored. In particular, the polarization signatures, which directly reflect the characteristic strongly dynamical evolution of magnetic field during reconnection, have not been carefully investigated. In this paper, we present a systematic study of radiation and polarization signatures arising from magnetic reconnection in blazars, based on combined PIC and polarized radiation transfer simulations with various physical parameters. We identify a harder-when-brighter trend in the spectral evolution. Moreover, higher-frequency bands (ultraviolet to X-ray) tend to flare earlier than lower-frequency bands (infrared to optical) in the synchrotron spectral component. Most importantly, polarization signatures appear more variable with higher frequencies. We find that the variation in temporal polarization depends strongly on the guide field strength. Specifically, reconnection with a significant guide field component leads to a very high polarization degree that contradicts typical blazar observations, while large polarization angle rotations are unique signatures of magnetic reconnection between nearly antiparallel magnetic field lines. These rotations are at least 90° and can extend to >180°, and they may be in either direction. Furthermore, these results imply that blazars that have shown large polarization angle rotations intrinsically have more nearly antiparallel magnetic field morphology.

79 ASTRONOMY AND ASTROPHYSICS↗

ASASSN-14ko is a Periodic Nuclear Transient in ESO 253-G003

We present the discovery that ASASSN-14ko is a periodically flaring active galactic nucleus at the center of the galaxy ESO 253-G003. At the time of its discovery by the All-Sky Automated Survey for Supernovae (ASAS-SN), it was classified as a supernova close to the nucleus. The subsequent 6 yr of V - and g -band ASAS-SN observations revealed that ASASSN-14ko has nuclear flares occurring at regular intervals. The 17 observed outbursts show evidence of a decreasing period over time, with a mean period of P 0 = 114.2 ± 0.4 days and a period derivative of . The most recent outburst in 2020 May, which took place as predicted, exhibited spectroscopic changes during the rise and had a UV bright, blackbody spectral energy distribution similar to tidal disruption events (TDEs). The X-ray flux decreased by a factor of 4 at the beginning of the outburst and then returned to its quiescent flux after ~8 days. The Transiting Exoplanet Survey Satellite observed an outburst during Sectors 4–6, revealing a rise time of 5.60 ± 0.05 days in the optical and a decline that is best fit with an exponential model. We discuss several possible scenarios to explain ASASSN-14ko’s periodic outbursts, but currently favor a repeated partial TDE. The next outbursts should peak in the optical on UT 2020 September 7.4±1.1 and UT 2020 December 26.5±1.4.

Astronomy & Astrophysics↗

A Possible Tidal Disruption Event Candidate in the Black Hole Binary System of OJ 287

The BL Lacertae OJ 287 is a supermassive black hole binary system with complex physics of its irregular flares. During 2016 October–2017 April, a surprising outburst in the X-ray, UV, and optical bands was detected, while no variability was seen in the γ-ray light curve. During the outburst, the X-ray light curves were dominated by the soft X-rays, whose peak in luminosity was ~10 46 erg s -1 —more than 10 times higher than the mean level before the outburst—and a "softer-when-brighter" phenomenon was exhibited. These above phenomena have been reported by some previous works. The hardness ratio showed negligible evolution with flare time and soft X-ray luminosity. Critically, the luminosity of the soft X-rays decayed following a power law of $t$ -5/3 , which occurs in most tidal disruption events (TDEs), and a similar trend can be seen in the UV and optical bands during the soft X-ray declining period. Helium and oxygen narrow emission lines are strengthened prominently in the optical spectra of post outburst epochs, that could be attributed to the surrounding gas caused by TDE. In conclusion, we discuss three possible origins of the event, including the jet's precession, the aftereffects of the black hole–disk impaction, and the TDE. Our results show that the TDE is the more likely scenario to explain the outburst.

79 ASTRONOMY AND ASTROPHYSICS↗

Slow Shock Formation Upstream of Reconnecting Current Sheets

The formation, development, and impact of slow shocks in the upstream regions of reconnecting current layers are explored. Slow shocks have been documented in the upstream regions of magnetohydrodynamic (MHD) simulations of magnetic reconnection as well as in similar simulations with the kglobal kinetic macroscale simulation model. They are therefore a candidate mechanism for preheating the plasma that is injected into the current layers that facilitate magnetic energy release in solar flares. Of particular interest is their potential role in producing the hot thermal component of electrons in flares. During multi-island reconnection, the formation and merging of flux ropes in the reconnecting current layer drives plasma flows and pressure disturbances in the upstream region. These pressure disturbances steepen into slow shocks that propagate along the reconnecting component of the magnetic field and satisfy the expected Rankine–Hugoniot jump conditions. Plasma heating arises from both compression across the shock and the parallel electric field that develops to maintain charge neutrality in a kinetic system. Shocks are weaker at lower plasma β, where shock steepening is slow. While these upstream slow shocks are intrinsic to the dynamics of multi-island reconnection, their contribution to electron heating remains relatively minor compared with that from Fermi reflection and the parallel electric fields that bound the reconnection outflow.

79 ASTRONOMY AND ASTROPHYSICS↗

Transient Radio Emission from Low-redshift Galaxies at z < 0.3 Revealed by the VLASS and FIRST Surveys

We present the discovery of a sample of 18 low-redshift (z < 0.3) galaxies with transient nuclear radio emission. These galaxies are not detected or are weakly detected in the Faint Images of the Radio Sky at Twenty cm survey, performed from 1993–2009, but have brightened significantly in radio flux (by a factor of ≳5) in the epoch I (2017–2019) observations of the Very Large Array Sky Survey (VLASS). All 18 galaxies have been detected in VLASS epoch II observations, from 2020–2021, from which the radio flux has been found to evolve slowly (with variability amplitudes of ≳40%) over a period of about 3 yr. 15 galaxies have been observed in the Rapid ASKAP Continuum Survey, and a flat or inverted spectral slope between 888 MHz and 3 GHz is found. Based on the Sloan Digital Sky Survey spectra taken before the radio brightening, 14 of the 18 galaxies can be classified as LINERs or normal galaxies with weak or no nuclear activity. Most galaxies are red and massive, with more than half having central black hole masses above 10 8 M ⊙ . We find that only one galaxy in our sample displays an optical flare lasting for at least two months, with a long decay in the infrared light curve that can be explained as the dust-heated echo emission of a central optical flare, such as a stellar tidal disruption event. We discuss several possibilities for the transient radio emission and conclude that it is likely associated with a newborn radio jet triggered by short sporadic fueling of a supermassive black hole. Such a scenario can be tested with further multifrequency radio observations of these sources, via measuring their radio flux variability and spectral evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Tidal Stripping of a White Dwarf by an Intermediate-mass Black Hole

During the inspiralling of a white dwarf (WD) into an intermediate-mass black hole (~10 2-5 M ⊙ ), both gravitational waves (GWs) and electromagnetic (EM) radiation are emitted. Once the eccentric orbit's pericenter radius approaches the tidal radius, the WD would be tidally stripped upon each pericenter passage. The accretion of this stripped mass would produce EM radiation. It is suspected that the recently discovered new types of transients, namely the quasiperiodic eruptions and the fast ultraluminous x-ray bursts, might originate from such systems. Modeling these flares requires a prediction of the amount of stripped mass from the WD and the details of the mass supply to the accretion disk. We run hydrodynamical simulations to study the orbital parameter dependence of the stripped mass. We find that our results match the analytical estimate that the stripped mass is proportional to z 5/2 , where z is the excess depth by which the WD overfills its instantaneous Roche lobe at the pericenter. The corresponding fallback rate of the stripped mass is calculated, which may be useful in interpreting the individual flaring light curve in candidate EM sources. We further calculate the long-term mass-loss evolution of a WD during its inspiral and the detectability of the GW and EM signals. The EM signal from the mass-loss stage can be easily detected: the limiting distance is ~320(M h /10 4 M ⊙ ) Mpc for the Einstein Probe. The GW signal, for space-borne detectors such as Laser Interferometer Space Antenna or TianQin, can be detected only within the Local Supercluster (~33 Mpc).

79 ASTRONOMY AND ASTROPHYSICS↗

Simulations of the Collision between a Debris Stream and an Outer Dusty Torus: A Possible Channel for Forming a Fast-rise and Long-delay Radio Outburst in Tidal Disruption Events

A geometrically thick dusty torus structure is believed to exist in the nuclear regions of galaxies (especially in active galactic nuclei). The debris stream from a tidal disruption event (TDE) will possibly collide with the dusty torus and produce a transient flare. We perform three-dimensional hydrodynamic simulations to model the dynamical evolution of the interaction between unbound debris and a dusty torus. During the continuous interaction, shocked material will be spilled out from the interaction region and form an outflow. We calculate the temporal evolution of synchrotron emission by assuming that the shock accelerates a fraction of electrons in the outflow into a nonthermal distribution. We find that radio emission from the debris–torus collision generates a steep-rise and slow-decline radio light curve due to the sharp edge and dense gas of the dusty torus, where the radio outburst is delayed relative to the main optical/X-ray outburst by several years or even several tens of years. We apply our model to a TDE that happened in a narrow-line Seyfert I galaxy (PS16dtm), where both the radio spectrum and the light curve can be roughly reproduced. Future high-sensitivity, wide-field-of-view radio surveys have the opportunity to detect more such radio flares.

79 ASTRONOMY AND ASTROPHYSICS↗

Ion-rich Acceleration during an Eruptive Flux Rope Event in a Multiple Null-point Configuration

We report on the γ-ray emission above 100 MeV from the GOES M3.3 flare SOL2012-06-03. The hard X-ray (HXR) and microwave emissions have typical time profiles with a fast rise to a well-defined peak followed by a slower decay. The >100 MeV emission during the prompt phase displayed a double-peaked temporal structure with the first peak following the HXR and microwaves, and the second one, about 3 times stronger, occurring 17 ± 2 s later. The time profiles seem to indicate two separate acceleration mechanisms at work, where the second γ-ray peak reveals a potentially pure or at least largely dominant ion acceleration. The Atmospheric Imaging Assembly imaging shows a bright elliptical ribbon and a transient brightening in the northwestern (NW) region. Nonlinear force-free extrapolations at the time of the impulsive peaks show closed field lines connecting the NW region to the southeastern part of the ribbon, and the magnetic topology revealed clusters of nulls. These observations suggest a spine-and-fan geometry, and based on these observations, we interpret the second γ-ray peak as being due to the predominant acceleration of ions in a region with multiple null points. The >100 MeV emission from this flare also exhibits a delayed phase with an exponential decay of roughly 350 s. We find that the delayed emission is consistent with ions being trapped in a closed flux tube with gradual escape via their loss cone to the chromosphere.

Pesce-Rollins, Melissa [Istituto Nazionale di Fisi↗

All-sky Neutrino Point-source Search with IceCube Combined Track and Cascade Data

Despite extensive efforts, discovery of high-energy astrophysical neutrino sources remains elusive. We present an event-level simultaneous maximum likelihood analysis of tracks and cascades using IceCube data collected from 2008 April 6 to 2022 May 23 to search the whole sky for neutrino sources, and using a source catalog, for coincidence of neutrino emission with gamma-ray emission. This is the first time a simultaneous fit of different detection channels is used to conduct a time-integrated all-sky scan with IceCube. Combining all-sky tracks, with superior pointing power and sensitivity in the northern sky, with all-sky cascades, with good energy resolution and sensitivity in the southern sky, we have developed the most sensitive point-source search to date by IceCube that targets the entire sky. The most significant point in the northern sky aligns with NGC 1068, a Seyfert II galaxy, which, from the catalog search, shows a 3.5σ excess over background after accounting for trials. The most significant point in the southern sky does not align with any source in the catalog and is not significant after accounting for trials. A search for the single most significant Gaussian flare at the locations of NGC 1068, PKS 1424+240, and the southern highest-significance point shows results consistent with expectations for steady emission. Notably, this is the first time that a flare shorter than four years has been excluded as being responsible for NGC 1068’s emergence as a neutrino source. Our results show that combining tracks and cascades when conducting neutrino source searches improves sensitivity and can lead to new discoveries.

Abbasi, R. [Loyola University, Chicago, IL (United↗

Millimeter Light Curves of Sagittarius A* Observed during the 2017 Event Horizon Telescope Campaign

The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5–11 in the 1.3 mm wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data sets, complementing the EHT very long baseline interferometry, are characterized by a cadence and signal-to-noise ratio previously unattainable for Sgr A* at millimeter wavelengths, and they allow for the investigation of source variability on timescales as short as a minute. While most of the light curves correspond to a low variability state of Sgr A*, the April 11 observations follow an X-ray flare and exhibit strongly enhanced variability. All of the light curves are consistent with a red-noise process, with a power spectral density (PSD) slope measured to be between −2 and −3 on timescales between 1 minute and several hours. Our results indicate a steepening of the PSD slope for timescales shorter than 0.3 hr. The spectral energy distribution is flat at 220 GHz, and there are no time lags between the 213 and 229 GHz frequency bands, suggesting low optical depth for the event horizon scale source. We characterize Sgr A*’s variability, highlighting the different behavior observed just after the X-ray flare, and use Gaussian process modeling to extract a decorrelation timescale and a PSD slope. We also investigate the systematic calibration uncertainties by analyzing data from independent data reduction pipelines.

79 ASTRONOMY AND ASTROPHYSICS↗

Findings on Subtask 3.1 - Bakken Rich Gas Enhanced Oil Recovery Project

Total in-place oil for the Bakken petroleum system (BPS) (which includes the Bakken and Three Forks Formations) has been estimated to be 600 billion barrels (bbl). However, BPS wells have decline rates as high as 85% over the first 3 years of their lives, and primary recovery factors typically range from 3% to 10% of original oil in place. Given the low initial recovery rates, even small incremental productivity improvements could dramatically increase technically recoverable oil in the BPS. One potential solution is enhanced oil recovery (EOR) using gas injection, such as carbon dioxide (CO2) or hydrocarbon (HC) gases. While commonly used in conventional reservoirs, CO2 EOR in unconventional tight oil reservoirs has been limited to pilot tests. EOR using rich gas (mixture of methane, ethane, and propane) has also been employed in numerous pilots in several unconventional plays and has recently been successfully applied in the Eagle Ford play. If successful, large-scale gas-based EOR in the BPS could dramatically increase oil productivity and recovery factors and extend the life of the play for decades. While CO2 may be a technically suitable working fluid for EOR in the BPS, supplies are limited and costs for using CO2 in EOR pilots are prohibitively high. Meanwhile, produced gas flaring has presented challenges for BPS operators in North Dakota. Analysis conducted by the North Dakota Pipeline Authority indicates that the current gas-gathering infrastructure in North Dakota is insufficient to accommodate all of the associated gas that is produced from the BPS. The geographically isolated location of North Dakota relative to large natural gas markets, combined with suppressed natural gas prices, has made it economically challenging for industry to invest capital in expanding gas-gathering infrastructure in the state. These circumstances led to a research program conducted by the Energy & Environmental Research Center (EERC) in partnership with Liberty Resources Management Company LLC (LR) to examine the potential to use rich gas injection for EOR and mitigate flaring. A rich gas EOR pilot test was designed and executed by LR at its Stomping Horse development area in Williams County, North Dakota. From July 2018 through May 2019, a total of 160 million standard cubic feet (MMscf) of rich produced gas was injected into the BPS using five different wells in a sequential injection strategy. LR’s Leon–Gohrick drill spacing unit (DSU) was used as the test site. Regulatory oversight was provided by the North Dakota Industrial Commission (NDIC). Technical support was provided by the EERC through a series of laboratory, modeling, and field-based activities, and additional post-pilot research activities incorporated learnings from the test, developed new laboratory data, improved fracture modeling methods, and developed machine learning and big data analytics. The results from the Stomping Horse rich gas EOR pilot activities indicate that developing an effective, economical EOR approach for the BPS will require more field tests. Another key lesson learned from the Stomping Horse tests is that detailed pre- and posttest data on reservoir conditions and fluids production are essential. Robust reservoir characterization provides information that is crucial to creating realistic geomodels and conducting valid dynamic simulations of potential EOR scenarios. A detailed understanding of the completions and production history of offset wells is also necessary for valid test result interpretations. This knowledge is essential to designing the operational parameters of injectivity tests and interpreting the results. A conformance control strategy is also essential to success. Laboratory-based examinations of rich gas interactions with reservoir fluids and rocks were conducted, with an emphasis on determining the ability to mobilize oil in the tight reservoir rocks and shales of the BPS. Injection fluid composition was shown to have a positive impact on reducing reservoir oil minimum miscibility pressure (MMP), reducing interfacial tension (IFT), and altering wettability. IFT and contact angle measurements demonstrated that wettability can be altered in the presence of rich gas, suggesting the potential to improve oil recovery. Iterative modeling of surface infrastructure and reservoir performance using data generated by the various project activities was conducted. A geologic model of the Stomping Horse area was built; history-matched oil, gas, and water production was used in simulations of various EOR scenarios. Early programmatic modeling results were used to support LR’s design and operation of the EOR pilot and to provide insight regarding optimization of future commercial-scale BPS EOR design and operations. Post-pilot modeling focused on alternative methods of understanding complex fracture networks and accelerating simulation time. These led to improved simulation run times and provide excellent history-matching results. Several of these iterative models were used as the bases for developing algorithms into machine learning and big data analytics. History matching in reservoir simulation is time-consuming and computer processing-intensive. Machine learning algorithms were created, and an automated history-matching tool was developed. A large set of synthetic reservoir simulations were created to generate well responses (oil, gas, and water production, well bottomhole pressure [BHP], and tracer or propane breakthrough) for a set of EOR operating parameters that included offset well status (open or closed), injectate (rich gas or propane), injection rate, and injection well BHP. A user interface was developed to provide real-time visualization. Machine learning-based models were developed to provide rapid forecasting of well performance given a set of user-defined EOR operating parameters. These predictive models allow the user to modify the offset well status, injection rate, and injection well BHP and rapidly forecast future production performance. The combination of real-time visualization tools with real-time forecasting tools provides a framework for real-time control—operational changes that the EOR site operator can enact (e.g., changing gas injection rates) to affect the observed performance and potentially improve the EOR outcome. There is great reason to be optimistic about the future of EOR in the Bakken. The results of the laboratory studies suggest significant potential for high rates of oil mobilization using produced field gas injection under the right conditions. The results of the lab studies, combined with rigorous statistical analysis of well production data and associated modeling efforts, confirm the notion that fluid mobility within the reservoir is controlled by fractures. As more knowledge is gained about the nature and distribution of fracture networks in the Bakken, the industry will be in a better position to predict and, ultimately, influence fluid mobility. New field tests are necessary to develop a more complete understanding of those conditions. Thoughtful and creatively engineered field tests within a well-characterized geologic setting will yield the fundamental knowledge needed to take Bakken oil production to the next level. This subtask was cofunded through the EERC–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission’s Oil and Gas Research Program and Computer Modelling Group.

04 OIL SHALES AND TAR SANDS↗

Development of PFLOTRAN Transport Capability for Use in the Waste Isolation Pilot Plant Performance Assessment - 20545

Waste Isolation Pilot Plant (WIPP) performance assessment (PA) calculations estimate the probability of radionuclide release from the repository to the land surface and across the land withdrawal boundary for a regulatory period of 10,000 years after facility closure. Simulations of flow and transport in the repository and the surrounding Salado Formation are foundational to the PA. Because proposed additional waste emplacement panels would result in an asymmetric repository layout, the US Department of Energy (DOE) is preparing to transition to use of a three-dimensional (3-D) model domain for simulation of flow and transport instead of the two-dimensional (2-D) flared grid domain currently used. DOE has charged Sandia National Laboratories with developing the capability necessary to simulate processes affecting flow and transport in the WIPP in PFLOTRAN, an open-source massively parallel multi-phase flow and reactive transport code. The new flow and transport capabilities developed in PFLOTRAN incorporate WIPP-specific process models and will replace the 2-D simulators (BRAGFLO and NUTS) that are currently utilized for Salado flow and transport calculations in WIPP PA. The focus of this paper is on the development of a new Nuclear Waste Transport (NWT) mode in PFLOTRAN that has all of the capabilities necessary for Salado transport simulations, including the ability to handle complete dry-out (100% gas saturation) of arbitrary cells in the model domain, radionuclide mass conservation at step changes in porosity associated with borehole intrusion, and the ability to calculate fluxes on a flared grid. The new PFLOTRAN transport capability and a suite of verification tests were designed around a list of functional requirements for WIPP PA calculations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

CubeSats for Gamma-Ray Astronomy [Book Chapter]

After many years of flying in space primarily for educational purposes, CubeSats – tiny satellites with form factors corresponding to arrangements of “1U” units, or cubes, each 10 cm on a side – have come into their own as valuable platforms for technology advancement and scientific investigations. CubeSats offer comparatively rapid, low-cost access to space for payloads that can be built, tested, and operated by relatively small teams, with substantial contributions from students and early career researchers. Continuing advances in compact, low-power detectors, readout electronics, and flight computers have now enabled X-ray and gamma-ray sensing payloads that can fit within the constraints of CubeSat missions, permitting in-orbit demonstrations of new techniques and innovative high-energy astronomy observations. Gamma-ray-sensing CubeSats are certain to make an important contribution in the new era of multi-messenger, time-domain astronomy by detecting and localizing bright transients such as gamma-ray bursts, solar flares, and terrestrial gamma-ray flashes; however, other astrophysical science areas requiring long observations in a low-background environment, including gamma-ray polarimetry, studies of nuclear lines, and measurement of diffuse backgrounds, will likely benefit as well. Here we present the primary benefits of CubeSats for high-energy astronomy, highlight the scientific areas currently or soon to be studied, and review the missions that are currently operating, under development, or proposed. A rich portfolio of CubeSats for gamma-ray astronomy already exists, and the potential for a broad range of creative and scientifically productive missions in the near future is very high.

79 ASTRONOMY AND ASTROPHYSICS↗