File Coverage

blib/lib/Physics/Psychrometry.pm

Criterion	Covered	Total	%
statement	64	65	98.4
branch	4	6	66.6
condition	2	3	66.6
subroutine	17	17	100.0
pod	12	14	85.7
total	99	105	94.2

# Psychrometry.pm
# Calculate psychrometric measures in moist air
package Physics::Psychrometry;

use 5.008007;
use strict;
use warnings;

require Exporter;

our @ISA = qw(Exporter);

# Items to export into callers namespace by default. Note: do not export
# names by default without a very good reason. Use EXPORT_OK instead.
# Do not simply export all your public functions/methods/constants.

# This allows declaration	use Physics::Psychrometry ':all';
# If you do not need this, moving things directly into @EXPORT or @EXPORT_OK
# will save memory.
our %EXPORT_TAGS = ( 'all' => [ qw(

) ] );

our @EXPORT_OK = ( @{ $EXPORT_TAGS{'all'} } );

our @EXPORT = qw(

);

our $VERSION = '1.2';

# based on http://www.bae.uky.edu/gates/psych/PTM/dbrh2al.c
#The dry-bulb temperature (C)       =? 24.6
#The relative humidity (%)          =? 68
#The air pressure (kPa)             =? 102
# e is 2.104209
#The dewpoint temperature (C)       = 18.324036
#The wet-bulb temperature (C)       = 20.326262
#The humidity ratio                 = 0.013101
#The enthalpy (kJ/kg)               = 58.095979
#The vapor pressure (kPa)           = 2.104209
#The degree of saturation (%)       = 67.325951
#The saturated vapor pressure (kPa) = 3.094426
#The specific air volume (m3/kg)    = 0.855627
#The density of moist air (kg/m^3)  = 1.184046
#The saturated humidity ratio       = 0.019460
#The specific humidity (ratio)      = 0.012932
#The absolute humidity (kg/m3)      = 0.015312

my $C14	= 6.54;
my $C15	= 14.526;
my $C16	= 0.7389;
my $C17	= 0.09486;
my $C18	= 0.4569;

sub dbdp2wb
{
    my ($db, $dp, $p) = @_;

    my $i = 0;
    my $e = t2es($dp);
    my $w = e2w($e, $p);
    my $wb1 = $dp;
    my $wb2 = $db;
    my $twb = $wb1 + ($wb2 - $wb1) / 2;
    while ((($wb2 - $wb1) > 0.001) && ($i < 1000))
    {
	my $tw = dbwb2w($db, $twb, $p);
	if($tw > $w)
	{		# /* overestimate the wetblb temperature */
	    $wb2 = $twb; #		/* the wb should be between wb1,twb */
	    $twb = $wb1 + ($wb2 - $wb1) / 2;
	}
	else
	{		#	/* underestimate the wb */
	    $wb1 = $twb;
	    $twb = $wb1 + ($wb2 - $wb1) / 2;
	}
    }
    my $wb = $twb;
    return $wb;
}

sub dbw2h
{
    my ($db, $w) = @_;
    return 1.006 * $db + $w * (2501.0 + 1.805 * $db);
}

sub e2w
{
    my ($e, $p) = @_;

    return 0.62198 * $e / ($p - $e);
}

sub dbwb2w
{
    my ($db, $wb, $p) = @_;


    my $estar = t2es($wb);
    my $wstar = e2w($estar, $p);
    my $t1 = (2501 - 2.381 * $wb) * $wstar -($db - $wb);
    my $t2 = 2501 + 1.805 * $db - 4.186 * $wb;
    return $t1 / $t2;
}

sub dbrh2e
{
    my ($db, $rh) = @_;

    return t2es($db) * $rh;
}

sub rhws2ds
{
    my ($rh, $ws) = @_;

    return $rh / (1 + (1 - $rh) * $ws / 0.62198);
}

sub t2es
{
    my ($t) = @_;

    $t += 273.15;
    if($t > 273.15)
    {
	return exp(-5800.2206 / $t + 1.3914993 - .048640239 * $t + (.41764768e-4) * pow($t, 2.0) - (.14452093e-7) * pow($t, 3.0) + 6.5459673 * log($t)) / 1000.0;
    }
    else
    {
      return exp(-5674.5359 / $t + 6.3925247 - (.9677843e-2) * $t + (.62215701e-6) * pow($t, 2.0) + (.20747825e-8) * pow($t, 3.0) - (.9484024e-12) * pow($t, 4.0) + 4.1635019 * log($t)) / 1000.0;
  }
}

sub e2dp
{
    my ($e) = @_;

    my $af = log($e);
    my $s1 = pow($af, 2.0);
    my $s2 = pow($af, 3.0);
    my $s3 = pow($e, 0.1984);
    my $t1 = $C14 + $C15 * $af + $C16 * $s1 + $C17 * $s2 + $C18 * $s3;
    my $t2 = 6.09 + 12.608 * $af + 0.4959 * $s1;
    return $t1 > 0 ? $t1 : $t2;
}

sub pow
{
    my ($n, $p) = @_;

    return $n ** $p;
}

sub dbw2v
{
    my ($db, $w, $p) = @_;

    $db += 273.15;
    return (287.055 * $db * ( 1 + 1.6078 * $w)) / ($p * 1000.0);
}

sub wv2da
{
    my ($w, $v) = @_;

    return (1 + $w) / $v;
}

sub es2ws
{
    my ($es, $p) = @_;

    return e2w($es, $p);
}

sub w2q
{
    my ($w) = @_;

    return $w / (1 + $w);
}

sub wv2X
{
    my ($w, $v) = @_;

    return $w / $v;
}

1;
__END__