Source author record

Radostin D. Simitev

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

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

11 published item(s)

preprint2022arXiv

Estimation of Parameters for an Archetypal Model of Cardiomyocyte Membrane Potentials

Contemporary realistic mathematical models of single-cell cardiac electrical excitation are immensely detailed. Model complexity leads to parameter uncertainty, high computational cost and barriers to mechanistic understanding. There is a need for reduced models that are conceptually and mathematically simple but physiologically accurate. To this end, we consider an archetypal model of single-cell cardiac excitation that replicates the phase-space geometry of detailed cardiac models, but at the same time has a simple piecewise-linear form and a relatively low-dimensional configuration space. In order to make this archetypal model practically applicable, we develop and report a robust method for estimation of its parameter values from the morphology of single-stimulus action potentials derived from detailed ionic current models and from experimental myocyte measurements. The procedure is applied to five significant test cases and an excellent agreement with target biomarkers is achieved. Action potential duration restitution curves are also computed and compared to those of the target test models and data, demonstrating conservation of dynamical pacing behaviour by the fine-tuned archetypal model. An archetypal model that accurately reproduces a variety of wet-lab and synthetic electrophysiology data offers a number of specific advantages such as computational efficiency, as also demonstrated in the study. Open-source numerical code of the models and methods used is provided.

preprint2022arXiv

Flute and kink instabilities in a dynamically twisted flux tube with anisotropic plasma viscosity

Magnetic flux tubes such as those in the solar corona are subject to a number of instabilities. Important among them is the kink instability which plays a central part in the nanoflare theory of coronal heating, and for this reason in numerical simulations it is usually induced by tightly-controlled perturbations and studied in isolation. In contrast, we find that when disturbances are introduced in our magnetohydrodynamic flux tube simulations by dynamic twisting of the flow at the boundaries fluting modes of instability are readily excited. We also find that the flute instability, which has been theorised but rarely observed in the coronal context, is strongly enhanced when plasma viscosity is assumed anisotropic. We proceed to investigate the co-existence and competition between flute and kink instabilities for a range of values of the resistivity and of the parameters of the anisotropic and isotropic models of viscosity. We conclude that while the flute instability cannot prevent the kink from ultimately dominating, it can significantly delay its development especially at strong viscous anisotropy induced by intense magnetic fields.

preprint2020arXiv

Onset of inertial magnetoconvection in rotating fluid spheres

The onset of convection in the form of magneto-inertial waves in a rotating fluid sphere permeated by a constant axial electric current is studied through a perturbation analysis. Explicit expressions for the dependence of the Rayleigh number on the azimuthal wavenumber are derived in the limit of high thermal diffusivity. Results for the case of thermally infinitely conducting boundaries and for the case of nearly thermally insulating boundaries are obtained.

preprint2020arXiv

Regimes of thermo-compositional convection and related dynamos in rotating spherical shells

Convection and magnetic field generation in the Earth and planetary interiors are driven by both thermal and compositional gradients. In this work numerical simulations of finite-amplitude double-diffusive convection and dynamo action in rapidly rotating spherical shells full of incompressible two-component electrically-conducting fluid are reported. Four distinct regimes of rotating double-diffusive convection identified in a recent linear analysis (Silva et al., 2019, Geophys. Astrophys. Fluid Dyn., doi:10.1080/03091929.2019.1640875) are found to persist significantly beyond the onset of instability while their regime transitions remain abrupt. In the semi-convecting and the fingering regimes characteristic flow velocities are small compared to those in the thermally- and compositionally-dominated overturning regimes, while zonal flows remain weak in all regimes apart from the thermally-dominated one. Compositionally-dominated overturning convection exhibits significantly narrower azimuthal structures compared to all other regimes while differential rotation becomes the dominant flow component in the thermally-dominated case as driving is increased. Dynamo action occurs in all regimes apart from the regime of fingering convection. While dynamos persist in the semi-convective regime they are very much impaired by small flow intensities and very weak differential rotation in this regime which makes poloidal to toroidal field conversion problematic. The dynamos in the thermally-dominated regime include oscillating dipolar, quadrupolar and multipolar cases similar to the ones known from earlier parameter studies. Dynamos in the compositionally-dominated regime exhibit subdued temporal variation and remain predominantly dipolar due to weak zonal flow in this regime. These results significantly enhance our understanding of the primary drivers of planetary core flows and magnetic fields.

preprint2019arXiv

The onset of thermo-compositional convection in rotating spherical shells

Double-diffusive convection driven by both thermal and compositional buoyancy in a rotating spherical shell can exhibit a rather large number of behaviours often distinct from that of the single diffusive system. In order to understand how the differences in thermal and compositional molecular diffusivities determine the dynamics of thermo-compositional convection we investigate numerically the linear onset of convective instability in a double-diffusive setup. We construct an alternative equivalent formulation of the non-dimensional equations where the linearised double-diffusive problem is described by an effective Rayleigh number, $\text{Ra}$, measuring the amplitude of the combined buoyancy driving, and a second parameter, $α$, measuring the mixing of the thermal and compositional contributions. This formulation is useful in that it allows for the analysis of several limiting cases and reveals dynamical similarities in the parameters space which are not obvious otherwise. We analyse the structure of the critical curves in this $\text{Ra}-α$ space, explaining asymptotic behaviours in $α$, transitions between inertial and diffusive regimes, and transitions between large scale (fast drift) and small scale (slow drift) convection. We perform this analysis for a variety of diffusivities, rotation rates and shell aspect ratios showing where and when new modes of convection take place.

preprint2015arXiv

Dynamo Effects Near The Transition from Solar to Anti-Solar Differential Rotation

Numerical MHD simulations play increasingly important role for understanding mechanisms of stellar magnetism. We present simulations of convection and dynamos in density-stratified rotating spherical fluid shells. We employ a new 3D simulation code for the solution of a physically consistent anelastic model of the process with a minimum number of parameters. The reported dynamo simulations extend into a "buoyancy-dominated" regime where the buoyancy forcing is dominant while the Coriolis force is no longer balanced by pressure gradients and strong anti-solar differential rotation develops as a result. We find that the self-generated magnetic fields, despite being relatively weak, are able to reverse the direction of differential rotation from anti-solar to solar-like. We also find that convection flows in this regime are significantly stronger in the polar regions than in the equatorial region, leading to non-oscillatory dipole-dominated dynamo solutions, and to concentration of magnetic field in the polar regions. We observe that convection has different morphology in the inner and at the outer part of the convection zone simultaneously such that organized geostrophic convection columns are hidden below a near-surface layer of well-mixed highly-chaotic convection. While we focus the attention on the buoyancy-dominated regime, we also demonstrate that conical differential rotation profiles and persistent regular dynamo oscillations can be obtained in the parameter space of the rotation-dominated regime even within this minimal model.

preprint2015arXiv

Mathematical Modelling of Heart Rate Changes in the Mouse

The CVS is composed of numerous interacting and dynamically regulated physiological subsystems which each generate measurable periodic components such that the CVS can itself be presented as a system of weakly coupled oscillators. The interactions between these oscillators generate a chaotic blood pressure waveform signal, where periods of apparent rhythmicity are punctuated by asynchronous behaviour. It is this variability which seems to characterise the normal state. We used a standard experimental data set for the purposes of analysis and modelling. Arterial blood pressure waveform data was collected from conscious mice instrumented with radiotelemetry devices over $24$ hours, at a $100$ Hz and $1$ kHz time base. During a $24$ hour period, these mice display diurnal variation leading to changes in the cardiovascular waveform. We undertook preliminary analysis of our data using Fourier transforms and subsequently applied a series of both linear and nonlinear mathematical approaches in parallel. We provide a minimalistic linear and nonlinear coupled oscillator model and employed spectral and Hilbert analysis as well as a phase plane analysis. This provides a route to a three way synergistic investigation of the original blood pressure data by a combination of physiological experiments, data analysis viz. Fourier and Hilbert transforms and attractor reconstructions, and numerical solutions of linear and nonlinear coupled oscillator models. We believe that a minimal model of coupled oscillator models that quantitatively describes the complex physiological data could be developed via such a method. Further investigations of each of these techniques will be explored in separate publications.

preprint2015arXiv

Modelling afferent nerve responses to bladder filling

A sensation of fullness in the bladder is a regular experience, yet the mechanisms that act to generate this sensation remain poorly understood. This is an important issue because of the clinical problems that can result when this system does not function properly. The aim of the study group activity was to develop mathematical models that describe the mechanics of bladder filling, and how stretch modulates the firing rate of afferent nerves. Several models were developed, which were qualitatively consistent with experimental data obtained from a mouse model.

preprint2014arXiv

Quasi-geostrophic approximation of anelastic convection

The onset of convection in a rotating cylindrical annulus with parallel ends filled with a compressible fluid is studied in the anelastic approximation. Thermal Rossby waves propagating in the azimuthal direction are found as solutions. The analogy to the case of Boussinesq convection in the presence of conical end surfaces of the annular region is emphasized. As in the latter case the results can be applied as an approximation for the description of the onset of anelastic convection in rotating spherical fluid shells. Reasonable agreement with three-dimensional numerical results published by Jones et al. (J. Fluid Mech., vol. 634, 2009, pp. 291-319) for the latter problem is found. As in those results the location of the onset of convection shifts outward from the tangent cylinder with increasing number $N_ρ$ of density scale heights until it reaches the equatorial boundary. A new result is that at a much higher number $N_ρ$ the onset location returns to the interior of the fluid shell.

preprint2013arXiv

Proceedings of the 2013 UK National Conference on Geophysical, Astrophysical and Industrial Magnetohydrodynamics at Glasgow

This document presents an accurate record of the scientific programme of the 2013 UK National Conference on Geophysical, Astrophysical and Industrial Magnetohydrodynamics held on 23rd and 24th of May 2013 at the School of Mathematics and Statistics of the University of Glasgow. The abstracts of presented contributions are listed in chronological order and a full list of participants and their affiliations is included.

preprint2011arXiv

Solar cycle properties described by simple convection-driven dynamos

Simple models of magnetic field generation by convection in rotating spherical shells exhibit properties resembling those observed on the sun. The {assumption of the Boussinesq approximation made in these models} prevents a realistic description of the solar cycle, but through a physically motivated change in the boundary condition for the differential rotation the propagation of dynamo waves towards higher latitudes can be reversed at least at low latitudes.