Multi-resolution isotropic strain limiting

In this paper we describe a fast strain-limiting method that allows stiff, incompliant materials to be simulated efficiently. Unlike prior approaches, which act on springs or individual strain components, this method acts on the strain tensors in a coordinate-invariant fashion allowing isotropic behavior. Our method applies to both two-and three-dimensional strains, and only requires computing the singular value decomposition of the deformation gradient, either a small 2x2 or 3x3 matrix, for each element. We demonstrate its use with triangular and tetrahedral linear-basis elements. For triangulated surfaces in three-dimensional space, we also describe a complementary edge-angle-limiting method to limit out-of-plane bending. All of the limits are enforced through an iterative, non-linear, Gauss-Seidel-like constraint procedure. To accelerate convergence, we propose a novel multi-resolution algorithm that enforces fitted limits at each level of a non-conforming hierarchy. Compared with other constraint-based techniques, our isotropic multi-resolution strain-limiting method is straightforward to implement, efficient to use, and applicable to a wide range of shell and solid materials. © 2010 ACM

On the non-local geometry of turbulence

A multi-scale methodology for the study of the non-local geometry of eddy structures in turbulence is developed. Starting from a given three-dimensional field, this consists of three main steps: extraction, characterization and classification of structures. The extraction step is done in two stages. First, a multi-scale decomposition based on the curvelet transform is applied to the full three-dimensional field, resulting in a finite set of component three-dimensional fields, one per scale. Second, by iso-contouring each component field at one or more iso-contour levels, a set of closed iso-surfaces is obtained that represents the structures at that scale. The characterization stage is based on the joint probability density function (p.d.f.), in terms of area coverage on each individual iso-surface, of two differential-geometry properties, the shape index and curvedness, plus the stretching parameter, a dimensionless global invariant of the surface. Taken together, this defines the geometrical signature of the iso-surface. The classification step is based on the construction of a finite set of parameters, obtained from algebraic functions of moments of the joint p.d.f. of each structure, that specify its location as a point in a multi-dimensional ‘feature space’. At each scale the set of points in feature space represents all structures at that scale, for the specified iso-contour value. This then allows the application, to the set, of clustering techniques that search for groups of structures with a common geometry. Results are presented of a first application of this technique to a passive scalar field obtained from 5123 direct numerical simulation of scalar mixing by forced, isotropic turbulence (Reλ = 265). These show transition, with decreasing scale, from blob-like structures in the larger scales to blob- and tube-like structures with small or moderate stretching in the inertial range of scales, and then toward tube and, predominantly, sheet-like structures with high level of stretching in the dissipation range of scales. Implications of these results for the dynamical behaviour of passive scalar stirring and mixing by turbulence are discussed

Geometry of enstrophy and dissipation, grid resolution effects and proximity issues in turbulence

We perform a multi-scale non-local geometrical analysis of the structures extracted from the enstrophy and kinetic energy dissipation-rate, instantaneous fields of a numerical database of incompressible homogeneous isotropic turbulence decaying in time obtained by DNS in a periodic box. Three different resolutions are considered: 256^3, 512^3 and 1024^3 grid points, with k_(max)η(overbar) approximately 1, 2 and 4, respectively, the same initial conditions and Re_λ ≈ 77. This allows a comparison of the geometry of the structures obtained for different resolutions. For the highest resolution, structures of enstrophy and dissipation evolve in a continuous distribution from blob-like and moderately stretched tube-like shapes at the large scales to highly stretched sheet-like structures at the small scales. The intermediate scales show a predominance of tube-like structures for both fields, much more pronounced for the enstrophy field. The dissipation field shows a tendency towards structures with lower curvedness than those of the enstrophy, for intermediate and small scales. The 256^3 grid resolution case (k_(max)η(overbar) ≈ 1) was unable to detect the predominance of highly stretched sheet-like structures at the smaller scales in both fields. The same non-local methodology for the study of the geometry of structures, but without the multi-scale decomposition, is applied to two scalar fields used by existing local criteria for the eduction of tube- and sheet-like structures in turbulence, Q and [A_ij]_+, respectively, obtained from invariants of the velocity-gradient tensor and alike in the 1024^3 case. This adds the non-local geometrical characterization and classification to those local criteria, assessing their validity in educing particular geometries. Finally, we introduce a new methodology for the study of proximity issues among structures of different fields, based on geometrical considerations and non-local analysis, by taking into account the spatial extent of the structures. We apply it to the four fields previously studied. Tube-like structures of Q are predominantly surrounded by sheet-like structures of [A_ij]_+, which appear at closer distances. For the enstrophy, tube-like structures at an intermediate scale are primarily surrounded by sheets of smaller scales of the enstrophy and structures of dissipation at the same and smaller scales. A secondary contribution results from tubes of enstrophy at smaller scales appearing at farther distances. Different configurations of composite structures are presented

Dislocation microstructures and strain-gradient plasticity with one active slip plane

We study dislocation networks in the plane using the vectorial phase-field model introduced by Ortiz and coworkers, in the limit of small lattice spacing. We show that, in a scaling regime where the total length of the dislocations is large, the phase field model reduces to a simpler model of the strain-gradient type. The limiting model contains a term describing the three-dimensional elastic energy and a strain-gradient term describing the energy of the geometrically necessary dislocations, characterized by the tangential gradient of the slip. The energy density appearing in the strain-gradient term is determined by the solution of a cell problem, which depends on the line tension energy of dislocations. In the case of cubic crystals with isotropic elasticity our model shows that complex microstructures may form, in which dislocations with different Burgers vector and orientation react with each other to reduce the total self energy

Nanocalorimetric Evidence for Nematic Superconductivity in the Doped Topological Insulator Sr0.1_{0.1}Bi2_{2}Se3_{3}

Spontaneous rotational-symmetry breaking in the superconducting state of doped $\mathrm{Bi}_2\mathrm{Se}_3$ has attracted significant attention as an indicator for topological superconductivity. In this paper, high-resolution calorimetry of the single-crystal $\mathrm{Sr}_{0.1}\mathrm{Bi}_2\mathrm{Se}_3$ provides unequivocal evidence of a two-fold rotational symmetry in the superconducting gap by a \emph{bulk thermodynamic} probe, a fingerprint of nematic superconductivity. The extremely small specific heat anomaly resolved with our high-sensitivity technique is consistent with the material's low carrier concentration proving bulk superconductivity. The large basal-plane anisotropy of $H_{c2}$ is attributed to a nematic phase of a two-component topological gap structure $\vec{\eta} = (\eta_{1}, \eta_{2})$ and caused by a symmetry-breaking energy term $\delta (|\eta_{1}|^{2} - |\eta_{2}|^{2}) T_{c}$. A quantitative analysis of our data excludes more conventional sources of this two-fold anisotropy and provides the first estimate for the symmetry-breaking strength $\delta \approx 0.1$, a value that points to an onset transition of the second order parameter component below 2K