-
Notifications
You must be signed in to change notification settings - Fork 1
/
tp1.in
153 lines (120 loc) · 4.53 KB
/
tp1.in
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
# ====================================================================== #
#
# Cálculo numérico avanzado, 2020, FIUBA.
#
# Grupo 4: Santiago Mosca, Santiago Pérez Raiden, Cristhian Zárate Evers.
#
# Trabajo Práctico n.° 1
#
# ====================================================================== #
#
# Archivo de input para el programa 'v2_main_tp1.py'
#
# Las líneas de ingreso deben escribirse como
# "<palabra_clave> = <valor>"
#
# Las unidades son siguientes:
# - longitudes y profundidad del dominio deben especificarse en metros
# - velocidad del río en m/s
# - descarga de contaminante en kg/día
# - decaimiento del contaminante en 1/día
# - tiempo total a simular en minutos
# - discretización espacial en metros
# - discretización temporal en segundos
#
# ====================================================================== #
# DATOS GEOMÉTRICOS
# -----------------
# Inicio y fin del dominio en coordenada 'x'
X_INI = 0.0
X_FIN = 2000.0
# Inicio y fin del dominio en coordenada 'y'
Y_INI = 0.0
Y_FIN = 350.0
# ====================================================================== #
# PARÁMETROS DE MODELADO
# ----------------------
# Velocidad de corriente en sentido 'x'
VEL = 0.8
# Profundidad del cauce
H = 10.0
# Factor de fricción [adimensional]
F = 0.04
# Coeficiente de difusión longitudinal [adimensional]
E_L = 30.0
# Coeficiente de difusión transversal [adim]
E_T = 0.3
# Descarga de contaminante
DESC_CONT = 60.0
# Constante de decaimiento del contaminante
C_DEC = 0.1
# Tiempo total a simular [min]
T_TOTAL = 120.0
# ====================================================================== #
# DISCRETIZACIÓN DEL PROBLEMA
# ---------------------------
DX = 10.0
DY = 5.0
# DISCRETIZACIÓN DEL PASO TEMPORAL. SE PUEDE ELEGIR POR UNA AUTOSELECCIÓN,
# QUE AJUSTA 'DT' SEGÚN CRITERIO DE ESTABILIDAD PARA MÉTODO EXPLÍCITO EN
# 2D (r <= 1/4). SI SE ACTIVA LA AUTOSLECCIÓN, LA ESPECIFICACIÓN DE 'DT'
# ES TODAVÍA NECESARIA PARA QUE EL PROGRAMA LEVANTE LAS VARIABLES SIN ERROR,
# AUNQUE NO SERÁ TENIDA EN CUENTA.
# AUTO_DT DEBE SER 'SI' O 'NO'
AUTO_DT = SI
DT = 2.0
# ====================================================================== #
# CONDICIONES DE BORDE
# --------------------
# Las condiciones de borde deben especificarse como 'DIR' (Dirichlet)
# o 'NEU' (Neumann).
# La condición de Dirichlet iguala el valor de la concentración a 0.
# La condición de Neumann iguala a 0 las derivadas primeras y segundas
# en la dirección normal al borde.
# Condición para el borde con coordenada 'X_INI'
CB_X_INI = DIR
# Condición para el borde con coordenada 'X_FIN'
CB_X_FIN = NEU
# Condición para el borde con coordenada 'Y_INI'
CB_Y_INI = NEU
# Condición para el borde con coordenada 'Y_FIN'
CB_Y_FIN = NEU
# ====================================================================== #
# POSICIÓN DE LA FORZANTE
# --------------------
# La forzante puede ubicarse en tres sitios predeterminados: al inicio, al
# medio o al fin de la longitud de la coordenada dada. Por simpleza de
# programación, en caso de ubicarse al inicio (o al fin) de la longitud el
# programa ubica la forzante en el nodo adyacente al borde
# correspondiente. Las opciones de selección son 'INI', 'MED' o 'FIN'.
# Ubicación de la forzante en coordenada 'x'
POS_X_FORZ = INI
# Ubicación de la forzante en coordenada 'y'
POS_Y_FORZ = INI
# ====================================================================== #
# SELECCIÓN DE MÉTODO DE SOLUCIÓN
# -------------------------------
# El valor de theta puede elegirse de los siguientes:
# 0.0 -----> Método explícito centrado
# 0.5 -----> Método de Crank-Nicolson
# 1.0 -----> Método fuertemente implícito
THETA = 0.0
# Discretización del término advectivo con 'upwinding'. <SI> o <NO>
UPWINDING = NO
# ====================================================================== #
# ESCRITURA A ARCHIVO
# -------------------
# Intervalo temporal, en minutos, para guardar la solución a archivo.
# Se utiliza también para la salida a gráfico de la solución.
T_SOL = 5
# ====================================================================== #
# SALIDA A GRÁFICO
# -------------------
# Dimensión de imagen. Si <1>, el gráfico utiliza una coordenada, 'x' o
# 'y', como variable independiente y la solución como dependiente.
# Si dimensión es <2>, entonces la salida es una vista superior del
# dominio, con la solución expresada según mapa de colores.
DIM_IMG = 2
# Nodo 'y' a lo largo del que se grafica la solución
Y_IMG = 1
# ====================================================================== #