A revised derivation of the discharge coefficient for flows over thin weirs and sills in the limits of wall overflow to a free overfall is given. Using dimensional analysis, we show that the discharge coefficient, $C_{d}$, in the classical weir-discharge equation is best understood as a weir Froude number, ${Fr}_{h}$, which accounts for the combined effects of inertia, contraction and viscous energy losses within the flow field. A comprehensive set of experimental data from historical studies is complimented by new data from the authors, featuring both laboratory flume experiments and three-dimensional numerical simulations of weir flows. Synthesis of these data elucidates the interaction between the coupled pressure and velocity fields, and the balance between inertial and contraction effects as ${Fr}_{h}$ varies. Analysis of the vertical pressure gradient reveals that the thickness of the nappe initially widens with increasing inertia, but then contracts again towards the free overfall limit due to diminishing flow separation at the base of the weir. These insights allow for a physical explanation of the transition between weir and sill flows using the channel Froude number. Practical limitations on predicting weir discharge and a description of characteristic flow regimes are also set forward.

We investigate the evolution of active galactic nucleus jets on kiloparsec-scales due to their interaction with the clumpy interstellar medium (ISM) of the host galaxy and, subsequently, the surrounding circumgalactic environment. Hydrodynamic simulations of this jet–environment interaction are presented for a range of jet kinetic powers, peak densities of the multiphase ISM, and scale radii of the larger-scale environment – characteristic of either a galaxy cluster or poor group. Synthetic radio images are generated by considering the combination of synchrotron radiation from the jet plasma and free-free absorption from the multiphase ISM. We find that jet propagation is slowed by interactions with a few very dense clouds in the host galaxy ISM, producing asymmetries in lobe length and brightness which persist to scales of tens of kpc for poor group environments. The classification of kiloparsec-scale jets is highly dependent on surface brightness sensitivity and resolution. Our simulations of young active sources can appear as restarted sources, showing double-double lobe morphology, high core prominence (CP $\gt 0.1$), and the expected radio spectra for both the inner- and outer-lobe components. We qualitatively reproduce the observed inverse correlation between peak frequency and source size and find that the peak frequency of the integrated radio spectrum depends on ISM density but not the jet power. Spectral turnover in resolved young radio sources therefore provides a new probe of the ISM.

Bow shocks generated by pulsars moving through weakly ionized interstellar medium (ISM) produce emission dominated by non-equilibrium atomic transitions. These bow shocks are primarily observed as H$\alpha$ nebulae. We developed a package, named Shu, that calculates non-LTE intensity maps in more than 150 spectral lines, taking into account geometrical properties of the pulsars’ motion and lines of sight. We argue here that atomic (C i, N i, O i) and ionic (S ii, N ii, O iii, Ne iv) transitions can be used as complementary and sensitive probes of ISM. We perform self-consistent 2D relativistic hydrodynamic calculations of the bow shock structure and generate non-LTE emissivity maps, combining global dynamics of relativistic flows, and detailed calculations of the non-equilibrium ionization states. We find that though typically $\text{H}_\alpha$ emission is dominant, spectral fluxes in [O iii], [S ii] and [N ii] may become comparable for relatively slowly moving pulsars. Overall, morphology of non-LTE emission, especially of the ionic species, is a sensitive probe of the density structures of the ISM.

Before a binary system enters into a common envelope (CE) phase, accretion from the primary star onto the companion star through Roche Lobe overflow (RLOF) will lead to the formation of an accretion disk, which may generate jets. Accretion before and during the CE may alter the outcome of the interaction. Previous studies have considered different aspects of this physical mechanism. Here we study the properties of an accretion disk formed via 3D hydrodynamic simulations of the RLOF mass transfer between a 7 M$_{\odot}$, red supergiant star and a 1.4 M$_{\odot}$, neutron star companion. We simulate only the volume around the companion for improved resolution. We use a 1D implicit mesa simulation of the evolution of the system during 30 000 yr between the on-set of the RLOF and the CE to guide the binary parameters and the mass-transfer rate, while we simulate only 21 yr of the last part of the RLOF in 3D using an ideal gas quasi-isothermal equation of state. We expect that a pre-CE disk under these parameters will have a mass of $\sim 5\times 10^{-3}$ M$_{\odot}$ and a radius of $\sim40\ R_\odot$ with a scale height of $\sim 5\ R$$_{\odot}$. The temperature profile of the disk is shallower than that predicted by the formalism of Shakura and Sunyaev, but more reasonable cooling physics would need to be included. We stress test these results with respect to a number of physical and numerical parameters, as well as simulation choices, and we expect them to be reasonable within a factor of a few for the mass and 15% for the radius. We also contextualise our results within those presented in the literature, in particular with respect to the dimensionality of simulations and the adiabatic index. We discuss the measured accretion rate in the context of the Shakura and Sunyaev formalism and debate the viscous mechanisms at play, finishing with a list of prospects for future work.

We introduce adaptive particle refinement for compressible smoothed particle hydrodynamics (SPH). SPH calculations have the natural advantage that resolution follows mass, but this is not always optimal. Our implementation allows the user to specify local regions of the simulation that can be more highly resolved. We test our implementation on practical applications including a circumbinary disc, a planet embedded in a disc, and a flyby. By comparing with equivalent globally high-resolution calculations, we show that our method is accurate and fast, with errors in the mass accreted onto sinks of less than 9% and speed ups of 1.07–6.62$\times$ for the examples shown. Our method is adaptable and easily extendable, for example, with multiple refinement regions or derefinement.

The innermost region of the Milky Way harbors the central molecular zone (CMZ). This region contains a large amount of molecular gas but a poor star formation rate considering the densities achieved by the gas in this region. We used the arepo code to perform a hydrodynamic and star formation simulation of the galaxy, where a Ferrers bar was adiabatically introduced. During the stage of bar imposition, the bar strength excites density waves close to the inner Lindblad resonance guiding material towards the inner galaxy, driving the formation of a ring that we qualitatively associate with the CMZ. During the simulation, we identified that the ring passes three main phases, namely: formation, instability, and quasi-stationary stages. During the whole evolution, and particularly in the quasi-stationary stage, we observe that the ring is associated with the x2 family of periodic orbits. Additionally, we found that most of the star formation occurs during the ring formation stage, while it drastically decreases in the instability stage. Finally, we found that when the gas has settled in a stable x2 orbit, the star formation takes place mostly after the dense gas passes the apocentre, triggering the conveyor-belt mechanism described in previous studies.

Odd Radio Circles (ORCs) are a class of low surface brightness, circular objects approximately one arcminute in diameter. ORCs were recently discovered in the Australian Square Kilometre Array Pathfinder (ASKAP) data and subsequently confirmed with follow-up observations on other instruments, yet their origins remain uncertain. In this paper, we suggest that ORCs could be remnant lobes of powerful radio galaxies, re-energised by the passage of a shock. Using relativistic hydrodynamic simulations with synchrotron emission calculated in post-processing, we show that buoyant evolution of remnant radio lobes is alone too slow to produce the observed ORC morphology. However, the passage of a shock can produce both filled and edge-brightnened ORC-like morphologies for a wide variety of shock and observing orientations. Circular ORCs are predicted to have host galaxies near the geometric centre of the radio emission, consistent with observations of these objects. Significantly offset hosts are possible for elliptical ORCs, potentially causing challenges for accurate host galaxy identification. Observed ORC number counts are broadly consistent with a paradigm in which moderately powerful radio galaxies are their progenitors.

An open problem of the derivation of the relativistic Vlasov equation for systems of charged particles moving with velocities up to the speed of light and creating the electromagnetic field in accordance with the full set of the Maxwell equations is considered. Moreover, the method of derivation is illustrated on the non-relativistic kinetic model. Independent derivation of the relativistic hydrodynamics is also demonstrated. The key role of these derivations of the hydrodynamic and kinetic equations includes the explicit operator of averaging on the physically infinitesimal volume suggested by L.S. Kuzmenkov.

Both observations and theoretical studies have convincingly shown that outflows (i.e., wind and jet) are common phenomena from black hole accretion systems with various accretion rates, although the physical driving mechanisms are not exactly same for different accretion modes. Outflows are not only important in the dynamics of black hole accretion, but also play an important role in AGN feedback; therefore it is crucial to investigate their main physical properties including mass flux and velocity. In this paper we summarize recent studies in investigating the properties and driving mechanisms from black hole accretion flows with various accretion rates.

Shallowly submerged oscillating structures may be found in wave energy devices or semi-submersible vessels. Predicting the force on such structures is critical for design purposes, but complicated due to nonlinear phenomena which can occur in shallow water, including wave breaking and bore formation. Such effects are particularly important around the first ‘resonance’ frequency of the fluid on top of the device, where linear theory predicts large flows on/off the cylinder and corresponding surface elevations and forces. In an effort to create a reliable and efficient model to predict the hydrodynamic force on a shallowly submerged truncated vertical cylinder, an axisymmetric nonlinear hybrid model is developed for forced heave oscillations. The flow above the cylinder is modelled using the nonlinear shallow water equations, and linear potential flow theory is used in the surrounding fluid. The model is compared with experimental results for forced heave oscillations and performs well for predicting the heave force. It is then used to examine linearised heave force for increasing amplitudes of (prescribed) harmonic heave motion. There is a significant reduction in the peaks of radiation damping and added mass coefficients with increasing amplitude, and associated shifts in the frequencies of the peaks.

Hydroinformatics is a technology that combines information and communications technologies together with various disciplinary optimization and simulation models that focus on the management of water. This paper reviews the historical development of hydroinformatics and summarizes the current state of this technology. It describes the range of modeling tools and applications currently described in hydroinformatics literature. The paper concludes with some speculations about possible future developments in hydroinformatics.

We present the Cosmological Double Radio Active Galactic Nuclei (CosmoDRAGoN) project: a large suite of simulated AGN jets in cosmological environments. These environments sample the intra-cluster media of galaxy clusters that form in cosmological smooth particle hydrodynamics (SPH) simulations, which we then use as inputs for grid-based hydrodynamic simulations of radio jets. Initially conical jets are injected with a range of jet powers, speeds (both relativistic and non-relativistic), and opening angles; we follow their collimation and propagation on scales of tens to hundreds of kiloparsecs, and calculate spatially resolved synthetic radio spectra in post-processing. In this paper, we present a technical overview of the project, and key early science results from six representative simulations which produce radio sources with both core- (Fanaroff-Riley Type I) and edge-brightened (Fanaroff-Riley Type II) radio morphologies. Our simulations highlight the importance of accurate representation of both jets and environments for radio morphology, radio spectra, and feedback the jets provide to their surroundings.