This paper focuses on the automatic recognition of map projection, its inverse and re-projection. Our analysis leads to the unconstrained optimization solved by the hybrid BFGS nonlinear least squares technique. The objective function is represented by the squared sum of the residuals. For the map re-projection the partial differential equations of the inverse transformation are derived. They can be applied to any map projection. Illustrative examples of the stereographic and globular Nicolosi projections...

The common goal of systems pharmacology, i.e. systems biology applied to the field of pharmacology, is to rely less on trial and error in designing an input-output systems, e.g. therapeutic schedules. In this paper we present, on the paradigmatic example of a regulatory network of drug-induced enzyme production, the further development of the study published by Duintjer Tebbens et al. (2019) in the Applications of Mathematics. Here, the key feature is that the nonlinear model in form of an ODE system...

Semi–smooth Newton methods are analyzed for a class of variational inequalities in infinite dimensions. It is shown that they are equivalent to certain active set strategies. Global and local super-linear convergence are proved. To overcome the phenomenon of finite speed of propagation of discretized problems a penalty version is used as the basis for a continuation procedure to speed up convergence. The choice of the penalty parameter can be made on the basis of an $L^{\infty }$ estimate for the penalized...

Semi–smooth Newton methods are analyzed for a class of variational inequalities in infinite dimensions. It is shown that they are equivalent to certain active set strategies. Global and local super-linear convergence are proved. To overcome the phenomenon of finite speed of propagation of discretized problems a penalty version is used as the basis for a continuation procedure to speed up convergence. The choice of the penalty parameter can be made on the basis of an L∞ estimate for the penalized...

A model shape optimal design in ${ℝ}^2$ is solved by means of the penalty method with extrapolation, which enables to obtain high order approximations of both the state function and the boundary flux, thus offering a reliable gradient for the sensitivity analysis. Convergence of the proposed method is proved for certain subsequences of approximate solutions.

Within the range of Prandtl-Reuss model of elasto-plasticity the following optimal design problem is solved. Given body forces and surface tractions, a part of the boundary, where the (two-dimensional) body is fixed, is to be found, so as to minimize an integral of the squared yield function. The state problem is formulated in terms of stresses by means of a time-dependent variational inequality. For approximate solutions piecewise linear approximations of the unknown boundary, piecewise constant...

The state problem of elasto-plasticity (for the model with strain-hardening) is formulated in terms of stresses and hardening parameters by means of a time-dependent variational inequality. The optimal design problem is to find the shape of a part of the boundary such that a given cost functional is minimized. For the approximate solutions piecewise linear approximations of the unknown boundary, piecewise constant triangular elements for the stress and the hardening parameter, and backward differences...

A minimization of a cost functional with respect to a part of a boundary is considered for an elasto-plastic axisymmetric body obeying Hencky's law. The principle of Haar-Kármán and piecewise linear stress approximations are used to solve the state problem. A convergence result and the existence of an optimal boundary is proved.

A minimization of a cost functional with respect to a part of the boundary, where the body is fixed, is considered. The criterion is defined by an integral of a yield function. The principle of Haar-Kármán and piecewise constant stress approximations are used to solve the state problem. A convergence result and the existence of an optimal boundary is proved.

The paper is devoted to the problem of verification of accuracy of approximate solutions obtained in computer simulations. This problem is strongly related to a posteriori error estimates, giving computable bounds for computational errors and detecting zones in the solution domain where such errors are too large and certain mesh refinements should be performed. A mathematical model embracing nonlinear elliptic variational problems is considered in this work. Based on functional type estimates developed...

A numerical technique for solving the classical brachistochrone problem in the calculus of variations is presented. The brachistochrone problem is first formulated as a nonlinear optimal control problem. Application of this method results in the transformation of differential and integral expressions into some algebraic equations to which Newton-type methods can be applied. The method is general, and yields accurate results.

The present paper deals with the numerical solution of 3D shape optimization problems in frictional contact mechanics. Mathematical modelling of the Coulomb friction problem leads to an implicit variational inequality which can be written as a fixed point problem. Furthermore, it is known that the discretized problem is uniquely solvable for small coefficients of friction. Since the considered problem is nonsmooth, we exploit the generalized Mordukhovich’s differential calculus to compute the needed...

Degenerate parabolic variational inequalities with convection are solved by means of a combined relaxation method and method of characteristics. The mathematical problem is motivated by Richard’s equation, modelling the unsaturated – saturated flow in porous media. By means of the relaxation method we control the degeneracy. The dominance of the convection is controlled by the method of characteristics.

Degenerate parabolic variational inequalities with convection are solved by means of a combined relaxation method and method of characteristics. The mathematical problem is motivated by Richard's equation, modelling the unsaturated – saturated flow in porous media. By means of the relaxation method we control the degeneracy. The dominance of the convection is controlled by the method of characteristics.

A class of variational data assimilation problems on reconstructing the initial-value functions is considered for the models governed by quasilinear evolution equations. The optimality system is reduced to the equation for the control function. The properties of the control equation are studied and the solvability theorems are proved for linear and quasilinear data assimilation problems. The iterative algorithms for solving the problem are formulated and justified.

A class of variational data assimilation problems on reconstructing the initial-value functions is considered for the models governed by quasilinear evolution equations. The optimality system is reduced to the equation for the control function. The properties of the control equation are studied and the solvability theorems are proved for linear and quasilinear data assimilation problems. The iterative algorithms for solving the problem are formulated and justified.

We prove the existence of a sequence $\left(x_k\right)$ satisfying $0\in f\left(x_k\right)+{\sum }_{i=1}^M{a}_i\nabla f(x_k+{\beta }_i(x_{k+1}-x_k))(x_{k+1}-x_k)+F\left(x_{k+1}\right)$, where f is a function whose second order Fréchet derivative ∇²f satifies a center-Hölder condition and F is a set-valued map from a Banach space X to the subsets of a Banach space Y. We show that the convergence of this method is superquadratic.