Towards near real- time data analysis and instrument calibration - - PowerPoint PPT Presentation

▶
towards near real time data analysis and instrument
SMART_READER_LITE
LIVE PREVIEW

Towards near real- time data analysis and instrument calibration - - PowerPoint PPT Presentation

Sep 27, 2022 456 likes •612 views

Towards near real- time data analysis and instrument calibration at the Necsa Neutron Strain Scanner D. Marais and A.M. Venter Radiation Science, Research and Development Division, The South African Nuclear Energy Corporation (Necsa) SOC

slide-1
SLIDE 1

Towards near real- time data analysis and instrument calibration at the Necsa Neutron Strain Scanner

  • D. Marais and A.M. Venter

Radiation Science, Research and Development Division, The South African Nuclear Energy Corporation (Necsa) SOC Limited

slide-2
SLIDE 2

Who we are and what we do

2

  • South African Fundamental

Atomic Research Installation

– A.k.a. SAFARI-1 Research Reactor

– 20MW Tank-in-pool – Commissioned in 1965 – Flux 4x1014 n.cm-2.s-1 – Medical isotopes – Irradiation services – Neutron activation analysis – Neutron beam line facilities

slide-3
SLIDE 3

Neutron Diffraction Facility (NDIFF)

3

Materials Probe for Internal Strain Investigations

MPISI

Powder Instrument for Transition in Structure Investigations

PITSI

Spotted hyena in Zulu Zebra in Sotho

slide-4
SLIDE 4

Neutron strain scanning

4

∆d

Monochromator Centre of rotation Beam apertures Sample Alignment

2.1 1.176 1.1765 1.177 1.1775 1.178

  • 60 -44 -28 -12 4

20 36 52 2.1 1.174 1.176 1.178 1.18 1.182

  • 60
  • 48
  • 36
  • 24
  • 12

12 24 36 48 60 2.1 1.175 1.176 1.177 1.178

  • 60
  • 48
  • 36
  • 24
  • 12

12 24 36 48 60

→

Q

slide-5
SLIDE 5

MPISI DAC

5

GumTree Gui Android remote control OpenInspire Camera Server Scanman Denex neutron detector ANSTO Histogram Server Beam monitor Monochromator – Si (331)

  • 6 Stages – Rotate, Tilt, X, Y,

Z, horizontal focus

  • λ = 1.67Å at 2θ_M = 83.5°

Incident beam aperture

  • 3 Stages – X, Y, Rotate
  • Hor: 0.3 – 5 mm, Ver: 0 – 20

mm Sample Table (Huber)

  • 5 Stages – X, Y, Z, Sample

Rotate, Detector Rotate

  • 250 kg and 250 mm travel

Diffracted beam aperture

  • 2 Stages – X, Y
  • Ver: 0 – 20 mm

Eularian Cradle 3 Stages – Omega, Phi, Chi Motion control cabinet Galil DMC2280 motor controllers Primary & Secondary Shutters Newport ¼ cradle 2 Stages – Phi, Chi

slide-6
SLIDE 6

Sample table Centre of Rotation

6

Focused sy translated som softzero set sy returned

OI Cam graphic sy

  • 1.499756

3.50 3.50

  • 1.499512

som 0.0 0.0 2.044 2.044 Edge positions 589.95, 871.28 647.82, 863.66 626.74, 842.16 596.55, 877.61 Centre 730.61 755.74 734.45 737.08

  • Align camera axes with

sample table axes

– Adjust focus for scaling factor

  • High definition camera

– 1392(h) x 1040(v) 6.45µm square pixels – Telecentric lens with distortion < 3µm

Rotation 45 90 All (aligned)

OI Cam graphic Edge positions 922.37, 1171.96 655.57, 880.58 272.92, 501.25 463.03, 729.16 Pin width (px) 249.59 225.01 228.33 266.13 Centre 1047.165 768.075 387.085 596.095 Corrected edges 891.48, 1202.85 612.39, 923.76 231.40, 542.77 440.41, 751.78

  • COR offset

calculation methods

– 3 Positions – 4 Positions

slide-7
SLIDE 7

Data acquisition strategies

7

Time based

User specifies

  • Fixed measurement time
  • Fixed sample positions

Moves to next fixed position Data acquisition Time reached ? More positions ? Peak fitting Data analysis No Yes Yes No User specifies

  • Max measurement time
  • Peak fit criteria
  • Fixed sample positions

Moves to next fixed position Data acquisition Fit criteria reached ? More positions ? Peak fitting Data analysis No Yes Yes No Max time reached ? No Yes

Statistics based

y = 0.0043x + 0.0002

0.005 0.01 0.015 0.02 0.025 0.03 1 2 3 4 5 6

Intensity / time [n/s] Gauge volume [mm3] y = 0.0063x - 0.0075

0.005 0.01 0.015 0.02 0.025 0.03 0.035 1 2 3 4 5 6

Peak / Background / Time [s-1] Gauge volume [mm3] y = 0.0045x + 0.0016

0.05 0.1 0.15 0.2 0.25 0.3 10 20 30 40 50 60

Intensity / time [n/s] Gauge volume [mm3] y = 0.0068ln(x) + 0.0005

0.005 0.01 0.015 0.02 0.025 0.03 10 20 30 40 50 60 Intensity / time [n/s]

Gauge volume [mm3]

Aluminium Titanium

slide-8
SLIDE 8

Statistics based DAQ implementation

8

ANSTO Histogram Server Histmem Source Fit Beam monitor Export

slide-9
SLIDE 9

Test case – Welded Iron sample

9

Const time Min intensity Stdev: 4e-5 Stdev: 1.5e-5 Max time 600 600 1200 1800 Min intensity to background N/A 3 3 3 Min intensity N/A 100 10 10 Min % certainty intensity N/A 95 95 95 Min % certainty position N/A 99.993 99.993 99.9987 Corresponding stdev d-value N/A 4.00E-05 4.00E-05 1.50E-05

5.20 2.10 8.31 11.41 14.52

60

  • 60

Weld collar Secondary aperture Primary aperture Welded iron sample Gauge volume 3x3x3

slide-10
SLIDE 10

Results

10

100 200 300 400 500 600 700

  • 60
  • 40
  • 20

20 40 60 500 1000 1500 2000 2500

  • 60
  • 40
  • 20

20 40 60 0.00E+00 2.00E-05 4.00E-05 6.00E-05 8.00E-05 1.00E-04

  • 60
  • 40
  • 20

20 40 60 StDev d-value 1.1789 1.1791 1.1793 1.1795 1.1797 1.1799 1.1801

  • 60
  • 40
  • 20

20 40 60 d-value

Const time

100 200 300 400 500 600 700

  • 60
  • 40
  • 20

20 40 60 20 40 60 80 100 120 140

  • 60
  • 40
  • 20

20 40 60

00 05 05 5 05 04

  • 60
  • 40
  • 20

20 40 60 StDev d-value 789 791 793 795 797 799 801

  • 60
  • 40
  • 20

20 40 60 d-value

Min intensity

200 400 600 800 1000 1200 1400

  • 60
  • 40
  • 20

20 40 60 20 40 60 80 100

  • 60
  • 40
  • 20

20 40 60

00 05 05 0E-05 05 04

  • 60
  • 40
  • 20

20 40 60 StDev d-value 789 791 793 795 797 799 801

  • 60
  • 40
  • 20

20 40 60 d-value

Stdev 4e-5

500 1000 1500 2000

  • 60
  • 40
  • 20

20 40 60 2.10 5.20 8.31 11.41 14.52 200 400 600 800 1000 1200 1400

  • 60
  • 40
  • 20

20 40 60 2.10 5.20 8.31 11.41 14.52

+00 05 05 05 05 04

  • 60
  • 40
  • 20

20 40 60 StDev d-value 2.10 5.20 8.31 11.41 14.52 789 791 793 795 797 799 801

  • 60
  • 40
  • 20

20 40 60 d-value 2.10 5.20 8.31 11.41 14.52

Stdev 1.5e-5 d-val Stdev d-val Intensity Time

60 2.10 5.20 8.31 11.41 14.52 60 2.10 5.20 8.31 11.41 14.52

2.10 5.20 8.31 11.41 14.52 2.10 5.20 8.31 11.41 14.52

slide-11
SLIDE 11

0.00E+00 1.00E-05 2.00E-05 3.00E-05 4.00E-05 5.00E-05 6.00E-05 7.00E-05

  • 60
  • 40
  • 20

20 40 60 stdev of d-values at 2.10 500 1000 1500 2000 2500

  • 60
  • 40
  • 20

20 40 60 Intensity at 2.10 100 200 300 400 500 600 700

  • 60
  • 40
  • 20

20 40 60 Time at 2.10

Comparison

11

Stdev d-val Intensity Time 2.10

C Mi s s

Time (h) Const time 5.17 Min intensity 0.36 stdev 4e-5 0.15 stdev 1.5e-5 1.41

0.00E+00 1.00E-05 2.00E-05 3.00E-05 4.00E-05 5.00E-05 6.00E-05 7.00E-05

  • 60
  • 40
  • 20

20 40 60 stdev of d-values at 8.31 200 400 600 800 1000

  • 60
  • 40
  • 20

20 40 60 500 1000 1500 2000

  • 60
  • 40
  • 20

20 40 60 Time at 8.31

8.31

C Mi s s

Time (h) Const time 5.17 Min intensity 2.09 stdev 4e-5 0.85 stdev 1.5e-5 9.75

50 100 150 200 250 300 350

  • 60
  • 40
  • 20

20 40 60 Intensity at 14.52 500 1000 1500 2000

  • 60
  • 40
  • 20

20 40 60 Time at 14.52 0.00E+00 2.00E-05 4.00E-05 6.00E-05 8.00E-05 1.00E-04

  • 60
  • 40
  • 20

20 40 60 stdev of d-values at 14.52

14.52

C Mi s s

Time (h) Const time 5.17 Min intensity 5.18 stdev 4e-5 5.71 stdev 1.5e-5 15.51

slide-12
SLIDE 12

Conclusions and future projects

12

  • Automated, ‘loggable’ COR alignment achieved
  • Statistics based DAQ can save time, but is not a

black box solution

  • Intelligent DAQ protocol

– Initial measurement on course mesh – Identify regions with high stress gradients – Re-measure these areas with finer mesh – Automate the process

slide-13
SLIDE 13

Acknowledgements

  • Vladimir Luzin for the entry curve software ‘Alihn’
  • Rainer Schneider and Stefan Flemming for

OpenInspire and telecentrical camera idea

  • ANSTO/PSI for the use of SICS and Gumtree
  • National Research Foundation of South Africa
  • Necsa SOC Limited
  • Nobugs for bringing the people together!

13

slide-14
SLIDE 14

Questions

14