5,666 research outputs found
Asymptotic integration algorithms for nonhomogeneous, nonlinear, first order, ordinary differential equations
New methods for integrating systems of stiff, nonlinear, first order, ordinary differential equations are developed by casting the differential equations into integral form. Nonlinear recursive relations are obtained that allow the solution to a system of equations at time t plus delta t to be obtained in terms of the solution at time t in explicit and implicit forms. Examples of accuracy obtained with the new technique are given by considering systems of nonlinear, first order equations which arise in the study of unified models of viscoplastic behaviors, the spread of the AIDS virus, and predator-prey populations. In general, the new implicit algorithm is unconditionally stable, and has a Jacobian of smaller dimension than that which is acquired by current implicit methods, such as the Euler backward difference algorithm; yet, it gives superior accuracy. The asymptotic explicit and implicit algorithms are suitable for solutions that are of the growing and decaying exponential kinds, respectively, whilst the implicit Euler-Maclaurin algorithm is superior when the solution oscillates, i.e., when there are regions in which both growing and decaying exponential solutions exist
A viscoplastic model with application to LiF-22 percent CaF2 hypereutectic salt
A viscoplastic model for class M (metal-like behavior) materials is presented. One novel feature is its use of internal variables to change the stress exponent of creep (where n is approximately = 5) to that of natural creep (where n = 3), in accordance with experimental observations. Another feature is the introduction of a coupling in the evolution equations of the kinematic and isotropic internal variables, making thermal recovery of the kinematic variable implicit. These features enable the viscoplastic model to reduce to that of steady-state creep in closed form. In addition, the hardening parameters associated with the two internal state variables (one scalar-valued, the other tensor-valued) are considered to be functions of state, instead of being taken as constant-valued. This feature enables each internal variable to represent a much wider spectrum of internal states for the material. The model is applied to a LiF-22 percent CaF2 hypereutectic salt, which is being considered as a thermal energy storage material for space-based solar dynamic power systems
Steady-state and transient Zener parameters in viscoplasticity: Drag strength versus yield strength
A hypothesis is put forth which enables the viscoplastician to formulate a theory of viscoplasticity that reduces, in closed form, to the classical theory of creep. This hypothesis is applied to a variety of drag and yield strength models. Because of two theoretical restrictions that are a consequence of this hypothesis, three different yield strength models and one drag strength model are shown to be theoretically admissible. One of these yield strength models is selected as being the most appropriate representation for isotropic hardening
Dynamics of Diblock Copolymers in Dilute Solutions
We consider the dynamics of freely translating and rotating diblock (A-B),
Gaussian copolymers, in dilute solutions. Using the multiple scattering
technique, we have computed the diffusion and the friction coefficients D_AB
and Zeta_AB, and the change Eta_AB in the viscosity of the solution as
functions of x = N_A/N and t = l_B/l_A, where N_A, N are the number of segments
of the A block and of the whole copolymer, respectively, and l_A, l_B are the
Kuhn lengths of the A and B blocks. Specific regimes that maximize the
efficiency of separation of copolymers with distinct "t" values, have been
identified.Comment: 20 pages Revtex, 7 eps figures, needs epsf.tex and amssymb.sty,
submitted to Macromolecule
Stress versus temperature dependent activation energies in creep
The activation energy for creep at low stresses and elevated temperatures is lattice diffusion, where the rate controlling mechanism for deformation is dislocation climb. At higher stresses and intermediate temperatures, the rate controlling mechanism changes from that of dislocation climb to one of obstacle-controlled dislocation glide. Along with this change, there occurs a change in the activation energy. It is shown that a temperature-dependent Gibbs free energy does a good job of correlating steady-state creep data, while a stress-dependent Gibbs free energy does a less desirable job of correlating the same data. Applications are made to copper and a LiF-22 mol. percent CaF2 hypereutectic salt
On the thermodynamics of stress rate in the evolution of back stress in viscoplasticity
A thermodynamic foundation using the concept of internal state variables is presented for the kinematic description of a viscoplastic material. Three different evolution equations for the back stress are considered. The first is that of classical, nonlinear, kinematic hardening. The other two include a contribution that is linear in stress rate. Choosing an appropriate change in variables can remove this stress rate dependence. As a result, one of these two models is shown to be equivalent to the classical, kinematic hardening model; while the other is a new model, one which seems to have favorable characteristics for representing ratchetting behavior. All three models are thermodynamically admissible
Critical behavior of long straight rigid rods on two-dimensional lattices: Theory and Monte Carlo simulations
The critical behavior of long straight rigid rods of length (-mers) on
square and triangular lattices at intermediate density has been studied. A
nematic phase, characterized by a big domain of parallel -mers, was found.
This ordered phase is separated from the isotropic state by a continuous
transition occurring at a intermediate density . Two analytical
techniques were combined with Monte Carlo simulations to predict the dependence
of on , being . The first involves
simple geometrical arguments, while the second is based on entropy
considerations. Our analysis allowed us also to determine the minimum value of
(), which allows the formation of a nematic phase on a
triangular lattice.Comment: 23 pages, 5 figures, to appear in The Journal of Chemical Physic
Geometric Aspects of D-branes and T-duality
We explore the differential geometry of T-duality and D-branes. Because
D-branes and RR-fields are properly described via K-theory, we discuss the
(differential) K-theoretic generalization of T-duality and its application to
the coupling of D-branes to RR-fields. This leads to a puzzle involving the
transformation of the A-roof genera in the coupling.Comment: 26 pages, JHEP format, uses dcpic.sty; v2: references added, v3:
minor change
Apex Exponents for Polymer--Probe Interactions
We consider self-avoiding polymers attached to the tip of an impenetrable
probe. The scaling exponents and , characterizing the
number of configurations for the attachment of the polymer by one end, or at
its midpoint, vary continuously with the tip's angle. These apex exponents are
calculated analytically by -expansion, and numerically by simulations
in three dimensions. We find that when the polymer can move through the
attachment point, it typically slides to one end; the apex exponents quantify
the entropic barrier to threading the eye of the probe
On the Quantum Invariant for the Brieskorn Homology Spheres
We study an exact asymptotic behavior of the Witten-Reshetikhin-Turaev
invariant for the Brieskorn homology spheres by use of
properties of the modular form following a method proposed by Lawrence and
Zagier. Key observation is that the invariant coincides with a limiting value
of the Eichler integral of the modular form with weight 3/2. We show that the
Casson invariant is related to the number of the Eichler integrals which do not
vanish in a limit . Correspondingly there is a
one-to-one correspondence between the non-vanishing Eichler integrals and the
irreducible representation of the fundamental group, and the Chern-Simons
invariant is given from the Eichler integral in this limit. It is also shown
that the Ohtsuki invariant follows from a nearly modular property of the
Eichler integral, and we give an explicit form in terms of the L-function.Comment: 26 pages, 2 figure
- …