lanl / lanl/THMC_verification_scripts

Question regarding coefficient of consolidation Cv values based on Table 1 parameter values.

Open
#1 2 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

Dominant language
Jupyter Notebook
Stars
2
Forks
0
PR merge metrics
No merged PRs in 30d

Description

Greetings,

I'm running your Jupyter notebook Terazaghi example, attempting to reproduce your results shown in Fig. 2 of your paper. With t set to 10 hours, pore pressure reaches 1e-19 Pa with the parameters given in the notebook. I changed P0 to 10e6 (10MPa), as specified in Table 1, which results in pore pressure reaching 1e-20 Pa. As shown in Fig. 2, pore pressure should be 8 MPa at 5m at 10 hr. This will result when the coefficient of consolidation Cv is O(2e-4) with no modification of P0. However, your Cv is 0.3941 [m^2/s] based on eq. 18 and the given values in Table 1. Can you shed some light on this discrepancy? What values of Cv yield the results shown in Fig. 2, and from what parameter values are these Cv values derived?

Thank you so much.

Chris

Contributor guide

No contributing guide indexed for this repository

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

Research direction

Start with the Terzaghi example Jupyter notebook and reproduce the reported run at t=10 hours and 5 m using the notebook values and P0=10 MPa. Compare the computed coefficient of consolidation from equation 18 and Table 1 with the Cv value needed to match Fig. 2, then document which parameter values account for the discrepancy.

Written by the indexing model from the issue text.

Assessment

Tech stack
jupyter-notebook
Domain
testing-qa
Issue type
Bug
Difficulty
4/5
Estimated time
3-5 days
Activity status
Stale
Clarity
Mostly clear
Newbie friendliness
35/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.