scholarly journals A stochastic controlled Schroedinger equation: convergence and robust stability for numerical solutions

2021 ◽  
pp. 178-184
Author(s):  
Cutberto Romero-Melendez ◽  
David Castillo-Fernandez

In this paper we study the stochastic stability of numerical solutions of a stochastic controlled Schr¨odinger equation. We investigate the boundedness in second moment, the convergence and the stability of the zero solution for this equation, using two new definitions of almost sure exponential robust stability and asymptotic stability, for the Euler-Maruyama numerical scheme. Considering that the diffusion term is controlled, by using the method of Lyapunov functions and the corresponding diffusion operator associated, we apply techniques of X. Mao and A. Tsoi for achieve our task. Finally, we illustrate this method with a problem in Nuclear Magnetic Resonance (NMR).

2014 ◽  
Vol 24 (3) ◽  
pp. 635-646 ◽  
Author(s):  
Deqiong Ding ◽  
Qiang Ma ◽  
Xiaohua Ding

Abstract In this paper, a NonStandard Finite Difference (NSFD) scheme is constructed, which can be used to determine numerical solutions for an epidemic model with vaccination. Here the NSFD method is employed to derive a set of difference equations for the epidemic model with vaccination. We show that difference equations have the same dynamics as the original differential system, such as the positivity of the solutions and the stability of the equilibria, without being restricted by the time step. Our proof of global stability utilizes the method of Lyapunov functions. Numerical simulation illustrates the effectiveness of our results


2001 ◽  
Vol 43 (2) ◽  
pp. 269-278 ◽  
Author(s):  
D. D. Bainov ◽  
I. M. Stamova

AbstractWe consider the stability of the zero solution of a system of impulsive functional-differential equations. By means of piecewise continuous functions, which are generalizations of classical Lyapunov functions, and using a technique due to Razumikhin, sufficient conditions are found for stability, uniform stability and asymptotical stability of the zero solution of these equations. Applications to impulsive population dynamics are also discussed.


Symmetry ◽  
2021 ◽  
Vol 13 (2) ◽  
pp. 353
Author(s):  
Ligia Munteanu ◽  
Dan Dumitriu ◽  
Cornel Brisan ◽  
Mircea Bara ◽  
Veturia Chiroiu ◽  
...  

The purpose of this paper is to study the sliding mode control as a Ricci flow process in the context of a three-story building structure subjected to seismic waves. The stability conditions result from two Lyapunov functions, the first associated with slipping in a finite period of time and the second with convergence of trajectories to the desired state. Simulation results show that the Ricci flow control leads to minimization of the displacements of the floors.


2002 ◽  
Vol 452 ◽  
pp. 163-187 ◽  
Author(s):  
C. L. BURCHAM ◽  
D. A. SAVILLE

A liquid bridge is a column of liquid, pinned at each end. Here we analyse the stability of a bridge pinned between planar electrodes held at different potentials and surrounded by a non-conducting, dielectric gas. In the absence of electric fields, surface tension destabilizes bridges with aspect ratios (length/diameter) greater than π. Here we describe how electrical forces counteract surface tension, using a linearized model. When the liquid is treated as an Ohmic conductor, the specific conductivity level is irrelevant and only the dielectric properties of the bridge and the surrounding gas are involved. Fourier series and a biharmonic, biorthogonal set of Papkovich–Fadle functions are used to formulate an eigenvalue problem. Numerical solutions disclose that the most unstable axisymmetric deformation is antisymmetric with respect to the bridge’s midplane. It is shown that whilst a bridge whose length exceeds its circumference may be unstable, a sufficiently strong axial field provides stability if the dielectric constant of the bridge exceeds that of the surrounding fluid. Conversely, a field destabilizes a bridge whose dielectric constant is lower than that of its surroundings, even when its aspect ratio is less than π. Bridge behaviour is sensitive to the presence of conduction along the surface and much higher fields are required for stability when surface transport is present. The theoretical results are compared with experimental work (Burcham & Saville 2000) that demonstrated how a field stabilizes an otherwise unstable configuration. According to the experiments, the bridge undergoes two asymmetric transitions (cylinder-to-amphora and pinch-off) as the field is reduced. Agreement between theory and experiment for the field strength at the pinch-off transition is excellent, but less so for the change from cylinder to amphora. Using surface conductivity as an adjustable parameter brings theory and experiment into agreement.


2011 ◽  
Vol 2011 ◽  
pp. 1-12 ◽  
Author(s):  
H. Saberi Najafi ◽  
A. Refahi Sheikhani ◽  
A. Ansari

We analyze the stability of three classes of distributed order fractional differential equations (DOFDEs) with respect to the nonnegative density function. In this sense, we discover a robust stability condition for these systems based on characteristic function and new inertia concept of a matrix with respect to the density function. Moreover, we check the stability of a distributed order fractional WINDMI system to illustrate the validity of proposed procedure.


2014 ◽  
Vol 2014 ◽  
pp. 1-6 ◽  
Author(s):  
Guiying Chen ◽  
Linshan Wang

The stability of a class of static interval neural networks with time delay in the leakage term is investigated. By using the method ofM-matrix and the technique of delay differential inequality, we obtain some sufficient conditions ensuring the global exponential robust stability of the networks. The results in this paper extend the corresponding conclusions without leakage delay. An example is given to illustrate the effectiveness of the obtained results.


2010 ◽  
Vol 2010 ◽  
pp. 1-23 ◽  
Author(s):  
Josef Diblík ◽  
Denys Ya. Khusainov ◽  
Irina V. Grytsay ◽  
Zdenĕk Šmarda

Many processes are mathematically simulated by systems of discrete equations with quadratic right-hand sides. Their stability is thought of as a very important characterization of the process. In this paper, the method of Lyapunov functions is used to derive classes of stable quadratic discrete autonomous systems in a critical case in the presence of a simple eigenvalueλ=1of the matrix of linear terms. In addition to the stability investigation, we also estimate stability domains.


Author(s):  
Marta J. Reith ◽  
Daniel Bachrathy ◽  
Gabor Stepan

Multi-cutter turning systems bear huge potential in increasing cutting performance. In this study we show that the stable parameter region can be extended by the optimal tuning of system parameters. The optimal parameter regions can be identified by means of stability charts. Since the stability boundaries are highly sensitive to the dynamical parameters of the machine tool, the reliable exploitation of the so-called stability pockets is limited. Still, the lower envelope of the stability lobes is an appropriate upper boundary function for optimization purposes with an objective function taken for maximal material removal rates. This lower envelope is computed by the Robust Stability Computation method presented in the paper. It is shown in this study, that according to theoretical results obtained for optimally tuned cutters, the safe stable machining parameter region can significantly be extended, which has also been validated by machining tests.


Author(s):  
Rama K. Yedavalli

This paper revisits the problem of checking the robust stability of matrix families generated by ‘Unstructured Convex Combinations’ of user supplied or externally supplied Vertex Matrices. A previous solution given by the author for this problem involved complete dependence on the quantitative (eigenvalue information) of a set of special matrices labeled the Kronecker Nonsingularity (KN) matrices. In this solution, the ‘convexity’ property is not explicit and transparent, to the extent that, unfortunately, the accuracy of the solution itself is being questioned and not embraced by the peer community. To erase this unforunate and unwarranted image of this author (in this specific problem) in the minds of the peer community, in this paper, the author treads a new path to find a solution that brings out the convexity property in an explicit and understandable way. In the new solution presented in this paper, we combine the qualitative (sign) as well as quantitative (magnitude) information of these KN matrices and present a vertex solution in which the convexity property of the solution is transparent making it more elegant and accepatble to the peer community, than the previous solution. The new solution clearly underscores the importance of using the sign structure of a matrix in assessing the stability of a matrix. This new solution is made possible by the new insight provided by the qualitative (sign) stability/instability derived from ecological principles. Examples are given which clearly demonstrate effectiveness of the new, convexity based algorithm. It is hoped that this new solution will be embraced by the peer community.


2018 ◽  
Vol 615 ◽  
pp. A67 ◽  
Author(s):  
P. A. González-Morales ◽  
E. Khomenko ◽  
T. P. Downes ◽  
A. de Vicente

The interaction of plasma with magnetic field in the partially ionised solar atmosphere is frequently modelled via a single-fluid approximation, which is valid for the case of a strongly coupled collisional media, such as solar photosphere and low chromosphere. Under the single-fluid formalism the main non-ideal effects are described by a series of extra terms in the generalised induction equation and in the energy conservation equation. These effects are: Ohmic diffusion, ambipolar diffusion, the Hall effect, and the Biermann battery effect. From the point of view of the numerical solution of the single-fluid equations, when ambipolar diffusion or Hall effects dominate can introduce severe restrictions on the integration time step and can compromise the stability of the numerical scheme. In this paper we introduce two numerical schemes to overcome those limitations. The first of them is known as super time-stepping (STS) and it is designed to overcome the limitations imposed when the ambipolar diffusion term is dominant. The second scheme is called the Hall diffusion scheme (HDS) and it is used when the Hall term becomes dominant. These two numerical techniques can be used together by applying Strang operator splitting. This paper describes the implementation of the STS and HDS schemes in the single-fluid code MANCHA3D. The validation for each of these schemes is provided by comparing the analytical solution with the numerical one for a suite of numerical tests.


Sign in / Sign up

Export Citation Format

Share Document