110 research outputs found
An algorithm for computing the 2D structure of fast rotating stars
Stars may be understood as self-gravitating masses of a compressible fluid
whose radiative cooling is compensated by nuclear reactions or gravitational
contraction. The understanding of their time evolution requires the use of
detailed models that account for a complex microphysics including that of
opacities, equation of state and nuclear reactions. The present stellar models
are essentially one-dimensional, namely spherically symmetric. However, the
interpretation of recent data like the surface abundances of elements or the
distribution of internal rotation have reached the limits of validity of
one-dimensional models because of their very simplified representation of
large-scale fluid flows. In this article, we describe the ESTER code, which is
the first code able to compute in a consistent way a two-dimensional model of a
fast rotating star including its large-scale flows. Compared to classical 1D
stellar evolution codes, many numerical innovations have been introduced to
deal with this complex problem. First, the spectral discretization based on
spherical harmonics and Chebyshev polynomials is used to represent the 2D
axisymmetric fields. A nonlinear mapping maps the spheroidal star and allows a
smooth spectral representation of the fields. The properties of Picard and
Newton iterations for solving the nonlinear partial differential equations of
the problem are discussed. It turns out that the Picard scheme is efficient on
the computation of the simple polytropic stars, but Newton algorithm is
unsurpassed when stellar models include complex microphysics. Finally, we
discuss the numerical efficiency of our solver of Newton iterations. This
linear solver combines the iterative Conjugate Gradient Squared algorithm
together with an LU-factorization serving as a preconditionner of the Jacobian
matrix.Comment: 40 pages, 12 figures, accepted in J. Comput. Physic
Harmonic density interpolation methods for high-order evaluation of Laplace layer potentials in 2D and 3D
We present an effective harmonic density interpolation method for the
numerical evaluation of singular and nearly singular Laplace boundary integral
operators and layer potentials in two and three spatial dimensions. The method
relies on the use of Green's third identity and local Taylor-like
interpolations of density functions in terms of harmonic polynomials. The
proposed technique effectively regularizes the singularities present in
boundary integral operators and layer potentials, and recasts the latter in
terms of integrands that are bounded or even more regular, depending on the
order of the density interpolation. The resulting boundary integrals can then
be easily, accurately, and inexpensively evaluated by means of standard
quadrature rules. A variety of numerical examples demonstrate the effectiveness
of the technique when used in conjunction with the classical trapezoidal rule
(to integrate over smooth curves) in two-dimensions, and with a Chebyshev-type
quadrature rule (to integrate over surfaces given as unions of non-overlapping
quadrilateral patches) in three-dimensions
Spectral tensor-train decomposition
The accurate approximation of high-dimensional functions is an essential task
in uncertainty quantification and many other fields. We propose a new function
approximation scheme based on a spectral extension of the tensor-train (TT)
decomposition. We first define a functional version of the TT decomposition and
analyze its properties. We obtain results on the convergence of the
decomposition, revealing links between the regularity of the function, the
dimension of the input space, and the TT ranks. We also show that the
regularity of the target function is preserved by the univariate functions
(i.e., the "cores") comprising the functional TT decomposition. This result
motivates an approximation scheme employing polynomial approximations of the
cores. For functions with appropriate regularity, the resulting
\textit{spectral tensor-train decomposition} combines the favorable
dimension-scaling of the TT decomposition with the spectral convergence rate of
polynomial approximations, yielding efficient and accurate surrogates for
high-dimensional functions. To construct these decompositions, we use the
sampling algorithm \texttt{TT-DMRG-cross} to obtain the TT decomposition of
tensors resulting from suitable discretizations of the target function. We
assess the performance of the method on a range of numerical examples: a
modifed set of Genz functions with dimension up to , and functions with
mixed Fourier modes or with local features. We observe significant improvements
in performance over an anisotropic adaptive Smolyak approach. The method is
also used to approximate the solution of an elliptic PDE with random input
data. The open source software and examples presented in this work are
available online.Comment: 33 pages, 19 figure
Hybrid collocation perturbation for PDEs with random domains
In this work we consider the problem of approximating the statistics of a
given Quantity of Interest (QoI) that depends on the solution of a linear
elliptic PDE defined over a random domain parameterized by random
variables. The random domain is split into large and small variations
contributions. The large variations are approximated by applying a sparse grid
stochastic collocation method. The small variations are approximated with a
stochastic collocation-perturbation method. Convergence rates for the variance
of the QoI are derived and compared to those obtained in numerical experiments.
Our approach significantly reduces the dimensionality of the stochastic
problem. The computational cost of this method increases at most quadratically
with respect to the number of dimensions of the small variations. Moreover, for
the case that the small and large variations are independent the cost increases
linearly
- âŠ