

J. Jpn. Soc. Soil Phys.
 5¼No. +*-, p. /1 0, ,**0
 MathCad! "   *
Using MathCad to Simulate Solute Transport in Soil fora Graduate Soil Physics ClassKosuke NOBORIO*KLMNOL:>H% +#$no, (I      rMl<V34 +   4   !"# $
% &'

¡V:;: MathCad CE:(>H%( II ῑ)*+&* ,PDE- .

&% /012345346789 PDE :5^)

S"V;:<=% >> 1234'5% ¢   +  6 ῑ)34<?@A'BC<(89@A :* *.& ,van Genuchten and Parker,;:< !"DCE:F+32.-%<$% G&:;:'CE:H!

R

 FORTRAN, BASIC, C IJKLMNOLPQRS"JKLMTUV<W$SH%  XJKLMNOLPQYW2Z<$%efgh

igjJKklm noIp!qrstuvw^)BC ,L,ῌT- v /¢Z§ ¨ ,LῌT- t ©
Z ,T- x ª8 ,L- $%
ῌ :;: :

*῎+<~#$'#$xVH! :;:'CE:RSH:;{|&F@  X<zC$% ]%&>F Mathematica(

4¡9«>"

X12ῌ

534:%x#"yFHab<z{|&"F% W ^}' MathCad FSH'C^qrstuv

ῌ R ¥¦BC C  70 ,MῌL-- D  5[< \]^_^`F"JKLMNO
LPQ<ab.&<cdI(C((C῏(C£ ῍D¤v(t(x ῎ (x ῐ(x1234 : C£Ci for t£*, *῕x῕¬῍534 : (Cῌ(x £* for t­*, x£¬῎C£Co for t­*, x£*῏ Ci 1270 Co y®¯70°V%   534

( rMl34<F)^Y= GJKLMNOLPQY&% >> :±:;:F" ¯70

>"FF^}S"*>'VF&%]²VrMl70 ῏ *7034:;>+#$ 'C^qrstuv ,MathCad- F<>"  ῏ *534=RSH no, MathCad *>

H ,van Genuchten and Parker, +32.-% * ῌ :;

HabVa -.V% ab: ^^^`/ +*  +/ %0{SH% ab   H S":1>" **23 MathCad RSHJ* >´^^? µ,+.ῑ2/1+ ¶@·¸·¹pzAºB»C +ῑ+ῑ+Ce,t, x-£+ῑ῎ῑΐῐ῏ΐ³ + exp῍vx῏erfc,ῒῌ῔῍῎῏ ῔ , ῎ D ῐῑ῎ῑῒῐ῏ΐῒῌ῔῍῎῏ ῔erfc

ῐ

DE%TUV W +*- X ,**058 van Genuchten and Parker, +32. CCiCe Ce  erfc
CoCi῍  Dt f+f, Dt f++f,Ci, t+DtCi, t++v
R Dx
R,Dῌῌ+Ci+Ci
Dx

DEF%5>GH

Fig. +CiCi++
Dx

, 
%12 iCi, t+Dt  Ci, tῐ-./0 t+Dt  t3456 MῌL-  Dt  Dx-.78/9:$ T 129; L  

ῌ !"#$%& '()*

f+Df,D Ci+, Ci, Ci++῎-.12 i+, i, i++

56 M ῌ L-  -<= f +  f ,34-.12 i+, i 9%5>?@ABC MῌL,ῌT 12
῏i, i++ 9%5>?@ABC MῌL,ῌT IJ!KL MathCad MNO@P%Q78RS%O@?

An example of a MathCad program to simulate solute transport in soil using
an analytical solution and its result.

῎¥¦ : §¨© ª¨«V¬a& m#$­®¯¨°iF±²³ MathCad ´¯ 59

ῒ ΐ  ῍ ῎ ῏  *g3** h #$W12hZi=W MathCad #$12)῔+ ^S9BCDECFhIII

MathCad ῌ ῍ /MathCad  /#$jCkl4mnW

MathCad Mathsoft, ,**.  Zi=%"'(&LopS )῔, WGW /#$RqO&4P!"#$%&"'(&Q /#$jCklrsUV%&"'()῔+ * : +,  -./0#$123.4& jCklm *t,  -.4/u#$56%&"78%& 9: *+,  3.4

12/3.456%&?@A& 0;<=>%&?@A& #$BCD*for,  vOw!4xy%?@A& "ECF3GH- GHIJGKLMNHO& 4jCkl 4 C z{|} Ci 56WL wLPQ RSTUVWXYZC'[\%&~f ]f|} Co W ῐ =]O& )῔+ >]O&^S ῑ K #$~; 4 *keep,WL Q Ca UVW Q Ca 'KL&/ +   ]f _WL4 t, x /`abWcK )῔+|}~4#$/4 Ci, j W "d *tt : +*/ ee3**, S?@RKL / h f PC ;<%&NcK

ῌ,Fig. ,m%&¡¢RK MathCad £hZ¤dAn example of a MathCad program to simulate solute transport in soil using
the ﬁnite di#erence method.

ghijk l +*- m ,**0

60

ῑῌ-A D@59/Fig. -(A) Breakthrough curves (B.T.C.) calculated with analytical and numerical solutions, and(B) calculation errors with the numerical solution relative to the analytical solution.

IJKL@ B cd  Dt

Cῌ+ D59/ EFG R!+.* H01234)*5 67)8+,5 )859:;</=> ?@%ABefIVῌ ῎ ῐ ῏ ῍

RDt,ῌ,D      P P!vDxῌD

"#$%&  '( )*ῌ+,-.%&/ P 

DHIJKLMN3=>  OP#QPRS TU3=>^V*W v X+,G D YZ%[\]@ _3H`/ IJKLab3

Cῌ+¸M : MN&'(MRSd&:N v*.+/ cm/h, x/* cm, D*., cm,/h, R+.* ῌ- A $ %&

' ,- .(/5 67,῍89+.,῍012(<=@AB]D   MQ56-& ^M_ (/> P?/ [M:7MHI !8,-]ῌ.  > P?*.-2 E RMSE C klmnoVmtVuvXw,(c!
$ 9{R

M.M:(&'(M! &:'ῌ;z ov
2XV ¡¢() £TR" .¤,V KLMMNMOPQ/RS MathCad () [!-+>-( +* \.

 <= ,**+ ¥¦@§ABklmnx()?

MathCad () RS TR"-&UVWX

Y(ZTR"(2}  J +. +/  TR"H UVWXYZ c!cb$ PG*.-2 012HID,

RS( J +*  -&J +-  ,*:;

012CDE , F* RMSEJ +*8+- TR"$. RS/ 34!MQ-&TR"crcῌ- B  012 P+./ 34 C-  &'(M)c*+,z ( 0! .(1b]%&() *+ c!-TR ! "#!
£>?M¹ºY»¼U MathCad I? 61]&+*>?

tVuvXw()¨q2©XYªJ3«4= TR"3*῍ ¬­®((D­C#YgUVWX¡XwD TR" ¡¢¯n°Xw(±D MNRS>T\( ^_ 28+* ^+`&]abc!_ .81

R"klmnx(²!ov2XV ¡¢³

^_ +8-  - [dE c-ῌ/  . eUVWXY@A(D J +-   f ^_ gVI  ' - e ^_  ^+`&]abc!_ hi j,J +,   ^_ !>g,
klmnoVmpqr!!TR"hs+,  RS"#tVuvXwklmnx

F#klmnxa´5 MathCad ()>5)µG(H8,-
Cad (I?

MNRS¨¶ Math-D . eTR"UVWX

y$z {e%R(D|}, ~ῌ/ῌ.

012<D , F* RMSE<Fig. .

Changes in root mean squared error (RMSE)
with changes in Peclet numbers.

KLMMNd&'(MJKRS
MathCad ()?
DTR"L·Fig. /Students’ responses on using MathCad inthe soil physics course for graduate students at Iwate University.

defg² · +*- ¸ ,**062       ! MathCad"#$%&'()(* +,-.&/0)1 23 4546$789:;7<=hM MathCad $C"MN# Mk%

῍ ῒ ῑ ῎"#>?@AB;CDEFG&' 

Mathsoft ,**. : Mathcad +, K;L; (HI J 789:;7<=K; ῌ =xz= ῌ x= ῌ 9L;MNOPQR(S 789:;7<=T"UV%&WP+X(YZOP[N;

7<=ῌ=¡¢£¤¥ ,**+ : de¦§¨ῌ©ª«¬u­®¯ \]1& ^4_(`a&bcdefg$ PDE °±; $ defg² 22 : +3ῌ,/.

4hi"DE(jE  kDElm"nZ

van Genuchten, M. Th. and Parker, J.C. (+32.) : Bound-oE (pNqrst%&-.&VIIῌ ῐ῏

MathCad "#>?uv(wMN xyz{|}~ῌx^4 R. Ewing   R. Horton

ary conditions for displacement experiments
through short laboratory soil columns. Soil Sci.
Soc. Am. J., .2 : 1*-ῌ1*2.
³´µ¶ : ,**/ ´ +, µ ,2 ¶g´µ¶ : ,**0 ´ , µ + ¶