<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
"<span style="display: inline !important; background-color: rgb(255, 255, 255);" class="ContentPasted0">So, to check densities in our simulation (we removed the liner btw) we compare with these papers by
</span><span style="display: inline !important; color: rgb(237, 92, 87); background-color: rgb(255, 255, 255);" class="ContentPasted0">multiplying</span><span style="display: inline !important; background-color: rgb(255, 255, 255);" class="ContentPasted0">
the FLASH density in the Z-pinch example with DD mass in grams (for the fill)"</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
I said multiplying but I meant dividing, sorry for the mistake.</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Alvaro Sanchez<br>
<b>Sent:</b> Wednesday, April 12, 2023 5:51 PM<br>
<b>To:</b> flash-users@flash.rochester.edu <flash-users@flash.rochester.edu><br>
<b>Subject:</b> Help with Z-pinch example.</font>
<div> </div>
</div>
<style type="text/css" style="display:none">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<div class="x_elementToProof" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
<div class="x_Ar x_Au x_Ao" style="display:block">
<div class="x_Am x_Al x_editable x_LW-avf x_tS-tW x_tS-tY x_ContentPasted0" tabindex="1" style="min-height:298px">
Hello dear folks from FLASH community, a pleasure to be here.
<div><br class="x_ContentPasted0">
</div>
<div class="x_ContentPasted0">I'm new in FLASH so I still probably have beginner mistakes and questions, but I really need help with a couple of things, so thank you very much in advance :)</div>
<div><br class="x_ContentPasted0">
</div>
<div class="x_ContentPasted0">We are working on a Z-pinch, using the example as a template, so we are trying to understand pretty much everything (or as much as possible) it has.</div>
<div><br class="x_ContentPasted0">
</div>
<div class="x_ContentPasted0">First is a very basic thing, density units. We suppose g/cm^3. The problem we have with the example is the high density of vacuum, even the density of the fill seems too high for us, as well as the density during the stagnation
period. We've been researching and, in almost every seen paper, number density is used.</div>
<div class="x_ContentPasted0">So, to check densities in our simulation (we removed the liner btw) we compare with these papers by multiplying the FLASH density in the Z-pinch example with DD mass in grams (for the fill). And we get much lower values than those
used in BC and simulated in the example (which could make sense, since our Z-pinch is a different one). We are also using lower values as initial conditions because we need, during our stagnation period, values of 1e17 cm^-3 for dens.</div>
<div><br class="x_ContentPasted0">
</div>
<div class="x_ContentPasted0">Are we making a mistake calculating the number density?</div>
<div class="x_ContentPasted0">We've seen, in ../HydroMain/unsplit/MHD_StaggeredMesh/<i><b class="x_ContentPasted0">hy_uhd_getFaceCurrents.F90</b></i></div>
<div><i><b><br class="x_ContentPasted0">
</b></i></div>
<div><i><b><br class="x_ContentPasted0">
</b></i></div>
<div><span style="background-color:rgb(0,0,255)"><font class="x_ContentPasted0" style="color:rgb(0,255,255)">!! Compute electron number density at (i,j,k)</font></span></div>
<div><span style="background-color:rgb(0,0,255)"><font class="x_ContentPasted0" style="color:rgb(0,255,255)">nele_rght = ye * hy_avogrado * U(DENS_VAR,i,j,k)</font></span></div>
<div><br class="x_ContentPasted0">
</div>
<div class="x_ContentPasted0">SO, we are worried about being wrong since the beginning.</div>
<div><br class="x_ContentPasted0">
</div>
<div><br class="x_ContentPasted0">
</div>
<div class="x_ContentPasted0">We found that .F90 in our search of the interpolation of the current from the Circuit source term. We don't understand how FLASH spatially interpolates the input current (we use the FileInput way) in a 2D cylindrical axisymmetric plane.</div>
<div class="x_ContentPasted0">Does someone know where in the internal code is this issue coded? or...</div>
<div class="x_ContentPasted0">Does someone know how it works? Does it goes through the higher density the higher the current density? Or am I completely wrong?</div>
<div><br class="x_ContentPasted0">
</div>
<div class="x_ContentPasted0">Thank you very very much if you read it, I really appreciate it, so thank you very much.</div>
<div><br class="x_ContentPasted0">
</div>
<span class="x_ContentPasted0">Any help is always nice, have a good day.</span></div>
</div>
<br>
</div>
</div>
</body>
</html>