modelElasticPenstock
Extends from Modelica.Icons.Example (Icon for runnable examples).
Information
Here is example of using the KP function to solve hyperbolic PDE (here, model for penstock with compressible water and elastic walls is used).
All calculation of the variables that is used for defining eigenvalues, source term S and vector F are implemented in additional function BasicEquation which used one time for centered values and then for boundary values.
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| Integer | N | 20 | |
| Geometry | |||
| SI.Height | H | 420 | Height over which water fall in the pipe |
| SI.Length | L | 600 | Length of the pipe |
| SI.Diameter | D | 3.3 | Diameter from the input side of the pipe |
| Initialization | |||
| SI.VolumeFlowRate | Vdot_0 | 20 | Initial flow rate in the pipe, m3/s |
Components
| Type | Name | Default | Description |
|---|---|---|---|
| Data | data | ||
| SI.Area | A_atm | D^2*pi/4 | |
| SI.Pressure[N] | p_p | ||
| SI.Pressure | dp | data.rho*data.g*H/N | |
| SI.Pressure | p_1 | 8e5 | |
| SI.Pressure | p_2 | 48e5 | |
| SI.Pressure[N,4] | p_ | ||
| SI.Length | dx | L/N | |
| SI.Length[N + 4] | B | zeros(N + 4) | |
| SI.MassFlowRate[N] | mdot | ||
| SI.MassFlowRate[N,4] | mdot_ | ||
| SI.MassFlowRate | mdot_R | Vdot_0*data.rho | |
| SI.MassFlowRate | mdot_V | Vdot_0*data.rho | |
| Real[8,N] | U_ | ||
| Real[2*N] | S_ | ||
| Real[2*N,4] | F_ | ||
| Real[N,4] | lam1 | ||
| Real[N,4] | lam2 | ||
| SI.VolumeFlowRate[N] | Vdot | ||
| Real | theta | 1.3 | |
| Real[2*N] | U | ||
| Real[N] | F_d | ||
| BasicEquation | basic | ||
| Functions.KP07.KPmethod | kP | ||
| BasicEquation | basicMid |