We present a nodal finite-element method that can be used to compute in parallel highly accurate solutions for 3-D controlled-source electromagnetic forward-modelling problems in anisotropic media. Secondary coupled-potential formulation of Maxwell's equations allows to avoid the singularities introduced by the sources, while completely unstructured tetrahedral meshes and mesh refinement support an accurate representation of geological and bathymetric complexity and improve the solution accuracy. Different complex iterative solvers and an efficient pre-conditioner based on the sparse approximate inverse are used for solving the resulting large sparse linear system of equations. Results are compared with the ones of other researchers to check the accuracy of the method. We demonstrate the performance of the code in large problems with tens and even hundreds of millions of degrees of freedom. Scalability tests on massively parallel computers show that our code is highly scalable.

The dynamics of unsteady flow in the human large airways during a rapid inhalation were investigated using highly detailed large-scale computational fluid dynamics on a subject-specific geometry. The simulations were performed to resolve all the spatial and temporal scales of the flow, thanks to the use of massive computational resources. A highly parallel finite element code was used, running on two supercomputers, solving the transient incompressible Navier-Stokes equations on unstructured meshes. Given that the finest mesh contained 350 million elements, the study sets a precedent for large-scale simulations of the respiratory system, proposing an analysis strategy for mean flow, fluctuations and wall shear stresses on a rapid and short inhalation (a so-called sniff). The geometry used encompasses the exterior face and the airways from the nasal cavity, through the trachea and up to the third lung bifurcation; it was derived from a contrast-enhanced computed tomography (CT) scan of a 48-year-old male. The transient inflow produces complex flows over a wide range of Reynolds numbers (Re). Thanks to the high fidelity simulations, many features involving the flow transition were observed, with the level of turbulence clearly higher in the throat than in the nose. Spectral analysis revealed turbulent characteristics persisting downstream of the glottis, and were captured even with a medium mesh resolution. However a fine mesh resolution was found necessary in the nasal cavity to observe transitional features. This work indicates the potential of large-scale simulations to further understanding of airway physiological mechanics, which is essential to guide clinical diagnosis; better understanding of the flow also has implications for the design of interventions such as aerosol drug delivery.

The human heart beats as a result of multiscale nonlinear dynamics coupling subcellular to whole organ processes, achieving electrophysiologically-driven mechanical contraction. Computational cardiac modelling and simulation have achieved a great degree of maturity, both in terms of mathematical models of underlying biophysical processes and the development of simulation software.

Background It is thought that aneurysms evolve as the result of progressive degradation of the wall in response to abnormal hemodynamics characterized by either high or low wall shear stress (WSS). Objective To investigate the effects of these two different hemodynamic pathways in a series of cerebral aneurysms with known rupture sites. Methods Nine aneurysms in which the rupture site could be identified in three-dimensional images were analyzed. The WSS distribution was obtained from computational fluid dynamics (CFD) simulations. Internal wall stresses were computed using structural wall models under hemodynamic loads determined by the CFD models. Wall properties (thickness and stiffness) were modulated with the WSS distribution (increased or decreased in regions of high or low WSS) to test possible wall degradation pathways. Rupture probability indices (RPI) were calculated to compare different wall models. Results Most rupture sites aligned with the intrasaccular flow stream and downstream of the primary impaction zone. The model that best explained the rupture site (produced higher RPI) in eight of the nine aneurysms (89%) had thinner and stiffer walls in regions of abnormally high WSS. The remaining case (11%) was best explained by a model with thinner and stiffer walls in regions of abnormally low WSS. Conclusions Aneurysm rupture seems to be caused by localized degradation and weakening of the wall in response to abnormal hemodynamics. Image-based computational models assuming wall thinning and stiffening in regions of abnormally high WSS were able to explain most of the observed rupture sites.