[FLASH-USERS] Problem of 2D cylindrical setup

[pchang]

Hi all,

Thank you for all the suggestions. I didn't have magnetic field in the 
simulations. Also, I did rotate the Windtunnel problem and have the 
lower y boundary as the axis of rotation. After some tests, I solved the 
problem by removing the line "REQUIRES 
Grid/GridBoundaryConditions/OneRow" in Config and it works now.

It seems like OneRow does not support axisymmetric.

Cheers,

Po-Yu


Lu, Yingchao 於 2022/5/24 下午 09:51 寫道:
>
> Do you have magnetic field in the simulations? If not with B-field, I 
> think "axisymmetric" and "reflect" do the same thing. If there is 
> B-field, as the user guide says "axisymmetric" is "like reflect, but 
> both normal and toroidal vector components change sign."
>
> Best,
>
> Yingchao
>
> *From: *flash-users <flash-users-bounces at flash.rochester.edu> on 
> behalf of Andy Sha Liao <andy at f.energy>
> *Date: *Tuesday, May 24, 2022 at 9:42 AM
> *To: *Po-Yu Chang <pchang at mail.ncku.edu.tw>
> *Cc: *flash-users at flash.rochester.edu <flash-users at flash.rochester.edu>
> *Subject: *Re: [FLASH-USERS] Problem of 2D cylindrical setup
>
> I think the cylindrical geometry only works if the lower y boundary is 
> the axis of rotation. So, you will have to manually rotate Windtunnel 
> to blow vertically (along y axis).
>
> Andy
>
> On Sun, May 22, 2022 at 1:37 PM Po-Yu Chang <pchang at mail.ncku.edu.tw> 
> wrote:
>
>     Hi all,
>
>     I am trying to simulate a hydrodynamic problem in 2D cylindrical
>     coordinate. Therefore, I started from the example "WindTunnel". This
>     example was originally in the cartesian coordinate. Therefore, I
>     had the
>     following for the setup:
>
>     ./setup -auto -2d +Cylindrical WindTunnel -objdir=WindTunnelCyl
>
>     I found that in flash.par, the geometry was still set as
>
>     "geometry = cartesian"
>
>     Therefore, I replaced that line with
>
>     "geometry = Cylindrical"
>
>     Then, I found that if I set the boundary:
>
>     "xl_boundary_type = axisymmetric"
>
>     It gave an error message: "Error message is unsupported boundary
>     condition on Lower Face"
>
>     If I replaced the boundary to
>
>     "xl_boundary_type = reflect"
>
>     Then, it ran without errors.
>
>     However, I believe the reflection boundary on the axis is not
>     correct.
>     It seems like the problem was still under cartesian coordinate even I
>     set it to cylindrical coordinate during the "./setup" and changed the
>     geometry to cylindrical in the flash.par. So, what did I miss for
>     setting the cylindrical coordinate?
>
>     Thank you.
>
>     Cheers,
>
>     Po-Yu
>
>     -- 
>     Po-Yu Chang
>     Assistant Professor
>     Institute of Space and Plasma Sciences (ISAPS)
>     National Cheng Kung University (NCKU)
>     No.1, University Road,
>     Tainan City 70101, Taiwan
>     (TEL) +886-6-2757575 ext. 65916
>     pchang at mail.ncku.edu.tw
>     http://capst.ncku.edu.tw/PGS/
>
>     _______________________________________________
>     flash-users mailing list
>     flash-users at flash.rochester.edu
>
>     For list info, including unsubscribe:
>     https://flash.rochester.edu/mailman/listinfo/flash-users
>
-- 
Po-Yu Chang
Assistant Professor
Institute of Space and Plasma Sciences (ISAPS)
National Cheng Kung University (NCKU)
No.1, University Road,
Tainan City 70101, Taiwan
(TEL) +886-6-2757575 ext. 65916
pchang at mail.ncku.edu.tw
http://capst.ncku.edu.tw/PGS/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20220524/13c26d78/attachment.htm>