14,583 research outputs found
The diffuse Nitsche method: Dirichlet constraints on phase-field boundaries
We explore diffuse formulations of Nitsche's method for consistently imposing Dirichlet boundary conditions on phase-field approximations of sharp domains. Leveraging the properties of the phase-field gradient, we derive the variational formulation of the diffuse Nitsche method by transferring all integrals associated with the Dirichlet boundary from a geometrically sharp surface format in the standard Nitsche method to a geometrically diffuse volumetric format. We also derive conditions for the stability of the discrete system and formulate a diffuse local eigenvalue problem, from which the stabilization parameter can be estimated automatically in each element. We advertise metastable phase-field solutions of the Allen-Cahn problem for transferring complex imaging data into diffuse geometric models. In particular, we discuss the use of mixed meshes, that is, an adaptively refined mesh for the phase-field in the diffuse boundary region and a uniform mesh for the representation of the physics-based solution fields. We illustrate accuracy and convergence properties of the diffuse Nitsche method and demonstrate its advantages over diffuse penalty-type methods. In the context of imaging based analysis, we show that the diffuse Nitsche method achieves the same accuracy as the standard Nitsche method with sharp surfaces, if the inherent length scales, i.e., the interface width of the phase-field, the voxel spacing and the mesh size, are properly related. We demonstrate the flexibility of the new method by analyzing stresses in a human vertebral body
A literature survey of low-rank tensor approximation techniques
During the last years, low-rank tensor approximation has been established as
a new tool in scientific computing to address large-scale linear and
multilinear algebra problems, which would be intractable by classical
techniques. This survey attempts to give a literature overview of current
developments in this area, with an emphasis on function-related tensors
Phase-field boundary conditions for the voxel finite cell method: surface-free stress analysis of CT-based bone structures
The voxel finite cell method employs unfitted finite element meshes and voxel quadrature rules to seamlessly
transfer CT data into patient-specific bone discretizations. The method, however, still requires the explicit
parametrization of boundary surfaces to impose traction and displacement boundary conditions, which
constitutes a potential roadblock to automation. We explore a phase-field based formulation for imposing
traction and displacement constraints in a diffuse sense. Its essential component is a diffuse geometry model
generated from metastable phase-field solutions of the Allen-Cahn problem that assumes the imaging data as
initial condition. Phase-field approximations of the boundary and its gradient are then employed to transfer
all boundary terms in the variational formulation into volumetric terms. We show that in the context of the
voxel finite cell method, diffuse boundary conditions achieve the same accuracy as boundary conditions
defined over explicit sharp surfaces, if the inherent length scales, i.e., the interface width of the phase-field,
the voxel spacing and the mesh size, are properly related. We demonstrate the flexibility of the new method
by analyzing stresses in a human femur and a vertebral body
Tensor Numerical Methods in Quantum Chemistry: from Hartree-Fock Energy to Excited States
We resume the recent successes of the grid-based tensor numerical methods and
discuss their prospects in real-space electronic structure calculations. These
methods, based on the low-rank representation of the multidimensional functions
and integral operators, led to entirely grid-based tensor-structured 3D
Hartree-Fock eigenvalue solver. It benefits from tensor calculation of the core
Hamiltonian and two-electron integrals (TEI) in complexity using
the rank-structured approximation of basis functions, electron densities and
convolution integral operators all represented on 3D
Cartesian grids. The algorithm for calculating TEI tensor in a form of the
Cholesky decomposition is based on multiple factorizations using algebraic 1D
``density fitting`` scheme. The basis functions are not restricted to separable
Gaussians, since the analytical integration is substituted by high-precision
tensor-structured numerical quadratures. The tensor approaches to
post-Hartree-Fock calculations for the MP2 energy correction and for the
Bethe-Salpeter excited states, based on using low-rank factorizations and the
reduced basis method, were recently introduced. Another direction is related to
the recent attempts to develop a tensor-based Hartree-Fock numerical scheme for
finite lattice-structured systems, where one of the numerical challenges is the
summation of electrostatic potentials of a large number of nuclei. The 3D
grid-based tensor method for calculation of a potential sum on a lattice manifests the linear in computational work, ,
instead of the usual scaling by the Ewald-type approaches
Asymptotic power of sphericity tests for high-dimensional data
This paper studies the asymptotic power of tests of sphericity against
perturbations in a single unknown direction as both the dimensionality of the
data and the number of observations go to infinity. We establish the
convergence, under the null hypothesis and contiguous alternatives, of the log
ratio of the joint densities of the sample covariance eigenvalues to a Gaussian
process indexed by the norm of the perturbation. When the perturbation norm is
larger than the phase transition threshold studied in Baik, Ben Arous and Peche
[Ann. Probab. 33 (2005) 1643-1697] the limiting process is degenerate, and
discrimination between the null and the alternative is asymptotically certain.
When the norm is below the threshold, the limiting process is nondegenerate,
and the joint eigenvalue densities under the null and alternative hypotheses
are mutually contiguous. Using the asymptotic theory of statistical
experiments, we obtain asymptotic power envelopes and derive the asymptotic
power for various sphericity tests in the contiguity region. In particular, we
show that the asymptotic power of the Tracy-Widom-type tests is trivial (i.e.,
equals the asymptotic size), whereas that of the eigenvalue-based likelihood
ratio test is strictly larger than the size, and close to the power envelope.Comment: Published in at http://dx.doi.org/10.1214/13-AOS1100 the Annals of
Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical
Statistics (http://www.imstat.org
Fluctuating volume-current formulation of electromagnetic fluctuations in inhomogeneous media: incandecence and luminescence in arbitrary geometries
We describe a fluctuating volume--current formulation of electromagnetic
fluctuations that extends our recent work on heat exchange and Casimir
interactions between arbitrarily shaped homogeneous bodies [Phys. Rev. B. 88,
054305] to situations involving incandescence and luminescence problems,
including thermal radiation, heat transfer, Casimir forces, spontaneous
emission, fluorescence, and Raman scattering, in inhomogeneous media. Unlike
previous scattering formulations based on field and/or surface unknowns, our
work exploits powerful techniques from the volume--integral equation (VIE)
method, in which electromagnetic scattering is described in terms of
volumetric, current unknowns throughout the bodies. The resulting trace
formulas (boxed equations) involve products of well-studied VIE matrices and
describe power and momentum transfer between objects with spatially varying
material properties and fluctuation characteristics. We demonstrate that thanks
to the low-rank properties of the associatedmatrices, these formulas are
susceptible to fast-trace computations based on iterative methods, making
practical calculations tractable. We apply our techniques to study thermal
radiation, heat transfer, and fluorescence in complicated geometries, checking
our method against established techniques best suited for homogeneous bodies as
well as applying it to obtain predictions of radiation from complex bodies with
spatially varying permittivities and/or temperature profiles
- …