Fast numerical methods for mixed, singular Helmholtz boundary value problems and Laplace eigenvalue problems - with applications to antenna design, sloshing, electromagnetic scattering and spectral geometry

This thesis presents a novel class of algorithms for the solution of scattering and eigenvalue problems on general two-dimensional domains under a variety of boundary conditions, including non-smooth domains and certain "Zaremba" boundary conditions - for which Dirichlet and Neumann conditions are specified on various portions of the domain boundary. The theoretical basis of the methods for the Zaremba problems on smooth domains concern detailed information, which is put forth for the first time in this thesis, about the singularity structure of solutions of the Laplace operator under boundary conditions of Zaremba type. The new methods, which are based on use of Green functions and integral equations, incorporate a number of algorithmic innovations, including a fast and robust eigenvalue-search algorithm, use of the Fourier Continuation method for regularization of all smooth-domain Zaremba singularities, and newly derived quadrature rules which give rise to high-order convergence even around singular points for the Zaremba problem. The resulting algorithms enjoy high-order convergence, and they can tackle a variety of elliptic problems under general boundary conditions, including, for example, eigenvalue problems, scattering problems, and, in particular, eigenfunction expansion for time-domain problems in non-separable physical domains with mixed boundary conditions.

Time-dependent monoenergetic neutron transport in two adjacent semi-infinite media

An exact solution to the monoenergetic Boltzmann equation is obtained for the case of a plane isotropic burst of neutrons introduced at the interface separating two adjacent, dissimilar, semi-infinite media. The method of solution used is to remove the time dependence by a Laplace transformation, solve the transformed equation by the normal mode expansion method, and then invert to recover the time dependence.

The general result is expressed as a sum of definite, multiple integrals, one of which contains the uncollided wave of neutrons originating at the source plane. It is possible to obtain a simplified form for the solution at the interface, and certain numerical calculations are made there.

The interface flux in two adjacent moderators is calculated and plotted as a function of time for several moderator materials. For each case it is found that the flux decay curve has an asymptotic slope given accurately by diffusion theory. Furthermore, the interface current is observed to change directions when the scattering and absorption cross sections of the two moderator materials are related in a certain manner. More specifically, the reflection process in two adjacent moderators appears to depend initially on the scattering properties and for long times on the absorption properties of the media.

This analysis contains both the single infinite and semi-infinite medium problems as special cases. The results in these two special cases provide a check on the accuracy of the general solution since they agree with solutions of these problems obtained by separate analyses.

Spectral density of first order Piecewise linear systems excited by white noise

The Fokker-Planck (FP) equation is used to develop a general method for finding the spectral density for a class of randomly excited first order systems. This class consists of systems satisfying stochastic differential equations of form ẋ + f(x) = m/Ʃ/j = 1 h_j(x)n_j(t) where f and the h_j are piecewise linear functions (not necessarily continuous), and the n_j are stationary Gaussian white noise. For such systems, it is shown how the Laplace-transformed FP equation can be solved for the transformed transition probability density. By manipulation of the FP equation and its adjoint, a formula is derived for the transformed autocorrelation function in terms of the transformed transition density. From this, the spectral density is readily obtained. The method generalizes that of Caughey and Dienes, J. Appl. Phys., 32.11.

This method is applied to 4 subclasses: (1) m = 1, h₁ = const. (forcing function excitation); (2) m = 1, h₁ = f (parametric excitation); (3) m = 2, h₁ = const., h₂ = f, n₁ and n₂ correlated; (4) the same, uncorrelated. Many special cases, especially in subclass (1), are worked through to obtain explicit formulas for the spectral density, most of which have not been obtained before. Some results are graphed.

Dealing with parametrically excited first order systems leads to two complications. There is some controversy concerning the form of the FP equation involved (see Gray and Caughey, J. Math. Phys., 44.3); and the conditions which apply at irregular points, where the second order coefficient of the FP equation vanishes, are not obvious but require use of the mathematical theory of diffusion processes developed by Feller and others. These points are discussed in the first chapter, relevant results from various sources being summarized and applied. Also discussed is the steady-state density (the limit of the transition density as t → ∞).

Recorrido y reflexiones en torno al pensamiento analítico de Georges Laplace: movimiento, interdependencia y arquetipos en la construcción de una arqueología científica

Homenaje a Georges Laplace, realizado en Vitoria-Gasteiz el 13,14 y 15 de noviembre de 2012. Edición a cargo de Aitor Calvo, Aitor Sánchez, Maite García-Rojas y Mónica Alonso-Eguíluz.

Arqueozoología analítica. Otro ejemplo práctico derivado de la obra de Georges Laplace

Homenaje a Georges Laplace, realizado en Vitoria-Gasteiz el 13,14 y 15 de noviembre de 2012. Edición a cargo de Aitor Calvo, Aitor Sánchez, Maite García-Rojas y Mónica Alonso-Eguíluz.

New doubly periodic waves of the (2+1)-dimensional double sine-Gordon equation

New exact solutions of the (2 + 1)-dimensional double sine-Gordon equation are studied by introducing the modified mapping relations between the cubic nonlinear Klein-Gordon system and double sine-Gordon equation. Two arbitrary functions are included into the Jacobi elliptic function solutions. New doubly periodic wave solutions are obtained and displayed graphically by proper selections of the arbitrary functions.

Novel equation to determine the hepatic triglyceride concentration in humans by MRI: diagnosis and monitoring of NAFLD in obese patients before and after bariatric surgery

Background: Non-alcoholic fatty liver disease (NAFLD) is caused by abnormal accumulation of lipids within liver cells. Its prevalence is increasing in developed countries in association with obesity, and it represents a risk factor for non-alcoholic steatohepatitis (NASH), cirrhosis and hepatocellular carcinoma. Since NAFLD is usually asymptomatic at diagnosis, new non-invasive approaches are needed to determine the hepatic lipid content in terms of diagnosis, treatment and control of disease progression. Here, we investigated the potential of magnetic resonance imaging (MRI) to quantitate and monitor the hepatic triglyceride concentration in humans. Methods: A prospective study of diagnostic accuracy was conducted among 129 consecutive adult patients (97 obesity and 32 non-obese) to compare multi-echo MRI fat fraction, grade of steatosis estimated by histopathology, and biochemical measurement of hepatic triglyceride concentration (that is, Folch value). Results: MRI fat fraction positively correlates with the grade of steatosis estimated on a 0 to 3 scale by histopathology. However, this correlation value was stronger when MRI fat fraction was linked to the Folch value, resulting in a novel equation to predict the hepatic triglyceride concentration (mg of triglycerides/g of liver tissue = 5.082 + (432.104 * multi-echo MRI fat fraction)). Validation of this formula in 31 additional patients (24 obese and 7 controls) resulted in robust correlation between the measured and estimated Folch values. Multivariate analysis showed that none of the variables investigated improves the Folch prediction capacity of the equation. Obese patients show increased steatosis compared to controls using MRI fat fraction and Folch value. Bariatric surgery improved MRI fat fraction values and the Folch value estimated in obese patients one year after surgery. Conclusions: Multi-echo MRI is an accurate approach to determine the hepatic lipid concentration by using our novel equation, representing an economic non-invasive method to diagnose and monitor steatosis in humans.