This preview shows pages 1–3. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: ODE  Medical Akzo Nobel problem II41 4 Medical Akzo Nobel problem 4.1 General information The problem consists of 2 partial di erential equations. Semidiscretization of this system yields a sti ODE. The parallelIVPalgorithm group of CWI contributed this problem to the test set in collaboration with R. van der Hout from Akzo Nobel Central Research. The software part of the problem is in the le medakzo.f available at [ MI03 ]. 4.2 Mathematical description of the problem The problem is of the form d y d t = f ( t; y ) ; y (0) = g; (II.4.1) with y 2 IR 2 N ; t 20 : Here, the integer N is a usersupplied parameter. The function f is given by f 2 j 1 = j y 2 j +1 y 2 j 3 2 + j y 2 j 3 2 y 2 j 1 + y 2 j +1 ( ) 2 k y 2 j 1 y 2 j ; f 2 j = k y 2 j y 2 j 1 ; where j = 2( j 1) 3 c 2 ; j = ( j 1) 4 c 2 : Here, j ranges from 1 to N , = 1 N , y 1 ( t ) = ( t ), y 2 N +1 = y 2 N 1 and g 2 IR 2 N is given by g = (0 ; v ; ; v ; : : : ; ; v ) T : The function is given by ( t ) = 2 for t 2 (0 ; 5] ; for t 2 (5 ; 20] : which means that f undergoes a discontinuity in time at t = 5. Suitable values for the parameters k , v and c are 100, 1 and 4, respectively. 4.3 Origin of the problem The Akzo Nobel research laboratories formulated this problem in their study of the penetration of radiolabeled antibodies into a tissue that has been infected by a tumor [ Hou94 ]. This study was carried out for diagnostic as well as therapeutic purposes. Let us consider a reaction di usion system in one spatial dimension: @u @t = @ 2 u @x 2 kuv; (II.4.2) @v @t = kuv; (II.4.3) II42 ODE  Medical Akzo Nobel problem which originates from the chemical reaction A + B k ! C: Here A , the radiolabeled antibody, reacts with substrate B , the tissue with the tumor, and k denotes the rate constant.the rate constant....
View Full
Document
 Spring '09
 Linda

Click to edit the document details