amrita-ebook Google Group

PDF and AMR_SOL's flowin command

2009-11-24T00:18:22Z

The attached PDF, myscript.pdf, is a proof-of-concept that it is
now possible for a document to issue AMR_SOL's flowin command
and access the raw data from a computational run. The attached simply
plots the computational grid when you hit the "flowin" button,
but one can imagine crafting all sorts of custom diagnostics.

Re: restart using previously computed fields and different BC and domains

2009-11-19T20:51:06Z

Matei,

It is possible, but not with AMR_SOL's existing idioms.
What you would need to do is use AMRITA's replace mechanism
on the keyword flowin providing a back-end that
does what you want. The back-end would basically lift
grids {G1,..,Glmax} from the target input and
and add them to the current G0. It would probably

restart using previously computed fields and different BC and domains

2009-11-19T20:31:10Z

Does anybody know if it is possible to combine computation results
from one run to the next.

For example, suppose I want to use a fully established 2D cellular
detonation profile as initial condition for subsequent studies where
the boundary conditions (e.g, diffraction, pertubations, etc...) would

Re: Grid adaption on curved surface

2009-11-19T17:04:57Z

Chris,

AMR_SOL's body-fitted machinery is not actively supported,
mainly for the reasons you've run into which are quite hard
to solve and package for general use. The alternative,
which works well for stylized shapes, like your cylinder,
is to employ a level-set approach that allows internal

RE: Grid adaption on curved surface

2009-11-19T17:07:46Z

Christian,
I have used the embedded boundary technique to represent the surface of
a cylinder in problems of interaction with shocks or detonation waves.
It behaved well and permitted to use a rectangular grid. If you want to
go that route, there are examples posted on the AMRITA e-book about how

Grid adaption on curved surface

2009-11-19T15:47:53Z

I have a curvilinear grid geometry with a curved surface in it
(circular). The flow currently shows artefacts on the wall due to the
N-sided nature of the surface. I could use a very fine mesh structure
for that wall but ideally, for the mesh adaption process, the elements
on the curved wall would need to be fitted to the exact computed wall

Re: Porous Wall Boundary Condition

2009-11-04T12:58:36Z

Gary,

I wasn't suggesting you lift a far-field bc verbatim, it was more a case
of suggesting you look at the associated characteristic analysis for some
insight as what to do.

At a practical level, the first thing you should do is to implement
your own reflecting wall bc using the template from the bdy manpage.

Re: Multiple inert gas EquationSet

2009-11-04T09:51:24Z

Christian, if you are just looking at having two ideal gas EoS's with
different gammas, there is a very simple, but robust method in
A Simple Method for Compressible Multifluid Flows
SIAM J. Sci. Comput. Volume 21, Issue 3, pp. 1115-1145 (1999)
by Saurel and Abgrall. This is also non-conserative but tests show

Re: Porous Wall Boundary Condition

2009-11-04T09:42:49Z

James,
hmmm...I'm no sure the 'far-field' bc is appropriate here
a) if you run a znd wave parallel to wall, the wave would
remain 1D so wouldn't have the desired the effect.
b) there would be no way to control the "porosity"
i.e. degree of absorbtion.
We are currently looking at specifing the normal velocity

Re: Multiple inert gas EquationSet

2009-11-03T20:42:51Z

Chris,

The attached script shows how a clue-in AMRITA user could
run a simulation with two inert gases; it nominally reproduces
the results of Quirk & Karni (1994). If you take a look
at the script using amrgi, you'll see it contains the
following procedures:

(i) BubbleEquations which defines a two-component EquationSet.

Multiple inert gas EquationSet

2009-11-03T19:05:21Z

I have a shock tube model which I'd like to convert to use two
distinct inert gases (driver/driven). Looking in stdlib/equations, I
found ShockBubbleEquations which looked promising. However, the
machinery required by BCG doesn't seem to be present, so the latter
chokes on its call with a 'BCG under construction' message. (ie, no

Re: Porous Wall Boundary Condition

2009-11-01T20:22:04Z

Gary,

So you really want an absorbing boundary condition.
These are often used in the aero community for far-field
boundary conditions, and you'll find a slew of references
in the literature.

If you think of Roe's schem, applied in the transverse
wave directon, then all you're doing is filling in a

Re: Porous Wall Boundary Condition

2009-11-01T17:23:33Z

Hi James,
in the experiments, a cellular detonation in a channel runs over a
section of "porous wall"
- in the experiments the "wall" consists of a network of criss-
crossing wires with
a given "pore density" I believe. Its role is to partially absorb or
destroy transverse waves
which impact on it. We want a boundary condition which mimics

Re: Porous Wall Boundary Condition

2009-11-01T15:50:45Z

Gary,

It's not difficult to add a transpiration boundary condition,
but without knowing the details of Matei's experimental setup
it's hard to advise on what formulation the transpiration
should take.

James

Porous Wall Boundary Condition

2009-10-30T17:34:54Z

Hi all,
Matei and I wanting to simulate his porous wall detonation
experiments.
Does anyone know of or have any ideas for a boundary condition which
mimics
the "loss" effects or non-perfect reflectivity of a porous wall?
My idea was to include a non-reactive layer of complient (i.e. a gas
of some non-infinite density)