<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
<title></title>
</head>
<body bgcolor="#ffffff" text="#000000">
In general, source term contributions should be mild, by some
standards. In your case, try limiting kinetic energy changes to perhaps
5-10% per time step.<br>
<br>
Also, it is best to call equation of state after all/any (depending on
what is being modified) source terms are applied to make sure the hydro
state is consistent.<br>
<br>
Driver_sourceTerms looks like the right place to modify hydro state
locally.<br>
<br>
Tomek<br>
--<br>
James Guillochon wrote:
<blockquote cite="mid:48D0FF0F.2040002@ucolick.org" type="cite">
<pre wrap="">OK, I will give that a try. Is there anything else I should be cautious
about? Is the Driver_sourceTerms routine the appropriate place to do
this sort of thing?
Thanks for your help!
Tomasz Plewa wrote:
</pre>
<blockquote type="cite">
<pre wrap="">James -
At the very least you need to update the total gas energy.
Tomek
--
James Guillochon wrote:
</pre>
<blockquote type="cite">
<pre wrap="">Hi everyone,
I am altering the velocity of a couple of grid squares in my
simulation in the Driver_sourceTerms function at each time step, and
after a couple hundred steps I end up with a grid cell or two with a
ridiculously high pressure/density, which causes the Helmholtz EOS
routine to not converge. I was wondering if there are any other
variables I should be altering when I give these grid cells their new
velocities? Or if there is a better place in Flash to be adjusting
the velocities than the Driver_sourceTerms routine?
Thanks!
- James
</pre>
</blockquote>
<pre wrap="">
!DSPAM:10135,48d0fd1a20044307041292!
</pre>
</blockquote>
<pre wrap=""><!---->
</pre>
</blockquote>
</body>
</html>