Double white dwarf (WD) binaries are increasingly being discovered at short orbital periods where strong tidal effects and significant tidal heating signatures may occur. We assume that the tidal potential of the companion excites outgoing gravity waves within the WD primary, the dissipation of which leads to an increase in the WD’s surface temperature. We compute the excitation and dissipation of the waves in cooling WD models in evolving
We describe new functionality in the GYRE stellar oscillation code for modeling tides in binary systems. Using a multipolar expansion in space and a Fourier-series expansion in time, we decompose the tidal potential into a superposition of partial tidal potentials. The equations governing the small-amplitude response of a spherical star to an individual partial potential are the linear, non-radial, nonadiabatic oscillation equations with an extra inhomogeneous forcing term. We introduce a new executable,
- Award ID(s):
- 1663696
- NSF-PAR ID:
- 10400480
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 945
- Issue:
- 1
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 43
- Size(s):
- ["Article No. 43"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract MESA binary simulations. Tidal heating is self-consistently computed and added to the models at every time step. As a binary inspirals to orbital periods less than ∼20 minutes, the WD’s behavior changes from cooling to heating, with temperature enhancements that can exceed 10,000 K compared with nontidally heated models. We compare a grid of tidally heated WD models to observed short-period systems with hot WD primaries. While tidal heating affects theirT eff, it is likely not the dominant luminosity. Instead, these WDs are probably intrinsically young and hot, implying that the binaries formed at short orbital periods. The binaries are consistent with undergoing common envelope evolution with a somewhat low efficiencyα CE. We delineate the parameter space where the traveling wave assumption is most valid, noting that it breaks down for WDs that cool sufficiently, where standing waves may instead be formed. -
more » « less
Abstract The salinity structure in an estuary is controlled by time‐dependent mixing processes. However, the locations and temporal variability of where significant mixing occurs is not well‐understood. Here we utilize a tracer variance approach to demonstrate the spatial and temporal structure of salinity mixing in the Hudson River Estuary. We run a 4‐month hydrodynamic simulation of the tides, currents, and salinity that captures the spring‐neap tidal variability as well as wind‐driven and freshwater flow events. On a spring‐neap time scale, salinity variance dissipation (mixing) occurs predominantly during the transition from neap to spring tides. On a tidal time scale, 60% of the salinity variance dissipation occurs during ebb tides and 40% during flood tides. Spatially, mixing during ebbs occurs primarily where lateral bottom salinity fronts intersect the bed at the transition from the main channel to adjacent shoals. During ebbs, these lateral fronts form seaward of constrictions located at multiple locations along the estuary. During floods, mixing is generated by a shear layer elevated in the water column at the top of the mixed bottom boundary layer, where variations in the along channel density gradients locally enhance the baroclinic pressure gradient leading to stronger vertical shear and more mixing. For both ebb and flood, the mixing occurs at the location of overlap of strong vertical stratification and eddy diffusivity, not at the maximum of either of those quantities. This understanding lends a new insight to the spatial and time dependence of the estuarine salinity structure.
-
more » « less
Abstract We examine the reasons for discrepancies between two alternative approaches to modeling small-amplitude tides in binary systems. The direct solution (DS) approach solves the governing differential equations and boundary conditions directly, while the modal decomposition (MD) approach relies on a normal-mode expansion. Applied to a model for the primary star in the heartbeat system KOI-54, the two approaches predict quite different behavior of the secular tidal torque. The MD approach exhibits the pseudosynchronization phenomenon, where the torque due to the equilibrium tide changes sign at a single, well-defined, and theoretically predicted stellar rotation rate. The DS approach instead shows “blurred” pseudosynchronization, where positive and negative torques intermingle over a range of rotation rates. We trace a major source of these differences to an incorrect damping coefficient in the profile functions describing the frequency dependence of the MD expansion coefficients. With this error corrected, some differences between the approaches remain; however, both are in agreement that pseudosynchronization is blurred in the KOI-54 system. Our findings generalize to any type of star for which the tidal damping depends explicitly or implicitly on the forcing frequency.
-
more » « less
Abstract The mechanisms and geographic distribution of global tidal dissipation in barotropic tidal models are examined using a high resolution unstructured mesh finite element model. Mesh resolution varies between 2 and 25 km and is especially focused on inner shelves and steep bathymetric gradients. Tidal response sensitivities to bathymetric changes are examined to put into context response sensitivities to frictional processes. We confirm that the Ronne Ice Shelf dramatically affects Atlantic tides but also find that bathymetry in the Hudson Bay system is a critical control. We follow a sequential frictional parameter optimization process and use TPXO9 data‐assimilated tidal elevations as a reference solution. From simulated velocities and depths, dissipation within the global model is estimated and allows us to pinpoint dissipation at high resolution. Boundary layer dissipation is extremely focused with 1.4% of the ocean accounting for 90% of the total. Internal tide friction is much more distributed with 16.7% of the ocean accounting for 90% of the total. Often highly regional dissipation can impact basin‐scale and even ocean wide tides. Optimized boundary layer friction parameters correlate very well with the physical characteristics of the locality with high friction factors associated with energetic tidal regions, deep ocean island chains, and ice covered areas. Global complex M2tide errors are 1.94 cm in deep waters. Total global boundary layer and internal tide dissipation are estimated, respectively, at 1.83 and 1.49 TW. This continues the trend in the literature toward attributing more dissipation to internal tides.
-
more » « less
ABSTRACT Many Sun-like stars are observed to host close-in super-Earths (SEs) as part of a multiplanetary system. In such a system, the spin of the SE evolves due to spin–orbit resonances and tidal dissipation. In the absence of tides, the planet’s obliquity can evolve chaotically to large values. However, for close-in SEs, tidal dissipation is significant and suppresses the chaos, instead driving the spin into various steady states. We find that the attracting steady states of the SE’s spin are more numerous than previously thought, due to the discovery of a new class of ‘mixed-mode’ high-obliquity equilibria. These new equilibria arise due to subharmonic responses of the parametrically driven planetary spin, an unusual phenomenon that arises in non-linear systems. Many SEs should therefore have significant obliquities, with potentially large impacts on the physical conditions of their surfaces and atmospheres.
