Each year, approximately 250,000 surgical procedures are performed to repair or to replace cardiac valves [1], and of these, approximately 50,000 are aortic valve replacement operations. Despite decades of development, many of the limitations of cardiac valve prostheses remain consequences of the fluid dynamics generated by the replacement valve [1]. To enable cross-validation studies that allow for detailed comparisons of experimental and computational results, we are developing fluid-structure interaction (FSI) models of the dynamics of aortic valve prostheses mounted in a ViVitro Systems, Inc. pulse duplicator apparatus. Such an experimentally validated computational model promises both to facilitate the design of novel valve prostheses, and also to assist in the regulatory process, by providing access to detailed spatially- and temporally-resolved flow data that are challenging to obtain via experimental approaches.Our numerical approach to FSI is based on the immersed boundary (IB) method [2]. The structural dynamics are described in Lagrangian form using a material coordinate system, whereas the momentum of the fluid-solid system and the viscosity and incompressibility of the fluid are described in Eulerian form using fixed Cartesian physical coordinates. Lagrangian and Eulerian variables are coupled by integral transforms with Dirac delta function kernels. When the equations are discretized, the singular delta function is replaced by a regularized version of the delta function. See Peskin [2] for details.To treat rigid body dynamics within the framework of the IB method, we adopt an approach similar to that recently developed by Kim and Peskin [3], in which the continuous equations are:(1)ρ(∂u∂t(x,t)+u(x,t)·∇u(x,t)) =-∇p(x,t)+μ∇2u(x,t)+f(x,t)(2)∇·u(x,t)=0(3)f(x,t)=∫UF(s,t)δ(x-X(s,t))ds(4)∂X∂t(s,t)=∫Ωu(x,t)δ(x-X(s,t))dx(5)∂Y∂t(s,t)=V(t)+W(t)×R(s,t)(6)F(s,t)=κ(Y(s,t)-X(s,t))in which x∈Ω are physical coordinates, s∈U are material coordinates, X(s,t) and Y(s,t) are time-dependent mappings from material coordinates to current coordinates, u(x,t) is the Eulerian velocity field, p(x,t) is the Eulerian pressure field, f(x,t) and F(s,t) are equivalent Eulerian and Lagrangian force densities, ρ is the fluid density, μ is the fluid viscosity, and δ(x)=δ(x)δ(y)δ(z) is the three-dimensional Dirac delta function. The force field F(s,t) acts to impose the rigidity constraint ∂X∂t(s,t)=U(s,t)=V(t)+W(t)×R(s,t)in which V(t) and W(t) are respectively the linear and angular velocity of the structure and R(s,t) is the radius vector. In the limit κ→∞, the constraint is imposed exactly; for finite κ>0, the constraint is only approximately satisfied. The dynamics of V(t) and W(t) are determined from the requirement that the time rate of change of the “excess” linear and angular momentum (i.e., in excess of the momentum accounted for by the momentum equation (1)) be proportional to the net force ∫Ω-F(s,t)ds and net torque ∫Ω-F(s,t)×R(s,t)ds acting on the body.Geometrical models of the aortic section of the ViVitro pulse duplicator and of a St. Jude Regent valve were created using SolidWorks (Dassault Systèmes SolidWorks Corp., Waltham, MA). These geometrical models were meshed using CUBIT (Sandia National Laboratory, Albuquerque, NM). To drive flow through the valve, we impose a physiological left ventricular pressure waveform, and characteristic downstream compliance and resistance to mimic the response of the systemic circulation. Simulations used the IBAMR software [4]. Initial simulation results are shown in Fig. 1, demonstrating that realistic flow rates can be obtained using this model under realistic driving and loading conditions. We are presently fine tuning this model in preparation for validation studies.We acknowledge discussions with A. P. S. Bhalla, A. Donev, D. M. McQueen, and C. S. Peskin, and funding from AHA award 10SDG4320049, NIH award HL117063, and NSF awards DMS 1016554 and OCI 1047734.