A sparse algebraic multigrid method is studied as a cheap and accurate way to compute approximations of Schur complements of matrices arising from the discretization of some symmetric and positive definite partial differential operators. The construction of such a multigrid is discussed and numerical experiments are used to verify the properties of the method.

Assuming an incident wave to be a field source, we calculate the field potential in a neighborhood of an inhomogeneous body. This problem which has been formulated in ${𝐑}^3$can be reduced to a bounded domain. Namely, a boundary condition for the potential is formulated on a sphere. Then the potential satisfies a well posed boundary value problem in a ball containing the body. A numerical approximation is suggested and its convergence is analyzed.

This contribution gives an overview of current research in applying object oriented programming to scientific computing at the computational mechanics laboratory (LABMEC) at the school of civil engineering – UNICAMP. The main goal of applying object oriented programming to scientific computing is to implement increasingly complex algorithms in a structured manner and to hide the complexity behind a simple user interface. The following areas are current topics of research and documented within the...

This contribution gives an overview of current research in applying object oriented programming to scientific computing at the computational mechanics laboratory (LABMEC) at the school of civil engineering – UNICAMP. The main goal of applying object oriented programming to scientific computing is to implement increasingly complex algorithms in a structured manner and to hide the complexity behind a simple user interface. The following areas are current topics of research and documented within the...

A priori estimates and strong solvability results in Sobolev space $W^{2,p}\left(\Omega \right)$, $1<p<\infty$ are proved for the regular oblique derivative problem $$\left\{\begin{array}{c}\sum _{i,j=1}^n{a}^{ij}\left(x\right)\frac{{\partial }^{2}u}{\partial x_i\partial x_j}=f\left(x\right)\text{a.e.}\Omega \hfill \\ \frac{\partial u}{\partial \ell }+\sigma \left(x\right)u=\varphi \left(x\right)\text{on}\partial \Omega \hfill \end{array}\right.$$ when the principal coefficients $a^{ij}$ are $VMO\cap L^{\infty }$ functions.

The aim of this paper is to extend the results of Calderón [1] and Kenig-Pipher [12] on solutions to the oblique derivative problem to the case where the data is assumed to be BMO or Hölder continuous.

In this paper, we make some observations on the work of Di Fazio concerning $W^{1,p}$ estimates, $1<p<\mathrm{\infty }$, for solutions of elliptic equations $\text{div}A\nabla u=\text{div}f$ , on a domain $\mathrm{\Omega }$ with Dirichlet data $0$ whenever $A\in VMO\left(\mathrm{\Omega }\right)$ and $f\in L^p\left(\mathrm{\Omega }\right)$. We weaken the assumptions allowing real and complex non-symmetric operators and $C^1$ boundary. We also consider the corresponding inhomogeneous Neumann problem for which we prove the similar result. The main tool is an appropriate representation for the Green (and Neumann) function on the upper half space. We propose...

This work is devoted to the study of a two-dimensional vector Poisson equation with the normal component of the unknown and the value of the divergence of the unknown prescribed simultaneously on the entire boundary. These two scalar boundary conditions appear prima facie alternative in a standard variational framework. An original variational formulation of this boundary value problem is proposed here. Furthermore, an uncoupled solution algorithm is introduced together with its finite element...

A Bernoulli free boundary problem with geometrical constraints is studied. The domain Ω is constrained to lie in the half space determined by x1 ≥ 0 and its boundary to contain a segment of the hyperplane  {x1 = 0}  where non-homogeneous Dirichlet conditions are imposed. We are then looking for the solution of a partial differential equation satisfying a Dirichlet and a Neumann boundary condition simultaneously on the free boundary. The existence and uniqueness of a solution have already been addressed...

A Bernoulli free boundary problem with geometrical constraints is studied. The domain Ω is constrained to lie in the half space determined by x1 ≥ 0 and its boundary to contain a segment of the hyperplane  {x1 = 0}  where non-homogeneous Dirichlet conditions are imposed. We are then looking for the solution of a partial differential equation satisfying a Dirichlet and a Neumann boundary condition simultaneously on the free boundary. The existence and uniqueness of a solution have already been addressed...