10,380 research outputs found
Solving large-scale dynamic systems using band Lanczos method in Rockwell NASTRAN on CRAY X-MP
The improved cost effectiveness using better models, more accurate and faster algorithms and large scale computing offers more representative dynamic analyses. The band Lanczos eigen-solution method was implemented in Rockwell's version of 1984 COSMIC-released NASTRAN finite element structural analysis computer program to effectively solve for structural vibration modes including those of large complex systems exceeding 10,000 degrees of freedom. The Lanczos vectors were re-orthogonalized locally using the Lanczos Method and globally using the modified Gram-Schmidt method for sweeping rigid-body modes and previously generated modes and Lanczos vectors. The truncated band matrix was solved for vibration frequencies and mode shapes using Givens rotations. Numerical examples are included to demonstrate the cost effectiveness and accuracy of the method as implemented in ROCKWELL NASTRAN. The CRAY version is based on RPK's COSMIC/NASTRAN. The band Lanczos method was more reliable and accurate and converged faster than the single vector Lanczos Method. The band Lanczos method was comparable to the subspace iteration method which was a block version of the inverse power method. However, the subspace matrix tended to be fully populated in the case of subspace iteration and not as sparse as a band matrix
Charge redistribution in the formation of one-dimensional lithium wires on Cu(001)
We describe the formation of one-dimensional lithium wires on a Cu(001) substrate, providing an atomic-scale description of the onset of metallization in this prototypical adsorption system. A combination of helium atom scattering and density-functional theory reveals pronounced changes in the electronic charge distribution on the formation of the c(5√2×√2)R45° Li/Cu(001) structure, as in-plane bonds are created. Charge donation from Li-substrate bonds is found to facilitate the formation of stable, bonded, and depolarized chains of Li adatoms that coexist with an interleaved phase of independent adatoms. The resultant overlayer has a commensurate charge distribution and lattice modulations but differs fundamentally from structurally similar charge-density wave systems
Ab Initio Liquid Hydrogen Muon Cooling Simulations with ELMS in ICOOL
This paper presents new theoretical results on the passage of muons through
liquid hydrogen which have been confirmed in a recent experiment. These are
used to demonstrate that muon bunches may be compressed by ionisation cooling
more effectively than suggested by previous calculations.
Muon cooling depends on the differential cross section for energy loss and
scattering of muons. We have calculated this cross section for liquid H2 from
first principles and atomic data, avoiding traditional assumptions. Thence, 2-D
probability maps of energy loss and scattering in mm-scale thicknesses are
derived by folding, and stored in a database. Large first-order correlations
between energy loss and scattering are found for H2, which are absent in other
simulations. This code is named ELMS, Energy Loss & Multiple Scattering. Single
particle trajectories may then be tracked by Monte Carlo sampling from this
database on a scale of 1 mm or less. This processor has been inserted into the
cooling code ICOOL. Significant improvements in 6-D muon cooling are predicted
compared with previous predictions based on GEANT. This is examined in various
geometries. The large correlation effect is found to have only a small effect
on cooling. The experimental scattering observed for liquid H2 in the MUSCAT
experiment has recently been reported to be in good agreement with the ELMS
prediction, but in poor agreement with GEANT simulation.Comment: 6 pages, 3 figure
Domino tilings and the six-vertex model at its free fermion point
At the free-fermion point, the six-vertex model with domain wall boundary
conditions (DWBC) can be related to the Aztec diamond, a domino tiling problem.
We study the mapping on the level of complete statistics for general domains
and boundary conditions. This is obtained by associating to both models a set
of non-intersecting lines in the Lindstroem-Gessel-Viennot (LGV) scheme. One of
the consequence for DWBC is that the boundaries of the ordered phases are
described by the Airy process in the thermodynamic limit.Comment: 14 pages, 8 figure
Minimal Brownian Ratchet: An Exactly Solvable Model
We develop an exactly-solvable three-state discrete-time minimal Brownian
ratchet (MBR), where the transition probabilities between states are
asymmetric. By solving the master equations we obtain the steady-state
probabilities. Generally the steady-state solution does not display detailed
balance, giving rise to an induced directional motion in the MBR. For a reduced
two-dimensional parameter space we find the null-curve on which the net current
vanishes and detailed balance holds. A system on this curve is said to be
balanced. On the null-curve, an additional source of external random noise is
introduced to show that a directional motion can be induced under the zero
overall driving force. We also indicate the off-balance behavior with biased
random noise.Comment: 4 pages, 4 figures, RevTex source, General solution added. To be
appeared in Phys. Rev. Let
Algebraic arctic curves in the domain-wall six-vertex model
The arctic curve, i.e. the spatial curve separating ordered (or `frozen') and
disordered (or `temperate) regions, of the six-vertex model with domain wall
boundary conditions is discussed for the root-of-unity vertex weights. In these
cases the curve is described by algebraic equations which can be worked out
explicitly from the parametric solution for this curve. Some interesting
examples are discussed in detail. The upper bound on the maximal degree of the
equation in a generic root-of-unity case is obtained.Comment: 15 pages, no figures; v2: metadata correcte
- …