<html>#2549: inlcude FLRWSolver in ET
<table style='border-spacing: 1ex 0pt; '>
<tr><td style='text-align:right'> Reporter:</td><td>Roland Haas</td></tr>
<tr><td style='text-align:right'> Status:</td><td>new</td></tr>
<tr><td style='text-align:right'>Milestone:</td><td>ET_2022_05</td></tr>
<tr><td style='text-align:right'> Version:</td><td>development version</td></tr>
<tr><td style='text-align:right'> Type:</td><td>enhancement</td></tr>
<tr><td style='text-align:right'> Priority:</td><td>major</td></tr>
<tr><td style='text-align:right'>Component:</td><td></td></tr>
</table>
<p>Comment (by Chi Tian):</p>
<p>Hi all, I have reviewed the documents and majority of the the codes. At least at physics level, everything looks fine to me. Only some tiny issues regarding the way FLRWSolver dealing with the extrinsic curvature:</p>
<ol>
<li>It seems there is a typo in Eq. (10) in the document, where minus sign is missing (at least for the case where the author mentions “Setting f =1/3 and n=2 in the above yields ∂tα=∂ta for EdS“).</li>
<li>The name of the variable <code>kvalue</code> in several <code>flrw_*.f90</code> files seems a little confusing. From its name, it seems this variable represents the extrinsic curvature K. However, its true value seems to be <code>3 * kvalue</code> if my understanding was correct.</li>
<li>By following the gauge condition <code>\dot{\alpha} = f * K</code>, effectively it is setting the coordinate time to something closed to the conformal time <code>\eta</code>, and the value of <code>\alpha</code> will grow roughly following the growth of the scale factor <code>a</code>. Could this possibly introduce numerical instabilities to <code>GRhydro</code> when starting from a very high redshift, and universe keeps expanding for quit long? Since I would guess most of astrophysical GR simulations may choose gauge conditions that make <code>\alpha</code> close to 1. </li>
</ol>
<p></p>
<p>--<br/>
Ticket URL: <a href='https://bitbucket.org/einsteintoolkit/tickets/issues/2549/inlcude-flrwsolver-in-et'>https://bitbucket.org/einsteintoolkit/tickets/issues/2549/inlcude-flrwsolver-in-et</a></p>
</html>