Researcher profile

William D. Henshaw

William D. Henshaw contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 15 - UnverifiedVerification L1Unclaimed author
3works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

3 published item(s)

preprint2022arXiv

Local Compatibility Boundary Conditions for High-Order Accurate Finite-Difference Approximations of PDEs

We describe a new approach to derive numerical approximations of boundary conditions for high-order accurate finite-difference approximations. The approach, called the Local Compatibility Boundary Condition (LCBC) method, uses boundary conditions and compatibility boundary conditions derived from the governing equations, as well as interior and boundary grid values, to construct a local polynomial, whose degree matches the order of accuracy of the interior scheme, centered at each boundary point. The local polynomial is then used to derive a discrete formula for each ghost point in terms of the data. This approach leads to centered approximations that are generally more accurate and stable than one-sided approximations. Moreover, the stencil approximations are local since they do not couple to neighboring ghost-point values which can occur with traditional compatibility conditions. The local polynomial is derived using continuous operators and derivatives which enables the automatic construction of stencil approximations at different orders of accuracy. The LCBC method is developed here for problems governed by second-order partial differential equations, and it is verified for a wide range of sample problems, both time-dependent and time-independent, in two space dimensions and for schemes up to sixth-order accuracy.

preprint2021arXiv

High-order accurate schemes for Maxwell's equations with nonlinear active media and material interfaces

We describe a fourth-order accurate finite-difference time-domain scheme for solving dispersive Maxwell's equations with nonlinear multi-level carrier kinetics models. The scheme is based on an efficient single-step three time-level modified equation approach for Maxwell's equations in second-order form for the electric field coupled to ODEs for the polarization vectors and population densities of the atomic levels. The resulting scheme has a large CFL-one time-step. Curved interfaces between different materials are accurately treated with curvilinear grids and compatibility conditions. A novel hierarchical modified equation approach leads to an explicit scheme that does not require any nonlinear iterations. The hierarchical approach at interfaces leads to local updates at the interface with no coupling in the tangential directions. Complex geometry is treated with overset grids. Numerical stability is maintained using high-order upwind dissipation designed for Maxwell's equations in second-order form. The scheme is carefully verified for a number of two and three-dimensional problems. The resulting numerical model with generalized dispersion and arbitrary nonlinear multi-level system can be used for many plasmonic applications such as for ab initio time domain modeling of nonlinear engineered materials for nanolasing applications, where nano-patterned plasmonic dispersive arrays are used to enhance otherwise weak nonlinearity in the active media.

preprint2020arXiv

Multigrid with Nonstandard Coarsening

We consider the numerical solution of Poisson's equation on structured grids using geometric multigrid with nonstandard coarse grids and coarse level operators. We are motivated by the problem of developing high-order accurate numerical solvers for elliptic boundary value problems on complex geometry using overset grids. Overset grids are typically dominated by large Cartesian background grids and thus fast solvers for Cartesian grids are highly desired. For flexibility in grid generation we would like to consider coarsening factors other than two, and lower-order accurate coarse-level approximations. We show that second-order accurate coarse-level approximations are very effective for fourth- or sixth-order accurate fine-level finite difference discretizations. We study the use of different Galerkin and non-Galerkin coarse-level operators. We use red-black smoothers with a relaxation parameter $ω$. Using local Fourier analysis we choose $ω$ and the coarse-level operators to optimize the overall multigrid convergence rate. Motivated by the use of red-black smoothers in one dimension that can result in a direct solver for the standard second-order accurate discretization to Poisson's equation, we show that this direct-solver property can be extended to two dimensions using a rotated grid that results from red-black coarsening. We evaluate the use of red-black coarsening in more general settings. We also study grid coarsening by a general factor and show that good convergence rates are retained for a range of coarsening factors near two. We ask the question of which coarsening factor leads to the most efficient algorithm.