No Arabic abstract
In Lipschitz two and three dimensional domains, we study the existence for the so--called Boussinesq model of thermally driven convection under singular forcing. By singular we mean that the heat source is allowed to belong to $H^{-1}(varpi,Omega)$, where $varpi$ is a weight in the Muckenhoupt class $A_2$ that is regular near the boundary. We propose a finite element scheme and, under the assumption that the domain is convex and $varpi^{-1} in A_1$, show its convergence. In the case that the thermal diffusion and viscosity are constants, we propose an a posteriori error estimator and show its reliability and local efficiency.
In two dimensions, we show existence of solutions to the stationary Navier Stokes equations on weighted spaces $mathbf{H}^1_0(omega,Omega) times L^2(omega,Omega)$, where the weight belongs to the Muckenhoupt class $A_2$. We show how this theory can be applied to obtain a priori error estimates for approximations of the solution to the Navier Stokes problem with singular sources.
We consider the periodic initial-value problem for the Serre equations of water-wave theory and its semidiscrete approximation in the space of smooth periodic polynomial splines. We prove that the semidiscrete problem is well posed, locally in time, and satisfies a discrete positivity property for the water depth.
A direct reconstruction algorithm based on Calderons linearization method for the reconstruction of isotropic conductivities is proposed for anisotropic conductivities in two-dimensions. To overcome the non-uniqueness of the anisotropic inverse conductivity problem, the entries of the unperturbed anisotropic tensors are assumed known emph{a priori}, and it remains to reconstruct the multiplicative scalar field. The quasi-conformal map in the plane facilitates the Calderon-based approach for anisotropic conductivities. The method is demonstrated on discontinuous radially symmetric conductivities of high and low contrast.
A Morley-Wang-Xu (MWX) element method with a simply modified right hand side is proposed for a fourth order elliptic singular perturbation problem, in which the discrete bilinear form is standard as usual nonconforming finite element methods. The sharp error analysis is given for this MWX element method. And the Nitsches technique is applied to the MXW element method to achieve the optimal convergence rate in the case of the boundary layers. An important feature of the MWX element method is solver-friendly. Based on a discrete Stokes complex in two dimensions, the MWX element method is decoupled into one Lagrange element method of Poisson equation, two Morley element methods of Poisson equation and one nonconforming $P_1$-$P_0$ element method of Brinkman problem, which implies efficient and robust solvers for the MWX element method. Some numerical examples are provided to verify the theoretical results.
In this paper we consider the numerical solution of Boussinesq-Peregrine type systems by the application of the Galerkin finite element method. The structure of the Boussinesq systems is explained and certain alternative nonlinear and dispersive terms are compared. A detailed study of the convergence properties of the standard Galerkin method, using various finite element spaces on unstructured triangular grids, is presented. Along with the study of the Peregrine system, a new Boussinesq system of BBM-BBM type is derived. The new system has the same structure in its momentum equation but differs slightly in the mass conservation equation compared to the Peregrine system. Further, the finite element method applied to the new system has better convergence properties, when used for its numerical approximation. Due to the lack of analytical formulas for solitary wave solutions for the systems under consideration, a Galerkin finite element method combined with the Petviashvili iteration is proposed for the numerical generation of accurate approximations of line solitary waves. Various numerical experiments related to the propagation of solitary and periodic waves over variable bottom topography and their interaction with the boundaries of the domains are presented. We conclude that both systems have similar accuracy when approximate long waves of small amplitude while the Galerkin finite element method is more effective when applied to BBM-BBM type systems.