<html><body><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; color: #000000"><div>Hello once again,</div><div><br data-mce-bogus="1"></div><div>To correct the previous mail, the setup was actually correct for running a perfect gas.<br></div><div>In a multiTemp model, the perfect gas equation I wrote down was innacurate. <!--StartFragment--><div><img src="cid:e52d6348a03126128a666ebb86f56d44e3ce4d12@zimbra" data-mce-src="cid:e52d6348a03126128a666ebb86f56d44e3ce4d12@zimbra" width="114" height="51"></div><div>The equation should take into account the contribution of both populations:<br data-mce-bogus="1"></div><!--EndFragment--> <br></div><div><img src="cid:9c785bb92e488c198569cf957f256384faafc526@zimbra" data-mce-src="cid:9c785bb92e488c198569cf957f256384faafc526@zimbra" width="206" height="53">(assuming that the effect of Trad is negligible)<br></div><div><br data-mce-bogus="1"></div><div>With this corrected formula, the output pressure and the calculated one match perfectly.<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Still, I would like to run my case with a 1-Temperature model.<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Sorry for the previous mail, no problem with the .par inputs for a perfect gas.</div><div><br data-mce-bogus="1"></div><div>Thibault G<br data-mce-bogus="1"></div><div><br></div><hr id="zwchr" data-marker="__DIVIDER__"><div data-marker="__HEADERS__"><b>De: </b>"Thibault Goudal" <thibault.goudal@u-bordeaux.fr><br><b>À: </b>"flash-users" <flash-users@flash.uchicago.edu><br><b>Envoyé: </b>Lundi 14 Décembre 2020 23:14:49<br><b>Objet: </b>Simple 1D LaserSlab</div><div><br></div><div data-marker="__QUOTED_TEXT__"><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; color: #000000;" data-mce-style="font-family: arial, helvetica, sans-serif; font-size: 12pt; color: #000000;"><div>Hi everyone,</div><div>I would like to benchmark with FLASH an analytical code I've just developed with some hypotheses. Here is the plot:</div><div>I would like to run a 1D cartesian LaserSlab example with a monoatomic perfect gas without any radiative effect and a 1-Temperature model (Ti=Te, very simple)</div><br><div>I tried to keep the "basic" laserSlab example:</div><br><div>I. <span style="text-decoration: underline;" data-mce-style="text-decoration: underline;">setup line</span>:</div><br><div><strong>./setup -auto LaserSlab -1d -nxb=16 +hdf5typeio species=cham,targ +laser +uhd3t +mtmmmt</strong></div><div>(I would like to keep two species so I think that +mtmmmt is essential even in 1T, +uhd3t is maybe not what I want to use in the 1Temp case)</div><br><br><br><div>II. <span style="text-decoration: underline;" data-mce-style="text-decoration: underline;">Specific parameters in .par</span>:</div><br><div>rt_useMGD = .false.<br>useOpacity = .false.</div><div>diff_useRadDiffusion = .false.<strong><br></strong></div><div><strong></strong><br></div><div>#EOS section<br>ms_targgamma = 1.666667<br>eos_targEosType = "eos_gam"<strong><br><br></strong></div><div><strong></strong><br></div><div>I compare the output pressure from the code and the one given by the perfect gas relation formed with the temperature:</div><br><div><img src="cid:e52d6348a03126128a666ebb86f56d44e3ce4d12@zimbra" data-mce-src="cid:e52d6348a03126128a666ebb86f56d44e3ce4d12@zimbra" width="114" height="51"></div><br><div>I cannot identify the two pressures.</div><div>I concluded that what I set as a perfect gas in the flash.par was maybe not that perfect...</div><br><div>Do you have any idea of such a behaviour?</div><br><div>Thanks</div><div><strong></strong><br></div><br><br><br><div>--</div><div><div style="text-align: left;" data-mce-style="text-align: left;"><span style="color: #000000;" data-mce-style="color: #000000;">Thibault Goudal</span></div><div style="text-align: left;" data-mce-style="text-align: left;"><span style="color: #000000;" data-mce-style="color: #000000;">PhD Student CEA</span></div><div style="text-align: left;" data-mce-style="text-align: left;"><span style="color: #000000;" data-mce-style="color: #000000;">CELIA CNRS UMR5107 laboratory</span></div><div style="text-align: left;" data-mce-style="text-align: left;"><span style="color: #000000;" data-mce-style="color: #000000;">+33664160842</span></div><div style="text-align: left;" data-mce-style="text-align: left;"><span style="color: #000000;" data-mce-style="color: #000000;">thibault.goudal@u-bordeaux.fr</span></div></div></div></div></div></body></html>