236 research outputs found
Measurement of Indeterminacy in Packings of Perfectly Rigid Disks
Static packings of perfectly rigid particles are investigated theoretically
and numerically. The problem of finding the contact forces in such packings is
formulated mathematically. Letting the values of the contact forces define a
vector in a high-dimensional space enable us to show that the set of all
possible contact forces is convex, facilitating its numerical exploration. It
is also found that the boundary of the set is connected with the presence of
sliding contacts, suggesting that a stable packing should not have more than
2M-3N sliding contacts in two dimensions, where M is the number of contacts and
N is the number of particles.
These results were used to analyze packings generated in different ways by
either molecular dynamics or contact dynamics simulations. The dimension of the
set of possible forces and the number of sliding contacts agrees with the
theoretical expectations. The indeterminacy of each component of the contact
forces are found, as well as the an estimate for the diameter of the set of
possible contact forces. We also show that contacts with high indeterminacy are
located on force chains. The question of whether the simulation methods can
represent a packing's memory of its formation is addressed.Comment: 12 pages, 13 figures, submitted to Phys Rev
Elastic behavior in Contact Dynamics of rigid particles
The systematic errors due to the practical implementation of the Contact
Dynamics method for simulation of dense granular media are examined. It is
shown that, using the usual iterative solver to simulate a chain of rigid
particles, effective elasticity and sound propagation with a finite velocity
occur. The characteristics of these phenomena are investigated analytically and
numerically in order to assess the limits of applicability of this simulation
method and to compare it with soft particle molecular dynamics.Comment: submitted to PRE, 7 pages, 6 figure
Motion of a driven tracer particle in a one-dimensional symmetric lattice gas
We study the dynamics of a tracer particle subject to a constant driving
force in a one-dimensional lattice gas of hard-core particles whose
transition rates are symmetric. We show that the mean displacement of the
driven tracer grows in time, , as , rather than the linear
time dependence found for driven diffusion in the bath of non-interacting
(ghost) particles. The prefactor is determined implicitly, as the
solution of a transcendental equation, for an arbitrary magnitude of the
driving force and an arbitrary concentration of the lattice gas particles. In
limiting cases the prefactor is obtained explicitly. Analytical predictions are
seen to be in a good agreement with the results of numerical simulations.Comment: 21 pages, LaTeX, 4 Postscript fugures, to be published in Phys. Rev.
E, (01Sep, 1996
Infrastructure for Detector Research and Development towards the International Linear Collider
The EUDET-project was launched to create an infrastructure for developing and
testing new and advanced detector technologies to be used at a future linear
collider. The aim was to make possible experimentation and analysis of data for
institutes, which otherwise could not be realized due to lack of resources. The
infrastructure comprised an analysis and software network, and instrumentation
infrastructures for tracking detectors as well as for calorimetry.Comment: 54 pages, 48 picture
Search for composite and exotic fermions at LEP 2
A search for unstable heavy fermions with the DELPHI detector at LEP is
reported. Sequential and non-canonical leptons, as well as excited leptons and
quarks, are considered. The data analysed correspond to an integrated
luminosity of about 48 pb^{-1} at an e^+e^- centre-of-mass energy of 183 GeV
and about 20 pb^{-1} equally shared between the centre-of-mass energies of 172
GeV and 161 GeV. The search for pair-produced new leptons establishes 95%
confidence level mass limits in the region between 70 GeV/c^2 and 90 GeV/c^2,
depending on the channel. The search for singly produced excited leptons and
quarks establishes upper limits on the ratio of the coupling of the excited
fermio
Search for lightest neutralino and stau pair production in light gravitino scenarios with stau NLSP
Promptly decaying lightest neutralinos and long-lived staus are searched for
in the context of light gravitino scenarios. It is assumed that the stau is the
next to lightest supersymmetric particle (NLSP) and that the lightest
neutralino is the next to NLSP (NNLSP). Data collected with the Delphi detector
at centre-of-mass energies from 161 to 183 \GeV are analysed. No evidence of
the production of these particles is found. Hence, lower mass limits for both
kinds of particles are set at 95% C.L.. The mass of gaugino-like neutralinos is
found to be greater than 71.5 GeV/c^2. In the search for long-lived stau,
masses less than 70.0 to 77.5 \GeVcc are excluded for gravitino masses from 10
to 150 \eVcc . Combining this search with the searches for stable heavy leptons
and Minimal Supersymmetric Standard Model staus a lower limit of 68.5 \GeVcc
may be set for the stau mas
Development and Validation of a Tokamak Skin Effect Transformer model
A control oriented, lumped parameter model for the tokamak transformer
including the slow flux penetration in the plasma (skin effect transformer
model) is presented. The model does not require detailed or explicit
information about plasma profiles or geometry. Instead, this information is
lumped in system variables, parameters and inputs. The model has an exact
mathematical structure built from energy and flux conservation theorems,
predicting the evolution and non linear interaction of the plasma current and
internal inductance as functions of the primary coil currents, plasma
resistance, non-inductive current drive and the loop voltage at a specific
location inside the plasma (equilibrium loop voltage). Loop voltage profile in
the plasma is substituted by a three-point discretization, and ordinary
differential equations are used to predict the equilibrium loop voltage as
function of the boundary and resistive loop voltages. This provides a model for
equilibrium loop voltage evolution, which is reminiscent of the skin effect.
The order and parameters of this differential equation are determined
empirically using system identification techniques. Fast plasma current
modulation experiments with Random Binary Signals (RBS) have been conducted in
the TCV tokamak to generate the required data for the analysis. Plasma current
was modulated in Ohmic conditions between 200kA and 300kA with 30ms rise time,
several times faster than its time constant L/R\approx200ms. The model explains
the most salient features of the plasma current transients without requiring
detailed or explicit information about resistivity profiles. This proves that
lumped parameter modeling approach can be used to predict the time evolution of
bulk plasma properties such as plasma inductance or current with reasonable
accuracy; at least in Ohmic conditions without external heating and current
drive sources
Exploration of the equilibrium operating space for NSTX-Upgrade
This paper explores a range of high-performance equilibrium scenarios available in the NSTX-Upgrade device [J.E. Menard, submitted for publication to Nuclear Fusion]. NSTX-Upgrade is a substantial upgrade to the existing NSTX device [M. Ono, et al., Nuclear Fusion 40, 557 (2000)], with significantly higher toroidal field and solenoid capabilities, and three additional neutral beam sources with significantly larger current drive efficiency. Equilibria are computed with freeboundary TRANSP, allowing a self consistent calculation of the non-inductive current drive sources, the plasma equilibrium, and poloidal field coil current, using the realistic device geometry. The thermal profiles are taken from a variety of existing NSTX discharges, and different assumptions for the thermal confinement scalings are utilized. The no-wall and idealwall n=1 stability limits are computed with the DCON code. The central and minimum safety factors are quite sensitive to many parameters: they generally increases with large outer plasmawall gaps and higher density, but can have either trend with the confinement enhancement factor. In scenarios with strong central beam current drive, the inclusion of non-classical fast ion diffusion raises qmin, decreases the pressure peaking, and generally improves the global stability, at the expense of a reduction in the non-inductive current drive fraction; cases with less beam current drive are largely insensitive to additional fast ion diffusion. The non-inductive current level is quite sensitive to the underlying confinement and profile assumptions. For instance, for BT=1.0 T and Pinj=12.6 MW, the non-inductive current level varies from 875 kA with ITER-98y,2 thermal confinement scaling and narrow thermal profiles to 1325 kA for an ST specific scaling expression and broad profiles. This sensitivity should facilitate the determination of the correct scaling of transport with current and field to use for future fully non-inductive ST devices. Scenarios are presented which can be sustained for 8-10 seconds, or (20-30)τCR, at βN=3.8-4.5, facilitating, for instance, the study of disruption avoidance for very long pulse. Scenarios have been documented which can operate with βT~25% and equilibrated qmin>1. The value of qmin can be controlled at either fixed non-inductive fraction of 100% or fixed plasma current, by varying which beam sources are used, opening the possibility for feedback qmin control. In terms of quantities like collisionality, neutron emission, non-inductive fraction, or stored energy, these scenarios represent a significant performance extension compared to NSTX and other present spherical torii
- …