scholarly journals Analytic Exploration of Safe Basins in a Benchmark Problem of Forced Escape

Author(s):  
Gleb Karmi ◽  
Pavel Kravetc ◽  
Oleg Gendelman

Abstract The paper presents an analytic approach for predicting the safe basins (SB) in a plane of initial conditions (IC) for escape of classical particle from the potential well under harmonic forcing. The solution is based on the approximation of isolated resonance, which reduces the dynamics to conservative flow on a two-dimensional resonance manifold (RM). Such a reduction allows easy distinction between escaping and non-escaping ICs. As a benchmark potential, we choose a common parabolic-quartic well with truncation at varying energy levels. The method allows accurate predictions of the SB boundaries for relatively low forcing amplitudes. The derived SBs demonstrate an unexpected set of properties, including decomposition into two disjoint zones in the IC plane for a certain range of parameters. The latter peculiarity stems from two qualitatively different escape mechanisms on the RM. For higher forcing values, the accuracy of the analytic predictions decreases to some extent due to the inaccuracies of the basic isolated resonance approximation, but mainly due to the erosion of the SB boundaries caused by the secondary resonances. Nevertheless, even in this case the analytic approximation can serve as a viable initial guess for subsequent numeric estimation of the SB boundaries.

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
E. V. Kirichenko ◽  
V. A. Stephanovich

AbstractWe study the joint effect of disorder and Coulomb interaction screening on the exciton spectra in two-dimensional (2D) structures. These can be van der Waals structures or heterostructures of organic (polymeric) semiconductors as well as inorganic substances like transition metal dichalcogenides. We consider 2D screened hydrogenic problem with Rytova–Keldysh interaction by means of so-called fractional Scrödinger equation. Our main finding is that above synergy between screening and disorder either destroys the exciton (strong screening) or promote the creation of a bound state, leading to its collapse in the extreme case. Our second finding is energy levels crossing, i.e. the degeneracy (with respect to index $$\mu $$ μ ) of the exciton eigenenergies at certain discrete value of screening radius. Latter effects may also be related to the quantum manifestations of chaotic exciton behavior in above 2D semiconductor structures. Hence, they should be considered in device applications, where the interplay between dielectric screening and disorder is important.


2013 ◽  
Vol 28 (16) ◽  
pp. 1350064 ◽  
Author(s):  
CATARINA BASTOS ◽  
ORFEU BERTOLAMI ◽  
NUNO COSTA DIAS ◽  
JOÃO NUNO PRATA

We consider a noncommutative description of graphene. This description consists of a Dirac equation for massless Dirac fermions plus noncommutative corrections, which are treated in the presence of an external magnetic field. We argue that, being a two-dimensional Dirac system, graphene is particularly interesting to test noncommutativity. We find that momentum noncommutativity affects the energy levels of graphene and we obtain a bound for the momentum noncommutative parameter.


1975 ◽  
Vol 69 (3) ◽  
pp. 417-443 ◽  
Author(s):  
Peter B. Rhines

Two-dimensional eddies in a homogeneous fluid at large Reynolds number, if closely packed, are known to evolve towards larger scales. In the presence of a restoring force, the geophysical beta-effect, this cascade produces a field of waves without loss of energy, and the turbulent migration of the dominant scale nearly ceases at a wavenumber kβ = (β/2U)½ independent of the initial conditions other than U, the r.m.s. particle speed, and β, the northward gradient of the Coriolis frequency.The conversion of turbulence into waves yields, in addition, more narrowly peaked wavenumber spectra and less fine-structure in the spatial maps, while smoothly distributing the energy about physical space.The theory is discussed, using known integral constraints and similarity solutions, model equations, weak-interaction wave theory (which provides the terminus for the cascade) and other linearized instability theory. Computer experiments with both finite-difference and spectral codes are reported. The central quantity is the cascade rate, defined as \[ T = 2\int_0^{\infty} kF(k)dk/U^3\langle k\rangle , \] where F is the nonlinear transfer spectrum and 〈k〉 the mean wavenumber of the energy spectrum. (In unforced inviscid flow T is simply U−1d〈k〉−1/dt, or the rate at which the dominant scale expands in time t.) T is shown to have a mean value of 3·0 × 10−2 for pure two-dimensional turbulence, but this decreases by a factor of five at the transition to wave motion. We infer from weak-interaction theory even smaller values for k [Lt ] kβ.After passing through a state of propagating waves, the homogeneous cascade tends towards a flow of alternating zonal jets which, we suggest, are almost perfectly steady. When the energy is intermittent in space, however, model equations show that the cascade is halted simply by the spreading of energy about space, and then the end state of a zonal flow is probably not achieved.The geophysical application is that the cascade of pure turbulence to large scales is defeated by wave propagation, helping to explain why the energy-containing eddies in the ocean and atmosphere, though significantly nonlinear, fail to reach the size of their respective domains, and are much smaller. For typical ocean flows, $k_{\beta}^{-1} = 70\,{\rm km} $, while for the atmosphere, $k_{\beta}^{-1} = 1000\,{\rm km}$. In addition the cascade generates, by itself, zonal flow (or more generally, flow along geostrophic contours).


2013 ◽  
Vol 2013 ◽  
pp. 1-13
Author(s):  
Motoaki Saruwatari ◽  
Masa-aki Hashimoto ◽  
Ryohei Fukuda ◽  
Shin-ichiro Fujimoto

We investigate the r-process nucleosynthesis during the magnetohydrodynamical (MHD) explosion of a supernova in a helium star of 3.3 M⊙, where effects of neutrinos are taken into account using the leakage scheme in the two-dimensional (2D) hydrodynamic code. Jet-like explosion due to the combined effects of differential rotation and magnetic field is able to erode the lower electron fraction matter from the inner layers. We find that the ejected material of low electron fraction responsible for the r-process comes out from just outside the neutrino sphere deep inside the Fe-core. It is found that heavy element nucleosynthesis depends on the initial conditions of rotational and magnetic fields. In particular, the third peak of the distribution is significantly overproduced relative to the solar system abundances, which would indicate a possible r-process site owing to MHD jets in supernovae.


2009 ◽  
Vol 2009 ◽  
pp. 1-22 ◽  
Author(s):  
Edson D. Leonel

A phase transition from integrability to nonintegrability in two-dimensional Hamiltonian mappings is described and characterized in terms of scaling arguments. The mappings considered produce a mixed structure in the phase space in the sense that, depending on the combination of the control parameters and initial conditions, KAM islands which are surrounded by chaotic seas that are limited by invariant tori are observed. Some dynamical properties for the largest component of the chaotic sea are obtained and described in terms of the control parameters. The average value and the deviation of the average value for chaotic components of a dynamical variable are described in terms of scaling laws, therefore critical exponents characterizing a scaling function that describes a phase transition are obtained and then classes of universality are characterized. The three models considered are: The Fermi-Ulam accelerator model, a periodically corrugate waveguide, and variant of the standard nontwist map.


Author(s):  
D.Y. Ivanov ◽  

Here we consider the initial-boundary value problems in a homogeneous cylindrical domain YI Ω ×+ ( Ω+ is an open two-dimensional bounded simply connected domain with a boundary 5 ∂Ω ∈C , 2 \ Ω≡ Ω − + R is the open exterior of the domain Ω+ , [0, ] YI ≡ Y is the height of the cylinder) on a time interval [0, ] TI ≡ T . The initial conditions and the boundary conditions on the bases of the cylinder are zero, and the boundary conditions on the lateral surface of the cylinder are given by the function 1 2 wx x yt ( , , ,) ( 1 2 (, ) x x ∈∂Ω , Y y ∈ I , T t I ∈ ). An approximate solution of such problems is obtained through the combined use of the Fourier method and the collocation boundary element method based on piecewise quadratic interpolation (PQI). The solution to the problem in the cylinder is expanded in a Fourier series in terms of eigenfunctions of the operator 2 By yy ≡ ∂ with the corresponding zero boundary conditions. The coefficients of such a Fourier series are solutions of problems for two-dimensional heat equations 2 2 t ∇ =∂ + u u ku . With a low smoothness of the functions w in the variable y, the weight of solutions at large values of k increases and the accuracy of solving the problem in the cylinder decreases. To maintain accuracy on a uniform grid, the step of discretization of the boundary function w with respect to the variable y is decreased by a factor of j. Here j is an averaged value of the quantity Y k π depending on the function w. In addition, the steps of discretization of functions ( ) 2 exp − τ k with respect to the variable τ in domains τ≤πT k are reduced by a factor of 2 2 k π . The steps in the remaining ranges of values τ and the steps by the other variables remain unchanged. The approximate solutions obtained on the basis of this procedure converge stably to exact solutions in the 2 ( ) LI I Y T × -norm with a cubic velocity uniformly with respect to sets of functions w, bounded by norm of functions with low smoothness in the variable y, uniformly along the length of the generatrix of the cylinder Y , and uniformly in the domain Ω . The latter is also associated with the use of PQI along the curve ∂Ω over the variable 2 2 ρ≡ − r d , which is carried out at small values of r ( d and r are the distances from the observed point of the domain Ω to the boundary ∂Ω and to the current point of integration along ∂Ω , respectively). The theoretical conclusions are confirmed by the results of the numerical solution of the problem in a circular cylinder, where the dependence of the boundary functions w on y is given by the normalized eigenfunctions of the differential operator By which vary in a sufficiently large range of values of k .


2021 ◽  
Vol 31 (15) ◽  
Author(s):  
Penghe Ge ◽  
Hongjun Cao

The existence of chaos in the Rulkov neuron model is proved based on Marotto’s theorem. Firstly, the stability conditions of the model are briefly renewed through analyzing the eigenvalues of the model, which are very important preconditions for the existence of a snap-back repeller. Secondly, the Rulkov neuron model is decomposed to a one-dimensional fast subsystem and a one-dimensional slow subsystem by the fast–slow dynamics technique, in which the fast subsystem has sensitive dependence on the initial conditions and its snap-back repeller and chaos can be verified by numerical methods, such as waveforms, Lyapunov exponents, and bifurcation diagrams. Thirdly, for the two-dimensional Rulkov neuron model, it is proved that there exists a snap-back repeller under two iterations by illustrating the existence of an intersection of three surfaces, which pave a new way to identify the existence of a snap-back repeller.


2021 ◽  
Vol 932 ◽  
Author(s):  
Guangzhao Zhou ◽  
Andrea Prosperetti

It is known that the dripping of a liquid film on the underside of a plate can be suppressed by tilting the plate so as to cause a sufficiently strong flow. This paper uses two-dimensional numerical simulations in a closed-flow framework to study several aspects of this phenomenon. It is shown that, in quasi-equilibrium conditions, the onset of dripping is closely associated with the curvature of the wave crests approaching a well-defined maximum value. When dynamic effects become significant, this connection between curvature and dripping weakens, although the critical curvature remains a useful reference point as it is intimately related to the short length scales promoted by the Rayleigh–Taylor instability. In the absence of flow, when the film is on the underside of a horizontal plate, the concept of a limit curvature is relevant only for small liquid volumes close to a critical value. Otherwise, the drops that form have a smaller curvature and a large volume. The paper also illustrates the peculiarly strong dependence of the dripping transition on the initial conditions of the simulations. This feature prevents the development of phase maps dependent only on the governing parameters (Reynolds number, Bond number, etc.) similar to those available for film flow on the upper side of an inclined plate.


2021 ◽  
Vol 928 ◽  
Author(s):  
Yu Liang ◽  
Lili Liu ◽  
Zhigang Zhai ◽  
Juchun Ding ◽  
Ting Si ◽  
...  

Shock-tube experiments on eight kinds of two-dimensional multi-mode air–SF $_6$ interface with controllable initial conditions are performed to examine the dependence of perturbation growth on initial spectra. We deduce and demonstrate experimentally that the amplitude development of each mode is influenced by the mode-competition effect from quasi-linear stages. It is confirmed that the mode-competition effect is closely related to initial spectra, including the wavenumber, the phase and the initial amplitude of constituent modes. By considering both the mode-competition effect and the high-order harmonics effect, a nonlinear model is established based on initial spectra to predict the amplitude growth of each individual mode. The nonlinear model is validated by the present experiments and data in the literature by considering diverse initial spectra, shock intensities and density ratios. Moreover, the nonlinear model is successfully extended based on the superposition principle to predict the growths of the total perturbation width and the bubble/spike width from quasi-linear to nonlinear stages.


Sign in / Sign up

Export Citation Format

Share Document