ODL-Sonde DLM1420-2200 55/2 Schaltung

[DL8BCN]

Das Thema Eigennulleffekt hatten wir an anderer Stelle schon diskutiert.
Leider ist der sehr schwierig zu bestimmen.
Weil selbst in einer sehr massiven Bleiburg immer noch kosmische Strahlung zu messen sein wird.
Bleibt nur der berühmte Felsenkeller in Dresden.

[Peter-1]

Jetzt frage ich mich wie die von mir gemessenen CPM bei 0 µSv/h bezeichnet werden.
Ist zwar nur eine Extrapolation, aber der Korrelationskoeffizient sieht ja nicht schlecht aus.

[Radiohörer]

    
...wie wäre es, wenn wir hier etwas mehr Realtität wagen würden!?
Seid mir nicht böse wg. der offenen Worte.
Ich finde, die, anscheinend für besorgte Bürger oder einem "ruhigen Sorgentelefon" geschuldetem, geschönten und bis zur Unkenntlichkeit gemittelten Messwerte des BfS, völlig für´n Arsch!
Jeder der sich ein 70031 in den Garten stellt und die Minuten oder 10 Minuten Werte in einer Kurve für einen Tag/Woche/Monate plotten läßt, kann es sehen! Um 100nSv/h springende Werte sind VÖLLIG NORMAL!
Schaut einfach mal bei der EU nach: https://remap.jrc.ec.europa.eu/Advanced.aspx , erweitert den Radius auf 50 km und probiert ein paar Sonden durch.

Und: wer stört sich jetzt noch daran, obs zwei oder drei Nachkommastellen hat???
Eigentlich ist die Einheit "µ" (Mikro) falsch gewählt: "n" (Nano) wäre die Richtigere...
Wahrscheinlich lullt ein Wert, der nach einer "0" vor dem "," steht mehr ein, als eine reelle Zahl

[Kermit]

Jetzt frage ich mich wie die von mir gemessenen CPM bei 0 µSv/h bezeichnet werden.

"intrinsischer Nulleffekt"

Ich würde mir diese Bezeichnung aus der NUK "leihen"...

[DL8BCN]

Mein Programmierer hat mir den Arduino Sketch nun angepasst, das ich auch über längere Zeiträume Impulse zählen kann
In Zeile 4 gibt es nun ,,Total Counts".
Da der Sketch nun um einiges verändert wurde, spricht wohl nichts mehr aus urheberrechtlicher Sicht dagegen, es hier zu posten:

Code Auswählen Erweitern

#include <LiquidCrystal_I2C.h>
#include <Wire.h>

unsigned long  counts;      //variable for GM Tube events
unsigned long  totalCounts; //Total counts since start
unsigned long previousMillis; //variable  for measuring time
float averageCPM;
float sdCPM;
int currentCPM;
float  calcCPM;
LiquidCrystal_I2C lcd(0x27, 20, 4);
float CPMArray[100];

#define  LOG_PERIOD 60000 // count rate (in milliseconds)

void setup() { //setup
  counts = 0;
  totalCounts = 0;
  currentCPM = 0;
  averageCPM = 0;
  sdCPM = 0;
  calcCPM  = 0;
  lcd.init();
  lcd.backlight();
  Serial.begin(9600);
  pinMode(2,  INPUT);
  attachInterrupt(digitalPinToInterrupt(2), impulse, FALLING); //define  external interrupts
}

void loop() { //main cycle
  lcd.setCursor(0,2);
  lcd.print("CPM Count: ");
  lcd.print(counts);
  lcd.setCursor(0,3);
  lcd.print("Total CPM Count: ");
  lcd.print(totalCounts);
  unsigned long currentMillis  = millis();
  if (currentMillis - previousMillis > LOG_PERIOD) {
    previousMillis  = currentMillis;
    CPMArray[currentCPM] = counts * 1;
    lcd.clear();
    lcd.setCursor(0,0);
    lcd.print("uSv/hr: ");
    lcd.print(outputSieverts(CPMArray[currentCPM]));
    counts = 0;
    averageCPM = 0;
    sdCPM = 0;
    //calc avg and sd
    for (int x=0;x<currentCPM+1;x++)  {
      averageCPM = averageCPM + CPMArray
[x];
    }
    averageCPM = averageCPM / (currentCPM + 1);
    for (int x=0;x<currentCPM+1;x++)  {
      sdCPM = sdCPM + sq(CPMArray
[x]- averageCPM);
    }
    sdCPM  = sqrt(sdCPM / currentCPM) / sqrt(currentCPM+1);

    Serial.println("Avg:  " + String(averageCPM) + " +/- " + String(sdCPM) + "  ArrayVal: " + String(CPMArray[currentCPM]));
    currentCPM = currentCPM + 1;
    displayAverageCPM();
  }
}

void  impulse() {
  counts++;
  totalCounts++;
}
void displayAverageCPM()  {
  lcd.setCursor(0,1);
  lcd.print("Avg: ");
  lcd.print(outputSieverts(averageCPM));
  lcd.print("+/-");
  lcd.print(outputSieverts(sdCPM));
}
float outputSieverts(float x)  {
  float y = (x-60) * 0.001284; //fuer Vacutec 70031
  return y;
}


[TritiumLicht]

Danke für das Sketch!
Weiter oben war doch die Frage, wie man die uSv mit 3 Nachkommastellen ausgibt (die bereits diskutierte Sinnhaftigkeit lasse ich außer Acht)
Laut Doku: Die Routine lcd.print gibt bei floats per Voreinstellung 2 Nachkommastellen aus.
Wenn man als Parameter die ", 3" mit angibt, sollten 3 Nachkommastellen ausgegeben werden.
Also die Zeile (ziemlich am Ende) wie folgt ändern: lcd.print(outputSieverts(averageCPM), 3);
(Bitte ausprobieren - ich habe selbst keinen Arduino, kann es daher nicht testen)
HTH!

[silfox]

Bitte den Ueberlauf vermeiden, das das Array nur 100 Elemente groß ist:

float CPMArray[100];

currentCPM = currentCPM + 1;

Man könnte einfügen:

if (currentCPM>=100) currentCPM=0;

Vielleicht aber wäre 60 zu bevorzugen, das wäre dann ein 10min Mittelwert.

[silfox]

Falsch: 10 x 1min = 10min

[DL8BCN]

Danke für die Info.
Ich habe selber sehr wenig Ahnung von Arduino Programmierung.
Finde aber, das man sich da reinfuchsen kann.
Zur Not frage ich meinen Sohn.
Der ist Fachinformatiker für Anwendungsentwicklung.
Aber so ein Forum ist natürlich super für den Austausch.
Schönen Sonntag!
Rainer

[DL8BCN]

@Tritium: Die Befehlszeile >lcd.print(outputSieverts(averageCPM), 3);< funktioniert. Der Avg-Wert wird nun mit 3 Stellen hinter dem Komma ausgegeben.
Beim Average Wert macht das auch eher Sinn als beim CPM-Wert, der nach einer Minute ausgegeben wird.

[DL8BCN]

Nun mal ein bisschen Praxis:
Ich war eben zu Vergleichsmessungen an der ODL-Sonde in Selsingen.
Dort wurde um 17:00 Uhr der Wert 0,071 µSv/h auf der ODL-Seite angezeigt.
Hier die Messungen:
17:00 Uhr bis 17:10 Uhr : 808 Imp.  =80,8 cpm entspricht 0,076 µSv/h
17:10 Uhr bis 17:20 Uhr : 1579 Imp. =77,1 cpm entspricht 0,071 µSv/h
17:20 Uhr bis 17:30 Uhr : 2377 Imp. =79,8 cpm entspricht 0,075 µSv/h
17:30 Uhr bis 17:40 Uhr : 3058 Imp. =68,1 cpm entspricht 0,057 µSv/h

Der Durchschnitt für diese 40 Minuten war also 3058/40 Min. =76,45 cpm entsprechend 0,069 µSv/h.
Meine Berechnungsformel lautet aktuell ODL=(cpm-30) x 0,00150 (Kehrwert = 666,66)

Am Ende meiner Messung zeigte die ODL im Internet für Selsingen nur noch 0,068 µSv/h an.
Wenn ich auch da den Mittelwert bilde, ergibt das (0,071 + 0,068)/2= 0,0695 µSv/h.

Eine viel bessere Übereinstimmung zu meiner Durchschnittsmessung (0,069 µSv/h) kann es wohl kaum geben.
Nun wird es interessant, ob sich diese Messgenauigkeit an einer anderen ODL-Station wiederholen läßt.

[NoLi]

...
Eine viel bessere Übereinstimmung zu meiner Durchschnittsmessung (0,069 µSv/h) kann es wohl kaum geben.
Nun wird es interessant, ob sich diese Messgenauigkeit an einer anderen ODL-Station wiederholen läßt.

Deine Messungen in allen Ehren, aber: wie hoch sind denn die tatsächlichen Messwert-Abweichungen (Energie- und Richtungsabhängigkeit, bezogen auf Cs-137) der einzelnen BfS-Sonden vom physikalischen Sollwert? Welche Spannbreiten haben die Toleranzgrenzen dieser Sonden?
Vielleicht kann @silfox etwas dazu sagen.

Norbert

[DL8BCN]

Meine Kalibrierformel beruht auf ,, Try and Error".
Den intrinsischen Nullwert von 30cpm habe ich laut Angaben der Fa. Vacutec gewählt.
Die hatten das im Felsenkeller in Dresden ermittelt und stellt einen durchschnittlichen Richtwert dar. Das hatte mir Vacutec damals auf Anfrage mitgeteilt. Dazu gibt es an anderer Stelle hier einen Thread.
Die 60 cpm aus dem Datenblatt bzw. aus früheren Posts hier erscheinen mir viel zu hoch, da ich dann ja über 50% der CPM Zählrate in Abzug bringe. Ebenso passt der Wert 840 bei mir nicht. Es muss eher 660 sein.
Den Multiplikator 0,00150 ( Kehrwert ist ca. 660) habe ich dann so gewählt, daß der Meßwert von ca. 0,07μSv/h passt.
Also ziemlich empirisch ermittelt.
Ob das eine erlaubte Methode ist, wird uns Peter vermutlich sagen, der das ganze mathematisch, wissenschaftlich angegangen ist.
Ob der Kalibrierwert halbwegs passt, sieht man vermutlich erst, wenn man mehrere ODL Sonden mit möglichst unterschiedlichen Messwerten besucht und Vergleichsmessungen macht.

[Peter-1]

Das sieht doch schon super aus ! 
Jetzt noch an andere Orte gehen und messen, dann wird es immer besser. Laß Dich nicht entmutigen, ich habe auch so angefangen. Je größer die unterschiedlichen ODL Werte an den Sonden sind um so besser wird Deine Einordnung werden.

Grüße
Peter