General treatment of essential boundary conditions in reduced order models for non-linear problems

Inhomogeneous essential boundary conditions must be carefully treated in the formulation of Reduced Order Models (ROMs) for non-linear problems. In order to investigate this issue, two methods are analysed: one in which the boundary conditions are imposed in an strong way, and a second one in which a weak imposition of boundary conditions is made. The ideas presented in this work apply to the big realm of a posteriori ROMs. Nevertheless, an a posteriori hyper-reduction method is specifically considered in order to deal with the cost associated to the non-linearity of the problems. Applications to nonlinear transient heat conduction problems with temperature dependent thermophysical properties and time dependent essential boundary conditions are studied. However, the strategies introduced in this work are of general application.Fil: Cosimo, Alejandro. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; ArgentinaFil: Cardona, Alberto. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; ArgentinaFil: Idelsohn, Sergio Rodolfo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; Argentina. Institució Catalana de Recerca i Estudis Avancats; España. International Center for Numerical Methods in Engineering; Españ

Spectral methods in general relativistic astrophysics

We present spectral methods developed in our group to solve three-dimensional partial differential equations. The emphasis is put on equations arising from astrophysical problems in the framework of general relativity.Comment: 51 pages, elsart (Elsevier Preprint), 19 PostScript figures, submitted to Journal of Computational & Applied Mathematic

Adaptive computation of gravitational waves from black hole interactions

We construct a class of linear partial differential equations describing general perturbations of non-rotating black holes in 3D Cartesian coordinates. In contrast to the usual approach, a single equation treats all radiative $\ell -m$ modes simultaneously, allowing the study of wave perturbations of black holes with arbitrary 3D structure, as would be present when studying the full set of nonlinear Einstein equations describing a perturbed black hole. This class of equations forms an excellent testbed to explore the computational issues of simulating black spacetimes using a three dimensional adaptive mesh refinement code. Using this code, we present results from the first fully resolved 3D solution of the equations describing perturbed black holes. We discuss both fixed and adaptive mesh refinement, refinement criteria, and the computational savings provided by adaptive techniques in 3D for such model problems of distorted black holes.Comment: 16 Pages, RevTeX, 13 figure

The Cauchy Problem for the Einstein Equations

Various aspects of the Cauchy problem for the Einstein equations are surveyed, with the emphasis on local solutions of the evolution equations. Particular attention is payed to giving a clear explanation of conceptual issues which arise in this context. The question of producing reduced systems of equations which are hyperbolic is examined in detail and some new results on that subject are presented. Relevant background from the theory of partial differential equations is also explained at some lengthComment: 98 page

Forward and inverse problems in fundamental and applied magnetohydrodynamics

This Minireview summarizes the recent efforts to solve forward and inverse problems as they occur in different branches of fundamental and applied magnetohydrodynamics. As for the forward problem, the main focus is on the numerical treatment of induction processes, including self-excitation of magnetic fields in non-spherical domains and/or under the influence of non-homogeneous material parameters. As an important application of the developed numerical schemes, the functioning of the von-K\'{a}rm\'{a}n-sodium (VKS) dynamo experiment is shown to depend crucially on the presence of soft-iron impellers. As for the inverse problem, the main focus is on the mathematical background and some first practical applications of the Contactless Inductive Flow Tomography (CIFT), in which flow induced magnetic field perturbations are utilized for the reconstruction of the velocity field. The promises of CIFT for flow field monitoring in the continuous casting of steel are substantiated by results obtained at a test rig with a low melting liquid metal. While CIFT is presently restricted to flows with low magnetic Reynolds numbers, some selected problems of non-linear inverse dynamo theory, with possible application to geo- and astrophysics, are also discussed.Comment: 15 pages, 4 figures, accepted for publication in European Physical Journal Special Topic

Asymptotics of conduction velocity restitution in models of electrical excitation in the heart

We extend a non-Tikhonov asymptotic embedding, proposed earlier, for calculation of conduction velocity restitution curves in ionic models of cardiac excitability. Conduction velocity restitution is the simplest non-trivial spatially extended problem in excitable media, and in the case of cardiac tissue it is an important tool for prediction of cardiac arrhythmias and fibrillation. An idealized conduction velocity restitution curve requires solving a non-linear eigenvalue problem with periodic boundary conditions, which in the cardiac case is very stiff and calls for the use of asymptotic methods. We compare asymptotics of restitution curves in four examples, two generic excitable media models, and two ionic cardiac models. The generic models include the classical FitzHugh–Nagumo model and its variation by Barkley. They are treated with standard singular perturbation techniques. The ionic models include a simplified “caricature” of Noble (J. Physiol. Lond. 160:317–352, 1962) model and Beeler and Reuter (J. Physiol. Lond. 268:177–210, 1977) model, which lead to non-Tikhonov problems where known asymptotic results do not apply. The Caricature Noble model is considered with particular care to demonstrate the well-posedness of the corresponding boundary-value problem. The developed method for calculation of conduction velocity restitution is then applied to the Beeler–Reuter model. We discuss new mathematical features appearing in cardiac ionic models and possible applications of the developed method