C#: Regne ut aktivitet til Radioaktiv kilde

[caid]

Jeg trenger hjelp til å finne den mest effektive metoden til å løse denne:

 

En industriell Iridium-192 isotop på 1000 GBq, har en halveringstid på ca 73,8 dager.

Dvs at styrken på isotopen etter 73,8 dager er halvert til 500 GBq.

 

Se mal under for hvordan jeg har tenkt funksjonen skal se ut.

 

const float halfLifeIr = 73.8F;

public float GetActivityIr(DateTime currentTime, DateTime newTime, float currentActivity )
{
 .....
 .....
 .....
 return newActivity;
}

 

Noen som greier å fullføre denne?

[dahwan]

Easy!

 

const float halfLifeIr = 73.8F;

public static float GetActivityIr(DateTime currentTime, DateTime newTime, float currentActivity)
{
   float daysPassed = (float)(newTime - currentTime).TotalDays;

   return currentActivity / (2 * halfLifeIr / daysPassed);
}

Endret 21. februar 2010 av dahwan

[caid]

Ja det var ganske EASY!!

Men din metode fungerer ikke

 

Men har nå funnet ut av det:

 

const float halfLifeIr = 73.8F;

public static float GetActivityIr(DateTime currentTime, DateTime newTime, float currentActivity)
{
 float daysPassed = (float)(newTime - currentTime).TotalDays;
 return (float)(currentActivity * (Math.Exp(-0.6931472 * daysPassed / halfLifeIr)));
 // 0,693 = log(2)
}

Endret 21. februar 2010 av caid

[caid]

Etter litt testing, ble den endelige metoden, som fungerer både frem og tilbake i tid, slik:

 

// Radioactive Decay Equation
public static double GetRadioactivity(DateTime initialTime, DateTime queryTime, double initialActivity, double halfLife)
{
 double elapsedTime = (queryTime - initialTime).TotalDays;
 return (initialActivity * (Math.Exp(-0.6931472 * elapsedTime / halfLife)));
 // Works also, but about 50% slower
 // return (initialActivity * Math.Pow(2.0, -elapsedTime / halfLife));
}

 

Takker dahwan for koden til elapsed time, elegant og enkelt

 

Om noen har en bedre løsning, så er dere velkommen til å poste den.

[dahwan]

definer "fungerte ikke"

 

fungerte perfekt under testene mine

[caid]

Skal prøve å forklare...

Setter vi f.eks elapsedTime til 147.6 (2 ganger halveringstiden)

skal mengden stråling fra isotopen halveres 2 ganger.

 

 

Eks med isotop på 1000 GBq:

Etter de første 73.8 dagene blir den halvert til 500 GBq,

etter 73.8 nye dager blir de resterende 500 GBq halvert til 250 GBq.

 

 

Bruker vi din formel, blir regnestykket slik:

1000 GBq / (2 * 73.8 / 147.6 ) = 1000 GBq

Altså det samme som utgangspunktet!!

[dahwan]

Til mitt forsvar kjørte jeg bare en test som fikk perfekt resultat

 

planla ikke å bruke helga mi på dette regnestykket