We calculate the motion of binary mass systems in gravity up to the sixth post--Newtonian order to the $G_N^3$ terms ab initio using momentum expansions within an effective field theory approach based on Feynman amplitudes in harmonic coordinates. For these contributions we construct a canonical transformation to isotropic and to EOB coordinates at 5PN and agree with the results in the literature \cite{Bern:2019nnu,Damour:2019lcq}. At 6PN we compare to the Hamiltonians in isotropic coordinates either given in \cite{Bern:2019nnu} or resulting from the scattering angle. We find a canonical transformation from our Hamiltonian in harmonic coordinates to \cite{Bern:2019nnu}, but not to \cite{Damour:2019lcq}. This implies that we also agree on all observables with \cite{Bern:2019nnu} to the sixth post--Newtonian order to $G_N^3$.

We calculate the motion of binary mass systems in gravity up to the fourth post--Newtonian order. We use momentum expansions within an effective field theory approach based on Feynman amplitudes in harmonic coordinates by applying dimensional regularization. We construct the canonical transformations to ADM coordinates and to effective one body theory (EOB) to compare with other approaches. We show that intermediate poles in the dimensional regularization parameter $\varepsilon$ vanish in the observables and the classical theory is not renormalized. The results are illustrated for a series of observables for which we agree with the literature.

The post-Newtonian and post-Minkowskian solutions for the motion of binary mass systems in gravity can be derived in terms of momentum expansions within effective field theory approaches. In the post-Minkowskian approach the expansion is performed in the ratio $G_N/r$, retaining all velocity terms completely, while in the post-Newtonian approach only those velocity terms are accounted for which are of the same order as the potential terms due to the virial theorem. We show that it is possible to obtain the complete post-Minkowskian expressions completely algorithmically, under most general purely mathematical conditions from a finite number of velocity terms and illustrate this up to the third post-Minkowskian order given in \cite{Bern:2019crd}.

FEARLESS (Fluid mEchanics with Adaptively Refined Large Eddy SimulationS) is a new numerical scheme arising from the combined use of subgrid scale (SGS) model for turbulence at the unresolved length scales and adaptive mesh refinement (AMR) for resolving the large scales. This tool is especially suitable for the study of turbulent flows in strongly clumped media. In this contribution, the main features of FEARLESS are briefly outlined. We then summarize the main results of FEARLESS cosmological simulations of galaxy cluster evolution. In clusters, the production of turbulence is closely correlated with merger events; for minor mergers, we find that turbulent dissipation affects the cluster energy budget only locally. The level of entropy in the cluster core is enhanced in FEARLESS simulations, in accord with a better modeling of the unresolved flow, and with its feedback on the resolved mixing in the ICM.

We present a numerical scheme for modelling unresolved turbulence in cosmological adaptive mesh refinement codes. As a first application, we study the evolution of turbulence in the intra-cluster medium and in the core of a galaxy cluster. Simulations with and without subgrid scale model are compared in detail. Since the flow in the ICM is subsonic, the global turbulent energy contribution at the unresolved length scales is smaller than 1% of the internal energy. We find that the production of turbulence is closely correlated with merger events occurring in the cluster environment, and its dissipation locally affects the cluster energy budget. Because of this additional source of dissipation, the core temperature is larger and the density is smaller in the presence of subgrid scale turbulence than in the standard adiabatic run, resulting in a higher entropy core value.

In the context of explosion models for Type Ia Supernovae, we present one- and two-dimensional simulations of fully resolved detonation fronts in degenerate C+O White Dwarf matter including clumps of previously burned material. The ability of detonations to survive the passage through sheets of nuclear ashes is tested as a function of the width and composition of the ash region. We show that detonation fronts are quenched by microscopically thin obstacles with little sensitivity to the exact ash composition. Front-tracking models for detonations in macroscopic explosion simulations need to include this effect in order to predict the amount of unburned material in delayed detonation scenarios.