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...

The goal of this article is to analyze the observability properties for a space semi-discrete approximation scheme derived from a mixed finite element method of the 1d wave equation on nonuniform meshes. More precisely, we prove that observability properties hold uniformly with respect to the mesh-size under some assumptions, which, roughly, measures the lack of uniformity of the meshes, thus extending the work [Castro and Micu, Numer. Math.102 (2006) 413–462] to nonuniform meshes. Our results...

Experimental evidence collected over the years shows that blood exhibits non-Newtonian characteristics such as shear-thinning, viscoelasticity, yield stress and thixotropic behaviour. Under certain conditions these characteristics become relevant and must be taken into consideration when modelling blood flow. In this work we deal with incompressible generalized Newtonian fluids, that account for the non-constant viscosity of blood, and present a new numerical method to handle fluid-rigid body interaction...

The paper is concerned with the study of a parabolic initial-boundary value problem with nonlinear Newton boundary condition considered in a two-dimensional domain. The goal is to prove the existence and uniqueness of a weak solution to the problem in the case when the nonlinearity in the Newton boundary condition does not satisfy any monotonicity condition and to analyze the finite element approximation.

The subject of the paper is the derivation and analysis of evolution Galerkin schemes for the two dimensional Maxwell and linearized Euler equations. The aim is to construct a method which takes into account better the infinitely many directions of propagation of waves. To do this the initial function is evolved using the characteristic cone and then projected onto a finite element space. We derive the divergence-free property and estimate the dispersion relation as well. We present some numerical...

In this paper the fluid-structure interaction problem is studied on a simplified model of the human vocal fold. The problem is mathematically described and the arbitrary Lagrangian-Eulerian method is applied in order to treat the time dependent computational domain. The viscous incompressible fluid flow and linear elasticity models are considered. The fluid flow and the motion of elastic body is approximated with the aid of finite element method. An attention is paid to the applied stabilization...

We consider a model for phase separation of a multi-component alloy with non-smooth free energy and a degenerate mobility matrix. In addition to showing well-posedness and stability bounds for our approximation, we prove convergence in one space dimension. Furthermore an iterative scheme for solving the resulting nonlinear discrete system is analysed. We discuss also how our approximation has to be modified in order to be applicable to a logarithmic free energy. Finally numerical experiments with...

We consider a model for phase separation of a multi-component alloy with non-smooth free energy and a degenerate mobility matrix. In addition to showing well-posedness and stability bounds for our approximation, we prove convergence in one space dimension. Furthermore an iterative scheme for solving the resulting nonlinear discrete system is analysed. We discuss also how our approximation has to be modified in order to be applicable to a logarithmic free energy. Finally numerical experiments...

The paper deals with numerical simulation of a compressible flow in time-dependent 2D domains with a special interest in medical applications to airflow in the human vocal tract. The mathematical model of this process is described by the compressible Navier-Stokes equations. For the treatment of the time-dependent domain, the arbitrary Lagrangian-Eulerian (ALE) method is used. The discontinuous Galerkin finite element method (DGFEM) is used for the space semidiscretization of the governing equations...

In this paper we consider the existence and asymptotic behavior of solutions of the following problem: $$u_{tt}{(t,x)-(\alpha +\beta \parallel \nabla u(t,x)\parallel }_2^2+{\beta \parallel \nabla v(t,x)\parallel }_2^2)\Delta u(t,x)+\delta |u_t{(t,x)|}^{p-1}{u}_t(t,x)={\mu \left|u(t,x)\right|}^{q-1}u(t,x),x\in \Omega ,t\ge 0,v_{tt}{(t,x)-(\alpha +\beta \parallel \nabla u(t,x)\parallel }_2^2+{\beta \parallel \nabla v(t,x)\parallel }_2^2)\Delta v(t,x)+\delta |v_t{(t,x)|}^{p-1}{v}_t(t,x)={\mu \left|v(t,x)\right|}^{q-1}v(t,x),x\in \Omega ,t\ge 0,u(0,x)=u_0\left(x\right),u_t(0,x)=u_1\left(x\right),x\in \Omega ,v(0,x)=v_0\left(x\right),v_t(0,x)=v_1\left(x\right),x\in \Omega ,u{{|}_{{}_{\partial \Omega }}=v|}_{{}_{\partial \Omega }}=0$$ where $q>1$, $p\ge 1$, $\delta >0$, $\alpha >0$, $\beta \ge 0$, $\mu \in ℝ$ and $\Delta$ is the Laplacian in ${ℝ}^N$.

A three-parameter family of Boussinesq type systems in two space dimensions is considered. These systems approximate the three-dimensional Euler equations, and consist of three nonlinear dispersive wave equations that describe two-way propagation of long surface waves of small amplitude in ideal fluids over a horizontal bottom. For a subset of these systems it is proved that their Cauchy problem is locally well-posed in suitable Sobolev classes. Further, a class of these systems is discretized...