Numerical Simulation of Wave Propagation Using the Shallow Water Equations
|
|
- Brandon Harper
- 6 years ago
- Views:
Transcription
1 umercal Smulaton of Wave Propagaton Usng the Shallow Water Equatons Junbo Par Harve udd College 6th Aprl 007 Abstract The shallow water equatons SWE were used to model water wave propagaton n one dmenson and two dmensons. SWE were appromated b usng fnte dfference method. One dmensonal SWE were tested usng varous ntal condtons. and SWE were used to model tsunam wave propagaton near coast lne. model correctl predcts the behavor of tsunam water wave. model develops numercal oscllaton once the wave reaches shallow area. A real basn data from Gulf of eco has been tested to smulate the wave propagatng n the area. The shallow water equatons descrbe propagaton of water wave whose wavelength s much longer than the depth of water. Tsunams are eamples of such waves. Wavelength of a tpcal tsunam can eceed 00 m. Ths s to be compared wth the average depth of ocean whch s appromatel 5 m []. Therefore tsunam wave propagaton can be modeled usng shallow water equatons. In ths report we use a numercal method to smulate wave propagaton n one and two dmensons. The same scheme s used to model tsunam wave propagaton near a shore. Shallow Water Equatons The dervaton of the shallow water equatons follows from the conservaton of mass and momentum. In dervng the equatons however we must mae an assumpton about the water that s beng modeled. We assume that there s no pressure varaton n the vertcal drecton []. Therefore t cannot correctl model the behavor of a water fall. Ths naturall forces the wavelength of the water wave to be larger than the depth of the water. Furthermore n ths report we neglected frcton between the flud and the basn. However frcton can actuall be accounted for Imamura et al. have successfull modeled the frctonal contrbuton [3]. One dmensonal Shallow Water Equatons The shallow water equatons n one dmenson are = 0 t g = 0. t.a.b
2 where g s the gravtatonal constant s the total thcness of water at and s a quantt defned as product of depth averaged veloct and water veloct n the drecton [4]. Equaton.a and.b are manfestaton of mass and momentum conservaton law respectvel. As mentoned earler frctonal terms have been gnored n our dscusson of SWE. 0.5 Undsturbed Water Surface h basn as Fg. Water level profle of an eample wave The shallow water equatons are coupled frst order dfferental equatons that can be uncoupled to produce two second-order dfferental equatons. t = g h h t = g h g t h where h s the depth of the water the relaton between and h s = h. The form of equaton.a and.b s ver smlar to that of a wave equaton. The speed of and propagaton at gven s therefore determned b total thcness of water c s g = gh..a.b umercal Scheme for Shallow Water Equatons To solve the shallow water equatons numercall we frst dscretzed space and tme. We used fnte dfferences to appromate the spatal and tme dervatves. The tme dervatves were appromated b usng the second-order centered dfferences followng the Cran-colson method. The spatal dervatves were appromated b usng the second-order centered dfferences. For the one dmensonal SWE we used a sem-mplct and mplct method to compare the numercal stablt of soluton. For all smulatons we used spatal grd sze = /8 m and tme step sze t = /3000 sec.
3 t Fg. Spatal and tme dscretzaton. Spatal doman s dscretzed usng ponts = L/. The boundar condton mposed on = 0 and = L are gven as = 0 = ; = 0. 3 where s an appromaton to t. Equaton 3 sets -veloct of water to be zero at boundar. As a consequence / = 0 at the boundar. Sem-Implct ethod Usng fnte dfference method we can appromate equaton as = / cg h = c / 4 c / 4 * * 4.a 4.b where c = t/ the rato between tme step and spatal step. s an appromaton to t. Gven an ntal condton at tme t = 0 ths method can be used to compute and at subsequent tme steps. However ths numercal scheme s not an eplct method but rather a sem-mplct method. To compute we use s from prevous tme step. To compute however we use and whch s obtaned b solvng equaton 4.a for frst. * * Fg 3. Stencl of Sem-Implct method for computng.
4 Implct ethod To ncrease the numercal stablt of the soluton we have mplemented an mplct method on the one dmensonal shallow water equatons. However equaton.b s non-lnear snce = h. In order to mplement an mplct scheme.b had to be separated nto lnear part and non-lnear part. Applng Cran-colson method and fnte dfference method we obtan followng equaton: / 4cgh c / 4 = = c / 4 / cg / 4cgh 5 Left hand sdes of equaton 5 are and terms at tme step. The rght hand sdes are terms at tme step. Values of and at tme step can be computed b solvng a lnear sstem of equatons. Test of Shallow Water Equatons The shallow water equatons n one dmenson were tested wth three dfferent ntal condtons. In all cases the ntal veloct of the water was set to be zero water was at rest at t = 0 and therefore = 0. We vared t = 0 to eamne the results of numercal smulatons. The frst row n Fgure 4 llustrates the three dfferent t = 0 s. The solutons at later tme t are plotted n the second row. The solutons obtaned b usng sem-mplct method and eplct method are plotted on the same graph for each plot. ote that the solutons almost overlap each other. We observe that the soluton to ntal condton n Fgure 4a and 4b develop numercal nstabltes. However these numercal nstabltes are not nherent to shallow water equaton. Equaton s a dfferental form of shallow water equaton whch assumes that the water profle s smooth. Yet the ntal condtons gven are not necessarl smooth. Ths can be understood f we eamne equaton. From equaton whch resembles a wave equaton we can deduce that the soluton to shallow water equaton can be appromated wth superposton of snusods. Therefore the ntal profle of can be thought of as a superposton of snusods. The step functon however can be epressed wth ver hgh frequenc snusods. Those hgh frequenc components develop nto numercal oscllatons. Yet the orgnal shallow water equaton allows for dscontnuous solutons such as a bore [5]. To solve problems wth non-smooth profle fnte volume method should be used.
5 The numercal oscllaton we observe does not arse from nstablt of sem-mplct method ether. Implct method does not mprove the results. The solutons are nearl eactl the same. The numercal nstablt arses from the lmts of fnte dfference method t=0 9.6 t=0 9.6 t= Soluton to Shallow water equaton 0 Soluton to Shallow water equaton 0 Soluton to Shallow water equaton t= t= t= a b c F g 4. Three dfferent ntal condtons: step functon slope snusod are plotted n the frst row. The second row s plots of at later tme step. Results from usng eplct method and mplct method have been plotted. The mamum and mnmum value of t = 0 are 8m and 0m respectvel. The basn depth s 0m for all three cases. a t = 0 and t = b t = 0 and t = c t = 0 and t =. Fnte dfference method assumes that dervatves est at all ponts. Therefore t cannot handle dscontnuous ntal condton such as step functon. It also cannot produce stable solutons for ntal condtons descrbed n fgure 4b snce the dervatve of the ntal condton s dscontnuous. As observed n fgure 4a and 4b numercal nstablt ensues f the ntal condton profle s not twce dfferentable. When gven a snusod as an ntal condton the soluton remans stable for much longer perod. Tsunam wave n The shallow water equatons were used to model wave propagaton near the shore. The basn profle has been appromated b a hperbolc tangent gven as 70 h = tanh 8 0 m 00 m 6 such that the depth would range from 5m to 95m. We made the depth of water to be postve at all
6 snce ver hgh numercal error occurs f we let depth profle be negatve ths would correspond to above sea level. The ntal condton was gven as a Gaussan for and 0 = 0.5ep 0 0 = 50ep 0 /8; /8 = 00 0; 7 was gven a non-zero value to mae the wave propagate to drecton a b Fg 5. a Plot of h the depth profle used n the smulaton. b Intal condton for and 0.0. The result of ths smulaton s plotted n Fgure 6. As the wave approaches the shore and enters the shallower area the wave slows down. Ths s what was epected snce the wave propagaton speed s appromatel c s = g see equaton. As the wave slows down the wdth of the Gaussan wave pacet s reduced. Ths results n ncreased ampltude of the water wave. Soluton to Shallow water equaton Soluton to Shallow water equaton Soluton to Shallow water equaton.5 t=0.5 t= t= depth profle/00m - depth profle/00m - depth profle/00m a b c Fg 6. at varous tme step. The depth profle has been scaled down b a factor of 00 a t = 0. b t =.37 a t = Two mensonal SWE Shallow water equatons n are etensons to equatons. The set of equatons are gven as
7 8.a here s product of. Equaton 8.a s derved om mass conservaton law. Equaton 8.b and 8.c are derved from momentum conservaton l SWE E was ver smlar to that of one a grd. The sze of the spatal doman t me dervatves were 9.a 9.b 9.c alculatng and frst usng equaton 9.b and 9.c. 8.b 8.c 0 g g 0 g g 0 = = = h t h t t w and depth averaged veloct n drecton [6] fr law. Frcton has been gnored n these equatons. umercal Scheme for Two mensona The numercal method used for two dmensonal SW dmensonal case. Space and tme have been dscretzed nto was determned to be a 00m b 00m. The spacng between each grd pont and was se to be m. The sze of the tme step was set to be t = /4000 sec. We used fnte dfference method to appromate SWE. The spatal dervatves were appromated b usng the second-order centered dfferences. The t appromated b usng the second-order centered dfferences followng the Cran-colson method. Equaton 8 therefore can be appromated b c / g c / g / c / c / g c / g / c / 4c / 4c c / h h * 4c / * 4 / = = = Agan the * n equaton 9.a ndcates that the values of and at tme step are from c The boundar condtons are smlar to that of case. and are zero on the
8 - boundares ths s also nown as no slp condton. Water cannot move nto or out of a boundar. It also cannot move along a boundar. Fg 7. Laout of dscretzed grd. The value of at the boundar ponts blue dot are alwas equated to the value of nsde the boundar red dot at each tme step. Smple odel: Hperbolc Tangent Basn The shallow water equatons were used to smulated water wave near a shore. The sea basn was modeled as a hperbolc tangent. The depth profle was gven as 50 h = tanh 0 The ntal condton was gven as 30 0 = 0.5ep = The result of numercal smulaton s shown n fgure 8. The ntal Gaussan profle pans out nto the - plane n all drecton. Bul of the water however travels toward drecton. There are features that are common to the and results. The waves slows down a b
9 a b c Fg 8. umercal smulaton of water propagaton near coast lne. a Soluton at t = 0. b Soluton at t = sec. c Soluton at t = 3 sec as t approaches the shore. ote that the wave front barel moved between t = and t = 3. Also the ampltude of the water ncrease as the wave pacet approaches the shore. Unle the case however we see more severe numercal oscllaton occurrng as the wave approaches the shore. umercal oscllaton s vsble n soluton at t = and t becomes more pronounced as the wave front approaches shallower parts of the shore at t = 3. Implct method has not been mplemented to chec f ths was an artfact of sem-mplct scheme that we used to solve two dmensonal SWE. However we now from the results of SWE smulatons that usng mplct method does not guarantee a stable soluton. One possble source of numercal oscllaton ma be the dscrepanc between analtcal doman of dependence and numercal doman of dependence. In a hperbolc equaton nformaton travels at the speed of wave propagaton. Fgure 9 llustrates the doman of dependence of two nd: analtc and numercal. umercal doman of dependence s determned b the grd spacng and the method used to appromate the dervatves. Analtc doman of dependence s determned b the speed of the wave propagaton c s and the tme step sze. When two doman of dependence dffer from another greatl numercal error s nduced. In hperbolc equatons ths leads to dspersve behavor of the numercal soluton [7]. Waves wth dfferent wavelength travel at dfferent speed.
10 t t c t F g 9. Schematc dagram of analtc doman of dependence and numercal doman of dependence. In ths case analtcal doman of dependence s narrower than the numercal doman of dependence. However as equaton suggests shallow water equatons have the form of a hperbolc equaton. Furthermore c pth of the water: c s s determned b the de s = g. Therefore as the waves approach the shallower regon the analtcal doman of dependence becomes narrower than numercal doman of dependence. Ths results n dspersve behavor whch maes waves wth dfferent wave length travel at dfferent speed. Snce a Gaussan wave pacet s a superposton of man snusods the wave breas down nto dscrete waves travelng at dfferent speed. The result s numercal oscllaton whch s a manfestaton of nterference between varous snusodal waves. Gulf of eco The two dmensonal shallow water equatons were further tested usng a real basn data whch can be obtaned from atonal Geophscal ata Center webste [8]. We obtaned the depth profle of sea basn near the Gulf of eco. The bathmetrc profle was sampled at dscrete ponts separated b appromatel 4000m. The sze of the basn area was 400 m b 400 m. The depth ranged from 50m to m. The raw bathmetr data however could not be used because the profle was not smooth enough. The agged terrans whch s a product of coarse resoluton of grd cell = 4 m was a maor source of numercal nstablt. In order to smooth out the terran and et retan the orgnal profle we appled convoluton to each data pont. h * = S l l where h * s the smoothed depth profle at = and =. S s gven as h l = ep. π 4 S 3
11 Equaton s a dscrete form of convoluton. The result of ths s shown n Fgure 0. a b Fg 0. a Raw depth profle. b Product obtaned after applng convoluton The ntal condton for ths smulaton was gven as = 0.5ep = The results of numercal smulaton are plotted n Fgure. a b
12 c Fg. Wave propagaton n Gulf of eco a Intal condton and sea basn b Soluton at t = 400 sec. c Soluton at t = 550 sec. The lateral rpples are not due to numercal oscllaton. The rpples arse from the shape of the basn. Compare ths wth fgure 8c. The lateral rpples observed n Fgure c are due to the shape of basn. Ths tpe of rpples does not develop when the depth profle s a smple hperbolc tangent gven n equaton 0. Such rpples are not observed n Fgure 8c. The curves observed n Fgure 8c s due to waves that reflected from the = 0 and = 00m boundar. The soluton however qucl becomes unstable as soon as the waves reflected from = 0 catches up wth the orgnal wave n Fgure b the orgnal wave s labeled wth red crcle and the reflected wave s labeled wth blue crcle. Interference between two waves produces numercal oscllatons. Acnowledgement I would le to than Professor Yong for provdng me wth resources and gudng me through out the proect. I than Sean eenehan for gudng me to fnd resources about bathmetr data. I than c Alger for hs dea of applng convoluton to raw bathmetrc data.
13 References [] [] Kevoran J. Partal fferental Equatons: Analtc Soluton Technques. Pacfc Grove Calf. c990. pp 0-9. [3] Imamura F. Yalcner A.C. Tsunam odelng anual raft pp [4] ImamuraF. et al. pp 9 [5] Kervoran pp. -3. [6] ImamuraF. et al. pp 6-7 [7] Yong. Prvate Communcaton Scentfc Computng Lecture 007. [8]
PART 8. Partial Differential Equations PDEs
he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationNew Method for Solving Poisson Equation. on Irregular Domains
Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationE91: Dynamics. E91: Dynamics. Numerical Integration & State Space Representation
E91: Dnamcs Numercal Integraton & State Space Representaton The Algorthm In steps of Δ f ( new ) f ( old ) df ( d old ) Δ Numercal Integraton of ODEs d d f() h Δ Intal value problem: Gven the ntal state
More informationLecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES
COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve
More informationHigher Order Wall Boundary Conditions for Incompressible Flow Simulations
THE 5 TH ASIAN COMPUTAITIONAL FLUID DYNAMICS BUSAN KOREA OCTOBER 7-30 003 Hgher Order Wall Boundary Condtons for Incompressble Flow Smulatons Hdetosh Nshda. Department of Mechancal and System Engneerng
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationChapter 4 The Wave Equation
Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a second-order lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as
More informationGrid Generation around a Cylinder by Complex Potential Functions
Research Journal of Appled Scences, Engneerng and Technolog 4(): 53-535, 0 ISSN: 040-7467 Mawell Scentfc Organzaton, 0 Submtted: December 0, 0 Accepted: Januar, 0 Publshed: June 0, 0 Grd Generaton around
More informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
More informationNumerical Transient Heat Conduction Experiment
Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use
More informationChapter 8. Potential Energy and Conservation of Energy
Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal
More informationAdditional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty
Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More information1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave
014 4th Internatonal Conference on Future nvronment and nergy IPCB vol.61 (014) (014) IACSIT Press, Sngapore I: 10.776/IPCB. 014. V61. 6 1-mensonal Advecton-ffuson Fnte fference Model ue to a Flow under
More informationChapter Newton s Method
Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve
More informationGravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11)
Gravtatonal Acceleraton: A case of constant acceleraton (approx. hr.) (6/7/11) Introducton The gravtatonal force s one of the fundamental forces of nature. Under the nfluence of ths force all objects havng
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationOne Dimensional Axial Deformations
One Dmensonal al Deformatons In ths secton, a specfc smple geometr s consdered, that of a long and thn straght component loaded n such a wa that t deforms n the aal drecton onl. The -as s taken as the
More informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More informationSeptic B-Spline Collocation Method for the Numerical Solution of the Modified Equal Width Wave Equation
Appled Mathematcs 79-749 do:.46/am..698 Publshed Onlne June (http://www.scrp.org/ournal/am) Septc B-Splne Collocaton Method for the umercal Soluton of the Modfed Equal Wdth Wave Equaton Abstract Turab
More informationNormally, in one phase reservoir simulation we would deal with one of the following fluid systems:
TPG4160 Reservor Smulaton 2017 page 1 of 9 ONE-DIMENSIONAL, ONE-PHASE RESERVOIR SIMULATION Flud systems The term sngle phase apples to any system wth only one phase present n the reservor In some cases
More informationNUMERICAL MODELLING OF TSUNAMI WAVE EQUATIONS MASTER OF SCIENCE IN MATHEMATICS
NUMERICAL MODELLING OF TSUNAMI WAVE EQUATIONS A thess submtted n partal fulfllment of the requrements for the award of the degree of MASTER OF SCIENCE IN MATHEMATICS submtted by HARI SHANKAR SHAW Roll
More information2.29 Numerical Fluid Mechanics
REVIEW Lecture 10: Sprng 2015 Lecture 11 Classfcaton of Partal Dfferental Equatons PDEs) and eamples wth fnte dfference dscretzatons Parabolc PDEs Ellptc PDEs Hyperbolc PDEs Error Types and Dscretzaton
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationANALYSIS OF ELECTROMAGNETIC FIELD USING THE CONSTRAINED INTERPOLATION PROFILE METHOD PHÂN TÍCH TRƯỜNG ĐIỆN TỪ SỬ DỤNG PHƯƠNG PHÁP CIP
ANALYSIS OF ELECTROMAGNETIC FIELD USING THE CONSTRAINED INTERPOLATION PROFILE METHOD PHÂN TÍCH TRƯỜNG ĐIỆN TỪ SỬ DỤNG PHƯƠNG PHÁP CIP LÊ VŨ HƯNG Cao đẳng kỹ thuật quốc ga Kushro, Nhật Bản Kushro Natonal
More informationFinite Difference Method
7/0/07 Instructor r. Ramond Rump (9) 747 698 rcrump@utep.edu EE 337 Computatonal Electromagnetcs (CEM) Lecture #0 Fnte erence Method Lecture 0 These notes ma contan coprghted materal obtaned under ar use
More informationOPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming
OPTIMIATION Introducton ngle Varable Unconstraned Optmsaton Multvarable Unconstraned Optmsaton Lnear Programmng Chapter Optmsaton /. Introducton In an engneerng analss, sometmes etremtes, ether mnmum or
More informationFour Bar Linkages in Two Dimensions. A link has fixed length and is joined to other links and also possibly to a fixed point.
Four bar lnkages 1 Four Bar Lnkages n Two Dmensons lnk has fed length and s oned to other lnks and also possbly to a fed pont. The relatve velocty of end B wth regard to s gven by V B = ω r y v B B = +y
More informatione i is a random error
Chapter - The Smple Lnear Regresson Model The lnear regresson equaton s: where + β + β e for,..., and are observable varables e s a random error How can an estmaton rule be constructed for the unknown
More informationAPPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS
Unversty of Oulu Student Laboratory n Physcs Laboratory Exercses n Physcs 1 1 APPEDIX FITTIG A STRAIGHT LIE TO OBSERVATIOS In the physcal measurements we often make a seres of measurements of the dependent
More informationDigital Signal Processing
Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over
More informationChapter 3 Differentiation and Integration
MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton
More informationFTCS Solution to the Heat Equation
FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use
More informationTHE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES
THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES W. C. Lao Department of Cvl Engneerng, Feng Cha Unverst 00 Wen Hwa Rd, Tachung, Tawan SUMMARY: The ndentaton etween clndrcal ndentor
More informationNumerical Solution of Boussinesq Equations as a Model of Interfacial-wave Propagation
BULLETIN of the Malaysan Mathematcal Scences Socety http://math.usm.my/bulletn Bull. Malays. Math. Sc. Soc. (2) 28(2) (2005), 163 172 Numercal Soluton of Boussnesq Equatons as a Model of Interfacal-wave
More informationRECEIVED. Negative Transverse Impedance
RECEVED SEP 2 3 996 OSTt > LS- 4 O C a f L W. Chou March 2, 989 (Rev. June 2, 9S9) Negatve Transverse mpedance ntroducton n Ref. ( we report an observaton that the horzontal and the vertcal loss factors
More informationPhysics 2A Chapter 3 HW Solutions
Phscs A Chapter 3 HW Solutons Chapter 3 Conceptual Queston: 4, 6, 8, Problems: 5,, 8, 7, 3, 44, 46, 69, 70, 73 Q3.4. Reason: (a) C = A+ B onl A and B are n the same drecton. Sze does not matter. (b) C
More informationWeek 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product
The Vector Product Week 11: Chapter 11 Angular Momentum There are nstances where the product of two vectors s another vector Earler we saw where the product of two vectors was a scalar Ths was called the
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More informationModelli Clamfim Equazione del Calore Lezione ottobre 2014
CLAMFIM Bologna Modell 1 @ Clamfm Equazone del Calore Lezone 17 15 ottobre 2014 professor Danele Rtell danele.rtell@unbo.t 1/24? Convoluton The convoluton of two functons g(t) and f(t) s the functon (g
More informationACTM State Calculus Competition Saturday April 30, 2011
ACTM State Calculus Competton Saturday Aprl 30, 2011 ACTM State Calculus Competton Sprng 2011 Page 1 Instructons: For questons 1 through 25, mark the best answer choce on the answer sheet provde Afterward
More informationApplication to Plane (rigid) frame structure
Advanced Computatonal echancs 18 Chapter 4 Applcaton to Plane rgd frame structure 1. Dscusson on degrees of freedom In case of truss structures, t was enough that the element force equaton provdes onl
More informationProcedia Computer Science
Avalable onlne at www.scencedrect.com Proceda Proceda Computer Computer Scence Scence 1 (01) 00 (009) 589 597 000 000 Proceda Computer Scence www.elsever.com/locate/proceda Internatonal Conference on Computatonal
More informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
More informationCHAPTER 4. Vector Spaces
man 2007/2/16 page 234 CHAPTER 4 Vector Spaces To crtcze mathematcs for ts abstracton s to mss the pont entrel. Abstracton s what makes mathematcs work. Ian Stewart The man am of ths tet s to stud lnear
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationDEMO #8 - GAUSSIAN ELIMINATION USING MATHEMATICA. 1. Matrices in Mathematica
demo8.nb 1 DEMO #8 - GAUSSIAN ELIMINATION USING MATHEMATICA Obectves: - defne matrces n Mathematca - format the output of matrces - appl lnear algebra to solve a real problem - Use Mathematca to perform
More information6 Supplementary Materials
6 Supplementar Materals 61 Proof of Theorem 31 Proof Let m Xt z 1:T : l m Xt X,z 1:t Wethenhave mxt z1:t ˆm HX Xt z 1:T mxt z1:t m HX Xt z 1:T + mxt z 1:T HX We consder each of the two terms n equaton
More informationONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00
ONE IMENSIONAL TRIANGULAR FIN EXPERIMENT Techncal Advsor: r..c. Look, Jr. Verson: /3/ 7. GENERAL OJECTIVES a) To understand a one-dmensonal epermental appromaton. b) To understand the art of epermental
More informationSolution Thermodynamics
Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationSIMULATION OF SOUND WAVE PROPAGATION IN TURBULENT FLOWS USING A LATTICE-BOLTZMANN SCHEME. Abstract
SIMULATION OF SOUND WAVE PROPAGATION IN TURBULENT FLOWS USING A LATTICE-BOLTZMANN SCHEME PACS REFERENCE: 43.20.Mv Andreas Wlde Fraunhofer Insttut für Integrerte Schaltungen, Außenstelle EAS Zeunerstr.
More informationSolutions to Homework 7, Mathematics 1. 1 x. (arccos x) (arccos x) 1
Solutons to Homework 7, Mathematcs 1 Problem 1: a Prove that arccos 1 1 for 1, 1. b* Startng from the defnton of the dervatve, prove that arccos + 1, arccos 1. Hnt: For arccos arccos π + 1, the defnton
More informationMACHINE APPLIED MACHINE LEARNING LEARNING. Gaussian Mixture Regression
11 MACHINE APPLIED MACHINE LEARNING LEARNING MACHINE LEARNING Gaussan Mture Regresson 22 MACHINE APPLIED MACHINE LEARNING LEARNING Bref summary of last week s lecture 33 MACHINE APPLIED MACHINE LEARNING
More informationTurbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH
Turbulence classfcaton of load data by the frequency and severty of wnd gusts Introducton Oscar Moñux, DEWI GmbH Kevn Blebler, DEWI GmbH Durng the wnd turbne developng process, one of the most mportant
More informationCopyright 2017 by Taylor Enterprises, Inc., All Rights Reserved. Adjusted Control Limits for P Charts. Dr. Wayne A. Taylor
Taylor Enterprses, Inc. Control Lmts for P Charts Copyrght 2017 by Taylor Enterprses, Inc., All Rghts Reserved. Control Lmts for P Charts Dr. Wayne A. Taylor Abstract: P charts are used for count data
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationChapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise.
Chapter - The Smple Lnear Regresson Model The lnear regresson equaton s: where y + = β + β e for =,..., y and are observable varables e s a random error How can an estmaton rule be constructed for the
More informationTools for large-eddy simulation
Center for Turbulence Research Proceedngs of the Summer Program 00 117 Tools for large-eddy smulaton By Davd A. Caughey AND Grdhar Jothprasad A computer code has been developed for solvng the ncompressble
More informationResource Allocation and Decision Analysis (ECON 8010) Spring 2014 Foundations of Regression Analysis
Resource Allocaton and Decson Analss (ECON 800) Sprng 04 Foundatons of Regresson Analss Readng: Regresson Analss (ECON 800 Coursepak, Page 3) Defntons and Concepts: Regresson Analss statstcal technques
More informationComputational Geodynamics: Advection equations and the art of numerical modeling
540 Geodynamcs Unversty of Southern Calforna, Los Angeles Computatonal Geodynamcs: Advecton equatons and the art of numercal modelng Sofar we manly focussed on dffuson equaton n a non-movng doman. Ths
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
More informationNumerical Methods. ME Mechanical Lab I. Mechanical Engineering ME Lab I
5 9 Mechancal Engneerng -.30 ME Lab I ME.30 Mechancal Lab I Numercal Methods Volt Sne Seres.5 0.5 SIN(X) 0 3 7 5 9 33 37 4 45 49 53 57 6 65 69 73 77 8 85 89 93 97 0-0.5 Normalzed Squared Functon - 0.07
More informationCIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M
CIS56: achne Learnng Lecture 3 (Sept 6, 003) Preparaton help: Xaoyng Huang Lnear Regresson Lnear regresson can be represented by a functonal form: f(; θ) = θ 0 0 +θ + + θ = θ = 0 ote: 0 s a dummy attrbute
More informationSimulation for Pedestrian Dynamics by Real-Coded Cellular Automata (RCA)
Smulaton for Pedestran Dynamcs by Real-Coded Cellular Automata (RCA) Kazuhro Yamamoto 1*, Satosh Kokubo 1, Katsuhro Nshnar 2 1 Dep. Mechancal Scence and Engneerng, Nagoya Unversty, Japan * kazuhro@mech.nagoya-u.ac.jp
More informationUnit 5: Quadratic Equations & Functions
Date Perod Unt 5: Quadratc Equatons & Functons DAY TOPIC 1 Modelng Data wth Quadratc Functons Factorng Quadratc Epressons 3 Solvng Quadratc Equatons 4 Comple Numbers Smplfcaton, Addton/Subtracton & Multplcaton
More informationThis model contains two bonds per unit cell (one along the x-direction and the other along y). So we can rewrite the Hamiltonian as:
1 Problem set #1 1.1. A one-band model on a square lattce Fg. 1 Consder a square lattce wth only nearest-neghbor hoppngs (as shown n the fgure above): H t, j a a j (1.1) where,j stands for nearest neghbors
More informationIV. Performance Optimization
IV. Performance Optmzaton A. Steepest descent algorthm defnton how to set up bounds on learnng rate mnmzaton n a lne (varyng learnng rate) momentum learnng examples B. Newton s method defnton Gauss-Newton
More informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
More informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
More informationSIMPLE LINEAR REGRESSION
Smple Lnear Regresson and Correlaton Introducton Prevousl, our attenton has been focused on one varable whch we desgnated b x. Frequentl, t s desrable to learn somethng about the relatonshp between two
More informationLine Drawing and Clipping Week 1, Lecture 2
CS 43 Computer Graphcs I Lne Drawng and Clppng Week, Lecture 2 Davd Breen, Wllam Regl and Maxm Peysakhov Geometrc and Intellgent Computng Laboratory Department of Computer Scence Drexel Unversty http://gcl.mcs.drexel.edu
More informationSolution Set #1
05-78-0 Soluton Set #. Fnd epressons and setch the results of the followng operatons: (a) COMB RECT The spacng of the elements of the COMB functon matches the wdth of the rectangle; we can do ths n ether
More informationMATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018
MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng
More informationKinematics in 2-Dimensions. Projectile Motion
Knematcs n -Dmensons Projectle Moton A medeval trebuchet b Kolderer, c1507 http://members.net.net.au/~rmne/ht/ht0.html#5 Readng Assgnment: Chapter 4, Sectons -6 Introducton: In medeval das, people had
More informationQuadratic speedup for unstructured search - Grover s Al-
Quadratc speedup for unstructured search - Grover s Al- CS 94- gorthm /8/07 Sprng 007 Lecture 11 001 Unstructured Search Here s the problem: You are gven a boolean functon f : {1,,} {0,1}, and are promsed
More informationHomework Assignment 3 Due in class, Thursday October 15
Homework Assgnment 3 Due n class, Thursday October 15 SDS 383C Statstcal Modelng I 1 Rdge regresson and Lasso 1. Get the Prostrate cancer data from http://statweb.stanford.edu/~tbs/elemstatlearn/ datasets/prostate.data.
More informationMath1110 (Spring 2009) Prelim 3 - Solutions
Math 1110 (Sprng 2009) Solutons to Prelm 3 (04/21/2009) 1 Queston 1. (16 ponts) Short answer. Math1110 (Sprng 2009) Prelm 3 - Solutons x a 1 (a) (4 ponts) Please evaluate lm, where a and b are postve numbers.
More information( ) [ ( k) ( k) ( x) ( ) ( ) ( ) [ ] ξ [ ] [ ] [ ] ( )( ) i ( ) ( )( ) 2! ( ) = ( ) 3 Interpolation. Polynomial Approximation.
3 Interpolaton {( y } Gven:,,,,,, [ ] Fnd: y for some Mn, Ma Polynomal Appromaton Theorem (Weerstrass Appromaton Theorem --- estence ε [ ab] f( P( , then there ests a polynomal
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationIrregular vibrations in multi-mass discrete-continuous systems torsionally deformed
(2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected
More informationχ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More informationCHAPTER 7 CONSTRAINED OPTIMIZATION 2: SQP AND GRG
Chapter 7: Constraned Optmzaton CHAPER 7 CONSRAINED OPIMIZAION : SQP AND GRG Introducton In the prevous chapter we eamned the necessary and suffcent condtons for a constraned optmum. We dd not, however,
More informationRectilinear motion. Lecture 2: Kinematics of Particles. External motion is known, find force. External forces are known, find motion
Lecture : Kneatcs of Partcles Rectlnear oton Straght-Lne oton [.1] Analtcal solutons for poston/veloct [.1] Solvng equatons of oton Analtcal solutons (1 D revew) [.1] Nuercal solutons [.1] Nuercal ntegraton
More informationTemperature. Chapter Heat Engine
Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the
More informationBinding energy of a Cooper pairs with non-zero center of mass momentum in d-wave superconductors
Bndng energ of a Cooper pars wth non-zero center of mass momentum n d-wave superconductors M.V. remn and I.. Lubn Kazan State Unverst Kremlevsaa 8 Kazan 420008 Russan Federaton -mal: gor606@rambler.ru
More informationDepartment of Electrical & Electronic Engineeing Imperial College London. E4.20 Digital IC Design. Median Filter Project Specification
Desgn Project Specfcaton Medan Flter Department of Electrcal & Electronc Engneeng Imperal College London E4.20 Dgtal IC Desgn Medan Flter Project Specfcaton A medan flter s used to remove nose from a sampled
More informationSCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors.
SCALARS AND ECTORS All phscal uanttes n engneerng mechancs are measured usng ether scalars or vectors. Scalar. A scalar s an postve or negatve phscal uantt that can be completel specfed b ts magntude.
More informationEffect of loading frequency on the settlement of granular layer
Effect of loadng frequency on the settlement of granular layer Akko KONO Ralway Techncal Research Insttute, Japan Takash Matsushma Tsukuba Unversty, Japan ABSTRACT: Cyclc loadng tests were performed both
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More information