Researcher profile

Louis-Alexandre Couston

Louis-Alexandre Couston contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 published item(s)

preprint2022arXiv

Competition between Rayleigh--Bénard and horizontal convection

We investigate the dynamics of a fluid layer subject to an imposed bottom heat flux and a top monotonically-increasing temperature profile driving horizontal convection. We use direct numerical simulations and consider a large range of flux-based Rayleigh numbers $10^6 \leq Ra_F \leq 10^9$ and imposed top horizontal to bottom vertical heat flux ratios $0 \leq Λ\leq 1$. The fluid domain is a closed two-dimensional box with aspect ratio $4\leq Γ\leq 16$ and we consider no-slip boundaries and adiabatic side walls. We demonstrate a regime transition from Rayleigh--Bénard convection (RB) to horizontal convection (HC) at $Λ\approx 10^{-2}$, which is independent of $Ra_F$ and $Γ$. At small $Λ$, the flow is organized in multiple overturning cells with approximately unit aspect ratio, while at large $Λ$ a single cell is obtained. The RB-relevant Nusselt number scaling with $Ra_F$ and the HC-relevant Nusselt number scaling with the horizontal Rayleigh number $Ra_L=Ra_FΛΓ^4$ are in good agreement with previous results from classical RB convection and HC studies in the limit $Λ\ll 10^{-2}$ and $Λ\gg 10^{-2}$, respectively. We demonstrate that the system is multi-stable near the transition $Λ\approx10^{-2}$, i.e. the exact number of cells not only depends on $Λ$ but also on the system's history. Our results suggest that subglacial lakes, which motivated this study, are likely to be dominated by RB convection, unless the slope of the ice-water interface, which controls the horizontal temperature gradient via the pressure-dependence of the freezing point, is greater than unity.

preprint2021arXiv

Dynamic flows create potentially habitable conditions in Antarctic subglacial lakes

Trapped beneath the Antarctic ice sheet lie over 400 subglacial lakes, which are considered to be extreme, isolated, yet viable habitats for microbial life. The physical conditions within subglacial lakes are critical to evaluating how and where life may best exist. Here, we propose that Earth&#39;s geothermal flux provides efficient stirring of Antarctic subglacial lake water. We demonstrate that most lakes are in a regime of vigorous turbulent vertical convection, enabling suspension of spherical particulates with diameters up to 36 micrometers. Thus, dynamic conditions support efficient mixing of nutrient- and oxygen-enriched meltwater derived from the overlying ice, which is essential for biome support within the water column. We caution that accreted ice analysis cannot always be used as a proxy for water sampling of lakes beneath a thin (<3.166 kilometers) ice cover, because a stable layer isolates the well-mixed bulk water from the ice-water interface where freezing may occur.

preprint2020arXiv

Aspect ratio affects iceberg melting

Iceberg meltwater is a critical freshwater flux from the cryosphere to the oceans. Global climate simulations therefore require simple and accurate parameterisations of iceberg melting. Iceberg shape is an important but often neglected aspect of iceberg melting. Icebergs have an enormous range of shapes and sizes, and distinct processes dominate basal and side melting. We show how different iceberg aspect ratios and relative ambient water velocities affect melting using a combined experimental and numerical study. The experimental results show significant variations in melting between different iceberg faces, as well as within each iceberg face. These findings are reproduced and explained with novel multiphysics numerical simulations. At high relative ambient velocities melting is largest on the side facing the flow, and mixing during vortex generation causes local increases in basal melt rates of over 50%. Double-diffusive buoyancy effects become significant when the relative ambient velocity is low. Existing melting parameterisations do not reproduce our findings. We propose several corrections to capture the influence of aspect ratio on iceberg melting.

preprint2020arXiv

Fabry-Perot resonance of water waves

We show that significant water wave amplification is obtained in a water resonator consisting of two spatially separated patches of small-amplitude sinusoidal corrugations on an otherwise flat seabed. The corrugations reflect the incident waves according to the so-called Bragg reflection mechanism, and the distance between the two sets controls whether the trapped reflected waves experience constructive or destructive interference within the resonator. The resulting amplification or suppression is enhanced with increasing number of ripples, and is most effective for specific resonator lengths and at the Bragg frequency, which is determined by the corrugation period. Our analysis draws on the analogous mechanism that occurs between two partially reflecting mirrors in optics, a phenomenon named after its discoverers Charles Fabry and Alfred Perot.

preprint2020arXiv

Landslide Tsunamis in Lakes

Landslides plunging into lakes and reservoirs can result in extreme wave runup at shores. This phenomenon has claimed lives and caused damage to near-shore properties. Landslide tsunamis in lakes are different from typical earthquake tsunamis in the open ocean in that (i) the affected areas are usually within the near-field of the source, (ii) the highest runup occurs within the time period of the geophysical event, and (iii) the enclosed geometry of a lake does not let the tsunami energy escape. To address the problem of transient landslide tsunami runup and to predict the resulting inundation, we utilize a nonlinear model equation in the Lagrangian frame of reference. The motivation for using such a scheme lies in the fact that the runup on an inclined boundary is directly and readily computed in the Lagrangian framework without the need to resort to approximations. In this work, we investigate the inundation patterns due to landslide tsunamis in a lake. We show by numerical computations that Airy&#39;s approximation of an irrotational theory using Lagrangian coordinates can legitimately predict runup of large amplitude. We also demonstrate that in a lake of finite size the highest runup may be magnified by constructive interference between edge-waves that are trapped along the shore and multiple reflections of outgoing waves from opposite shores, and may occur somewhat later after the first inundation.

preprint2020arXiv

Shape and size of large-scale vortices: A generic fluid pattern in geophysical fluid dynamics

Planetary rotation organizes fluid motions into coherent, long-lived swirls, known as large-scale vortices (LSVs), which play an important role in the dynamics and long-term evolution of geophysical and astrophysical fluids. Here, using direct numerical simulations, we show that LSVs in rapidly rotating mixed convective and stably stratified fluids, which approximates the two-layer, turbulent-stratified dynamics of many geophysical and astrophysical fluids, have a generic shape and that their size can be predicted. We show that LSVs emerge in the convection zone from upscale energy transfers and can penetrate into the stratified layer. At the convective-stratified interface, the LSV cores have a positive buoyancy anomaly. Due to the thermal wind constraint, this buoyancy anomaly leads to winds in the stratified layer that decay over a characteristic vertical length scale. Thus LSVs take the shape of a depth-invariant cylinder with a finite-size radius in the turbulent layer and of a penetrating half dome in the stratified layer. Importantly, we demonstrate that when LSVs penetrate all the way through the stratified layer and reach a boundary that is no-slip, they saturate by boundary friction. We provide a prediction for the penetration depth and maximum radius of LSVs as a function of the LSV vorticity, the stratified layer depth, and the stratification. Our results, which apply for cyclonic LSVs, suggest that LSVs in slowly rotating stars and Earth&#39;s liquid core are confined to the convective layer, while in Earth&#39;s atmosphere and oceans they can penetrate far into the stratified layer.

preprint2020arXiv

The Science Case for a Titan Flagship-class Orbiter with Probes

We outline a flagship-class mission concept focused on studying Titan as a global system, with particular emphasis on the polar regions. Investigating Titan from the unique standpoint of a polar orbit would enable comprehensive global maps to uncover the physics and chemistry of the atmosphere, and the topography and geophysical environment of the surface and subsurface. The mission includes two key elements: (1) an orbiter spacecraft, which also acts as a data relay, and (2) one or more small probes to directly investigate Titan&#39;s seas and make the first direct measurements of their liquid composition and physical environment. The orbiter would carry a sophisticated remote sensing payload, including a novel topographic lidar, a long-wavelength surface-penetrating radar, a sub-millimeter sounder for winds and for mesospheric/thermospheric composition, and a camera and near-infrared spectrometer. An instrument suite to analyze particles and fields would include a mass spectrometer to focus on the interactions between Titan&#39;s escaping upper atmosphere and the solar wind and Saturnian magnetosphere. The orbiter would enter a stable polar orbit around 1500 to 1800 km, from which vantage point it would make global maps of the atmosphere and surface. One or more probes, released from the orbiter, would investigate Titan&#39;s seas in situ, including possible differences in composition between higher and lower latitude seas, as well as the atmosphere during the parachute descent. The number of probes, as well as the instrument complement on the orbiter and probe, remain to be finalized during a mission study that we recommend to NASA as part of the NRC Decadal Survey for Planetary Science now underway, with the goal of an overall mission cost in the &#34;small flagship&#34; category of ~$2 bn. International partnerships, similar to Cassini-Huygens, may also be included for consideration.

preprint2020arXiv

Topography generation by melting and freezing in a turbulent shear flow

We report an idealized numerical study of a melting and freezing solid adjacent to a turbulent, buoyancy-affected shear flow, in order to improve our understanding of topography generation by phase changes in the environment. We use the phase-field method to dynamically couple the heat equation for the solid with the Navier-Stokes equations for the fluid. We investigate the evolution of an initially-flat and horizontal solid boundary overlying a pressure-driven turbulent flow. We assume a linear equation of state for the fluid and change the sign of the thermal expansion coefficient, such that the background density stratification is either stable, neutral or unstable. We find that channels aligned with the direction of the mean flow are generated spontaneously by phase changes at the fluid-solid interface. Streamwise vortices in the fluid, the interface topography and the temperature field in the solid influence each other and adjust until a statistical steady state is obtained. The crest-to-trough amplitude of the channels are larger than about 10$δ_ν$ in all cases, with $δ_ν$ the viscous length scale, but are much larger and more persistent for an unstable stratification than for a neutral or stable stratification. This happens because a stable stratification makes the cool melt fluid buoyant such that it shields the channel from further melting, whereas an unstable stratification makes the cool melt fluid sink, inducing further melting by rising hot plumes. The statistics of flow velocities and melt rates are investigated, and we find that channels and keels emerging in our simulations do not significantly change the mean drag coefficient.

preprint2019arXiv

Low-frequency variability in massive stars: Core generation or surface phenomenon?

Bowman et al. (2019) reported low-frequency photometric variability in 164 O- and B-type stars observed with K2 and TESS. They interpret these motions as internal gravity waves, which could be excited stochastically by convection in the cores of these stars. The detection of internal gravity waves in massive stars would help distinguish between massive stars with convective or radiative cores, determine core size, and would provide important constraints on massive star structure and evolution. In this work, we study the observational signature of internal gravity waves generated by core convection. We calculate the \textit{wave transfer function}, which links the internal gravity wave amplitude at the base of the radiative zone to the surface luminosity variation. This transfer function varies by many orders of magnitude for frequencies $\lesssim 1 \, {\rm d}^{-1}$, and has regularly-spaced peaks near $1 \, {\rm d}^{-1}$ due to standing modes. This is inconsistent with the observed spectra which have smooth ``red noise&#39;&#39; profiles, without the predicted regularly-spaced peaks. The wave transfer function is only meaningful if the waves stay predominately linear. We next show that this is the case: low frequency traveling waves do not break unless their luminosity exceeds the radiative luminosity of the star; and, the observed luminosity fluctuations at high frequencies are so small that standing modes would be stable to nonlinear instability. These simple calculations suggest that the observed low-frequency photometric variability in massive stars is not due to internal gravity waves generated in the core of these stars. We finish with a discussion of (sub)surface convection, which produces low-frequency variability in low-mass stars, very similar to that observed in Bowman et al. (2019) in higher mass stars.