14,071 research outputs found
Classical and all-floating FETI methods for the simulation of arterial tissues
High-resolution and anatomically realistic computer models of biological soft
tissues play a significant role in the understanding of the function of
cardiovascular components in health and disease. However, the computational
effort to handle fine grids to resolve the geometries as well as sophisticated
tissue models is very challenging. One possibility to derive a strongly
scalable parallel solution algorithm is to consider finite element tearing and
interconnecting (FETI) methods. In this study we propose and investigate the
application of FETI methods to simulate the elastic behavior of biological soft
tissues. As one particular example we choose the artery which is - as most
other biological tissues - characterized by anisotropic and nonlinear material
properties. We compare two specific approaches of FETI methods, classical and
all-floating, and investigate the numerical behavior of different
preconditioning techniques. In comparison to classical FETI, the all-floating
approach has not only advantages concerning the implementation but in many
cases also concerning the convergence of the global iterative solution method.
This behavior is illustrated with numerical examples. We present results of
linear elastic simulations to show convergence rates, as expected from the
theory, and results from the more sophisticated nonlinear case where we apply a
well-known anisotropic model to the realistic geometry of an artery. Although
the FETI methods have a great applicability on artery simulations we will also
discuss some limitations concerning the dependence on material parameters.Comment: 29 page
Long time dynamics and coherent states in nonlinear wave equations
We discuss recent progress in finding all coherent states supported by
nonlinear wave equations, their stability and the long time behavior of nearby
solutions.Comment: bases on the authors presentation at 2015 AMMCS-CAIMS Congress, to
appear in Fields Institute Communications: Advances in Applied Mathematics,
Modeling, and Computational Science 201
Nonlinear multigrid based on local spectral coarsening for heterogeneous diffusion problems
This work develops a nonlinear multigrid method for diffusion problems
discretized by cell-centered finite volume methods on general unstructured
grids. The multigrid hierarchy is constructed algebraically using aggregation
of degrees of freedom and spectral decomposition of reference linear operators
associated with the aggregates. For rapid convergence, it is important that the
resulting coarse spaces have good approximation properties. In our approach,
the approximation quality can be directly improved by including more spectral
degrees of freedom in the coarsening process. Further, by exploiting local
coarsening and a piecewise-constant approximation when evaluating the nonlinear
component, the coarse level problems are assembled and solved without ever
re-visiting the fine level, an essential element for multigrid algorithms to
achieve optimal scalability. Numerical examples comparing relative performance
of the proposed nonlinear multigrid solvers with standard single-level
approaches -- Picard's and Newton's methods -- are presented. Results show that
the proposed solver consistently outperforms the single-level methods, both in
efficiency and robustness
Characteristic Evolution and Matching
I review the development of numerical evolution codes for general relativity
based upon the characteristic initial value problem. Progress in characteristic
evolution is traced from the early stage of 1D feasibility studies to 2D
axisymmetric codes that accurately simulate the oscillations and gravitational
collapse of relativistic stars and to current 3D codes that provide pieces of a
binary black hole spacetime. Cauchy codes have now been successful at
simulating all aspects of the binary black hole problem inside an artificially
constructed outer boundary. A prime application of characteristic evolution is
to extend such simulations to null infinity where the waveform from the binary
inspiral and merger can be unambiguously computed. This has now been
accomplished by Cauchy-characteristic extraction, where data for the
characteristic evolution is supplied by Cauchy data on an extraction worldtube
inside the artificial outer boundary. The ultimate application of
characteristic evolution is to eliminate the role of this outer boundary by
constructing a global solution via Cauchy-characteristic matching. Progress in
this direction is discussed.Comment: New version to appear in Living Reviews 2012. arXiv admin note:
updated version of arXiv:gr-qc/050809
Towards absorbing outer boundaries in General Relativity
We construct exact solutions to the Bianchi equations on a flat spacetime
background. When the constraints are satisfied, these solutions represent in-
and outgoing linearized gravitational radiation. We then consider the Bianchi
equations on a subset of flat spacetime of the form [0,T] x B_R, where B_R is a
ball of radius R, and analyze different kinds of boundary conditions on
\partial B_R. Our main results are: i) We give an explicit analytic example
showing that boundary conditions obtained from freezing the incoming
characteristic fields to their initial values are not compatible with the
constraints. ii) With the help of the exact solutions constructed, we determine
the amount of artificial reflection of gravitational radiation from
constraint-preserving boundary conditions which freeze the Weyl scalar Psi_0 to
its initial value. For monochromatic radiation with wave number k and arbitrary
angular momentum number l >= 2, the amount of reflection decays as 1/(kR)^4 for
large kR. iii) For each L >= 2, we construct new local constraint-preserving
boundary conditions which perfectly absorb linearized radiation with l <= L.
(iv) We generalize our analysis to a weakly curved background of mass M, and
compute first order corrections in M/R to the reflection coefficients for
quadrupolar odd-parity radiation. For our new boundary condition with L=2, the
reflection coefficient is smaller than the one for the freezing Psi_0 boundary
condition by a factor of M/R for kR > 1.04. Implications of these results for
numerical simulations of binary black holes on finite domains are discussed.Comment: minor revisions, 30 pages, 6 figure
- âŠ