Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency
|
|
- Theodore McKenzie
- 5 years ago
- Views:
Transcription
1 Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com
2 Abstract In this paper the influence of asorption an trapping processes as escribe by TLM are investigate as a function of scaling of the spatial iscretisation. It is observe that there is a continuum between the macro-scale an the micro-scale that can be moelle by TLM. The resolution of some initial ifficulties in interpreting the effect of scaling of the 'capacity ratio' has provie useful insights into the nature of that parameter.
3 1. Introuction The first paper in this series [1] introuce the process of constructing TLM analogues for iffusion with asorption/trapping. The electromagnetic network can be solve exactly, so the only approximation in this approach is how well the construct escribes the physical process being moelle. In general we get reasonable agreement between theory an experiment at the macro-level. This paper is concerne with consistency; whether the unerlying PDE which is escribe by a particular scheme remains unchange as we scale own a moel from macro level to micro level. In many moelling scenarios these two levels of abstraction are consiere as istinct an the intermeiate cases ten to be avoie, which is unfortunate. It has been shown that simple iffusion processes which are escribe by Fick's secon law reuce to a Telegraphers' equation as Δx an Δt ten towar zero when the ratio Δx/Δt is hel constant [,3]. So long as the correct space/time scaling is use it appears that there is inee a continuum between the macro- an micro-levels of moelling. Here we will exten this work an emonstrate that the network moels for iffusion with asorption an trapping are inee consistent with the physical processes as space an time iscretisations are scale with constant velocity. It will be seen that although the parameters which escribe scattering in the presence of asorption/trapping appear complicate, they possess an unerlying simplicity which can be relate to more funamental moels. The paper starts by consiering the relationship between the effective iffusion constant an the TLM scattering parameters. This is followe by an analysis of the scattering processes as moelle by TLM. We are then in a position to consier what happens when such network moels are subjecte to space/time scaling. It will be shown that although the starting point may appear more complicate the results are not out of line with what has been previously observe [1,,3,].. The effective iffusion constant an the TLM scattering parameters In part A of this paper we ha an effective iffusion constant, which for the Langmuir isotherm is D eff D / (1 + K, where K was the constant or proportionality relating the concentration of immobilise iffusant to the concentration of free iffusant. This can be relate to the TLM moel for an iniviual pulse arriving at a noe which we can now consier as a scattering/asorbing site. We use the same efinitions as before [1] for the reflection, transmission an retention coefficients, which for ease of reference are reprouce in Appenix I. We can express the conservation of mass flux in a lossy network which inclues stubs as where ρ ( ρ LL I τ ( τ LL I an σ ( τ LS I ρ + τ + σ 1 (1 Thus σ is the fraction that is absorbe an ρ + τ represents fraction of free iffusant (1 σ For the Langmuir isotherm the factor 1 the fraction of free-iffusant (1 σ (1+ K
4 Thus the effective iffusion constant which has been erive strictly from a TLM formulation D eff D(1 σ is consistent with Eqn ( in [1], D eff D / (1+ K, which was erive from conventional calculus consierations. 3. An insight into the physical processes as moelle by TLM We can rearrange the expression for the one-imensional line-line reflection coefficient in an informative way Since ρ LL ( S + Z S + Z ( + Z( + Z + S We rewrite this as ρ LL ( S + Z S ( + Z( + Z + S + ( + Z( Z ( + Z( + Z + S or ρ LL ( S + Z S ( + Z( + Z + S + ( + Z ( + Z( + Z + S Z ( + Z + S The first two terms can be amalgamate an reference to the efinition of τ LS gives ρ LL ( S + ( + Z ( + Z( + Z + S τ LS ( which reuces to the one-imensional reflection coefficient in the absence of iffusion with an aitional term ρ LL ( + Z τ LS ( + Z σ (3a Using an ientical line of argument we can show that one imensional transmission coefficient can be represente as τ LL Z ( + Z τ LS Z ( + Z σ (3b These results o not require that ρ SS 0 an hence are applicable to both asorption an trapping processes. We have effectively prove that the one-imensional scattering coefficients in the presence of asorption/trapping is exactly what it woul be in the absence of asorption/trapping, but with half of the absorption component apportione to reflection an half to transmission, inepenent of the relative magnitues of these two coefficients.
5 . The scaling of the lossy TLM scattering parameters in the presence of asorption e Cogan [] has shown how the reflection coefficient (in the absence of asorption can be scale an how this provies a continuum of results between macro an micro-scale moels. This means that if we start with the reflection coefficient in the absence of absorption ρ LL ( + Z 1 1+ D Δt (a Δx an if we choose a unit propagation velocity which we assume remains constant uner scaling. Then, as Δx Δx/f an Δt Δt/f, Eqn (a becomes ρ LL 1 1+ Df (b The reflection coefficient scales as shown, an as f it tens towars zero:- the problem moves from a iffusion which can be expresse by a parabolic equation, through the telegraphers' equation an ultimately to a purely wave equation. The succeeing sections of this paper investigate how the various scattering terms which are involve in the process of iffusion plus asorption/trapping behave as we move from the macro- to the micro-scale. The basic noe of the transmission line network analogue is shown in figure 1 where an Z represent the components which contribute to iffusion while s an Z s both contribute to trapping an storage. In our treatment we will move from these values to istribute values (base on resistance or capacitance per unit length. These are summarise below Figure 1 A one-imensional lossy noe (length Δx with an open-circuit stub of length Δx/ an impeance Z S, which is in series with a resistor, S. Z Δt C Δx Z s Δt Δx s Δx Δx (because there are two resistors within a noe (because the stub-length is Δx/ (because the stub-length is Δx/
6 .1 The transmission line to transmission line reflection coefficient ρ LL ( S + Z S + Z ( + Z( + Z + S (5a When the above substitutions are mae in the expression for the line-to-line reflection coefficient an when the scalings Δx Δx/f an Δt Δt/f are inclue, then the ivien (top part of eqn (5a becomes while the ivisor becomes f D f (5b D + + D f + D D + C f (5c ( ρ LL I D + 8 f D f + D f + D D + C f (5 This is consistent in that if either is infinite or if is zero (Z s is infinite then eqn (5 reverts to eqn (b. However, because of the nature of the ivien the reflection coefficient can assume negative values an we nee to investigate this an what it means. As a first step we can consier the case of infinite scaling when Δx an Δt ten towars zero (f. In the extreme ( ρ LL I D D + Z Z + Z s (5e The relationship on the right above is obtaine by multiplying the mile expression above an below by Δx/Δt an using D 1/( C, Z Δt/C an Z s Δt/( Δx/. If stub capacitance is zero (stub impeance is infinite then the reflection coefficient is zero, which is consistent with the case in eqn (b when f. Using D 1/( C, Z Δt/C an Z s Δt/( Δx/ simple rearrangement an multiplication above an below by Δx/Δt yiels the result on the right of (5e, namely that in the extreme as f, an s play no part. The negative sign in (5e is then the result of a signal travelling on a lossless line Z encountering a pair of lossless lines, one of which has impeance Z an the other Z s.
7 . The scaling of the line-to-line transmission coefficient ( τ LL I Z( S + Z S ( + Z( + Z + S (6a The ivien is D f +16D 1 f (6b while the ivisor is ientical with eqn (5c Thus ( τ LL I 1+ + D + D f + 16D 1 f + D f + D D + C f (6c This is consistent in that if either is infinite or if is zero (Z s is infinite then the line-to-line transmission coefficient becomes unity. When the scaling becomes very large most terms in (6c isappear an we are left with a situation where s plays no part. ( τ LL I (6 D + If stub capacitance is zero (stub impeance is infinite then the coefficient is unity. Simple rearrangement an multiplication above an below by Δx/Δt yiels (6e, which shows that also plays no part. We simply have lossless line-to-line transmission at the iscontinuity where the stub is joine to the main network. ( τ LL I Z s Z + Z s (6e.3 The scaling of the line-to-stub transmission coefficient ( τ LS I Z ( + Z + S The ivisor is ientical to (5c, while the ivien is Df+8D f (7a (7b so ( τ LS I D + Df + 8D f + D f + D D + C f (7c It is clear that if either is infinite or if is zero (Z s is infinite then the line-to-stub transmission coefficient becomes zero. However if is infinite an/or if is non-zero then as
8 f we have a situation which is inepenent of the stub resistance, but epens only on the ratio of the line an stub istribute capacitances. D ( τ LS I D + This can be translate into network impeances as above an we have Z ( τ LS I Z + Z z (7 (7e. The scaling of the stub-to-stub reflection coefficient ( ρ SS I + Z + Z S S + Z + S (8a The ivien is 1+ + D - + D f + D D - C f (8b so ( ρ SS I D D + + D + D f + D D - C f f + D D + C f (8c This is perhaps the most variable of the scattering coefficients in terms of the way in which the various parameters can influence its polarity. However, as we scale to large values of f then we once again get a situation where s plays no part D - ( ρ SS I (8 D + In fact, plays no part at the extreme where f an this can be seen by transforming (8 back into impeances Z ( ρ SS I - Z s (8e Z + Z s A signal travelling along the stub an approaching the iscontinuity sees two impeances (each of magnitue, Z in parallel.
9 .5 The scaling of the stub-to-line transmission coefficient ( τ SL I Z S + Z + S (9a 1 The ivien is 8 f + 16D 1 f (9b 1 8 f + 16D 1 f so ( τ SL I 1+ + D + + D f + D D + C f Which is positive sign for all network parameters. As f the above reuces to (9c 1 ( τ SL I (9 D + which is ientical to (6, the expression for the line-to-line transmission coefficient. Similarly, the transformation to impeances is ientical to (6e. Z ( τ SL I s (9e Z + Z s.6 The scaling of eqn (1 Having evelope scaleable expressions for sum is always unity, inepenent of scaling. ( ρ LL I,( τ LL I an ( τ LS I it is easy to show that their.7 Scaling of the capacity ratio In orer to scale this quantity correctly we must recognise that is represents the ratio between the concentration of trappe species to the concentration of free species. Sections have been concerne with the scaling of flux-relate parameters which are not applicable here. Instea we must consier the voltage scattering parameters which are reprouce in Appenix I. We consier the single event where species enters the noe. Some is transmitte, some reflecte an some trappe. The ratio of trappe to free is then K ( τ LS V ( ρ LL V + τ LL ( V (10a
10 If we use the same approach to scaling as use in the previous sections we obtain for this ratio K( f 8 f f Df (10b In the extreme, as f this reuces to K Z s Z s Z (10c Which is exactly what we woul get if we were to take the expressions for the parameters in (10a an set both an s 0. As elsewhere, the lossy components are reuce as the problems is scale own an tens towars zero as f.8 Discussion on scaling It is known that in the extreme of infinite scaling even a simple iffusion formulation reuces to a lossless wave-propagation treatment. Of course we are never at quite this extreme, but rather at the level where the physical processes can be moelle by the telegraph equation. It has been emonstrate here that all of the scattering parameters which are involve in iffusion with asorption/trapping/recombination o inee scale an in the extreme, they too behave like propagation on a lossless network. We have emonstrate that matter is conserve at all levels of scaling an although it has come as a surprise that the capacity ratio oes not remain constant with scaling, we have been able to make sense of what has been obtaine. In spite of what has been sai above we live in a finite worl. We o not expect real physical processes to operate at the level where Δx an Δt are actually zero. However we woul like to be able to treat them where the iscretisations are at an atomic level an the question of an appropriate moification to the telegraph equation must now be consiere. 5. TLM moels for iffusion with asorption an the telegraphers' equation It is easiest to evelop an appropriate set of equation by using a parallel arrangement of two stubs of length Δx, instea of the one of length Δx/ which we have been using up to now. This equivalence allows us to use a metho similar to that given by Christopoulos [5]. A pseuo twoimensional representation is shown in figure where the arms in the y-irection represent the one-imensional line an the arms in the x-irection represent the lossy stub-line which is terminate at both ens by infinite (open-circuit impeances. We will assume that all transmission lines have the same istribute inuctance, L an that ifferences in impeance are expresse in the ifferent capacitances, C for Z an C S for Z S. In all treatments up to now we have use s which is given by s Δx However, since the stub length here is Δx we use s ' Δx
11 Figure : A simplifie network for one-imensional iffusion with asorption, where all lengths are reuce to the same iscretisations, Δx. In the x-irection we have V x L i x t i x an i x x C ' i x t (11a where C ' C + C S Δx In the y-irection we have V y L i y t i y an i y y C ' i y t (11b These can be manipulate so that the currents (fluxes are eliminate to give V x L ' C V t + ' C V t an V y L ' C V t + ' C V t (1 The sum is V x + V y L ( C + V t + ( + ( C + V t (13
12 Now since we are measuring the voltage (concentration at the centre of the noe in figure an since Δx Δy we have the telegraph equation for one-imensional iffusion with asorption/trapping V x L ( C + V t + ( + ( C + V t (1 This equation is entirely consistent for the stanar lossy telegraph equation if s 0 an C s 0. V x L C V t + C V t (15 We see that the inclusion of a provision for asorption has two effects. In the first case the iffusion constant is reuce, which is generally accepte to be the case, but at the microscopic level we also see that the wave-like character is slowe own ue to loaing by the capacitance in the stub. 6. Conclusions It has been emonstrate here that all of the TLM scattering parameters which are involve in iffusion with asorption/trapping/recombination o inee scale in a way which is at all times consistent with the unerlying physics. It has been particularly interesting to see how the associate TLM networks can provie fresh insights into the physical processes which they are being use to moel. This inclues the 'effective iffusion constant' which is applicable in chromatography/electrophoresis. This in turn utilises a parameter calle the capacity ratio an the way in which this scales was not obvious at first sight. This eficiency probably has more to o with the inconsistencies of efinition of the parameter itself, which effectively fuses ieas of surface concentration (trappe species an volume concentration (free species, which only become resolve as scaling tens towars infinity. While this paper has commente on the fact that much research elsewhere iscusse either macro scale or micro scale processes while ignoring the intermeiate cases, in the en it has been guilty of the same fault. However, expressions for the scattering parameters an unerlying equations as a function of scaling have been evelope, but their investigation must be part of another paper.
13 eferences 1. D. e Cogan Transmission Line Matrix (TLM network analogues of reversible trapping processes (Part A: introuction an methoology D. e Cogan, The relationship between parabolic an hyperbolic transmission line matrix moels for heat-flow Microelectronics Journal 30 ( D. e Cogan, Some observations on spurious effects in Numerical Moels International Journal of Mathematical Algorithms 1( ( D. e Cogan, X. Gui an M. ak, Some Observations on the TLM Numerical Solution of the Laplace Equation, Journal of Mathematical Moelling an Algorithms, vol. 8, no., pp , Dec C. Christopoulos, The Transmission Line Moelling Metho IEEE/OUP 1995, p. 171
14 Appenix I The various one-imensional TLM current/flux scattering coefficients which were erive as Eqn (5 in [1] are reprouce here for ease of reference ( ρ LL I ( τ LL I ( τ LS I ( S + Z S + Z ( + Z( + Z + S Z( S + Z S ( + Z( + Z + S Z ( + Z + S line-to-line reflection line-to-line transmission line-to-stub transmission (asorption ( ρ SS I + Z + S Z S + Z + S ( τ SL I Z S + Z + S stub-to-stub reflection (retention stub-to-line transmission (release The voltage/concentration scattering parameters are ( ρ LL V ( ρ LL I line-to-line reflection ( τ LL V ( τ LL I line-to-line transmission ( τ LS V Z S ( + Z ( + Z( + Z + S line-to-stub transmission (asorption ( ρ SS V ( ρ SS I stub-to-stub reflection (retention ( τ SL V Z( + Z ( ( + Z + Z + S stub-to-line transmission (release note that some of these parameters have been ivie above an below by ( + Z in orer to clarify some of the erivations in the text.
Chapter 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 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 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 informationHomework 7 Due 18 November at 6:00 pm
Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
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 informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
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 informationConsider for simplicity a 3rd-order IIR filter with a transfer function. where
Basic IIR Digital Filter The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient ifference equation
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationBasic IIR Digital Filter Structures
Basic IIR Digital Filter Structures The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient
More information. Using a multinomial model gives us the following equation for P d. , with respect to same length term sequences.
S 63 Lecture 8 2/2/26 Lecturer Lillian Lee Scribes Peter Babinski, Davi Lin Basic Language Moeling Approach I. Special ase of LM-base Approach a. Recap of Formulas an Terms b. Fixing θ? c. About that Multinomial
More informationTHE USE OF KIRCHOFF S CURRENT LAW AND CUT-SET EQUATIONS IN THE ANALYSIS OF BRIDGES AND TRUSSES
Session TH US O KIRCHO S CURRNT LAW AND CUT-ST QUATIONS IN TH ANALYSIS O BRIDGS AND TRUSSS Ravi P. Ramachanran an V. Ramachanran. Department of lectrical an Computer ngineering, Rowan University, Glassboro,
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
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 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 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 informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
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 informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS
ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an
More informationCenter of Gravity and Center of Mass
Center of Gravity an Center of Mass 1 Introuction. Center of mass an center of gravity closely parallel each other: they both work the same way. Center of mass is the more important, but center of gravity
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 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 informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More information1 The Derivative of ln(x)
Monay, December 3, 2007 The Derivative of ln() 1 The Derivative of ln() The first term or semester of most calculus courses will inclue the it efinition of the erivative an will work out, long han, a number
More informationPhysics 2212 K Quiz #2 Solutions Summer 2016
Physics 1 K Quiz # Solutions Summer 016 I. (18 points) A positron has the same mass as an electron, but has opposite charge. Consier a positron an an electron at rest, separate by a istance = 1.0 nm. What
More informationStability of solution of the diffusion equation using various formulation of finite difference method
Stability of solution of the diffusion equation using various formulation of finite difference method Michael J. Malachowski * Abstract A probabilistic approach has been used to analyze the stability of
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 observation on developments in lossy TLM
Some observation on developments in lossy TLM Donard de Cogan School of Information Systems, UEA, Norwich N4 7TJ (UK) Xiang Gui Department of Electrical and Computer Engineering, The University of Calgary,
More informationDynamics of the synchronous machine
ELEC0047 - Power system ynamics, control an stability Dynamics of the synchronous machine Thierry Van Cutsem t.vancutsem@ulg.ac.be www.montefiore.ulg.ac.be/~vct These slies follow those presente in course
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
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 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 informationEfficient Macro-Micro Scale Coupled Modeling of Batteries
A00 Journal of The Electrochemical Society, 15 10 A00-A008 005 0013-651/005/1510/A00/7/$7.00 The Electrochemical Society, Inc. Efficient Macro-Micro Scale Couple Moeling of Batteries Venkat. Subramanian,*,z
More informationCalculus in the AP Physics C Course The Derivative
Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.
More informationA Second Time Dimension, Hidden in Plain Sight
A Secon Time Dimension, Hien in Plain Sight Brett A Collins. In this paper I postulate the existence of a secon time imension, making five imensions, three space imensions an two time imensions. I will
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 information6 Wave equation in spherical polar coordinates
6 Wave equation in spherical polar coorinates We now look at solving problems involving the Laplacian in spherical polar coorinates. The angular epenence of the solutions will be escribe by spherical harmonics.
More informationPCCP PAPER. 1 Introduction. A. Nenning,* A. K. Opitz, T. M. Huber and J. Fleig. View Article Online View Journal View Issue
PAPER View Article Online View Journal View Issue Cite this: Phys. Chem. Chem. Phys., 2014, 16, 22321 Receive 4th June 2014, Accepte 3r September 2014 DOI: 10.1039/c4cp02467b www.rsc.org/pccp 1 Introuction
More informationHarmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method
1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse
More informationDesigning Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations
EECS 6B Designing Information Devices an Systems II Fall 07 Note 3 Secon Orer Differential Equations Secon orer ifferential equations appear everywhere in the real worl. In this note, we will walk through
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationModeling time-varying storage components in PSpice
Moeling time-varying storage components in PSpice Dalibor Biolek, Zenek Kolka, Viera Biolkova Dept. of EE, FMT, University of Defence Brno, Czech Republic Dept. of Microelectronics/Raioelectronics, FEEC,
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationAngles-Only Orbit Determination Copyright 2006 Michel Santos Page 1
Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More informationStrength Analysis of CFRP Composite Material Considering Multiple Fracture Modes
5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation
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 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 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 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 informationError Floors in LDPC Codes: Fast Simulation, Bounds and Hardware Emulation
Error Floors in LDPC Coes: Fast Simulation, Bouns an Harware Emulation Pamela Lee, Lara Dolecek, Zhengya Zhang, Venkat Anantharam, Borivoje Nikolic, an Martin J. Wainwright EECS Department University of
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 informationUnified kinetic model of dopant segregation during vapor-phase growth
PHYSICAL REVIEW B 72, 195419 2005 Unifie kinetic moel of opant segregation uring vapor-phase growth Craig B. Arnol 1, * an Michael J. Aziz 2 1 Department of Mechanical an Aerospace Engineering an Princeton
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
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 informationMODELLING DEPENDENCE IN INSURANCE CLAIMS PROCESSES WITH LÉVY COPULAS ABSTRACT KEYWORDS
MODELLING DEPENDENCE IN INSURANCE CLAIMS PROCESSES WITH LÉVY COPULAS BY BENJAMIN AVANZI, LUKE C. CASSAR AND BERNARD WONG ABSTRACT In this paper we investigate the potential of Lévy copulas as a tool for
More informationConstruction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems
Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu
More informationARCH 614 Note Set 5 S2012abn. Moments & Supports
RCH 614 Note Set 5 S2012abn Moments & Supports Notation: = perpenicular istance to a force from a point = name for force vectors or magnitue of a force, as is P, Q, R x = force component in the x irection
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 informationLyapunov Functions. V. J. Venkataramanan and Xiaojun Lin. Center for Wireless Systems and Applications. School of Electrical and Computer Engineering,
On the Queue-Overflow Probability of Wireless Systems : A New Approach Combining Large Deviations with Lyapunov Functions V. J. Venkataramanan an Xiaojun Lin Center for Wireless Systems an Applications
More informationMULTI-SCALE METHODS FOR THE SOLUTION OF THE RADIATIVE TRANSFER EQUATION
MULTI-SCALE METHODS FOR THE SOLUTION OF THE RADIATIVE TRANSFER EQUATION Pero J. Coelho *,, Nicolas Crouseilles **, Pero Pereira * an Maxime Roger *** * LAETA, IDMEC, Dept. of Mechanical Engineering, Instituto
More informationWe G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies
We G15 5 Moel Reuction Approaches for Solution of Wave Equations for Multiple Frequencies M.Y. Zaslavsky (Schlumberger-Doll Research Center), R.F. Remis* (Delft University) & V.L. Druskin (Schlumberger-Doll
More informationWhat s in an Attribute? Consequences for the Least Common Subsumer
What s in an Attribute? Consequences for the Least Common Subsumer Ralf Küsters LuFG Theoretical Computer Science RWTH Aachen Ahornstraße 55 52074 Aachen Germany kuesters@informatik.rwth-aachen.e Alex
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 information1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity
AP Physics Multiple Choice Practice Electrostatics 1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity. A soli conucting sphere is given a positive charge Q.
More informationPerfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs
Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs Ashish Goel Michael Kapralov Sanjeev Khanna Abstract We consier the well-stuie problem of fining a perfect matching in -regular bipartite
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationPRACTICE 4. CHARGING AND DISCHARGING A CAPACITOR
PRACTICE 4. CHARGING AND DISCHARGING A CAPACITOR. THE PARALLEL-PLATE CAPACITOR. The Parallel plate capacitor is a evice mae up by two conuctor parallel plates with total influence between them (the surface
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 informationDot trajectories in the superposition of random screens: analysis and synthesis
1472 J. Opt. Soc. Am. A/ Vol. 21, No. 8/ August 2004 Isaac Amiror Dot trajectories in the superposition of ranom screens: analysis an synthesis Isaac Amiror Laboratoire e Systèmes Périphériques, Ecole
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 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 informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationAverage value of position for the anharmonic oscillator: Classical versus quantum results
verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive
More informationAN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E.
AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathca Release 14 Khyruin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering an Applie Science Gonzaga University SDC
More informationfv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n
Chapter 11 Rossby waves Supplemental reaing: Pelosky 1 (1979), sections 3.1 3 11.1 Shallow water equations When consiering the general problem of linearize oscillations in a static, arbitrarily stratifie
More informationthere is no special reason why the value of y should be fixed at y = 0.3. Any y such that
25. More on bivariate functions: partial erivatives integrals Although we sai that the graph of photosynthesis versus temperature in Lecture 16 is like a hill, in the real worl hills are three-imensional
More informationState-Space Model for a Multi-Machine System
State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal
More informationQuantum mechanical approaches to the virial
Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationq = F If we integrate this equation over all the mass in a star, we have q dm = F (M) F (0)
Astronomy 112: The Physics of Stars Class 4 Notes: Energy an Chemical Balance in Stars In the last class we introuce the iea of hyrostatic balance in stars, an showe that we coul use this concept to erive
More informationA simple model for the small-strain behaviour of soils
A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:
More informationCalculating the effective delayed neutron fraction using Monte Carlo techniques
Calculating the effective elaye neutron fraction using Monte Carlo techniques 1. Introuction Steven C. van er Marck * an Robin Klein Meulekamp NRG, Petten, the Netherlans We present a true Monte Carlo
More informationA Review of Multiple Try MCMC algorithms for Signal Processing
A Review of Multiple Try MCMC algorithms for Signal Processing Luca Martino Image Processing Lab., Universitat e València (Spain) Universia Carlos III e Mari, Leganes (Spain) Abstract Many applications
More informationDiophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations
Diophantine Approximations: Examining the Farey Process an its Metho on Proucing Best Approximations Kelly Bowen Introuction When a person hears the phrase irrational number, one oes not think of anything
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 informationOrdinary Differential Equations: Homework 1
Orinary Differential Equations: Homework 1 M. Gameiro, J.-P. Lessar, J.D. Mireles James, K. Mischaikow January 12, 2017 2 Chapter 1 Motivation 1.1 Exercises Exercise 1.1.1. (Frictionless spring) Consier
More informationABDELSHAFY, OTHMAN, OSHMARIN, AL-MUTAWA, CAPOLINO: EPD IN CTLS UC IRVINE, SEP 2018
arxiv: [physics.app-ph] 13 SEP 2018 ABDELSHAFY, OTHMAN, OSHMARIN, AL-MUTAWA, CAPOLINO: EPD IN CTLS UC IRVINE, SEP 2018 Exceptional Points of Degeneracy in Perioically- Couple Waveguies an the Interplay
More informationLecture 6: Control of Three-Phase Inverters
Yoash Levron The Anrew an Erna Viterbi Faculty of Electrical Engineering, Technion Israel Institute of Technology, Haifa 323, Israel yoashl@ee.technion.ac.il Juri Belikov Department of Computer Systems,
More informationModel for Dopant and Impurity Segregation During Vapor Phase Growth
Mat. Res. Soc. Symp. Proc. Vol. 648, P3.11.1-7 2001 Materials Research Society Moel for Dopant an Impurity Segregation During Vapor Phase Growth Craig B. Arnol an Michael J. Aziz Division of Engineering
More informationThe Press-Schechter mass function
The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for
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 informationIntroduction to Markov Processes
Introuction to Markov Processes Connexions moule m44014 Zzis law Gustav) Meglicki, Jr Office of the VP for Information Technology Iniana University RCS: Section-2.tex,v 1.24 2012/12/21 18:03:08 gustav
More informationMark J. Machina CARDINAL PROPERTIES OF "LOCAL UTILITY FUNCTIONS"
Mark J. Machina CARDINAL PROPERTIES OF "LOCAL UTILITY FUNCTIONS" This paper outlines the carinal properties of "local utility functions" of the type use by Allen [1985], Chew [1983], Chew an MacCrimmon
More informationIntegration Review. May 11, 2013
Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More information