Convergence of a discrete scheme in a regularization method for the quasi-steady Maxwell system in a nonhomogeneous conducting medium

Convergence of the discrete solution to the solution of a regularized system of the Maxwell equations written in terms of a vector magnetic potential with a special calibration of medium’s conduction is considered. The problem is discretized by a Nedelec vector finite element method in space and by an implicit Euler scheme in time. An optimal theoretical energy estimate of the approximate solution error in 3D Lipschitz polyhedral domains is obtained.     

A regularization method for stationary maxwell equations in an inhomogeneous conducting medium

We consider the problem of defining a vector potential of a magnetic field with a nonstandard calibration in an inhomogeneous conducting medium. The problem in question is one with constraints on the right-hand side and on the solution itself. A generalized and regularized statement of this problem without constraints is proposed and substantiated. Such a problem statement is equivalent to the original generalized problem with constraints.     

多维广义SRLW方程的Chebyshev拟谱方法分析

考虑了一类多维的广义对称正则长波(SRLW)方程的齐次初边值问题Chebyshev拟谱逼近,构造了全离散的Chebyshev拟谱格式,给出了这种格式近似解的收敛性和最优误差估计。     

Mixed finite element analysis for dissipative SRLW equations with damping term

In this paper a mixed finite element (MFE) formulation is proposed for the initial-boundary value problem of dissipative symmetric regularized long wave (SRLW) equations with damping. Existence and uniqueness of its generalized solution and of the fully discrete mixed finite element solution are proved. Error estimates based on energy methods are given. Numerical experiments verify the theoretical analysis.     

用3个向量对构造梁振动系统的刚度矩阵

针对梁的离散化模型的刚度矩阵是五对角矩阵,梁振动反问题的实质是实对称五对角矩阵的特征值反问题.该文利用向量对、Moore-Penrose广义逆给出了实对称五对角矩阵向量对反问题存在唯一解的条件,并结合矩阵分块讨论了双对称五对角矩阵向量对反问题解存在唯一的条件,进而计算了次对角线位置元素为负,其它位置元素均为正的实对称五对角矩阵特征值反问题.由于构造梁的离散模型需要的数据可由测试得到,故而其结果适合于模态分析、系统结构的分析与设计等方面应用.最后给出了数值算例,通过数值讨论说明方法的有效性.     

Regularization method for solving the quasi-stationary Maxwell equations in an inhomogeneous conducting medium

Nedelec vector finite elements are used for the numerical solution of a regularized version of the quasi-stationary Maxwell equations written in terms of a scalar and a vector magnetic potential with special calibration taking into account the conductivity of the medium. An optimal energy estimate for the error of the approximate solution in Lipschitz polyhedral domains is established. Numerical results are presented that demonstrate the stability of the method.     

A generalized global Arnoldi method for ill-posed matrix equations

This paper discusses the solution of large-scale linear discrete ill-posed problems with a noise-contaminated right-hand side. Tikhonov regularization is used to reduce the influence of the noise on the computed approximate solution. We consider problems in which the coefficient matrix is the sum of Kronecker products of matrices and present a generalized global Arnoldi method, that respects the structure of the equation, for the solution of the regularized problem. Theoretical properties of the method are shown and applications to image deblurring are described.     

Optimal control of Maxwell’s equations with regularized state constraints

This paper is devoted to an optimal control problem of Maxwell??s equations in the presence of pointwise state constraints. The control is given by a divergence-free three-dimensional vector function representing an applied current density. To cope with the divergence-free constraint on the control, we consider a vector potential ansatz. Due to the lack of regularity of the control-to-state mapping, existence of Lagrange multipliers cannot be guaranteed. We regularize the optimal control problem by penalizing the pointwise state constraints. Optimality conditions for the regularized problem can be derived straightforwardly. It also turns out that the solution of the regularized problem enjoys higher regularity which then allows us to establish its convergence towards the solution of the unregularized problem. The second part of the paper focuses on the numerical analysis of the regularized optimal control problem. Here the state and the control are discretized by Nédélec??s curl-conforming edge elements. Employing the higher regularity property of the optimal control, we establish an a priori error estimate for the discretization error in the $\boldsymbol{H}(\bold{curl})$ -norm. The paper ends by numerical results including a numerical verification of our theoretical results.     

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.     

Prox-Regularization Methods for Generalized Fractional Programming

If a fractional program does not have a unique solution or the feasible set is unbounded, numerical difficulties can occur. By using a prox-regularization method that generates a sequence of auxiliary problems with unique solutions, these difficulties are avoided. Two regularization methods are introduced here. They are based on Dinkelbach-type algorithms for generalized fractional programming, but use a regularized parametric auxiliary problem. Convergence results and numerical examples are presented.     

The unconditional stable and convergent difference methods with intrinsic parallelism for quasilinear parabolic systems

A kind of the general finite difference schemes with intrinsic parallelism forthe boundary value problem of the quasilinear parabolic system is studied without assum-ing heuristically that the original boundary value problem has the unique smooth vectorsolution. By the method of a priori estimation of the discrete solutions of the nonlineardifference systems, and the interpolation formulas of the various norms of the discretefunctions and the fixed-point technique in finite dimensional Euclidean space, the exis-tence and uniqueness of the discrete vector solutions of the nonlinear difference systemwith intrinsic parallelism are proved. Moreover the unconditional stability of the generalfinite difference schemes with intrinsic parallelism is justified in the sense of the continu-ous dependence of the discrete vector solution of the difference schemes on the discretedata of the original problems in the discrete w_2~(2,1) norms. Finally the convergence of thediscrete vector solutions of the certain differe     

General difference schemes with intrinsic parallelism for nonlinear parabolic systems

The boundary value problem for nonlinear parabolic system is solved by the finite difference method with intrinsic parallelism. The existence of the discrete vector solution for the general finite difference schemes with intrinsic parallelism is proved by the fixed-point technique in finite-dimensional Euclidean space. The convergence and stability theorems of the discrete vector solutions of the nonlinear difference system with intrinsic parallelism are proved. The limitation vector function is just the unique generalized solution of the original problem for the parabolic system.     

Approximation of an Inverse Initial Problem for a Biparabolic Equation

In this paper, we consider the problem of finding the initial distribution for the linear inhomogeneous and nonlinear biparabolic equation. The problem is severely ill-posed in the sense of Hadamard. First, we apply a general filter method to regularize the linear nonhomogeneous problem. Then, we also give a regularized solution and consider the convergence between the regularized solution and the sought solution. Under the a priori assumption on the exact solution belonging to a Gevrey space, we consider a generalized nonlinear problem by using the Fourier truncation method to obtain rigorous convergence estimates in the norms on Hilbert space and Hilbert scale space.     

Regularization of Discrete Ill-Posed Problems

The paper concerns conditioning aspects of finite-dimensional problems arising when the Tikhonov regularization is applied to discrete ill-posed problems. A relation between the regularization parameter and the sensitivity of the regularized solution is investigated. The main conclusion is that the condition number can be decreased only to the square root of that for the nonregularized problem. The convergence of solutions of regularized discrete problems to the exact generalized solution is analyzed just in the case when the regularization corresponds to the minimal condition number. The convergence theorem is proved under the assumption of the suitable relation between the discretization level and the data error. As an example the method of truncated singular value decomposition with regularization is considered. This revised version was published online in July 2006 with corrections to the Cover Date.     

UNCONDITIONAL STABLE DIFFERENCE METHODS WITH INTRINSIC PARALLELISM FOR SEMILINEAR PARABOLIC SYSTEMS OF DIVERGENCE TYPE

The general finite difference schemes with intrinsic parallelism for the boundary value problem of the semilinear parabolic system of divergence type with bounded measurable coefficients is studied. By the approach of the discrete functional analysis, the existence and uniqueness of the discrete vector solutions of the nonlinear difference system with intrinsic parallelism are proved. Moreover the unconditional stability of the general difference schemes with intrinsic parallelism justified in the sense of the continuous dependence of the discrete vector solution of the difference schemes on the discrete initial data of the original problems in the discrete W_2~(2,1) (Q△) norms. Finally the convergence of the discrete vector solutions of the certain difference schemes with intrinsic parallelism to the unique generalized solution of the original semilinear parabolic problem is proved.     

Finite elements for parabolic equations backwards in time

Summary Several regularization methods for parabolic equations backwards in time together with the usual additional constraints for their solution are considered. The error of the regularization is estimated from above and below. For a boundary value problem in time-method, finite elements as well as a time discretization are introduced and the error with respect to the regularized solution is estimated, thus giving an overall error of the discrete regularized problem. The algorithm is tested in simple numerical examples.     

The regularizing properties of anisotropic radial basis functions

In the present work we consider the problem of interpolating scattered data using radial basis functions (RBF). In general, it is well known that this leads to a discrete linear inverse problem that needs to be regularized in order to provide a meaningful solution. The work focuses on a metric-regularization approach, based on a new class of RBF, called anisotropic RBF. The work provides theoretical justifications for the regularization approach and it considers a suitable proposal for the metric, supporting it by numerical examples.     

Pointwise error estimates for a generalized Oseen problem and an application to an optimal control problem

We derive pointwise error estimates for a generalized Oseen when it is approximated by a low order Taylor‐Hood finite element scheme in two dimensions. The analysis is based on estimates for regularized Green's functions associated with a generalized Oseen problem on weighted Sobolev spaces and weighted interpolation results. We apply the maximum norm results to obtain convergence in an optimal control problem governed by a generalized Oseen equation and present a numerical example that allows us to show the behavior of the error.     

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     

Fully discrete T‐ψ finite element method to solve a nonlinear induction hardening problem

We study an induction hardening model described by Maxwell's equations coupled with a heat equation. The magnetic induction field is assumed a nonlinear constitutional relation and the electric conductivity is temperature‐dependent. The T ‐ψ method is to transform Maxwell's equations to the vector–scalar potential formulations and to solve the potentials by means of the finite element method. In this article, we present a fully discrete T ‐ψ finite element scheme for this nonlinear coupled problem and discuss its solvability. We prove that the discrete solution converges to a weak solution of the continuous problem. Finally, we conclude with two numerical experiments for the coupled system.