Diffusion equation linearization

V = vapor rate, Ib/hr ako, used as molecular volume for component in diffusivity equation or, linear velocity, ft/sec. r... 

The strategies discussed in the previous chapter are generally applicable to convection-diffusion equations such as Eq. (32). If the function O is a component of the velocity field, the incompressible Navier-Stokes equation, a non-linear partial differential equation, is obtained. This stands in contrast to O representing a temperature or concentration field. In these cases the velocity field is assumed as given, and only a linear partial differential equation has to be solved. The non-linear nature of the Navier-Stokes equation introduces some additional problems, for which special solution strategies exist. Corresponding numerical techniques are the subject of this section. 

In all considered above models, the equilibrium morphology is chosen from the set of possible candidates, which makes these approaches unsuitable for discovery of new unknown structures. However, the SCFT equation can be solved in the real space without any assumptions about the phase symmetry [130], The box under the periodic boundary conditions in considered. The initial quest for uy(r) is produced by a random number generator. Equations (42)-(44) are used to produce density distributions T(r) and pressure field ,(r). The diffusion equations are numerically integrated to obtain q and for 0 < s < 1. The right-hand size of Eq. (47) is evaluated to obtain new density profiles. The volume fractions at the next iterations are obtained by a linear mixing of new and old solutions. The iterations are performed repeatedly until the free-energy change... 

In the particular case dealt with now (fully labile complexation), due to the linearity of a combined diffusion equation for DmCm + DmlL ml, the flux in equation (65) can still be seen as the sum of the independent diffusional fluxes of metal and complex, each contribution depending on the difference between the surface and bulk concentration value of each species. But equation (66) warns against using just a rescaling factor for the total metal or for the free metal alone. In general, if the diffusion is coupled with some nonlinear process, the resulting flux is not proportional to bulk-to-surface differences, and this complicates the use of mass transfer coefficients (see ref. [II] or Chapter 3 in this volume). 

As discussed in Section 4.3, the linear-eddy model solves a one-dimensional reaction-diffusion equation for all length scales. Inertial-range fluid-particle interactions are accounted for by a random rearrangement process. This leads to significant computational inefficiency since step (3) is not the rate-controlling step. Simplifications have thus been introduced to avoid this problem (Baldyga and Bourne 1989). 

The observations above can be rapidly turned into a semi-quantitative theory for star-polymer stress-relaxation [24] which is amenable to more quantitative refinement [25]. The key observation is that the diffusion equation for stress-re-lease, which arises in linear polymers via the passage of free ends out of deformed tube segment, is now modified in star polymers by the potential of Eq. (16). Apart from small displacements of the end, the diffusion to any position s along the arm will now need to be activated and so is exponentially suppressed. Each position along the arm, s, will possess its own characteristic stress relaxation time T(s) given approximately by... 

The diffusion coefficient D can be thought of as the velocity of the analyte as it moves from the bulk of the solution towards the electrode just prior to the electron-transfer reaction. Because D is a velocity, larger values of D relate to a faster motion of analyte through the solution, while smaller values relate to slower motion. It is assumed, when deriving equation (6.3), that diffusion is linear... 

Transient is a C-program for solving systems of generally non-linear, parabolic partial differential equations in two variables (that is, space and time), in particular, reaction-diffusion equations within the generalized Crank-Nicolson Finite Difference Method. 

Notice that the right-hand side of Eq. (34) is equal to the ratio of the transformed concentration at the second measurement point to the transformed concentration at the first measurement point. In the terminology of control engineering, this quantity is the transfer function of the system between Xo and Xm- The Laplace-transform method is possible because the diffusion equation is a linear differential equation. Thus, the right-hand side of Eq. (34) could in principle be used in a control-system analysis of an axial-dispersion process. 

Prom the diffusion equation, it can be seen that if d Cldx = 0, then dCldt = 0, meaning that the concentration at the position would not vary with time. Hence, if the initial concentration is uniform, the concentration would not change with time. If the initial concentration profile is linear and the concentrations at the two ends are not changed from linear distribution, because... 

The diffusion equation with constant diffusivity (Equation 3-8) is said to be linear, which means that if fand g are solutions to the equation, then any linear combination of f and g, i.e., u = af+ bg, where a and b are constants, is also a solution. To show this, we can write... 

The above derivation has not made use of the initial and boundary conditions yet, and shows only that A may take any constant value. The value of A can be constrained by boundary conditions to be discrete Ai, A2,..., as can be seen in the specific problem below. Because each function corresponding to given A is a solution to the diffusion equation, based on the principle of superposition, any linear combination of these functions is also a solution. Hence, the general solution for the given boundary conditions is... 

The superposition principle can be used to combine solutions for linear partial differential equations, like the diffusion equation. It is stated as follows ... 

A linear PDF, such as the diffusion equation, can only have first-order and zero-order source or sink rates. But what if the source or sink term is of higher order An example would be the generalized reaction... 

Rate equations There are two basic types of kinetic rate expressions. The first and simpler is the case of linear diffusion equations or linear driving forces (LDF) and the second and more rigorous is the case of classic Fickian differential equations. 

Linear diffusion equations This is the simplest case and is used extensively in the related literature (Perry and Green, 1999 Hashimoto et al., 1977 Cooney, 1990, 1993). The equations are the following ... 

Equation (1.57a) implies that in the locally electro-neutral ambipolar diffusion concentration of both ions evolves according to a single linear diffusion equation with an effective diffusivity given by (1.57b). Physically, the role of the electric field, determined from the elliptic current continuity equation... 

For illustration, we present in Fig. 3.3.1a,b the results of a numerical solution of the original system (3.3.19) (Curve 1) for e = 10-2, 10-3, 7=1 together with a plot of the leading term (3.3.48) (Curve 2). We also present for comparison a plot of 1 erfcx (Curve 3), the similarity solution for the linear diffusion equation with the boundary and initial conditions analogous... 

The local stability of a given stationary-state profile can be determined by the same sort of test applied to the solutions for a CSTR. Of course now, when we substitute in a = ass + Aa etc., we have the added complexity that the profile is a function of position, as may be the perturbation. Stability and instability again are distinguished by the decay or growth of these small perturbations, and except for special circumstances the governing reaction-diffusion equation for SAa/dr will be a linear second-order partial differential equation. Thus the time dependence of Aa will be governed by an infinite series of exponential terms ... 

Collective Modes and Time-Correlation Functions. Our linear equations (6.15) and (6.16) describe two characteristic kinds of collective modes in gels the longitudinal part of u obeys Tanaka s equation (4.16) and g = V u is governed by the diffusion equation (4.18), while the transverse parts of u and are coupled to form a slow transverse sound at small wave numbers. By assuming the space-time dependence as exp(i[Pg.99]

Fick s second law (Eq. 18-14) is a second-order linear partial differential equation. Generally, its solutions are exponential functions or integrals of exponential functions such as the error function. They depend on the boundary conditions and on the initial conditions, that is, the concentration at a given time which is conveniently chosen as t = 0. The boundary conditions come in different forms. For instance, the concentration may be kept fixed at a wall located atx0. Alternatively, the wall may be impermeable for the substance, thus the flux at x0 is zero. According to Eq. 18-6, this is equivalent to keeping dC/dx = 0 at x0. Often it is assumed that the system is unbounded (i.e., that it extends from x = - °o to + °°). For this case we have to make sure that the solution C(x,t) remains finite when x -a °°. In many cases, solutions are found only by numerical approximations. For simple boundary conditions, the mathematical techniques for the solution of the diffusion equation (such as the Laplace transformation) are extensively discussed in Crank (1975) and Carslaw and Jaeger (1959). 

To avoid the difficulties associated with the spherical diffusion equation, a useful hypothesis is the linear-driving-force concept. This arises when a parabolic concentration profile within the spherical particles is supposed - which is a good approximation in cases where there is a Thiele modulus of a maximum volume of 2-5 (that is, with some intra--particle resistance [50]). In these conditions, the volume-averaged intra-particle concentration is defined as ... 

Methods to solve the diffusion equation for specific boundary and initial conditions are presented in Chapter 5. Many analytic solutions exist for the special case that D is uniform. This is generally not the case for interdiffusivity D (Eq. 3.25). If D does not vary rapidly with composition, it can be replaced by successive approximations of a uniform diffusivity and results in a linearization of the diffusion equation. The... 

Then, because the diffusion equation is a linear second-order differential equation, r(x, ) = p(x,t) + q(x,t) is a solution for the boundary conditions and the initial condition ... 

The anisotropic form of Fick s law would seem to complicate the diffusion equation greatly. However, in many cases, a simple method for treating anisotropic diffusion allows the diffusion equation to keep its simple form corresponding to isotropic diffusion. Because Dtj is symmetric, it is always possible to find a linear coordinate transformation that will make the Dij diagonal with real components (the eigenvalues of D). Let elements of such a transformed system be identified by a hat. Then... 

The general treatment for multicomponent diffusion results in linear systems of diffusion equations. A linear transformation of the concentrations produces a simplified system of uncoupled linear diffusion equations for which general solutions can be obtained by methods presented in Chapter 5. 

---