Diferencia entre revisiones de «Nivel piezométrico - Grupo 19»
(→Método de Éuler explícito) |
(→Método de Éuler explícito) |
||
| Línea 42: | Línea 42: | ||
===Método de Éuler explícito=== | ===Método de Éuler explícito=== | ||
| − | + | ||
%Datos | %Datos | ||
L=20;T=50;D=0.01;hp=35;ho=40; | L=20;T=50;D=0.01;hp=35;ho=40; | ||
| Línea 72: | Línea 72: | ||
[xx,tt]=meshgrid(x,t); | [xx,tt]=meshgrid(x,t); | ||
surf(xx,tt,sol) | surf(xx,tt,sol) | ||
| − | |||
===Método de Éuler implícito=== | ===Método de Éuler implícito=== | ||
Revisión del 16:00 19 may 2014
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Nivel piezométrico. Grupo 19-B |
| Asignatura | Ecuaciones Diferenciales |
| Curso | Curso 2013-14 |
| Autores | Javier Abad, Jesús Castaño, Ignacio Embid, Javier Pérez |
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
Se define el concepto de nivel piezométrico como la altura de la superficie libre de agua sobre el nivel del mar, en los acuíferos libres. En los confinados, es la altura que alcanzaría el agua en el interior de un sondeo hasta equilibrarse con la presión atmosférica. Para poder conocer la variación del nivel piezométrico se utiliza la ecuación de la conservación de la masa y la ley de Darcy.
La Ley de Darcy describe la relación entre la cantidad o la velocidad de flujo del agua, la permeabilidad del acuífero y el gradiente piezométrico (o gradiente hidráulico).
Ecuación de la conservación de la masa: [math] S ·\frac{ \partial h }{ \partial t } + div q = 0[/math]
Ley de Darcy particularizada para nuestro problema: :[math]{q}=-{K}\cdot\nabla(h(ρ))[/math]
Contenido
- 1 Cálculo del Laplaciano y ecuación diferencial en coordenadas polares
- 2 Condiciones de contorno y sistema de ecuaciones
- 3 Resolución del problema
- 4 Gráfica del comportamiento del nivel piezométrico
- 5 Resolución del problema
- 6 Gráfica del nivel piezométrico para tiempos grandes
- 7 Valor estacionario
- 8 Estimación de la capacidad de recuperación de acuífero
- 9 Método de Fourier
1 Cálculo del Laplaciano y ecuación diferencial en coordenadas polares
Cálculo del Laplaciano y ecuación diferencial en coordenadas polares: Divergencia: [math] \nabla\cdot\vec F = \frac{\partial F_x}{\partial x}+ \frac{\partial F_y}{\partial y}+ \frac{\partial F_z}{\partial z} [/math]
2 Condiciones de contorno y sistema de ecuaciones
3 Resolución del problema
3.1 Método de las diferencias finitas
3.2 Método del trapecio
4 Gráfica del comportamiento del nivel piezométrico
5 Resolución del problema
5.1 Método de Éuler explícito
%Datos
L=20;T=50;D=0.01;hp=35;ho=40;
%Discretizacion
dx=1/10;dt=dx;
%Vector x
x0=1;
x=x0:dx:L;
xint=x0+dx:dx:L-dx;
%Vector tiempo
t=0:dt:T;
%Aplicar euler explicito
N=length(x)-2;
K=(1/dx^2)*(diag(2*ones(1,N))-diag(ones(1,N-1),1)-diag(ones(1,N-1),-1));
A=(1/(dx*2))*(diag(ones(1,N-1),1)-diag(ones(1,N-1),-1));
A=diag(1./xint)*A;
h0=(40*ones(N,1));
M=length(t)-1;
F=(zeros(N,1));
F(1)=D*((hp/dx^2)-(hp/(2*dx)*1/xint(1)));
F(N)=D*((ho/dx^2)+ho/(dx*2*xint(N)));
sol(1,:)=[40;h0;40]';
h=h0;
for n=1:M
h=(((eye(N)-(dt*D)*(K-A)))*h+(dt)*(F));
sol(n+1,:)=[35 h' 40];
end
figure(1)
[xx,tt]=meshgrid(x,t);
surf(xx,tt,sol)
