A Coupled Spatial-Network Model: A Mathematical Framework for Applications in Epidemiology.

There is extensive evidence that network structure (e.g., air transport, rivers, or roads) may significantly enhance the spread of epidemics into the surrounding geographical area. A new compartmental modeling framework is proposed which couples well-mixed (ODE in time) population centers at the vertices, 1D travel routes on the graph's edges, and a 2D continuum containing the rest of the population to simulate how an infection spreads through a population. The edge equations are coupled to the vertex ODEs through junction conditions, while the domain equations are coupled to the edges through boundary conditions. A numerical method based on spatial finite differences for the edges and finite elements in the 2D domain is described to approximate the model, and numerical verification of the method is provided. The model is illustrated on two simple and one complex example geometries, and a parameter study example is performed. The observed solutions exhibit exponential decay after a certain time has passed, and the cumulative infected population over the vertices, edges, and domain tends to a constant in time but varying in space, i.e., a steady state solution.

Finite-difference time-domain simulation of spacetime cloak.

In this work, we present a numerical method that remedies the instabilities of the conventional FDTD approach for solving Maxwell's equations in a space-time dependent magneto-electric medium with direct application to the simulation of the recently proposed spacetime cloak. We utilize a dual grid FDTD method overlapped in the time domain to provide a stable approach for the simulation of a magneto-electric medium with time and space varying permittivity, permeability and coupling coefficient. The developed method can be applied to explore other new physical possibilities offered by spacetime cloaking, metamaterials, and transformation optics.

Subpixel smoothing finite-difference time-domain method for material interface between dielectric and dispersive media.

In this Letter, we have shown that the subpixel smoothing technique that eliminates the staircasing error in the finite-difference time-domain method can be extended to material interface between dielectric and dispersive media by local coordinate rotation. First, we show our method is equivalent to the subpixel smoothing method for dielectric interface, then we extend it to a more general case where dispersive/dielectric interface is present. Finally, we provide a numerical example on a scattering problem to demonstrate that we were able to significantly improve the accuracy.

Kink shape solutions of the Maxwell-Lorentz system.

In the limit of high amplitude oscillating electromagnetic fields, a sequence of kink antikink shaped optical waves has been found in the Maxwell's equations coupled to a single Lorentz oscillator and with Kerr nonlinearity. The individual kinks and antikinks result from a traveling wave assumption and their stability has been assessed by numerical simulations. For typical physical parameter values the kink width is of the order of tens of femtoseconds.