skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Optimal limiters for the spectral element method.

Authors:
; ;
Publication Date:
Research Org.:
Sandia National Laboratories
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1072624
Report Number(s):
SAND2013-0057J
DOE Contract Number:
AC04-94AL85000
Resource Type:
Journal Article
Resource Relation:
Journal Name: Proposed for publication in Journal of Computational Physics.
Country of Publication:
United States
Language:
English

Citation Formats

Guba, Oksana, Taylor, Mark A., and St-Cyr, Amik. Optimal limiters for the spectral element method.. United States: N. p., 2013. Web.
Guba, Oksana, Taylor, Mark A., & St-Cyr, Amik. Optimal limiters for the spectral element method.. United States.
Guba, Oksana, Taylor, Mark A., and St-Cyr, Amik. Tue . "Optimal limiters for the spectral element method.". United States. doi:.
@article{osti_1072624,
title = {Optimal limiters for the spectral element method.},
author = {Guba, Oksana and Taylor, Mark A. and St-Cyr, Amik},
abstractNote = {},
doi = {},
journal = {Proposed for publication in Journal of Computational Physics.},
number = ,
volume = ,
place = {United States},
year = {Tue Jan 01 00:00:00 EST 2013},
month = {Tue Jan 01 00:00:00 EST 2013}
}
  • We have solved the problem of layer-by-layer laser-light dosimetry in biological tissues and of selecting an individual therapeutic dose in laser therapy. A method is proposed for real-time monitoring of the radiation density in tissue layers in vivo, concentrations of its endogenous (natural) and exogenous (specially administered) chromophores, as well as in-depth distributions of the spectrum of light action on these chromophores. As the background information use is made of the spectrum of diffuse light reflected from a patient's tissue, measured by a fibre-optic spectrophotometer. The measured spectrum is quantitatively analysed by the method of approximating functions for fluxes ofmore » light multiply scattered in tissue and by a semi-analytical method for calculating the in-depth distribution of the light flux in a multi-layered medium. We have shown the possibility of employing the developed method for monitoring photosensitizer and oxyhaemoglobin concentrations in tissue, light power absorbed by chromophores in tissue layers at different depths and laser-induced changes in the tissue morphology (vascular volume content and ratios of various forms of haemoglobin) during photodynamic therapy. (biophotonics)« less
  • Design sensitivity analysis, along with the shape optimization of heat diffusion problems using the boundary element method (BEM), is presented in this paper. The present approach utilizes the implicit differentiation of discretized boundary integral equations with respect to the design variables to yield the sensitivity equations. A technique based on the response of an object to a constant boundary temperature is presented for the evaluation of singular terms in the thermal sensitivity kernels. A procedure for the design sensitivity analysis of a reduced system of equations obtained via substructuring and condensation is also presented. The BEM formulations are implemented formore » both two-dimensional and axisymmetric objects. A number of sample problems are solved to demonstrate the accuracy of the present sensitivity formulation and to obtain optimal configurations of some mechanical components of practical interest, which are subjected to different thermal environments.« less
  • Uniform estimates for the error in the finite element method are derived for a model problem on a general triangular mesh in two dimensions. These are optimal if the degree of the piecewise polynomials is greater than one. Similar estimates of the error are also derived in L/sup p/. As an intermediate step, an L/sup 1/ estimate of the gradient of the error in the finite element approximation of the Green's function is proved that is optimal for all degrees.
  • Highlights: •New strategies to improve the accuracy of the reconstruction through mesh and finite element parameterization. •Use of gradient filtering through an alternative inner product within the adjoint method. •An integral form of the cost function is used to make the reconstruction compatible with all finite element formulations, continuous and discontinuous. •Gradient-based algorithm with the adjoint method is used for the reconstruction. -- Abstract: Optical tomography is mathematically treated as a non-linear inverse problem where the optical properties of the probed medium are recovered through the minimization of the errors between the experimental measurements and their predictions with a numericalmore » model at the locations of the detectors. According to the ill-posed behavior of the inverse problem, some regularization tools must be performed and the Tikhonov penalization type is the most commonly used in optical tomography applications. This paper introduces an optimized approach for optical tomography reconstruction with the finite element method. An integral form of the cost function is used to take into account the surfaces of the detectors and make the reconstruction compatible with all finite element formulations, continuous and discontinuous. Through a gradient-based algorithm where the adjoint method is used to compute the gradient of the cost function, an alternative inner product is employed for preconditioning the reconstruction algorithm. Moreover, appropriate re-parameterization of the optical properties is performed. These regularization strategies are compared with the classical Tikhonov penalization one. It is shown that both the re-parameterization and the use of the Sobolev cost function gradient are efficient for solving such an ill-posed inverse problem.« less
  • In this paper, we study a coupling of spectral methods and the [rho]-version finite element methods for elliptic boundary value problems containing singularities. The method of auxiliary mapping, which is a recent development to deal with domain singularities in the [rho]-version of the finite element method, is employed to remove the pollution effect caused by singularities. An iterative interfacial coupling between spectral methods and the [rho]-version of the finite element method is used and investigated numerically. The advantages of such an approach are demonstrated by the high accuracy of spectral methods for the smooth part of solutions and the flexibilitymore » of the [rho]-version of the finite element method for dealing with singularities and irregular domains. The efficiency of the coupling method is also evaluated by comparing results obtained by this method with those obtained by the full finite element algorithm. 26 refs., 9 figs., 6 tabs.« less