title n.
Download
Skip this Video
Loading SlideShow in 5 Seconds..
Title PowerPoint Presentation
Download Presentation
Title

Loading in 2 Seconds...

play fullscreen
1 / 17

Title - PowerPoint PPT Presentation


  • 142 Views
  • Uploaded on

Koji Niita: RIST, Japan Hiroshi Iwase: GSI, Germany Tatsuhiko Sato, JAEA, Japan Yousuke Iwamoto, Norihito Matsuda, Yukio Sakamoto, Hiroshi Nakashima: JAEA, Japan Davide Mancusi, Lembit Sihver: Chalmers, Sweden. Title. P HI T S.

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about 'Title' - koen


An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
title

Koji Niita:RIST, Japan

Hiroshi Iwase: GSI, Germany

Tatsuhiko Sato, JAEA, Japan

Yousuke Iwamoto, Norihito Matsuda,

Yukio Sakamoto, Hiroshi Nakashima: JAEA, Japan

Davide Mancusi, Lembit Sihver: Chalmers, Sweden

Title

PHITS

Multi-Purpose Particle and Heavy Ion Transport code System

overview 1

1951 1983 1997 2001 2003

NMTC NMTC/JAERI NMTC/JAERI97 NMTC/JAM PHITS

(ORNL) High Energy Fission CG geometry JAM, GG JQMD, MCNP

Magnetic Field Gravity

PHITS = MCNP + JAM + JQMD

Transport Particle and Energy

Proton

Neutron

Meson

Barion

Nucleus

Photon

Electron

0 ~ 200 GeV

10-5 eV ~ 200 GeV

0 ~ 200 GeV

0 ~ 200 GeV

0 ~ 100 GeV/u

1 keV ~ 1 GeV

1 keV ~ 1 GeV

MCNP

Neutron, Photon, Electron

Transport by Nuclear Data

JAM

Hadron-Nucleus Collisions

up to 200 GeV

JQMD

Nucleus-Nucleus Collisions

by Molecular Dynamics

Geometry: CG and GG

Tally, Mesh and Graphic

Tally: Track, Cross, Heat, Star,

Time, DPA, Product, LET

Mesh: cell, r-z, xyz

Counter:

Graphic: ANGEL (PS generator)

Language and Parallelism

FORTRAN 77

MPI

Overviewof PHITS(Particle and Heavy Ion Transport code System) (1)

Overview (1)

External Field:Magnetic Field, Gravity

Optical and Mechanical

devices

slide3

(JQMD)

(Bertini)

(JQMD)

(Bertini)

(JQMD)

MCNP

with

nuclear data

Event Generator

MCNP

with

nuclear data

Overviewof PHITS(Particle and Heavy Ion Transport code System) (2)

photons

electrons

hadrons

neutrons

protons

nucleus

200 GeV

200 GeV

200 GeV

100 GeV/u

100 GeV

 JAM, Hadron cascade model 

JQMD

In progress

 GEM, Evaporation process 

1 GeV

 SPAR, ATIMA, Ionization process 

20 MeV

10 MeV/u

1 MeV

1 MeV

only transport

with dE/dx (SPAR, ATIMA)

1 keV

0 MeV

0 MeV/u

thermal

0 MeV

slide4

JAM code for Hadron Nucleus Collisions up to 200 GeV

JAM

Introducing JAM(Jet AA Microscopic Transport Model) Y. Nara et.al. Phys. Rev.C61 (2000) 024901

JAM is a Hadronic Cascade Model, which explicitly treats all established hadronic states including resonances with explicit spin and isospin as well as their anti-particles.

We have parameterized all Hadron-Hadron Cross Sections, based on Resonance Model and String Model by fitting the available experimental data.

Au+Au 200GeV/u in cm

jam 1

DDX of elementary process, p(12GeV/c) + p , by JAM

JAM (1)

JAM: Y. Nara et.al. Phys. Rev.C61 (2000) 024901

rapidity distribution

transverse momentum distribution

Data : Nucl. Phys. B69 (1974) 454

jqmd 1

JQMD code for Nucleus-Nucleus Collisions up to 100 GeV/u

JQMD-1

JQMD(Jaeri Quantum Molecular Dynamics) forSimulation of Nucleus-Nucleus Collisions

K. Niiita et.al. Phys. Rev. C52 (1995) 2620http://hadron31.tokai.jaeri.go.jp/jqmd/

56Fe 800 MeV/u on 208Pb

Analysis of Nucleus-Nucleus Collisions by JQMD

jqmd 2

Results of PHITS

JQMD code inPHITS

JQMD-2

Introducing JQMD in PHITS :H. Iwase et.al. J. Nucl. Sci.Technol.39 (2002) 1142

Neutron Spectra from Thick Target

new feature 1

Only one-body observables

Event Generator

Any observables

Event Generator Modeof PHITS

What are we doing in Monte Carlo calculations for particle transport ?

New Feature (1)
  • Solving one-body Boltzmann equation

by using the evaluated nuclear data.

  • MCNP type code

energy is conserved in average.

no correlations

  • Simulating real phenomena

by using event generators.

  • PHITS for high energy by JAM, JQMD.

treat all ejectiles of collisions.

energy and momentum are conserved

in each collision.

Observables beyond one-body quantities are often required.

Deposit energy distribution in a cell (pulse height tally).

LET distribution in a micro dosimetry,

Coincident experiments, ………

new feature 2

Event Generator Mode for low energy neutrons in PHITS

New Feature (2)

Neutron data + Special Evaporation Model

We use the channel cross sections and neutron energy spectrum

of the first neutron and assume the binary decay of recoiled nucleus.

capture

charged particle and photon decay

elastic

final state is uniquely determined

Neutron

channels

(n,n’)

charged particle and photon decay

after the first neutron emission

(n,Nn’)

all particle and photon decay

after the first neutron emission

By this model, we can determine all ejectiles (neutrons, charged particles,

nucleus and photons) with keeping energy and momentum conservation.

PHITS can transport all charged particle and nucleus down to zero energy

and estimate deposit energy without local approximation (kerma factor).

event generator mode 2
Event Generator Mode (2)

An example of Event Generator Mode

Neutron-induced semiconductor soft error (SEU: single event upset)

Deposit Energy Distribution by PHITS

Leakage of

recoil nucleus

from Si

Neutron 19 MeV

Kerma

Si

slide12

Benchmark test of HEAT : compared with KEK experiment

Heat

Exp. and Cal. by H. Ohnishi et.al.

Nucl. Instr. and Meth. A545 (2005) 88

12 GeV

MCNPX: calculated by N. Matsuda

slide13

forward

45°

90°

90 degree

120°

MCNPX

PHITS

Ep=12GeV Cu:1mmt

backward

DDX of Secondary Particles

space radiation

Neutron spectrum in atmosphere

Space Radiation

T. Sato and K. Niita; Radiat. Res. 166 (2006) 544

CREME96

PHITS with JENDL-HEfile

slide15

EXPACS :Excel-based Program for calculating

Atmospheric Cosmic-ray Spectrum

by T. Sato

new features

New Featuresof PHITS(Recent Developments, released by ver.2.12)

New Features

(1) Duct Source option

for Neutron Long Beam Line calculations

(2) Super Mirror function

for Neutron Optics: reflection on the wall

(3) Time Dependent Materials

for moving material, T0 chopper, …

(4) Time Dependent Magnetic Fields

Pulse Magnet for neutron optics,

Wobbler Magnet (AC magnet)

(5) Dumpall option

for re-calculation of whole process

(history tape)

(6) Event Generator Mode

a full correlated transport

for all particles and energies

(7) LET tally

dose and track length

as a function of LET (dE/dx)

(8) Deposit2 tally

correlation of the deposit energy

distribution of two regions

(9) Timer (clock of particle)

simulation of TOF with dE correlation

distribution

PHITSver.2.13 and ANGEL ver.4.35

have been released from JAEA.

PHITSver.2.13 and ANGEL ver.4.35

have been released from JAEA.

PHITSver.2.13 and ANGEL ver.4.35

have been released from JAEA.

PHITSver.2.13 and ANGEL ver.4.35

have been released from JAEA.

niita@tokai.rist.or.jp

RIST. Code Center: http://www.rist.or.jp

niita@tokai.rist.or.jp

RIST. Code Center: http://www.rist.or.jp

niita@tokai.rist.or.jp

RIST. Code Center: http://www.rist.or.jp

niita@tokai.rist.or.jp

RIST. Code Center: http://www.rist.or.jp

distribution

Distribution of PHITS