Comments (8)
cbuahin
commented on July 21, 2024
@brycemihal and @btneilson, I have added the friction heat flux to the code(see 2e064c4). You turn it on by setting the new flag called FLUID_FRICTION_HEAT to YES in the OPTIONS section. Please pull it, compile, and replace the CSHComponent in your existing installation. Remember to compile in release mode. Let me know if you have any trouble compiling.
from cshcomponent.
brycemihal
commented on July 21, 2024
Thanks Caleb for working on this. I was able to pull and compile the new
CSH component but when I turn the FLUID_FRICTION_HEAT option to YES the
model crashes. When I have it set to NO, it runs fine. I'm not comfortable
going through QT to debug the code and find the problem. Perhaps this is
something you could walk me through sometime.
The model crashes on the first CSH time step and I get the output of:
…
[CVODE WARNING] CVode
Internal t = 0 and h = 0 are such that t + h = t on the next step. The
solver will continue anyway.
[CVODE ERROR] CVode
At t = 0 and h = 0, the corrector convergence test failed repeatedly or
with |h| = hmin.
CSH Solver failed
On Sat, Jun 8, 2019 at 8:47 PM Caleb Buahin ***@***.***> wrote:
@brycemihal <https://github.com/brycemihal> and @btneilson
<https://github.com/btneilson>, I have added the friction heat flux to
the code(see 2e064c4
<2e064c4>).
You turn it on by setting the new flag called FLUID_FRICTION_HEAT to YES in
the OPTIONS section. Please pull it, compile, and replace the CSHComponent
in your existing installation. Remember to compile in release mode. Let me
know if you have any trouble compiling.
—
You are receiving this because you were mentioned.
Reply to this email directly, view it on GitHub
<#1?email_source=notifications&email_token=ACX5IHPF4PCOI6FITZNENN3PZRVNVA5CNFSM4HVZ4VV2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGODXICK7A#issuecomment-500180348>,
or mute the thread
<https://github.com/notifications/unsubscribe-auth/ACX5IHNMHZFVT66B4DFR7GTPZRVNVANCNFSM4HVZ4VVQ>
.
from cshcomponent.
cbuahin
commented on July 21, 2024
@brycemihal Sorry about that. I think I forgot to initialize the friction heat variable to zero. I am going to fix it shortly and recommit it.
from cshcomponent.
cbuahin
commented on July 21, 2024
I actually checked and I am initializing things correctly. Check the friction equation in element.cpp in the function Element::computeFluidFrictionHeat() to make sure the formulation is right.
from cshcomponent.
brycemihal
commented on July 21, 2024
Caleb, I did find an error in your equation: *fluidFrictionHeatFlux =
9805.0 * flow.value / (width * slope)*;. Slope should not be in the
denominator. See the attached screenshot from Theurer et al. 1985. I don't
think the formulation is throwing the error though. I looked into my
element widths in previous model runs and every element is zero for the
first time step. This would cause a NaN value when computed. It looks like
I have the initial XSECTION_AREA, WIDTH, SLOPE, FLOW, and DISPERSION_COEFF
all set to zero in my CSH .inp file. I'm guessing that if I set my WIDTH
value to something other than zero the problem will be fixed. I will test
this and let you know if it works.
[image: Screen Shot 2019-06-10 at 1.56.33 PM.png]
…On Mon, Jun 10, 2019 at 1:47 PM Caleb Buahin ***@***.***> wrote:
I actually checked and I am initializing things correctly. Check the
friction equation in element.cpp in the function
Element::computeFluidFrictionHeat() to make sure the formulation is right.
—
You are receiving this because you were mentioned.
Reply to this email directly, view it on GitHub
<#1?email_source=notifications&email_token=ACX5IHLUSUTJBA4POZOTAYLPZ2VT5A5CNFSM4HVZ4VV2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGODXLA2WA#issuecomment-500567384>,
or mute the thread
<https://github.com/notifications/unsubscribe-auth/ACX5IHOATPNATCK2C3TUFGTPZ2VT5ANCNFSM4HVZ4VVQ>
.
from cshcomponent.
cbuahin
commented on July 21, 2024
@brycemihal that sounds good. Also, you can modify the code to remove my check for non zero slope since the slope is in the numerator. Please commit any changes you make. In the message just type #1 to reference this issue.
from cshcomponent.
cbuahin
commented on July 21, 2024
@brycemihal I have actually fixed it here so nvm about committing your changes.
from cshcomponent.
brycemihal
commented on July 21, 2024
Caleb, I just committed my changes but under the Bryce branch. The model
appears to be working now with the change in initial width.
…On Mon, Jun 10, 2019 at 2:21 PM Caleb Buahin ***@***.***> wrote:
@brycemihal <https://github.com/brycemihal> I have actually fixed it here
so nvm about committing your changes.
—
You are receiving this because you were mentioned.
Reply to this email directly, view it on GitHub
<#1?email_source=notifications&email_token=ACX5IHLSSAMBXBBZQ5TEPE3PZ2ZUNA5CNFSM4HVZ4VV2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGODXLDXWI#issuecomment-500579289>,
or mute the thread
<https://github.com/notifications/unsubscribe-auth/ACX5IHNHXFIH6OO6UP2TJVDPZ2ZUNANCNFSM4HVZ4VVQ>
.
from cshcomponent.
Related Issues (2)
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from cshcomponent.