|
syms x
D( E9 X2 n( i+ T, gK = 5;%常量
& @5 ^' K, L' Z! _9 s! I- ^h = 2;%常量 多孔储液介质厚度
4 D3 |' Z2 s3 U% D/ k) Z/ R# hD = 0.012;%半径
7 ?0 O5 d3 @1 X) Z! ~( IH = 0.2;%压缩后的高度
7 x: u6 C, `& {! U7 {" G/ ML = @( x ) 0.2.*x;
6 i! O; \$ f- P9 X+ q+ \% ?0 o% % % % 计算与孔隙率相关的brinkman数8 @" r# p: |. a9 |
c = @( t )( ( t.^3.*D.^2 )./( 36.*( 1-t ).^2.*K ) ).^2;%渗透率求解, \% W7 C. Q. f! U) u: g
t1 = 0.0033;%孔隙率8 v; ~/ m b# B* G: m
a =@( x ) L( x )./c( t1 );%a为brinkman数
7 Y1 Q; o, b- W6 `. ]0 R; ?b = H./c( t1 );, |8 D+ q5 y- H$ z# B
% % % % 计算方程中的f,控制此处f不与x一起变化( K g; \0 T3 m& a* d, U# R) R
f = ( cos( h.*b.*H )-1 )/( b.*sin( h.*b.*H ) );
3 L; l" [5 }% v% % % % 计算方程中常数余项C,此处x只与a有关9 h ^. _. C3 C* }7 k$ @) _# F+ ` J. K
o = @( x )( a( x ).^2.*( H-f )./( 2.*f-H ) );
% F: l. Y: x, [" {6 @p = @( x )( a( x ).^2./( 2.*f-H ) );
7 G( O: N2 a9 k% ?O = int( o( x ), 0, 1 );
2 ^, a$ q5 V) t! x& V( SP = int( p( x ), 0, 1 ); |
|