The Steady-State Problem

We will discuss briefly the steady-state limit of the normal problem in three dimensions. This was analyzed in detail for the plane case in Sect. 3.11. The formulae developed in that section go over to three dimensions with minimal and obvious changes essentially, these amount to substituting the three-dimensional elastic pressure distribution for the two dimensional form - for example (5.2.32) with /e = /o, replacing (3.11.46). 

However, the depth of penetration was not discussed, since it is indeterminate in the plane problem. We will indicate the formulae determining this quantity and draw certain conclusions. As in Sect. 3.11, we focus on a period [A, Ai] around the time when the contact area is a minimum. The values assigned to zdi, ZI2, at the end of Sect. 3.11 for a load given by (3.11.58) also apply here. During the decreasing phase, before Iq, D(t) is determined by (5.2.17), which becomes [see (2.4.18)]  

In the case of the standard linear model, 77 (0 is given by (3.11.17), while 

Let us now agree that t to, and write out the two equations more explicitly. Equation (5.2.37) becomes, on transforming the integral to [/, + 

Note that we deduce from (5.2.40), in a manner similar to the derivation of (3.11.55), that 

---

Equations (3.77) and (3.78) give the solutions to the steady-state problem, given... 

Find the analytical solution to the steady-state problem in Example 4.2. 

A radically different approach to the steady-state problem was investigated by Hsing (H6). In this approach the steady-state flow problem was formulated as the following constrained minimization problem ... 

The steady-state problem yields a system of simultaneous linear algebraic equations that can be solved by Gaussian elimination and back substitution. I shall turn now to calculating the time evolution of this system, starting from a phosphate distribution that is not in steady state. In this calculation, assume that the phosphate concentration is initially the same in all reservoirs and equal to the value in river water, 10 I 3 mole P/m3. How do the concentrations evolve from this starting value to the steady-state values just calculated ... 

If available, use Twopnt or other user-oriented boundary-value software to solve the steady-state problem. 

In the steady stagnation-flow formulation the thermodymanic pressure may be assumed to be constant and treated as a specified parameter. The small pressure variations in the axial direction, which may be determined from the axial momentum equaiton, become decoupled from the system of governing equations (Section 6.2). The small radial pressure variations associated with the pressure-curvature eigenvalue A are also presumed to be negligible. While this formulation works very well for the steady-state problem, it can lead to significant numerical difficulties in the transient case. A compressible formulation that retains the pressure as a dependent variable (not a fixed parameter) relieves the problem [323],... 

I. The steady-state problem. Diffusion and reaction, Chem. Eng. Sci. 1957, 6, 262-268. 

We apply the finite difference scheme to the first integral of the above equation and treat the second integral in the same way as we did in the steady-state problem of the previous section. This results in... 

The resultant single differential equation can then readily be transformed into Eq. (18). Thus, it becomes clear that the O Toole (1965) formulation of the steady-state problem is exactly equivalent to that of Smith and Ewart (1948), notwithstanding that the approach to the problem is somewhat different. 

The periodic problem (85) differs from the steady state problem (59) by one term, namely i Q [R, Z]. This term change the solution to give the amplimde and phase lag of temperature oscillation in a point with coordinates R and Z. 

Static bifurcation can be studied by means of singular perturbation methodp 1.6. taking c - 0. The problem of finding static bifurcation points turns to the steady state problem similar to the CSTR where is expressed by (16)-From this and from the monotonicity of the fvmction G follows that for a given value of bifurcation parameter Da the maximal number of static bifurcation points corresponds to that of the systems with = 1 and Da <0 Da >. 

To illustrate the features of our proposed algorithm, we apply it to the case of a tubular reactor with axial dispersion, where an elementary first order irreversible exothermic reaction takes place A —[8]. The steady state problem is described two nonlinear partial differential equations (PDEs), which in dimensionless form are ... 

A sketch of the steady-state problem that we will consider is shown in Fig. 3 13. The tube radius is denoted as a, and we utilize the standard cylindrical coordinates (r, (f>, z).18 We assume that the wall of the tube is insulated so that the heat flux is zero for positions z < 0. On the other hand, beginning at z = 0 and for all z > 0 there is a constant positive heat flux q through the wall of the tube. A reasonable approximation to this condition can be realized if the tube is uniformly wrapped with a heat tape or wire resistance heater beginning at z = 0. It may be noted that the problem would be identical from a mathematical point of view if there were a negative heat flux prescribed so that the temperature of the fluid in the tube decreased rather than increased for z > 0. We assume that the entry temperature of the fluid into the heated portion of the tube is Qo and that viscous dissipation can be neglected. 

One point that has not been emphasized is that all of the preceding analysis and discussion pertains only to the steady-state problem. From this type of analysis, we cannot deduce anything about the stability of the spherical (Hadamard Rybczynski) shape. In particular, if a drop or bubble is initially nonspherical or is perturbed to a nonspherical shape, we cannot ascertain whether the drop will evolve toward a steady, spherical shape. The answer to this question requires additional analysis that is not given here. The result of this analysis26 is that the spherical shape is stable to infinitesimal perturbations of shape for all finite capillary numbers but is unstable in the limit Ca = oo (y = 0). In the latter case, a drop that is initially elongated in the direction of motion is predicted to develop a tail. A drop that is initially flattened in the direction of motion, on the other hand, is predicted to develop an indentation at the rear. Further analysis is required to determine whether the magnitude of the shape perturbation is a factor in the stability of the spherical shape for arbitrary, finite Ca.21 Again, the details are not presented here. The result is that finite deformation can lead to instability even for finite Ca. Once unstable, the drop behavior for finite Ca is qualitatively similar to that predicted for infinitesimal perturbations of shape at Ca = oo that is, oblate drops form an indentation at the rear, and prolate drops form a tail. 

Statement of the Blasius problem. We consider the steady-state problem on the longitudinal zero-pressure-gradient flow (VP = 0) past a half-infinite flat plate (0 < X < oo). We assume that the coordinates X and Y are directed along the plate and transverse to the plate, respectively, and the origin is placed at the front edge of the plate. The velocity of the incoming flow is U. ... 

Suppose that on a vertical wall whose temperature is constant and equal to Ts, stagnant dry saturated vapor is condensing. Let us consider the steady-state problem under the assumption that we have laminar waveless flow in the condensate film. According to [200], we make the following assumptions the film motion is determined by gravity and viscosity forces the heat transfer is only across the film due to heat conduction there is no dynamic interaction between the liquid and vapor phases the temperature on the outer surface of the condensate film is constant and equal to the saturation temperature Tg the physical parameters of the condensate are independent of temperature and the vapor density is small compared with the condensate density the surface tension on the free surface of the film does not affect the flow. 

R. Aris, Shape Factors for Irregular Particles. 1. The Steady-State Problem. Diffusion and Re-... 

First solve the steady state problem for the heat exchange area A for normal operation T= 358 K. 

Interfacial Limitation The solution of the steady-state problem for only interfacial limitation is given by (12.91). The aqueous-phase reaction rate is then... 

Finally, we often assume that the diffusivity, thermal conductivity and partial molar enthalpies are independent of temperature and composition to produce the following coupled mass and energy balances for the steady-state problem... 

Since it was not within our ability to solve the time-dependent equahon, we naturally solved the steady-state problem so that the accumulation term went to zero, which left only the spatial derivative ... 

Fig. 4.6. Distribution of ice-like clusters in water (a) at a temperature above o °C, (6) at a temperature below o °C if an equilibrium distribution is assumed, (c) at a temperature below o °C for the steady-state problem.

In his first studies of the Fourier equation for temperature distribution in a reactive system, Frank-Kamenetskii restricted his attention to three shapes, the infinite slab, the infinite cylinder, and the sphere. For these three geometries (class A) the Laplacian operator can be expressed in terms of a single co-ordinate, and the steady-state problem is reduced to solving the ordinary differential equation... 

For most actual problems the source of neutrons is at very high velocities compared to thermal velocities. The convergence of the eigenfunction solution (11) or (12) will be extremely slow under these conditions, and the direct solution in terms of eigenfunctions will not be very useful. For the steady state problem, the solution for N(v) in the region far above thermal velocities, but far below source velocities, is well known [10] to go as Ijv, It would be interesting to show that the solution (12) has this property. We recall that the singular solutions to Equation (10) have this asymptotic... 

The eigenvalue approach to the steady state problem has the advantage that the solution (12) explicitly exhibits the dependence of the neutron velocity distribution on absorber concentration. In the limit of weak absorption, a Ai < A2 , the solution becomes... 

The study undertaken in this section is motivated by a so-called method of harmonics associated with criticality determinations in a thermal reactor. The harmonics method for the steady state problem has been studied by Goertzel and Garabedian [13] and Edlund and Noderer [14]. In this method, functions representing neutron fluxes, absorptions, productions, and therma-lizations are expanded in spatial eigenfunctions which vanish at the outer boundary of the total reactor system. In this section the domain of applicability of the method is extended to the kinetics problem described above. 

Aris R. On shape factors for irregular particles. 1. The steady state problem of diffusion and reaction. Chemical Engineering Science 1957 6(6) 262-268. 

As has been mentioned in Sect. 7.3, the continuous injection-molding operation results in a cyclic heat transfer behavior in the mold, after a short transient period. The cycle-averaged temperature can be represented by a steady state heat conduction equation, i.e., Eq. 7.10. The mold cooling analysis can be greatly simplihed by solving the steady state problem. The boundary integral equation of ( 7.10) is... 

The solution of the steady state problem described above was performed using the commercial software COMSOL Multiphysics v 3.5a. The numerical technique used by that software is the Finite Element Method (FEM). The shape functions, chosen for the simulation, are Lagrange quadratic shape functions. 

The main idea in the development of the algorithm for the nonlinear back-off synthesis method is to use the steady state formulation in order to generate promising control structures that are then evaluated under dynamic conditions. Since the dynamics of the plant can only further restrict the feasible space, when compared to the steady state, dynamic economics caimot be better than the steady state economics (see also Fig. 1). Any structure that is feasible under dynamic conditions (ie[0,T/]) is also feasible at steady state (t=0) since the latter is a subset of the former. However, the converse is not true. That is, a structure that is feasible at steady state can be dynamically infeasible. Those structures that correspond to steady state feasible but dynamically infeasible solutions have to be excluded fi-om the feasible space of the steady state problem in order to obtain the feasible space of the dynamic... 

---