ON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
|
|
- Jane Black
- 5 years ago
- Views:
Transcription
1 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 In this paper, a problem of shape esign for a uct with the flow governe by the one-imensional Euler equations is analyze. The flow is assume to be transonic, in the sense that we have a shock embee in the flow. The presence of the shock introuces analytical an numerical ifficulties that are overcome by introucing the shock location as an explicit state variable. We justify the use of ajoint variables by establishing conitions to ensure that the linearize constraint is surjective. These ajoint variables are use to state necessary optimality conitions an to compute graients. In aition, our characterization of the ajoint variables is of interest for comparison to ajoints associate with shock-capturing formulations.. Introuction In this paper we apply the shock-fitting framework evelope in Cliff et al 3 to an optimal uct esign problem governe by the one imensional Euler equations. The objective of this research is the rigorous mathematical justification of ifferentiability of the problem, surjectivity of the linearize constraints, an the accurate formulation of the ajoint equations for optimality. The purposes of this analysis are two-fol: first, as the founation for the efficient numerical solution of the problem by an all at once metho, an; secon, for comparison to the behavior of ajoints arising from more formal treatments, incluing those arising from shock-capturing proceures. Reynols Metal Professor, Aerospace Engineering Department, Associate Fellow AIAA y Assistant Professor, Mathematics Department z Grauate Stuent, Aerospace Engineering Department, Stuent Member AIAA The following notation will be use: = ensity u = velocity e = specific total energy h = height of the channel q = h x =h logarithmic erivative of channel height = specific heat ratio p = ( ; )(e ; u ) pressure: A quasi one imensional invisci steay state flow can be escribe by the steay state Euler equations in conservation form F(U(x)) where F(U) x q(x) Q(U(x)) = x () U U U U 3 A u e A u p u A Q(U) (e p)u u u A : (e p)u See e.g. Anerson. It shoul be note that, that F an Q are homogeneous of egree one, i.e. F(tU) =tf(u) an Q(tU) = tq(u) for all t IR. One can easily euce that F(U) =A(U) U an Q(U) =S(U) U () Copyright c996 by the authors. Publishe by the, Inc. with permission.
2 where an A(U) = F(U) U ;3 u (3 ; ) u ; ;ue ( ; ) u 3 e ; 3 ; u u S(U) = Q(U) U ;u u ;ue ( ; ) u 3 e ; 3 ; u u A A : It can be shown that A(U) has eigenvalues u ; a, u, u a, where a is the spee of soun. Bounary conitions for the system will be specifie below. Subsequently we will often omit the argument x. The problem we are intereste in is the esign of an area function h generating a flow that best approximates a esire pressure istribution p in the least squares sense. We assume that the area function h is fixe at some point, say x =. This allows us to use the logarithmic erivative q of h as the esign parameter. For given q(x) =h (x)=h(x) the area function can be compute from h(x) = exp ln(h()) Z x (3) (4) q(t) t : (5) The mathematical formulation of our esign problem is as follows: Minimize J(U q)= (p(x) ; p (x)) x (6) subject to the steay state Euler equations F(U) q Q(U) = x (7) x The flow is assume to be transonic, with the shock implicitly efine by the Rankine Hugoniot conition, F(U(x s ;)) = F(U(x s )): (8) We have bounary conitions at the inlet B i (U()) u e A in (u) in (e) in A = (9) an bounary conitions at the outlet B o (U) = p ; p out = () = ( ; ) U 3 ; U ; p U out : In aition, we impose bouns on the logarithmic erivative of the area function q low q(x) q upp x [ ]: () Iollo et al 6 have stuie this problem using an ajoint metho to formulate the optimality conitions. The ajoint system erive by them iffers from that in the current paper. The bounary conitions impose on the ajoint variables at the shock by Iollo et al 6 seem to be incorrect, except in a special case. This will be iscusse in greater etail in Section 4. An extension of their results to the two imensional case is consiere in Iollo an Salas. 7 Jameson 8 has also use an ajoint metho to formulate a similar class of optimal esign problems. A similar problem has been stuie by Frank an Shubin. 5 They point out that solutions of this simple moel exhibit phenomena quite similar to those in two imensional invisci flow over an airfoil. They have shown that uner given assumptions, the continuity an energy equations can be integrate out reucing the Euler system to a single ODE. In their version of the problem, a esire velocity profile is set as a target instea of a pressure istribution. They iscuss numerical solutions for this problem, employing stanar iscretization schemes, an using both a black box metho an a Sequential Quaratic Programming (SQP) metho. They have reporte that the SQP scheme is not very robust, as compare to the black box scheme, though when it oes converge, it proves to be the superior metho. Similar results, using the SQP metho, have been reporte in Sachikova et al. The failure of the SQP metho stems from the existence of the shock in the flow. In particular, when Gounov type schemes with low continuity properties are use to capture the shock, we fin that the SQP metho, which is esigne for smooth problems, performs poorly. In previous research 3 we iscusse an approach to overcome this problem by treating the shock location as an explicit state variable an enforcing the shock conition as an aitional state constraint. This was one in the context of the single ODE moel escribe by Frank an Shubin. 5 We foun that this approach yiels better smoothness properties for the constraints. We provie a rigorous mathematical analysis of the existence an
3 nature of optimal solutions for the infinite imensional problem, an stuie the relation between the infinite imensional an the iscrete problem, which reveale valuable insight into improving the performance of the optimization algorithm. Other aspects of this esign problem an other methos for its solution have been stuie by Borggaar, by Shenoy an Cliff an by Narucci, Grossman, an Haftka. In Borggaar the flow variables are viewe as functions of the esign variables an a sensitivity equation approach is use to compute the graient. An optimal control approach is use by Shenoy an Cliff. They assume that the shock location is known. The area an the flow velocity are the state variables. The esign problem is then formulate as an optimal control problem governe by a system of ODEs an solve using a multiple shooting metho for the two point bounary value problem which is obtaine from the necessary optimality conitions. In Narucci et al sensitivities for various (semi ) iscretizations are stuie. Although the shock location is treate implicitly, it enters the iscretization scheme for the objective function. The esign problem is solve numerically using the black box approach with a steepest escent metho.. The Design Problem in the Transforme Domain For a rigorous mathematical treatment of the optimal control problem we transform the esign problem presente in the introuction. In this transformation, the shock location is introuce as an explicit variable an the regions left an right of the shock are transforme onto the fixe interval [ ]: an [ x s ] ;! [ ] U(x) ;! U(x s ) [ x s ] ;! [ ] U(x) ;! U( ; ( ; x s )): we fin that the Euler equations (7), the Rankine Hugoniot shock conitions (8) an the bounary conitions (9), () are transforme into the following system of equations for U l an U r : F(U l ) x sq l Q(U l ) = [ ] F(U r ) (3) (x s ; )q r Q(U r ) = [ ] (4) with Rankine Hugoniot shock conition F(U l ()) = F(U r ()) (5) an bounary conitions B i (U l ()) = (6) at the inlet an with bounary conitions B o (U r ()) = (7) at the outlet. For convenience, we summarize these conitions in a composite constraint function C(U l U r x s q l q r )= F(U l ) x sq l Q(U l ) F(U r ) (x s ; )q r Q(U r ) F(U l ()) ; F(U r ()) B i (U l ()) B o (U r ()) : C A (8) With this formulation the system (3 7) can be written as: C(U l U r x s q l q r )=: The systems (7 ) an (3 7) are equivalent. If U is a solution of the original system (7 ) with shock location x s, then the solution of the transforme system (3 7) can be compute from (). On the other han, if (U l U r x s ) is a solution to the transforme system (3 7), then the solution U to the original system can be compute from 8 < x U l U(x) = : ;x U r x s x [ x s ) ;x s x [x s ]: (9) With To formulate the control problem for the transforme problem, we efine the esire pressures U l () = U(x s ) U r () = U( ; ( ; x s )) () p l ( x s) = p (x s ) p r ( x s) = p ( ; ( ; x s )): () 3
4 With the transformations (), () the objective function (6) is given by (p(x) ; p (x)) x = Z xs (p(x) ; p (x)) x x s (p(x) ; p (x)) x = x s (p l () ; p l ( x s)) ( ; x s ) (p r () ; p r ( x s)) : The control problem we have to solve is given as follows: Minimize J(U l U r x s q l q r ) = x s ; x s subject to the equality constraints (p l () ; p l ( x s)) () (p r () ; p r ( x s)) C(U l U r x s q l q r )= () an to the inequality constraints x s (3) q low q l () q r () q upp [ ]: (4) The states are given by the triple (U l U r x s ) an the controls are (q l q r ). The equations (3), (4), (5), (6), (7) are the state equations. The flow solutions (U l U r ) an the esign functions (q l q r ) o not epen on the shock location explicitly; they are couple to shock location x s only through the state equations (3), (4), (5). For a complete escription of the control problem, we have to specify the control space an the state space. It will be seen that an appropriate choice for the control space is Q = L ( ) L ( ): (5) The set of amissible parameters is efine by Q a = f(q l q r ) Qjq low q l () q r () q upp g (6) for [ ]. Recall that q low. A typical control vector is enote by q (q l q r ): (7) The state space is chosen to be V = W ( ) W ( ) IR (8) where W ( ) = (W ( )) 3 an a typical vector is enote as Defining the space v (U l U r x s ): (9) C = L ( ) L ( ) IR 7 (3) with L ( ) = (L ( )) 3, the constraint, (8), is viewe as a map from VQinto C: C : VQ!C: For later purposes we also efine the ajoint space = W ( ) W ( ) IR 7 : (3) As usual, for a function f :[ ]! IR we set kfk = ess sup [ ] jf ()j kfk = kfk kf k : For a vector value function F we efine kfk kfk analogously. 3. Differentiability an Surjectivity The transformation of the esign problem performe in Section allows us to establish the Fréchet ifferentiability of the constraint function C an of the objective function J. The proofs for Fréchet ifferentiability are very similar to those presente in Cliff et al 3 for the Frank an Shubin moel problem. In the following we make use of the notation (3), (4) an (9). Theorem 3. The nonlinear operator C : VQ!C is Fréchet ifferentiable at all (v q) VQ a with l r >. The partial Fréchet erivatives are given by C v (v q)(ev) = [A(U l )eu l ] x sq l S(U l )eu l ex sq l Q(U l ) [A(U r)eu r] (x s;)q rs(u r)eu rex sq rq(u r) A(U l ())eu l ();A(U r())eu r() U Bi (U l ())eu l () U Bo (U r())eu r() C A (3) 4
5 where U Bi (U l ()) = I 33 is equal to the ientity, ev = ( e U l e U r ex s ) an an U Bo (U r ()) = p r() U = ; ; u r () ;( ; )u r() ; C q (v q)(eq) = where eq =(eq l eq r ). x s Q(U l )eq l (x s ; )Q(U r )eq r (33) (34) C A Proof: The proof of this result is similar to the corresponing one in Cliff et al 3 an is therefore omitte. Theorem 3. If the esire pressure p is ifferentiable with absolutely continuous erivative, then the objective function J is Fréchet ifferentiable. The partial Fréchet erivatives are given by J Ul (v q) e U l = x s (p l ; p l ) p l U l eu l (35) J Ur (v q) e U r = (36) where p i (33), an ( ; x s ) U i =( ; (p r ; p r ) p r U r eu r u i ;( ; )u i ; ), i = l r, cf. J xs (v q) (37) = ; ; (p l() ; p l ( x s)) x s (p l () ; p l ( x s)) (p l ) x() (p r() ; p r ( x s)) ( ; x s )(p r () ; p r ( x s)) (p r ) x() : In the following, we interpret p r() U as a row or a column, as appropriate. The form shoul be apparent from the context. Proof: The assertion follows from the efinition of J using stanar estimates. The proof is therefore omitte. The final result of this section concerns the surjectivity of the linearize constraints. This property is important for the formulation of optimality conitions, see e.g. Maurer an Zowe. 9 Due to the partitioning of variables into states an controls, the surjectivity of the linearization of C is guarantee if C v (v q) is surjective. Recall that a l an a r enote the spee of soun in the regions left an right of the shock, respectively. Theorem 3.3 Let (U l U r x s ) be given with u l () >, u l () 6= a l (), u r () >, u r () 6= a r () for all [ ] an x s [ ]. If where A(U r ()) ;T p r() U T [T r q r Q(U r ) ;T r ()T l () ; T l q l Q(U l ) 6= T l () = exp T r() = exp ;R x s q l (t) S(U l (t))a(u l (t)) ; t ;R (x s;) q r(t) S(U ; r(t))a(u r(t)) t (38) then for all r =(r l r r r s r in r out ) Cthe linear system C v (v q)(ev) =r (39) amits a unique solution ev. Proof: Uner the assumptions of the theorem, the matrices A(U l ) an A(U r ) are invertible for all [ ]. Therefore, the first equation in (39) is equivalent to [A(U l ) e U l ] x s q l S(U l )A(U l ) ; A(U l ) e U l ex s q l Q(U l )=r l : Using the integrating factor T l the solution can shown to be eu l () =A(U l ()) ; T l () ; A(U l ()) e U l () Z T l () rl () ; ex s q l ()Q(U l ()) : (4) Similarly, the solution to the secon equation [A(U r ) U e r ] (x s ; )q r S(U r ) U e r ex s q r Q(U r )=r r 5
6 can be compute using the integrating factor T r an is given by eu r () =A(U r ()) ; T r () ; A(U r ()) e U r () Z T r () rr () ; ex s q r ()Q(U r ()) : (4) Inserting (4) an (4) into the thir equation A(U l ()) U e l () ; A(U r ()) U e r () = r s yiels ex s T r () ; T r ()q r ()Q(U r ()) ;T l () ; T l ()q l ()Q(U l ()) = r s ; T l () ; A(U l ()) e U l () T l ()r l () T r () ; A(U r ()) e U r () T r ()r r () : (4) Next, we have to incorporate the equations arising from the bounary conitions. The equation U Bi (U l ()) U e l () = r in implies U e l () = r in. If w w are linearly inepenent vectors orthogonal to p r() U, then U Bo (U l ()) U e r () = r out implies that eu r () = r out p r() k p r() U! w! w : U k Notice that, since 6=, the inequality k p r() U k > hols true. Inserting these expressions for U e l () an eu r () into (4) shows that the existence of a unique solution of (39), i.e. the existence of a unique solution (ex s!! ), can be guarantee if the 3 3 matrix [T r () ; T r q r Q(U r ) ;T l () ; T l q l Q(U l )] ;T r () ; A(U r ())w ;T r () ; A(U r ())w is nonsingular. This is the case if an only if A(U r ()) ; [T r q r Q(U r ) ;T r ()T l () ; T l q l Q(U l )] ;w ;w is nonsingular. Since w w are linearly inepenent vectors orthogonal to p r(), the latter conition is U equivalent to (38). We efine the vector an the Lagrange function L(v ) = x s 4. Optimality Conitions ( l r s i o ) ; x s (p l ; p l ) (p r ; p r ) T l [(F(U l)) x s q l Q(U l )] T r [(F(U r)) (x s ; ) q r Q(U r )] T s [F(U l()) ; F(U r ())] T i B i (U l ()) o B o (U r ()): (43) If the shock location at the optimum obeys x s ( ), then the first orer necessary optimality conitions are = L v (v q )(ev) (44) = L (v q )(e) (45) L q (v q )(eq) (46) for all ev C, for all e, an for all eq with (q eq) Q a. The secon equation (45) in the optimality system yiels the state system (3 7). Using integration by parts an appropriate test functions U e l U e r ex s we fin that the first equation (44) in the optimality system yiels the ajoint equations T A(U l ) ( l ) 6
7 h T = x s q l S(U l ) l x s (p l ; p pl l ) U l T A(U r ) ( r ) = (x s ; )q r S(U r ) T r h ( ; x s )(p r ; p pr i r ) U r with conitions i (47) (48) l () = ; s r () = s (49) at the shock, the bounary conitions T i o U Bi (U l ()) = l () T A(U l ()) U Bo (U r ()) = r () T A(U r ()) (5) an the conition q l T l Q(U l )q r T r Q(U r ) J xs (v q) =: (5) Using the efinitions (9) an () one can see that the first equation in (5) is equivalent to T A(U l ()) l () = i (5) which implies that the Lagrange multipliers l are free at the left en point =. If we enote the components of r by () r () r (3) r, then the secon equation in (5) is equivalent to ;3 u r ()() r () ;ur ()e r () ( ; ) u 3 r () (3) r () = o ; u r () () r ()(3; ) u r() () r () er () ;3 ; u r () (3) r () = ; o( ; )u r () ( ; ) () r () u r() (3) r () = o ( ; ) (53) The thir equation (46) in the optimality system is equivalent to 8 < if q l () =q low x s l () T Q(U l ) = if q l () (q low q upp ) : if q l () =q upp (54) an (x s ;) r () T Q(U r ) 8 < : if q r () =q low = if q r () (q low q upp ): if q r () =q upp : (55) Note the shock conition (49) for the co state. From the examination of the other equations it can be seen that l () = ; r () = can be guarantee only if the esire pressures can be matche exactly, i.e. the source terms in (47), (48), an the term J xs (U l U r x s q l q r ) in (5) vanish. This homogeneity of the co states at the shock was assume generally in Iollo et al. 6 However, the rigorous presentation shows that the ajoints ( l r ) will in general be nonzero at the shock with opposite signs left an right of the shock. The fact that the co states are nonzero at the shock if the optimal value of the objective function is nonzero can also be observe in the numerical experiments for the Frank an Shubin problem reporte in Cliff et al Conclusions an Outlook We have provie a careful erivation of an optimality system for a shape-esign problem base on -D Euler flows. The analysis shows that the ajoints are smooth (inee in W ) in regions where the flow is smooth, but just as with the flow variables, there is a shock in the ajoint variables. The conitions erive in Section 4 can be use for numerical implementation in several ways. The statemoel (3 7) efines a bounary-value problem to be solve for the flow variables an the shock-location given the area-istribution. With the state/control in-han, the ajoint system (47 49, 5, 53) efines the ajoint variables ( l r ). Finally, the optimality conitions are given as (54,55). A simple graient metho woul procee by solving the first of these for the state, the secon for the ajoint an then using an increment in the control as erive from the optimality conitions. As a globalization strategy, one might use the control increment as a search irection in the control space Q an perform a line-search for improvement of the cost function. This woul require solving the system () for the variables (U l U r x s ) for each trial control increment. As an alternative approach one might use an SQP framework, wherein the states an controls are treate as inepenent variables an the state moel () is treate as an explicit constraint. In this all-at-once framework iterates procee towar optimality an feasibility at the same time. The rigorous analysis of ifferentiability an 7
8 existence of Lagrange multipliers is expecte to be important here, as it was in Cliff et al. 3 The stuy in Frank an Shubin 5 showe that formal evelopment of an SQP algorithm i not lea to a satisfactory computational proceure. This is because the shock-capturing flow-moel use oes not possess the necessary smoothness properties to justify Lagrangian formulation. In Cliff et al 3 we utilize the same shock-fitting formulation employe here, but for the single ODE moel of the flow. 5 In the present paper we have applie this approach to the full -D Euler system. Aitionally, it is worth noting that care is neee in the formulation of the finite-imensional approximation to preserve consistency with the continous problem. In Cliff et al 3 we observe that a moification of the finite-imensional state/ajoint inner-prouct ha a substantial impact on numerical efficiency of the overall proceure. The analysis leas to a characterization of optimal shapes an can be numerically implemente in several ways. The flow formulation can be interprete as a shock-fitting scheme an so is of questionable practical utiltity as a computational proceure for 3-imensional flows. On the other han, the analysis clearly reveals the nature of a shock in the ajoint system an this is of value for comparison to shock-capturing methos. There is no reason to believe that a iscretization that leas to accurate shock moeling for the flow system is necessarily appropriate for moeling the shock in the ajoint system. The results in this paper can be use to construct test examples since we have a characterization of the ajoint system without numerical artifacts, such as artificial viscosity. Acknowlegement This research was supporte by the Air Force Office of Scientific Research uner grant F an by the NSF uner Grant DMS References [] J. D. Anerson. Moern Compressible Flow with Historical Prespective. McGraw Hill Series in Aeronautical an Aerospace Engineering. McGraw Hill, New York, 99. [] Borggaar, J. T., The Sensitivity Equation Metho for Optimal Design, PhD thesis, Virginia Polytechnic Institute an State University, Department of Mathematics, Blacksburg, VA 46 3, USA, 994. [3] Cliff, E.M., Heinkenschloss, M., an Shenoy, A., Optimal Control for Flows with Discontinuities, Technical Report 95 9, Interisciplinary Center for Applie Mathematics, Virginia Polytechnic Institute an State University, Blacksburg, VA 46, 995, also Journal of Optimization, Theory an Application, to appear. [4] Dennis, J. E., Heinkenschloss, M., an Vicente, L. N., Trust region interior point algorithms for a class of nonlinear programming problems, Technical Report 94, Interisciplinary Center for Applie Mathematics, Virginia Polytechnic Institute an State University, Blacksburg, VA 46, 994. [5] Frank, P. D. an Shubin, G. Y., A comparison of optimization base approaches for a moel computational aeroynamics esign problem, J. Comput. Physics, vol. 98, 99, pp [6] Iollo, A., Salas, M. D. an Ta asan, S., Shape Optimization governe by the Euler Equations using an Ajoint Metho, ICASE Report No , also NASA CR 9555, November 993. [7] Iollo, A. an Salas, M. D., Contribution to the Optimal Shape Design of Two-Dimensional Internal Flows with Embee Shocks, ICASE Report No. 95-, also NASA CR 956, March 995. [8] Jameson, A., Optimum Aeroynamic Design Using CFD an Control Theory, In AIAA Paper 95-79, presente at the th Computational Flui Dynamics Conference, San Diego, California, June 995. [9] Maurer, H. an Zowe, J. First an Secon orer necessary an sufficient optimality conitions for infinite imensional programming problems, Mathematical Programming, vol. 6, 979, pp 98. [] Narucci, R., Grossman, B., an Haftka, R. T., Sensitivity algorithms for an inverse esign problem involving a shock wave, In AIAA Paper 94-96, presente at the 3n Aerospace Sciences Meeting & Exhibit, Reno, Nevaa, January -3, 994. [] Sachikova, E., Shenoy, A. an Cliff, E.M., Computational Issues in Optimization Base Design, In Proceeings from the 34th IEEE Conference on Decision an Control, New Orleans, December 995, pp
9 [] Shenoy, A. an Cliff, E.M., An optimal control formulation for a flow matching problem, In Proceeings from the 5th AIAA/USAF/NASA/ISSMO Symposium On Multiisciplinray Analysis An Optimization, Panama City Beach, September 7-9, 994, pp
Problems 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 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 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 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 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 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 informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More 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 informationTMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments
Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary
More informationCalculus of Variations
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 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 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 information7.1 Support Vector Machine
67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to
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 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 informationCOUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS
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
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 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 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 informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationTRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS
TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS Francesco Bullo Richar M. Murray Control an Dynamical Systems California Institute of Technology Pasaena, CA 91125 Fax : + 1-818-796-8914 email
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 informationOptimal Variable-Structure Control Tracking of Spacecraft Maneuvers
Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign
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 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 informationLower bounds on Locality Sensitive Hashing
Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More 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 informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationDiscrete Hamilton Jacobi Theory and Discrete Optimal Control
49th IEEE Conference on Decision an Control December 15-17, 2010 Hilton Atlanta Hotel, Atlanta, GA, USA Discrete Hamilton Jacobi Theory an Discrete Optimal Control Tomoi Ohsawa, Anthony M. Bloch, an Melvin
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 informationSINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
Communications on Stochastic Analysis Vol. 2, No. 2 (28) 289-36 Serials Publications www.serialspublications.com SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
More informationON ISENTROPIC APPROXIMATIONS FOR COMPRESSIBLE EULER EQUATIONS
ON ISENTROPIC APPROXIMATIONS FOR COMPRESSILE EULER EQUATIONS JUNXIONG JIA AND RONGHUA PAN Abstract. In this paper, we first generalize the classical results on Cauchy problem for positive symmetric quasilinear
More informationApplications of the Wronskian to ordinary linear differential equations
Physics 116C Fall 2011 Applications of the Wronskian to orinary linear ifferential equations Consier a of n continuous functions y i (x) [i = 1,2,3,...,n], each of which is ifferentiable at least n times.
More informationOptimization of a point-mass walking model using direct collocation and sequential quadratic programming
Optimization of a point-mass walking moel using irect collocation an sequential quaratic programming Chris Dembia June 5, 5 Telescoping actuator y Stance leg Point-mass boy m (x,y) Swing leg x Leg uring
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 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 informationProblem set 2: Solutions Math 207B, Winter 2016
Problem set : Solutions Math 07B, Winter 016 1. A particle of mass m with position x(t) at time t has potential energy V ( x) an kinetic energy T = 1 m x t. The action of the particle over times t t 1
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 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 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 informationSYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is
SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. Uniqueness for solutions of ifferential equations. We consier the system of ifferential equations given by x = v( x), () t with a given initial conition
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 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 informationCONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS
BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN
More 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 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 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 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 informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
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 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 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 informationLower Bounds for the Smoothed Number of Pareto optimal Solutions
Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.
More information2 GUANGYU LI AND FABIO A. MILNER The coefficient a will be assume to be positive, boune, boune away from zero, an inepenent of t; c will be assume con
A MIXED FINITE ELEMENT METHOD FOR A THIRD ORDER PARTIAL DIFFERENTIAL EQUATION G. Li 1 an F. A. Milner 2 A mixe finite element metho is escribe for a thir orer partial ifferential equation. The metho can
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 informationConservation Laws. Chapter Conservation of Energy
20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action
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 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 informationCalculus of Variations
Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,
More informationA. Incorrect! The letter t does not appear in the expression of the given integral
AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question
More 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 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 informationOptimization of Geometries by Energy Minimization
Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.
More informationA Universal Model for Bingham Fluids with Two Characteristic Yield Stresses
A Universal Moel for Bingham Fluis with Two Characteristic Yiel Stresses NGDomostroeva 1 an NNTrunov DIMeneleyev Institute for Metrology Russia, StPetersburg 195 Moskovsky pr 19 February, 4, 9 Abstract:
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 informationLecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012
CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration
More informationAn Optimal Algorithm for Bandit and Zero-Order Convex Optimization with Two-Point Feedback
Journal of Machine Learning Research 8 07) - Submitte /6; Publishe 5/7 An Optimal Algorithm for Banit an Zero-Orer Convex Optimization with wo-point Feeback Oha Shamir Department of Computer Science an
More informationStructure preserving finite volume methods for the shallow water equations
Structure preserving finite volume methos for the shallow water equations by Ulrik Skre Fjorholm THESIS for the egree of MASTER OF SCIENCE Master i Anvent matematikk og mekanikk) Department of Mathematics
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 informationTractability results for weighted Banach spaces of smooth functions
Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More 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 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 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 informationEnergy behaviour of the Boris method for charged-particle dynamics
Version of 25 April 218 Energy behaviour of the Boris metho for charge-particle ynamics Ernst Hairer 1, Christian Lubich 2 Abstract The Boris algorithm is a wiely use numerical integrator for the motion
More informationCalculus and optimization
Calculus an optimization These notes essentially correspon to mathematical appenix 2 in the text. 1 Functions of a single variable Now that we have e ne functions we turn our attention to calculus. A function
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationGeneralized Tractability for Multivariate Problems
Generalize Tractability for Multivariate Problems Part II: Linear Tensor Prouct Problems, Linear Information, an Unrestricte Tractability Michael Gnewuch Department of Computer Science, University of Kiel,
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 informationA Weak First Digit Law for a Class of Sequences
International Mathematical Forum, Vol. 11, 2016, no. 15, 67-702 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1288/imf.2016.6562 A Weak First Digit Law for a Class of Sequences M. A. Nyblom School of
More informationFLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction
FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS ALINA BUCUR, CHANTAL DAVID, BROOKE FEIGON, MATILDE LALÍN 1 Introuction In this note, we stuy the fluctuations in the number
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 informationNested Saturation with Guaranteed Real Poles 1
Neste Saturation with Guarantee Real Poles Eric N Johnson 2 an Suresh K Kannan 3 School of Aerospace Engineering Georgia Institute of Technology, Atlanta, GA 3332 Abstract The global stabilization of asymptotically
More informationProof of SPNs as Mixture of Trees
A Proof of SPNs as Mixture of Trees Theorem 1. If T is an inuce SPN from a complete an ecomposable SPN S, then T is a tree that is complete an ecomposable. Proof. Argue by contraiction that T is not a
More 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 informationarxiv:hep-th/ v1 3 Feb 1993
NBI-HE-9-89 PAR LPTHE 9-49 FTUAM 9-44 November 99 Matrix moel calculations beyon the spherical limit arxiv:hep-th/93004v 3 Feb 993 J. Ambjørn The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø,
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 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 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 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 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 information3 The variational formulation of elliptic PDEs
Chapter 3 The variational formulation of elliptic PDEs We now begin the theoretical stuy of elliptic partial ifferential equations an bounary value problems. We will focus on one approach, which is calle
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 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 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 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 informationChapter 9 Method of Weighted Residuals
Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe
More information