High order generalized upwind schemes and numerical solution of singular perturbation problems

High even order generalizations of the traditional upwind method are introduced to solve second order ODE-BVPs without recasting the problem as a first order system. Both theoretical analysis and numerical comparison with central difference schemes of the same order show that these new methods may avoid typical oscillations and achieve high accuracy. Singular perturbation problems are taken into account to emphasize the main features of the proposed methods. AMS subject classification (2000)  65L10, 65L12, 65L50     

Mixed Elastico-Plasticity Problems with Partially Unknown Boundaries

In this paper,we study mixed elastico-plasticity problems in which part of the boundary is known,while the other part of the boundary is unknown and is a free boundary.Under certain conditions,this problemcan be transformed into a Riemann-Hilbert boundary value problem for analytic functions and a mixed boundaryvalue problem for complex equations.Using the theory of generalized analytic functions,the solvability of theproblem is discussed.     

The statement and numerical solution of an optimization problem in X-ray tomography

A nonclassical problem is considered for the transport equation with coefficients depending on the energy of radiation. The task is to find the discontinuity surfaces for the coefficients of the equation from measurements of the radiation flux leaving the medium. For this tomography problem, an optimization problem is stated and numerically analyzed. The latter consists in determining the radiation energy that ensures the best reconstruction of the unknown medium. A simplified optimization problem is solved analytically.     

Exceptional set of a representation with fractional powers

We prove that for any given c, 1 < c < 17/11, almost all natural numbers are representable in the form [x ^c] + [p ^c], where x is a natural number and p is a prime.     

Parametric Quintic-Spline Approach to the Solution of a System of Fourth-Order Boundary-Value Problems

In this paper, we use parametric quintic splines to derive some consistency relations which are then used to develop a numerical method for computing the solution of a system of fourth-order boundary-value problems associated with obstacle, unilateral, and contact problems. It is known that a class of variational inequalities related to contact problems in elastostatics can be characterized by a sequence of variational inequations, which are solved using some numerical method. Numerical evidence is presented to show the applicability and superiority of the new method over other collocation, finite difference, and spline methods.     

Mass‐lumping or not mass‐lumping for eigenvalue problems

In this article we analyze the effect of mass‐lumping in the linear triangular finite element approximation of second‐order elliptic eigenvalue problems. We prove that the eigenvalue obtained by using mass‐lumping is always below the one obtained with exact integration. For singular eigenfunctions, as those arising in non convex polygons, we prove that the eigenvalue obtained with mass‐lumping is above the exact eigenvalue when the mesh size is small enough. So, we conclude that the use of mass‐lumping is convenient in the singular case. When the eigenfunction is smooth several numerical experiments suggest that the eigenvalue computed with mass‐lumping is below the exact one if the mesh is not too coarse. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 653–664, 2003     

Shape Optimization and Fictitious Domain Approach for Solving Free Boundary Problems of Bernoulli Type

This contribution deals with an efficient method for the numerical realization of the exterior and interior Bernoulli free boundary problems. It is based on a shape optimization approach. The state problems are solved by a fictitious domain solver using boundary Lagrange multipliers.     

Iterative regularization methods for atmospheric remote sensing

In this paper we present different inversion algorithms for nonlinear ill-posed problems arising in atmosphere remote sensing. The proposed methods are Landweber's method (LwM), the iteratively regularized Gauss-Newton method, and the conventional and regularizing Levenberg-Marquardt method. In addition, some accelerated LwMs and a technique for smoothing the Levenberg-Marquardt solution are proposed. The numerical performance of the methods is studied by means of simulations. Results are presented for an inverse problem in atmospheric remote sensing, i.e., temperature sounding with an airborne uplooking high-resolution far-infrared spectrometer.     

Creation, Dissipation and Recycling of Resources in Non-Equilibrium Systems

In this article, we define stochastic dynamics for a system coupled to reservoirs. The rules for forward and backward transitions are related by a generalized detailed balance identity involving the system and its reservoirs. We compare the variation of information and of entropy. We define the Carnot dissipation and prove that it can be expressed in terms of cyclic transformations. Lower bounds for partial dissipations are also studied, as well as the effect of switching off certain reservoirs. We also study the near degeneracy of the stochastic matrix, relate it to phase transitions and we show that the reduced dynamics on the set of phases is a permutation. Finally, we relate these concepts to heat, work and more generally to the dissipation and creation of resources, in general systems.     

On quasi-solution to infeasible linear complementarity problem obtained by Lemke’s method

For a linear complementarity problem with inconsistent system of constraints a notion of quasi-solution of Tschebyshev type is introduced. It’s shown that this solution can be obtained automatically by Lemke’s method if the constraint matrix of the original problem is copositive plus or belongs to the intersection of matrix classes P ₀ and Q ₀.