Scribd
Upload
70 views

Computational Physics Test

This document provides instructions for a mid-semester exam on finite difference methods for solving partial differential equations. Students are asked to solve two problems analytically and numerically. The first problem involves solving the Laplace equation to model hydraulic potential in a drainage basin, and comparing analytical and numerical solutions with varying grid sizes. The second problem involves solving the water hammer equation to model pressure and velocity in a pipe, and similarly comparing analytical and numerical solutions. Students must submit their work by May 22, 2014.

Uploaded by

naufal Riyandi
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
70 views

Computational Physics Test

This document provides instructions for a mid-semester exam on finite difference methods for solving partial differential equations. Students are asked to solve two problems analytically and numerically. The first problem involves solving the Laplace equation to model hydraulic potential in a drainage basin, and comparing analytical and numerical solutions with varying grid sizes. The second problem involves solving the water hammer equation to model pressure and velocity in a pipe, and similarly comparing analytical and numerical solutions. Students must submit their work by May 22, 2014.

Uploaded by

naufal Riyandi
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

Universitas Padjadjaran Fisika Komputasi Dr.

Irwan Ary Dharmawan


Program Studi Fisika Semester Genap 2013/2014 [email protected]

Ujian Tengah Semester


Metoda Beda Hingga untuk Laplace, Poisson, Wave Equation - Deadline 22 Mei 2014

1. The hydraulic potential within a small drainage basin that lies in a shallow valley can
be written in two-dimensional Laplace Equation

∂ 2u ∂ 2u
+ = 0, 0 < x < L, 0 < y < z0
∂x2 ∂y 2
with the boundary conditions

u(x, z0 ) = gz0 + gcx

∂u(0, y) ∂u(L, y)
= = 0,
∂x ∂x
and
∂u(x, 0)
=0
∂y
where u(x, y) is the hydraulic potential, g is the acceleration due to gravity, and c gives
the slope of the topography. The condition ∂u(L,y)
∂x
= 0, and ∂u(x,0)
∂y
= 0 specify a no-flow
condition through the bottom and sides of the aquifer. The condition ∂u(0,y)
∂x
= 0 ensures
symmetry about the x = 0 line. The boundary condition u(x, z0 ) = gz0 + gcx gives the
fluid potential at the water table, where z0 is the elevation of the water table above the
standard datum. The term gcx expresses the increase of the potential from the valley
bottom toward the water divide. On average it closely follows the topography.
From above problem, show that the analytical complete solution is
∞
gcL 4gcL X cos [(2m − 1)πx/L] cosh [(2m − 1)πy/L]
u(x, y) = gz0 + − 2
2 π m=1 (2m − 1)2 cosh [(2m − 1)πz0 /L]

Develop the numerical solution based on the Finite Difference method, compare the
results between analytical and numerical results, analyse the error for different grid
space.

http://geophys.unpad.ac.id/~iad/lect/fiskom/uts.pdf
Universitas Padjadjaran Fisika Komputasi Dr. Irwan Ary Dharmawan
Program Studi Fisika Semester Genap 2013/2014 [email protected]

2. The pressure and velocity distribution from water hammer in a pipe without friction are
given by the equations
∂p ∂u
= −ρc2
∂t ∂x
and
∂u 1 ∂p
=−
∂t ρ ∂x
where p(x, t) denotes the pressure perturbation, u(x, t) is the velocity perturbation, c is
the speed of sound in water, and ρ is the density of water. These two first-order partial
differential equations can be combined to yield

∂ 2p 2
2∂ p
=c
∂t2 ∂x2
Find the analytical solution to this partial differential equation if p(0, t) = p0 , and
u(L, t) = 0, and the initial conditions are p(x, 0) = p0 , and u(x, 0) = u0 . Develop the
numerical solution based on the Finite Difference method, compare the results between
analytical and numerical results, analyse the error for different grid space.

http://geophys.unpad.ac.id/~iad/lect/fiskom/uts.pdf

You might also like