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

[FELIPE ORELLANA ROVIROSA]

 Hi Elbridge,

       If the periodic BC you refer to are reflective about the wall (have
wavelength equal to the box horizontal length), that implies that the walls
have zero net heat flux.

       Check the cookbook 2D convection example. As far as I remember, it
is set exactly for that case. Homogeneous Newman null on the side walls,
and homogeneous Dirichlet on the base and top. I dont remember exactly, but
check it out.

   If you mean a Temp that is periodic inhomogeneous function of 'y' over
the walls, that will need the command function.

cheers,

Felipe

On Tue, Jun 14, 2016 at 5:26 PM, Elbridge Gerry Puckett <
egp at math.ucdavis.edu> wrote:

> 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
>
>
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20160614/edbd31da/attachment.html>