[FLASH-USERS] BC: flux direct prescription

[Dongwook Lee]

Hi Massimo,

One easy way to do so is to call Grid_getCellCoords to get the location
of the boundary where you want to inject fluxes. Once you detect the
boundary  in the routine FLASH computes fluxes (e.g.,
hy_uhd_getFaceFluxes.F90 in the unsplit solver, and
rieman.F90 or hydro_1d.F90), you can easily modify fluxes for
injection.

Hope this helps,

Dongwook

----------------------------------------------------
Dongwook Lee, Ph.D., Research Scientist
ASC Center for Astrophysical Thermonuclear Flashes
University of Chicago
5747 S. Ellis, Rm 319
(773) 834-6830

On Fri, 7 Jan 2011, Massimo Gaspari wrote:

>
> Hi all,
> I am trying to set a particular, but simple, boundary condition:I want to directly set the fluxes at the boundary intercells of oneside of the domain.
> I want in fact to avoid to prescribe fictious data values at the BC ghost cells (the usual centered values), because first I have physically only the fluxes and, second, it is probably impossibleto give a centered value to guess the right fluxes (PPM goes throughprimitive interpolation, correction and using exact Riemann problem,in order to find the "right" fluxes at interfaces).
> I think there should be a subroutine or a "spot" to call my own fluxesat boundary interfaces somewhere... but I don't know exactly where.
> Can you please help me?
> Thank you very much!
>
> Max