非对称分数阶对流弥散的数值模拟和参数反演.pdfVIP

2013年 12月 高 等 学 校 计 算 数 学 学 报 第 35卷第 4期 非对称分数阶对流弥散的 数值模拟及参数反演 李新洁 李功胜 贾现正 f山东理工大学理学院应用数学所，淄博 255049) NUM ERICAL SIM ULATION AND PARAM ETERS INVERSION FOR NON．SYM M ETRIC FRACTIONAL ADVECTION DISPERSION EQUATION LiXinjie LiGongsheng JiaXianzheng (InstituteofAppliedMathematics，ShandongUniversityofTechnologY, Zibo255049) Abstract Thispaperdealswith numericalsolution forone—dimensionalnon- svmmetricfractionaladvectioll-dispersionequation(FADEinshort)inafint。 domain andparametersinversion．A differenceschemeispresentedbasedonthe Grtln、)lrald—Letnikovdefinitionofthefractionalderivative，andsensitivityofthe solution tothefractionalorderandtheskew parameterisanalyzedfordifferent inputsources．Furthermore，aninverseproblem ofdeterminingthefractionalo卜 deranddispersioncoefficientsimultaneouslyisinvestigatedbyfinalobservations· Numericalinversionsarecarriedoutsuccessfullybyapplyinganoptimalperturba— tionregularizationalgorithm，andimpactofthefractionalorderonthe nverSon algorithm isdiscussed． 国家 自然科学基金资助项目10926194)和山东省自然科学基金资助项目 (ZR2011AQ014)． 收稿 日期：2011—02—07． 310 · 李新洁等：非对称分数阶对流弥散的数值模拟及参数反演 第 4期 Keywords fractionaladvection—dispersionequation，Griinwald—Letnikovderiva- tive，skewnessparameter，finitedifferencescheme，parameterinversion，numerical simulation AMS(2000)subjectclassifications 35Rll，35R30，65M06 中图法分类号 O175，O241，0242 1 引 言 溶质在多孔介质中的运移一般以对流与机械弥散为主，基于Fick扩散定律的对流弥 散方程 (ADE)是描述溶质运移行为模式的经典模型 [18，241．然而，上世纪90年代 以来， 有关研究 [1,2，6，7，14，19，26]表明，用ADE模型来描述溶质的运移过程时，有时并不能很好 的模拟穿透 曲线的提前穿透或拖尾现象，同时也不能很好地解释溶质长距离传输中表现 出的弥散系数的尺度效应等问题．事实上，ADE模型揭示了溶质运移的瞬时扩散性与对称 性，但多孔介质中的溶质运移有时表现出明显的非对称性等非