<< Chapter < Page | Chapter >> Page > |
The pressure relative to some datum can be determined by the square of the magnitude of velocity. This is easily calculated from the complex potential.
There are some theorems that facilitate the integration of pressure around bodies in the complex plane, but they will not be discussed here. The pressure and tangential velocity profiles for the inviscid flow around an object are needed for calculation of the viscous flow in the boundary layer between the solid boundary and the inviscid outer flow.
% Calculate pressure distribution from square of velocity
(your code here to calculate pressure field)pcolor(x,y,p)
colormap(hot)shading flat
axis image
The conservation equations for material, momentum, and energy reduce to first order PDE in the absence of diffusivity, dispersion, viscosity, and heat conduction. In thin films, viscosity may be a dominant effect in the velocity profile normal to the surface but the continuity equation integrated over the film thickness will have only first order spatial derivatives unless the effects of interfacial curvature become important. In one dimension, the system of first order partial differential equations can be calculated by the method of characteristics [A. Jeffrey (1976), H.-K. Rhee, R. Aris, and N. R. Amundson (1986, 1987)]. Here we will only consider the case of a single dependent variable with constant initial and boundary conditions. Denote the dependent variable as S and the independent variables as and . The differential equation with its initial and boundary conditions are as follows.
The dependent variable can be normalized such that the initial condition is equal to zero and the boundary condition is equal to unity. Thus, henceforth it is assumed such a transformation has been made. The dependent variable will be called 'saturation' and the flux called 'fractional flow' to use the nomenclature for multiphase flow in porous media. However, the dependent variable could be film thickness in film drainage or height of a free surface as in water waves. The PDE, IC, and BC are rewritten as follows.
The differential, is easily calculated since there is only one independent saturation. If there were three or more phases this differential would be a Jacobian matrix. The locus of constant saturation will be sought by taking the total differential of .
Notification Switch
Would you like to follow the 'Transport phenomena' conversation and receive update notifications?