Study of Georgian Cladoniaceae

Through updating of the identifications of 1306 specimens housed at the largest herbarium collection in Tbilisi (TBI) and some new collections made in the summer of 2015 along the Greater Caucasus Range in Georgia, a comprehensive list and a key of forty three Georgian Cladonia species are presented. Cladonia acuminata, C.bacilliformis, C. borealis, C. conista, C. cyanipes, C. cyathomorpha, C. cf. libifera, C. macrophyllodes, C. sulphurina, and C. symphycarpa are reported as new to Georgia. Thirteen species of Cladonia were deleted from the earlier checklists.Peer reviewe

Properties of a method of fundamental solutions for the parabolic heat equation

We show that a set of fundamental solutions to the parabolic heat equation, with each element in the set corresponding to a point source located on a given surface with the number of source points being dense on this surface, constitute a linearly independent and dense set with respect to the standard inner product of square integrable functions, both on lateral- and time-boundaries. This result leads naturally to a method of numerically approximating solutions to the parabolic heat equation denoted a method of fundamental solutions (MFS). A discussion around convergence of such an approximation is included

Integral Equations for Heat Kernel in Compound Media

By making use of the potentials of the heat conduction equation the integral equations are derived which determine the heat kernel for the Laplace operator $-a^2\Delta$ in the case of compound media. In each of the media the parameter $a^2$ acquires a certain constant value. At the interface of the media the conditions are imposed which demand the continuity of the `temperature' and the `heat flows'. The integration in the equations is spread out only over the interface of the media. As a result the dimension of the initial problem is reduced by 1. The perturbation series for the integral equations derived are nothing else as the multiple scattering expansions for the relevant heat kernels. Thus a rigorous derivation of these expansions is given. In the one dimensional case the integral equations at hand are solved explicitly (Abel equations) and the exact expressions for the regarding heat kernels are obtained for diverse matching conditions. Derivation of the asymptotic expansion of the integrated heat kernel for a compound media is considered by making use of the perturbation series for the integral equations obtained. The method proposed is also applicable to the configurations when the same medium is divided, by a smooth compact surface, into internal and external regions, or when only the region inside (or outside) this surface is considered with appropriate boundary conditions.Comment: 26 pages, no figures, no tables, REVTeX4; two items are added into the Reference List; a new section is added, a version that will be published in J. Math. Phy

Three-dimensional desingularized boundary integral methods for potential problems

The concept of desingularization in three-dimensional boundary integral computations is re-examined. The boundary integral equation is desingularized by moving the singular points away from the boundary and outside the problem domain. We show that the desingularization gives better solutions to several problems. As a result of desingularization, the surface integrals can be evaluated by simpler techniques, speeding up the computation. The effects of the desingularization distance on the solution and the condition of the resulting system of algebraic equations are studied for both direct and indirect versions of the boundary integral method. Computations show that a broad range of desingularization distances gives accurate solutions with significant savings in the computation time. The desingularization distance must be carefully linked to the mesh size to avoid problems with uniqueness and ill-conditioning. As an example, the desingularized indirect approach is tested on unsteady non-linear three-dimensional gravity waves generated by a moving submerged disturbance; minimal computational difficulties are encountered at the truncated boundary.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/50203/1/1650120807_ftp.pd

Potentials of Gaussians and approximate wavelets

We derive new formulas for harmonic, diffraction, elastic, and hydrodynamic potentials acting on anisotropic Gaussians and approximate wavelets. These formulas can be used to construct accurate cubature formulas for these potentials

The method of fundamental solutions for three-dimensional inverse geometric elasticity problems

We investigate the numerical reconstruction of smooth star-shaped voids (rigid inclusions and cavities) which are compactly contained in a three-dimensional isotropic linear elastic medium from a single set of Cauchy data (i.e. nondestructive boundary displacement and traction measurements) on the accessible outer boundary. This inverse geometric problem in three-dimensional elasticity is approximated using the method of fundamental solutions (MFS). The parameters describing the boundary of the unknown void, its centre, and the contraction and dilation factors employed for selecting the fictitious surfaces where the MFS sources are to be positioned, are taken as unknowns of the problem. In this way, the original inverse geometric problem is reduced to finding the minimum of a nonlinear least-squares functional that measures the difference between the given and computed data, penalized with respect to both the MFS constants and the derivative of the radial coordinates describing the position of the star-shaped void. The interior source points are anchored and move with the void during the iterative reconstruction procedure. The feasibility of this new method is illustrated in several numerical examples

A new fast multi-domain BEM to model seismic wave propagation and amplification in 3D geological structures

International audienceThe analysis of seismic wave propagation and amplification in complex geological structures raises the need for efficient and accurate numerical methods. The solution of the elastodynamic equations using traditional boundary element methods (BEMs) is greatly hindered by the fully-populated nature of the matrix equations arising from the discretization. In a previous study limited to homogeneous media, the present authors have established that the Fast Multipole (FM) method reduces the complexity of a 3-D elastodynamic BEM to $N \log N$ per GMRES iteration and demonstrated its effectiveness on 3-D canyon configurations. In this article, the frequency-domain FM-BEM methodology is extented to 3-D elastic wave propagation in piecewise-homogeneous domains in the form of a FM-accelerated multi-region BE-BE coupling approach. This new method considerably enhances the capability of the BEM for studying the propagation of seismic waves in 3-D alluvial basins of arbitrary geometry embedded in semi-infinite media. Several fully 3-D examples (oblique SV-waves) representative of such configurations validate and demonstrate the capabilities of the multi-domain fast multipole approach. They include comparisons with available (low-frequency) results for various types of incident wavefields, and time-domain results obtained by means of Fourier synthesis