Collin E. Haese

L. River Spencer, Reagan A. Cardoza, Vijay K. Dubey et al.

Ultrasound imaging tasks such as calibration, inverse parameter estimation, and acquisition design require models that are physically grounded, efficient, and differentiable with respect to meaningful material and system parameters. While full-wave solvers offer high fidelity, they are often too expensive for iterative optimization, and existing ray-based methods have mostly been limited to forward simulation. In this work, we present a fully differentiable end-to-end ultrasound simulation framework based on full-path Monte Carlo ray tracing. Building on UltraRay, the method propagates gradients from image-space losses back through acoustic transport, beamforming, and post-processing, enabling gradient-based optimization over scene and acquisition parameters. The framework combines differentiable ray transport in Mitsuba 3/Dr.Jit with a custom differentiable bridge through the ultrasound image-formation pipeline. Forward examples reproduce expected geometric image features and capture more complex anatomical structures. In inverse problems, the method recovers known parameters in a simulated-reference setting and identifies effective parameters that improve agreement between simulated and experimental B-mode images in a simulation-to-real setting. Finite-difference comparisons further support the consistency of the computed gradients. Overall, this work provides a practical foundation for differentiable, physics-based ultrasound simulation and optimization.

Vijay K. Dubey, Collin E. Haese, Osman Gültekin et al.

Surrogate models for the rapid inference of nonlinear boundary value problems in mechanics are helpful in a broad range of engineering applications. However, effective surrogate modeling of applications involving the contact of deformable bodies, especially in the context of varying geometries, is still an open issue. In particular, existing methods are confined to rigid body contact or, at best, contact between rigid and soft objects with well-defined contact planes. Furthermore, they employ contact or collision detection filters that serve as a rapid test but use only the necessary and not sufficient conditions for detection. In this work, we present a graph neural network architecture that utilizes continuous collision detection and, for the first time, incorporates sufficient conditions designed for contact between soft deformable bodies. We test its performance on two benchmarks, including a problem in soft tissue mechanics of predicting the closed state of a bioprosthetic aortic valve. We find a regularizing effect on adding additional contact terms to the loss function, leading to better generalization of the network. These benefits hold for simple contact at similar planes and element normal angles, and complex contact at differing planes and element normal angles. We also demonstrate that the framework can handle varying reference geometries. However, such benefits come with high computational costs during training, resulting in a trade-off that may not always be favorable. We quantify the training cost and the resulting inference speedups on various hardware architectures. Importantly, our graph neural network implementation results in up to a thousand-fold speedup for our benchmark problems at inference.