To appear in International Journal of Numerical Methods in Fluids in Stability analysis of numerical interface conditions in uid-structure therm
|
|
- MargaretMargaret Chambers
- 5 years ago
- Views:
Transcription
1 To appear in International Journal of Numerical Method in Fluid in 997. Stability analyi of numerical interface condition in uid-tructure thermal analyi M. B. Gile Oxford Univerity Computing Laboratory Numerical Analyi Group Thi paper analye the numerical tability of coupling procedure in modelling the thermal diuion in a olid and uid with continuity of temperature and heat ux at the interface. A imple onedimenional model i employed with uniform material propertie and grid denity in each domain. A number of dierent explicit and implicit algorithm are conidered for both the interior equation and the boundary condition. The analyi how that, in general, thee are table provided Dirichlet boundary condition are impoed on the uid and Neumann boundary condition are impoed on the olid; in each cae, the impoed value are obtained from the other domain. Oxford Univerity Computing Laboratory Numerical Analyi Group Wolfon Building Park Road Oxford, England OX 3QD April, 997
2 Introduction Thi analyi i motivated by interet in numerical procedure for coupling eparate computation of thermal diuion in a olid and a uid. A typical example application i the computation of heat tranfer to a blade in a ga turbine. The urrounding air in a high preure turbine i on average at a much higher temperature and therefore there i a ignicant heat ux from the uid into the turbine blade. In teady-tate, thi i matched by a correponding heat tranfer from the blade to relatively cold air owing through internal cooling paage. One approach to the numerical approximation of the above ituation would be the ue of a ingle conitent, fully-coupled dicretiation modelling both the olid and the uid, plu the boundary condition at the interface [9]. However, for the olid it i the calar unteady parabolic p.d.e. which decribe the thermal diuion, while for the uid the appropriate equation are the Navier- Stoke equation, with uitable turbulence modelling. Therefore, the production of a ingle fully-coupled code for the combined diuion application can be a much work a writing the individual program for the eparate olid and uid application. Since there are often exiting code which accurately and eciently olve thee individual problem, a more practical approach in many circumtance i to ue thee together to analye coupled problem [, 4, 8, 7, 3, ]. Both CFD code and thermal analyi code uually have the capability to pecify either the temperature or the heat ux at boundarie. A natural choice therefore for coupling thee code i to pecify the urface temperature at the interface in one code, taking the value from the other code, and pecify the boundary heat ux in the econd code, taking it value from the rt code [, 4]. A concern i whether there i any poibility that thi coupling procedure will introduce a numerical intability which doe not exit for the uncoupled problem. Thi i the iue that i addreed in thi tudy. The general theory for the analyi of numerical interface or boundary condition intabilitie i well-etablihed [6, ] but can be very complicated to apply in practice. In D and 3D computation of engineering interet, one cla of error mode which might be untable are thoe whoe variation i purely in the direction normal to the interface between the olid and uid. The nite difference or nite volume equation for thi cla of error mode reduce to being one-dimenional, and therefore in thi paper we implify the analyi for thi diuion problem by retricting attention to a imple D model problem with a uniform grid on either ide of the interface. Since there i no velocity component normal to the olid boundary in D and 3D ow, it i appropriate in the D model problem to omit any convection term. Stability for thi D model problem i a neceary condition for the tability of the real D and 3D computation. It may, or may not, be a ucient condition for tability but undertanding the nature of poible D intabilitie clearly give inight into the potential intabilitie
3 3 in D and 3D computation. Analytic problem The parabolic p.d.e. decribing unteady thermal diuion @x q = (.) Here T (x; t) i the temperature, q(x; t) i the heat ux, c(x) i the heat capacity and k(x) i the conductivity. Thee equation are valid for arbitrary, piecewie continuou poitive function c(x); k(x). The nite volume algorithm to be analyed are all baed on the integral verion of thi equation, d dt Z x x c T dx = [q] x x : (.) At any interface at which c and/or k are dicontinuou, the equation are augmented by the requirement that T and q mut be continuou. The boundary condition a x! are that the temperature aymptote to a contant value, T, and o the heat ux tend to zero. Dening T max (t) = max x T (x; t); T min (t) = min x T (x; t); (.3) an important property of olution of the unteady diuion equation i that for any non-uniform initial condition T (x; 0) and for all t > 0, dt max dt 0; dt min dt 0: (.4) Furthermore, if T = T + then T min (t)! T max (t) a t!. The behaviour of the maximum and minimum temperature will be important in dening the numerical tability of the coupled ytem. Although the above theory i given for general c(x); k(x), in thi paper we will now retrict attention to a ingle interface at x = 0 with c and k having uniform value c ; k for x <0, and c + ; k + for x >0. 3 Fully-coupled dicretiation In thi ection we examine the tability of fully-coupled dicretiation of the model problem. The theory for thi i well etablihed ince it i imply a pecial cae of the more general problem of the dicretiation of a parabolic p.d.e. with patially varying diuivity [5, 0]. There are everal reaon for doing thi analyi even though it i believed that the fully-coupled approach i not the mot
4 4 practical approach to real coupled application. The rt i to how that in a good fully-coupled dicretiation there are no intabilitie aociated with the interface treatment, that tability of the dicretiation on the uniform meh on either ide of the interface i a neceary and ucient for tability of the fullycoupled dicretiation. The econd i to have a benchmark againt which to compare the `weakly-coupled' dicretiation in the next ection. Thee will be hown to have interface intabilitie under certain condition, and it i informative to ee how thee are related to dierence in the interface treatment relative to the fully-coupled dicretiation. Uing a computational grid with uniform pacing x for x <0 and uniform pacing x + for x >0, the location of grid node i given by x = ( x ; 0 x + ; 0 (3.) Aociated with each grid node i the dicrete temperature variable T n to approximate the analytic olution T (x; t) at x=x ; t=nt. which i 3. An explicit algorithm Uing forward Euler time dierencing and conervative patial dierencing baed on the integral form of the unteady diuion equation on the interval x xx + give the following explicit algorithm, C (T n+ T n ) = (q n + q n ); q n + = K + (T+ T n n ) (3.) where 8 >< C = c x t c x t ; < 0 + c +x + ; = 0 t (3.3) and >: K + c + x + ; > 0 t = 8 >< >: k x ; + < 0 k + x + ; + > 0 (3.4) Note that the equation for = 0 involve the conductivity and heat capacity on both ide of the interface. In particular, t C 0 i the heat capacity of the whole nite volume computational cell extending from x to x +.
5 5 For 6=0, the dierence equation reduce to T n+ = T n + d T n + T n + T n (3.5) where d = k t cx : (3.6) Standard Fourier tability analyi on either ide of the interface how that a dicrete Fourier mode i table provided d. We will now prove that if the requirement of Fourier tability are atied on each ide of the interface, then the fully-coupled dicretiation i table in the ene that T n+ max T n max ; T n+ min T n min ; (3.7) where T n max max T n ; We begin by noting that if d T n min min T n : (3.8) then for any poitive value r d + rd + ( + r) =) d + rd + + r : (3.9) We will ue thi reult with r dened a the ratio of the heat capacitie of the computational cell on either ide of the interface, where r = c +x + c x : (3.0) The next tep i to re-write the full dierence equation a T n+ = ( a b ) T n + a T n + + b T n (3.) a = b = d ; < 0 a 0 = d + r ; b 0 = rd + + r ; (3.) a = b = d + ; > 0 0 < d o for all, a ; b and a b are poitive quantitie and thu T n+ i a poitive weighted average of T n + ; T n ; T n. Hence, T n min min(t n + ; T n ; T n ) T n+ max(t n + ; T n ; T n ) T n max (3.3) Thi i true for all, and o taking the maximum over all, and the minimum over all, give the deired reult, Equation (3.7).
6 6 3. An implicit algorithm Replacing the forward Euler time dierencing with backward Euler time dierencing give the following implicit algorithm. C (T n+ T n ) = (q n+ + q n+ ) q n+ + = K + (T n+ + T n+ ); (3.4) with C ; K + a dened before. Fourier tability analyi of the dicretiation on either ide of the interface how it to be unconditionally table. The fully-coupled dicretiation i alo unconditionally table in the ame ene a before. To prove thi, the dierence equation i re-written a where T n+ = ( a b ) T n + a T n+ + + b T n+ ; (3.5) a = b = K + K + + K K + C ; ; (3.6) K + + K + C a b = C K + + K + C : It i clear that a ; b and a b are poitive quantitie, for all. We now chooe J uch that T n+ J = Tmax. n+ Subtracting Tmax n+ from both ide of the dierence equation give ( a J b J )(T n J T n+ max) + a J (T n+ J+ T n+ max) + b J (T n+ J T n+ max) = 0: (3.7) Becaue T n+ J ; T n+ J+ Tmax n+ and a J ; b J ; a J b J are all poitive, either TJ n > T max n+ or TJ n = T n+ J+ = T n+ J+ = Tmax. n+ In the rt cae, we immediately get the reult that Tmax n > T n+ max. In the econd cae, we can repeat the argument with = J. By further repetition if neceary, we conclude that either Tmax n > T max n+ or T n = T n+ = Tmax, n+ for all, in which cae Tmax n = Tmax. n+ Exactly the ame argument can be ued to prove that Tmin n T n+ min equality occurring only in the trivial cae in which T n i contant. with 4 Looely-coupled dicretiation In the looely-coupled dicretiation, each half of the domain i olved eparately with boundary condition containing information from the other. The natural boundary condition for a diuion problem are either Dirichlet (the pecication
7 7 of the boundary temperature) or Neumann (the pecication of the boundary heat ux). Therefore we will conider a looely-coupled procedure in which the calculation for x 0 ue Dirichlet data obtained from the olution for x 0, while the calculation for x 0 ue Neumann data obtained from the olution for x An explicit algorithm Given exiting olution at time level n in both halve of the domain, the implet and mot natural explicit numerical algorithm for determining T n+ for 0 i c x t c x t (T n+ T n ) = k x (T n + T n +T n ); < 0 (T n+ 0 T n 0 ) = q w k x (T n 0 T n ); (4.) where q w i the heat ux pecied a the interface boundary condition. Uing a nite volume derivation, the equation for < 0 correpond to the control volume [x ; x + ] of width x, wherea the equation for =0 correpond to the control volume [x ; 0] of width x. The implet conitent equation for determining the heat ux at the interface from the data in 0 i q w = k + x + (T n T n 0 ): (4.) Thi one-ided approximation to the temperature gradient at the urface i only rt order accurate during unteady tranient. However, it i typical of the numerical method ued for practical computation [8, 7]. The correponding explicit numerical algorithm for imultaneouly determining T n+ for >0 i c + x + t (T n+ T n ) = k + x + (T n + T n +T n ): (4.3) The equation for = require the variable T n 0 and thi i et by the Dirichlet boundary condition T n 0 = T w ; (4.4) where T w i the interface temperature. The obviou value for thi i imply T n 0 from the computation for 0. To ummarie the communication between the two calculation for 0 and 0, at each timetep there i an exchange of data, with the program or ubroutine performing the calculation for 0 upplying the value of T w to the other program or ubroutine performing the calculation for 0, while the latter
8 8 end q w to the former. It i then poible that the computation for the two halve could proceed in parallel (perhap uing eparate procee on eparate worktation) until they again exchange data before the next timetep. By comparing Equation (4.,4.3) with Equation (3.), it can be een that the only dierence i the omiion of the term c + x + t in the equation for = 0. If c + x + c x, then thi omitted term i negligible compared to the retained term c x t and o it eem likely that no intability will be introduced by it omiion. On the other hand, if c + x + c x, then the omitted term may be very ignicant. Thi indicate very imply that a key parameter in the following analyi will be the variable r, dened earlier in Equation (3.0) a the ratio of thee two quantitie. For the purpoe of analyi it i more convenient to conolidate and implify the equation into the following form, T n+ = T n + d T n + T n + T n ; < 0 T n+ 0 = T n 0 d T n 0 T n + rd+ (T n T n 0 ) ; (4.5) T n+ = T n + d + T n + T n + T n ; > 0 where d and r are a dened previouly. In applying the tability theory of Godunov and Ryabenkii [6, ], the tak i to invetigate the exitence of eparable normal mode of the form T n = z n f : (4.6) The dicretiation i untable if the dierence equation admit uch olution which atify the far-eld boundary condition, f! 0 a!, and have z >, giving exponential growth in time. The form of the olution i very imilar to the aumed Fourier mode, except that the amplitude of the patial ocillation decay exponentially with away from the interface. For thi application the normal mode mut be of the form 8 < T n = : z n ; 0 z n +; 0 : (4.7)
9 9 The dierence equation, Equation (4.5) are atied provided the three variable z; ; + atify the following equation. z = + d ( + ) z = + d ( ) + rd + ( + ) (4.8) z = + d + ( ) Solving the rt of thee equation to obtain give 0 = 4d d z A : (4.9) To atify the far-eld boundary condition a! it i neceary to chooe the negative quare root when the argument i real and poitive. When it i complex, the choice of root i dened by the requirement that <. Similarly, olving the third of the equation give + = z d + 4d + A : (4.0) z Subtituting thee into the econd equation give the following nonlinear equation for z. 0 4d 4d + A = 0 (4.) z z There i no imple cloed form olution to thi, giving z a an explicit function of the parameter d ; d + ; r. Intead, we conider aymptotic olution under dierent aumption. When d ; d +, the quare root term can be expanded to give the following approximate equation and olution. d z rd + z 0; =) z d rd + : (4.) The requirement for tability i z <. The olution z(r) lie inide z = for uciently mall value of r, but then croe it at z = when r = d +. Thu for d ; d + the tability requirement i r < d +. Expanding the analyi to conider arbitrary value for d ; d +, we begin by conidering the aymptotic behaviour when r and r. When r, the econd term in Equation (4.) i relatively mall, and the approximate olution i 4d z 0; =) z 4d : (4.3)
10 0 Since d mut atify 0 < d for the dicretiation to be table according to Fourier tability analyi, it follow that z. Thu, there i no coupled intability when r. When r, the rt term in Equation (4.) i relatively mall, and o to a rt approximation the olution i 4d + z 0; =) 4d + z 0 =) z : (4.4) To get a more accurate approximate olution, the rt term i approximated uing z to obtain 4d + z r =) 4d + z r =) z rd + : (4.5) Thu for xed d + and uciently large r, there i an intability with z being large, real and negative. The correponding value of and + will be mall, real and negative, o the intability will appear a a `awtooth' ocillation mode, both patially and in time, with an amplitude which decay exponentially away from the interface, but grow exponentially in time. Since the looely coupled ytem i table for r and untable for r the remaining quetion i the value of r at which the intability begin. Thi correpond to the lowet poitive real value of r for which z =. Becaue of the requirement that r i real, it can be hown from Equation (4.) that thi again require z =, in which cae p d r = p (4.6) d + Thu the condition for tability i r < p d p d + (4.7) A typical calculation with a timetep cloe to the Fourier tability limit might have d = d + = 3, for which the coupled tability limit i r <. The key to 8 obtaining tability in practical computation i the correct choice of which half of the domain ue the Dirichlet boundary condition and which half ue the Neumann boundary condition. The uual practice for the coupled blade/air computation dicued in the Introduction i to ue Neumann boundary condition for the olid computation, and Dirichlet boundary condition for the uid computation. For thi choice, the correponding value of r i given by r = c uid x uid c olid x : (4.8) olid
11 Given typical value for the parameter involved, r i uually very mall and o thi i table. If, on the other hand, one were to ue Dirichlet boundary condition for the olid computation and Neumann boundary condition for the uid computation, then the appropriate value for r would be the invere of the above quantity, which would be very large. In thi cae the coupled calculation would be untable unle one ued an extremely mall timetep. Uing the approximate olution for r in Equation (4.5), the timetep tability limit i given by d + r ; (4.9) o tability of the coupled ytem would require the ue of a timetep very much maller than that needed for Fourier tability. Thi analyi i upported by the numerical reult preented in Figure and. The computation ue the nite domain , initial condition T 0 = for < 0 and T 0 = for 0 and boundary condition T n 000 = ; T n 000 =. In addition, all of the computation ue d = d + = 3 for 8 which the analyi above predict the coupled ytem to be table only for r <. Figure how two et of reult with T n plotted for the rt 0 iteration in each cae. In a), r = and the olution i clearly table, with an initial tranient at the interface decaying very quickly, while in b), r = and the olution i very untable. Figure how another two et of reult with T n plotted every 5 iteration. In a), r = 0:99 and the olution appear to be table, although with the interface tranient decaying more lowly in thi cae, while in b), r = :0 and the olution i clearly untable. 4. A hybrid algorithm The next algorithm to conider i a hybrid one, in which the computation i unaltered for >0, but the algorithm for 0 i replaced by the correponding implicit method baed on a backward Euler time dicretiation. c x t c x t (T n+ T n ) = k n+ (T+ T x n+ +T n+ ); < 0 (T n+ 0 T n 0 ) = q w k n+ (T0 T x n+ ): (4.0) The boundary heat ux q w i again dened explicitly by q w = k + x + (T n T n 0 ): (4.) The dierence equation for >0 are unchanged, a i the communication of data between the calculation for 0 and >0.
12 The conolidated, implied form of the equation i T n+ = T n + d T n+ + T n+ T n+ 0 = T n 0 d T n+ 0 T n+ T n+ = T n + d + T n + T n + T n and the normal mode i again of the form 8 < T n = : + T n+ ; < 0 z n ; 0 + rd+ (T n T n 0 ) ; (4.) ; > 0 z n +; 0 : (4.3) The dierence equation, Equation (4.) are atied provided the three variable z; ; + atify the following equation. = z + d ( + ) = z + d ( ) + rd + z ( + ) (4.4) z = + d + ( ) The third of thee equation require that + depend on z in exactly the ame way a for the purely explicit algorithm. Solving the rt of thee equation ubect to the far-eld boundary condition give = + z d Subtituting thee into the econd equation give + 4d z r When r, the aymptotic olution i and o the dicretiation i table for all value of d. When r, the aymptotic olution i + 4d A : (4.5) z 4d + A = 0 (4.6) z z = ( + 4d ) + O(r); (4.7) z rd + p + 4d + O(r ); (4.8) and o the coupled dicretiation i till untable for uciently large value of r.
13 3 The cro-over from tability to intability again occur when z =, giving p + d r = p (4.9) d + Thu the condition for tability i r < p + d p d + (4.30) Comparing thi reult with the correponding reult for the purely explicit algorithm, it can be een that the new tability region i greater except when d. Thi ha a phyical interpretation; when d i not mall, the trong implicit coupling of the computational cell for 0 increae the eective thermal capacity of the cell aected in one timetep by the interface heat ux. Numerical experiment were performed on the ame domain and with the ame initial and boundary condition a before, and with d = 4; d + = 3 8. The analyi above predict tability provided r < 6, and thi i upported by Figure 3 which how two et of reult with T n plotted every 5 iteration. In a), r = 5:95 and the olution i table with a lowly decaying interface tranient, while in b), r =6:05 and the olution i clearly untable. 4.3 An implicit algorithm We now conider an algorithm which i implicit on each ide of the interface, but with explicit updating of the data ued for the interface boundary condition. The implicit numerical algorithm for 0 i again c x t c x t with q w dened explicitly by (T n+ T n ) = k n+ (T+ T x n+ +T n+ ); < 0 (T0 n+ T0 ) n = q w k n+ (T x 0 T n+ ): (4.3) q w = k + x + (T n T n 0+): (4.3) An important point in the above equation i the ditinction between T0, n the value of T n at =0 a calculated for the domain 0, and T0+, n the value of T n at = 0 for the domain 0. In the previou dicretiation thee two value have been identical but thi will not be true in thi cae. The correponding implicit numerical algorithm for imultaneouly determining T n+ for >0 i c + x + t (T n+ T n ) = k + x + (T n+ + T n+ +T n+ ): (4.33)
14 4 The equation for = require the variable T n+ 0+ and thi i et by the Dirichlet boundary condition T n+ 0+ = T w ; (4.34) where T w i the interface temperature. Uing explicit updating of boundary data, T w = T n 0 ; (4.35) p T 0+ lag T 0 by one iteration. The pattern of communication between the calculation for 0 and 0 i exactly the ame a for the explicit algorithm. They exchange the value of T w and q w at the beginning of the timetep, perform the timetep calculation independently (poibly in parallel on eparate worktation) and then repeat the proce for the next timetep. For the purpoe of analyi it i again more convenient to conolidate and implify the equation into the following form, T n+ = T n + d T n+ + T n+ + T n+ T0 n+ = T0 n d T n+ 0 T n+ T n+ = T n + d + T n+ + T n+ + T n+ T n+ 0+ = T n 0 : ; < 0 + rd+ T n T n 0+ ; > 0 The form of the normal mode olution for thi cae i 8 < T n = : z n ; = 0 ; ; ; 3; : : : z n +; = 0+; ; ; 3; : : : ; (4.36) : (4.37) The fourth equation in Equation (4.36) i automatically atied by the above choice of normal mode. The other three equation require that the variable z; ; + atify the following equation. = z + d ( + ) = z + d ( ) + rd + z ( + ) (4.38) = z + d + ( ) Solution of the rt and third of thee equation, ubect to the far-eld boundary condition, give = + z d + = + z d d z A ; + 4d + A : (4.39) z
15 5 Subtituting thee into the econd equation give + 4d z + rz 0 When r, the aymptotic olution i and o the dicretiation i table for all value of d. When r, the aymptotic olution + 4d + A = 0: (4.40) z z = ( + 4d ) + O(r); (4.4) z i p r p! + 4d+ p + O(): (4.4) + 4d Thu for xed d ; d + and uciently large r, the coupled ytem i untable. It i not poible for general value of d ; d + to determine explicitly the value of r above which the olution procedure i untable. It i poible however to obtain an aymptotic olution under the aumption d ; d +. Thi i a reaonable aumption ince the motivation in uing implicit method i to ue much larger timetep than would be table uing explicit method. Under the aumption d ; d +, Equation (4.40) reduce to qd + rz q d + 0; =) z i p r d + d! 4 : (4.43) Hence, under thee condition the approximate tability limit i which can alo be re-expreed a r < c 3 +x 4 + k + d d + ; (4.44) < c 3 x 4 k : (4.45) Provided, a before, that the correct choice i made a to which domain ue the Neumann b.c.' and which ue the Dirichlet b.c.', then r hould be uciently mall that practical computation will be table. Numerical experiment were performed on the ame domain and with the ame initial and boundary condition a before, and with d = d + = 50. The approximate tability limit for thee value i r <, and thi i upported by Figure 4 which how two et of reult with T n plotted every 5 iteration. In a), r = : and the olution i table with a lowly decaying interface tranient, while in b), r =: and the olution i clearly untable.
16 6 5 Concluding remark The tability analyi in thi paper ha hown the viability of a looely-coupled approach to computing the temperature and heat ux in coupled uid/tructure interaction. The key point to achieving numerical tability i the ue of Neumann boundary condition for the tructural calculation and Dirichlet boundary condition for the uid calculation. Although the analyi wa performed here for the D model diuion equation, the reult are believed to be applicable to the real ituation in which the 3D diuion equation i ued to model the heat ux in the tructure and the 3D Navier-Stoke equation are ued to model the behaviour of the uid. Thi i upported by the practical experience of 3D computation performed uing thi coupling procedure [, 4]. The analyi alo aumed a time-accurate modelling of the uid/tructure interaction. In practical computation, the point of engineering interet i often the teady-tate temperature and heat ux ditribution. In uch cae, the computation in the tructure and uid can both proceed with dierent timetep given by their repective Fourier tability limit. The coupled normal mode analye remain valid uing the value of d ; d + baed on the timetep t ; t + ued in the two domain. Acknowledgement Thi reearch wa upported by Roll-Royce plc and ha benetted from dicuion with Dr. Peter Stow of Roll-Royce plc and Dr. Mehmet Imregun of Imperial College. Reference [] R.S. Amano, K.D. Wang, and V. Pavelic. A tudy of rotor cavitie and heat tranfer in a cooling proce in a ga turbine. Journal of Turbomachinery, 6:333{338, 994. [] D. Bohn, G. Lang, H. Schonenborn, and B. Bonho. Determination of thermal tre and train baed on a combined aerodynamic and thermal analyi for a turbine nozzle guide vane. ASME Paper 95-CTP-89, 995. [3] D. Bohn, H. Schonenborn, B. Bonho, and H. Wilhelmi. Prediction of the lm-cooling eectivene in ga turbine blade uing a numerical model for the coupled imulation of uid ow and diabatic wall. ISABE Conference Paper , 995.
17 7 [4] J. Chew, I.J. Taylor, and J.J. Bonell. CFD development for turbine blade heat tranfer. In 3rd International Conference on Reciprocating Engine and Ga Turbine, I. Mech E., London, number C , 994. [5] J. Crank. The Mathematic of Diuion. Clarendon Pre, nd edition, 975. [6] S.K. Godunov and V.S. Ryabenkii. The Theory of Dierence Scheme{An Introduction. North Holland, Amterdam, 964. [7] A. Heelhau and D.T. Vogel. Numerical imulation of turbine blade cooling with repect to blade heat condution and inlet temperature prole. AIAA Paper , 995. [8] A. Heelhau, D.T. Vogel, and H. Krain. Coupling of 3D Navier-Stoke external ow calculation and internal 3D heat conduction calculation for cooled turbine blade. In AGARD CP-57, Heat Tranfer and Cooling in Ga Turbine, 99. [9] J. Moore, J. G. Moore, G. S. Henry, and U. Chaudry. Flow and heat tranfer in turbine tip gap. Journal of Turbomachinery, :30{309, July 989. [0] K.W. Morton and D.F. Mayer. Numerical Solution of Partial Dierential Equation { an Introduction. Cambridge Univerity Pre, Cambridge, 994. [] R.D. Richtmyer and K.W. Morton. Dierence Method for Initial-Value Problem. Wiley-Intercience, nd edition, 967. Reprint edn (994) Krieger Publihing Company, Malabar.
18 8 :5 a) r = 0:5 T 0:5 : : b) r = 75: T 0: 75: Figure : Explicit algorithm with reult every iteration
19 9 :5 a) r = 0:99 0:5 T 0:5 : : b) r = :0 4: T 0: 4: Figure : Explicit algorithm with reult every 5 iteration
20 0 :5 a) r = 5:95 0:5 T 0:5 : : b) r = 6:05 : T : 6: Figure 3: Hybrid algorithm with reult every 5 iteration
21 :5 a) r = : 0:5 T 0:5 : : b) r = : 5: T 5: 5: Figure 4: Implicit algorithm with reult every 5 iteration
Bogoliubov Transformation in Classical Mechanics
Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How
More informationCHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS
CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3
More informationPRESSURE WORK EFFECTS IN UNSTEADY CONVECTIVELY DRIVEN FLOW ALONG A VERTICAL PLATE
Proceeding of 3ICCHMT 3 rd International Conference on Computational Heat and Ma Tranfer May 6 3, 3, Banff, CANADA Paper Number 87 PRESSURE WORK EFFECTS IN UNSTEADY CONVECTIVELY DRIVEN FLOW ALONG A VERTICAL
More information7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281
72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition
More informationSocial Studies 201 Notes for November 14, 2003
1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationControl Systems Analysis and Design by the Root-Locus Method
6 Control Sytem Analyi and Deign by the Root-Locu Method 6 1 INTRODUCTION The baic characteritic of the tranient repone of a cloed-loop ytem i cloely related to the location of the cloed-loop pole. If
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. VIII Decoupling Control - M. Fikar
DECOUPLING CONTROL M. Fikar Department of Proce Control, Faculty of Chemical and Food Technology, Slovak Univerity of Technology in Bratilava, Radlinkého 9, SK-812 37 Bratilava, Slovakia Keyword: Decoupling:
More informationLecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)
Lecture 7: Analytic Function and Integral (See Chapter 4 in Boa) Thi i a good point to take a brief detour and expand on our previou dicuion of complex variable and complex function of complex variable.
More informationin a circular cylindrical cavity K. Kakazu Department of Physics, University of the Ryukyus, Okinawa , Japan Y. S. Kim
Quantization of electromagnetic eld in a circular cylindrical cavity K. Kakazu Department of Phyic, Univerity of the Ryukyu, Okinawa 903-0, Japan Y. S. Kim Department of Phyic, Univerity of Maryland, College
More informationLecture 10 Filtering: Applied Concepts
Lecture Filtering: Applied Concept In the previou two lecture, you have learned about finite-impule-repone (FIR) and infinite-impule-repone (IIR) filter. In thee lecture, we introduced the concept of filtering
More informationCalculation of the temperature of boundary layer beside wall with time-dependent heat transfer coefficient
Ŕ periodica polytechnica Mechanical Engineering 54/1 21 15 2 doi: 1.3311/pp.me.21-1.3 web: http:// www.pp.bme.hu/ me c Periodica Polytechnica 21 RESERCH RTICLE Calculation of the temperature of boundary
More informationPDF hosted at the Radboud Repository of the Radboud University Nijmegen
PDF hoted at the Radboud Repoitory of the Radboud Univerity Nijmegen The following full text i an author' verion which may differ from the publiher' verion. For additional information about thi publication
More informationON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang
Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang
More informationBernoulli s equation may be developed as a special form of the momentum or energy equation.
BERNOULLI S EQUATION Bernoulli equation may be developed a a pecial form of the momentum or energy equation. Here, we will develop it a pecial cae of momentum equation. Conider a teady incompreible flow
More informationOne Class of Splitting Iterative Schemes
One Cla of Splitting Iterative Scheme v Ciegi and V. Pakalnytė Vilniu Gedimina Technical Univerity Saulėtekio al. 11, 2054, Vilniu, Lithuania rc@fm.vtu.lt Abtract. Thi paper deal with the tability analyi
More informationJump condition at the boundary between a porous catalyst and a homogeneous fluid
From the SelectedWork of Francico J. Valde-Parada 2005 Jump condition at the boundary between a porou catalyt and a homogeneou fluid Francico J. Valde-Parada J. Alberto Ochoa-Tapia Available at: http://work.bepre.com/francico_j_valde_parada/12/
More informationSocial Studies 201 Notes for March 18, 2005
1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationSolutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam
BSc - Sample Examination Digital Control Sytem (5-588-) Prof. L. Guzzella Solution Exam Duration: Number of Quetion: Rating: Permitted aid: minute examination time + 5 minute reading time at the beginning
More informationOnline supplementary information
Electronic Supplementary Material (ESI) for Soft Matter. Thi journal i The Royal Society of Chemitry 15 Online upplementary information Governing Equation For the vicou flow, we aume that the liquid thickne
More informationBUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL
J. Autral. Math. Soc. Ser. B 4(999), 332 349 BUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL J. LEE and J.-M. VANDEN-BROECK 2 (Received 22 April 995; revied 23 April 996)
More informationChapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog
Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou
More informationFluid-structure coupling analysis and simulation of viscosity effect. on Coriolis mass flowmeter
APCOM & ISCM 11-14 th December, 2013, Singapore luid-tructure coupling analyi and imulation of vicoity effect on Corioli ma flowmeter *Luo Rongmo, and Wu Jian National Metrology Centre, A*STAR, 1 Science
More informationGiven the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is
EE 4G Note: Chapter 6 Intructor: Cheung More about ZSR and ZIR. Finding unknown initial condition: Given the following circuit with unknown initial capacitor voltage v0: F v0/ / Input xt 0Ω Output yt -
More informationLecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004
18.997 Topic in Combinatorial Optimization April 29th, 2004 Lecture 21 Lecturer: Michel X. Goeman Scribe: Mohammad Mahdian 1 The Lovaz plitting-off lemma Lovaz plitting-off lemma tate the following. Theorem
More informationonline learning Unit Workbook 4 RLC Transients
online learning Pearon BTC Higher National in lectrical and lectronic ngineering (QCF) Unit 5: lectrical & lectronic Principle Unit Workbook 4 in a erie of 4 for thi unit Learning Outcome: RLC Tranient
More informationEXTENDED STABILITY MARGINS ON CONTROLLER DESIGN FOR NONLINEAR INPUT DELAY SYSTEMS. Otto J. Roesch, Hubert Roth, Asif Iqbal
EXTENDED STABILITY MARGINS ON CONTROLLER DESIGN FOR NONLINEAR INPUT DELAY SYSTEMS Otto J. Roech, Hubert Roth, Aif Iqbal Intitute of Automatic Control Engineering Univerity Siegen, Germany {otto.roech,
More informationMarch 18, 2014 Academic Year 2013/14
POLITONG - SHANGHAI BASIC AUTOMATIC CONTROL Exam grade March 8, 4 Academic Year 3/4 NAME (Pinyin/Italian)... STUDENT ID Ue only thee page (including the back) for anwer. Do not ue additional heet. Ue of
More informationConvex Hulls of Curves Sam Burton
Convex Hull of Curve Sam Burton 1 Introduction Thi paper will primarily be concerned with determining the face of convex hull of curve of the form C = {(t, t a, t b ) t [ 1, 1]}, a < b N in R 3. We hall
More informationinto a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get
Lecture 25 Introduction to Some Matlab c2d Code in Relation to Sampled Sytem here are many way to convert a continuou time function, { h( t) ; t [0, )} into a dicrete time function { h ( k) ; k {0,,, }}
More informationX R. U x U x B. U x U x B X R. U x U x B. U x U x B. one solution. solution. two solutions no solution. one. R two solutions no solution.
CLASSIFICATION OF CODIMENSION-ONE RIEMANN SOLUTIONS STEPHEN SCHECTER, BRADLEY J. PLOHR, AND DAN MARCHESIN Abtract. We invetigate olution of Riemann problem for ytem of two conervation law in one patial
More informationDigital Control System
Digital Control Sytem - A D D A Micro ADC DAC Proceor Correction Element Proce Clock Meaurement A: Analog D: Digital Continuou Controller and Digital Control Rt - c Plant yt Continuou Controller Digital
More informationGain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays
Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,
More informationPreemptive scheduling on a small number of hierarchical machines
Available online at www.ciencedirect.com Information and Computation 06 (008) 60 619 www.elevier.com/locate/ic Preemptive cheduling on a mall number of hierarchical machine György Dóa a, Leah Eptein b,
More informationSingular perturbation theory
Singular perturbation theory Marc R. Rouel June 21, 2004 1 Introduction When we apply the teady-tate approximation (SSA) in chemical kinetic, we typically argue that ome of the intermediate are highly
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction
More informationNew bounds for Morse clusters
New bound for More cluter Tamá Vinkó Advanced Concept Team, European Space Agency, ESTEC Keplerlaan 1, 2201 AZ Noordwijk, The Netherland Tama.Vinko@ea.int and Arnold Neumaier Fakultät für Mathematik, Univerität
More informationPhysics 741 Graduate Quantum Mechanics 1 Solutions to Final Exam, Fall 2014
Phyic 7 Graduate Quantum Mechanic Solution to inal Eam all 0 Each quetion i worth 5 point with point for each part marked eparately Some poibly ueful formula appear at the end of the tet In four dimenion
More informationCONGESTION control is a key functionality in modern
IEEE TRANSACTIONS ON INFORMATION TEORY, VOL. X, NO. X, XXXXXXX 2008 On the Connection-Level Stability of Congetion-Controlled Communication Network Xiaojun Lin, Member, IEEE, Ne B. Shroff, Fellow, IEEE,
More informationIEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation
IEOR 316: Fall 213, Profeor Whitt Topic for Dicuion: Tueday, November 19 Alternating Renewal Procee and The Renewal Equation 1 Alternating Renewal Procee An alternating renewal proce alternate between
More informationChapter 4. The Laplace Transform Method
Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination
More informationDYNAMIC REDESIGN OF A FLOW CONTROL SERVO-VALVE USING A PRESSURE CONTROL PILOT
Proceeding of IMECE ASME International Mechanical Engineering Congre & Exhibition November -6,, New York, New York, USA IMECE/DSC-B- DYNAMIC REDESIGN OF A FLOW CONTROL SERVO-VALVE USING A PRESSURE CONTROL
More informationS_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS
S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS by Michelle Gretzinger, Daniel Zyngier and Thoma Marlin INTRODUCTION One of the challenge to the engineer learning proce control i relating theoretical
More informationarxiv: v2 [nucl-th] 3 May 2018
DAMTP-207-44 An Alpha Particle Model for Carbon-2 J. I. Rawlinon arxiv:72.05658v2 [nucl-th] 3 May 208 Department of Applied Mathematic and Theoretical Phyic, Univerity of Cambridge, Wilberforce Road, Cambridge
More informationPractice Problems - Week #7 Laplace - Step Functions, DE Solutions Solutions
For Quetion -6, rewrite the piecewie function uing tep function, ketch their graph, and find F () = Lf(t). 0 0 < t < 2. f(t) = (t 2 4) 2 < t In tep-function form, f(t) = u 2 (t 2 4) The graph i the olid
More informationSMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD
SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD S.P. Teeuwen, I. Erlich U. Bachmann Univerity of Duiburg, Germany Department of Electrical Power Sytem
More information[Saxena, 2(9): September, 2013] ISSN: Impact Factor: INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY
[Saena, (9): September, 0] ISSN: 77-9655 Impact Factor:.85 IJESRT INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY Contant Stre Accelerated Life Teting Uing Rayleigh Geometric Proce
More informationCHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL
98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i
More informationAMS 212B Perturbation Methods Lecture 20 Part 1 Copyright by Hongyun Wang, UCSC. is the kinematic viscosity and ˆp = p ρ 0
Lecture Part 1 Copyright by Hongyun Wang, UCSC Prandtl boundary layer Navier-Stoke equation: Conervation of ma: ρ t + ( ρ u) = Balance of momentum: u ρ t + u = p+ µδ u + ( λ + µ ) u where µ i the firt
More informationStudy of a Freely Falling Ellipse with a Variety of Aspect Ratios and Initial Angles
Study of a Freely Falling Ellipe with a Variety of Apect Ratio and Initial Angle Dedy Zulhidayat Noor*, Ming-Jyh Chern*, Tzyy-Leng Horng** *Department of Mechanical Engineering, National Taiwan Univerity
More informationOn the Curv ature of Space ²
Gener al Relativity and Gravitation, Vol. 31, No. 1, 1999 On the Curv ature of Space ² By A. Friedman in Peterburg * With one gure. Received on 9. June 19 1. 1. In their well-known work on general comologic
More informationFactor Analysis with Poisson Output
Factor Analyi with Poion Output Gopal Santhanam Byron Yu Krihna V. Shenoy, Department of Electrical Engineering, Neurocience Program Stanford Univerity Stanford, CA 94305, USA {gopal,byronyu,henoy}@tanford.edu
More informationGreen-Kubo formulas with symmetrized correlation functions for quantum systems in steady states: the shear viscosity of a fluid in a steady shear flow
Green-Kubo formula with ymmetrized correlation function for quantum ytem in teady tate: the hear vicoity of a fluid in a teady hear flow Hirohi Matuoa Department of Phyic, Illinoi State Univerity, Normal,
More informationFinite Element Analysis of a Fiber Bragg Grating Accelerometer for Performance Optimization
Finite Element Analyi of a Fiber Bragg Grating Accelerometer for Performance Optimization N. Baumallick*, P. Biwa, K. Dagupta and S. Bandyopadhyay Fiber Optic Laboratory, Central Gla and Ceramic Reearch
More informationA Single Particle Thermal Model for Lithium Ion Batteries
A Single Particle Thermal Model for Lithium Ion Batterie R. Painter* 1, B. Berryhill 1, L. Sharpe 2 and S. Keith Hargrove 2 1 Civil Engineering, Tenneee State Univerity, Nahville, TN, USA 2 Mechanical
More informationFaculty of Environmental Sciences, Institute of Waste Management and Contaminated Site Treatment. The Simulation Software.
Faculty of Environmental Science, Intitute of Wate Management and Contaminated Site Treatment The Simulation Software PCSiWaPro Overview 1. Modelling in the unaturated oil zone 2. The oftware PCSiWaPro
More informationCake ltration analysis the eect of the relationship between the pore liquid pressure and the cake compressive stress
Chemical Engineering Science 56 (21) 5361 5369 www.elevier.com/locate/ce Cake ltration analyi the eect of the relationhip between the pore liquid preure and the cake compreive tre C. Tien, S. K. Teoh,
More informationA BATCH-ARRIVAL QUEUE WITH MULTIPLE SERVERS AND FUZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH
Mathematical and Computational Application Vol. 11 No. pp. 181-191 006. Aociation for Scientific Reearch A BATCH-ARRIVA QEE WITH MTIPE SERVERS AND FZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH Jau-Chuan
More informationSTOCHASTIC DIFFERENTIAL GAMES:THE LINEAR QUADRATIC ZERO SUM CASE
Sankhyā : he Indian Journal of Statitic 1995, Volume 57, Serie A, Pt. 1, pp.161 165 SOCHASIC DIFFERENIAL GAMES:HE LINEAR QUADRAIC ZERO SUM CASE By R. ARDANUY Univeridad de Salamanca SUMMARY. hi paper conider
More informationSimple Food Chain in a Chemostat with Distinct Removal Rates
Journal of Mathematical Analyi and Application 242, 7592 Ž 2000. doi:0.006jmaa.999.6655, available online at http:www.idealibrary.com on Simple Food Chain in a Chemotat with Ditinct Removal Rate Bingtuan
More informationRecent progress in fire-structure analysis
EJSE Special Iue: Selected Key Note paper from MDCMS 1 1t International Conference on Modern Deign, Contruction and Maintenance of Structure - Hanoi, Vietnam, December 2007 Recent progre in fire-tructure
More informationTHE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER
Proceeding of IMAC XXXI Conference & Expoition on Structural Dynamic February -4 Garden Grove CA USA THE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER Yung-Sheng Hu Neil S Ferguon
More informationModeling of Transport and Reaction in a Catalytic Bed Using a Catalyst Particle Model.
Excerpt from the Proceeding of the COMSOL Conference 2010 Boton Modeling of Tranport and Reaction in a Catalytic Bed Uing a Catalyt Particle Model. F. Allain *,1, A.G. Dixon 1 1 Worceter Polytechnic Intitute
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Laplace Tranform Paul Dawkin Table of Content Preface... Laplace Tranform... Introduction... The Definition... 5 Laplace Tranform... 9 Invere Laplace Tranform... Step Function...4
More informationResearch Article A New Kind of Weak Solution of Non-Newtonian Fluid Equation
Hindawi Function Space Volume 2017, Article ID 7916730, 8 page http://doi.org/10.1155/2017/7916730 Reearch Article A New Kind of Weak Solution of Non-Newtonian Fluid Equation Huahui Zhan 1 and Bifen Xu
More informationBeta Burr XII OR Five Parameter Beta Lomax Distribution: Remarks and Characterizations
Marquette Univerity e-publication@marquette Mathematic, Statitic and Computer Science Faculty Reearch and Publication Mathematic, Statitic and Computer Science, Department of 6-1-2014 Beta Burr XII OR
More informationA Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems
A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement
More informationOnline Appendix for Managerial Attention and Worker Performance by Marina Halac and Andrea Prat
Online Appendix for Managerial Attention and Worker Performance by Marina Halac and Andrea Prat Thi Online Appendix contain the proof of our reult for the undicounted limit dicued in Section 2 of the paper,
More information84 ZHANG Jing-Shang Vol. 39 of which would emit 5 He rather than 3 He. 5 He i untable and eparated into n + pontaneouly, which can alo be treated a if
Commun. Theor. Phy. (Beijing, China) 39 (003) pp. 83{88 c International Academic Publiher Vol. 39, No. 1, January 15, 003 Theoretical Analyi of Neutron Double-Dierential Cro Section of n+ 11 B at 14. MeV
More informationSERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lectures 41-48)
Chapter 5 SERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lecture 41-48) 5.1 Introduction Power ytem hould enure good quality of electric power upply, which mean voltage and current waveform hould
More informationMATEMATIK Datum: Tid: eftermiddag. A.Heintz Telefonvakt: Anders Martinsson Tel.:
MATEMATIK Datum: 20-08-25 Tid: eftermiddag GU, Chalmer Hjälpmedel: inga A.Heintz Telefonvakt: Ander Martinon Tel.: 073-07926. Löningar till tenta i ODE och matematik modellering, MMG5, MVE6. Define what
More informationOn the Normal-Mode Frequency Spectrum of Kinetic Magnetohydrodynamics. J.J. Ramos. December, 2014
PSFC/JA-14-33 On the Normal-Mode Frequency Spectrum of Kinetic Magnetohydrodynamic J.J. Ramo December, 14 Plama Science and Fuion Center Maachuett Intitute of Technology Cambridge MA 139, U.S.A. Thi work
More informationIII.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SUBSTANCES
III.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SBSTANCES. Work purpoe The analyi of the behaviour of a ferroelectric ubtance placed in an eternal electric field; the dependence of the electrical polariation
More informationEP225 Note No. 5 Mechanical Waves
EP5 Note No. 5 Mechanical Wave 5. Introduction Cacade connection of many ma-pring unit conitute a medium for mechanical wave which require that medium tore both kinetic energy aociated with inertia (ma)
More informationCHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS
Copyright 22 IFAC 5th Triennial World Congre, Barcelona, Spain CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS Tritan Pérez Graham C. Goodwin Maria M. Serón Department of Electrical
More informationAn Inequality for Nonnegative Matrices and the Inverse Eigenvalue Problem
An Inequality for Nonnegative Matrice and the Invere Eigenvalue Problem Robert Ream Program in Mathematical Science The Univerity of Texa at Dalla Box 83688, Richardon, Texa 7583-688 Abtract We preent
More informationThermal Stress in a Half-Space with Mixed Boundary Conditions due to Time Dependent Heat Source
IOSR Journal of Mathematic (IOSR-JM) e-issn: 2278-5728, p-issn: 239-765X Volume, Iue 6 Ver V (Nov - Dec 205), PP 9-25 wwwiorjournalorg Thermal Stre in a Half-Space with Mixed Boundary Condition due to
More informationConvergence criteria and optimization techniques for beam moments
Pure Appl. Opt. 7 (1998) 1221 1230. Printed in the UK PII: S0963-9659(98)90684-5 Convergence criteria and optimization technique for beam moment G Gbur and P S Carney Department of Phyic and Atronomy and
More informationTHE PARAMETERIZATION OF ALL TWO-DEGREES-OF-FREEDOM SEMISTRONGLY STABILIZING CONTROLLERS. Tatsuya Hoshikawa, Kou Yamada and Yuko Tatsumi
International Journal of Innovative Computing, Information Control ICIC International c 206 ISSN 349-498 Volume 2, Number 2, April 206 pp. 357 370 THE PARAMETERIZATION OF ALL TWO-DEGREES-OF-FREEDOM SEMISTRONGLY
More informationStandard Guide for Conducting Ruggedness Tests 1
Deignation: E 69 89 (Reapproved 996) Standard Guide for Conducting Ruggedne Tet AMERICA SOCIETY FOR TESTIG AD MATERIALS 00 Barr Harbor Dr., Wet Conhohocken, PA 948 Reprinted from the Annual Book of ASTM
More informationCharacterization of the heat transfer in open-cell metal foam
Characterization of the heat tranfer in open-cell metal foam C. Briano-Calcagno, J. Fontánez-Delgado & N. Dukhan Department of Mechanical Engineering, Univerity of Puerto Rico Mayagüez, Mayagüez, P.R.,
More informationConvective Heat Transfer
Convective Heat Tranfer Example 1. Melt Spinning of Polymer fiber 2. Heat tranfer in a Condener 3. Temperature control of a Re-entry vehicle Fiber pinning The fiber pinning proce preent a unique engineering
More informationUnified Correlation between SPT-N and Shear Wave Velocity for all Soil Types
6 th International Conference on Earthquake Geotechnical Engineering 1-4 ovember 15 Chritchurch, ew Zealand Unified Correlation between SPT- and Shear Wave Velocity for all Soil Type C.-C. Tai 1 and T.
More informationMulticolor Sunflowers
Multicolor Sunflower Dhruv Mubayi Lujia Wang October 19, 2017 Abtract A unflower i a collection of ditinct et uch that the interection of any two of them i the ame a the common interection C of all of
More informationCONTROL SYSTEMS. Chapter 5 : Root Locus Diagram. GATE Objective & Numerical Type Solutions. The transfer function of a closed loop system is
CONTROL SYSTEMS Chapter 5 : Root Locu Diagram GATE Objective & Numerical Type Solution Quetion 1 [Work Book] [GATE EC 199 IISc-Bangalore : Mark] The tranfer function of a cloed loop ytem i T () where i
More informationME 375 EXAM #1 Tuesday February 21, 2006
ME 375 EXAM #1 Tueday February 1, 006 Diviion Adam 11:30 / Savran :30 (circle one) Name Intruction (1) Thi i a cloed book examination, but you are allowed one 8.5x11 crib heet. () You have one hour to
More informationThen C pid (s) S h -stabilizes G(s) if and only if Ĉpid(ŝ) S 0 - stabilizes Ĝ(ŝ). For any ρ R +, an RCF of Ĉ pid (ŝ) is given by
9 American Control Conference Hyatt Regency Riverfront, St. Loui, MO, USA June -, 9 WeC5.5 PID Controller Synthei with Shifted Axi Pole Aignment for a Cla of MIMO Sytem A. N. Gündeş and T. S. Chang Abtract
More informationOn the Stability Region of Congestion Control
On the Stability Region of Congetion Control Xiaojun Lin and Ne B. Shroff School of Electrical and Computer Engineering Purdue Univerity, Wet Lafayette, IN 47906 {linx,hroff}@ecn.purdue.edu Abtract It
More informationNONLINEAR CONTROLLER DESIGN FOR A SHELL AND TUBE HEAT EXCHANGER AN EXPERIMENTATION APPROACH
International Journal of Electrical, Electronic and Data Communication, ISSN: 232-284 Volume-3, Iue-8, Aug.-25 NONLINEAR CONTROLLER DESIGN FOR A SHELL AND TUBE HEAT EXCHANGER AN EXPERIMENTATION APPROACH
More informationHyperbolic Partial Differential Equations
Hyperbolic Partial Differential Equation Evolution equation aociated with irreverible phyical procee like diffuion heat conduction lead to parabolic partial differential equation. When the equation i a
More informationEvolutionary Algorithms Based Fixed Order Robust Controller Design and Robustness Performance Analysis
Proceeding of 01 4th International Conference on Machine Learning and Computing IPCSIT vol. 5 (01) (01) IACSIT Pre, Singapore Evolutionary Algorithm Baed Fixed Order Robut Controller Deign and Robutne
More informationIntroduction to Laplace Transform Techniques in Circuit Analysis
Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found
More informationModeling the scalar wave equation with Nyström methods
GEOPHYSICS, VOL. 71, NO. 5 SEPTEMBER-OCTOBER 200 ; P. T151 T158, 7 FIGS. 10.1190/1.2335505 Modeling the calar wave equation with Nytröm method Jing-Bo Chen 1 ABSTRACT High-accuracy numerical cheme for
More informationList coloring hypergraphs
Lit coloring hypergraph Penny Haxell Jacque Vertraete Department of Combinatoric and Optimization Univerity of Waterloo Waterloo, Ontario, Canada pehaxell@uwaterloo.ca Department of Mathematic Univerity
More informationDEVELOPMENT OF A STRUCTURED THERMOCLINE THERMAL ENERGY STORAGE SYSTEM
DEVELOPMENT OF A STRUCTURED THERMOCLINE THERMAL ENERGY STORAGE SYSTEM Brad M. Brown Matt N. Straer R. Paneer Selvam Univerity of Arkana Department of Civil Engineering 4190 Bell Engineering Center Fayetteville,
More informationORE Open Research Exeter
ORE Open Reearch Exeter TITLE The onet of convection in rotating circular cylinder with experimental boundary condition AUTHORS Zhang, Keke; Liao, X. JOURNAL Journal of Fluid Mechanic DEPOSITED IN ORE
More informationThe Impact of Imperfect Scheduling on Cross-Layer Rate. Control in Multihop Wireless Networks
The mpact of mperfect Scheduling on Cro-Layer Rate Control in Multihop Wirele Network Xiaojun Lin and Ne B. Shroff Center for Wirele Sytem and Application (CWSA) School of Electrical and Computer Engineering,
More informationFermi Distribution Function. n(e) T = 0 T > 0 E F
LECTURE 3 Maxwell{Boltzmann, Fermi, and Boe Statitic Suppoe we have a ga of N identical point particle in a box ofvolume V. When we ay \ga", we mean that the particle are not interacting with one another.
More informationHELICAL TUBES TOUCHING ONE ANOTHER OR THEMSELVES
15 TH INTERNATIONAL CONFERENCE ON GEOMETRY AND GRAPHICS 0 ISGG 1-5 AUGUST, 0, MONTREAL, CANADA HELICAL TUBES TOUCHING ONE ANOTHER OR THEMSELVES Peter MAYRHOFER and Dominic WALTER The Univerity of Innbruck,
More informationHybrid Projective Dislocated Synchronization of Liu Chaotic System Based on Parameters Identification
www.ccenet.org/ma Modern Applied Science Vol. 6, No. ; February Hybrid Projective Dilocated Synchronization of Liu Chaotic Sytem Baed on Parameter Identification Yanfei Chen College of Science, Guilin
More informationRELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS
www.arpapre.com/volume/vol29iue1/ijrras_29_1_01.pdf RELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS Sevcan Demir Atalay 1,* & Özge Elmataş Gültekin
More information