We present a reservoir simulation framework for coupled thermal-compositional-mechanics processes. We use finite-volume methods to discretize the mass and energy conservation equations and finite-element methods for the mechanics problem. We use the first-order backward Euler for time. We solve the resulting set of nonlinear algebraic equations using fully implicit (FI) and sequential-implicit (SI) solution schemes. The FI approach is attractive for general-purpose simulation due to its unconditional stability. However, the FI method requires the development of a complex thermo-compositional-mechanics framework for the nonlinear problems of interest, and that includes the construction of the full Jacobian matrix for the coupled multi-physics discrete system of equations. On the other hand, SI-based solution schemes allow for relatively fast development because different simulation modules can be coupled more easily. The challenge with SI schemes is that the nonlinear convergence rate depends strongly on the coupling strength across the physical mechanisms and on the details of the sequential updating strategy across the different physics modules. The flexible automatic differentiation-based framework described here allows for detailed assessment of the robustness and computational efficiency of different coupling schemes for a wide range of multi-physics subsurface problems.

The first-known tidal disruption event (TDE) with strong evidence for a relativistic jet - based on extensive multiwavelength campaigns - is Swift J1644+5734. In order to directly measure the apparent speed of the radio jet, we performed very long baseline interferometry (VLBI) observations with the European VLBI network (EVN) at 5 GHz. Our observing strategy was to identify a very nearby and compact radio source with the real-time e-EVN, and then utilize this source as a stationary astrometry reference point in the later five deep EVN observations.With respect to the in-beam source FIRST J1644+5736, we have achieved a statistical astrometric precision about 12 μas (68 per cent confidence level) per epoch. This is one of the best phase-referencing measurements available to date. No proper motion has been detected in the Swift J1644+5734 radio ejecta. We conclude that the apparent average ejection speed between 2012.2 and 2015.2 was less than 0.3c with a confidence level of 99 per cent. This tight limit is direct observational evidence for either a very small viewing angle or a strong jet deceleration due to interactions with a dense circum-nuclear medium, in agreement with some recent theoretical studies.

Few high-redshift, radio-loud quasars are known to date. The extremely luminous, radio-bright quasar, SDSS J013127.34-032100.1 was recently discovered at a redshift of z = 5.18. We observed the source with high resolution very long baseline interferometry (VLBI) at 1.7 GHz with the European VLBI Network (EVN) and found a single compact radio component. We estimated a lower limit to the brightness temperature of the detected radio component, T_B ~10¹¹ K. Additionaly, when compared to archival radio data, the source showed significant flux density variation. These two findings are indicative of the blazar nature of the source.

Steam Assistant Gravity Drainage (SAGD) is widely used to recover heavy oil and bitumen reservoirs. Typical SAGD operations involve a pair of horizontal wells separated vertically. Steam, or a steam-solvent mixture (e.g., Expanding- Solvent SAGD), is injected into the upper well to form a steam chamber and mobilize oil, which drains to the lower production well. Significant mechanical stresses associated with this process can increase the risk of fracturing the reservoir, or cap-rock. We perform a fully-coupled thermal-compositional-mechanical numerical simulation of SAGD and ES-SAGD processes for a typical bitumen reservoir in the Fort McMurray region of Alberta, Canada. A mixed finite-volume approximation for the flow and a Galerkin finite-element approximation for the mechanics are used, and the resulting set of nonlinear equations are solved using a fully implicit formulation. The two discretizations share the same unstructured grid. We demonstrate that thermo-mechanical effects can be quite significant in SAGD operations. The sharper gradients associated with the standard SAGD process increase the risk of damaging of the cap-rock. On the other hand, ES-SAGD operations lead to more dispersed temperature and pressure distributions, which decreases the possibility of damaging the cap-rock.