
#1
Sep111, 08:40 AM

P: 41

1. The problem statement, all variables and given/known data
I'm going to model the velocity of a river stream v(y,z) as a function of the location (y,z) in the crosssection of the river. Here y represents the distance from one of the shores and z the depth of the river. I have done measurements of the depth z on different locations y, so I can find an approximation for a function z(y) and hence find the profile of the river. I have also done measurements of the velocity on different points in the river, so now I wish to find a function v(y,z) that describes the speed at every point in the crosssection. This will give me a way to find the average speed by integrating over the function for the crossarea. 2. Relevant equations Basically, through a number of assumptions and approximations, the problem can be reduced to a set of two differential equations: [itex]\eta_{y}\frac{\partial^{2} v(y,z)}{\partial y^2} + \eta_{z}\frac{\partial^{2} v(y,z)}{\partial z^2} = \rho g Sin(\alpha)[/itex] By assuming by turns [itex]\eta_{y} = 0[/itex] and [itex]\eta_{z} = 0[/itex] I can find functions v(y) and v(z) that describe a kind of an average velocity for any given y or given z. The problem now is the integration of these two into one function v(y,z). Given boundary conditions: * v(0,z) = 0 * v(W,z) = 0 > W represents max width * v(y,0) = 0 > (y,0) is a point on the bottom * [itex]\frac{\partial v}{\partial z}[H] = 0 [/itex] > change of speed is zero in for the maximum height H. 3. The attempt at a solution I've found solved the equations for the assumptions [itex]\eta_{y} = 0[/itex] and [itex]\eta_{z} = 0[/itex] and using the boundary conditions found: [itex]v(y) = \frac{\rho g Sin[\alpha]}{2\eta_{y}} (Wy)y[/itex] [itex]v(z) = \frac{\rho g Sin[\alpha]}{2\eta_{z}} (2Hz)z[/itex] The question now is whether there is an easier way to find v(y,z) by directly solving it mathematically, possibly by combining solutions for v(y) and v(z)? I can see that a linear combination of the two would work, but I can't see how I can find the general solution of the function. I have tried to find a solution using Mathematica, but by filling in the boundary conditions it reduced to the function v(y) that is independent of z. Here is what my code in Mathematica more or less looks like:




#2
Sep111, 04:42 PM

P: 1,666

Have not gone through it in detail but the first thing I notice is that you seem to be misinterpreting the line:
[tex]\left\{\left\{f\to \text{Function}\left[\{y,z\},\frac{g y^2 \rho \text{Sin}[\alpha ]}{2 a}+C[1]\left[\frac{\sqrt{a b} y}{a}+z\right]+C[2]\left[\frac{\sqrt{a b} y}{a}+z\right]\right]\right\}\right\}[/tex] The C[1] and C[2] are function names and not constants c1 and c2 multiplied by those values in brackets. You realize that or no? So the first one is an arbritrary function of a single variable named C1 at the value given in brackets. Or am I insulting your intelligence by saying that? Sorry if so. 



#3
Sep211, 02:59 AM

P: 41

Ah, thanks, I actually didn't know that! Indeed, now the solution seems to make more sense to me. It's identical to the wave function except that it is nonhomogeneous. Now I presume I can find the functions C[1] and C[2] in the same way as the wave function and then add the factor [itex] \frac{gy^2\rho Sin[\alpha]}{2\eta_{y}} [/itex] to it?
EDIT: I have simply no idea how to solve for the functions C[1] and C[2] using the boundary conditions. I looked up the solution to the wave equation but the approach there is quite different since it is possible to use the boundary conditions of amplitude and speed at t=0 for x=0. I lack knowledge of differential equations to solve for the functions. Also, if I try to insert the boundary conditions in DSolve in Mathematica, it doesn't give me any solution. 



#4
Sep211, 06:45 AM

P: 1,666

Model involving SOD of two variables[tex]\frac{\partial^2 v}{\partial y^2}+\frac{\partial^2 v}{\partial z^2}=1/2,\quad 0\leq y\leq W,\quad 0\leq z\leq H[/tex] [tex]\begin{array}{cc} v(0,z)=0 & v(W,z)=0 \\ v(y,0)=0 & \frac{\partial v}{\partial z}\biggr_{z=H}=0 \end{array} [/tex] Or something simple like that. So I'll set up a grid say for now 10x10, do all the difference equations on it with a spacing of one, then solve for all the internal points so that's what 64 linear equations in 64 variables. Alright, I did this quick so not entirely confident of the code but the backsubstitution on the results using a sample point in the grid are close to 1/2 which is the first numeric result and the second result is the derivative with respect to z along the top edge at a select point. It's 0.07. Not zero but pretty close to what we want. So then maybe the code is workable so if you want, you can start adding complexity to the PDE and make the grid spacing smaller. When you're done with that, it's not then hard to integrate the interpolation function as you indicated above.
Plot3D[thev[x,y],{x,1,10},{y,1,10}] which is below. Note at the top edge with respect to z, it's almost flat which I assume is the first partial set to zero there. 



#5
Sep211, 05:10 PM

P: 41

Thanks a lot. The graph indeed shows the expected result, with the velocity highest midway between the shores and decreasing in the neighbourhood of the shores and bottom.
I'm going to complete it by adding extending the formula and making the grid larger. Also, since I have real data for a number of points (around 20 I believe), I could insert those in my grid and then try to find an Interpolation through all the points. I'll show some results if I have some more. 


Register to reply 
Related Discussions  
A binomial problem involving 2 different random variables.  Set Theory, Logic, Probability, Statistics  6  
how to solve two equations involving complex variables??  Linear & Abstract Algebra  11  
proof involving FTC to two variables  Calculus  1  
[Mathematica Help] An integration involving many vector variables  Math & Science Software  1  
A Probability Problem Involving 6 Random Variables  Calculus & Beyond Homework  1 