ray tracing

16 views
Skip to first unread message

Hans Hornung

unread,
Feb 27, 2012, 3:57:52 PM2/27/12
to amrita-ebook
One of my graduate students is doing some interesting work on acoustic
ray tracing in hypersonic flows. In particular, the region near the
stagnation point of a blunt body is of interest. So I did some
computations using Amrita. What he needs is X, Y, U, V and speed of
sound in the coarse and refined grids. Since the command flowout
stores this information and much more, there must be a way of
extracting all of what he needs from one of the files .dstr or .solv.
We did it so far with just the coarse grid by using probe, but finer
resolution is needed near the shock for the ray tracing to behave
properly. Could anyone please help us with that?

Perhaps someone has already produced a script for acoustic ray
tracing.

Hans

James Quirk

unread,
Feb 27, 2012, 4:17:58 PM2/27/12
to Hans Hornung, amrita-ebook
Hans,

On Mon, 27 Feb 2012, Hans Hornung wrote:

> One of my graduate students is doing some interesting work on acoustic
> ray tracing in hypersonic flows. In particular, the region near the
> stagnation point of a blunt body is of interest. So I did some
> computations using Amrita. What he needs is X, Y, U, V and speed of
> sound in the coarse and refined grids. Since the command flowout
> stores this information and much more, there must be a way of
> extracting all of what he needs from one of the files .dstr or .solv.

The attached script, run_get_flowfield, shows how to grab
the information you're after. It does this by compiling and running
a small FORTRAN fragment to loop over all the cells in AMR_SOL's
current solution field. Thus there is no need to access the .dstr and
.solv files, directly. You will still need to connect the dots,
as it were, but the attached will at least get ypu started.

James

run_get_flowfield

Matei Radulescu

unread,
Feb 27, 2012, 4:36:34 PM2/27/12
to amrita...@googlegroups.com
Hans,
Gary Sharpe, Carlos Leung and myself have done something simpler for
1D detonation dynamics, where we integrated the location of the C+, C- and
C0 characteristics. Gary implemented this by integrating the
characteristic relations using a first order Euler in time and first order
interpolation in space.
I have attached our script, see the fold close to line 189.
Matei
________________________________________
From: amrita...@googlegroups.com [amrita...@googlegroups.com] on behalf of Hans Hornung [hgho...@gmail.com]
Sent: Monday, February 27, 2012 3:57 PM
To: amrita-ebook
Subject: ray tracing

Hans

--
You received this message because you are subscribed to the Google Groups "amrita-ebook" group.
To post to this group, send email to amrita...@googlegroups.com.
To unsubscribe from this group, send email to amrita-ebook...@googlegroups.com.
For more options, visit this group at http://groups.google.com/group/amrita-ebook?hl=en.

2stepflowin
Reply all
Reply to author
Forward
0 new messages