Computer program for predicting creep behavior of bodies of revolution

Computer program, CRAB, uses finite-element method to calculate creep behavior and predict steady-state stresses in an arbitrary body of revolution subjected to a time-dependent axisymmetric load. Creep strains follow a time hardening law and a Prandtl-Reuss stress-strain relationship

A Lanczos Method for Approximating Composite Functions

We seek to approximate a composite function h(x) = g(f(x)) with a global polynomial. The standard approach chooses points x in the domain of f and computes h(x) at each point, which requires an evaluation of f and an evaluation of g. We present a Lanczos-based procedure that implicitly approximates g with a polynomial of f. By constructing a quadrature rule for the density function of f, we can approximate h(x) using many fewer evaluations of g. The savings is particularly dramatic when g is much more expensive than f or the dimension of x is large. We demonstrate this procedure with two numerical examples: (i) an exponential function composed with a rational function and (ii) a Navier-Stokes model of fluid flow with a scalar input parameter that depends on multiple physical quantities

Spin diffusion of correlated two-spin states in a dielectric crystal

Reciprocal space measurements of spin diffusion in a single crystal of calcium fluoride (CaF$_2$) have been extended to dipolar ordered states. The experimental results for the component of the spin diffusion parallel with the external field are $D_{D}^{||}=29 \pm 3 \times 10^{-12}$ cm$^{2}$/s for the [001] direction and $D_{D}^{||}=33 \pm 4 \times 10^{-12}$ cm$^{2}$/s for the [111] direction. The diffusion rates for dipolar order are significantly faster than those for Zeeman order and are considerably faster than predicted by simple theoretical models. It is suggested that constructive interference in the transport of the two spin state is responsible for this enhancement. As expected the anisotropy in the diffusion rates is observed to be significantly less for dipolar order compared to the Zeeman case.Comment: 4 pages, 2 figures. Resubmitted to PRL - new figure added / discussion expande

A computer program for plotting stress-strain data from compression, tension, and torsion tests of materials

A computer program for plotting stress-strain curves obtained from compression and tension tests on rectangular (flat) specimens and circular-cross-section specimens (rods and tubes) and both stress-strain and torque-twist curves obtained from torsion tests on tubes is presented in detail. The program is written in FORTRAN 4 language for the Control Data 6000 series digital computer with the SCOPE 3.0 operating system and requires approximately 110000 octal locations of core storage. The program has the capability of plotting individual strain-gage outputs and/or the average output of several strain gages and the capability of computing the slope of a straight line which provides a least-squares fit to a specified section of the plotted curve. In addition, the program can compute the slope of the stress-strain curve at any point along the curve. The computer program input and output for three sample problems are presented

Coming to America: Multiple Origins of New World Geckos

Geckos in the Western Hemisphere provide an excellent model to study faunal assembly at a continental scale. We generated a time-calibrated phylogeny, including exemplars of all New World gecko genera, to produce a biogeographic scenario for the New World geckos. Patterns of New World gecko origins are consistent with almost every biogeographic scenario utilized by a terrestrial vertebrate with different New World lineages showing evidence of vicariance, dispersal via temporary land bridge, overseas dispersal, or anthropogenic introductions. We also recovered a strong relationship between clade age and species diversity, with older New World lineages having more species than more recently arrived lineages. Our data provide the first phylogenetic hypothesis for all New World geckos and highlight the intricate origins and ongoing organization of continental faunas. The phylogenetic and biogeographical hypotheses presented here provide an historical framework to further pursue research on the diversification and assembly of the New World herpetofauna

Hydrodynamic approach to coherent nuclear spin transport

We develop a linear response formalism for nuclear spin diffusion in a dipolar coupled solid. The theory applies to the high-temperature, long-wavelength regime studied in the recent experiments of Boutis et al. [Phys. Rev. Lett. 92, 137201 (2004)], which provided direct measurement of interspin energy diffusion in such a system. A systematic expansion of Kubo's formula in the flip-flop term of the Hamiltonian is used to calculate the diffusion coefficients. We show that this approach is equivalent to the method of Lowe and Gade [Phys. Rev. 156, 817 (1967)] and Kaplan [Phys. Rev. B 2, 4578 (1970)], but has several calculational and conceptual advantages. Although the lowest orders in this expansion agree with the experimental results for magnetization diffusion, this is not the case for energy diffusion. Possible reasons for this disparity are suggested.Comment: 7 pages, REVTeX4; Published Versio

Penalized Orthogonal Iteration for Sparse Estimation of Generalized Eigenvalue Problem

We propose a new algorithm for sparse estimation of eigenvectors in generalized eigenvalue problems (GEP). The GEP arises in a number of modern data-analytic situations and statistical methods, including principal component analysis (PCA), multiclass linear discriminant analysis (LDA), canonical correlation analysis (CCA), sufficient dimension reduction (SDR) and invariant co-ordinate selection. We propose to modify the standard generalized orthogonal iteration with a sparsity-inducing penalty for the eigenvectors. To achieve this goal, we generalize the equation-solving step of orthogonal iteration to a penalized convex optimization problem. The resulting algorithm, called penalized orthogonal iteration, provides accurate estimation of the true eigenspace, when it is sparse. Also proposed is a computationally more efficient alternative, which works well for PCA and LDA problems. Numerical studies reveal that the proposed algorithms are competitive, and that our tuning procedure works well. We demonstrate applications of the proposed algorithm to obtain sparse estimates for PCA, multiclass LDA, CCA and SDR. Supplementary materials are available online