<div dir="ltr"><div><span style="font-family:arial,sans-serif;font-size:12.800000190734863px">There are numerous reasons why the current implementation of the multigroup radiation solver is unfit for use in astrophysical problems. The radiation solver was designed for HEDP experiments, not astrophysics. It currently requires computational timesteps that are prohibitively short for astrophysical problems and does not correctly calculate the pressure gradient term in the equations of radiation hydrodynamics.</span><br>
</div><div dir="ltr"><div dir="ltr" style="font-family:arial,sans-serif;font-size:12.800000190734863px"><div><br></div><div>I am currently working on fixing both of these issues, and anticipate completing this in the next few months.</div>
<div><br></div><div>As for your segmentation fault, I didn't look at your output, but I do recall there was a segfault that would occur in the computeDt module in the vanilla solver. That may be where your issue lies. Regardless, I wouldn't bother with the RHD solver for astrophysics problems just yet.</div>
<div><br></div><div>Chris Lindner</div><div><a href="mailto:christopherclindner@gmail.com" target="_blank">christopherclindner@gmail.com</a></div><div>PhD Candidate</div><div>Department of Astronomy</div><div>The University of Texas at Austin</div>
</div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jun 25, 2014 at 12:00 PM, <span dir="ltr"><<a href="mailto:flash-users-request@flash.uchicago.edu" target="_blank">flash-users-request@flash.uchicago.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send flash-users mailing list submissions to<br>
<a href="mailto:flash-users@flash.uchicago.edu">flash-users@flash.uchicago.edu</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://flash.uchicago.edu/mailman/listinfo/flash-users" target="_blank">http://flash.uchicago.edu/mailman/listinfo/flash-users</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:flash-users-request@flash.uchicago.edu">flash-users-request@flash.uchicago.edu</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:flash-users-owner@flash.uchicago.edu">flash-users-owner@flash.uchicago.edu</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of flash-users digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Help with multigroup radiation for CCSN (Joseph Smidt)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Tue, 24 Jun 2014 13:34:41 -0600<br>
From: Joseph Smidt <<a href="mailto:josephsmidt@gmail.com">josephsmidt@gmail.com</a>><br>
Subject: [FLASH-USERS] Help with multigroup radiation for CCSN<br>
To: <a href="mailto:flash-users@flash.uchicago.edu">flash-users@flash.uchicago.edu</a><br>
Message-ID:<br>
<CACfQXCqjdsuykGgjXPGzTq_JBj7UgmSLxv0xf=<a href="mailto:BApR2AkVpsww@mail.gmail.com">BApR2AkVpsww@mail.gmail.com</a>><br>
Content-Type: text/plain; charset=UTF-8<br>
<br>
Hey everyone,<br>
<br>
I downloaded the new FLASH4.2.1 which contains a core collapse<br>
supernova (CCSN) test problem by Sean Couch. I wanted to take this<br>
problem and turn on radiation. (As well as supply my own initial<br>
profile)<br>
<br>
Unfortunately I am getting segfaults and I am new enough to the<br>
code that none of my debugging has helped. Let me now walk through<br>
the process for what I did.<br>
<br>
I took the vanilla CCSN test problem that comes with flash,<br>
renamed it A140 and added the recommended initial conditions lines<br>
shown here [1] to the Simulation_initBlack.F90 as follows: [2]. (See<br>
lines 46-51 for the exact changes) I also changed the Config to this<br>
[3].<br>
<br>
I then setup the code using this command:<br>
<br>
./setup A140 -auto -1d +spherical -objdir a140 -nxb=16<br>
threadBlockList=False +pm4dev threadWithinBlock=False +newMpole +mgd<br>
mgd_meshgroups=1<br>
<br>
I then compiled and ran with this flash.par [4]. You will note I<br>
am just trying to run with one group with constant opacities. The CCSN<br>
test requires an initial profile file which is here [5]. When the<br>
code runs I get this error [6] and gdb gives this backtrace [7].<br>
<br>
Sorry for all the pastebins, but I hope they detail all my steps.<br>
I am just trying to turn on mgd radiation for the CCSN test problem<br>
provided by Flash and I am stuck so any help is appreciated.<br>
<br>
<br>
<br>
[1] <a href="http://flash.uchicago.edu/site/flashcode/user_support/flash4_ug_4p22/node29.html" target="_blank">http://flash.uchicago.edu/site/flashcode/user_support/flash4_ug_4p22/node29.html</a><br>
<br>
[2] <a href="http://pastebin.com/vS8F67Bz" target="_blank">http://pastebin.com/vS8F67Bz</a><br>
<br>
[3] <a href="http://pastebin.com/rGwmZMBq" target="_blank">http://pastebin.com/rGwmZMBq</a><br>
<br>
[4] <a href="http://pastebin.com/NY2mAzx7" target="_blank">http://pastebin.com/NY2mAzx7</a><br>
<br>
[5] <a href="http://pastebin.com/DKaUSZLU" target="_blank">http://pastebin.com/DKaUSZLU</a><br>
<br>
[6] <a href="http://pastebin.com/e8TXDUwT" target="_blank">http://pastebin.com/e8TXDUwT</a><br>
<br>
[7] <a href="http://pastebin.com/FdYEPjsB" target="_blank">http://pastebin.com/FdYEPjsB</a><br>
<br>
<br>
--<br>
------------------------------------------------------------------------<br>
Joseph Smidt <<a href="mailto:josephsmidt@gmail.com">josephsmidt@gmail.com</a>><br>
<br>
Theoretical Division<br>
P.O. Box 1663, Mail Stop B283<br>
Los Alamos, NM 87545<br>
Office: <a href="tel:505-665-9752" value="+15056659752">505-665-9752</a><br>
Fax: <a href="tel:505-667-1931" value="+15056671931">505-667-1931</a><br>
<br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
flash-users mailing list<br>
<a href="mailto:flash-users@flash.uchicago.edu">flash-users@flash.uchicago.edu</a><br>
<a href="http://flash.uchicago.edu/mailman/listinfo/flash-users" target="_blank">http://flash.uchicago.edu/mailman/listinfo/flash-users</a><br>
<br>
<br>
End of flash-users Digest, Vol 80, Issue 10<br>
*******************************************<br>
</blockquote></div><br></div></div></div>