Search results for "Computational Mathematic"
showing 10 items of 987 documents
Wavelet-like bases for thin-wire integral equations in electromagnetics
2005
AbstractIn this paper, wavelets are used in solving, by the method of moments, a modified version of the thin-wire electric field integral equation, in frequency domain. The time domain electromagnetic quantities, are obtained by using the inverse discrete fast Fourier transform. The retarded scalar electric and vector magnetic potentials are employed in order to obtain the integral formulation. The discretized model generated by applying the direct method of moments via point-matching procedure, results in a linear system with a dense matrix which have to be solved for each frequency of the Fourier spectrum of the time domain impressed source. Therefore, orthogonal wavelet-like basis trans…
Invariant rotational curves in Sitnikov's Problem
1993
The Sitnikov's Problem is a Restricted Three-Body Problem of Celestial Mechanics depending on a parameter, the eccentricity,e. The Hamiltonian,H(z, v, t, e), does not depend ont ife=0 and we have an integrable system; ife is small the KAM Theory proves the existence of invariant rotational curves, IRC. For larger eccentricities, we show that there exist two complementary sequences of intervals of values ofe that accumulate to the maximum admissible value of the eccentricity, 1, and such that, for one of the sequences IRC around a fixed point persist. Moreover, they shrink to the planez=0 ase tends to 1.
Comparison among three boundary element methods for torsion problems: CPM, CVBEM, LEM
2011
This paper provides solutions for De Saint-Venant torsion problem on a beam with arbitrary and uniform cross-section. In particular three methods framed into complex analysis have been considered: Complex Polynomial Method (CPM), Complex Variable Boundary Element Method (CVBEM) and Line Element-less Method (LEM), recently proposed. CPM involves the expansion of a complex potential in Taylor series, computing the unknown coefficients by means of collocation points on the boundary. CVBEM takes advantage of Cauchy’s integral formula that returns the solution of Laplace equation when mixed boundary conditions on both real and imaginary parts of the complex potential are known. LEM introduces th…
The Method of Fundamental Solutions in Solving Coupled Boundary Value Problems for M/EEG
2015
The estimation of neuronal activity in the human brain from electroencephalography (EEG) and magnetoencephalography (MEG) signals is a typical inverse problem whose solution pro- cess requires an accurate and fast forward solver. In this paper the method of fundamental solutions is, for the first time, proposed as a meshfree, boundary-type, and easy-to-implement alternative to the boundary element method (BEM) for solving the M/EEG forward problem. The solution of the forward problem is obtained by numerically solving a set of coupled boundary value problems for the three-dimensional Laplace equation. Numerical accuracy, convergence, and computational load are investigated. The proposed met…
Modelling of Systems with a Dispersed Phase: “Measuring” Small Sets in the Presence of Elliptic Operators
2016
When modelling systems with a dispersed phase involving elliptic operators, as is the case of the Stokes or Navier-Stokes problem or the heat equation in a bounded domain, the geometrical structure of the space occupied by the dispersed phase enters in the homogenization process through its capacity, a quantity which can be used to define the equivalence classes in \(H^1\). We shall review the relationship between capacity and homogenization terms in the limit when the number of inclusions becomes large, focusing in particular on the situation where the distribution of inclusions is not necessarily too regular (i.e. it is not periodic).
Polarization tensors of planar domains as functions of the admittivity contrast
2014
(Electric) polarization tensors describe part of the leading order term of asymptotic voltage perturbations caused by low volume fraction inhomogeneities of the electrical properties of a medium. They depend on the geometry of the support of the inhomogeneities and on their admittivity contrast. Corresponding asymptotic formulas are of particular interest in the design of reconstruction algorithms for determining the locations and the material properties of inhomogeneities inside a body from measurements of current flows and associated voltage potentials on the body's surface. In this work we consider the two-dimensional case only and provide an analytic representation of the polarization t…
Explicit Algorithms for a New Time Dependent Model Based on Level Set Motion for Nonlinear Deblurring and Noise Removal
2000
In this paper we formulate a time dependent model to approximate the solution to the nonlinear total variation optimization problem for deblurring and noise removal introduced by Rudin and Osher [ Total variation based image restoration with free local constraints, in Proceedings IEEE Internat. Conf. Imag. Proc., IEEE Press, Piscataway, NJ, (1994), pp. 31--35] and Rudin, Osher, and Fatemi [ Phys. D, 60 (1992), pp. 259--268], respectively. Our model is based on level set motion whose steady state is quickly reached by means of an explicit procedure based on Roe's scheme [ J. Comput. Phys., 43 (1981), pp. 357--372], used in fluid dynamics. We show numerical evidence of the speed of resolution…
SIMPRE: A software package to calculate crystal field parameters, energy levels, and magnetic properties on mononuclear lanthanoid complexes based on…
2013
This work presents a fortran77 code based on an effective electrostatic model of point charges around a rare earth ion. The program calculates the full set of crystal field parameters, energy levels spectrum, and wave functions, as well as the magnetic properties such as the magnetization, the temperature dependence of the magnetic susceptibility, and the Schottky contribution to the specific heat. It is designed for real systems that need not bear ideal symmetry and it is able to determine the easy axis of magnetization. Its systematic application to different coordination environments allows magneto-structural studies. The package has already been successfully applied to several mononucle…
The minimal model of Hahn for the Calvin cycle.
2018
There are many models of the Calvin cycle of photosynthesis in the literature. When investigating the dynamics of these models one strategy is to look at the simplest possible models in order to get the most detailed insights. We investigate a minimal model of the Calvin cycle introduced by Hahn while he was pursuing this strategy. In a variant of the model not including photorespiration it is shown that there exists exactly one positive steady state and that this steady state is unstable. For generic initial data either all concentrations tend to infinity at lates times or all concentrations tend to zero at late times. In a variant including photorespiration it is shown that for suitable v…
A Nonlinear Primal-Dual Method for Total Variation-Based Image Restoration
1999
We present a new method for solving total variation (TV) minimization problems in image restoration. The main idea is to remove some of the singularity caused by the nondifferentiability of the quantity $|\nabla u|$ in the definition of the TV-norm before we apply a linearization technique such as Newton's method. This is accomplished by introducing an additional variable for the flux quantity appearing in the gradient of the objective function, which can be interpreted as the normal vector to the level sets of the image u. Our method can be viewed as a primal-dual method as proposed by Conn and Overton [ A Primal-Dual Interior Point Method for Minimizing a Sum of Euclidean Norms, preprint,…