Cooldt Code

/*
 * Cooldt cool-down time calculator for multi-layers pipe
 *
 * Copyright 2010, 2013, 2014, Benjamin DEGLO DE BESSES. All rights reserved.
 *
 * Redistribution and use in source and binary forms, with or without modification, are
 * permitted provided that the following conditions are met:
 *
 *    1. Redistributions of source code must retain the above copyright notice, this list of
 *       conditions and the following disclaimer.
 *
 *    2. Redistributions in binary form must reproduce the above copyright notice, this list
 *       of conditions and the following disclaimer in the documentation and/or other materials
 *       provided with the distribution.
 *
 * THIS SOFTWARE IS PROVIDED BY Benjamin DEGLO DE BESSES "AS IS" AND ANY EXPRESS OR IMPLIED
 * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND
 * FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL Benjamin DEGLO DE BESSES OR
 * CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
 * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON
 * ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
 * NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF
 * ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 *
 */

#include "cppinterface.h"

/*
 * C function declaration
 */
extern "C"
{
#include "cooldt.h"
}

double CoolDownTime(double Tinitial, double Tambient, double Tfinal, double Hout, const MyArray& Diam,
                    const MyArray& Conduc, const MyArray& Dens, const MyArray& Cp)
{
    /*
     * Check if Diam, Conduc, Dens & Cp have the same size
     */
    if ( ( Diam.size() != Conduc.size() ) || ( Diam.size() != Dens.size() ) || ( Diam.size() != Cp.size() ) )
    {
        THROW_XLW("#" << C_ERROR_NB_NOT_EQUALS) ;
    }

    PipeProblem* PP = (PipeProblem*)malloc(sizeof(PipeProblem)) ;

    /*
     * Copy problem parameters
     */
    PP->Tinitial = Tinitial ;
    PP->Tambient = Tambient ;
    PP->Tfinal = Tfinal ;
    PP->Hout = Hout ;

    /*
     * Size the PP structures pointers
     */
    PP->N = Diam.size() - 1 ;
    SizePipeProblem(PP) ;

    /*
     * Copy layers parameters
     */
    int i = PP->N +1 ;
    while(i--)
    {
        PP->r[i] = 0.5 * Diam[i] ;
        PP->k[i] = Conduc[i] ;
        PP->d[i] = Dens[i] ;
        PP->Cp[i] = Cp[i] ;
    }

    /*
     * Check if input data values are OK
     */
    char ErrorMessage[CharBufferSize] ;
    FILE *fp = fopen( "NUL", "w") ;
    if (  CheckInputData(fp, ErrorMessage, PP) != 0)
    {
        fclose(fp) ;
        FreePipeProblem(PP) ;
        free(PP) ;

        THROW_XLW("#" << ErrorMessage) ; /* return error messages */
    }

    /*
     * The data are OK: Calculate CDT
     * Calculate the time step and cells number to reach the accuracy
     */
    aaControl(fp, PP) ;

    /*
     * Calculate the solution
     * then convert it in hours
     */
    double cdt_out = SolutionSimpleOutput(PP) / 3600.0 ;

    /*
     * Free data structures and pointers
     */
    fclose(fp) ;
    FreePipeProblem(PP) ;
    free(PP) ;

    /*
     * return the cdt
     */
    return cdt_out ;
}