Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Stokes problem description #343

Open
dgm2 opened this issue Sep 12, 2024 · 1 comment
Open

Stokes problem description #343

dgm2 opened this issue Sep 12, 2024 · 1 comment
Labels
help wanted Extra attention is needed question Further information is requested

Comments

@dgm2
Copy link

dgm2 commented Sep 12, 2024

Hello
could you add a brief description of the Stoke problem similar to the other tutorials?
https://github.com/mathLab/PINA/blob/fefba810341379b51b9a49e59d45fd82d668b1bb/examples/problems/stokes.py
referred to the run_stokes file

Best regards

@dgm2 dgm2 added the help wanted Extra attention is needed label Sep 12, 2024
@dario-coscia dario-coscia added the question Further information is requested label Sep 12, 2024
@dario-coscia
Copy link
Collaborator

Hi @dgm2 ! Unfortunately on examples we do not put the mathematical formulations as you noticed. The problem formulation is the following:

$$ \begin{cases} -\Delta \mathbf{u}(x,y) + \nabla p(x,y) = 0\quad (x,y)\in\Omega \quad \texttt{momentum} \\ \nabla\cdot\mathbf{u}(x,y)= 0\quad (x,y)\in\Omega \quad \texttt{continuity} \\ u_x(x,y) - 2(1-y^2) = 0 \quad (x,y)\in\Gamma_{in} \quad \texttt{inlet}\\ p(x,y) = 0 \quad (x,y)\in\Gamma_{out} \quad \texttt{outlet}\\ u(x,y) = 0 \quad (x,y)\in\Gamma_{top}\cup\Gamma_{bot} \quad \texttt{wall}\\ \end{cases} $$

Where $u(x,y)=[u_x(x,y), u_y(x,y)]$ is a vectorial function, and $p$ is a scalar function. The domain is a square domain $\Omega = [-2,2]\times[-1,1]$, with $\Gamma_{in}={-2}\times[-1,1]$, $\Gamma_{out}={2}\times[-1,1]$, $\Gamma_{top}=[-2,2] \times 1$, $\Gamma_{bot}=[-2,2] \times -1$.

I hope this will help😄

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
help wanted Extra attention is needed question Further information is requested
Projects
None yet
Development

No branches or pull requests

2 participants