Search

We study a two-dimensional isotropic rotating system and obtain both theoretically and numerically a $K^{-2}$ energy spectrum under the rapidly rotating condition ($Ro\ll 1$), which was initially obtained by Zeman (1994) and Zhou (1995). In rotating turbulence, the $K^{-2}$ energy spectrum was proposed under the assumption of isotropy, however, the direction selectivity of rotation breaks isotropy, making this $K^{-2}$ spectrum not easily observable. To fill the gap between theoretical assumptions and realizability, we study the turbulence of inertial waves in an artificial two-dimensional isotropic rotating turbulence system. In the limit of a small Rossby number, we asymptotically derive a nonlinear amplitude equation for inertial waves, which gives the $K^{-2}$ spectrum using a strong turbulence argument. This scaling is justified by numerical simulations of both the amplitude equation and the original system.

With the rapid growth of Large Language Models (LLMs) across various domains, numerous new LLMs have emerged, each possessing domain-specific expertise. This proliferation has highlighted the need for quick, high-quality, and cost-effective LLM query response methods. Yet, no single LLM exists to efficiently balance this trilemma. Some models are powerful but extremely costly, while others are fast and inexpensive but qualitatively inferior. To address this challenge, we present TO-Router, a non-monolithic LLM querying system that seamlessly integrates various LLM experts into a single query interface and dynamically routes incoming queries to the most high-performant expert based on query's requirements. Through extensive experiments, we demonstrate that when compared to standalone expert models, TO-Router improves query efficiency by up to 40\%, and leads to significant cost reductions of up to 30%, while maintaining or enhancing model performance by up to 10%., Comment: 14 pages, 7 figures, 2 tables

Large language models (LLMs) have surged in popularity and are extensively used in commercial applications, where the efficiency of model serving is crucial for the user experience. Most current research focuses on optimizing individual sub-procedures, e.g. local inference and communication, however, there is no comprehensive framework that provides a holistic system view for optimizing LLM serving in an end-to-end manner. In this work, we conduct a detailed analysis to identify major bottlenecks that impact end-to-end latency in LLM serving systems. Our analysis reveals that a comprehensive LLM serving endpoint must address a series of efficiency bottlenecks that extend beyond LLM inference. We then propose ScaleLLM, an optimized system for resource-efficient LLM serving. Our extensive experiments reveal that with 64 concurrent requests, ScaleLLM achieves a 4.3x speed up over vLLM and outperforms state-of-the-arts with 1.5x higher throughput.

In this paper, we systematically investigate the general dark matter-electron interactions within the framework of effective field theories (EFT). We consider both the non-relativistic (NR) EFT and the relativistic EFT descriptions of the interactions with the spin of dark matter (DM) up to one, i.e., the scalar ($\phi$), fermion ($\chi$), and vector $(X)$ DM scenarios. We first collect the leading-order NR EFT operators describing the DM-electron interactions, and construct especially the NR operators for the vector DM case. Next, we consider all possible leading-order relativistic EFT operators including those with a photon field and perform the NR reduction to match them onto the NR EFT. Then we rederive the DM-bound-electron scattering rate within the NR EFT framework and find that the matrix element squared, which is the key input that encodes the DM and atomic information, can be compactly decomposed into three terms. Each term is a product of a DM response function $(a_{0,1,2})$, which is essentially a factor of Wilson coefficients squared, and its corresponding generalized atomic response function ($\widetilde W_{0,1,2}$). Lastly, we employ the electron recoil data from the DM direct detection experiments (including XENON10, XENON1T, and PandaX-4T) to constrain all the non-relativistic and relativistic operators in all three DM scenarios. We set strong bounds on the DM-electron interactions in the sub-GeV region. Particularly, we find that the latest PandaX-4T S2-only data provide stringent constraints on dark matter with a mass greater than approximately 20 MeV, surpassing those from the previous XENON10 and XENON1T experiments., Comment: 37 pages, 8 figures, minor corrections for several typos, and added one reference. Accepted by JHEP

In this letter we revisit general dark matter (DM)-bound-electron interactions studied previously in the influential work of [Catena et al., Phys. Rev. Res. 2, 033195 (2020)]. We derive the DM-electron response functions and find a crucial minus sign was missed for the second atomic response function $W_2$ defined in that work. The minus sign has significant phenomenological consequences when explaining experimental bounds on specific DM scenarios. Furthermore, for the most general DM-electron nonrelativistic or relativistic interactions for DM with spin up to one, we find there are three DM response functions ($a_{0,1,2}$) whose corresponding atomic response functions ($\widetilde W_{0,1,2}$) are linear combinations of the four response functions ($W_{1,2,3,4}$) given in that work, $$ \widetilde W_0 = W_1, \, \widetilde W_2 = |\mathbf{v}_0^\perp|^2 W_1+ W_3 - 2 {m_e\, \mathbf{q}\cdot \mathbf{v}_0^\perp \over \mathbf{q}^2} W_2,\, \widetilde W_3 = { (\mathbf{q}\cdot \mathbf{v}_0^\perp)^2 \over \mathbf{q}^2} W_1 + {m_e^2 \over \mathbf{q}^2}W_4 - 2 {m_e\, \mathbf{q}\cdot \mathbf{v}_0^\perp \over \mathbf{q}^2} W_2. $$ Due to the minus sign correction for $W_2$, there can be significant cancellations between the $W_2$ and $W_{3,4}$ terms, so that $\widetilde W_{2,3}$ are dominated by the usual response function $W_1$ in some cases. Ignoring the sign could thus result in misinterpretation of the experimental data in some DM scenarios. As an example, we show that the recent XENON1T constraint on the fermionic DM anapole moment is weakened by a factor of 2 or so. Many DM scenarios involving DM or electron axial-vector current can yield $W_2$ and thus are potentially affected by the sign., Comment: 6+2 pages, 4 figures

Helical magnetohydrodynamic turbulence with Hall effects is ubiquitous in heliophysics and plasma physics, such as star formation and solar activities, and its intrinsic mechanisms are still not clearly explained. Direct numerical simulations reveal that when the forcing scale is comparable to the ion inertial scale, Hall effects induce remarkable cross helicity. It then suppresses the inverse cascade efficiency, leading to the accumulation of large-scale magnetic energy and helicity. The process is accompanied by the breaking of current sheets via filaments along magnetic fields. Using the Ulysses data, the numerical findings are separately confirmed. These results suggest a novel mechanism wherein small-scale Hall effects could strongly affect large-scale magnetic fields through cross helicity., Comment: 5 figures, 6 pages

Effective field theory (EFT) provides a model-independent framework for interpreting the results of dark matter (DM) direct detection experiments. In this study, we demonstrate that the two fermionic DM-quark tensor operators, $(\bar{\chi} i\sigma^{\mu\nu} \gamma^5 \chi) (\bar{q} \sigma_{\mu\nu}q)$ and $(\bar{\chi} \sigma^{\mu\nu} \chi) (\bar{q} \sigma_{\mu\nu} q)$, can contribute to the DM electric and magnetic dipole moments via nonperturbative QCD effects, in addition to the well-studied contact DM-nucleon operators. We then investigate the constraints on these two operators by considering both the contact and the dipole contributions using the XENON1T nuclear recoil and Migdal effect data. We also recast other existing bounds on the DM dipole operators, derived from electron and nuclear recoil measurements in various direct detection experiments, as constraints on the two tensor operators. For $m_\chi \lesssim 1\,\rm GeV$, our results significantly extend the reach of constraints on the DM-quark tensor operators to masses as low as $5\,\rm MeV$, with the bound exceeding that obtained by the Migdal effect with only contact interactions by an order of magnitude or so. In particular, for the operator $(\bar{\chi} \sigma^{\mu\nu}i\gamma_5 \chi) (\bar{q} \sigma_{\mu\nu}q)$ with DM mass $m_\chi \gtrsim 10\,\rm GeV$, the latest PandaX constraint on the DM electric dipole moment puts more stringent bounds than the previous direct detection limit., Comment: v1: 19 pages, 4 figures; v2: We have updated discussion on nucleon form factors by using more recent numbers, without modifying our final numerical results; v3: Add comparisons with other constraints and provide a UV model for the generation of tensor operators from tree-level matching; To be published in CPC

The deposition and intercalation of metal atoms can induce superconductivity in monolayer and bilayer graphenes. For example, it has been experimentally proved that Li-deposited graphene is a superconductor with critical temperature $T_{c}$ of 5.9 K, Ca-intercalated bilayer graphene C$_{6}$CaC$_{6}$ and K-intercalated epitaxial bilayer graphene C$_{8}$KC$_{8}$ are superconductors with $T_{c}$ of 2-4 K and 3.6 K, respectively. However, the $T_{c}$ of them are relatively low. To obtain higher $T_{c}$ in graphene-based superconductors, here we predict a new Ca-intercalated bilayer graphene C$_{2}$CaC$_{2}$, which shows higher Ca concentration than the C$_{6}$CaC$_{6}$. It is proved to be thermodynamically and dynamically stable. The electronic structure, electron-phonon coupling (EPC) and superconductivity of C$_{2}$CaC$_{2}$ are investigated based on first-principles calculations. The EPC of C$_{2}$CaC$_{2}$ mainly comes from the coupling between the electrons of C-$p_{z}$ orbital and the high- and low-frequency vibration modes of C atoms. The calculated EPC constant $\lambda$ of C$_{2}$CaC$_{2}$ is 0.75, and the superconducting $T_{c}$ is 18.9 K, which is much higher than other metal-intercalated bilayer graphenes. By further applying -4\% biaxial compressive strain to C$_{2}$CaC$_{2}$, the $T_{c}$ can be boosted to 26.6 K. Thus, the predicted C$_{2}$CaC$_{2}$ provides a new platform for realizing superconductivity with the highest $T_{c}$ in bilayer graphenes., Comment: 29 pages,11 figures

We show that the rotating Rayleigh-Benard convection, where a rotating fluid is heated from below, exhibits non-Hermitian topological states. Recently, Favier and Knobloch (JFM 2020) hypothesized that the robust wall modes in rapidly rotating convection are topologically protected. We study the linear problem around the conduction profile, and by considering a Berry curvature defined in the complex wavenumber space, particularly, by introducing a complex vertical wavenumber, we find that these modes can be characterized by a non-zero integer Chern number, indicating their topological nature. The eigenvalue problem is intrinsically non-Hermitian, therefore the definition of Berry curvature generalizes that of the stably stratified problem. Moreover, the three-dimensional setup naturally regularizes the eigenvector at the infinite horizontal wavenumber. Under the hydrostatic approximation, it recovers a two-dimensional analogue of the one which explains the topological origin of the equatorial Kelvin and Yanai waves. The existence of the tenacious wall modes relies only on rotation when the fluid is stratified, no matter whether it is stable or unstable. However, the neutrally stratified system does not support a topological edge state. In addition, we define a winding number to visualize the topological nature of the fluid., Comment: 16 pages, 3 figures

We introduce the effective field theory of two different light dark particles interacting with the standard model (SM) light states in a single vertex, termed dark sector effective field theory (DSEFT). We focus on the new light particles with spin up to 1 and being real in essence, namely, new real scalars $\phi$ and $S$, Majorana fermions $\chi$ and $\psi$, and real vectors $X_\mu$ and $V_\mu$. In the framework of low energy effective field theory with QED and QCD symmetry, the DSEFT can be classified into six categories, including the scalar-scalar-SM ($\phi S$-SM), fermion-fermion-SM ($\chi\psi$-SM), vector-vector-SM ($X V$-SM), scalar-fermion-SM ($\phi \chi$-SM), scalar-vector-SM ($\phi X$-SM), and fermion-vector-SM ($\chi X$-SM) cases. For each case, we construct the effective operator basis up to canonical dimension 7, which will cover most interesting phenomenology at low energy. As a phenomenological example, we investigate the longstanding neutron lifetime anomaly through the neutron dark decay modes $n \to \chi \phi \text{ or } \chi X$ from the effective interactions in the fermion-scalar-SM or fermion-vector-SM case. When treating the light fermion as a dark matter candidate, we also explore the constraints from DM-neutron annihilation signal at Super-Kamiokande. We find the neutron dark decay in each scenario can accommodate the anomaly, at the same time, without contradicting with the Super-Kamiokande limit., Comment: 33 pages, 4 figures, typos are corrected and several new references are included. To appear in JHEP

We study the statistically steady states of the forced dissipative three-dimensional homogeneous isotropic turbulence at scales larger than the forcing scale in real separation space. The probability density functions (PDFs) of longitudinal velocity difference at large separations are close to but deviate from Gaussian, measured by their non-zero odd parts. Under the assumption that forcing controls the large-scale dynamics, we propose a conjugate regime to Kolmogorov's inertial range, independent of the forcing scale, to capture the odd parts of PDFs. The analytical expressions of the third-order longitudinal structure functions derived from the K\'arm\'an-Howarth-Monin equation prove that the odd-part PDFs of velocity differences at large separations are small but non-zero, and show that the odd-order longitudinal structure functions have a universal power-law decay with exponent $-2$ as the separation tends to infinity regardless of the particular forcing form, implying a significant coupling between large and small scales. Thus, dynamics of large scales depart from the absolute equilibrium, and we can partially recover small-scale information without explicitly resolving small-scale dynamics. The departure from the statistical equilibrium is quantified and found to be viscosity independent. Even though this departure is small, it is significant and should be considered when studying the large scales of the forced three-dimensional homogeneous isotropic turbulence.

This paper introduces FedSecurity, an end-to-end benchmark that serves as a supplementary component of the FedML library for simulating adversarial attacks and corresponding defense mechanisms in Federated Learning (FL). FedSecurity eliminates the need for implementing the fundamental FL procedures, e.g., FL training and data loading, from scratch, thus enables users to focus on developing their own attack and defense strategies. It contains two key components, including FedAttacker that conducts a variety of attacks during FL training, and FedDefender that implements defensive mechanisms to counteract these attacks. FedSecurity has the following features: i) It offers extensive customization options to accommodate a broad range of machine learning models (e.g., Logistic Regression, ResNet, and GAN) and FL optimizers (e.g., FedAVG, FedOPT, and FedNOVA); ii) it enables exploring the effectiveness of attacks and defenses across different datasets and models; and iii) it supports flexible configuration and customization through a configuration file and some APIs. We further demonstrate FedSecurity's utility and adaptability through federated training of Large Language Models (LLMs) to showcase its potential on a wide range of complex applications.

Solar filament eruptions, flares and coronal mass ejections (CMEs) are manifestations of drastic release of energy in the magnetic field, which are related to many eruptive phenomena from the Earth magnetosphere to black hole accretion disks. With the availability of high-resolution magnetograms on the solar surface, observational data-based modelling is a promising way to quantitatively study the underlying physical mechanisms behind observations. By incorporating thermal conduction and radiation losses in the energy equation, we develop a new data-driven radiative magnetohydrodynamic (MHD) model, which has the capability to capture the thermodynamic evolution compared to our previous zero-\b{eta} model. Our numerical results reproduce major observational characteristics of the X1.0 flare on 2021 October 28 in NOAA active region (AR) 12887, including the morphology of the eruption, kinematic of flare ribbons, extreme-ultraviolet (EUV) radiations, and two components of the EUV waves predicted by the magnetic stretching model, i.e., a fast-mode shock wave and a slower apparent wave due to successive stretching of magnetic field lines. Moreover, some intriguing phenomena are revealed in the simulation. We find that flare ribbons separate initially and ultimately stop at the outer stationary quasi-separatrix layers (QSLs). Such outer QSLs correspond to the border of the filament channel and determine the final positions of flare ribbons, which can be used to predict the size and the lifetime of a flare before it occurs. In addition, the side view of the synthesized EUV and white-light images exhibit typical three-part structures of CMEs, where the bright leading front is roughly cospatial with the non-wave component of the EUV wave, reinforcing the magnetic stretching model for the slow component of EUV waves., Comment: 46 pages, 15 figures, accepted for publication in The Astrophysical Journal Supplement

We explore the spatial-temporal structure functions of Burgers turbulence driven by a temporal Ornstein-Uhlenbeck (OU) process, where the characteristic time scale of the OU process is much larger than that of the energy flux across spatial scales. Based on the K\'arm\'an-Howarth-Monin equation, we obtain an expression for the third-order spatial-temporal structure function away from the dissipation scale. This expression combines Kolmogorov's exact result of spatial structure function and the exponential temporal decay of the external force. We numerically justify this expression and find that the high-order structure functions also decay exponentially, however, the dependence of decay rates on order is different for the odd- and even-order structure functions. Comparing the OU-driven Burgers turbulence with that driven by temporal white noise, their spatial structure functions are identical when the energy injection rates are the same, which justifies Kolmogorov's theory, but these two systems' temporal structure functions differ. Also, the velocity pdf in the OU-driven Burgers turbulence shows a bimodal distribution, contradicting the near-Gaussian distribution in white-noise-driven turbulence., Comment: 18 pages, 9 figures

The Ocean's turbulent energy cycle has a paradox; large-scale eddies under the control of Earth's rotation primarily transfer kinetic energy (KE) to larger scales via an inverse cascade, while a transfer to smaller scales is needed to accomplish dissipation. It has been argued, using numerical simulations, that fronts, waves and other turbulent structures can produce a forward cascade of KE toward dissipation scales. However, this forward cascade and its coexistence with known inverse cascade were not confirmed in observations. Here we present the first evidence of a dual KE cascade in the Ocean by analyzing velocity measurements from surface drifters released in the Gulf of Mexico. Our results show that KE is injected at two dominant scales and transferred to both large and small scales, with the downscale flux dominating at scales smaller than ~1-10km. The cascade rates are modulated seasonally, with stronger KE injection and forward transfer during winter., Comment: No comments

Since the famous work by Kolmogorov on incompressible turbulence, the structure-function theory has been a key foundation of modern turbulence study. Due to the simplicity of Burgers turbulence, structure functions are calculated to arbitrary orders, which provides numerous implications for other compressible turbulent systems. We present the derivation of exact forcing-scale resolving expressions for high-order structure functions of the burgers turbulence. Compared with the previous theories where the structure functions are calculated in the inertial range based on the statistics of shocks, our expressions link high-order structure functions in different orders without extra information on the flow structure and are valid beyond the inertial range, therefore they are easily checked by numerical simulations., Comment: 11 pages 4 figures

In this paper, we present an autonomous navigation system for goal-driven exploration of unknown environments through deep reinforcement learning (DRL). Points of interest (POI) for possible navigation directions are obtained from the environment and an optimal waypoint is selected, based on the available data. Following the waypoints, the robot is guided towards the global goal and the local optimum problem of reactive navigation is mitigated. Then, a motion policy for local navigation is learned through a DRL framework in a simulation. We develop a navigation system where this learned policy is integrated into a motion planning stack as the local navigation layer to move the robot between waypoints towards a global goal. The fully autonomous navigation is performed without any prior knowledge while a map is recorded as the robot moves through the environment. Experiments show that the proposed method has an advantage over similar exploration methods, without reliance on a map or prior information in complex static as well as dynamic environments.

Wind forcing injects energy into the mesoscale eddies and near-inertial waves (NIWs) in the ocean, and the NIW is believed to solve the puzzle of mesoscale energy budget by absorbing energy from mesoscale eddies followed by a forward cascade of NIW energy which finally dissipates at the ocean interior. This work studies the turbulent energy transfer in the NIW--quasigeostrophic mean mesoscale eddy coupled system based on a previously derived two-dimensional model which has a Hamiltonian structure and inherits conserved quantities in the Boussinesq equations (Xie \& Vanneste, \textit{J. Fluid Mech.}, vol. 774, 2015, pp. 147--169). Based on the conservation of energy, potential enstrophy and wave action, we propose a heuristic argument predicting the existence of phase transition with changing the relative strength between NIW and mean flow. By running forced-dissipative numerical simulations with varying parameter $R$, the ratio of the magnitude of NIW and mean-flow forcing, we justify the existence of phase transition, which is found to be second-order, around critical value $R_c$. When $0R_c$, energy transfers downscale, wave action transfers bidirectionally, and vortex filaments are dominant. We find the catalytic wave induction (CWI) mechanism where the NIW induces a downscale energy flux of the mean flow. The CWI mechanism differs from the stimulated loss of balance by the absence of energy conversion from the mesoscale eddy to NIW, and it is found to be effective in the toy-model study, making it potentially important for ocean energetics., Comment: 17 pages 7 figures

Estimating accurate depth from a single image is challenging because it is an ill-posed problem as infinitely many 3D scenes can be projected to the same 2D scene. However, recent works based on deep convolutional neural networks show great progress with plausible results. The convolutional neural networks are generally composed of two parts: an encoder for dense feature extraction and a decoder for predicting the desired depth. In the encoder-decoder schemes, repeated strided convolution and spatial pooling layers lower the spatial resolution of transitional outputs, and several techniques such as skip connections or multi-layer deconvolutional networks are adopted to recover the original resolution for effective dense prediction. In this paper, for more effective guidance of densely encoded features to the desired depth prediction, we propose a network architecture that utilizes novel local planar guidance layers located at multiple stages in the decoding phase. We show that the proposed method outperforms the state-of-the-art works with significant margin evaluating on challenging benchmarks. We also provide results from an ablation study to validate the effectiveness of the proposed method.

We derive and test a new heuristic theory for third-order structure functions that resolve the forcing scale in the scenario of simultaneous spectral energy transfer to both small and large scales, which can occur naturally in rotating stratified turbulence or magnetohydrodynamical~(MHD) turbulence, for example. The theory has three parameters, namely the upscale/downscale energy transfer rates and the forcing scale, and it includes the classic inertial range theories as local limits. When applied to measured data, our global-in-scale theory can deduce the energy transfer rates using the full range of data, therefore it has broader applications compared with the local theories, especially in the situations where the data is imperfect. In addition, because of the resolution of forcing scales, the new theory can detect the scales of energy input, which was impossible before. We test our new theory with a two-dimensional simulation of MHD turbulence., Comment: 10 pages, 2 figures

Motivated by the recent discovery of subsurface oceans on planetary moons and the interest they have generated, we explore convective flows in shallow spherical shells of dimensionless gap width $\varepsilon^2\ll 1$ in the rapid rotation limit $\mathrm{E}\ll1$, where $\mathrm{E}$ is the Ekman number. We employ direct numerical simulations (DNS) of the Boussinesq equations to compute the local heat flux $\mathrm{Nu}(\lambda)$ as a function of the latitude $\lambda$ and use the results to characterize the trapping of convection at low latitudes, around the equator. We show that these results are quantitatively reproduced by an asymptotically exact nonhydrostatic equatorial $\beta$-plane convection model at a much more modest computational cost than DNS. We identify the trapping parameter $\beta=\varepsilon \mathrm{E}^{-1}$ as the key parameter that controls the vigor and latitudinal extent of convection for moderate thermal forcing when $\mathrm{E}\sim\varepsilon$ and $\varepsilon\downarrow 0$. This model provides a new theoretical paradigm for nonlinear investigations., Comment: 26 pages, 10 figures. Accepted for publication in Physical Review Fluids

A simple model of nonlinear salt-finger convection in two dimensions is derived and studied. The model is valid in the limit of small solute to heat diffusivity ratio and large density ratio, which is relevant to both oceanographic and astrophysical applications. Two limits distinguished by the magnitude of the Schmidt number are found. For order one Schmidt numbers, appropriate for astrophysical applications, a modified Rayleigh-B\'enard system with large-scale damping due to a stabilizing temperature is obtained. For large Schmidt numbers, appropriate for the oceanic setting, the model combines a prognostic equation for the solute field and a diagnostic equation for inertia-free momentum dynamics. Two distinct saturation regimes are identified for the second model: The weakly driven regime is characterized by a large-scale flow associated with a balance between advection and linear instability, while the strongly driven regime produces multiscale structures, resulting in a balance between the energy input through linear instability and the energy transfer between scales. For both regimes, we analytically predict and numerically confirm the dependence of the kinetic energy and salinity fluxes on the ratio between solute and heat Rayleigh numbers. The spectra and probability density functions are also computed., Comment: 25 pages, 16 figures

This paper proposes a novel model-free screening procedure for ultrahigh dimensional data analysis. By utilizing slicing technique which has been successfully ap- plied to continuous variables, we construct a new index called the fused mean-variance for feature screening. This method has the following merits: (i) it is model-free, i.e., without specifying regression form of predictors and response variable; (ii) it can be used to analyze various types of variables including discrete, categorical and continuous vari- ables; (iii) it still works well even when the covariates/random errors are heavy-tailed or the predictors are strongly dependent. Under some regularity conditions, we establish the sure screening and rank consistency. Simulation studies are conducted to assess the performance of the proposed approach. A real data is used to illustrate the proposed method.

Mountain-generated inertia-gravity waves (IGWs) affect the dynamics of both the atmosphere and the ocean through the mean force they exert as they interact with the flow. A key to this interaction is the presence of critical-level singularities or, when planetary rotation is taken into account, inertial-level singularities, where the Doppler-shifted wave frequency matches the local Coriolis frequency. We examine the role of the latter singularities by studying the steady wavepacket generated by a multiscale mountain in a rotating linear shear flow at low Rossby number. Using a combination of WKB and saddle-point approximations, we provide an explicit description of the form of the wavepacket, of the mean forcing it induces, and of the mean-flow response. We identify two distinguished regimes of wave propagation: Regime I applies far enough from a dominant inertial level for the standard ray-tracing approximation to be valid; Regime II applies to a thin region where the wavepacket structure is controlled by the inertial-level singularities. The wave--mean-flow interaction is governed by the change in Eliassen--Palm (or pseudomomentum) flux. This change is localised in a thin inertial layer where the wavepacket takes a limiting form of that found in Regime II. We solve a quasi-geostrophic potential-vorticity equation forced by the divergence of the Eliassen--Palm flux to compute the wave-induced mean flow. Our results, obtained in an inviscid limit, show that the wavepacket reaches a large-but-finite distance downstream of the mountain (specifically, a distance of order $k_*^{1/2} \Delta^{3/2}$, where $k_*^{-1}$ and $\Delta$ measure the wave and envelope scales of the mountain) and extends horizontally over a similar scale., Comment: 26 pages, 8 figures

Wind forcing of the ocean generates a spectrum of inertia-gravity waves that is sharply peaked near the local inertial (or Coriolis) frequency. The corresponding near-inertial waves (NIWs) are highly energetic and play a significant role in the slow, large-scale dynamics of the ocean. To analyse this role, we develop a new model of the nondissipative interactions between NIWs and balanced motion. The model is derived using the generalised-Lagrangian-mean (GLM) framework (specifically, the glm variant of Soward & Roberts (2010)), taking advantage of the time-scale separation between the two types of motion to average over the short NIW period. We combine Salmon's (2013) variational formulation of GLM with Whitham averaging to obtain a system of equations governing the joint evolution of NIWs and mean flow. Assuming that the mean flow is geostrophically balanced reduces this system to a simple model coupling Young & Ben Jelloul's (1997) equation for NIWs with a modified quasi-geostrophic equation. In this coupled model, the mean flow affects the NIWs through advection and refraction; conversely, the NIWs affect the mean flow by modifying the potential-vorticity inversion - the relation between advected potential vorticity and advecting mean velocity - through a quadratic wave term, consistent with the GLM results of Buhler & McIntyre (1998). The coupled model is Hamiltonian and its conservation laws, for wave action and energy in particular, prove illuminating: on their basis, we identify a new interaction mechanism whereby NIWs forced at large scales extract energy from the balanced flow as their horizontal scale is reduced by differential advection and refraction so that their potential energy increases. A rough estimate suggests that this mechanism could provide a significant sink of energy for mesoscale motion and play a part in the global energetics of the ocean., Comment: 25 pages, 4 figures

A rigid spherical particle in an acoustic wave field oscillates at the wave period but has also a mean motion on a longer time scale. The dynamics of this mean motion is crucial for numerous applications of acoustic microfluidics, including particle manipulation and flow visualisation. It is controlled by four physical effects: acoustic (radiation) pressure, streaming, inertia and viscous drag. In this paper, we carry out a systematic multiscale analysis of the problem in order to assess the relative importance of these effects depending on the parameters of the system that include wave amplitude, wavelength, sound speed, sphere radius, and viscosity. We identify two distinguished regimes characterised by a balance among three of the four effects, and we derive the equations that govern the mean particle motion in each regime. This recovers and organises classical results by King, Gor'kov and Doinikov, clarifies the range of validity of these results, and reveals a new nonlinear dynamical regime. In this regime, the mean motion of the particle remains intimately coupled to that of the surrounding fluid, and while viscosity affects the fluid motion, it plays no part in the acoustic pressure. Simplified equations, valid when only two physical effects control the particle motion, are also derived. They are used to obtain sufficient conditions for the particle to behave as a passive tracer of the Lagrangian-mean fluid motion., Comment: 31 pages, 3 figures

In microfluidic applications involving high-frequency acoustic waves over a solid boundary, the Stokes boundary-layer thickness $\delta$ is so small that some non-negligible slip may occur at the fluid-solid interface. This paper assesses the impact of the slip by revisiting the classical problem of steady acoustic streaming over a flat boundary with the Navier boundary condition $u|_{y=0} = L_\mathrm{s} \partial_y u|_{y=0}$, where $u$ is the velocity tangent to the boundary $y=0$, and the parameter $L_\mathrm{s}$ is the slip length. The limit outside the boundary layer provide an effective slip velocity. A general expression is obtained for the streaming velocity outside the boundary layer as a function of the dimensionless parameter $L_\mathrm{s}/\delta$. Particularising to travelling and standing waves shows that the boundary slip respectively increases and decreases the streaming velocity., Comment: 11 pages, 3 figures

The greatest difficulty that is encountered by students in Thermodynamics classes is to find relationships between variables and solve a total differential equation that relates one thermodynamic state variable to two mutually-independent state variables. Rules of differentiation, including the total differential and the cyclic rule, are fundamental to solving thermodynamics problems. However, many students do not understand such basic mathematics principles, because ordinary calculus and thermodynamics textbooks do not provide specific geometric explanations for the differential rules. Spatial visualization allows students to develop a new line of reasoning. In this investigation, the total differential rule and the cyclic rule are interpreted geometrically., Comment: The paper has been withdrawn by the author due to the request from the Journal of Chemical Education according to the copyright agreement