  Extreme scattering events (ESEs) are observed as dramatic ($>50\%$) drops in
flux density that occur over an extended period of weeks to months. Discrete
plasma lensing structures are theorized to scatter the radio waves produced by
distant sources such as pulsars, causing the signature decrease in flux density
and characteristic caustic spikes in ESE light curves. While plasma lens models
in the extant literature have reproduced key features of ESE light curves, they
have all faced the problem of being highly over-dense and over-pressured
relative to the surrounding interstellar medium (ISM) by orders of magnitude.
We model ESEs by numerically ray-tracing through analytic, volumetric plasma
lens models by solving the eikonal equation. Delaunay triangulation connecting
the rays approximates the wavefront, generating a mapping from the observer
plane to the source plane to account for multiple-imaging. This eikonal method
of ray-tracing is tested against known analytic solutions and is then applied
to a three-dimensional Gaussian-distributed electron volume density lens, and a
filament model inspired by Grafton et al. (2023). We find convergence of our
numerical results with established analytic solutions validating our numerical
method, and reproduce ESE-like light curves. Our numerical ray-tracing method
lends itself well to exploring the lensing effects of volumetric turbulence as
well as sheet-like lenses, which is currently in progress.
