par_contour = {'cmap':'cool', 'levels':20}
#
# Datos del problema
#
longitud_x = 1.0 # meters
longitud_y = 1.0
k = 0.01 #0.001 # W/m.K
Nx = 41 # Número de nodos 81 x 81
Ny = 41
dt = 0.0001
Tmax = 400 # 3000
#
# Definición del dominio y condiciones de frontera
#
placa = Rectangle(longitud_x, longitud_y)
placa.boundaryConditions(dirichlet = {'LEFT':-0.0, 'RIGHT':0.0}, neumman ={'BOTTOM':0, 'TOP':0})
#
# Creamos la malla y obtenemos datos importantes
#
malla = placa.constructMesh(Nx, Ny)
ivx, ivy, _ = malla.bounds(bi = 1, ei = Nx-1,
bj = 1, ej = Ny-1)
nx = malla.nx # Número de nodos
ny = malla.ny # Número de nodos
nvx = malla.vx # Número de volúmenes
nvy = malla.vy # Número de volúmenes
dx = malla.dx # Tamaño de los volúmenes
dy = malla.dy # Tamaño de los volúmenes
xv, yv, _ = malla.coordinatesMeshFVM()
xd, yd, _ = malla.coordinatesMeshFDM()
xg, yg = np.meshgrid(xv, yv, sparse=False)
longitud_x = 1.0 # meters
longitud_y = 2.0
TL = 0 # °C
TR = 0 # °C
TB = 0 # °C
TT = 1 # °C
k = 1 # W/m.K
Nx = 11 # Número de nodos
Ny = 21
#
# Definición del dominio y condiciones de frontera
#
placa = Rectangle(longitud_x, longitud_y)
placa.boundaryConditions(dirichlet={
'LEFT': TL,
'RIGHT': TR,
'BOTTOM': TB,
'TOP': TT
})
#
# Creamos la malla y obtenemos datos importantes
#
malla = placa.constructMesh(Nx, Ny)
ivx, ivy, _ = malla.bounds(bi=1, ei=Nx - 1, bj=1, ej=Ny - 1)
nx = malla.nx # Número de nodos
ny = malla.ny # Número de nodos
nvx = malla.vx # Número de volúmenes
nvy = malla.vy # Número de volúmenes
dx = malla.dx # Tamaño de los volúmenes
dy = malla.dy # Tamaño de los volúmenes
#
# Imprimimos los datos del problema (nicely)
longitud_y = 2.0
TL = 100 # °C
FR = 0 # Flujo a la derecha
FB = 0 # Flujo abajo
TT = 50 # °C
k = 1 # W/m.K
Nx = 11 # Número de nodos
Ny = 21
#
# Definición del dominio y condiciones de frontera
#
placa = Rectangle(longitud_x, longitud_y)
placa.boundaryConditions(dirichlet={
'LEFT': TL,
'TOP': TT
},
neumman={
'BOTTOM': FB,
'RIGHT': FR
})
#
# Creamos la malla y obtenemos datos importantes
#
malla = placa.constructMesh(Nx, Ny)
ivx, ivy, _ = malla.bounds(bi=1, ei=Nx - 1, bj=1, ej=Ny - 1)
nx = malla.nx # Número de nodos
ny = malla.ny # Número de nodos
nvx = malla.vx # Número de volúmenes
nvy = malla.vy # Número de volúmenes
dx = malla.dx # Tamaño de los volúmenes
dy = malla.dy # Tamaño de los volúmenes
#