放射性废料的处理问题

放射性废料的处理问题美国原子能委员会以往处理浓缩的放射性废料的方法，一直是把它们装入密封的圆桶里，然后扔到水深为90多米的海底。生态学家和科学家们表示担心，怕圆桶下沉到海底时与海底碰撞而发生破裂，从而造成核污染。原子能委员会分辨说这是不可能的。为此工程师们进行了碰撞实验。发现当圆桶下沉速度超过12.2m/s与海底相撞时，圆桶就可能发生碰裂。这样为避免圆桶碰裂，需要计算一下圆桶沉到海底时速度是多少?这时已知圆桶重量为239.46kg，体积为0.2058m3，海水密度为1035.71kg/m3，如果圆桶速度小于12.2m/s就说明这种方法是安全可靠的，否则就要禁止使用这种方法来处理放射性废料。假设水的阻力与速度大小成正比例，其正比例常数6.0k。现要求建立合理的数学模型，解决如下实际问题：1.判断这种处理废料的方法是否合理?2.一般情况下，v大，k也大；v小，k也小。当v很大时，常用kv来代替k，那么这时速度与时间关系如何?并求出当速度不超过12.2m/s，圆桶的运动时间和位移应不超过多少?(k的值仍设为0.6)3.1问题分析与建立模型的建立(1)首先要找出圆桶的运动规律，由于圆桶在运动过程中受到本身的重力以及水的浮力F和水的阻力f的作用，所以根据牛顿运动定律得到以下的方程:fFGF合(1)又因为gVFmgGdtsdmdtdvmmaF,,22合以及dtdskkvf，可得到圆桶的位移和速度分别满足下面的微分方程:)3()2(22kvgVmgdtdvmdtdskgVmgdtsdm(2)由题设这时圆桶受到的阻力应改为22)(dtdskkvf，类似上面，可得到这时圆桶的速度应满足如下的微分方程:2kvgVmgdtdvm(4)3.2模型求解(1)根据方程(2)和(3)，加上初始条件0,0000tttsvdtds以及题设的初始数据。通过Matlab就可以求出圆桶的位移合速度的方程。Matlab源程序如下:m=239.46;w=0.2058;g=9.8;p=1035.71;k=0.6;dsolve('m*D2s=m*g-p*g*w-k*Ds','t')C3/exp((k*t)/m)-(C2-g*m*t+g*p*t*w)/k-(g*m*(m-p*w))/k^2]]],[},0]0[,0]0[],[****][*[{[;6.0;71.1035;8.9;2058.0;46.239''''ttssstskwgpgmtsmDsolveChopkpgwmans=(2362544856843450433*t)/5497558138880000-(5*C2)/3+C3/exp((10*t)/3991)-9428916523662210678103/54975581388800000]],[},0]0[],[****][*[{'ttvvtvkwgpgmtvmDSolve得出位移方程为:tetts00250564.0171511744.429171511)((5)速度方程为tetv00250564.0744.429744.429)((6)通过方程(5)及mts90)(，利用下面Mathematica程序:}]13,{,171511]00250564.0[/171511744.42990[ttExptFindRoot求出时间st002.13，再把它带入方程(6)求出圆桶的速度为smv/7729.13.显然此圆桶的速度已超过sm/2.12，可以得出这种处理废料的方法不合理。因此，美国原子能委员会已经禁止用这种方法来处理放射性废料。(2)根据方程(4)，加上初始条件00tv，利用下面的Mathematica程序:]],[},0]0[,2])^[(****][*[{'ttvvtvkVgpgmtvmDSolve求出圆桶的速度:).()(tmFkTanhkFtv若把题设中的kF,(仍为0.6的话)和m的值带入上式，可得到圆桶的速度为)0519436.0(7303.20)(tTanhtv这时若该速度要小于sm/2.12，那么利用Mathematica可得圆桶的运动时间就不能超过s13，位移不能超过84.8m。对应的Mathematica源程序:}]12,{,2.12]0519426.0[*7303.20[ttTanhFindRoot(求时间)}]13,0,{],0519426.0[*7303.20[ttTanhIntegrate(求位移)从这个模型，我们也可以得到原来处理核废料的方法是不合理的.3.3结果分析(1)由于在实际中k与v的关系很难确定，所以上面的模型有它的局限性，而且对不同的介质比如在水中与在空气中k与v的关系也不同.如果按照(1)的假设k为常数的话，那么水中的这个k就比在空气中对应的k要大一些.在一般情况下，k应是v的函数，即k=k(v)，至于是什么样的函数，这个问题至今还没有能解决.(2)这个模型还可以推广到其他方面，比如说一个物体从高空落向地面的道理也是一样的.尽管物体越高，落到地面的速度越大，但决不会无限大.3.4相关的思考题(1)若在其他介质如在空气中物体的运行速度与它所受的空气阻力仍成正比，那么这个正比例系数与在水中的正比例系数是否一样?(2)把所建立的数学模型应用到求火车的速度问题上，并说明为何火车的速度不能无限增大?