Por Daniel Duque
Las condiciones de contorno (BC) son vitales en una simulación:
Corresponden a distintas situaciones físicas
Afectan la propia existencia de soluciones (lo que se traduce en "inestabilidad" en los códigos = cuelgues)
Lo sutil es que las condiciones de velocidad y presión están relacionadas. P.e. las típicas no-slip:
0/p:
pared
{
type zeroGradient;
}
0/U:
pared
{
type noSlip;
}
En general:
0/p:
pared
{
type zeroGradient;
}
0/U:
pared
{
type fixedValue;
value uniform ( u v w );
}
En un caso sencillo (Poiseuille, flujo "libre" ):
0/p:
entrada
{
type fixedValue;
value uniform p1;
}
0/U:
entrada
{
type zeroGradient;
}
Alterar nuestra caja genérica para este caso
Comenzar con este blockMeshDict
Para el cilindro, nos gustaría especificar directamente la velocidad de entrada. Entonces:
0/p:
entrada
{
type zeroGradient;
}
salida
{
type fixedValue;
value uniform 0;
}
0/U:
entrada
{
type fixedValue;
value uniform (1 0 0);
}
salida
{
type zeroGradient;
}
Comenzar con este blockMeshDict
0/U: ```c++ salida { type inletOutlet; inletValue uniform (0 0 0); value uniform (1 0 0); }