Mathematical modelling of the cardiovascular system

In this paper we will address the problem of developing mathematical models for the numerical simulation of the human circulatory system. In particular, we will focus our attention on the problem of haemodynamics in large human arteries

Domain decomposition methods for systems of conservation laws: Spectral collocation approximations

Hyperbolic systems of conversation laws are considered which are discretized in space by spectral collocation methods and advanced in time by finite difference schemes. At any time-level a domain deposition method based on an iteration by subdomain procedure was introduced yielding at each step a sequence of independent subproblems (one for each subdomain) that can be solved simultaneously. The method is set for a general nonlinear problem in several space variables. The convergence analysis, however, is carried out only for a linear one-dimensional system with continuous solutions. A precise form of the error reduction factor at each iteration is derived. Although the method is applied here to the case of spectral collocation approximation only, the idea is fairly general and can be used in a different context as well. For instance, its application to space discretization by finite differences is straight forward

Domain decomposition preconditioners for the spectral collocation method

Several block iteration preconditioners are proposed and analyzed for the solution of elliptic problems by spectral collocation methods in a region partitioned into several rectangles. It is shown that convergence is achieved with a rate which does not depend on the polynomial degree of the spectral solution. The iterative methods here presented can be effectively implemented on multiprocessor systems due to their high degree of parallelism

On the coupling of hyperbolic and parabolic systems: Analytical and numerical approach

The coupling of hyperbolic and parabolic systems is discussed in a domain Omega divided into two distinct subdomains omega(+) and omega(-). The main concern is to find the proper interface conditions to be fulfilled at the surface separating the two domains. Next, they are used in the numerical approximation of the problem. The justification of the interface conditions is based on a singular perturbation analysis, i.e., the hyperbolic system is rendered parabolic by adding a small artifical viscosity. As this goes to zero, the coupled parabolic-parabolic problem degenerates into the original one, yielding some conditions at the interface. These are taken as interface conditions for the hyperbolic-parabolic problem. Actually, two alternative sets of interface conditions are discussed according to whether the regularization procedure is variational or nonvariational. It is shown how these conditions can be used in the frame of a numerical approximation to the given problem. Furthermore, a method of resolution is discussed which alternates the resolution of the hyperbolic problem within omega(-) and of the parabolic one within omega(+). The spectral collocation method is proposed, as an example of space discretization (different methods could be used as well); both explicit and implicit time-advancing schemes are considered. The present study is a preliminary step toward the analysis of the coupling between Euler and Navier-Stokes equations for compressible flows

Preconditioned Minimal Residual Methods for Chebyshev Spectral Caluclations

The problem of preconditioning the pseudospectral Chebyshev approximation of an elliptic operator is considered. The numerical sensitiveness to variations of the coefficients of the operator are investigated for two classes of preconditioning matrices: one arising from finite differences, the other from finite elements. The preconditioned system is solved by a conjugate gradient type method, and by a DuFort-Frankel method with dynamical parameters. The methods are compared on some test problems with the Richardson method and with the minimal residual Richardson method

On the boundary treatment in spectral methods for hyperbolic systems

Spectral methods were successfully applied to the simulation of slow transients in gas transportation networks. Implicit time advancing techniques are naturally suggested by the nature of the problem. The correct treatment of the boundary conditions are clarified in order to avoid any stability restriction originated by the boundaries. The Beam and Warming and the Lerat schemes are unconditionally linearly stable when used with a Chebyshev pseudospectral method. Engineering accuracy for a gas transportation problem is achieved at Courant numbers up to 100

Homoclinic snaking of localized states in doubly diffusive convection

Numerical continuation is used to investigate stationary spatially localized states in two-dimensional thermosolutal convection in a plane horizontal layer with no-slip boundary conditions at top and bottom. Convectons in the form of 1-pulse and 2-pulse states of both odd and even parity exhibit homoclinic snaking in a common Rayleigh number regime. In contrast to similar states in binary fluid convection, odd parity convectons do not pump concentration horizontally. Stable but time-dependent localized structures are present for Rayleigh numbers below the snaking region for stationary convectons. The computations are carried out for (inverse) Lewis number \tau = 1/15 and Prandtl numbers Pr = 1 and Pr >> 1

An Open-Domain Dialog Act Taxonomy

This document defines the taxonomy of dialog acts that are necessary to encode domain-independent dialog moves in the context of a task-oriented, open-domain dialog. Such taxonomy is formulated to satisfy two complementary requirements: on the one hand, domain independence, i.e. the power to cover all the range of possible interactions in any type of conversation (particularly conversation oriented to the performance of tasks). On the other hand, the ability to instantiate a concrete set of tasks as defined by a specific knowledge base (such as an ontology of domain concepts and actions) and within a particular language. For the modeling of dialog acts, inspiration is taken from several well-known dialog annotation schemes, such as DAMSL (Core & Allen, 1997), TRAINS (Traum, 1996) and VERBMOBIL (Alexandersson et al., 1997)

Modeling dimensionally-heterogeneous problems: analysis, approximation and applications

In the present work a general theoretical framework for coupled dimensionally-heterogeneous partial differential equations is developed. This is done by recasting the variational formulation in terms of coupling interface variables. In such a general setting we analyze existence and uniqueness of solutions for both the continuous problem and its finite dimensional approximation. This approach also allows the development of different iterative substructuring solution methodologies involving dimensionally-homogeneous subproblems. Numerical experiments are carried out to test our theoretical result

Computational Reduction for Parametrized PDEs: Strategies and Applications

In this paper we present a compact review on the mostly used techniques for computational reduction in numerical approximation of partial differential equations. We highlight the common features of these techniques and provide a detailed presentation of the reduced basis method, focusing on greedy algorithms for the construction of the reduced spaces. An alternative family of reduction techniques based on surrogate response surface models is briefly recalled too. Then, a simple example dealing with inviscid flows is presented, showing the reliability of the reduced basis method and a comparison between this technique and some surrogate model