System with Dynamical Force Elements

The equations of motion presented so far are second order ODEs with or without constraints. In some cases there might be also additional first order differential equations. Consider, for example, the system shown in Fig. 1.6. If the mass m2 is taken to be zero the mass mi is suspended by an interconnection consisting of a spring and a damper in series. 

To establish the equations of motion of that system under the influence of gravitation we will flrst give m2 a nonzero mass. This yields 

If in the second equation m2 is set to zero, we obtain a mixed system consisting in a second and a first order ODE 

Mixed systems of this type arise also when mechanical systems are influenced by non mechanical dynamic forces, like electro-magnetic forces in a magnetically levitated (maglev) vehicle. Also, multibody systems with additional control devices, hydraulic components etc. may result in such a mixed form. 

---

As already noted, the present study of dynamic fuel cell behavior involves the analysis of systems with capacitive elements. These elements control the rate at which process parameters change due to net forces imposed by other coupled process parameters. A general dynamic equation showing capacitance behavior is ... 

Perhaps the most common computer simulation method for nonequilibrium systems is the nonequilibrium molecular dynamics (NEMD) method [53, 88]. This typically consists of Hamilton s equations of motion augmented with an artificial force designed to mimic particular nonequilibrium fluxes, and a constraint force or thermostat designed to keep the kinetic energy or temperature constant. Here is given a brief derivation and critique of the main elements of that method. 

We will introduce basic kinetic concepts that are frequently used and illustrate them with pertinent examples. One of those concepts is the idea of dynamic equilibrium, as opposed to static (mechanical) equilibrium. Dynamic equilibrium at a phase boundary, for example, means that equal fluxes of particles are continuously crossing the boundary in both directions so that the (macroscopic) net flux is always zero. This concept enables us to understand the non-equilibrium state of a system as a monotonic deviation from the equilibrium state. Driven by the deviations from equilibrium of certain functions of state, a change in time for such a system can then be understood as the return to equilibrium. We can select these functions of state according to the imposed constraints. If the deviations from equilibrium are sufficiently small, the result falls within a linear theory of process rates. As long as the kinetic coefficients can be explained in terms of the dynamic equilibrium properties, the reaction rates are directly proportional to the deviations. The thermodynamic equilibrium state is chosen as the reference state in which the driving forces X, vanish, but not the random thermal motions of structure elements i. Therefore, systems which we wish to study kinetically must first be understood at equilibrium, where the SE fluxes vanish individually both in the interior of all phases and across phase boundaries. This concept will be worked out in Section 4.2.1 after fluxes of matter, charge, etc. have been introduced through the formalism of irreversible thermodynamics. 

The sample s dynamic stiffness is obtained from the measured force and motion amplitudes. The formula for this dynamic stiffness is derived below with reference to the Figure 2 measurement system model. The support impedance element in this model represents the combined effects of the flexibility of the downward force train into the shaker field assembly, the mass of the field assembly, and the oscillating magnetic reaction force exerted on the field assembly from the armature. 

If it can be supposed that the process is essentially controlled by the dynamics over a potential surface element dS of dimentions which are small when compared with the dimensions of the system considered, //, n, R, V and J can be regarded as nearly conserved over dS. The second term in Equation 14 // v /2 thus gives the kinetic energy of the motion perpendicular to dS which is available for the barrier crossing. The first and the third term in Equation 14 are nearly conserved on dS. They become exactly conserved in the limit when the range of force goes to zero. 

The quantum mechanical many-body nature of the interatomic forces is taken into account naturally through the Hellmann-Feynman theorem. Since the scheme usually uses a minimal basis set for the electronic structure calculation and the Hamiltonian matrix elements are parametrized, large numbers of atoms can be tackled within the present computer capabilities. One of the distinctive features of this scheme in comparison with other empirical schemes is that all the parameters in the model can be obtained theoretically. It is therefore very useful for studying novel materials where experimental data are not readily available. The scheme has been demonstrated to be a powerful method for studying various structural, dynamical, and electronic properties of covalent systems. 

---