Numerical Analysis of Nonstationary Thermistor Problem

1.AMathematicalModelandaDiscreteSchemeThemodelofanonstationaxythermistorproblemisderivedfromtheconservationlawsofcurrentandenergy(see[l][2]l3]):Findapair{W,u}suchthatV-(a(u)VW)=oinQT=flX(O,T),(1.l)W=Woonoflx(O,T),(1.2)ut-bu=a(U)IVWl'inQT,(1'3)u=oonoflx(O,T),(1.5)u(x,o)=uo(x)infl(1.5)whereflCR"(N21)isaboundeddomain,occupiedbythethermistor;W=yt(x,t),u=u(x,t)aredistributionsoftheelectricalpotentialandthetemperatureinfl,respectively;J(u)isthetemperaturedependentelectricalconductivity;a…     

热敏电阻问题的有限差分解法及其收敛性分析

The thermistor problem is an initial-boundary value problem of coupled nonlinear differential equations.The nonlinear PDEs consist of a heat equation with the Joule heating as a source and a current conservation equation with temperature-dependent electrical conductivity.This problem has important applications in industry.In this paper,A new finite difference scheme is proposed on nonuniform rectangular partition for the thermistor problem.In the theoretical analyses,the second-order error estimates are obtained for electrical potential in discrete L2 and H1 norms,and for the temperature in L2 norm.In order to get these second-order error estimates,the Joule heating source is used in a changed equivalent form.     

The mathematical modeling of the electric field in the media with anisotropic objects

We present a numerical scheme for modeling the electric field in the media with tensor conductivity. This scheme is based on vector finite element method in frequency domain. The numerical computations of the electric field in the anisotropic medium are done. The conductivity of the anisotropic medium is positive defined dense tensor in general case. We consider the electric field from anisotropic layer, inclined anisotropic layer and some anisotropic objects in isotropic half-space.     

A mathematical model for impulse resistance welding

We present a mathematical model of impulse resistance welding. It accounts for electrical, thermal and mechanical effects, which are non‐linearly coupled by the balance laws, constitutive equations and boundary conditions. The electrical effects of the weld machine are incorporated by a discrete oscillator circuit which is coupled to the field equations by a boundary condition. We prove the existence of weak solutions for a slightly simplified model which however still covers most of its essential features, e.g. the quadratic Joule heat term and a quadratic term due to non‐elastic energy dissipation. We discuss the numerical implementation in a 2D setting, present some numerical results and conclude with some remarks on future research. Copyright © 2003 John Wiley & Sons, Ltd.     

Analysis of an optimal control problem for the tridomain model in cardiac electrophysiology

In the present paper, an optimal control problem constrained by the tridomain equations in electrocardiology is investigated. The state equations consisting in a coupled reaction–diffusion system modeling the propagation of the intracellular and extracellular electrical potentials, and ionic currents, are extended to further consider the effect of an external bathing medium. The existence and uniqueness of solution for the tridomain problem and the related control problem is assessed, and the primal and dual problems are discretized using a finite volume method which is proved to converge to the corresponding weak solution. In order to illustrate the control of the electrophysiological dynamics, we present some preliminary numerical experiments using an efficient implementation of the proposed scheme.     

求解带有非线性边界条件的涡流方程的A-φ解耦有限元格式

本文研究边界条件符合幂指数型非线性关系H × n = n × (|E × n|^α-1E × n)(0 < α ≤ 1)的涡流方程.使用A-φ耦合有限元格式数值求解这类问题具有较高精度,但计算开销大. A-φ解耦有限元计算格式能够在每个时间步上分别求解矢量A和标量φ,以此降低计算规模,提高计算效率.我们证明了解耦格式中解的存在唯一性,并且给出了它的误差估计.最后给出的数值实验证明了本文所提供的解耦算法是稳定和有效的.     

Numerical method for solving a two-dimensional electrical impedance tomography problem in the case of measurements on part of the outer boundary

The two-dimensional electrical impedance tomography problem is considered in the case of a piecewise constant electrical conductivity. The task is to determine the unknown boundary separating the regions with different conductivity values, which are known. Input information is the electric field measured on a portion of the outer boundary of the domain. A numerical method for solving the problem is proposed, and numerical results are presented.     

On the compact difference scheme for the two-dimensional coupled nonlinear Schrödinger equations

A high-order finite difference method for the two-dimensional coupled nonlinear Schrödinger equations is considered. The proposed scheme is proved to preserve the total mass and energy in a discrete sense and the solvability of the scheme is shown by using a fixed point theorem. By converting the scheme in the point-wise form into a matrix–vector form, we use the standard energy method to establish the optimal error estimate of the proposed scheme in the discrete L²-norm. The convergence order is proved to be of a fourth-order in space and a second-order in time, respectively. Finally, some numerical examples are given in order to confirm our theoretical results for the numerical method. The numerical results are compared with exact solutions and other existing method. The comparison between our numerical results and those of Sun and Wangreveals that our method improves the accuracy of space and time directions.     

Convergence of the mixed finite element method for Maxwell's equations with nonlinear conductivity

In this paper, we study a numerical scheme to solve coupled Maxwell's equations with a nonlinear conductivity. This model plays an important role in the study of type‐II superconductors. The approximation scheme is based on backward Euler discretization in time and mixed conforming finite elements in space. We will prove convergence of this scheme to the unique weak solution of the problem and develop the corresponding error estimates. As a next step, we study the stability of the scheme in the quasi‐static limit ? → 0 and present the corresponding convergence rate. Finally, we support the theory by several numerical experiments. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical simulation of thermal problems coupled with magnetohydrodynamic effects in aluminium cell

A system of partial differential equations describing the thermal behavior of aluminium cell coupled with magnetohydrodynamic effects is numerically solved. The thermal model is considered as a two-phases Stefan problem which consists of a non-linear convection–diffusion heat equation with Joule effect as a source. The magnetohydrodynamic fields are governed by Navier–Stokes and by static Maxwell equations. A pseudo-evolutionary scheme (Chernoff) is used to obtain the stationary solution giving the temperature and the frozen layer profile for the simulation of the ledges in the cell. A numerical approximation using a finite element method is formulated to obtain the fluid velocity, electrical potential, magnetic induction and temperature. An iterative algorithm and 3-D numerical results are presented.     

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.     

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.     

Weak solutions of the forward problem in EEG for different conductivity values

The potential distribution on the scalp produced by current sources in the brain can be measured by an EEG recorder. The relationship between these sources and the scalp potential distribution may be described by a well-known mathematical model where some simplifications are usually introduced. The head is modeled as a multicompartment nested set and the conductivity of the different tissues is approximated by a positive piecewise constant function. This simplified model is used to solve the forward problem (FP), i.e., to calculate the scalp potential for a current source configuration. In this work, we prove that the weak solutions of the FP are continuous with respect to the conductivity values, that is, the difference between the scalp potentials is small if the conductivity values are closed enough. We present numerical examples that illustrates this property.     

On explicit representation and approximations of Dirichlet-to-Neumann semigroup

In his book (Functional Analysis, Wiley, New York, 2002), P. Lax constructs an explicit representation of the Dirichlet-to-Neumann semigroup, when the matrix of electrical conductivity is the identity matrix and the domain of the problem in question is the unit ball in ? ⁿ . We investigate some representations of Dirichlet-to-Neumann semigroup for a bounded domain. We show that such a nice explicit representation as in Lax book, is not possible for any domain except Euclidean balls. It is interesting that the treatment in dimension 2 is completely different than other dimensions. Finally, we present a natural and probably the simplest numerical scheme to calculate this semigroup in full generality by using Chernoff’s theorem.     

A numerical method for mass conservative coupling between fluid flow and solute transport

We present a new coupled discretization approach for species transport in an incompressible fluid. The Navier-Stokes equations for the flow are discretized by the divergence-free Scott-Vogelius element on barycentrically refined meshes guaranteeing LBB stability. The convection-diffusion equation for species transport is discretized by the Voronoi finite volume method. In accordance to the continuous setting, due to the exact integration of the normal component of the flow through the Voronoi surfaces, the species concentration fulfills discrete global and local maximum principles. Besides of the numerical scheme itself, we present important aspects of its implementation. Further, for the case of homogeneous Dirichlet boundary conditions, we give a convergence proof for the coupled scheme. We report results of the application of the scheme to the interpretation of limiting current measurements in an electrochemical flow cell with cylindrical shape.     

Unsteady convective boundary layer flow of a viscous fluid at a vertical surface with variable fluid properties

In this paper we present numerical solutions to the unsteady convective boundary layer flow of a viscous fluid at a vertical stretching surface with variable transport properties and thermal radiation. Both assisting and opposing buoyant flow situations are considered. Using a similarity transformation, the governing time-dependent partial differential equations are first transformed into coupled, non-linear ordinary differential equations with variable coefficients. Numerical solutions to these equations subject to appropriate boundary conditions are obtained by a second order finite difference scheme known as the Keller-Box method. The numerical results thus obtained are analyzed for the effects of the pertinent parameters namely, the unsteady parameter, the free convection parameter, the suction/injection parameter, the Prandtl number, the thermal conductivity parameter and the thermal radiation parameter on the flow and heat transfer characteristics. It is worth mentioning that the momentum and thermal boundary layer thicknesses decrease with an increase in the unsteady parameter.     

A finite element approach for analysis and computational modelling of coupled reaction diffusion models

In this article, we establish the existence and uniqueness of solutions to the coupled reaction–diffusion models using Banach fixed point theorem. The Galerkin finite element method is used for the approximation of solutions, and an a priori error estimate is derived for such approximations. A scheme is proposed by combining the Crank–Nicolson and the predictor–corrector methods for the time discretization. Some numerical examples are considered to illustrate the accuracy and efficiency of the proposed scheme. It is found that the scheme is second‐order convergent. In addition, nonuniform grids are used in some cases to enhance the accuracy of the scheme.     

Global impulsive exponential synchronization of time-delayed coupled chaotic systems

In this paper, the impulsive exponential synchronization problem for time-delayed coupled chaotic systems is investigated. By establishing an impulsive differential delay inequality and using the property of P-cone, some simple conditions of impulsive exponential synchronization of two coupled chaotic systems are derived. To illustrate the effectiveness of the new scheme, some numerical examples are given.     

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.     

A Decomposition Method for a Unilateral Contact Problem with Tresca Friction Arising in Electro-elastostatics

This article is concerned with the numerical modeling of unilateral contact problems in an electro-elastic material with Tresca friction law and electrical conductivity condition. First, we prove the existence and uniqueness of the weak solution of the model. Rather than deriving a solution method for the full coupled problem, we present and study a successive iterative (decomposition) method. The idea is to solve successively a displacement subproblem and an electric potential subproblem in block Gauss-Seidel fashion. The displacement subproblem leads to a constraint non-differentiable (convex) minimization problem for which we propose an augmented Lagrangian algorithm. The electric potential unknown is computed explicitly using the Riesz's representation theorem. The convergence of the iterative decomposition method is proved. Some numerical experiments are carried out to illustrate the performances of the proposed algorithm.