[FLASH-USERS] Implementing Jet via
Ryan Farber
rjfarber at umich.edu
Tue Jul 23 16:56:54 EDT 2019
Hi Alex,
I agree that the WindTunnel's Grid_applyBCEdge is a good starting point.
In Grid_applyBCEdge, the user boundary condition is triggered if bcType ==
USER_DEFINED where USER_DEFINED is #defined to be -38 in my constants.h
Grid_init grabs xl_boundary_type from flash.par and maps that string to an
integer. So you could check to make sure that initial part is proceeding
correctly (i.e., if you print out gr_domainBC(LOW,:) do you get a value
matching USER_DEFINED in your constants.h?
If so, you can then forward trace gr_domainBC until you get to where your
abort is that Ernesto pointed out (or backtrace from there).
In any case, I would guess that using +cylindrical geometry is the issue.
Do you get the abort message without that (i.e., for cartesian)?
Also, what version of FLASH are you using?
Best,
--------
Ryan
On Tue, Jul 23, 2019 at 8:47 AM Ernesto Zurbriggen <ezurbriggen at unc.edu.ar>
wrote:
> Hi Alex!
>
> Did you have a look at the
> routine /Grid/GridBoundaryConditions/gr_bcApplyToOneFace.F90, specifically
> at line 180?
> It's probable you are setting a boundary condition type that is not being
> handled.
>
> Best!
>
>
> El mar., 23 jul. 2019 a las 11:00, Alex Robinson - UKRI STFC (<
> alex.robinson at stfc.ac.uk>) escribió:
>
>> Hello flash users,
>> I would have thought that this has already been done, however I
>> would like to set up a simulation where a jet is injected by implementing a
>> boundary condition for this. I had thought that using ‘Wind Tunnel’ as a
>> basis for this would be a good idea, and using the
>> ‘GridBoundaryConditions/OneRow’ approach where I have my own
>> Grid_applyBCEdge.F90 file in the simulation directory. This compiles fine
>> (setup line reads : ./setup Jet -auto -2d +cylindrical +ug -nxb=200
>> -nyb=200) . However when I try to run this, the simulation aborts with the
>> message : DRIVER_ABORT : unsupported boundary condition on Lower Face.
>> This is after I set the yl_boundary_type to be either ‘user’ or
>> ‘user-defined’. I have specified REQUIRES
>> Grid/GridBoundaryConditions/OneRow in the config file.
>>
>> This message implies that despite my Grid_applyBCEdge.F90 file being in
>> the simulation directory that it hasn’t overridden the default method as no
>> ‘user’ implementation appears in the standard file in the source. However
>> I can’t see what I haven’t done that has lead to the overriding not
>> happening. Any advice?
>>
>> Thanks,
>> Alex
>>
>
>
> --
> *PhD Ernesto Zurbriggen*
>
> *Email: ezurbriggen at unc.edu.ar <ezurbriggen at unc.edu.ar>*
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20190723/86211656/attachment-0001.htm>
More information about the flash-users
mailing list