Skip to content

Fix input_t.prm in the nonlinear channel flow benchmark#6851

Merged
alarshi merged 1 commit intogeodynamics:mainfrom
YiminJin:input_t_prm
Feb 10, 2026
Merged

Fix input_t.prm in the nonlinear channel flow benchmark#6851
alarshi merged 1 commit intogeodynamics:mainfrom
YiminJin:input_t_prm

Conversation

@YiminJin
Copy link
Contributor

@YiminJin YiminJin commented Feb 7, 2026

This PR is a follow-up of the latest issue displayed in #6471 . The template input file input_t.prm in benchmarks/newton_solver_benchmark_set/nonlinear_channel_flow uses include input_v.prm for simplicity, but it makes the bash script run.sh produce the same result for different cases, because commands like -e "s/set Output directory .*/set Output directory = results\/$dirname_clean/g" does not work. This PR replaces include input_v.prm by the actual parameter settings, and the problem is mended. The results with traction boundary conditions are show below:
figure_t

Copy link
Contributor

@alarshi alarshi left a comment

Choose a reason for hiding this comment

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

Thank you for noticing this bug and fixing it!
What you did is correct but I was wondering if we can achieve the same changes by just copying the Solver parameters subsection (ln 19- ln 35) and including the input_v.prm file as before, since everything else is not changed by the run.sh script.
Besides one small change in the comment, it looks good to me if the tester is happy.

# Like input_v.prm and based on the nonlinear channel
# flow benchmark. This is used to test the traction boundary conditions
# Like the input_v.prm test and based on the nonlinear channel
# flow benchmark. This is used to test the velocity boundary conditions
Copy link
Contributor

Choose a reason for hiding this comment

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

I think you accidentally changed it. Please change it back to traction.

@alarshi
Copy link
Contributor

alarshi commented Feb 8, 2026

/rebuild

Copy link
Member

@MFraters MFraters left a comment

Choose a reason for hiding this comment

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

Looks good to me, besides the comment from @alarshi and a small comment I have.

set Variable names = x,y

# For velocity boundary conditions both are used, for pressure boundary conditions only the first (x) component
set Function expression = 0;(1e-37/(n+1))*((1e9/8e3)^n)*(((5e3)^(n+1))-(abs(x-(5e3))^(n+1)));
Copy link
Member

Choose a reason for hiding this comment

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

I think the second part can be zero as well now as it is no y component. It should not matter, but I think it looks more clean than a complicated function.

@YiminJin
Copy link
Contributor Author

Thanks @alarshi @MFraters ! I think your advice is smart and I have followed it. It is true that with fewer modifications there is lower possibility to make mistakes, and the total size of the software becomes smaller.

Copy link
Contributor

@alarshi alarshi left a comment

Choose a reason for hiding this comment

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

Looks good, thank you!

@alarshi alarshi merged commit ae0bab6 into geodynamics:main Feb 10, 2026
8 of 9 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants