Dynamic Behaviour of Process Systems.pptVIP

Modelling Simulation of Chemical Engineering Systems Department of Chemical Engineering King Saud University LECTURE #9Numerical Solution of Boundary-Value Differential Equations BVP Examples Diffusion with chemical reaction in a catalyst. Diffusion or condition in a tubular reactor Heat mass transfer in boundary layer problems. Disecretization of nonlinear PDEs Initial Boundary Conditions Numerical Methods BVP Formulation Canonical Form Original Form Shooting Method An initial-value technique Use successive approximation to solve for the missing initial conditions. It can utilize highly efficient initial value procedure. Shooting can be tried from either end of the BVP. Shooting Method Algorithm Problem: Guess initial condition for y2 : Integrate the system forward, calculated y2 at xf is Define the function Solution : Example MATLAB Code Finite Difference Method Replace the derivative with finite difference approximation ODEs become algebraic equations Example of FD method Example: solve the BVP using FD technique MATLAB command BVP4C Solve two-point boundary value problems for ODEs by collocation.  SOL = BVP4C(ODEFUN,BCFUN,SOLINIT) integrates a system of ordinary differential equations of the form y = f(x,y) on the interval [a,b], Example solinit = bvpinit([0 1 2 3 4],[1 0]); sol = bvp4c(@twoode,@twobc,solinit); solve a BVP on the interval [0,4] with differential equations and boundary conditions computed by functions twoode and twobc, respectively. This example uses [0 1 2 3 4] as an initial mesh, and [1 0] as an initial approximation of the solution components at the mesh points. xint = linspace(0,4); yint = deval(sol,xint); evaluate the solution at 100 equally spaced points in [0 4]. The first component of the solution is then plotted with plot(xint,yint(1,:)); HW Solve the problem of flow of a Newtonian flow fluid in a pipe described