A Fast Algorithm to Solve Nonlinear Hypersingular Integral Equations Arising in a Crack Problem

A fast algorithm related to the generalized minimal residual algorithm (GMRES) is proposed to approximate solution of a nonlinear hypersingular integral equation arising in a crack problem. At first, a collocation method is proposed and developed in weighted Sobolev space. Then, the Newton-Kantorovjch method is used for solving the obtained system of nonlinear equations.

Endothelial shear stress from large-scale blood flow simulations

We discuss the optimal evaluation of endothelial shear stress for real-life case studies based on anatomic data acquisition. The fluid dynamic simulations require smoothing of the geometric dataset to avoid major artefacts in the flow patterns, especially in the proximity of bifurcations. A systematic series of simulations at different corrugation levels shows that, below a smoothing length of about 0.5 mm, the numerical data are insensitive to further smoothing. © 2011 The Royal Society.

Multiscale hemodynamics using GPU clusters

The parallel implementation of MUPHY, a concurrent multiscale code for large-scale hemodynamic simulations in anatomically realistic geometries, for multi-GPU platforms is presented. Performance tests show excellent results, with a nearly linear parallel speed-up on up to 32GPUs and a more than tenfold GPU/CPU acceleration, all across the range of GPUs. The basic MUPHY scheme combines a hydrokinetic (Lattice Boltzmann) representation of the blood plasma, with a Particle Dynamics treatment of suspended biological bodies, such as red blood cells.

Petaflop biofluidics simulations on a two million-core system

We present a computational framework for multi-scale simulations of real-life biofluidic problems. The framework allows to simulate suspensions composed by hundreds of millions of bodies interacting with each other and with a surrounding fluid in complex geometries. We apply the methodology to the simulation of blood flow through the human coronary arteries with a spatial resolution comparable with the size of red blood cells, and physiological levels of hematocrit (the red blood cell volume fraction).

Quantitative morphotectonics of the Pliocene to Quaternary Auletta basin, southern Italy

The geomorphological evolution of the Pliocene-Quaternary Auletta basin, a wide fault-bounded depression of the southern Apennines axial zone, Italy, was reconstructed using both DEM-based morphometric analysis and classical morphotectonic investigations. Morphotectonic analyses have been integrated with geological, structural and paleomagnetic data in order to reconstruct the Quaternary evolution of the area. The Auletta basin coincides with the lower valley of the Tanagro River and is filled by Pliocene to Pleistocene marine and continental sediments.

Runge-Kutta Discretizations of Infinite Horizon Optimal Control Problems with Steady-State Invariance

Direct numerical approximation of a continuous-time infinite horizon control problem, requires to recast the model as a discrete-time, finite-horizon control model. The quality of the optimization results can be heavily degraded if the discretization process does not take into account features of the original model to be preserved.

Three-Dimensional Lattice Pseudo-Potentials for Multiphase Flow Simulations at High Density Ratios

It is shown that the combination of generalized Van der Waals equations of state with high-order discrete velocity lattices, permits to simulate the dynamics of liquid droplets at air-water density ratios, with very moderate levels of spurious currents near the droplet interface. Satisfactory agreement with experimental data on droplet collisions at density ratios of order thousand is reported.