Particletracing method is usually used to simulate the solute transport in fracture media.In this method,the concentration at one point is proportional to the number of particles visiting this point.However,this method is rather inefficient at the points with small concentration.Few particles visit these points,which leads to violent oscillation or gives zero value of concentration.In this paper,we proposed a particletracing method with variable weights.The concentration at one point is proportional to the sum of the weights of the particles visiting it.By adjusting weights,the number of visiting particles distributes evenly in the whole range.Through this variable weights scheme,we can eliminate the violent oscillation and increase the accuracy of orders of magnitudes.
Jiang Jianguo*,Wu Jichun.
Modeling solute transport in fracture network by particle-tracing method with variable weights [J]. Journal of Nanjing University(Natural Sciences), 2016, 52(3): 464-469
{{custom_sec.title}}
{{custom_sec.title}}
{{custom_sec.content}}
参考文献
[1] Berkowitz B,Scher H.Anomalous transport in random fracture networks.Physical Review Letters,1997,79(20):4038. [2] Edery Y,Guadagnini A,Scher H,et al.Origins of anomalous transport in heterogeneous media:Structural and dynamic controls.Water Resoures Research,2014,50(2):1490-1505. [3] Tsang Y W,Tsang C F.A particletracking method for advective transport in fractures with diffusion into finite matrix blocks.Water Resources Research,2001,37(3),831-835. [4] 薛禹群,谢春红.地下水数值模拟.北京:科学出版社,2007,189-193.(Xue Y Q,Xie C H.Numerical simulation for groundwater.Beijing:Science Press,2007,189-193.) [5] Bijeljic B,Raeini A,Mostaghimi P,et al.Predictions of nonfickian solute transport in different classes of porous media using direct simulation on porescale images.Physical Review E,2013,87(1):013011. [6] Zhang X,Lv M.Persistence of anomalous dispersion in uniform porous medium demonstrated by porescale simulations.Water Resources Research,2007,43:W07437. [7] Kang P K,LeBorgne T,Dentz M,et al.Impact of velocity correlation and distribution on transport in fractured media:Field evidence and theoretical model.Water Resources Research,2015,51(2):940-959. [8] Makedonska N,Painter S L,Bui Q M,et al.Particle tracking approach for transport in threedimensional discrete fracture networks.Computational Geosciences,2015,19(5):1123-1137. [9] Dentz M,Kang P K,Le Borgne T.Continuous time random walks for nonlocal radial solute transport.Advances in Water Resources,2015,82:16-26. [10] LaBolle E M,Quastel J,Fogg G E.Random walk simulation of transport in heterogeneous porous media:Local mass conservation problem and implementation methods.Water Resources Research,1996,32:583-593. [11] Srinivasan S,Tartakovsky D M,Dentz M,et al.Random walk particle tracking simulations of nonFickian transport in heterogeneous media.Journal of Computational Physics,2010,229:4304-4314. [12] Zhang Y,Benson D,Meerschaert M,et al.On using random walks to solve the spacefractional advectiondispersion equations.Journal of Statistical Physics,2006,123:89-110. [13] Grassberger P.Prunedenriched rosenbluth method:Simulations of θ polymers of chain length up to 1000000.Physical Review E,1997,56(3):3682. [14] Prellberg T,Krawczyk J.Flat histogram version of the pruned and enriched rosenbluth method.Physical Review Letters,2004,92(12):120602. [15] Jiang J,Wu J.Simulation of continuoustime random walks by the pruned enriched method.Physical Review E,2011,84(3):036710. [16] Bijeljic B,Muggeridge A H,Blunt M J.Porescale modeling of longitudinal dispersion.Water Resources Research,2004,40(11):W11501.