DOCSLIB.ORG

Mark Boady, Pavel Grinfeld, Jeremy Johnson

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Mark Boady, Pavel Grinfeld, Jeremy Johnson Approximation of Fourier Series to Determine Laplace-Dirichlet Eigenvalues Mark Boady, Pavel Grinfeld, Jeremy Johnson MOTIVATION EVALUATION RESULTS In 2004, Grinfeld and Strang [3] posed the following boundary problem. The CMS provides expressions for the derivatives of λ but evaluation is Evaluating more terms increases accuracy but also the time and number of What is the series in 1/N for the simple Laplace eigenvalues λN on a regular dependent on the surface. terms that need to be simplified. polygon with N sides under the Dirichlet boundary condition? The first four Z terms in the series were computed by hand in 2012 [4]. Our goal is to automate i λ1 = − Criur udS: (7) the process and derive more terms in the series. S Z The Calculus of Moving Surfaces (CMS), an extension of tensor calculus 2 α i δC i i @u λ2 = (C Bα r uriu − r uriu − 2Cr riu to deforming manifold, provides an analytic method to describe the change S δτ @t caused by the deformation of a surface. The Laplace eigenvalues on the unit 2 m i − 2C N riur rmu)dS (8) circle under the Dirichlet boundary condition are known. We have applied the CMS to analyze the change in the eigenvalues between the unit circle and N- The surface deformation is encapsulated in the surface velocity, C, and the sided regular polygon[2]. It remains to evaluate these expressions and form a @u partial derivatives of u, @t . The terms only need to be evaluated at t = 0. Taylor series for the eigenvalues on the N-sided polygon. This problem is interesting in applied mathematics[5], physics[1] and pat- − cos(θ) + cos(π=N) Cjt=0 = (9) We believe that all values will converge to integer multiple of the ζ function tern recognition[6]. It also presents a number of analytical challenges associ- cos(θ) iζ(j)ρk in the form N j for integers i,j, and k. ated with the singularities present in the problem. This problem has important To work with this expression we represent it as a Fourier series, implications for the analysis of numerical methods in which a smooth bound- P1 ikNθ c0(k)e where ary is replaced by an approximated shape. k=1 ( π2 − 3N 2 k = 0 c0(k) = (−1)k (−1)kπ2 5(−1)k (10) SURFACE DEFORMATION N 2k2 + 3N 4k4 − N 4k4 + ··· otherwise The Laplace-Dirichlet Eigenvalue problem is described by the following Next, we evaluate λ1 in terms of the Fourier coefficients. set of equations for a boundary S embedded in space Ω. 2 2 2 2π ρ 4u = −λu (1) λ1 = −2c0(0)ρ = (11) 3N 2 ujS = 0 (2) For λ2, the result includes convolutions of coefficients. 1 Z Our current Taylor series up to N 5 is u2dΩ = 1 (3) Ω p 2 2 δu1 π2ρ2 ζ (3) ρ2 14 π4ρ2 λ2 = 2conv(c0; c0)(0)ρ − 2c1(0)ρ + 4 πconv c0; jr=1 (0) (12) λ(t = 1) = ρ2 + 2=3 + 4 + The unknowns, u and λ, can be found when S is the unit circle and Ω is de- δr N 2 N 3 45 N 4 scribed in polar coordinates. 2 2 @u 1 ρ −45360 ζ (5) ρ + 453600 ζ (5) Where u1(k), is the k-th coefficient of the Fourier series for @t . The ex- + (17) J (ρr) 22680 N 5 u(r; θ) = p 0 (4) pressions for ui are derived from the CMS and will also include convolutions πJ1(ρ) of the ci coefficients. 2 CHALLENGES λ = ρ (5) ( c (0)(ρrJ (ρr)−J (ρr)) 0 p1 0 k = 0 πJ1(ρ) Although the results are simple, the expression swells rapidly during evalu- u1(k) = c (k)ρJ (ρr) (13) In this expression, Jm is the m-th Bessel Function and ρ is the n-th root of J0. 0 kNp ation. The highest we can currently determine λ4 is with convolutions of range J (ρ) π otherwise The surface will start at t = 0 with the unit circle and deform into an N-sided kN m = −32 ··· 32. For this sum, we need to evaluate 22,494 terms which took regular polygon at time t = 1. The area will be constant for the deformation. By approximating these convolutions, we can determine the solution to λ2. over an hour. Our goal is to optimize the calculations using methods such as 2 4 2 parallel processing to solve λ4 for range m = −256 ··· 256 in minutes. 8ζ(3)ρ 32π ρ −5 λ2 = + + O(N ) (14) N 3 45N 4 REFERENCES [1] P. Amore. Spectroscopy of drums and quantum billiards: perturbative and non-perturbative CONVOLUTIONS results. J. Math. Phys., 51:052105, 2010. Convolutions are approximated using a custom Maple library. [2] M. Boady, P. Grinfeld, and J. Johnson. A Term Rewriting System for the Calculus of Moving Surfaces. International Symposium on Symbolic and Algebraic Computation (ISSAC), 2013. m X [3] P. Grinfeld and G. Strang. Laplace eigenvalues on polygons. Computers and Mathematics conv(a; b)(i) = a(k)b(i − k) (15) with Applications 48:1121–1133, 2004. k=−m [4] P. Grinfeld and G. Strang. Laplace eigenvalues on regular polygons: A series in 1/N. Journal We will use a Taylor series to describe λ(t = 1) in terms of the i-th deriva- of Mathematical Analysis and Applications 385(1):135–149, 2012. Deeper convolutions require more terms approximate. [5] P. Guidotti and J. V. Lambers. Eigenvalue characterization and computation for the laplacian tives of λi(t = 0) = λi. on general 2-d domains. Numerical Functional Analysis and Optimization. 29:507-531, m m ! 1 1 1 X X 2008. λ(t = 1) = λ + λ + λ + λ + λ + ··· (6) conv(c; conv(a; b))(i) = c(j) a(k)b((i − j) − k) (16) 1 2! 2 3! 3 4! 4 [6] M. Khabou, L. Hermi, and M. Rhouma. Shape recognition using eigenvalues of the dirichlet j=−m k=−m laplacian. Pattern Recognition, 40(1):141-153, 2007..
Load more

Recommended publications
  • ED16 Abstracts 73
    ED16 Abstracts 73 IP1 industrial problems. She will give examples of some indus- Mathematical Modeling with Elementary School- trial research problems that student teams have tackled, Aged Students discuss some of the skills that are needed to be successful working with and in industry, discuss the challenges that Modeling, a cyclic process by which mathematicians de- one may face when working with corporate partners, and velop and use mathematical tools to represent, understand, present a summary of some of the lessons that have been and solve real-world problems, provides important learning learned over the years in these collaborations. opportunities for school students. Modeling opportunities in secondary schools are apparent, but what about in the Suzanne L. Weekes younger grades? Two questions are critical in mathemati- Worcester Polytechnic Institute (WPI) cal modeling in K-5 settings. (1) How should opportunities [email protected] for modeling in K-5 settings be constructed and carried out? (2) What are the tasks of teaching when engaging el- ementary students in mathematical modeling? In this talk IP5 I will present a framework for teaching mathematical mod- Title Not Available at Time of Publication eling in elementary classrooms and provide illustrations of its use by elementary grades teachers. Abstract not available at time of publication. Elizabeth A. Burroughs Philip Uri Treisman Montana State University The University of Texas at Austin [email protected] [email protected] IP2 CP1 Graduate Student Education in Computational Regime Switching Models and the Mental Account- Mathematics and Scientific Computing ing Framework Abstract not available at time of publication.
    [Show full text]
  • 2006 - 2007 Annual Report Department of Mathematics College of Arts & Sciences Drexel University Photo by R
    2006 - 2007 Annual Report Department of Mathematics College of Arts & Sciences Drexel University Photo by R. Andrew Hicks Andrew R. by Photo TABLE OF CONTENTS PEOPLE 1 Message From the Department Head 2 Tenure Track Faculty 4 Auxiliary Faculty 5 Visiting Faculty & Fellows 5 Emeritus Faculty 6 Staff 6 Teaching Assistants & Research Assistants 7 New Faculty Profiles AWARDS / HIGHLIGHTS 8 Faculty Awards 9 Faculty Grants 10 Faculty Appointments / Conference Organizations 11 Faculty Publications / Patents 13 Faculty Presentations 17 Faculty New Course Offerings 18 Honors Day Awards 20 Graduate Student Day Awards 21 Albert Herr Teaching Assistant Award 21 Student Presentations 22 Bachelor of Science Degrees Awarded 22 Master of Science Degrees Awarded 22 Doctor of Philosophy Degree Awarded DEPARTMENT 23 Colloquium Photo by R. Andrew Hicks 25 Analysis Seminars 27 Dean’s Seminar 28 Departmental Committees 30 Mathematics Resource Center 31 Departmental Outreach 33 Financial Mathematics Advisory Board EVENTS 35 Student Activities 36 SIAM Chapter / Math Bytes 37 Departmental Social Events GIVING 38 Support the Mathematics Department MESSAGE FROM THE DEPARTMENT HEAD Dear Alumni and Friends, It has been another great year for our department. It is wonderful to have award winning faculty in the department. Dr. Pavel Grinfeld received the College of Arts and Sciences Excellence Award for both his research and teaching accomplishments, and Ms. Marna Mozeff Hartmann received the Barbara G. Hornum Award for Excellence in Teaching for her dedication to teaching and advising our undergraduates. Also several of our students have been recognized. Teaching Assistants Emek Kose and Jason Aran have won two of Drexel’s Teaching Assistant Excellence Awards.
    [Show full text]
  • Philosophical Magazine A
    This article was downloaded by:[Massachusetts Institute of Technology] On: 7 January 2008 Access Details: [subscription number 768485848] Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer House, 37-41 Mortimer Street, London W1T 3JH, UK Philosophical Magazine A Publication details, including instructions for authors and subscription information: http://www.informaworld.com/smpp/title~content=t713396797 Towards thermodynamics of elastic electric conductors M. Grinfeld a; P. Grinfeld b a The Educational Testing Service, Princeton, New Jersey, USA b Department of Mathematics, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA Online Publication Date: 01 May 2001 To cite this Article: Grinfeld, M. and Grinfeld, P. (2001) 'Towards thermodynamics of elastic electric conductors', Philosophical Magazine A, 81:5, 1341 - 1354 To link to this article: DOI: 10.1080/01418610108214445 URL: http://dx.doi.org/10.1080/01418610108214445 PLEASE SCROLL DOWN FOR ARTICLE Full terms and conditions of use: http://www.informaworld.com/terms-and-conditions-of-access.pdf This article maybe used for research, teaching and private study purposes. Any substantial or systematic reproduction, re-distribution, re-selling, loan or sub-licensing, systematic supply or distribution in any form to anyone is expressly forbidden. The publisher does not give any warranty express or implied or make any representation that the contents will be complete or accurate or up to date. The accuracy of any instructions, formulae and drug doses should be independently verified with primary sources. The publisher shall not be liable for any loss, actions, claims, proceedings, demand or costs or damages whatsoever or howsoever caused arising directly or indirectly in connection with or arising out of the use of this material.
    [Show full text]
  • Applications of Symbolic Computation to the Calculus of Moving Surfaces a Thesis Submitted to the Faculty of Drexel University B
    Applications of Symbolic Computation to the Calculus of Moving Surfaces A Thesis Submitted to the Faculty of Drexel University by Mark W. Boady in partial fulfillment of the requirements for the degree of Doctor of Philosophy June 2016 c Copyright 2016 Mark W. Boady. All Rights Reserved. i Dedications To my lovely wife, Petrina. Without your endless kindness, love, and patience, this thesis never would have been completed. To my parents, Bill and Charlene, you have loved and supported me in every endeavor. I am eternally grateful to be your son. To my in-laws, Jerry and Annelie, your endless optimism and assistance helped keep me motivated throughout this project. ii Acknowledgments I would like to express my sincere gratitude to my supervisors Dr. Jeremy Johnson and Dr. Pavel Grinfeld. Dr. Johnson's patience and insight provided a safe harbor during the seemingly endless days of debugging and testing code. I cannot thank him enough for improvements in my writing, coding, and presenting skills. This thesis could not have existed without Dr. Grinfeld's pioneering research. His vast knowledge and endless enthusiasm propelled us forward at every roadblock. A special thanks also goes out to the members of my thesis committee: Dr. Bruce Char, Dr. Dario Salvucci, and Dr. Hoon Hong. Without their valuable input, this work would not have been fully realized. They consistently helped to keep this project focused while I was buried in the details of development. All the members of the Drexel Computer Science Department have provided endless help in my development as a teacher and researcher.
    [Show full text]
  • Lecture Notes: Mathematical Methods
    LECTURE NOTES IMPERIAL COLLEGE LONDON DEPARTMENT OF COMPUTING Mathematical Methods Instructors: Marc Deisenroth, Mahdi Cheraghchi Version: November 30, 2016 Contents 1 Course Overview3 2 Sequences6 2.1 The Convergence Definition.......................6 2.2 Illustration of Convergence.......................7 2.3 Common Converging Sequences.....................8 2.4 Combinations of Sequences.......................9 2.5 Sandwich Theorem............................9 2.6 Ratio Tests for Sequences........................ 11 2.7 Proof of Ratio Tests............................ 12 2.8 Useful Techniques for Manipulating Absolute Values......... 13 2.9 Properties of Real Numbers....................... 14 3 Series 15 3.1 Geometric Series............................. 16 3.2 Harmonic Series............................. 16 3.3 Series of Inverse Squares......................... 17 3.4 Common Series and Convergence.................... 18 3.5 Convergence Tests............................ 19 3.6 Absolute Convergence.......................... 22 3.7 Power Series and the Radius of Convergence.............. 23 3.8 Proofs of Ratio Tests........................... 24 4 Power Series 27 4.1 Basics of Power Series.......................... 27 4.2 Maclaurin Series............................. 27 4.3 Taylor Series............................... 31 4.4 Taylor Series Error Term......................... 33 4.5 Deriving the Cauchy Error Term..................... 34 4.6 Power Series Solution of ODEs..................... 34 5 Linear Algebra 36 5.1 Linear Equation Systems........................
    [Show full text]
  • GILBERT STRANG Education Positions Held Awards and Duties
    CURRICULUM VITAE : GILBERT STRANG Work: Department of Mathematics, M.I.T. Cambridge MA 02139 Phone: 617{253{4383 Fax: 617{253{4358 e-mail: [email protected] URL: http://math.mit.edu/ gs Home: 7 Southgate Road Wellesley MA 02482 Phone: 781{235{9537 Education 1. 1952{1955 William Barton Rogers Scholar, M.I.T. (S.B. 1955) 2. 1955{1957 Rhodes Scholar, Oxford University (B.A., M.A. 1957) 3. 1957{1959 NSF Fellow, UCLA (Ph.D. 1959) Positions Held 1. 1959{1961 C.L.E. Moore Instructor, M.I.T. 2. 1961{1962 NATO Postdoctoral Fellow, Oxford University 3. 1962{1964 Assistant Professor of Mathematics, M.I.T. 4. 1964{1970 Associate Professor of Mathematics, M.I.T. 5. 1970{ Professor of Mathematics, M.I.T. Awards and Duties 1. Alfred P. Sloan Fellow (1966{1967) 2. Chairman, M.I.T. Committee on Pure Mathematics (1975{1979) 3. Chauvenet Prize, Mathematical Association of America (1976) 4. Council, Society for Industrial and Applied Mathematics (1977{1982) 5. NSF Advisory Panel on Mathematics (1977{1980) (Chairman 1979{1980) 6. CUPM Subcommittee on Calculus (1979{1981) 7. Fairchild Scholar, California Institute of Technology (1980{1981) 8. Honorary Professor, Xian Jiaotong University, China (1980) 9. American Academy of Arts and Sciences (1985) 10. Taft Lecturer, University of Cincinnati (1977) 11. Gergen Lecturer, Duke University (1983) 12. Lonseth Lecturer, Oregon State University (1987) 13. Magnus Lecturer, Colorado State University (2000) 14. Blumberg Lecturer, University of Texas (2001) 15. AMS-SIAM Committee on Applied Mathematics (1990{1992) 16. Vice President for Education, SIAM (1991{1996) 17.
    [Show full text]
  • Electroconducting Sphere Inside Unbounded Isotropic Matrix for ALEGRA Verification by Michael Grinfeld, Pavel Grinfeld, John Niederhaus, and Angel Rodriguez
    ARL-TR-8994 ● JULY 2020 Electroconducting Sphere inside Unbounded Isotropic Matrix for ALEGRA Verification by Michael Grinfeld, Pavel Grinfeld, John Niederhaus, and Angel Rodriguez Approved for public release; distribution is unlimited. NOTICES Disclaimers The findings in this report are not to be construed as an official Department of the Army position unless so designated by other authorized documents. Citation of manufacturer’s or trade names does not constitute an official endorsement or approval of the use thereof. Destroy this report when it is no longer needed. Do not return it to the originator. ARL-TR-8994 ● JULY 2020 Electroconducting Sphere inside Unbounded Isotropic Matrix for ALEGRA Verification Michael Grinfeld Weapons and Materials Research Directorate, CCDC Army Research Laboratory Pavel Grinfeld Drexel University John Niederhaus and Angel Rodriguez Sandia National Laboratories Approved for public release; distribution is unlimited. Form Approved REPORT DOCUMENTATION PAGE OMB No. 0704-0188 Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing the burden, to Department of Defense, Washington Headquarters Services, Directorate for Information Operations and Reports (0704-0188), 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to any penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number.
    [Show full text]
  • Closed, Two Dimensional Surface Dynamics
    Closed, Two Dimensional Surface Dynamics David V. Svintradze∗ Department of Physics, Tbilisi State University, Chavchavadze Ave. 03, 0179 Tbilisi, Georgia (Dated: February 21, 2018) We present dynamic equations for two dimensional closed surfaces and analytically solve it for some simplified cases. We derive final equations for surface normal motions by two different ways. The solution of the equations of motions in normal direction indicates that any closed, two di- mensional, homogeneous surface with time invariable surface energy density adopts constant mean curvature shape when it comes in equilibrium with environment. As an example, we apply the formalism to analyze equilibrium shapes of micelles and explain why they adopt spherical, lamellar and cylindrical shapes. We show that theoretical calculation for micellar optimal radius is in good agreement with all atom simulations and experiments. I. INTRODUCTION balance equations [11, 12]. Specific dynamical equations accounting for bending as well as electrodynamic effects Biological systems exhibit a variety of morphologies have also been reported [13{15]. Furthermore, active and experience large shape deformations during a mo- membrane theories have extended our understanding of tion. Such 'choreography' of shape motility is charac- passive membranes. Active membrane theories include teristic not only for all living organisms and cells [1] external forces [16{19] and provide a framework for the but also for proteins, nucleic acids and to all biomacro- study of active biological or chemical processes at sur- molecules in general. Shape motility, which is a motion faces, such as the cell cortex, the mechanics of epithelial of two-dimensional surfaces, may be a result of active tissues, or reconstituted active systems on surfaces [16].
    [Show full text]
  • Mathematical Methods
    LECTURE NOTES IMPERIAL COLLEGE LONDON DEPARTMENT OF COMPUTING Mathematical Methods Marc Deisenroth, Mahdi Cheraghchi Version: December 1, 2017 Contents 1 Course Overview3 2 Sequences6 2.1 The Convergence Definition.......................6 2.2 Illustration of Convergence.......................7 2.3 Common Converging Sequences.....................8 2.4 Combinations of Sequences.......................9 2.5 Sandwich Theorem............................9 2.6 Ratio Tests for Sequences........................ 11 2.7 Proof of Ratio Tests............................ 12 2.8 Useful Techniques for Manipulating Absolute Values......... 13 2.9 Properties of Real Numbers....................... 14 3 Series 15 3.1 Geometric Series............................. 16 3.2 Harmonic Series............................. 16 3.3 Series of Inverse Squares......................... 17 3.4 Common Series and Convergence.................... 18 3.5 Convergence Tests............................ 19 3.6 Absolute Convergence.......................... 22 3.7 Power Series and the Radius of Convergence.............. 23 3.8 Proofs of Ratio Tests........................... 24 4 Power Series 27 4.1 Basics of Power Series.......................... 27 4.2 Maclaurin Series............................. 27 4.3 Taylor Series............................... 31 4.4 Taylor Series Error Term......................... 33 4.5 Deriving the Cauchy Error Term..................... 34 4.6 Power Series Solution of ODEs..................... 34 5 Linear Algebra 36 5.1 Linear Equation Systems........................
    [Show full text]
  • PD09 Abstracts 45
    PD09 Abstracts 45 IP1 functions. Free Boundary Regularity Yann Brenier We discuss level surfaces of solutions u to singular semilin- CNRS ear elliptic equations such as Δu = δ(u). The level surface Universit´ede Nice Sophia-Antipolis u = 0 can be interpreted as a solution to a free boundary [email protected] problem, the problem of finding the optimal shape of in- sulating material. A variant of this equation, the Prandtl- Batchelor equation, describes the wake created by a boat. IP4 We will prove that stable free boundaries are smooth in The Importance of PDEs for Modeling and Simu- three-dimensional space. Numerical methods are needed lation on Peta and Exascale Systems to explore counterexamples in higher dimensions and to establish a more robust regularity theory with numerically As scientists and engineers tackle more complex problems effective bounds. involving multiphysics and multiscales and then seek to solve these problems through computation, the partial dif- David Jerison ferential equation formulations which incorporate and more Massachusetts Institute of Technology closely approximate the complex physical or biological phe- [email protected] nomena being modeled will also be important to unlocking the approach taken on Peta and Exascale systems. Scien- tists and engineers working on problems in such areas as IP2 aerospace, biology, climate modeling, energy and other ar- Fluid Biomembranes: Modeling and Computation eas demand ever increasing compute power for their prob- lems. For Petascale and Exascale systems to be useful mas- We study two models for biomembranes. The first one is sive parallelism at the chip level is not sufficient.
    [Show full text]
  • Introduction to Tensor Analysis and the Calculus of Moving Surfaces
    Introduction to Tensor Analysis and the Calculus of Moving Surfaces Pavel Grinfeld Introduction to Tensor Analysis and the Calculus of Moving Surfaces 123 Pavel Grinfeld Department of Mathematics Drexel University Philadelphia, PA, USA ISBN 978-1-4614-7866-9 ISBN 978-1-4614-7867-6 (eBook) DOI 10.1007/978-1-4614-7867-6 Springer New York Heidelberg Dordrecht London Library of Congress Control Number: 2013947474 Mathematics Subject Classifications (2010): 4901, 11C20, 15A69, 35R37, 58A05, 51N20, 51M05, 53A05, 53A04 © Springer Science+Business Media New York 2013 This work is subject to copyright. All rights are reserved by the Publisher, whether the whole or part of the material is concerned, specifically the rights of translation, reprinting, reuse of illustrations, recitation, broadcasting, reproduction on microfilms or in any other physical way, and transmission or information storage and retrieval, electronic adaptation, computer software, or by similar or dissimilar methodology now known or hereafter developed. Exempted from this legal reservation are brief excerpts in connection with reviews or scholarly analysis or material supplied specifically for the purpose of being entered and executed on a computer system, for exclusive use by the purchaser of the work. Duplication of this publication or parts thereof is permitted only under the provisions of the Copyright Law of the Publisher’s location, in its current version, and permission for use must always be obtained from Springer. Permissions for use may be obtained through RightsLink at the Copyright Clearance Center. Violations are liable to prosecution under the respective Copyright Law. The use of general descriptive names, registered names, trademarks, service marks, etc.
    [Show full text]
  • Notes on the Nonlinear Dependence of a Multiscale Coupled System with Respect to the Interface
    NOTES ON THE NONLINEAR DEPENDENCE OF A MULTISCALE COUPLED SYSTEM WITH RESPECT TO THE INTERFACE Fernando A Morales Escuela de Matem´aticas,Universidad Nacional de Colombia, Sede Medell´ın. Calle 59 A No 63-20, Of 43-106. Medell´ın,Colombia. Abstract. This work studies the dependence of the solution with respect to interface geometric perturbations in a multiscaled coupled Darcy flow system in direct variational formulation. A set of admissible perturbation functions and a sense of convergence are presented, as well as sufficient conditions on the forcing terms, in order to conclude strong convergence statements. For the rate of convergence of the solutions we start solving completely the one dimensional case using orthogonal decompositions on appropriate subspaces. Finally, the rate of convergence question is analyzed in a simple multiple dimen- sional setting, studying the nonlinear operators introduced by the geometric perturbations. 1. Introduction The study of saturated flow in geological porous media frequently presents nat- ural structures with a dense network of fissures nested in the rock matrix [21]. It is also frequent to observe vuggy porous media [1], which have the presence of cavities in the rock matrix significantly larger than the average pore size of the medium. This is a multiple scale physics phenomenon, because there are regions of the medium where the flow velocity is significantly larger than the velocity on the other ones. The modeling of the interface between regions is subject to very active research: first, fluid transmission conditions across the interface are of great importance, see [19,3] for discussion of the governing laws; see [7,2] for a numerical point of view; see [1, 12,5, 16] for the analytic approach and [20,4] for a more general perspec- tive.
    [Show full text]