Scale Perturbation Research Articles

Permafrost modelling with OpenFOAM®: New advancements of the permaFoam solver

Permafrost, i.e., soil that is year-round frozen in depth, is covering a quarter of the continents of the northern hemisphere. It currently experiences fast changes due to climate change at global scale and technogenic perturbations at local scale, and the assessment and anticipation of these changes are of primary importance for many environmental and engineering applications in cold regions. To these ends, permafrost modelling is required, while the strong couplings and non-linearities involved in the physics at stake make it highly challenging, especially from a computational point of view. In this work we present a new solver for permafrost hydrology developed in the framework of OpenFOAM®, allowing to benefit from its advanced high-performance computing capabilities. The solver is tested for realistic, field-based cases, and its parallel performances are characterized up to ∼16000 cores on IRENE supercomputer (TGCC, CEA). Program summaryProgram Title: permaFoamCPC Library link to program files:https://doi.org/10.17632/swp88cvpwb.1Developer's repository link:https://develop.openfoam.com/Community/hydrology/Licensing provisions: GPLv3Programming language: C++Nature of problem: This software solves the coupled equations that govern water flow and heat transfer in variably saturated and variably frozen porous media, for transient problems in Three-dimensional, heterogeneous domains. The equation for water flow is Richards equation, which is a very popular model for water transfer in variably saturated porous media (e.g.: soils), and the equation for heat transfer is a Fourier equation including advection and the freeze/thaw of the pore water. The solver is designed to take advantage of the massively parallel computing performance of OpenFOAM®. The goal is to be able to model natural hydrosystems of cold regions on large temporal and spatial scales.Solution method: For each equation a mixed implicit (FVM for Finite Volume Method in the object oriented OpenFOAM framework) and explicit (FVC for Finite Volume Calculus in the object oriented OpenFOAM® framework) discretization with backward time scheme is embedded in an iterative linearization procedure (Picard algorithm). The coupling between the two equations is performed through an operator splitting approach. The implementation has been carried out with a concern for robustness and parallel efficiency.Additional comments including restrictions and unusual features: This version of permaFoam has been tested with OpenFOAM_v1912, v2106, v2112 and v2206, thus everything might not work with other (especially older) versions of OpenFOAM®. When using permaFoam, one should be careful to use fine enough spatial and temporal discretisations where and when steep fronts (freeze/thaw fronts, imbibition/drainage fronts) occur, otherwise numerical stability problems might arise.

Operator complexity for quantum scalar fields and cosmological perturbations

We calculate the operator complexity for the displacement, squeeze and rotation operators of a quantum harmonic oscillator, assuming equal computational cost for the corresponding fundamental gates. The complexity of the time-dependent displacement operator is constant, equal to the magnitude of the coherent state parameter, while the complexity of unitary evolution by a generic quadratic Hamiltonian is proportional to the amount of squeezing and is sensitive to the time-dependent phase of the unitary operator. We apply these results to study the complexity of a free massive scalar field, finding that the complexity has a period of rapid linear growth followed by a saturation determined by the UV cutoff and the number of spatial dimensions. We also study the complexity of the unitary evolution of quantum cosmological perturbations in de Sitter space, which can be written as time-dependent squeezing and rotation operators on individual Fourier mode pairs. The complexity of a single mode pair at late times grows linearly with the number of $e$-folds, while the complexity at early times oscillates rapidly due to the sensitivity of operator complexity to the phase of unitary time evolution. Integrating over all modes, the total complexity of cosmological perturbations scales as the square root of the (exponentially) growing volume of de Sitter space, suggesting that inflation leads to an explosive growth in complexity of the Universe.

Propagating Alfvén waves in open structures with random structuring

ABSTRACT We consider the behaviour of Alfvén waves propagating in a medium with random density perturbations. The imposed density perturbations have a broad-band spectrum and their characteristic spatial scale may be defined according to the peak in the spectrum. The interaction of the boundary driven Alfvén waves with the medium generates reflections most efficiently when their wavelength is comparable to the spatial scale of the density perturbations. For our monotonic driver, this leads to the generation of quasi-periodic oscillations. The periods of oscillation of the propagating Alfvén waves is no longer only associated with the driver. Additional periodicities may be associated with one or more characteristic spatial scales in the density profile, or with beating between other spectral components. Multiple wave reflections cause oscillatory power to be retained at low altitudes, increasing opportunities to contribute to heating at those locations.

Effective heat conductivity of a composite with hexagonal lattice of perfectly conducting circular inclusions: An analytical solution

AbstractPaper is devoted to the effective stationary heat conductivity for the fibre composite materials. We are aimed on getting on analytical expression for effective thermal conductivity coefficient. Asymptotic homogenization approach, based on the multiple scale perturbation method, is used. This allows to reduce the original boundary value problem in multiply connected domain to the sequence of boundary value problems in simply connected domains. These problems include: the local problem for the periodically repeated cell and homogenized problem with effective coefficient. It is shown that for densely packed and high contrast fibre composites, the cell problem can be solved analytically. For this aim, lubrication approach (asymptotics of thin layer) has been employed. We also generalise obtained solution to the case of medium‐sized inclusions in the framework of the Padé approximants.

A nonlinear time delay control influence in simulation aircraft tail contain a vertical cantilever beam

Nonlinear time delay saturation controller (NTDSC) is used to minimize the vibrations of a vertical cantilever beam simulated as an airplane tail within harmonic excitation force. In the case of primary and 1:2 internal resonance as a simultaneous resonance case, an analytic solution is obtained by employing the multiple time scale perturbation technique to solve the nonlinear differential equations and modeling the system with NTDSC. Time histories and Poincare portraits were created to demonstrate the system without and with controller. The time-history response, as well as the impacts of the parameters on the system and controller, were simulated numerically using the MATLAB program. A good match is found when the approximate solutions are compared to the numerical simulations via the Runge–Kutta method fourth-order (RK4). The stability investigation of the steady state solution in the studied simultaneous resonance case is reviewed and analyzed via the Routh-Huriwitz approach. The impact of controller parameters and time delays on system response curves are examined. Also, the operation’s safe zone is established from the outcomes of the effect via various values of time delays on the system after adding the control scheme. The numerical and analytical solutions at time-history were compared using the MATLAB program to confirm their comparability. The outcomes of this study were used to develop a theoretical foundation for the system’s design and construction. A comparison is made with recently released papers to appear the difference between this study and the others.

Parton physics of the large- Nc mesons

We initiate the studies on the structural physics of the tower of stable large-$N_c$ mesons through a first computation of the collinear quark-structure of a large-$N_c$ pion using lattice Monte-Carlo methods. We adapt the large-$N_c$ continuum reduction for the determination of meson correlation functions involving the spatially-extended quasi-PDF operators as a perfect strategy to concentrate only on the short perturbative length scales. We find the internal structures of pion in the large-$N_c$ and $N_c=3$ theories to be quite similar. Interestingly, we find hints that even the observed differences could arise to a large extent via the different perturbative QCD evolution in the two theories from similar initial conditions at low factorization scales.

Turbulence generation by shock interaction with a highly nonuniform medium.

An initially planar shock wave propagating into a medium of nonuniform density will be perturbed, leading to the generation of postshock velocity perturbations. Using numerical simulations we study this phenomenon in the case of highly nonuniform density (order-unity normalized variance, σ_{ρ}/ρ[over ¯]∼1) and strong shocks (shock Mach numbers M[over ¯]_{s}≳10). This leads to a highly disrupted shock and a turbulent postshock flow. We simulate this interaction for a range of shock drives and initial density configurations meant to mimic those which might be presently achieved in experiments. Theoretical considerations lead to scaling relations, which are found to reasonably predict the postshock turbulence properties. The turbulent velocity dispersion and turbulent Mach number are found to depend on the preshock density dispersion and shock speed in a manner consistent with the linear Richtymer-Meshkov instability prediction. We also show a dependence of the turbulence generation on the scale of density perturbations. The postshock pressure and density, which can be substantially reduced relative to the unperturbed case, are found to be reasonably predicted by a simplified analysis that treats the extended shock transition region as a single normal shock.

Abstract 197: A systemic model derivation platform for generating 3D neuroendocrine tumor cell spheroids to accelerate cancer research

Abstract The Broad Cancer Cell Line Factory (CCLF) aims to increase the number and representation of in vitro/ex vivo cell models for common and rare cancer types. Neuroendocrine tumor (NET) cell model derivation is one of the CCLFs focus because it lacks well-characterized, publicly available models.Two major barriers existed in deriving NET cell models, including how to collect sufficient patient tumor tissue samples for the model derivation pilot and how to systematically iterate model derivation strategies since there is no prior knowledge for NET cell model generation success. Thus, we partnered with the MD Anderson Cancer Center, the Dana-Farber Cancer Institute, and the Rare Cancer Research Foundation to collect patient tissue samples. All patient’s NET tissues were sequenced with a targeted Pan-Cancer panel to ensure high tumor content. To reduce fibroblast outgrowth, we combined an empirical rich media matrix (HYBRID technology) with a 3D spheroid culture system to initiate one sample in 16-64 conditions. The growing cultures at passage 3-5 were genomically credentialed to ensure the driver events matched with the original patient tissue. So far, we have received more than 70 NET samples. While several derived models are still under culture, we successfully generated 5 genomically verified NET tumor models, including small intestinal, pancreas, and liver subtypes. To phenotypically characterize these NET models, neuroendocrine biomarkers such as chromogranin A, synaptophysin, SSTR2, and VMAT 1/2 were also evaluated using qRT-PCR and ELISA. We observed that these NET spheroid models display long doubling times (2-4 weeks) at later passages which limits their utility for large scale perturbation experiments and model sharing capability with the research community. While we are currently working on several strategies to improve the propagation ability in these models, 1 (out of 5) model has reached passage 15 with a 3 day doubling time. Genomic studies, such as RNAseq, will be performed to address the model transcriptome changes after overcoming growth plateaus. Here we showed that it is feasible to derive NET models from patient biospecimens using our HYBRID strategy. As we expand our NET cohort, we will further refine disease-specific model generation protocols for different NET types. Our goal is to share our model generation experience and make these tumor cell models publicly available to the research community in order to accelerate cancer research. Citation Format: Adel Attari, Madison Liistro, Barbara Van Hare, Jennifer Chan, Emma Coleman, Tim Heffernan, Bianca Amador, Matthew Meyerson, Jesse Boehm, William Sellers, Yuen-Yi (Moony) Tseng. A systemic model derivation platform for generating 3D neuroendocrine tumor cell spheroids to accelerate cancer research [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 197.

The golden number seen in a mechanical oscillator

A seemingly ubiquitous irrational number often appearing in nature and in man-made things like structures, paintings, and physical systems, is the golden number. Here, we show that this astonishing number appears in the periodic solutions of an underactuated mass-spring oscillator driven by a nonlinear self-excitation. Specifically, by using the two-time scale perturbation method, it is analytically demonstrated that the golden number appears in the ratio of amplitudes, as well as in the oscillation frequency of the periodic solution, which is referred to as golden solution and, by applying the Poincaré method, it is demonstrated that this solution is asymptotically stable. Additionally, the analytic results are illustrated by means of numerical simulations and also, an experimental study is conducted.

Eliminate the Nonlinear Oscillations of the Modified duffing Equation by using the Nonlinear Integrated Positive Position Feedback

In this article, the nonlinear integrated positive position feedback (NIPPF) control adds to a nonlinear dynamical system modeled as the well known Duffng oscillators. This control is proposed to mitigate system nonlinear vibrations. The whole system mathematical model is analyzed by applying the multiple time scales perturbation method. The slow- ow modulation equations that govern the oscillation amplitudes of both the main system and controller are derived. The stability of the steady-state solution is presented and studied applying frequency response equations near the simultaneous primary and internal resonance cases. Before and after (NIPPF) control the nonlinear systems' steady-state amplitude are examined, the comparison is made to validate the closeness between the numerical solution and the analytical perturbative one at time-history and frequency response curves.

Quantum and Relativistic Effects on the KdV and Envelope Solitons in Ion-Plasma Waves

In this article, we investigate the linear and nonlinear behavior of ion-plasma waves in a two-component plasma containing ions and dust particles. We have studied the linear dispersion characteristics, nonlinear KdV solitons, and amplitude modulated electrostatic waves. Such ion-plasma wave modes in dust-ion plasma have been theoretically investigated by incorporating quantum diffraction and relativistic effects. Numerical analysis of the linear dispersion relation with variations of different parameters has been carried out. By employing the reductive perturbation technique the KdV equation describing the small amplitude solitary profile has been studied. Using the standard multiple scale perturbation method, a nonlinear Schrödinger (NLS) equation has been derived by including quantum and relativistic effects, which describes the nonlinear evolution of the ion-plasma waves in dust-ion plasma. The wave is found to become modulationally unstable beyond certain critical wavenumber. The critical wavenumber and the growth rate of modulational instability are shown to depend significantly on the relativistic motion of plasma particles and the quantum diffraction effect. The results presented in the article are expected to be helpful for understanding the propagation of finite amplitude ion-plasma waves in some laboratory and astrophysical plasma systems. <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">PACS</i> —52.59.Hq, 67.10.Jn, 52.20 −j, 52.25 −b, 52.30 −d, 52.35 Mw, 52.65 Vv

Gravitational Wave Signatures of Black Hole Quasinormal Mode Instability.

Black hole (BH) spectroscopy has emerged as a powerful approach to extracting spacetime information from gravitational wave (GW) observed signals. Yet, quasinormal mode (QNM) spectral instability under small scale perturbations has been recently shown to be a common classical general relativistic phenomenon [J. L. Jaramillo et al., Phys. Rev. X 11, 031003 (2021)PRXHAE2160-330810.1103/PhysRevX.11.031003]. This requires assessing its impact on the BH QNM spectrum, in particular on BH QNM overtone frequencies. We conclude (i)perturbed BH QNM overtones are indeed potentially observable in the GW waveform, providing information on small-scale environment BH physics, and (ii)their detection poses a challenging data analysis problem of singular interest for LISA astrophysics. We adopt a twofold approach, combining theoretical results from scattering theory with a fine-tuned data analysis on a highly accurate numerical GW ringdown signal. The former introduces a set of effective parameters (partially relying on a BH Weyl law) to characterize QNM instability physics. The latter provides a proof of principle demonstrating that the QNM spectral instability is indeed accessible in the time-domain GW waveform, though certainly requiring large signal-to-noise ratios. Particular attention is devoted to discussing the patterns of isospectrality loss under QNM instability, since the disentanglement between axial and polar GW parities may already occur within the near-future detection range.

Global sensitivity analysis in hydrodynamic modeling and flood inundation mapping

Flood simulation with two-dimensional hydrodynamic models is subject to different sources of uncertainties in model configuration, boundary conditions, and model parametrization. HEC-RAS 2D is a widely used hydrodynamic model. Here, we assess the sensitivity of the HEC-RAS 2D to its configuration factors and parameters. We evaluate the impacts of different model configuration factors, including the floodplain and channel roughness coefficients, terrain and mesh size, as well as river boundary conditions on the dynamic of water levels, maximum water level, and flood extent, and determine the importance of these contributing factors for reliable flood inundation modeling using both variogram and variance-based analyses. For our case study, we found that depending on the perturbation scale, the performance measure, and the predicted output, the sensitivity to input factors changes. While for simulation of water level dynamics, the DEM and mesh resolutions are the most important factors, for flood extent mapping the floodplain roughness coefficient and the upstream boundary condition are the key factors. In addition, our analyses indicated that for a reliable simulation of maximum flood water level, modelers need to spend resources on the calibration of floodplain roughness coefficient while using fine DEM and mesh resolutions. Here, we also investigate the role of computational time interval and mesh resolution on the model's run time and accuracy. Our results suggested that reducing the computational time interval has minimum impact, as it increases the model run time without much improving the model accuracy.

Impacts of filler loading and particle size on the transition to linear-nonlinear dichotomy in the rheological responses of particle-filled polymer solutions

In this study, the nonlinear behavior of carbon black-filled polybutadiene solutions under large-amplitude oscillatory shear is investigated. The results show that in the nonlinear regime, the third harmonic intensity, as measured by the ratio of the third to the first harmonics I3/I1, decreases significantly above a critical concentration ϕc of the polymer in the matrix, which results in the amplitude stress deviating strongly from the linear dependence of strain, while the time dependence of stress remains sinusoidal. Increasing the filler particle size significantly decreases the critical ϕc. However, increasing the filler loading basically has no effect on the transition to linear-nonlinear dichotomy. This transition happens when the mesh size ξ of the entangled polymer network in the matrix becomes smaller than the primary filler particle size. Above ϕc, the topological hindrance of the entangled polymer chains apparently considerably slows down the recovery speed of the broken filler network in the material. Hence, the quasisinusoidal response in the system that has a strain-dependent modulus is probably due to the restoration of the broken filler network requiring longer than the time scale of a typical dynamic perturbation.

The Beaufort Sea acoustic duct's variability and its impact on acoustic propagation using the mode interaction parameter

The Beaufort duct is a subsurface sound channel formed by cold Pacific Winter Water sandwiched between warmer Pacific Summer Water and Atlantic Water in the Western Arctic Ocean. This duct traps sound waves and allows them to travel long distances without losing energy to lossy interactions with sea ice and surface waves. This study quantifies Beaufort duct variability based on Canada Basin Acoustic Propagation Experiment (CANAPE) and Coordinated Arctic Acoustic Thermometry Experiment (CAATEX) oceanographic observations. Deterministic ocean features induce coupling between acoustic modes confined to the Beaufort duct and non-ducted modes by weakening the duct or causing it to take on an asymmetric form. A non-dimensional mode interaction parameter (MIP) can be defined using the acoustic frequency and the vertical and horizontal scales of sound speed perturbation to characterize coupling strength. It identifies three important wave propagation regimes: sudden approximation, adiabatic approximation, and maximum interaction regime (Colosi and Zinicola-Lapin, 2021). When the MIP between ducted and lossy modes is more and less than 1, strong and weak acoustic variability is predicted, respectively. Variability is high when the MIP between two ducted modes surpasses 1, but modei–mode interference patterns grow more complicated. Acoustic numerical simulations are used to demonstrate various effects.

Dispersion in lattices with patterns of hardening and softening stiffness nonlinearity

Recent focus has been given to analytically predicting amplitude-dependent dispersion in periodic structures with weak stiffness nonlinearity. These dispersion relationships may inform devices such as amplitude-dependent filters, waveguides, and diodes. However, attention is generally restricted to either spatially uniform nonlinearity or mutually exclusive distributions of hardening or softening stiffness. This study investigates dispersion in lattices with spatial modulations in stiffness nonlinearity. Special attention is given to modulations consisting of both hardening and softening nonlinear terms. A multiple scales perturbation analysis reveals that patterns of hardening and softening stiffness enable both lifting and lowering of a passband’s frequencies. In such cases, passbands contain discrete frequency and wavenumber pairs that prevent amplitude-dependent dispersion shifting. Numerical integration of the lattice equations of motion are carried out to confirm the analytically predicted dispersion behavior. A design of experiment is proposed in which strength and sign of nonlinearity can be tuned with the initial angle of additively manufactured grounding springs. [Work supported by the Office of Naval Research.]

Vorticity Budget and Formation Mechanisms of a Mesoscale Convective Vortex in a Heavy-Rainstorm Episode

Mesoscale convective vortices (MCVs) often cause rainstorms. To deepen our understanding of MCV formation mechanisms, reanalysis data from the National Centers for Environmental Prediction and the Weather Research and Forecasting model were used to simulate MCV activity in East China in August 2009. The simulations could reproduce the MCV and associated convective activities well. The vorticity budget and MCV formation mechanisms were then analyzed. The results show that the planetary vorticity advection is much smaller than other terms of the vorticity equation. The MCV initiates in the convective precipitation region and below 800 hPa. When the MCV initiates, there are vorticity-variation couplets within the vortex, and the MCV moves towards the positive vorticity-variation direction. In positive vorticity-variation areas, the divergence term and the tilting term are the vorticity source. The equilibrium response to diabatic heating is one of the forming mechanisms of this MCV. The latent-heating level is relatively low in this MCV case, and the MCV-forming level is also relatively low. Another forming mechanism of this MCV is the tilting of the horizontal vortex tube caused by the upward motion. At the MCV initiation, the perturbation scale of the vortex is found to be larger than the Rossby deformation radius, and thus the MCV could have a long duration.

The effect of stochastically perturbed parametrisation tendencies (SPPT) on rapidly ascending air streams

AbstractThe stochastically perturbed parametrisation tendency (SPPT) scheme is a well‐established technique in ensemble forecasting to address model uncertainty by introducing perturbations into the tendencies provided by the physics parametrisations. The magnitude of the perturbations scales with the local net parametrisation tendency, resulting in large perturbations where diabatic processes are active. Rapidly ascending air streams, such as warm conveyor belts (WCBs) and organized tropical convection, are often driven by cloud diabatic processes and are therefore prone to such perturbations. This study investigates the effects of SPPT and initial condition perturbations on rapidly ascending air streams by computing trajectories in sensitivity experiments with the European Centre for Medium‐Range Weather Forecasts (ECMWF) ensemble prediction system, which are set up to disentangle the effects of initial conditions and physics perturbations. The results demonstrate that SPPT systematically increases the frequency of rapidly ascending air streams. The effect is observed globally, but is enhanced in regions where the latent heating along the trajectories is larger. Despite the frequency changes, there are only minor modifications to the physical properties of the trajectories due to SPPT. In contrast to SPPT, initial condition perturbations do not affect WCBs and tropical convection systematically. An Eulerian perspective on vertical velocities reveals that SPPT increases the frequency of strong upward motions compared with experiments with unperturbed model physics. Consistent with the altered vertical motions, precipitation rates are also affected by the model physics perturbations. The unperturbed control member shows the same characteristics as the experiments without SPPT regarding rapidly ascending air streams. We make use of this to corroborate the findings from the sensitivity experiments by analyzing the differences between perturbed and unperturbed members in operational ensemble forecasts of ECMWF. Finally, we give an explanation of how symmetric, zero‐mean perturbations can lead to a unidirectional response when applied in a nonlinear system.

Dark energy model with very large-scale inhomogeneity

We consider a dynamical model for dark energy based on an ultralight mass scalar field with very large-scale inhomogeneities. This model may cause observable impacts on the anisotropic properties of the cosmic microwave background (CMB) intensity and luminosity distance. We formulate the model as the cosmological perturbations of the superhorizon scales, focusing on the local region of our universe. Moreover, we investigated the characteristic properties of the late-time evolution of inhomogeneous dark energy. Our numerical solutions show that the model can mimic the standard $\Lambda$CDM cosmology while including spatially dependent dark energy with flexible ranges of the model parameters. We put a constraint on the amplitude of these inhomogeneities of the dark energy on very large scales with the observations of the CMB anisotropies. We also discuss their influence on the estimation of the luminosity distance.

Nonperturbative running of the quark mass for Nf=3 QCD from the chirally rotated Schrödinger functional

We study the renormalization group (RG) running of the quark mass, for ${N}_{f}=3$ QCD with Wilson fermions in a mixed action setup, with standard Schr\"odinger functional (SF) boundary conditions for sea quarks and chirally rotated Schr\"odinger functional ($\ensuremath{\chi}\mathrm{SF}$) boundary conditions for valence quarks. This necessitates the tuning of the boundary factor ${z}_{\mathrm{f}}({g}_{0}^{2})$ of the $\ensuremath{\chi}\mathrm{SF}$ valence action, in order to ensure that QCD symmetries are fully recovered in the continuum. The properties of this novel setup are monitored through the ratios ${Z}_{S}/{Z}_{P}$ and ${\mathrm{\ensuremath{\Sigma}}}_{S}/{\mathrm{\ensuremath{\Sigma}}}_{P}$ of the renormalization parameters and step scaling functions of the scalar and pseudoscalar densities. Where comparison is possible, our ${Z}_{S}/{Z}_{P}$ results are found to agree with previous determinations, based on a mass ratio method [G. M. de Divitiis et al. (ALPHA Collaboration), Eur. Phys. J. C 79, 797 (2019)] and Ward identities [ J. Heitger et al. (ALPHA Collaboration), Eur. Phys. J. C 80, 765 (2020); J. Heitger et al., Eur. Phys. J. C 81, 606 (2021)], with Schr\"odinger functional boundary conditions. The behavior of ${\mathrm{\ensuremath{\Sigma}}}_{S}/{\mathrm{\ensuremath{\Sigma}}}_{P}$ confirms the theoretical expectations of $\ensuremath{\chi}\mathrm{SF}$ QCD, related to the restoration of the theory's symmetries in the continuum limit. From the step-scaling function of the pseudoscalar density we obtain the quark mass RG-running function from hadronic to perturbative energy scales. This is fully compatible with the earlier result obtained in a similar setup for Wilson quarks with Schr\"odinger functional boundary conditions [I. Campos et al., Eur. Phys. J. C 78, 387 (2018)], and provides a strong universality test for the two lattice setups.