Fast modal extraction in NASTRAN via the FEER computer program

A new eigensolution routine, FEER (Fast Eigensolution Extraction Routine), used in conjunction with NASTRAN at Israel Aircraft Industries is described. The FEER program is based on an automatic matrix reduction scheme whereby the lower modes of structures with many degrees of freedom can be accurately extracted from a tridiagonal eigenvalue problem whose size is of the same order of magnitude as the number of required modes. The process is effected without arbitrary lumping of masses at selected node points or selection of nodes to be retained in the analysis set. The results of computational efficiency studies are presented, showing major arithmetic operation counts and actual computer run times of FEER as compared to other methods of eigenvalue extraction, including those available in the NASTRAN READ module. It is concluded that the tridiagonal reduction method used in FEER would serve as a valuable addition to NASTRAN for highly increased efficiency in obtaining structural vibration modes

Placing three-dimensional isoparametric elements into NASTRAN

Linear (8 node), parabolic (20 node), cubic (32 node) and mixed (some edges linear, some parabolic and some cubic) have been inserted into NASTRAN, level 15.1. First the dummy element feature was used to check out the stiffness matrix generation routines for the linear element in NASTRAN. Then, the necessary modules of NASTRAN were modified to include the new family of elements. The matrix assembly was changed so that the stiffness matrix of each isoparametric element is only generated once as the time to generate these higher order elements tends to be much longer than the other elements in NASTRAN. This paper presents some of the experiences and difficulties of inserting a new element or family of elements into NASTRAN

Blocking and Persistence in the Zero-Temperature Dynamics of Homogeneous and Disordered Ising Models

A ``persistence'' exponent theta has been extensively used to describe the nonequilibrium dynamics of spin systems following a deep quench: for zero-temperature homogeneous Ising models on the d-dimensional cubic lattice, the fraction p(t) of spins not flipped by time t decays to zero like t^[-theta(d)] for low d; for high d, p(t) may decay to p(infinity)>0, because of ``blocking'' (but perhaps still like a power). What are the effects of disorder or changes of lattice? We show that these can quite generally lead to blocking (and convergence to a metastable configuration) even for low d, and then present two examples --- one disordered and one homogeneous --- where p(t) decays exponentially to p(infinity).Comment: 8 pages (LaTeX); to appear in Physical Review Letter

Non-equilibrium Dynamics of Finite Interfaces

We present an exact solution to an interface model representing the dynamics of a domain wall in a two-phase Ising system. The model is microscopically motivated, yet we find that in the scaling regime our results are consistent with those obtained previously from a phenomenological, coarse-grained Langevin approach.Comment: 12 pages LATEX (figures available on request), Oxford preprint OUTP-94-07

Spinning BTZ Black Hole versus Kerr Black Hole : A Closer Look

By applying Newman's algorithm, the AdS_3 rotating black hole solution is ``derived'' from the nonrotating black hole solution of Banados, Teitelboim, and Zanelli (BTZ). The rotating BTZ solution derived in this fashion is given in ``Boyer-Lindquist-type'' coordinates whereas the form of the solution originally given by BTZ is given in a kind of an ``unfamiliar'' coordinates which are related to each other by a transformation of time coordinate alone. The relative physical meaning between these two time coordinates is carefully studied. Since the Kerr-type and Boyer-Lindquist-type coordinates for rotating BTZ solution are newly found via Newman's algorithm, next, the transformation to Kerr-Schild-type coordinates is looked for. Indeed, such transformation is found to exist. And in this Kerr-Schild-type coordinates, truely maximal extension of its global structure by analytically continuing to ``antigravity universe'' region is carried out.Comment: 17 pages, 1 figure, Revtex, Accepted for publication in Phys. Rev.

Spectral densities of scale-free networks

The spectral densities of the weighted Laplacian, random walk and weighted adjacency matrices associated with a random complex network are studied using the replica method. The link weights are parametrized by a weight exponent $\beta$. Explicit results are obtained for scale-free networks in the limit of large mean degree after the thermodynamic limit, for arbitrary degree exponent and $\beta$.Comment: 14 pages, two figure

Random acyclic networks

Directed acyclic graphs are a fundamental class of networks that includes citation networks, food webs, and family trees, among others. Here we define a random graph model for directed acyclic graphs and give solutions for a number of the model's properties, including connection probabilities and component sizes, as well as a fast algorithm for simulating the model on a computer. We compare the predictions of the model to a real-world network of citations between physics papers and find surprisingly good agreement, suggesting that the structure of the real network may be quite well described by the random graph.Comment: 4 pages, 2 figure

Measurement of the interaction strength in a Bose-Fermi mixture with 87Rb and 40K

A quantum degenerate, dilute gas mixture of bosonic and fermionic atoms was produced using 87Rb and 40K. The onset of degeneracy was confirmed by observing the spatial distribution of the gases after time-of-flight expansion. Further, the magnitude of the interspecies scattering length between the doubly spin polarized states of 87Rb and 40K, |a_RbK|, was determined from cross-dimensional thermal relaxation. The uncertainty in this collision measurement was greatly reduced by taking the ratio of interspecies and intraspecies relaxation rates, yielding |a_RbK| = 250 +/- 30 a_0, which is a lower value than what was reported in [M. Modugno et al., Phys. Rev. A 68, 043626 (2003)]. Using the value for |a_RbK| reported here, current T=0 theory would predict a threshold for mechanical instability that is inconsistent with the experimentally observed onset for sudden loss of fermions in [G. Modugno et al., Science 297, 2240 (2002)].Comment: RevTeX4 + 4 eps figures; Replaced with published versio