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 $k$ ($k$-mers) on square and triangular lattices at intermediate density has been studied. A nematic phase, characterized by a big domain of parallel $k$-mers, was found. This ordered phase is separated from the isotropic state by a continuous transition occurring at a intermediate density $\theta_c$. Two analytical techniques were combined with Monte Carlo simulations to predict the dependence of $\theta_c$ on $k$, being $\theta_c(k) \propto k^{-1}$. 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 $k$ ($k_{min}=7$), 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 $\gamma_1$ and $\gamma_2$, 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 $\epsilon$-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 $\Sigma(p_1,p_2,p_3)$ 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 $\tau\to N \in \mathbb{Z}$. 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