363 lines
32 KiB
C#
363 lines
32 KiB
C#
using System.Runtime.InteropServices;
|
||
|
||
namespace CapMachine.Wpf.PPCalculation
|
||
{
|
||
//https://github.com/usnistgov/REFPROP-wrappers/tree/master 64位的看这个
|
||
|
||
/// <summary>
|
||
/// This is C# Wrapper Class for REFPROP.
|
||
/// REFPROP.DLL 是32程序的Dll文件，如果应用程序是32位的，需要使用REFPROP.DLL，修改为： [DllImport(@".\PPCalculation\REFPROP\REFPROP.DLL")]
|
||
/// REFPRP64.DLL 是64位程序的Dll文件，如果应用程序是64位的，需要使用REFPRP64.DLL 修改为： [DllImport(@".\PPCalculation\REFPROP\REFPRP64.DLL")]
|
||
/// 还请搭配流体文件的数据（请注意路径的长度）
|
||
/// </summary>
|
||
public static class IRefProp64
|
||
{
|
||
/**************************************************************
|
||
* Set the path where the fluid files are located.
|
||
* The path does not need to contain the ending "\" and it can
|
||
* point directly to the location where the DLL is stored if a
|
||
* fluids subdirectory (with the corresponding fluid files) is
|
||
* located there.
|
||
*************************************************************/
|
||
[DllImport(@".\PPCalculation\REFPROP\REFPRP64.DLL")]
|
||
public static extern void SETPATHdll(string htype, ref long ln);
|
||
|
||
/******************************************************************************
|
||
* Define models and initialize arrays. A call to this subroutine is required
|
||
*******************************************************************************/
|
||
[DllImport(@".\PPCalculation\REFPROP\REFPRP64.DLL", CharSet = CharSet.Ansi)]
|
||
public static extern void SETUPdll
|
||
(
|
||
ref long nComps, // (INPUT) number of components (1 for pure fluid) [integer]
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string hfld, // (INPUT) list of file names specifying fluid/mixture components. Separated by pipes "|"
|
||
// e.g., METHANE|AMMONIA|ARGON
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string hfmix, // (INPUT) mixture coefficients. File name containing coefficients for mixture model, if applicable
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string hrf, // (INPUT) reference state for thermodynamic calculations [3 character string]
|
||
// 'DEF' : default reference state as specified in fluid file is applied to each pure component
|
||
// 'NBP' : h,s = 0 at pure component normal boiling point(s)
|
||
// 'ASH' : h,s = 0 for sat liquid at -40 C (ASHRAE convention)
|
||
// 'IIR' : h= 200, s = 1.0 for sat liq at 0 C (IIR convention)
|
||
// other choices are possible, but these require a separate call to SETREF
|
||
ref long ierr, // (OUTPUT) error flag:
|
||
// 0 = successful
|
||
// 101 = error in opening file
|
||
// 102 = error in file or premature end of file
|
||
// -103 = unknown model encountered in file
|
||
// 104 = error in setup of model
|
||
// 105 = specified model not found
|
||
// 111 = error in opening mixture file
|
||
// 112 = mixture file of wrong type
|
||
// 114 = nc<>nc from setmod
|
||
// -117 = binary pair not found, all parameters will be estimated
|
||
// 117 = no mixture data are available, mixture is outside the range of the model and calculations will not be made
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string herr, // (OUTPUT) error string
|
||
ref long hfldLen,
|
||
ref long hfmixLen,
|
||
ref long hrfLen,
|
||
ref long herrLen);
|
||
|
||
|
||
/***************************************************************************************
|
||
* General Flash subroutine taking a temperature and pressure.
|
||
* This routine accepts both single-phase and two-phase states as input.
|
||
* If the the phase is know, the specialized routines are faster.
|
||
**************************************************************************************/
|
||
[DllImport(@".\PPCalculation\REFPROP\REFPRP64.DLL", CharSet = CharSet.Ansi)]
|
||
public static extern void TPFLSHdll
|
||
(
|
||
ref double t, // (INPUT) temperature [K]
|
||
ref double p, // (INPUT) pressure [kPa]
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] x, // (INPUT) overall (bulk) composition [array of mol frac]
|
||
ref double d, // (OUTPUT) overall (bulk) molar density [mol/L]
|
||
ref double Dl, // (OUTPUT) molar density of liquid phase [mol/L]
|
||
ref double Dv, // (OUTPUT) molar density of the vapor phase [mol/L]
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] xliq, // (OUTPUT) composition of liquid phase [array of mol frac]
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] xvap, // (OUTPUT) composition of vapor phase [array of mol frac]
|
||
ref double q, // (OUTPUT) vapor quality on a MOLAR basis [moles vapor/total moles]
|
||
// q < 0 indicates subcooled (compressed) liquid
|
||
// q = 0 indicates saturated liquid
|
||
// q = 1 indicates saturated vapor
|
||
// q > 1 indicates superheated vapor
|
||
// q = 998 superheated vapor, but quality not defined (in most situations, t > Tc)
|
||
// q = 999 indicates supercritical state (t > Tc) and (p > Pc)
|
||
ref double ee, // (OUTPUT) overall (bulk) internal energy [J/mol]
|
||
ref double h, // (OUTPUT) overall (bulk) enthalpy [J/mol]
|
||
ref double ss, // (OUTPUT) overall (bulk) entropy [J/mol-K]
|
||
ref double Cv, // (OUTPUT) isochoric (constant V) heat capacity [J/mol-K]
|
||
ref double Cp, // (OUTPUT) isobaric (constatnt p) heat capacity [J/mol-K] (not defined for 2-phase)
|
||
ref double w, // (OUTPUT) speed of sound [m/s] (not defined for 2-phase)
|
||
ref long ierr, // (OUTPUT) error flag
|
||
// 0 = all inputs within limits
|
||
// <> 0 = one or more inputs outside limits
|
||
// -1 = 1.5*tmax > t > tmax
|
||
// 1 = t < tmin or t > 1.5*tmax
|
||
// 2 = D > Dmax or D < 0
|
||
// -4 = 2*pmax > p > pmax
|
||
// 4 = p < 0 or p > 2*pmax
|
||
// 8 = component composition < 0 or > 1 and/or composition sum < 0 or > 1
|
||
// 16 = p > pmelt
|
||
// -16 = t < ttrp (important for water)
|
||
// If multiple inputs are outside limits, ierr = abs[sum(ierr)] with the sign determined by the most severe excursion
|
||
// (ierr < 0 indicates a warning -- results may be questionable, ierr > 0 indicates an error -- calculations impossible)
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string herr, // (OUTPUT) error string
|
||
ref long ln
|
||
);
|
||
|
||
|
||
|
||
/****************************************************************************************
|
||
* Open a mixture file (e.g., R410A.mix) and read constituents and mole fractions
|
||
****************************************************************************************/
|
||
[DllImport(@".\PPCalculation\REFPROP\REFPRP64.DLL", CharSet = CharSet.Ansi)]
|
||
public static extern void SETMIXdll
|
||
(
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string hmxnme, // (INPUT) mixture file name to be read in
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string hfmix, // (INPUT) mixture coefficients. File name containing coefficients for mixture model
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string hrf, // (INPUT) reference state for thermodynamic calculations (see same variable in SETUPdll)
|
||
ref long nccp, // (OUTPUT) number of fluids in mixture
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string hfiles, // (OUTPUT) array of file names specifying mixture components that were used to call setup
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] x, // (OUTPUT) array of mole fractions for the specified mixture
|
||
ref long ierr, // (OUTPUT) error flag
|
||
// 0 = successful
|
||
// 101 = error in opening file
|
||
// -102 = mixture file contains mixing parameters
|
||
// -103 = composition not equal to 1
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string herr, // (OUTPUT) error string
|
||
ref long hfldLen,
|
||
ref long hfmixLen,
|
||
ref long hrfLen,
|
||
ref long hfilesLen,
|
||
ref long herrLen);
|
||
|
||
/**********************************************************************************************
|
||
* iterate for saturated liquid and vapor states given temperature and the composition of one phase
|
||
***********************************************************************************************/
|
||
[DllImport(@".\PPCalculation\REFPROP\REFPRP64.DLL", CharSet = CharSet.Ansi)]
|
||
public static extern void SATTdll //!< a member function.
|
||
(
|
||
ref double t, // (INPUT) temperature [K]
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] x, // (INPUT) composition [array of mol frac]
|
||
ref long kph, // (INPUT) phase flag
|
||
// 1 = input x is liquid composition (bubble point)
|
||
// 2 = input x is vapor composition (dew point)
|
||
// 3 = input x is liquid composition (freezing point)
|
||
// 4 = input x is vapor composition (sublimation point)
|
||
ref double p, // (OUTPUT) pressure [kPa]
|
||
ref double Dl, // (OUTPUT) molar density of saturated liquid [mol/L]
|
||
ref double Dv, // (OUTPUT) molar density of saturated vapor [mol/L]
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] xliq, // (OUTPUT) liquid phase composition [array of mol frac]
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] xvap, // (OUTPUT) vapor phase composition [array of mol frac]
|
||
ref long ierr, // (OUTPUT) error flag
|
||
// 0 = successful
|
||
// 1 = T < Tmin
|
||
// 8 = x out of range
|
||
// 9 = T and x out of range
|
||
// 120 = CRITP did not converge
|
||
// 121 = T > Tcrit
|
||
// 122 = TPRHO-liquid did not converge (pure fluid)
|
||
// 123 = TPRHO-vapor did not converge (pure fluid)
|
||
// 124 = pure fluid iteration did not converge
|
||
// following 3 error codes are advisory -- iteration will either converge on later guess or error out
|
||
// -125 = TPRHO did not converge for parent ph (mix)
|
||
// -126 = TPRHO did not converge for incipient (mix)
|
||
// -127 = composition iteration did not converge
|
||
// 128 = mixture iteration did not converge (error out from previous 3)
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string herr, // (OUTPUT) error string
|
||
ref long herrLen
|
||
);
|
||
|
||
|
||
/*******************************************************************************************************************************
|
||
* calculates the phase boundary of a mixture at a given composition, and the critical point, cricondentherm, and cricondenbar
|
||
******************************************************************************************************************************/
|
||
[DllImport(@".\PPCalculation\REFPROP\REFPRP64.DLL", CharSet = CharSet.Ansi)]
|
||
public static extern void SATSPLNdll
|
||
(
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] x, // (INPUT) composition [array of mol frac]
|
||
ref long ierr, // (OUTPUT) error flag: 0 = successful
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string herr, // (OUTPUT) error string
|
||
ref long herrLen
|
||
);
|
||
|
||
|
||
/*************************************************************************
|
||
* converts quality and composition on a mole basis to a mass basis
|
||
*************************************************************************/
|
||
[DllImport(@".\PPCalculation\REFPROP\REFPRP64.DLL", CharSet = CharSet.Ansi)]
|
||
public static extern void QMASSdll
|
||
(
|
||
ref double qmol, // (INPUT) molar quality [moles vapor/total moles] See parameter 'q' in TPFLSHdll
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] xl, // (INPUT) composition of liquid phase [array of mol frac]
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] xv, // (INPUT) composition of vapor phase [array of mol frac]
|
||
ref double qkg, // (OUTPUT) quality on mass basis [mass of vapor/total mass]
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] xlkg, // (OUTPUT) mass composition of liquid phase [array of mass frac]
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] xvkg, // (OUTPUT) mass composition of vapor phase [array of mass frac]
|
||
ref double wliq, // (OUTPUT) molecular weight of liquid phase [g/mol]
|
||
ref double wvap, // (OUTPUT) molecular weight of vapor phase [g/mol]
|
||
ref long ierr, // (OUTPUT) error flag: 0 = all inputs within limits, -19 = q < 0 or > 1
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string herr,
|
||
ref long herrLen
|
||
);
|
||
|
||
|
||
/*************************************************************************************
|
||
* compute thermal quantities as a function of temperature, density, and compositions using core
|
||
* functions (Helmholtz free energy, ideal gas heat capacity and various derivatives and integrals)
|
||
************************************************************************************/
|
||
[DllImport(@".\PPCalculation\REFPROP\REFPRP64.DLL", CharSet = CharSet.Ansi)]
|
||
public static extern void THERMdll
|
||
(
|
||
ref double t, // temperature [K]
|
||
ref double Dl, // molar density [mol/L]
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] x, // composition [array of mol frac]
|
||
ref double p, // (OUTPUT) pressure [kPa]
|
||
ref double el, // (OUTPUT) internal energy [J/mol]
|
||
ref double hl, // (OUTPUT) enthalpy [J/mol]
|
||
ref double sl, // (OUTPUT) entropy [J/mol-K]
|
||
ref double cvl, // (OUTPUT) isochoric heat capacity [J/mol-K]
|
||
ref double cpl, // (OUTPUT) isobaric heat capacity [J/mol-K]
|
||
ref double wl, // (OUTPUT) speed of sound [m/s]
|
||
ref double hjt // (OUTPUT) isenthalpic Joule-Thomson coefficient [K/kPa]
|
||
);
|
||
|
||
|
||
/********************************************************************************
|
||
* Solve for density given T and P using the most robust iterative routine.
|
||
* This mirrors the LabVIEW TPRHO DLL call path to obtain phase-specific density.
|
||
********************************************************************************/
|
||
[DllImport(@".\PPCalculation\REFPROP\REFPRP64.DLL", CharSet = CharSet.Ansi)]
|
||
public static extern void TPRHOdll
|
||
(
|
||
ref double t, // (INPUT) temperature [K]
|
||
ref double p, // (INPUT) pressure [kPa]
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] x, // (INPUT) composition [array of mol frac]
|
||
ref long kph, // (INPUT) phase flag: 1=liquid, 2=vapor
|
||
ref long kguess, // (INPUT) initial guess: 0=let routine choose
|
||
ref double D, // (OUTPUT) molar density [mol/L]
|
||
ref long ierr, // (OUTPUT) error flag: 0=success
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string herr, // (OUTPUT) error string
|
||
ref long herrLen
|
||
);
|
||
|
||
|
||
|
||
/******************************************************************
|
||
* compute the transport properties of thermal conductivity and
|
||
* viscosity as functions of temperature, density, and composition
|
||
*******************************************************************/
|
||
[DllImport(@".\PPCalculation\REFPROP\REFPRP64.DLL", CharSet = CharSet.Ansi)]
|
||
public static extern void TRNPRPdll
|
||
(
|
||
ref double t, // (INPUT) temperature [K]
|
||
ref double d, // (INPUT) molar density [mol/L]
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] x, // (INPUT) composition array [mol frac]
|
||
ref double eta, // (OUTPUT) viscosity [uPa.s]
|
||
ref double tcx, // (OUTPUT) thermal conductivity [W/m.K]
|
||
ref long ierr, // (OUTPUT) error flag:
|
||
// 0 = successful
|
||
// -31 = temperature out of range for conductivity
|
||
// -32 = density out of range for conductivity
|
||
// -33 = T and D out of range for conductivity
|
||
// -41 = temperature out of range for viscosity
|
||
// -42 = density out of range for viscosity
|
||
// -43 = T and D out of range for viscosity
|
||
// -51 = T out of range for both visc and t.c.
|
||
// -52 = D out of range for both visc and t.c.
|
||
// -53 = T and/or D out of range for visc and t.c.
|
||
// 39 = model not found for thermal conductivity
|
||
// 40 = model not found for thermal conductivity or viscosity
|
||
// 49 = model not found for viscosity
|
||
// 50 = ammonia/water mixture (no properties calculated)
|
||
// 51 = exactly at t, rhoc for a pure fluid; k is infinite
|
||
// -58,-59 = ECS model did not converge
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string herr,
|
||
ref long herrLen
|
||
);
|
||
|
||
|
||
|
||
/***************************************************
|
||
* critical parameters as a function of composition
|
||
***************************************************/
|
||
[DllImport(@".\PPCalculation\REFPROP\REFPRP64.DLL", CharSet = CharSet.Ansi)]
|
||
public static extern void CRITPdll
|
||
(
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] x, // (INPUT) composition [array of mol frac]
|
||
ref double tc, // (OUTPUT) critical temperature [K]
|
||
ref double pc, // (OUTPUT) critical pressure [kPa]
|
||
ref double Dc, // (OUTPUT) critical density [mol/L]
|
||
ref long ierr, // (OUTPUT) error flag: 0=successful, 1=did not converge
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string herr,
|
||
ref long herrLen
|
||
);
|
||
|
||
|
||
/*************************************************************************************************
|
||
* iterate for saturated liquid and vapor states given pressure and the composition of one phase
|
||
************************************************************************************************/
|
||
[DllImport(@".\PPCalculation\REFPROP\REFPRP64.DLL", CharSet = CharSet.Ansi)]
|
||
public static extern void SATPdll
|
||
(
|
||
ref double p, // (INPUT) pressure [kPa]
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] x, // (INPUT) composition array [array of mol frac]
|
||
ref long kph, // (INPUT) phase flag. see same parameter in SATPdll
|
||
ref double t, // (OUTPUT) temperature [K]
|
||
ref double Dl, // (OUTPUT) density of liquid [mol/L]
|
||
ref double Dv, // (OUTPUT) density of vapor [mol/L]
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] xliq, // (OUTPUT) liquid composition [array of mol frac]
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] xvap, // (OUTPUT) vapor composition [array of mol frac]
|
||
ref long ierr, // (OUTPUT) error flag:
|
||
// 0 = successful
|
||
// 2 = P<Ptp
|
||
// 4 = P< 0
|
||
// 8 = x out of range
|
||
// 12 = P and x out of range
|
||
// 140 = CRITP did not converge
|
||
// 141 = P> Pcrit
|
||
// 142 = TPRHO-liquid did not converge (pure fluid)
|
||
// 143 = TPRHO-vapor did not converge (pure fluid)
|
||
// 144 = pure fluid iteration did not converge
|
||
// following 3 error codes are advisory--iteration will either converge on later guess or error out (ierr = 148)
|
||
// -144 = Raoult's law (mixture initial guess) did not converge
|
||
// -145 = TPRHO did not converge for parent ph(mix)
|
||
// -146 = TPRHO did not converge for incipient(mix)
|
||
// -147 = composition iteration did not converge
|
||
// 148 = mixture iteration did not converge (error out from previous 3)
|
||
[MarshalAs(UnmanagedType.VBByRefStr)] ref string herr, // (OUTPUT) error string
|
||
ref long herrLen
|
||
);
|
||
|
||
|
||
/***********************************************************
|
||
* molecular weight for a mixture of specified composition
|
||
***********************************************************/
|
||
[DllImport(@".\PPCalculation\REFPROP\REFPRP64.DLL", CharSet = CharSet.Ansi)]
|
||
public static extern void WMOLdll
|
||
(
|
||
[MarshalAs(UnmanagedType.LPArray, SizeParamIndex = 0)] double[] x, // (INPUT) composition array [array of mol frac]
|
||
ref double wm // (OUTPUT) molar mass [g/mol]
|
||
);
|
||
|
||
|
||
/******************************************************
|
||
* Provides fluid constants for specified component
|
||
******************************************************/
|
||
[DllImport(@".\PPCalculation\REFPROP\REFPRP64.DLL", CharSet = CharSet.Ansi)]
|
||
public static extern void INFOdll
|
||
(
|
||
ref long icomp, // (INPUT) Component number in mixture. 1 for pure fluid
|
||
ref double wmm, // (OUTPUT) molecular weight [g/mol]
|
||
ref double ttrp, // (OUTPUT) triple point temperature [K]
|
||
ref double tnbpt, // (OUTPUT) normal boiling point temperature [K]
|
||
ref double tc, // (OUTPUT) critical temperature [K]
|
||
ref double pc, // (OUTPUT) critical pressure [kPa]
|
||
ref double Dc, // (OUTPUT) critical density [mol/L]
|
||
ref double Zc, // (OUTPUT) compressibility at critical point [pc/(Rgas*Tc*Dc)]
|
||
ref double acf, // (OUTPUT) acentric factor [dimensionless]
|
||
ref double dip, // (OUTPUT) dipole moment [debeye]
|
||
ref double Rgas // (OUTPUT) gas constant [J/mol-K]
|
||
);
|
||
|
||
}
|
||
}
|