Search

The turbulent wake behind a square-back Ahmed body in close proximity to the ground exhibits bimodal switching. This manifests as the center of the wake switching between one of two asymmetric positions, either horizontally or vertically. Switches occur over random timescales, with the wake recovering symmetry in the long time average. The present work employs wall-resolved large eddy simulations to investigate feedback control for suppressing horizontal (lateral) wake bimodality of a square-back Ahmed body at Reynolds number, Re H ∼ 3.3 × 10 4 based on the body height. Base-mounted pressure sensors are used to estimate the position of the wake as an input signal for the controller, while actuation targets the near-wake region via synthetic jets emanating from a gap around the perimeter of the Ahmed body base. A nonlinear feedback controller based on a Langevin model of the wake dynamics is synthesized. This successfully suppresses the wake lateral bimodal switching. However, this switching is replaced by a time-periodic streamwise motion of the large coherent structure occupying the near-wake region, leading to amplification of the higher frequency dynamical wake modes. The action of feedback control also leads to base pressure recovery and a reduction in pressure drag. Upon varying the controller parameters, a trade-off between the degree of bimodality suppression and drag reduction is observed. A maximum drag reduction of 7.4 % is achieved for a semisymmetrized wake, with a fully symmetrized wake achieving 2.5 % reduction. Bimodality suppression is proposed to have an indirect link to drag reduction through the effect of the wake state on the separated free shear layers and the upstream boundary layers.

High spatial resolution information on urban air pollution levels is unavailable in many areas globally, partially due to the high input data needs of existing estimation approaches. We introduced a computer vision method to estimate annual means for air pollution levels from street-level images. We used annual mean estimates of NO2 and PM2.5 concentrations from locally calibrated models as labels from London, New York, and Vancouver to allow for compilation of a sufficiently large dataset (similar to 250 k images for each city). Our experimental setup is designed to quantify intra- and intercity transferability of image-based model estimates. Performances were high and comparable to traditional land-use regression (LUR) and dispersion models when training and testing images from the same city (R-2 values between 0.51 and 0.95 when validated on data from ground monitoring stations). Similar to LUR models, transferability of models between cities in different geographies is more difficult. Specifically, transferability between the three cities (London, New York, and Vancouver), which have similar pollution source profiles, was moderately successful (R-2 values between zero and 0.67). Comparatively, performances when transferring models trained on cities with very different source profiles, such as Accra in Ghana and Hong Kong, were lower (R-2 between zero and 0.21). This suggests a need for local calibration, using additional measurement data from cities that share similar source profiles., Remote Sensing, 14 (14), ISSN:2072-4292

This paper proposes a cloud-based mangrove monitoring framework that uses Google Collaboratory and Google Earth Engine to classify mangroves in Southeast Asia (SEA) using satellite remote sensing imagery (SRSI). Three multi-class classification convolutional neural network (CNN) models were generated, showing F1-score values as high as 0.9 in only six epochs of training. Mangrove forests are tropical and subtropical environments that provide essential ecosystem services to local biota and coastal communities and are considered the most efficient vegetative carbon stock globally. Despite their importance, mangrove forest cover continues to decline worldwide, especially in SEA. Scientists have produced monitoring tools based on SRSI and CNNs to identify deforestation hotspots and drive targeted interventions. Nevertheless, although CNNs excel in distinguishing between different landcover types, their greatest limitation remains the need for significant computing power to operate. This may not always be feasible, especially in developing countries. The proposed framework is believed to provide a robust, low-cost, cloud-based, near-real-time monitoring tool that could serve governments, environmental agencies, and researchers, to help map mangroves in SEA.

Cilia and flagella are organelles that play central roles in unicellular locomotion, embryonic development, and fluid transport around tissues. In these examples, multiple cilia are often found in close proximity and exhibit coordinated motion. Inspired by the flagellar motion of biflagellate cells, we examine the synchrony exhibited by a filament pair surrounded by a viscous fluid and tethered to a rigid planar surface. A geometrically-switching base moment drives filament motion, and we characterize how the stability of synchonized states depends of the base torque magnitude. In particular, we study the emergence of bistability that occurs when the anti-phase, breast-stroke branch becomes unstable. Using a bisection algorithm, we find the unstable edge-state that exists between the two basins of attraction when the system exhibits bistability. We establish a bifurcation diagram, study the nature of the bifurcation points, and find that the observed dynamical system can be captured by a modified version of Adler’s equation. The bifurcation diagram and presence of bistability reveal a simple mechanism by which the anti-phase breast stroke can be modulated, or switched entirely to in-phase undulations through the variation of a single bifurcation parameter.

PurposeStereotactic Ablative Radiotherapy (SABR) has recently emerged as a favourable treatment option for prostate cancer patients. With higher doses delivered over fewer fractions, motion adaptation is a requirement for accurate delivery of SABR. This study compared the efficacy of multileaf collimator (MLC) tracking vs. gating as a real-time motion adaptation strategy for prostate SABR patients enrolled in a clinical trial.MethodsForty-four prostate cancer patients treated over five fractions in the TROG 15.01 SPARK trial were analysed in this study. Forty-nine fractions were treated using MLC tracking and 166 fractions were treated using beam gating and couch shifts. A time-resolved motion-encoded dose reconstruction method was used to evaluate the dose delivered using each motion adaptation strategy and compared to an estimation of what would have been delivered with no motion adaptation strategy implemented.ResultsMLC tracking and gating both delivered doses closer to the plan compared to when no motion adaptation strategy was used. Differences between MLC tracking and gating were small with differences in the mean discrepancy from the plan of -0.3% (CTV D98%), 1.4% (CTV D2%), 0.4% (PTV D95%), 0.2% (rectum V30Gy) and 0.0% (bladder V30Gy). On average, 0.5 couch shifts were required per gated fractions with a mean interruption duration of 1.8 ± 2.6 min per fraction treated using gating.ConclusionBoth MLC tracking and gating were effective strategies at improving the accuracy of the dose delivered to the target and organs at risk. While dosimetric performance was comparable, gating resulted in interruptions to treatment.Clinical trial registration numberNCT02397317.

Turbulent magnetic reconnection in a quasi-parallel shock under parameters relevant to the Earth's bow shock is investigated by means of a two-dimensional particle-in-cell simulation. The addressed aspects include the reconnection electric field, the reconnection rate, and the electron and the ion outflow speeds. In the shock transition region, many current sheets are generated in shock-driven turbulence, and electron-only reconnection and reconnection where both ions and electrons are involved can occur in those current sheets. The electron outflow speed in electron-only reconnection shows a positive correlation with the theoretical speed, which is close to the local electron Alfvén speed, and a strong convection electric field is generated by the large electron outflow. As a result, the reconnection electric field becomes much larger than those in the standard magnetopause or magnetotail reconnection. In shock-driven reconnection that involves ion dynamics, both electron outflows and ion outflows can reach of the order of 10 times the Alfvén speed in the X-line rest frame, leading to a reconnection electric field the same order as that in electron-only reconnection. An electron-only reconnection event observed by the magnetospheric multiscale mission downstream of a quasi-parallel shock is qualitatively similar to those in the simulation and shows that the outflow speed reaches approximately half the local electron Alfvén speed, supporting the simulation prediction.

An optical model is established to investigate the effects of refractive index changes on the measurement of interfacial tension by the pendant drop method with axisymmetric drop shape analysis. In such measurements, light passes from the pendant drop through a surrounding bulk phase, an optical window and air to reach the lens of the camera system. The relation between object and image size is typically determined by calibration and, if the refractive indices of any of the materials in the optical path change between calibration and measurement, a correction should be made. The simple model derived in this paper allows corrections to be calculated along with the corresponding contribution to the overall uncertainty of the interfacial tension. The model was verified by measurements of the interfacial tension between decane and water under two different calibration conditions. Neglect of the correction was shown to cause errors of up to 6 % when the bulk phase changed from air (during calibration) to water (during measurements) and of about 9 % when the system was calibrated without the optical window used for the final measurements. The refraction changes due to high pressures and supercritical fluid states can also lead to measurement errors. The proposed model can facilitate more accurate interfacial tension measurements and reduce the amount of repetitive calibration work required.

We combine extended proper orthogonal decomposition (EPOD) together with the Fukagata-Iwamoto-Kasagi (FIK) identity to quantify the role of individual coherent structures on the wall mass transfer in a turbulent pipe flow. Direct numerical simulation at a Reynolds number of 5300 (based on bulk velocity) is performed with the passive scalar released at the pipe inlet. The proper orthogonal decomposition (POD) eigenvalues show that the scalar field can be described by a more compact set of modes compared to the velocity field, and that these modes are skewed towards higher azimuthal wave numbers. POD modes for the scalar and EPOD modes for the velocity are visualized in the cross-stream plane to infer the capacity of each mode to transport scalar to and from the wall. A form of the FIK identity is derived for the wall mass transfer coefficient (Sherwood number, Sh) and employed to separate the contributions of the mean and fluctuating velocity and scalar fields. The FIK decomposition shows that the turbulent velocity/scalar correlations account for up to 65.8% of the total Sh. The contribution of each POD and EPOD mode to the Sh number is also computed; it is found that, using azimuthal wave numbers m=1–15 and POD modes n=1–10, it is possible to reconstruct 49% of the turbulent component of Sh, with the velocity modes containing only 31% of the turbulent kinetic energy. Quadrant analysis shows that these modes are related to ejection and sweep events near the wall, with the ejection events dominating.

Turbulent plasmas generate a multitude of thin current structures that can be sites for magnetic reconnection. The Magnetospheric Multiscale (MMS) mission has recently enabled the detailed examination of such turbulent current structures in Earth's magnetosheath and revealed that a novel type of reconnection, known as electron-only reconnection, can occur. In electron-only reconnection, ions do not have enough space to couple to the newly reconnected magnetic fields, suppressing ion jet formation and resulting in thinner sub-proton-scale current structures with faster super-Alfvénic electron jets. In this study, MMS observations are used to examine how the magnetic correlation length ( λC) of the turbulence, which characterizes the size of the large-scale magnetic structures and constrains the length of the current sheets formed, influences the nature of turbulence-driven reconnection. We systematically identify 256 reconnection events across 60 intervals of magnetosheath turbulence. Most events do not appear to have ion jets; however, 18 events are identified with ion jets that are at least partially coupled to the reconnected magnetic field. The current sheet thickness and electron jet speed have a weak anti-correlation, with faster electron jets at thinner current sheets. When [Formula: see text] ion inertial lengths, as is typical near the sub-solar magnetosheath, a tendency for thinner current sheets and potentially faster electron jets is present. The results are consistent with electron-only reconnection being more prevalent for turbulent plasmas with relatively short λC and may be relevant to the nonlinear dynamics and energy dissipation in turbulent plasmas.

In their seminal monograph ’Molecular Theory of Capillarity’, Rowlinson and Widom describe different possible shapes of density profiles at the vapor-liquid interface of mixtures. They postulated that in some instances, density profiles could possibly be non-monotonic, exhibiting either a maximum and/ or a minimum. This contribution revisits this statement in the light of four decades of posterior research. We summarize the distinct morphologies at the vapor-liquid interface suggested in the literature recognizing that the condition of a single minimum in the profile has not yet been reported. Interfacial density profiles with a single maximum as well as fully monotonic density profiles have been observed and reported extensively. The case of a simultaneous maximum and minimum is more controversial, as it has only been predicted using theoretical approaches such as density gradient theory (DGT). This ambiguity is further investigated in this work using the example of the vapor-liquid interface of cyclohexane + butanol. Both DGT in combination with several distinct equations of state and molecular dynamics simulations are used. The results from the two methods are found to be contradictory: while the DGT results predict a maximum/minimum structure, the computer experiment results indicate only a single maximum in the density profiles. This work thereby emphasizes that results from DGT for highly non-ideal mixtures should not be taken for granted.

The turbulent axisymmetric bluff body wake is studied with a large eddy simulation (LES). The effect of pulsed jet forcing on the wake is investigated. For completeness, the pulsed jet actuator is included in the computational domain. Spectral proper orthogonal decomposition (SPOD) is applied to analyze the wake and the effect of pulsed jet forcing. The numerical results show good agreement with the experimental results and successfully reproduce the dominant modes of the axisymmetric wake. SPOD identifies the axisymmetric bubble pumping mode ( St D = 0.06 ) , axisymmetric-breaking vortex shedding mode ( St D = 0.22 ) , and an axisymmetric-breaking mode close to the subharmonic of the shedding mode ( St D = 0.1 ) . When the high-frequency pulsed jet ( St θ = 0.107 ) is applied to the wake, the pressure fluctuation on the base and the azimuthal modes are globally weakened without any mode selection. The mean flow properties show that the inlet flow is driven towards the wake centerline by the high-frequency jet, and the shear layer deviates, forming a concave separation streamline. The entrainment in the wake is suppressed, and the wake is narrowed, leading to a global pressure rise. The high-frequency pulsed jet is concentric and greatly enhances the entrainment in the vicinity of the separation point. A low-pressure region is generated there, which drives the inlet flow towards the wake. Low-frequency forcing ( St θ = 0.029 ) generates a diffusive pulsed jet which enhances the mixing and entrainment in its trajectory. Low-frequency forcing reduces the wake length and accelerates recirculating flow near the base, leading to a reduction in base pressure, and the vortex shedding mode ( St D = 0.2 ) is amplified.

At the Earth's low-latitude magnetopause, the Kelvin–Helmholtz instability (KHI), driven by the velocity shear between the magnetosheath and magnetosphere, has been frequently observed during northward interplanetary magnetic field (IMF) periods. However, the signatures of the KHI have been much less frequently observed during southward IMF periods, and how the KHI develops under southward IMF has been less explored. Here, we performed a series of realistic 2D and 3D fully kinetic simulations of a KH wave event observed by the Magnetospheric Multiscale (MMS) mission at the dusk-flank magnetopause during southward IMF on September 23, 2017. The simulations demonstrate that the primary KHI bends the magnetopause current layer and excites the Rayleigh–Taylor instability (RTI), leading to penetration of high-density arms into the magnetospheric side. This arm penetration disturbs the structures of the vortex layer and produces intermittent and irregular variations of the surface waves which significantly reduces the observational probability of the periodic KH waves. The simulations further demonstrate that in the non-linear growth phase of the primary KHI, the lower-hybrid drift instability (LHDI) is induced near the edge of the primary vortices and contributes to an efficient plasma mixing across the magnetopause. The signatures of the large-scale surface waves by the KHI/RTI and the small-scale fluctuations by the LHDI are reasonably consistent with the MMS observations. These results indicate that the multi-scale evolution of the magnetopause KH waves and the resulting plasma transport and mixing as seen in the simulations may occur during southward IMF

We study properties of the turbulent/nonturbulent interface (TNTI) between two layers of stratified fluids through direct numerical simulations (DNSs). Zero mean shear forcing creates moderate turbulence in one of the layers with the Taylor microscale Reynolds numbers in the mixed region of Re λ = 35 , 44 . We focus on the similarities and differences of the properties of stratified TNTIs due to two distinct types of forcing: (a) convection due to a boundary heat source and (b) agitation resembling a vertically oscillating grid experiment. Similarly to other stratified flows, the small scale dynamics of the TNTI in the present DNSs differ from what would be expected in comparable yet unstratified TNTIs. The interface cannot be indeed uniquely identified by the commonly used vorticity ω . Instead, the potential enstrophy Π 2 is shown to be the most appropriate marker in these flow cases. It is emphasized that the Kolmogorov lengthscale η K ∼ √ ν / ω is not representative of the small scale dynamics of the interface. Hence, an alternative lengthscale, η Π , is defined, in analogy to the Kolmogorov scale, based on the potential enstrophy, η Π = ( ν 3 / Π ∗ ) 1 / 6 , being Π ∗ = | g / ρ 0 Π | . The conditionally averaged profiles of potential enstrophy Π 2 , enstrophy ω 2 , and turbulent kinetic energy dissipation ε of the two distinctly different turbulence forcing cases collapsed when scaled by η Π at different time instants in each simulation. This implies not only the self-similarity of the small scale statistics of the TNTI in either of the two cases, but also the similarity between the statistics of the two different turbulent flows in the proximity of TNTI.

Introducing an axis of reflectional symmetry in a quasicrystal leads to the creation of localised edge modes that can be used to build waveguides. We develop theory that characterises reflection-induced localised modes in materials that are formed by recursive tiling rules. This general theory treats a one-dimensional continuous differential model and describes a broad class of both quasicrystalline and periodic materials. We present an analysis of a material based on the Fibonacci sequence, which has previously been shown to have exotic, Cantor-like spectra with very wide spectral gaps. Our approach provides a way to create localised edge modes at frequencies within these spectral gaps, giving strong and stable wave localisation. We also use our general framework to make a comparison with reflection-induced modes in periodic materials. These comparisons show that while quasicrystalline waveguides enjoy enhanced robustness over periodic materials in certain settings, the benefits are less clear if the decay rates are matched. This shows the need to carefully consider equivalent structures when making robustness comparisons and to draw conclusions on a case-by-case basis, depending on the specific application.

Analysis of ground movement rates along the coastline and upper sections of the Ventnor landslide complex was carried out utilizing Persistent Scatterer Interferometric Synthetic Aperture Radar methods using Sentinel-1 SAR data from 2015 to 2019 (four years). Results were compared with rainfall data, historical ground investigation records and monitoring surveys carried out at Ventnor to relate observations to geology, geomorphology and rainfall. Decomposition of InSAR viewing geometries to vertical and horizontal aligned well with previous ground-based studies. Subsidence of −9.8 mm a−1 at the Lowtherville Graben and heave of +8.5 mm a−1 along the coastline south of Ventnor Park were observed. Decomposition to east-west geometry results showed an eastward displacement of approximately 12.4 mm a−1 along the coastline south of Ventnor Park, and a westward displacement of −3.7 mm a−1 throughout built up sections of Ventnor town, indicating the landslide was displacing more in an eastern direction than vertically. The cause of this movement was investigated by using publicly available intrusive boreholes paired with Persistent Scatterer Interferometry, and a new ground model spanning east-west parallel to the coastline was presented. No evidence of significant ground movement was observed along heavily protected sections of the coastline, suggesting coastal defences comprised of concrete aprons and rip rap appear to be an effective coastal management/landslide stabilisation tool when compared to rip rap alone. The mechanism of this increased stability is likely due to the combination of toe weighting and reduced toe erosion. A lag of approximately 13–20 days was observed between high rainfall events and subsequent peaks in ground displacement, which was shorter than a 29 day lag observed in a previous study. Similar observations of prolonged rainfall resulting in prolonged displacements were also observed. The study demonstrates the capabilities of the PSI methodology in identifying the same ground movements that conventional methods provide. By providing detailed analysis of ground deformation of the Ventnor landslide, we demonstrate small ground movements, validated with existing ground movement surveys. Similar methodology can be applied to coastal landslides in urban environments worldwide, providing a relatively cheap and rapid resource for coastal landslide monitoring.

We report the speeds of sound in liquid toluene (methylbenzene) measured using double-path pulse-echo apparatus independently at The University of Western Australia (UWA) and Imperial College London (ICL). The UWA data were measured at temperatures between (306 and 423) K and at pressures up to 65 MPa with standard uncertainties of between (0.02 and 0.04)%. At ICL, measurements were made at temperatures between (283.15 and 473.15) K and at pressures up to 390 MPa with standard uncertainty of 0.06%. By means of thermodynamic integration, the measured sound-speed data were combined with initial density and isobaric heat capacity values obtained from extrapolated experimental data to derive a comprehensive set of thermodynamic properties of liquid toluene over the full measurement range. Extensive uncertainty analysis was performed by studying the response of derived properties to constant and dynamic perturbations of the sound-speed surface, as well as the initial density and heat capacity values. The relative expanded uncertainties at 95% confidence of derived density, isobaric heat capacity, isobaric expansivity, isochoric heat capacity, isothermal compressibility, isentropic compressibility, thermal pressure coefficient and internal pressure were estimated to be (0.2, 2.2, 1.0, 2.6, 0.6, 0.2, 1.0 and 2.7)%, respectively. Due to their low uncertainty, these data and derived properties should be well suited for developing a new and improved fundamental Helmholtz equation of state for toluene.

Counter-flow configurations, whereby two streams of fluid are brought together from opposite directions, are highly efficient mixers due to the high turbulence intensities that can be maintained. In this paper, a simplified version of the problem is introduced that is amenable to direct numerical simulation. The resulting turbulent flow problem is confined between two walls, with one non-zero mean velocity component varying in the space direction normal to the wall, corresponding to a simple shear flow. Compared to conventional channel flows, the mean flow is inflectional and the maximum turbulence intensity relative to the maximum mean velocity is nearly an order of magnitude higher. The numerical requirements and turbulence properties of this configuration are first determined. The Reynolds shear stress is required to vary linearly by the imposed forcing, with a peak at the channel centreline. A similar behaviour is observed for the streamwise Reynolds stress, the budget of which shows an approximately uniform distribution of dissipation, with large contributions from production, pressure-strain and turbulent diffusion. A viscous sublayer is obtained near the walls and with increasing Reynolds number small-scale streaks in the streamwise momentum are observed, superimposed on the large-scale structures that buffet this region. When the peak local mean Mach number reaches 0.55, turbulent Mach numbers of 0.6 are obtained, indicating that this flow configuration can be useful to study compressibility effects on turbulence.

The Rayleigh-Taylor (RT) instability occurs at the interface between two fluids of different densities when the lighter fluid pushes against the heavier. For example, RT instability fingers are evident in the Crab Nebula. It plays a critical role in all known forms of fusion. The variable density model (VDM), which is comprised of a set of partial differential equations, best models the flow in the mixing layer. This paper shows that the VDM can be exactly transformed into the inhomogeneous, incompressible, forced Navier-Stokes equations.

Biermann battery magnetic field generation driven by high power laser-solid interactions is explored in experiments performed with the OMEGA EP laser system. Proton deflectometry captures changes to the strength, spatial profile, and temporal dynamics of the self-generated magnetic fields as the target material or laser intensity is varied. Measurements of the magnetic flux during the interaction are used to help validate extended magnetohydrodynamic (MHD) simulations. Results suggest that kinetic effects cause suppression of the Biermann battery mechanism in laser-plasma interactions relevant to both direct and indirect-drive inertial confinement fusion. Experiments also find that more magnetic flux is generated as the target atomic number is increased, which is counter to a standard MHD understanding., 9 pages, 5 figures

In this work, we present a model, based on rough hard-sphere theory, for the tracer diffusion coefficients of gaseous solutes in non-polar liquids. This work extends an earlier model developed specifically for carbon dioxide in hydrocarbon liquids and establishes a general correlation for gaseous solutes in non-polar liquids. The solutes considered were light hydrocarbons, carbon dioxide, nitrogen and argon, while the solvents were all hydrocarbon liquids. Application of the model requires knowledge of the temperature-dependent molar core volumes of the solute and solvent, which can be determined from pure-component viscosity data, and a temperature-independent roughness factor which can be determined from a single diffusion coefficient measurement in the system of interest. The new model was found to correlate the experimental data with an average absolute relative deviation of 2.7 %. The model also successfully represents computer-simulation data for tracer diffusion coefficients of hard-sphere mixtures and reduces to the expected form for self-diffusion when the solute and solvent become identical.

This work is framed within AIChE's 10th Industrial Fluid Properties Simulation Challenge, with the aim of assessing the capability of molecular simulation methods and force fields to accurately predict the pressure dependence of the shear viscosity of 2,2,4-trimethylhexane at 293.15 K (20 °C) at pressures up to 1 GPa. In our entry for the challenge, we employ coarse-grained intermolecular models parametrized via a top-down technique where an accurate equation of state is used to link the experimentally-observed macroscopic volumetric properties of fluids to the force-field parameters. The state-of-the-art version of the statistical associating fluid theory (SAFT) for potentials of variable range as reformulated in the Mie incarnation is employed here. The potentials are used as predicted by the theory, with no fitting to viscosity data. Viscosities are calculated by molecular dynamics (MD) employing two independent methods; an equilibrium-based procedure based on the analysis of the pressure fluctuations through a Green-Kubo formulation and a non-equilibrium method where periodic perturbations of the boundary conditions are employed to simulate experimental shear stress conditions. There is an indication that, at higher pressures, the model predicts a solid phase (freezing) which we believe to be an artefact of the simplified molecular geometry used in the modelling. A comparison (made after disclosure of the experimental data) show that the model consistently underpredicts the viscosity by about 30%, but follows the pressure dependency accurately.

1-component Extended Hard-Sphere (1-cEHS) model has been developed recently to predict the viscosity of liquid, n-alkane mixtures. It represents a mixture by a single pseudo-component characterized by an appropriate molecular weight and calculates the viscosity by means of the modified, extended hard-sphere model (EHS) that makes use of a universal function relating reduced viscosity to reduced volume. In this work we have extended the model to also predict the viscosity of mixtures containing alkylbenzene and cycloalkane species. Furthermore, we have developed a new 3-component Extended Hard-Sphere (3-cEHS) model which requires only a knowledge of the overall composition of n-alkane, alkylbenzene and cycloalkane species. Extensive comparison with the available experimental data indicates that both models (1-cEHS and 3-cEHS) predict the viscosity of binary and multicomponent mixtures containing n-alkane, alkylbenzene and cycloalkane species with uncertainty of 5–10%. The proposed models are a precursor of a new family of models that do not require a knowledge of the detailed composition of the mixture, but still take advantage of the underlying molecular description.

The paper evaluates the interaction between atomization quality and atomization efficiency, which has not been understood, although it is commonly observed that the atomization quality of different atomizers does not improve linearly with addition of energy. The results quantify the energy exchange between the air and liquid streams of a twin water impinging jets atomizer and its consequences on atomization characteristics and explain the behavior of quality and efficiency. The liquid jet breakup length, liquid jets separation distance at the breakup region and spray angles were measured with high speed photography and the droplet characteristics, such as spatial distributions of mean droplet velocities and diameters and normalized liquid volume flux with Phase Doppler Anemometry (PDA). The results show that the breakup length decreased and the separation distance of the interacting liquid jets at the geometrical ‘impingement’ region increased rapidly as Air-to-Liquid Momentum Ratio (ALMR) increased and then remained constant for ALMR>9. Spray angles were different on different planes through the spray and generally decreased with increasing ALMR and were insensitive to the liquid jets impingement angle. The spatial distributions of average droplet size, velocity and normalized liquid volume flux in the sprays became elongated normal to the plane of the two liquid jets for larger liquid flow rates, in agreement with the spray angle in the near nozzle region. The spatially-averaged Sauter Mean Diameter (SMD) of the sprays quantified uniquely the atomization quality and showed, for the first time, that it did not depend on liquid jets impingement angle. The average SMD

Parameter identification for wildfire forecasting models often relies on case-by-case tuning or posterior diagnosis/analysis, which can be computationally expensive due to the complexity of the forward prediction model. In this paper, we introduce an efficient parameter flexible fire prediction algorithm based on machine learning and reduced order modelling techniques. Using a training dataset generated by physics-based fire simulations, the method forecasts burned area at different time steps with a low computational cost. We then address the bottleneck of efficient parameter estimation by developing a novel inverse approach relying on data assimilation techniques (latent assimilation) in the reduced order space. The forward and the inverse modellings are tested on two recent large wildfire events in California. Satellite observations are used to validate the forward prediction approach and identify the model parameters. By combining these forward and inverse approaches, the system manages to integrate real-time observations for parameter adjustment, leading to more accurate future predictions.

An analytical solution is derived for the flow generated by a self-propelling two-dimensional Marangoni boat driven by reactive insoluble surfactant on a deep layer of fluid of viscosity μ at zero Reynolds number, capillary number, and surface Péclet number. In the model, surfactant emitted from the edges of the boat causes a surface tension disparity across the boat. Once emitted, the surfactant diffuses along the interface and sublimates to the upper gas phase. A linear equation of state relates the surface tension to the surfactant concentration. The propulsion speed of the boat is shown to be U0=Δσ(2πμ)−1e√DaK0(√Da) where Da is a Damköhler number measuring the reaction rate of the surfactant to its surface diffusion, Δσ is the surface tension disparity between the front and rear of the boat, and K0 is the order-zero modified Bessel function. Explicit expressions for the stream function associated with the Stokes flow beneath the boat are found facilitating ready examination of the Marangoni-induced streamlines. An integral formula, derived using the reciprocal theorem, is also given for the propulsion speed of the boat in response to a more general Marangoni stress distribution.

In this paper, the turbulence fields generated by a group of modified fractal grids, referred to as the rectangular fractal grids (RFGs), are documented and discussed. The experiments were carried out using hot-wire anemometry in three facilities at Imperial College London and the Laboratory of Fluid Mechanics in Lille, France. Due to the large Reynolds number of the resulting turbulence, several data processing methods for turbulence properties are carefully evaluated. Two spectral models were adopted, respectively, to correct the large and small wave-number ranges of the measured spectrum. After the technical discussion, the measurement results are presented in terms of one-point statistics, length scales, homogeneity, isotropy, and dissipation. The main conclusions are twofold. First, the location of maximum turbulence intensity xpeak is shown to be independent of the inlet Reynolds number but dependent on the ratio between the lengths of the largest grid bars in the transverse and vertical directions. This is crucial to the production of prescribed features of turbulent flows in laboratory. Second, these RFG-generated turbulent flows are shown to be quasihomogeneous in the decay region for x/xpeak>1.5, but the isotropy is poorer than that of the previous studied fractal square grid-generated turbulence. In the beginning of the decay region, a decreasing pattern of the integral length scale Lu and Taylor microscale λ was observed, yet the ratio Lu/λ remained roughly constant along the centerline, so Cε∼Re−1λ, complying with the nonequilibrium scaling relation reported in previous studies for various turbulent flows.

Run-and-tumble (RnT) motion is an example of active motility where particles move at constant speed and change direction at random times. In this work we study RnT motion with diffusion in a harmonic potential in one dimension via a path integral approach. We derive a Doi-Peliti field theory and use it to calculate the entropy production and other observables in closed form. All our results are exact.

To investigate the influence of the orifice geometry on near-field coherent structures in a jet, Fourier proper orthogonal decomposition (Fourier-POD) is applied. Velocity and vorticity snapshots obtained from tomographic particle image velocimetry at the downstream distance of two equivalent orifice diameters are analyzed. Jets issuing from a circular orifice and from a fractal orifice are examined, where the fractal geometry is obtained from a repeating fractal pattern applied to a base square shape. While in the round jet energy is mostly contained at wave number m=0, associated to the characteristic Kelvin-Helmholtz vortex rings, in the fractal jet modal structures at the fundamental azimuthal wave number m=4 capture the largest amount of energy. In addition, energy is scattered across a wider range of wave numbers than in the round jet. The radial Fourier-POD profiles, however, are nearly insensitive to the orifice geometry, and collapse to a universal distribution when scaled with a characteristic radial length. A similar collapse was recently observed in POD analysis of turbulent structures in pipe flow. However, unlike in pipe flow, the azimuthal-to-radial aspect ratio of the Fourier-POD structures is not constant and varies greatly with the wave number. The second part of the paper focuses on the relationship between streamwise vorticity and streamwise velocity, to characterize the role of the orifice geometry on the lift-up mechanism recently found to be active in turbulent jets [P. Nogueira, A. Cavalieri, P. Jordan, and V. Jaunet, Large-scale streaky structures in turbulent jets, J. Fluid Mech. 873, 211 (2019)JFLSA70022-112010.1017/jfm.2019.365]. The averaging of the streamwise vorticity conditioned on intense positive fluctuations of streamwise velocity reveals a pair of vorticity structures of opposite sign flanking the conditioning point, inducing a radial flow towards the jet periphery. This pair of structures is observed in both jets, even if the azimuthal extent of this pattern is 30% larger in the jet issuing from the circular orifice. The coupling between streamwise vorticity and velocity motions is also examined using Fourier-POD. The analysis reveals that in the jet with a circular orifice lower wave-number modes, corresponding to structures at larger scales, capture a larger fraction of the vorticity-velocity coupling. This evidences that the orifice geometry directly influences the interaction between velocity and vorticity.

The characteristics of the internal layers of intense shear are examined in a mixing layer and in a jet, in the range of Reynolds numbers 134

A minimal form of quasilinear approximation (QLA), recently proposed with a stochastic forcing and proper orthogonal decomposition modes [Hwang and Ekchardt, J. Fluid Mech. 894, A23 (2020)], has been extended by employing a resolvent framework. A particular effort is made to reach an extremely high Reynolds number by carefully controlling the approximation without loss of the general scaling properties in the spectra, while setting out the main limitations and accuracy of the proposed QLA with possibility of further improvement. The QLA is subsequently applied to turbulent channel flow up to Reτ=106 (Reτ is the friction Reynolds number). While confirming that the logarithmic wall-normal dependence in streamwise and spanwise turbulence intensities robustly appears, it reveals some nontrivial difference from the scaling of the classical attached eddy model based on inviscid flow assumption. First, the spanwise wave number spectra do not show any clearly visible inverse-law behavior due to the viscous wall effect prevailing in a significant portion of the lower part of the logarithmic layer. Second, the near-wall peak streamwise and spanwise turbulence intensities are found to deviate from lnReτ scaling for Reτ≳104. Importantly, the near-wall streamwise turbulence intensity is inversely proportional to 1/U+cl (U+cl is the inner-scaled channel centreline velocity), consistent with the scaling obtained from an asymptotic analysis of the Navier-Stokes equations [Monkewitz and Nagib, J. Fluid Mech. 783, 474 (2015)]. The same behavior was also observed for the streamwise turbulence intensity in the logarithmic region, as was predicted with the asymptotic analysis. Finally, the streamwise turbulence intensity in the logarithmic region is found to become greater than the near-wall one at Reτ≃O(105). It is shown that this behavior originates from the near-wall spectra associated with large-scale inactive motions, the intensity of which gradually decays as Reτ→∞.

The Orr mechanism is revisited to understand its precise role in the transition of plane Couette flow. By considering homogeneous shear flow and plane Couette flow, it is identified that the Orr mechanism induces a lift-up effect which significantly amplifies spanwise velocity. An optimal perturbation analysis for an individual velocity component reveals that the amplification of spanwise velocity is most active at the streamwise length comparable to the given spanwise length of the perturbation. The relevance of this mechanism to transition is subsequently examined in plane Couette flow. To this end, a set of initial conditions, which combines the optimal perturbation for spanwise velocity with the one for all the velocity components, is considered by varying their amplitudes. Two representative transition scenarios are found: oblique and streak transitions. In the former, the spanwise velocity perturbation amplified with the Orr mechanism initiates both streak amplification and breakdown, whereas in the latter, its role is limited only to the streak breakdown at the late stage of transition. As such, the oblique transition offers a route to turbulence energetically more efficient than the streak transition, at least for the cases examined in the present paper. Finally, the oblique transition is found to share many physical similarities with the transition by the minimal seed.

MIF approaches, such as the MagLIF experiment, use magnetic fields in dense plasma to suppress cross-field thermal conduction, attempting to reduce heat losses and trap alpha particles to achieve ignition. However, the magnetic field can introduce other transport effects, some of which are deleterious. An understanding of these processes is thus crucial for accurate modelling of MIF. We generalise past work exploiting self-similar solutions to describe transport processes in planar geometry and compare the model to the radiation-magnetohydrodynamics code Chimera. We solve the 1D extended MHD equations under pressure balance, making no assumptions about the ratio of magnetic and thermal pressures in the plasma. The resulting ODE boundary value problem is solved using a shooting method, combining an implicit ODE solver and a Newton-Raphson root finder. We show that the Nernst effect dominates over resistive diffusion in high $\beta$ plasma, but its significance is reduced as the $\beta$ decreases. On the other hand, we find that Ettingshausen and Ohmic heating effects are dominant in low $\beta$ plasma, and can be observable in even order unity $\beta$ plasma, though in the presence of a strong temperature gradient heat conduction remains dominant. We then present a test problem for the Ohmic heating and Ettingshausen effects which will be useful to validate codes modelling these effects. We also observe that the Ettingshausen effect plays a role in preventing temperature separation when Ohmic heating is strong. Neglecting this term may lead to overestimates for the electron temperature at a vacuum-plasma interface, such as at the edge of a z-pinch. The model developed can be used to provide test problems with arbitrary boundary conditions for magnetohydrodynamics codes, with the ability to freely switch on terms to compare their individual implementations., Comment: 32 pages, 7 figures

We report Magnetospheric Multiscale four-spacecraft observations of a thin reconnecting current sheet with weakly asymmetric inflow conditions and a guide field of approximately twice the reconnecting magnetic field. The event was observed at the interface of interlinked magnetic field lines at the flank magnetopause when the maximum spacecraft separation was 370 km and the spacecraft covered & SIM;1.7 ion inertial lengths (d(i)) in the reconnection outflow direction. The ion-scale spacecraft separation made it possible to observe the transition from electron-only super ion-Alfvenic outflow near the electron diffusion region (EDR) to the emergence of sub-Alfvenic ion outflow in the ion diffusion region (IDR). The EDR to IDR evolution over a distance less than 2 d(i) also shows the transition from a near-linear reconnecting magnetic field reversal to a more bifurcated current sheet as well as significant decreases in the parallel electric field and dissipation. Both the ion and electron heating in this diffusion region event were similar to the previously reported heating in the far downstream exhausts. The dimensionless reconnection rate, obtained four different ways, was in the range of 0.13-0.27. This event reveals the rapid spatial evolution of the plasma and electromagnetic fields through the EDR to IDR transition region.& nbsp;(C) 2021 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Industrial inspection protocols are qualified using mock-ups manufactured according to the same procedure as the plant part. For coarse-grained castings, known for their low inspectability, relying on mock-ups becomes particularly challenging owing to the variability of grain properties among components. Consequently, there is a keen interest in the capability to verify whether the grain size of the component under test matches the qualification specification in-situ. This paper investigates the potential of an attenuation measurement for assessing the ultrasonic inspectability of coarse-grained components using qualified procedures in a practical setting. The experimental part of the study focuses on an industrial Inconel 600 mock-up with spatially varying attenuation, measured across the entire sample in an immersion tank. Three zones with distinctly different attenuations were examined using metallography, which allowed for calculating classical grain size histograms and two-point correlation functions. For one of the zones, we synthesised the microstructure with the same statistical properties numerically and simulated the propagation of ultrasound using a grain-scale finite element model. The results showed good agreement with the experiment, and lead to several suggestions for the reasons for the discrepancy, the varying grain size statistics being the most likely. A parametric study, which followed, depicted the effect of the mean and standard deviation-to-mean ratio of the log-normal grain size distribution on the attenuation of ultrasound and its frequency dependence. Most notably, we demonstrated the known non-uniqueness of the relationship between the log-normal grain size distribution parameters and the attenuation. We suggested that the correlation length calculated from a single exponential fit to the two-point correlation function is a more robust metric describing grain statistics for this context, which can be obtained from attenuation. The correlation lengths estimated from measured attenuation using the second-order approximation model for the three zones of the studied mock-up yielded results of acceptable accuracy. We concluded that this metric could replace the average grain size in practical settings, as it retains more statistical information than the mean grain size and allows for linking measurements to the established theoretical attenuation models which this paper demonstrates.

Paper-based devices have a wide range of applications in point-of-care diagnostics, environmental analysis, and food monitoring. Paper-based devices can be deployed to resource-limited countries and remote settings in developed countries. Paper-based point-of-care devices can provide access to diagnostic assays without significant user training to perform the tests accurately and timely. The market penetration of paper-based assays requires decreased device fabrication costs, including larger packing density of assays (i.e., closely packed features) and minimization of assay reagents. In this review, we discuss fabrication methods that allow for increasing packing density and generating closely packed features in paper-based devices. To ensure that the paper-based device is low-cost, advanced fabrication methods have been developed for the mass production of closely packed assays. These emerging methods will enable minimizing the volume of required samples (e.g., liquid biopsies) and reagents in paper-based microfluidic devices.

We present an evaluation of Aerosol Optical Depth (AOD) retrievals from the Moderate Resolution Imaging Spectroradiometer (MODIS) over Australia covering the period 2001–2020. We focus on retrievals from the Deep Blue (DB) and Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithms, showing how these compare to one another in time and space. We further employ speciated AOD estimates from Copernicus Atmospheric Monitoring Service (CAMS) reanalyses to help diagnose aerosol types and hence sources. Considering Australia as a whole, monthly mean AODs show similar temporal behaviour, with a well-defined seasonal peak in the Austral summer. However, excepting periods of intense biomass burning activity, MAIAC values are systematically higher than their DB counterparts by, on average, 50%. Decomposing into seasonal maps, the patterns of behaviour show distinct differences, with DB showing a larger dynamic range in AOD, with markedly higher AODs (ΔAOD∼0.1) in northern and southeastern regions during Austral winter and summer. This is counter-balanced by typically smaller DB values across the Australian interior. Site level comparisons with all available level 2 AOD data from Australian Aerosol Robotic Network (AERONET) sites operational during the study period show that MAIAC tends to marginally outperform DB in terms of correlation (RMAIAC = 0.71, RDB = 0.65) and root-mean-square error (RMSEMAIAC = 0.065, RMSEDB = 0.072). To probe this behaviour further, we classify the sites according to the predominant surface type within a 25 km radius. This analysis shows that MAIAC’s advantage is retained across all surface types for R and all but one for RMSE. For this surface type (Bare, comprising just 1.2% of Australia) the performance of both algorithms is relatively poor, (RMAIAC = 0.403, RDB = 0.332).

Interfacial tension (IFT) between reservoir fluids is an important property in enhanced oil recovery (EOR) and carbon geological storage (CGS). Quantitative knowledge of IFT is needed to support and assist the interpretation of multiphase flow and wetting behaviour in porous media and to facilitate numerical reservoir simulation. Iododecane and iodide-containing brines are common contrast agents in visualisation of multiphase flow in porous media by X-ray CT imaging. The effect of the introduced contrast agents on the IFT was studied in this work by means of pendant-drop experiments and modelling with the density-gradient theory. We report experimental IFTs between N2, 7 mass% KI (aq), and decane-iododecane mixtures with various iododecane mass fractions at temperatures from 298 K to 353 K and pressures from 1 MPa to 30 MPa. The IFTs between N2 and the liquid phases decrease with the increase of either pressure or temperature and increase with the increasing KI molality or iododecane mass fraction. The IFTs between H2O and decane-iododecane mixtures decrease with temperature or iododecane mass fraction and increase slightly with increasing pressure. The IFT data were modelled by means of the density-gradient theory coupled with the volume-translated Peng-Robinson equation of state. Empirical equations were also developed to correlate all of the measured data. A workflow was proposed for estimating the IFTs between gas, brine and the doped hydrocarbon systems based on the experimental and modelling work.