Numerical Solution of Some Diffusion Problems in 3-Layered 3D Domain

In this paper we consider averaging and finite difference methods for solving the 3-D boundary-value problem in a multilayered domain. We consider the metal concentration in the 3 layered peat blocks. Using experimental data the mathematical model for calculating the concentration of metal at different points in peat layers is developed. A specific feature of these problems is that it is necessary to solve the 3-D boundary-value problems for the partial differential equations （PDEs） of the elliptic type of second order with piece-wise diffusion coefficients in the three layer domain. We develop here a finite-difference method for solving a problem of the above type with the periodical boundary condition in x direction. This procedure allows reducing the 3-D problem to a system of 2-D problems by using a circulant matrix.

A COUPLED CONTINUOUS-DISCONTINUOUS FEM APPROACH FOR CONVECTION DIFFUSION EQUATIONS

In this article, we introduce a coupled approach of local discontinuous Calerkin and standard finite element method for solving convection diffusion problems. The whole domain is divided into two disjoint subdomains. The discontinuous Galerkin method is adopted in the subdomain where the solution varies rapidly, while the standard finite element method is used in the other subdomain due to its lower computational cost. The stability and a priori error estimate are established. We prove that the coupled method has O（ε1/2 ＋ h1/2）hk） convergence rate in an associated norm, where ε is the diffusion coefficient, h is the mesh size and k is the degree of polynomial. The numerical results verify our theoretical results. Moreover, 2k-order superconvergence of the numerical traces at the nodes, and the optimal convergence of the errors under L2 norm are observed numerically on the uniform mesh. The numerical results also indicate that the coupled method has the same convergence order and almost the same errors as the purely LDG method.

On the Robustness of the xy-Zebra-Gauss-Seidel Smoother on an Anisotropic Diffusion Problem

Studies of problems involving physical anisotropy are applied in sciences and engineering,for instance,when the thermal conductivity depends on the direction.In this study,the multigrid method was used in order to accelerate the convergence of the iterative methods used to solve this type of problem.The asymptotic convergence factor of the multigrid was determined empirically(computer aided)and also by employing local Fourier analysis(LFA).The mathematical model studied was the 2D anisotropic diffusion equation,in whichε>0 was the coefficient of a nisotropy.The equation was discretized by the Finite Difference Method(FDM)and Central Differencing Scheme(CDS).Correction Scheme(CS),pointwise Gauss-Seidel smoothers(Lexicographic and Red-Black ordering),and line Gauss-Seidel smoothers(Lexicographic and Zebra ordering)in x and y directions were used for building the multigrid.The best asymptotic convergence factor was obtained by the Gauss-Seidel method in the direction x for 0<ε<<1 and in the direction y forε>>1.In this sense,an xy-zebra-GS smoother was proposed,which proved to be efficient and robust for the different anisotropy coefficients.Moreover,the convergence factors calculated empirically and by LFA are in agreement.

ANALYSIS OF BOUNDARY LAYER SINGULARITYIN A NONLINEAR DIFFUSION PROBLEM

In this paper, the author analyzes the singularity of a boundary layer in a nonlinear diffusion problem. Results show when the limiting solution satisfies the boundary condition, there is no boundary singularity. Otherwise, the boundary layer exists, and its thickness is proportional to epsilon(1/2), here epsilon is a small positive real parameter.

Localization of Solutions of a Nonlinear Diffusion Problem

This paper concerns with properties of solutions of a nonlinear diffusion problem in non-divergence form. By constructing proper test functions, it is proved that solutions of the problem possess the property of localization.

An Approximate Riemann Solver for Advection-Diffusion Based on the Generalized Riemann Problem

We construct an approximate Riemann solver for scalar advection-diffusion equations with piecewise polynomial initial data.The objective is to handle advection and diffusion simultaneously to reduce the inherent numerical diffusion produced by the usual advection flux calculations.The approximate solution is based on the weak formulation of the Riemann problem and is solved within a space-time discontinuous Galerkin approach with two subregions.The novel generalized Riemann solver produces piecewise polynomial solutions of the Riemann problem.In conjunction with a recovery polynomial,the Riemann solver is then applied to define the numerical flux within a finite volume method.Numerical results for a piecewise linear and a piecewise parabolic approximation are shown.These results indicate a reduction in numerical dissipation compared with the conventional separated flux calculation of advection and diffusion.Also,it is shown that using the proposed solver only in the vicinity of discontinuities gives way to an accurate and efficient finite volume scheme.

Effects of environmentally friendly agricultural land protection programs: Evidence from the Lake Seyfe area of Turkey

fully supported by the Turkish Scientific Research Council （TUBITAK in Turkish, 110O747）

Computable extensions of generalized fractional kinetic equations in astrophysics

Fractional calculus and special functions have contributed a lot to mathematical physics and its various branches. The great use of mathematical physics in distinguished astrophysical problems has attracted astronomers and physicists to pay more attention to available mathematical tools that can be widely used in solving several problems of astrophysics/physics. In view of the great importance and usefulness of kinetic equations in certain astrophysical problems, the authors derive a generalized fractional kinetic equation involving the Lorenzo-Hartley function, a generalized function for fractional calculus. The fractional kinetic equation discussed here can be used to investigate a wide class of known （and possibly also new） fractional kinetic equations, hitherto scattered in the literature. A compact and easily computable solution is established in terms of the Lorenzo-Hartley function. Special cases, involving the generalized Mittag-Leffler function and the R-function, are considered. The obtained results imply the known results more precisely.

DiscreteMaximumPrinciple Based on Repair Technique for Finite Element Scheme of Anisotropic Diffusion Problems

In this paper,we construct a global repair technique for the finite element scheme of anisotropic diffusion equations to enforce the repaired solutions satisfying the discrete maximum principle.It is an extension of the existing local repair technique.Both of the repair techniques preserve the total energy and are easy to be implemented.The numerical experiments show that these repair techniques do not destroy the accuracy of the finite element scheme,and the computational cost of the global repair technique is cheaper than the local repair technique when the diffusion tensors are highly anisotropic.

New Splitting Methods for Convection-Dominated Diffusion Problems and Navier-Stokes Equations

We present a new splitting method for time-dependent convention-dominated diffusion problems.The original convention diffusion system is split into two sub-systems:a pure convection system and a diffusion system.At each time step,a convection problem and a diffusion problem are solved successively.A few important features of the scheme lie in the facts that the convection subproblem is solved explicitly and multistep techniques can be used to essentially enlarge the stability region so that the resulting scheme behaves like an unconditionally stable scheme;while the diffusion subproblem is always self-adjoint and coercive so that they can be solved efficiently using many existing optimal preconditioned iterative solvers.The scheme can be extended for solving the Navier-Stokes equations,where the nonlinearity is resolved by a linear explicit multistep scheme at the convection step,while only a generalized Stokes problem is needed to solve at the diffusion step and the major stiffness matrix stays invariant in the time marching process.Numerical simulations are presented to demonstrate the stability,convergence and performance of the single-step and multistep variants of the new scheme.

Local Stability for an Inverse Coefficient Problem of a Fractional Diffusion Equation

Time-fractional diffusion equations are of great interest and importance on describing the power law decay for diffusion in porous media. In this paper, to identify the diffusion rate, i.e., the heterogeneity of medium, the authors consider an inverse coefficient problem utilizing finite measurements and obtain a local HSlder type conditional stability based upon two Carleman estimates for the corresponding differential equations of integer orders.

ANALYSIS ON A NUMERICAL SCHEME WITH SECOND-ORDER TIME ACCURACY FOR NONLINEAR DIFFUSION EQUATIONS

A nonlinear fully implicit finite difference scheme with second-order time evolution for nonlinear diffusion problem is studied.The scheme is constructed with two-layer coupled discretization(TLCD)at each time step.It does not stir numerical oscillation,while permits large time step length,and produces more accurate numerical solutions than the other two well-known second-order time evolution nonlinear schemes,the Crank-Nicolson(CN)scheme and the backward difference formula second-order(BDF2)scheme.By developing a new reasoning technique,we overcome the difficulties caused by the coupled nonlinear discrete diffusion operators at different time layers,and prove rigorously the TLCD scheme is uniquely solvable,unconditionally stable,and has second-order convergence in both s-pace and time.Numerical tests verify the theoretical results,and illustrate its superiority over the CN and BDF2 schemes.

A MFE method combined with L1-approximation for a nonlinear time-fractional coupled diffusion system

In this paper,a nonlinear time-fractional coupled diffusion system is solved by using a mixed finite element(MFE)method in space combined with L1-approximation and implicit second-order backward difference scheme in time.The stability for nonlinear fully discrete finite element scheme is analyzed and a priori error estimates are derived.Finally,some numerical tests are shown to verify our theoretical analysis.

Asymptotic-Preserving Discrete Schemes for Non-Equilibrium Radiation Diffusion Problem in Spherical and Cylindrical Symmetrical Geometries

We study the asymptotic-preserving fully discrete schemes for nonequilibrium radiation diffusion problem in spherical and cylindrical symmetric geometry.The research is based on two-temperature models with Larsen’s flux-limited diffusion operators.Finite volume spatially discrete schemes are developed to circumvent the singularity at the origin and the polar axis and assure local conservation.Asymmetric second order accurate spatial approximation is utilized instead of the traditional first order one for boundary flux-limiters to consummate the schemes with higher order global consistency errors.The harmonic average approach in spherical geometry is analyzed,and its second order accuracy is demonstrated.By formal analysis,we prove these schemes and their corresponding fully discrete schemes with implicitly balanced and linearly implicit time evolutions have first order asymptoticpreserving properties.By designing associated manufactured solutions and reference solutions,we verify the desired performance of the fully discrete schemes with numerical tests,which illustrates quantitatively they are first order asymptotic-preserving and basically second order accurate,hence competent for simulations of both equilibrium and non-equilibrium radiation diffusion problems.

BEHAVIOR OF SOLUTIONS TO A DEGENERATE DIFFUSION PROBLEM

The goal of this work is to study the behavior of solutions to a degenerate diffusion problem. The main interests center on the extinction properties and evolution of the supports of solutions for the equations considered.

AN EXPANDED CHARACTERISTIC-MIXED FINITE ELEMENT METHOD FOR A CONVECTION-DOMINATED TRANSPORT PROBLEM

In this paper, we propose an Expanded Characteristic-mixed Finite Element Method for approximating the solution to a convection dominated transport problem. The method is a combination of characteristic approximation to handle the convection part in time and an expanded mixed finite element spatial approximation to deal with the diffusion part. The scheme is stable since fluid is transported along the approximate characteristics on the discrete level. At the same time it expands the standard mixed finite element method in the sense that three variables are explicitly treated： the scalar unknown, its gradient, and its flux. Our analysis shows the method approximates the scalar unknown, its gradient, and its flux optimally and simultaneously. We also show this scheme has much smaller time-truncation errors than those of standard methods. A numerical example is presented to show that the scheme is of high performance.

Fast Numerical Simulation of Two-Phase Transport Model in the Cathode of a Polymer Electrolyte Fuel Cell

In this paper,we apply streamline-diffusion and Galerkin-least-squares fi-nite element methods for 2D steady-state two-phase model in the cathode of polymer electrolyte fuel cell(PEFC)that contains a gas channel and a gas diffusion layer(GDL).This two-phase PEFC model is typically modeled by a modified Navier-Stokes equation for the mass and momentum,with Darcy’s drag as an additional source term in momentum for flows through GDL,and a discontinuous and degenerate convectiondiffusion equation for water concentration.Based on the mixed finite element method for the modified Navier-Stokes equation and standard finite element method for water equation,we design streamline-diffusion and Galerkin-least-squares to overcome the dominant convection arising from the gas channel.Meanwhile,we employ Kirchhoff transformation to deal with the discontinuous and degenerate diffusivity in water concentration.Numerical experiments demonstrate that our finite element methods,together with these numerical techniques,are able to get accurate physical solutions with fast convergence.

Near-Optimal Controls of Differential Systems with Switching and Random Jumps Subject to Fast Switching and Wideband Noise Perturbation

This work develops near-optimal controls for systems given by differential equations with wideband noise and random switching.The random switching is modeled by a continuous-time,time-inhomogeneous Markov chain.Under broad conditions,it is shown that there is an associated limit problem,which is a switching jump diffusion.Using near-optimal controls of the limit system,we then build controls for the original systems.It is shown that such constructed controls are nearly optimal.