Source author record

Nader Saffari

Nader Saffari appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.med-ph math.NA Numerical Analysis physics.comp-ph physics.flu-dyn

Catalog footprint

What is connected

3works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

preprint2022arXiv

Frequency-robust preconditioning of boundary integral equations for acoustic transmission

The scattering and transmission of harmonic acoustic waves at a penetrable material are commonly modelled by a set of Helmholtz equations. This system of partial differential equations can be rewritten into boundary integral equations defined at the surface of the objects and solved with the boundary element method (BEM). High frequencies or geometrical details require a fine surface mesh, which increases the number of degrees of freedom in the weak formulation. Then, matrix compression techniques need to be combined with iterative linear solvers to limit the computational footprint. Moreover, the convergence of the iterative linear solvers often depends on the frequency of the wave field and the objects' characteristic size. Here, the robust PMCHWT formulation is used to solve the acoustic transmission problem. An operator preconditioner based on on-surface radiation conditions (OSRC) is designed that yields frequency-robust convergence characteristics. Computational benchmarks compare the performance of this novel preconditioned formulation with other preconditioners and boundary integral formulations. The OSRC preconditioned PMCHWT formulation effectively simulates large-scale problems of engineering interest, such as focused ultrasound treatment of osteoid osteoma.

preprint2015arXiv

A numerical model for the study of photoacoustic imaging of brain tumours

Photoacoustic imaging has shown great promise for medical imaging, where optical energy absorption by blood haemoglobin is used as the contrast mechanism. A numerical method was developed for the in-silico assessment of the photoacoustic image reconstruction of the brain. Image segmentation techniques were used to prepare a digital phantom from MR images. Light transport through brain tissue was modelled using a Finite Element approach. The resulting acoustic pressure was then estimated by pulsed photoacoustics considerations. The forward acoustic wave propagation was modelled by the linearized coupled first order wave equations and solved by an acoustic k-space method. Since skull bone is an elastic solid and strongly attenuates ultrasound (due to both scattering and absorption), a k-space method was developed for elastic media. To model scattering effects, a new approach was applied based on propagation in random media. In addition, absorption effects were incorporated using a power law. Finally, the acoustic pressure was reconstructed using the k-space time reversal technique. The simulations were ran in 3D to produce the photoacoustic tomogram of a brain tumour. The results demonstrate the convergence of the models, and their suitability for investigating the photoacoustic imaging process.

preprint2009arXiv

The effects of nonlinear wave propagation on the stability of inertial cavitation

In the context of forecasting temperature and pressure fields in high-intensity focussed ultrasound, the accuracy of predictive models is critical for the safety and efficacy of treatment. In such fields inertial cavitation is often observed. Classically, estimations of cavitation thresholds have been based on the assumption that the incident wave at the surface of a bubble was the same as in the far-field, neglecting the effect of nonlinear wave propagation. By modelling the incident wave as a solution to Burgers' equation using weak shock theory, the effects of nonlinear wave propagation on inertial cavitation are investigated using both numerical and analytical techniques. From radius-time curves for a single bubble, it is observed that there is a reduction in the maximum size of a bubble undergoing inertial cavitation and that the inertial collapse occurs earlier in contrast with the classical case. Corresponding stability thresholds for a bubble whose initial radius is slightly below the critical Blake radius are calculated. Bifurcation diagrams and frequency-response curves are presented associated with the loss of stability. The consequences and physical implications of the results are discussed with respect to the classical results.