[FLASH-USERS] Adding new source terms to FLASH

Robert Fisher rfisher1 at umassd.edu
Sat Mar 25 14:47:08 EDT 2017


Hi Simon :

   Adding your own source terms to FLASH is straightforward due to the
clean organization of the code. You can read about the process in very
broad strokes in the manual:

http://flash.uchicago.edu/site/flashcode/user_support/flash4_ug_4p4/node204.html

which includes a stub example module and instructions on how to introduce
your new source term to the Driver module. However, I'd recommend using one
of the actual source modules (Heat, Ionize, etc.) close to the physics of
your intended source term, and working backward from there.

  The process above will allow you access to integral timestep state data,
and should suffice for a first-order accurate source term. Accessing the
half-step data for higher-order accuracy is more subtle, and requires
adding the source terms directly to the hydro solver. For an example, you
can look at

source/physics/Hydro/HydroMain/split/PPM/PPMKernel/hydro_1d.F90

and the routines calling it to see how the gravitational accelerations are
applied to the momentum and energy equations at the half-timestep intervals.

  Best wishes,

  Bob

On Sat, Mar 25, 2017 at 1:59 PM, Simon Daley-Yates <
sdaley at star.sr.bham.ac.uk> wrote:

> Dear all
>
> I'm new to flash and just starting out.
>
> My question is: how easy is it to add new source terms to FLASH? Is this
> someting trivial or does it require intimate knolage of the way FLASH is
> written? Is there anyone that has done this and has any advice on how to go
> about achiving this?
>
> The calculation of the new source term requires access to the hole
> solution array every timestep/half timestep. Is this possible?
>
> Many thanks in advance,
> Simon
>



-- 
Dr. Robert Fisher
Associate Professor / Graduate Program Director
University of Massachusetts/Dartmouth
Department of Physics
285 Old Westport Road
North Dartmouth, Massachusetts 02747
robert.fisher at umassd.edu
http://www.novastella.org
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20170325/b46c953b/attachment.htm>


More information about the flash-users mailing list