Numerical Simulation of a Nonlinear Problem Arising in Heat Transfer and Magnetostatics

The unknown u can be temperature, magnetic potential, etc., | · | denotes the Euclidean norm and ∇ and ∇· denote, respectively, the usual gradient and divergence operators. The functions k, f and g are data, and we assume that k : Ω× [0,+∞) → R satisfies the following bound, for two positive constants k1 and k2: k1 ≤ k(x, s) + s ∂ ∂s k(x, s) ≤ k2, ∀(x, s) ∈ Ω× [0,+∞). Problem (1) allows to predict a wide range of physical phenomena, such as magnetostatics and heat transfer. Due to the nonlinearity, we need to resort to numerical methods. We apply an adaptive stabilized mixed finite element method and present some results from our numerical experiments.


Introduction
We consider the mathematical model from [1]. The model comprises a second order nonlinear elliptic partial differential equation over a bounded domain Ω in R 2 , provided with Dirichlet boundary conditions over the contour Γ:    −∇ · (k(·, |∇u|)∇u) = f in Ω , u = g on Γ. (1) The unknown u can be temperature, magnetic potential, etc., | · | denotes the Euclidean norm and ∇ and ∇· denote, respectively, the usual gradient and divergence operators. The functions k, f and g are data, and we assume that k : Ω × [0, +∞) → R satisfies the following bound, for two positive constants k 1 and k 2 : k 1 ≤ k(x, s) + s ∂ ∂s k(x, s) ≤ k 2 , ∀(x, s) ∈ Ω × [0, +∞). Problem (1) allows to predict a wide range of physical phenomena, such as magnetostatics and heat transfer. Due to the nonlinearity, we need to resort to numerical methods. We apply an adaptive stabilized mixed finite element method and present some results from our numerical experiments.

Adaptive Augmented Mixed Finite Element Method
We introduce two additional unknowns, t = ∇u and σ = k(·, |t|)t, and define the function space In order to approximate the triplet (t, σ, u), the domain is divided into triangular cells that conform a mesh (see Figure 1a). In each of those cells an approximation of the unknowns is calculated. The corresponding function space is X h := P 0 × RT 0 × P 1 , with P k being the space of polynomials of degree ≤ k and RT 0 is the lowest order Raviart Thomas space. We consider the following discrete problem: with κ i being positive constants. The nonlinear form A(·, ·) is strongly monotone, Lipschitz-continuous and bounded for appropriate values of the parameters κ i so that problem (2) has a unique solution. We consider an adaptive algorithm, which estimates the local error in each cell, and refines the mesh over those regions with the highest error. On each triangle, we introduce the following novel error indicator, which is reliable and locally efficient [2]: where E(T) denotes the set of edges of T and E(Γ) the set of edges in the boundary.

Results and Discussion
The adaptive algorithm was tested on a singular example over the unit square with a singularity in its upper right corner. The exact solution is u(x, y) = (2.1 − x − y) −1/3 . In Figure 1b, we can see that the algorithm refines the mesh near the singularity. Figure 1c shows a graph that compares, for both uniform and adaptive refinements, the exact error and the defined indicator vs. degrees of freedom (dof). With the adaptive algorithm the total error is lower, and a good agreement between error and indicator is seen. In fact, the efficiency index, which is the quotient between indicator and total error, tends to one as the mesh size decreases.
These simulations were carried out in FreeFEM++, following the mesh refinement algorithm proposed in [3], and the nonlinear system was solved with Newton's method with a tolerance of 10 −9 , needing four iterations.

Conflicts of Interest:
The authors declare no conflict of interest.