A numerical assessment of the random walk particle tracking method for heterogeneous aquifers
|
|
- Simon Baldwin
- 5 years ago
- Views:
Transcription
1 288 Calibration and Reliability in Groundwater Modelling: From Uncertainty to Decision Making (Proceedings of ModelCARE 2005, The Hague, The Netherlands, June 2005). IAHS Publ. 304, A numerical assessment of the random walk article tracking method for heterogeneous aquifers P. SALAMON, D. FERNÀNDEZ-GARCIA & J. J. GÓMEZ-HERNÁNDEZ Universidad Politécnica de Valencia, Deartamento de Ingeniería Hidráulica y Medio Ambiente, Camino de Vera, s/n., E Valencia, Sain salamon@hma.uv.es Abstract The random walk article tracking method has become an imortant tool for the uncertainty assessment of solute transort models due to its comutational efficiency and the non-existence of numerical disersion. Yet, in highly heterogeneous aquifers the smoothness assumtion of the velocity field, requisite for a correct solution of the random walk equations, may not hold anymore. Three different numerical methods to overcome this roblem are numerically evaluated. The results demonstrate that in an aquifer with an isotroic satial correlation, and following a multi-gaussian random function with a small connectivity of extreme values, all techniques are able to model solute transort correctly, even for a hydraulic conductivity variance of σ 2 lnk = 4. Keywords heterogeneity; local mass conservation; numerical imlementation; random walk article tracking; solute transort INTRODUCTION Eulerian transort models are often lagued by numerical disersion or artificial oscillations, esecially for advection-dominated roblems. To reduce these roblems, a higher grid resolution and smaller time stes may be alied, resulting in large comutational times. One alternative to solve transort in heterogeneous orous media is the Lagrangian aroach. In articular, the random walk article tracking (RWPT) method, treats the transort of a solute by artitioning the solute mass into a large number of articles. It moves each article through the orous medium using the velocity field obtained from the solution of the flow equation to simulate advection, and adds a Brownian random dislacement to simulate disersion. This aroach avoids solving the transort equation directly and therefore is virtually free of numerical disersion and artificial oscillations. Furthermore, comutational times, even for models with a high grid resolution and characterized by strong heterogeneities, are significantly smaller than using the traditional Eulerian, mixed Eulerian-Lagrangian, or total variation diminishing (TVD) schemes (e.g. Tomson & Gelhar, 1990; Zheng & Bennett, 2002). These advantages have made the RWPT method a oular choice for comlex, highresolution transort roblems, inverse modelling, and uncertainty assessment of contaminant transort.
2 Numerical assessment of the random walk article tracking method for heterogeneous aquifers 289 BASIC PRINCIPLES Random walk article tracking is a method from statistical hysics which has been used in the analysis of disersion and diffusion rocesses in orous media. It is based on a similarity between the Fokker-Planck equation and the advection disersion equation. Using the similarity between these equations, solute concentrations are reresented by the density function of articles that are moving following the Itô- Taylor integration scheme (Gardiner, 1990): ( t + Δt) = ( t) + A(, t) Δt + B(, t) Δw Δt (1) T with A = u + D and 2 D = B B where (t) is the osition of a article at time t; A is a drift vector; B is a tensor defining the strength of disersion which is related to the disersion tensor D as shown above; u is the groundwater velocity; and Δw, an uncorrelated stochastic force, is a vector of indeendent normally distributed random variables with zero mean and unit variance. However, using Stratonovich s interretation of a stochastic integral leads to the following equation, by which the density function of articles also fulfils the advection disersion equation: Δt Δt ( t + Δt) = ( t) + A, t + Δt + B, t + Δw Δt (2) T T with: A = u + B ( B ) and 2 D = B B 2 It can be seen that the drift vector and the disersion tensor are evaluated at a moment halfway along the time ste to the next article osition. Equation (2) is rarely used for the RWPT method, as this rocedure requires an additional iterative scheme within each time ste. Nevertheless, the modified velocity contains a derivative term reduced by a factor of 1/2. THE PROBLEM OF LOCAL SOLUTE MASS CONSERVATION Numerical imlementation of the random walk equations is relatively simle, with one excetion. When solving the flow equation using numerical methods the resulting hydraulic heads and the associated velocity field are usually comuted at discrete oints. Yet, simulation of solute transort by the random walk methodology requires continuous information about the velocity field. Therefore, a ma of velocities from this discrete information has to be generated. This velocity ma should fulfil the local fluid mass balance at any location and the local solute mass conservation at any gridcell interface. In general, there is not a simle solution to this roblem, but several aroaches have been roosed in the literature. The interolation method The interolation method is certainly the aroach most commonly used to address this roblem. It uses linear interolation to obtain the groundwater velocity at any
3 290 P. Salamon et al. oint and bilinear interolation to calculate the gradient of the disersion term and the disersion tensor. By means of this so-called hybrid scheme, the interolation method fulfils both the fluid mass balance and the local solute mass balance. Nevertheless, using bilinear interolation introduces a greater smoothing, which can lead to errors esecially for highly heterogeneous aquifers. LaBolle et al. (1996) resented a detailed analysis of this method. The reflection rincile The reflection rincile was first resented by Uffink (1985) and is based on the idea of diffusion across a boundary with different diffusion coefficients. He suggested that art of a set of articles crossing from one domain into another with different diffusion coefficients must be reflected in order to obtain the correct concentration distribution across the boundary between the domains. This is done in RWPT by assigning a certain robability for a article of crossing the domain interface, which deends on the diffusion coefficients on either side. The local flow mass balance is satisfied by using linear interolation to calculate groundwater velocity. Various authors (Ackerer, 1987; Cordes & Rouvé, 1991; Semra et al., 1993; Hoteit et al., 2002) have used this aroach and suggested slight modifications. Generalized stochastic differential equations Strictly seaking, stochastic theory only alies when transort roerties are smooth functions in sace. LaBolle et al. (2000) resented generalized stochastic differential equations (GSDE) for the case of discontinuous transort roerties and their numerical integration for RWPT. In ractical terms this means that, for the case of isotroic disersion and homogeneous orosity, a two-ste rocedure is used in order to account for the local solute mass balance. First, the velocity at the articles osition, shifted only by the disersive art, is evaluated. Then, this velocity is used to calculate the final article ste. The local fluid mass balance is fulfilled by using linear interolation as in the other two aroaches. NUMERICAL ANALYSIS The aroaches were numerically imlemented into the random walk transort model RW3D (Fernàndez-Garcia et al., 2005). Solute transort in the following synthetic cases is conservative, two-dimensional and ore-scale disersion is assumed to be isotroic (α = 0.01 m). A constant dislacement scheme is used as it is comutationally more efficient than the constant time ste scheme (Wen & Gómez-Hernández, 1996). MODFLOW (McDonald & Harbaugh, 1988) was used to solve the flow equation and to comute the cell-interface velocities. Two heterogeneous hydraulic conductivity fields with a discretization of Δx = Δy = 0.3 m and a total domain size of x = 60 m and y = 21 m were generated using stochastic
4 Numerical assessment of the random walk article tracking method for heterogeneous aquifers 291 simulation techniques as shown in Fig. 1. Whereas for model A a Gaussian random function is used, model B was created using indicator simulation. Satial correlation is isotroic for model A (λ x = λ y = 1.2 m) and strongly anisotroic for model B (λ x = 20 m, λ y = 1.2 m). The satial correlation for both fields is modelled using an exonential semivariogram and the indicator simulation was erformed using a mosaic model with nine thresholds. The variance is increased from σ 2 lnk = 0.5 to 4 in order to simulate increasing heterogeneity and the mean hydraulic conductivity is ket constant at lnk = m s -1. In this article, the interolation method, using the Itô or the Stratonovich interretation of a stochastic integral, and the GSDE method in highly heterogeneous aquifers are evaluated. The third order TVD (ULTIMATE) scheme, included in the MT3DMS model (Zheng & Wang, 1999), was used for comarison with these methods. This scheme offers an attractive alternative to the traditional mixed Eulerian-Lagrangian schemes (e.g. method of characteristics, MOC) due to its mass conservation roerty and the ability of reducing numerical disersion (Zheng & Wang, 1999; Zheng & Bennett, 2002). Nevertheless, as the resent roblem is strongly advection-dominated, the domain was refined to a discretization of Δx = Δy = 0.1 m to minimize any ossible numerical disersion. One model run for model B with σ 2 lnk = 2 was reeated with a discretization of Δx = Δy = m and the two results obtained with the TVD scheme were comared. Neither artificial oscillations nor numercial disersion were observed. The differences between the different numerical schemes are investigated by comaring the time-deendence of corresonding solute lume satial moments. The satial moments are calculated following Tomson & Gelhar (1990) Model A Model B Y Y Fig. 1 Models generated using sequential Gaussian simulation (Model A) and indicator simulation (Model B) with different satial correlation. Model dimensions are given in metres; hydraulic conductivity is as lnk RESULTS All methods erformed well for model A, even when resenting an extreme heterogeneity of σ 2 lnk = 4 (Fig. 2). For model B, all the techniques illustrate good results u to a hydraulic conductivity variance of σ 2 lnk = 2. Increasing the heterogeneity for this model, only the interolation method in combination with the Itô- or Stratonovich- Fokker-Planck equation is able to reroduce the average velocity and the macrodisersivity correctly (Fig. 3). The GSDE method underestimates the average velocity and overestimates the macrodisersivity as deicted in Fig. 3. With resect to the numerical imlementation, the interolation method in combination with the Itô-Fokker-Planck equation is the easiest method. All other
5 292 P. Salamon et al Average Velocity U(t) [m/s] Longitudinal Macrodisersivity A L (t) [m] Transient Variation Diminishing (TVD) Scheme Interolation Method (Ito-Fokker-Planck) Interolation Method (Stratonovich-Fokker-Planck) Generalized Stochastic Differential Equations Fig. 2 Average velocity and longitudinal macrodisersivity using model A with σ 2 lnk = 4 (where G (t) denotes the location of the centre of mass at time t). Average Velocity U(t) [m/s] Longitudinal Macrodisersivity A L (t) [m] Transient Variation Diminishing (TVD) Scheme Interolation Method (Ito-Fokker-Planck) Interolation Method (Stratonovich-Fokker-Planck) Generalized Stochastic Differential Equations Fig. 3 Average velocity and longitudinal macrodisersivity using model B with σ 2 lnk = 2 (where G (t) denotes the location of the centre of mass at time t). aroaches require either an additional iterative scheme (Stratonovich interretation), or are comlex to imlement for three-dimensional, heterogeneous systems (reflection rincile, GSDE). Finally, it should be noted that comutation times, esecially for the strongly heterogeneous cases, are significantly lower for the RWPT method than for the TVD scheme. Whereas the TVD scheme required from aroximately 30 minutes (for model A with σ 2 lnk = 0.5) to several hours (for model B with σ 2 lnk = 4) to solve the transort equation, RWPT, using 2000 articles, required aroximately 15 minutes, clearly demonstrating the comutational efficiency.
6 Numerical assessment of the random walk article tracking method for heterogeneous aquifers 293 Acknowledgements This work was funded by the Sanish Ministry of Science & Technology, Project ESCOPA REN REFERENCES Ackerer, Ph. (1987) Random-walk method to simulate ollutant transort in alluvial aquifers or fractured rocks. In: Advances in Analytical and Numerical Groundwater Flow and Quality Modelling (ed. by E. Custodio, A. Gurgui & J. P. Lobo Ferreira), NATO ASI Series C, vol Cordes, C. H. D. & Rouvé,G. (1991) A new very efficient algorithm for article tracking in layered aquifers. In: Comuter Methods in Water Resources II, vol. 1, Groundwater Modelling and Pressure Flow (ed. by D. B. Sari et al.), Sringer Verlag, New York, USA. Fernàndez-Garcia, D., Illangasekare, T. H. & Rajaram, H. (2005) Differences in the scale deendence of disersivity estimated from temoral and satial moments in chemically and hysically heterogeneous orous media. Adv. Water Resour. 28(7), Gardiner, C. W. (1990) Handbook of Stochastic Methods for Physics, Chemistry, and the Natural Sciences. Sringer Verlag, New York, USA. Hoteit, H., Mosé, R., Younes, A., Lehmann, F. & Ackerer, Ph. (2002) Three-dimensional modeling of mass transfer in orous media using the mixed hybrid finite elements and the random walk methods. Math. Geology (34)4, LaBolle, E. M., Fogg, G. E. & Tomson, A. F. B. (1996) Random-walk simulation of transort in heterogeneous orous media: Local mass-conservation roblem and imlementation methods. Water Resour. Res. 32(3), LaBolle, E. M., Quastel, J., Fogg, G. E. & Gravner, J. (2000) Diffusion rocesses in comosite orous media and their numerical integration by random walks: generalized stochastic differential equations with discontinuous coefficients. Water Resour. Res. 36(3), McDonald, M. G. & Harbaugh, A. W. (1988) A Modular Three-dimensional Finite-difference Ground-water Flow Model. US Geol. Survey Techniques of Water Resources Investigations, Book 6. Semra, K., Ackerer, P. & Mosé, R. (1993) Three-dimensional groundwater quality modelling in heterogeneous media. In: Water Pollution II: Modeling, Measuring and Predicition, Comut. Mech., Billerica, Massachusetts, USA. Tomson, A. F. B. & Gelhar, L. W. (1990) Numerical simulation of solute transort in three-dimensional, randomly heterogeneous orous media. Water Resour. Res. 26(10), Uffink, G. J. M. (1985) A random walk method for the simulation of macrodisersion in a stratified aquifer. In: Relation of Groundwater Quality and Quantity, IAHS Publ IAHS Press, Wallingford, UK. Wen,. H. & Gómez-Hernández, J. J. (1996) The constant dislacement scheme for tracking articles in heterogeneous aquifers. Groundwater 34(1), Zheng, C. & Wang, P. P. (1999) MT3DMS, a modular three-dimensional multisecies transort model for simulation of advection, disersion, and chemical reactions of contaminants in groundwater systems. Reort to the US Army Cors of Engineers, Washington, DC, USA. Zheng, C. & Bennett, G. D. (2002) Alied Contaminant Transort Modeling, 2nd edn. John Wiley and Sons, Inc., New York, USA.
7
Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationslow velocity fast velocity
Hydrodynamic disersion Disersive transort of contaminants results from the heterogeneous distribution of water flow velocities within and between different soil ores (Figure 1, left). Disersion can be
More informationMATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK
Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment
More informationRadial Basis Function Networks: Algorithms
Radial Basis Function Networks: Algorithms Introduction to Neural Networks : Lecture 13 John A. Bullinaria, 2004 1. The RBF Maing 2. The RBF Network Architecture 3. Comutational Power of RBF Networks 4.
More informationEffect of geometry on flow structure and pressure drop in pneumatic conveying of solids along horizontal ducts
Journal of Scientific LAÍN & Industrial SOMMERFELD Research: PNEUMATIC CONVEYING OF SOLIDS ALONG HORIZONTAL DUCTS Vol. 70, February 011,. 19-134 19 Effect of geometry on flow structure and ressure dro
More informationPaper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation
Paer C Exact Volume Balance Versus Exact Mass Balance in Comositional Reservoir Simulation Submitted to Comutational Geosciences, December 2005. Exact Volume Balance Versus Exact Mass Balance in Comositional
More informationVISCOELASTIC PROPERTIES OF INHOMOGENEOUS NANOCOMPOSITES
VISCOELASTIC PROPERTIES OF INHOMOGENEOUS NANOCOMPOSITES V. V. Novikov ), K.W. Wojciechowski ) ) Odessa National Polytechnical University, Shevchenko Prosekt, 6544 Odessa, Ukraine; e-mail: novikov@te.net.ua
More informationEstimating function analysis for a class of Tweedie regression models
Title Estimating function analysis for a class of Tweedie regression models Author Wagner Hugo Bonat Deartamento de Estatística - DEST, Laboratório de Estatística e Geoinformação - LEG, Universidade Federal
More informationDeriving Indicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V.
Deriving ndicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V. Deutsch Centre for Comutational Geostatistics Deartment of Civil &
More informationOn the Fluid Dependence of Rock Compressibility: Biot-Gassmann Refined
Downloaded 0/9/3 to 99.86.4.8. Redistribution subject to SEG license or coyright; see Terms of Use at htt://library.seg.org/ On the luid Deendence of Rock Comressibility: Biot-Gassmann Refined Leon Thomsen,
More informationA Model for Randomly Correlated Deposition
A Model for Randomly Correlated Deosition B. Karadjov and A. Proykova Faculty of Physics, University of Sofia, 5 J. Bourchier Blvd. Sofia-116, Bulgaria ana@hys.uni-sofia.bg Abstract: A simle, discrete,
More information8 STOCHASTIC PROCESSES
8 STOCHASTIC PROCESSES The word stochastic is derived from the Greek στoχαστικoς, meaning to aim at a target. Stochastic rocesses involve state which changes in a random way. A Markov rocess is a articular
More informationModeling Chlorine Decay in Pipes using Two- State Random Walk Approach
International Journal of Engineering and Advanced Technology (IJEAT) ISSN: 2249 8958, Volume-4 Issue-3, February 2015 Modeling Chlorine Decay in Pies using Two- State Random Walk Aroach Y. M. Mahrous,
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationFeedback-error control
Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller
More informationAnalysis of Pressure Transient Response for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia
roceedings World Geothermal Congress 00 Bali, Indonesia, 5-9 Aril 00 Analysis of ressure Transient Resonse for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia Jorge A.
More informationarxiv: v1 [physics.data-an] 26 Oct 2012
Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch
More informationUncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning
TNN-2009-P-1186.R2 1 Uncorrelated Multilinear Princial Comonent Analysis for Unsuervised Multilinear Subsace Learning Haiing Lu, K. N. Plataniotis and A. N. Venetsanooulos The Edward S. Rogers Sr. Deartment
More informationSIMULATION OF DIFFUSION PROCESSES IN LABYRINTHIC DOMAINS BY USING CELLULAR AUTOMATA
SIMULATION OF DIFFUSION PROCESSES IN LABYRINTHIC DOMAINS BY USING CELLULAR AUTOMATA Udo Buschmann and Thorsten Rankel and Wolfgang Wiechert Deartment of Simulation University of Siegen Paul-Bonatz-Str.
More informationP043 Anisotropic 2.5D - 3C Finite-difference Modeling
P04 Anisotroic.5D - C Finite-difference Modeling A. Kostyukevych* (esseral echnologies Inc.), N. Marmalevskyi (Ukrainian State Geological Prosecting Institute), Y. Roganov (Ukrainian State Geological Prosecting
More informationEffective conductivity in a lattice model for binary disordered media with complex distributions of grain sizes
hys. stat. sol. b 36, 65-633 003 Effective conductivity in a lattice model for binary disordered media with comlex distributions of grain sizes R. PIASECKI Institute of Chemistry, University of Oole, Oleska
More informationON THE GRID REFINEMENT RATIO FOR ONE-DIMENSIONAL ADVECTIVE PROBLEMS WITH NONUNIFORM GRIDS
Proceedings of COBEM 2005 Coyright 2005 by ABCM 18th International Congress of Mechanical Engineering November 6-11, 2005, Ouro Preto, MG ON THE GRID REFINEMENT RATIO FOR ONE-DIMENSIONAL ADVECTIVE PROBLEMS
More informationAn Investigation on the Numerical Ill-conditioning of Hybrid State Estimators
An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators S. K. Mallik, Student Member, IEEE, S. Chakrabarti, Senior Member, IEEE, S. N. Singh, Senior Member, IEEE Deartment of Electrical
More informationTime Domain Calculation of Vortex Induced Vibration of Long-Span Bridges by Using a Reduced-order Modeling Technique
2017 2nd International Conference on Industrial Aerodynamics (ICIA 2017) ISBN: 978-1-60595-481-3 Time Domain Calculation of Vortex Induced Vibration of Long-San Bridges by Using a Reduced-order Modeling
More informationCHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules
CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules. Introduction: The is widely used in industry to monitor the number of fraction nonconforming units. A nonconforming unit is
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationScaling Multiple Point Statistics for Non-Stationary Geostatistical Modeling
Scaling Multile Point Statistics or Non-Stationary Geostatistical Modeling Julián M. Ortiz, Steven Lyster and Clayton V. Deutsch Centre or Comutational Geostatistics Deartment o Civil & Environmental Engineering
More informationNumerical Methods for Particle Tracing in Vector Fields
On-Line Visualization Notes Numerical Methods for Particle Tracing in Vector Fields Kenneth I. Joy Visualization and Grahics Research Laboratory Deartment of Comuter Science University of California, Davis
More information5. PRESSURE AND VELOCITY SPRING Each component of momentum satisfies its own scalar-transport equation. For one cell:
5. PRESSURE AND VELOCITY SPRING 2019 5.1 The momentum equation 5.2 Pressure-velocity couling 5.3 Pressure-correction methods Summary References Examles 5.1 The Momentum Equation Each comonent of momentum
More informationEuler-Lagrange modeling of a gas-liquid stirred reactor with consideration. of bubble breakage and coalescence
Euler-Lagrange modeling of a gas-liquid stirred reactor with consideration of bubble breakage and coalescence R. Sungkorn Institute for Process and Particle Engineering, Graz University of Technology,
More informationLow field mobility in Si and GaAs
EE30 - Solid State Electronics Low field mobility in Si and GaAs In doed samles, at low T, ionized imurity scattering dominates: τ( E) ------ -------------- m N D πe 4 ln( + γ ) ------------- + γ γ E 3
More information1 University of Edinburgh, 2 British Geological Survey, 3 China University of Petroleum
Estimation of fluid mobility from frequency deendent azimuthal AVO a synthetic model study Yingrui Ren 1*, Xiaoyang Wu 2, Mark Chaman 1 and Xiangyang Li 2,3 1 University of Edinburgh, 2 British Geological
More informationProbability Estimates for Multi-class Classification by Pairwise Coupling
Probability Estimates for Multi-class Classification by Pairwise Couling Ting-Fan Wu Chih-Jen Lin Deartment of Comuter Science National Taiwan University Taiei 06, Taiwan Ruby C. Weng Deartment of Statistics
More informationEvaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models
Evaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models Ketan N. Patel, Igor L. Markov and John P. Hayes University of Michigan, Ann Arbor 48109-2122 {knatel,imarkov,jhayes}@eecs.umich.edu
More informationTowards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK
Towards understanding the Lorenz curve using the Uniform distribution Chris J. Stehens Newcastle City Council, Newcastle uon Tyne, UK (For the Gini-Lorenz Conference, University of Siena, Italy, May 2005)
More informationProbability of dike failure due to uplifting and piping
Probability of dike failure due to lifting and iing J.M. van Noortwijk KV Consultants, P.O. Bo 0, 803 AC Lelystad, The Netherlands A.C.W.M. Vrouwenvelder TNO Building and Construction Research, P.O. Bo
More informationRecursive Estimation of the Preisach Density function for a Smart Actuator
Recursive Estimation of the Preisach Density function for a Smart Actuator Ram V. Iyer Deartment of Mathematics and Statistics, Texas Tech University, Lubbock, TX 7949-142. ABSTRACT The Preisach oerator
More informationUniformly best wavenumber approximations by spatial central difference operators: An initial investigation
Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations
More informationLogistics Optimization Using Hybrid Metaheuristic Approach under Very Realistic Conditions
17 th Euroean Symosium on Comuter Aided Process Engineering ESCAPE17 V. Plesu and P.S. Agachi (Editors) 2007 Elsevier B.V. All rights reserved. 1 Logistics Otimization Using Hybrid Metaheuristic Aroach
More informationLower Confidence Bound for Process-Yield Index S pk with Autocorrelated Process Data
Quality Technology & Quantitative Management Vol. 1, No.,. 51-65, 15 QTQM IAQM 15 Lower onfidence Bound for Process-Yield Index with Autocorrelated Process Data Fu-Kwun Wang * and Yeneneh Tamirat Deartment
More informationImplementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis
CST0 191 October, 011, Krabi Imlementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis Chakrit Suvanjumrat and Ekachai Chaichanasiri* Deartment of Mechanical Engineering, Faculty
More informationOn Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm
On Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm Gabriel Noriega, José Restreo, Víctor Guzmán, Maribel Giménez and José Aller Universidad Simón Bolívar Valle de Sartenejas,
More informationA Bound on the Error of Cross Validation Using the Approximation and Estimation Rates, with Consequences for the Training-Test Split
A Bound on the Error of Cross Validation Using the Aroximation and Estimation Rates, with Consequences for the Training-Test Slit Michael Kearns AT&T Bell Laboratories Murray Hill, NJ 7974 mkearns@research.att.com
More informationPreconditioning techniques for Newton s method for the incompressible Navier Stokes equations
Preconditioning techniques for Newton s method for the incomressible Navier Stokes equations H. C. ELMAN 1, D. LOGHIN 2 and A. J. WATHEN 3 1 Deartment of Comuter Science, University of Maryland, College
More informationarxiv: v2 [cond-mat.stat-mech] 10 Feb 2016
Surface roerties and scaling behavior of a generalized ballistic deosition model in (1 + 1)-dimension arxiv:141.66v [cond-mat.stat-mech] 1 Feb 16 Baisakhi Mal, 1,, Subhankar Ray, 1, and J. Shamanna, 1
More informationUsing the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process
Using the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process P. Mantalos a1, K. Mattheou b, A. Karagrigoriou b a.deartment of Statistics University of Lund
More informationLower bound solutions for bearing capacity of jointed rock
Comuters and Geotechnics 31 (2004) 23 36 www.elsevier.com/locate/comgeo Lower bound solutions for bearing caacity of jointed rock D.J. Sutcliffe a, H.S. Yu b, *, S.W. Sloan c a Deartment of Civil, Surveying
More informationUnderstanding DPMFoam/MPPICFoam
Understanding DPMFoam/MPPICFoam Jeroen Hofman March 18, 2015 In this document I intend to clarify the flow solver and at a later stage, the article-fluid and article-article interaction forces as imlemented
More informationMaximum Entropy and the Stress Distribution in Soft Disk Packings Above Jamming
Maximum Entroy and the Stress Distribution in Soft Disk Packings Above Jamming Yegang Wu and S. Teitel Deartment of Physics and Astronomy, University of ochester, ochester, New York 467, USA (Dated: August
More informationOn split sample and randomized confidence intervals for binomial proportions
On slit samle and randomized confidence intervals for binomial roortions Måns Thulin Deartment of Mathematics, Usala University arxiv:1402.6536v1 [stat.me] 26 Feb 2014 Abstract Slit samle methods have
More informationHomogeneous and Inhomogeneous Model for Flow and Heat Transfer in Porous Materials as High Temperature Solar Air Receivers
Excert from the roceedings of the COMSOL Conference 1 aris Homogeneous and Inhomogeneous Model for Flow and Heat ransfer in orous Materials as High emerature Solar Air Receivers Olena Smirnova 1 *, homas
More informationChurilova Maria Saint-Petersburg State Polytechnical University Department of Applied Mathematics
Churilova Maria Saint-Petersburg State Polytechnical University Deartment of Alied Mathematics Technology of EHIS (staming) alied to roduction of automotive arts The roblem described in this reort originated
More informationOn Wald-Type Optimal Stopping for Brownian Motion
J Al Probab Vol 34, No 1, 1997, (66-73) Prerint Ser No 1, 1994, Math Inst Aarhus On Wald-Tye Otimal Stoing for Brownian Motion S RAVRSN and PSKIR The solution is resented to all otimal stoing roblems of
More informationEfficient & Robust LK for Mobile Vision
Efficient & Robust LK for Mobile Vision Instructor - Simon Lucey 16-623 - Designing Comuter Vision As Direct Method (ours) Indirect Method (ORB+RANSAC) H. Alismail, B. Browning, S. Lucey Bit-Planes: Dense
More informationDeveloping A Deterioration Probabilistic Model for Rail Wear
International Journal of Traffic and Transortation Engineering 2012, 1(2): 13-18 DOI: 10.5923/j.ijtte.20120102.02 Develoing A Deterioration Probabilistic Model for Rail Wear Jabbar-Ali Zakeri *, Shahrbanoo
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell October 25, 2009 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationEstimating Laterally Loaded Pile Response
Estimating Laterally Loaded Pile Resonse C. Y. Lee, PhD Deartment of Civil Engineering College of Engineering Universiti Tenaga Nasional ABSTRACT A simlified and ractical aroach for estimating the laterally
More informationCFL Conditions for Runge-Kutta Discontinuous Galerkin Methods on Triangular Grids
CFL Conditions for Runge-Kutta Discontinuous Galerkin Methods on Triangular Grids T. Toulorge a,, W. Desmet a a K.U. Leuven, Det. of Mechanical Engineering, Celestijnenlaan 3, B-31 Heverlee, Belgium Abstract
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell February 10, 2010 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationONE. The Earth-atmosphere system CHAPTER
CHAPTER ONE The Earth-atmoshere system 1.1 INTRODUCTION The Earth s atmoshere is the gaseous enveloe surrounding the lanet. Like other lanetary atmosheres, it figures centrally in transfers of energy between
More informationRisk Analysis in Water Quality Problems. Souza, Raimundo 1 Chagas, Patrícia 2 1,2 Departamento de Engenharia Hidráulica e Ambiental
Risk Analysis in Water Quality Problems. Downloaded from aselibrary.org by Uf - Universidade Federal Do Ceara on 1/29/14. Coyright ASCE. For ersonal use only; all rights reserved. Souza, Raimundo 1 Chagas,
More informationPrinciples of Computed Tomography (CT)
Page 298 Princiles of Comuted Tomograhy (CT) The theoretical foundation of CT dates back to Johann Radon, a mathematician from Vienna who derived a method in 1907 for rojecting a 2-D object along arallel
More informationEfficient Approximations for Call Admission Control Performance Evaluations in Multi-Service Networks
Efficient Aroximations for Call Admission Control Performance Evaluations in Multi-Service Networks Emre A. Yavuz, and Victor C. M. Leung Deartment of Electrical and Comuter Engineering The University
More informationEigenanalysis of Finite Element 3D Flow Models by Parallel Jacobi Davidson
Eigenanalysis of Finite Element 3D Flow Models by Parallel Jacobi Davidson Luca Bergamaschi 1, Angeles Martinez 1, Giorgio Pini 1, and Flavio Sartoretto 2 1 Diartimento di Metodi e Modelli Matematici er
More informationChapter 2 Introductory Concepts of Wave Propagation Analysis in Structures
Chater 2 Introductory Concets of Wave Proagation Analysis in Structures Wave roagation is a transient dynamic henomenon resulting from short duration loading. Such transient loadings have high frequency
More informationDistributed Rule-Based Inference in the Presence of Redundant Information
istribution Statement : roved for ublic release; distribution is unlimited. istributed Rule-ased Inference in the Presence of Redundant Information June 8, 004 William J. Farrell III Lockheed Martin dvanced
More informationarxiv: v2 [cond-mat.stat-mech] 17 Nov 2011
Quasi-stationary simulations of the contact rocess on quenched networks Silvio C. Ferreira, 1, Ronan S. Ferreira, 1 Claudio Castellano, 2, 3 and Romualdo Pastor-Satorras 4 1 Deartamento de Física, Universidade
More informationFE FORMULATIONS FOR PLASTICITY
G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND
More informationJournal of Computational Physics
Journal of Comutational Physics 30 (011) 6916 6957 Contents lists available at ScienceDirect Journal of Comutational Physics ournal homeage: www.elsevier.com/locate/c Numerical imlementation of mixing
More informationA SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE
THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE K.W. Gan*, M.R. Wisnom, S.R. Hallett, G. Allegri Advanced Comosites
More informationSTA 250: Statistics. Notes 7. Bayesian Approach to Statistics. Book chapters: 7.2
STA 25: Statistics Notes 7. Bayesian Aroach to Statistics Book chaters: 7.2 1 From calibrating a rocedure to quantifying uncertainty We saw that the central idea of classical testing is to rovide a rigorous
More informationTHREE-DIMENSIONAL POROUS MEDIA MODEL OF A HORIZONTAL STEAM GENERATOR. T.J.H. Pättikangas, J. Niemi and V. Hovi
THREE-DIMENSIONAL POROUS MEDIA MODEL OF A HORIZONTAL STEAM GENERATOR T.J.H. Pättikangas, J. Niemi and V. Hovi VTT Technical Research Centre of Finland, P.O.B. 000, FI 0044 VTT, Finland T. Toila and T.
More informationEstimation of the large covariance matrix with two-step monotone missing data
Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo
More informationDomain Dynamics in a Ferroelastic Epilayer on a Paraelastic Substrate
Y. F. Gao Z. Suo Mechanical and Aerosace Engineering Deartment and Princeton Materials Institute, Princeton University, Princeton, NJ 08544 Domain Dynamics in a Ferroelastic Eilayer on a Paraelastic Substrate
More informationA Two-State Random Walk Approach for Modeling Chlorine Decay in Water Distribution Network
International Journal of Civil & Environmental Engineering IJCEE-IJENS Vol: 15 No: 05 39 A Two-State Random Walk Aroach for Modeling Chlorine Decay in Water Distribution Network Y. M. Mahrous, Abdullah
More informationCharacterizing the Behavior of a Probabilistic CMOS Switch Through Analytical Models and Its Verification Through Simulations
Characterizing the Behavior of a Probabilistic CMOS Switch Through Analytical Models and Its Verification Through Simulations PINAR KORKMAZ, BILGE E. S. AKGUL and KRISHNA V. PALEM Georgia Institute of
More informationApplied Mathematics and Computation
Alied Mathematics and Comutation 217 (2010) 1887 1895 Contents lists available at ScienceDirect Alied Mathematics and Comutation journal homeage: www.elsevier.com/locate/amc Derivative free two-oint methods
More informationAn Analysis of Reliable Classifiers through ROC Isometrics
An Analysis of Reliable Classifiers through ROC Isometrics Stijn Vanderlooy s.vanderlooy@cs.unimaas.nl Ida G. Srinkhuizen-Kuyer kuyer@cs.unimaas.nl Evgueni N. Smirnov smirnov@cs.unimaas.nl MICC-IKAT, Universiteit
More informationChapter 1 Fundamentals
Chater Fundamentals. Overview of Thermodynamics Industrial Revolution brought in large scale automation of many tedious tasks which were earlier being erformed through manual or animal labour. Inventors
More informationStochastic integration II: the Itô integral
13 Stochastic integration II: the Itô integral We have seen in Lecture 6 how to integrate functions Φ : (, ) L (H, E) with resect to an H-cylindrical Brownian motion W H. In this lecture we address the
More informationKeywords: pile, liquefaction, lateral spreading, analysis ABSTRACT
Key arameters in seudo-static analysis of iles in liquefying sand Misko Cubrinovski Deartment of Civil Engineering, University of Canterbury, Christchurch 814, New Zealand Keywords: ile, liquefaction,
More informationAnalysis of M/M/n/K Queue with Multiple Priorities
Analysis of M/M/n/K Queue with Multile Priorities Coyright, Sanjay K. Bose For a P-riority system, class P of highest riority Indeendent, Poisson arrival rocesses for each class with i as average arrival
More informationModeling Volumetric Coupling of the Dispersed Phase using the Eulerian-Lagrangian Approach
ILASS Americas 21st Annual Conference on Liquid Atomization and Sray Systems, Orlando, FL, May 28 Modeling Volumetric Couling of the Disersed Phase using the Eulerian-Lagrangian Aroach Ehsan Shams and
More informationSYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY
SYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY FEDERICA PASQUOTTO 1. Descrition of the roosed research 1.1. Introduction. Symlectic structures made their first aearance in
More informationFig. 4. Example of Predicted Workers. Fig. 3. A Framework for Tackling the MQA Problem.
217 IEEE 33rd International Conference on Data Engineering Prediction-Based Task Assignment in Satial Crowdsourcing Peng Cheng #, Xiang Lian, Lei Chen #, Cyrus Shahabi # Hong Kong University of Science
More informationSupplementary Materials for Robust Estimation of the False Discovery Rate
Sulementary Materials for Robust Estimation of the False Discovery Rate Stan Pounds and Cheng Cheng This sulemental contains roofs regarding theoretical roerties of the roosed method (Section S1), rovides
More informationMeshless Methods for Scientific Computing Final Project
Meshless Methods for Scientific Comuting Final Project D0051008 洪啟耀 Introduction Floating island becomes an imortant study in recent years, because the lands we can use are limit, so eole start thinking
More informationUncertainty Modeling with Interval Type-2 Fuzzy Logic Systems in Mobile Robotics
Uncertainty Modeling with Interval Tye-2 Fuzzy Logic Systems in Mobile Robotics Ondrej Linda, Student Member, IEEE, Milos Manic, Senior Member, IEEE bstract Interval Tye-2 Fuzzy Logic Systems (IT2 FLSs)
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
More informationConvex Optimization methods for Computing Channel Capacity
Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem
More informationSystem Reliability Estimation and Confidence Regions from Subsystem and Full System Tests
009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 0-, 009 FrB4. System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests James C. Sall Abstract
More informationDETC2003/DAC AN EFFICIENT ALGORITHM FOR CONSTRUCTING OPTIMAL DESIGN OF COMPUTER EXPERIMENTS
Proceedings of DETC 03 ASME 003 Design Engineering Technical Conferences and Comuters and Information in Engineering Conference Chicago, Illinois USA, Setember -6, 003 DETC003/DAC-48760 AN EFFICIENT ALGORITHM
More informationBrownian Motion and Random Prime Factorization
Brownian Motion and Random Prime Factorization Kendrick Tang June 4, 202 Contents Introduction 2 2 Brownian Motion 2 2. Develoing Brownian Motion.................... 2 2.. Measure Saces and Borel Sigma-Algebras.........
More informationPreliminary Uncertainty Estimation of the Pressure Distortion Coefficient of a Pressure. Balance by FEM Calculations
Preliminary Uncertainty Estimation of the Pressure Distortion Coefficient of a Pressure Balance by FEM Calculations G. Molinar*, M. Bergoglio*, G. Mosso*,G. Buonanno**, M. Dell Isola** * Istituto di Metrologia
More informationUniform Law on the Unit Sphere of a Banach Space
Uniform Law on the Unit Shere of a Banach Sace by Bernard Beauzamy Société de Calcul Mathématique SA Faubourg Saint Honoré 75008 Paris France Setember 008 Abstract We investigate the construction of a
More informationAn Improved Calibration Method for a Chopped Pyrgeometer
96 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 17 An Imroved Calibration Method for a Choed Pyrgeometer FRIEDRICH FERGG OtoLab, Ingenieurbüro, Munich, Germany PETER WENDLING Deutsches Forschungszentrum
More informationA Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition
A Qualitative Event-based Aroach to Multile Fault Diagnosis in Continuous Systems using Structural Model Decomosition Matthew J. Daigle a,,, Anibal Bregon b,, Xenofon Koutsoukos c, Gautam Biswas c, Belarmino
More informationYixi Shi. Jose Blanchet. IEOR Department Columbia University New York, NY 10027, USA. IEOR Department Columbia University New York, NY 10027, USA
Proceedings of the 2011 Winter Simulation Conference S. Jain, R. R. Creasey, J. Himmelsach, K. P. White, and M. Fu, eds. EFFICIENT RARE EVENT SIMULATION FOR HEAVY-TAILED SYSTEMS VIA CROSS ENTROPY Jose
More information