Share an example only consists of one compressor and one motor written by me

[Wang-Yihu]

From the MOOSE code IdealGasFluidProperties.C, ADShaftConnectedCompressor1PhaseUserObject.C, the equations for compressor are:
(only for the steady)

############################################
#Eq 1: s_1 = cv * ln( T_1^γ / p_1^(γ-1) )
#Eq 2: p_1 * v_1 = R/M * T_1
#Eq 3: u_1 = m_dot * v_1 / A_1
#Eq 4: h0_1 = cv * T_1 + p_1 * v_1 + 0.5 * u_1^2
#Eq 5: p0_1 = (h0_1 / γ / cv)^(γ / (γ-1)) * exp(-s_1 / (γ-1) /cv)
#Eq 6: v0_1 = R / M / p0_1 * exp( (s_1+cv * (γ-1) * lnp0_1) / γ / cv )
#Eq 7: T0_1 = p0_1 * v0_1 * M / R
#Eq 8: c0_1 = sqrt(γ * R * T0_1 / M)
#Eq 9: p0_2 = p0_1 * Rp( m_dot * v0_1 / c0_1 / (rated), ω / c0_1 / (rated) )
#Eq 10: v0_2_is = R / M / p0_2 * exp( (s_1 + cv * (γ-1) * lnp0_2) / γ / cv )
#Eq 11: T0_2_is = p0_2 * v0_2_is * M / R
#Eq 12: h0_2_is = cv * T0_2_is + p0_2 * v0_2_is
#Eq 13: h0_2 = (h0_2_is - h0_1 +η * h0_1) / η
#Eq 14: τ_c = -m_dot / ω * (h0_2 - h0_1)
#Eq 15: (u_1^2 / v_1 + p_1) * A_1 - (u_2^2 / v_2 + p_2) * A_2 + (p0_2 - p0_1) *A_ref = 0
#Eq 16: h0_2 = cv * T_2 + p_2 * v_2 + 0.5 * u_2^2
#Eq 17: p_2 * v_2 = R/M * T_2
#Eq 18: u_2 = m_dot * v_2/A_2
######################################################

There are 18 equations and 22 variables:

s_1, T_1, p_1, v_1, u_1
m_dot, h0_1, p0_1, v0_1, T0_1
c0_1, p0_2, v0_2_is, T0_2_is, h0_2_is
h0_2, τ_c, ω, p_2, u_2
T_2, v_2

In MOOSE, we always know p0_1, T0_1, τ, p_2
I write a python code to solve this equation group

#https://mooseframework.inl.gov/source/components/ShaftConnectedCompressor1Phase.html
import math
from scipy.optimize import fsolve
from scipy.interpolate import interp1d
############################################
#Eq 1: s_1 = cv*ln( T_1^γ/p_1^(γ-1) )
#Eq 2: p_1*v_1 = R/M*T_1
#Eq 3: u_1 = m_dot*v_1/A_1
#Eq 4: h0_1 = cv*T_1 + p_1*v_1 + 0.5*u_1^2
#Eq 5: p0_1 = (h0_1/γ/cv)^(γ/(γ-1)) * exp(-s_1/(γ-1)/cv)
#Eq 6: v0_1 = R/M/p0_1 * exp( (s_1+cv*(γ-1)*lnp0_1)/γ/cv )
#Eq 7: T0_1 = p0_1*v0_1*M/R
#Eq 8: c0_1 = sqrt(γ*R*T0_1/M)
#Eq 9: p0_2 = p0_1*Rp( m_dot*v0_1/c0_1/(rated), ω/c0_1/(rated) )
#Eq 10: v0_2_is = R/M/p0_2 * exp( (s_1+cv*(γ-1)*lnp0_2)/γ/cv )
#Eq 11: T0_2_is = p0_2*v0_2_is*M/R
#Eq 12: h0_2_is = cv*T0_2_is + p0_2*v0_2_is
#Eq 13: h0_2 = (h0_2_is - h0_1 +η*h0_1)/η
#Eq 14: τ_c = -m_dot/ω * (h0_2 - h0_1)
#Eq 15: (u_1^2/v_1 + p_1)*A_1 - (u_2^2/v_2 + p_2)*A_2 + (p0_2 - p0_1)*A_ref = 0
#Eq 16: h0_2 = cv*T_2 + p_2*v_2 + 0.5*u_2^2
#Eq 17: p_2*v_2 = R/M*T_2
#Eq 18: u_2 = m_dot*v_2/A_2
######################################################
def interpolate_2d_comp(x_input, y_input):
    # Table data: y: ([x values], [z values])
    # y for omega, x for mdot, z for Rp
    table = {
        0.5208: ([0.2868480, 0.3084640, 0.3300760, 0.3516920, 0.3733040, 0.3949200, 0.4165200, 0.4381600, 0.4597600, 0.4813600, 0.5030000, 0.5187200],
                 [1.5153000, 1.5088000, 1.5032000, 1.4984000, 1.4914000, 1.4863000, 1.4846000, 1.4756000, 1.4540000, 1.4237000, 1.3856000, 1.3506000]),

        0.6250: ([0.3674120, 0.3890240, 0.4106400, 0.4322400, 0.4523200, 0.4754800, 0.4970800, 0.5177200, 0.5442400, 0.5658800, 0.5874800, 0.6090800, 0.6307200, 0.6484000, 0.6621600, 0.6759200, 0.6886800, 0.6994800, 0.7093200, 0.7181600, 0.7253600, 0.7338800, 0.7407600, 0.7466400, 0.7517600, 0.7595600, 0.7633600],
                 [1.7741000, 1.7626000, 1.7550000, 1.7536000, 1.7497000, 1.7455000, 1.7376000, 1.7276000, 1.7026000, 1.6780000, 1.6405000, 1.5914000, 1.5233000, 1.4564000, 1.3890000, 1.3187000, 1.2318000, 1.1631000, 1.0851000, 1.0143000, 0.9425300, 0.8601900, 0.7844300, 0.7160400, 0.6450700, 0.5733400, 0.5341100]),

        0.7292: ([0.4067200, 0.4283200, 0.4499200, 0.4715600, 0.4931600, 0.5147600, 0.5364000, 0.5580000, 0.5796400, 0.6012400, 0.6228400, 0.6444800, 0.6660800, 0.6872000, 0.7093200, 0.7299200, 0.7476400, 0.7643200, 0.7800400, 0.7928000, 0.8036400, 0.8134400, 0.8222800, 0.8301600, 0.8380000, 0.8448800, 0.8507600, 0.8566800, 0.8635600, 0.8704400, 0.8782800],
                 [2.1535000, 2.1352000, 2.1251000, 2.1192000, 2.1135000, 2.1088000, 2.1071000, 2.1068000, 2.1043000, 2.0975000, 2.0877000, 2.0742000, 2.0539000, 2.0272000, 1.9813000, 1.9245000, 1.8643000, 1.7959000, 1.7227000, 1.6551000, 1.5896000, 1.5110000, 1.4403000, 1.3714000, 1.2969000, 1.2326000, 1.1655000, 1.0985000, 1.0177000, 0.9385600, 0.8578700]),

        0.8333: ([0.4872800, 0.5088800, 0.5304800, 0.5521200, 0.5737200, 0.5953600, 0.6169600, 0.6385600, 0.6602000, 0.6818000, 0.7034000, 0.7250400, 0.7466400, 0.7682400, 0.7898800, 0.8114800, 0.8330800, 0.8547200, 0.8753600, 0.8940000, 0.9097200, 0.9224800, 0.9342800, 0.9450800, 0.9568800],
                 [2.6634000, 2.6485000, 2.6366000, 2.6248000, 2.6184000, 2.6124000, 2.6110000, 2.6102000, 2.6110000, 2.6102000, 2.6110000, 2.6099000, 2.6060000, 2.5989000, 2.5874000, 2.5642000, 2.5309000, 2.4824000, 2.4202000, 2.3500000, 2.2832000, 2.2195000, 2.1437000, 2.0706000, 1.9947000]),

        0.9375: ([0.5658800, 0.5874800, 0.6090800, 0.6307200, 0.6523200, 0.6739200, 0.6955600, 0.7171600, 0.7387600, 0.7604000, 0.7820000, 0.8036400, 0.8252400, 0.8468400, 0.8684800, 0.8900800, 0.9116800, 0.9333200, 0.9549200, 0.9765200, 0.9981600, 1.0197600, 1.0413600, 1.0600400, 1.0757600, 1.0904800, 1.1032800, 1.1150800, 1.1248800, 1.1327600, 1.1406000, 1.1484800, 1.1563200, 1.1632000, 1.1690800, 1.1750000, 1.1808800, 1.1867600, 1.1916800, 1.1976000],
                 [3.4718000, 3.4428000, 3.4189000, 3.3970000, 3.3728000, 3.3455000, 3.3249000, 3.3106000, 3.2976000, 3.2886000, 3.2819000, 3.2799000, 3.2844000, 3.2892000, 3.2898000, 3.2892000, 3.2889000, 3.2830000, 3.2703000, 3.2551000, 3.2281000, 3.1873000, 3.1291000, 3.0660000, 2.9999000, 2.9281000, 2.8556000, 2.7823000, 2.7101000, 2.6451000, 2.5778000, 2.5099000, 2.4280000, 2.3562000, 2.2850000, 2.2128000, 2.1344000, 2.0580000, 1.9920000, 1.9198000]),
    }

    y_list = sorted(table.keys())
    z_y_interp = []

    for y in y_list:
        x_vals, z_vals = table[y]
        f = interp1d(x_vals, z_vals, kind='linear', fill_value='extrapolate')
        z_y_interp.append(f(x_input))

    f_y = interp1d(y_list, z_y_interp, kind='linear', fill_value='extrapolate')
    return float(f_y(y_input))

# Example usage
z = interpolate_2d_comp(0.3890240, 0.6250)  # Extrapolates if needed
print(f"Interpolated z = {z:.4f}")
######################################################
#Input parameters:
p0_1 = 129155; T0_1 = 319.56; p_2 = 360000; tau = 49.45
eff = 0.79
M = 0.029; R = 8.314; gamma = 1.4; cv = 1/(gamma-1)*R/M; cp = gamma*cv
A_1 = math.pi/4*0.15**2; A_2 = math.pi/4*0.15**2; A_ref = (A_1 + A_2)/2
c0_rated = 358.13321019260536
rho0_rated = 1.146881112
m_dot_rated_c = 2.9319973795574246
omega_rated = 10053
###########################################
#These variables can be got without assumption:
#Eq 7, 6, 5, 8
v0_1 = T0_1*R/p0_1/M; rho_0_1 = 1/v0_1
s_1 = (math.log(v0_1) - math.log(R/M/p0_1)) * gamma*cv - cv*(gamma-1)*math.log(p0_1)
h0_1 = gamma*cv*(  p0_1/math.exp(-s_1/(gamma-1)/cv)  )**((gamma-1)/gamma)
c0_1 = math.sqrt(gamma*R*T0_1/M)

def equation(vars):
    m_dot, omega, p_1 = vars
#Combine Eq 1, 4, 3, 2
    T_1 = math.exp(  (s_1/cv + (gamma-1)*math.log( abs(p_1) ))/gamma  )
    u_1 = math.sqrt(  2*abs(h0_1 - cp*T_1)  )
    v_1 = u_1*A_1/m_dot
    eq1 = p_1*v_1 - R/M*T_1

    rel_corrected_flow = (m_dot/rho_0_1/c0_1)/(m_dot_rated_c/rho0_rated/c0_rated)
    rel_corrected_speed = (omega/c0_1)/(omega_rated/c0_rated)

# Combine Eq 9, 10, 11, 12, 13, 14
    p0_2 = p0_1 * interpolate_2d_comp(rel_corrected_flow, rel_corrected_speed)
    v0_2_is = R / M / p0_2 * math.exp((s_1 + cv * (gamma - 1) * math.log( abs(p0_2) )) / gamma / cv)
    T0_2_is = p0_2 * v0_2_is * M / R
    h0_2_is = cv * T0_2_is + p0_2 * v0_2_is
    h0_2 = (h0_2_is - h0_1 + eff * h0_1) / eff
    eq2 = tau - m_dot/omega * (h0_2 - h0_1)

#Eq 16, and replace u_2 -> Eq 18,and 17
    A = 0.5*m_dot**2*R**2/A_2**2/M**2/p_2**2; B = cp; C = -h0_2
    Delta = B**2 - 4*A*C
    T_2 = (  -B + math.sqrt(Delta)  )/2/A
    u_2 = math.sqrt(2 * (h0_2 - cp * T_2))
    v_2 = u_2 * A_2 / m_dot
# Eq 15
    eq3 = (u_1**2/v_1 + p_1)*A_1 - (u_2**2/v_2 + p_2)*A_2 + (p0_2 - p0_1)*A_ref

    return [eq1, eq2, eq3]

result =  fsolve(equation, [4.0, 8000, 101325])
print(result)  #

And MOOSE input file

# (I)-----(C)-----(O)
#      1       2    
#

omega_rated = 10053
c0_rated = 358.13321019260536
rho0_rated = 1.146881112 
m_dot_rated_c = 2.9319973795574246
#m_dot_in = 4.34
#T_in = 298.15
Diameter = 0.15
x1 = 0.0
L1 = 10.0
L2 = 10.0
n_elems1 = 100
n_elems2 = 100
p_ambient = 100000
T_ambient = 298.15
I_comp = 1.0
I_motor = 1.0
const_eff = 0.79

[GlobalParams]
  #For ShaftConnectedCompressor1Phase
  initial_p = ${p_ambient}
  initial_T = ${T_ambient}       
  initial_vel_x = 0
  initial_vel_y = 0
  initial_vel_z = 0
  scaling_factor_rhoV = 1            # Scaling factor for rho*V
  scaling_factor_rhouV = 1           # Scaling factor for rho*u*V 
  scaling_factor_rhovV = 1           # Scaling factor for rho*v*V 
  scaling_factor_rhowV = 1           # Scaling factor for rho*w*V 
  scaling_factor_rhoEV = 1e-5        # Scaling factor for rho*E*V

  #For FlowChannel1Phase
  fp = fp_air                        # Fluid properties user object
  closures = closures                # Closures type
  f = 0                              # Wall friction factor
  initial_vel = 0                    # Initial velocity in the flow channel [m/s]
  scaling_factor_1phase = '1 1 1e-5' # Scaling factors for each single phase variable (rhoA, rhouA, rhoEA)
  rdg_slope_reconstruction = none    # Slope reconstruction type for rDG spatial discretization 
  orientation = '1 0 0'
  gravity_vector = '0 0 0'
[]

[FluidProperties]
  [fp_air]
    type = IdealGasFluidProperties
  []
[]

[Closures]
  [closures]
    type = Closures1PhaseSimple
  []
[]

[Components]
  [shaft]
    type = Shaft
    connected_components = 'compressor motor'
    initial_speed = 0.0
  []

  [inlet]
    type = InletStagnationPressureTemperature1Phase
    input = 'pipe1:in'
    p0 = 129155  
    T0 = 319.56  
  []
  [pipe1]
    type = FlowChannel1Phase
    position = '${x1} 0 0'
    length = ${L1}
    n_elems = ${n_elems1}
    A = ${fparse pi/4*Diameter^2}
  []

  [compressor]
    type = ShaftConnectedCompressor1Phase
    position = '${fparse x1 + L1} 0 0'
    inlet = 'pipe1:out'
    outlet = 'pipe2:in'
    A_ref = ${fparse pi/4*Diameter^2}
    volume = ${fparse pi/4*Diameter^2 * 1.0}

    omega_rated = ${omega_rated}
    mdot_rated = ${m_dot_rated_c}
    c0_rated = ${c0_rated}
    rho0_rated = ${rho0_rated}

    speeds = '0.5208 0.6250 0.7292 0.8333 0.9375'
    Rp_functions =  'rp_comp1 rp_comp2 rp_comp3 rp_comp4 rp_comp5'
    eff_functions = 'eff_comp1 eff_comp1 eff_comp1 eff_comp1 eff_comp1'

    min_pressure_ratio = 1.0

    speed_cr_I = 0
    inertia_const = ${I_comp}
    inertia_coeff = '${I_comp} 0 0 0'

    # assume no shaft friction
    speed_cr_fr = 0
    tau_fr_const = 0
    tau_fr_coeff = '0 0 0 0'
  []

  [motor]
    type = ShaftConnectedMotor
    inertia = ${I_motor}
    torque = 49.45
  []
  
  [pipe2]
    type = FlowChannel1Phase
    position = '${fparse x1+L1} 0 0'
    length = ${L2}
    n_elems = ${n_elems2}
    A = ${fparse pi/4*Diameter^2}
  []  

  [outlet]
    type = Outlet1Phase
    input = 'pipe2:out'
    p = 360000
  []

[]

[Executioner]
  type = Transient
  scheme = 'bdf2'

  end_time = 10000000000
  [TimeStepper]
    type = IterationAdaptiveDT
    dt = 0.01
    growth_factor = 1.1
    cutback_factor = 0.9
  []
  dtmin = 1e-5
  dtmax = 1000
  steady_state_detection = true
  steady_state_start_time = 200000
  steady_state_tolerance = 1e-15

  solve_type = NEWTON
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-8
  nl_max_its = 15

  l_tol = 1e-4
  l_max_its = 10

  petsc_options_iname  = '-pc_type'
  petsc_options_value  = ' lu     '
[]

[Postprocessors]
  [shaft_speed]
    type = ScalarVariable
    variable = 'shaft:omega'
    execute_on = 'INITIAL TIMESTEP_END'
  []
  [motor_torque]
    type = RealComponentParameterValuePostprocessor
    component = motor
    parameter = torque
    execute_on = 'INITIAL TIMESTEP_END'
  []
#For compressor inlet
  [a_m_dot_1]
    type = ADFlowJunctionFlux1Phase
    boundary = pipe1:out
    connection_index = 0
    equation = mass
    junction = compressor
    outputs = csv_2
  []
  [a_vel1]
    type = SideAverageValue
    variable = vel_x
    boundary = pipe1:out
    outputs = csv_2
  []
  [a_vel2]
    type = SideAverageValue
    variable = vel_y
    boundary = pipe1:out
    outputs = csv_2
  []
  [a_vel3]
    type = SideAverageValue
    variable = vel_z
    boundary = pipe1:out
    outputs = csv_2
  []
  [a_e1]
    type = SideAverageValue
    variable = e
    boundary = pipe1:out
    outputs = csv_2
  []
  [a_v1]
    type = SideAverageValue
    variable = v
    boundary = pipe1:out
    outputs = csv_2
  []
  [a_p1]
    type = SideAverageValue
    variable = p
    boundary = pipe1:out
    outputs = csv_2
  []
  [a_T1]
    type = SideAverageValue
    variable = T
    boundary = pipe1:out
    outputs = csv_2
  []
  [a_H1]
    type = SideAverageValue
    variable = H
    boundary = pipe1:out
    outputs = csv_2
  []   
#For compressor outlet
  [a_m_dot_2]
    type = ADFlowJunctionFlux1Phase
    boundary = pipe2:in
    connection_index = 0
    equation = mass
    junction = compressor
    outputs = csv_2
  []
  [a_vel1_2]
    type = SideAverageValue
    variable = vel_x
    boundary = pipe2:in
    outputs = csv_2
  []
  [a_vel2_2]
    type = SideAverageValue
    variable = vel_y
    boundary = pipe2:in
    outputs = csv_2
  []
  [a_vel3_2]
    type = SideAverageValue
    variable = vel_z
    boundary = pipe2:in
    outputs = csv_2
  []
  [a_e2]
    type = SideAverageValue
    variable = e
    boundary = pipe2:in
    outputs = csv_2
  []
  [a_v2]
    type = SideAverageValue
    variable = v
    boundary = pipe2:in
    outputs = csv_2
  []
  [a_p2]
    type = SideAverageValue
    variable = p
    boundary = pipe2:in
    outputs = csv_2
  []
  [a_T2]
    type = SideAverageValue
    variable = T
    boundary = pipe2:in
    outputs = csv_2
  []
  [a_H2]
    type = SideAverageValue
    variable = H
    boundary = pipe2:in
    outputs = csv_2
  []
[]


[Outputs]
  exodus = true
  [csv]
    type = CSV
    file_base = 'zhiyouyige_compressor'
    execute_vector_postprocessors_on = 'INITIAL'
  []
  [csv_2]
    type = CSV
    file_base = 'zhiyouyige_compressor_2'
  []
[]

[Functions]
  [rp_comp1]
    type = PiecewiseLinear
    x = '0.2868480  0.3084640  0.3300760  0.3516920  0.3733040  0.3949200  0.4165200  0.4381600  0.4597600  0.4813600  0.5030000  0.5187200'
    y = '1.5153000  1.5088000  1.5032000  1.4984000  1.4914000  1.4863000  1.4846000  1.4756000  1.4540000  1.4237000  1.3856000  1.3506000'
    extrap = true
  []
  [rp_comp2]
    type = PiecewiseLinear
    x = '0.3674120  0.3890240  0.4106400  0.4322400  0.4523200  0.4754800  0.4970800  0.5177200  0.5442400  0.5658800  0.5874800  0.6090800  0.6307200  0.6484000  0.6621600  0.6759200  0.6886800  0.6994800  0.7093200  0.7181600  0.7253600  0.7338800  0.7407600  0.7466400  0.7517600  0.7595600  0.7633600'
    y = '1.7741000  1.7626000  1.7550000  1.7536000  1.7497000  1.7455000  1.7376000  1.7276000  1.7026000  1.6780000  1.6405000  1.5914000  1.5233000  1.4564000  1.3890000  1.3187000  1.2318000  1.1631000  1.0851000  1.0143000  0.9425300  0.8601900  0.7844300  0.7160400  0.6450700  0.5733400  0.5341100'
    extrap = true
  []
  [rp_comp3]
    type = PiecewiseLinear
    x = '0.4067200  0.4283200  0.4499200  0.4715600  0.4931600  0.5147600  0.5364000  0.5580000  0.5796400  0.6012400  0.6228400  0.6444800  0.6660800  0.6872000  0.7093200  0.7299200  0.7476400  0.7643200  0.7800400  0.7928000  0.8036400  0.8134400  0.8222800  0.8301600  0.8380000  0.8448800  0.8507600  0.8566800  0.8635600  0.8704400  0.8782800'
    y = '2.1535000  2.1352000  2.1251000  2.1192000  2.1135000  2.1088000  2.1071000  2.1068000  2.1043000  2.0975000  2.0877000  2.0742000  2.0539000  2.0272000  1.9813000  1.9245000  1.8643000  1.7959000  1.7227000  1.6551000  1.5896000  1.5110000  1.4403000  1.3714000  1.2969000  1.2326000  1.1655000  1.0985000  1.0177000  0.9385600  0.8578700'
    extrap = true
  []
  [rp_comp4]
    type = PiecewiseLinear
    x = '0.4872800  0.5088800  0.5304800  0.5521200  0.5737200  0.5953600  0.6169600  0.6385600  0.6602000  0.6818000  0.7034000  0.7250400  0.7466400  0.7682400  0.7898800  0.8114800  0.8330800  0.8547200  0.8753600  0.8940000  0.9097200  0.9224800  0.9342800  0.9450800  0.9568800'
    y = '2.6634000  2.6485000  2.6366000  2.6248000  2.6184000  2.6124000  2.6110000  2.6102000  2.6110000  2.6102000  2.6110000  2.6099000  2.6060000  2.5989000  2.5874000  2.5642000  2.5309000  2.4824000  2.4202000  2.3500000  2.2832000  2.2195000  2.1437000  2.0706000  1.9947000'
    extrap = true
  []
  [rp_comp5]
    type = PiecewiseLinear
    x = '0.5658800  0.5874800  0.6090800  0.6307200  0.6523200  0.6739200  0.6955600  0.7171600  0.7387600  0.7604000  0.7820000  0.8036400  0.8252400  0.8468400  0.8684800  0.8900800  0.9116800  0.9333200  0.9549200  0.9765200  0.9981600  1.0197600  1.0413600  1.0600400  1.0757600  1.0904800  1.1032800  1.1150800  1.1248800  1.1327600  1.1406000  1.1484800  1.1563200  1.1632000  1.1690800  1.1750000  1.1808800  1.1867600  1.1916800  1.1976000'
    y = '3.4718000  3.4428000  3.4189000  3.3970000  3.3728000  3.3455000  3.3249000  3.3106000  3.2976000  3.2886000  3.2819000  3.2799000  3.2844000  3.2892000  3.2898000  3.2892000  3.2889000  3.2830000  3.2703000  3.2551000  3.2281000  3.1873000  3.1291000  3.0660000  2.9999000  2.9281000  2.8556000  2.7823000  2.7101000  2.6451000  2.5778000  2.5099000  2.4280000  2.3562000  2.2850000  2.2128000  2.1344000  2.0580000  1.9920000  1.9198000'
    extrap = true
  []
  [eff_comp1]
    type = ConstantFunction
    value = ${const_eff}
  []
[]

And you will find that the final result from MOOSE is similar to that from this Python

[GiudGiud]

Thanks for sharing!
@licharlot @joshuahansel