Works of A.A. Samarskii on Computational Mathematics

2009 ◽  
Vol 9 (1) ◽  
pp. 5-36
Author(s):  
P.N. Vabishchevich

Abstract This is a review of the main results in computational mathematics that were obtained by the eminent Russian mathematician Alexander Andreevich Samarskii (February 19, 1919 – February 11, 2008). His outstanding research output addresses all the main questions that arise in the construction and justification of algorithms for the numerical solution of problems from mathematical physics. The remarkable works of A.A. Samarskii include statements of the main principles re- quired in the construction of difference schemes, rigorous mathematical proofs of the stability and convergence of these schemes, and also investigations of their algorith- mic implementation. A.A. Samarskii and his collaborators constructed and applied in practical calculations a large number of algorithms for solving various problems from mathematical physics, including thermal physics, gas dynamics, magnetic gas dynam- ics, plasma physics, ecology and other important models from the natural sciences.

1998 ◽  
Vol 08 (06) ◽  
pp. 1055-1070 ◽  
Author(s):  
ALEXANDER A. SAMARSKII ◽  
PETR P. MATUS ◽  
PETR N. VABISHCHEVICH

Nowadays the general theory of operator-difference schemes with operators acting in Hilbert spaces has been created for investigating the stability of the difference schemes that approximate linear problems of mathematical physics. In most cases a priori estimates which are uniform with respect to the t norms are usually considered. In the investigation of accuracy for evolutionary problems, special attention should be given to estimation of the difference solution in grid analogs of integral with respect to the time norms. In this paper a priori estimates in such norms have been obtained for two-level operator-difference schemes. Use of that estimates is illustrated by convergence investigation for schemes with weights for parabolic equation with the solution belonging to [Formula: see text].


2020 ◽  
Vol 72 (4) ◽  
pp. 73-77
Author(s):  
L.M. Tukenova ◽  

Mathematical models of oceanology are equations of the Navier-Stokes type, the construction of stable effective algorithms for their solution is associated with certain difficulties due to the well-known problems of setting boundary conditions, the presence of integro-differential relations, etc. In practice, when solving problems of oceanology, finitedifference methods are widely used, but there are no works in the literature devoted to theoretical studies of the stability and convergence of the algorithms used. In most cases, stability and convergence tests are established through computational experiments. Therefore, we believe that the development and mathematical substantiation of converging methods for solving the system of oceanology equations are urgent problems of computational mathematics. The paper studies variants of the fictitious domain method for a nonlinear ocean model. An existence theorem for the convergence of solutions to approximate models obtained using the fictitious domain method is investigated. An unimprovable estimate of the rate of convergence of the solution of the fictitious domain method is derived.


2021 ◽  
Vol 3 (2) ◽  
Author(s):  
Abdullah Al-Mamun ◽  
S. M. Arifuzzaman ◽  
Sk. Reza-E-Rabbi ◽  
Umme Sara Alam ◽  
Saiful Islam ◽  
...  

AbstractThe perspective of this paper is to characterize a Casson type of Non-Newtonian fluid flow through heat as well as mass conduction towards a stretching surface with thermophoresis and radiation absorption impacts in association with periodic hydromagnetic effect. Here heat absorption is also integrated with the heat absorbing parameter. A time dependent fundamental set of equations, i.e. momentum, energy and concentration have been established to discuss the fluid flow system. Explicit finite difference technique is occupied here by executing a procedure in Compaq Visual Fortran 6.6a to elucidate the mathematical model of liquid flow. The stability and convergence inspection has been accomplished. It has observed that the present work converged at, Pr ≥ 0.447 indicates the value of Prandtl number and Le ≥ 0.163 indicates the value of Lewis number. Impact of useful physical parameters has been illustrated graphically on various flow fields. It has inspected that the periodic magnetic field has helped to increase the interaction of the nanoparticles in the velocity field significantly. The field has been depicted in a vibrating form which is also done newly in this work. Subsequently, the Lorentz force has also represented a great impact in the updated visualization (streamlines and isotherms) of the flow field. The respective fields appeared with more wave for the larger values of magnetic parameter. These results help to visualize a theoretical idea of the effect of modern electromagnetic induction use in industry instead of traditional energy sources. Moreover, it has a great application in lung and prostate cancer therapy.


2009 ◽  
Vol 131 (3) ◽  
Author(s):  
Mario Durán ◽  
Jean-Claude Nédélec ◽  
Sebastián Ossandón

An efficient numerical method, using integral equations, is developed to calculate precisely the acoustic eigenfrequencies and their associated eigenvectors, located in a given high frequency interval. It is currently known that the real symmetric matrices are well adapted to numerical treatment. However, we show that this is not the case when using integral representations to determine with high accuracy the spectrum of elliptic, and other related operators. Functions are evaluated only in the boundary of the domain, so very fine discretizations may be chosen to obtain high eigenfrequencies. We discuss the stability and convergence of the proposed method. Finally we show some examples.


Author(s):  
David Roos Launchbury ◽  
Luca Mangani ◽  
Ernesto Casartelli ◽  
Francesco Del Citto

Abstract In the industrial simulation of flow phenomena, turbulence modeling is of prime importance. Due to their low computational cost, Reynolds-averaged methods (RANS) are predominantly used for this purpose. However, eddy viscosity RANS models are often unable to adequately capture important flow physics, specifically when strongly anisotropic turbulence and vortex structures are present. In such cases the more costly 7-equation Reynolds stress models often lead to significantly better results. Unfortunately, these models are not widely used in the industry. The reason for this is not mainly the increased computational cost, but the stability and convergence issues such models usually exhibit. In this paper we present a robust implementation of a Reynolds stress model that is solved in a coupled manner, increasing stability and convergence speed significantly compared to segregated implementations. In addition, the decoupling of the velocity and Reynolds stress fields is addressed for the coupled equation formulation. A special wall function is presented that conserves the anisotropic properties of the model near the walls on coarser meshes. The presented Reynolds stress model is validated on a series of semi-academic test cases and then applied to two industrially relevant situations, namely the tip vortex of a NACA0012 profile and the Aachen Radiver radial compressor case.


Author(s):  
Farhad Aghili

A heavy payload attached to the wrist force/moment (F/M) sensor of a manipulator can cause the conventional impedance controller to fail in establishing the desired impedance due to the noncontact components of the force measurement, i.e., the inertial and gravitational forces of the payload. This paper proposes an impedance control scheme for such a manipulator to accurately shape its force-response without needing any acceleration measurement. Therefore, no wrist accelerometer or a dynamic estimator for compensating the inertial load forces is required. The impedance controller is further developed using an inner/outer loop feedback approach that not only overcomes the robot dynamics uncertainty, but also allows the specification of the target impedance model in a general form, e.g., a nonlinear model. The stability and convergence of the impedance controller are analytically investigated, and the results show that the control input remains bounded provided that the desired inertia is selected to be different from the payload inertia. Experimental results demonstrate that the proposed impedance controller is able to accurately shape the impedance of a manipulator carrying a relatively heavy load according to the desired impedance model.


2012 ◽  
Vol 433-440 ◽  
pp. 6795-6801
Author(s):  
Xue Gui Zhu ◽  
Zhi Hong Fu ◽  
Xing Zhe Hou

Some modifications of conventional neuron proportional-integral-differential controller (NPID) are presented in this paper to prevent its slow dynamic response and loss of control in case of improper learning factors. The quasi-step signal replaces the step signal as the reference signal to improve the dynamic characteristics. The control output of NPID is modified every step by multiplying a penalty factor called senior teacher signal to suppress further the overshoot and compress the settling time. The steady-state error from the modified NPID (MNPID) is reduced or removed by adjusting dynamically reference input signal while excluding the pseudo steady state. Lots of simulation experiments are done to prove the stability and convergence of the MNPID control algorithm.


2013 ◽  
Vol 2013 ◽  
pp. 1-8 ◽  
Author(s):  
Abdon Atangana ◽  
Dumitru Baleanu

A kind of parabolic equation was extended to the concept of fractional calculus. The resulting equation is, however, difficult to handle analytically. Therefore, we presented the numerical solution via the explicit and the implicit schemes. We presented together the stability and convergence of this time-fractional parabolic equation with two difference schemes. The explicit and the implicit schemes in this case are stable under some conditions.


Sign in / Sign up

Export Citation Format

Share Document