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

Correct the Poisson equation in the cavity flow example #2308

Merged
merged 1 commit into from
Feb 13, 2024

Conversation

rafael-fuente
Copy link
Contributor

There is a sign mistake in the Poisson equation. It doesn't significantly affect the solution of this particular example, but it can be important when solving Navier Stokes for other problems.

Copy link

Check out this pull request on  ReviewNB

See visual diffs & provide feedback on Jupyter Notebooks.


Powered by ReviewNB

@mloubout mloubout added the examples examples label Feb 12, 2024
@mloubout
Copy link
Contributor

Good catch, thanks!

Copy link

codecov bot commented Feb 12, 2024

Codecov Report

All modified and coverable lines are covered by tests ✅

Comparison is base (8c86bc2) 86.71% compared to head (f68c098) 86.71%.

Additional details and impacted files
@@           Coverage Diff           @@
##           master    #2308   +/-   ##
=======================================
  Coverage   86.71%   86.71%           
=======================================
  Files         229      229           
  Lines       42958    42958           
  Branches     7958     7958           
=======================================
+ Hits        37251    37253    +2     
+ Misses       5023     5021    -2     
  Partials      684      684           

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

@@ -474,7 +474,7 @@
"# First order derivatives will be handled with p.dxc\n",
"eq_u =Eq(u.dt + u*u.dx + v*u.dy, -1./rho * p.dxc + nu*(u.laplace), subdomain=grid.interior)\n",
"eq_v =Eq(v.dt + u*v.dx + v*v.dy, -1./rho * p.dyc + nu*(v.laplace), subdomain=grid.interior)\n",
"eq_p =Eq(p.laplace,rho*(1./dt*(u.dxc+v.dyc)-(u.dxc*u.dxc)+2*(u.dyc*v.dxc)+(v.dyc*v.dyc)), subdomain=grid.interior)\n",
"eq_p =Eq(p.laplace,rho*(1./dt*(u.dxc+v.dyc)-(u.dxc*u.dxc)-2*(u.dyc*v.dxc)-(v.dyc*v.dyc)), subdomain=grid.interior)\n",
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

There's also a missing space around the equals here right (unless it's rendering weird on my laptop)? Would it be possible to tag that onto the PR?

@mloubout mloubout merged commit bb7491d into devitocodes:master Feb 13, 2024
31 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
examples examples
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants