Message: Re: Concentric cylinders and energy deposition Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Concentric cylinders and energy deposition 

Forum: Geometry
Re: Question Concentric cylinders and energy deposition (Gabriel Sawakuchi)
Date: 18 Oct, 2005
From: michel maire <michel maire>

User Gabriel Sawakuchi wrote:

>> Hello,
>> 
>> I am a new user and I am not sure if this is the right forum for my
>> question. Anyway, I would like to get the total energy deposited by
>> delta rays in concentric cylinders when an ion pass through the axis of
>> the inner cylinder.
>> 
>> I have modified the Geometry of example Testm7 but I think I am doing
>> something wrong.
>> 
>> Case 1: The energy deposited
>> in the inner cylinder (tally1 0) 
>> (RMin=0, RMax=1000 nm,Height = 4000 nm) 
>> is bigger than that deposited in the
>> outer cylinder (tally2 0)
>> (RMin=0, RMax=1000 nm, Height = 4000 nm). 
>> It should not be correct since the volume of
>> the outer cylinder is larger and therefore it is expected to
>> happen more energy deposition events.

  Tally1 and Tally2 have the same mother: lAbsor. 
  They MUST NOT overlap (see user guide). 
  In your case, Tally2 must be from R1max to R2max
>>              __________
>>             |          |
>>             |__________| 
>>             |(x1,y1,z1)|
>>     --ion-->|(x2,y2,z2)|
>>             |__________|
>>             |          |
>>             |__________|
>> 
>> ----------------------------------------------------------
>> 
>> Cumulated Doses : Edep; Edep/Ebeam; Dose; Mass
>> 
>> tally1 0: 888.22 keV; 0.0088822 %; 4.8624 Gy; 2.9267e-08 mg
>> 
>> tally2 0: 31.704 keV; 0.00031704 %; 0.04339 Gy; 1.1707e-07 mg
>> 
>> ----------------------------------------------------------
>> 
>> Case 2: When I set two identical cylinders 
>> (RMin=0, RMax=1000 nm,Height = 4000 nm)
>> (same dimension and position) I get 0 eV for energy 
>> deposition in one of them.
>>              __________ 
>>             |(x1,y1,z1)|
>>     --ion-->|(x2,y2,z2)|
>>             |__________|

see above

>> 
>> ----------------------------------------------------------
>> 
>> Cumulated Doses : Edep; Edep/Ebeam; Dose; Mass
>> 
>> tally1 0: 899.09 keV; 0.0089909 %; 4.9219 Gy; 2.9267e-08 mg
>> 
>> tally2 0: 0 eV; 0 %; 0 Gy; 2.9267e-08 mg
>> 
>> -----------------------------------------------------------
>> 
>> Case 3:. Another strange thing happen when I define only one cylinder
>> (RMin=0, RMax=1000 nm,Height = 4000 nm).
>> The energy deposited is different from the case when I have
>> two cylinders with the same dimensions (Case 2).
>>              __________ 
>>             |          |
>>     --ion-->|(x1,y1,z1)|
>>             |__________|

    this case is correct
>> 
>> ------------------------------------------------------------
>> 
>> Cumulated Doses : Edep; Edep/Ebeam; Dose; Mass
>> 
>> tally1 0: 971.71 keV; 0.0097171 %; 5.3195 Gy; 2.9267e-08 mg
>> 
>> ------------------------------------------------------------
>> 
>> Case 4: On the other hand, when I have one cylinder in front of the other
>> through the path of the ion I get reasonable results but different
>> from the result I get in Case 3 where the first cylinder has the same
>> dimensions and position.
>>              __________ __________
>>             |          |          |
>>     --ion-->|(x1,y1,z1)|(x2,y2,z2)|
>>             |__________|__________|
>> 
    this case is correct too, but the geometry may be different of case 3.
    What is the tally material ?

How many events in a run ?

>> -------------------------------------------------------------
>> 
>> Cumulated Doses : Edep; Edep/Ebeam; Dose; Mass
>> 
>> tally1 0: 941.18 keV; 0.0094118 %; 5.1523 Gy; 2.9267e-08 mg
>> 
>> tally2 0: 946.59 keV; 0.0094659 %; 5.1819 Gy; 2.9267e-08 mg
>> 

  By the way, TestEm2 is also close to your problem.

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

 Add Message Add Message
to: "Re: Concentric cylinders and energy deposition"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]