This work investigates heat transport in rotating internally heated convection, for a horizontally periodic fluid between parallel plates under no-slip and isothermal boundary conditions. The main results are the proof of lower bounds on the mean temperature, $\overline {{\langle {T} \rangle }}$, and the heat flux out of the bottom boundary, ${\mathcal {F}}_B$, at infinite Prandtl number, where the Prandtl number is the non-dimensional ratio of viscous to thermal diffusion. The lower bounds are functions of the Rayleigh number quantifying the ratio of internal heating to diffusion and the Ekman number, $E$, which quantifies the ratio of viscous diffusion to rotation. We utilise two different estimates on the vertical velocity, $w$, one pointwise in the domain (Yan, J. Math. Phys., vol. 45(7), 2004, pp. 2718–2743) and the other an integral estimate over the domain (Constantin et al., Phys. D: Non. Phen., vol. 125, 1999, pp. 275–284), resulting in bounds valid for different regions of buoyancy-to-rotation dominated convection. Furthermore, we demonstrate that similar to rotating Rayleigh–Bénard convection, for small $E$, the critical Rayleigh number for the onset of convection asymptotically scales as $E^{-4/3}$.

This paper derives the upper bound on heat transport in supergravitational turbulent thermal convection analytically and numerically. Using a piecewise background profile, the functional inequality analysis delivers a suboptimal bound $Nu\lesssim -({3\sqrt {3}}/{2})({\eta \ln (\eta )}/({1-\eta ^2}))\,Ra^{1/2}$ as $Ra\rightarrow \infty$, where $Nu$ is the Nusselt number, $\eta$ is the radius ratio of the inner cylinder to the outer cylinder ($0<\eta <1$), and $Ra$ is the Rayleigh number. A variational problem yielded from Doering–Constantin–Hopf formalism is solved asymptotically and numerically, which delivers a better upper bound than the piecewise background profile. The asymptotic analysis and numerical data indicate that the current best bound is given by $Nu\leqslant -0.106({\eta \ln (\eta )}/({(1-\eta )(1+\sqrt {\eta })^2}))\,Ra^{1/2}$. Both analytical and numerical results demonstrate that the upper bound can be significantly reduced by the curvature effect. Unlike the traditional Rayleigh–Bénard turbulence, in which the optimal perturbation yielded from the variational problem is always two-dimensional, the present study shows that three-dimensional perturbations, annular perturbations and axisymmetric perturbations can be induced by the curvature effect simultaneously. However, we show that the bound yielded from the three-dimensional variational problem is very close to the axisymmetric situation as $\eta$ increases and $Ra$ increases.

In this article, we study the following Schrödinger equation

\begin{align*}\begin{cases}-\Delta u -\frac{\mu}{|x|^2} u+\lambda u =f(u), &\text{in}~ \mathbb{R}^N\backslash\{0\},\\\int_{\mathbb{R}^{N}}|u|^{2}\mathrm{d} x=a, & u\in H^1(\mathbb{R}^{N}),\end{cases}\end{align*}

where $N\geq 3$, a > 0, and $\mu \lt \frac{(N-2)^2}{4}$. Here $\frac{1}{|x|^2} $ represents the Hardy potential (or ‘inverse-square potential’), λ is a Lagrange multiplier, and the nonlinearity function f satisfies the general Sobolev critical growth condition. Our main goal is to demonstrate the existence of normalized ground state solutions for this equation when $0 \lt \mu \lt \frac{(N-2)^2}{4}$. We also analyse the behaviour of solutions as $\mu\to0^+$ and derive the existence of normalized ground state solutions for the limiting case where µ = 0. Finally, we investigate the existence of normalized solutions when µ < 0 and analyse the asymptotic behaviour of solutions as $\mu\to 0^-$.

We consider Euler flows on two-dimensional (2-D) periodic domain and are interested in the stability, both linear and nonlinear, of a simple equilibrium given by the 2-D Taylor–Green vortex. As the first main result, numerical evidence is provided for the fact that such flows possess unstable eigenvalues embedded in the band of the essential spectrum of the linearized operator. However, the unstable eigenfunction is discontinuous at the hyperbolic stagnation points of the base flow and its regularity is consistent with the prediction of Lin (Intl Math. Res. Not., vol. 2004, issue 41, 2004, pp. 2147–2178). This eigenfunction gives rise to an exponential transient growth with the rate given by the real part of the eigenvalue followed by passage to a nonlinear instability. As the second main result, we illustrate a fundamentally different, non-modal, growth mechanism involving a continuous family of uncorrelated functions, instead of an eigenfunction of the linearized operator. Constructed by solving a suitable partial differential equation (PDE) optimization problem, the resulting flows saturate the known estimates on the growth of the semigroup related to the essential spectrum of the linearized Euler operator as the numerical resolution is refined. These findings are contrasted with the results of earlier studies of a similar problem conducted in a slightly viscous setting where only the modal growth of instabilities was observed. This highlights the special stability properties of equilibria in inviscid flows.

The present study offers a twofold contribution on counter-gradient transport (CGT) of turbulent scalar flux. First, by examining turbulent scalar mixing through synchronized particle image velocimetry and planar laser-induced fluorescence on an inclined jet in cross-flow, we clarify the previously unexplained phenomenon of CGT, revealing key flow structures, their spatial distribution and modelling implications. Statistical analysis identifies two distinct CGT regions: local cross-gradient transport in the windward shear layer and non-local effects near the wall after injection. These behaviours are driven by specific flow structures, namely Kelvin–Helmholtz vortices (local) and wake vortices (non-local), suggesting that scalar flux can be decomposed into a gradient-type term for gradient diffusion and a term for large-eddy stirring. Second, we propose a new approach for reconstruction of turbulent mean flow and scalar fields using continuous adjoint data assimilation (DA). By rectifying model-form errors through anisotropic correction under observational constraints, our DA model minimizes discrepancies between experimental measurements and numerical predictions. As expected, the introduced forcing term effectively identifies regions where traditional models fall short, particularly in the jet centreline and near-wall regions, thereby enhancing the accuracy of the mean scalar field. These enhancements occur not only within the observation region but also in unseen regions, underscoring present DA approach's reliability and practicality for reproducing mean flow behaviours from limited data. These findings lay a solid foundation for adjoint-based model-consistent data-driven methods, offering promising potential for accurately predicting complex flow scenarios like film cooling.

We consider an internally heated fluid between parallel plates with fixed thermal fluxes. For a large class of heat sources that vary in the direction of gravity, we prove that $\smash { \smash {{\langle {\delta T} \rangle _h}} } \geq \sigma R^{-1/3} - \mu$, where $\smash { \smash {{\langle {\delta T} \rangle _h}} }$ is the average temperature difference between the bottom and top plates, $R$ is a ‘flux’ Rayleigh number and the constants $\sigma,\mu >0$ depend on the geometric properties of the internal heating. This result implies that mean downward conduction (for which $\smash { \smash {{\langle {\delta T} \rangle _h}} }< 0$) is impossible for a range of Rayleigh numbers smaller than a critical value $R_0:=(\sigma /\mu )^{3}$. The bound demonstrates that $R_0$ depends on the heating distribution and can be made arbitrarily large by concentrating the heating near the bottom plate. However, for any given fixed heating profile of the class we consider, the corresponding value of $R_0$ is always finite. This points to a fundamental difference between internally heated convection and its limiting case of Rayleigh–Bénard convection with fixed-flux boundary conditions, for which $\smash {{\langle {\delta T} \rangle _h}}$ is known to be positive for all $R$.

The Choquard equation is a partial differential equation that has gained significant interest and attention in recent decades. It is a nonlinear equation that combines elements of both the Laplace and Schrödinger operators, and it arises frequently in the study of numerous physical phenomena, from condensed matter physics to nonlinear optics.

In particular, the steady states of the Choquard equation were thoroughly investigated using a variational functional acting on the wave functions.

In this article, we introduce a dual formulation for the variational functional in terms of the potential induced by the wave function, and use it to explore the existence of steady states of a multi-state version the Choquard equation in critical and sub-critical cases.

We use the entropy method to analyse the nonlinear dynamics and stability of a continuum kinetic model of an active nematic suspension. From the time evolution of the relative entropy, an energy-like quantity in the kinetic model, we derive a variational bound on relative entropy fluctuations that can be expressed in terms of orientational order parameters. From this bound we show isotropic suspensions are nonlinearly stable for sufficiently low activity, and derive upper bounds on spatiotemporal averages in the unstable regime that are consistent with fully nonlinear simulations. This work highlights the self-organising role of activity in particle suspensions, and places limits on how organised such systems can be.

Applying a variational analysis, a minimum-enstrophy vortex in three-dimensional (3-D) fluids with continuous stratification is found, under the quasi-geostrophic hypothesis. The buoyancy frequency is held constant. This vortex is an ideal limiting state in a flow with an enstrophy decay while energy and generalized angular momentum remain fixed. The variational method used to obtain two-dimensional (2-D) minimum-enstrophy vortices is applied here to 3-D integral quantities. The solution from the first-order variation is expanded on a basis of orthogonal spherical Bessel functions. By computing second-order variations, the solution is found to be a true minimum in enstrophy. This solution is weakly unstable when inserted in a numerical code of the quasi-geostrophic equations. After a stage of linear instability, nonlinear wave interaction leads to the reorganization of this vortex into a tripolar vortex. Further work will relate our solution with maximal entropy 3-D vortices.

Successive drops of coloured ink mixed with surfactant are deposited onto a thin film of water to create marbling patterns in the Japanese art technique of Suminagashi. To understand the physics behind this and other applications where surfactant transports adsorbed passive matter at gas–liquid interfaces, we investigate the Lagrangian trajectories of material particles on the surface of a thin film of a confined viscous liquid under Marangoni-driven spreading by an insoluble surfactant. We study a model problem in which several deposits of exogenous surfactant simultaneously spread on a bounded rectangular surface containing a pre-existing endogenous surfactant. We derive Eulerian and Lagrangian formulations of the equations governing the Marangoni-driven surface flow. Both descriptions show how confinement can induce drift and flow reversal during spreading. The Lagrangian formulation captures trajectories without the need to calculate surfactant concentrations; however, concentrations can still be inferred from the Jacobian of the map from initial to current particle position. We explore a link between thin-film surfactant dynamics and optimal transport theory to find the approximate equilibrium locations of material particles for any given initial condition by solving a Monge–Ampère equation. We find that as the endogenous surfactant concentration $\delta$ vanishes, the equilibrium shapes of deposits using the Monge–Ampère approximation approach polygons with corners curving in a self-similar manner over lengths scaling as $\delta ^{1/2}$. We explore how Suminagashi patterns may be produced by using computationally efficient successive solutions of the Monge–Ampère equation.

A variational principle is proposed to derive the governing equations for the problem of ocean wave interactions with a floating ice shelf, where the ice shelf is modelled by the full linear equations of elasticity and has an Archimedean draught. The variational principle is used to form a thin-plate approximation for the ice shelf, which includes water–ice coupling at the shelf front and extensional waves in the shelf, in contrast to the benchmark thin-plate approximation for ocean wave interactions with an ice shelf. The thin-plate approximation is combined with a single-mode approximation in the water, where the vertical motion is constrained to the eigenfunction that supports propagating waves. The new terms in the approximation are shown to have a major impact on predictions of ice shelf strains for wave periods in the swell regime.

This paper presents a novel Hamiltonian formulation of the isotropic Navier–Stokes problem based on a minimum-action principle derived from the principle of least squares. This formulation uses the velocities $u_{i}(x_{j},t)$ and pressure $p(x_{j},t)$ as the field quantities to be varied, along with canonically conjugate momenta deduced from the analysis. From these, a conserved Hamiltonian functional $H^{*}$ satisfying Hamilton's canonical equations is constructed, and the associated Hamilton–Jacobi equation is formulated for both compressible and incompressible flows. This Hamilton–Jacobi equation reduces the problem of finding four separate field quantities ($u_{i}$,$p$) to that of finding a single scalar functional in those fields – Hamilton's principal functional ${S}^{*}[u_{i},p,t]$. Moreover, the transformation theory of Hamilton and Jacobi now provides a prescribed recipe for solving the Navier–Stokes problem: find ${S}^{*}$. If an analytical expression for ${S}^{*}$ can be obtained, it will lead via canonical transformation to a new set of fields which are simply equal to their initial values, giving analytical expressions for the original velocity and pressure fields. Failing that, if one can only show that a complete solution to this Hamilton–Jacobi equation does or does not exist, that will also resolve the question of existence of solutions. The method employed here is not specific to the Navier–Stokes problem or even to classical mechanics, and can be applied to any traditionally non-Hamiltonian problem.

Cavitation inception originates from nuclei in a liquid. This paper proposes a Gibbs free energy approach that provides a smooth transition from homogeneous to heterogeneous nucleation when gas is present. The impact of gas content on nucleation is explored. It is found that the gas content stabilises nuclei, a phenomenon not present in pure liquid–vapour systems. This reduces the energy barrier over that required to nucleate a vapour bubble. Different gas saturation levels are studied. Gas content can significantly reduce the energy barrier required for nucleation, and under certain circumstances eliminate it. An analytic solution for the critical radius and activation energy is obtained that accounts for gas content. The classical Blake radius is recovered as a limiting case. The hysteresis between incipience and desinence is explained using the asymmetry observed in the critical radii. The solution is used to obtain the initial bubble radius, given a critical pressure condition in cavitation susceptibility meter experiments. The relationship between initial bubble diameter and critical pressure is described by an analytic solution that accounts for gas content. A model for the derivative of the cumulative nuclei histogram with respect to bubble diameter is proposed. An analytic expression is obtained that shows good agreement with decades worth of experimental data compiled by Khoo et al. (Exp. Fluids, vol. 61, issue 2, 2020, pp. 1–20) from ocean to water tunnels. The expression recovers the $-4$ power law that is observed experimentally.

Capillary folding is the process of folding planar objects into three-dimensional (3-D) structures using capillary force. It has many important industrial applications, e.g. the fabrication of microelectromechanical systems. In this work, we propose a 3-D model for the capillary folding of thin elastic sheets with pinned contact lines. The energy of the system consists of interfacial energies between the different phases and the elastic energy of the sheet. The latter is described by the nonlinear Koiter's model, which can accommodate large deformations of the sheet. From the energy, we derive the governing equations for the static system using a variational approach. We then develop a numerical method to find equilibrium solutions via a relaxation dynamics. At each time step, we evolve the sheet by using the subdivision element method, and update the droplet surface by minimizing a squared area functional using linear finite elements. Qualitatively, numerical solutions for the equilibrium configurations of the sheet–droplet system agree well with those observed in experiments. Furthermore, we identify the critical bendabilities and present bifurcation diagrams for the folding of a triangular sheet. The results exhibit rich and fully 3-D behaviours that were not captured by previous two-dimensional models. Our results provide new insights into the nonlinear process of capillary folding, and may contribute to the advancement of microfabrication techniques.

The super-temporal-resolution (STR) reconstruction of turbulent flows is an important data augmentation application for increasing the data reach in measurement techniques and understanding turbulence dynamics. This paper proposes a data assimilation (DA) strategy based on weak-constraint four-dimensional variation to conduct an STR reconstruction in a turbulent jet beyond the Nyquist limit from given low-sampling-rate observations. Highly resolved large-eddy simulation (LES) data are used to produce synthetic measurements, which are used as observations and for validation. A segregated assimilation procedure is realised to assimilate the initial condition, inflow boundary condition and model error separately. Different types of observational data are tested. The first type is down-sampled LES data containing many small-scale turbulence structures with or without synthetic noise. The DA results show that the temporal variation of the small-scale structures is well recovered even with noise in the observations. The spectra are resolved to a frequency approximately one order of magnitude higher than what can be captured within the Nyquist limit. The second type of observation is low-sampling-rate tomographic particle image velocimetry (tomo-PIV) data with or without the injection of small-scale structures. The modulation between the large-scale structures contained in the tomo-PIV fields and the small scales injected from the observations is improved. The resultant small scales in the STR reconstruction have the characteristics of authentic turbulence to a considerable extent. Additionally, DA yields much smaller errors in the prediction of particle positions when compared with the Wiener filter, demonstrating the great potential for Lagrangian particle tracking in measurement techniques.

Invariant solutions of the Navier–Stokes equations play an important role in the spatiotemporally chaotic dynamics of turbulent shear flows. Despite the significance of these solutions, their identification remains a computational challenge, rendering many solutions inaccessible and thus hindering progress towards a dynamical description of turbulence in terms of invariant solutions. We compute equilibria of three-dimensional wall-bounded shear flows using an adjoint-based matrix-free variational approach. To address the challenge of computing pressure in the presence of solid walls, we develop a formulation that circumvents the explicit construction of pressure and instead employs the influence matrix method. Together with a data-driven convergence acceleration technique based on dynamic mode decomposition, this yields a practically feasible alternative to state-of-the-art Newton methods for converging equilibrium solutions. We compute multiple equilibria of plane Couette flow starting from inaccurate guesses extracted from a turbulent time series. The variational method outperforms Newton(-hookstep) iterations in converging successfully from poor initial guesses, suggesting a larger convergence radius.

This paper first studies the multiplicity of normalized solutions to the non-autonomous Schrödinger equation with mixed nonlinearities

\begin{equation*}\begin{cases}-\Delta u=\lambda u+h(\epsilon x)|u|^{q-2}u+\eta |u|^{p-2}u,\quad x\in \mathbb{R}^N, \\\int_{\mathbb{R}^N}|u|^2\,\textrm{d}x=a^2,\end{cases}\end{equation*}

$a, \epsilon, \eta \gt 0$

$\lambda\in \mathbb{R}$

$p=2N/(N-2)$

The Schrödinger–Poisson system describes standing waves for the nonlinear Schrödinger equation interacting with the electrostatic field. In this paper, we are concerned with the existence of positive ground states to the planar Schrödinger–Poisson system with a nonlinearity having either a subcritical or a critical exponential growth in the sense of Trudinger–Moser. A feature of this paper is that neither the finite steep potential nor the reaction satisfies any symmetry or periodicity hypotheses. The analysis developed in this paper seems to be the first attempt in the study of planar Schrödinger–Poisson systems with lack of symmetry.

In this note we provide an analytical proof of the zero helicity condition for systems governed by the Gross–Pitaevskii equation (GPE). The proof is based on the hydrodynamic interpretation of the GPE, and the direct use of Noether's theorem by applying Kleinert's multi-valued gauge theory. As a by-product we also demonstrate the conservation and quantization of the circulation for the GPE.

We investigate the equilibrium of an axisymmetric system consisting of sessile and pendent drops on pre-stretched nonlinear elastic membranes. The membrane experiences large deformations due to the drop's weight and interfacial interactions. We first show that force balance alone leads to non-unique equilibrium solutions. Identifying the system's equilibrium with the minimum of its free energy, we then demonstrate that the equilibrium solution is made unique by requiring the continuity of meridional stretches across the three-phase contact circle. For a special class of nonlinear elastic materials – $I_2$ materials – we then compute the equilibrium configurations of the drop–membrane system for a range of drop volumes and membrane pre-tensions. Finally, the present work facilitates two important applications: (a) the membrane's pre-tension and current tension are related exactly to help in utilizing the system as an elastocapillary probe for membrane pre-tension and (b) we suggest an experimental protocol for measuring the membrane's surface properties.