2,786 research outputs found
Compressive Space-Time Galerkin Discretizations of Parabolic Partial Differential Equations
We study linear parabolic initial-value problems in a space-time variational
formulation based on fractional calculus. This formulation uses "time
derivatives of order one half" on the bi-infinite time axis. We show that for
linear, parabolic initial-boundary value problems on , the
corresponding bilinear form admits an inf-sup condition with sparse tensor
product trial and test function spaces. We deduce optimality of compressive,
space-time Galerkin discretizations, where stability of Galerkin approximations
is implied by the well-posedness of the parabolic operator equation. The
variational setting adopted here admits more general Riesz bases than previous
work; in particular, no stability in negative order Sobolev spaces on the
spatial or temporal domains is required of the Riesz bases accommodated by the
present formulation. The trial and test spaces are based on Sobolev spaces of
equal order with respect to the temporal variable. Sparse tensor products
of multi-level decompositions of the spatial and temporal spaces in Galerkin
discretizations lead to large, non-symmetric linear systems of equations. We
prove that their condition numbers are uniformly bounded with respect to the
discretization level. In terms of the total number of degrees of freedom, the
convergence orders equal, up to logarithmic terms, those of best -term
approximations of solutions of the corresponding elliptic problems.Comment: 26 page
Un método Wavelet-Galerkin para ecuaciones diferenciales parciales parabólicas
In this paper an Adaptive Wavelet-Galerkin method for the solution ofparabolic partial differential equations modeling physical problems withdifferent spatial and temporal scales is developed. A semi-implicit timedifference scheme is applied andB-spline multiresolution structure on theinterval is used. As in many cases these solutions are known to presentlocalized sharp gradients, local error estimators are designed and an ef-ficient adaptive strategy to choose the appropriate scale for each time isdeveloped. Finally, experiments were performed to illustrate the applica-bility and efficiency of the proposed method.En este trabajo se desarrolla un método Wavelet-Galerkin Adaptativopara la resolución de ecuaciones diferenciales parabólicas que modelanproblemas físicos, con diferentes escalas en el espacio y en el tiempo. Seutiliza un esquema semi-implícito en diferencias temporales y la estructuramultirresolución de las B-splines sobre intervalo.Como es sabido que enmuchos casos las soluciones presentan gradientes localmente altos, se handiseñado estimadores locales de error y una estrategia adaptativa eficientepara elegir la escala apropiada en cada tiempo. Finalmente, se realizaronexperimentos que ilustran la aplicabilidad y la eficiencia del método pro-puestoFil: Vampa, Victoria Cristina. Universidad Nacional de La Plata. Facultad de Ingeniería; ArgentinaFil: Martín, María Teresa. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. Universidad Nacional de La Plata. Facultad de Ingeniería; Argentin
Analytic Regularity and GPC Approximation for Control Problems Constrained by Linear Parametric Elliptic and Parabolic PDEs
This paper deals with linear-quadratic optimal control problems constrained by a parametric or stochastic elliptic or parabolic PDE. We address the (difficult) case that the state equation depends on a countable number of parameters i.e., on with , and that the PDE operator may depend non-affinely on the parameters. We consider tracking-type functionals and distributed as well as boundary controls. Building on recent results in [CDS1, CDS2], we show that the state and the control are analytic as functions depending on these parameters . We
establish sparsity of generalized polynomial chaos (gpc) expansions of both, state and control, in terms of the stochastic coordinate sequence of the random inputs, and prove convergence rates of best -term truncations of these expansions. Such truncations are the key for subsequent computations since they do {\em not} assume that the stochastic input data has a finite expansion. In the follow-up paper [KS2], we explain two methods how such best -term truncations can practically be computed, by greedy-type algorithms
as in [SG, Gi1], or by multilevel Monte-Carlo methods as in
[KSS]. The sparsity result allows in conjunction with adaptive wavelet Galerkin schemes for sparse, adaptive tensor discretizations of control problems constrained by linear elliptic and parabolic PDEs developed in [DK, GK, K], see [KS2]
A multiresolution space-time adaptive scheme for the bidomain model in electrocardiology
This work deals with the numerical solution of the monodomain and bidomain
models of electrical activity of myocardial tissue. The bidomain model is a
system consisting of a possibly degenerate parabolic PDE coupled with an
elliptic PDE for the transmembrane and extracellular potentials, respectively.
This system of two scalar PDEs is supplemented by a time-dependent ODE modeling
the evolution of the so-called gating variable. In the simpler sub-case of the
monodomain model, the elliptic PDE reduces to an algebraic equation. Two simple
models for the membrane and ionic currents are considered, the
Mitchell-Schaeffer model and the simpler FitzHugh-Nagumo model. Since typical
solutions of the bidomain and monodomain models exhibit wavefronts with steep
gradients, we propose a finite volume scheme enriched by a fully adaptive
multiresolution method, whose basic purpose is to concentrate computational
effort on zones of strong variation of the solution. Time adaptivity is
achieved by two alternative devices, namely locally varying time stepping and a
Runge-Kutta-Fehlberg-type adaptive time integration. A series of numerical
examples demonstrates thatthese methods are efficient and sufficiently accurate
to simulate the electrical activity in myocardial tissue with affordable
effort. In addition, an optimalthreshold for discarding non-significant
information in the multiresolution representation of the solution is derived,
and the numerical efficiency and accuracy of the method is measured in terms of
CPU time speed-up, memory compression, and errors in different norms.Comment: 25 pages, 41 figure
Adaptive Galerkin approximation algorithms for partial differential equations in infinite dimensions
Space-time variational formulations of infinite-dimensional Fokker-Planck (FP) and Ornstein-Uhlenbeck (OU) equations for functions on a separable Hilbert space are developed. The well-posedness of these equations in the Hilbert space of functions on , which are square-integrable with respect to a Gaussian measure on , is proved. Specifically, for the infinite-dimensional FP equation, adaptive space-time Galerkin discretizations, based on a tensorized Riesz basis, built from biorthogonal piecewise polynomial wavelet bases in time and the Hermite polynomial chaos in the Wiener-Itô decomposition of , are introduced and are shown to converge quasioptimally with respect to the nonlinear, best -term approximation benchmark. As a consequence, the proposed adaptive Galerkin solution algorithms perform quasioptimally with respect to the best -term approximation in the finite-dimensional case, in particular. All constants in our error and complexity bounds are shown to be independent of the number of "active" coordinates identified by the proposed adaptive Galerkin approximation algorithms
Adaptive multiresolution computations applied to detonations
A space-time adaptive method is presented for the reactive Euler equations
describing chemically reacting gas flow where a two species model is used for
the chemistry. The governing equations are discretized with a finite volume
method and dynamic space adaptivity is introduced using multiresolution
analysis. A time splitting method of Strang is applied to be able to consider
stiff problems while keeping the method explicit. For time adaptivity an
improved Runge--Kutta--Fehlberg scheme is used. Applications deal with
detonation problems in one and two space dimensions. A comparison of the
adaptive scheme with reference computations on a regular grid allow to assess
the accuracy and the computational efficiency, in terms of CPU time and memory
requirements.Comment: Zeitschrift f\"ur Physicalische Chemie, accepte
- …