scholarly journals Axisymmetric inertial modes in a spherical shell at low Ekman numbers

2018 ◽  
Vol 844 ◽  
pp. 597-634 ◽  
Author(s):  
M. Rieutord ◽  
L. Valdettaro
Keyword(s):  
Spherical Shell ◽  
Numerical Solutions ◽  
Asymptotic Properties ◽  
Analytical Formula ◽  
Shear Layers ◽  
Periodic Forcing ◽  
Oscillatory Flows ◽  
Ekman Number ◽  
Critical Latitude

We investigate the asymptotic properties of axisymmetric inertial modes propagating in a spherical shell when viscosity tends to zero. We identify three kinds of eigenmodes whose eigenvalues follow very different laws as the Ekman number $E$ becomes very small. First are modes associated with attractors of characteristics that are made of thin shear layers closely following the periodic orbit traced by the characteristic attractor. Second are modes made of shear layers that connect the critical latitude singularities of the two hemispheres of the inner boundary of the spherical shell. Third are quasi-regular modes associated with the frequency of neutral periodic orbits of characteristics. We thoroughly analyse a subset of attractor modes for which numerical solutions point to an asymptotic law governing the eigenvalues. We show that three length scales proportional to $E^{1/6}$, $E^{1/4}$ and $E^{1/3}$ control the shape of the shear layers that are associated with these modes. These scales point out the key role of the small parameter $E^{1/12}$ in these oscillatory flows. With a simplified model of the viscous Poincaré equation, we can give an approximate analytical formula that reproduces the velocity field in such shear layers. Finally, we also present an analysis of the quasi-regular modes whose frequencies are close to $\sin (\unicode[STIX]{x03C0}/4)$ and explain why a fluid inside a spherical shell cannot respond to any periodic forcing at this frequency when viscosity vanishes.

scholarly journals Internal shear layers from librating objects

2017 ◽  
Vol 826 ◽  
pp. 653-675 ◽  
Author(s):  
Stéphane Le Dizès ◽  
Michael Le Bars
Keyword(s):  
Asymptotic Solution ◽  
Shear Layer ◽  
Numerical Solutions ◽  
Shear Layers ◽  
Similarity Solutions ◽  
Curvature Radius ◽  
Ekman Number ◽  
Layer Structures ◽  
Two Parameters ◽  
Critical Latitude

In this work, we analyse the internal shear layer structures generated by the libration of an axisymmetric object in an unbounded fluid rotating at a rotation rate $\unicode[STIX]{x1D6FA}^{\ast }$ using direct numerical simulation and small Ekman number asymptotic analysis. We consider weak libration amplitude and libration frequency $\unicode[STIX]{x1D714}^{\ast }$ within the inertial wave interval $(0,2\unicode[STIX]{x1D6FA}^{\ast })$ such that the fluid dynamics is mainly described by a linear axisymmetric harmonic solution. The internal shear layer structures appear along the characteristic cones of angle $\unicode[STIX]{x1D703}_{c}=\text{acos}(\unicode[STIX]{x1D714}^{\ast }/(2\unicode[STIX]{x1D6FA}^{\ast }))$ which are tangent to the librating object at so-called critical latitudes. These layers correspond to thin viscous regions where the singularities of the inviscid solution are smoothed. We assume that the velocity field in these layers is described by the class of similarity solutions introduced by Moore & Saffman (Phil. Trans. R. Soc. Lond. A, vol. 264, 1969, pp. 597–634). These solutions are characterized by two parameters only: a real parameter $m$, which measures the strength of the underlying singularity, and a complex amplitude coefficient $C_{0}$. We first analyse the case of a disk for which a general asymptotic solution for small Ekman numbers is known when the disk is in a solid plane. We demonstrate that the numerical solutions obtained for a free disk and for a disk in a solid plane are both well described by the asymptotic solution and by its similarity form within the internal shear layers. For the disk, we obtain a parameter $m=1$ corresponding to a Dirac source at the edge of the disk and a coefficient $C_{0}\propto E^{1/6}$ where $E$ is the Ekman number. The case of a smoothed librating object such as a spheroid is found to be different. By asymptotically matching the boundary layer solution to similarity solutions close to a critical latitude on the surface, we show that the adequate parameter $m$ for the similarity solution is $m=5/4$, leading to a coefficient $C_{0}\propto E^{1/12}$, that is larger than for the case of a disk for small Ekman numbers. A simple general expression for $C_{0}$ valid for any axisymmetric object is obtained as a function of the local curvature radius at the critical latitude in agreement with this change of scaling. This result is tested and validated against direct numerical simulations.

scholarly journals Inertial waves in a rotating spherical shell: attractors and asymptotic spectrum

2001 ◽  
Vol 435 ◽  
pp. 103-144 ◽  
Author(s):  
M. RIEUTORD ◽  
B. GEORGEOT ◽  
L. VALDETTARO
Keyword(s):  
Lyapunov Exponent ◽  
Velocity Field ◽  
Spherical Shell ◽  
Asymptotic Properties ◽  
Analytical Form ◽  
Shear Layers ◽  
Three Dimensions ◽  
Geometrical Approach ◽  
Steady Flows ◽  
Whole Space

We investigate the asymptotic properties of inertial modes confined in a spherical shell when viscosity tends to zero. We first consider the mapping made by the characteristics of the hyperbolic equation (Poincaré's equation) satisfied by inviscid solutions. Characteristics are straight lines in a meridional section of the shell, and the mapping shows that, generically, these lines converge towards a periodic orbit which acts like an attractor (the associated Lyapunov exponent is always negative or zero). We show that these attractors exist in bands of frequencies the size of which decreases with the number of reflection points of the attractor. At the bounding frequencies the associated Lyapunov exponent is generically either zero or minus infinity. We further show that for a given frequency the number of coexisting attractors is finite.We then examine the relation between this characteristic path and eigensolutions of the inviscid problem and show that in a purely two-dimensional problem, convergence towards an attractor means that the associated velocity field is not square-integrable. We give arguments which generalize this result to three dimensions. Then, using a sphere immersed in a fluid filling the whole space, we study the critical latitude singularity and show that the velocity field diverges as 1/√d, d being the distance to the characteristic grazing the inner sphere.We then consider the viscous problem and show how viscosity transforms singularities into internal shear layers which in general reveal an attractor expected at the eigenfrequency of the mode. Investigating the structure of these shear layers, we find that they are nested layers, the thinnest and most internal layer scaling with E1/3, E being the Ekman number; for this latter layer, we give its analytical form and show its similarity to vertical 1/3-shear layers of steady flows. Using an inertial wave packet travelling around an attractor, we give a lower bound on the thickness of shear layers and show how eigenfrequencies can be computed in principle. Finally, we show that as viscosity decreases, eigenfrequencies tend towards a set of values which is not dense in [0, 2Ω], contrary to the case of the full sphere (Ω is the angular velocity of the system).Hence, our geometrical approach opens the possibility of describing the eigenmodes and eigenvalues for astrophysical/geophysical Ekman numbers (10−10–10−20), which are out of reach numerically, and this for a wide class of containers.

Inertial waves in a rotating spherical shell

1997 ◽  
Vol 341 ◽  
pp. 77-99 ◽  
Author(s):  
M. RIEUTORD ◽  
L. VALDETTARO
Keyword(s):  
Spherical Shell ◽  
Cross Section ◽  
Shear Layers ◽  
Inertial Waves ◽  
Ekman Number ◽  
Distribution Of Eigenvalues ◽  
Simple Rules ◽  
Meridional Cross Section ◽  
Zero Viscosity ◽  
The Web

The structure and spectrum of inertial waves of an incompressible viscous fluid inside a spherical shell are investigated numerically. These modes appear to be strongly featured by a web of rays which reflect on the boundaries. Kinetic energy and dissipation are indeed concentrated on thin conical sheets, the meridional cross-section of which forms the web of rays. The thickness of the rays is in general independent of the Ekman number E but a few cases show a scaling with E1/4 and statistical properties of eigenvalues indicate that high-wavenumber modes have rays of width O(E1/3). Such scalings are typical of Stewartson shear layers. It is also shown that the web of rays depends on the Ekman number and shows bifurcations as this number is decreased.This behaviour also implies that eigenvalues do not evolve smoothly with viscosity. We infer that only the statistical distribution of eigenvalues may follow some simple rules in the asymptotic limit of zero viscosity.

Analysis of singular inertial modes in a spherical shell: the slender toroidal shell model

2002 ◽  
Vol 463 ◽  
pp. 345-360 ◽  
Author(s):  
M. RIEUTORD ◽  
L. VALDETTARO ◽  
B. GEORGEOT
Keyword(s):  
Spherical Shell ◽  
Numerical Solutions ◽  
Three Dimensional ◽  
Two Dimensions ◽  
Toroidal Shell ◽  
Perturbative Approach ◽  
Quadratic Potential ◽  
Ekman Number ◽  
Damping Rate ◽  
Thin Spherical Shell

We derive the asymptotic spectrum (as the Ekman number E → 0) of axisymmetric inertial modes when the problem is restricted to two dimensions. We show that the damping rate of such modes scales with the square root of the Ekman number and that the width of the shear layers of the eigenfunctions scales with E1/4. The eigenfunctions obey a Schrödinger equation with a quadratic potential; we provide the analytical expression for eigenvalues (frequency and damping rate). These results validate the picture that attractors act like a potential well, trapping inertial waves which resist confinement owing to viscosity. Using three-dimensional numerical solutions, we show that the results can be applied to equatorially trapped modes in a thin spherical shell; in fact, these two-dimensional solutions give the first step (the zeroth order) of a perturbative approach to three-dimensional solutions in a spherical shell. Our method is applicable in a straightforward way to any other container where bi-dimensionality dominates.

scholarly journals SIR Model for Dengue Disease with Effect of Dengue Vaccination

2018 ◽  
Vol 2018 ◽  
pp. 1-14 ◽  
Author(s):  
Pratchaya Chanprasopchai ◽  
I. Ming Tang ◽  
Puntani Pongsumpun
Keyword(s):  
Mathematical Model ◽  
Mortality Rate ◽  
Dengue Virus ◽  
Numerical Solutions ◽  
Specific Treatment ◽  
Sir Model ◽  
Dengue Vaccine ◽  
Dengue Disease ◽  
Modelling Methods

The dengue disease is caused by dengue virus, and there is no specific treatment. The medical care by experienced physicians and nurses will save life and will lower the mortality rate. A dengue vaccine to control the disease is available in Thailand since late 2016. A mathematical model would be an important way to analyze the effects of the vaccination on the transmission of the disease. We have formulated an SIR (susceptible-infected-recovered) model of the transmission of the disease which includes the effect of vaccination and used standard dynamical modelling methods to analyze the effects. The equilibrium states and their stabilities are investigated. The trajectories of the numerical solutions plotted into the 2D planes and 3D spaces are presented. The main contribution is determining the role of dengue vaccination in the model. From the analysis, we find that there is a significant reduction in the total hospitalization time needed to treat the illness.

scholarly journals Extended logistic growth model for heterogeneous populations

2017 ◽  
Author(s):  
Wang Jin ◽  
Scott W McCue ◽  
Matthew J Simpson
Keyword(s):  
Cell Proliferation ◽  
Growth Model ◽  
Numerical Solutions ◽  
Cellular Level ◽  
Logistic Growth ◽  
Logistic Growth Model ◽  
Living Tissues ◽  
Discrete Mathematical Model ◽  
The Continuum

AbstractCell proliferation is the most important cellular-level mechanism responsible for regulating cell population dynamics in living tissues. Modern experimental procedures show that the proliferation rates of individual cells can vary significantly within the same cell line. However, in the mathematical biology literature, cell proliferation is typically modelled using a classical logistic equation which neglects variations in the proliferation rate. In this work, we consider a discrete mathematical model of cell migration and cell proliferation, modulated by volume exclusion (crowding) effects, with variable rates of proliferation across the total population. We refer to this variability as heterogeneity. Constructing the continuum limit of the discrete model leads to a generalisation of the classical logistic growth model. Comparing numerical solutions of the model to averaged data from discrete simulations shows that the new model captures the key features of the discrete process. Applying the extended logistic model to simulate a proliferation assay using rates from recent experimental literature shows that neglecting the role of heterogeneity can, at times, lead to misleading results.

scholarly journals The Potential Role of the ZIKV NS5 Nuclear Spherical-Shell Structures in Cell Type-Specific Host Immune Modulation during ZIKV Infection

Cells ◽  
2019 ◽  
Vol 8 (12) ◽  
pp. 1519 ◽  
Author(s):  
Min Jie Alvin Tan ◽  
Kitti Wing Ki Chan ◽  
Ivan H. W. Ng ◽  
Sean Yao Zu Kong ◽  
Chin Piaw Gwee ◽  
...  
Keyword(s):  
Spherical Shell ◽  
Immune Modulation ◽  
Transcriptional Response ◽  
Type I ◽  
Primary Role ◽  
Structural Protein ◽  
Cell Type ◽  
Zikv Infection ◽  
Cell Type Specific

The Zika virus (ZIKV) non-structural protein 5 (NS5) plays multiple viral and cellular roles during infection, with its primary role in virus RNA replication taking place in the cytoplasm. However, immunofluorescence assay studies have detected the presence of ZIKV NS5 in unique spherical shell-like structures in the nuclei of infected cells, suggesting potentially important cellular roles of ZIKV NS5 in the nucleus. Hence ZIKV NS5′s subcellular distribution and localization must be tightly regulated during ZIKV infection. Both ZIKV NS5 expression or ZIKV infection antagonizes type I interferon signaling, and induces a pro-inflammatory transcriptional response in a cell type-specific manner, but the mechanisms involved and the role of nuclear ZIKV NS5 in these cellular functions has not been elucidated. Intriguingly, these cells originate from the brain and placenta, which are also organs that exhibit a pro-inflammatory signature and are known sites of pathogenesis during ZIKV infection in animal models and humans. Here, we discuss the regulation of the subcellular localization of the ZIKV NS5 protein, and its putative role in the induction of an inflammatory response and the occurrence of pathology in specific organs during ZIKV infection.

scholarly journals Generalized Independence in the q-Voter Model: How Do Parameters Influence the Phase Transition?

Entropy ◽  
2020 ◽  
Vol 22 (1) ◽  
pp. 120 ◽  
Author(s):  
Angelika Abramiuk ◽  
Katarzyna Sznajd-Weron
Keyword(s):  
Phase Transition ◽  
Critical Point ◽  
Broad Spectrum ◽  
Order Parameter ◽  
Analytical Formula ◽  
Voter Model ◽  
Pair Approximation ◽  
Scaling Relation ◽  
Independent Behavior

We study the q-voter model with flexibility, which allows for describing a broad spectrum of independence from zealots, inflexibility, or stubbornness through noisy voters to self-anticonformity. Analyzing the model within the pair approximation allows us to derive the analytical formula for the critical point, below which an ordered (agreement) phase is stable. We determine the role of flexibility, which can be understood as an amount of variability associated with an independent behavior, as well as the role of the average network degree in shaping the character of the phase transition. We check the existence of the scaling relation, which previously was derived for the Sznajd model. We show that the scaling is universal, in a sense that it does not depend neither on the size of the group of influence nor on the average network degree. Analyzing the model in terms of the rescaled parameter, we determine the critical point, the jump of the order parameter, as well as the width of the hysteresis as a function of the average network degree ⟨ k ⟩ and the size of the group of influence q.

scholarly journals Dynamic Crushing Analysis of a Three-Dimensional Re-Entrant Auxetic Cellular Structure

Materials ◽  
2019 ◽  
Vol 12 (3) ◽  
pp. 460 ◽  
Author(s):  
Tao Wang ◽  
Zhen Li ◽  
Liangmo Wang ◽  
Zhengdong Ma ◽  
Gregory Hulbert
Keyword(s):  
Energy Dissipation ◽  
Relative Density ◽  
Cellular Structure ◽  
Three Dimensional ◽  
Analytical Formula ◽  
Absorption Capacity ◽  
Plastic Energy ◽  
Theoretical Predictions ◽  
Auxetic Structure

Dynamic behaviors of the three-dimensional re-entrant auxetic cellular structure have been investigated by performing beam-based crushing simulation. Detailed deformation process subjected to various crushing velocities has been described, where three specific crushing modes have been identified with respect to the crushing velocity and the relative density. The crushing strength of the 3D re-entrant auxetic structure reveals to increase with increasing crushing velocity and relative density. Moreover, an analytical formula of dynamic plateau stress has been deduced, which has been validated to present theoretical predictions agreeing well with simulation results. By establishing an analytical model, the role of relative density on the energy absorption capacity of the 3D re-entrant auxetic structure has been further studied. The results indicate that the specific plastic energy dissipation is increased by increasing the relative density, while the normalized plastic energy dissipation has an opposite sensitivity to the relative density when the crushing velocity exceeds the critical transition velocity. The study in this work can provide insights into the dynamic property of the 3D re-entrant auxetic structure and provides an extensive reference for the crushing resistance design of the auxetic structure.

Sign in / Sign up

Export Citation Format

Share Document