How Do You Numerically Solve Coupled PDEs with Continuity Constraints?

  • Thread starter Thread starter keyns
  • Start date Start date
  • Tags Tags
    Coupled
keyns
Messages
3
Reaction score
0
I'm new, hi all.

I have two coupled equations, one of which is continuity. Basically, my problem comes down to the following system:

(1) u=f(v) (similarly, v=g(u). Here, u and v are the components of a vector field, ie u=u(x,y) and v=v(x,y).

(2) Continuity: \nabla \cdot \textbf{u} = 0 or u_{x}+v_{y}=0

From here, I can find the following expressions

u_{x} = -g_{y} \left( u \right)
v_{y} = -f_{x} \left( v \right)

Which I think leaves an equation of the form G \left( u,u_{x},u_{y} \right)=0 and F \left( v,v_{x},v_{y} \right)=0

It seems to me my original problem has two variables and I have two equations. I think this should be solvable, but I don't know how. Any help please? Thanks in advance!

--edit-- p.s. I'm looking for a numerical (discrete) solution.
 
Last edited:
Physics news on Phys.org
keyns said:
I'm new, hi all.

I have two coupled equations, one of which is continuity. Basically, my problem comes down to the following system:

(1) u=f(v) (similarly, v=g(u). Here, u and v are the components of a vector field, ie u=u(x,y) and v=v(x,y).

If u = f(v) and v = g(u), don't you just have a (possibly nonlinear) system of equations to solve for v and u? What is the continuity condition (2) for? Do you get more than one solution by solving (1) and you need (2) to choose a solution of interest?
 
Mute said:
If u = f(v) and v = g(u), don't you just have a (possibly nonlinear) system of equations to solve for v and u? What is the continuity condition (2) for? Do you get more than one solution by solving (1) and you need (2) to choose a solution of interest?

Actually I have only one relation for u and v that I can write it as u(v) or v(u). Sorry for the confusion. Otherwise you would be right. To clarify my equations:

(1) A relation for u and v (if I have u, I have v and vice versa)
(2) A relation for u_{x} and v_{y} (if I have u_{x}, I have v_{y} and vice versa

--edit-- Which then, after some rewriting, leads to two ODE's G and F as stated before. I just don't know how to solve those.
 
Last edited:
There is the following linear Volterra equation of the second kind $$ y(x)+\int_{0}^{x} K(x-s) y(s)\,{\rm d}s = 1 $$ with kernel $$ K(x-s) = 1 - 4 \sum_{n=1}^{\infty} \dfrac{1}{\lambda_n^2} e^{-\beta \lambda_n^2 (x-s)} $$ where $y(0)=1$, $\beta>0$ and $\lambda_n$ is the $n$-th positive root of the equation $J_0(x)=0$ (here $n$ is a natural number that numbers these positive roots in the order of increasing their values), $J_0(x)$ is the Bessel function of the first kind of zero order. I...
Are there any good visualization tutorials, written or video, that show graphically how separation of variables works? I particularly have the time-independent Schrodinger Equation in mind. There are hundreds of demonstrations out there which essentially distill to copies of one another. However I am trying to visualize in my mind how this process looks graphically - for example plotting t on one axis and x on the other for f(x,t). I have seen other good visual representations of...
Back
Top