To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
We propose a computational framework for simulating the self-similar regime of turbulent Rayleigh–Taylor (RT) mixing layers in a statistically stationary manner. By leveraging the anticipated self-similar behaviour of RT mixing layers, a transformation of the vertical coordinate and velocities is applied to the Navier–Stokes equations (NSE), yielding modified equations that resemble the original NSE but include two sets of additional terms. Solving these equations, a statistically stationary RT (SRT) flow is achieved. Unlike temporally growing Rayleigh–Taylor (TRT) flow, SRT flow is independent of initial conditions and can be simulated over infinite simulation time without escalating resolution requirements, hence guaranteeing statistical convergence. Direct numerical simulations (DNS) are performed at an Atwood number of 0.5 and unity Schmidt number. By varying the ratio of the mixing layer height to the domain width, a minimal flow unit of aspect ratio 1.5 is found to approximate TRT turbulence in the self-similar mode-coupling regime. The SRT minimal flow unit has one-sixteenth the number of grid points required by the equivalent TRT simulation of the same Reynolds number and grid resolution. The resultant flow corresponds to a theoretical limit where self-similarity is observed in all fields and across the entire spatial domain – a late-time state that existing experiments and DNS of TRT flow have difficulties attaining. Simulations of the SRT minimal flow unit span TRT-equivalent Reynolds numbers (based on mixing layer height) ranging from 500 to 10 800. The SRT results are validated against TRT data from this study as well as from Cabot & Cook (Nat. Phys., vol. 2, 2006, pp. 562–568).
There is growing evidence that the broadband radio spectral energy distributions (SEDs) of star-forming galaxies (SFGs) contain a wealth of complex physics. In this paper we aim to determine the physical emission and loss processes causing radio SED curvature and steepening to see what observed global astrophysical properties, if any, are correlated with radio SED complexity. To do this, we have acquired radio continuum data between 70 MHz and 17 GHz for a sample of 19 southern local ($z \lt 0.04$) SFGs. Of this sample 11 are selected to contain low-frequency ($ \lt $300 MHz) turnovers (LFTOs) in their SEDs and eight are control galaxies with similar global properties. We model the radio SEDs for our sample using a Bayesian framework whereby radio emission (synchrotron and free-free) and absorption or loss processes are included modularly. We find that without the inclusion of higher frequency data ($ \gt $17 GHz) single synchrotron power-law based models are always preferred for our sample; however, additional processes including free-free absorption (FFA) and synchrotron losses are often required to accurately model radio SED complexity in SFGs. The fitted synchrotron spectral indices range from $-0.45$ to $-1.07$ and are strongly anticorrelated with stellar mass suggesting that synchrotron losses are the dominant mechanism acting to steepen the spectral index in larger/more massive nearby SFGs. We find that LFTOs in the radio SED are independent from the inclination of SFGs; however, higher inclination galaxies tend to have steeper fitted spectral indices indicating losses to diffusion of cosmic ray electrons into the galactic halo. Four of five of the merging systems in our SFG sample have elevated specific star formation rates and flatter fitted spectral indices with unconstrained LFTOs. Lastly, we find no significant separation in global properties between SFGs with or without modelled LFTOs. Overall these results suggest that LFTOs are likely caused by a combination of FFA and ionisation losses in individual recent starburst regions with specific orientations and interstellar medium properties that, when averaged over the entire galaxy, do not correlate with global astrophysical properties.
Fluid dynamics systems driven by dominant, near-periodic dynamics are common across wakes, jets, rotating machinery and high-speed flows. Traditional modal decomposition techniques have been used to gain insight into these flows, but can require many modes to represent key physical processes. With the aim of generating modes that intuitively convey the underlying physical mechanisms, we propose an intrinsic phase-based proper orthogonal decomposition (IPhaB POD) method, which creates energetically ranked modes that evolve along a characteristic cycle of the dominant near-periodic dynamics. Our proposed formulation is set in the time domain, which is particularly useful in cases where the cyclical content is imperfectly periodic. We formally derive IPhaB POD within a POD framework that therefore inherits the energetically ranked decomposition and optimal low-rank representation inherent to POD. As part of this derivation, a dynamical systems representation is utilized, facilitating a definition of phase within the system's near-periodic cycle in the time domain. An expectation operator and inner product are also constructed relative to this definition of phase in a manner that allows for the various cycles within the data to demonstrate imperfect periodicity. The formulation is tested on two sample problems: a simple, low Reynolds number aerofoil wake, and a complex, high-speed pulsating shock wave problem. The method is compared to space-only POD, spectral POD (SPOD) and cyclostationary SPOD. The method is shown to better isolate the dominant, near-periodic global dynamics in a time-varying IPhaB mean, and isolate the tethered, local-in-phase dynamics in a series of time-varying modes.
Flow field in the near wake of a short-finite circular cylinder at $L/D=1.0$ with an angle of attack between 0$^\circ$–15$^\circ$, where the transition from the non-reattaching flow to the reattaching flow appears, is investigated in wind tunnel tests with a supportless condition. Stereo particle image velocimetry measurements were applied to the experiments at the Reynolds number of $3.46\times 10^4$, and velocity fields in the near wake were obtained. The data was mainly analysed using spectral proper orthogonal decomposition. Characteristic large-scale wake structures of recirculation bubble pumping and large-scale vortex shedding were identified in the near wake of the cylinder regardless of the angle of attack. The phase difference of expansion and contraction of the recirculation flow appears in the recirculation bubble pumping at $\alpha \neq 0^\circ$. On the other hand, the eigenfunctions of velocity fluctuations at the vortex shedding frequency show a similar spatial pattern regardless of $\alpha$. Frequency analyses of wake position calculated from the reconstructed velocity field clarified that peak frequency is different between two in-plane directions when $\alpha \neq 0^\circ$. In addition, three vortex shedding patterns (anticlockwise/clockwise circular and flapping) are identified not only at $\alpha =0^\circ$ but also $\alpha \neq 0^\circ$. The feature of wake position in the radial direction for each pattern is observed regardless of the angle of attack. The relationship between the recirculation bubble pumping and the wake position in the radial direction is apparent in the non-reattaching flow but is weaker with $\alpha$ in the reattaching flow.
Observations of Galactic supernova remnants (SNRs) are crucial to understanding supernova explosion mechanisms and their impact on our Galaxy’s evolution. SNRs are usually identified by searching for extended, circular structures in all-sky surveys. However, the resolution and sensitivity of any given survey results in selection biases related to the brightness and angular scale of a subset of the total SNR population. As a result, we have only identified 1/3 of the expected number of SNRs in our Galaxy. We used data collected by the Murchison Widefield Array (MWA) to perform a visual search for SNR candidates over $ 285^{\circ} \lt l \lt 70^{\circ}$ and $|b| \lt 16^{\circ}$. We then used the Widefield Infrared Survey Explorer to eliminate likely Hii regions from our SNR candidate sample. By exploiting the resolution and sensitivity of MWA data, we have successfully detected 10 new candidates using our proposed method. In addition, our method has also enabled us to detect and verify 10 previously known but unconfirmed candidates. The 20 SNR candidates described in the paper will increase the known SNR population in the Galaxy by 7%.
Super-resolution of turbulence is a term used to describe the prediction of high-resolution snapshots of a flow from coarse-grained observations. This is typically accomplished with a deep neural network and training usually requires a dataset of high-resolution images. An approach is presented here in which robust super-resolution can be performed without access to high-resolution reference data, as might be expected in an experiment. The training procedure is similar to data assimilation, wherein the model learns to predict an initial condition that leads to accurate coarse-grained predictions at later times, while only being shown coarse-grained observations. Implementation of the approach requires the use of a fully differentiable flow solver in the training loop to allow for time-marching of predictions. A range of models are trained on data generated from forced, two-dimensional turbulence. The networks have reconstruction errors which are similar to those obtained with ‘standard’ super-resolution approaches using high-resolution data. Furthermore, the methods are comparable to the performance of standard data assimilation for state estimation on individual trajectories, outperforming these variational approaches at initial time and remaining robust when unrolled in time where performance of the standard data-assimilation algorithm improves.
Using the ONEDFEL code we perform free electron laser simulations in the astrophysically important guide-field dominated regime. For wigglers’ (Alfvén waves) wavelengths of tens of kilometres and beam Lorentz factor ${\sim }10^3$, the resulting coherently emitted waves are in the centimetre range. Our simulations show a growth of the wave intensity over fourteen orders of magnitude, over the astrophysically relevant scale of approximately a few kilometres. The signal grows from noise (unseeded). The resulting spectrum shows fine spectral substructures, reminiscent of those observed in fast radio bursts.
The regulation of electron heat transport in high-$\beta$, weakly collisional, magnetized plasma is investigated. A temperature gradient oriented along a mean magnetic field can induce a kinetic heat-flux-driven whistler instability (HWI), which back-reacts on the transport by scattering electrons and impeding their flow. Previous analytical and numerical studies have shown that the heat flux for the saturated HWI scales as $\beta _e^{-1}$. These numerical studies, however, had limited scale separation and consequently large fluctuation amplitudes, which calls into question their relevance at astrophysical scales. To this end, we perform a series of particle-in-cell simulations of the HWI across a range of $\beta _e$ and temperature-gradient length scales under two different physical set-ups. The saturated heat flux in all of our simulations follows the expected $\beta _e^{-1}$ scaling, supporting the robustness of the result. We also use our simulation results to develop and implement several methods to construct an effective collision operator for whistler turbulence. The results point to an issue with the standard quasi-linear explanation of HWI saturation, which is analogous to the well-known $90^{\circ }$ scattering problem in the cosmic-ray community. Despite this limitation, the methods developed here can serve as a blueprint for future work seeking to characterize the effective collisionality caused by kinetic instabilities.
Ion cyclotron resonance heating is a versatile heating method that has been demonstrated to be able to efficiently couple power directly to the ions via the fast magnetosonic wave. However, at temperatures relevant for reactor grade devices such as DEMO, electron damping becomes increasingly important. To reduce electron damping, it is possible to use an antenna with a power spectrum dominated by low parallel wavenumbers. Moreover, using an antenna with a unidirectional spectrum, such as a travelling wave array antenna, the parallel wavenumber can be downshifted by mounting the antenna in an elevated position relative to the equatorial plane. This downshift can potentially enhance ion heating as well as fast wave current drive efficiency. Thus, such a system could benefit ion heating during the ramp-up phase and be used for current drive during flat-top operation. To test this principle, both ion heating and current drive have been simulated in a DEMO-like plasma for a few different mounting positions of the antenna using the FEMIC code. We find that moving the antenna off the equatorial plane makes ion heating more efficient for all considered plasma temperatures at the expense of on-axis heating. Moreover, although current drive efficiency is enhanced, electron damping is reduced for lower mode numbers, thus reducing the driven current in this part of the spectrum.
We apply a continuation method to recently optimized stellarator equilibria with excellent quasi-axisymmetry to generate new equilibria with a wide range of rotational transform profiles. Using these equilibria, we investigate how the rotational transform affects fast-particle confinement, the maximum coil–plasma distance, the maximum growth rate in linear gyrokinetic ion-temperature gradient simulations and the ion heat flux in corresponding nonlinear simulations. We find values of two-term quasi-symmetry error comparable to or lower than those of the similar Landreman–Paul (Phys. Rev. Lett., vol. 128, 2022, 035001) configuration for values of the mean rotational transform $\bar {\iota }$ between $0.12$ and $0.75$. The fast-particle confinement improves with $\bar {\iota }$ until $\bar {\iota } = 0.73$, at which point the degradation in quasi-symmetry outweighs the benefits of further increasing $\bar {\iota }$. The required coil–plasma distance only varies by about ${\pm }10\,\%$ for the configurations under consideration, and is between $2.8$ and $3.3\ \mathrm {m}$ when the configuration is scaled up to reactor size (minor radius $a=1.7\ \mathrm {m}$ and volume-averaged magnetic field strength of $5.86\ \mathrm {T}$). The maximum growth rate from linear gyrokinetic simulations increases with $\bar {\iota }$, but also shifts towards higher $k_y$ values. The maximum linear growth rate is sensitive to the choice of flux tube at rational $\bar {\iota }$, but this can be compensated for by taking the maximum over several flux tubes. The corresponding ion heat fluxes from nonlinear simulations display a non-monotonic relation to $\bar {\iota }$. Sufficiently large positive shear is destabilizing. This is reflected in both linear growth rates and nonlinear heat fluxes.
We present an innovative design for a two-head, gas-cooled multi-slab high-energy, high-repetition-rate amplifier aimed at mitigating thermally induced depolarization in a wide-bandwidth neodymium-doped glass gain medium. This architecture employs two quartz rotators (QRs) with opposite-handedness, strategically positioned within each multi-slab amplifier head, to enhance depolarization compensation. Theoretical modeling of this amplifier configuration demonstrates a 20× reduction in depolarization losses for a 70 mm beam operating at the central wavelength, compared to conventional approaches that utilize a single QR positioned between the amplifier heads. In addition, for a wide bandwidth source, the integration of QRs with opposite-handedness yields a 9× improvement in depolarization losses at the spectral extremes compared to the use of two QRs exhibiting the same optical handedness in both amplifier heads.
The multi-colour complete light curves and low-resolution spectra of two short period eclipsing Am binaries V404 Aur and GW Gem are presented. The stellar atmospheric parameters of the primary stars were derived through the spectra fitting. The observed and TESS-based light curves of them were analysed by using the Wilson-Devinney code. The photometric solutions suggest that both V404 Aur and GW Gem are semi-detached systems with the secondary component filling its critical Roche Lobe, while the former should be a marginal contact binary. The $O-C$ analysis found that the period of V404 Aur is decreasing at a rate of $dP/dt=-1.06(\pm0.01)\times 10^{-7}\,\mathrm{d}\,\mathrm{ yr}^{-1}$, while the period of GW Gem is increasing at $dP/dt=+2.41(\pm0.01)\times 10^{-8} \mathrm{d}\,\mathrm{yr}^{-1}$. The period decrease of V404 Aur may mainly be caused by the combined effects of the angular momentum loss (AML) via an enhanced stellar wind of the more evolved secondary star and mass transfer between two components. The period increase of GW Gem supports the mass transfer from the secondary to the primary. Both targets may be in the broken contact stage predicted by the thermal relaxation oscillations theory and will eventually evolve to the contact stage. We have collected about 54 well-known eclipsing Am binaries with absolute parameters from the literature. The relations of these parameters are summarised. There are some components that have a higher degree of evolution. The majority of their hydrogen shell may have been stripped away and the stellar internal layer exposed. The accretion processes from such evolved components may be very important for the formation of Am peculiarity in binaries.
The stability and dynamics of solitary waves propagating along the surface of an inviscid ferrofluid jet in the absence of gravity are investigated analytically and numerically. For the axisymmetric geometry, the problem is shown to be a conservative system with total energy as the Hamiltonian; however, one of the canonical variables differs from those in the classic water-wave problem in the Cartesian coordinate system. The Dirichlet–Neumann operator appearing in the kinetic energy is then expanded as a Taylor series, described in homogeneous powers of the surface displacement. Based on the further analysis of the Dirichlet–Neumann operator, a systematic procedure is proposed to derive reduced model equations of multiple scales in various asymptotic limits from the full Euler equations in the Hamiltonian/Lagrangian framework. Particularly, a fully dispersive model arising from retaining terms valid up to the quartic order in the series expansion of the kinetic energy, which results in quadratic and cubic algebraic nonlinearities in Hamilton's equations and henceforth is abbreviated as the cubic full-dispersion model, is proposed. By comparing bifurcation curves and wave profiles of various types of axisymmetric solitary waves among different model equations, the cubic full-dispersion model is found to agree well with the full Euler equations, even for waves of considerably large amplitudes. The stability properties of axisymmetric solitary waves subjected to longitudinal disturbances are verified with the newly proposed model. Our analytical results, consistent with Saffman's theory, indicate that in the axisymmetric cylindrical system, the stability exchange subjected to superharmonic perturbations also occurs at the stationary point of the speed-energy bifurcation curve. A series of numerical experiments for the stability and dynamics of solitary waves are performed via the numerical time integration of the model equation, and collision interactions between stable solitary waves show non-elastic features.
We demonstrate a low-cost radio telescope using a 2.4-m satellite dish, an inexpensive printer circuit board (PCB)-based dual-pole antenna and commodity-off-the-shelf components. Open-source, radio acquisition and professional pulsar processing tools are used to successfully monitor on a daily basis the Southern Hemisphere pulsar J0835−4510 (the Vela pulsar) at 820 MHz as it transits the meridian. The system successfully detected a ‘glitch’ in the Vela pulsar in real-time at Woodchester Observatory in South Australia. Woodchester represents a good balance of bandwidth, observation time, and cost to achieve scientifically interesting results on pulsar timing for amateur radio astronomy and/or STEM outreach projects.
Granular column collapse is a simple but important problem to the granular material community, due to its links to dynamics of natural hazards, such as landslides and pyroclastic flows, and many industrial situations, as well as its potential of analysing transient and non-local rheology of granular flows. This article proposes a new dimensionless number to describe the run-out behaviour of granular columns on inclined planes based on both previous experimental data and dimensional analysis. With the assistance of the sphero-polyhedral discrete element method (DEM), we simulate inclined granular column collapses with different initial aspect ratios, particle contact properties and initial solid fractions on inclined planes with different inclination angles ($2.5^{\circ }\unicode{x2013}20.0^{\circ }$) to verify the proposed dimensional analysis. Detailed analyses are further provided for better understanding of the influence of different initial conditions and boundary conditions, and to help unify the description of the run-out scaling of systems with different inclination angles. This work determines the similarity and unity between granular column collapses on inclined planes and those on horizontal planes, and helps investigate the transient rheological behaviour of granular flows, which has direct relevance to various natural and engineering systems.
Two major areas of modern radio astronomy, namely, explosive astrophysical transient phenomena and observations of cosmological structures, are driving the design of aperture arrays towards large numbers of low-cost elements consisting of multiple spatial scales spanning the dimensions of individual elements, the size of stations (groupings of individual elements), and the spacing between stations. Such multi-scale, hierarchical aperture arrays require a combination of data processing architectures – pre-correlation beamformer, generic version of fast Fourier transform (FFT)-based direct imager, post-correlation beamformer, and post-correlation FFT imager – operating on different ranges of spatial scales to obtain optimal performance in imaging the entire field of view. Adopting a computational cost metric based on the number of floating point operations, its distribution over the dimensions of discovery space, namely, field of view, angular resolution, polarisation, frequency, and time is examined to determine the most efficient hybrid architectures over the parameter space of hierarchical aperture array layouts. Nominal parameters of specific upcoming and planned arrays – the SKA at low frequencies (SKA-low), SKA-low-core, a proposed long baseline extension to SKA-low (LAMBDA-I), Compact All-Sky Phased Array (CASPA), and a lunar array (FarView-core) – are used to determine the most optimal architecture hierarchy for each from a computational standpoint and provide a guide for designing hybrid architectures for multi-scale aperture arrays. For large, dense-packed layouts, a FFT-based direct imager is most efficient for most cadence intervals, and for other layouts that have relatively lesser number of elements or greater sparsity in distribution, the best architecture is more sensitive to the cadence interval, which in turn is determined by the science goals.
The flexible delivery of single-frequency lasers is far more challenging than that of conventional lasers due to the onset of stimulated Brillouin scattering (SBS). Here we present the successful delivery of 100 W single-frequency laser power through 100 m of anti-resonant hollow-core fiber (AR-HCF) in an all-fiber configuration, with the absence of SBS. By employing a custom-designed AR-HCF with a mode-field diameter matching that of a large-mode-area panda fiber, the system achieves high coupling efficiency without the need for free-space components or fiber post-processing. The AR-HCF attains a transmission efficiency of 92%, delivering an output power of 100.3 W with a beam quality factor (M2) of 1.22. The absence of SBS is confirmed through monitoring backward light, which shows no increase in intensity. This all-fiber architecture ensures high stability, compactness and efficiency, potentially expanding the application scope of single-frequency lasers in high-precision metrology, optical communication, light detection and ranging systems, gravitational wave detection and other advanced applications.
Elemental abundances, particularly the C/O ratio, are seen as a way to connect the composition of planetary atmospheres with planet formation scenario and the disc chemical environment. We model the chemical composition of gas and ices in a self-gravitating disc on timescales of 0.5 Myr since its formation to study the evolution of C/O ratio due to dust dynamics and growth and phase transitions of the volatile species. We use the thin-disc hydrodynamic code FEOSAD, which includes disc self-gravity, thermal balance, dust evolution, and turbulent diffusion, and treats dust as a dynamically different and evolving component interacting with the gas. It also describes freeze-out, sublimation, and advection of four most abundant volatile species: H$_2$O, CO$_2$, CH$_4$, and CO. We demonstrate the effect of gas and dust substructures such as spirals and rings on the distribution of volatiles and C/O ratios, including the formation of multiple snowlines of one species, and point out the anticorrelation between dust-to-gas ratio and total C/O ratio emerging due to the contribution of oxygen-rich ice mantles. We identify time and spatial locations where two distinct trigger mechanisms for planet formation are operating and differentiate them by C/O ratio range: wide range of the C/O ratios of $0-1.4$ for streaming instability, and a much narrower range $0.3-0.6$ for gravitational instability (with the initial value of 0.34). This conclusion is corroborated by observations, showing that transiting exoplanets, which possibly experienced migration through a variety of disc conditions, have significantly larger spread of C/O in comparison with directly imaged exoplanets likely formed in gravitationally unstable outer disk regions. We show that the ice-phase $\textrm{C/O}\approx$0.2–0.3 between the CO, CO$_2$, and CH$_4$ snowlines corresponds to the composition of the Solar system comets, that represent primordial planetesimals.
Understanding the physics of electromagnetic pulse (EMP) emission and nozzle damage is critical for the long-term operation of laser experiments with gas targets, particularly at facilities looking to produce stable sources of radiation at high repetition rates. We present a theoretical model of plasma formation and electrostatic charging when high-power lasers are focused inside gases. The model can be used to estimate the amplitude of gigahertz EMPs produced by the laser and the extent of damage to the gas jet nozzle. Looking at a range of laser and target properties relevant to existing high-power laser systems, we find that EMP fields of tens to hundreds of kV/m can be generated several metres from the gas jet. Model predictions are compared with measurements of EMPs, plasma formation and nozzle damage from two experiments on the VEGA-3 laser and one experiment on the Vulcan Petawatt laser.