Final Report on COOMET Key Comparison of National Pressure Standards in the Range 100 Pa to 5 kpa of Gauge Pressure (COOMET.M.

Size: px
Start display at page:

Download "Final Report on COOMET Key Comparison of National Pressure Standards in the Range 100 Pa to 5 kpa of Gauge Pressure (COOMET.M."

Transcription

1 Verson 1 of Fnal Report on COOMET Key Comparson of Natonal Pressure Standards n the Range 100 Pa to 5 kpa of Gauge Pressure (COOMET.M.P-K14) H. Ahrendt 1, J. Könemann 1, W. Sabuga 1, Y. Kselev 2, O. Vtkovsky 2, D. Pražák 3, Z. Krajcek 3, K. Dapkevčenė 4, ABSTRACT Ths report descrbes a COOMET key comparson of pneumatc gauge pressure standards of four Natonal Metrology Insttutes, lsted n the chronologcal order of ther measurements, that was carred out n the perod from November 2008 to March 2010 n order to determne ther degree of equvalence n the range of 100 Pa to 5 kpa of gauge pressure. The plot laboratory was PTB. The reference pressure standards of the partcpants were of dfferent desgn. The transfer standard was a pston gauge model V1600 of the company Pressurements. The quantty under comparson was the effectve area of the transfer standard at dfferent pressure values reported together wth uncertanty contrbutons and the conclusve combned uncertanty of measurement. All partcpants results agree wth the key comparson reference values wthn the expanded uncertantes calculated wth a coverage factor 2, all but one results even wthn the standard uncertantes. For the partcpants results compared n pars, all of totally 48 pars show agreement wthn the expanded and 46 pars wthn the standard uncertantes. The results of the comparson demonstrate equvalence of the laboratory standards and support ther measurement capabltes stated n the KCDB of BIPM. 1 Physkalsch-Technsche Bundesanstalt (PTB), Braunschweg, Germany 2 D.I. Mendeleyev Insttute for Metrology (VNIIM), St. Petersburg, Russa 3 Czech Metrology Insttute (CMI), Brno, Czech Republc 4 Vlnus Metrology Centre (VMC), Vlnus, Lthuana

2 CONTENTS 1. INTRODUCTION PARTICIPANTS PTB pressure standards and measurement method VNIIM pressure standard and measurement method CMI pressure standard and measurement method VMC pressure standard and measurement method TRANSFER STANDARD ORGANIZATION AND CHRONOLOGY OF THE COMPARISON MEASUREMENT PROCEDURES RESULTS Stablty of the transfer standard Results of the Partcpants Reference value calculaton Degree of equvalence DISCUSSION CONCLUSIONS...20 REFERENCES...20 APPENDIX A. Relatve devatons of the partcpants results from the reference value ((A p, - A p,ref )/A p,ref ) wth the expanded uncertantes (k = 2) of these devatons at pressures 100 Pa to 5 kpa...21 APPENDIX B. Relatve dfferences between the partcpants results (d j ) and ther expanded uncertantes (U(d j )) (k = 2) at pressures 100 Pa to 5 kpa...24 APPENDIX C. Stablty of TS

3 1. INTRODUCTION At the COOMET TCM meetng held at KazInMetr, Astana, on 8 October 2008, t was decded to carry out a key comparson (KC) n the range of 5 kpa of pneumatc gauge pressure. Pont of am of ths comparson was to state the equvalence of the natonal pressure standards of the partcpatng natonal metrologcal laboratores n the 5 kpa range. Two partcpants n ths comparson, the Physkalsch-Technsche Bundesanstalt (PTB), Germany, and the Czech Metrology Insttute (CMI), Czech Republc, are also partcpants n the respectve KCs of EURAMET, EURAMET.M.P-K4.2010, beng n progress, and a CCM KC beng stll n preparaton. Ths wll allow a lnk of the current COOMET KC to the correspondng EURAMET and CCM KCs. PTB was agreed to be a plot laboratory n ths KC and provded a transfer standard (TS) for ths KC. The comparson was carred out n accordance wth ts Techncal Protocol whch specfed the procedures to be followed n the comparson and had been prepared n accordance wth the Gudelnes for CIPM Key Comparsons, 1 March PARTICIPANTS The laboratory standards (LS) used n ths comparson were of dfferent types, namely a dvng bell manometer, deadweght and force-controlled pston gauges, and a water column mcromanometer. Dfferent methods were appled by the partcpants to compare ther standards wth the TS. All uncertantes relatng to the partcpants standards and ther results n ths report are standard ones f not explctly ndcated as expanded. 2.1 PTB pressure standards and measurement method Two dfferent laboratory standards were used by PTB n ths KC. For the pressure values up to nomnally 4 kpa, the dvng bell known as TGM 5 wth the mass set dentfed by g01ma2 was compared to TS. Ths bell manometer was developed at the ASMW [1] and was transferred to PTB n 1990 actng there as a reference standard for small gauge and dfferental pressures up to 4 kpa. In order to cover ths pressure range t was equpped wth Fluornert nstead of nonane as a workng flud, and, moreover, the load mechansm was drectly attached to the bell. It was used n a comparson of European dfferental pressure standards from 3 Pa to 1 kpa [2]. The effectve area of the dvng bell s traceable to a prmary mercury manometer of PTB. Table 1 summarzes propertes of the dvng bell manometers and the related measurement condtons. 3

4 Table 1. PTB dvng bell manometer for pressures up to 4 kpa and measurement condtons Type of devce Dvng bell manometer Manufacturer ASMW, Germany Measurement range n Pa, mode (1 4000), gauge Materal of dvng bell or pston brass Workng lqud (for dvng bell manometer) Fluornert Pressure-transmttng gas Ar or ntrogen Zero-pressure effectve area (A 0 ) at reference temperature n m Relatve uncertanty of A 0 n Uncertanty of mass (m) peces or of mass measurement 125 mg m Lnear thermal expanson coeffcent of dvng bell manometer (α p ) n C Reference temperature (t 0 ) n C 20 Local gravty (g) n m/s Relatve uncertanty of g n Heght dfference between laboratory standard (LS) and TS (h, postve f LS s hgher than 0 TS) n mm Uncertanty of h n mm 5 For pressure values 4800 and 5000 Pa a gas-operated pressure balance TL1568 from Ruska company was used together wth the mass set p04ma. The propertes of ths pressure balance are presented n Table 2. 4

5 Table 2. PTB pressure balance for pressures of 4.8 and 5 kpa Type of devce Manufacturer Measurement range n kpa, mode Materal of dvng bell or pston Materal of cylnder (for pressure balance) Pressure-transmttng gas Zero-pressure effectve area (A 0 ) at reference temperature n m Relatve uncertanty of A 0 n Uncertanty of mass peces at p max, mg 3 Lnear thermal expanson coeffcent of dvng bell or pston (α p ) n C Lnear thermal expanson coeffcent of cylnder (for pressure balance, α c ) n C Reference temperature (t 0 ) n C 20 Local gravty (g) n m/s Relatve uncertanty of g n Heght dfference between laboratory standard (LS) and TS (h, postve f LS s hgher than TS) n mm TL1568 Ruska , gauge and absolute steel Tungsten carbde Dry gas/ntrogen Uncertanty of h n mm VNIIM pressure standard and measurement method The VNIIM pressure standard s a water column mcromanometer whch s a part of the natonal specal standard of unt of pressure GET (GOST «State system for ensurng the unformty of measurements. The State specal standard and All-Unon verfcaton schedule for means of measurements of the dfference of pressures up to Pa») havng desgnaton "MKШ". Physcally, t s a water-based compensated mcromanometer workng n a pressure range up to 5 kpa wth a traceablty of ts propertes to VNIIM standards. Wth the TS used n the KC and under the condtons of the experment, the maxmum pressure of the laboratory standard was lmted by 4.8 kpa. Hence, devatng from the Techncal Protocol, the hghest measured pressure was 4.8 kpa. TS was transported to VNIIM by the plot laboratory on 21 November 2008 and measured there between 24 and 26 November It was returned to PTB on 28 November The propertes of the VNIIM pressure standard and the measurement condtons are gven n Table 3. 5

6 Table 3. VNIIM water column mcromanometer and measurement condtons Type of devce Water column mcromanometer Manufacturer VNIIM, Russa Pressure (p) measurement range n Pa Relatve uncertanty of pressure for 100 Pa: for 5000 Pa: Workng lqud Dstlled water Densty of workng lqud (ρ l ) kg/m Relatve uncertanty of ρ l n Reference temperature of workng lqud (t 0 ) n C 20 Uncertanty of t 0 n C 0.1 Pressure-transmttng gas Ntrogen Local gravty (g) n m/s Relatve uncertanty of g n Heght dfference between laboratory standard (LS) and TS (h, postve f LS s hgher than 20 TS) n mm Uncertanty of h n mm CMI pressure standard and measurement method The CMI pressure standard used n ths comparson s a force-balanced pston gauge (FPG) manufactured by DH Instruments, USA, model FPG 8601, whose descrpton of the physcal prncple s descrbed n [3] and the commercal realsaton n [4]. The setup of the CMI FPG was modfed wth respect to the commercal one n such a way that a turbo-molecular pump was added to the pumpng system of the reference sde n order to reduce the resdual pressure on the reference sde and, therefore, ts uncertanty contrbuton. A thorough evaluaton of the CMI system s presented n [5]. The effectve area was evaluated both by a dmensonal measurement of the pstoncylnder geometry and by cross-floatng technques comparng t aganst the CMI PG 7601 pressure balance. These comparsons were performed wth the CMI standards durng the year An ntercomparson wth the Slovak natonal metrology nsttute (SMU) was performed n December 2002 from 2 to 15 kpa n both gauge and absolute pressure mode wth a Bell and Howell pressure balance as a transfer standard. Another ntercomparson wth the Fnnsh natonal metrology nsttute, (MIKES), was realzed n July 2003 from 1 Pa to 15 kpa n gauge pressure mode and from 6 Pa to 15 kpa n absolute pressure mode, the MIKES standard beng another FPG A thrd blateral comparson, EURAMET.M.P-S2, was carred out wth PTB that used an FRS5 force-compensated pressure balance manufactured by Furness Controls, UK. For the actual measurements an MKS Baratron dfferental pressure cell (DPC) was used as a zero pressure ndcator nstalled between the CMI FPG and the TS. Wth ths DPC CMI could automate the measurements usng software Compass of the FPG. The techncal detals of the CMI standard and the measurement condtons are gven n Table 4. 6

7 The transfer standard was transported to CMI by the plot laboratory on 11 May The measurements were performed n the perod from 12 to 13 May 2009, accordng to the Techncal Protocol. Afterwards, TS was brought back to PTB on 15 May Table 4. CMI force-balanced pston gauge and measurement condtons Type of devce FPG 8601 Manufacturer DH Instruments Measurement range n Pa , gauge and absolute Materal of dvng bell or pston tungsten-carbde Materal of cylnder (for pressure balance) tungsten-carbde Pressure-transmttng gas ntrogen Zero-pressure effectve area (A 0 ) at reference temperature n m Relatve uncertanty of A 0 n Uncertanty of mass (m) peces or of mass measurement 0.5 mg m Lnear thermal expanson coeffcent of dvng bell or pston (α p ) n C Lnear thermal expanson coeffcent of cylnder (for pressure balance, α c ) n C Reference temperature (t 0 ) n C 20 Local gravty (g) n m/s Relatve uncertanty of g n Heght dfference between laboratory standard (LS) and TS (h, postve f LS s hgher than 124 TS) n mm Uncertanty of h n mm 2/ VMC pressure standard and measurement method The VMC laboratory standard was a pston gauge Metran-505, Vozduch I, manufactured by Metran, Russa, seral number 094, made n The devce was calbrated at PTB (Germany) n August Thus, ts pressure value s traceable to the German natonal pressure standards. The pressure values realsed wth pstons M and B n combnaton wth weghts and the uncertantes of that pressure values were determned at PTB and reported n Calbraton certfcate Nr PTB Fnally, the effectve area A(p) of pston M assembly and pston B assembly were calculated at VMC. The techncal detals of the VMC standard and the measurement condtons are gven n Table 5. 7

8 Table 5. VMC pston gauge and measurement condtons Type of devce Metran-505, Vozduch I, seral Nr. 094 Manufacturer Metran Measurement range, mode 20 Pa to 25 kpa: gauge pressure 5 Pa to 25 kpa: dfferental pressure Materal of pstons up to 125 Pa - alumnum 160 Pa to 25 kpa - stanless steel Materal of cylnder (for pressure balance) stanless steel Pressure-transmttng gas Ar/ntrogen Zero-pressure effectve area (A 0 ) at reference temperature n m 2 3 pstons Relatve uncertanty of A 0 n Uncertanty of mass (m) peces or of mass measurement, mg 7.1 Lnear thermal expanson coeffcent of pston and cylnder (α p ) n C -1 (23±2) 10-6 Reference temperature (t 0 ) n C 20 Local gravty (g) n m/s Relatve uncertanty of g n Heght dfference between LS and TS (h, postve f LS s hgher than TS) n mm -6 Uncertanty of h n mm 5 TS and the laboratory standard were nterconnected va a DPC usng ntrogen as a pressure-transmttng medum. TS and the laboratory standard were loaded wth masses correspondng to the nomnal pressures as specfed n the Techncal Protocol and the resdual pressure dfferences between them were measured wth the DPC. Hence, the uncertanty of the DPC, whch accordng to the plot laboratory s Uncertanty of zero ndcator 0.02 Pa was added to the uncertanty of the VMC standard. TS was brought to VMC on 23 November 2009 and measured between 24 and 27 November, before t was returned to PTB on 27 November TRANSFER STANDARD The transfer standard was descrbed n detal n the techncal protocol. It was a pston gauge V1600/1D manufactured by Pressurements Ltd., Bedfordshre, England, n 1998 and dentfed by seral number The assembly conssted of a cylnder dentfed by seral number B 192, whch s engraved on the cylnder face, and three dfferent pstons all carryng number N 192 on the lower pston face and dentfed by markngs "1 mbar", "2 mbar" and "L". The nomnal effectve area of the assembly was A 0,nom = 4.91 cm 2. The cylnder s materal was stanless steel, whereas pston "1 mbar" was made of alumnum and pstons "2 mbar" and "L" were fabrcated agan from stanless steel. 8

9 Followng the manufacturer nformaton, the thermal expanson coeffcent of the assembly of cylnder and any of the pstons (ether of steel or of alumnum) was taken as (22 ± 2) 10-6 C -1 wth ts standard uncertanty based both on the manufacturer nformaton and the experence of the plot laboratory. Table 6. Masses and denstes of pstons Pston True mass n g Densty n kg/m 3 "1 mbar" ± ± 25 "2 mbar" ± ± 25 "L" ± ± 25 The weght carrer was marked wth "L" and made of alumnum. Table 7. Masses and denstes of the weght carrer Weght carrer True mass n g Densty n kg/m ± ± 25 Seven weghts were dentfed by markngs "1 mbar", "5 mbar", "10 mbar", "20_1 mbar", "20_2 mbar", and addtonally by seral number N192 on each pece. Table 8. Masses and denstes of the weghts Weght True mass n g Densty n kg/m 3 "1 mbar" ± ± 25 "2_1 mbar" ± ± 25 "2_2 mbar" ± ± 25 "5 mbar" ± ± 25 "10 mbar" ± ± 25 "20_1 mbar" ± ± 25 "20_2 mbar" ± ± 25 In combnaton wth the pston gauge, a zero pressure ndcator FC 014 was delvered to the partcpants, whch had been manufactured by Furness Controls and had seral number Ths ndcator had a scale nterval of 0.02 Pa and a maxmum dfferental pressure range of 100 Pa at a maxmum acceptable lne pressure of 100 kpa and a maxmum overload dfferental pressure of 100 kpa. 4. ORGANIZATION AND CHRONOLOGY OF THE COMPARISON The measurements were performed n the order and tmes gven n Table 9. 9

10 Table 9. Chronology of measurements Insttute Measurement start date Measurement end date PTB, ntal nvestgaton of TS 4 Nov Nov 2008 VNIIM 24 Nov Nov 2008 PTB, ntermedate check of TS 30 Mar Apr 2009 CMI 12 May May 2009 PTB, ntermedate check of TS 12 Nov Nov 2009 VMC 24 Nov Nov 2009 PTB, fnal check of TS 4 Mar Mar 2010 The plot laboratory performed three measurement cycles and one half-cycle measurement (November 2009) of TS. The measurement n March 2010 was taken as a PTB contrbuton to ths KC. 5. MEASUREMENT PROCEDURES The transfer standard had to be handled and the pston-cylnder assembly mounted n accordance wth the nstructons gven n the User s Manual Reference of the V1600 provded to the partcpants. TS was operated throughout all measurements by one member of the plot laboratory (Mrs. Ahrendt). The TS was recommended to be located close to the laboratory s reference standard to keep the pressure lne between the two nstruments as short as possble. The pston gauge V1600/ID of TS was equpped wth two cylnders. The cylnder dentfed by "B" was used n the comparson measurements. The reference level of TS was the upper face of cylnder B. The horzontalty of TS expected to be better than 0.1 mm/m was checked wth a sprt level placed on the upper face of cylnder B. The temperature of TS was measured wth a thermometer of the partcpant. It was attached to cylnder B of TS. The reference temperature of the comparson was 20 C. For measurements performed at a temperature devatng from 20 C, the effectve area of the TS was referred to 20 C usng the pston-cylnder thermal expanson coeffcent. The zero-pressure ndcator was swtched n the pressure lne between the gauge of TS and the laboratory standard to control equalty of the pressures generated by the two pressure standards and to avod pressure gradents along the lne. A bypass lne wth a valve was connectng both sdes of the zero ndcator to set ts zero readng. The workng gas of the assembly was ether dry ar or ntrogen. The pston gauge of TS dd not have any electroncs to be warmed up, but the zero-pressure ndcator had to be swtched on at least 5 mnute before the measurements. The pston had to be n a self-centred stable floatng poston. The tme between a pressure level change and the acquston of the data correspondng to the equlbrum of the laboratory standard and TS had to be not shorter than 3 mnutes. The measurements ncluded four cycles each wth nomnal pressures generated n the followng order (100, 200, 500, 1000, 2000, , 4000, 5000, 5000, 4000, 3000, 2000, 1000, 500, 200, 100) Pa. Thus, 64 measurements were performed n total. 10

11 Table 10. Load of TS to generate the nomnal pressures Nomnal pressure Loads (pstons, weght carrer, weghts) n Pa 100 Pston "1 mbar" 200 Pston "2 mbar" 500 Pston "L" + weght carrer "L"+ "1 mbar" 1000 Pston "L" + weght carrer "L" + "1 mbar" + "5 mbar" 2000 Pston "L" + weght carrer "L" + "1 mbar" + "5 mbar" + "10 mbar" 3000 Pston "L" + weght carrer "L" + "1 mbar" + "5 mbar" + "20_1 mbar" 4000 Pston "L" + weght carrer "L" + "1 mbar" + "5 mbar" + "10 mbar" + "20_1 mbar" 5000 Pston "L" + weght carrer "L" + "1 mbar" + "5 mbar" + "20_1 mbar" +"20_2 mbar" Due to measurement set-up of the actual comparson, VNIIM was only able to measure a nomnal pressure value of 4800 Pa as the hghest achevable, see secton 2.2. The correspondng combnaton of mass peces for ths pont was the followng: Table 11. Load of TS to generate the nomnal pressure of 4800 Pa Nomnal pressure n Pa 4800 Loads (pstons, weght carrer, weghts) Pston "L" + weght carrer "L" + "2_1 mbar" +"2_2 mbar" + "20_1 mbar" +"20_2 mbar" Comparablty of the results obtaned at ths pressure wth the results of other partcpants obtaned at 5 kpa was determned by measurements of the plot laboratory performed at both nomnal pressures, 4.8 and 5 kpa. No addtonal/dfferent loads were appled to the pstons of TS. The equlbrum between TS and the laboratory standard, whch was controlled wth the help of the zero-pressure ndcator, was acheved by adjustng the pressure n the laboratory standard. In the case of VMC usng a smlar devce to TS as a laboratory standard, the zero ndcator readng was used as a drect measure of the pressure dfference between TS and LS. The effectve area of the TS determned for a partcular measurement (A p ) referred to 20 C was calculated wth the equaton A p = g m ( 1 ρa / ρ ) p [1 + ( α + β )( t 20 C) ], (1) where m are true masses of the pston, the weght carrer and the mass peces placed on the weght carrer of TS; ρ are denstes of the parts wth masses m ; 11

12 ρ a s ar densty; g s local gravty acceleraton; p s pressure generated by the laboratory standard at the TS reference level; (α+β) s thermal expanson coeffcent of the pston-cylnder; t s temperature of TS. The values of p and ρ a were calculated and t was measured by the partcpatng laboratory. The condtons of the measurements are gven n Annexes A1-A4. 6. RESULTS 6.1 Stablty of the transfer standard Pror to the comparson, the long-term stablty of TS has been evaluated on the bass of calbraton results obtaned n 1999, 2001 and 2007 (Fgs. 1 and 2). green 2001 blue 1999 magenta Pa Gauge pressure n Pa Fgure 1. Devatons from reference pressure, pstons "1 mbar" and "2 mbar" 12

13 green 2001 purple 1999 magenta Pa Gauge pressure n Pa Fgure 2. Devatons from reference pressure, pston "L" As a result, no systematc changes could be seen at pressures up to 200 Pa. At the maxmum pressure of 5 kpa, the dfferences dd not exceed 1 Pa. Durng the comparson, the TS stablty was measured by the plot laboratory at the begnnng, n the mddle, and at the end of the comparson. The results for the one year perod of the actual KC are presented n Fgs. 3 and 4. Durng the tme of the KC, the plot laboratory performed four measurements n March 2010, November 2009, March 2009 and November 2008, see table 8, although the measurements of November 2009 ncluded only one measurement at each pressure. Fgure 3. Stablty of the transfer standard. Mean effectve areas and ther standard uncertantes for pstons "1 mbar" and "2 mbar" measured by PTB n 2008, 2009 and

14 Fgure 4. Stablty of the transfer standard. Mean effectve areas and ther standard uncertantes for pston L measured by PTB n 2008, 2009 and 2010 From an analyss of these results, see Appendx C, t could be concluded that the TS dd not underlay any systematc tme drft durng ths KC. Therefore, no tme dependent correcton was appled to the partcpants' results, and the nstablty of TS, u nstab (A p ), was expressed by the standard devaton of the mean values of the effectve area measured by the plot lab for each nomnal pressure, A p,plot, mean : N Ap, plot, mean = Ap, plot, N, (2) = 1 N 2 ( A ) ( A A ) ( N ) u nstab p = p, plot, p, plot, mean 1, (3) = 1 where N s number of plot laboratory measurements over the tme of the comparson, N = 4. The uncertanty of the TS determned n ths way s presented n Table 12. Table 12. Uncertanty contrbuton due to nstablty of TS, u nstab (A p ) Nomnal pressure n Pa u nstab (A p ) / mm 2 u nstab (A p )/A p Fgures 5 and 6 show relatve devatons of the effectve areas measured by the plot four tmes from ther mean values. It s obvous that no systematc tme drft occurred. Hence, the partcpants values of the effectve areas dd not have to be corrected for the comparson

15 Fgure 5. Instablty of TS: Relatve devatons of the sngle effectve areas from the mean ones of pstons "1 mbar" and "2 mbar" measured by the plot laboratory together wth the relatve standard devatons of the mean values (uncertanty bars) extracted from Table 12 Fgure 6. Instablty of TS: Relatve devatons of the sngle effectve areas from the mean ones of pston "L" measured by the plot laboratory together wth the relatve standard devatons of the mean values (uncertanty bars) extracted from Table Results of the Partcpants The mean effectve areas, the standard devatons of the effectve areas at each pressure and the combned relatve standard uncertantes of the effectve areas reported by the partcpants are presented n Table 13. The graphcal vsualzaton of the summary of all measurement cycles and partcpants s gven n Fgures 7 and 8. 15

16 Table 13. Mean effectve areas (<A p >), ther relatve standard devatons (s(a p )/A p ), relatve standard uncertantes of pressures (u(p)/p), standard uncertanty of the temperature of TS (u(t)) and combned standard uncertantes of the mean effectve areas (u(a p )/A p ) Nomnal pressure n Pa <A p> s(a p)/a p u(p)/p u(t) / mm 2 / C PTB VNIIM CMI VMC u(a p)/a p <A p> s(a p)/a p u(p)/p u(t) / mm 2 / C u(a p)/a p <A p> s(a p)/a p u(p)/p u(t) / mm 2 / C u(a p)/a p <A p> s(a p)/a p u(p)/p u(t) / mm 2 / C u(a p)/a p 16

17 A p / mm PTB 2010 CMI VNIIM VMC p / Pa Fgure 7. Mean effectve areas and ther uncertantes of TS wth pstons "1 mbar" and "2 mbar" determned by the partcpants A p / mm PTB 2010 CMI VNIIM VMC p / Pa Fgure 8. Mean effectve areas and ther uncertantes of TS wth pston "L" determned by the partcpants 17

18 6.3 Reference value calculaton The key comparson reference value (KCRV), A p,ref, was calculated at each pressure pont of the KC as a weghted mean, see eq. (4), where N s the number of the ndependent partcpants, PTB, CMI and VNIIM, N = 3, A p, and u(a p, )/A p, are the mean effectve area and ts uncertanty, respectvely, of partcpant at pressure p. Ths method s recommended n [6] and has been used for evaluaton of numerous KCs. The uncertanty of A p,ref, u(a p,ref ), was calculated by equaton (5). The weghted mean method s applcable only when the results of the partcpants succeed the consstency check based on the ch-squared test. Results are consdered as consstent f the observed ch-squared value χ 2 obs calculated by (6) s smaller than the value of the ch-square dstrbuton calculated for degree of freedom ν = N 1 at probablty Pr = 0.05, χ 2 (ν, Pr). A p,ref = N u A 2 = 1 p, N 1 p, ( ) 2 A u ( A ) 1 (,ref ) = N Ap 2 1 u ( Ap, ) N ( ) 2 Ap, Ap,ref obs = 2 u ( A ) = 1 p, 0.5 = u (5) = 1 p, 2 χ (6) Note, that although four NMIs partcpated n ths KC, VMC was excluded from the evaluaton of KCRV snce t s traceable to PTB (calbraton certfcate Nr PTB 2009, see secton 2.4) and, hence, s statstcally not ndependent. Takng nto account the nstablty of TS as defned by equaton (3), the combned standard uncertanty of KCRV, u*(a p,ref ), s defned as u*(a p,ref ) = [u 2 (A p,ref ) + u 2 nstab(a p )] 0.5. (7) Results of the KCRV evaluaton are summarsed n Table 14. Table 14. Key comparson reference values (A p,ref ), ther relatve standard uncertantes (u(a p,ref )/A p,ref ) and uncertantes combned wth the TS nstablty (u*(a p,ref )/A p,ref ) as well as parameters of the ch-squared test χ 2 obs and χ 2 (ν, 0.05) Nomnal pressure n Pa A p,ref / u(a p,ref ) / u*(a p,ref ) / mm 2 A p,ref A p,ref χ 2 obs χ 2 (ν, 0.05), ν = N (4) 18

19 1) Ths value of A p,ref was obtaned by combnng the PTB and CMI results obtaned at 5 kpa wth the result of VNIIM obtaned at 4.8 kpa. Such a combnaton s justfed by the measurements of the plot laboratory performed at pressures of 4.8 and 5 kpa whch showed no dfference n the effectve area. The consstency check supports the choce of the KCRV calculaton as a weghted mean. 6.4 Degree of equvalence The degrees of equvalence of the partcpants n relaton to KCRV consdered at each pressure are expressed n terms of relatve devatons of the partcpants' results from KCRV ( A p, / A p,ref ) and relatve expanded uncertantes of these devatons (U( A p, ) / A p,ref ), the latter beng calculated as: ( Ap, ) Ap,ref = 2 [ u ( Ap, ) u ( Ap,ref ) unstab ( Ap )] Ap, ref U +. (8) Numercal data for the devatons and the uncertantes at all pressures are lsted n Table 15, a graphcal presentaton s gven n Appendx A. Table 15. Relatve devatons of the partcpants results from the reference values ( A p /A p ) and ther expanded uncertantes (U( A p /A p )) Nomnal pressure (Pa) Ap /Ap PTB VNIIM CMI VMC U( Ap/Ap) Ap /Ap U( Ap/Ap) Ap /Ap U( Ap/Ap) Ap /Ap U( Ap/Ap) The degrees of equvalence between the laboratores are expressed by relatve dfferences between them (d j ) and relatve expanded uncertantes of these dfferences (U(d j )) calculated as: d ( Ap, Ap, j ) Ap, ref ( d ) [ ( ) ( ) ( )] 1 j 2 u Ap, + u Ap, j unstab Ap Ap, ref j =, (9) U = +, (10) whch are presented n the tables n Appendx B. 19

20 7. DISCUSSION From Table 13, the approprate performance of TS can be deduced. At the hghest pressure of 5 kpa, the typcal relatve standard devatons of A p of the partcpants range from 5 ppm to 80 ppm. A comparson of s(a p ) wth u(a p ) ndcates that for most NMIs n ths KC a major contrbuton to the uncertanty of A p was from the uncertanty of the laboratory standards. Comparson of typcal s(a p )/A p values from laboratory to laboratory clearly show dfferent performance of pressure measurements of the partcpatng NMIs. Among the laboratores, all values of the reported A p agree wthn ther expanded uncertantes (k = 2). All NMIs agree wth the reference value of KC wthn ther expanded uncertantes (k = 2), (Table 15 and Appendx A). The most results agree wth the reference value of the KC even wthn ther standard uncertantes (k = 1) except for PTB at 5 kpa. A comparson n pars demonstrates that the results of all partcpants agree wth each other wthn ther expanded uncertantes (k = 2) at all pressures (Appendx B). Wthn 48 compared pars of results, n no case there s a dsagreement on the level of the expanded uncertantes (k = 2). Only n two cases there are dfferences between the laboratores whch are bgger than the standard uncertantes of these dfferences. 8. CONCLUSIONS The transfer standard was stable wthn Pa p n terms of pressure n the perod of the KC. For all laboratores all the results are equvalent wth the key comparson reference values wthn the expanded uncertantes (k = 2), all but one results even wthn the standard uncertantes. For the NMIs results compared n pars, all of totally 48 pars show agreement wthn the expanded (k = 2) and 46 pars wthn the standard uncertantes. Wth the results of ths comparson PTB, VNIIM, CMI and VMC have supported ther measurement capabltes stated n the KCDB of BIPM. REFERENCES [1] Pppg E., Uhthoff E., Tauchglocken-Manometer als Gruppen-Prmärnormal des ASMW zur Darstellung der Druckskala m Überdruckberech bs 1600 Pa, Fengerätetechnk, 24 (1975) 20 [2] Perkn M. et al., Comparson of European dfferental pressure standards n the range 3 Pa to 1000 Pa, Metrologa, 36 (1999) 1 [3] Oowa A., Novel Nonrotatonal Pston Gauge wth Weght Balance Mechansm for the Measurement of Small Dfferental Pressures, Metrologa, 30 (1994) 607 [4] Delajoud P., Grard M., Pston gauge wth electronc balance for calbratons of very low pressures, Vakuum n Forschung und Praxs, 15 (2003) 24 [5] Tesar J., Repa P., Prazak D., Krajcek Z., Peksa L., The new method of traceablty of a force balanced pston gauge used as prmary vacuum standard, Vacuum, 76 ( 2004) 491 [6] Cox M.G., The evaluaton of key comparson data, Metrologa, 39 (2002)

21 APPENDIX A. Relatve devatons of the partcpants results from the reference value ((Ap, - Ap,ref)/Ap,ref) wth the expanded uncertantes (k = 2) of these devatons at pressures 100 Pa to 5 kpa 100 Pa: 200 Pa: 21

22 500 Pa: 1000 Pa: 2000 Pa: 22

23 3000 Pa: 4000 Pa: 5000 Pa: 23

24 APPENDIX B. Relatve dfferences between the partcpants results () and ther expanded uncertantes () (k = 2) at pressures 100 Pa to 5 kpa j p = 100 Pa CMI VNIIM PTB VMC CMI VNIIM PTB VMC j p = 200 Pa CMI VNIIM PTB VMC CMI VNIIM PTB VMC j p = 500 Pa CMI VNIIM PTB VMC CMI VNIIM PTB VMC j p = 1 kpa CMI VNIIM PTB VMC CMI VNIIM PTB VMC

25 j p = 2 kpa CMI VNIIM PTB VMC CMI VNIIM PTB VMC j p =3 kpa CMI VNIIM PTB VMC CMI VNIIM PTB VMC j p = 4 kpa CMI VNIIM PTB VMC CMI VNIIM PTB VMC j p = 5 kpa CMI VNIIM PTB VMC CMI VNIIM PTB VMC

26 APPENDIX C. Stablty of TS The followng dscusson focuses by chance on the nomnal pressure pont of 3000 Pa as a typcal data pont, compare to the whole dataset of TS measurements at the plot laboratory n Fg. 12. It could be shown that a lnear ft wthout slope contanng all measurements wthout weghng procedure yelded a satsfyng result ndcatng the stablty of TS. Performng a weghted lnear ft y = ax + b the followng values for a and b were obtaned b = ± mm 2 and a = (2.36 ± 3.57) 10-5 mm 2 /day, see Fgure C1. Fgure C1. Weghted lnear ft of the mean effectve areas wth correspondng combned uncertantes at p = 3000 Pa measured n March 2010 (0 days), March 2009 and November 2008, assumng a tme drft. Due to the uncertantes of the ft, a constant zero slope was deduced as a consequence,.e. t was assumed that bascally the TS had no systematc drft durng the tme of the comparson. Fgure C2. Weghted lnear ft of the mean effectve areas wth correspondng combned uncertantes at p = 3000 Pa measured n March 2010 (0 days), March 2009 and November 2008, assumng no tme drft. On the other hand, the TS s assocated wth a tme-constant uncertanty to be evaluated. Performng a weghted lnear ft y = ax + b assumng a = 0 we obtaned 26

27 then b = ± mm 2 (Fgure C2). Ths analyss was compared wth a nonweghted lnear ft wth zero-slope, see Fgure C3. Fgure C3. Non-weghted lnear ft of the mean effectve areas wth correspondng combned uncertantes at p = 3000 Pa measured n March 2010 (0 days), March 2009 and November 2008, assumng no tme drft. We ganed b = ± mm 2. Fnally, we analysed all data ponts avalable,.e. we consdered the sngle pressure-value measurements n November 2009 for a non-weghted lnear ft wth zero-slope, see Fgure C4. Fgure C4. Non-weghted lnear ft of the mean effectve areas wth correspondng combned uncertantes at p = 3000 Pa measured n March 2010 (0 days), November 2009, March 2009 and November 2008, assumng no tme drft. The result s then b = ± mm 2. Hence, we realzed that the dfferent analyses wth zero-slope bascally yelded smlar results. Consequently, we performed non-weghted lnear fts wth the complete data set to estmate the uncertanty of TS due to ts nstablty. 27

EURAMET.M.D-S2 Final Report Final report

EURAMET.M.D-S2 Final Report Final report Fnal report on ERAMET blateral comparson on volume of mass standards Project number: 1356 (ERAMET.M.D-S2) Volume of mass standards of 10g, 20 g, 200 g, 1 kg Zoltan Zelenka 1 ; Stuart Davdson 2 ; Cslla

More information

Measurement of the Activity Concentration of the Radionuclide Am-241 in a Solution COOMET PROJECT 359/RU/06

Measurement of the Activity Concentration of the Radionuclide Am-241 in a Solution COOMET PROJECT 359/RU/06 Measurement of the Actvty Concentraton of the Radonuclde Am-4 n a Soluton COOMET PROJECT 359/RU/06 I.A. Khartonov, A.V. Zanevsky ), V. Mlevsk, A. Ivanukovch ), P. Oropesa and Y. Moreno 3) ), ) BelGIM,

More information

«Key comparisons in the field of multiples and submultiples of the kilogram»

«Key comparisons in the field of multiples and submultiples of the kilogram» FINAL REPORT: СООМЕТ.М.М-К5 «Key comparsons n the feld of multples and submultples of the klogram» (CООМЕТ 258/RU/2) V.S. Snegov, VNIIM, St.-Petersburg December 213 1. Introducton Ths report summarzes

More information

Final report. Absolute gravimeter Intercomparison

Final report. Absolute gravimeter Intercomparison Federal Department of Justce and Polce FDJP Federal Offce of Metrology METAS Baumann Henr 16.04.010 Fnal report Absolute gravmeter Intercomparson EURAMET Project no. 1093 Coordnator of the comparson Henr

More information

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

More information

1 Introduction. Final Report

1 Introduction. Final Report Fnal Report Fnal Report on EURAMET key comparson (EURAMET.M.M-K4. and EURAMET.M.M-K.) of 1 kg and submultples of the klogram standards n stanless steel (project code: EURAMET 110) 1 Introducton Ths report

More information

BIPM comparison BIPM.RI(II)-K1.Eu-155 of the activity measurements of the radionuclide 155 Eu. G. Ratel and C. Michotte BIPM

BIPM comparison BIPM.RI(II)-K1.Eu-155 of the activity measurements of the radionuclide 155 Eu. G. Ratel and C. Michotte BIPM BIPM comparson BIPM.RI(II)-K1.Eu-155 of the actvty measurements of the radonuclde 155 Eu G. Ratel and C. Mchotte BIPM Abstract In 1993, a natonal metrology nsttute, the NPL (UK), submtted a sample of known

More information

Experiment 1 Mass, volume and density

Experiment 1 Mass, volume and density Experment 1 Mass, volume and densty Purpose 1. Famlarze wth basc measurement tools such as verner calper, mcrometer, and laboratory balance. 2. Learn how to use the concepts of sgnfcant fgures, expermental

More information

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient 58:080 Expermental Engneerng 1 OBJECTIVE Lab 2e Thermal System Response and Effectve Heat Transfer Coeffcent Warnng: though the experment has educatonal objectves (to learn about bolng heat transfer, etc.),

More information

Chapter 13: Multiple Regression

Chapter 13: Multiple Regression Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to

More information

LNG CARGO TRANSFER CALCULATION METHODS AND ROUNDING-OFFS

LNG CARGO TRANSFER CALCULATION METHODS AND ROUNDING-OFFS CARGO TRANSFER CALCULATION METHODS AND ROUNDING-OFFS CONTENTS 1. Method for determnng transferred energy durng cargo transfer. Calculatng the transferred energy.1 Calculatng the gross transferred energy.1.1

More information

SPANC -- SPlitpole ANalysis Code User Manual

SPANC -- SPlitpole ANalysis Code User Manual Functonal Descrpton of Code SPANC -- SPltpole ANalyss Code User Manual Author: Dale Vsser Date: 14 January 00 Spanc s a code created by Dale Vsser for easer calbratons of poston spectra from magnetc spectrometer

More information

Regulation No. 117 (Tyres rolling noise and wet grip adhesion) Proposal for amendments to ECE/TRANS/WP.29/GRB/2010/3

Regulation No. 117 (Tyres rolling noise and wet grip adhesion) Proposal for amendments to ECE/TRANS/WP.29/GRB/2010/3 Transmtted by the expert from France Informal Document No. GRB-51-14 (67 th GRB, 15 17 February 2010, agenda tem 7) Regulaton No. 117 (Tyres rollng nose and wet grp adheson) Proposal for amendments to

More information

Evaluation of inconsistent data of key comparisons of measurement standards

Evaluation of inconsistent data of key comparisons of measurement standards ACTA IMEKO December 013, Volume, Number, 8 33 www.meko.org Evaluaton of nconsstent data of key comparsons of measurement standards A. Chunovkna, N. Zvagn, N. Burmstrova D.I.Mendeleyev Insttute for Metrology,

More information

Second Order Analysis

Second Order Analysis Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

Uncertainty and auto-correlation in. Measurement

Uncertainty and auto-correlation in. Measurement Uncertanty and auto-correlaton n arxv:1707.03276v2 [physcs.data-an] 30 Dec 2017 Measurement Markus Schebl Federal Offce of Metrology and Surveyng (BEV), 1160 Venna, Austra E-mal: markus.schebl@bev.gv.at

More information

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00 ONE IMENSIONAL TRIANGULAR FIN EXPERIMENT Techncal Advsor: r..c. Look, Jr. Verson: /3/ 7. GENERAL OJECTIVES a) To understand a one-dmensonal epermental appromaton. b) To understand the art of epermental

More information

Bilateral comparison. between. the FG5X-216 and the FG5-242

Bilateral comparison. between. the FG5X-216 and the FG5-242 Olver Francs 15.06.015 Internal Report Blateral comparson between the FG5X-16 and the FG5-4 Walferdange Underground Laboratory for Geodynamcs, Luxembourg Olver Francs - Unversty of Luxembourg (UL) - Luxembourg

More information

FUZZY FINITE ELEMENT METHOD

FUZZY FINITE ELEMENT METHOD FUZZY FINITE ELEMENT METHOD RELIABILITY TRUCTURE ANALYI UING PROBABILITY 3.. Maxmum Normal tress Internal force s the shear force, V has a magntude equal to the load P and bendng moment, M. Bendng moments

More information

Comparison of Regression Lines

Comparison of Regression Lines STATGRAPHICS Rev. 9/13/2013 Comparson of Regresson Lnes Summary... 1 Data Input... 3 Analyss Summary... 4 Plot of Ftted Model... 6 Condtonal Sums of Squares... 6 Analyss Optons... 7 Forecasts... 8 Confdence

More information

/ n ) are compared. The logic is: if the two

/ n ) are compared. The logic is: if the two STAT C141, Sprng 2005 Lecture 13 Two sample tests One sample tests: examples of goodness of ft tests, where we are testng whether our data supports predctons. Two sample tests: called as tests of ndependence

More information

Analytical Chemistry Calibration Curve Handout

Analytical Chemistry Calibration Curve Handout I. Quck-and Drty Excel Tutoral Analytcal Chemstry Calbraton Curve Handout For those of you wth lttle experence wth Excel, I ve provded some key technques that should help you use the program both for problem

More information

Statistics Chapter 4

Statistics Chapter 4 Statstcs Chapter 4 "There are three knds of les: les, damned les, and statstcs." Benjamn Dsrael, 1895 (Brtsh statesman) Gaussan Dstrbuton, 4-1 If a measurement s repeated many tmes a statstcal treatment

More information

PHYS 450 Spring semester Lecture 02: Dealing with Experimental Uncertainties. Ron Reifenberger Birck Nanotechnology Center Purdue University

PHYS 450 Spring semester Lecture 02: Dealing with Experimental Uncertainties. Ron Reifenberger Birck Nanotechnology Center Purdue University PHYS 45 Sprng semester 7 Lecture : Dealng wth Expermental Uncertantes Ron Refenberger Brck anotechnology Center Purdue Unversty Lecture Introductory Comments Expermental errors (really expermental uncertantes)

More information

x = , so that calculated

x = , so that calculated Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to

More information

Negative Binomial Regression

Negative Binomial Regression STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

CHAPTER IV RESEARCH FINDING AND DISCUSSIONS

CHAPTER IV RESEARCH FINDING AND DISCUSSIONS CHAPTER IV RESEARCH FINDING AND DISCUSSIONS A. Descrpton of Research Fndng. The Implementaton of Learnng Havng ganed the whole needed data, the researcher then dd analyss whch refers to the statstcal data

More information

HYDROMETER CALIBRATION: LINKING SIM NMIs TO THE EURAMET KEY COMPARISON REFERENCE VALUE OF EURAMET.M.D-K4

HYDROMETER CALIBRATION: LINKING SIM NMIs TO THE EURAMET KEY COMPARISON REFERENCE VALUE OF EURAMET.M.D-K4 Smposo de Metrología 00 7 al 9 de Octubre HYROMETER CALIBRATION: LINKING SIM NMIs TO THE EURAMET KEY COMPARISON REFERENCE VALUE OF EURAMET.M.-K4 Becerra L O, Lorefce S, Pennecch F Centro Naconal de Metrología

More information

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA 4 Analyss of Varance (ANOVA) 5 ANOVA 51 Introducton ANOVA ANOVA s a way to estmate and test the means of multple populatons We wll start wth one-way ANOVA If the populatons ncluded n the study are selected

More information

is the calculated value of the dependent variable at point i. The best parameters have values that minimize the squares of the errors

is the calculated value of the dependent variable at point i. The best parameters have values that minimize the squares of the errors Multple Lnear and Polynomal Regresson wth Statstcal Analyss Gven a set of data of measured (or observed) values of a dependent varable: y versus n ndependent varables x 1, x, x n, multple lnear regresson

More information

A LINEAR PROGRAM TO COMPARE MULTIPLE GROSS CREDIT LOSS FORECASTS. Dr. Derald E. Wentzien, Wesley College, (302) ,

A LINEAR PROGRAM TO COMPARE MULTIPLE GROSS CREDIT LOSS FORECASTS. Dr. Derald E. Wentzien, Wesley College, (302) , A LINEAR PROGRAM TO COMPARE MULTIPLE GROSS CREDIT LOSS FORECASTS Dr. Derald E. Wentzen, Wesley College, (302) 736-2574, wentzde@wesley.edu ABSTRACT A lnear programmng model s developed and used to compare

More information

Final Report on Bilateral Supplementary Comparison APMP.M.P-S5 in Hydraulic Gauge Pressure from 1 MPa to 10 MPa

Final Report on Bilateral Supplementary Comparison APMP.M.P-S5 in Hydraulic Gauge Pressure from 1 MPa to 10 MPa ASIA-PACIFIC METROLOGY PROGRAMME 10 MPa HYDRAULIC GAUGE PRESSURE BILATERAL SUPPLEMENTARY COMPARISON Comarson Identfer: APMP.M.P-S5 Fnal Reort on Blateral Sulementary Comarson APMP.M.P-S5 n Hydraulc Gauge

More information

Chapter 11: Simple Linear Regression and Correlation

Chapter 11: Simple Linear Regression and Correlation Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests

More information

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS Unversty of Oulu Student Laboratory n Physcs Laboratory Exercses n Physcs 1 1 APPEDIX FITTIG A STRAIGHT LIE TO OBSERVATIOS In the physcal measurements we often make a seres of measurements of the dependent

More information

Composite Hypotheses testing

Composite Hypotheses testing Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter

More information

N.G. Domostroeva, A.G. Chunovkina. with contribution by:

N.G. Domostroeva, A.G. Chunovkina. with contribution by: СOOMET.M.V-K1 key ntercomparson of lqud vscosty measurements N.G. Domostroeva, A.G. Chunovkna wth contrbuton by: Natala Domostroeva: D.I. Mendeleyev Insttute for Metrology, Russan Federaton Kalna Dmtrova:

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

Final Report of the SIM 60 Co Absorbed- Dose-to-Water Comparison

Final Report of the SIM 60 Co Absorbed- Dose-to-Water Comparison Fnal Report of the SIM 60 Co Absorbed- Dose-to-Water Comparson (KCDB Entry SIM.RI(I)-K4) May 2008 C. K. Ross and K. R. Shortt (NRC, Canada; plot laboratory) M. Sarav (CNEA, Argentna) A. Meghzfene (IAEA)

More information

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

The optimal delay of the second test is therefore approximately 210 hours earlier than =2.

The optimal delay of the second test is therefore approximately 210 hours earlier than =2. THE IEC 61508 FORMULAS 223 The optmal delay of the second test s therefore approxmately 210 hours earler than =2. 8.4 The IEC 61508 Formulas IEC 61508-6 provdes approxmaton formulas for the PF for smple

More information

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth

More information

ITTC - Recommended Procedures and Guidelines

ITTC - Recommended Procedures and Guidelines Page of 9 able of Contents UCERAIY AALYSIS - EXAMPLE FOR WAERJE PROPULSIO ES... 2. PURPOSE OF PROCEDURE... 2 2. EXAMPLE FOR WAERJE PROPULSIO ES... 2 2. est Desgn... 2 2.2 Measurement Systems and Procedure...

More information

Temperature. Chapter Heat Engine

Temperature. Chapter Heat Engine Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

Limited Dependent Variables

Limited Dependent Variables Lmted Dependent Varables. What f the left-hand sde varable s not a contnuous thng spread from mnus nfnty to plus nfnty? That s, gven a model = f (, β, ε, where a. s bounded below at zero, such as wages

More information

DETERMINATION OF UNCERTAINTY ASSOCIATED WITH QUANTIZATION ERRORS USING THE BAYESIAN APPROACH

DETERMINATION OF UNCERTAINTY ASSOCIATED WITH QUANTIZATION ERRORS USING THE BAYESIAN APPROACH Proceedngs, XVII IMEKO World Congress, June 7, 3, Dubrovn, Croata Proceedngs, XVII IMEKO World Congress, June 7, 3, Dubrovn, Croata TC XVII IMEKO World Congress Metrology n the 3rd Mllennum June 7, 3,

More information

Report of Key Comparison SIM.QM-K91

Report of Key Comparison SIM.QM-K91 Report of Key Comparson SIM.QM-K91 ph of phthalate buffer (nomnal ph ~4.01 at 25 ºC) Fnal Report F.B. Gonzaga, J.C. Das, K.W. Pratt, J. Waters, L. Dmtrova, M. Delgado, G.T. Canaza, R.O. Crstancho, L.A.

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

2016 Wiley. Study Session 2: Ethical and Professional Standards Application

2016 Wiley. Study Session 2: Ethical and Professional Standards Application 6 Wley Study Sesson : Ethcal and Professonal Standards Applcaton LESSON : CORRECTION ANALYSIS Readng 9: Correlaton and Regresson LOS 9a: Calculate and nterpret a sample covarance and a sample correlaton

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

International Comparison EUROMET.QM-K1c

International Comparison EUROMET.QM-K1c Internatonal Comparson EUROMET.QM-K1c Fnal Report A.M.H. van der Veen 1, G. Neuwenkamp 1, R. Oudwater 1, R.M. Wessel 1, J. Novak, J-F. Perrochet 3, A. Ackermann 3, A. Rakowska 4, L. Cortez 5, F. Das 5,

More information

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11)

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11) Gravtatonal Acceleraton: A case of constant acceleraton (approx. hr.) (6/7/11) Introducton The gravtatonal force s one of the fundamental forces of nature. Under the nfluence of ths force all objects havng

More information

UNIVERSITY OF TORONTO Faculty of Arts and Science. December 2005 Examinations STA437H1F/STA1005HF. Duration - 3 hours

UNIVERSITY OF TORONTO Faculty of Arts and Science. December 2005 Examinations STA437H1F/STA1005HF. Duration - 3 hours UNIVERSITY OF TORONTO Faculty of Arts and Scence December 005 Examnatons STA47HF/STA005HF Duraton - hours AIDS ALLOWED: (to be suppled by the student) Non-programmable calculator One handwrtten 8.5'' x

More information

The Geometry of Logit and Probit

The Geometry of Logit and Probit The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.

More information

BOOTSTRAP METHOD FOR TESTING OF EQUALITY OF SEVERAL MEANS. M. Krishna Reddy, B. Naveen Kumar and Y. Ramu

BOOTSTRAP METHOD FOR TESTING OF EQUALITY OF SEVERAL MEANS. M. Krishna Reddy, B. Naveen Kumar and Y. Ramu BOOTSTRAP METHOD FOR TESTING OF EQUALITY OF SEVERAL MEANS M. Krshna Reddy, B. Naveen Kumar and Y. Ramu Department of Statstcs, Osmana Unversty, Hyderabad -500 007, Inda. nanbyrozu@gmal.com, ramu0@gmal.com

More information

Laboratory 1c: Method of Least Squares

Laboratory 1c: Method of Least Squares Lab 1c, Least Squares Laboratory 1c: Method of Least Squares Introducton Consder the graph of expermental data n Fgure 1. In ths experment x s the ndependent varable and y the dependent varable. Clearly

More information

Ph 219a/CS 219a. Exercises Due: Wednesday 23 October 2013

Ph 219a/CS 219a. Exercises Due: Wednesday 23 October 2013 1 Ph 219a/CS 219a Exercses Due: Wednesday 23 October 2013 1.1 How far apart are two quantum states? Consder two quantum states descrbed by densty operators ρ and ρ n an N-dmensonal Hlbert space, and consder

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

Linear Regression Analysis: Terminology and Notation

Linear Regression Analysis: Terminology and Notation ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented

More information

Supplemental document

Supplemental document Electronc Supplementary Materal (ESI) for Physcal Chemstry Chemcal Physcs. Ths journal s the Owner Socetes 01 Supplemental document Behnam Nkoobakht School of Chemstry, The Unversty of Sydney, Sydney,

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

A PAPER CLOCK MODEL FOR THE CESIUM CLOCK ENSEMBLE OF TL

A PAPER CLOCK MODEL FOR THE CESIUM CLOCK ENSEMBLE OF TL A PAPER CLOCK MODEL FOR THE CESIUM CLOCK ESEMBLE OF TL Shnn Yan Ln and Hsn Mn Peng atonal Standard Tme and Frequency Lab., TL, Chunghwa Telecom Co., Ltd. o. Lane 55, Mn-Tsu Rd. Sec. 5, YangMe, TaoYuan,

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Chapter 8. Potential Energy and Conservation of Energy

Chapter 8. Potential Energy and Conservation of Energy Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal

More information

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law: CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and

More information

This column is a continuation of our previous column

This column is a continuation of our previous column Comparson of Goodness of Ft Statstcs for Lnear Regresson, Part II The authors contnue ther dscusson of the correlaton coeffcent n developng a calbraton for quanttatve analyss. Jerome Workman Jr. and Howard

More information

MEASUREMENT OF MOMENT OF INERTIA

MEASUREMENT OF MOMENT OF INERTIA 1. measurement MESUREMENT OF MOMENT OF INERTI The am of ths measurement s to determne the moment of nerta of the rotor of an electrc motor. 1. General relatons Rotatng moton and moment of nerta Let us

More information

THE NEED FOR MEASUREMENT

THE NEED FOR MEASUREMENT AUSTRALIAN CURRICULUM PHYSICS GETTING STARTED WITH PHYSICS THE NEED FOR MEASUREMENT The heart and sole of any scence s measurement and expermentaton. In scence, attempts are made to make sense of our physcal

More information

Here is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y)

Here is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y) Secton 1.5 Correlaton In the prevous sectons, we looked at regresson and the value r was a measurement of how much of the varaton n y can be attrbuted to the lnear relatonshp between y and x. In ths secton,

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

On the correction of the h-index for career length

On the correction of the h-index for career length 1 On the correcton of the h-ndex for career length by L. Egghe Unverstet Hasselt (UHasselt), Campus Depenbeek, Agoralaan, B-3590 Depenbeek, Belgum 1 and Unverstet Antwerpen (UA), IBW, Stadscampus, Venusstraat

More information

CinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Model. formulation. procedure

CinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Model. formulation. procedure nhe roblem-solvng Strategy hapter 4 Transformaton rocess onceptual Model formulaton procedure Mathematcal Model The mathematcal model s an abstracton that represents the engneerng phenomena occurrng n

More information

U-Pb Geochronology Practical: Background

U-Pb Geochronology Practical: Background U-Pb Geochronology Practcal: Background Basc Concepts: accuracy: measure of the dfference between an expermental measurement and the true value precson: measure of the reproducblty of the expermental result

More information

Laboratory 3: Method of Least Squares

Laboratory 3: Method of Least Squares Laboratory 3: Method of Least Squares Introducton Consder the graph of expermental data n Fgure 1. In ths experment x s the ndependent varable and y the dependent varable. Clearly they are correlated wth

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

See Book Chapter 11 2 nd Edition (Chapter 10 1 st Edition)

See Book Chapter 11 2 nd Edition (Chapter 10 1 st Edition) Count Data Models See Book Chapter 11 2 nd Edton (Chapter 10 1 st Edton) Count data consst of non-negatve nteger values Examples: number of drver route changes per week, the number of trp departure changes

More information

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850) hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In

More information

Statistical Inference. 2.3 Summary Statistics Measures of Center and Spread. parameters ( population characteristics )

Statistical Inference. 2.3 Summary Statistics Measures of Center and Spread. parameters ( population characteristics ) Ismor Fscher, 8//008 Stat 54 / -8.3 Summary Statstcs Measures of Center and Spread Dstrbuton of dscrete contnuous POPULATION Random Varable, numercal True center =??? True spread =???? parameters ( populaton

More information

Strength Requirements for Fore Deck Fittings and Equipment

Strength Requirements for Fore Deck Fittings and Equipment (Nov 2002) (Rev.1 March 2003) (Corr.1 July 2003) (Rev.2 Nov 2003) (Rev.3 July 2004) (Rev.4 Nov 2004) (Rev.5 May 2010) Strength Requrements for Fore Deck Fttngs and Equpment 1. General 1.1 Ths UR S 27 provdes

More information

Chapter 9: Statistical Inference and the Relationship between Two Variables

Chapter 9: Statistical Inference and the Relationship between Two Variables Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

More information

COOMET Supplementary Comparison on Fiber Optic Power Responsivity

COOMET Supplementary Comparison on Fiber Optic Power Responsivity COOMET Supplementary Comparson on Fber Optc Power Responsvty COOMET.PR-S6.01 Fnal Report Alexey Svetlchny, Vladmr Kravtsov All-Russan Research Insttute for Optcal and Physcal Measurements (VNIIOFI), Ozernaya

More information

Numerical Transient Heat Conduction Experiment

Numerical Transient Heat Conduction Experiment Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use

More information

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

Department of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution

Department of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution Department of Statstcs Unversty of Toronto STA35HS / HS Desgn and Analyss of Experments Term Test - Wnter - Soluton February, Last Name: Frst Name: Student Number: Instructons: Tme: hours. Ads: a non-programmable

More information

Problem Points Score Total 100

Problem Points Score Total 100 Physcs 450 Solutons of Sample Exam I Problem Ponts Score 1 8 15 3 17 4 0 5 0 Total 100 All wor must be shown n order to receve full credt. Wor must be legble and comprehensble wth answers clearly ndcated.

More information

Torsion Stiffness of Thin-walled Steel Beams with Web Holes

Torsion Stiffness of Thin-walled Steel Beams with Web Holes Torson Stffness of Thn-walled Steel Beams wth Web Holes MARTN HORÁČEK, JNDŘCH MELCHER Department of Metal and Tmber Structures Brno Unversty of Technology, Faculty of Cvl Engneerng Veveří 331/95, 62 Brno

More information

SINGLE EVENTS, TIME SERIES ANALYSIS, AND PLANETARY MOTION

SINGLE EVENTS, TIME SERIES ANALYSIS, AND PLANETARY MOTION SINGLE EVENTS, TIME SERIES ANALYSIS, AND PLANETARY MOTION John N. Harrs INTRODUCTION The advent of modern computng devces and ther applcaton to tme-seres analyses permts the nvestgaton of mathematcal and

More information

VERIFICATION OF THE CONVENTIONAL MEASURING UNCERTAINTY EVALUATION MODEL WITH MONTE CARLO SIMULATION

VERIFICATION OF THE CONVENTIONAL MEASURING UNCERTAINTY EVALUATION MODEL WITH MONTE CARLO SIMULATION ISSN 1726-4529 Int. j. smul. model. 4 (2005) 2, 76-84 Orgnal scentfc paper VERIFICATION OF THE CONVENTIONAL MEASURING UNCERTAINTY EVALUATION MODEL WITH MONTE CARLO SIMULATION Acko, B. & Godna, A. Unversty

More information

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS

More information

Color Rendering Uncertainty

Color Rendering Uncertainty Australan Journal of Basc and Appled Scences 4(10): 4601-4608 010 ISSN 1991-8178 Color Renderng Uncertanty 1 A.el Bally M.M. El-Ganany 3 A. Al-amel 1 Physcs Department Photometry department- NIS Abstract:

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information