Momentum equations integral formulation

The Lippmann—Schwinger equations (6.73) are written formally in terms of a discrete notation i) for the complete set of target states, which includes the ionisation continuum. For a numerical solution it is necessary to have a finite set of coupled integral equations. We formulate the coupled-channels-optical equations that describe reactions in a channel subspace, called P space. This is projected from the chaimel space by an operator P that includes only a finite set of target states. The entrance channel 0ko) is included in P space. The method was first discussed by Feshbach (1962). Its application to the momentum-space formulation of electron—atom scattering was introduced by McCarthy and Stelbovics... 

The derivation of the given whole field formulation, introducing the Dirac delta function (d/) into the surface tension force relation to maintain the discontinuous (singular) nature of this term, is to a certain extent based on physical intuition rather than first principles (i.e., in mathematical terms this approach is strictly not characterized as a continuum formulation on the differential form). Chandrasekhar [31] (pp 430-433) derived a similar model formulation and argued that to some extent the whole field momentum equation can be obtained by a formal mathematical procedure. However, the fact that the equation involves /-functions means that to interpret the equation correctly at a point of discontinuity, we must integrate the equation, across the interface, over an infinitesimal volume element including the discontinuity. 

In order to properly solve (17.5), sharp changes in the properties as well as pressure forces due to surface tension effects have to be resolved. In particular, surface tension results in a jump in pressure across a curved interface. The pressure jump is discontinuous and located only at the interface. This singularity creates difficulties when deriving a continuum formulation of the momentum equation. The interfacial conditions should be embedded in the field equations as source terms. Once the equations are discretized in a finite-thickness interfacial zone, the fiow properties are allowed to change smoothly. It is therefore necessary to create a continuum surface force (CSF) equal to the surface tension at the interface, or in a transitional region, and zero elsewhere. Therefore, the surface integral term in (17.5) could be rewritten into an appropriate volume integral... 

For slow flames, Eq. (4.79) may be uncoupled from the remainder of the calculation (as has been done so far), and Eq. (2.7b) may be used to determine the steady-state pressure profile at the end of the integration. For much faster flames where there are appreciable gasdynamic effects and associated density changes, the momentum equation must be coupled directly into the system, and the energy equations (2.19), (2.20) or (2.20q) must be used in place of Eq. (2.20b). In the finite-difference formulation discussed in Section 4.2, it then also becomes necessary to modify Eq. (4.44) to include the effect of variable pressure on the density and to introduce the condition... 

The formulation outlined above is in configuration space, but several authors, notably Ghosh and his collaborators (Chaudhury, Ghosh and Sil, 1974) and Mitroy (1993), also Mitroy, Berge and Stelbovics (1994) and Mitroy and Ratnavelu (1995), have preferred to work in momentum space with a set of coupled integral equations rather than the coupled integro-differential equations (3.31) and (3.32). 

Next, we apply Galerkin s weighted residual method and reduce the order of integration of the various terms in the above equations using the Green-Gauss Theorem (9.1.2) for each element. For a simpler presentation we will deal with each term in the above equations separately. The terms of the x-component (eqn. (9.95)) of the penalty formulation momentum balance become... 

The nodal unknows a are to be chosen so as to satisfy the governing equations in an integral sense this can be done by using a Galerkin weighted residual formulation of the conservation equations for momentum and energy transport ... 

It is desirable to construct formulations and numerical methods which exactly (i.e. up to rounding error) preserve the total momentum from step to step. One obvious approaeh to this problem is to simply project the momenta onto the linear momenrnm constraint at the end of each step (or after some number of steps). Such a projection introduces potential issues in terms of convergence order and would certainly complicate the analyses presented thus far in this book. Moreover, the optimal choice of projection is unclear and it is easy to define poor schemes (for example, modifying always the momentum of just the first particle in order to balance all the remaining components) which are likely to introduce artifacts (bias) in simulation. For this reason, there is interest in building in momentum conservation into the equations of motion (and indeed the integrator). Ideally this should be done in a localized and homogeneous way so that momentum is not transferred by a nonphysical mechanism between distant particles. 

It should be recalled that, because of the presence of the external potential and the nonlocal form of Ep given by Eq. (11.11), all operators resulting from these unitary transformations are well defined only in momentum space (compare the discussion of the square-root operator in the context of the Klein-Gordon equation in chapter 5 and the momentum-space formulation of the Dirac equation in section 6.10). Whereas So acts as a simple multiplicative operator, all higher-order terms containing the potential V are integral operators and completely described by specifying their kernel. For example, the... 

A collision between the two bodies is known as an impact during which forces are created that act and disappear over a short period of time. The duration of the contact period governs the choice of the method used to analyze the impact. The methods for predicting the impact responses can primarily be classified into two groups. In one, the impact is treated as a discontinuous event. Momentum-balance/impulse equations are usually formulated by integrating the acceleration-based form or the canonical form of the governing equations of motion. The solution to these equations gives the jump in the... 

---