Fast maximum intensity projections of large medical data sets by exploiting hierarchical memory architectures.

Maximum intensity projections (MIPs) are an important visualization technique for angiographic data sets. Efficient data inspection requires frame rates of at least five frames per second at preserved image quality. Despite the advances in computer technology, this task remains a challenge. On the one hand, the sizes of computed tomography and magnetic resonance images are increasing rapidly. On the other hand, rendering algorithms do not automatically benefit from the advances in processor technology, especially for large data sets. This is due to the faster evolving processing power and the slower evolving memory access speed, which is bridged by hierarchical cache memory architectures. In this paper, we investigate memory access optimization methods and use them for generating MIPs on general-purpose central processing units (CPUs) and graphics processing units (GPUs), respectively. These methods can work on any level of the memory hierarchy, and we show that properly combined methods can optimize memory access on multiple levels of the hierarchy at the same time. We present performance measurements to compare different algorithm variants and illustrate the influence of the respective techniques. On current hardware, the efficient handling of the memory hierarchy for CPUs improves the rendering performance by a factor of 3 to 4. On GPUs, we observed that the effect is even larger, especially for large data sets. The methods can easily be adjusted to different hardware specifics, although their impact can vary considerably. They can also be used for other rendering techniques than MIPs, and their use for more general image processing task could be investigated in the future.     

Phosphoinositide 3-kinase induces scattering and tubulogenesis in epithelial cells through a novel pathway

Hepatocyte growth factor/scatter factor (HGF/SF) treatment of the Madin-Darby canine kidney epithelial cell line causes scattering of cells grown in monolayer culture and the formation of branching tubules by cells grown in collagen gels. HGF/SF causes prolonged activation of both the mitogen-activated protein (MAP) kinase extracellular signal-regulated kinase 2 (ERK2) and the phosphoinositide 3-OH kinase (PI 3-kinase) target protein kinase B (PKB)/Akt; inhibition of either the MAP kinase pathway by the MAP kinase/ERK kinase inhibitor PD98059 or the PI 3-kinase pathway by LY294002 blocks HGF/SF induction of scattering, although in morphologically distinct ways. Expression of constitutively activated PI 3-kinase, Ras, or R-Ras will cause scattering, but activated Raf will not, indicating that activation of the MAP kinase pathway is not sufficient for this response. Downstream of PI 3-kinase, activated PKB/Akt and Rac are both unable to induce scattering, implicating a novel pathway. Scattering induced by Ras or PI 3-kinase is sensitive to PD98059, as well as to LY294002, suggesting that basal MAP kinase activity is required, but not sufficient, for the scattering response. Induction of MDCK cell tubulogenesis in collagen gels by HGF/SF is inhibited by PD98059; expression of activated Ras and Raf causes disorganized growth in this system, but activated PI 3-kinase or R-Ras causes branching tubule formation similar to that seen with HGF/SF treatment. These data indicate that multiple signaling pathways acting downstream of Met and Ras are needed for these morphological effects; scattering is induced primarily by the PI 3-kinase pathway, which acts through effectors other than PKB/Akt or Rac and requires at least basal MAP kinase function. Elevated PI 3-kinase activity induces tubulogenesis, but total inhibition and excess activation of the MAP kinase pathway both oppose this effect.     

A new monitoring system for piping systems in fossil fired power plants

A CEC-funded project has been performed to tackle the problem of producing an advanced Life Monitoring System (LMS) which would calculate the creep and fatigue damage experienced by high temperature pipework components. Four areas were identified where existing Life Monitoring System technology could be improved:

1. 1. the inclusion of creep relaxation

2. 2. the inclusion of external loads on components

3. 3. a more accurate method of calculating thermal stresses due to temperature transients

4. 4. the inclusion of high cycle fatigue terms.

The creep relaxation problem was solved using stress reduction factors in an analytical in-elastic stress calculation. The stress reduction factors were produced for a number of common geometries and materials by means of non-linear finite element analysis. External loads were catered for by producing influence coefficients from in-elastic analysis of the particular piping system and using them to calculate bending moments at critical positions on the pipework from load and displacement measurements made at the convenient points at the pipework. The thermal stress problem was solved by producing a completely new solution based on Green's Function and Fast Fourier transforms. This allowed the thermal stress in a complex component to be calculated from simple non-intrusive thermocouple measurements made on the outside of the component. The high-cycle fatigue problem was dealt with precalculating the fatigue damage associated with standard transients and adding this damage to cumulative total when a transient occurred.

The site testing provided good practical experience and showed up problems which would not otherwise have been detected.     

Simulation of the submerged energy nozzle-mold water model system using laser-optical and computational fluid dynamics methods

The present work describes quantitative digital particle image velocimetry measurements of a full-scale water model of a thin slab mold. Different casting speeds and two submerged entry nozzles with one and two outlet ports have been investigated. The flow pattern of the single-port nozzle shows a counterclockwise-rotating double vortex that is nearly steady-state but leads to high stationary surface waves. The flow jets out of the two-port nozzle oscillate and produce a transient flow pattern with low wave amplitudes. The amplitudes for the one-port nozzle show a linear variation with the volumetric flow rate. The experimental results lead to a good interpretation of the flow phenomena and are used to validate steady-state numerical simulations with the commercial program, CFX, on the basis of the Reynolds equations. To describe anisotropic turbulence effects, the Reynolds stress model (RSM) is used for the flat single-port nozzle and the standard k-ɛ model for the mold flow. The calculated mean velocities and wave amplitudes, predicted from pressure distribution at the water surface, are generally in the consensus of the experimental data. An erratum to this article is available at .     

Three-dimensional MHD high-resolution computations with CWENO employing adaptive mesh refinement

Until recently, numerical simulations of discontinuities in highly super-Alfvénic plasmas have been severely limited by comparatively crude resolution and accuracy. Significant progress in the numerical simulation of such plasmas was achieved with the recently implemented Central Weighted Essentially Non-Oscillatory (CWENO) scheme. Combining this technique with that of adaptive mesh refinement (AMR), we have developed a third-order numerical scheme, which is able to efficiently capture strong gradients on spatial scales being small compared to the overall scale of the plasma system considered. Here, we first describe important algorithmic aspects of the scheme as well as the physics included in it. Second, we present the results of various performance tests. And, third, we illustrate its application to ‘real world problems’ using the example of the dynamics of a Sedov-type explosion.