2025-06-19 | | Total: 21
The present work consists of topics covered through a course currently taught by the author at SIMIS.
The shallow water equations often assume a constant velocity profile along the vertical axis. However, this assumption does not hold in many practical applications. To better approximate the vertical velocity distribution, models such as the shallow water moment expansion models have been proposed. Nevertheless, under non-slip bottom boundary conditions, both the standard shallow water equation and its moment-enhanced models struggle to accurately capture the vertical velocity profile due to the stiff source terms. In this work, we propose modified shallow water equations and corresponding moment-enhanced models that perform well under both non-slip and slip boundary conditions. The primary difference between the modified and original models lies in the treatment of the source term, which allows our modified moment expansion models to be readily generalized, while maintaining compatibility with our previous analysis on the hyperbolicity of the model. To assess the performance of both the standard and modified moment expansion models, we conduct a comprehensive numerical comparison with the incompressible Navier--Stokes equations -- a comparison that is absent from existing literature.
We present a novel artificial diffusion method to circumvent the instabilities associated with the standard finite element approximation of convection-diffusion equations. Motivated by the micromorphic approach, we introduce an auxiliary variable, which is related to the gradient of the field of interest, and which leads to a coupled problem. Conditions for well-posedness of the resulting formulation are established. We carry out a comprehensive numerical study to compare the proposed methodology against some well-established approaches in one- and two-dimensional settings. The proposed method outperforms established approaches in general in approximating accurately the solutions to pertinent and challenging problems.
In this article, we study the stability of solutions to a nonlinear viscoelastic plate problem with frictional damping of a memory on a part of the boundary, and a logarithmic source in a bounded domain $\Omega \subset \mathbb{R}^2.$ In this problem the relaxation function $r$ satisfies $r^{\prime}\left( t\right)\leq -h\left( t\right) \mathcal{G}\left( r\left( t\right) \right)$ for all $t\geq 0$, where $h$ is a nonincreasing positive function. This work extends previous works with viscoelastic plate problems and improves earlier results in the
We study coarse-graining methods for stochastic differential equations. In particular we consider averaging and a type of projection operator method, sometimes referred to as effective dynamic via conditional expectations. The projection method (PM) we consider is related to the ``mimicking marginals'' coarse graining approach proposed by Gyöngy. The first contribution of this paper is to provide further theoretical background for the PM and a rigorous link to the Gyöngy method. Moreover, we compare PM and averaging. While averaging applies to systems with time scale separation, the PM can in principle be applied irrespective of this. However it is often assumed that the two methods coincide in presence of scale separation. The second contribution of this paper is to make this statement precise, provide sufficient conditions under which these two methods coincide and then show -- via examples and counterexamples -- that this needs not be the case.
There exist elegant methods of aligning point clouds in $\mathbb R^3$. Unfortunately, these methods rely on the positive definite property of the Euclidean metric, and do not easily extend to the indefinite Minkowski metric. In this paper, we propose two solutions to the following problem: given inertial reference frames $A$ and $B$, and given (possibly noisy) measurements of a set of 4-vectors $\{v_i\}$ made in those reference frames with components $\{v_{A,i}\}$ and $\{v_{B,i}\}$, find the optimal Lorentz transformation $\Lambda$ such that $\Lambda v_{A,i}=v_{B,i}$. The method we outline is conceptually simple and easily extends to alignment problems in other matrix Lie groups.
Ion-acoustic waves in a dusty plasma are investigated where it is assumed that the ions follow a Cairns distribution and the electrons are Boltzmann distributed. Two theoretical methods are applied: Sagdeev pseudopotential analysis (SPA) and reductive perturbation theory (RPT). Since SPA incorporates all nonlinearities of the model it is the most accurate but deriving soliton profiles requires numerical integration of Poisson's equation. By contrast, RPT is a perturbation method which at second order yields the Gardner equation incorporating both the quadratic nonlinearity of the KdV equation with the cubic nonlinearity of the modified KdV equation. For consistency with the perturbation scheme the coefficient of the quadratic term needs to be at least an order of magnitude smaller than the coefficient of the cubic term. Solving the Gardner equation yields an analytic expression of the soliton profile. Selecting an appropriate set of compositional parameters, the soliton solutions obtained from SPA and RPT are analyzed and compared.
We relate graph complexes, Calabi-Yau $A_\infty$-categories and Kontsevich's cocycle construction. Our main result produces a commutative square of shifted Poisson algebras; one of its edges is the Loday-Quillen-Tsygan map, generalized to $A_\infty$-categories. We describe a quantized version via Beilinson-Drinfeld algebras. The larger context is to provide categorical methods which relate enumerative geometry (as in mirror symmetry) and large $N$ gauge theories.
We present an elementary approach to observe frequency cascade on forced nonlinear Schrödinger equations. The forcing term consists of a constant term, perturbed by a modulated Gaussian well. Algebraic computations provide an explicit frequency cascade when time and space derivatives are discarded from the nonlinear Schrödinger equation. We provide stability results, showing that when derivatives are incorporated in the model, the initial algebraic solution may be little affected, possibly over long time intervals. Numerical simulations are provided, which support the analysis.
This paper is a summary of discussions at the recent ITEP-JINR-YerPhI workshop on Vogel theory in Dubna. We consider relation between Vogel divisor(s) and the old Dynkin classification of simple Lie algebras. We consider application to knot theory and the hidden role of Jacobi identities in the definition/invariance of Kontsevich integral, which is the knot polynomial with the values in diagrams, capable of revealing all Vassiliev invariants -- including the ones, not visible in other approaches. Finally we comment on the possible breakdown of Vogel universality after the Jack/Macdonald deformation. Generalizations to affine, Yangian and DIM algebras are also mentioned. Especially interesting could be the search for universality in ordinary Yang-Mills theory and its interference with confinement phenomena.
We establish a unified theoretical framework that connects classical orthogonal polynomial systems to matrix Lyapunov equations through the fundamental physics of energy dissipation in stochastic dynamical systems. Starting from the energy balance principle in infinite-dimensional Hilbert spaces, we derive a master integral representation that naturally encompasses both spectral geometry and covariance dynamics. The theory reveals that established orthogonal polynomials (Zernike, Hermite, spherical harmonics) and matrix Lyapunov equations are dual manifestations of the same underlying energy dissipation structure. We provide rigorous mathematical foundations showing how finite-dimensional projections of infinite-dimensional energy integrals reproduce classical matrix equations, with specific structure determined by the symmetries of noise processes. The framework demonstrates that adding uniform dissipation to classical differential operators preserves their polynomial eigenfunction structure while ensuring the energy balance conditions required for physical consistency.
In the paper we study the impact of the boundary vorticity distribution on the dynamics of enstrophy for flows around streamlined body. A new energy identity is derived in the article, which includes the boundary values of the vortex function. For the Stokes system the dissipativity of enstrophy is proved. For the Navier-Stokes system a new equation of the enstrophy dynamics is obtained.
In this paper we consider a particular class of solutions of the Rayleigh-Boltzmann equation, known in the nonlinear setting as homoenergetic solutions, which have the form $g\left( x,v,t \right) =f\left( v-L\left( t\right)x,t\right)$ where the matrix $L(t)$ describes a shear flow deformation. We began this analysis in [22] where we rigorously proved the existence of a stationary non-equilibrium solution and established the different behaviour of the solutions for small and large values of the shear parameter, for cut-off collision kernels with homogeneity parameter $0\leq \gamma <1$, including Maxwell molecules and hard potentials. In this paper, we concentrate in the case where the deformation term dominates the collision term for large times (hyperbolic-dominated regime). This occurs for collision kernels with $\gamma < 0$ and in particular we focus on $\gamma \in (-1,0)$. In such a hyperbolic-dominated regime, it appears challenging to provide a clear description of the long-term asymptotics of the solutions. Here we present a formal analysis of the long-time asymptotics for the distribution of velocities and provide the explicit form for the asymptotic profile. Additionally, we discuss the different asymptotic behaviour expected in the case of homogeneity $\gamma < -1$. Furthermore, we provide a probabilistic interpretation describing a stochastic process consisting in a combination of collisions and shear flows. The tagged particle velocity $\{v(t)\}_{t\geq 0}$ is a Markov process that arises from the combination of free flights in a shear flow along with random jumps caused by collisions.
We study the discrete Bak-Sneppen model introduced by Barbay and Kenyon (2001) "On the discrete Bak-Sneppen model of self-organized criticality". We extend their results as well as the non-triviality result of Meester and Znamenskiy (2002) for a finite segment of $\mathbb{Z}^1$ with the periodic boundary condition to a large class of graphs, by using coupling between the Bak-Sneppen model and the oriented percolation in a quadrant. This allows us to avoid dealing with the so-called avalanches, thus simplifying many arguments.
Randomisation is widely used in quantum algorithms to reduce the number of quantum gates and ancillary qubits required. A range of randomised algorithms, including eigenstate property estimation by spectral filters, Hamiltonian simulation, and perturbative quantum simulation, though motivated and designed for different applications, share common features in the use of unitary decomposition and Hadamard-test-based implementation. In this work, we start by analysing the role of randomised linear-combination-of-unitaries (LCU) in quantum simulations, and present several quantum circuits that realise the randomised composite LCU. A caveat of randomisation, however, is that the resulting state cannot be deterministically prepared, which often takes an unphysical form $U \rho V^\dagger$ with unitaries $U$ and $V$. Therefore, randomised LCU algorithms are typically restricted to only estimating the expectation value of a single Pauli operator. To address this, we introduce a quantum instrument that can realise a non-completely-positive map, whose feature of frequent measurement and reset on the ancilla makes it particularly suitable in the fault-tolerant regime. We then show how to construct an unbiased estimator of the effective (unphysical) state $U \rho V^\dagger$ and its generalisation. Moreover, we demonstrate how to effectively realise the state prepared by applying an operator that admits a composite LCU form. Our results reveal a natural connection between randomised LCU algorithms and shadow tomography, thereby allowing simultaneous estimation of many observables efficiently. As a concrete example, we construct the estimators and present the simulation complexity for three use cases of randomised LCU in Hamiltonian simulation and eigenstate preparation tasks.
The work reported in ``Granular micromechanics-based identification of isotropic strain gradient parameters for elastic geometrically nonlinear deformations" misidentified key terms in the grain-pair objective relative displacement when accounting for the second gradient of placement. In this paper, we correct that oversight by deriving a revised expression for the grain-pair objective relative displacement within the granular micromechanics framework. The amended terms, which resemble Christoffel symbols expressed in terms of strain gradients, modify the contributions of both the normal and tangential components to the strain energy and, consequently, alter the identified strain gradient elastic parameters. Importantly, the identification of the standard (first gradient) elastic tensor remains unchanged. This brief paper presents the corrected derivation, the resulting stiffness tensors for anisotropic strain gradient elasticity, and updated analytical expressions for the material parameters in both 2D and 3D isotropic settings.
The classical Floquet theory allows to map a time-periodic system of linear differential equations into an autonomous one. By looking at it in a geometrical way, we extend the theory to a class of non-autonomous non-periodic equations. This is obtained by considering a change of variables which depends on time in a non-trivial way, i.e. introducing gauge transformations, well known in fundamental Physics and Field Theory -- but which seems to have received little attention in this context.
The challenges posed by the development of field theories, both classical and quantum, force us to question their most basic and foundational ideas like the role and origin of space-time, the meaning of physical states, etc. Among them the notion of ``field'' itself is notoriously difficult to address. These notes aim to analyze such notion from the perspective offered by the groupoid description of quantum mechanics inspired by Schwinger's picture of quantum mechanics. Then, a natural interpretation of the notion of physical fields as functors among appropriate groupoids will emerge. The domain of a field in this new picture is a groupoid that describes ``test particles'', and its codomain is a groupoid that describes the intrinsic nature of the system being probed. Such a space of functors carries some natural structures, which are best described in a categorical language. Some illustrative examples will be presented that could help clarify the various abstract notions discussed in the text.
We study the energy landscape of the Random Energy model (REM) integrated along trajectories of the simple random walk on the hypercube. We show that the quenched cumulant generating function of the time integral of the REM energy undergoes phase transitions in the large $N$ limit for trajectories of any time extent, and identify phases distinguished by the activity and value of the time integral. This is achieved by relating the dynamical behavior to the spectral properties of Hamiltonians associated with the Quantum Random Energy Model (QREM). Of independent interest are deterministic $ \ell^p $-properties of the resolvents of such Hamiltonians, which we establish.
We consider real tensors of order $D$, that is $D$-dimensional arrays of real numbers $T_{a^1a^2 \dots a^D}$, where each index $a^c$ can take $N$ values. The tensor entries $T_{a^1a^2 \dots a^D}$ have no symmetry properties under permutations of the indices. The invariant polynomials built out of the tensor entries are called trace invariants. We prove that for a Gaussian random tensor with $D\ge 3$ indices (that is such that the entries $T_{a^1a^2 \dots a^D}$ are independent identically distributed Gaussian random variables) the cumulant, or connected expectation, of a product of trace invariants is not always suppressed in scaling in $N$ with respect to the product of the expectations of the individual invariants. Said otherwise, not all the multi-trace expectations factor at large $N$ in terms of the single-trace ones and the Gaussian scaling is not subadditive on the connected components. This is in stark contrast to the $D=2$ case of random matrices in which the multi-trace expectations always factor at large $N$. The best one can do for $D\ge 3$ is to identify restricted families of invariants for which the large $N$ factorization holds and we check that this indeed happens when restricting to the family of melonic observables, the dominant family in the large $N$ limit.
This is an introduction to measure theory, integration and function spaces, with all the needed preliminaries included, and with some applications included as well. We first discuss some basic motivations, coming from discrete probability, that we develop in detail, as a preliminary to general measure theory. Then we discuss measure theory, integration and function spaces, all developed in a standard way, and with emphasis on the explicit computation of various integrals. Finally, we come back to probability, discrete and continuous, with a more advanced discussion, of quantum flavor.