Researcher profile

Rupert Klein

Rupert Klein contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
12topics
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)

preprint2025arXiv

Affine Invariant Langevin Dynamics for rare-event sampling

We introduce an affine invariant Langevin dynamics (ALDI) framework for the efficient estimation of rare events in nonlinear dynamical systems. Rare events are formulated as Bayesian inverse problems through a nonsmooth limit-state function whose zero level set characterises the event of interest. To overcome the nondifferentiability of this function, we propose a smooth approximation that preserves the failure set and yields a posterior distribution satisfying the small-noise limit. The resulting potential is sampled by ALDI, a (derivative-free) interacting particle system whose affine invariance allows it to adapt to the local anisotropy of the posterior. We demonstrate the performance of the method across a hierarchy of benchmarks, namely two low-dimensional examples (an algebraic problem with convex geometry and a dynamical problem of saddle-type instability) and a point-vortex model for atmospheric blockings. In all cases, ALDI concentrates near the relevant near-critical sets and provides accurate proposal distributions for self-normalised importance sampling. The framework is computationally robust, potentially gradient-free, and well-suited for complex forward models with strong geometric anisotropy. These results highlight ALDI as a promising tool for rare-event estimation in unstable regimes of dynamical systems.

preprint2023arXiv

On the boundary layer arising from fast internal waves dynamics

In this paper, we investigate the boundary layer arising from the fast internal waves in the Boussinesq equations with the Brunt-Vaisälä frequency of order $ \mathcal O(1/\varepsilon) $. For the homogeneous-in-height stratification, previous work by \emph{Desjardins, Lannes, Saut, 3(1):153--192, Water Waves, 2021} establishes uniform-in-$ε$ estimates locally in time, with additional constraints on the boundary data initially, which essentially restricts the dynamics in the spatially periodic domain. Removing such constraints, our goal is to investigate the general near-boundary behavior. We observe that the fast internal waves will give rise to large growth of the spatial derivatives in the normal direction of the solutions in the vicinity of the boundary. To capture this phenomenon, we introduce an inviscid boundary layer using a natural scaling. In addition, we investigate the well-posedness of such a boundary layer system in the space of analytic functions.

preprint2022arXiv

Chemical diffusion master equation: formulations of reaction--diffusion processes on the molecular level

The chemical diffusion master equation (CDME) describes the probabilistic dynamics of reaction--diffusion systems at the molecular level [del Razo et al., Lett. Math. Phys. 112:49, 2022]; it can be considered the master equation for reaction--diffusion processes. The CDME consists of an infinite ordered family of Fokker--Planck equations, where each level of the ordered family corresponds to a certain number of particles and each particle represents a molecule. The equations at each level describe the spatial diffusion of the corresponding set of particles, and they are coupled to each other via reaction operators --linear operators representing chemical reactions. These operators change the number of particles in the system, and thus transport probability between different levels in the family. In this work, we present three approaches to formulate the CDME and show the relations between them. We further deduce the non-trivial combinatorial factors contained in the reaction operators, and we elucidate the relation to the original formulation of the CDME, which is based on creation and annihilation operators acting on many-particle probability density functions. Finally we discuss applications to multiscale simulations of biochemical systems among other future prospects.

preprint2022arXiv

The soundproof model of an acoustic--internal waves system with low stratification

This work is devoted to investigating a compressible fluid system with low stratification, which is driven by fast acoustic waves and internal waves. The approximation using a soundproof model is justified. More precisely, the soundproof model captures the dynamics of both the non-oscillating mean flows and the oscillating internal waves, while filters out the fast acoustic waves, of the compressible system with or without initial acoustic waves. Moreover, the fast-slow oscillation structure is investigated.

preprint2021arXiv

QG-DL: Dynamics of a diabatic layer in the quasi-geostrophic framework

Quasi-geostrophic (QG) theory describes the dynamics of synoptic scale flows in the trophosphere that are balanced with respect to both acoustic and internal gravity waves. Within this framework, effects of (turbulent) friction near the ground are usually represented by Ekman Layer theory. The troposphere covers roughly the lowest ten kilometers of the atmosphere while Ekman layer heights are typically just a few hundred meters. However, this two-layer asymptotic theory does not explicitly account for substantial changes of the potential temperature stratification due to diabatic heating associated with cloud formation or with radiative and turbulent heat fluxes, which, in the middle latitudes, can be particularly important in about the lowest three kilometers. To address this deficiency, this paper extends the classical QG-Ekman layer model by introducing an intermediate, dynamically and thermodynamically active layer, called the "diabatic layer" (DL) from here on. The flow in this layer is also in acoustic, hydrostatic, and geostrophic balance but, in contrast to QG flow, variations of potential temperature are not restricted to small deviations from a stable and time independent background stratification. Instead, within the diabatic layer, diabatic processes are allowed to affect the leading-order stratification. As a consequence, the diabatic layer modifies the pressure field at the top of the Ekman layer, and with it the intensity of Ekman pumping seen by the quasi-geostrophic bulk flow. The result is the proposed extended quasi-geostrophic three-layer QG-DL-Ekman model for mid-latitude (dry and moist) dynamics.

preprint2020arXiv

Asymptotics for moist deep convection I: Refined scalings and self-sustaining updrafts

Moist processes are among the most important drivers of atmospheric dynamics,and scale analysis and asymptotics are cornerstones of theoretical meteorology. Accounting for moist processes in systematic scale analyses therefore seems of considerable importance for the field. Klein & Majda (TCFD, 20, 525--552, (2006)) proposed a scaling regime for the incorporation of moist bulk microphysics closures in multiscale asymptotic analyses of tropical deep convection. This regime is refined here to allow for mixtures of ideal gases and to establish consistency with a more general multiple scales modelling framework for atmospheric flows. Deep narrow updrafts, so-called "hot towers", constitute principal building blocks of larger scale storm systems. They are analysed here in a sample application of the new scaling regime. A single quasi-onedimensional columnar cloud is considered on the vertical advective (or tower life cycle) time scale. The refined asymptotic scaling regime is essential for this example as it reveals a new mechanism for the self-sustainance of such updrafts. Even for strongly positive convectively available potential energy (CAPE), a vertical balance of buoyancy forces is found in the presence of precipitation. This balance induces a diagnostic equation for the vertical velocity and it is responsible for the generation of self-sustained balanced updrafts. The time dependent updraft structure is encoded in a Hamilton-Jacobi equation for the precipitation mixing ratio. Numerical solutions of this equation suggest that the self-sustained updrafts may strongly enhance hot tower life cycles.

preprint2020arXiv

Global Well-posedness for the Primitive Equations Coupled to Nonlinear Moisture Dynamics with Phase Changes

In this work we study the global solvability of the primitive equations for the atmosphere coupled to moisture dynamics with phase changes for warm clouds, where water is present in the form of water vapor and in the liquid state as cloud water and rain water. This moisture model contains closures for the phase changes condensation and evaporation, as well as the processes of autoconversion of cloud water into rainwater and the collection of cloud water by the falling rain droplets. It has been used by Klein and Majda in \cite{KM} and corresponds to a basic form of the bulk microphysics closure in the spirit of Kessler \cite{Ke} and Grabowski and Smolarkiewicz \cite{GS}. The moisture balances are strongly coupled to the thermodynamic equation via the latent heat associated to the phase changes. In \cite{HKLT} we assumed the velocity field to be given and proved rigorously the global existence and uniqueness of uniformly bounded solutions of the moisture balances coupled to the thermodynamic equation. In this paper we present the solvability of a full moist atmospheric flow model, where the moisture model is coupled to the primitive equations of atmospherical dynamics governing the velocity field. For the derivation of a priori estimates for the velocity field we thereby use the ideas of Cao and Titi \cite{CT}, who succeeded in proving the global solvability of the primitive equations.

preprint2020arXiv

Open Systems out of Equilibrium: Theory and Simulation

We consider the theoretical model of Bergmann and Lebowitz for open systems out of equilibrium and translate its principles in the adaptive resolution molecular dynamics technique (AdResS). We simulate Lennard-Jones fluids with open boundaries in a thermal gradient and find excellent agreement of the stationary responses with results obtained from the simulation of a larger, locally forced closed system. The encouraging results pave the way for a computational treatment of open systems far from equilibrium framed in a well-established theoretical model that avoids possible numerical artifacts and physical misinterpretations.

preprint2019arXiv

A semi-implicit compressible model for atmospheric flows with seamless access to soundproof and hydrostatic dynamics

We introduce a second-order numerical scheme for compressible atmospheric motions at small to planetary scales. The collocated finite volume method treats the advection of mass, momentum, and mass-weighted potential temperature in conservation form while relying on Exner pressure for the pressure gradient term. It discretises the rotating compressible equations by evolving full variables rather than perturbations around a background state, and operates with time steps constrained by the advection speed only. Perturbation variables are only used as auxiliary quantities in the formulation of the elliptic problem. Borrowing ideas on forward-in-time differencing, the algorithm reframes the authors' previously proposed schemes into a sequence of implicit midpoint, advection, and implicit trapezoidal steps that allows for a time integration unconstrained by the internal gravity wave speed. Compared with existing approaches, results on a range of benchmarks of nonhydrostatic- and hydrostatic-scale dynamics are competitive. The test suite includes a new planetary-scale inertia-gravity wave test highlighting the properties of the scheme and its large time step capabilities. In the hydrostatic-scale cases the model is run in pseudo-incompressible and hydrostatic mode with simple switching within a uniform discretization framework. The differences with the compressible runs return expected relative magnitudes. By providing seamless access to soundproof and hydrostatic dynamics, the developments represent a necessary step towards an all-scale blended multimodel solver.