On the solution of the Helmholtz equation on regions with corners.

In this paper we solve several boundary value problems for the Helmholtz equation on polygonal domains. We observe that when the problems are formulated as the boundary integral equations of potential theory, the solutions are representable by series of appropriately chosen Bessel functions. In addition to being analytically perspicuous, the resulting expressions lend themselves to the construction of accurate and efficient numerical algorithms. The results are illustrated by a number of numerical examples.

Optimal schedules of light exposure for rapidly correcting circadian misalignment.

Jet lag arises from a misalignment of circadian biological timing with the timing of human activity, and is caused by rapid transmeridian travel. Jet lag's symptoms, such as depressed cognitive alertness, also arise from work and social schedules misaligned with the timing of the circadian clock. Using experimentally validated mathematical models, we develop a new methodology to find mathematically optimal schedules of light exposure and avoidance for rapidly re-entraining the human circadian system. In simulations, our schedules are found to significantly outperform other recently proposed schedules. Moreover, our schedules appear to be significantly more robust to both noise in light and to inter-individual variations in endogenous circadian period than other proposed schedules. By comparing the optimal schedules for thousands of different situations, and by using general mathematical arguments, we are also able to translate our findings into general principles of optimal circadian re-entrainment. These principles include: 1) a class of schedules where circadian amplitude is only slightly perturbed, optimal for dim light and for small shifts 2) another class of schedules where shifting occurs along the shortest path in phase-space, optimal for bright light and for large shifts 3) the determination that short light pulses are less effective than sustained light if the goal is to re-entrain quickly, and 4) the determination that length of daytime should be significantly shorter when delaying the clock than when advancing it.

On the efficient evaluation of the azimuthal Fourier components of the Green's function for Helmholtz's equation in cylindrical coordinates.

In this paper, we develop an efficient algorithm to evaluate the azimuthal Fourier components of the Green's function for the Helmholtz equation in cylindrical coordinates. A computationally efficient algorithm for this modal Green's function is essential for solvers for electromagnetic scattering from bodies of revolution (e.g., radar cross sections, antennas). Current algorithms to evaluate this modal Green's function become computationally intractable when the source and target are close or when the wavenumber is large or complex. Furthermore, most state-of-the-art methods cannot be easily parallelized. In this paper, we present an algorithm for evaluating the modal Green's function that has performance independent of both source-to-target proximity and wavenumber, and whose cost grows as O(m), where m is the Fourier mode. Our algorithm's performance is independent of whether the wavenumber is real or complex. Furthermore, our algorithm is embarrassingly parallelizable.