<div dir="ltr">Hi Paul,<div><br></div><div>Thanks so much for the help and quick response. The code is running fine. I do realize that this is not a Womersley inlet. In fact I have access to patient specific Doppler measured velocity profiles (V of T in this case) that I am implementing with the assumption of a fully developed flow. This is justified to some extent since the vein in reality, before the spot that we consider the inlet for our study, is long enough to develop a parabolic profile in the tube. </div>
<div><br></div><div>Am still waiting for any ideas if the Sponge outlet condition can be used with stress formulation for non-Newtonian flows. </div><div><br></div><div>Cheers,</div><div>Javid </div></div><div class="gmail_extra">
<br><br><div class="gmail_quote">On Tue, Aug 6, 2013 at 12:00 PM, <span dir="ltr"><<a href="mailto:nek5000-users-request@lists.mcs.anl.gov" target="_blank">nek5000-users-request@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send Nek5000-users mailing list submissions to<br>
<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:nek5000-users-request@lists.mcs.anl.gov">nek5000-users-request@lists.mcs.anl.gov</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:nek5000-users-owner@lists.mcs.anl.gov">nek5000-users-owner@lists.mcs.anl.gov</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of Nek5000-users digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Pulsatile Flow in NEK5000 (<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>)<br>
2. Re: Pulsatile Flow in NEK5000 (<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Mon, 5 Aug 2013 16:19:37 -0500<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>, Cassel Kevin <<a href="mailto:cassel@iit.edu">cassel@iit.edu</a>><br>
Subject: [Nek5000-users] Pulsatile Flow in NEK5000<br>
Message-ID:<br>
<<a href="mailto:mailman.12039.1375737604.15097.nek5000-users@lists.mcs.anl.gov">mailman.12039.1375737604.15097.nek5000-users@lists.mcs.anl.gov</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Hi There NEKs,<br>
<br>
I am currently trying to run a case with pulsatile inlet condition. The<br>
context is pulsatile stenotic flows. I have two Question regrading the set<br>
up:<br>
<br>
1st Question) How do I access the "Time" variable in NEK5000? I need to<br>
define a time dependent inlet velocity condition in the user file. V as a<br>
function of Y and T. I already have the function which is a patient<br>
specific velocity inlet profile. Please note that this should be a<br>
pulsatile flow. This, I assume, is different than a periodic flow. I'll be<br>
investigation the effect of pulsatility on some 2D geometries.<br>
It would be wonderful to know how to implement time dependent inlet<br>
velocity condition in user file. Previously I have only used steady<br>
and parabolic inlet velocity.<br>
<br>
2nd Question) I will be using a simple Dirichlet BC on pressure for the<br>
outlet. I have also used the Nozzle outlet condition but, as it turns out,<br>
I cannot use the Nozzle condition with stress formulation. I wanted to know<br>
if the Sponge condition for outlet is proper and if it can be coupled with<br>
stress formulation? I am considering a non-Newtonian fluid with variable<br>
viscosity for my simulation.<br>
<br>
Any assistance is greatly appreciated.<br>
<br>
Best of luck to everyone,<br>
Javid<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20130805/8178cef4/attachment-0001.html" target="_blank">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20130805/8178cef4/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Mon, 5 Aug 2013 17:24:32 -0500 (CDT)<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
Subject: Re: [Nek5000-users] Pulsatile Flow in NEK5000<br>
Message-ID:<br>
<<a href="mailto:mailman.12049.1375741473.15097.nek5000-users@lists.mcs.anl.gov">mailman.12049.1375741473.15097.nek5000-users@lists.mcs.anl.gov</a>><br>
Content-Type: text/plain; charset=utf-8<br>
<br>
<br>
Hi Javid,<br>
<br>
You simply need to use the "time" variable, e.g.,<br>
<br>
<br>
ux=1-y*y<br>
<br>
gives the steady state solution, whereas<br>
<br>
eps = 0.1<br>
period = 5.0<br>
arg = 2*pi*time/period<br>
ux = (1-y*y)*(1+eps*sin(arg))<br>
<br>
would give a sinusoidal fluctuation. (Note that this is<br>
_not_ the Womersely solution that is the exact inlet<br>
condition for fully developed pulsatile flow...)<br>
<br>
Regarding outflow, I don't know why the standard nozzle<br>
approach would not work with the stress formulation.<br>
Am on travel at the moment and don't have time to dig<br>
into this, but perhaps someone else can assist.<br>
<br>
Best, Paul<br>
<br>
<br>
<br>
----- Original Message -----<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>, "Cassel Kevin" <<a href="mailto:cassel@iit.edu">cassel@iit.edu</a>><br>
Sent: Monday, August 5, 2013 4:19:37 PM<br>
Subject: [Nek5000-users] Pulsatile Flow in NEK5000<br>
<br>
<br>
<br>
Hi There NEKs,<br>
<br>
<br>
I am currently trying to run a case with pulsatile inlet condition. The context is pulsatile stenotic flows. I have two Question regrading the set up:<br>
<br>
<br>
1st Question) How do I access the "Time" variable in NEK5000? I need to define a time dependent inlet velocity condition in the user file. V as a function of Y and T. I already have the function which is a patient specific velocity inlet profile. Please note that this should be a pulsatile flow. This, I assume, is different than a periodic flow. I'll be investigation the effect of pulsatility on some 2D geometries.<br>
It would be wonderful to know how to implement time dependent inlet velocity condition in user file. Previously I have only used steady and parabolic inlet velocity.<br>
<br>
<br>
2nd Question) I will be using a simple Dirichlet BC on pressure for the outlet. I have also used the Nozzle outlet condition but, as it turns out, I cannot use the Nozzle condition with stress formulation. I wanted to know if the Sponge condition for outlet is proper and if it can be coupled with stress formulation? I am considering a non-Newtonian fluid with variable viscosity for my simulation.<br>
<br>
<br>
Any assistance is greatly appreciated.<br>
<br>
<br>
Best of luck to everyone,<br>
Javid<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
<br>
End of Nek5000-users Digest, Vol 54, Issue 1<br>
********************************************<br>
</blockquote></div><br></div>