Компонент стандартной библиотеки Modelica AbruptAdaptor

Question

Протестировано с:

Версии MSL:

• 3.2.1, 2013-08-14, сборка 2 (2013-08-14 08: 44: 41Z)
• 3.2.1, 2013-08-14, сборка 4 (2015-09-30 09:15:00Z)

Ниже представлена ​​модель с 4 компонентами:

1. Источник массового расхода (M_bound)
2. AbruptAdapter (abruptAdapter)
3. Динамическая труба (pipeSmall)
4. Источник давления (P_bound)

Два случая подключения модели также содержатся в equation раздел:

Случай № 1: M_bound -> abruptAdaptor -> pipeSmall -> P_bound

Случай № 2: M_bound -> pipeSmall -> abruptAdaptor -> P_bound

Резюме: Случай № 1 имитирует нормально, но Случай № 2 дает сбой. Сгенерированная ошибка показана ниже:

The following error was detected at time: 0

Model error - division by zero: (data.zeta1) / ((if data.zeta1_at_a then rho_a_des*A_a^2 else rho_b_des*A_b^2)) = (0.402964) / (0)

The stack of functions is:
Modelica.Fluid.Fittings.BaseClasses.QuadraticTurbulent.pressureLoss_m_flow_totalPressure
Modelica.Fluid.Fittings.BaseClasses.QuadraticTurbulent.pressureLoss_m_flow_totalPressure(
 -pipeBig.port_b.m_flow, 
abruptAdaptor.state_a.d, 
abruptAdaptor.state_b_des.d, 
abruptAdaptor.state_b.d, 
abruptAdaptor.state_a_nondes.d, 
abruptAdaptor.data, 
abruptAdaptor.m_flow_small)
First evaluation failed for non-linear solver.

Вопрос: Почему произошел сбой в случае № 2 и рекомендации по его устранению? Любое понимание было бы здорово. Это может быть даже ошибка в AbruptAdaptor составная часть.

Код модели

model AbruptAdaptor_1Pipe

  replaceable package Medium = Modelica.Media.Water.StandardWater;

  inner Modelica.Fluid.System system;

  Modelica.Fluid.Sources.Boundary_pT P_bound(
    redeclare package Medium = Medium,
    nPorts=1,
    p=system.p_ambient,
    T=M_bound.T);

  Modelica.Fluid.Sources.MassFlowSource_T M_bound(
    redeclare package Medium = Medium,
    m_flow=1,
    T=300,
    nPorts=1);

  Modelica.Fluid.Fittings.AbruptAdaptor abruptAdaptor(
    redeclare package Medium = Medium,
    diameter_a=1,
    diameter_b=0.5,
    m_flow_nominal=M_bound.m_flow,
    m_flow_start=M_bound.m_flow);

  Modelica.Fluid.Pipes.DynamicPipe pipeSmall(
    redeclare package Medium = Medium,
    T_start=M_bound.T,
    m_flow_start=M_bound.m_flow,
    length=1,
    diameter=0.5,
    modelStructure=Modelica.Fluid.Types.ModelStructure.av_b); // switch to av_vb with case #2

equation 

  // Case #1
  connect(M_bound.ports[1],abruptAdaptor.port_a);
  connect(abruptAdaptor.port_b,pipeSmall.port_a);
  connect(pipeSmall.port_b,P_bound.ports[1]);

  // Case #2
  //connect(M_bound.ports[1],pipeSmall.port_b);
  //connect(abruptAdaptor.port_b,pipeSmall.port_a);
  //connect(abruptAdaptor.port_a,P_bound.ports[1]);

  annotation (uses(Modelica(version="3.2.1")));
end AbruptAdaptor_1Pipe;

Answer 1

Ваша смещенная сетка в порядке, вы устанавливаете объемы там, где должны были. Не волнуйтесь, это не ваша ошибка, это отладка, я пытался использовать ее несколько раз, и она всегда дает сбой одним из способов (я не знаю положение вашего адаптера, но это не имеет значения). Вы всегда можете использовать отверстия, где вы определяете отношения dp(Areas_rate), избегая потоков реверсирования в обоих случаях, поэтому соединяйте два объема до и после отверстий (например, два клапана параллельно). Это должно работать.