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

Pseudo-Steady Solver Logic #50

Merged
merged 6 commits into from
Jul 29, 2021
Merged

Pseudo-Steady Solver Logic #50

merged 6 commits into from
Jul 29, 2021

Conversation

eyoung55
Copy link
Collaborator

@eyoung55 eyoung55 commented Jul 8, 2021

This PR allows the user to specify final_time = none in the yaml to trigger a pseudo-steady run of the unsteady solver, where the solver will continue iterating until the average velocity field converges.

  • Add stabilization coefficient to the list of default parameters (allows stabilization to be turned off with stability_eps: 0)
  • Adds mixing length model to unsteady solver for comparison to steady
  • Adds new logic for computing average velocity and testing convergence

It most likely won't pass GHA because it also uses the fast/optimized solver+preconditioner pairs which leads to slightly different results.

@eyoung55 eyoung55 changed the base branch from master to dev July 8, 2021 17:48
Comment on lines +1000 to +1003
# FIXME: We should be consistent about which way the turbine
# forces are oriented.
if self.farm.turbine_method != "alm":
self.tf *= -1.0
Copy link
Collaborator Author

Choose a reason for hiding this comment

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

@jefalon what do you think about standardizing the direction the turbine force points? I'd like to swap everything so that turbine forces point opposite the incoming flow, but recognize that ALM is the odd one out and it would be easier to bring it into line with the others.

Copy link
Collaborator

Choose a reason for hiding this comment

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

I think that's a good idea and having it opposite to flow is logical since it's suppose to be forcing it to slow down. Thinks to check to make sure this sign flip doesn't break anything is:

  • the sign of the tf term in the system functional for both the stabilized and Taylor hood steady solver
  • the sign on the power and 2d_power objective functions
  • probably some other place I can't think of

Copy link
Collaborator

Choose a reason for hiding this comment

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

We need to fix the opt output for plotting

@jefalon jefalon merged commit 4d974cc into NREL:dev Jul 29, 2021
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

2 participants