A Hybrid Differential Transforms and Finite Difference Method to Numerical Solution of Convection–Diffusion Equation

In this work, we discuss a hybrid-based method on differential transforms and a finite difference method to numerical solution of convection–diffusion equation with Dirichlet’s type boundary conditions. The developed method is tested on various problems and the numerical results are reported in tabular and figure form. This method can be easily extended to handle non-linear convection–diffusion partial differential equations.

Implicit-explicit multistep formulations for finite element discretisations using continuous interior penalty

We consider a finite element method with symmetric stabilisation for the discretisation of the transient convection--diffusion equation. For the time-discretisation we consider either the second order backwards differentiation formula or the Crank-Nicolson method. Both the convection term and the associated stabilisation are treated explicitly using an extrapolated approximate solution. We prove stability of the method and the $\tau^2 + h^{p+{\frac12}}$ error estimates for the $L^2$-norm under either the standard hyperbolic CFL condition, when piecewise affine ($p=1$) approximation is used, or in the case of finite element approximation of order $p \ge 1$, a stronger, so-called $4/3$-CFL, i.e. $\tau \leq C h^{4/3}$. The theory is illustrated with some numerical examples.

Shallow water model for pollutant distribution in the Ypacarai Lake

In this paper, we analyze the distribution of a non-reactive contaminant in Ypacarai Lake. We propose a shallow-water model that considers wind-induced currents, inflow and outflow conditions in the tributaries, and bottom effects due to the lakebed. The hydrodynamic is based on the depth-averaged Navier-Stokes equations considering wind stresses as force terms which are functions of the wind velocity. Bed (bottom) stress is based on Manning's equation, the lakebed characteristics, and wind velocities. The contaminant transportation is modeled by a 2D convection-diffusion equation taking into consideration water level. Comparisons between the simulation of the model, analytical solutions, and laboratory results confirm that the model captures the complex dynamic phenomenology of the lake. In the simulations, one can see the regions with the highest risk of accumulation of contaminants. It is observed the effect of each term and how it can be used them to mitigate the impact of the pollutants.