数学建模社区-数学中国

标题: 【转】过桥米线里的肉片热传导模型 [打印本页]
作者: 我要吃章鱼丸子    时间: 2016-3-31 14:22
标题: 【转】过桥米线里的肉片热传导模型
过桥米线里的肉片热传导模型
; ]  p- [0 W4 e/ y$ F' F' q/ G  N7 o
作者：学习部  来源：桂工数模  查看：198  评论：0
1 z0 K0 s8 \( A. p, z5 u  _6 Z" l8 D+ x( D$ S

8 g! c0 X- `8 e! R9 ]7 q1 c5 d; |7 ]+ r! k3 W) w
$ ]. ?& U4 h4 T) `4 X$ t7 B! ~
$ y) y, N' D; O9 u

. K/ E6 f$ m8 P" W, ^8 L( p' k" s. H( S5 c, _9 n
     近来，网上关于云南过桥米线和旋毛虫的问题有很多的讨论。咱也凑个热闹，说说过桥米线里面的物理问题。讨论的一个关键，在于肉片放到热汤里面之后，多长时间之内肉片中心会达到什么样的温度，能否达到灭杀旋毛虫幼虫所需的温度。
" ?  a% r: y. P/ K" D. A. \+ E
5 V* A* `! V0 A6 D: u
一维热传导：建模与求解- o3 n8 t# X" m1 c4 Y2 s

( h1 A: y  w+ G    首先让我们来简化一下这个问题。由于一般来说，热汤的量远多于肉片的量，可近似认为加入肉片不会造成汤温变化。也就是说，肉片的大小是没有关系的，肉片的厚度才是更关键的因素。而若非肉片叠在一起，肉的两边都是和汤接触的，因此我们可以简单地把这个问题考虑成一个一维的热传导问题。0 C, r. c9 u  w; [

6 {" R! b" H9 l; k

/ L6 I" H% j) r: q4 e+ b2 d( n) y% w0 D( p
    诚然，肉片的厚度可能会不均匀，肉片在汤里面可能是有弯折的。这些因素对肉片中心的温度都会造成影响。但是为了得到一个初步的了解，我们可以先不考虑这些因素。这样，我们就得到了一个单纯的一维热传导问题。在肉片放入汤中的时候，整个体系的温度如下图所示：
) c4 K7 T$ a7 H# l7 f' W# A: R& v. s, a" p4 K8 m

0 ~& \) \, n* @9 e0 m1 ^/ I7 f( A% _/ f6 V
    这里我们把坐标x的原点定在肉片的中心，设肉片的厚度为L，则肉片从x=-L/2到L/2。假设汤的温度恒定为Tsoup不变，肉的温度在开始 (t=0) 时为Tmeat，这些就是这个问题里面的边界条件。现在我们需要的，就是利用物理和数学的方法去在这个体系里面求解一维热传导方程[1] ：
/ A: B" R9 a, g2 G5 w: Z
% q" Q: s' x' e/ k% P* H
: Z7 [- l, ~/ \  d
& _" s( |! y3 {( B& A  R$ W
    其中，α是肉片里面的热扩散系数，数值越大说明热在物体里面越容易传播，我们可以用水的值来代替，α=1.4×10-7m2/s [2] 。3 [& u$ y$ U& v( T

, A+ ]. [( l7 H- b8 X    具体的求解需要利用分离变量法和余弦函数的正交展开，仔细推导过程见文章最末，在此先略过不提。6 C$ n: [7 O! |" \
- o( x, r' q1 w5 d
    对于任意的Tsoup和Tmeat我们可以得到如下的解：4 n4 z% n" S. _1 }

. P; M& s; D, P2 O' I
[点击查看原图]
9 m) Y7 V! ?' v$ p% i4 k; n- S4 r6 Q6 l: ^
9 T4 F0 F0 A1 g' Y8 J. }

# C  \- v& O5 |. A7 x7 r- g5 f( H薄肉片，导热快
# r2 T4 v+ |; ]  t4 _) f* e6 m( L: ]4 E: O6 E) n( S* c" l
    根据这个解，如果知道汤和肉片的初始温度以及肉的厚度，我们就可以给出肉内部温度随时间变化的曲线。如果以汤的最低温度 80℃，肉的初始温度是 0℃ （不妨就设肉品刚刚解冻），肉的厚度是1毫米为初始条件，那么我们可以画出当时间 t 为 0 秒， 1 秒， 2 秒， 3 秒的时候肉内部的温度的变化图：
; T6 X7 g: r1 P, E2 `& t7 P8 Z; Q6 }* @9 X! y0 m% s+ A( A& V
1 `, c+ j) c/ L. {2 B/ R! C
' W5 P0 s% {4 @
    根据上图，3 秒之后肉片中心温度大约是 78℃，基本接近汤的温度了。0 D) q" j: ?& L! R  L- N( I; Y

' a0 B* [- s, Q    另一方面，我们也可以画出肉片中心的温度随时间变化的曲线，如下图所示，从左到右依次为厚度为 1 毫米、 2 毫米和 10 毫米的肉片，汤的温度是 80℃，肉片开始时是 0℃。
7 s" J9 h1 M% u/ F: Q+ }3 `1 i1 a4 h4 {
2 j$ R' }) w% k& O" T' Y
    我们可以看到，比较薄的肉片放到热汤里，几秒钟的时间就会达到和汤几乎一致的温度，而比较厚的肉片中心需要很长的时间才能够热起来，对于米线来说，汤可能早就凉了。这个时间可以从前面给出的解里面估计出来，如果肉片的厚度L是以毫米为单位的话，那么需要大概数个L2秒才能把肉片的中心加热到汤的温度。涮过羊肉的读者应该在生活里面做过实验的验证了。
" D7 M! i$ ^( n* o0 O. d" E4 [
/ F4 ]' z# r2 P0 O' g  q: _    从以上的分析来说，比较薄的肉片（毫米厚度）在数秒之内就可以加热到汤的温度，而这个时间内汤并不会变凉，因此，只要汤的温度比杀死旋毛虫幼虫所需的62.2摄氏度高一些（参阅果壳网《旋毛虫风波，过桥米线还能吃吗？》），则完全不需要担心。这在一般的过桥米线里面应该不是问题。
* q6 C# a+ x  L6 y7 g
% s% `$ J) s  m9 D& d    但是如果肉片比较厚，或者肉片互相叠放（相当于增加了肉片的厚度），结果就是需要很长的时间才能将中心加热到能够杀死旋毛虫幼虫的温度。万一肉本身带有旋毛虫幼虫的话，食用起来恐怕就会有感染的风险。这里也建议大家把肉片放到汤里面的时候把肉片分开一些，不要叠放在一起。
7 u7 w" V! q& A; u6 E3 o# t! S: ?/ Z) E. B! r$ `8 T
模型的局限与应用
1 D% ^. c: N. M, k7 d
( K0 c# F, n. ^5 Q! s( k# a    需要注意的是，在这个分析里面，我们做了一些近似。8 T9 r5 |; V3 T7 q% }4 f, S

' ]+ u- \" v+ c    由于汤比肉多得多，我们假设在肉片升温的过程中汤的温度是不变的。和肉片接触的汤在把热传递给肉片时温度会降低，会在汤里面形成复杂的流体运动（对流）；而且随着时间的变化，汤是会变凉的。但是如果我们考虑的时间比汤变量所需的小的多，这个假设就没有问题。因此，对于比较薄的肉片，这个假设是合理的，对于比较厚的肉片加热需要的时间会比较长，实际的情况会比这个分析更差一些。
, ~9 ?4 F# b. F# ]9 `: s  Q  A& h4 d) c; ~
    在分析里面，我们用了水的热扩散系数。实际测量得出的猪肉的热扩散系数和水的确比较接近[4]。肉里面绝大部分都是水，含有的蛋白质和脂肪等有机物会影响肉片的传导系数，但是改变不会太大也不会产生太大的影响。哈佛大学应用数学和应用物理系的Michael Brenner教授曾经在一个科普讲座里面用金枪鱼块测量了水的热扩散系数，得出了很准确的数字[3] 。而且，即使这个估计不准确也没有关系，比如说，如果把肉片换成导热性更差的木头，其热扩散系数略大于水的一半[2] ，这样求出来的时间会是原来的约2倍，而不会改变结果的数量级。也就是说，即使是1mm厚的木片插到汤里面的话， 5 ~ 6 秒钟后木片中心就应该跟汤的温度基本一致了。9 k+ H3 ?* ]) V" c# t
- z8 {$ A3 p+ D$ f4 G  \
    模型告诉我们：肉片的厚度会导致内部温度上升时间的极大延长。但因为这只是理论上的计算，而实际生活中存在很多变量，所以最后还请以相关标准为指导。例如，根据《云南省旋毛虫病公共卫生应急预案（试行）》的建议，在吃过桥米线时最好能将肉片在 85℃ 的汤中浸烫 1 - 2 分钟。
1 ~. V: W" E2 }" d5 A7 s6 L4 @% Q
% T$ n( g: l! t参考资料：
& M; U+ Z% }9 c/ J4 G8 Z1、维基百科，“heat equation”，http://en.wikipedia.org/wiki/Heat_equation
& i. j# q1 H7 H' ?. C$ o1 `$ h: q& P# S
2、维基百科，“Thermal diffusivity”，http://en.wikipedia.org/wiki/Thermal_diffusivity% j, c5 f) k6 d5 x* {

" q! v2 i8 x5 b# @3、Boulder Summer School 2011, Public lecture “Science and Cooking”./ [) Z9 O8 y9 n$ b0 S
/ |. {4 M# P/ N1 A
4、Maria Elena Sosa-Morales et. al., Journal of Food Engineering 77, 731-738 (2006).
2 O, n2 e+ g- z
/ b. E, c8 G: e& u# H- V7 Y附推导过程：
( J1 e7 `: [2 r5 ?6 J0 E2 i
0 @7 A0 Z$ P6 T
[点击查看原图]; A0 O  o% ]0 R( [% z/ {& \/ w

5 S# G; b6 g2 f+ B- A* l1 W
" j" \, m4 E( W8 E7 S: ^9 q: C  ?) d: g1 i6 Y  v/ X
（本文转自果壳网）
1 X8 }- w' Z6 g8 w% N# Y) j5 _  v% o
  ^' r- }: m& ~& v+ n

9 Q8 y( W# C; q5 H! q1 r: ?5 {



欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5