Immersed Boundary Method for generalised finite volume and finite difference Navier-Stokes solvers

In Immersed Boundary Methods (IBM) the effect of complex geometries is introduced through the forces added in the Navier-Stokes solver at the grid points in the vicinity of the immersed boundaries. Most of the methods in the literature have been used with Cartesian grids. Moreover many of the methods developed in the literature do not satisfy some basic conservation properties (the conservation of torque, for instance) on non-uniform meshes. In this paper we will follow the RKPM method originated by Liu et al. [1] to build locally regularized functions that verify a number of integral conditions. These local approximants will be used both for interpolating the velocity field and for spreading the singular force field in the framework of a pressure correction scheme for the incompressible Navier-Stokes equations. We will also demonstrate the robustness and effectiveness of the scheme through various examples. Copyright © 2010 by ASME.

Shape vs. volume: Invariant shape descriptors for 3D region of interest characterization in MRI

We propose a new approach for quantifying regions of interest (ROIs) in medical image data. Rotationally invariant shape descriptors (ISDs) were applied to 3D brain regions extracted from MRI scans of 5 Parkinson's patients and 10 control subjects. We concentrated on the thalamus and the caudate nucleus since prior studies have suggested they are affected in Parkinson's disease (PD). In the caudate, both the ISD and volumetric analyses found significant differences between control and PD subjects. The ISD analysis however revealed additional differences between the left and right caudate nuclei in both control and PD subjects. In the thalamus, the volumetric analysis showed significant differences between PD and control subjects, while ISD analysis found significant differences between the left and right thalami in control subjects but not in PD patients, implying disease-induced shape changes. These results suggest that employing ISDs for ROI characterization both complements and extends traditional volumetric analyses. © 2006 IEEE.

On the transition from finite-volume negatively buoyant releases to continuous fountains

An experimental investigation to identify the source conditions that distinguish finite-volume negatively buoyant fluid projectile behaviour from fountain behaviour in quiescent environments of uniform density is described. Finite-volume releases are governed by their source Froude number Fr D and the aspect ratio L/D of the release, where L denotes the length of the column of fluid dispensed vertically from the nozzle of diameter D. We establish the influence of L/D on the peak rise heights of a release formed by dispensing saline solution into fresh water for 0volume or continuous flux. The critical aspect ratio (L/D) f, for a given Fr D, which when exceeded no longer influenced release behaviour, led to the determination of Fr D, (L/D) f paired source conditions that give rise to solely Froude-number-dependent, i.e. fountain-like, behaviour. As such, we make the link between finite-volume releases and continuous fountains. The Fr D(L/D) f pairs led us directly to the classification of a Fr D, L/D space from which source conditions giving rise to either negatively buoyant projectiles or fountains may be readily identified. The variation of (L/D) f with Fr D corresponds closely to established fountain regimes of very weak, weak and forced fountains. Moreover, our results indicate that the formation or otherwise of a primary vortex, as fluid is ejected, has a profound influence on the length of the dispensed fluid column that is necessary to achieve rise heights equal to fountain rise heights. © 2012 Cambridge University Press.

Semi-implicit second order finite volume jet stream computations

Semi-implicit, second order temporal and spatial finite volume computations of the flow in a differentially heated rotating annulus are presented. For the regime considered, three cyclones and anticyclones separated by a relatively fast moving jet of fluid or "jet stream" are predicted. Two second order methods are compared with, first order spatial predictions, and experimental measurements. Velocity vector plots are used to illustrate the predicted flow structure. Computations made using second order central differences are shown to agree best with experimental measurements, and to be stable for integrations over long time periods (> 1000s). No periodic smoothing is required to prevent divergence.

Model reductions for inference: generality of pairwise, binary, and planar factor graphs.

We offer a solution to the problem of efficiently translating algorithms between different types of discrete statistical model. We investigate the expressive power of three classes of model-those with binary variables, with pairwise factors, and with planar topology-as well as their four intersections. We formalize a notion of "simple reduction" for the problem of inferring marginal probabilities and consider whether it is possible to "simply reduce" marginal inference from general discrete factor graphs to factor graphs in each of these seven subclasses. We characterize the reducibility of each class, showing in particular that the class of binary pairwise factor graphs is able to simply reduce only positive models. We also exhibit a continuous "spectral reduction" based on polynomial interpolation, which overcomes this limitation. Experiments assess the performance of standard approximate inference algorithms on the outputs of our reductions.

Pulsed-field magnetization of drilled bulk high-temperature superconductors: flux front propagation in the volume and on the surface

We present a method for characterizing the propagation of the magnetic flux in an artificially drilled bulk high-temperature superconductor (HTS) during a pulsed-field magnetization. As the magnetic pulse penetrates the cylindrical sample, the magnetic flux density is measured simultaneously in 16 holes by means of microcoils that are placed across the median plane, i.e. at an equal distance from the top and bottom surfaces, and close to the surface of the sample. We discuss the time evolution of the magnetic flux density in the holes during a pulse and measure the time taken by the external magnetic flux to reach each hole. Our data show that the flux front moves faster in the median plane than on the surface when penetrating the sample edge; it then proceeds faster along the surface than in the bulk as it penetrates the sample further. Once the pulse is over, the trapped flux density inside the central hole is found to be about twice as large in the median plane than on the surface. This ratio is confirmed by modelling.

Study by Hall probe mapping of the trapped flux modification produced by local heating in YBCO HTS bulks for different surface/volume ratios

The aim of this report is to compare the trapped field distribution under a local heating created at the sample edge for different sample morphologies. Hall probe mappings of the magnetic induction trapped in YBCO bulk samples maintained out of thermal equilibrium were performed on YBCO bulk single domains, YBCO single domains with regularly spaced hole arrays, and YBCO superconducting foams. The capability of heat draining was quantified by two criteria: the average induction decay and the size of the thermally affected zone caused by a local heating of the sample. Among the three investigated sample shapes, the drilled single domain displays a trapped induction which is weakly affected by the local heating while displaying a high trapped field. Finally, a simple numerical modelling of the heat flux spreading into a drilled sample is used to suggest some design rules about the hole configuration and their size. © 2005 IOP Publishing Ltd.