Fractional-RWHet: An enhanced solver for solute transport with both spatiotemporal memory and conditioning on local aquifer properties
|
|
- Hortense Walsh
- 6 years ago
- Views:
Transcription
1 Fractional-RWHet: An enhanced soler for solute transport with both spatiotemporal memory and conditioning on local aquifer properties Zhang Yong 1, Daid A. Benson 2, Eric M. LaBolle 3, Donald M. Reees 4 1 Desert Research Institute, yong.zhang@dri.edu, Las Vegas, NV, USA 2 Colorado School of Mines, yzhang@mines.edu, dbenson@mines.edu, Golden, CO, USA 3 Uniersity of California, emlabolle@ucdais.edu, Dais, CA, USA 4 Desert Research Institute, mreees@dri.edu, Reno, NV, USA ABSRAC racer transport through natural geologic formations is affected by physical and/or chemical heterogeneity, which is present across multiple scales. Local random walk transport models based on the 2 nd -order adection-dispersion equation, such as RWHet, can (at least theoretically) incorporate the site characterization information as fine as the measurement scale. Emerging nonlocal transport models, such as the fractional-order adection-dispersion equation, are demonstrated to capture coneniently the influence of fine scale heterogeneity (under the measurement scale) on solute transport. A combination of local and nonlocal models seems ery attractie and has drawn serious attention recently in the hydrology community. In this study, we combine fractional dynamics and random walks to derie an enhanced soler, named Fractional-RWHet, based on the nonlinear Langein approach. Fractional- RWHet can capture the sub-grid, non-fickian dispersion, including the solute retention process due to relatiely low-permeable sediments and the fast-moing plume front due to the fine-scale preferential flow paths. It also maintains the adantage of local approach, where the obserable heterogeneity (such as hydrofacies from driller s logs) and the measurement of aquifer properties (i.e., elocity from obseration wells) can be incorporated. he Fractional-RWHet therefore enables us to efficiently simulate solute transport with both spatiotemporal memory and conditioning on local aquifer properties. In addition, Fractional-RWHet is built on the solid footing of a fully Lagrangian framework, and can interface directly with the popular USGS MODFLOW code. INRODUCION Increasing concerns of the fate and longeity of contaminant in groundwater necessitates modeling of complex small-scale transport process oer ery large space and time. Many, if not most, natural geologic formations exhibit fractal structure [Cushman, 199] with multiscale physical/chemical heterogeneities. Recent studies [such as Berkowitz et al., 26] further show that the solute transport through heterogeneous media at arious scales can be non-fickan, which cannot be captured efficiently (if at all) by the classical simulator based on the 2 nd -order adection-dispersion equation (ADE). New computational theories/models are therefore needed to extend the 2 nd order ADE based simulator to 1) simulate non-fickian dispersion at all scales, 2) propagate transport information between arious scales, and 3) incorporate measurements of aquifer properties (i.e., elocity) at any practical resolutions. his study proposes an enhanced soler by combining the local simulator deeloped by LaBolle [26] and the nonlocal transport theory proposed by Benson [1998] to capture contaminant transport through aquifer/aquitard systems with physical/chemical heterogeneities, where the influences of sub-grid heterogeneity, local aquifer properties, and spatiotemporal nonlocal dependence on solute dynamics can be characterized efficiently. A preious effort to capture the sub-grid elocity can be found in Feehley et al. [2], where a dual-domain mass transfer (i.e., single-rate mass transfer) methodology is embedded into M3D [Zheng, 199]. See also the oeriew of multiscaling techniques for the study of microfluidics and nanofluidics by Hu and Li [27]. Here the local and nonlocal theories follow the definitions gien by Cushman [199, page 3]: If the constitutie ariable depends on what is happening at a point in spacetime or a ery small neighborhood of the point, the ariable is said to be local and deried from a local theory. On the other hand, if information is needed to define the constitutie ariables from regions of space and time distinct from a neighborhood of the space-time point where ealuation of the ariables is to be made, then the theories and constitutie ariables are said to be non-local in character. he
2 classical soler based on the 2 nd -order ADE is a local approach, where the ariation of solute concentration depends on a ery small neighborhood (i.e., a REV-scale grid) in space and time. he fractional-order adection-dispersion equation (fade) model deeloped by Benson [1998] is a nonlocal transport theory, where the local ariation of solute mass flux is affected by a wide range of upstream domains (i.e., multiple grids) and a long historical release of contaminants at the local point. he influence of sub-grid heterogeneity (i.e., below the REV or measurement scale) on solute transport can also be captured by the fade model [Benson, 1998; Benson et al., 28]. Also note that ariable nonlocal approaches hae been deeloped, such as the multi-rate mass transfer theory [Haggerty et al., 2] and the continuous time random walk (CRW) framework [Berkowitz et al., 26]. Cushman and Moroni [21] showed that the arious nonlocal theories could find a home in the nonlocal dispersie constitutie theory, where the space and time memory kernels define the nonlocality. Our soler can incorporate other forms of memory functions (such as the tame Léy flights) coneniently, so it is actually a generalized soler for arious memory kernels with physical/hydrological definitions. We name this enhanced soler the fractional random walk and the resulting code Fractional-RWHet. he random walk is a fully Lagrangian scheme, with encouraging adantages including minimal numerical dispersion, sub-flow-grid scale resolution of concentrations, easy implementation and near linear scaling on parallel computers due to computational expense that depends on independent particle numbers rather than grid dimensions [Green, 22]. It is also a promising tool to obiate the grid-aeraging problem when used for the local scale reaction with sharp concentration fronts. In addition, Fractional- RWHet can interface directly with the popular USGS MODFLOW code (note that other flow solutions can also be incorporated). his soler allows us to simulate non-fickian dispersion without the need to explicitly define local-scale heterogeneity in a numerical model. In the following we introduce the methodology and examples, including simulations of solute transport through arious heterogeneous media generated by a Marko Chain based geostatistical tool [Carle, 1998]. MEHODOLOGY DEVELOPMEN We proposed the following noel multiscaling fade to capture the local ariation of transport strength, the possibly heay-tailed particle retention processes, and the possibly preferential leading fronts β C ( x t) Θ +, = ( x t) M( dθ ) D( x t) C ( x t) + f ( x) δ ( t) t t, H 1 I,,, (1) where is the elocity ector, D is a scalar dispersion coefficient, ( < 2 ) is the order of the time fractional deriatie (Caputo type), β is the capacity coefficient for < 1 (note 1+ β is the retardation factor for = 1 ) and the parameter defining the waiting time distribution for 1 < 2, the sign function Θ = 1 for < 1 and -1 otherwise, I is the identity matrix, H -1 is the inerse of the scaling matrix is a function uniquely determined by the requirement of (1) haing a mild solution (see Baeumer et al. [25] for the definition of mild solution). C is the solute concentration in the total (mobile+immobile) phase. proiding the order and direction of the fractional deriaties, and M(dθ) is the mixing measure. f ( x ) When β = and and D are constant, (1) reduces to the multiscaling fade proposed by Meerschaert et al. [21]. When deriing (1), we selected the Langein equation (see eq. (3)), where the net flux due to nonlocal dispersion is the same form as the one in the Fokker-Planck equation (FPE). Similar treatment (to the ADE) was proposed by Cortis et al. [24]. his option allows the ariable dispersion coefficient, and it can also reduce significantly the computational burden and thus make the following particle tracking much more conenient. Additional experiments are needed to ealuate the FPE type of fractional dispersion, such as a detailed analysis of the influence of the gradient of D on particle moement. We deeloped a non-markoian particle tracking scheme to approximate the complex model (1). his scheme includes two steps. First, the indiidual jump ector of each particle at time step i, which is essentially a multiscaling compound Poisson process, can be calculated by [ tm/ dt] [ tm/ dt ] H Z() t = x i = R r i θ i, (2) i= 1 i= 1
3 where Z(t ) denotes the particle location at time t, n = [t / dt ] is the number of random jumps by time t using time step size dt, R i is the random length of the i-th jump, the jump direction θ i is a random unit ector drawn from the distribution of the mixing measure M(dθ), and t m is the operational time that particles spend in motion by the clock time t. he jump length of the particle along the eigenector belonging to the k-th eigenalue 1 /α k of H can be calculated by the Langein approach R1 /α k = (x k, t ) dt + D(x k, t )1 /α k dl α k (t ), k = 1, L d,, (3) where dl α is the random noise underlying an α -order Léy motion, k represents the direction of the k-th eigenector of H, and d is the dimension of H. (1) and (3) are also alid for the transient flow. -1 (a) ime nonlocal (α=2).16 (b) Space nonlocal (a=) C/C =.1 =.5 = α = 1.3 α = 1.7 α = (c) Conergence (=.5, α=2) FF RW with nt = 5 RW with nt = 2 RW with nt = RW with nt = Figure 1. Snapshots using the random walk approach (symbols) ersus the FF solution (lines) for the time fade (a) and the space fade (b). he dashed line is the ADE result. In (a), =1, D=.1, =3, α =2, =.1~.9, β=.1, and the initial instantaneous point source is at location. In (b), =1, D=1, =4, α=1.3~1.9, β=, and the source is at location 3. (c) Conergence of the random walk approximation. nt denotes the number of time steps used in the numerical experiments. In the second step, the operational time can be simulated as the number of renewals by time Streamline t > for a gien waiting time distribution with power-law probability tails. Here the particle motion is not instantaneous, and thus we can distinguish the status (mobile or immobile) for 4 each particle at any gien time. his distinction is critical in modeling field-measured plumes, since the sampling process tends to collect mobile solutes preferentially. Also note that the 2 indiidual operational time is an exponentially distributed random ariable, which can be drawn by generating a uniform ( 1) number corresponding to the CDF of the exponential function (with a mean corresponding to the mean motion time). In addition, when 1 < 2, scaling Source location Particle plume X (Horizontal) limit of the operational time process is split into two parts, the mean wait and the deiation from Figure 2. RWHet plume (left) s. Fractionalthe mean, to presere the near and intermediate RWHet plume (right). Here α=1.8, =5, β=. he time behaior. he number of jumps is streamline is exaggerated for ision purpose. associated with the max-process [Baeumer et al., 25; Benson et al., 28], since deiations from the mean waiting time may be negatie. Physically this means that the particle pauses until the next time the accumulated deiations become positie, so that the max increases. 2nd-order ADE Y (Longitudinal) 6 Space fade
4 We erified the aboe random walk scheme extensiely. A few examples are shown in Figure 1. he random walk solutions match well the fast Fourier transform (FF) solutions (Figures 1a, b). We also checked the conergence speed. he solution of (1) is equal to the density of the stochastic process obtained as a scaling limit of the CRW. herefore, the probability density of the CRW will conerge to the solution of (1) after long time or large time steps. Numerical experiments (Figure 1c) reeal that the conergence can be achieed after affordable time steps, such as. SIMULAION RESULS We then ealuated the particle plumes goerned by the time fade, underlying the fine-resolution elocity shown in Figure 2. Results (Figure 3) show that the index controls the slope of the late-time tail of BC, while β affects the mass decay rate. We also analyzed the particle plume distribution. In Figure 4, the background image shows the conductiity distribution. he longitudinal and horizontal dispersiity is.1 and.1 m, respectiely, and the molecular dispersion coefficient is m2/day. An instantaneous line source was used along x=185~2m and y=5m. he plume exhibits both leading and trailing edges, just as expected. We then tested the influence of mixing measure on particle plumes. Results (Figure 5) show clearly the moement of particles along the two preferential directions defined by the mixing measure. o drawn a conclusion, the aboe numerical tests show the applicability of Fractional-RWHet. In the next study, we will apply this enhanced soler to simulate contaminant transport through real aquifers, including the non-fickian, 3-d, plumes obsered at the well-studied MADE test site. Longitudinal Distance (m) C/C C/C We first applied the Fractional-RWHet to simulating solute transport through a heterogeneous porous medium, where the goerning equation of solute transport is the space fade. A 2-d heterogeneity model was generated first (see RWHet result =1.3, β= the background image (adction only) shown in Figure 2) [Carle, =1.5, β= =.3, β= ], and the steady-2 =.5, β=.1 state flow field was =1.7, β= -4 calculated by USGS MODFLOW code. We -4 =.7, β=.1 then released particles at -6 an upstream point, to track the moement of particles =1.5, β=3 along the streamline =1.5, β= -2 =.5, β=.2 (adection only). A new -2 subroutine, called =.5, β=.1 =1.5, β= MOVEP5, was added into -4 Fraction-RWHet, to -4 =.5, β=.1 calculate the fractional -6 dispersion with local ime ariation of transport strength. he fractional Figure 3. Breakthrough cures simulated by Fractional-RWHet dispersion was projected ersus RWHet (with adection only, so the BC is a ertical line). along the streamlines. he resultant particle plume has a heaier leading edge comparing to the plume goerned by the 2nd-order ADE (Figure 2) X (m) 2 22 Figure 4. Particle plumes goerned by the ADE (light dots) and the time fade (dark dots), with =1.5, β=2, and =5.
5 (b) 2nd-order ADE (a) Operator stable parameters (c) Space fade 35 α 1 =1.8 M 1 =.7 3 M 2 =.3 3o 2o Y (m) α 2 = Source Source Figure 5. Particle plume captured the 2nd-order ADE (b) and the space (multiscaling) fade (c). he multiscaling fade used for (c) has pre-defined operator stable parameters shown in (a). he direction of in (b) and (c) shows the mean flow direction., particles were released at the point showing with a diamond in (b) and (c). Plumes in (b) and (c) were simulated by RWHet [LaBolle, 26], and FractionalRWHet, respectiely X (m) ACKNOWLEDGMEN his work was supported by the National Science Foundation (NSF) under EAR , EAR , and EAR his paper does not necessary reflect the iew of the NSF. REFERENCES Baeumer, B., D. A. Benson, and M. M. Meerschaert, 25. Adection and dispersion in time and space, Phys. A, 35(2-4), Benson, D. A., he fractional adection-dispersion equation: Deelopment and application. Ph.D. hesis, Uniersity of Neada, Reno, 157 pp. Benson, D. A., M. M. Meerschaert, Y. Zhang, B. Baeumer, and R. Schumer, 28. Stochastic model for multi-rate mobile/immobile contaminant transport, Water Resour. Res., in reiew. Berkowitz, B., A. Cortis, M. Dentz, and H. Scher, 26. Modeling non-fickian transport on geological formations as a continuous time random walk, Re. Geophy., 44(2), RG23. Carle, S., PROGS: ransition Probability Geostatistical Software, 2.1. Uniersity of California, Dais, 84 pp. Cortis, A., C. Gallo, H. Scher, and B. Berkowitz, 24. Numerical simulation of non-fickian transport in geological formations with multiple-scale heterogeneities, Water Resour. Res., 4, W429. Cushman, J. H., 199. Dynamics of Fluids in Hierarchical Porous Media, Academic Press, San Diego, CA, 55 pp. Cushman, J. H., and M. M. Moroni, 21. Statistical mechanics with three-dimensional particle tracking elocimetry experiments in the study of anomalous dispersion. I. heory, Phys. Fluids, 13, Feehley, C. E., C. Zheng, and F. J. Molz, 2. A dual-domain mass transfer approach for modeling solute transport in heterogeneous porous media, application to the MADE site, Water Resour. Res., 36(9), Green, C.., 22. Effects of heterogeneity on reactie transport in geologic media. Ph.D. hesis, Uniersity of California, Dais, 42 pp. Haggerty, R., S. A. McKenna, and L. C. Meigs, 2. On the late-time behaiour of tracer test breakthrough cures, Water Resour. Res., 36(12), Hu, G. P., and D. Q. Li, 27. Multiscale phenomena in microfluidics and nanofluidics, Chem. Eng. Sci., 62, LaBolle, E. M., 26. RWHet: Random Walk Particle Model for Simulating ransport n Heterogeneous Permeable Media,.3.2. User s Manual and Program Documentation, 27 pp. Meerschaert, M. M., D. A. Benson, and B. Baeumer, 21. Operator Léy motion and multiscaling anomalous diffusion, Phys. Re. E, 63, Zheng, C., 199. M3D: A modular three-dimensional transport model for simulation of adection, dispersion and chemical reactions of contaminants in groundwater systems, S.S. Papadoulos and Associates, Inc., Rockile, MD, 163 pp.
Lagrangian characterization of contaminant transport through multidimensional heterogeneous media with limited heterogeneity information
Lagrangian characterization of contaminant transport through multidimensional heterogeneous media with limited heterogeneit information Zhang Yong 1, David A. Benson 1, Mark M. Meerschaert 2, Eric M. LaBolle
More informationA tempered multiscaling stable model to simulate transport in regional scale fractured media
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 37,, doi:10.1029/2010gl043609, 2010 A tempered multiscaling stable model to simulate transport in regional scale fractured media Yong Zhang,
More informationASSESSING TEMPORAL MOMENT CHARACTERISTICS OF SOLUTE TRANSPORT PHENOMENA IN HETEROGENEOUS AQUIFERS USING RANDOM WALK PARTICLE TRACKING
International Journal of GEOMATE, May, 17, Vol.1, Issue 33, pp. 114-11 Geotec., Const. Mat. & En., ISSN:186-99, Japan, DOI: http://dx.doi.org/1.166/17.33.641 ASSESSING TEMPORAL MOMENT CHARACTERISTICS OF
More informationWATER RESOURCES RESEARCH, VOL. 43, W05439, doi: /2006wr004912, 2007
Click Here for Full Article WATER RESOURCES RESEARCH, VOL. 43,, doi:0.029/2006wr00492, 2007 Space-fractional advection-dispersion equations with variable parameters: Diverse formulas, numerical solutions,
More informationPredicting the Tails of Breakthrough Curves in Regional-Scale Alluvial Systems
Predicting the Tails of Breakthrough Curves in Regional-Scale Alluvial Systems by Yong Zhang 1,DavidA.Benson 2,andBorisBaeumer 3 Abstract The late tail of the breakthrough curve (BTC) of a conservative
More information0 a 3 a 2 a 3 0 a 1 a 2 a 1 0
Chapter Flow kinematics Vector and tensor formulae This introductory section presents a brief account of different definitions of ector and tensor analysis that will be used in the following chapters.
More informationVISUAL PHYSICS ONLINE RECTLINEAR MOTION: UNIFORM ACCELERATION
VISUAL PHYSICS ONLINE RECTLINEAR MOTION: UNIFORM ACCELERATION Predict Obsere Explain Exercise 1 Take an A4 sheet of paper and a heay object (cricket ball, basketball, brick, book, etc). Predict what will
More informationUNDERSTAND MOTION IN ONE AND TWO DIMENSIONS
SUBAREA I. COMPETENCY 1.0 UNDERSTAND MOTION IN ONE AND TWO DIMENSIONS MECHANICS Skill 1.1 Calculating displacement, aerage elocity, instantaneous elocity, and acceleration in a gien frame of reference
More informationDomenico Solution Is It Valid?
Domenico Solution Is It Valid? by V. Sriniasan 1,T.P.Clement, and K.K. Lee 3 Abstract The Domenico solution is widely used in seeral analytical models for simulating ground water contaminant transport
More informationDO PHYSICS ONLINE. WEB activity: Use the web to find out more about: Aristotle, Copernicus, Kepler, Galileo and Newton.
DO PHYSICS ONLINE DISPLACEMENT VELOCITY ACCELERATION The objects that make up space are in motion, we moe, soccer balls moe, the Earth moes, electrons moe, - - -. Motion implies change. The study of the
More informationChapter 2: 1D Kinematics Tuesday January 13th
Chapter : D Kinematics Tuesday January 3th Motion in a straight line (D Kinematics) Aerage elocity and aerage speed Instantaneous elocity and speed Acceleration Short summary Constant acceleration a special
More informationMOTION OF FALLING OBJECTS WITH RESISTANCE
DOING PHYSICS WIH MALAB MECHANICS MOION OF FALLING OBJECS WIH RESISANCE Ian Cooper School of Physics, Uniersity of Sydney ian.cooper@sydney.edu.au DOWNLOAD DIRECORY FOR MALAB SCRIPS mec_fr_mg_b.m Computation
More informationParameters Identification of Equivalent Circuit Diagrams for Li-Ion Batteries
Parameters Identification of Equialent Circuit Diagrams for Li-Ion eries Ahmad ahmoun, Helmuth Biechl Uniersity of Applied ciences Kempten Ahmad.ahmoun@stud.fh-empten.de, biechl@fh-empten.de Abstract-eries
More informationarxiv: v1 [cond-mat.stat-mech] 5 Oct 2018
On the canonical distributions of a thermal particle in the weakly confining potential of special type Tatsuaki Wada Region of Electrical and Electronic Systems Engineering, Ibaraki Uniersity, Nakanarusawa-cho,
More informationSUPPLEMENTARY MATERIAL. Authors: Alan A. Stocker (1) and Eero P. Simoncelli (2)
SUPPLEMENTARY MATERIAL Authors: Alan A. Stocker () and Eero P. Simoncelli () Affiliations: () Dept. of Psychology, Uniersity of Pennsylania 34 Walnut Street 33C Philadelphia, PA 94-68 U.S.A. () Howard
More informationIntroduction to Thermodynamic Cycles Part 1 1 st Law of Thermodynamics and Gas Power Cycles
Introduction to Thermodynamic Cycles Part 1 1 st Law of Thermodynamics and Gas Power Cycles by James Doane, PhD, PE Contents 1.0 Course Oeriew... 4.0 Basic Concepts of Thermodynamics... 4.1 Temperature
More informationSEG/San Antonio 2007 Annual Meeting
Characterization and remoal of errors due to magnetic anomalies in directional drilling Nathan Hancock * and Yaoguo Li Center for Graity, Electrical, and Magnetic studies, Department of Geophysics, Colorado
More informationResidual migration in VTI media using anisotropy continuation
Stanford Exploration Project, Report SERGEY, Noember 9, 2000, pages 671?? Residual migration in VTI media using anisotropy continuation Tariq Alkhalifah Sergey Fomel 1 ABSTRACT We introduce anisotropy
More informationFeb 6, 2013 PHYSICS I Lecture 5
95.141 Feb 6, 213 PHYSICS I Lecture 5 Course website: faculty.uml.edu/pchowdhury/95.141/ www.masteringphysics.com Course: UML95141SPRING213 Lecture Capture h"p://echo36.uml.edu/chowdhury213/physics1spring.html
More informationReal Gas Thermodynamics. and the isentropic behavior of substances. P. Nederstigt
Real Gas Thermodynamics and the isentropic behaior of substances. Nederstigt ii Real Gas Thermodynamics and the isentropic behaior of substances by. Nederstigt in partial fulfillment of the requirements
More informationTailing of the breakthrough curve in aquifer contaminant transport: equivalent longitudinal macrodispersivity and occurrence of anomalous transport
5 GQ7: Securing Groundwater Quality in Urban and Industrial Environments (Proc. 6th International Groundwater Quality Conference held in Fremantle, Western Australia, 2 7 December 27). IAHS Publ. no. XXX,
More informationStresses in Sedimentary Strata, including Coals, and the Effects of Fluid Withdrawal on Effective Stress and Permeability
Uniersity of Wollongong Research Online Coal Operators' Conference Faculty of Engineering and Information Sciences 2011 Stresses in Sedimentary Strata, including Coals, and the Effects of Fluid Withdrawal
More informationE : Ground-penetrating radar (GPR)
Geophysics 3 March 009 E : Ground-penetrating radar (GPR) The EM methods in section D use low frequency signals that trael in the Earth by diffusion. These methods can image resistiity of the Earth on
More informationNotes on Linear Minimum Mean Square Error Estimators
Notes on Linear Minimum Mean Square Error Estimators Ça gatay Candan January, 0 Abstract Some connections between linear minimum mean square error estimators, maximum output SNR filters and the least square
More informationAn Introduction to Three-Dimensional, Rigid Body Dynamics. James W. Kamman, PhD. Volume II: Kinetics. Unit 3
Summary An Introduction to hree-dimensional, igid ody Dynamics James W. Kamman, PhD Volume II: Kinetics Unit Degrees of Freedom, Partial Velocities and Generalized Forces his unit defines the concepts
More informationPurpose of the experiment
Impulse and Momentum PES 116 Adanced Physics Lab I Purpose of the experiment Measure a cart s momentum change and compare to the impulse it receies. Compare aerage and peak forces in impulses. To put the
More informationAstrometric Errors Correlated Strongly Across Multiple SIRTF Images
Astrometric Errors Correlated Strongly Across Multiple SIRTF Images John Fowler 28 March 23 The possibility exists that after pointing transfer has been performed for each BCD (i.e. a calibrated image
More informationNote: the net distance along the path is a scalar quantity its direction is not important so the average speed is also a scalar.
PHY 309 K. Solutions for the first mid-term test /13/014). Problem #1: By definition, aerage speed net distance along the path of motion time. 1) ote: the net distance along the path is a scalar quantity
More informationarxiv: v3 [physics.data-an] 9 Oct 2013
Uncoering wind turbine properties through two-dimensional stochastic modeling of wind dynamics arxi:1210.71613 [physics.data-an] 9 Oct 2013 Frank Raischel, 1, 2 Teresa Scholz, 3, 4 Vitor V. Lopes, 3 and
More informationLecture 21: Physical Brownian Motion II
Lecture 21: Physical Brownian Motion II Scribe: Ken Kamrin Department of Mathematics, MIT May 3, 25 Resources An instructie applet illustrating physical Brownian motion can be found at: http://www.phy.ntnu.edu.tw/jaa/gas2d/gas2d.html
More informationRelativity in Classical Mechanics: Momentum, Energy and the Third Law
Relatiity in Classical Mechanics: Momentum, Energy and the Third Law Roberto Assumpção, PUC-Minas, Poços de Caldas- MG 37701-355, Brasil assumpcao@pucpcaldas.br Abstract Most of the logical objections
More informationClassical dynamics on graphs
Physical Reiew E 63 (21) 66215 (22 pages) Classical dynamics on graphs F. Barra and P. Gaspard Center for Nonlinear Phenomena and Complex Systems, Uniersité Libre de Bruxelles, Campus Plaine C.P. 231,
More informationNoise constrained least mean absolute third algorithm
Noise constrained least mean absolute third algorithm Sihai GUAN 1 Zhi LI 1 Abstract: he learning speed of an adaptie algorithm can be improed by properly constraining the cost function of the adaptie
More informationWhy does Saturn have many tiny rings?
2004 Thierry De Mees hy does Saturn hae many tiny rings? or Cassini-Huygens Mission: New eidence for the Graitational Theory with Dual Vector Field T. De Mees - thierrydemees @ pandora.be Abstract This
More informationLect-19. In this lecture...
19 1 In this lecture... Helmholtz and Gibb s functions Legendre transformations Thermodynamic potentials The Maxwell relations The ideal gas equation of state Compressibility factor Other equations of
More informationState-space Modelling of Hysteresis-based Control Schemes
European Control Conference (ECC) July 7-9,, Zürich, Switzerland. State-space Modelling of Hysteresis-based Control Schemes Soumya Kundu Ian A. Hiskens Abstract The paper deelops a state-space model for
More informationESCI 485 Air/sea Interaction Lesson 3 The Surface Layer
ESCI 485 Air/sea Interaction Lesson 3 he Surface Layer References: Air-sea Interaction: Laws and Mechanisms, Csanady Structure of the Atmospheric Boundary Layer, Sorbjan HE PLANEARY BOUNDARY LAYER he atmospheric
More informationContinuous Time Random Walks for Non-Local Radial Solute Transport
Continuous Time Random Walks for Non-Local Radial Solute Transport Marco Dentz Institute of Environmental Assessment and Water Research (IDÆA), Spanish National Research Council (CSIC), 08034 Barcelona,
More informationClassical Mechanics NEWTONIAN SYSTEM OF PARTICLES MISN NEWTONIAN SYSTEM OF PARTICLES by C. P. Frahm
MISN-0-494 NEWTONIAN SYSTEM OF PARTICLES Classical Mechanics NEWTONIAN SYSTEM OF PARTICLES by C. P. Frahm 1. Introduction.............................................. 1 2. Procedures................................................
More informationare applied to ensure that physical principles are not iolated in the definition of the discrete transition model. The oerall goal is to use this fram
A Comprehensie Methodology for Building Hybrid Models of Physical Systems Pieter J. Mosterman Λ Institute of Robotics and System Dynamics DLR Oberpfaffenhofen P.O. Box 1116 D-8223 Wessling Germany Gautam
More informationSupplementary Information Microfluidic quadrupole and floating concentration gradient Mohammad A. Qasaimeh, Thomas Gervais, and David Juncker
Mohammad A. Qasaimeh, Thomas Gerais, and Daid Juncker Supplementary Figure S1 The microfluidic quadrupole (MQ is modeled as two source (Q inj and two drain (Q asp points arranged in the classical quardupolar
More informationAn Elastic Contact Problem with Normal Compliance and Memory Term
Machine Dynamics Research 2012, Vol. 36, No 1, 15 25 Abstract An Elastic Contact Problem with Normal Compliance and Memory Term Mikäel Barboteu, Ahmad Ramadan, Mircea Sofonea Uniersité de Perpignan, France
More informationGeneral Lorentz Boost Transformations, Acting on Some Important Physical Quantities
General Lorentz Boost Transformations, Acting on Some Important Physical Quantities We are interested in transforming measurements made in a reference frame O into measurements of the same quantities as
More informationThe Production of Interactive Software for Supporting the Kinematics Study on Linear Motion and Swing Pendulum
The Production of Interactie Software for Supporting the Kinematics Study on Linear Motion and Swing Pendulum Liliana, Kartika Gunadi, Yonathan Rindayanto Ongko Informatics Engineering, Petra Christian
More informationTrajectory Estimation for Tactical Ballistic Missiles in Terminal Phase Using On-line Input Estimator
Proc. Natl. Sci. Counc. ROC(A) Vol. 23, No. 5, 1999. pp. 644-653 Trajectory Estimation for Tactical Ballistic Missiles in Terminal Phase Using On-line Input Estimator SOU-CHEN LEE, YU-CHAO HUANG, AND CHENG-YU
More informationStatus: Unit 2, Chapter 3
1 Status: Unit, Chapter 3 Vectors and Scalars Addition of Vectors Graphical Methods Subtraction of Vectors, and Multiplication by a Scalar Adding Vectors by Components Unit Vectors Vector Kinematics Projectile
More informationA possible mechanism to explain wave-particle duality L D HOWE No current affiliation PACS Numbers: r, w, k
A possible mechanism to explain wae-particle duality L D HOWE No current affiliation PACS Numbers: 0.50.-r, 03.65.-w, 05.60.-k Abstract The relationship between light speed energy and the kinetic energy
More informationv v Downloaded 01/11/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
The pseudo-analytical method: application of pseudo-laplacians to acoustic and acoustic anisotropic wae propagation John T. Etgen* and Serre Brandsberg-Dahl Summary We generalize the pseudo-spectral method
More informationarxiv: v1 [physics.comp-ph] 17 Jan 2014
An efficient method for soling a correlated multi-item inentory system Chang-Yong Lee and Dongu Lee The Department of Industrial & Systems Engineering, Kongu National Uniersity, Kongu 34-70 South Korea
More informationOn the Connections Between Laplace and ELzaki Transforms
Adances in Theoretical and Applied Mathematics. ISSN 973-4554 Volume 6, Number (), pp. Research India Publications http://www.ripublication.com/atam.htm On the Connections Between Laplace and ELzaki Transforms
More informationSIMULATIONS OF CHARACTERISTICS OF TUNED LIQUID COLUMN DAMPER USING AN ELLIPTICAL FLOW PATH ESTIMATION METHOD
October -7, 008, Beijing, China SIMULATIONS OF CHARACTERISTICS OF TUNED LIQUID COLUMN DAMPER USING AN ELLIPTICAL FLOW PATH ESTIMATION METHOD P. Chaiiriyawong, S. Limkatanyu and T. Pinkaew 3 Lecturer, Dept.
More informationME224 Lab 5 - Thermal Diffusion
ME4 Lab 5 ME4 Lab 5 - hermal Diffusion (his lab is adapted from IBM-PC in the laboratory by B G homson & A F Kuckes, Chapter 5) 1. Introduction he experiments which you will be called upon to do in this
More informationMagnetic Fields Part 3: Electromagnetic Induction
Magnetic Fields Part 3: Electromagnetic Induction Last modified: 15/12/2017 Contents Links Electromagnetic Induction Induced EMF Induced Current Induction & Magnetic Flux Magnetic Flux Change in Flux Faraday
More informationAnalysis of cylindrical heat pipes incorporating the e ects of liquid±vapor coupling and non-darcian transportða closed form solution
International Journal of Heat and Mass Transfer 42 (1999) 3405±341 Analysis of cylindrical heat pipes incorporating the e ects of liquid±apor coupling and non-darcian transportða closed form solution N.
More informationSubordinated advection-dispersion equation for contaminant transport
WATER RESOURCES RESEARCH, VO. 37, NO. 6, PAGES 1543 1550, JUNE 2001 Subordinated advection-dispersion equation for contaminant transport Boris Baeumer, 1 David A. Benson, 2 Mark M. Meerschaert, 3 and Stephen
More informationThe Kinetic Theory of Gases
978-1-107-1788-3 Classical and Quantum Thermal Physics The Kinetic Theory of Gases CHAPTER 1 1.0 Kinetic Theory, Classical and Quantum Thermodynamics Two important components of the unierse are: the matter
More informationA continuous time random walk approach to transient flow in heterogeneous porous media
WATER RESOURCES RESEARCH, VOL. 42,, doi:10.1029/2006wr005227, 2006 A continuous time random walk approach to transient flow in heterogeneous porous media Andrea Cortis 1 and Christen Knudby 2 Received
More informationReversal in time order of interactive events: Collision of inclined rods
Reersal in time order of interactie eents: Collision of inclined rods Published in The European Journal of Physics Eur. J. Phys. 27 819-824 http://www.iop.org/ej/abstract/0143-0807/27/4/013 Chandru Iyer
More informationChapter 1: Kinematics of Particles
Chapter 1: Kinematics of Particles 1.1 INTRODUCTION Mechanics the state of rest of motion of bodies subjected to the action of forces Static equilibrium of a body that is either at rest or moes with constant
More informationCOMPARISON OF ANALYTICAL SOLUTIONS FOR CMSMPR CRYSTALLIZER WITH QMOM POPULATION BALANCE MODELING IN FLUENT
CHINA PARTICUOLOGY Vol. 3, No. 4, 13-18, 5 COMPARISON OF ANALYTICAL SOLUTIONS FOR CMSMPR CRYSTALLIZER WITH QMOM POPULATION BALANCE MODELING IN FLUENT Bin Wan 1, Terry A. Ring 1, *, Kumar M. Dhanasekharan
More informationTransmission Line Transients
8 5 Transmission Line Transients CHAPTER OBJECTIES After reading this chapter, you should be able to: Proide an analysis of traelling waes on transmission lines Derie a wae equation Understand the effect
More informationDynamics ( 동역학 ) Ch.2 Motion of Translating Bodies (2.1 & 2.2)
Dynamics ( 동역학 ) Ch. Motion of Translating Bodies (. &.) Motion of Translating Bodies This chapter is usually referred to as Kinematics of Particles. Particles: In dynamics, a particle is a body without
More informationInformation Propagation Speed in Bidirectional Vehicular Delay Tolerant Networks
Information Propagation Speed in Bidirectional Vehicular Delay Tolerant Networks Emmanuel Baccelli, Philippe Jacquet, Bernard Mans, Georgios Rodolakis INRIA, France, {firstname.lastname}@inria.fr Macquarie
More informationTransmission lines using a distributed equivalent circuit
Cambridge Uniersity Press 978-1-107-02600-1 - Transmission Lines Equialent Circuits, Electromagnetic Theory, and Photons Part 1 Transmission lines using a distributed equialent circuit in this web serice
More informationOn Using Random Walks to Solve the Space-Fractional Advection-Dispersion Equations
Journal of Statistical Physics, Vol. 123, No. 1, April 2006 ( C 2006 ) DOI: 10.1007/s10955-006-9042-x On Using Random Walks to Solve the Space-Fractional Advection-Dispersion Equations Zhang Yong, 1,2
More informationSolution to 1-D consolidation of non-homogeneous soft clay *
Xie et al. / J Zhejiang Uni SCI 25 6A(Suppl. I):29-34 29 Journal of Zhejiang Uniersity SCIENCE ISSN 19-395 http://www.zju.edu.cn/jzus E-mail: jzus@zju.edu.cn Solution to 1-D consolidation of non-homogeneous
More informationLECTURE NOTE THERMODYNAMICS (GEC 221)
LETURE NOTE ON THERMODYNAMIS (GE ) Thermodynamics is the branch of science that treats the arious phenomena of energy and related properties of matter especially the relationship between heat, work and
More informationERAD THE SEVENTH EUROPEAN CONFERENCE ON RADAR IN METEOROLOGY AND HYDROLOGY
Multi-beam raindrop size distribution retrieals on the oppler spectra Christine Unal Geoscience and Remote Sensing, TU-elft Climate Institute, Steinweg 1, 68 CN elft, Netherlands, c.m.h.unal@tudelft.nl
More informationarxiv: v1 [stat.ml] 15 Feb 2018
1 : A New Algorithm for Streaming PCA arxi:1802.054471 [stat.ml] 15 Feb 2018 Puyudi Yang, Cho-Jui Hsieh, Jane-Ling Wang Uniersity of California, Dais pydyang, chohsieh, janelwang@ucdais.edu Abstract In
More informationDynamic Vehicle Routing with Moving Demands Part II: High speed demands or low arrival rates
ACC 9, Submitted St. Louis, MO Dynamic Vehicle Routing with Moing Demands Part II: High speed demands or low arrial rates Stephen L. Smith Shaunak D. Bopardikar Francesco Bullo João P. Hespanha Abstract
More information4-vectors. Chapter Definition of 4-vectors
Chapter 12 4-ectors Copyright 2004 by Daid Morin, morin@physics.harard.edu We now come to a ery powerful concept in relatiity, namely that of 4-ectors. Although it is possible to derie eerything in special
More informationSemi-implicit Treatment of the Hall Effect in NIMROD Simulations
Semi-implicit Treatment of the Hall Effect in NIMROD Simulations H. Tian and C. R. Soinec Department of Engineering Physics, Uniersity of Wisconsin-Madison Madison, WI 5376 Presented at the 45th Annual
More informationAnalytical solutions for compartmental models of contaminant transport in enclosed. spaces
ESGI8 Analytical solutions for compartmental models of contaminant transport in enclosed Problem presented by Simon Parker DSTL Executie Summary Understanding the transport of hazardous airborne materials
More informationA Regularization Framework for Learning from Graph Data
A Regularization Framework for Learning from Graph Data Dengyong Zhou Max Planck Institute for Biological Cybernetics Spemannstr. 38, 7076 Tuebingen, Germany Bernhard Schölkopf Max Planck Institute for
More informationWATER RESOURCES RESEARCH, VOL. 44, W05410, doi: /2008wr006858, 2008
Click Here for Full Article WATER RESOURCES RESEARCH, VOL. 44, W05410, doi:10.1029/2008wr006858, 2008 Transport of conservative solutes in simulated fracture networks: 2. Ensemble solute transport and
More informationTo string together six theorems of physics by Pythagoras theorem
To string together six theorems of physics by Pythagoras theorem H. Y. Cui Department of Applied Physics Beijing Uniersity of Aeronautics and Astronautics Beijing, 00083, China ( May, 8, 2002 ) Abstract
More informationModeling Hydraulic Accumulators for use in Wind Turbines
The 13th Scandinaian International Conference on Fluid Power, SICFP213, une 3-5, 213, Linköping, Sweden Modeling Hydraulic Accumulators for use in Wind Turbines Henrik Brun Hansen and Peter Windfeld Rasmussen
More informationDynamics of a Bertrand Game with Heterogeneous Players and Different Delay Structures. 1 Introduction
ISSN 749-889 (print), 749-897 (online) International Journal of Nonlinear Science Vol.5(8) No., pp.9-8 Dynamics of a Bertrand Game with Heterogeneous Players and Different Delay Structures Huanhuan Liu,
More informationSemi-Empirical Formulas of Drag/Lift Coefficients for High Speed Rigid Body Manoeuvring in Water Column
Semi-Empirical Formulas of Drag/Lift Coefficients for High Speed Rigid Body Manoeuring in Water Column Peter C. Chu and Chenwu Fan Naal Ocean Analysis and Prediction Laboratory Naal Postgraduate School,
More informationA Guidance Law for a Mobile Robot for Coverage Applications: A Limited Information Approach
Proceedings of the Third International Conference on Adances in Control and Optimization of Dynamical Systems, Indian Institute of Technology anpur, India, March 1-15, 014 Frcsh1.6 A Guidance Law for a
More informationTarget Trajectory Estimation within a Sensor Network
Target Trajectory Estimation within a Sensor Network Adrien Ickowicz IRISA/CNRS, 354, Rennes, J-Pierre Le Cadre, IRISA/CNRS,354, Rennes,France Abstract This paper deals with the estimation of the trajectory
More informationLesson 3: Free fall, Vectors, Motion in a plane (sections )
Lesson 3: Free fall, Vectors, Motion in a plane (sections.6-3.5) Last time we looked at position s. time and acceleration s. time graphs. Since the instantaneous elocit is lim t 0 t the (instantaneous)
More informationTools for Investigation of Dynamics of DC-DC Converters within Matlab/Simulink
Tools for Inestigation of Dynamics of DD onerters within Matlab/Simulink Riga Technical Uniersity, Riga, Latia Email: pikulin03@inbox.l Dmitry Pikulin Abstract: In this paper the study of complex phenomenon
More informationInternational Journal of Scientific & Engineering Research, Volume 7, Issue 11, November-2016 ISSN Nullity of Expanded Smith graphs
85 Nullity of Expanded Smith graphs Usha Sharma and Renu Naresh Department of Mathematics and Statistics, Banasthali Uniersity, Banasthali Rajasthan : usha.sharma94@yahoo.com, : renunaresh@gmail.com Abstract
More information4. A Physical Model for an Electron with Angular Momentum. An Electron in a Bohr Orbit. The Quantum Magnet Resulting from Orbital Motion.
4. A Physical Model for an Electron with Angular Momentum. An Electron in a Bohr Orbit. The Quantum Magnet Resulting from Orbital Motion. We now hae deeloped a ector model that allows the ready isualization
More informationThe dispersion of a light solid particle in high-reynolds number homogeneous stationary turbulence: LES approach with stochastic sub-grid model
Computational Methods in Multiphase Flow III 65 The dispersion of a light solid particle in high-reynolds number homogeneous stationary turbulence: ES approach with stochastic sub-grid model M. Gorokhoski
More informationNUMERICAL SIMULATION OF HYDRODYNAMIC FIELD FROM PUMP-TURBINE RUNNER
UMERICAL SIMULATIO OF YDRODYAMIC FIELD FROM PUMP-TURBIE RUER A. Iosif and I. Sârbu Department of Building Serices, Politehnica Uniersity of Timisoara, Romania E-Mail: ioan.sarbu@ct.upt.ro STRACT One of
More informationEE243 Advanced Electromagnetic Theory Lec # 15 Plasmons. Reading: These PPT notes and Harrington on corrugated surface 4.8.
EE Applied EM Fall 6, Neureuther Lecture #5 Ver //6 EE43 Adanced Electromagnetic Theory Lec # 5 Plasmons Fixed Slides 6, (Group and Phase Velocity for Guided Waes) (Sign Reersals for Propagation Direction
More informationA014 Uncertainty Analysis of Velocity to Resistivity Transforms for Near Field Exploration
A014 Uncertainty Analysis of Velocity to Resistiity Transforms for Near Field Exploration D. Werthmüller* (Uniersity of Edinburgh), A.M. Ziolkowski (Uniersity of Edinburgh) & D.A. Wright (Uniersity of
More informationProbabilistic Engineering Design
Probabilistic Engineering Design Chapter One Introduction 1 Introduction This chapter discusses the basics of probabilistic engineering design. Its tutorial-style is designed to allow the reader to quickly
More informationInelastic Collapse in One Dimensional Systems with Ordered Collisions
Lake Forest College Lake Forest College Publications Senior Theses Student Publications 4-4-05 Inelastic Collapse in One Dimensional Systems with Ordered Collisions Brandon R. Bauerly Lake Forest College,
More informationExperimental link of coarsening rate and volume distribution in dry foam
August 2012 EPL, 99 (2012 48003 doi:.1209/0295-5075/99/48003 www.epljournal.org Eperimental link of coarsening rate and olume distribution in dry foam J. Lambert 1, F. Graner 2,3, R. Delannay 1 and I.
More informationReduction of friction by normal oscillations. I. Influence of contact stiffness.
Reduction of friction by normal oscillations. I. Influence of contact stiffness. M. Popo,,3, V. L. Popo,,3, N. V. Popo Berlin Uniersity of Technology, 63 Berlin, Germany Tomsk Polytechnic Uniersity, 6345
More informationDamage detection in a sandwich composite beam using wavelet transforms
Damage detection in a sandwich composite beam using waelet transforms T. A. Dawood *a, R. A. Shenoi a, S. M. Veres a, M. Sahin a, M. J Gunning b a School of Engineering Sciences, Uniersity of Southampton,
More informationAn Alternative Characterization of Hidden Regular Variation in Joint Tail Modeling
An Alternatie Characterization of Hidden Regular Variation in Joint Tail Modeling Grant B. Weller Daniel S. Cooley Department of Statistics, Colorado State Uniersity, Fort Collins, CO USA May 22, 212 Abstract
More informationPlasma electron hole kinematics. I. Momentum conservation
Plasma electron hole kinematics. I. Momentum conseration The MIT Faculty has made this article openly aailable. Please share how this access benefits you. Your story matters. Citation As Published Publisher
More informationStationary states and energy cascades in inelastic gases
Los Alamos National Laboratory From the SelectedWorks of Eli Ben-Naim May, 005 Stationary states and energy cascades in inelastic gases E. Ben-Naim, Los Alamos National Laboratory J. Machta, Uniersity
More informationEVALUATION OF MOISTURE DIFFUSION THEORIES IN POROUS MATERIALS. By Juan C. Alvarez. MASTER OF SCIENCE in MECHANICAL ENGINEERING APPROVED:
EVALUAION OF OISURE DIFFUSION HEORIES IN POROUS AERIALS By Juan C. Alarez hesis submitted to the Faculty of the Virginia Polytechnic Institute and State Uniersity In partial fulfillment of the requirements
More informationSimulations of bulk phases. Periodic boundaries. Cubic boxes
Simulations of bulk phases ChE210D Today's lecture: considerations for setting up and running simulations of bulk, isotropic phases (e.g., liquids and gases) Periodic boundaries Cubic boxes In simulations
More informationPhysics 111. Help sessions meet Sunday, 6:30-7:30 pm in CLIR Wednesday, 8-9 pm in NSC 098/099
ics Announcements day, ember 7, 2007 Ch 2: graphing - elocity s time graphs - acceleration s time graphs motion diagrams - acceleration Free Fall Kinematic Equations Structured Approach to Problem Soling
More information