A note on groundwater flow along a hillslope
|
|
- Douglas Gallagher
- 5 years ago
- Views:
Transcription
1 WATER RESOURCES RESEARCH, VOL. 40,, oi: /2003wr002438, 2004 A note on grounwater flow along a hillslope Eoaro Daly Dipartimento i Iraulica Trasporti e Infrastrutture Civili, Politecnico i Torino, Torino, Italy Amilcare Porporato Department of Civil an Environmental Engineering, Duke University, Durham, North Carolina, USA Receive 28 June 2003; revise 2 October 2003; accepte 22 October 2003; publishe 6 January [1] Particular cases of grounwater flow along a hillslope are stuie using the grounwater hyraulic theory (e.g., Dupuit approximation) escribe by the Boussinesq equation. Analytical solutions are foun as simple transformations of known similarity solutions. In particular, the expressions for phreatic surface level an volume flux an the corresponing loop-rating curve for the flow of a grounwater moun as well as the travelling wave solutions of the Boussinesq equation are escribe in etail. INDEX TERMS: 1829 Hyrology: Grounwater hyrology; 1899 Hyrology: General or miscellaneous; 3220 Mathematical Geophysics: Nonlinear ynamics; KEYWORDS: Boussinesq equation, Dupuit approximation, similarity solution Citation: Daly, E., an A. Porporato (2004), A note on grounwater flow along a hillslope, Water Resour. Res., 40,, oi: /2003wr Introuction [2] Unconfine grounwater flow in a sloping aquifer may be moele using the nonlinear Boussinesq equation, base on the classical Dupuit approximation. Solutions of such an equation are of interest in grounwater catchment hyrology [Troch et al., 1993; Sanfor et al., 1993; Szilagyi an Parlange, 1998; Mizimura, 2002], coastal grounwater hyraulics [e.g., Li et al., 2000a, 2000b], an several other nonlinear iffusion problems [e.g., Peletier, 1971; Aronson, 1986; Gratton an Minotti, 1990; Diez et al., 1992]. Although exact analytical solutions of the Boussinesq equation on a sloping be o not seem to have been given before, many approximate solutions have been propose in the literature [e.g., Lockington et al., 2000; Telyakovskiy et al., 2002], among which the solutions of the linearize equation [e.g., Brutsaert, 1994] an of the kinematic wave equation [e.g., Beven, 1981; Fan an Bras, 1998; Troch et al., 2002] must also be quote. [3] This paper presents some analytical solutions of the Boussinesq equation escribing the grounwater flow along a hillslope. By means of a simple travelling wave coorinate transformation, the Boussinesq equation is written as if the flow occurre on a horizontal impermeable be, for which exact similarity solutions are well known [e.g., Polubarinova- Kochina, 1962; Aravin an Numerov, 1965; Bear, 1972; Gratton an Minotti, 1990]. The practical applicability of the approach is somewhat limite by the fact that the same transformation must also apply to the initial an bounary conitions; however, the analysis may be useful to clarify some mathematical an physical aspects of the problem an to furnish benchmarks for the valiation of numerical simulations. 2. Moel Description [4] The water flow within a shallow saturate zone of a soil mantle overlying an impermeable berock of slope q can be stuie following the classical approach of Dupuit [e.g., Polubarinova-Kochina, 1962; Bear, 1972], which neglects the curvature of the flow streamlines, which are assume to be parallel to the impermeable be. As iscusse by Dagan [1967], such an approximation is equivalent to the shallow-water approximation in open channel flows [e.g., Stoker, 1957]. [5] For isotropic an homogeneous soils, the continuity equation an the equation of motion are cos q þ g þ 1 ¼ sin q ð2þ where t is time, x is the spatial coorinate, n is the average soil porosity, h is the water level, Q is the time volume flux per unit with, i.e., Q = hu, u is the seepage velocity, g is the gravity constant, sinq is the be slope, an j is the friction slope. Neglecting the inertial terms an assuming j = u/k s, where K s is the average soil hyraulic conuctivity, the equation of motion becomes u ¼ K s sin q K s cos q; ð3þ which is the well-known Darcy s law. The introuction of equation (3) in equation (1) leas to ð1þ Copyright 2004 by the American Geophysical Union /04/2003WR ¼ K s sin q þ K h cos q; ð4þ 1of5
2 DALY AND PORPORATO: TECHNICAL NOTE that is usually referre to as the Boussinesq equation for the general case of an unconfine aquifer on a sloping be. A few exact similarity solutions of equation (4) have been erive for horizontal be, i.e., sin q = 0 [e.g., Polubarinova- Kochina, 1962; Aravin an Numerov, 1965; Bear, 1972; Gratton an Minotti, 1990]. [6] A common simplification of equation (2) is to assume that the hyraulic graient is equal to the be slope [e.g., Beven, 1981]; with this hypothesis, equation (4) is simplifie to a linear kinematic wave equation, whose solutions are waves moving ownhill at a constant spee with unaltere shape. In fact, the presence of linear convection suggests the use of a travelling wave coorinate also to transform equation (4). Thus introucing a new system of coorinates whose origin moves with the velocity c 0 = K s sin q/n, i.e., equation (4) becomes x ¼ x c 0 t; t ¼ t; ¼ h ; ð6þ where K is equal to K s cos q/n. Obviously, the initial an bounary conitions associate to equation (4) must be transforme too. It is interesting to note that c 0 represents the effective mean velocity in the soil pores when / =0, i.e., c 0 = u/n, so that the origin of the new coorinate system follows the movement of the water along the hillslope. 3. Grounwater Moun [7] The evolution in time an space of a given water volume per unit with, V 0, initially concentrate in the section x = 0 of a completely ry soil can be escribe by equation (4), with the conition Z þ1 hx; ð tþx ¼ V 0 =n: ð7þ Clearly, thanks to the transformation (5), such a problem is equivalent to the classical case of a grounwater moun spreaing on a horizontal be. Following Barenblatt [1996], the similarity solution can be elegantly foun by imensional analysis. The level of the phreatic surface h epens on the governing parameters, t, x, V 0 /n, an K. Using the characteristic length in the irection orthogonal to the be, H, the length in the flow irection, L, an a characteristic timescale, T, the imensions of h an of the governing parameters can be expresse as ½hŠ ¼ H; ½tŠ ¼ T; ½Š¼L; x ½V 0 =nš ¼ HL; ½KŠ ¼ L 2 T H : ð8þ Dimensional analysis leas to the relation c = (h), where h c ¼ ðv 0 =nþ 2=3 K =3 t =3 h ¼ x ðv 0 =nþ 1=3 K 1=3 t 1=3 : ð9þ The substitution of (9) in (6) leas to the orinary ifferential equation h h þ h ð10þ 3 which, integrate twice (the first integration constant is zero by symmetry), yiels c ¼ ðhþ ¼C 1 h2 6 ; ð11þ where C 1 is equal to (3/32) 1/3 for the conition (7). [8] Returning to the initial coorinate system, the equation escribing the evolution of the phreatic surface along the hillslope of the initial water volume V 0 thus reas hx; ð tþ ¼ 3 1=3 V 2= ðktþ ð x c 0tÞ 2 1=3 6Kt ; ð12þ which escribes a parabola that moves ownhill with velocity c 0 expaning its base an reucing its peak. As the be slope, q, increases so oes c 0, while the spreaing of the moun base, which is proportional to (cos q) 1/3, becomes slower; as a consequence, the water volume being constant, the peak of the moun also ecays more slowly, i.e., proportionally to (cos q) /3. It shoul be notice, moreover, that the previous solution is not limite to small be slopes, but only by the Dupuit approximation (i.e., small curvatures of the streamlines). [9] The corresponing expression of the ischarge is Qx; ð tþ ¼ " # 3 1=3 V 2= ðktþ ðx c 2 0tÞ 1=3 6Kt nx þ 2K s sin q t : ð13þ 3t As appears from Figure 1a, uring the initial perio the water front moves uphill, leaing to a negative flux in the sections upstream x = 0. From equation (12), the position of the upper front moves in time accoring to " x u ¼ 1 2n 2K s sin q t 6 2=3 n V # 0K s cos q t 1=3 n 2 : ð14þ The uppermost point, x u,m, is touche by the front at pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi t u ¼ nv 0 tan q sin q pffiffi ; 6 Ks ð15þ as shown in the example of Figure 2. [10] The combine analytical expressions of h an Q (equations (12) an (13)) escribe the so-calle loop-rating curve [e.g., Henerson, 1966], which summarizes the temporal ynamics of the grounwater wave in a given section (Figure 3). [11] The extension to the case with cylinrical symmetry is also of interest. Assuming an isotropic an homogeneous soil, the two-imensional continuity equation an the two equations of motion along a hillslope are ð hu xþ hu ð16þ 2of5
3 DALY AND PORPORATO: TECHNICAL NOTE Using cylinrical coorinates, x = r cos f an z = r sin f, equation (22) may be written ¼ 1 2 h 2 þ 1 2 þ 2 h 2 r 2 ; ð23þ where the last term on the right han sie is zero by raial symmetry. [12] Being now [V 0 /n] =HL 2, imensional analysis leas to c = (r), where h c ¼ ðv 0 =nþ 1=2 K =2 t =2 r r ¼ ðv 0 =nþ 1=4 K 1=4 t ; 1=4 ð24þ which substitute in equation (22) gives the orinary ifferential equation whose solution is r 2 r r þ r2 2 ð25þ Figure 1. Flow of a grounwater moun along a hillslope: (a) phreatic surface an (b) water flux per unit with. Parameters: n = 0.4, K s =0.5m, q =15, V 0 =1m 2 (initial time 0.2 ; final time 3.8 ; time step 0.4 ). u x ¼ K s sin q K s cos q; ð17þ u y ¼ K s cos ð18þ where u x an u y are the seepage velocity along respectively x an y. Substituting the two expressions of the velocities in equation (16), the two-imensional Boussinesq equation is obtaine as ¼ K s sin q þ K s cos : ð19þ ¼ C 2 r2 8 ; ð26þ p where C 2 is equal to 1/(2 ffiffiffi p ), as follows from equation (20). [13] Therefore, in terms of the original coorinates system, the evolution of the phreatic surface, hx; ð y; tþ ¼ V " # 1=2 0=n 1 p Kt 2 ffiffiffi ðx c 0tÞ 2 þ y 2 ; ð27þ p 8ðV 0 =nþ 1=2 ðktþ 1=2 represents a paraboloi that subtens a constant volume equal to V 0 /n, an that moves ownhill expaning its circular base an reucing its peak. 4. Traveling Waves [14] Since equation (4) is translationally invariant, it amits travelling wave solutions of type [e.g., Gratton an Minotti, 1990] h ¼ hðþ; x x ¼ x ct; ð28þ Introucing a water volume V 0 concentrate in the origin of a ry soil, such that Z þ1 Z þ1 hx; ð y; tþxy ¼ V 0 =n; ð20þ the substitution yiels x ¼ x c 0 t; z ¼ y; t ¼ t; ¼ : ð22þ 3of5 Figure 2. Space-time evolution of the upper front of a grounwater moun flowing along a hillslope. The two terms of equation (14) are also shown. Parameters as in Figure 1.
4 DALY AND PORPORATO: TECHNICAL NOTE Figure 3. Loop-rating curve showing the temporal evolution of the grounwater wave at section x =4 m. Parameters as in Figure 1. where c is a constant etermine by the bounary conitions. Such solutions are closely connecte with selfsimilarity [Barenblatt, 1996], an represent phreatic surface profiles which move along the hillslope without changing their shape. [15] Using equation (28), equation (4) may be written as ðc c 0 Þh þ Kh h x x which after a first integration leas to ð29þ ðc c 0 Þh þ Kh h x ¼ C 3: ð30þ The solution corresponing to C 3 = 0 is simply h ¼ h 0 þ c 0 c x; ð31þ K which escribes a straight line that avances with spee c with a front locate at x =(h 0 K)/(c c 0 ), where h 0 is the water level at x = 0 (Figure 4a). The flow might be imagine as prouce by a piston, whose axis is parallel to the be, pushing at constant spee a given volume of water [e.g., Gratton an Minotti, 1990]. [16] A ifferent solution is obtaine when C 3 6¼ 0. Letting C 3 =(c c 0 )h 0,itis x þ C 4 ¼ K ½ c 0 c h þ h 0 lnðh h 0 ÞŠ; ð32þ Figure 4. Traveling wave solutions efine by equation (30). (a) C 3 = 0 an (b) C 3 6¼ 0. Parameters: n = 0.4, K s = 0.5 m, q =45, h 0 =1m,C 4 chosen such that h(x =0) =3m. 4of5
5 DALY AND PORPORATO: TECHNICAL NOTE where C 4 is etermine choosing a value of h for a certain x. This case might represent the asymptotics of a flow prouce by a piston [Gratton an Minotti, 1990] pushing a layer of flui of initial thickness h 0 (Figure 4b) or, equivalently, a bounary conition h = h 0 moving with constant spee c. All the profiles for c < c 0 ten to h 0 when x!, while they are asymptotic to Kx/(c 0 c) for positive x (i.e., horizonal for c = 0). The situation is reverse for c > c 0. References Aravin, V. I., an S. N. Numerov (1965), Theory of Flui Flow in Uneformable Porous Meia, Isr. Program for Sci. Transl., Jerusalem. Aronson, D. G. (1986), The porous meium equation, in Some Problems in Nonlinear Diffusion, eite by A. Fasano an M. Primicerio, Lect. Notes Math., 1224, Barenblatt, G. I. (1996), Scaling, Self-Similarity, an Intermeiate Asimptotics, Cambrige Univ. Press, New York. Bear, J. (1972), Dynamics of Fluis in Porous Meia, Elsevier Sci., New York. Beven, K. (1981), Kinematic Subsurface Stormflow, Water Resour. Res., 17(5), Brutsaert, W. (1994), The unit response of grounwater outflow from a hillslope, Water Resour. Res., 30(10), Dagan, G. (1967), Secon orer theory of shallow free surface flow in porous meia, Q. J. Mech. Appl. Math., 20(4), Diez, J. A., R. Gratton, an J. Gratton (1992), Self-similar solution of the secon kin for a convergent viscous gravity current, Phys. Fluis, 4(6), Fan, Y., an R. L. Bras (1998), Analytical solutions to hillslope subsurface storm flow an saturation overlan flow, Water Resour. Res., 34(4), Gratton, J., an F. Minotti (1990), Self-similar viscous gravity currents: Phase-plane formalism, J. Flui Mech., 210, Henerson, F. M. (1966), Open Channel Flow, Macmillan, New York. Li, L., D. A. Barry, F. Stagnitti, J.-Y. Parlange, an D.-S. Jeng (2000a), Beach water table fluctuations ue to spring-neap ties: Moving bounary effects, Av. Water Resour., 23, Li, L., D. A. Barry, C. Cunningham, F. Stagnitti, an J.-Y. Parlange (2000b), A two imensional analytical solution of grounwater responses to tial loaing in an estuary an ocean, Av. Water Resour., 23, Lockington, D. A., J.-Y. Parlange, M. B. Parlange, an J. Selker (2000), Similarity solution of the Boussinesq equation, Av. Water Resour., 23, Mizimura, K. (2002), Drought flow from hillslope, J. Hyrol. Eng., 7(2), Peletier, L. A. (1971), Asymptotic behavior of solutions of the porous meia equation, SIAM J. Appl. Math., 21(4), Polubarinova-Kochina, P. Y. (1962), Theory of Grounwater Movement, Princeton Univ. Press, Princeton, N. J. Sanfor, W. E., J.-Y. Parlange, an T. S. Steenhuis (1993), Hillslope rainage with suen ryown: Close form solution an laboratory experiments, Water Resour. Res., 29(7), Stoker, J. J. (1957), Water Waves, Wiley-Interscience, Hoboken, N. J. Szilagyi, J., an M. B. Parlange (1998), Baseflow separation base on analytical solutions of the Boussinesq equation, J. Hyrol., 204, Telyakovskiy, A. S., G. A. Braga, an F. Furtao (2002), Approximate similarity solutions of the Boussinesq equation, Av. Water Resour., 25, Troch, P. A., F. P. De Troch, an W. Brutsaert (1993), Effective water table epth to escribe initial conitions prior to storm rainfall in humi regions, Water Resour. Res., 29(2), Troch, P., E. van Loon, an A. Hilberts (2002), Analytical solutions to hillslope kinematic wave equation for subsurface flow, Av. Water Resour., 25, E. Daly, Dipartimento i Iraulica Trasporti e Infrastrutture Civili, Politecnico i Torino, Torino, Italy. A. Porporato, Department of Civil an Environmental Engineering, Duke University, 127 Huson Hall, Durham, NC 27708, USA. (amilcare@ uke.eu.) 5of5
Chapter 2 Governing Equations
Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationA Short Note on Self-Similar Solution to Unconfined Flow in an Aquifer with Accretion
Open Journal o Flui Dynamics, 5, 5, 5-57 Publishe Online March 5 in SciRes. http://www.scirp.org/journal/oj http://x.oi.org/.46/oj.5.57 A Short Note on Sel-Similar Solution to Unconine Flow in an Aquier
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationDrainage of a horizontal Boussinesq aquifer with a power law hydraulic conductivity profile
WATER RESOURCES RESEARCH, VOL. 41, W11422, doi:1.129/25wr4241, 25 Drainage of a horizontal Boussinesq aquifer with a power law hydraulic conductivity profile David E. Rupp and John S. Selker Department
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationComparative Approaches of Calculation of the Back Water Curves in a Trapezoidal Channel with Weak Slope
Proceeings of the Worl Congress on Engineering Vol WCE, July 6-8,, Lonon, U.K. Comparative Approaches of Calculation of the Back Water Curves in a Trapezoial Channel with Weak Slope Fourar Ali, Chiremsel
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
More informationCompletely passive natural convection
Early View publication on wileyonlinelibrary.com (issue an page numbers not yet assigne; citable using Digital Object Ientifier DOI) ZAMM Z. Angew. Math. Mech., 1 6 (2011) / DOI 10.1002/zamm.201000030
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More informationEXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL
Tome XVI [18] Fascicule 3 [August] 1. Charles Chinwuba IKE EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL 1. Department of Civil Engineering, Enugu State University
More informationCan we distinguish Richards and Boussinesq s equations for hillslopes?: The Coweeta experiment revisited
WATER RESOURCES RESEARCH, VOL. 35, NO. 2, PAGES 589 593, FEBRUARY 1999 Can we distinguish Richards and Boussinesq s equations for hillslopes?: The Coweeta experiment revisited T. S. Steenhuis, 1 J.-Y.
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationTOWARDS THERMOELASTICITY OF FRACTAL MEDIA
ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationAn analytical investigation into filmwise condensation on a horizontal tube in a porous medium with suction at the tube surface
Heat Mass Transfer (29) 45:355 361 DOI 1.17/s231-8-436-y ORIGINAL An analytical investigation into filmwise conensation on a horizontal tube in a porous meium with suction at the tube surface Tong Bou
More informationfv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n
Chapter 11 Rossby waves Supplemental reaing: Pelosky 1 (1979), sections 3.1 3 11.1 Shallow water equations When consiering the general problem of linearize oscillations in a static, arbitrarily stratifie
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationConvective heat transfer
CHAPTER VIII Convective heat transfer The previous two chapters on issipative fluis were evote to flows ominate either by viscous effects (Chap. VI) or by convective motion (Chap. VII). In either case,
More informationApplications of First Order Equations
Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more
More informationHomework 7 Due 18 November at 6:00 pm
Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine
More informationMulti-component bi-hamiltonian Dirac integrable equations
Chaos, Solitons an Fractals 9 (009) 8 8 www.elsevier.com/locate/chaos Multi-component bi-hamiltonian Dirac integrable equations Wen-Xiu Ma * Department of Mathematics an Statistics, University of South
More informationarxiv: v1 [math-ph] 2 May 2016
NONLINEAR HEAT CONDUCTION EQUATIONS WITH MEMORY: PHYSICAL MEANING AND ANALYTICAL RESULTS PIETRO ARTALE HARRIS 1 AND ROBERTO GARRA arxiv:165.576v1 math-ph] May 16 Abstract. We stuy nonlinear heat conuction
More informationA simple model for the small-strain behaviour of soils
A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:
More informationu t v t v t c a u t b a v t u t v t b a
Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More informationSudden drawdown and drainage of a horizontal aquifer
WATER RESOURCES RESEARCH, VOL. 37, NO. 8, PAGES 2097 20, AUGUST 200 Sudden drawdown and drainage of a horizontal aquifer J.-Y. Parlange, F. Stagnitti, 2 A. Heilig, J. Szilagyi, 3 M. B. Parlange, 4 T. S.
More informationLecture 2 - First order linear PDEs and PDEs from physics
18.15 - Introuction to PEs, Fall 004 Prof. Gigliola Staffilani Lecture - First orer linear PEs an PEs from physics I mentione in the first class some basic PEs of first an secon orer. Toay we illustrate
More informationModelling the Zero-Inertia, Horizontal Viscous Dam-Break Problem
r WSEAS International Conference on APPLIED an TEORETICAL MECANICS, Spain, December 4-6, 7 8 Moelling the Zero-Inertia, orizontal Viscous Dam-Break Problem BLAISE NSOM, WILFRIED NDONG AND BLAISE RAVELO
More informationLagrangian and Hamiltonian Dynamics
Lagrangian an Hamiltonian Dynamics Volker Perlick (Lancaster University) Lecture 1 The Passage from Newtonian to Lagrangian Dynamics (Cockcroft Institute, 22 February 2010) Subjects covere Lecture 2: Discussion
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationarxiv: v1 [cond-mat.stat-mech] 9 Jan 2012
arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,
More informationPoroelasticity and tidal loading
Chapter 6: Poroelasticity an tial loaing 63 Chapter 6: Poroelasticity an tial loaing 6. Introuction It is establishe above, in Chapter 4, that there are many known examples of the tial moulation of seafloor
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationDusty Plasma Void Dynamics in Unmoving and Moving Flows
7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND SPACE SCIENCES (EUCASS) Dusty Plasma Voi Dynamics in Unmoving an Moving Flows O.V. Kravchenko*, O.A. Azarova**, an T.A. Lapushkina*** *Scientific an Technological
More informationAnalytic Scaling Formulas for Crossed Laser Acceleration in Vacuum
October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945
More informationOptimization of a point-mass walking model using direct collocation and sequential quadratic programming
Optimization of a point-mass walking moel using irect collocation an sequential quaratic programming Chris Dembia June 5, 5 Telescoping actuator y Stance leg Point-mass boy m (x,y) Swing leg x Leg uring
More informationDynamics of magmatic intrusions in the upper crust: Theory and applications to laccoliths on Earth and the Moon
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116,, oi:1.129/21jb818, 211 Dynamics of magmatic intrusions in the upper crust: Theory an applications to laccoliths on Earth an the Moon Chloé Michaut 1 Receive 18
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationApplication of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate
Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite
More informationRelation between the propagator matrix of geodesic deviation and the second-order derivatives of the characteristic function
Journal of Electromagnetic Waves an Applications 203 Vol. 27 No. 3 589 60 http://x.oi.org/0.080/0920507.203.808595 Relation between the propagator matrix of geoesic eviation an the secon-orer erivatives
More informationCouncil for Innovative Research
ISSN: 347-3487 A solution of Fractional Laplace's equation y Moifie separation of variales ABSTRACT Amir Pishkoo,, Maslina Darus, Fatemeh Tamizi 3 Physics an Accelerators Research School (NSTRI) P.O. Box
More informationPhysics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1
Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson
JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises
More informationTravel time approach to kinetically sorbing solute by diverging radial flows through heterogeneous porous formations
WATER RESOURCES RESEARCH, VOL. 4, W157, oi:1.19/1wr16, 1 Travel time approach to kinetically sorbing solute by iverging raial flows through heterogeneous porous formations Geraro Severino, 1 Samuele De
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationAnalytical accuracy of the one dimensional heat transfer in geometry with logarithmic various surfaces
Cent. Eur. J. Eng. 4(4) 014 341-351 DOI: 10.478/s13531-013-0176-8 Central European Journal of Engineering Analytical accuracy of the one imensional heat transfer in geometry with logarithmic various surfaces
More informationP. A. Martin b) Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom
Time-harmonic torsional waves in a composite cyliner with an imperfect interface J. R. Berger a) Division of Engineering, Colorao School of Mines, Golen, Colorao 80401 P. A. Martin b) Department of Mathematics,
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationState-Space Model for a Multi-Machine System
State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationCapillary effect on water table fluctuations in unconfined aquifers
WATER RESOURCES RESEARCH, VOL. 49, 3064 3069, doi:10.1002/wrcr.20237, 2013 Capillary effect on water table fluctuations in unconfined aquifers Jun Kong, 1 Cheng-Ji Shen, 2 Pei Xin, 2 Zhiyao Song, 3 Ling
More informationThe Principle of Least Action and Designing Fiber Optics
University of Southampton Department of Physics & Astronomy Year 2 Theory Labs The Principle of Least Action an Designing Fiber Optics 1 Purpose of this Moule We will be intereste in esigning fiber optic
More informationEXACT TRAVELING WAVE SOLUTIONS FOR A NEW NON-LINEAR HEAT TRANSFER EQUATION
THERMAL SCIENCE, Year 017, Vol. 1, No. 4, pp. 1833-1838 1833 EXACT TRAVELING WAVE SOLUTIONS FOR A NEW NON-LINEAR HEAT TRANSFER EQUATION by Feng GAO a,b, Xiao-Jun YANG a,b,* c, an Yu-Feng ZHANG a School
More informationNumerical Integrator. Graphics
1 Introuction CS229 Dynamics Hanout The question of the week is how owe write a ynamic simulator for particles, rigi boies, or an articulate character such as a human figure?" In their SIGGRPH course notes,
More informationTraveling wave solution of the Boussinesq equation for groundwater flow in horizontal aquifers
WATER RESOURCES RESEARCH, VOL. 49, 1668 1679, doi:1.1/wrcr.168, 13 Traveling wave solution of the Boussinesq equation for groundwater flow in horizontal aquifers H. A. Basha 1 Received 14 August 1; revised
More informationEffect of Rotation on Thermosolutal Convection. in a Rivlin-Ericksen Fluid Permeated with. Suspended Particles in Porous Medium
Av. Theor. Appl. Mech., Vol. 3,, no. 4, 77-88 Effect of Rotation on Thermosolutal Convection in a Rivlin-Ericksen Flui Permeate with Suspene Particles in Porous Meium A. K. Aggarwal Department of Mathematics
More information6. Friction and viscosity in gasses
IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner
More informationwater adding dye partial mixing homogenization time
iffusion iffusion is a process of mass transport that involves the movement of one atomic species into another. It occurs by ranom atomic jumps from one position to another an takes place in the gaseous,
More informationEfficient Macro-Micro Scale Coupled Modeling of Batteries
A00 Journal of The Electrochemical Society, 15 10 A00-A008 005 0013-651/005/1510/A00/7/$7.00 The Electrochemical Society, Inc. Efficient Macro-Micro Scale Couple Moeling of Batteries Venkat. Subramanian,*,z
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationSensors & Transducers 2015 by IFSA Publishing, S. L.
Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity
More information3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects
IOP Conference Series: Materials Science an Engineering 3-D FEM Moeling of fiber/matrix interface eboning in UD composites incluing surface effects To cite this article: A Pupurs an J Varna 2012 IOP Conf.
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 14 (2009) 3901 3913 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Dynamics of a ring network
More informationThe Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis
The Kepler Problem For the Newtonian 1/r force law, a miracle occurs all of the solutions are perioic instea of just quasiperioic. To put it another way, the two-imensional tori are further ecompose into
More informationarxiv: v1 [math-ph] 5 May 2014
DIFFERENTIAL-ALGEBRAIC SOLUTIONS OF THE HEAT EQUATION VICTOR M. BUCHSTABER, ELENA YU. NETAY arxiv:1405.0926v1 [math-ph] 5 May 2014 Abstract. In this work we introuce the notion of ifferential-algebraic
More informationThermal Modulation of Rayleigh-Benard Convection
Thermal Moulation of Rayleigh-Benar Convection B. S. Bhaauria Department of Mathematics an Statistics, Jai Narain Vyas University, Johpur, Inia-3400 Reprint requests to Dr. B. S.; E-mail: bsbhaauria@reiffmail.com
More informationStability of travelling waves
Stability of travelling waves Björn Sanstee Department of Mathematics Ohio State University 231 West 18th Avenue Columbus, OH 43210, USA E-mail: sanstee.1@osu.eu Abstract An overview of various aspects
More informationStorage-dependent drainable porosity for complex hillslopes
WATER RESOURCES RESEARCH, VOL. 41,, doi:10.1029/2004wr003725, 2005 Storage-dependent drainable porosity for complex hillslopes A. G. J. Hilberts and P. A. Troch Hydrology and Quantitative Water Management
More informationA. Incorrect! The letter t does not appear in the expression of the given integral
AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question
More informationNUMERICAL STUDY OF THERMAL RADIATIONS AND THERMAL STRATIFICATION MECHANISMS IN MHD CASSON FLUID FLOW. and Sardar Muhammad BILAL c
NUMERICAL STUDY OF THERMAL RADIATIONS AND THERMAL STRATIFICATION MECHANISMS IN MHD CASSON FLUID FLOW Khalil Ur REHMAN b c * Noor Ul SABA b Iffat ZEHRA c Muhamma Yousaf MALIK ab an Sarar Muhamma BILAL c
More informationSemiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom
PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,
More informationOn the number of isolated eigenvalues of a pair of particles in a quantum wire
On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in
More informationObjective: To introduce the equations of motion and describe the forces that act upon the Atmosphere
Objective: To introuce the equations of motion an escribe the forces that act upon the Atmosphere Reaing: Rea pp 18 6 in Chapter 1 of Houghton & Hakim Problems: Work 1.1, 1.8, an 1.9 on pp. 6 & 7 at the
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationThe Three-dimensional Schödinger Equation
The Three-imensional Schöinger Equation R. L. Herman November 7, 016 Schröinger Equation in Spherical Coorinates We seek to solve the Schröinger equation with spherical symmetry using the metho of separation
More information12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes
Differentiation of vectors 12.5 Introuction The area known as vector calculus is use to moel mathematically a vast range of engineering phenomena incluing electrostatics, electromagnetic fiels, air flow
More informationBasic Differentiation Rules and Rates of Change. The Constant Rule
460_00.q //04 4:04 PM Page 07 SECTION. Basic Differentiation Rules an Rates of Change 07 Section. The slope of a horizontal line is 0. Basic Differentiation Rules an Rates of Change Fin the erivative of
More informationMAE 210A FINAL EXAM SOLUTIONS
1 MAE 21A FINAL EXAM OLUTION PROBLEM 1: Dimensional analysis of the foling of paper (2 points) (a) We wish to simplify the relation between the fol length l f an the other variables: The imensional matrix
More informationOn Using Unstable Electrohydraulic Valves for Control
Kailash Krishnaswamy Perry Y. Li Department of Mechanical Engineering, University of Minnesota, 111 Church St. SE, Minneapolis, MN 55455 e-mail: kk,pli @me.umn.eu On Using Unstable Electrohyraulic Valves
More informationPAijpam.eu RELATIVE HEAT LOSS REDUCTION FORMULA FOR WINDOWS WITH MULTIPLE PANES Cassandra Reed 1, Jean Michelet Jean-Michel 2
International Journal of Pure an Applie Mathematics Volume 97 No. 4 2014 543-549 ISSN: 1311-8080 (printe version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu oi: http://x.oi.org/10.12732/ijpam.v97i4.13
More informationOPTIMAL CONTROL OF A PRODUCTION SYSTEM WITH INVENTORY-LEVEL-DEPENDENT DEMAND
Applie Mathematics E-Notes, 5(005), 36-43 c ISSN 1607-510 Available free at mirror sites of http://www.math.nthu.eu.tw/ amen/ OPTIMAL CONTROL OF A PRODUCTION SYSTEM WITH INVENTORY-LEVEL-DEPENDENT DEMAND
More informationA Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges
Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key
More informationThe dynamics of the simple pendulum
.,, 9 G. Voyatzis, ept. of Physics, University of hessaloniki he ynamics of the simple penulum Analytic methos of Mechanics + Computations with Mathematica Outline. he mathematical escription of the moel.
More informationA new identification method of the supply hole discharge coefficient of gas bearings
Tribology an Design 95 A new ientification metho of the supply hole ischarge coefficient of gas bearings G. Belforte, F. Colombo, T. Raparelli, A. Trivella & V. Viktorov Department of Mechanics, Politecnico
More informationSOLUTIONS for Homework #3
SOLUTIONS for Hoework #3 1. In the potential of given for there is no unboun states. Boun states have positive energies E n labele by an integer n. For each energy level E, two syetrically locate classical
More informationThe proper definition of the added mass for the water entry problem
The proper efinition of the ae mass for the water entry problem Leonaro Casetta lecasetta@ig.com.br Celso P. Pesce ceppesce@usp.br LIE&MO lui-structure Interaction an Offshore Mechanics Laboratory Mechanical
More informationStability Analysis and Analytical Solution of a Nonlinear Model for Controlled Drug Release: Travelling Wave Fronts
Stability Analysis an Analytical Solution of a Nonlinear Moel for Controlle Drug Release: Travelling Wave Fronts Chontita Rattanaul, an Yongwimon Lenbury Abstract In this paper, the process of rug issolution
More informationA Second Time Dimension, Hidden in Plain Sight
A Secon Time Dimension, Hien in Plain Sight Brett A Collins. In this paper I postulate the existence of a secon time imension, making five imensions, three space imensions an two time imensions. I will
More informationCE2253- APPLIED HYDRAULIC ENGINEERING (FOR IV SEMESTER)
CE5-APPLIED HYDRAULIC ENGINEERING/UNIT-II/UNIFORM FLOW CE5- APPLIED HYDRAULIC ENGINEERING (FOR IV SEMESTER) UNIT II- UNIFORM FLOW CE5-APPLIED HYDRAULIC ENGINEERING/UNIT-II/UNIFORM FLOW CE5- APPLIED HYDRAULIC
More informationarxiv:nlin/ v1 [nlin.cd] 21 Mar 2002
Entropy prouction of iffusion in spatially perioic eterministic systems arxiv:nlin/0203046v [nlin.cd] 2 Mar 2002 J. R. Dorfman, P. Gaspar, 2 an T. Gilbert 3 Department of Physics an Institute for Physical
More informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More information1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.
Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency
More information[Kumar*, 5(2): February, 2016] ISSN: (I2OR), Publication Impact Factor: 3.785
[Kumar*, 5(): February, 6] ISSN: 77-9655 (IOR), Publication Impact Factor: 3.785 IJESRT INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY THERMOSOLUTAL CONVECTION IN A HETEROGENEOUS VISCO-ELASTIC
More information