A neutron star in a compact binary is expected to be well-approximated by a barotropic flow during the inspiral phase. During the merger phase, where tidal disruption and shock-heating occur, a baroclinic description is needed instead. In the barotropic case, a Hamiltonian formulation potentially offers unique benefits for numerical relativity simulations of the inspiral phase, including highly accurate conservation of circulation and superconvergence of the fluid variables, and is actively being explored. In this work, we investigate the viability of a Hamiltonian formulation in the baroclinic case. At odds with the barotropic case, this formulation is non-conservative, yet it can be treated well with approximate Riemann solver algorithms since the non-conservative terms vanish across genuinely nonlinear fields. Nonetheless, using numerical 1-dimensional shock tube tests we find that the weak solutions of the Hamiltonian system differ from the standard ones obtained by enforcing conservation of rest mass density, momentum density, and energy density across discontinuities. We also show that barotropic Hamiltonian formulations can admit shockwaves at fluid-vacuum interfaces, which may be related to the unstable behavior of stellar surfaces observed in past numerical tests. In light of the unphysical weak solutions, we expect that in future implementations of the Hamiltonian formulation of hydrodynamics in numerical relativity it will be necessary to use an explicitly barotropic formulation during the inspiral phase, and then switch to a robust baroclinic formulation prior to merger.

We write down and apply the linearized fluid and gravitational equations consistent with pseudo-Newtonian simulations, whereby Newtonian hydrodynamics is used with a pseudo-Newtonian monopole and standard Newtonian gravity for higher multipoles. We thereby eliminate the need to use mode function matching to identify the active non-radial modes in pseudo-Newtonian core-collapse supernova simulations, in favor of the less complex and less costly mode frequency matching method. In doing so, we are able to measure appropriate boundary conditions for a mode calculation.

The dynamics of self-gravitating fluid bodies is described by the Euler-Einstein system of partial differential equations. The break-down of well-posedness on the fluid-vacuum interface remains a challenging open problem, which is manifested in simulations of oscillating or inspiraling binary neutron-stars. We formulate and implement a well-posed canonical hydrodynamic scheme, suitable for neutron-star simulations in numerical general relativity. The scheme uses a variational principle by Carter-Lichnerowicz stating that barotropic fluid motions are conformally geodesic and Helmholtz's third theorem stating that initially irrotational flows remain irrotational. We apply this scheme in 3+1 numerical general relativity to evolve the canonical momentum of a fluid element via the Hamilton-Jacobi equation. We explore a regularization scheme for the Euler equations, that uses a fiducial atmosphere in hydrostatic equilibrium and allows the pressure to vanish, while preserving strong hyperbolicity on the vacuum boundary. The new regularization scheme resolves a larger number of radial oscillation modes compared to standard, non-equilibrium atmosphere treatments.

We investigate correlated gravitational wave and neutrino signals from rotating core-collapse supernovae with simulations. Using an improved mode identification procedure based on mode function matching, we show that a linear quadrupolar mode of the core produces a dual imprint on gravitational waves and neutrinos in the early post-bounce phase of the supernova. The angular harmonics of the neutrino emission are consistent with the mode energy around the neutrinospheres, which points to a mechanism for the imprint on neutrinos. Thus, neutrinos carry information about the mode amplitude in the outer region of the core, whereas gravitational waves probe deeper in. We also find that the best-fit mode function has a frequency bounded above by $\sim 420$ Hz, and yet the mode's frequency in our simulations is $\sim 15\%$ higher, due to the use of Newtonian hydrodynamics and a widely used pseudo-Newtonian gravity approximation. This overestimation is particularly important for the analysis of gravitational wave detectability and asteroseismology, pointing to limitations of pseudo-Newtonian approaches for these purposes, possibly even resulting in excitation of incorrect modes. In addition, mode frequency matching (as opposed to mode function matching) could be resulting in mode misidentification in recent work. Lastly, we evaluate the prospects of a multimessenger detection of the mode using current technology. The detection of the imprint on neutrinos is most challenging, with a maximum detection distance of $\sim\!1$ kpc using the IceCube Neutrino Observatory. The maximum distance for detecting the complementary gravitational wave imprint is $\sim\!5$ kpc using Advanced LIGO at design sensitivity.

Part IA: We present numerical measurements of relativistic scaling relations in $(2+1)$-dimensional conformal fluid turbulence, which perform favourably over their non-relativistic versions. As seen with incompressible turbulence in past studies, we find that the energy spectrum exhibits $k^{-2}$ scaling rather than the Kolmogorov/Kraichnan expectation of $k^{-5/3}$. Part IB: We compute the fractal dimension $D$ of a turbulent anti-deSitter black brane reconstructed from boundary fluid data using the fluid-gravity duality. Our value of $D=2.584(1)$ is consistent with the upper bound $D\leq 3$, resolving a recent claim that $D=3+1/3$. We describe how to covariantly define the fractal dimension of spatial sections of the horizon, and we speculate on assigning a `bootstrapped' value to the entire horizon. Part II: We report progress implementing a fluid code with post-Newtonian (PN) gravity in spherical symmetry. The PN formalism couples a fluid, its self-gravity, and a black hole via elliptic equations. This eliminates radiative modes, allowing larger time steps, which is helpful for studying systems with very long time scales, eg. tidal disruption events. Part III: Asteroseismology of rotating core-collapse supernovae is possible with a multimessenger strategy. We show an $l=2$, $m=0$, $n\gtrsim 2$, $f\lesssim 280$ Hz mode of the core is responsible for emission in gravitational waves and neutrinos. The angular harmonics of the neutrino emission is consistent with the mode energy around the neutrinospheres, where $r\sim 70$ km. Thus, neutrinos carry information about the mode in the outer region of the core, whereas gravitational waves probe the deep inner core $r \lesssim 30$ km.

We derive exact scaling relations for two-dimensional relativistic hydrodynamic turbulence in the inertial range of scales. We consider both the energy cascade towards large scales and the enstrophy cascade towards small scales. We illustrate these relations by numerical simulations of turbulent weakly compressible flows. Intriguingly, the fluid-gravity correspondence implies that the gravitational field in spacetimes with anti-de Sitter asymptotics should exhibit similar scaling relations.