<div dir="ltr">Hi Wang,<div><br></div><div>This revealed a small bug in the HYPRE solver for this problem. I'm attaching a patch here. </div><div><br></div><div>You can copy it to the "FLASH4.7" directory and apply it with</div><div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">path -p1 < MagDiff.patch</blockquote><div><br></div><div>Thanks for reporting the issue! </div><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr" style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px">*************************************************************************</div><div dir="ltr" style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px">Adam Reyes<div><img src="https://flash.rochester.edu/site/common/images/email_signature/FLASH.jpg"><br></div><div>Code Group Leader, Flash Center for Computational Science <span> </span><br></div><div>Research Scientist, Dept. of Physics and Astronomy</div><div>University of Rochester</div><div>River Campus: Bausch and Lomb Hall, 369 <span> </span><br></div><div>500 Wilson Blvd. PO Box 270171, Rochester, NY 14627</div><div dir="ltr">Email<span> </span><a href="mailto:adam.reyes@rochester.edu" target="_blank">adam.reyes@rochester.edu</a></div><div>Web<span> </span><a href="https://flash.rochester.edu/" target="_blank">https://flash.rochester.edu</a><br></div><div> (he / him / his)</div><div><img src="https://flash.rochester.edu/site/common/images/email_signature/FLASH-pride-sml.png"><br></div>*************************************************************************</div></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Apr 28, 2023 at 6:58 AM 王东升 <<a href="mailto:wangdswork@stu.xjtu.edu.cn">wangdswork@stu.xjtu.edu.cn</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">Hello Eddie,<br>
Thank you very much for your reply.<br>
I have checked the drawing of species variables (fill, line, vacu) and they were all still present in the simulation process. The resistivity of the vacuum has also been checked which remained 1e9(cm2/s) in the vacuum region.<br>
In fact, I didn’t make a lot of changes with Z-Pinch file. Here are some of the changes I've made:<br>
1. I set the fluid velocity to zero in the "Simulation_adjustEvolution.F90" file to rule out the magnetic freezing effect.<br>
2. The current was set to 1MA in the "current.bat" file.<br>
3. The thermal conductivity and radiation transport module was set to false in the "flash.par" file.<br>
4. The density of liner was set to 3e-3(g/cm3) as the plasma region, and the vacuum resistivity was set to 1e9(cm2/s).<br>
5. I set 0 to 0.8cm for species "line" as the plasma region and 0.8 to 1cm for species "vacu" as the vacuum region in the "Simulation_initBlock.F90" file.<br>
I attached the image of the species variable and the resistivity, and I also attached the "flash.par".<br>
Thank you again and look forward to your reply.<br>
Best wishes,<br>
Wang<br>
<br>
<br>
在 2023-04-28 03:04:33，"Hansen, Eddie" <<a href="mailto:ehansen@pas.rochester.edu" target="_blank">ehansen@pas.rochester.edu</a>> 写道：<br>
<br>
>Hi Wang,<br>
><br>
> <br>
><br>
>Have you plotted the species variables (fill, line, vacu) to make sure you still have vacuum in the domain? It would also be good to check the resistivity (res2).<br>
><br>
> <br>
><br>
>It’s hard to tell what might be going wrong without more info. Can you send the par file?<br>
><br>
> <br>
><br>
>--<br>
><br>
>Eddie Hansen<br>
><br>
>Research Scientist<br>
><br>
>Flash Center for Computational Science<br>
><br>
> <br>
><br>
> <br>
><br>
>From: flash-users <<a href="mailto:flash-users-bounces@flash.rochester.edu" target="_blank">flash-users-bounces@flash.rochester.edu</a>> on behalf of 王东升 <<a href="mailto:wangdswork@stu.xjtu.edu.cn" target="_blank">wangdswork@stu.xjtu.edu.cn</a>><br>
>Date: Thursday, April 27, 2023 at 1:19 PM<br>
>To: flash-users <<a href="mailto:flash-users@flash.rochester.edu" target="_blank">flash-users@flash.rochester.edu</a>><br>
>Subject: [FLASH-USERS] Diffusion problem of angular magnetic field in one dimension<br>
><br>
>Dear FLASH user & developers,<br>
>    I have a question about implicit magnetic diffusion in FLASH4.7. I wanted to observe the diffusion of angular magnetic field in one dimension, so I set the fluid velocity to zero in the "Simulation_adjustEvolution.F90" file to rule out the magnetic freezing effect. However, I found that the magnetic field diffused slowly in vacuum, which is quite different from the result of a paper I saw. I have attached the comparison diagram of the two results.<br>
>    The vacuum resistivity I set is 1e9(cm2/s), and the vacuum region is between 0.8 and 1cm. I modified it on the Z-Pinch file, and the current was set to 1MA.<br>
>    I would appreciate any help.<br>
>    Thanks!<br>
>Best,<br>
>Wang Dongsheng<br>
_______________________________________________<br>
flash-users mailing list<br>
<a href="mailto:flash-users@flash.rochester.edu" target="_blank">flash-users@flash.rochester.edu</a><br>
<br>
For list info, including unsubscribe:<br>
<a href="https://flash.rochester.edu/mailman/listinfo/flash-users" rel="noreferrer" target="_blank">https://flash.rochester.edu/mailman/listinfo/flash-users</a><br>
</blockquote></div>