Convergence analysis on Browder-Tikhonov regularization for second-order evolution hemivariational inequality

This paper studies the Browder-Tikhonov regularization of a second-order evolution hemivariational inequality (SOEHVI) with non-coercive operators. With duality mapping, the regularized formulations and a derived first-order evolution hemivariational inequality (FOEHVI) for the problem considered are presented. By applying the Browder-Tikhonov regularization method to the derived FOEHVI, a sequence of regularized solutions to the regularized SOEHVI is constructed, and the strong convergence of the whole sequence of regularized solutions to a solution to the problem is proved.     

An inverse problem formulation of the immersed-boundary method

We formulate the immersed-boundary method (IBM) as an inverse problem. A control variable is introduced on the boundary of a larger domain that encompasses the target domain. The optimal control is the one that minimizes the mismatch between the state and the desired boundary value along the immersed target-domain boundary. We begin by investigating a naïve problem formulation that we show is ill-posed: in the case of the Laplace equation, we prove that the solution is unique, but it fails to depend continuously on the data; for the linear advection equation, even solution uniqueness fails to hold. These issues are addressed by two complimentary strategies. The first strategy is to ensure that the enclosing domain tends to the true domain, as the mesh is refined. The second strategy is to include a specialized parameter-free regularization that is based on penalizing the difference between the control and the state on the boundary. The proposed inverse IBM is applied to the diffusion, advection, and advection-diffusion equations using a high-order discontinuous Galerkin discretization. The numerical experiments demonstrate that the regularized scheme achieves optimal rates of convergence and that the reduced Hessian of the optimization problem has a bounded condition number, as the mesh is refined.     

Regularization in 3D for anisotropic elastodynamic crack and obstacle problems

We propose, in this paper, a unified method of generating a regularized integral equation in the double layer potential approach for 3D anisotropic elastodynamics. Our regularization preserves the causality in the time-domain. The method is based on a special decomposition of the hypersingular kernel which appears in the integral representation of the stress tensor.     

A new positive‐definite regularization of incompressible Navier–Stokes equations discretized with Q1/P0 finite element

A new regularization method is proposed for the Galerkin approximation of the incompressible Navier–Stokes equations with Q1/P0 element, by newly introducing a square‐type linear form into the variational divergence‐free constraint regularized with the global pressure jump (GPJ) method. The addition of the square‐type linear form is intended to eliminate the hydrostatic pressure mode appearing in confined flows, and to make the discretized matrix positive definite and then non‐singular without the pressure pegging trick. Effects of the free parameters for the regularization on the solutions are numerically examined with a 2‐D driven cavity flow problem. Furthermore, the convergences in the conjugate gradient iteration for the solution of the pressure Poisson equation are compared among the mixed method, the GPJ method and the present method for both leaky and non‐leaky 3‐D driven cavity flows. Finally, the non‐leaky 3‐D cavity flows at different Re numbers are solved to compare with the literature data and to demonstrate the accuracy of the proposed method. Copyright © 2003 John Wiley & Sons, Ltd.     

Initial-Boundary value problems for the Boltzmann equation: Global existence of weak solutions

For an open set of ³ bounded or not, we consider initial-boundary value problems for the Boltzmann equation. For general gas-surface interaction laws and for hard potentials, we prove a global existence result for weak solutions. The proof uses the regularization of the collision operator and the renormalization method for the regularized problem. By using weak compactness in L¹ and averaged stability ofQ(f,f), we prove the existence of weak solutions of our problem.Dedicated to the Memory of Ronald DiPerna     

Regularization of Neutral Delay Differential Equations with Several Delays

For neutral delay differential equations the right-hand side can be multi-valued, when one or several delayed arguments cross a breaking point. This article studies a regularization via a singularly perturbed problem, which smooths the vector field and removes the discontinuities in the derivative of the solution. A low-dimensional dynamical system is presented, which characterizes the kind of generalized solution that is approximated. For the case that the solution of the regularized problem has high frequency oscillations around a codimension-2 weak solution of the original problem, a new stabilizing regularization is proposed and analyzed.     

Damage detection using a new regularization method with variable parameter

In this research, a sensitivity approach to finite element model updating is used to determine stiffness reduction factors from measured structural response. The used method causes a set of nonlinear ill-conditioned equations that need to be linearized and regularized in order to find the solution. A new approach to solve the problem is presented using variable regularization parameter. Utilization of variable regularization parameter eliminates dependency on the number of iterations and prevents the loss of regularization effect due to iterations. A new stopping criteria is used which is based on the difference between mean and variance of last iterations. Furthermore the results show that using wavelet transform to update the model yields better results than modal parameters. Expedient performance of the proposed method is shown through a numerical simulation.     

网络RTK参考站间模糊度动态解算中病态方程的一种解算方法

在网络RTK参考站间的模糊度估计中,若误差方程严重病态,将导致模糊度解与其准确值偏差较大或整周模糊度无法固定,因此提出了一种适于网络RTK模糊度动态解算的新方案：1）法方程病态性的判断;2）Tikhonov正则化解算病态方程;3）LAMBDA方法搜索固定整周模糊度。同时,深入研究了Tikhonov正则化矩阵的构造方法和正则化参数的选取准则。最后以实例验证了采用此方案解算病态方程是可行的,通过选取合适的正则化参数可以解得准确的整周模糊度;详细讨论了选择不同的正则化参数对模糊度解算结果的影响。     

Some results concerning approximation of regularized compressible flow

We summarize here some theoretical results for fictitious gas regularization of compressible flow and give error estimates for the finite element approximation to the regularized problem.     

Identification of a moving boundary for a heat conduction problem in a multilayer medium

In this paper, we give a uniqueness theorem for the moving boundary of a heat problem in a composite medium. Through solving the Cauchy problem of heat equation in each subdomain, we finally find an approximation to the moving boundary for one-dimensional heat conduction problem in a multilayer medium. The numerical scheme is based on the use of the method of fundamental solutions and a discrete Tikhonov regularization technique with the generalized cross-validation choice rule for a regularization parameter. Numerical experiments for five examples show that our proposed method is effective and stable.     

On one class of differential equations of fractional order

We consider the Darboux problem for a differential equation of fractional order that contains a regularized mixed derivative. Sufficient conditions for the existence and uniqueness of a solution of this problem are obtained in the class of continuous functions. We also propose a method for finding an approximate solution of this problem and prove the convergence of this method.     

A Space-Time PGD-DIC Algorithm:

The aim of this study is to develop a new regularized Digital Image Correlation (DIC) method for time dependent measurements. The correlation problem is written as a minimization problem over the space-time domain in a general formulation including 2D-DIC and Stereo DIC (SDIC). The unknown time-resolved displacement field is found as a sum of products of space and time functions, similarly to the Proper Generalized Decomposition in computational mechanics. It is shown that the space fields are less sensitive to noise as time regularity acts as a physical regularization of the space fields. The proposed method is illustrated by vibration measurement under harmonic excitation in 2D-DIC and SDIC.     

Degenerate Regularization of Forward–Backward Parabolic Equations: The Regularized Problem

We study a quasilinear parabolic equation of forward–backward type in one space dimension, under assumptions on the nonlinearity which hold for a number of important mathematical models (for example, the one-dimensional Perona–Malik equation), using a degenerate pseudoparabolic regularization proposed in Barenblatt et al. (SIAM J Math Anal 24:1414–1439, 1993), which takes time delay effects into account. We prove existence and uniqueness of positive solutions of the regularized problem in a space of Radon measures. We also study qualitative properties of such solutions, in particular concerning their decomposition into an absolutely continuous part and a singular part with respect to the Lebesgue measure. In this respect, the existence of a family of viscous entropy inequalities plays an important role.     

Distributional and regularized radiation fields of non-uniformly moving straight dislocations,and elastodynamic Tamm problem

This work introduces original explicit solutions for the elastic fields radiated by non-uniformly moving, straight, screw or edge dislocations in an isotropic medium, in the form of time-integral representations in which acceleration-dependent contributions are explicitly separated out. These solutions are obtained by applying an isotropic regularization procedure to distributional expressions of the elastodynamic fields built on the Green tensor of the Navier equation. The obtained regularized field expressions are singularity-free, and depend on the dislocation density rather than on the plastic eigenstrain. They cover non-uniform motion at arbitrary speeds, including faster-than-wave ones. A numerical method of computation is discussed, that rests on discretizing motion along an arbitrary path in the plane transverse to the dislocation, into a succession of time intervals of constant velocity vector over which time-integrated contributions can be obtained in closed form. As a simple illustration, it is applied to the elastodynamic equivalent of the Tamm problem, where fields induced by a dislocation accelerated from rest beyond the longitudinal wave speed, and thereafter put to rest again, are computed. As expected, the proposed expressions produce Mach cones, the dynamic build-up and decay of which is illustrated by means of full-field calculations.     

Numerical Flow Simulation for Bingham Plastics in a Single-Screw Extruder

Numerical simulations have been performed concerning the operation of a single-screw extruder, pumping a Bingham plastic under isothermal, developed flow conditions. Under the assumption of sufficiently low Reynolds numbers, inertia effects are neglected. The singular rheological behavior of the Bingham plastic is considered as the limiting case within a class of generalized Newtonian liquids with smooth constitutive equations. The validation of this regularization process is shown for a related flow problem where the Bingham solution is known analytically. A mixed finite-element method is applied to the flow in the screw-extruder to reduce the equations of motion, the continuity equation, and the regularized constitutive equation to a set of nonlinear algebraic equations, which are solved using a Newton method. In particular, the pumping characteristics of a given screw geometry are extracted from the finite-element calculations, i.e., the dependence of the volumetric flow rate and of the power requirement on the axial pressure drop, on the screw speed, and on the rheological parameters. Calculated flow fields clearly show the size and position of regions in the extruder channel where the Bingham plastic behaves like a solid. Received: 12 December 1995 and accepted 12 November 1996     

The concept of material forces in phase transition problems within the level-set framework

The dynamics of a phase transition front in solids using the level set method is examined in this paper. Introducing an implicit representation of singular surfaces, a regularized version of the sharp interface model arises. The interface transforms into a thin transition layer of nonzero thickness where all quantities take inhomogeneous expressions within the body. It is proved that the existence of an inhomogeneous energy of the material predicts inhomogeneity forces that drive the singularity. The driving force is a material force entering the canonical momentum equation (pseudo-momentum) in a natural way. The evolution problem requires a kinetic relation that determines the velocity of the phase transition as a function of the driving force. Here, the kinetic relation is produced by invoking relations that can be considered as the regularized versions of the Rankine–Hugoniot jump conditions. The effectiveness of the method is illustrated in a shape memory alloy bar.     

A new method for detecting non-linear damping and restoring forces in non-linear oscillation systems from transient data

The purpose of this study is to recover the functional form of both non-linear damping and non-linear restoring forces in the non-linear oscillatory motions of an autonomous system. Using two sets of measured motion response data of the system, an inverse problem is formulated for recovering (or identification): the differential equation of motion is transformed into an equivalent integral equation of motion. The identification, which is non-linear, is shown to be one-to-one. However, the inverse problem formulated herein is concerned with the Volterra-type of non-linear integral equation of the first kind. This leads to numerical instability: solutions of the inverse problem lack stability properties. In order to overcome the difficulty, a regularization method is applied to the identification process. In addition, an L-curve criterion, combined with regularization, is introduced to find an optimal choice for the regularization parameter (i.e., the number of iterations), in the presence of noisy data. The workability of the identification is investigated for simultaneously recovering the functional form of the non-linear damping and the non-linear restoring forces through a numerical experiment.     

半空间SH波方程非线性井间双参数反演

以半空间的ＳＨ波方程出发,采用Ｂｏｒｎ迭代法求解半空间弹性介质中密度和剪切模量分布的非线性反演问题。首先,采用矩量法和正则化方法,给出井间反演积分方程的离散形式,然后应用Ｂｒｏｎ迭代法求解非线性反演问题。     

Regularized formulation of the variational brittle fracture with unilateral contact: Numerical experiments

This paper presents a modified regularized formulation of the Ambrosio-Tortorelli type to introduce the crack non-interpenetration condition in the variational approach to fracture mechanics proposed by Francfort and Marigo [1998. Revisiting brittle fracture as an energy minimization problem. J. Mech. Phys. Solids 46 (8), 1319-1342]. We focus on the linear elastic case where the contact condition appears as a local unilateral constraint on the displacement jump at the crack surfaces. The regularized model is obtained by splitting the strain energy in a spherical and a deviatoric parts and accounting for the sign of the local volume change. The numerical implementation is based on a standard finite element discretization and on the adaptation of an alternate minimization algorithm used in previous works. The new regularization avoids crack interpenetration and predicts asymmetric results in traction and in compression. Even though we do not exhibit any gamma-convergence proof toward the desired limit behavior, we illustrate through several numerical case studies the pertinence of the new model in comparison to other approaches.     

Relaxation time spectra from short frequency range small-angle dynamic rheometry

In this work, the influence of the experimental frequency range over the relaxation time spectrum is studied. The relaxation time spectra were calculated from dynamic moduli, a well-known ill-posed problem, using a regularization method. The method solves the ill-posed problem by simultaneous minimization of the regularized standard deviation and a restriction function. The solution was validated using a simulated spectrum. Truncated moduli data generated from simulated spectra were used to evaluate the method for smaller frequency range data. Finally, experimental data of a wormlike micellar system mixed in aqueous solution with a zwitterionic copolymer were used to validate the method. It was possible to obtain relaxation time spectra from short frequency range data if the relaxation time range is allowed to be higher than the inverse of the highest and lowest experimental frequencies. These spectra can be used qualitatively to describe complex systems when no time-temperature superposition experiments are feasible.