The finite element method explicit scheme for a solution of one problem of surface and ground water combined movement

Similar documents
ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

The Finite Element Method: A Short Introduction

On a nonlinear compactness lemma in L p (0, T ; B).

Numerical Simulation of One-Dimensional Wave Equation by Non-Polynomial Quintic Spline

Adaptive Kernel Estimation of the Conditional Quantiles

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Convexity preserving interpolation by splines of arbitrary degree

Solution for singularly perturbed problems via cubic spline in tension

Sharp integral inequalities involving high-order partial derivatives. Journal Of Inequalities And Applications, 2008, v. 2008, article no.

TR/95 February Splines G. H. BEHFOROOZ* & N. PAPAMICHAEL

Stanford University CS254: Computational Complexity Notes 7 Luca Trevisan January 29, Notes for Lecture 7

Shuai Dong. Isaac Newton. Gottfried Leibniz

5 The Laplace Equation in a convex polygon

Appendix B. The Finite Difference Scheme

TR/28. OCTOBER CUBIC SPLINE INTERPOLATION OF HARMONIC FUNCTIONS BY N. PAPAMICHAEL and J.R. WHITEMAN.

Numerical Heat and Mass Transfer

An identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites

More metrics on cartesian products

General viscosity iterative method for a sequence of quasi-nonexpansive mappings

Multivariate Ratio Estimator of the Population Total under Stratified Random Sampling

APPENDIX A Some Linear Algebra

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

Ballot Paths Avoiding Depth Zero Patterns

A boundary element method with analytical integration for deformation of inhomogeneous elastic materials

The Order Relation and Trace Inequalities for. Hermitian Operators

Homogenization of reaction-diffusion processes in a two-component porous medium with a non-linear flux-condition on the interface

Pathan Mahabub Basha *, Vembu Shanthi

Perfect Competition and the Nash Bargaining Solution

arxiv: v1 [math.co] 12 Sep 2014

Existence results for a fourth order multipoint boundary value problem at resonance

Y. Guo. A. Liu, T. Liu, Q. Ma UDC

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Solving Singularly Perturbed Differential Difference Equations via Fitted Method

Modelli Clamfim Equazioni differenziali 7 ottobre 2013

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity

A Discrete Approach to Continuous Second-Order Boundary Value Problems via Monotone Iterative Techniques

2 Finite difference basics

Continuous Time Markov Chain

Finite Element Modelling of truss/cable structures

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

arxiv: v1 [math.co] 1 Mar 2014

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices

CENTROID (AĞIRLIK MERKEZİ )

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

On Pfaff s solution of the Pfaff problem

ACTM State Calculus Competition Saturday April 30, 2011

A Spline based computational simulations for solving selfadjoint singularly perturbed two-point boundary value problems

Quantum Particle Motion in Physical Space

ON THE JACOBIAN CONJECTURE

A Linearized Finite Difference Scheme with Non-uniform Meshes for Semi-linear Parabolic Equation

CENTROID (AĞIRLIK MERKEZİ )

Numerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

Exercise Solutions to Real Analysis

A new Approach for Solving Linear Ordinary Differential Equations

Report on Image warping

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Binomial transforms of the modified k-fibonacci-like sequence

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION

Research Article A Generalized Sum-Difference Inequality and Applications to Partial Difference Equations

Fixed point method and its improvement for the system of Volterra-Fredholm integral equations of the second kind

TR/01/89 February An O(h 6 ) cubic spline interpolating procedure for harmonic functions. N. Papamichael and Maria Joana Soares*

2nd International Conference on Electronics, Network and Computer Engineering (ICENCE 2016)

UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 2017/2018 FINITE ELEMENT AND DIFFERENCE SOLUTIONS

PHYS 705: Classical Mechanics. Canonical Transformation II

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

The Minimum Universal Cost Flow in an Infeasible Flow Network

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

NP-Completeness : Proofs

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

Multigrid Methods and Applications in CFD

Appendix B. Criterion of Riemann-Stieltjes Integrability

Curvature and isoperimetric inequality

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

Two point fuzzy boundary value problem with eigenvalue parameter contained in the boundary condition

Kernel Methods and SVMs Extension

High resolution entropy stable scheme for shallow water equations

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

arxiv: v1 [math.ho] 18 May 2008

Research Article Numerov s Method for a Class of Nonlinear Multipoint Boundary Value Problems

3 Basic boundary value problems for analytic function in the upper half plane

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Positivity Preserving Interpolation by Using

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper

Inductance Calculation for Conductors of Arbitrary Shape

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

SL n (F ) Equals its Own Derived Group

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

Maximizing the number of nonnegative subsets

International Conference on Advanced Computer Science and Electronics Information (ICACSEI 2013) equation. E. M. E. Zayed and S. A.

Formal solvers of the RT equation

Lecture 21: Numerical methods for pricing American type derivatives

Case Study of Markov Chains Ray-Knight Compactification

Math 702 Midterm Exam Solutions

A Pursuit Problem Described by Infinite System of Differential Equations with Coordinate-Wise Integral Constraints on Control Functions

Transcription:

IOP Conference Seres: Materals Scence and Engneerng PAPER OPEN ACCESS e fnte element metod explct sceme for a soluton of one problem of surface and ground water combned movement o cte ts artcle: L L Glazyrna and M F Pavlova 6 IOP Conf. Ser.: Mater. Sc. Eng. 58 4 Related content - On surface electromagnetc waves V N Datsko and A A Kopylov - An Explct Sceme for te KdV Equaton Wang Hu-Png Wang Yu-Sun and Hu Yng-Yng - A new parallelzaton algortm of ocean model wt explct sceme X D Fu Vew te artcle onlne for updates and enancements. s content was downloaded from IP address 37.44.93. on 7//7 at 3:5

t Internatonal Conference on "Mes metods for boundary-value problems and applcatons" IOP Publsng IOP Conf. Seres: Materals Scence and Engneerng 58 (6) 4 do:.88/757-899x/58//4 e fnte element metod explct sceme for a soluton of one problem of surface and ground water combned movement L L Glazyrna and M F Pavlova Kazan Federal Unversty 8 Kremlyovskaya str. Kazan 48 Russa E-mal: glazyrna-ludmla@ya.ru Abstract. We nvestgate te ntal-boundary value problem for te system of two nonlnear degenerate parabolc combned equatons wen te one of equatons s set n te doman we consder and te second s set on te conducted secton of ts doman. Wt elp of te semdscretzaton metods wt respect to tme varable and fnte element metod (FEM) wt respect to space varable we construct an explct dfference sceme. We obtaned condtons wc provde us convergence of constructed metod. e result of convergence was obtaned wt te mnmum condtons for smootness of te ntal data.. Statement of te problem s paper s devoted to te convergence nvestgaton of an approxmate soluton metod for te followng ntal-boundary problem: ( u) ( x t) ( k ( x u( x t) u( x t) )) f ( x t) x \ t ( ) t x () ( u) u ( x t ) k x u ( x t ) ( x t ) t s s ( ( )) [ k( x u( x t) u( x t) ) cos( n x) ] f( x t) () [ u] x t ( ) u( x) u ( x) x; u( x t) x t ( ). (3) Here s a bounded doman n space R s ts boundary s a conducted secton of and dvdng t nto two connected domans [ ] s a dscontnuty of a functon wen ts passng troug te secton n s a normal to. Problem () (3) descrbes (see []) te process of surface-water and groundwater jont moton n addton te secton s a bed rver or cannel te requred functon u defnes te egt of te water face above te mpenetrable base. We assume tat functons ( ) are absolutely contnuous strongly ncreasng and satsfy te followng condtons for arbtrary R Content from ts work may be used under te terms of te Creatve Commons Attrbuton 3. lcence. Any furter dstrbuton of ts work must mantan attrbuton to te autor(s) and te ttle of te work journal ctaton and DOI. Publsed under lcence by IOP Publsng Ltd

t Internatonal Conference on "Mes metods for boundary-value problems and applcatons" IOP Publsng IOP Conf. Seres: Materals Scence and Engneerng 58 (6) 4 do:.88/757-899x/58//4 ere ' 3 b b ( ) ( t) t dt b b (4) ( ) b b (5) 4 5 ( '( ) )' (6) b b b b b b. 3 4 5 We assume functons k ( x )( ) k( x ) are contnuous wt respect to x are measurable wt respect to x and satsfy te followng condtons for arbtrary x x R R p j j k ( x ) M M M M p (7) p k x M M3 M M3 ( ) (8) ( k ( x ) k ( x ))( ) (9) k ( x ) M M M M p () p k ( x ) M M M M () p 3 3 ( k ( x ) k ( x ))( ). () rougout te followng to defne te generalzed soluton for te problem () (3) we use normed spaces V V ( ) W( ) obtaned by closng C( ) and C ( ; C ( )) n te followng norms u u u V W ( ) W ( ) p p u u u V Lp W p Lp W p ( ) ( ; ( )) ( ; ( )) u u u u W ( ) V ( ) L( ; L ( )) L ( ; L ( )) and we denote te functonal J suc tat ts value on element v can be defned as ere z V ( ) d J( z( t)) v ( z( t)) v( x) dx ( z( t)) v( s) ds Fv s a value of functonal F In papers [] [3] we proved tat space V dt V ( ) (3) ( V ) on te element. concdes wt te set of functons from W p ( ) and ter ncrement on fromw p ( ). Also n papers [] [3] te smlar result was obtaned for V( ).

t Internatonal Conference on "Mes metods for boundary-value problems and applcatons" IOP Publsng IOP Conf. Seres: Materals Scence and Engneerng 58 (6) 4 do:.88/757-899x/58//4 Defnton. Functon u W ( ) suc tat u( x) u ( x) a.e. n and J( u) dt( V ( )) (4) s called generalzed soluton of te problem () (3) f te ntegral dentty v J ( u) v dt k ( x u u) dxdt x u v k x u dsdt f v dt f v dt ( ) s s (5) olds for arbtrary v W ( ) ; ere f v ( f v ) s te value of functonal f from te space L ( ; W ( )) ( and from L p p ( ; W ( )) respectvely) on element v from W( ). p p s paper s a contnuaton of te nvestgatons wc ave been begun n papers [4] [5] were was gven te generalzed statement of te consderng problem te teorems of exstence and unqueness of generalzed soluton ave been proved.. Descrpton of te approxmate metod. e convergence nvestgaton. For problem () (3) we construct te approxmate soluton metod wt respect to varable t and FEM wt respect to space varables. o ts end on [] we defne a unform grd { N } were s te grd ncrement. Under te assumpton tat s a convex doman we defne te polygon nscrbed n wt boundary. s polyedron satsfes te followng condtons: ) for any pont x tere exsts a pont lyng from t at a dstance not more tan ( s te grd ncrement wt respect to space varables); ) ntersecton ponts of and are te vertexes of polygon. e trangulaton of doman s carred out te followng way. For eac of te subdomans and on wc te secton dvdes te doman te trangulaton s performed autonomously suc tat te sets of mes nodes constructed on and and lyng on concde. Suc fragmentaton n negborood of secton provdes us trangles wt one curvlnear sde. Let te V l be a set of te functons from V suc tat ter restrcton on eac fnte element of te trangulaton s an mage of a polynomal of degree l on te basc element (see [6] p.3 [7] p. 47). Defnton. A functon y() t l s called semdscrete problem soluton f for arbtrary z l and t{ } te followng condtons are satsfed V V z y t z x k x y y x dx ( t( ( )) ( ) ( ) ( )) x ( ( ( )) ( ) ( y ) z ( )) t y t z s k x y s ds (6) s s 3

t Internatonal Conference on "Mes metods for boundary-value problems and applcatons" IOP Publsng IOP Conf. Seres: Materals Scence and Engneerng 58 (6) 4 do:.88/757-899x/58//4 Here t f ( x t) f( x ) d. t f ( x t) z( x) dx f ( s t) z( s) ds e followng statements were proved. Lemma. Let te condtons (4) () be satsfed. In addton ncrements were y( x) u ( x) a.e.n and. (7) f L ( ; W ( )) f L ( ; W ( )) u V L ( ) p p p p and are pcked te way tat te followng nequalty olds { } c (8) max p ( p ) p max { p p } p ( p ) max { p } p p p max p p p. { } en te followng pror estmates old for te soluton of problem (6) (7) t y( t) y( t) c (9) L ( ) ( ) L [ ] t p p t W p( ) W p( ) y( t) y( t) c () t t t y ( t) c () t t L ( ) y ( t) c t () t L ( ) k ( ( y( t k )) ( y( t)) )( y( t k ) y( t) ) dx k c (3) k N. Lemma. Let y be a sequence of functons for wc a pror estmates (9) (3) are vald. en tere exsts te functon u W ( ) followng lmt relatons as : and ncrement sequences and tat satsfy te y( t) u n V ( ) (4) 4

t Internatonal Conference on "Mes metods for boundary-value problems and applcatons" IOP Publsng IOP Conf. Seres: Materals Scence and Engneerng 58 (6) 4 do:.88/757-899x/58//4 y( t) u -weakly n L ( ; L ( )) (5) y( t) u -weakly n L ( ; L ( )) (6) were y( t) u a.e. n Q and (7) y( x t) ( x t) y( x t) ( xt ) ( \ ) z z are te pecewse constant extensons of te mes functon z. ey are equal to value zt ( ) on te sets [ t t ][ t t ] respectvely t. Let s consder some proof steps of Lemma. A pror estmates (9) () provde us boundedness of sequences { yt ( )} n V( ) L ( ; L ( )) and n L ( ; L ( )). at s wy tere exst functon u from W( ) and subsequence of sequences { }{ } suc tat as and relatons (4) (6) old. en we take subsequences { }{ } suc tat a pror estmates (9) (3) te compactness teorem (see [8] lemmas.8.9) and (4) (6) provde us te lmt relaton (7). And relatons (8) (9) follow from te estmate () and nequaltes (7) () k ( x y y) K n L ( Q ) (8) p y k( x y ) K n Lp ( ). s e next we proved tat te lmt functon u satsfes () (3). For tat we perform te passage of te lmt n te followng equaton ( y( t)) v dxdt ( u ) v () dx t ( y( t)) v dsdt ( u ) v () ds t v v k x y y dxdt k x y dsdt y ( ) ( ) x s s f v dt f v dt wc s obtaned from (6) as z v ( v s an nterpolant from l constructed for te functon v C( ) C ( ): ( ) ) after te transformaton of te frst and te trd summands wt t elp of te summaton by parts formula and te use of extenson wt respect to varable t and ntegraton over []. As a result we obtan te followng equaton as d v ( u) v dxdt ( u ) v() dx K dxdt x dt V (9) 5

t Internatonal Conference on "Mes metods for boundary-value problems and applcatons" IOP Publsng IOP Conf. Seres: Materals Scence and Engneerng 58 (6) 4 do:.88/757-899x/58//4 d v ( u) v dsdt ( u ) v() ds K dsdt s (3) dt f v dt f v dt. In consderaton of C ( ) we rewrte (3) n te form ( u) vdx ( u) vds dt f v { } { t v v K dx f v K ds dt G( t) dt. } (3) x s From te equaton (3) follows tat functon Gt () s a generalzed dervatve from L( ) of te functon from t follows tat (see (3)) ( u) vdx ( u) vds v v J( u( t)) v f v K dx f v K ds s x almost everywere n []. Next we use te monotone metod to determne tat By usng t we prove tat as suc tat (3) v v { } x s K dx K ds dt v u v { } x s s k ( x u u) dx k ( x u ) ds dt. t t ( ˆ y( t) y( t)) k ( x y y) dxdt x ( ˆ y( t) y( t)) k ( x y y) dsdt s were s from lemma y ˆ( t) y( t ). e man result of ts paper s eorem. Let be a convex doman n te space condtons (9) (3) and let / (33) R functons f L ( ; W ( )) f L ( ; W ( )) u V L ( ) p p p p k ( ) satsfy k 6

t Internatonal Conference on "Mes metods for boundary-value problems and applcatons" IOP Publsng IOP Conf. Seres: Materals Scence and Engneerng 58 (6) 4 do:.88/757-899x/58//4 and condton (9) s satsfed. en te pecewse constant extensons subsequence of sceme soluton (6) (7) wc s satsfyng relatons (4) (9) converges to generalzed soluton of te problem () (3). In condtons of soluton unqueness (see [5]) te wole sequence as ts property. Acknowledgments e researc was supported by te Russan Foundaton for researc (projects 5--5686 5-4- 35). References [] Antontsev S N and Meyrmanov A M 979 Matematceske model covmestnogo dvgenya povernostny I podzemny vod (Novosbrsk: Novosbrsk State Unversty publsng) [] merbaev M R 4 Spaces wt a grap norm and strengtened Sobolev spaces I Russan Matmatcs 5 55 65 [3] merbaev M R 4 Spaces wt a grap norm and strengtened Sobolev spaces II Russan Matmatcs 9 46 53 [4] Glazyrna L L and Pavlova M F 994 On solvablty of certan problem on surface and ground water jont moton Russan Matmatcs 9 5 6 [5] Glazyrna L L and Pavlova M F A unqueness teorem for te soluton of a problem n te teory of te jon moton of cannel and underground Russan Matmatcs [6] Dautov R Z and Karcevsk M M Vvedene d teoryu metoda konecny elementov (Kazan: Kazan State Unversty publsng) [7] Korneev V G 977 Semy metoda konecny elementov vysok poryadkov tocnost (Lenngrad: Lengrad State Unversty publsng) [8] Alt H.W. and Luckaus S. 983 Quaslnear ellptc-parabolc dfferental equaton Mat.Z. 8 3 4 7