modelDynamicTwoPhaseFlowRiser

Dynamic two-phase flow riser

Extends from ThermoSysPro.Fluid.Interfaces.IconColors (Colors of the static, singular and dynamic model icons).

Information

## Copyright © EDF 2002 - 2026   
## ThermoSysPro Version 4.2

Parameters

Type	Name	Default	Description
Units.SI.Length	L	10.	Pipe length
Units.SI.Diameter	D	0.02	Hydraulic diameter
Units.SI.Length	rugosrel	0.0007	Pipe relative roughness
Integer	ntubes	1	Number of pipes in parallel
Units.SI.Position	z1	0	Pipe inlet altitude
Units.SI.Position	z2	0	Pipe outlet altitude
Real	rgliss	1	Phase slip coefficient
Integer	a	4200	Phase pressure loss coefficient
Real	dpfCorr	1.00	Corrective term for the friction pressure loss (dpf) for each node
Real	hcCorr	1.00	Corrective term for the heat exchange coefficient (hc) for each node
Integer	Ns	10	Number of segments
Units.SI.Temperature[Ns]	T0	fill(300, Ns)	Initial fluid temperature (active if steady_state=false and option_temperature=true)
Units.SI.SpecificEnthalpy[Ns]	h0	fill(1e5, Ns)	Initial fluid specific enthalpy (active if steady_state=false and option_temperature=false)
Boolean	inertia	true	true: momentum balance equation with inertia - false: without inertia
Boolean	advection	false	true: momentum balance equation with advection term - false: without advection term
Boolean	dynamic_energy_balance	true	true: dynamic energy balance equation - false: static energy balance equation
Boolean	dynamic_mass_balance	true	true: dynamic mass balance equation - false: static mass balance equation (active if dynamic_energy_balance=true)
Boolean	simplified_dynamic_energy_balance	true	true: simplified dynamic energy balance equation - false: full dynamic energy balance equation (active if dynamic_energy_balance=true and dynamic_mass_balance=true)
Boolean	steady_state	true	true: start from steady state - false: start from T0 (if option_temperature = true) or h0 (if option_temperature=false)(active if dynamic_energy_balance=true)
Boolean	option_temperature	true	true: initial temperature is fixed - false: initial specific enthalpy is fixed (active if steady_state=false)
Boolean	continuous_flow_reversal	false	true: continuous flow reversal - false: discontinuous flow reversal
Boolean	diffusion	false	true: energy balance equation with diffusion - false: energy balance equation without diffusion
Fluid › Fluid properties
IF97Region	region	IF97Region.All_regions	IF97 region (active for IF97 water/steam only)

Connectors

Type	Name	Default	Description
Interfaces.Connectors.FluidInlet	C1
ThermoSysPro.Thermal.Connectors.ThermalPort[Ns]	CTh1
Interfaces.Connectors.FluidOutlet	C2
ThermoSysPro.Thermal.Connectors.ThermalPort[Ns]	CTh2

Components

Type	Name	Default	Description
Units.SI.AbsolutePressure[N + 1]	P		Fluid pressure in node i
Units.SI.MassFlowRate[N]	Q		Mass flow rate in node i
Units.SI.SpecificEnthalpy[N + 1]	h		Fluid specific enthalpy in node i
Units.SI.SpecificEnthalpy[N]	hb		Fluid specific enthalpy at the boundary of node i
Units.SI.AbsolutePressure[N + 1]	Pb		Bounded fluid pressure in node i
Units.SI.Density[N - 1]	rho1		Fluid density in thermal node i
Units.SI.Density[N]	rho2		Fluid density in hydraulic node i
Units.SI.Density[N + 1]	rhoc		Fluid density at the boudary of node i
Units.SI.MassFlowRate[N - 1]	BQ		Right hand side of the mass balance equation for thermal node i
Units.SI.Power[N - 1]	BH		Right hand side of the energy balance equation for thermal node i
Units.SI.Power[N - 1]	dW1_conv		Convection thermal power exchanged on the liquid side for node i
Units.SI.Power[N - 1]	dW1_diff		Diffusion thermal power exchanged on the liquid side for node i
Units.SI.Power[N - 1]	dW1		Thermal power exchanged on the liquid side for node i
Units.SI.Power[N - 1]	dW2_conv		Convection thermal power exchanged on the liquid side for node i
Units.SI.Power[N - 1]	dW2_diff		Diffusion thermal power exchanged on the liquid side for node i
Units.SI.Power[N - 1]	dW2		Thermal power exchanged on the liquid side for node i
Units.SI.Power	W1t		Total power exchanged on the liquid side
Units.SI.Temperature[N - 1]	Tp1		Wall temperature in node i
Units.SI.Temperature[N - 1]	Tp2		Wall temperature in node i
Units.SI.CoefficientOfHeatTransfer[N - 1]	hi		Fluid heat exchange coefficient in node i
Units.SI.CoefficientOfHeatTransfer[N - 1]	hcl		Fluid heat exchange coefficient in node i for the liquid fraction
Units.SI.CoefficientOfHeatTransfer[N - 1]	hcv		Fluid heat exchange coefficient in node i for the vapor fraction
Real[N - 1]	S		Corrective terme correctif for nucleation removal
Real[N - 1]	E		Corrective term for hcl
Units.SI.CoefficientOfHeatTransfer[N - 1]	heb		Fluid heat exchange coefficient for vaporization in thermal node i
Units.SI.ReynoldsNumber[N - 1]	Rel1		Reynolds number in thermal node i for the liquid
Units.SI.ReynoldsNumber[N]	Rel2		Reynolds number in hydraulic node i for the liquid
Units.SI.ReynoldsNumber[N - 1]	Rev1		Reynolds number in thermal node i for the vapor
Units.SI.ReynoldsNumber[N]	Rev2		Reynolds number in hydraulic node i for the vapor
Real[N - 1]	Prl		Fluid Prandtl number in node i for the liquid
Real[N - 1]	Prv		Fluid Prandtl number in node i for the vapor
Units.SI.ThermalConductivity[N - 1]	kl1		Thermal conductivity in thermal node i for the liquid
Units.SI.ThermalConductivity[N]	kl2		Thermal conductivity in hydraulic node i for the liquid
Units.SI.ThermalConductivity[N - 1]	kv1		Thermal conductivity in thermal node i for the vapor
Units.SI.ThermalConductivity[N]	kv2		Thermal conductivity in hydraulic node i for the vapor
Real[N - 1]	xv1		Vapor mass fraction in thermal node i
Real[N]	xv2		Vapor mass fraction in hydraulic node i
Real[N - 1]	xbs		Bounded upper value for the vapor mass fraction
Real[N - 1]	xbi		Bounded lower value for the vapor mass fraction
Units.SI.DynamicViscosity[N - 1]	mul1		Dynamic viscosity in thermal node i for the liquid
Units.SI.DynamicViscosity[N]	mul2		Dynamic viscosity in hydraulic node i for the liquid
Units.SI.DynamicViscosity[N - 1]	muv1		Dynamic viscosity in thermal node i for the vapor
Units.SI.DynamicViscosity[N]	muv2		Dynamic viscosity in hydraulic node i for the vapor
Units.SI.SpecificHeatCapacity[N - 1]	cpl1		Specific heat capacity in thermal node i for the liquid
Units.SI.SpecificHeatCapacity[N]	cpl2		Specific heat capacity in hydraulic node i for the liquid
Units.SI.SpecificHeatCapacity[N - 1]	cpv1		Specific heat capacity in thermal node i for the vapor
Units.SI.SpecificHeatCapacity[N]	cpv2		Specific heat capacity in hydraulic node i for the vapor
Real[N - 1]	Bo		Boiling number
Real[N - 1]	Xtt		Martinelli number
Units.SI.SpecificEnthalpy[N - 1]	lv		Specific enthalpy for vaporisation
Units.SI.Density[N - 1]	rhol1		Fluid density in thermal node i for the liquid
Units.SI.Density[N]	rhol2		Fluid density in hydraulic node i for the liquid
Units.SI.Density[N - 1]	rhov1		Fluid density in thermal node i for the vapor
Units.SI.Density[N]	rhov2		Fluid density in hydraulic node i for the vapor
Units.SI.Temperature[N - 1]	T1		Fluid temperature in thermal node i
Units.SI.Temperature[N]	T2		Fluid temperature in hydraulic node i
ThermoSysPro.Units.SI.PressureDifference[N]	dpa		Advection term for the mass balance equation in node i
ThermoSysPro.Units.SI.PressureDifference[N]	dpf		Friction pressure loss in node i
ThermoSysPro.Units.SI.PressureDifference[N]	dpg		Gravity pressure loss in node i
Real[N]	khi		Hydraulic pressure loss coefficient in node i
Real[N]	lambdal		Friction pressure loss coefficient in node i for the liquid
Real[N]	lambdav		Friction pressure loss coefficient in node i for the vapor)
Real[N]	filo		Pressure loss coefficient for two-phase flow
FluidType	ftype		Fluid type
Integer	fluid	Integer(ftype)	Fluid number
ThermoSysPro.Units.SI.MassFraction	Xco2		CO2 mass fraction
ThermoSysPro.Units.SI.MassFraction	Xh2o		H2O mass fraction
ThermoSysPro.Units.SI.MassFraction	Xo2		O2 mass fraction
ThermoSysPro.Units.SI.MassFraction	Xso2		SO2 mass fraction
Real[N]	diff_res_l		Diffusion resistance in hydraulic node i for the liquid
Real[N]	diff_res_v		Diffusion resistance in hydraulic node i for the vapor
Real[N]	diff_res		Total diffusion resistance in hydraulic node i
Real	diff_res_t		Total diffusion resistance in the pipe
Real[N - 1]	diff_res_e		Diffusion resistance at inlet of thermal node i
Real[N - 1]	diff_res_s		Diffusion resistance at outlet of thermal node i
Units.SI.MassFlowRate[N]	gamma		Total diffusion conductance in hydraulic node i
Units.SI.MassFlowRate[N - 1]	gamma_e		Diffusion conductance at inlet of thermal node i
Units.SI.MassFlowRate[N - 1]	gamma_s		Diffusion conductance at outlet of thermal node i
Units.SI.Power[N - 1]	Je		Thermal power diffusion from inlet of thermal node i
Units.SI.Power[N - 1]	Js		Thermal power diffusion from outlet of thermal node i
Units.SI.Power[N - 1]	J		Total thermal power diffusion of thermal node i
Real[N - 1]	re		Value of r(Q/gamma) for inlet of thermal node i
Real[N - 1]	rs		Value of r(Q/gamma) for outlet of thermal node i
ThermoSysPro.Properties.WaterSteam.Common.ThermoProperties_ph[N - 1]	pro1
ThermoSysPro.Properties.WaterSteam.Common.ThermoProperties_ph[2]	proc
ThermoSysPro.Properties.WaterSteam.Common.ThermoProperties_ph[N]	pro2
ThermoSysPro.Properties.WaterSteam.Common.PropThermoSat[N]	vsat2
ThermoSysPro.Properties.WaterSteam.Common.PropThermoSat[N]	lsat2
ThermoSysPro.Properties.WaterSteam.Common.PropThermoSat[N - 1]	vsat1
ThermoSysPro.Properties.WaterSteam.Common.PropThermoSat[N - 1]	lsat1

Revisions

Authors  

Baligh El Hefni  
Daniel Bouskela