Multigrid Solution of a Lavrentiev-Regularized State-Constrained Parabolic Control Problem

A mesh-independent, robust, and accurate multigrid scheme to solve a linear state-constrained parabolic optimal control problem is presented. We first consider a Lavrentiev regularization of the state-constrained optimization problem. Then, a multigrid scheme is designed for the numerical solution of the regularized optimality system. Central to this scheme is the construction of an iterative pointwise smoother which can be formulated as a local semismooth Newton iteration. Results of numerical experiments and theoretical two-grid local Fourier analysis estimates demonstrate that the proposed scheme is able to solve parabolic state-constrained optimality systems with textbook multigrid efficiency.     

p-Adic Analogue of the Wave Equation

In this paper, a p-adic analogue of the wave equation with Lipschitz source is considered. Since it is hard to arrive the solution of the problem, we propose a regularized method to solve the problem from a modified p-adic integral equation. Moreover, we give an iterative scheme for numerical computation of the regularlized solution.

     

Solution of a regularized problem for a stationary magnetic field in a nonhomogeneous conducting medium by a finite element method

The use of a vector finite element method for solving a regularized stationary magnetic field problem, which is formulated in terms of a vector magnetic potential, is validated. A generalized solutionis approximated using first-order Nedelec vector elements of the second kind on tetrahedrons. The existence and uniqueness of the solution to a discrete regularized problem and its convergence to the generalized solution for the case of an inhomogeneous (in the electromagnetic properties) domain are justified. Some issues of the numerical solution to the discrete regularized problem are discussed. Approaches to optimize the algorithms are shown on a series of numerical experiments.     

A Tikhonov-based projection iteration for nonlinear Ill-posed problems with sparsity constraints

In this paper, we consider nonlinear inverse problems where the solution is assumed to have a sparse expansion with respect to a preassigned basis or frame. We develop a scheme which allows to minimize a Tikhonov functional where the usual quadratic regularization term is replaced by a one-homogeneous (typically weighted ℓ _p ) penalty on the coefficients (or isometrically transformed coefficients) of such expansions. For (p < 2), the regularized solution will have a sparser expansion with respect to the basis or frame under consideration. The computation of the regularized solution amounts in our setting to a Landweber-fixed-point iteration with a projection applied in each fixed-point iteration step. The performance of the resulting numerical scheme is demonstrated by solving the nonlinear inverse single photon emission computerized tomography (SPECT) problem.     

A dynamic viscoelastic contact problem with normal compliance,finite penetration and nonmonotone slip rate dependent friction

We consider a mathematical model which describes the dynamic evolution of a viscoelastic body in frictional contact with an obstacle. The contact is modelled with normal compliance and unilateral constraint, associated to a rate slip-dependent version of Coulomb’s law of dry friction. In order to approximate the contact conditions, we consider a regularized problem wherein the contact is modelled by a standard normal compliance condition without finite penetrations. For each problem, we derive a variational formulation and an existence result of the weak solution of the regularized problem is obtained. Next, we prove the convergence of the weak solution of the regularized problem to the weak solution of the initial nonregularized problem. Then, we introduce a fully discrete approximation of the variational problem based on a finite element method and on a second order time integration scheme. The solution of the resulting nonsmooth and nonconvex frictional contact problems is presented, based on approximation by a sequence of nonsmooth convex programming problems. Finally, some numerical simulations are provided in order to illustrate both the behaviour of the solution related to the frictional contact conditions and the convergence result.     

Regularized state-constrained boundary optimal control of the Navier-Stokes equations

The numerical solution of the Dirichlet boundary optimal control problem of the Navier-Stokes equations in presence of pointwise state constraints is investigated. Two different regularization techniques are considered. First, a Moreau-Yosida regularization of the problem is studied. Optimality conditions are derived and the convergence of the regularized solutions towards the original one is proved. A source representation of the control combined with a Lavrentiev type regularization strategy is also presented. The analysis concerning optimality conditions and convergence of the regularized solutions is carried out. In the last part of the paper numerical experiments are presented. For the numerical solution of each regularized problem a semi-smooth Newton method is applied.     

Semi-smooth Newton methods for the Signorini problem

Semi-smooth Newton methods are analyzed for the Signorini problem. A proper regularization is introduced which guarantees that the semi-smooth Newton method is superlinearly convergent for each regularized problem. Utilizing a shift motivated by an augmented Lagrangian framework, to the regularization term, the solution to each regularized problem is feasible. Convergence of the regularized problems is shown and a report on numerical experiments is given.     

A Regularized Newton Method in Inverse Scattering Problem from Cavities

We consider the inverse problem to determine the shape of a open cavity embedded in the infinite ground plane from knowledge of the far-field pattern of the scattering of TM polarization.For its approximate solution we propose a regularized Newton iteration scheme.For a foundation of Newton type methods we establish the Fréchet differentiability of solution to the scattering problem with respect to the boundary of the cavity.Some numerical examples of the feasibility of the method are presented.     

On the resolution of singularities of ordinary differential systems

We show how some differential geometric ideas help to resolve some singularities of ordinary differential systems. Hence a singular problem is replaced by a regular one, which facilitates further analysis of the system. The methods employed are constructive and the regularized systems can also be used for numerical computations. This revised version was published online in June 2006 with corrections to the Cover Date.     

A combined BDF‐semismooth Newton approach for time‐dependent Bingham flow

This article is devoted to the numerical simulation of time‐dependent convective Bingham flow in cavities. Motivated by a primal‐dual regularization of the stationary model, a family of regularized time‐dependent problems is introduced. Well posedness of the regularized problems is proved, and convergence of the regularized solutions to a solution of the original multiplier system is verified. For the numerical solution of each regularized multiplier system, a fully discrete approach is studied. A stable finite element approximation in space together with a second‐order backward differentiation formula for the time discretization are proposed. The discretization scheme yields a system of Newton differentiable nonlinear equations in each time step, for which a semismooth Newton algorithm is used. We present two numerical experiments to verify the main properties of the proposed approach. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

On a terminal value problem for a generalization of the fractional diffusion equation with hyper-Bessel operator

In this paper, we consider an inverse problem of recovering the initial value for a generalization of time-fractional diffusion equation, where the time derivative is replaced by a regularized hyper-Bessel operator. First, we investigate the existence and regularity of our terminal value problem. Then we show that the backward problem is ill-posed, and we propose a regularizing scheme using a fractional Tikhonov regularization method. We also present error estimates between the regularized solution and the exact solution using two parameter choice rules.     

Control analysis for water quality problem

The optimal control of water quality model is considered. We analyze a regularized version of the pointwise control. A numerical algorithm is presented. The solution approach is a combination of a TVD scheme for the shallow water equations and a purley deterministic particle method for the pollutant problem. We present some computational results. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

解洞穴逆散射问题的一种正则化Newton方法

We consider the inverse problem to determine the shape of a open cavity embedded in the infinite ground plane from knowledge of the far-field pattern of the scattering of TM polarization.For its approximate solution we propose a regularized Newton iteration scheme.For a foundation of Newton type methods we establish the Fréchet differentiability of solution to the scattering problem with respect to the boundary of the cavity.Some numerical examples of the feasibility of the method are presented.     

Modeling and analysis of the unilateral contact of a piezoelectric body with a conductive support

We consider a mathematical model which describes the quasistatic process of contact between a piezoelectric body and an electrically conductive support, the so-called foundation. We model the material's behavior with a nonlinear electro-viscoelastic constitutive law; the contact is frictionless and is described with the Signorini condition and a regularized electrical conductivity condition. We derive a variational formulation for the problem and then we prove the existence of a unique weak solution to the model. The proof is based on arguments of nonlinear equations with multivalued maximal monotone operators and fixed point. Then we introduce a fully discrete scheme, based on the finite element method to approximate the spatial variable and the backward Euler scheme to discretize the time derivatives. We treat the unilateral contact conditions by using an augmented Lagrangian approach. We implement this scheme in a numerical code then we present numerical simulations in the study of two-dimensional test problems, together with various comments and interpretations.     

Lagrange multipliers, (exact) regularization and error bounds for monotone variational inequalities

We examine two central regularization strategies for monotone variational inequalities, the first a direct regularization of the operative monotone mapping, and the second via regularization of the associated dual gap function. A key link in the relationship between the solution sets to these various regularized problems is the idea of exact regularization, which, in turn, is fundamentally associated with the existence of Lagrange multipliers for the regularized variational inequality. A regularization is said to be exact if a solution to the regularized problem is a solution to the unregularized problem for all parameters beyond a certain value. The Lagrange multipliers corresponding to a particular regularization of a variational inequality, on the other hand, are defined via the dual gap function. Our analysis suggests various conceptual, iteratively regularized numerical schemes, for which we provide error bounds, and hence stopping criteria, under the additional assumption that the solution set to the unregularized problem is what we call weakly sharp of order greater than one.     

Stability and consistency of the semi-implicit co-volume scheme for regularized mean curvature flow equation in level set formulation

We show stability and consistency of the linear semi-implicit complementary volume numerical scheme for solving the regularized, in the sense of Evans and Spruck, mean curvature flow equation in the level set formulation. The numerical method is based on the finite volume methodology using the so-called complementary volumes to a finite element triangulation. The scheme gives the solution in an efficient and unconditionally stable way.     

Solvability of a dynamic contact problem between a piezoelectric body and a conductive foundation

We consider a mathematical model which describes the dynamic process of contact between a piezoelectric body and an electrically conductive foundation. We model the material’s behavior with a nonlinear electro-viscoelastic constitutive law; the contact is frictionless and is described with the normal compliance condition and a regularized electrical conductivity condition. We derive a variational formulation for the problem and then, under a smallness assumption on the data, we prove the existence of a unique weak solution to the model. We also investigate the behavior of the solution with respect the electric data on the contact surface and prove a continuous dependence result. Then, we introduce a fully discrete scheme, based on the finite element method to approximate the spatial variable and the backward Euler scheme to discretize the time derivatives. We treat the contact by using a penalized approach and a version of Newton’s method. We implement this scheme in a numerical code and, in order to verify its accuracy, we present numerical simulations in the study of two-dimensional test problems. These simulations provide a numerical validation of our continuous dependence result and illustrate the effects of the conductivity of the foundation, as well.     

Numerical scheme for multilayer shallow-water model in the low-Froude number regime

The aim of this note is to present a multi-dimensional numerical scheme approximating the solutions to the multilayer shallow-water model in the low-Froude-number regime. The proposed strategy is based on a regularized model where the advection velocity is modified with a pressure gradient in both mass and momentum equations. The numerical solution satisfies the dissipation of energy, which acts for mathematical entropy, and the main physical properties required for simulations within oceanic flows.     

An inverse boundary value problem for the Oseen equation

Based on the two‐dimensional stationary Oseen equation we consider the problem to determine the shape of a cylindrical obstacle immersed in a fluid flow from a knowledge of the fluid velocity on some arc outside the obstacle. First, we obtain a uniqueness result for this ill‐posed and non‐linear inverse problem. Then, for the approximate solution we propose a regularized Newton iteration scheme based on a boundary integral equation of the first kind. For a foundation of Newton‐type methods we establish the Fréchet differentiability of the solution to the Dirichlet problem for the Oseen equation with respect to the boundary and investigate the injectivity of the linearized mapping. Some numerical examples for the feasibility of the method are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

THE QUASI-BOUNDARY VALUE METHOD FOR IDENTIFYING THE INITIAL VALUE OF THE SPACE-TIME FRACTIONAL DIFFUSION EQUATION

In this article, we consider to solve the inverse initial value problem for an inhomogeneous space-time fractional diffusion equation. This problem is ill-posed and the quasi-boundary value method is proposed to deal with this inverse problem and obtain the series expression of the regularized solution for the inverse initial value problem. We prove the error estimates between the regularization solution and the exact solution by using an a priori regularization parameter and an a posteriori regularization parameter choice rule. Some numerical results in one-dimensional case and two-dimensional case show that our method is effcient and stable.