Exact solution of the diffusionconvection equation in 在对流扩散方程精确解.pptVIP

NS, LPPD Exact solution of the diffusion-convection equation in cylindrical geometry Oleksandr Ivanchenko, Nikhil Sindhwani, and Andreas A. Linninger Laboratory for Product and Process Design, M/C 063 University of Illinois at Chicago Chicago, IL 60607-7000, USA * NS, LPPD LPPD seminar: 2nd September 2009 Problem Formulation NS, LPPD * C(r,t)=CS/C0 q= Qin/h, or the radial conductive velocity field. NS, LPPD * Boundary conditions: Convection-diffusion equation in radial co-ordinates: Initial conditions: General solution, Separation of variables. NS, LPPD * Substituting this in the convection-diffusion equation. Temporal solution. λi used because there could be more than one λ that satisfy the equation above. λi has to be a positive number for the solution to be stable. NS, LPPD * Radial solution: where, μ=1-V0/D, also, By using Forbenius power series method of solution, This can also be written by using the definition of a Bessel function of the First Kind and order λ. NS, LPPD * Using , we can find the roots of the bessel function equation at r=L. Where, si, are the roots of the bessel function. Now, λ is known and we can write the final analytical solution NS, LPPD * By using the initial condition, we can find out Ai, by fourier-bessel decomposition: NOTE: This solution fails when f(r)=0. This is the case when the domain is empty initially. Advanced solution: by decomposing steady state and dynamic parts. NS, LPPD * Steady state part and dynamic part of the solution are separated. Boundary Conditions: Initial Condition: NS, LPPD * Steady state solution: At t=0: This gives the initial condition for dynamic part of the solution, which can be replaced in the general solution to give: NS, LPPD * Final Solution: NS, LPPD * Solution for D=0.01 and γ=1 NS, LPPD * N=20 N=60 Effect of number of terms used, more terms, smoother solution. NS, LPPD * Solution in Non-dimensional form: τ=Dt/L2 , ρ=r/L, P=Vo/D Steady state solution NS, LPPD * Effec