BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Keith Julien

Keith Julien appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynphysics.geo-phastro-ph.EPastro-phastro-ph.SRmath-phmath.MPmath.NAnlin.CDNumerical Analysis

Catalog footprint

What is connected

13works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

A gyroscopic polynomial basis in the sphere

Standard spectral codes for full sphere dynamics utilize a combination of spherical harmonics and a suitableradial basis to represent fluid variables. These basis functions have a rotational invariance not present ingeophysical flows. Gyroscopic alignment - alignment of dynamics along the axis of rotation - is ahallmark of geophysical fluids in the rapidly rotating regime. The Taylor-Proudman theorem, resultingfrom a dominant balance of the Coriolis force and the pressure gradient force, yields nearly invariant flows along this axial direction.In this paper we tailor a coordinate system to the cylindrical structures found in rotating spherical flows.This "spherindrical" coordinate system yields a natural hierarchy of basis functions, composed of Jacobi polynomialsin the radial and vertical direction, regular throughout the ball.We expand fluid variables using this basis and utilize sparse Jacobi polynomial algebra to implement all operatorsrelevant for partial differential equations in the spherical setting. We demonstrate the representation power ofthe basis in three eigenvalue problems for rotating fluids.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Experimental observation of the geostrophic turbulence regime of rapidly rotating convection

The competition between turbulent convection and global rotation in planetary and stellar interiors governs the transport of heat and tracers, as well as magnetic-field generation. These objects operate in dynamical regimes ranging from weakly rotating convection to the `geostrophic turbulence' regime of rapidly rotating convection. However, the latter regime has remained elusive in the laboratory, despite a worldwide effort to design ever-taller rotating convection cells over the last decade. Building on a recent experimental approach where convection is driven radiatively, we report heat transport measurements in quantitative agreement with this scaling regime, the experimental scaling-law being validated against direct numerical simulations (DNS) of the idealized setup. The scaling exponent from both experiments and DNS agrees well with the geostrophic turbulence prediction. The prefactor of the scaling-law is greater than the one diagnosed in previous idealized numerical studies, pointing to an unexpected sensitivity of the heat transport efficiency to the precise distribution of heat sources and sinks, which greatly varies from planets to stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Connections between Non-Rotating, Slowly Rotating, and Rapidly Rotating Turbulent Convection Transport Scalings

Buoyancy-driven convection is likely the dominant driver of turbulent motions in the universe, and thus, is widely studied by physicists, engineers, geophysicists and astrophysicists. Maybe unsurprisingly, these different communities discuss the gross convective behaviors in different ways, often without significant cross-talk existing between them. Here, we seek to draw connections between these communities. We do so by carrying out a set of basic scale estimations for how heat and fluid momentum transport should behave in non-rotating, slowly rotating and rapidly rotating buoyancy-driven convective environments. We find that slowly and rapidly rotating scalings can be inter-related via one parameter, the so-called convective Rossby number $\RoC$, a dissipation-free parameter measuring the importance of buoyancy driving relative to rotation. Further, we map between non-flux-based and the flux-based, buoyancy-driven scalings used by different groups. In doing so, these scalings show that there are clean connections between the different communities' approaches and that a number of the seemingly different scalings are actually synonymous with one another.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A nonlinear model for rotationally constrained convection with Ekman pumping

It is a well established result of linear theory that the influence of differing mechanical boundary conditions, i.e., stress-free or no-slip, on the primary instability in rotating convection becomes asymptotically small in the limit of rapid rotation. This is accounted for by the diminishing impact of the viscous stresses exerted within Ekman boundary layers and the associated vertical momentum transport by Ekman pumping. By contrast, in the nonlinear regime recent experiments and supporting simulations are now providing evidence that the efficiency of heat transport remains strongly influenced by Ekman pumping in the rapidly rotating limit. In this paper, a reduced model is developed for the case of low Rossby number convection in a plane layer geometry with no-slip upper and lower boundaries held at fixed temperatures. A complete description of the dynamics requires the existence of three distinct regions within the fluid layer: a geostrophically balanced interior where fluid motions are predominately aligned with the axis of rotation, Ekman boundary layers immediately adjacent to the bounding plates, and thermal wind layers driven by Ekman pumping in between. The reduced model uses a classical Ekman pumping parameterization to alleviate the need for spatially resolving the Ekman boundary layers. Results are presented for both linear stability theory and a special class of nonlinear solutions described by a single horizontal spatial wavenumber. It is shown that Ekman pumping allows for significant enhancement in the heat transport relative to that observed in simulations with stress-free boundaries. Without the intermediate thermal wind layer the nonlinear feedback from Ekman pumping would be able to generate a heat transport that diverges to infinity. This layer arrests this blowup resulting in finite heat transport at a significantly enhanced value.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convection-driven kinematic dynamos at low Rossby and magnetic Prandtl numbers

Most large-scale planetary magnetic fields are thought to be driven by low Rossby number convection of a low magnetic Prandtl number fluid. Here kinematic dynamo action is investigated with an asymptotic, rapidly rotating dynamo model for the plane layer geometry that is intrinsically low magnetic Prandtl number. The thermal Prandtl number and Rayleigh number are varied to illustrate fundamental changes in flow regime, ranging from laminar cellular convection to geostrophic turbulence in which an inverse energy cascade is present. A decrease in the efficiency of the convection to generate a dynamo, as determined by an increase in the critical magnetic Reynolds number, is observed as the buoyancy forcing is increased. This decreased efficiency may result from both the loss of correlations associated with the increasingly disordered states of flow that are generated, and boundary layer behavior that enhances magnetic diffusion locally. We find that the spatial characteristics of $α$, and thus the large-scale magnetic field, is dependent only weakly on changes in flow behavior. However, our results are limited to the linear, kinematic dynamo regime, and future simulations including the Lorentz force are therefore necessary to assess the robustness of this result. In contrast to the large-scale magnetic field, the behavior of the small-scale magnetic field is directly dependent on, and therefore shows significant variations with, the small-scale convective flow field.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convection-driven kinematic dynamos at low Rossby and magnetic Prandtl numbers: single mode solutions

The onset of dynamo action is investigated within the context of a newly developed low Rossby, low magnetic Prandtl number, convection-driven dynamo model. This multiscale model represents an asymptotically exact form of an $α^2$ mean field dynamo model in which the small-scale convection is represented explicitly by finite amplitude, single mode solutions. Both steady and oscillatory convection are considered for a variety of horizontal planforms. The kinetic helicity is observed to be a monotonically increasing function of the Rayleigh number. As a result, very small magnetic Prandtl number dynamos can be found for sufficiently large Rayleigh numbers. All dynamos are found to be oscillatory with an oscillation frequency that increases as the strength of the convection is increased and the magnetic Prandtl number is reduced. Kinematic dynamo action is strongly controlled by the profile of the helicity; single mode solutions which exhibit boundary layer behavior in the helicity show a decrease in the efficiency of dynamo action due to the enhancement of magnetic diffusion in the boundary layer regions. For a given value of the Rayleigh number, lower magnetic Prandtl number dynamos are excited for the case of oscillatory convection in comparison to steady convection. With regards to planetary dynamos, these results suggest that the low magnetic Prandtl number dynamos typical of liquid metals are more easily driven by thermal convection than by compositional convection.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The effects of Ekman pumping on quasi-geostrophic Rayleigh-Benard convection

Numerical simulations of 3D, rapidly rotating Rayleigh-Benard convection are performed using an asymptotic quasi-geostrophic model that incorporates the effects of no-slip boundaries through (i) parameterized Ekman pumping boundary conditions, and (ii) a thermal wind boundary layer that regularizes the enhanced thermal fluctuations induced by pumping. The fidelity of the model, obtained by an asymptotic reduction of the Navier-Stokes equations that implicitly enforces a pointwise geostrophic balance, is explored for the first time by comparisons of simulations against the findings of direct numerical simulations and laboratory experiments. Results from these methods have established Ekman pumping as the mechanism responsible for significantly enhancing the vertical heat transport. This asymptotic model demonstrates excellent agreement over a range of thermal forcing for Pr ~1 when compared with results from experiments and DNS at maximal values of their attainable rotation rates, as measured by the Ekman number (E ~ 10^{-7}); good qualitative agreement is achieved for Pr > 1. Similar to studies with stress-free boundaries, four spatially distinct flow morphologies exists. Despite the presence of frictional drag at the upper and/or lower boundaries, a strong non-local inverse cascade of barotropic (i.e., depth-independent) kinetic energy persists in the final regime of geostrophic turbulence and is dominant at large scales. For mixed no-slip/stress-free and no-slip/no-slip boundaries, Ekman friction is found to attenuate the efficiency of the upscale energy transport and, unlike the case of stress-free boundaries, rapidly saturates the barotropic kinetic energy. For no-slip/no-slip boundaries, Ekman friction is strong enough to prevent the development of a coherent dipole vortex condensate. Instead vortex pairs are found to be intermittent, varying in both time and strength.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Multiscale Dynamo Model Driven by Quasi-geostrophic Convection

A convection-driven multiscale dynamo model is developed in the limit of low Rossby number for the plane layer geometry in which the gravity and rotation vectors are aligned. The small-scale fluctuating dynamics are described by a magnetically-modified quasi-geostrophic equation set, and the large-scale mean dynamics are governed by a diagnostic thermal wind balance. The model utilizes three timescales that respectively characterize the convective timescale, the large-scale magnetic evolution timescale, and the large-scale thermal evolution timescale. Distinct equations are derived for the cases of order one and low magnetic Prandtl number. It is shown that the low magnetic Prandtl number model is characterized by a magnetic to kinetic energy ratio that is asymptotically large, with ohmic dissipation dominating viscous dissipation on the large-scales. For the order one magnetic Prandtl number model the magnetic and kinetic energies are equipartitioned and both ohmic and viscous dissipation are weak on the large-scales; large-scale ohmic dissipation occurs in thin magnetic boundary layers adjacent to the solid boundaries. For both magnetic Prandtl number cases the Elsasser number is small since the Lorentz force does not enter the leading order force balance. The new models can be considered fully nonlinear, generalized versions of the dynamo model originally developed by Childress and Soward [Phys. Rev. Lett., \textbf{29}, p.837, 1972]. These models may be useful for understanding the dynamics of convection-driven dynamos in regimes that are only just becoming accessible to direct numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exact coherent structures in a reduced model of parallel shear flow

A reduced description of shear flows consistent with the Reynolds number scaling of lower-branch exact coherent states in plane Couette flow [J. Wang et al., Phys. Rev. Lett. 98, 204501 (2007)] is constructed. Exact time-independent nonlinear solutions of the reduced equations corresponding to both lower and upper branch states are found for Waleffe flow [F. Waleffe, Phys. Fluids 9, 883--900 (1997)]. The lower branch solution is characterized by fluctuations that vary slowly along the critical layer while the upper branch solutions display a bimodal structure and are more strongly focused on the critical layer. The reduced model provides a rational framework for investigations of subcritical spatiotemporal patterns in parallel shear flows.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Reduced description of exact coherent states in parallel shear flows

Exact coherent states of a linearly stable, plane parallel shear flow confined between stationary stress-free walls and driven by a sinusoidal body force (a flow first introduced by F. Waleffe, Phys. Fluids 9, 883 (1997)) are computed using equations obtained from a large Reynolds-number asymptotic reduction of the Navier-Stokes equations. The reduced equations employ a decomposition into streamwise-averaged (mean) and streamwise-varying (fluctuation) components and are characterized by an effective order one Reynolds number in the mean equations along with a formally higher-order diffusive regularization of the fluctuation equations. A robust numerical algorithm for computing exact coherent states is introduced. Numerical continuation of the lower branch states to lower Reynolds numbers reveals the presence of a saddle-node; the saddle-node allows access to upper branch states that, like the lower branch states, appear to be self-consistently described by the reduced equations. Both lower and upper branch states are characterized in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Breakdown of the Anelastic Approximation in Rotating Compressible Convection: Implications for Astrophysical Systems

The linear theory for rotating compressible convection in a plane layer geometry is presented for the astrophysically-relevant case of low Prandtl number gases. When the rotation rate of the system is large, the flow remains geostrophically balanced for all stratification levels investigated and the classical (i.e., incompressible) asymptotic scaling laws for the critical parameters are recovered. For sufficiently small Prandtl numbers, increasing stratification tends to further destabilise the fluid layer, decrease the critical wavenumber and increase the oscillation frequency of the convective instability. In combination, these effects increase the relative magnitude of the time derivative of the density perturbation contained in the conservation of mass equation to non-negligible levels; the resulting convective instabilities occur in the form of compressional quasi-geostrophic oscillations. We find that the anelastic equations, which neglect this term, cannot capture these instabilities and possess spuriously-growing eigenmodes in the rapidly rotating, low Prandtl number regime. It is shown that the Mach number for rapidly rotating compressible convection is intrinsically small for all background states, regardless of the departure from adiabaticity.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Chaotic Advection and the Emergence of Tori in the Küppers-Lortz State

Motivated by the roll-switching behavior observed in rotating Rayleigh-Bénard convection, we define a Küppers-Lortz (K-L) state as a volume-preserving flow with periodic roll switching. For an individual roll state, the Lagrangian particle trajectories are periodic. In a system with roll-switching, the particles can exhibit three-dimensional, chaotic motion. We study a simple phenomenological map that models the Lagrangian dynamics in a K-L state. When the roll axes differ by $120^{\circ}$ in the plane of rotation, we show that the phase space is dominated by invariant tori if the ratio of switching time to roll turnover time is small. When this parameter approaches zero these tori limit onto the classical hexagonal convection patterns, and, as it gets large, the dynamics becomes fully chaotic and well-mixed. For intermediate values, there are interlinked toroidal and poloidal structures separated by chaotic regions. We also compute the exit time distributions and show that the unbounded chaotic orbits are normally diffusive. Although the map presumes instantaneous switching between roll states, we show that the qualitative features of the flow persist when the model has smooth, overlapping time-dependence for the roll amplitudes (the Busse-Heikes model).

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Baroclinic Vorticity Production in Protoplanetary Disks; Part I: Vortex Formation

The formation of vortices in protoplanetary disks is explored via pseudo-spectral numerical simulations of an anelastic-gas model. This model is a coupled set of equations for vorticity and temperature in two dimensions which includes baroclinic vorticity production and radiative cooling. Vortex formation is unambiguously shown to be caused by baroclinicity because (1) these simulations have zero initial perturbation vorticity and a nonzero initial temperature distribution; and (2) turning off the baroclinic term halts vortex formation, as shown by an immediate drop in kinetic energy and vorticity. Vortex strength increases with: larger background temperature gradients; warmer background temperatures; larger initial temperature perturbations; higher Reynolds number; and higher resolution. In the simulations presented here vortices form when the background temperatures are $\sim 200K$ and vary radially as $r^{-0.25}$, the initial vorticity perturbations are zero, the initial temperature perturbations are 5% of the background, and the Reynolds number is $10^9$. A sensitivity study consisting of 74 simulations showed that as resolution and Reynolds number increase, vortices can form with smaller initial temperature perturbations, lower background temperatures, and smaller background temperature gradients. For the parameter ranges of these simulations, the disk is shown to be convectively stable by the Solberg-Høiland criteria.

0 citations0 reviewsNo code linkedOpen access