Ver el código de «El Tanque de Tormentas de Arroyofresno (Grupo 30)»

{{TrabajoED 
| El Tanque de Tormentas de Arroyofresno. Grupo 30 
| [[:Categoría:Teoría de Campos|Teoría de Campos]]
| [[:Categoría:TC25/26|2025-26]] 
| David Bretaña Blanco <br />
Alba García Martín <br />
Lorena Fernandez Soler <br />
Carlos Álvarez Gajero <br />
Jesus Gil Gutiérrez <br />
}}

<div style="font-size:120%; width:350px; padding:1em; line-height:1.6; margin-bottom:2em;">
__TOC__
</div>
==  Funcionamiento del Tanque de Arroyofresno ==
<div style="float:right; width:300px; margin-left:20px; clear:right;">
[[Archivo:xdkk.jpg|300px|thumb]]
[[Archivo:xd1.jpg|300px|thumb]]
</div>

Es un gran tanque de tormentas subterráneo situado bajo el Club de Campo de
Madrid. Diseñado para retener las “primeras aguas de lluvia” y evitar sobrecargas
al sistema de saneamiento y vertidos al Manzanares. Su volumen útil de
almacenamiento es aproximadamente 400 000 m³.<br />
Funciona como regulador hidráulico: el agua de lluvia llega por un colector de
varios km de largo y se deriva al tanque cuando el caudal de la red excede la
capacidad de desagüe o cuando el sistema operativo decide derivar para proteger
la Estación Depuradora de Aguas Residuales (EDAR).<br />
El agua no permanece un tiempo fijo dentro del tanque: depende de la intensidad
de la lluvia, de la configuración del tanque y de la capacidad de la EDAR aguas
abajo. En la práctica, el agua se mantiene hasta que la depuradora puede recibirla,
es decir, desde horas hasta pocos días según el episodio y la capacidad de
tratamiento. Arroyofresno puede llenarse completamente en 7–13 horas según
intensidad de lluvia y que solo se ha llegado al llenado máximo en contadas
ocasiones.<br />
El hormigon al estar expuesto a aguas contaminadas puede sufrir problemas de durabilidad.Los procesos responsables del deterioro son bien conocidos:
- Corrosión de armaduras por cloruros (ataque por iones Cl⁻).
Los cloruros penetran, rompen la pasivación del acero y provocan corrosión
localizada generando una pérdida de sección y aparición de fisuras por
expansión.
- Carbonatación.
Reacción del CO₂ con hidróxidos del cemento reduce el pH del hormigón
desprotege la armadura y facilita una posterior corrosión.
- Ataque por sulfatos y corrosión biológica.
Al oxidarse en la atmósfera de la cámara superior generan ácido sulfúrico, que
ataca al cemento generando una descalcificación en el material. Es uno de los
mecanismos más destructivos en estructuras de aguas residuales.
- Ataques por sales y ciclos de humedecerse/secarse
Formación de productos expansivos (ettringita u otros) que generan fisuración
y exfoliación.<br />
Las soluciones que se pueden emplear para estos problemas son variadas como por ejemplo aplicar hormigón de baja permeabilidad y mayor recubrimiento de armaduras para
retrasar la llegada de cloruros y la carbonatación o uso de armaduras inoxidables o revestidas con epoxi en ambientes muy
agresivos.<br />

Se recomiendan planes activos de conservación como: Inspecciones
periódicas, reparación precoz de revestimientos, limpieza de material
biológico, monitorización de potencial de corrosión de armaduras y contenido
de cloruros. <br />

A lo largo de España podemos encontrar otros tanques similares:<br />
- Butarque (Madrid). Otro depósito de aproximadamente 400 000 m³.
Construido con el mismo objetivo que Arroyofresno. Éste y Butarque son
dos de los mayores tanques de tormentas en Europa.<br />

- Red de tanques de Madrid. La ciudad tiene decenas de tanques (unos 38 en
el término municipal y cerca de 70 en la Comunidad de Madrid) con
capacidad total en la región de aproximadamente 1 400m³. Se usan como
sistema distributivo para evitar vertidos y regular el paso de aguas hacia la
EDAR.<br />

- Galindo (Bizkaia). Tanque/túnel de tormentas en construcción (proyecto
grande, túnel singular) para reducir vertidos en Ría de Bilbao (Capacidad
desconocida debido a construcción en curso). Ejemplo de proyecto
europeo moderno de gran escala.

== Presión y fuerzas estructurales sobre columnas ==
<big>'''2.1 Representación del campo de presiones</big><br />

En este apartado se describe cómo se representa el campo de presión hidrostática sobre la superficie lateral de una columna del tanque cuando éste se encuentra completamente lleno. El objetivo es visualizar cómo varía la presión en función de la profundidad mediante un mapa de colores aplicado sobre la superficie cilíndrica.<br />


1) Campo de presión hidrostática<br />
La presión en el interior del tanque depende únicamente de la coordenada vertical ''z'', aumentando linealmente con la profundidad:<br />

<math>P(z) = P_0 - \rho g\, z,\qquad z \in [-H,0]</math><br />

donde:<br />
* <math>P_0 = 101325\ \text{Pa}</math> es la presión atmosférica en la superficie del agua,<br />
* <math>\rho = 1000\ \text{kg/m}^3</math> es la densidad del agua,<br />
* <math>g = 9.81\ \text{m/s}^2</math> es la aceleración de la gravedad,<br />
* <math>H = 22\ \text{m}</math> es la profundidad total del tanque.<br />

Dado que ''z'' toma valores negativos hacia el fondo, el término <math>-\rho g z</math> hace que la presión aumente progresivamente al descender.<br />


2) Superficie cilíndrica de la columna<br />

La columna se considera un cilindro circular recto de radio:<br />

<math>R_c = 0.75\ \text{m}</math><br />

Para estudiar la distribución de presión sobre su superficie exterior, se parametriza el cilindro mediante:<br />

<math>
\begin{cases}
x = R_c \cos\theta, \\
y = R_c \sin\theta, \\
z = z,
\end{cases}
</math><br />

donde:
* <math>\theta \in [0,2\pi]</math> recorre toda la circunferencia,<br />
* <math>z \in [-H,0]</math> recorre la altura total sumergida.<br />


3) Representación mediante mapa de color<br />

Una vez definido el campo escalar <math>P(z)</math> y la superficie cilíndrica, se utiliza un mapa de colores para representar visualmente la presión:<br />

* Las zonas cercanas al fondo muestran colores asociados a presiones más elevadas.<br />
* Las zonas próximas a la superficie aparecen con colores correspondientes a presiones menores.<br />


4) Código MATLAB<br/>
[[Archivo:2_1.png|thumb|500px|Campo de presión hidrostática]]
{{matlab|codigo=

clear;clc;
%Datos
H=22;
rc=0.75;
Po=101325;
Da=1000;
g=9.81;

%Mallado en cilindricas
theta=linspace(0,2*pi,250);
z=linspace(-H,0,250);
[theta,z]=meshgrid(theta,z);

%Coordenadas cartesianas superficie
R=rc;
X=R*cos(theta);
Y=R*sin(theta);
Z=z;

%Campo de presion hidrostatica
P=Po-Da*g*Z;

%Representacion superficie con mapa de color
figure;
surf(X,Y,Z,P,'EdgeColor','none');
colorbar;
xlabel('X(m)');
ylabel('Y(m)');
zlabel('Z(m)');
title('Campo de presión hidrostatica');
axis equal;
}}<br />

<big>'''2.2 Gradiente de presión y representación del campo vectorial'''</big><br />

En este apartado se analiza el gradiente del campo de presión hidrostática dentro del tanque y se representa dicho gradiente sobre la superficie lateral de una de sus columnas. El objetivo es mostrar, mediante vectores, cómo cambia la presión y hacia dónde aumenta cuando el tanque está lleno.<br />


1) Gradiente del campo de presión<br />

Según la fórmula referida en el punto anterior y puesto que la presión sólo varía con la profundidad, la presión sólo cambia en la dirección vertical.<br />

El gradiente indica hacia dónde cambia más un valor. Como en este caso la presión solo depende de la profundidad ''z'', su gradiente es:<br />

<math>
\nabla P = \left( \frac{\partial P}{\partial x},\; \frac{\partial P}{\partial y},\; \frac{\partial P}{\partial z} \right)
= (0,\;0,\; -\rho_{\text{agua}}\, g).
</math><br />

Esto implica que:<br />
el vector gradiente apunta constantemente hacia abajo, con magnitud proporcional al peso específico del agua.<br />

El gradiente es por tanto un campo vectorial uniforme dirigido hacia el fondo del tanque.<br />


2) Representación del gradiente como campo vectorial<br />

Al superponer el gradiente sobre la superficie cilíndrica se obtiene un campo de vectores que muestra:<br />

* dirección constante hacia abajo, indicando que la presión aumenta a medida que descendemos en la columna,<br />
* magnitud uniforme, puesto que el gradiente es constante en todo el volumen del agua.<br />


3) Código de matlab:
[[Archivo:2_2.png|thumb|500px|Gradiente de presión]]
{{matlab|codigo=
clear;clc;
%Datos
H=22;
rc=0.75;
Po=101325;
Da=1000;
g=9.81;

%Mallado en cilindricas
theta=linspace(0,2*pi,25);
z=linspace(-H,0,25);
[theta,z]=meshgrid(theta,z);

%Coordenadas cartesianas superficie
R=rc;
X=R*cos(theta);
Y=R*sin(theta);
Z=z;

%Gradiente de presion
dPx=zeros(size(Z));
dPy=zeros(size(Z));
dPz=-Da*g*ones(size(Z));

%Representacion
figure;
surf(X,Y,Z,'FaceAlpha',0.2,'EdgeColor','none');
hold on;
quiver3(X,Y,Z,dPx,dPy,dPz,'AutoScaleFactor',0.5);
xlabel('X(m)');
ylabel('Y(m)');
zlabel('Z(m)');
title('Gradiente de presion sobre la superficie de la columna');
axis equal;
grid on;
hold off;
}}

<big>'''2.3 Fuerza total que el agua ejerce sobre una columna'''</big><br />

1) Parametrización de la superficie<br />
Para calcular la fuerza total que el liquido(agua)ejerce sobre la columna primero analizaremos la superficie de la columna<br />

Radio de la columna:<br />
<math>R_c=0.75 </math><br />

La parametrizacion de la superficie lateral será entonces:<br />

<math>
\mathbf{r}(\theta, z) =
\begin{pmatrix}
R_c \cos\theta \\
R_c \sin\theta \\
z
\end{pmatrix}
, \theta\in[0,2\pi] z\in[-H,0]
</math><br />

El vector normal exterior es:<br />
<math>\mathbf{n}(\theta, z)=\begin{pmatrix}
\cos\theta \\
\sin\theta \\
0
\end{pmatrix}
</math><br />

El diferencial de superficie sera entonces:<br />
<math>\quad dS = R_c\, d\theta\, dz,\quad </math><br />


2) El campo vectorial de la fuerza de presión es:<br />

<math>
\mathbf{F}(\theta, z) = -P(z)\,\mathbf{n}(\theta, z) = -(P_0 - \rho g z)
\begin{pmatrix}
\cos\theta \\
\sin\theta \\
0
\end{pmatrix}.
</math><br />


3) La fuerza total:<br />

<math>
 F_{\mathrm{tot}}
= \iint_{S}  F(\theta,z)\, dS
= -R_c
\int_{z=-H}^{0}
\int_{\theta=0}^{2\pi}
\left(P_0 - \rho g z \right)
\begin{pmatrix}
\cos\theta \\[4pt]
\sin\theta \\[4pt]
0
\end{pmatrix}
\, d\theta\, dz
= -R_c
\int_{z=-H}^{0}
\left(P_0 - \rho g z \right)
\left[
\int_{\theta=0}^{2\pi}
\begin{pmatrix}
\cos\theta \\[4pt]
\sin\theta \\[4pt]
0
\end{pmatrix}
\, d\theta
\right]
dz
</math><br />
Si resolvemos la segunda integral esta da 0 por lo que la <math>F_{\mathrm{tot}}=0</math>


4) Resultado final<br />

<math>F_{\mathrm{tot}}=0</math><br />
El resultado de la fuerza total ejercida es 0 ya que cada fuerza horizontal ejercida sobre la columna se anula con su simétrica, el agua ejerce la misma fuerza a la misma altura asique al sumar todas las fuerzas el conjunto resulta ser nulo.

== Infiltración y corrosión en columnas ==
'''
<big>'''3.1 Representación del campo de concentración infiltrado'''</big><br />

En este apartado se describe cómo se representa el campo de concentración de contaminantes que ha penetrado en el interior de una columna .El objetivo es visualizar de forma clara cómo los contaminantes presentes en el agua penetran en las columnas, cómo su concentración aumenta con la profundidad del depósito y cómo disminuye progresivamente desde la superficie hacia el interior.<br />


1) Concentración de contaminantes en el agua del tanque<br />
Durante los episodios de tormenta, el agua que entra en el tanque arrastra contaminantes, los cuales tienden a sedimentar hacia el fondo. Por ello, la concentración en el agua varía con la profundidad y aumenta conforme se desciende:<br />

<math>C_{\text{agua}}(z)=C_0\left(1+\alpha\,\frac{|z|}{H}\right)</math><br />

donde:<br />

* <math>C_0 = 0.3\,\text{kg/m}^3</math> concentración en superficie,<br />
* <math>\alpha=3</math> incremento de contaminación hacia el fondo,<br />
* <math>z\in[-H,0]</math> coordenada vertical, con <math>z=-H</math>.<br />

Esto refleja que el agua cercana al fondo contiene más partículas sedimentadas.<br />


2) Infiltración radial en las columnas de hormigón<br />
Las columnas cilíndricas que soportan la estructura están sometidas al contacto directo con el agua contaminada. Los contaminantes penetran en su interior siguiendo una pérdida rápida de concentración conforme se penetra en el hormigón.<br />

<math>C_{\text{col}}(\rho,z)= C_{\text{agua}}(z)\,e^{-\lambda(\rho_c-\rho)}</math><br />

donde:<br />

* <math>\rho</math> es la distancia al eje de la columna,<br />
* <math>\rho_c = 0.75\,\text{m}</math> es el radio exterior de la columna,<br />
* <math>\lambda = 10\,\text{m}^{-1}</math> es el coeficiente de atenuación.<br />

La concentración es máxima en la superficie de la columna (donde <math>\rho=\rho_c</math>) y disminuye rápidamente hacia el interior.<br />


3) Posición de la armadura<br />
En el interior de cada columna se encuentra la armadura de acero, situada a un radio:<br />

<math>\rho_a = 0.65\,\text{m}</math><br />


4) Representación en una sección vertical<br />
Para estudiar la variación vertical de la concentración infiltrada, se considera una sección vertical que pasa por el eje de la columna.<br />

En esta representación:<br />

* El eje horizontal corresponde al radio <math>\rho</math>, desde el centro hasta la superficie de la columna.<br />
* El eje vertical corresponde a la profundidad <math>z</math>, desde la superficie del agua (<math>z=0</math>) hasta el fondo (<math>z=-H</math>).<br />

El mapa de colores permite visualizar:

* El aumento de concentración hacia el fondo.<br />
* La disminución radial hacia el interior.<br />
* La localización de la armadura en <math>\rho=\rho_a</math><br />


Código MATLAB: 
[[Archivo:Punto5_1.png|thumb|500px|Sección vertical]]
{{matlab|codigo=
clear; clc; close all;

% Datos
H = 22;        
rc = 0.75;     
ra = 0.65;     
Co = 0.3;      
alpha = 3;
lambda = 10;   
%Mallado
r = linspace(0, rc, 300);
z = linspace(-H, 0, 300);
[R, Z] = meshgrid(r, z); 

% Concentración agua
Cagua =@(z) Co .* (1 + alpha .* abs(z) ./ H);

% Concentración infiltrada 
Ccol = Cagua(Z) .* exp(-lambda .* (rc - R));

%Sección vertical 
figure;
surf(R, Z, Ccol,'EdgeColor','none');
view(2);
colorbar;
xlabel('R(m)');
ylabel('Z(m)');
title('Concentración infiltrada-sección vertical');
axis tight;
hold on;
%armadura 
zline = max(Ccol(:)) * 1.05;
plot3([ra ra], [-H 0],[zline zline], 'k--', 'LineWidth', 3);
text(ra-0.08, -H+2,zline, 'Armadura', 'Color','k','FontWeight','bold','FontSize',16);
hold off;
}} <br />


5) Representación en una sección transversal<br />
Se estudia también una sección transversal de la columna en un plano horizontal situado a mitad de la altura del tanque (<math>z=-H/2</math>).<br />

En esta otra representacion:<br />

* La concentración depende únicamente del radio <math>\rho</math>.<br />

El mapa de colores muestra:

* El máximo de concentración en el borde exterior.
* Una bajada fuerte conforme se va hacia el interior.
* La circunferencia interior correspondiente a la armadura (<math>\rho=\rho_a</math>).<br />


Código MATLAB: 
[[Archivo:Punto5_2.png|thumb|500px|Sección transversal]]
{{matlab|codigo=
clear; clc; close all;

%Datos
H = 22;        
rc = 0.75;     
ra = 0.65;     
Co = 0.3;      
alpha = 3;
lambda = 10;   
%Mallado
r = linspace(0, rc, 250);
z = linspace(0, 2*pi, 250);
[R, Z] = meshgrid(r, z);

%Coordenadas cartesianas
X = R .* cos(Z);
Y = R .* sin(Z);

%Concentración mitad de la altura
z0 = -H/2;
Cagua = @(z) Co .* (1 + alpha .* abs(z) ./ H);
Ccol = Cagua(z0) .* exp(-lambda .* (rc - R));

% Seccion transversal 
figure;
surf(X,Y,Ccol,'EdgeColor','none');
view(2);
colorbar;
xlabel('X(m)');
ylabel('Y(m)');
title('Concentración infiltrada-seccion transversal');
axis tight;

hold on;
%Armadura
th=linspace(0,2*pi,200);
Xa = ra * cos(th);
Ya = ra * sin(th);
zline = max(Ccol(:));
plot3(Xa,Ya,zline* ones(size(Xa)), 'k--', 'LineWidth', 3);
text(Xa(1)*0.45, Ya(1)*1.05,zline, 'Armadura', 'Color','k','FontWeight','bold','FontSize',16);
hold off;
}}<br />

6) Interpretación física y estructural<br />
Las mayores concentraciones se encuentran en las zonas profundas del tanque.<br />


<big>'''3.2 Representación del campo del gradiente de concentración infiltrada'''</big><br />

En este apartado se analiza el gradiente del campo de concentración de contaminantes en el interior de una columna. El objetivo es identificar cómo y en qué dirección varía más intensamente la concentración dentro del material.<br />


1) Sección vertical de la columna y dominio de estudio<br />
Para analizar la penetración de contaminantes, se considera una sección vertical que pasa por el eje de la columna;<br />

* el dominio radial: <math>\rho \in [0,\rho_c]</math>,<br />
* la profundidad: <math>z \in [-H,0]</math>, con <math>H=22\ \text{m}</math>.<br />

Esta sección permite estudiar simultáneamente la variación radial  y la variación vertical.<br />


Usando las fórmulas ya indicadas en el punto anterior:<br />


2) Campo del gradiente de concentración<br />
Definido el campo escalar <math>C_{\text{col}}(\rho,z)</math>, se calcula su gradiente:<br />

<math>
\nabla C_{\text{col}} =
\left(
\frac{\partial C_{\text{col}}}{\partial \rho},
\frac{\partial C_{\text{col}}}{\partial z}
\right)
</math><br />

El gradiente representa:<br />
* la dirección en la que la concentración aumenta más rápidamente,<br />
* la intensidad del cambio local en cada punto.<br />

El componente radial muestra cómo la concentración baja desde la superficie hacia el centro de la columna, mientras que el componente vertical indica que en las zonas más profundas la concentración es mayor.<br />


4) Código MATLAB
[[Archivo:Punto6.png|thumb|500px|Sección vertical]]
{{matlab|codigo=
clear,clc;

% Datos 
H = 22;
rc = 0.75;
Co = 0.3;
alpha = 3;
lambda = 10;

%Mallado 
r = linspace(0, rc, 200);
z = linspace(0, -H, 200);
[R, Z] = meshgrid(r, z);

% Concentración agua
Cagua =@(z) Co .* (1 + alpha .* abs(z) ./ H);

% Concentración infiltrada 
R1 = abs(R);
Ccol = Cagua(Z) .* exp(-lambda .* (rc - R1));
Ccol(R1 > rc) = NaN;

%Gradiente de concentracion

dr = r(2) - r(1);
dz = z(2) - z(1);

dCdr = zeros(size(Ccol));
dCdz = zeros(size(Ccol));

dCdr(:,1:end-1) = Ccol(:,2:end) - Ccol(:,1:end-1);
dCdr(:,end) = 0;
dCdr = dCdr / dr;

dCdz(1:end-1,:) = Ccol(2:end,:) - Ccol(1:end-1,:);
dCdz(end,:) = 0;
dCdz = dCdz / dz;



%Representacion flechas gradiente
skip = 5; 
quiver(R(1:skip:end,1:skip:end), Z(1:skip:end,1:skip:end), ...
       dCdr(1:skip:end,1:skip:end), dCdz(1:skip:end,1:skip:end), 'g');

xlabel('R (m)');
ylabel('Z (m)');
title('Gradiente concentración infiltrada-sección vertical');
axis tight;
hold off;
}}<br />

* Las flechas indican la dirección de máximo incremento de concentración.<br />
* La longitud representa la magnitud del gradiente en cada punto.<br />
* Cerca de la superficie exterior aparecen vectores más intensos debido al rápido descenso de concentración hacia el interior.<br />
* En el fondo del tanque los vectores son mayores, ya que la concentración del agua es más alta.<br />


<big>'''3.3 Representación de la función de concentración en la armadura en función de la profundidad'''</big><br />

En este apartado se representa cómo varía la concentración de contaminantes que alcanza la armadura de acero de una columna en función de la profundidad. El objetivo es identificar las zonas en las que la concentración supera el umbral crítico capaz de iniciar la corrosión del acero.<br />


2) Infiltración hacia la armadura<br />
Evaluando la expresión <math>C_{\text{col}}(\rho,z) = C_{\text{agua}}(z)\,e^{-\lambda(\rho_c - \rho)}</math> en la posición del acero (<math>\rho = \rho_a</math>) se obtiene la concentración que alcanza la armadura:<br />

<math>C_{\text{acero}}(z) = C_{\text{agua}}(z)\,e^{-\lambda(\rho_c - \rho_a)}</math><br />


3) Identificación del umbral crítico de corrosión<br />
La corrosión de la armadura comienza cuando la concentración de contaminantes supera el valor mínimo necesario para activar el proceso electroquímico. Este umbral se fija en:<br />

<math>C_{\text{crit}} = 0.05\ \text{kg/m}^3</math><br />


4) Código MATLAB:<br />
[[Archivo:Punto7.png|thumb|500px|Concentración en la armadura]]
{{matlab|codigo=
clear, clc;

%Datos
H = 22;        
rc = 0.75;     
ra = 0.65;     
Co = 0.3;      
alpha = 3;
lambda = 10;
Ccrit=0.05; %Umbral critico

% vector z
z=linspace(-H,0,200);

% Concentración acero
Cagua = @(z) Co .* (1 + alpha .* abs(z) ./ H);
Ccol = Cagua(z) .* exp(-lambda .* (rc - ra));
Cacero=Ccol;

%Grafica

figure;
plot(Cacero,z, 'b', 'LineWidth', 2);  
hold on;
yline(Ccrit, 'r--', 'LineWidth', 2);    
xlabel('C_{acero}');
ylabel('Z(m)');
title('Concentración en la armadura y Umbral crítico');
legend('C_{acero}(z)', 'Umbral crítico');
set(gca,'YDir','normal')


axis padded;
hold off;
}}<br />
* La curva <math>C_{\text{acero}}(z)</math> aumenta hacia el fondo debido a la mayor concentración en las zonas más profundas del agua.<br />
* La intersección con la línea del umbral crítico permite identificar los estratos donde existe riesgo de corrosión.<br />
Las zonas donde se cumple <math>C_{\text{acero}}(z) > C_{\text{crit}}</math> corresponden a niveles de profundidad en los que la armadura podría iniciar procesos de degradación.<br />


<big>'''3.4 Representación de superficies de isoconcentración dentro de una columna'''</big><br />

En este apartado se muestra cómo se representan las superficies tridimensionales de isoconcentración de contaminantes dentro de una columna. El objetivo es visualizar cómo se distribuyen los contaminantes que se infiltran desde el agua hacia el interior del material.<br />


2) Geometría cilíndrica de la columna<br />
La columna se modela como un sólido cilíndrico recto definido por coordenadas cilíndricas <math>(\rho,\theta,z)</math>:

<math>
\begin{cases}
x=\rho\cos\theta,\\
y=\rho\sin\theta,\\
z=z,
\end{cases}
</math><br />

con:<br />

* <math>\rho\in[0,\rho_c]</math>: desde el eje hasta la superficie exterior,<br />
* <math>\theta\in[0,2\pi]</math>: <br />
* <math>z\in[-H,0]</math>: altura total sumergida.<br />

En este volumen se define el campo de concentración <math>C_{\text{col}}(\rho,z)</math>.<br />


3) Superficies de isoconcentración:<br />
Las isoconcentraciones son superficies tridimensionales en las que la concentración del contaminante toma un valor constante.<br />


En este caso se representan las superficies correspondientes a:<br />

<math>
C_{\text{col}} = 0.1,\ 0.2,\ 0.3,\ 0.4\ \text{kg/m}^3.
</math><br />


4) Representación del cilindro de la armadura<br />

Dentro de la columna existe un anillo interior donde se encuentra el acero, ubicado a un radio:<br />

<math>\rho_a = 0.65\ \text{m}</math>,<br />


Geométricamente, la armadura se visualiza como un cilindro coaxial más pequeño:<br />

<math>
\begin{cases}
x=\rho_a\cos\theta,\\
y=\rho_a\sin\theta,\\
z=z.
\end{cases}
</math><br />



5) Código MATLAB<br />
[[Archivo:Punto8.png|thumb|500px|Superficies de isoconcentración]]
{{matlab|codigo=
clear,clc; 

%Datos
H=22;
rc=0.75;
ra = 0.65;   
Co = 0.3;      
alpha = 3;
lambda = 10;
Ccrit=0.05;

%Mallado en cilindricas
r = linspace(0, rc, 90);
z = linspace(-H, 0, 90);
theta=linspace(0,2*pi,90);
[R,TH,Z]=meshgrid(r,theta,z);

%Coordenadas cartesianas superficie
X=R.*cos(TH);
Y=R.*sin(TH);

% Concentración 
Cagua = @(z) Co .* (1 + alpha .* abs(Z) ./ H);
Ccol = Cagua(z) .* exp(-lambda .* (rc - R));


%Isosuperficies

IsoS = [0.1 0.2 0.3 0.4];
colores = [0 0 1; 0 1 0; 1 1 0; 1 0 0];

figure; 
hold on;

for k = 1:length(IsoS)
    p = patch(isosurface(X, Y, Z, Ccol, IsoS(k)));
    set(p, 'FaceColor', colores(k,:), 'EdgeColor','none', 'FaceAlpha',0.35);
end


% Armadura 
theta2 = linspace(0, 2*pi, 60); 
z2 = linspace(0, H, 2);         
[TH_2, Z_2] = meshgrid(theta2, z2);  

X_a = ra * cos(TH_2);  
Y_a = ra * sin(TH_2); 

surf(X_a, Y_a, Z_2, 'FaceColor',[0.5 0.5 0.5], 'EdgeColor','none', 'FaceAlpha',0.5)
axis equal
xlabel('X'); ylabel('Y'); zlabel('Z');
title('Armadura interna')

% Estética de la figura
axis equal;
xlabel('X (m)'); ylabel('Y (m)'); zlabel('Z (m)');
title('Superficies de isoconcentracion');
view(3);
camlight headlight;
lighting gouraud;

hold off;
}}<br />
* Las isosuperficies más altas aparecen cerca del fondo y del perímetro exterior, reflejando la mayor concentración del agua profunda.<br />
* Las isosuperficies más bajas se sitúan más cerca del eje, indicando regiones menos expuestas.<br />


<big>'''3.5 Vida útil de la columna en la zona critica '''</big><br />

Queremos calcular cuanto tiempo tardara en corroerse el 25% de la barra.<br />


1) Datos<br />
Diámetro de la barra=20 mm<br />
Reducción de la barra=25%<br />
Velocidad de corrosión=0.1 mm/año<br />


2) Calculo del tiempo que tarda en reducirse<br />
Primero calcularemos cuanto se va a reducir el diametro de la columna:<br /> 
<math>
D_r = 20\left(\frac{25}{100}\right) = 5\ \text{mm}
</math><br />
Ahora sabemos que 5 mm es lo que se reduce el diametro, aplicando un calculo sencillo a continuacion sabremos entonces cuanto tiempo tarda en reducirse.<br />
<math>
t=\frac{D_r}{v_{\text{corrosión}}}
=\frac{5\ \text{mm}}{0.1\ \text{mm/año}}
=50\ \text{años}
</math><br />
Es decir, necesitamos unos 5 años aproximadamente para que se produzca esta pérdida del 25%.



<big>'''3.6 Masa total de contaminantes infiltrados en una columna'''</big><br />

En este apartado queremos calcular la masa total de contaminantes infiltrados en una columna para saber cuánto contaminante entra en el hormigón, evaluar el riesgo de corrosión y degradación, y comparar cuánta contaminación queda dentro de la columna frente al agua exterior para estimar sus necesidades de mantenimiento.<br />


1) Cálculo de la masa total de contaminantes infiltrados en una columna 
La concentración infiltrada en el hormigón es:

<math>
C_{\text{col}}(\rho,z)=C_{\text{agua}}(z)\,e^{-\lambda(\rho_c-\rho)}.
</math>


La masa total infiltrada se obtiene con la integral triple:

<math>
M=\int_{-H}^{0}\int_{0}^{\rho_c}\int_{0}^{2\pi}
C_{\text{col}}(\rho,z)\,\rho\,d\theta\,d\rho\,dz.
</math>


<math>
\int_0^{2\pi} d\theta = 2\pi.
</math>


<math>
\int_{0}^{\rho_c} \rho\, e^{-\lambda(\rho_c-\rho)} d\rho
= \frac{\lambda\rho_c - 1 + e^{-\lambda\rho_c}}{\lambda^2}.
</math>


La concentración del agua es:

<math>
C_{\text{agua}}(z)=C_0\left(1+\alpha\frac{|z|}{H}\right).
</math>


La integral en z queda:

<math>
\int_{-H}^{0} C_{\text{agua}}(z)\,dz
= C_0\,H\left(1+\frac{\alpha}{2}\right).
</math>


2) Resultado final 

<math>
M = 2\pi\,C_0\,H\left(1+\frac{\alpha}{2}\right)
\frac{\lambda\rho_c - 1 + e^{-\lambda\rho_c}}{\lambda^2}
\approx 6.7\ \text{kg}.
</math>


3) Comparación con la masa de contaminantes en un volumen equivalente de agua 
Volumen del cilindro de agua:

<math>
V = \pi\rho_c^2 H.
</math>


Con concentración media:

<math>
\overline{C}_{\text{agua}} = C_0\left(1+\frac{\alpha}{2}\right),
</math>


La masa sería:

<math>
M_{\text{agua}}
= V\,\overline{C}_{\text{agua}}
\approx 29\ \text{kg}.
</math>


4) Conclusión <br />

<math>
\text{Fracción infiltrada} =
\frac{M_{\text{col}}}{M_{\text{agua}}}
= \frac{6.7}{29}
\approx 0.23
= 23%.
</math><br />

La columna infiltra solo aproximadamente el 23% del contaminante.<br />

== Decantación y distribución de contaminantes ==

<big>'''4.1 Representación del campo de concentración de contaminantes en el agua'''</big><br />

En este apartado se describe la representación del campo de concentración de contaminantes presente en el agua, cuya variación depende de la profundidad debido a los procesos de sedimentación. El objetivo es visualizar esta distribución mediante la gráfica de la función que define la concentración y un mapa de colores aplicado sobre una pared interna del tanque.


1) Representación gráfica de la función <math>C_{\text{agua}}(z)</math>:<br />
Dado que el campo de concentración depende únicamente de la coordenada vertical <math>z</math>, se representa inicialmente como una curva que muestra su variación desde la superficie del agua (<math>z=0</math>) hasta el fondo del tanque (<math>z=-H</math>).<br />


Código MATLAB:
[[Archivo: Punto11_1.png|thumb|500px|Función <math>C_{\text{agua}}(z)</math>]]
{{matlab|codigo=
clear; close all; clc;

% Datos
Co = 0.3;       
alpha = 3;
H = 22;         
L = 290;        

% Grafico funcion Cagua(z)
z = linspace(-H, 0, 200);
Cagua = Co * (1 + alpha * abs(z)/H);

figure;
plot(Cagua,z, 'LineWidth', 2);
ylabel('Z (m)');
xlabel('C_{agua}(z) (kg/m^3)');
title('Concentración de contaminantes en el agua');
grid on;
}}<br />

Esta gráfica permite reconocer:<br />
* la concentración mínima en la superficie,<br />
* el incremento gradual hacia el fondo,<br />


2) Mapa de colores sobre una pared del tanque:<br />
El mismo campo escalar <math>C_{\text{agua}}(z)</math> se proyecta sobre una pared vertical del tanque, como un plano rectangular de longitud <math>L = 290 \ \text{m}</math> y altura <math>H = 22 \ \text{m}</math>.<br />


Codigo MATLAB:
[[Archivo: Punto11_2.png|thumb|500px|Mapa de colores sobre una pared]]
{{matlab|codigo=
clear; close all; clc;

% Datos
Co = 0.3;       
alpha = 3;
H = 22;         
L = 290;    

% Mapa de colores sobre pared del tanque
x = linspace(0, L, 60);
z2 = linspace(-H, 0, 60);  
[X, Z] = meshgrid(x, z2);

Cagua = Co * (1 + alpha * abs(Z)/H);

figure;
pcolor(X, Z, Cagua);
shading interp;
colorbar;
xlabel('X (m)');
ylabel('Z (m)');
title('Mapa de colores concentración C_{agua}(z)');
}}<br />
* Los colores intensos cerca del fondo indican concentraciones elevadas.<br />
* Los tonos suaves próximos a la superficie representan concentraciones bajas.<br />


<big>'''4.2 Gradiente de concentración de contaminantes en el agua'''</big><br />

En este apartado se describe la representación del gradiente del campo de concentración de contaminantes en el agua. El objetivo es visualizar la dirección y la magnitud del flujo asociado a la variación vertical de la concentración.<br />


1) Gradiente del campo de concentración:<br />

Como la concentración solo varía con <math>z</math>, el gradiente se expresa como:<br />

<math>
\nabla C_{\text{agua}} = \left(0,\ 0,\ \frac{dC_{\text{agua}}}{dz}\right),
</math><br />

donde su componente vertical es constante:<br />

<math>
\frac{dC_{\text{agua}}}{dz} = -\,\frac{C_0\,\alpha}{H}.
</math><br />

El signo negativo indica que la concentración aumenta al descender (hacia valores más negativos de <math>z</math>). Representa el flujo que tendería a desplazarse hacia la superficie.


2) Proyección sobre una pared vertical del tanque:<br />

Para visualizar espacialmente el gradiente, se proyecta este campo vectorial sobre una pared interna del tanque como un plano rectangular definido por:<br />

* longitud <math>L = 290\ \text{m}</math>,<br />
* altura <math>H = 22\ \text{m}</math>,<br />
* coordenadas <math>(x,z)</math>.<br />

Dado que el campo depende únicamente de <math>z</math>, el gradiente mantiene el mismo valor en todos los puntos de la pared.<br />


3) Código MATLAB:<br />
[[Archivo: Punto12.png|thumb|500px|Gradiente de concentración de contaminante en el agua]]
{{matlab|codigo=
clear; close all; clc;

%Datos
Co = 0.3;     
alpha = 3;
H = 22;       
L = 290;
 

% Malla pared
x = linspace(0, L, 60);
z = linspace(-H, 0, 50);   
[X, Z] = meshgrid(x, z);

% Concentración
Cagua = Co * (1 + alpha*abs(Z)/H);

% Gradiente analítico 
dCdz =  -Co*alpha/H;  

% Componentes del gradiente
U = zeros(size(X));           
V = dCdz * ones(size(X)); 

% Representación
figure;
pcolor(X, Z, Cagua);
shading interp;
colorbar;
hold on;

skip = 3;

quiver(X(1:skip:end,1:skip:end), ...
       Z(1:skip:end,1:skip:end), ...
       U(1:skip:end,1:skip:end), ...
       V(1:skip:end,1:skip:end), ...
       'k', 'AutoScale', 'on');

xlabel('x (m)');
ylabel('z (m)');
title('Gradiente de concentración de contaminante en el agua');
set(gca,'YDir','normal');

hold off;
}}<br />

* Las zonas cercanas al fondo muestran colores correspondientes a concentraciones elevadas.<br />
* Las zonas próximas a la superficie presentan colores asociados a valores bajos.<br />
* Todos los vectores apuntan verticalmente hacia abajo, indicando la dirección en la que la concentración aumenta.<br />


<big>''' 4.3 Superficie de isoconcentración en el agua y representación 3D'''</big><br />

En este apartado se describe cómo se representan las superficies de isoconcentración del contaminante en el agua. El objetivo es visualizar cómo se distribuyen los contaminantes en el volumen de agua y cómo aumenta su concentración con la profundidad debido al proceso de sedimentación.<br />


2) Superficies de isoconcentración en el interior del tanque:<br />
Puesto que <math>z = 0</math> corresponde a la superficie y <math>z = -H</math> al fondo, el término <math>\frac{|z|}{H}</math> crece linealmente con la profundidad. De este modo, el valor de <math>C_{\text{agua}}(z)</math> es mínimo en la superficie e incrementa progresivamente conforme se desciende.<br />


Para visualizar este campo tridimensional, se representan isosuperficies:<br />

<math>
C_{\text{agua}}(x,y,z) = C^\*
</math><br />

donde <math>C^\*</math> es un nivel de concentración concreto.<br /> 

En este caso se consideran los valores:<br />

<math>
C^\* = 0.4,\; 0.6,\; 0.8,\; 1.0\ \text{kg/m}^3
</math><br />


3) Código MATLAB:<br />
[[Archivo: Punto13.png|thumb|500px|Superficies de isoconcentración en el agua]]
{{matlab|codigo=
clear; close all; clc;

%Datos
L = 290;   
W = 140;  
H = 22;   
Co = 0.3;      
alpha = 3;     

% Mallado 3D 
x = linspace(-L/2, L/2, 40);
y = linspace(-W/2, W/2, 40);
z = linspace(-H, 0, 40);   
[X, Y, Z] = meshgrid(x, y, z);

% Concentración
Cagua = Co .* (1 + alpha .* abs(Z) ./ H);

% Valores de isosuperficie
IsoS = [0.4 0.6 0.8 1.0];

% Figura
figure; 
hold on;
title('Superficies de isoconcentración');
xlabel('x (m)'); ylabel('y (m)'); zlabel('z (m)');

colors = lines(length(IsoS));

for k = 1:length(IsoS)
    p = patch(isosurface(X, Y, Z, Cagua, IsoS(k)));
    set(p,'FaceColor',colors(k,:), ...
          'EdgeColor','none', ...
          'FaceAlpha',0.35);
end

colormap turbo;
colorbar;  
grid on;
axis equal;
view(3);
camlight; 
lighting gouraud;

hold off;
}}<br />

* Cada isosuperficie es el conjunto de puntos del tanque donde el agua presenta exactamente la misma concentración de contaminantes.<br />
* Como la concentración depende únicamente de ''z'', cada superficie de isoconcentración adopta la forma de una capa horizontal cerrada.
* Cuanto mayor es el valor de <math>C^\*</math>, más profunda se encuentra la superficie correspondiente, ya que las concentraciones más altas aparecen cerca del fondo.<br />
* Cada superficie mantiene la geometría del tanque, pero a diferentes niveles ''z''.<br />


<big>'''4.4 Masa total de contaminantes en el agua'''</big><br />

La masa total de contaminantes presentes en el agua del tanque se determina integrando el campo de concentración del contaminante sobre el volumen ocupado por el agua. Dado que la concentración depende únicamente de la profundidad se puede simplificar la integral. Este procedimiento permite estimar la cantidad total de contaminantes disueltos en el tanque.<br />


1) Volumen de agua en el tanque<br />
El tanque se modela como un prisma rectangular. El área de su sección horizontal es:<br />

<math>
A = L \times W
</math><br />

con:<br />

* <math>L = 290\,\text{m}</math><br />
* <math>W = 140\,\text{m}</math><br />

De donde:<br />

<math>
A = 290\,\text{m} \times 140\,\text{m} = 40\,600\,\text{m}^2
</math><br />

Si el fondo del tanque está en <math>z = -22\,\text{m}</math> y la superficie del agua en <math>z = -7\,\text{m}</math>, la altura de agua acumulada es:<br />

<math>
H_{\text{agua}} = |-7 - (-22)| = 15\,\text{m}
</math><br />

y el volumen total de agua:<br />

<math>
V_{\text{total}} = A \times H_{\text{agua}}
                   = 40\,600\,\text{m}^2 \times 15\,\text{m}
                   = 609\,000\,\text{m}^3
</math><br />


2) Integración para la masa total de contaminantes:<br />
La masa total de contaminantes <math>M_{\text{cont}}</math> en el volumen de agua se obtiene mediante:<br />

<math>
M_{\text{cont}}
= \iiint_{V_{\text{agua}}} C_{\text{agua}}(z)\, dV
= \int_{-H}^{-7} C_{\text{agua}}(z)\, A\, dz
</math><br />

Sustituyendo el campo de concentración:<br />

<math>
M_{\text{cont}}
= A\, C_0
\int_{-22}^{-7}
\left(1 - \alpha\frac{z}{H}\right)\, dz
</math><br />

Con los valores numéricos:<br />

<math>
M_{\text{cont}}
= (40\,600\,\text{m}^2)\,
(0.3\,\text{kg/m}^3)\,
\int_{-22}^{-7}
\left(1 - 3\frac{z}{22}\right)\, dz
</math><br />

La integral indefinida resulta:<br />

<math>
\int \left(1 - \frac{3}{22}z\right)\, dz
= z - \frac{3}{44}z^2
</math><br />

Evaluando entre <math>z=-22</math> y <math>z=-7</math>:<br />

<math>
\left[z - \frac{3}{44} z^2\right]_{-22}^{-7}
= (-7 - \tfrac{147}{44}) - (-22 - 33)
\approx 44.6591\,\text{m}
</math><br />

Finalmente, multiplicando por el coeficiente constante:<br />

<math>
M_{\text{cont}}
= 12\,180\,\text{kg} \times 44.6591
</math><br />

<math>
\boxed{
M_{\text{cont}} \approx 5.44\times 10^{5}\,\text{kg}
}
</math><br />


<big>'''4.5 Pasadas 24 Horas'''</big><br />

1) Masa total de contaminantes en el tanque:<br />
El tanque se modela como un prisma rectangular con:<br />
* Longitud: ''L'' = 290 m  <br />
* Anchura: ''W'' = 140 m  <br />
* Profundidad: ''H'' = 22 m  <br />

La concentración de contaminantes en el agua depende únicamente de la profundidad y aumenta hacia el fondo por sedimentación, entonces la masa total se obtiene integrando el campo de concentración en el agua a lo largo de la profundidad y multiplicando por el área efectiva del fondo del tanque.<br />


2) Área efectiva del fondo:<br />
Las columnas estructurales tienen radio:<br />

<math>\rho_c = 0.75\ \text{m}</math><br />

y están dispuestas en una retícula con separación de 15 m. <br />

El número aproximado de columnas es:<br />

<math>N \approx \frac{L}{15}\,\frac{W}{15}</math><br />

La superficie ocupada por todas las columnas es:<br />

<math>A_{\text{col}} = N\,\pi \rho_c^2</math><br />

El área libre disponible para la acumulación de lodo es:<br />

<math>A_{\text{ef}} = L\,W - A_{\text{col}}</math><br />


3) Masa de contaminantes depositada tras 24 horas:<br />
Después de un periodo de decantación de 24 h, el 60 % de la masa total de contaminantes se deposita en el fondo<br />

<math>M_{\text{dep}} = 0.6\,M_{\text{tot}}</math><br />


4) Volumen del lodo formado:
El material depositado forma en el fondo un lodo con densidad aparente<br />

<math>\rho_{\text{lodo}} = 1400\ \text{kg/m}^3</math><br />

El volumen total de lodo es<br />

<math>V_{\text{lodo}} = \frac{M_{\text{dep}}}{\rho_{\text{lodo}}}</math><br />


5) Espesor promedio de la capa de lodo:<br />
El espesor medio de la capa de lodo se obtiene dividiendo el volumen depositado entre el área efectiva<br />

<math>e = \frac{V_{\text{lodo}}}{A_{\text{ef}}}</math><br />

Este valor representa la altura media del sedimento acumulado en un único episodio de lluvia.<br />


6) Acumulación anual y criterio de limpieza:<br />
El tanque opera en promedio 8 veces al año. Si no se realiza limpieza, el espesor acumulado tras un año es<br />

<math>e_{\text{anual}} = 8\,e</math><br />

El criterio de mantenimiento establece que la limpieza debe realizarse cuando el espesor supera:<br />

<math>e_{\text{crit}} = 0.30\ \text{m}</math><br />

La comparación entre ambos valores permite determinar si la limpieza anual es suficiente o si se requiere una frecuencia mayor.<br />


7) Código de MATLAB:<br />
{{matlab|codigo=
niv = input('nivel maximo de agua: ');
dens = input('introduce la densidad del lodo: ');

% Integral de concentracion
Ct = integral(concentracion, -22, niv);
fprintf('Masa por unidad de superficie = %.4f kg/m^2\n', Ct);

% Geometria del estanque 
L = 290; W = 140;
A_total = L * W;   

% Columnas
r = 0.75;
n_col = (L/15) * (W/15);
A_col = n_col * pi * r^2;

A_ef = A_total - A_col;  

% Masa total 
Mtot = Ct * A_ef;
fprintf('Masa total de contaminante = %.2f kg\n', Mtot);

% Masa depositada (60%)
Mdep = 0.6 * Mtot;
fprintf('Masa depositada (60%%) = %.2f kg\n', Mdep);

% Volumen del lodo   
Vdep = Mdep / dens;
fprintf('Volumen del lodo = %.2f m^3\n', Vdep);

% Espesor medio 
espesor = Vdep / A_ef;
fprintf('Espesor medio del lodo = %.4f m\n', espesor);

% Comprobacion limpieza

lodo_acumulado = espesor * 8;
if lodo_acumulado < (30/100)
    disp ('es suficiente una limpiez anual.')
else
    disp('Se requiere una limpieza más frecuente.');
end
}}<br />
{{matlab|codigo=
% Concentracion
function [Cz] = concentracion (z)
Cz= @(z) 0.3*(1+3*(abs(z)/22));
}}<br />

== Póster para la presentación ==
[[Archivo:El Tanque de Tormentas de ArroyoFresno (Grupo 30) A1 DEFINITIVO.png|thumb|700px| Póster para la presentación]]


[[Categoría:Teoría de Campos]]
[[Categoría:TC25/26]]