Class Leachman

java.lang.Object
neqsim.thermo.util.leachman.Leachman
public class Leachman extends Object
The Leachman class provides methods to calculate the density, pressure, and various thermodynamic properties of hydrogen using the Leachman method. It supports different types of hydrogen: normal, para, and ortho.

The class includes methods for setting up the parameters for different hydrogen types, calculating the density from temperature and pressure, and calculating various thermodynamic properties.

The main methods include:

The class also includes helper methods for calculating the ideal gas and real gas Helmholtz energy and their derivatives.

Example usage: 


Leachman leachman = new Leachman();
leachman.SetupLeachman("normal");
doubleW density = new doubleW(30);
doubleW pressure = new doubleW(500.0);
intW ierr = new intW(0);
StringW herr = new StringW("");
leachman.DensityLeachman(0, 500, pressure.val, density, ierr, herr);
System.out.println("Density: " + density.val);

Note: The class uses the doubleW and intW classes from the org.netlib.util package to handle output parameters.

Author:
victor
See Also:
  • doubleW
  • intW
  • StringW

  • Field Summary

    Fields
    Modifier and Type
    Field
    Description
    (package private) double[]
    a0k
     
    (package private) double[]
    a0k_normal
     
    (package private) double[]
    a0k_ortho
     
    (package private) double[]
    a0k_para
     
    (package private) double[]
    b0k
     
    (package private) double[]
    b0k_normal
     
    (package private) double[]
    b0k_ortho
     
    (package private) double[]
    b0k_para
     
    (package private) double[]
    beta_i
     
    (package private) double[]
    beta_i_normal
     
    (package private) double[]
    beta_i_ortho
     
    (package private) double[]
    beta_i_para
     
    (package private) double[]
    d_i
     
    (package private) double[]
    D_i
     
    (package private) double[]
    d_i_normal
     
    (package private) double[]
    D_i_normal
     
    (package private) double[]
    d_i_ortho
     
    (package private) double[]
    D_i_ortho
     
    (package private) double[]
    d_i_para
     
    (package private) double[]
    D_i_para
     
    (package private) double
    Dc
     
    (package private) double
    dPdDsave
     
    (package private) double
    epsilon
     
    (package private) double[]
    gamma_i
     
    (package private) double[]
    gamma_i_normal
     
    (package private) double[]
    gamma_i_ortho
     
    (package private) double[]
    gamma_i_para
     
    private String
    hydrogenType
     
    (package private) int
    K
     
    (package private) int
    K_normal
     
    (package private) int
    K_ortho
     
    (package private) int
    K_para
     
    (package private) int
    l
     
    (package private) int
    m
     
    (package private) double
    M_L
     
    (package private) int
    n
     
    (package private) double[]
    N_i
     
    (package private) double[]
    N_i_normal
     
    (package private) double[]
    N_i_ortho
     
    (package private) double[]
    N_i_para
     
    (package private) double[]
    p_i
     
    (package private) double[]
    p_i_normal
     
    (package private) double[]
    p_i_ortho
     
    (package private) double[]
    p_i_para
     
    (package private) double[]
    phi_i
     
    (package private) double[]
    phi_i_normal
     
    (package private) double[]
    phi_i_ortho
     
    (package private) double[]
    phi_i_para
     
    (package private) double
    R_L
     
    (package private) double[]
    t_i
     
    (package private) double[]
    t_i_normal
     
    (package private) double[]
    t_i_ortho
     
    (package private) double[]
    t_i_para
     
    (package private) double
    Tc
     

  • Constructor Summary

    Constructors
    Constructor
    Description
    Leachman()
     

  • Method Summary

    Modifier and Type
    Method
    Description
    (package private) void
    Alpha0Leachman(double T, double D, org.netlib.util.doubleW[] a0)
    Calculate the ideal gas Helmholtz energy and its derivatives with respect to tau and delta.
    (package private) void
    AlpharLeachman(int itau, int idelta, double T, double D, org.netlib.util.doubleW[][] ar)
    Calculate the residual Helmholtz energy and its derivatives with respect to tau and delta.
    void
    DensityLeachman(int iFlag, double T, double P, org.netlib.util.doubleW D, org.netlib.util.intW ierr, org.netlib.util.StringW herr)
    Calculate density as a function of temperature and pressure using the Leachman method.
    static void
    main(String[] args)
    main.
    void
    PressureLeachman(double T, double D, org.netlib.util.doubleW P, org.netlib.util.doubleW Z)
    Calculate the pressure and compressibility factor using the Leachman method.
    void
    propertiesLeachman(double T, double D, org.netlib.util.doubleW P, org.netlib.util.doubleW Z, org.netlib.util.doubleW dPdD, org.netlib.util.doubleW d2PdD2, org.netlib.util.doubleW d2PdTD, org.netlib.util.doubleW dPdT, org.netlib.util.doubleW U, org.netlib.util.doubleW H, org.netlib.util.doubleW S, org.netlib.util.doubleW Cv, org.netlib.util.doubleW Cp, org.netlib.util.doubleW W, org.netlib.util.doubleW G, org.netlib.util.doubleW JT, org.netlib.util.doubleW Kappa, org.netlib.util.doubleW A)
    PropertiesGERG.
    (package private) void
    PseudoCriticalPointLeachman(org.netlib.util.doubleW Tcx, org.netlib.util.doubleW Dcx)
     
    void
    SetupLeachman()
    Initialize the constants and parameters for the default hydrogen type (normal) using the Leachman method.
    void
    SetupLeachman(String hydrogenType)
    SetupGERG.

    Methods inherited from class Object

    clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

  • Field Details

    • hydrogenType

      private String hydrogenType

    • R_L

      double R_L

    • M_L

      double M_L

    • Tc

      double Tc

    • Dc

      double Dc

    • dPdDsave

      double dPdDsave

    • K

      int K

    • a0k

      double[] a0k

    • b0k

      double[] b0k

    • K_normal

      int K_normal

    • a0k_normal

      double[] a0k_normal

    • b0k_normal

      double[] b0k_normal

    • K_para

      int K_para

    • a0k_para

      double[] a0k_para

    • b0k_para

      double[] b0k_para

    • K_ortho

      int K_ortho

    • a0k_ortho

      double[] a0k_ortho

    • b0k_ortho

      double[] b0k_ortho

    • l

      int l

    • m

      int m

    • n

      int n

    • N_i

      double[] N_i

    • t_i

      double[] t_i

    • d_i

      double[] d_i

    • p_i

      double[] p_i

    • phi_i

      double[] phi_i

    • beta_i

      double[] beta_i

    • gamma_i

      double[] gamma_i

    • D_i

      double[] D_i

    • N_i_normal

      double[] N_i_normal

    • t_i_normal

      double[] t_i_normal

    • d_i_normal

      double[] d_i_normal

    • p_i_normal

      double[] p_i_normal

    • phi_i_normal

      double[] phi_i_normal

    • beta_i_normal

      double[] beta_i_normal

    • gamma_i_normal

      double[] gamma_i_normal

    • D_i_normal

      double[] D_i_normal

    • N_i_para

      double[] N_i_para

    • t_i_para

      double[] t_i_para

    • d_i_para

      double[] d_i_para

    • p_i_para

      double[] p_i_para

    • phi_i_para

      double[] phi_i_para

    • beta_i_para

      double[] beta_i_para

    • gamma_i_para

      double[] gamma_i_para

    • D_i_para

      double[] D_i_para

    • N_i_ortho

      double[] N_i_ortho

    • t_i_ortho

      double[] t_i_ortho

    • d_i_ortho

      double[] d_i_ortho

    • p_i_ortho

      double[] p_i_ortho

    • phi_i_ortho

      double[] phi_i_ortho

    • beta_i_ortho

      double[] beta_i_ortho

    • gamma_i_ortho

      double[] gamma_i_ortho

    • D_i_ortho

      double[] D_i_ortho

    • epsilon

      double epsilon

  • Constructor Details

    • Leachman

      public Leachman()

  • Method Details

    • DensityLeachman

      public void DensityLeachman(int iFlag, double T, double P, org.netlib.util.doubleW D, org.netlib.util.intW ierr, org.netlib.util.StringW herr)
      Calculate density as a function of temperature and pressure using the Leachman method. This is an iterative routine that calls PressureLeachman to find the correct state point. Generally, only 6 iterations at most are required. If the iteration fails to converge, the ideal gas density and an error message are returned. No checks are made to determine the phase boundary, which would have guaranteed that the output is in the gas phase (or liquid phase when iFlag=2). It is up to the user to locate the phase boundary, and thus identify the phase of the T and P inputs. If the state point is 2-phase, the output density will represent a metastable state.
      Parameters:
      iFlag - Set to 0 for strict pressure solver in the gas phase without checks (fastest mode, but output state may not be stable single phase). Set to 1 to make checks for possible 2-phase states (result may still not be stable single phase, but many unstable states will be identified). Set to 2 to search for liquid phase (and make the same checks when iFlag=1).
      T - Temperature (K)
      P - Pressure (kPa)
      D - Density (mol/l). For the liquid phase, an initial value can be sent to the routine to avoid a solution in the metastable or gas phases. The initial value should be sent as a negative number.
      ierr - Error number (0 indicates no error)
      herr - Error message if ierr is not equal to zero

    • PressureLeachman

      public void PressureLeachman(double T, double D, org.netlib.util.doubleW P, org.netlib.util.doubleW Z)
      Calculate the pressure and compressibility factor using the Leachman method.
      Parameters:
      T - Temperature (K)
      D - Density (mol/l)
      P - Output parameter for pressure (kPa)
      Z - Output parameter for compressibility factor

    • AlpharLeachman

      void AlpharLeachman(int itau, int idelta, double T, double D, org.netlib.util.doubleW[][] ar)
      Calculate the residual Helmholtz energy and its derivatives with respect to tau and delta. Outputs:
      • ar(0,0) - Residual Helmholtz energy (dimensionless, =a/RT)
      • ar(0,1) - delta*partial (ar)/partial(delta)
      • ar(0,2) - delta^2*partial^2(ar)/partial(delta)^2
      • ar(0,3) - delta^3*partial^3(ar)/partial(delta)^3
      • ar(1,0) - tau*partial (ar)/partial(tau)
      • ar(1,1) - tau*delta*partial^2(ar)/partial(tau)/partial(delta)
      • ar(2,0) - tau^2*partial^2(ar)/partial(tau)^2
      Parameters:
      itau - Order of derivative with respect to tau
      idelta - Order of derivative with respect to delta
      T - Temperature (K)
      D - Density (mol/l)
      ar - Output array for Helmholtz energy and its derivatives

    • Alpha0Leachman

      void Alpha0Leachman(double T, double D, org.netlib.util.doubleW[] a0)
      Calculate the ideal gas Helmholtz energy and its derivatives with respect to tau and delta.
      Parameters:
      T - Temperature (K)
      D - Density (mol/l)
      a0 - Output array for ideal gas Helmholtz energy and its derivatives

    • propertiesLeachman

      public void propertiesLeachman(double T, double D, org.netlib.util.doubleW P, org.netlib.util.doubleW Z, org.netlib.util.doubleW dPdD, org.netlib.util.doubleW d2PdD2, org.netlib.util.doubleW d2PdTD, org.netlib.util.doubleW dPdT, org.netlib.util.doubleW U, org.netlib.util.doubleW H, org.netlib.util.doubleW S, org.netlib.util.doubleW Cv, org.netlib.util.doubleW Cp, org.netlib.util.doubleW W, org.netlib.util.doubleW G, org.netlib.util.doubleW JT, org.netlib.util.doubleW Kappa, org.netlib.util.doubleW A)

      PropertiesGERG.

      Parameters:
      T - a double
      D - a double
      P - a doubleW object
      Z - a doubleW object
      dPdD - a doubleW object
      d2PdD2 - a doubleW object
      d2PdTD - a doubleW object
      dPdT - a doubleW object
      U - a doubleW object
      H - a doubleW object
      S - a doubleW object
      Cv - a doubleW object
      Cp - a doubleW object
      W - a doubleW object
      G - a doubleW object
      JT - a doubleW object
      Kappa - a doubleW object
      A - a doubleW object

    • PseudoCriticalPointLeachman

      void PseudoCriticalPointLeachman(org.netlib.util.doubleW Tcx, org.netlib.util.doubleW Dcx)
      Parameters:
      Tcx - temperature in Kelvin
      Dcx - density

    • SetupLeachman

      public void SetupLeachman(String hydrogenType)

      SetupGERG.

      Parameters:
      hydrogenType - String

    • SetupLeachman

      public void SetupLeachman()
      Initialize the constants and parameters for the default hydrogen type (normal) using the Leachman method.

    • main

      public static void main(String[] args)

      main.

      Parameters:
      args - an array of String objects