Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical Bifurcation Analysis of Marine Ice Sheet Models
63
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The climate variability associated with the Pleistocene Ice Ages is one of the most fascinating puzzles in the Earth Sciences still awaiting a satisfactory explanation. In particular, the explanation of the dominant 100 kyr period of the glacial cycles over the last million years is a long-standing problem. Based on bifurcation analyses of low-order models, many theories have been suggested to explain these cycles and their frequency. The new aspect in this contribution is that, for the first time, numerical bifurcation analysis is applied to a two-dimensional marine ice sheet model with a dynamic grounding line. In this model, we find Hopf bifurcations with an oscillation period of about 100 kyr which may be relevant to glacial cycles.
rate research
Read More
Using various techniques from dynamical systems theory, we rigorously study an experimentally validated model by [Barkley et al., Nature, 526:550-553, 2015], which describes the rise of turbulent pipe flow via a PDE system of reduced complexity. The fast evolution of turbulence is governed by reaction-diffusion dynamics coupled to the centerline velocity, which evolves with advection of Burgers type and a slow relaminarization term. Applying to this model a spatial dynamics ansatz and geometric singular perturbation theory, we prove the existence of a heteroclinic loop between a turbulent and a laminar steady state and establish a cascade of bifurcations of various traveling waves mediating the transition to turbulence. The most complicated behaviour can be found in an intermediate Reynolds number regime, where the traveling waves exhibit arbitrarily long periodic-like dynamics indicating the onset of chaos. Our analysis provides a systematic mathematical approach to identifying the transition to spatio-temporal turbulent structures that may also be applicable to other models arising in fluid dynamics.
We perform a bifurcation analysis of the steady state solutions of Rayleigh--Benard convection with no-slip boundary conditions in two dimensions using a numerical method called deflated continuation. By combining this method with an initialisation strategy based on the eigenmodes of the conducting state, we are able to discover multiple solutions to this non-linear problem, including disconnected branches of the bifurcation diagram, without the need of any prior knowledge of the dynamics. One of the disconnected branches we find contains a s-shape bifurcation with hysteresis, which is the origin of the flow pattern that may be related to the dynamics of flow reversals in the turbulent regime. Linear stability analysis is also performed to analyse the steady and unsteady regimes of the solutions in the parameter space and to characterise the type of instabilities.
We study a self-similar solution of the kinetic equation describing weak wave turbulence in Bose-Einstein condensates. This solution presumably corresponds to an asymptotic behavior of a spectrum evolving from a broad class of initial data, and it features a non-equilibrium finite-time condensation of the wave spectrum $n(omega)$ at the zero frequency $omega$. The self-similar solution is of the second kind, and it satisfies boundary conditions corresponding to a nonzero constant spectrum (with all its derivative being zero) at $omega=0$ and a power-law asymptotic $n(omega) to omega^{-x}$ at $omega to infty ;; xin mathbb{R}^+$. Finding it amounts to solving a nonlinear eigenvalue problem, i.e. finding the value $x^*$ of the exponent $x$ for which these two boundary conditions can be satisfied simultaneously. To solve this problem we develop a new high-precision algorithm based on Chebyshev approximations and double exponential formulas for evaluating the collision integral, as well as the iterative techniques for solving the integro-differential equation for the self-similar shape function. This procedures allow to achieve a solution with accuracy $approx 4.7 %$ which is realized for $x^* approx 1.22$.
The dryland vegetation model proposed by Rietkerk and collaborators has been explored from a bifurcation perspective in several previous studies. Our aim here is to explore in some detail the bifurcation phenomena present when the coefficients of the model are allowed to vary in a wide range of parameters. In addition to the primary bifurcation parameter, the precipitation, we allow the two infiltration rate parameters to vary as well. We find that these two parameters control the size and stability of nonhomogeneous biomass states in a way that can be predicted. Further, they control when certain homogeneous and inhomogeneous (in space) periodic (in time) orbits exist. Finally, we show that the model possesses infinitely many unphysical steady state branches. We then present a modification of the model which eliminates these unphysical solutions, and briefly explore this new model for a fixed set of parameters.
An electrohydrodynamic (EHD) flow in a point-to-ring corona configuration is investigated experimentally, analytically and via a multiphysics numerical model. The interaction between the accelerated ions and the neutral gas molecules is modeled as an external body force in the Navier-Stokes equation (NSE). The gas flow characteristics are solved from conservation principles with spectral methods. The analytical and numerical simulation results are compared against experimental measurements of the cathode voltage, ion concentration, and velocity profiles. A nondimensional parameter, X, is formulated as the ratio of the local electric force to the inertial term in the NSE. In the region of X > 1, the electric force dominates the flow dynamics, while in the X << 1 region, the balance of viscous and inertial terms yields traditional pipe flow characteristics.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
