Dirichlet 3D–0D genBC Coupling Diverges in CFD

Description

When the 3D model is coupled to a 0D genBC using a Dirichlet boundary condition in a CFD simulation, the simulation diverges immediately.

Here is the usage I have in xml file:

   <Add_BC name="lumen_inlet" >
      <Type> Dir </Type>
      <Time_dependence> Coupled </Time_dependence>
      <Profile> Parabolic </Profile>
      <Impose_flux> true </Impose_flux>
      <Effective_direction> (0, 0, 1) </Effective_direction>
   </Add_BC>

Simulation diverges immediately:

---------------------------------------------------------------------
 Eq     N-i     T       dB  Ri/R1   Ri/R0    R/Ri     lsIt   dB  %t
---------------------------------------------------------------------
 NS 1-1  5.100e-01  [0 1.000e+00 1.000e+00 5.551e-04]  [4 -16 55]
 NS 1-2  1.017e+00  [1 1.139e+00 1.139e+00 2.022e-04]  [5 -17 62]
 NS 1-3  1.467e+00  [16 6.686e+00 6.686e+00 8.851e-04]  [4 -18 58]
 NS 1-4  1.920e+00  !41 1.163e+02 1.163e+02 7.885e-04!  [4 -17 58]
 NS 1-5  2.737e+00  !58 8.396e+02 8.396e+02 6.908e-04!  [8 -7 76]
 Error! NAN encountered..Dir
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_OVERFLOW_FLAG
 NS 1-6  3.903e+00  !45 1.984e+02 1.984e+02 4.206e-03!  !9 -3 83!  WARNING: The linear system solution has not converged
 Error! NAN encountered..Dir
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_OVERFLOW_FLAG
 NS 1-7  5.023e+00  !33 4.891e+01 4.891e+01 2.628e-03!  !9 -4 83!  WARNING: The linear system solution has not converged
 Error! NAN encountered..Dir
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_OVERFLOW_FLAG
 NS 1-8  6.217e+00  [19 9.818e+00 9.818e+00 3.852e-03]  !9 -3 84!  WARNING: The linear system solution has not converged
 Error! NAN encountered..Dir
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_OVERFLOW_FLAG
 NS 1-9  7.390e+00  [8 2.556e+00 2.556e+00 3.756e-03]  !9 -3 84!  WARNING: The linear system solution has not converged
 Error! NAN encountered..Dir
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_OVERFLOW_FLAG

I’m going to set up a simple test case with a prescribed flow rate in genBC and check whether this issue also occurs in svFSI.

Reproduction

Attached is a simulation folder that couples the inlet and outlet to a closed loop lumped network.

Conduit_genBC.zip

Expected behavior

The simulation should converge and produce reasonable results.

Additional context

No response

Code of Conduct

• I agree to follow this project's Code of Conduct and Contributing Guidelines