Researcher profile

W. D. Henshaw

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

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - Baseline
5works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

5 published item(s)

preprint2016arXiv

Computed Flow and Fluorescence Over the Ocular Surface

Fluorescein is perhaps the most commonly used substance to visualize tear film thickness and dynamics; better understanding of this process aids understanding of dry eye syndrome which afflicts millions of people. We study a mathematical model for tear film flow, evaporation, solutal transport and fluorescence over the exposed ocular surface during the interblink. Transport of the fluorescein ion by fluid flow in the tear film affects the intensity of fluorescence via changes in concentration and tear film thickness. Evaporation causes increased osmolarity and potential irritation over the ocular surface; it also alters fluorescein concentration and thus fluorescence. Using thinning rates from in vivo measurements together with thin film equations for flow and transport of multiple solutes, we compute dynamic results for tear film quantities of interest. We compare our computed intensity distributions with in vivo observations. A number of experimental features are recovered by the model.

preprint2015arXiv

A stable partitioned FSI algorithm for incompressible flow and deforming beams

An added-mass partitioned (AMP) algorithm is described for solving fluid-structure interaction (FSI) problems coupling incompressible flows with thin elastic structures undergoing finite deformations. The new AMP scheme is fully second-order accurate and stable, without sub-time-step iterations, even for very light structures when added-mass effects are strong. The fluid, governed by the incompressible Navier-Stokes equations, is solved in velocity-pressure form using a fractional-step method; large deformations are treated with a mixed Eulerian-Lagrangian approach on deforming composite grids. The motion of the thin structure is governed by a generalized Euler-Bernoulli beam model, and these equations are solved in a Lagrangian frame using two approaches, one based on finite differences and the other on finite elements. Special treatment of the AMP condition is required to couple the finite-element beam solver with the finite-difference-based fluid solver, and two coupling approaches are described. A normal-mode stability analysis is performed for a linearized model problem involving a beam separating two fluid domains, and it is shown that the AMP scheme is stable independent of the ratio of the mass of the fluid to that of the structure. A traditional partitioned (TP) scheme using a Dirichlet-Neumann coupling for the same model problem is shown to be unconditionally unstable if the added mass of the fluid is too large. A series of benchmark problems of increasing complexity are considered to illustrate the behavior of the AMP algorithm, and to compare the behavior with that of the TP scheme. The results of all these benchmark problems verify the stability and accuracy of the AMP scheme. Results for one benchmark problem modeling blood flow in a deforming artery are also compared with corresponding results available in the literature.

preprint2013arXiv

A stable FSI algorithm for light rigid bodies in compressible flow

In this article we describe a stable partitioned algorithm that overcomes the added mass instability arising in fluid-structure interactions of light rigid bodies and inviscid compressible flow. The new algorithm is stable even for bodies with zero mass and zero moments of inertia. The approach is based on a local characteristic projection of the force on the rigid body and is a natural extension of the recently developed algorithm for coupling compressible flow and deformable bodies. Normal mode analysis is used to prove the stability of the approximation for a one-dimensional model problem and numerical computations confirm these results. In multiple space dimensions the approach naturally reveals the form of the added mass tensors in the equations governing the motion of the rigid body. These tensors, which depend on certain surface integrals of the fluid impedance, couple the translational and angular velocities of the body. Numerical results in two space dimensions, based on the use of moving overlapping grids and adaptive mesh refinement, demonstrate the behavior and efficacy of the new scheme. These results include the simulation of the difficult problem of a shock impacting an ellipse of zero mass.

preprint2013arXiv

A stable partitioned FSI algorithm for incompressible flow and elastic solids

A stable partitioned algorithm for coupling incompressible flows with compressible elastic solids is described. This added-mass partitioned (AMP) scheme requires no sub-iterations, can be made fully second- or higher-order accurate, and remains stable even in the presence of strong added-mass effects. The approach extends the scheme of Banks et al. [1,2] for compressible flow, and uses Robin (mixed) boundary conditions with the fluid and solid solvers at the interface. The AMP Robin conditions are derived from a local characteristic decomposition in the solid at the interface. Two forms of the Robin conditions are derived depending on whether the fluid equations are advanced with a fractional-step method or not. A normal mode analysis for a discretization of an FSI model problem is performed to show that the new AMP algorithm is stable for any ratio of the solid and fluid densities, including the case of very light solids when the added-mass effects are large. In contrast, it is shown that a traditional partitioned algorithm involving a Dirichlet-Neumann coupling for the same FSI problem is formally unconditionally unstable for any ratio of densities. Exact traveling wave solutions are derived for three FSI model problems of increasing complexity, and these solutions are used to verify the stability and accuracy of the corresponding numerical results obtained from the AMP algorithm for the cases of light, medium and heavy solids.

preprint2013arXiv

A stable partitioned FSI algorithm for incompressible flow and structural shells

A stable partitioned algorithm for fluid-structure interaction (FSI) problems that couple viscous incompressible flow with structural shells or beams is described. This added-mass partitioned (AMP) scheme uses Robin (mixed) interface conditions for the pressure and velocity in the fluid that are derived directly from the governing equations. The AMP scheme is stable even for very light structures, requires no sub-iterations, and can be made fully second-order, or higher-order, accurate. The new scheme is shown to be stable through the analysis of a model problem. Exact traveling wave solutions for three FSI model problems are derived. Numerical results for a linearized FSI problem in two-dimensions, using these exact solutions, demonstrate that the scheme is stable and accurate, even for very light structures.