Re: [fluka-discuss]: low density question

From: Francesco Cerutti <Francesco.Cerutti_at_cern.ch>
Date: Sun, 24 Nov 2013 21:14:12 +0100

Dear Alec,

your diagnosis is pretty correct, the crash is due to the terribly low
density value, as used in the Sternheimer parameterization of the density
effect correction to stopping power. But please note that values below
10^(-10) g/cm^3, far from working well, are actually meaningless, since
the code in this case just resets the respective material to vacuum
(issuing a "material changed" message - by the way incomplete - in the
output file).

Kind regards

Francesco

**************************************************
Francesco Cerutti
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
tel. ++41 22 7678962
fax ++41 22 7668854

On Fri, 22 Nov 2013, Alexander MacKinnon wrote:

> Dear colleagues,
>
> I've been experimenting with FLUKA as a tool for particle transport in an astrophysical setting. I've followed to some extent the example of the cards for studying cosmic rays in the terrestrial atmosphere but have run up against an interesting problem. Our lengthscales are very large and our densities very small, by most people's standards, and I think it's from density that the problem arises.
>
> The attached .inp file attempts to implement a target with successive layers of monotonically increasing density, according to one of the models for solar atmospheric structure. With this .inp file FLUKA fails without producing any .out, .log or .err files. The temporary versions of the .out and .log files do exist, however, in the fluka_XXXX folder that doesn't get deleted when there's a crash. The .log file from there tells us there's a floating point exception in a function called dpdx.f which I guess calculates the energy loss rate.
>
> By experiment we discovered that this only happens if one of the layers has a density in a well-defined range, between 10^(-10) and 1.82 10^(-10) gm.cm^(-3) . As long as all densities are greater or less than this value FLUKA runs without crashing.
>
> There are some ways round this problem. Only the total column density in each of the layers really matters for most purposes, so we could easily generate a different density structure that would tell us what we need to know, avoiding the problematic density range, and then re-apply the results in the real atmospheric structure. For some purposes we only need thick target yields or secondaries emergent at the top of the atmosphere and again the detailed density structure would not be important. Etc. But we'd prefer to understand what's happening in case this highlights a bigger issue we should be aware of.
>
> Thanks in advance for advice.
>
> Regards,
>
> Alec MacKinnon
>
> PS let me anticipate some more general comments. The higher layers of our atmosphere are of course not neutral but there's very little column depth there. The solar atmosphere is almost neutral throughout a large fraction of the stopping depth of the particles of interest. Results from a code like FLUKA give a valuable reference case even if there are significant plasma physics complications in transport.
>
>
> Dr Alexander L MacKinnon
> Centre for Open Studies
> University of Glasgow
> 11 Eldon Street
> GLASGOW G3 6NH
> UK
> (+44)(0)141 330 1857
> http://www.gla.ac.uk/services/centreforopenstudies/staff/alexandermackinnon/
Received on Sun Nov 24 2013 - 22:14:26 CET

This archive was generated by hypermail 2.3.0 : Sun Nov 24 2013 - 22:14:26 CET