L.P. Csernai Bergen

(1)

BCPL Presentation -Csernai 1

L.P. Csernai

Bergen

Computational Physics Lab.

(2)

Founding BCPL

• 1999, Joint application by the Theoretical and Computational Physics Section and the

Para//ab High Performance Computing Lab to EU for establishing a Research

Infrastructure

• Basis: Theoretical Reaction Modeling

Capacity lags behind experimental resources (CERN, BNL, FNAL) hindering research advance

• 2000, The EC supports BCPL as one of the 11 new Research Infrastructures!

• 2001, EC awards 2

^nd

RI grant to BCPL for

(3)

€ 813 000.-

(w/ MCTS)

(4)

Parallab

Unifob

(5)

Spring 2002: 30 visitors

Since 2000: 44 Research Projects  100 visitors

(6)

Countries of the users' home institutions

0 10 20 30 40 50 60 70 80

Home institution

Number of users

EU member states EU Assoc. States Other

Total

(7)

(8)

Extension Proposal to EU,

February 2001

(9)

(10)

(11)

(12)

Highlights of BCPL achievements in 2001

Up to now BCPL has accepted 44 European projects for support from the beginning of yr. 2000. Near to 100 researchers visited BCPL. The capacity to receive researchers from the EU was increased by 30% (extra offices and workstations) compared to the one offered in the contract in the beginning of 2001. Although we did not receive additional financing from the EU to cover the excess demand, we accommodate 10-20% more visitors, than the contract allows, so we might run out of EU financing before the end of the contract period.

In this report only a small fraction of the projects can be mentioned with important and recent results. The projects are grouped in thematic classes, as some of them are strongly related.

(13)

NUCLEAR WASTE MANAGEMENT

In the project "Transmutation of Heavy Nuclei" (project no. 3) from Calabria, Italy, led by L. Jenkovszky has completed its activity. It made an important step forward in "burning" radioactive nuclear waste to short lived and non- radioactive isotopes. Several core configurations were tested on our

supercomputer, exploiting parallel processing. As a side effect with come core configurations 1000 times the invested energy could be recovered from the forced decay of radioactive waste nuclei.

Another project, addressing topics related to nuclear waste, "Optimization of the Th-U Fuel Cycle for Molten-Salt Reactor" (project no. 30) from CEA, Saclay, France led by D. Ridikas, worked also on optimization of reactor core and fuel configurations. The aim of the project was to check the advantages of parallel processing in these Monte-Carlo codes. The project was successful, most codes could be converted by using local help, and most test runs

completed successfully.

(14)

SUBATOMIC REACTIONS

Project no.7 & 20 "Event Generator for A+A Collisions at RHIC"

from Ins. fuer Theoretische Phys., Uni. Frankfurt, Gremany led by M. Gorenstein had the ambitious goal to unify or match parton cascade, hydrodynamic and hadron cascade models for the description of A+A reactions. At the present stage physical

assumptions on how to match different models were tested on the computer, to see what are the consequences of different assumptions andwhich can be realized in a computational model at all. Most

assumptions lead to excessive computational task, and it is not trivial

to find an accurate but computationally well realizable model.

(15)

Project: BCPL07 (Univ. of Frankfurt):

Pb+Pb 160 GeV / nucleon, UrQMD model

Heavy Ion Reactions at CERN-SPS

[Stöcker]

(16)

(17)

(18)

Project no. 29 "Hadronization and Freeze Out in Fluid Dynamics" from Uni. Lisbon, Portugal, led by V.K. Magas, had similar aim of coupling different modules of reaction

models. The model was successful in working out the initial stage and its coupling to the CFD module. The publication reporting this work is submitted to Phys. Rev. D for publication already. The final module and the problem of Entropy production at this stage is still under intensive study.

Project no. 16 "Strangeness production in heavy ion collisions" from Univ. of Oulu, Finland led by A. Keranen, works exactly on the same part of the problem. A highly parallelized hadronization code was worked out, which can be used in connection with a preceding CFD module. Several publications are expected soon based on these results, and this will lead to the evaluation of measurable hadron data in a detailed and accurate reaction model.

A specific problem related to the same topic is studied in project no. 15 "Lambda (1520) production in NN collisions" from Uni. Marburg, Germany led by F. Puehlhofer. This particle has the same energetic properties as most other hadrons but its formation cross section is much smaller. Consequently it can make a difference between coalescence type of reaction models and thermal equilibrium models. The project works intensively on evaluating the consequences of the different model assumptions in different reaction models and

compare the model results to the experimental data measured by their collaboration. This project illustrates one of the tightest collaborations between theory and experiment.

(19)

3-dim Hydro for RHIC Energies

Au+Au E

_CM

=65 GeV/nucl. b=0.1 b

_max

A

_σ

=0.08 => σ~10 GeV/fm

n / n

₀

[ 1 ]

e [ GeV / fm

³

]

T= 0.0 fm/c n

_max

= 8.67 e

_max

=32.46 GeV / fm

³

L

_x,y

= 1.45 fm L

_z

=0.145 fm

. .

Bergen - Frankfurt - Los Alamos Nat. Lab. [Strottman]

(20)

3-dim Hydro for RHIC Energies

Au+Au E

_CM

=65 GeV/nucl. b=0.1 b

_max

A

_σ

=0.08 => σ~10 GeV/fm

n / n

₀

[ 1 ]

e [ GeV / fm

³

]

T=1.9 fm/c n

_max

= 8.66 e

_max

= 31.82 GeV / fm

³

L

_x,y

= 1.45 fm L

_z

=0.145 fm

. .

(21)

3-dim Hydro for RHIC Energies

Au+Au E

_CM

=65 GeV/nucl. b=0.1 b

_max

A

_σ

=0.08 => σ~10 GeV/fm

n / n

₀

[ 1 ]

e [ GeV / fm

³

]

T= 3.8 fm/c n

_max

= 7.77 e

_max

= 27.22 GeV / fm

³

L

_x,y

= 1.45 fm L

_z

=0.145 fm .

.

(22)

3-dim Hydro for RHIC Energies

Au+Au E

_CM

=65 GeV/nucl. b=0.1 b

_max

A

_σ

=0.08 => σ~10 GeV/fm

n / n

₀

[ 1 ]

e [ GeV / fm

³

]

T= 5.7 fm/c n

_max

= 6.36 e

_max

= 26.31 GeV / fm

³

L

_x,y

= 1.45 fm L

_z

=0.145 fm

. .

(23)

3-dim Hydro for RHIC Energies

Au+Au E

_CM

=65 GeV/nucl. b=0.1 b

_max

A

_σ

=0.08 => σ~10 GeV/fm

n / n

₀

[ 1 ]

e [ GeV / fm

³

]

T= 7.6 fm/c n

_max

= 5.22 e

_max

= 37.16 GeV / fm

³

L

_x,y

= 1.45 fm L

_z

=0.145 fm

. .

(24)

3-dim Hydro for RHIC Energies

Au+Au E

_CM

=65 GeV/nucl. b=0.1 b

_max

A

_σ

=0.08 => σ~10 GeV/fm

n / n

₀

[ 1 ]

e [ GeV / fm

³

]

T= 9.5 fm/c n

_max

= 4.45 e

_max

= 32.86 GeV / fm

³

L

_x,y

= 1.45 fm L

_z

=0.145 fm

. .

(25)

SUBATOMIC REACTIONS and ASTROPHYSICS

Project no. 6 "Internal structure of neutron stars", from Opava, Czech Republic, led by Z. Stuchlik, studied of the consequences of QGP formation, as these

influence the properties of neutron stars. Going beyond usual simplified

studies, the group numerically evaluated the properties of rotating compact

stars and their dependence on the QGP Equation of State. This way the

astrophysical and accelerator studies can mutually complement each other

for better understanding the high energy properties of matter.

(26)

Project 42:

Opava, CZ

[Z.Stuchlik,S.Hledik]

Cooling of extremely

compact stars

(27)

Radiation damage in solids [ Project: BCPL09]

TU Darmstadt

[Balogh]

(28)

Project: BCPL09

(29)

ATOMIC and STATISTICAL PHYSICS

Project no. 33: Dynamics of Atomic Bose Condensates, from Belfast, UK, led by J.F. McCann, aims for developing efficient high performance codes for modeling ultra-cold weakly interacting gases in a time dependent Hartree equation approach. The group developed a new method for studying complex deformation dynamics in irregular trap shapes such as asymmetric ellipsoids.

The findings of this work are prepared for publication in Phys. Review.

In addition the group has identified gyroscopic precession of vortex lines in ellipsoidal traps. This finding is now under further study.

(30)

Dynamics of Atomic Bose Condensates,BCPL33

[McCann, Belfast]

(31)

Belfast, BCPL33

(32)

Project no. 13 "Fluctuations and correlations - - analogies between fully

developed turbulence and relativistic heavy-ion collisions" from TU Dresden, Germany lead by M. Greiner, studied FD similarities on different scales,

reaching interesting and novel observations. The results are published in Phys.

Lett. A and cond-mat/0106347.

Project no. 27: Cusp Conditions on Wavelet Density Operators, from TU of Budapest, Hungary, led by J. Pipek, has found superstructures in electron density distributions in some mesoscopic systems. These superstructures are detected by calculating the structural entropy of distributions. These

calculations have required massive computing power.

(33)

Example for scaling

s ₀ (r)  s ₁ (r)

Daubechies 4 scaling function

J. Pipek pr. # 27:

(34)

Iterative inclusion of details into rough scaled distributions

The animation is the reverse play of the following process:

A two-dimensional wavelet transform (Mallat's pyramidal

algorithm, 1989) is applied to the source image, resulting in a

smoothed version of the original that can be further attacked in exactly the same way as the original image to obtain the wavelet coefficients at the

second highest resolution level.

This process repeats over and

(35)

Hordaland District Government – Regional Industry

• Rolls-Royce Marine AS – Ulstein Turbine AS

• Rogaland Forskning – Bergen: Multi-phase flow

• Naturgas West / Shell Int. – Techn. Applications: Fuel cells

• Rubitech AS: Dypvannsteknologi, Hydroacoustics

• Naturgas West: Fluidized beds

(36)

Hordaland Fylkeskommune - UiBed

Future Plans:

(37)

The Eurodyn’s Norwegian Heritage...

• Mid 1960s, Kongsberg’s KG2

(38)

The Eurodyn’s Norwegian Heritage...

• 21st Century, Eurodyn’s dual entry compressor and radial turbine

• Unique

– all radial technology gas generator

– not an aero-derivative

(39)

CONCLUSIONS

The selected 10 projects give a taste of the variability of BCPL's research projects. The other 3/4th of the projects not mentioned in this brief report are similarly important and effective. Most projects were fully successful.

In case of two projects we have registered problems due to the not perfect human

assistance. (The researchers of these projects were in Bergen during the Easter period last year where most of the BCPL and Parallab staff is on vacation.) To, improve human

assistance particularly in Computational Physics, Parallab employed in full position from Jan. 1, 2001 a Ph.D. Physicist, who also strengthens BCPL staff of project hosts. We hope this will further improve our service to EU users.

This also indicates that human help and assistance is crucial in the operation of a Research Infrastructure, so BCPL puts maximum emphasis on this.

The offered computational power to EU-RI user groups is increased by a factor of 12.8 from the beginning of 2002 by installing an IBM e-server supercomputer (For the same access charge to the EU as before!). At the time of the installation it was the 6th most powerful supercomputer in Europe.

Bergen, Feb. 28, 2002