Template-based quadrilateral meshing

Generating quadrilateral meshes is a highly non-trivial task, as design decisions are frequently driven by specific application demands. Automatic techniques can optimize objective quality metrics, such as mesh regularity, orthogonality, alignment and adaptivity; however, they cannot make subjective design decisions. There are a few quad meshing approaches that offer some mechanisms to include the user in the mesh generation process; however, these techniques either require a large amount of user interaction or do not provide necessary or easy to use inputs. Here, we propose a template-based approach for generating quad-only meshes from triangle surfaces. Our approach offers a flexible mechanism to allow external input, through the definition of alignment features that are respected during the mesh generation process. While allowing user inputs to support subjective design decisions, our approach also takes into account objective quality metrics to produce semi-regular, quad-only meshes that align well to desired surface features. Published by Elsevier Ltd.

Efficient solutions to robust, semi-implicit discretizations of the immersed boundary method

The immersed boundary method is a versatile tool for the investigation of flow-structure interaction. In a large number of applications, the immersed boundaries or structures are very stiff and strong tangential forces on these interfaces induce a well-known, severe time-step restriction for explicit discretizations. This excessive stability constraint can be removed with fully implicit or suitable semi-implicit schemes but at a seemingly prohibitive computational cost. While economical alternatives have been proposed recently for some special cases, there is a practical need for a computationally efficient approach that can be applied more broadly. In this context, we revisit a robust semi-implicit discretization introduced by Peskin in the late 1970s which has received renewed attention recently. This discretization, in which the spreading and interpolation operators are lagged. leads to a linear system of equations for the inter-face configuration at the future time, when the interfacial force is linear. However, this linear system is large and dense and thus it is challenging to streamline its solution. Moreover, while the same linear system or one of similar structure could potentially be used in Newton-type iterations, nonlinear and highly stiff immersed structures pose additional challenges to iterative methods. In this work, we address these problems and propose cost-effective computational strategies for solving Peskin`s lagged-operators type of discretization. We do this by first constructing a sufficiently accurate approximation to the system`s matrix and we obtain a rigorous estimate for this approximation. This matrix is expeditiously computed by using a combination of pre-calculated values and interpolation. The availability of a matrix allows for more efficient matrix-vector products and facilitates the design of effective iterative schemes. We propose efficient iterative approaches to deal with both linear and nonlinear interfacial forces and simple or complex immersed structures with tethered or untethered points. One of these iterative approaches employs a splitting in which we first solve a linear problem for the interfacial force and then we use a nonlinear iteration to find the interface configuration corresponding to this force. We demonstrate that the proposed approach is several orders of magnitude more efficient than the standard explicit method. In addition to considering the standard elliptical drop test case, we show both the robustness and efficacy of the proposed methodology with a 2D model of a heart valve. (C) 2009 Elsevier Inc. All rights reserved.

Spaces of compact operators on C(2(m) circle plus [0, alpha]) spaces

We classify up to isomorphism the spaces of compact operators K(E, F), where E and F are Banach spaces of all continuous functions defined on the compact spaces 2(m) circle plus [0, alpha], the topological sum of Cantor cubes 2(m) and the intervals of ordinal numbers [0, alpha]. More precisely, we prove that if 2(m) and aleph(gamma) are not real-valued measurable cardinals and n >= aleph(0) is not sequential cardinal, then for every ordinals xi, eta, lambda and mu with xi >= omega(1), eta >= omega(1), lambda = mu < omega or lambda, mu is an element of [omega(gamma), omega(gamma+1)[, the following statements are equivalent: (a) K(C(2(m) circle plus [0, lambda]), C(2(n) circle plus [0, xi])) and K(C(2(m) circle plus [0, mu]), C(2(n) circle plus [0, eta]) are isomorphic. (b) Either C([0, xi]) is isomorphic to C([0, eta] or C([0, xi]) is isomorphic to C([0, alpha p]) and C([0, eta]) is isomorphic to C([0,alpha q]) for some regular cardinal alpha and finite ordinals p not equal q. Thus, it is relatively consistent with ZFC that this result furnishes a complete isomorphic classification of these spaces of compact operators. (C) 2010 Elsevier Inc. All rights reserved.