16.12.2012 Views

Computer Algebra Recipes

Computer Algebra Recipes

Computer Algebra Recipes

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

266 CHAPTER 6. LINEAR PDE MODELS. PART 2<br />

6.2.3 Radioactive Contamination<br />

The unexamined life is not worth living.<br />

Socrates, Greek philosopher (470{399 BC)<br />

When coauthor Richard was a student, he worked one summer as a chemistry<br />

lab assistant at a uranium mine on Great Bear Lake, which straddles the<br />

Arctic Circle in Northern Canada. Because it was a summer job, there was a<br />

period of several weeks when the sun never set, so evening baseball games were<br />

never called o® because of darkness. The baseball diamond was located on the<br />

leveled mine tailings, the only relatively °at spot in the small northern min-<br />

0<br />

x<br />

Figure 6.1: The radioactive disposal site.<br />

ing town of Port Radium, which is perched on almost treeless primordial rock<br />

formations. The author has always wondered what long-term health problems<br />

eventually arose among the permanent workers because of working in the mine<br />

and playing baseball on the radioactive mine tailings. Motivated by these reminiscences,<br />

we could not resist including a related example involving radioactive<br />

contamination.<br />

A radioactive gas is di®using at a steady rate into the atmosphere from<br />

a leveled contaminated disposal site. The ground and the atmosphere will<br />

be taken to be semi-in¯nite media with X = 0 at the boundary as shown in<br />

Figure 6.1. The concentration C(X; T) of radioactive gas in the atmosphere<br />

obeys the concentration equation (a modi¯ed di®usion equation)<br />

@C(X; T)<br />

@T<br />

= d @2 C(X; T)<br />

@X 2 ¡ ¸C(X; T); (6.8)<br />

where d is the di®usion constant and ¸ is the decay rate of the radioactive gas.<br />

The boundary condition at X =0isgivenbyFick's law,<br />

¡d @C(0;T)<br />

= K; (6.9)<br />

@X

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!