Skip to content

Heat transport process added#1328

Closed
grubbymoon wants to merge 22 commits intoufz:masterfrom
grubbymoon:HeatTransport_Pull
Closed

Heat transport process added#1328
grubbymoon wants to merge 22 commits intoufz:masterfrom
grubbymoon:HeatTransport_Pull

Conversation

@grubbymoon
Copy link
Copy Markdown
Contributor

validated with 1D heat conduction analytical solution with dirichlet boundary

@grubbymoon grubbymoon force-pushed the HeatTransport_Pull branch 2 times, most recently from cd2b89c to 053c3c1 Compare August 1, 2016 15:53
@ogsbot
Copy link
Copy Markdown
Member

ogsbot commented Aug 1, 2016

Jenkins: OGS-6/Linux-PRs failed: https://svn.ufz.de:8443/job/OGS-6/job/Linux-PRs/2442/

Comment thread ProcessLib/CMakeLists.txt Outdated
endif()
if(TARGET Boost)
add_dependencies(ProcessLib Boost)
add_dependencies(ProcessLib Boost)
Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

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

If possible, please keep the original formatting of the lines you don't change.

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

Thanks, corrected

@endJunction
Copy link
Copy Markdown
Member

I would not call this process a "heat transport" process, which, for me, implies non-linear advective transport; better to call it "linear heat diffusion", or "linear heat equation", or "conductive heat transfer", or "heat conduction".

Comment thread Applications/CLI/Tests.cmake Outdated
temperature_analytical.vtu line_60_heat_pcs_0_ts_81_t_31640625.000000.vtu Temperature_Analytical_1year temperature
)

AddTest(
Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

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

Indentations in this AddTest needs four spaces.

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

corrected

@ogsbot
Copy link
Copy Markdown
Member

ogsbot commented Aug 2, 2016

Jenkins: OGS-6/Mac-PRs failed: https://svn.ufz.de:8443/job/OGS-6/job/Mac-PRs/2310/

@ogsbot
Copy link
Copy Markdown
Member

ogsbot commented Aug 2, 2016

Jenkins: OGS-6/Linux-PRs-dynamic failed: https://svn.ufz.de:8443/job/OGS-6/job/Linux-PRs-dynamic/668/

@ogsbot
Copy link
Copy Markdown
Member

ogsbot commented Aug 2, 2016

Jenkins: OGS-6/Linux-PRs failed: https://svn.ufz.de:8443/job/OGS-6/job/Linux-PRs/2445/

@ogsbot
Copy link
Copy Markdown
Member

ogsbot commented Aug 2, 2016

Jenkins: OGS-6/Gui/Gui-Linux-PRs failed: https://svn.ufz.de:8443/job/OGS-6/job/Gui/job/Gui-Linux-PRs/2392/

@chleh
Copy link
Copy Markdown
Collaborator

chleh commented Aug 2, 2016

Jenkins, test this please.

@ogsbot
Copy link
Copy Markdown
Member

ogsbot commented Aug 2, 2016

Jenkins: OGS-6/Linux-PRs failed: https://svn.ufz.de:8443/job/OGS-6/job/Linux-PRs/2446/

@ogsbot
Copy link
Copy Markdown
Member

ogsbot commented Aug 2, 2016

Jenkins: OGS-6/Linux-PRs-dynamic failed: https://svn.ufz.de:8443/job/OGS-6/job/Linux-PRs-dynamic/669/

@ogsbot
Copy link
Copy Markdown
Member

ogsbot commented Aug 2, 2016

Jenkins: OGS-6/Gui/Gui-Linux-PRs failed: https://svn.ufz.de:8443/job/OGS-6/job/Gui/job/Gui-Linux-PRs/2394/

@wenqing
Copy link
Copy Markdown
Member

wenqing commented Aug 2, 2016

@endJunction As for the name of process, I think HEAT_TRANSPORT is adequate because that we only need one object to handle the simulation of heat transport process, the PDE of which could be linear or non-linear, or which is coupled to other physical processes. The advection term, which is a coupling term to fluid process, is not calculated in this PR so far.

@wenqing
Copy link
Copy Markdown
Member

wenqing commented Aug 2, 2016

@endJunction The advection term could be added in future in the implementation of coupled TH process.

Comment thread Applications/CLI/Tests.cmake Outdated
EXECUTABLE_ARGS line_60_heat.prj
WRAPPER time
TESTER vtkdiff
ABSTOL 1e-2 RELTOL 1e-2
Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

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

I'd say that in absence of any singularities in the solution, the tolerances are a little to sloppy. Why is it so? Can we improve it?

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

i have improved the tolerances from 1e-2 to 1e-4. The error between analytical solution and numerical solution is normally a little bit sloppy compared with the error between two numerical solutions

@ogsbot
Copy link
Copy Markdown
Member

ogsbot commented Aug 4, 2016

Jenkins: OGS-6/Linux-PRs-dynamic failed: https://svn.ufz.de:8443/job/OGS-6/job/Linux-PRs-dynamic/685/

@grubbymoon
Copy link
Copy Markdown
Contributor Author

@endJunction @wenqing @chleh Please review again !

@@ -0,0 +1,380 @@
/**
Copy link
Copy Markdown
Collaborator

@chleh chleh Sep 1, 2016

Choose a reason for hiding this comment

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

This file an artefact of your rebasing. Please remove it.
And please also remove the file .merge_file_R8GkSe

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

removed them.

@ogsbot
Copy link
Copy Markdown
Member

ogsbot commented Sep 1, 2016

Jenkins: OGS-6/Linux-PRs failed: https://svn.ufz.de:8443/job/OGS-6/job/Linux-PRs/2781/

{
class HeatConductionProcess final : public Process
{
using Base = Process;
Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

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

This line is probably unneeded. Please remove it.

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

corrected

@chleh
Copy link
Copy Markdown
Collaborator

chleh commented Sep 1, 2016

Only small things from my side. Afterwards and after data cleanup 👍

@ogsbot
Copy link
Copy Markdown
Member

ogsbot commented Sep 1, 2016

Jenkins: OGS-6/Linux-PRs failed: https://svn.ufz.de:8443/job/OGS-6/job/Linux-PRs/2782/

_localM.noalias() += sm.N.transpose() * density * heat_capacity *
sm.N * sm.detJ * wp.getWeight();
// heat flux only computed for output.
const NodalVectorType heat_flux =
Copy link
Copy Markdown
Member

@endJunction endJunction Sep 2, 2016

Choose a reason for hiding this comment

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

Needs checking of the correct size. And remove extra parentheses.

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

corrected

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

@endJunction I have tested the error threshold. 1e-4 is required for neumann boundary and 1e-5 is fine with dirichlet boundary ! I have updated the threshold.

endJunction added a commit that referenced this pull request Sep 4, 2016
@endJunction
Copy link
Copy Markdown
Member

Merged in master

Had to update Tests/Data due to unrelated files. Also squashed the fixes for clarity.

@endJunction endJunction closed this Sep 4, 2016
@ogsbot
Copy link
Copy Markdown
Member

ogsbot commented Jun 19, 2020

OpenGeoSys development has been moved to GitLab.

See this pull request on GitLab.

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.

5 participants