考虑包气带通过一个孔道经过孔隙水的扩散向潜水传送组分的传质过程。与瞬时源不同的是,这时扩散组分通过在孔隙水表层形成一个稳定的饱和层,在气体相平衡溶解的作用下成为浓度恒定且持续的源,如图8—6所示。
图8—3 扩散距离z=2cm处浓度的历时曲线
图8—4 扩散过程不同时刻浓度的分布曲线
图8—5 瞬时源扩散方程曲面图
a—单向扩散曲面图;b—双向扩散曲面图
连续扩散源可表示成无限多的瞬时源(单个脉冲)连续组成的。扩散系统空间位置P点扩散组分浓度是由无限多个瞬时源的连续累加形成的,其浓度可用叠加的形式表达。
图8—6 连续恒定源单向一维分子扩散
连续恒定源扩散问题的初始条件为:c=0(z>0,t≥0);边界条件:无限远边界条件:c(∞,t)=0(t>0);若以恒定浓度常量c0表达dW,使得瞬时源强度dW=c0dξ(dξ为微小的距离),则用δ函数δ(ξ)表示的初始边界条件为c(0,t)=dW·δ(ξ)=c0dξ·δ(ξ)(z=0,t≥0),δ(ξ)的量纲为[L—1]。
坐标ξ(为负值)的一个瞬时源扩散作用到P(z)点,即为一个dW瞬时源的扩散传质方程在z位置产生的浓度值,依据式(8—7a)或式(8—7b)积分,得式(8—9)中的z—ξ为瞬时源到P(z)点的位置。
水文地球化学基础
连续恒定源单向扩散至z位置产生的总浓度值c(z,t),等于持续恒定源中每个微分量扩散浓度值dc(z,t)的叠加,符合叠加原理,表达为c(z,t)=
水文地球化学基础
将参变量改写为ξ=z—ξ,则积分上下限ξ=[—∞ 0]变为ξ=[∞ z],dξ=—dξ,式(8—11)可仍以可取代号作变量写为
水文地球化学基础
式(8—10)或式(8—11)右边的积分可用误差函数erf或余误差函数erfc表示。
由误差函数的定义:
水文地球化学基础
误差函数常用性质有:erf(0)=0,erf(—x)+erf(x)=0,erf(∞)=1。
余误差函数erfc的定义:
erfc(x)=1—erf(x)
在应用误差函数作数量计算时,可利用软件计算其函数值,一般精度的计算可查图8—7得到误差函值。
对照erf(x)定义,取z=
图8—7 误差函数曲线
利用误差函数的性质处理后得到的式(8—12),是以余误差函数erfc()表达的连续恒定源单向扩散问题解:
水文地球化学基础
依据方程式(8—12),式中的参数c0=10mg/L,D=1.3×10—4m2/d,绘制图8—8。式中
式(8—10)~式(8—12)是持续恒定源单向分子扩散模型的解析解形式,图8—8所示曲面是解的三维图形解形式。
图8—8 持续源一维扩散浓度变化曲面
参数取值:c0=10mg/L,D=1.3×10—4m2/d
作图范围:0<z<0.5,0<t<300
为了直观表现扩散组分状态的动态变化,可以利用非稳态问题的解,作连续多个不同时刻t下的z—c曲线,制作动画,即为动画形式的解。因此,解析解、数值解是最常用的两种形式,有时其他的表达形式可能更适合应用的需要。
【例题8—3】某水井内径25cm,渗流涌入1.0m深的井水,井水温度10℃,表层水饱和溶解氧浓度可持续达10mg/L,井水中无产氧和耗氧过程,水井附近含水层地下水中的溶解氧浓度为0。试计算300d的平均扩散距离
解:t=300d,D=1.3×10—4m2/d,由式(8—5)得
按300d时平均扩散距离为0.27m,可认为扩散浓度波前锋可达到
依据式(8—8),在Microsoft Excel表单中可利用Excel内置函数erfc计算。数表解见表8—3。
表8—3 解的数据表
再将数值结果作图(图8—9),以利于观察浓度变化趋势和得到更多位置上的浓度值。
图8—9 持续源一维单向扩散浓度变化曲线
图8—8为本题单向持续恒定源分子扩散模型的图形解。图8—9为图8—8中的一个剖面线,是在一定时间(t=300 d)时持续源一维单向扩散浓度变化曲线。
扩散源扩散组分的来源,有多种形式。如前所述的瞬时源是最基本的源,其他的源可以用瞬时源表示。除了无限连续源、无限持续源以外,还存在有限连续或持续源;持续源的源强度可以如例8—3为恒定源,也可是强度随时间变化的源。
本文地址: http://www.goggeous.com/f/1/57324
文章来源:天狐定制
版权声明:除非特别标注,否则均为本站原创文章,转载时请以链接形式注明文章出处。
2025-01-08职业培训
2025-01-08职业培训
2025-01-08职业培训
2025-01-08职业培训
2025-01-08职业培训
2025-01-08职业培训
2025-01-08职业培训
2025-01-08职业培训
2025-01-08职业培训
2025-01-08职业培训
2024-11-26 19:19:00职业培训
2024-11-26 19:18:59职业培训
2024-11-26 19:18:49职业培训
2024-11-26 19:18:48职业培训
2024-11-26 19:18:47职业培训
2024-11-26 19:18:46职业培训
2024-11-26 19:18:46职业培训
2024-11-26 19:18:45职业培训
2024-11-26 19:18:44职业培训
2024-11-26 19:18:43职业培训
2025-01-01 10:26职业培训
2024-12-11 10:22职业培训
2024-12-22 03:07职业培训
2025-01-02 23:57职业培训
2024-12-14 10:51职业培训
2024-12-28 01:04职业培训
2024-12-01 22:04职业培训
2025-01-01 09:02职业培训
2024-12-11 07:43职业培训
2024-12-22 00:27职业培训
扫码二维码
获取最新动态