One Dimensional Convection: Interpolation Models for CFD
|
|
- Augustus Bradford
- 6 years ago
- Views:
Transcription
1 One Dimensional Convection: Interpolation Moels for CFD ME 448/548 Notes Geral Recktenwal Portlan State University Department of Mechanical Engineering ME 448/548: D Convection-Di usion Equation Overview Demonstrate the problem arising from central i erence approimation of the convection term.. Introuce upwin concept an provie a Matlab implementation. Introuce non-uniform mesh behavior, an show that the local Peclet number matters. Demonstrate numerical behavior with Matlab computations. ME 448/548: D Convection-Di usion Equation page
2 Moel Problem The one-imensional convection-i usion equation is (u ) S =0 0apple apple L () For S =0an the bounary conitions the eact solution to Equation () is (0) = 0 (L) = L, (2) 0 L 0 = ep(u/ ) ep(pe L ) (3) where Pe L = ul (4) ME 448/548: D Convection-Di usion Equation page 2 Eact Solution to the Moel Problem The parameter is Pe L = ul φ ME 448/548: D Convection-Di usion Equation page 3
3 Finite Volume Mesh Two-imensional mesh δ w δ Δy W P E i- i i+ y w Δ One-imensional mesh, with emphasis on cell sizes near the bounary. Δ 2 i = Δ m m 2 m m δ w,2 δ,2 δ w,m δ,m ME 448/548: D Convection-Di usion Equation page 4 Integrate the Moel Equation over a Control Volume Obtain a numerical moel of the convection-i usion equation (u ) S =0 by integrating with respect to over the control volume. δ w δ Z w (u ) Z w Z w S=0 (5) y Δy i- W w i P Δ E i+ ME 448/548: D Convection-Di usion Equation page 5
4 Integrate the Di usion Term Integrate i usion term once Z w (6) = e w y Δy i- W δ w w i P Δ δ E i+ = flu of across the east face, out of the C.V. e = flu of across the west face, into the C.V. w ME 448/548: D Convection-Di usion Equation page 6 Integrate the Di usion Term Replace the i usion flues with finite i erence approimations an e w D e = e E P = D e ( E P ) w P W w = D w ( P W ) e D w = w (7) w y Δy i- W δ w w i P Δ δ E i+ = E P w = P W. (8) W, P,an E are values at the noes, i.e.,thesearetheunknowns(theiscrete epenent variables) in the numerical moel. ME 448/548: D Convection-Di usion Equation page 7
5 Integrate the Di usion Term Assume that is uniform: e = w =. Z w = e w (9) D e ( E P ) D w ( P W ) (0) ME 448/548: D Convection-Di usion Equation page 8 Integrate the Source Term Integrate the source term, assuming that S is uniform over the control volume. Z w S S P p. () y Δy i- W δ w w i P Δ δ E i+ ME 448/548: D Convection-Di usion Equation page 9
6 Convection Term Overview. Integrate once: What is the value of at the interface? 2. Linear interpolation leas to central i erence 3. Central i erence causes unacceptable wiggles on a coarse mesh 4. Upwin interpolation suppresses the wiggles, but reuces the accuracy ME 448/548: D Convection-Di usion Equation page 0 Convection Term First Try The first step is easy Z w (u ) =(u ) e (u ) w (2) How o we efine (u ) e an (u ) w? (u ) e = prouct of u an at the east interface (u ) w = prouct of u an at the west interface ME 448/548: D Convection-Di usion Equation page
7 Use Linear Interpolation at the Interface Use linear interpolation to estimate e from P an E, antoestimate w from W an P where e = e E + ( e) P (3) w = w W + ( w) P (4) P e = an w = P w. E P P W ϕ W ϕ P W w P E ϕe ME 448/548: D Convection-Di usion Equation page 2 Reconstitute the Convection Terms Substituting Equation (3) an Equation (4) into Equation (2) an rearranging gives Z w (u ) = u e e( E P ) u w w ( W P )+u e P u w P (5) The last two terms in the preceing equation cancel because u is a uniform parameter, i.e. u e = u w. Z w (u ) = u e e( E P ) u w w ( W P ) (6) ME 448/548: D Convection-Di usion Equation page 3
8 Final Form of the Discrete Convection Di usion Equation Substituting Equation (0), Equation () an Equation (6) into Equation (5) an simplifying gives a E E + a P P a W W = b (7) where a E = P (D e u e e ) (8) a W = P (D w + u w w ) (9) a P = a E + a W (20) b = S P (2) ME 448/548: D Convection-Di usion Equation page 4 Another Triiagonal System of Equations a P, a E, a W,2 a P,2... a E, a W,i a P,i a E,i i. 3 2 = b b 2. b i (22) a W,m a P,m m b m Coe cients in the first an last equations are ajuste to enforce the bounary conitions ( = 0) = 0 =) a P, =, a E, =0, b = 0 ( = L) = L =) a P,m =, a W,m =0, b m = L ME 448/548: D Convection-Di usion Equation page 5
9 Numerical Solution with the Central Di erence Scheme Central i erence solutions to Equation () for Pe L = 50, Pe =5. This is the steay solution. The oscillations are not an instability: the solution oes not grow as the result of iterations or time steps..5 φ 0.5 CDS solution eact PeL = 50.0, Pe = 5.0 CDS scheme, Marror = ME 448/548: D Convection-Di usion Equation page 6 Numerical Solution with the Central Di erence Scheme Central i erence solutions to Equation () for Pe L = 50, Pe =.7..5 CDS solution eact Increasing the number of noes has eliminate the wil oscillations in the numerical solution. φ 0.5 PeL = 50.0, Pe =.7 CDS scheme, Marror = ME 448/548: D Convection-Di usion Equation page 7
10 Consequences of Negative Coe cients Consier the possible values taken by a E where the interpolation coe cient e is Rearrange Equation (8) as a E = P (D e u e e ) (8) e = P /2 a E = P (D e u e e )= P = P e e P /2 u e Pe 2 ME 448/548: D Convection-Di usion Equation page 8 Consequences of Negative Coe cients We have a E = P e Pe 2 Since P e > 0 always, we observe the a E becomes negative when or Pe < 0 2 Pe > 2 ME 448/548: D Convection-Di usion Equation page 9
11 Consequences of Negative Coe cients When Pe > 2 the value of a E becomes negative an the eigenvalues of the coe an the numerical solution oscillates. cients become comple Note that we can keep Pe = u e e P < 2 by making P small enough. That is why reucing the mesh size (increasing the number of noes) eliminate the error in the sample calculation. ME 448/548: D Convection-Di usion Equation page 20 Upwin Di erencing: A cure with a cost. Negative a E can be avoie with juicious choice of e. Recallthat 0.9 an e = e E + ( e) P 0.8 φ Increasing u a E = P (D e u e e ) Physical reasoning suggests that e shoul be influence by u e. 3. Take inspiration from the analytical solution to Equation () As the velocity increases, the value of at the mipoint gets closer to the value of the upstream bounary. The Upwin Di erence Scheme gives a strong bias to the upstream value of. ME 448/548: D Convection-Di usion Equation page 2
12 Upwin Di erencing: A cure with a cost Strategy: Ignore the value of at the ownwin noe: ϕ W u w > 0 u e > 0 u w < 0 u e < 0 ϕ W ϕ P ϕ E ϕ P ϕ E W w P E W w P E u w > 0, u e > 0: u w < 0, u e < 0: w = e =0 w = W e = P w =0 e = w = P e = E ME 448/548: D Convection-Di usion Equation page 22 Upwin Di erencing: A cure with a cost ϕ W u w > 0 u e > 0 u w < 0 u e < 0 ϕ W ϕ P ϕ E ϕ P ϕ E W w P E W w P E Restate the rules e = ( 0 if ue 0 if u e < 0 w = ( if uw 0 0 if u w < 0 ME 448/548: D Convection-Di usion Equation page 23
13 Compare Upwin Solution with Central Di erence Solution Compare UDS an CDS for m = 2, Pe =5. No oscillations for the UDS scheme 2.5 Central Upwin Eact φ 0.5 Pe L = 50.0, Pe = e c =.583, e u = ME 448/548: D Convection-Di usion Equation page 24 Compare Upwin Solution with Central Di erence Solution Compare UDS an CDS for m = 32, Pe =.7. No oscillations for either scheme. 2.5 Central Upwin Eact φ 0.5 Pe L = 50.0, Pe =.7 0 e c = 0.268, e u = ME 448/548: D Convection-Di usion Equation page 25
14 Convergence stuy for Upwin an Central Di erence Schemes Solve the moel problem on a series of finer meshes Upwin Central Di erence m Pe e u e u ratio e c e c ratio ME 448/548: D Convection-Di usion Equation page 26 Convergence stuy for Upwin an Central Di erence Schemes m e u ratio e c ratio Marror Central Upwin ME 448/548: D Convection-Di usion Equation page 27
15 Improving Accuracy with a Non-uniform Mesh A non-uniform mesh can improve accuracy of a numerical solution. Noes shoul be close together in regions where the graient is higher. Demonstrate the iea with a stretche mesh L Δ Δ 2 Δ 3 Δ n ME 448/548: D Convection-Di usion Equation page 28 Stretche Mesh () A growing or shrinking cell size, an hence noe spacing, can be achieve by specifying where r is a fie constant. i+ i = r (23) r> causes noe spacing to increase as i increases. r< causes noe spacing to ecrease as i increases. L Δ Δ 2 Δ 3 Δ n ME 448/548: D Convection-Di usion Equation page 29
16 Stretche Mesh (2). The total length L of the region to be subivie is given. 2. Choose the stretching ratio r, an the number of control volumes n. Constraint: The withs of the control volumes must a up to L. L = nx i= i = + r + r r n = +r + r r n The sum of terms can be rewritten. S =+r + r r n = rn r (24) Thus, L = S or = L S. (25) ME 448/548: D Convection-Di usion Equation page 30 Stretche Mesh (3) To create a mesh with a geometric progression of sizes. Specify L, r, ann. 2. Compute from Equation (25). 3. Compute the remaining cell withs with a loop. Note that the values of r an n interact in etermining the actual spacing. It is possible to get istorte meshes by changing r with n hel constant, or changing n with r hel constant. ME 448/548: D Convection-Di usion Equation page 3
17 Stretche Mesh (3) Eamples: Use L =, n = 0 an vary r r =.05 n = 0 r = 0.95 n = 0 r =.0 n = 0 r = 0.90 n = 0 r =.20 n = 0 r = 0.80 n = 0 r =.30 n = 0 r = 0.70 n = 0 r =.40 n = 0 r = 0.60 n = 0 ME 448/548: D Convection-Di usion Equation page 32 Stretche Mesh (4) Eamples: Use L =, an vary r r =.05 r = 0.95 r =.0 r = 0.90 r =.20 r = 0.80 r =.30 r = 0.70 r =.40 r = 0.60 For large n, avoi etreme istortion by using r. ME 448/548: D Convection-Di usion Equation page 33
18 Stretche Mesh (4) Eamples: Use L =, r =0.90, r =0.95 an vary n r = 0.90 n = 0 r = 0.95 n = 0 r = 0.90 r = 0.95 r = 0.90 n = 30 r = 0.95 n = 30 r = 0.90 n = 40 r = 0.95 n = 40 For large n, avoi etreme istortion by using r. ME 448/548: D Convection-Di usion Equation page 34 Stretche Solution 2.5 Central Upwin Eact n = 24, r =0.90 n = 24, r = Central Upwin Eact φ φ 0.5 Pe L = 50.0, Pe = Pe L = 50.0, Pe = e c = 0.042, e u = e c = 0.007, e u = Observations: Careful choice of stretching parameters can yiel goo accuracy on a moest mesh. No obvious way of preicting the accuracy when the eact solution is not available. ME 448/548: D Convection-Di usion Equation page 35
19 Summary () We showe that Conventional finite volume schemes store epenent variables at cell centers. Avection terms require a seconary scheme for interpolating epenent variables at cell faces. The central i erence scheme can cause oscillations The upwin i erence scheme suppresses oscillations, but introuces aitional truncation error. ME 448/548: D Convection-Di usion Equation page 36 Summary (2) We i not inclue the following etails No proofs were given for truncation error estimates. Many other schemes eist. QUICK (quaratic upstream interpolation for convective kinematics),. ENO (essentially non-oscillatory),. TVD (total variational iminishing),. Secon-orer upwining,. Streamline upwining,. an others. ME 448/548: D Convection-Di usion Equation page 37
20 Summary (3) Users nee to be aware of trae-o s Oscillatory higher orer schemes cause solutions to fail, i.e., not converge Venors may use a efault scheme that is stable, but inaccurate, instea of a more accurate scheme Users may choose a stable an inaccurate scheme just to get the coe to converge. That s not a goo iea, especially if the accuracy of the solution is important. Don t suppress the wiggles They re trying to tell you something! P. Gresho an R. Lee, 98 ME 448/548: D Convection-Di usion Equation page 38
One Dimensional Convection: Interpolation Models for CFD
One Dimensional Convection: Interpolation Models for CFD ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@pdx.edu ME 448/548: 1D Convection-Diffusion
More informationChapter 9 Method of Weighted Residuals
Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationSection 7.1: Integration by Parts
Section 7.1: Integration by Parts 1. Introuction to Integration Techniques Unlike ifferentiation where there are a large number of rules which allow you (in principle) to ifferentiate any function, the
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More informationTHE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS
THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS
More informationDissipative numerical methods for the Hunter-Saxton equation
Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a
More information1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.
Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency
More informationTMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments
Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary
More informationCalculus of Variations
Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,
More informationUnit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule
Unit # - Families of Functions, Taylor Polynomials, l Hopital s Rule Some problems an solutions selecte or aapte from Hughes-Hallett Calculus. Critical Points. Consier the function f) = 54 +. b) a) Fin
More informationLecture contents. Metal-semiconductor contact
1 Lecture contents Metal-semiconuctor contact Electrostatics: Full epletion approimation Electrostatics: Eact electrostatic solution Current Methos for barrier measurement Junctions: general approaches,
More informationSummary: Differentiation
Techniques of Differentiation. Inverse Trigonometric functions The basic formulas (available in MF5 are: Summary: Differentiation ( sin ( cos The basic formula can be generalize as follows: Note: ( sin
More informationProof of SPNs as Mixture of Trees
A Proof of SPNs as Mixture of Trees Theorem 1. If T is an inuce SPN from a complete an ecomposable SPN S, then T is a tree that is complete an ecomposable. Proof. Argue by contraiction that T is not a
More informationWeb-Based Technical Appendix: Multi-Product Firms and Trade Liberalization
Web-Base Technical Appeni: Multi-Prouct Firms an Trae Liberalization Anrew B. Bernar Tuck School of Business at Dartmouth & NBER Stephen J. Reing LSE, Yale School of Management & CEPR Peter K. Schott Yale
More informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationMath 1271 Solutions for Fall 2005 Final Exam
Math 7 Solutions for Fall 5 Final Eam ) Since the equation + y = e y cannot be rearrange algebraically in orer to write y as an eplicit function of, we must instea ifferentiate this relation implicitly
More informationDesigning Information Devices and Systems I Spring 2018 Lecture Notes Note 16
EECS 16A Designing Information Devices an Systems I Spring 218 Lecture Notes Note 16 16.1 Touchscreen Revisite We ve seen how a resistive touchscreen works by using the concept of voltage iviers. Essentially,
More informationThe Exact Form and General Integrating Factors
7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily
More informationSome New Thoughts on the Multipoint Method for Reactor Physics Applications. Sandra Dulla, Piero Ravetto, Paolo Saracco,
Jeju, Korea, April 16-20, 2017, on USB 2017 Some New Thoughts on the Multipoint Metho for Reactor Physics Applications Sanra Dulla, Piero Ravetto, Paolo Saracco, Politecnico i Torino, Dipartimento Energia,
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationFirst Order Linear Differential Equations
LECTURE 8 First Orer Linear Differential Equations We now turn our attention to the problem of constructing analytic solutions of ifferential equations; that is to say,solutions that can be epresse in
More information5.4 Fundamental Theorem of Calculus Calculus. Do you remember the Fundamental Theorem of Algebra? Just thought I'd ask
5.4 FUNDAMENTAL THEOREM OF CALCULUS Do you remember the Funamental Theorem of Algebra? Just thought I' ask The Funamental Theorem of Calculus has two parts. These two parts tie together the concept of
More informationSources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences
Sources an Sinks of Available Potential Energy in a Moist Atmosphere Olivier Pauluis 1 Courant Institute of Mathematical Sciences New York University Submitte to the Journal of the Atmospheric Sciences
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More information18 EVEN MORE CALCULUS
8 EVEN MORE CALCULUS Chapter 8 Even More Calculus Objectives After stuing this chapter you shoul be able to ifferentiate an integrate basic trigonometric functions; unerstan how to calculate rates of change;
More informationx = c of N if the limit of f (x) = L and the right-handed limit lim f ( x)
Limit We say the limit of f () as approaches c equals L an write, lim L. One-Sie Limits (Left an Right-Hane Limits) Suppose a function f is efine near but not necessarily at We say that f has a left-hane
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson
JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises
More informationII. First variation of functionals
II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationA Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation
A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationNonlinear Dielectric Response of Periodic Composite Materials
onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics
More informationCMSC 313 Preview Slides
CMSC 33 Preview Slies These are raft slies. The actual slies presente in lecture may be ifferent ue to last minute changes, scheule slippage,... UMBC, CMSC33, Richar Chang CMSC 33 Lecture
More informationLearning in Monopolies with Delayed Price Information
Learning in Monopolies with Delaye Price Information Akio Matsumoto y Chuo University Ferenc Sziarovszky z University of Pécs February 28, 2013 Abstract We call the intercept of the price function with
More informationIERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210
IERCU Institute of Economic Research, Chuo University 50th Anniversary Special Issues Discussion Paper No.210 Discrete an Continuous Dynamics in Nonlinear Monopolies Akio Matsumoto Chuo University Ferenc
More informationProblem 3.84 of Bergman. Consider one-dimensional conduction in a plane composite wall. The outer surfaces are exposed to a fluid at T
1/10 bergman3-84.xmc Problem 3.84 of Bergman. Consier one-imensional conuction in a plane composite wall. The outer surfaces are expose to a flui at T 5 C an a convection heat transfer coefficient of h1000
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationLecture 6: Calculus. In Song Kim. September 7, 2011
Lecture 6: Calculus In Song Kim September 7, 20 Introuction to Differential Calculus In our previous lecture we came up with several ways to analyze functions. We saw previously that the slope of a linear
More informationCode_Aster. Detection of the singularities and computation of a card of size of elements
Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : Josselin DLMAS Clé : R4.0.04 Révision : 9755 Detection of the singularities an computation of a car of size
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationA. Incorrect! The letter t does not appear in the expression of the given integral
AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question
More informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More informationECE 422 Power System Operations & Planning 7 Transient Stability
ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3 Transient Stability The ability of
More informationLecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.
b. Partial erivatives Lecture b Differential operators an orthogonal coorinates Recall from our calculus courses that the erivative of a function can be efine as f ()=lim 0 or using the central ifference
More informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More information3.2 Differentiability
Section 3 Differentiability 09 3 Differentiability What you will learn about How f (a) Might Fail to Eist Differentiability Implies Local Linearity Numerical Derivatives on a Calculator Differentiability
More informationx f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0.
Engineering Mathematics 2 26 February 2014 Limits of functions Consier the function 1 f() = 1. The omain of this function is R + \ {1}. The function is not efine at 1. What happens when is close to 1?
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationMonopoly Part III: Multiple Local Equilibria and Adaptive Search
FH-Kiel University of Applie Sciences Prof Dr Anreas Thiemer, 00 e-mail: anreasthiemer@fh-kiele Monopoly Part III: Multiple Local Equilibria an Aaptive Search Summary: The information about market eman
More informationAn inductance lookup table application for analysis of reluctance stepper motor model
ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationFundamental Laws of Motion for Particles, Material Volumes, and Control Volumes
Funamental Laws of Motion for Particles, Material Volumes, an Control Volumes Ain A. Sonin Department of Mechanical Engineering Massachusetts Institute of Technology Cambrige, MA 02139, USA August 2001
More informationAdvanced Partial Differential Equations with Applications
MIT OpenCourseWare http://ocw.mit.eu 18.306 Avance Partial Differential Equations with Applications Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.eu/terms.
More informationSolution Methods. Steady convection-diffusion equation. Lecture 05
Solution Methods Steady convection-diffusion equation Lecture 05 1 Navier-Stokes equation Suggested reading: Gauss divergence theorem Integral form The key step of the finite volume method is to integrate
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationSection 2.1 The Derivative and the Tangent Line Problem
Chapter 2 Differentiation Course Number Section 2.1 The Derivative an the Tangent Line Problem Objective: In this lesson you learne how to fin the erivative of a function using the limit efinition an unerstan
More informationChapter 11: Feedback and PID Control Theory
Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by
More informationChapter 1 Overview: Review of Derivatives
Chapter Overview: Review of Derivatives The purpose of this chapter is to review the how of ifferentiation. We will review all the erivative rules learne last year in PreCalculus. In the net several chapters,
More informationON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationElectromagnet Gripping in Iron Foundry Automation Part II: Simulation
www.ijcsi.org 238 Electromagnet Gripping in Iron Founry Automation Part II: Simulation Rhythm-Suren Wahwa Department of Prouction an Quality Engineering, NTNU Tronheim, 7051, Norway Abstract This paper
More information3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes
Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we
More informationOptimization of Geometries by Energy Minimization
Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.
More informationCONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS
BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN
More informationCS9840 Learning and Computer Vision Prof. Olga Veksler. Lecture 2. Some Concepts from Computer Vision Curse of Dimensionality PCA
CS9840 Learning an Computer Vision Prof. Olga Veksler Lecture Some Concepts from Computer Vision Curse of Dimensionality PCA Some Slies are from Cornelia, Fermüller, Mubarak Shah, Gary Braski, Sebastian
More informationUNIFYING PCA AND MULTISCALE APPROACHES TO FAULT DETECTION AND ISOLATION
UNIFYING AND MULISCALE APPROACHES O FAUL DEECION AND ISOLAION Seongkyu Yoon an John F. MacGregor Dept. Chemical Engineering, McMaster University, Hamilton Ontario Canaa L8S 4L7 yoons@mcmaster.ca macgreg@mcmaster.ca
More informationSturm-Liouville Theory
LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory
More informationA Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges
Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationA SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS
International Journal on Engineering Performance-Base Fire Coes, Volume 4, Number 3, p.95-3, A SIMPLE ENGINEERING MOEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PROCTS V. Novozhilov School of Mechanical
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationDusty Plasma Void Dynamics in Unmoving and Moving Flows
7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND SPACE SCIENCES (EUCASS) Dusty Plasma Voi Dynamics in Unmoving an Moving Flows O.V. Kravchenko*, O.A. Azarova**, an T.A. Lapushkina*** *Scientific an Technological
More informationPlanar sheath and presheath
5/11/1 Flui-Poisson System Planar sheath an presheath 1 Planar sheath an presheath A plasma between plane parallel walls evelops a positive potential which equalizes the rate of loss of electrons an ions.
More informationIPMSM Inductances Calculation Using FEA
X International Symposium on Inustrial Electronics INDEL 24, Banja Luka, November 68, 24 IPMSM Inuctances Calculation Using FEA Dejan G. Jerkan, Marko A. Gecić an Darko P. Marčetić Department for Power,
More informationMathematics 116 HWK 25a Solutions 8.6 p610
Mathematics 6 HWK 5a Solutions 8.6 p6 Problem, 8.6, p6 Fin a power series representation for the function f() = etermine the interval of convergence. an Solution. Begin with the geometric series = + +
More informationCode_Aster. Detection of the singularities and calculation of a map of size of elements
Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : DLMAS Josselin Clé : R4.0.04 Révision : Detection of the singularities an calculation of a map of size of
More informationTransmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency
Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com
More informationSemiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom
PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,
More informationQuantum Search on the Spatial Grid
Quantum Search on the Spatial Gri Matthew D. Falk MIT 2012, 550 Memorial Drive, Cambrige, MA 02139 (Date: December 11, 2012) This paper explores Quantum Search on the two imensional spatial gri. Recent
More informationarxiv: v1 [hep-lat] 19 Nov 2013
HU-EP-13/69 SFB/CPP-13-98 DESY 13-225 Applicability of Quasi-Monte Carlo for lattice systems arxiv:1311.4726v1 [hep-lat] 19 ov 2013, a,b Tobias Hartung, c Karl Jansen, b Hernan Leovey, Anreas Griewank
More informationChapter 11: Feedback and PID Control Theory
Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationBalancing Expected and Worst-Case Utility in Contracting Models with Asymmetric Information and Pooling
Balancing Expecte an Worst-Case Utility in Contracting Moels with Asymmetric Information an Pooling R.B.O. erkkamp & W. van en Heuvel & A.P.M. Wagelmans Econometric Institute Report EI2018-01 9th January
More informationd dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1
Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of
More information2-7. Fitting a Model to Data I. A Model of Direct Variation. Lesson. Mental Math
Lesson 2-7 Fitting a Moel to Data I BIG IDEA If you etermine from a particular set of ata that y varies irectly or inversely as, you can graph the ata to see what relationship is reasonable. Using that
More informationKinematics of Rotations: A Summary
A Kinematics of Rotations: A Summary The purpose of this appenix is to outline proofs of some results in the realm of kinematics of rotations that were invoke in the preceing chapters. Further etails are
More informationPhysics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2
Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationChapter 11: Feedback and PID Control Theory
Chapter 11: Feeback an ID Control Theory Chapter 11: Feeback an ID Control Theory I. Introuction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works
More informationSurvey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013
Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing
More informationRevisiting the Charge Concept in HBT/BJT Models
evisiting the Charge Concept in HBT/BJT Moels Zoltan Huska an Ehrenfrie Seebacher austriamicrosystems AG 23r Bipolar Arbeitkeis BipAK Meeting at STM Crolles, France, 5 October 2 Outline recalling the junction
More informationShort Intro to Coordinate Transformation
Short Intro to Coorinate Transformation 1 A Vector A vector can basically be seen as an arrow in space pointing in a specific irection with a specific length. The following problem arises: How o we represent
More informationSemi-continuous challenge
Maurice Clerc Maurice.Clerc@WriteMe.com Last upate: 004-0-7. Why? Semi-continuous challenge The aim is to improve my parameter free (aaptive) particle swarm optimizer Tribes [CLE 03]. In orer to o that,
More informationAn Analytical Expression of the Probability of Error for Relaying with Decode-and-forward
An Analytical Expression of the Probability of Error for Relaying with Decoe-an-forwar Alexanre Graell i Amat an Ingmar Lan Department of Electronics, Institut TELECOM-TELECOM Bretagne, Brest, France Email:
More informationarxiv:cond-mat/ v1 12 Oct 2005
Creation an estruction of a spin gap in weakly couple quarter-fille laers B. Eegger,, H.G. Evertz, an R.M. Noack Institut für Theoretische Physik, Technische Universität Graz, A-8 Graz, Austria Institut
More informationConstraint Reformulation and a Lagrangian Relaxation based Solution Algorithm for a Least Expected Time Path Problem Abstract 1.
Constraint Reformulation an a Lagrangian Relaation base Solution Algorithm for a Least Epecte Time Path Problem Liing Yang State Key Laboratory of Rail Traffic Control an Safety, Being Jiaotong University,
More information