# Navier code

**URL:** <https://discourse.julialang.org/t/navier-code/86898>\
**Category:** New to Julia\
**Created:** [September 7, 2022, 3:48pm UTC](https://discourse.julialang.org/t/navier-code/86898 "2022-09-07T15:48:51Z")\
**Posts on this page:** 6\
**Page:** 1

<div class="post-metadata">

**Author:** ![simonluzara](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/simonluzara/32/205260_2.png) [@simonluzara](https://discourse.julialang.org/u/simonluzara)\
**Post date:** [September 7, 2022, 3:48pm UTC](https://discourse.julialang.org/t/navier-code/86898/1 "2022-09-07T15:48:51Z")

</div>

the site suggested I should send you a pm. Hope it’s fine with you.  
Why did you say the code it wrong? Do you have any comment looking at it?

It is a simple annular flow with harmonic compressible boundary conditions (a fundamental Hydraulic Damper).

---

<div class="post-metadata">

**Author:** ![Oscar\_Smith](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/oscar_smith/32/25343_2.png) [@Oscar\_Smith](https://discourse.julialang.org/u/Oscar_Smith)\
**Post date:** [September 7, 2022, 3:54pm UTC](https://discourse.julialang.org/t/navier-code/86898/2 "2022-09-07T15:54:27Z")

</div>

yeah. rad is a scalar, so you aren’t integrating the thing you meant to.

---

<div class="post-metadata">

**Author:** ![simonluzara](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/simonluzara/32/205260_2.png) [@simonluzara](https://discourse.julialang.org/u/simonluzara)\
**Post date:** [September 7, 2022, 3:56pm UTC](https://discourse.julialang.org/t/navier-code/86898/3 "2022-09-07T15:56:01Z")

</div>

yeah but velocity is a vector. its like doing a r1u1+r2u2…to get the volume flow rate.

---

<div class="post-metadata">

**Author:** ![Oscar\_Smith](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/oscar_smith/32/25343_2.png) [@Oscar\_Smith](https://discourse.julialang.org/u/Oscar_Smith)\
**Post date:** [September 7, 2022, 3:56pm UTC](https://discourse.julialang.org/t/navier-code/86898/4 "2022-09-07T15:56:49Z")

</div>

but rad is a scalar, so it’s rad_u1+rad_u2…

---

<div class="post-metadata">

**Author:** ![simonluzara](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/simonluzara/32/205260_2.png) [@simonluzara](https://discourse.julialang.org/u/simonluzara)\
**Post date:** [September 7, 2022, 3:58pm UTC](https://discourse.julialang.org/t/navier-code/86898/5 "2022-09-07T15:58:19Z")

</div>

rad is a vector, I have defined it above as range(0.01,0.017,length=N) from r1 to r1+delta with 101 point divisions making it 100 numerical spaces.

---

<div class="post-metadata">

**Author:** ![Oscar\_Smith](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/oscar_smith/32/25343_2.png) [@Oscar\_Smith](https://discourse.julialang.org/u/Oscar_Smith)\
**Post date:** [September 7, 2022, 3:59pm UTC](https://discourse.julialang.org/t/navier-code/86898/6 "2022-09-07T15:59:16Z")

</div>

sorry, I missed that.