[aspect-devel] How does one set the Temperature BC in a 2D Periodic Box?

[Elbridge Gerry Puckett]

Hi!

We're running a 2D computation in a 3 X 1 rectangle with period BCs on 
the left and right walls, but NOT on the top and bottom.

How do I set the temperature BC on the left and right walls to be periodic?

I'm including the relevant parts of my parameter file below.

Here is the Geometry model:

    subsection Geometry model
       set Model name = box
       subsection Box
         set X extent = 3.0
         set X repetitions = 192
         set Y extent = 1
         set Y repetitions = 64
         set X periodic = true
       end
    end

Here is the Boundary temperature model

    subsection Boundary temperature model

       set Model name = box

       subsection Box
         set Bottom temperature = 1
         set Left temperature   = 0
         set Right temperature  = 0
         set Top temperature    = 0
       end
    end

Here are the Model settings

    subsection Model settings

       set Remove nullspace = linear x momentum

       set Fixed temperature boundary indicators   = bottom, top # left,
    right

       set Zero velocity boundary indicators       = top, bottom
       set Prescribed velocity boundary indicators =
       set Tangential velocity boundary indicators = left,
    right                <----------------------------- THIS SEEMS
    INCORRECT, IS IT?

       set Include adiabatic heating               = false
       set Include shear heating                   = false
    end


Thanks!

- Gerry & Harsha

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20160614/3f18b647/attachment.html>