Computational Models of Material Interfaces for the Study of Extracorporeal Shock Wave Therapy

Extracorporeal Shock Wave Therapy (ESWT) is a noninvasive treatment for a variety of musculoskeletal ailments. A shock wave is generated in water and then focused using an acoustic lens or reflector so the energy of the wave is concentrated in a small treatment region where mechanical stimulation enhances healing. In this work we have computationally investigated shock wave propagation in ESWT by solving a Lagrangian form of the isentropic Euler equations in the fluid and linear elasticity in the bone using high-resolution finite volume methods. We solve a full three-dimensional system of equations and use adaptive mesh refinement to concentrate grid cells near the propagating shock. We can model complex bone geometries, the reflection and mode conversion at interfaces, and the the propagation of the resulting shear stresses generated within the bone. We discuss the validity of our simplified model and present results validating this approach

A simple and efficient BEM implementation of quasistatic linear visco-elasticity

A simple, yet efficient procedure to solve quasistatic problems of special linear visco-elastic solids at small strains with equal rheological response in all tensorial components, utilizing boundary element method (BEM), is introduced. This procedure is based on the implicit discretisation in time (the so-called Rothe method) combined with a simple "algebraic" transformation of variables, leading to a numerically stable procedure (proved explicitly by discrete energy estimates), which can be easily implemented in a BEM code to solve initial-boundary value visco-elastic problems by using the Kelvin elastostatic fundamental solution only. It is worth mentioning that no inverse Laplace transform is required here. The formulation is straightforward for both 2D and 3D problems involving unilateral frictionless contact. Although the focus is to the simplest Kelvin-Voigt rheology, a generalization to Maxwell, Boltzmann, Jeffreys, and Burgers rheologies is proposed, discussed, and implemented in the BEM code too. A few 2D and 3D initial-boundary value problems, one of them with unilateral frictionless contact, are solved numerically

Three-dimensional CFD simulations with large displacement of the geometries using a connectivity-change moving mesh approach

This paper deals with three-dimensional (3D) numerical simulations involving 3D moving geometries with large displacements on unstructured meshes. Such simulations are of great value to industry, but remain very time-consuming. A robust moving mesh algorithm coupling an elasticity-like mesh deformation solution and mesh optimizations was proposed in previous works, which removes the need for global remeshing when performing large displacements. The optimizations, and in particular generalized edge/face swapping, preserve the initial quality of the mesh throughout the simulation. We propose to integrate an Arbitrary Lagrangian Eulerian compressible flow solver into this process to demonstrate its capabilities in a full CFD computation context. This solver relies on a local enforcement of the discrete geometric conservation law to preserve the order of accuracy of the time integration. The displacement of the geometries is either imposed, or driven by fluid–structure interaction (FSI). In the latter case, the six degrees of freedom approach for rigid bodies is considered. Finally, several 3D imposed-motion and FSI examples are given to validate the proposed approach, both in academic and industrial configurations

Extended Smoothed Boundary Method for Solving Partial Differential Equations with General Boundary Conditions on Complex Boundaries

In this article, we describe an approach for solving partial differential equations with general boundary conditions imposed on arbitrarily shaped boundaries. A continuous function, the domain parameter, is used to modify the original differential equations such that the equations are solved in the region where a domain parameter takes a specified value while boundary conditions are imposed on the region where the value of the domain parameter varies smoothly across a short distance. The mathematical derivations are straightforward and generically applicable to a wide variety of partial differential equations. To demonstrate the general applicability of the approach, we provide four examples herein: (1) the diffusion equation with both Neumann and Dirichlet boundary conditions; (2) the diffusion equation with both surface diffusion and reaction; (3) the mechanical equilibrium equation; and (4) the equation for phase transformation with the presence of additional boundaries. The solutions for several of these cases are validated against corresponding analytical and semi-analytical solutions. The potential of the approach is demonstrated with five applications: surface-reaction-diffusion kinetics with a complex geometry, Kirkendall-effect-induced deformation, thermal stress in a complex geometry, phase transformations affected by substrate surfaces, and a self-propelled droplet.Comment: This document is the revised version of arXiv:0912.1288v

Smoothed Boundary Method for Solving Partial Differential Equations with General Boundary Conditions on Complex Boundaries

In this article, we describe an approach for solving partial differential equations with general boundary conditions imposed on arbitrarily shaped boundaries. A function that has a prescribed value on the domain in which a differential equation is valid and smoothly but rapidly varying values on the boundary where boundary conditions are imposed is used to modify the original differential equations. The mathematical derivations are straight forward, and generically applicable to a wide variety of partial differential equations. To demonstrate the general applicability of the approach, we provide four examples: (1) the diffusion equation with both Neumann and Dirichlet boundary conditions, (2) the diffusion equation with surface diffusion, (3) the mechanical equilibrium equation, and (4) the equation for phase transformation with additional boundaries. The solutions for a few of these cases are validated against corresponding analytical and semi-analytical solutions. The potential of the approach is demonstrated with five applications: surface-reaction diffusion kinetics with a complex geometry, Kirkendall-effect-induced deformation, thermal stress in a complex geometry, phase transformations affected by substrate surfaces, and a self-propelling droplet.Comment: A better smooth algorithm has been developed and tested, will soon replace Eq. 58 in page 16. We have also developed a level-set moving boundary SBM method, and it will replace the Navier-Stokes-Cahn-Hilliard type domain parameter tracking method in Section 5.