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 "cooldt.h"

/*
 * Function to export
 */
double __stdcall cooldt(double Tinitial, double Tambient, double Tfinal, double Hout, int NbLayer,
                        const double * Diam, const double * Conduc, const double * Dens, const double * Cp,
                        char * ErrorMessage)
{
    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 = NbLayer ;
    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
     */
    FILE *fp = fopen( "NUL", "w") ;
    if ( CheckInputData(fp, ErrorMessage, PP) != 0)
    {
        fclose(fp) ;
        FreePipeProblem(PP) ;
        free(PP) ;
        return NAN ;
    }

    /*
     * 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 ;
}