Journal of Computational Physics 149, (1999) Article ID jcph , available online at
|
|
- Rudolf Malone
- 5 years ago
- Views:
Transcription
1 Joural of Computatioal Physics 149, (1999) Article ID jcph , available olie at o NOTE Defiig Wave Amplitude i Characteristic Boudary Coditios Key Words: Euler compressible equatios; characteristic boudary coditios; oreflectig coditios; iitial coditios. Characteristic treatmet of boudary coditios for the Euler equatios relies o determiig the stregth of the waves eterig the computatioal domai as a fuctio of the stregth of the outgoig waves ad the physical boudary coditios. The purpose of this ote is to demostrate how critical the defiitio chose for the wave amplitudes ca be. The 2D Euler equatios may be expressed i quasi-liear form as V + A V x + B V = 0. (1) y Here V = (ρ, u,v,p) T is the vector of primitive variables ad each of the matrices A ad B has its ow complete set of real eigevalues ad right ad left eigevectors. The matrix E defied as A x + B y ca be itroduced, where is chose as the outward ormal to the boudary uder cosideratio. By diagoalizig E the eigevalue matrix Λ = L E L 1 = diag ( λ 1,λ2,λ3,λ4 ) =diag(u, u, u + c, u c), (2) is obtaied, where u = u ad c is the speed of soud. The matrices L (L 1 ) with left (right) eigevectors as rows (colums) relate variatios i the characteristic variables W to variatios i the primitive vector V through the relatios δw = L δv, δv = L 1 δw. (3) I 2D, the four characteristic variables satisfy a set of covectio equatios with the speed of propagatio give by Eq. (2), with source terms related to pressure ad velocity variatios i the s-directio, where s forms a orthoormal basis (, s) with. These equatios are /99 $30.00 Copyright c 1999 by Academic Press All rights of reproductio i ay form reserved. 418
2 WAVE AMPLITUDES 419 obtaied by premultiplyig Eq. (1) by L. The fourth equatio reads + (u c) W 4 + cs W 2 = 0. (4) Applyig a explicit Euler time discretizatio to Eq. (1), the update of primitive variables ca be writte as [ V = V +1 V = tr= t A V ] x +B V. (5) y For a give boudary with ormal, the full residual R of Eq. (5) ca be split ito a ormal compoet R (ivolvig oly ormal derivatives) ad a tagetial compoet R s (ivolvig oly derivatives alog s). Let us defie V as the boudary value at time level, ad δv P the predicted boudary update from the iterior scheme, prior to applicatio of the boudary coditio. Defie also δvw P as the compoet of δvp to which the boudary coditio will be applied, ad δv U = δv P δvw P as the part of the boudary update which is ot affected by the characteristic boudary coditio. Typically, a characteristic based boudary treatmet is applied as follows: 1. Choose the part of the residual (δvw P ) to which the boudary coditios are to be applied. If it is the complete residual, δv U = From Eq. (3), decompose δvw P, ito characteristic variatios δwi,p igoig ad outgoig waves, with correspodig primitive variatios δv i,p w 3. Modify the amplitude of the icomig wave(s) δw i,p ad δw out due to ad δv out w. accordig to the physical. Retai the requiremets at the boudary. This produces the corrected amplitudes, δw i,c outgoig waves δw out or δv out w as they are. 4. Combie the waves δw i,c primitive variables. This gives δv C w ad δw out, ad usig L 1, Eq. (3), trasform back to. The boudary poit is the updated as V +1 = V + δv U + δv C w = V + δv U + δv i,c + δv out. The decompositio of the Euler equatios ito a set of waves travelig ormally to the boudary provides a theoretical basis to derive proper boudary coditio treatmets, followig steps above. However, such theory gives o idicatio of the best defiitio for the part of the update to which the boudary coditios are to be applied. This is the pricipal reaso why so may differet formulatios are discussed i the literature [1 8]. Let us defie a approach, we shall call the full residual approach, as a boudary treatmet such that i step 1 above δvw P = tr. Followig steps 2 4, i the case of a 2D subsoic outlet, this leads to the oreflectig boudary coditio = 0 = u + 1 P ρc which is equivalet to that proposed i [1 4]. Similarly, we shall defie the ormal approach as the boudary treatmet such that for step 1, δv P w = tr. This leads to the oreflectig coditio (u c) [ = 0 = (u c) u + 1 ρc ] P (6) (7)
3 420 F. NICOUD which is equivalet to the forms i [5, 6], although preseted i a completely differet formalism. Hirsch [7] argues that the oreflectig coditio has to be applied to the advectio terms of the bicharacteristic equatios (u c) + u s = 0. (8) Followig Giles [8], the aalysis of the liearized Euler equatios based o a Fourier decompositio of the solutio at the boudary gives W 2 = u u s ; (9) Eqs. (6), (7), (8), ad (9) are all oreflectig boudary coditios based o characteristic aalysis. I 1D, they all reduce to oe of the forms (u c) ( /) = 0or / = 0 which are equivalet sice the last characteristic equatio is simply / + (u c) ( /) i this case. However, these boudary treatmets are ot equivalet i 2D. It traspires that uder certai circumstaces they ca eve produce completely differet results. For example, cosider the computatio domai defied spatially by 0 < x < 1 ad 0 < y < 1. The iitial coditio is uiform for the desity ad the static pressure, ad zero for the velocity i the y-directio. For the streamwise velocity, we impose u(x, y) = U 0 (1.5 + tah(10(y 0.5))) for x = 0 ad u(y) = 0 elsewhere. U 0 is chose such that the flow is subsoic everywhere. The full residual approach is used at the ilet to impose the velocity compoets ad the temperature while a oreflectig coditio is tested at the outlet. The ormal oreflectig characteristic coditio is used for both y = 0 ad y = 1 to allow acoustic disturbaces i the y-directio to leave the domai. The velocity profile is expected to propagate dowstream durig the computatio. The steady solutio is obviously u(x, y) = u(0, y) for all x. Typical velocity profiles obtaied after covergece with both the full residual ad ormal formulatios at the outlet boudary are show i Fig. 1. Clearly the full residual outlet coditio prevets the give velocity profile from propagatig alog the x-directio. Istead, the u-velocity teds to be uiform ear the exit. O the other had, the use of the ormal approach leads to the correct velocity profiles. Both the Hirsh ad the Giles formulatios allow the hyperbolic taget profile to propagate as FIG. 1. Velocity profiles at differet abscissa for the full residual (left) ad the ormal (right) o-reflectig outlet boudary coditio.
4 WAVE AMPLITUDES 421 FIG. 2. Time evolutios of the streamwise (left) ad the ormal (right) velocity compoet at a poit at the outlet boudary for the Giles (dashed lie) ad ormal (solid lie) formulatios. expected, but the relaxatio time to the steady state is loger with the Giles coditio (see Fig. 2). These dramatic differeces may be explaied by formulatig all the boudary coditios i the same framework. Ay (oreflectig) boudary coditio ca be writte either i terms of time derivatives (temporal form) or i terms of ormal derivatives (spatial form). These two forms are liked through the compatibility relatio, Eq. (4), ad imposig a boudary coditio o the time derivative ca be traslated ito a coditio o the ormal derivative ad vice versa. A overview of the various coditios cosidered, writte i their two equivalet forms (temporal ad spatial), is give i Table I. This table provides a formal compariso of these boudary coditios. Of course, the results of a computatio depeds oly o the choice of the boudary (the rows i the table) ad ot o the form uder which it is writte (the colums i the table). The full residual formulatio imposes W 4 / = 0 ad thus forces the temporal evolutios of the streamwise velocity ad the pressure to remai early proportioal. At the iitial time, both quatities are uiform at the exit, so that their profiles keep the same shape durig the computatio if 1/ρc does ot deped o y. This feature of the solutio at the boudary is well predicted by the computatio (ot show) but is ot compatible with the preset physical cofiguratio. Actually the full residual approach is likely to give icorrect results as soo as the iitial coditios at the boudary are icosistet with the actual flow structure. The temporal forms of the other three boudary treatmets (see Table I) show that the temporal evolutio TABLE I Correspodece betwee the Temporal ad the Spatial Form for Some Noreflectig Boudary Coditios; 2D Case Name Temporal form Spatial form Thompso [4] Poisot [6] Hirsh [7] Giles [8] = 0, Eq. (6) W = u 4 s = c W2 = u W 2 u s, Eq. (9) = 1 u u c s W4 +c W2 =0, Eq. (7) = us (u c) = W 2 +c W2,Eq. (8)
5 422 F. NICOUD of the velocity ad the pressure are o loger proportioal. Three differet terms appear i the right-had side of the temporal forms, amely A = u s ( /), B = c( W 2 /) ad C = u ( W 2 /). At least oe term is eeded to esure that the give boudary coditio ca reach the correct steady state. For Hirsh s coditio which leads to results that are equivalet to the ormal approach, the term deoted above as A is ot critical (this term is the differece betwee two boudary coditios which give the same results; see Table I). Thus the term B is resposible for the success of the computatios with those coditios. Oe observes also that the term C i the Giles treatmet is othig but a factor of u /c smaller tha B i the preset subsoic test case. Accordigly, with the Giles treatmet the relaxatio time to the steady state has bee foud to be loer tha for the other coditios (see Fig. 2). More details of the preset study are available i [9], icludig the 3D versio of Table I ad implemetatio details of the boudary treatmets for a flow solver based o hybrid meshes. Some prelimiary rus were performed by Dr. G. Heradez. All the computatios were doe with the Fortra library AVBP/COUPL developed at CERFACS. REFERENCES 1. S. Chakravarthy, Euler equatios implicit scheme ad boudary coditio, AIAA J. 21, 699 (1983). 2. M. Hayder ad E. Turkel, Noreflectig boudary coditios for jet flow computatios, AIAA J. 33(12), 2264 (1995). 3. M. Hayder ad E. Turkel, High Order Accurate Solutios of Viscous Problems, AIAA Paper , K. W. Thompso, Time depedet boudary coditios for hyperbolic systems. II. J. Comput. Phys. 89, 439 (1990). 5. K. W. Thompso, Time depedet boudary coditios for hyperbolic systems, J. Comput. Phys. 68, 1 (1987). 6. T. J. Poisot ad S. K. Lele, Boudary coditios for direct simulatios of compressible viscous flows, J. Comput. Phys. 101, 104 (1991). 7. C. Hirsh, Numerical Computatio of Iteral ad Exteral Flow (Wiley, New York, 1990), Vol M. Giles, No-reflectig boudary coditios for euler equatio calculatio, AIAA J. 28(12), 2050 (1990). 9. F. Nicoud, O the Amplitude of the Waves i the Characteristic Boudary Coditios, Techical Report TR/CFD/98/21, CERFACS, Received Jauary 20, 1998; revised August 10, 1998 F. Nicoud 1 CERFACS Cetre Europée de Recherche et de Formatio Avacée e Calcul Scietifique 42, Aveue Gaspard Coriolis Toulouse cedex, Frace icoud@cerfacs.fr 1 Preset address: Ceter for Turbulece Research, Staford Uiversity, Staford, CA icoud@ ctr.staford.edu.
TMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods
TMA4205 Numerical Liear Algebra The Poisso problem i R 2 : diagoalizatio methods September 3, 2007 c Eiar M Røquist Departmet of Mathematical Scieces NTNU, N-749 Trodheim, Norway All rights reserved A
More informationPrinciple Of Superposition
ecture 5: PREIMINRY CONCEP O RUCUR NYI Priciple Of uperpositio Mathematically, the priciple of superpositio is stated as ( a ) G( a ) G( ) G a a or for a liear structural system, the respose at a give
More informationMachine Learning for Data Science (CS 4786)
Machie Learig for Data Sciece CS 4786) Lecture 9: Pricipal Compoet Aalysis The text i black outlies mai ideas to retai from the lecture. The text i blue give a deeper uderstadig of how we derive or get
More informationAppendix: The Laplace Transform
Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio
More informationAlgebra of Least Squares
October 19, 2018 Algebra of Least Squares Geometry of Least Squares Recall that out data is like a table [Y X] where Y collects observatios o the depedet variable Y ad X collects observatios o the k-dimesioal
More informationMachine Learning for Data Science (CS 4786)
Machie Learig for Data Sciece CS 4786) Lecture & 3: Pricipal Compoet Aalysis The text i black outlies high level ideas. The text i blue provides simple mathematical details to derive or get to the algorithm
More informationLecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods
Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace
More informationQuestion 1: The magnetic case
September 6, 018 Corell Uiversity, Departmet of Physics PHYS 337, Advace E&M, HW # 4, due: 9/19/018, 11:15 AM Questio 1: The magetic case I class, we skipped over some details, so here you are asked to
More informationSimilarity between quantum mechanics and thermodynamics: Entropy, temperature, and Carnot cycle
Similarity betwee quatum mechaics ad thermodyamics: Etropy, temperature, ad Carot cycle Sumiyoshi Abe 1,,3 ad Shiji Okuyama 1 1 Departmet of Physical Egieerig, Mie Uiversity, Mie 514-8507, Japa Istitut
More information6. Kalman filter implementation for linear algebraic equations. Karhunen-Loeve decomposition
6. Kalma filter implemetatio for liear algebraic equatios. Karhue-Loeve decompositio 6.1. Solvable liear algebraic systems. Probabilistic iterpretatio. Let A be a quadratic matrix (ot obligatory osigular.
More informationTrue Nature of Potential Energy of a Hydrogen Atom
True Nature of Potetial Eergy of a Hydroge Atom Koshu Suto Key words: Bohr Radius, Potetial Eergy, Rest Mass Eergy, Classical Electro Radius PACS codes: 365Sq, 365-w, 33+p Abstract I cosiderig the potetial
More informationFirst, note that the LS residuals are orthogonal to the regressors. X Xb X y = 0 ( normal equations ; (k 1) ) So,
0 2. OLS Part II The OLS residuals are orthogoal to the regressors. If the model icludes a itercept, the orthogoality of the residuals ad regressors gives rise to three results, which have limited practical
More informationFrequency Response of FIR Filters
EEL335: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we itroduce the idea of the frequecy respose of LTI systems, ad focus specifically o the frequecy respose of FIR filters.. Steady-state
More informationGeometry of LS. LECTURE 3 GEOMETRY OF LS, PROPERTIES OF σ 2, PARTITIONED REGRESSION, GOODNESS OF FIT
OCTOBER 7, 2016 LECTURE 3 GEOMETRY OF LS, PROPERTIES OF σ 2, PARTITIONED REGRESSION, GOODNESS OF FIT Geometry of LS We ca thik of y ad the colums of X as members of the -dimesioal Euclidea space R Oe ca
More informationA STUDY ON MHD BOUNDARY LAYER FLOW OVER A NONLINEAR STRETCHING SHEET USING IMPLICIT FINITE DIFFERENCE METHOD
IRET: Iteratioal oural of Research i Egieerig ad Techology eissn: 39-63 pissn: 3-7308 A STUDY ON MHD BOUNDARY LAYER FLOW OVER A NONLINEAR STRETCHING SHEET USING IMPLICIT FINITE DIFFERENCE METHOD Satish
More informationPhysics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.
Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;
More informationCO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS
CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS C.PRAX ad H.SADAT Laboratoire d'etudes Thermiques,URA CNRS 403 40, Aveue du Recteur Pieau 86022 Poitiers Cedex,
More informationSimilarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall
Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig
More informationEstimation of Backward Perturbation Bounds For Linear Least Squares Problem
dvaced Sciece ad Techology Letters Vol.53 (ITS 4), pp.47-476 http://dx.doi.org/.457/astl.4.53.96 Estimatio of Bacward Perturbatio Bouds For Liear Least Squares Problem Xixiu Li School of Natural Scieces,
More informationMatsubara-Green s Functions
Matsubara-Gree s Fuctios Time Orderig : Cosider the followig operator If H = H the we ca trivially factorise this as, E(s = e s(h+ E(s = e sh e s I geeral this is ot true. However for practical applicatio
More informationMath 257: Finite difference methods
Math 257: Fiite differece methods 1 Fiite Differeces Remember the defiitio of a derivative f f(x + ) f(x) (x) = lim 0 Also recall Taylor s formula: (1) f(x + ) = f(x) + f (x) + 2 f (x) + 3 f (3) (x) +...
More informationPAPER : IIT-JAM 2010
MATHEMATICS-MA (CODE A) Q.-Q.5: Oly oe optio is correct for each questio. Each questio carries (+6) marks for correct aswer ad ( ) marks for icorrect aswer.. Which of the followig coditios does NOT esure
More informationFluid Physics 8.292J/12.330J % (1)
Fluid Physics 89J/133J Problem Set 5 Solutios 1 Cosider the flow of a Euler fluid i the x directio give by for y > d U = U y 1 d for y d U + y 1 d for y < This flow does ot vary i x or i z Determie the
More informationFilter banks. Separately, the lowpass and highpass filters are not invertible. removes the highest frequency 1/ 2and
Filter bas Separately, the lowpass ad highpass filters are ot ivertible T removes the highest frequecy / ad removes the lowest frequecy Together these filters separate the sigal ito low-frequecy ad high-frequecy
More informationDeterminants of order 2 and 3 were defined in Chapter 2 by the formulae (5.1)
5. Determiats 5.. Itroductio 5.2. Motivatio for the Choice of Axioms for a Determiat Fuctios 5.3. A Set of Axioms for a Determiat Fuctio 5.4. The Determiat of a Diagoal Matrix 5.5. The Determiat of a Upper
More informationAE/ME 339 Computational Fluid Dynamics (CFD)
AE/ME 339 Computatioal Fluid Dyamics (CFD 0//004 Topic0_PresCorr_ Computatioal Fluid Dyamics (AE/ME 339 Pressure Correctio Method The pressure correctio formula (6.8.4 Calculatio of p. Coservatio form
More informationAssignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.
PHYSICS 34 QUANTUM PHYSICS II (25) Assigmet 2 Solutios 1. With respect to a pair of orthoormal vectors ϕ 1 ad ϕ 2 that spa the Hilbert space H of a certai system, the operator  is defied by its actio
More informationPAijpam.eu ON TENSOR PRODUCT DECOMPOSITION
Iteratioal Joural of Pure ad Applied Mathematics Volume 103 No 3 2015, 537-545 ISSN: 1311-8080 (prited versio); ISSN: 1314-3395 (o-lie versio) url: http://wwwijpameu doi: http://dxdoiorg/1012732/ijpamv103i314
More information5.1 Review of Singular Value Decomposition (SVD)
MGMT 69000: Topics i High-dimesioal Data Aalysis Falll 06 Lecture 5: Spectral Clusterig: Overview (cotd) ad Aalysis Lecturer: Jiamig Xu Scribe: Adarsh Barik, Taotao He, September 3, 06 Outlie Review of
More informationPOWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS
Joural of Applied Mathematics ad Computatioal Mechaics 4 3(3) 3-8 POWER SERIES SOLUION OF FIRS ORDER MARIX DIFFERENIAL EQUAIONS Staisław Kukla Izabela Zamorska Istitute of Mathematics Czestochowa Uiversity
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationHydro code II. Simone Recchi
Hydro code II Simoe Recchi 16.06.2011 Literature E.F. Toro: Riema Solvers ad Numerical Methods for Fluid Dyamics, Spriger R.J. LeVeque: Fiite Volume Methods for Hyperbolic Problems, Cambridge Uiversity
More informationSequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018
CSE 353 Discrete Computatioal Structures Sprig 08 Sequeces, Mathematical Iductio, ad Recursio (Chapter 5, Epp) Note: some course slides adopted from publisher-provided material Overview May mathematical
More informationSection 5.5. Infinite Series: The Ratio Test
Differece Equatios to Differetial Equatios Sectio 5.5 Ifiite Series: The Ratio Test I the last sectio we saw that we could demostrate the covergece of a series a, where a 0 for all, by showig that a approaches
More informationTaylor expansion: Show that the TE of f(x)= sin(x) around. sin(x) = x - + 3! 5! L 7 & 8: MHD/ZAH
Taylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. A ay poit i the eighbourhood of 0, the fuctio ƒ() ca be represeted by a power series of the followig form: X 0 f(a) f() f() ( ) f( ) ( )
More informationRay Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET
Ray Optics Theory ad Mode Theory Dr. Mohammad Faisal Dept. of, BUT Optical Fiber WG For light to be trasmitted through fiber core, i.e., for total iteral reflectio i medium, > Ray Theory Trasmissio Ray
More informationSequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece 1, 1, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet
More informationL 5 & 6: RelHydro/Basel. f(x)= ( ) f( ) ( ) ( ) ( ) n! 1! 2! 3! If the TE of f(x)= sin(x) around x 0 is: sin(x) = x - 3! 5!
aylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. At ay poit i the eighbourhood of =0, the fuctio ca be represeted as a power series of the followig form: X 0 f(a) f() ƒ() f()= ( ) f( ) (
More informationModified Decomposition Method by Adomian and. Rach for Solving Nonlinear Volterra Integro- Differential Equations
Noliear Aalysis ad Differetial Equatios, Vol. 5, 27, o. 4, 57-7 HIKARI Ltd, www.m-hikari.com https://doi.org/.2988/ade.27.62 Modified Decompositio Method by Adomia ad Rach for Solvig Noliear Volterra Itegro-
More informationIP Reference guide for integer programming formulations.
IP Referece guide for iteger programmig formulatios. by James B. Orli for 15.053 ad 15.058 This documet is iteded as a compact (or relatively compact) guide to the formulatio of iteger programs. For more
More informationDiscrete Mathematics for CS Spring 2008 David Wagner Note 22
CS 70 Discrete Mathematics for CS Sprig 2008 David Wager Note 22 I.I.D. Radom Variables Estimatig the bias of a coi Questio: We wat to estimate the proportio p of Democrats i the US populatio, by takig
More informationNotes for Lecture 11
U.C. Berkeley CS78: Computatioal Complexity Hadout N Professor Luca Trevisa 3/4/008 Notes for Lecture Eigevalues, Expasio, ad Radom Walks As usual by ow, let G = (V, E) be a udirected d-regular graph with
More information1 Inferential Methods for Correlation and Regression Analysis
1 Iferetial Methods for Correlatio ad Regressio Aalysis I the chapter o Correlatio ad Regressio Aalysis tools for describig bivariate cotiuous data were itroduced. The sample Pearso Correlatio Coefficiet
More informationSalmon: Lectures on partial differential equations. 3. First-order linear equations as the limiting case of second-order equations
3. First-order liear equatios as the limitig case of secod-order equatios We cosider the advectio-diffusio equatio (1) v = 2 o a bouded domai, with boudary coditios of prescribed. The coefficiets ( ) (2)
More informationR is a scalar defined as follows:
Math 8. Notes o Dot Product, Cross Product, Plaes, Area, ad Volumes This lecture focuses primarily o the dot product ad its may applicatios, especially i the measuremet of agles ad scalar projectio ad
More informationNumerical Methods for Ordinary Differential Equations
Numerical Methods for Ordiary Differetial Equatios Braislav K. Nikolić Departmet of Physics ad Astroomy, Uiversity of Delaware, U.S.A. PHYS 460/660: Computatioal Methods of Physics http://www.physics.udel.edu/~bikolic/teachig/phys660/phys660.html
More informationMeasurement uncertainty of the sound absorption
Measuremet ucertaity of the soud absorptio coefficiet Aa Izewska Buildig Research Istitute, Filtrowa Str., 00-6 Warsaw, Polad a.izewska@itb.pl 6887 The stadard ISO/IEC 705:005 o the competece of testig
More informationStreamfunction-Vorticity Formulation
Streamfuctio-Vorticity Formulatio A. Salih Departmet of Aerospace Egieerig Idia Istitute of Space Sciece ad Techology, Thiruvaathapuram March 2013 The streamfuctio-vorticity formulatio was amog the first
More informationLyapunov Stability Analysis for Feedback Control Design
Copyright F.L. Lewis 008 All rights reserved Updated: uesday, November, 008 Lyapuov Stability Aalysis for Feedbac Cotrol Desig Lyapuov heorems Lyapuov Aalysis allows oe to aalyze the stability of cotiuous-time
More informationChapter 7: The z-transform. Chih-Wei Liu
Chapter 7: The -Trasform Chih-Wei Liu Outlie Itroductio The -Trasform Properties of the Regio of Covergece Properties of the -Trasform Iversio of the -Trasform The Trasfer Fuctio Causality ad Stability
More informationExplicit Group Methods in the Solution of the 2-D Convection-Diffusion Equations
Proceedigs of the World Cogress o Egieerig 00 Vol III WCE 00 Jue 0 - July 00 Lodo U.K. Explicit Group Methods i the Solutio of the -D Covectio-Diffusio Equatios a Kah Bee orhashidah Hj. M. Ali ad Choi-Hog
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2016 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More information1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.
.3 Covergece Theorems of Fourier Series I this sectio, we preset the covergece of Fourier series. A ifiite sum is, by defiitio, a limit of partial sums, that is, a cos( kx) b si( kx) lim a cos( kx) b si(
More informationSingular value decomposition. Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine
Lecture 11 Sigular value decompositio Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaie V1.2 07/12/2018 1 Sigular value decompositio (SVD) at a glace Motivatio: the image of the uit sphere S
More informationInfinite Sequences and Series
Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet
More informationFFTs in Graphics and Vision. The Fast Fourier Transform
FFTs i Graphics ad Visio The Fast Fourier Trasform 1 Outlie The FFT Algorithm Applicatios i 1D Multi-Dimesioal FFTs More Applicatios Real FFTs 2 Computatioal Complexity To compute the movig dot-product
More informationOn forward improvement iteration for stopping problems
O forward improvemet iteratio for stoppig problems Mathematical Istitute, Uiversity of Kiel, Ludewig-Mey-Str. 4, D-24098 Kiel, Germay irle@math.ui-iel.de Albrecht Irle Abstract. We cosider the optimal
More informationA sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece,, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet as
More informationLesson 10: Limits and Continuity
www.scimsacademy.com Lesso 10: Limits ad Cotiuity SCIMS Academy 1 Limit of a fuctio The cocept of limit of a fuctio is cetral to all other cocepts i calculus (like cotiuity, derivative, defiite itegrals
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 9 Multicolliearity Dr Shalabh Departmet of Mathematics ad Statistics Idia Istitute of Techology Kapur Multicolliearity diagostics A importat questio that
More informationTHE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.
THE SOLUTION OF NONLINEAR EQUATIONS f( ) = 0. Noliear Equatio Solvers Bracketig. Graphical. Aalytical Ope Methods Bisectio False Positio (Regula-Falsi) Fied poit iteratio Newto Raphso Secat The root of
More informationA proposed discrete distribution for the statistical modeling of
It. Statistical Ist.: Proc. 58th World Statistical Cogress, 0, Dubli (Sessio CPS047) p.5059 A proposed discrete distributio for the statistical modelig of Likert data Kidd, Marti Cetre for Statistical
More informationBoundary layer problem on conveyor belt. Gabriella Bognár University of Miskolc 3515 Miskolc-Egyetemváros, Hungary
Boudary layer problem o coveyor belt Gabriella Bogár Uiversity of Miskolc 355 Miskolc-Egyetemváros, Hugary e-mail: matvbg@ui-miskolc.hu Abstract: A techologically importat source of the boudary layer pheomeo
More informationMETHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS
Please cite this article as: Staisław Kula, Method of fudametal solutios for Helmholtz eigevalue problems i elliptical domais, Scietific Research of the Istitute of Mathematics ad Computer Sciece, 009,
More informationDiscrete-Time Systems, LTI Systems, and Discrete-Time Convolution
EEL5: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we begi our mathematical treatmet of discrete-time s. As show i Figure, a discrete-time operates or trasforms some iput sequece x [
More informationSpecial Modeling Techniques
Colorado School of Mies CHEN43 Secial Modelig Techiques Secial Modelig Techiques Summary of Toics Deviatio Variables No-Liear Differetial Equatios 3 Liearizatio of ODEs for Aroximate Solutios 4 Coversio
More informationPSYCHOLOGICAL RESEARCH (PYC 304-C) Lecture 9
Hypothesis testig PSYCHOLOGICAL RESEARCH (PYC 34-C Lecture 9 Statistical iferece is that brach of Statistics i which oe typically makes a statemet about a populatio based upo the results of a sample. I
More informationAxioms of Measure Theory
MATH 532 Axioms of Measure Theory Dr. Neal, WKU I. The Space Throughout the course, we shall let X deote a geeric o-empty set. I geeral, we shall ot assume that ay algebraic structure exists o X so that
More informationNumerical Astrophysics: hydrodynamics
Numerical Astrophysics: hydrodyamics Part 1: Numerical solutios to the Euler Equatios Outlie Numerical eperimets are a valuable tool to study astrophysical objects (where we ca rarely do direct eperimets).
More information1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y
Questio (a) A square matrix A= A is called positive defiite if the quadratic form waw > 0 for every o-zero vector w [Note: Here (.) deotes the traspose of a matrix or a vector]. Let 0 A = 0 = show that:
More informationOn Nonsingularity of Saddle Point Matrices. with Vectors of Ones
Iteratioal Joural of Algebra, Vol. 2, 2008, o. 4, 197-204 O Nosigularity of Saddle Poit Matrices with Vectors of Oes Tadeusz Ostrowski Istitute of Maagemet The State Vocatioal Uiversity -400 Gorzów, Polad
More informationStability of fractional positive nonlinear systems
Archives of Cotrol Scieces Volume 5(LXI), 15 No. 4, pages 491 496 Stability of fractioal positive oliear systems TADEUSZ KACZOREK The coditios for positivity ad stability of a class of fractioal oliear
More informationNICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =
AN INTRODUCTION TO SCHRÖDER AND UNKNOWN NUMBERS NICK DUFRESNE Abstract. I this article we will itroduce two types of lattice paths, Schröder paths ad Ukow paths. We will examie differet properties of each,
More informationA CONFINEMENT MODEL OF HIGH STRENGTH CONCRETE
3 th World Coferece o Earthquake Egieerig Vacouver, B.C., Caada August -6, 24 Paper No. 873 A CONFINEMENT MODEL OF HIGH STRENGTH CONCRETE Nobutaka NAKAZAWA, Kazuhiko KAWASHIMA 2, Gakuho WATANABE 3, Ju-ichi
More informationDiscrete Orthogonal Moment Features Using Chebyshev Polynomials
Discrete Orthogoal Momet Features Usig Chebyshev Polyomials R. Mukuda, 1 S.H.Og ad P.A. Lee 3 1 Faculty of Iformatio Sciece ad Techology, Multimedia Uiversity 75450 Malacca, Malaysia. Istitute of Mathematical
More informationIntroduction to Machine Learning DIS10
CS 189 Fall 017 Itroductio to Machie Learig DIS10 1 Fu with Lagrage Multipliers (a) Miimize the fuctio such that f (x,y) = x + y x + y = 3. Solutio: The Lagragia is: L(x,y,λ) = x + y + λ(x + y 3) Takig
More informationDecoupling Zeros of Positive Discrete-Time Linear Systems*
Circuits ad Systems,,, 4-48 doi:.436/cs..7 Published Olie October (http://www.scirp.org/oural/cs) Decouplig Zeros of Positive Discrete-Time Liear Systems* bstract Tadeusz Kaczorek Faculty of Electrical
More informationTaylor polynomial solution of difference equation with constant coefficients via time scales calculus
TMSCI 3, o 3, 129-135 (2015) 129 ew Treds i Mathematical Scieces http://wwwtmscicom Taylor polyomial solutio of differece equatio with costat coefficiets via time scales calculus Veysel Fuat Hatipoglu
More informationChapter 1. Complex Numbers. Dr. Pulak Sahoo
Chapter 1 Complex Numbers BY Dr. Pulak Sahoo Assistat Professor Departmet of Mathematics Uiversity Of Kalyai West Begal, Idia E-mail : sahoopulak1@gmail.com 1 Module-2: Stereographic Projectio 1 Euler
More informationMicroscopic traffic flow modeling
Chapter 34 Microscopic traffic flow modelig 34.1 Overview Macroscopic modelig looks at traffic flow from a global perspective, whereas microscopic modelig, as the term suggests, gives attetio to the details
More informationA note on the modified Hermitian and skew-hermitian splitting methods for non-hermitian positive definite linear systems
A ote o the modified Hermitia ad skew-hermitia splittig methods for o-hermitia positive defiite liear systems Shi-Liag Wu Tig-Zhu Huag School of Applied Mathematics Uiversity of Electroic Sciece ad Techology
More informationQuantum Annealing for Heisenberg Spin Chains
LA-UR # - Quatum Aealig for Heiseberg Spi Chais G.P. Berma, V.N. Gorshkov,, ad V.I.Tsifriovich Theoretical Divisio, Los Alamos Natioal Laboratory, Los Alamos, NM Istitute of Physics, Natioal Academy of
More informationA Block Cipher Using Linear Congruences
Joural of Computer Sciece 3 (7): 556-560, 2007 ISSN 1549-3636 2007 Sciece Publicatios A Block Cipher Usig Liear Cogrueces 1 V.U.K. Sastry ad 2 V. Jaaki 1 Academic Affairs, Sreeidhi Istitute of Sciece &
More information(b) What is the probability that a particle reaches the upper boundary n before the lower boundary m?
MATH 529 The Boudary Problem The drukard s walk (or boudary problem) is oe of the most famous problems i the theory of radom walks. Oe versio of the problem is described as follows: Suppose a particle
More informationAn Introduction to Randomized Algorithms
A Itroductio to Radomized Algorithms The focus of this lecture is to study a radomized algorithm for quick sort, aalyze it usig probabilistic recurrece relatios, ad also provide more geeral tools for aalysis
More informationOptics. n n. sin. 1. law of rectilinear propagation 2. law of reflection = 3. law of refraction
Optics What is light? Visible electromagetic radiatio Geometrical optics (model) Light-ray: extremely thi parallel light beam Usig this model, the explaatio of several optical pheomea ca be give as the
More informationCONTENTS. Course Goals. Course Materials Lecture Notes:
INTRODUCTION Ho Chi Mih City OF Uiversity ENVIRONMENTAL of Techology DESIGN Faculty Chapter of Civil 1: Orietatio. Egieerig Evaluatio Departmet of mathematical of Water Resources skill Egieerig & Maagemet
More informationNumerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION
alaysia Joural of athematical Scieces 3(1): 83-93 (9) umerical Coformal appig via a Fredholm Itegral Equatio usig Fourier ethod 1 Ali Hassa ohamed urid ad Teh Yua Yig 1, Departmet of athematics, Faculty
More informationLainiotis filter implementation. via Chandrasekhar type algorithm
Joural of Computatios & Modellig, vol.1, o.1, 2011, 115-130 ISSN: 1792-7625 prit, 1792-8850 olie Iteratioal Scietific Press, 2011 Laiiotis filter implemetatio via Chadrasehar type algorithm Nicholas Assimais
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationSolutions to Homework 7
Solutios to Homework 7 Due Wedesday, August 4, 004. Chapter 4.1) 3, 4, 9, 0, 7, 30. Chapter 4.) 4, 9, 10, 11, 1. Chapter 4.1. Solutio to problem 3. The sum has the form a 1 a + a 3 with a k = 1/k. Sice
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationThe z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j
The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.
More informationTHE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS
R775 Philips Res. Repts 26,414-423, 1971' THE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS by H. W. HANNEMAN Abstract Usig the law of propagatio of errors, approximated
More informationSimple Linear Regression
Chapter 2 Simple Liear Regressio 2.1 Simple liear model The simple liear regressio model shows how oe kow depedet variable is determied by a sigle explaatory variable (regressor). Is is writte as: Y i
More informationAnalysis of a Numerical Scheme An Example
http://www.d.edu/~gtryggva/cfd-course/ Computatioal Fluid Dyamics Lecture 3 Jauary 5, 7 Aalysis of a Numerical Scheme A Example Grétar Tryggvaso Numerical Aalysis Example Use the leap-frog method (cetered
More informationCALCULATION OF FIBONACCI VECTORS
CALCULATION OF FIBONACCI VECTORS Stuart D. Aderso Departmet of Physics, Ithaca College 953 Daby Road, Ithaca NY 14850, USA email: saderso@ithaca.edu ad Dai Novak Departmet of Mathematics, Ithaca College
More information2D DSP Basics: 2D Systems
- Digital Image Processig ad Compressio D DSP Basics: D Systems D Systems T[ ] y = T [ ] Liearity Additivity: If T y = T [ ] The + T y = y + y Homogeeity: If The T y = T [ ] a T y = ay = at [ ] Liearity
More informationLecture III-2: Light propagation in nonmagnetic
A. La Rosa Lecture Notes ALIED OTIC Lecture III2: Light propagatio i omagetic materials 2.1 urface ( ), volume ( ), ad curret ( j ) desities produced by arizatio charges The objective i this sectio is
More informationDETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.
ICSV4 Cairs Australia 9- July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON- UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics
More information