- 3
- 0
- 约2.26千字
- 约 3页
- 2023-12-04 发布于辽宁
- 举报
一个圆柱形的容器,内装350升均匀混合的盐水溶液。如果纯水以每秒14 升的速度从容器顶部流入,同时,容器内的混合的盐水以每秒10.5升的速度从 容器底部流出。开始时,容器内盐的含量为7千克。求经过时间t后容器内盐的 含量。
一.问题分析
已知:水的密度为1kg/L,盐溶解度为36g(在假设(1)中)
由题意可以知道,此题中的t时间时容器内的盐含量的问题即为关于t的微 分方程的求解问题。在流入到流出的过程中,由于混合在水中的盐含量是不同的, 所以溶解于水中的盐的量每一时刻都是不同的。而所求的这个瞬时量即为微分方 程的解。而由题意知,可计算出7kg盐所需要的溶剂为194L水。因此,由混合 液体积即可知开始时刻的7kg盐是完全溶于水中的,并且没有饱和。所以,整个 过程为食盐水被再次稀释的过程,则不会出现有盐析出现象。即加入纯水盐水就 被稀释,同时向外流出含有盐的混合液。因此,整体变量只与时间t有关。则可 根据时间t列出微分方程求解。
二假设
1) 由资料知,温度对盐在水中的溶解度变化影响不大,我们可以假设盐溶解度 恒为室温条件下的溶解度为36g,又由条件知该过程为稀释过程,则不存在盐析 出结品的过程,即不考虑结晶度问题。
2) 任意时刻容器内混合的、流出的盐水都均匀(即任意时刻容器内及流出的盐 水密度恒定包含初始状态350升混合盐水)
3) 水流入及盐水流出的速度均为匀速的。
4) 假设注水时间为t、t时刻时容器内含盐量为S(t)、纯水流入容器的速度为v1, 混合液流出的速度为v2。
模型建立
计算过程如下:
1) t时刻容器内含盐量可表示为S(t)
2) 纯水以v1=14L/s的速度流入容器,混合液以v2=10.5L/s的速度流出容器,
则根据溶解度公式:溶解度=巴溶质*i00g,可推出7kg食盐是完全溶解在水中 m溶剂
的
r时刻混合液中盐含量t时刻混合液量3)根据假设1)、
r时刻混合液中盐含量t时刻混合液量
* △ t时刻流出的混合液量=△ t时刻混合液中盐的变化量
(假设其中At足够小)
4)根据上述比例式即可得到相应的方程为: 一*10.5 * At = P(r)- (t+ At)
350 + 3.5 * t
而使At - 0可得微分方程:* = - 10.5 * P(t)。
dt 350 + 3.5 * t
5)求解该微分方程当0 p(t) 7时的曲线图像
模型求解与检验
设溶液中盐的浓度为C,时间变量为t
当 t=0 时,C = C0, t=t 时,C = C
当由t变化到t+dt, C变化到C + dC
由物料守恒,
输入一输出=积累+反应
0—VoCdt=(350+(Vin-Vo)t)*C-(350+(Vin-Vo)*(t+dt))*(C+dC)
化简此方程,并从上面两个条件积分,得到
C=C0/(1+0.01t)人4
从而,含量为:
M=C*(350+3.5t)
数值解
dsolve(Dy=-10.5*y/(350+3.5*x) ,x)
ans =
C2/(x + 100)人3
又由于初始x=0,时溶液中盐含量为7kg,所以可以计算出C2= 7000000
图形解
则带入C2= 7000000则可求出解的图像:
x=linspace(0,300,10000);
y= 7000000./ (x + 100).人3;
plot(x,y) ,xlabel(时间’),ylabel(混合中盐含量),title(混合液中盐含 量曲线)
将仁0代入y (t)可以得到y=7,与初始状态的含盐量为7相符。显然可以得到 y (t)在区间[0,8)上函数是连续的。使用MATLAB求函数的导数可以看出, y(t)在区间[0,8)上导数恒小于0,所以y(t)在[0,8)为减函数。使用MATLAB 求函数右极限,可以得出t越大时y是接近于0的,但是不等于0,与真实情况 相符。
六.模型评估
模型模型模型模型评估评估评估评估::::该模型在水流速度较小的情况 下能较为真实的反应容器内的含盐量,但是水流速度如果较大的时,达到盐在不 同浓度的水中的扩散至均匀的时间不能忽略的时候将不能很好的反应容器内的 含盐量。
结论及应用
溶液混合问题的更一般的提法是□设有一容器装有某种质量浓度的溶液,以流量
1V注入质量浓度为1C的溶液口指同一种类溶液,只是质量浓度不同口,假定溶 液立即被搅匀,并以2V的流量流出这种混合溶液,试建立容器中质量浓度与时间 的数学模型.该模型不仅适用于液体的混合,而且还适用于讨论气体的混合.
求盐水混合物含盐量的问题不仅仅局限于数学问题的求解,在实际中也有广 泛的应用:
在农田的施肥过程中,施化肥的过程相当于注入较高浓度的盐水,而在农田 的实际管理过程中,灌溉用水得从农田流过,因此,农田里化肥水浓度会发生改 变,而农作物
原创力文档

文档评论(0)