Search results for "numerical approximation"
showing 6 items of 16 documents
The average over a sphere
1980
Abstract The N points ri and the N segments ΔΩi of the unit sphere used in the numerical approximation of the average over the sphere are optimized to approximate the average of the set of spherical harmonics {;Yl,m;l = 0, 1, 2, …, L}; up to L = 18. The symmetry of f( r ) can be taken into acount by using only a distinct subquantity of the N point {; r i , ΔΩ i }; . Sets for N = 48n (n = 1, 2, …, 6) are tabulated. The advantage of the method is shown by the calculation of a powder Mossbauer spectrum including electric and magnetic hyperfine interactions.
A brief overview on the numerical behavior of an implicit meshless method and an outlook to future challenges
2015
In this paper recent results on a leapfrog ADI meshless formulation are reported and some future challenges are addressed. The method benefits from the elimination of the meshing task from the pre-processing stage in space and it is unconditionally stable in time. Further improvements come from the ease of implementation, which makes computer codes very flexible in contrast to mesh based solver ones. The method requires only nodes at scattered locations and a function and its derivatives are approximated by means of a kernel representation. A perceived obstacle in the implicit formulation is in the second order differentiations which sometimes are eccesively sensitive to the node configurat…
A novel numerical meshless approach for electric potential estimation in transcranial stimulation
2015
In this paper, a first application of the method of fundamental solutions in estimating the electric potential and the spatial current density distribution in the brain due to transcranial stimulation, is presented. The coupled boundary value p roblems for the electric potential are solved in a meshless way, so avoiding the use of grid based numerical methods. A multi-spherical geometry is considered and numerical results are discussed.
Accretion shock on CTTSs and its X-ray emission
2009
High spectral resolution X-ray observations of classical T Tauri stars (CTTSs) demonstrate the presence of plasma at T~2-3×10^6 K and ne~10^11-10^13 cm-3. Stationary models suggest that this emission is due to shock-heated accreting material. We address this issue by a 1-D hydrodynamic model of the impact of the accretion flow onto a chromosphere of a CTTS with the aim of investigating the stability of accretion shock and the role of the chromosphere. Our simulations include the effects of gravity, radiative losses from optically thin plasma, the thermal conduction and a detailed modeling of the stellar chromosphere. Here we present the results of a simulation based on the parameters of the…
Discontinuous Galerkin semi-Lagrangian method for Vlasov-Poisson
2011
We present a discontinuous Galerkin scheme for the numerical approximation of the one-dimensional periodic Vlasov-Poisson equation. The scheme is based on a Galerkin-characteristics method in which the distribution function is projected onto a space of discontinuous functions. We present comparisons with a semi-Lagrangian method to emphasize the good behavior of this scheme when applied to Vlasov-Poisson test cases. Une méthode de Galerkin discontinu est proposée pour l’approximation numérique de l’équation de Vlasov-Poisson 1D. L’approche est basée sur une méthode Galerkin-caractéristiques où la fonction de distribution est projetée sur un espace de fonctions discontinues. En particulier, …
Higher-order Hamilton–Jacobi perturbation theory for anisotropic heterogeneous media: dynamic ray tracing in Cartesian coordinates
2018
With a Hamilton–Jacobi equation in Cartesian coordinates as a starting point, it is common to use a system of ordinary differential equations describing the continuation of first-order derivatives of phase-space perturbations along a reference ray. Such derivatives can be exploited for calculating geometrical spreading on the reference ray and for establishing a framework for second-order extrapolation of traveltime to points outside the reference ray. The continuation of first-order derivatives of phase-space perturbations has historically been referred to as dynamic ray tracing. The reason for this is its importance in the process of calculating amplitudes along the reference ray. We exte…