Edgar Knobloch
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Edgar Knobloch has a total of 82 co-authors. Here are the top co-authors by number of shared publications:
- Keith Julien
University of Colorado
10 shared publications - Adrian van Kan
University of California at Berkeley
8 shared publications - Uwe Thiele
Institute of Theoretical Physics, University of Münster, Münster, Germany
8 shared publications - Cédric Beaume
Department of Aeronautics, Imperial College London, London SW7 2AZ, United Kingdom
7 shared publications - Chang Liu
University of Connecticut
6 shared publications - Hsien-Ching Kao
Wolfram Research Inc., Champaign, Illinois 61820, USA
5 shared publications
Latest Publications
- 13 November 2026
- Published by American Physical Society (APS) in Physical Review E
- Vol. 114 (5)
Accurate simulation of pulled and pushed fronts in the nonautonomous Fisher-Kolmogorov-Petrovsky-Piskunov equation
- 13 November 2026
- Published by American Physical Society (APS) in Physical Review E
- Vol. 114 (5)
Accurate simulation of pulled and pushed fronts in the nonautonomous Fisher-Kolmogorov-Petrovsky-Piskunov equation
- Open Access
- 4 September 2026
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1042
Quasi-geostrophic Rayleigh–Bénard convection on the tilted bold italic f $\boldsymbol{f}$ -planeRapidly rotating Rayleigh–Bénard convection on an f $f$ -plane at colatitude theta Subscript f $\vartheta _{\!f}$ is investigated numerically using an asymptotically reduced equation set valid in the limit of very rapid rotation. The equations provide a non-hydrostatic but quasi-geostrophic description in a non-orthogonal coordinate system. The tilt changes the structure of the large-scale barotropic condensate from large-scale vortices to zonal flows as the colatitude of the f $f$ -plane increases, with bistable states present for certain parameter ranges, extending prior work to a geophysically significant parameter regime. This behaviour is understood through the impact of broken rotation symmetry on the barotropic source terms resulting from baroclinic vorticity stresses and baroclinic torque. As the tilt angle theta Subscript f $\vartheta _{\!f}$ increases, global heat and momentum transport is reduced relative to upright-polar convection, a result that is explained through linear theory and nonlinear power maps both of which demonstrate increased attenuation of the domain of dynamically active spatial scales as the convective modes depart from a north–south alignment in the horizontal plane. A key finding is that the predominance of lateral thermal mixing allows for the maintenance of a persistent unstable mean temperature gradient that saturates at increasing forcing levels and remains insensitive to the colatitude.- 4 September 2026
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1042
Open AccessQuasi-geostrophic Rayleigh–Bénard convection on the tilted bold italic f $\boldsymbol{f}$ -planeRapidly rotating Rayleigh–Bénard convection on an f $f$ -plane at colatitude theta Subscript f $\vartheta _{\!f}$ is investigated numerically using an asymptotically reduced equation set valid in the limit of very rapid rotation. The equations provide a non-hydrostatic but quasi-geostrophic description in a non-orthogonal coordinate system. The tilt changes the structure of the large-scale barotropic condensate from large-scale vortices to zonal flows as the colatitude of the f $f$ -plane increases, with bistable states present for certain parameter ranges, extending prior work to a geophysically significant parameter regime. This behaviour is understood through the impact of broken rotation symmetry on the barotropic source terms resulting from baroclinic vorticity stresses and baroclinic torque. As the tilt angle theta Subscript f $\vartheta _{\!f}$ increases, global heat and momentum transport is reduced relative to upright-polar convection, a result that is explained through linear theory and nonlinear power maps both of which demonstrate increased attenuation of the domain of dynamically active spatial scales as the convective modes depart from a north–south alignment in the horizontal plane. A key finding is that the predominance of lateral thermal mixing allows for the maintenance of a persistent unstable mean temperature gradient that saturates at increasing forcing levels and remains insensitive to the colatitude. - Open Access
- 30 July 2026
- Published by The Royal Society in Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
- Vol. 384 (2325)
Two-dimensional turbulent condensates without bottom drag
The extent to which equilibrium statistical theory is applicable to driven dissipative dynamics remains an important open question in many systems. We use extensive direct numerical simulations (DNS) of the incompressible two-dimensional (2D) Navier–Stokes equation to examine the steady state of large-scale condensates in 2D turbulence at finite Reynolds number Re in the absence of bottom drag. Large-scale condensates appear above a critical Reynolds number Rec≈4.19. Close to this onset, we find a power-law scaling of the energy with Re-Rec, with the energy spectrum at large scales following the absolute equilibrium form proposed by Kraichnan. At larger Re, the energy spectrum deviates from this form, displaying a steep power-law range at low wavenumbers with exponent -5, with most of the energy dissipation occurring within the condensate at large scales. We show that this spectral exponent is consistent with the logarithmic radial vorticity profile of the condensate vortices predicted by quasi-linear theory for a viscously saturated condensate. Our findings shed new light on the classical problem of large-scale turbulent condensation in forced dissipative 2D flows in finite domains, showing that the large scales are close to equilibrium dynamics in weakly turbulent flows but not in the strong condensate regime with Re≫1. This article is part of the theme issue ‘Frontiers of turbulence and statistical physics’.
- 30 July 2026
- Published by The Royal Society in Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
- Vol. 384 (2325)
Open AccessTwo-dimensional turbulent condensates without bottom drag
The extent to which equilibrium statistical theory is applicable to driven dissipative dynamics remains an important open question in many systems. We use extensive direct numerical simulations (DNS) of the incompressible two-dimensional (2D) Navier–Stokes equation to examine the steady state of large-scale condensates in 2D turbulence at finite Reynolds number Re in the absence of bottom drag. Large-scale condensates appear above a critical Reynolds number Rec≈4.19. Close to this onset, we find a power-law scaling of the energy with Re-Rec, with the energy spectrum at large scales following the absolute equilibrium form proposed by Kraichnan. At larger Re, the energy spectrum deviates from this form, displaying a steep power-law range at low wavenumbers with exponent -5, with most of the energy dissipation occurring within the condensate at large scales. We show that this spectral exponent is consistent with the logarithmic radial vorticity profile of the condensate vortices predicted by quasi-linear theory for a viscously saturated condensate. Our findings shed new light on the classical problem of large-scale turbulent condensation in forced dissipative 2D flows in finite domains, showing that the large scales are close to equilibrium dynamics in weakly turbulent flows but not in the strong condensate regime with Re≫1. This article is part of the theme issue ‘Frontiers of turbulence and statistical physics’.
- 1 July 2026
- Published by AIP Publishing in Chaos: An Interdisciplinary Journal of Nonlinear Science
- Vol. 36 (7)
Localized dynamic patterns in the complex cubic–quintic Swift–Hohenberg equation in one and two dimensions
We use bifurcation analysis and numerical methods to study dynamic patterns in a complex Swift–Hohenberg equation focusing on structures present on a 1D interval with periodic boundary conditions and on a 2D disk with Neumann boundary conditions. In 1D, the equation features Hopf bifurcations from the trivial state u ≡ 0 at finite wave number k ≠ 0 (also known as wave bifurcations), leading to a bifurcation problem with O(2)×S1 symmetry, i.e., symmetry under translations and reflections in space, and phase rotations. The O(2) symmetry implies the simultaneous appearance of traveling waves and standing waves, and we consider the case where both bifurcations are subcritical and yield secondary bifurcations to spatially localized structures in the form of modulated traveling waves and localized standing waves, respectively. We use numerical continuation to show that some of the solutions thus obtained exhibit homoclinic snaking associated with their gradual growth in spatial extent. We exploit gauge symmetry to compute these dynamic patterns as relative equilibria, which allows us to locate efficiently tertiary bifurcations to two-frequency localized standing waves and to localized drifting waves. The intricate behavior identified in 1D is shown to provide a “road map” for the organization of wall-attached states on 2D disks in the form of rotating and standing waves at the wall, which may further localize in angle, yielding rotating spots or stationary breathing spots at the wall, and tertiary bifurcations from these. Some direct numerical simulations are used to identify further structures, including dynamics present in the disk bulk.
- 1 July 2026
- Published by AIP Publishing in Chaos: An Interdisciplinary Journal of Nonlinear Science
- Vol. 36 (7)
Localized dynamic patterns in the complex cubic–quintic Swift–Hohenberg equation in one and two dimensions
We use bifurcation analysis and numerical methods to study dynamic patterns in a complex Swift–Hohenberg equation focusing on structures present on a 1D interval with periodic boundary conditions and on a 2D disk with Neumann boundary conditions. In 1D, the equation features Hopf bifurcations from the trivial state u ≡ 0 at finite wave number k ≠ 0 (also known as wave bifurcations), leading to a bifurcation problem with O(2)×S1 symmetry, i.e., symmetry under translations and reflections in space, and phase rotations. The O(2) symmetry implies the simultaneous appearance of traveling waves and standing waves, and we consider the case where both bifurcations are subcritical and yield secondary bifurcations to spatially localized structures in the form of modulated traveling waves and localized standing waves, respectively. We use numerical continuation to show that some of the solutions thus obtained exhibit homoclinic snaking associated with their gradual growth in spatial extent. We exploit gauge symmetry to compute these dynamic patterns as relative equilibria, which allows us to locate efficiently tertiary bifurcations to two-frequency localized standing waves and to localized drifting waves. The intricate behavior identified in 1D is shown to provide a “road map” for the organization of wall-attached states on 2D disks in the form of rotating and standing waves at the wall, which may further localize in angle, yielding rotating spots or stationary breathing spots at the wall, and tertiary bifurcations from these. Some direct numerical simulations are used to identify further structures, including dynamics present in the disk bulk.
- Open Access
- 25 September 2025
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1020, R1
Spontaneous generation of helical flows by salt fingers
We study the dynamics of salt fingers in the regime of slow salinity diffusion (small inverse Lewis number) and strong stratification (large density ratio), focusing on regimes relevant to Earth’s oceans. Using three-dimensional direct numerical simulations in periodic domains, we show that salt fingers exhibit rich, multiscale dynamics in this regime, with vertically elongated fingers that are twisted into helical shapes at large scales by mean flows and disrupted at small scales by isotropic eddies. We use a multiscale asymptotic analysis to motivate a reduced set of partial differential equations that filters internal gravity waves and removes inertia from all parts of the momentum equation except for the Reynolds stress that drives the helical mean flow. When simulated numerically, the reduced equations capture the same dynamics and fluxes as the full equations in the appropriate regime. The reduced equations enforce zero helicity in all fluctuations about the mean flow, implying that the symmetry-breaking helical flow is generated spontaneously by strictly non-helical fluctuations.
- 25 September 2025
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1020, R1
Open AccessSpontaneous generation of helical flows by salt fingers
We study the dynamics of salt fingers in the regime of slow salinity diffusion (small inverse Lewis number) and strong stratification (large density ratio), focusing on regimes relevant to Earth’s oceans. Using three-dimensional direct numerical simulations in periodic domains, we show that salt fingers exhibit rich, multiscale dynamics in this regime, with vertically elongated fingers that are twisted into helical shapes at large scales by mean flows and disrupted at small scales by isotropic eddies. We use a multiscale asymptotic analysis to motivate a reduced set of partial differential equations that filters internal gravity waves and removes inertia from all parts of the momentum equation except for the Reynolds stress that drives the helical mean flow. When simulated numerically, the reduced equations capture the same dynamics and fluxes as the full equations in the appropriate regime. The reduced equations enforce zero helicity in all fluctuations about the mean flow, implying that the symmetry-breaking helical flow is generated spontaneously by strictly non-helical fluctuations.
- 6 August 2025
- Published by Elsevier in Journal of Computational Physics
- Vol. 541
Rescaled equations for well-conditioned direct numerical simulations of rapidly rotating convection
- 6 August 2025
- Published by Elsevier in Journal of Computational Physics
- Vol. 541
Rescaled equations for well-conditioned direct numerical simulations of rapidly rotating convection
- 9 May 2025
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1010, A42
Bridging the Rossby number gap in rapidly rotating thermal convection
Geophysical and astrophysical fluid flows are typically driven by buoyancy and strongly constrained at large scales by planetary rotation. Rapidly rotating Rayleigh–Bénard convection (RRRBC) provides a paradigm for experiments and direct numerical simulations (DNS) of such flows, but the accessible parameter space remains restricted to moderately fast rotation rates (Ekman numbers ${ {Ek}} \gtrsim 10^{-8}$ ), while realistic ${Ek}$ for geo- and astrophysical applications are orders of magnitude smaller. On the other hand, previously derived reduced equations of motion describing the leading-order behaviour in the limit of very rapid rotation ( $ {Ek}\to 0$ ) cannot capture finite rotation effects, and the physically most relevant part of parameter space with small but finite ${Ek}$ has remained elusive. Here, we employ the rescaled rapidly rotating incompressible Navier–Stokes equations (RRRiNSE) – a reformulation of the Navier–Stokes–Boussinesq equations informed by the scalings valid for ${Ek}\to 0$ , recently introduced by Julien et al. (2024) – to provide full DNS of RRRBC at unprecedented rotation strengths down to $ {Ek}=10^{-15}$ and below, revealing the disappearance of cyclone–anticyclone asymmetry at previously unattainable Ekman numbers ( ${Ek}\approx 10^{-9}$ ). We also identify an overshoot in the heat transport as ${Ek}$ is varied at fixed $\widetilde { {Ra}} \equiv {Ra}{Ek}^{4/3}$ , where $Ra$ is the Rayleigh number, associated with dissipation due to ageostrophic motions in the boundary layers. The simulations validate theoretical predictions based on thermal boundary layer theory for RRRBC and show that the solutions of RRRiNSE agree with the reduced equations at very small ${Ek}$ . These results represent a first foray into the vast, largely unexplored parameter space of very rapidly rotating convection rendered accessible by RRRiNSE.- 9 May 2025
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1010, A42
Bridging the Rossby number gap in rapidly rotating thermal convection
Geophysical and astrophysical fluid flows are typically driven by buoyancy and strongly constrained at large scales by planetary rotation. Rapidly rotating Rayleigh–Bénard convection (RRRBC) provides a paradigm for experiments and direct numerical simulations (DNS) of such flows, but the accessible parameter space remains restricted to moderately fast rotation rates (Ekman numbers ${ {Ek}} \gtrsim 10^{-8}$ ), while realistic ${Ek}$ for geo- and astrophysical applications are orders of magnitude smaller. On the other hand, previously derived reduced equations of motion describing the leading-order behaviour in the limit of very rapid rotation ( $ {Ek}\to 0$ ) cannot capture finite rotation effects, and the physically most relevant part of parameter space with small but finite ${Ek}$ has remained elusive. Here, we employ the rescaled rapidly rotating incompressible Navier–Stokes equations (RRRiNSE) – a reformulation of the Navier–Stokes–Boussinesq equations informed by the scalings valid for ${Ek}\to 0$ , recently introduced by Julien et al. (2024) – to provide full DNS of RRRBC at unprecedented rotation strengths down to $ {Ek}=10^{-15}$ and below, revealing the disappearance of cyclone–anticyclone asymmetry at previously unattainable Ekman numbers ( ${Ek}\approx 10^{-9}$ ). We also identify an overshoot in the heat transport as ${Ek}$ is varied at fixed $\widetilde { {Ra}} \equiv {Ra}{Ek}^{4/3}$ , where $Ra$ is the Rayleigh number, associated with dissipation due to ageostrophic motions in the boundary layers. The simulations validate theoretical predictions based on thermal boundary layer theory for RRRBC and show that the solutions of RRRiNSE agree with the reduced equations at very small ${Ek}$ . These results represent a first foray into the vast, largely unexplored parameter space of very rapidly rotating convection rendered accessible by RRRiNSE.- 17 March 2025
- Published by American Physical Society (APS) in Physical Review Fluids
- Vol. 10 (3)
Traveling spatially localized convective structures in an inclined porous medium
- 17 March 2025
- Published by American Physical Society (APS) in Physical Review Fluids
- Vol. 10 (3)
Traveling spatially localized convective structures in an inclined porous medium
- 27 June 2024
- Published by American Physical Society (APS) in Physical Review Fluids
- Vol. 9 (6) , 064605
Fluctuation-induced transitions in anisotropic two-dimensional turbulence
Two-dimensional (2D) turbulence features an inverse energy cascade that produces large-scale flow structures such as large-scale vortices (LSVs) and unidirectional jets. We investigate the dynamics of such large-scale structures using extensive direct numerical simulations (DNS) of stochastically forced, viscously damped 2D turbulence within a periodic rectangular (Cartesian) domain . LSVs form and dominate the system when the domain aspect ratio , while unidirectional jets predominate at . At intermediate values of , both structures are metastable, and fluctuation-induced transitions between LSVs and jets are observed. Based on large-scale energy balance in the condensate, we derive and verify predictions for the dependence of the total kinetic energy and the flow polarity on the nondimensional control parameters. We further collect detailed statistics on the lifetimes of LSVs and jets from DNS runs of up to viscous diffusive time units in length. The distribution of the lifetimes is consistent with that of a memoryless Poisson process. The data are compatible with an exponential dependence of the mean lifetime on the aspect ratio . In addition, the mean lifetimes depend sensitively on the Reynolds number : As increases, the energy gap between LSV (lower energy) and jet states (higher energy) arising from anisotropic dissipation increases, leading to an increase in lifetimes that is approximately exponential in for both LSVs and jets. Similarly, as the ratio of the forcing scale to the domain size increases, the transition rates increase sharply, confirming earlier findings. We investigate the transition dynamics in terms of kinetic energy, flow polarity, modal amplitude, and 2D phase-space diagrams, revealing that the transitions occur in two stages: In the initial stage, an efficient redistribution of kinetic energy by nonlinear triadic interactions facilitates a rapid transition from LSVs to jets and vice versa. In the second stage, the kinetic energy of the newly formed structure slowly adjusts to its associated (higher or lower) equilibrium value on a longer, viscous timescale, leading to a time delay that results in hysteretic transition behavior. Fluctuation-induced transitions may also occur between different numbers of jets. Our findings shed new light on the dynamics of coherent large-scale structures in anisotropic turbulence.- 27 June 2024
- Published by American Physical Society (APS) in Physical Review Fluids
- Vol. 9 (6) , 064605
Fluctuation-induced transitions in anisotropic two-dimensional turbulence
Two-dimensional (2D) turbulence features an inverse energy cascade that produces large-scale flow structures such as large-scale vortices (LSVs) and unidirectional jets. We investigate the dynamics of such large-scale structures using extensive direct numerical simulations (DNS) of stochastically forced, viscously damped 2D turbulence within a periodic rectangular (Cartesian) domain . LSVs form and dominate the system when the domain aspect ratio , while unidirectional jets predominate at . At intermediate values of , both structures are metastable, and fluctuation-induced transitions between LSVs and jets are observed. Based on large-scale energy balance in the condensate, we derive and verify predictions for the dependence of the total kinetic energy and the flow polarity on the nondimensional control parameters. We further collect detailed statistics on the lifetimes of LSVs and jets from DNS runs of up to viscous diffusive time units in length. The distribution of the lifetimes is consistent with that of a memoryless Poisson process. The data are compatible with an exponential dependence of the mean lifetime on the aspect ratio . In addition, the mean lifetimes depend sensitively on the Reynolds number : As increases, the energy gap between LSV (lower energy) and jet states (higher energy) arising from anisotropic dissipation increases, leading to an increase in lifetimes that is approximately exponential in for both LSVs and jets. Similarly, as the ratio of the forcing scale to the domain size increases, the transition rates increase sharply, confirming earlier findings. We investigate the transition dynamics in terms of kinetic energy, flow polarity, modal amplitude, and 2D phase-space diagrams, revealing that the transitions occur in two stages: In the initial stage, an efficient redistribution of kinetic energy by nonlinear triadic interactions facilitates a rapid transition from LSVs to jets and vice versa. In the second stage, the kinetic energy of the newly formed structure slowly adjusts to its associated (higher or lower) equilibrium value on a longer, viscous timescale, leading to a time delay that results in hysteretic transition behavior. Fluctuation-induced transitions may also occur between different numbers of jets. Our findings shed new light on the dynamics of coherent large-scale structures in anisotropic turbulence.- 1 May 2024
- Published by AIP Publishing in Chaos: An Interdisciplinary Journal of Nonlinear Science
- Vol. 34 (5)
Emergence of rogue-like waves in a reaction-diffusion system: Stochastic output from deterministic dissipative dynamics
Rogue waves are an intriguing nonlinear phenomenon arising across different scales, ranging from ocean waves through optics to Bose–Einstein condensates. We describe the emergence of rogue wave-like dynamics in a reaction-diffusion system that arise as a result of a subcritical Turing instability. This state is present in a regime where all time-independent states are unstable and consists of intermittent excitation of spatially localized spikes, followed by collapse to an unstable state and subsequent regrowth. We characterize the spatiotemporal organization of spikes and show that in sufficiently large domains the dynamics are consistent with a memoryless process.
- 1 May 2024
- Published by AIP Publishing in Chaos: An Interdisciplinary Journal of Nonlinear Science
- Vol. 34 (5)
Emergence of rogue-like waves in a reaction-diffusion system: Stochastic output from deterministic dissipative dynamics
Rogue waves are an intriguing nonlinear phenomenon arising across different scales, ranging from ocean waves through optics to Bose–Einstein condensates. We describe the emergence of rogue wave-like dynamics in a reaction-diffusion system that arise as a result of a subcritical Turing instability. This state is present in a regime where all time-independent states are unstable and consists of intermittent excitation of spatially localized spikes, followed by collapse to an unstable state and subsequent regrowth. We characterize the spatiotemporal organization of spikes and show that in sufficiently large domains the dynamics are consistent with a memoryless process.