COUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS
|
|
- Maximillian Gardner
- 5 years ago
- Views:
Transcription
1 VI International Conference on Computational Methos for Couple Problems in Science an Engineering COUPLED PROBLEMS 15 B. Schrefler, E. Oñate an M. Paparakakis(Es) COUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS Jan Norström, Fatemeh Ghasemi Department of Mathematics, Computational Mathematics, Linköping University, SE Linköping, Sween. web page: janno11/ Department of Mathematics, Computational Mathematics, Linköping University, SE Linköping, Sween - Web page: fatgh43/ Key wors: multi-physics problem, well pose problems, stability, coupling proceure, high orer finite ifferences, summation-by-parts operatorsm, weak interface conitions Abstract. We iscuss well-poseness an stability of multi-physics problems by stuying a moel problem. By applying the energy metho, bounary an interface conitions are erive such that the continuous an semi-iscrete problem are well-pose an stable. The numerical scheme is implemente using high orer finite ifference operators on summation-by-parts (SBP) form an weakly impose bounary an interface conitions. Numerical experiments involving a spectral analysis corroborate the theoretical finings. 1 INTRODUCTION Roughly speaking, a well pose initial bounary value problem require that a unique solution that can be estimate in terms of the ata, exist. The most common proceure for showing well poseness is the so calle energy-metho where one multiply the governing partial ifferential equations (PDEs) with the solution, integrate by parts an impose bounary conitions [1]. The same general knowlege is not wie-sprea when it comes to the mathematical coupling of multi-physics problems. The reason for that is the more complex an to some extent more unclear nature of coupling conitions compare to imposing bounary conitions. Firstly, accuracy relations must exist such that combinations of variables for one set of PDEs at the interface is equal to combinations of variables for the other set. Seconly, the number of accuracy relations must fit both problems. Too many conitions ruin existence an too few ruin uniqueness. If the number of accuracy relations are too few, aitional conitions requiring external ata must be ae. Thirly, the accuracy relations must be such that no artificial growth or ecay is generate. 1
2 We will investigate the problems mentione above an generalize the investigation in [, 3] where we erive the coupling conitions by only emaning a well pose problem. Coupling of hyperbolic PDEs of ifferent size at the interface will be our primary focus. Once the coupling conitions are known for the continuous multi-physics problem we will iscretize using high orer finite ifferences on summation-by-parts form an inclue the coupling conitions weakly using the SAT technique [4, 5]. THE MODEL PROBLEM We will consier the following system, an the scalar equation u t + Au x =, 1 x, t >, (1) u(x, ) = f(x), v t + bv x =, x 1, t >, () v(x, ) = g(x). In (1), u = (u 1, u ) T is a vector of unknowns, f(x) = (f 1 (x), f (x)) T is a vector of given ata an for simplicity we choose [ ] a A =, a >. (3) a Two bounary/interface conitions are neee for the system (1) while equation () nees one bounary/interface conition..1 The interface conitions We apply the energy metho to both equations an a them together to get t ( u + α v ) = u T Au x= + αbv x= = w T Ew, (4) where α is positive free weight, w = [u 1, u, v] T an E is a E = a. (5) αb In (4), the bounary terms at the outer bounaries x = ±1 are ignore. The eigenvalues of E are {a, a, αb}. If b <, one of the eigenvalues is positive an we nee one conition at x =, otherwise we nee two conitions.
3 In orer to couple the problems we nee at least one accuracy conition. Let The relation (6) inserte in equation (4) leas to v = C T u, C = [c 1, c ] T. (6) t ( u + α v ) = u T (, t)du(, t), (7) where D = (αbcc T A). The characteristic polynomial relate to the eigenvalues λ of D is λ αb(c 1 + c )λ + αabc 1 c a. (8) To simplify the following iscussion we let s 1 = αb(c 1 + c ) an s = αabc 1 c a, which yiels the roots λ 1, = s 1 ± s 1 s. (9) First we consier b <. This leas to a positive s 1. If c 1 c a/bα, then s an both roots of the characteristic polynomial are negative, which means that D is negative efinite. This means that if c 1 an c have opposite sign, the couple problems satisfy an energy estimate for all choices of α. But if c 1 an c have the same sign, the energy estimate is not satisfie for any value of α. Consequently the couple problems with the interface conition v = C T u satisfy an energy estimate for b < if an only if c 1 an c have opposite signs. Next, consier b >. This leas to negative s 1 an at least one of the eigenvalues must be positive, which means that we nee an aitional conition. As mentione above, two conitions are neee at x =. One of them is an interface conition an the other one must be such that the right-han sie of (7) is negative semi-efinite. We will refer to this aitional conition as a bounary conition. If c 1 c a/bα, then s an one of the eigenvalues of D is positive (λ + ) an the other one is negative (λ ). Let D = Y ΛY T an rewrite (7) as t ( u + α v ) = u T (, t)(y ΛY T )u(, t), (1) where Λ = iag{λ +, λ } an Y is the matrix of eigenvectors to D. Let Λ = Λ + + Λ, where Λ + = iag{λ +, } an Λ = iag{, λ }. Furthermore we have D = D + + D where D + = Y Λ + Y T an D = Y Λ Y T. Then (1) leas to t ( u + α v ) = (Y T u(, t)) T (Λ + + Λ )(Y T u(, t)), (11) 3
4 The most general conition base on (11) is (Y T + R r Y T )u(, t) = h(t), x =, (1) where Y + an Y are the eigenvectors relate to the positive an negative eigenvalues, respectively. Letting h(t) = an inserting (1) into (11) leas to t ( u + α v ) = (λ + R rλ + )(Y T u(, t)). (13) If λ + R rλ +, then the right-han sie of (13) is boune an we have a well-pose coupling. Note that with R r = we have the so calle characteristic bounary conitions. Consequently an energy estimate is obtaine if c 1 an c are chosen such that c 1 c a/bα. This means that c 1 an c must be less than an arbitrary positive number that we can choose. In short: all values of c 1 an c lea to a well-pose problem if b >.. The semi-iscrete problem Let A be an M N matrix an B a P R matrix. The Kronecker prouct of these matrices is efine as a 11 B a 1N B A B =. (14) a M1 B a MN B First, we consier b <. The semi-iscrete SBP-SAT formulations of (1) an () are, u t + (D u A)u = (P 1 u E u N Σ)(C T ũ N v )e u N, v t + bd v v = P 1 v σ(v C T ũ N )e v. (15) In (15), the outer bounary conitions are ignore as in the continuous case, D u,v = Pu,v 1 Q u,v are the ifference operators, P u,v are positive efinite matrices an Q u,v satisfy Q u,v + Q T u,v = iag[ 1,, 1]. The iscrete gri functions, relate to the gri vectors x u = (x = 1,, x N = ) an x v = (y =,, y M = 1) are u = (u 1, u, u 1N, u N ), v = (v,, v M ). (16) The vectors e u N = (,,, 1, 1)T an e v = (1,, ) T are N 1 an M 1, respectively. EN u = iag[,, 1] an Ev = iag[1,, ] are N N an M M, respectively. The penalty matrix Σ is given by [ ] σ1 σ Σ =, (17) σ 3 σ 4 4
5 σ is penalty parameter an ũ N = [u 1N, u N ] T. Next, consier b >. The semi-iscrete SBP-SAT formulations of (1) an () are, u t + (D u A)u = (Pu 1 EN u Σ)(C T ũ N v )e u N + (Pu 1 EN u Ξ)((I N H)u e u N h), v t + bd v v = Pv 1 σ(v C T ũ N )e v. (18) where the penalty matrix Ξ an H are given by [ ] [ χ1 χ Ξ = 1 Rr, H = χ 3 χ 4 ] Y T. (19) The bounary ata h is efine as h = [, h] T. Also, in the following analysis we will use the iscrete norms u P u I = u T (P u I)u, v P v = v T P v v. ()..1 Stability conitions at the interface First we consier b <. The iscrete energy metho is applie to (15) by multiplying the two equations with u T (P u I) an v T P v, respectively. The SBP properties of D u,v yiels t ( u P u I + α v P v ) = ũ T NAũ N + α bv + ũ T NΣH + α σv (v C T ũ N ). (1) In(1), α is a positive weight (not necessarily the same as in the continuous case) an H = [C T ũ N v, C T ũ N v ] T. In orer to mimic the continuous case, we choose σ = c 1 αb/ an σ 4 = c αb/ an σ 1 = σ 3 =. The final penalty matrix in block form becomes Σ = αb/ [ C ]. By inserting that into (1) we get t ( u P u I + α v P v ) =ũ T NDũ N + α v (b + σ) σv C T ũ N (αb + α σ). () If we choose σ = αb/α, for α α the right-han sie of () will be boune ue to the continuous result above. Next, we consier b > an let h(t) =. Multiplying (18) by u T (P u I) an v T P v leas to t ( u P u I + α v P v ) ũ T N(D + Ξ H + (Ξ H) T )ũ N, (3) where we have chosen Σ an σ as for the case b <. By using Y Y T = I, we can rewrite the right-han sie of (3) as ũ T N(D + Ξ H + (Ξ H) T )ũ N =(Y T ũ N ) T (Λ + (Y T Ξ HY ) + (Y T Ξ HY ) T )(Y T ũ N ). (4) 5
6 Let Ξ = Y T Ξ an choose Ξ such that Ξ = iag( χ 1, χ ). We also use the following split, [ ] (Y (Y T T ũ N ) = + ũ N ) (Y T. (5) ũ N ) Now, we can rewrite (4) as, [ ũ T N(D + Ξ H + (Ξ H) (Y T T )ũ N = + ũ N ) (Y T ũ N ) ] T [ λ + + χ 1 R r χ 1 R r χ 1 λ ] [ (Y T + ũ N ) (Y T ũ N ) ]. (6) By the choice χ 1 = λ +, the right-han sie of (6) can be rewritten as (λ + R rλ + )(Y T ũ N ) λ + ((Y T + ũ N ) R r (Y T ũ N )), (7) which is negative ue to the continuous result. Consequently if we choose Σ = αb/ [ C ] [ ] λ, σ = αb/α, Ξ = Y +, then for α α, (18) is stable... Stability conitions at the left bounary In orer to have a well-pose problem, we nee one conition at x = 1. We consier the homogeneous bounary conition with R l < 1. The SAT term at x = 1 is (Pu 1 [ ] π1 π Π =, Ĥ = (X T + R l X T )u( 1, t) =, (8) π 3 π 4 E u Π)(I N Ĥ)u, where [ 1 Rl ] X T, (9) an X is the matrix of eigenvectors to A. It can be shown that an energy estimate is obtaine if π 1 = a/(4(r l + 1)), π =, π 3 = a/(4(r l 1)), π 4 =. (3)..3 Stability conitions at the right bounary For the case b <, one conition at x = 1 is also neee. We choose the homogeneuos v(1, t) =. The SAT term at x = 1 is Pv 1 where θ satisfies θv N ev N θ b/. (31) 6
7 3 THE SPECTRUM In this section, we consier the continuous an iscrete spectrum for our problem. 3.1 The spectrum for the continuous problem By applying the Laplace transform to (1) an () we get the following system of orinary ifferential equations sû + Aû x =, 1 x, sˆv + bˆv x =, x 1. (3) We have ignore the initial conitions, since they o not influence the spectra an make the ansatz û = e kx ψ an ˆv = e k3x ψ 3. This leas to [ ] A (si + Bk)Ψ =, B =, (33) b where Ψ = [ψ 1, ψ, ψ 3 ] T. This system of equations have a non-trivial solution only when et(si + Bk) =, which leas to k 1 = s, k a = s an k a 3 = s. b The general solution incluing the eigenvectors is ŵ = α 1 e s a x α e s a x α 3 e s b x 1, (34) where ŵ = [û 1, û, ˆv] T. The unknowns α 1, α an α 3, will be etermine by the bounary an interface conitions. First we consier b <, with the conitions (X T + R l X T )u( 1, t) =, C T u(, t) v(, t) =, (35) v(1, t) =, an R l 1. The interface an bounary conitions are such that the couple problem is well-pose. By applying these conitions to (34), we obtain c 1 + c c 1 c 1 Eα =, E = e s a R l e s a, (36) e s b where α = [α 1, α, α 3 ] T. A non-trivial solution, require et(e) = e s b (e s a (c c 1 ) R l e s a (c1 + c )) =. (37) 7
8 The zeros of et(e) which form the spectrum of (3) are s = { a ln( R l(c 1 +c ) c c 1 ) + naπi, n Z, if R l(c 1 +c ) c c 1 >, a l(c 1 +c ) ) + naπi + aπi, n Z, if R l(c 1 +c ) c c 1 <. c c 1 (38) The real part of s is negative if R l(c 1 +c ) c c 1 < 1. It is easy to verify that this hols for R l < 1 an c 1, c with opposite sign. This means that if we choose c 1, c an R l such that the couple problem leas to an energy estimate, then the real part of s will be negative. Recall that this require that c 1 an c must have opposite signs. Next, consier b >. In this case we have the conitions (X T + R l X T )u( 1, t) =, C T u(, t) v(, t) =, (39) (Y T + R r Y T )u(, t) =, an λ + Rrλ + an Rl 1. The couple problems with the conitions (39) satisfy an energy estimate. By applying (39) to (34) leas to the following system of equations c 1 + c c 1 c 1 Eα =, E = e s a R l e s a, (4) y 1 R r y R r y 1 R r y R r where α = [α 1, α, α 3 ] T. The zeros of et(e) in this case are { a s = ln( R l(y 1 R ry R r) y 1 +R ry R r ) + naπi, n Z, if R l(y 1 R ry R r) y 1 +R ry R r >, a ln( R l(y 1 R ry R r) ) + naπi + aπi, n Z, if R l(y 1 R ry R r) y 1 +R ry R r <. y 1 +R ry R r The real part of s is negative if (41) R l(y 1 R r y R r ) y 1 + R r y R r < 1. (4) Note that the eterminant of E is inepenent of c 1 an c. This means that if (4) hols, then the real part of s is negative for all c 1 an c. Recall that for all values of c 1 an c, suitable choices of α such that the couple problems satisfy an energy estimate coul be foun. This implies that there is no limitation on c 1 an c in both the energy an spectral analysis. However, recall that this require an aitional bounary conition. 3. The semi-iscrete spectrum Consier b <. The SBP-SAT approximation of (1) an (), incluing (35) is u t + (D u A)u = (Pu 1 E u 1 Π)(I N Ĥ)u + (Pu EN u Σ)(C T ũ N v )e u N, v t + bd v v = P 1 v σ(v C T ũ N )e v + P 1 v θv N e N. (43) 8
9 In orer to etermine the semi-iscrete spectrum, we follow [6] an rewrite (43) in matrix form as where W = [u 1, u, v] T an [ (Qu A) H i = bq v W t = P 1 (H i + H c )W, (44) ] [ P, P 1 1 u I = Pv 1 ]. (45) The penalty matrix H c which is zero except at the bounaries an interface has the structure ΠĤ... EC H c = T E σc T, (46) σ... θ where E = [σ, σ 4 ] T. Next, consier b >. The SBP-SAT approximation of (1) an (), with conitions (39) is u t + (D u A)u = (Pu 1 E u Π)(I N 1 Ĥ)u + (Pu EN u Σ)(C T ũ N v )e u N + (Pu 1 EN u Ξ)(I N H)u, v t + bd v v = Pv 1 σ(v C T ũ N )e v. (47) The approximation (47) can be written on the the form (44) where in this case ΠĤ... H c = EC T + Ξ H E σc T σ, (48)... while H i is the same as before an given above. The eigenvalues of the matrix P 1 (H i +H c ) form the iscrete spectrum of (43) an (47). 4 NUMERICAL RESULTS In this section, we use the metho of manufacture solution in orer to test the accuracy of the approximations. RK3 is use to iscretize time. We also iscuss the relation between the continuous an semi-iscrete spectrum. 9
10 SBP 1 SBP 4 SBP 63 SBP 84 N error rate error rate error rate error rate e- - e-3 - e-3-1e-3-4 6e e e e e e e e e e e e e-4.4 4e e e e e e e Table 1: error an rate q u for b <. SBP 1 SBP 4 SBP 63 SBP 84 N error rate error rate error rate error rate e-1-3e- - 3e- - 4e-3-4 4e e e e e-.5 4e e e e e e e e-4.5 6e e e e-4.1 7e e e Table : error an rate q v for b <. 4.1 Accuracy The analytical solution that we use in the metho of manufacture solution is q u = ln u 1 (x, t) = u (x, t) = cos(π(x t)), v(x, t) = sin(3π(x bt)). (49) The rate of convergence is calculate as ( (u N 1 1, u N 1 ) (u 1, u ) Pu I (u N 1, u N ) (u 1, u ) Pu I ) / ln ( N1 N ), q v = ln ( v N 1 v Pv v N v Pv ) / ln ( N1 N ),(5) where u 1, u an v are the analytical solutions an u N i 1, u N i an v N i are the corresponing numerical solutions with N i gri points. First, we consier b <. The choosen coefficients are α = α = 1, a = 1, b = 1. To have a well-pose problem, we choose R l < 1 an c 1, c such that c 1 c 1/. Let R l =.5 an c 1 = 1, c =. Tables 1 an show the error an convergence rate q u an q v, respectively, for SBP operators with th, 3 th, 4 th an 5 th orer. Next, we consier b >. We again choose α = α = 1, a = 1, b = 1 an take R l =.5, R r =.5, c 1 = 1 an c = 1 in orer to have a well-pose problem. Tables 3 an 4 show the error an convergence rates for q u an q v, respectively. Clearly, the esign orer of accuracy is obtaine. 1
11 SBP 1 SBP 4 SBP 63 SBP 84 N error rate error rate error rate error rate e- - 5e-3-3e-3 - e-3-4 6e e e e e-3.6 8e e e e e e e e-5.5 1e e e e e e e Table 3: error an rate q u for b >. SBP 1 SBP 4 SBP 63 SBP 84 N error rate error rate error rate error rate 3e- - 8e-3-7e-3-3e-3-4 7e-3. 1e e e e-3.3 1e e e e-4. e e e e-4. e e e e-5. 3e e e Table 4: error an rate q v for b >. 4. The spectrum of the continuous an semi-iscrete operators Figures 1-3 show the iscrete an continuous spectrum for ifferent gris using the SBP4 operator. One can clearly see the convergence of the iscrete spectrum to the continuous one as the gris are refine. This convergence hol both for positive an negative b an show that the solutions of the semi-iscrete scheme converge to the continuous one. 5 SUMMARY AND CONCLUSIONS We have iscusse well-poseness an stability of multi-physics problems by analyzing a moel problem. It was shown that for ceartain wave spees, only interface conitions were require, while in other cases aitional information in the form of bounary conitions must be supplie. By applying the energy metho, we erive bounary an interface conitions such that the continuous an semi-iscrete problem are well-pose an stable. The numerical scheme was implemente using high orer finite ifference operators on SBP form an weakly impose bounary an interface conitions using the SAT technique. It was shown that we obtaine esign orer of accuracy, an that the spectrum of the 11
12 6 4 N=8 N=16 N=3 N=64 Continuous N=8 N=16 N=3 N=64 Continuous Figure 1: Global view: the iscrete an continuous spectrum, b < (left) an b > (right) Imaginary part 5 N=8 N=16 N=3 N=64 Continuous 5 5 N=8 N=16 N=3 N=64 Continuous Real part Figure : Meium view: the iscrete an continuous spectrum, b < (left) an b > (right) N=8 N=16 N=3 N=64 Continuous N=8 N=16 N=3 N=64 Continuous Figure 3: Zoome view: the iscrete an continuous spectrum, b < (left) an b > (right). 1
13 iscrete operator converge to the spectrum of the continuous operator. The numerical experiments in combination with the theoretical erivations showe what type of analysis that is require to obtain accurate numerical simulations of multi-physics problems. Future work will inclue a generalization of this investigation for hyperbolic problems, an an extension to coupling of incompletely parabolic problems such as the compressible Navier-Stoke s equations. REFERENCES [1] Norström, J. an Svär, M. Well-pose bounary conitions for the Navier-Stokes equations. SIAM J. Numer. Anal. (5) 43: [] Norström, J. an Eriksson, S. Flui structure interaction problems: the necessity of a well pose, stable an accurate formulation. Commun. comput. phys(cicp). (1) 8: [3] Linström, J. an Norström, J. A stable an high orer accurate conjugate heat transfer problem. J. Comput. Phys. (1) 9: [4] Norström, J. an Berg, J. Conjugate heat transfer for the unsteay compressible Navier-Stokes equations using a multi-block coupling. Comput Fluis. (13) 7:- 9. [5] Norström, J., Gong, J., Van er Weie, E. an Svär, M. A stable an conservative high orer multi-block metho for the compressible Navier-Stokes equations. J. Comput. Phys. (9) 8: [6] Norström, J. an Carpenter, M. H. High-orer finite ifference methos, multiimensional linear problems an curvilinear coorinates. J. Comput. Phys. (1) 173:
Stable 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 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 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 informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
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 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 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 informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More information1.4.3 Elementary solutions to Laplace s equation in the spherical coordinates (Axially symmetric cases) (Griffiths 3.3.2)
1.4.3 Elementary solutions to Laplace s equation in the spherical coorinates (Axially symmetric cases) (Griffiths 3.3.) In the spherical coorinates (r, θ, φ), the Laplace s equation takes the following
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 informationWell-posedness, stability and conservation for a discontinuous interface problem: an initial investigation.
Well-posedness, stability and conservation for a discontinuous interface problem: an initial investigation. Cristina La Cognata and Jan Nordström Abstract A robust interface treatment for the discontinuous
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 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 informationAbstract A nonlinear partial differential equation of the following form is considered:
M P E J Mathematical Physics Electronic Journal ISSN 86-6655 Volume 2, 26 Paper 5 Receive: May 3, 25, Revise: Sep, 26, Accepte: Oct 6, 26 Eitor: C.E. Wayne A Nonlinear Heat Equation with Temperature-Depenent
More informationThe effect of dissipation on solutions of the complex KdV equation
Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,
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 informationLINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form
LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1 We consier ifferential equations of the form y + a()y = b(), (1) y( 0 ) = y 0, where a() an b() are functions. Observe that this class of equations inclues equations
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 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 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 informationProblems Governed by PDE. Shlomo Ta'asan. Carnegie Mellon University. and. Abstract
Pseuo-Time Methos for Constraine Optimization Problems Governe by PDE Shlomo Ta'asan Carnegie Mellon University an Institute for Computer Applications in Science an Engineering Abstract In this paper we
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 informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
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 informationFirst Order Linear Differential Equations
LECTURE 6 First Orer Linear Differential Equations A linear first orer orinary ifferential equation is a ifferential equation of the form ( a(xy + b(xy = c(x. Here y represents the unknown function, y
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 informationMath 300 Winter 2011 Advanced Boundary Value Problems I. Bessel s Equation and Bessel Functions
Math 3 Winter 2 Avance Bounary Value Problems I Bessel s Equation an Bessel Functions Department of Mathematical an Statistical Sciences University of Alberta Bessel s Equation an Bessel Functions We use
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationSpectral analysis of the incompressible Navier-Stokes equations with different boundary conditions
Spectral analysis of the incompressible Navier-Stokes equations with different boundary conditions Cristina La Cognata, Jan Nordström Department of Mathematics, Computational Mathematics, Linköping University,
More informationWell-posedness, Stability and Conservation for a Discontinuous Interface Problem
Department of Mathematics Well-posedness, Stability and Conservation for a Discontinuous Interface Problem Cristina La Cognata and Jan Nordström LiTH-MAT-R--214/16--SE Department of Mathematics Linköping
More informationChaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena
Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos
More informationDifferentiability, Computing Derivatives, Trig Review. Goals:
Secants vs. Derivatives - Unit #3 : Goals: Differentiability, Computing Derivatives, Trig Review Determine when a function is ifferentiable at a point Relate the erivative graph to the the graph of an
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 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 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 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 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 informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
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 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 informationThe Three-dimensional Schödinger Equation
The Three-imensional Schöinger Equation R. L. Herman November 7, 016 Schröinger Equation in Spherical Coorinates We seek to solve the Schröinger equation with spherical symmetry using the metho of separation
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 informationIPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy
IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation
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 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 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 informationTwo formulas for the Euler ϕ-function
Two formulas for the Euler ϕ-function Robert Frieman A multiplication formula for ϕ(n) The first formula we want to prove is the following: Theorem 1. If n 1 an n 2 are relatively prime positive integers,
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 informationPartial Differential Equations
Chapter Partial Differential Equations. Introuction Have solve orinary ifferential equations, i.e. ones where there is one inepenent an one epenent variable. Only orinary ifferentiation is therefore involve.
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 informationMany problems in physics, engineering, and chemistry fall in a general class of equations of the form. d dx. d dx
Math 53 Notes on turm-liouville equations Many problems in physics, engineering, an chemistry fall in a general class of equations of the form w(x)p(x) u ] + (q(x) λ) u = w(x) on an interval a, b], plus
More informationSection 2.7 Derivatives of powers of functions
Section 2.7 Derivatives of powers of functions (3/19/08) Overview: In this section we iscuss the Chain Rule formula for the erivatives of composite functions that are forme by taking powers of other functions.
More informationMathcad Lecture #5 In-class Worksheet Plotting and Calculus
Mathca Lecture #5 In-class Worksheet Plotting an Calculus At the en of this lecture, you shoul be able to: graph expressions, functions, an matrices of ata format graphs with titles, legens, log scales,
More informationMathematical Review Problems
Fall 6 Louis Scuiero Mathematical Review Problems I. Polynomial Equations an Graphs (Barrante--Chap. ). First egree equation an graph y f() x mx b where m is the slope of the line an b is the line's intercept
More information2.5 SOME APPLICATIONS OF THE CHAIN RULE
2.5 SOME APPLICATIONS OF THE CHAIN RULE The Chain Rule will help us etermine the erivatives of logarithms an exponential functions a x. We will also use it to answer some applie questions an to fin slopes
More informationVIBRATIONS OF A CIRCULAR MEMBRANE
VIBRATIONS OF A CIRCULAR MEMBRANE RAM EKSTROM. Solving the wave equation on the isk The ynamics of vibrations of a two-imensional isk D are given by the wave equation..) c 2 u = u tt, together with the
More informationParameter estimation: A new approach to weighting a priori information
Parameter estimation: A new approach to weighting a priori information J.L. Mea Department of Mathematics, Boise State University, Boise, ID 83725-555 E-mail: jmea@boisestate.eu Abstract. We propose a
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 informationImplicit Differentiation
Implicit Differentiation Implicit Differentiation Using the Chain Rule In the previous section we focuse on the erivatives of composites an saw that THEOREM 20 (Chain Rule) Suppose that u = g(x) is ifferentiable
More informationarxiv: v1 [math-ph] 5 May 2014
DIFFERENTIAL-ALGEBRAIC SOLUTIONS OF THE HEAT EQUATION VICTOR M. BUCHSTABER, ELENA YU. NETAY arxiv:1405.0926v1 [math-ph] 5 May 2014 Abstract. In this work we introuce the notion of ifferential-algebraic
More informationMath 115 Section 018 Course Note
Course Note 1 General Functions Definition 1.1. A function is a rule that takes certain numbers as inputs an assigns to each a efinite output number. The set of all input numbers is calle the omain of
More informationGLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS
Electronic Journal of Differential Equations, Vol. 015 015), No. 99, pp. 1 14. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu GLOBAL SOLUTIONS FOR D COUPLED
More informationSpace-time Linear Dispersion Using Coordinate Interleaving
Space-time Linear Dispersion Using Coorinate Interleaving Jinsong Wu an Steven D Blostein Department of Electrical an Computer Engineering Queen s University, Kingston, Ontario, Canaa, K7L3N6 Email: wujs@ieeeorg
More informationDifferentiability, Computing Derivatives, Trig Review
Unit #3 : Differentiability, Computing Derivatives, Trig Review Goals: Determine when a function is ifferentiable at a point Relate the erivative graph to the the graph of an original function Compute
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 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 informationANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS
ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics
More informationLecture 2 - First order linear PDEs and PDEs from physics
18.15 - Introuction to PEs, Fall 004 Prof. Gigliola Staffilani Lecture - First orer linear PEs an PEs from physics I mentione in the first class some basic PEs of first an secon orer. Toay we illustrate
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 informationINVERSE PROBLEM OF A HYPERBOLIC EQUATION WITH AN INTEGRAL OVERDETERMINATION CONDITION
Electronic Journal of Differential Equations, Vol. 216 (216), No. 138, pp. 1 7. ISSN: 172-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu INVERSE PROBLEM OF A HYPERBOLIC EQUATION WITH AN
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
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 informationStable finite difference schemes for the magnetic induction equation with Hall effect
Eigenössische Technische Hochschule Zürich Ecole polytechnique féérale e Zurich Politecnico feerale i Zurigo Swiss Feeral Institute of Technology Zurich Stable finite ifference schemes for the magnetic
More informationThe Generalized Incompressible Navier-Stokes Equations in Besov Spaces
Dynamics of PDE, Vol1, No4, 381-400, 2004 The Generalize Incompressible Navier-Stokes Equations in Besov Spaces Jiahong Wu Communicate by Charles Li, receive July 21, 2004 Abstract This paper is concerne
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More information1 Heisenberg Representation
1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.
More informationReview of Differentiation and Integration for Ordinary Differential Equations
Schreyer Fall 208 Review of Differentiation an Integration for Orinary Differential Equations In this course you will be expecte to be able to ifferentiate an integrate quickly an accurately. Many stuents
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 informationRobust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k
A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine
More informationStable Poiseuille Flow Transfer for a Navier-Stokes System
Proceeings of the 26 American Control Conference Minneapolis, Minnesota, USA, June 14-16, 26 WeB2.2 Stable Poiseuille Flow Transfer for a Navier-Stokes System Rafael Vázquez, Emmanuel Trélat an Jean-Michel
More informationFinal Exam Study Guide and Practice Problems Solutions
Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making
More informationCharacterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations
Characterizing Real-Value Multivariate Complex Polynomials an Their Symmetric Tensor Representations Bo JIANG Zhening LI Shuzhong ZHANG December 31, 2014 Abstract In this paper we stuy multivariate polynomial
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
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 informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationarxiv: v1 [cs.lg] 22 Mar 2014
CUR lgorithm with Incomplete Matrix Observation Rong Jin an Shenghuo Zhu Dept. of Computer Science an Engineering, Michigan State University, rongjin@msu.eu NEC Laboratories merica, Inc., zsh@nec-labs.com
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 informationDiscrete Operators in Canonical Domains
Discrete Operators in Canonical Domains VLADIMIR VASILYEV Belgoro National Research University Chair of Differential Equations Stuencheskaya 14/1, 308007 Belgoro RUSSIA vlaimir.b.vasilyev@gmail.com Abstract:
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 informationStudents need encouragement. So if a student gets an answer right, tell them it was a lucky guess. That way, they develop a good, lucky feeling.
Chapter 8 Analytic Functions Stuents nee encouragement. So if a stuent gets an answer right, tell them it was a lucky guess. That way, they evelop a goo, lucky feeling. 1 8.1 Complex Derivatives -Jack
More informationApplications of First Order Equations
Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more
More informationarxiv: v1 [cond-mat.stat-mech] 9 Jan 2012
arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,
More informationBEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi
BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica
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 information4. Important theorems in quantum mechanics
TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this
More informationUNDERSTANDING INTEGRATION
UNDERSTANDING INTEGRATION Dear Reaer The concept of Integration, mathematically speaking, is the "Inverse" of the concept of result, the integration of, woul give us back the function f(). This, in a way,
More informationOn the number of isolated eigenvalues of a pair of particles in a quantum wire
On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in
More informationQubit channels that achieve capacity with two states
Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March
More informationSome vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10
Some vector algebra an the generalize chain rule Ross Bannister Data Assimilation Research Centre University of Reaing UK Last upate 10/06/10 1. Introuction an notation As we shall see in these notes the
More information