A Statistical Shape Model of the Morphological Variation of the Infrarenal Abdominal Aortic Aneurysm Neck

Hostile aortic neck characteristics, such as short length and large diameter, have been associated with type Ia endoleaks and reintervention after endovascular aneurysm repair (EVAR). However, such characteristics partially describe the complex aortic neck morphology. A more comprehensive quantitative description of 3D neck shape might lead to new insights into the relationship between aortic neck morphology and EVAR outcomes in individual patients. This study identifies the 3D morphological shape components that describe the infrarenal aortic neck through a statistical shape model (SSM). Pre-EVAR CT scans of 97 patients were used to develop the SSM. Parameterization of the morphology was based on the center lumen line reconstruction, a triangular surface mesh of the aortic lumen, 3D coordinates of the renal arteries, and the distal end of the aortic neck. A principal component analysis of the parametrization of the aortic neck coordinates was used as input for the SSM. The SSM consisted of 96 principal components (PCs) that each described a unique shape feature. The first five PCs represented 95% of the total morphological variation in the dataset. The SSM is an objective model that provides a quantitative description of the neck morphology of an individual patient

Axially symmetric membranes with polar tethers

Axially symmetric equilibrium configurations of the conformally invariant Willmore energy are shown to satisfy an equation that is two orders lower in derivatives of the embedding functions than the equilibrium shape equation, not one as would be expected on the basis of axial symmetry. Modulo a translation along the axis, this equation involves a single free parameter c.If c\ne 0, a geometry with spherical topology will possess curvature singularities at its poles. The physical origin of the singularity is identified by examining the Noether charge associated with the translational invariance of the energy; it is consistent with an external axial force acting at the poles. A one-parameter family of exact solutions displaying a discocyte to stomatocyte transition is described.Comment: 13 pages, extended and revised version of Non-local sine-Gordon equation for the shape of axi-symmetric membrane

Multi-scale 3-D Surface Description: Open and Closed Surfaces

A novel technique for multi-scale smoothing of a free-form 3-D surface is presented. Complete triangulated models of 3-D objects are constructed automatically and using a local parametrization technique, are then smoothed using a 2-D Gaussian filter. Our method for local parametrization makes use of semigeodesic coordinates as a natural and efficient way of sampling the local surface shape. The smoothing eliminates the surface noise together with high curvature regions such as sharp edges, therefore, sharp corners become rounded as the object is smoothed iteratively. Our technique for free-form 3-D multi-scale surface smoothing is independent of the underlying triangulation. It is also argued that the proposed technique is preferrable to volumetric smoothing or level set methods since it is applicable to incomplete surface data which occurs during occlusion. Our technique was applied to closed as well as open 3-D surfaces and the results are presented here

Multi-Scale Free-Form Surface Description and Curvature Estimation

A novel technique for multi-scale smoothing of a free-form 3-D surface is presented. Complete triangulated models of 3-D objects are constructed at our center [4] and using a local parametrization technique, are then smoothed using a 2-D Gaussian filter. Our method for local parametrization makes use of semigeodesic coordinates as a natural and efficient way of sampling the local surface shape. The smoothing eliminates the surface noise together with high curvature regions such as sharp edges, therefore, sharp corners become rounded as the object is smoothed iteratively. Our technique for free-form 3-D multi-scale surface smoothing is independent of the underlying triangulation. It is also argued that the proposed technique is preferrable to volumetric smoothing or level set methods since it is applicable to incomplete surface data which occurs during occlusion. The technique was applied to simple and complex 3-D objects and the results are presented here

Arbitrary Lagrangian–Eulerian finite element method for curved and deforming surfaces: I. General theory and application to fluid interfaces

An arbitrary Lagrangian–Eulerian (ALE) finite element method for arbitrarily curved and deforming two-dimensional materials and interfaces is presented here. An ALE theory is developed by endowing the surface with a mesh whose in-plane velocity need not depend on the in-plane material velocity, and can be specified arbitrarily. A finite element implementation of the theory is formulated and applied to curved and deforming surfaces with in-plane incompressible flows. Numerical inf–sup instabilities associated with in-plane incompressibility are removed by locally projecting the surface tension onto a discontinuous space of piecewise linear functions. The general isoparametric finite element method, based on an arbitrary surface parametrization with curvilinear coordinates, is tested and validated against several numerical benchmarks. A new physical insight is obtained by applying the ALE developments to cylindrical fluid films, which are computationally and analytically found to be stable to non-axisymmetric perturbations, and unstable with respect to long-wavelength axisymmetric perturbations when their length exceeds their circumference. A Lagrangian scheme is attained as a special case of the ALE formulation. Though unable to model fluid films with sustained shear flows, the Lagrangian scheme is validated by reproducing the cylindrical instability. However, relative to the ALE results, the Lagrangian simulations are found to have spatially unresolved regions with few nodes, and thus larger errors

A many-body interatomic potential for ionic systems: application to MgO

An analytic representation of the short-range repulsion energy in ionic systems is described that allows for the fact that ions may change their size and shape depending on their environment. This function is extremely efficient to evaluate relative to previous methods of modeling the same physical effects. Using a well-defined parametrization procedure we have obtained parameter sets for this energy function that reproduce closely the density functional theory potential energy surface of bulk MgO. We show how excellent agreement can be obtained with experimental measurements of phonon frequencies and temperature and pressure dependences of the density by using this effective potential in conjunction with ab initio parametrization.Comment: To appear in Journal of Chemical Physics (Oct 15th 2003