src/FDM/UIUCModel/uiuc_ice.cpp

   1 /**********************************************************************
   2
   3  FILENAME:     uiuc_ice.cpp
   4
   5 ----------------------------------------------------------------------
   6
   7  DESCRIPTION:  reads in clean coefficient and icing severity
   8                parameters and returns iced coefficient
   9
  10 ----------------------------------------------------------------------
  11
  12  STATUS:       alpha version
  13
  14 ----------------------------------------------------------------------
  15
  16  REFERENCES:
  17
  18 ----------------------------------------------------------------------
  19
  20  HISTORY:      02/22/2000   initial release
  21                04/25/2000   (JS) added uiuc_ice_eta function
  22                             (removed from uiuc_coefficients)
  23
  24 ----------------------------------------------------------------------
  25
  26  AUTHOR(S):    Jeff Scott         <jscott@mail.com>
  27
  28 ----------------------------------------------------------------------
  29
  30  VARIABLES:
  31
  32 ----------------------------------------------------------------------
  33
  34  INPUTS:       uiuc_ice_eta:
  35                -Simtime
  36                -icing times
  37                -final icing severity (eta_ice_final)
  38
  39                uiuc_ice_filter:
  40                -clean aero coefficient
  41                -icing parameter for that coefficient (kC)
  42                -icing severity (eta_ice)
  43
  44 ----------------------------------------------------------------------
  45
  46  OUTPUTS:      uiuc_ice_eta:
  47                -icing severity (eta_ice)
  48
  49                uiuc_ice_filter:
  50                -iced aero coefficient
  51
  52 ----------------------------------------------------------------------
  53
  54  CALLED BY:    uiuc_coefficients
  55                uiuc_coef_drag
  56                uiuc_coef_lift
  57                uiuc_coef_pitch
  58                uiuc_coef_sideforce
  59                uiuc_coef_roll
  60                uiuc_coef_yaw
  61
  62 ----------------------------------------------------------------------
  63
  64  CALLS TO:     none
  65
  66 ----------------------------------------------------------------------
  67
  68  COPYRIGHT:    (C) 2000 by Michael Selig
  69
  70  This program is free software; you can redistribute it and/or
  71  modify it under the terms of the GNU General Public License
  72  as published by the Free Software Foundation.
  73
  74  This program is distributed in the hope that it will be useful,
  75  but WITHOUT ANY WARRANTY; without even the implied warranty of
  76  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  77  GNU General Public License for more details.
  78
  79  You should have received a copy of the GNU General Public License
  80  along with this program; if not, write to the Free Software
  81  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
  82
  83 **********************************************************************/
  84
  85 #include "uiuc_ice.h"
  86
  87
  88 void uiuc_ice_eta()
  89 {
  90   double slope = 0;
  91
  92   if (Simtime >= iceTime)
  93     {
  94       // set ice_on flag
  95       ice_on = true;
  96
  97       // slowly increase icing severity over period of transientTime
  98       if (Simtime < (iceTime + transientTime))
  99         {
 100           slope = eta_ice_final / transientTime;
 101           eta_ice = slope * (Simtime - iceTime);
 102         }
 103       else
 104         {
 105           eta_ice = eta_ice_final;
 106         }
 107     }
 108   return;
 109 }
 110
 111
 112 double uiuc_ice_filter( double Ca_clean, double kCa )
 113 {
 114   double Ca_iced = 0;
 115
 116   //cout << "Ice Model Engaged" << endl;
 117
 118   Ca_iced = Ca_clean * (1 + kCa * eta_ice);
 119
 120   return Ca_iced;
 121 }
 122
 123 // end uiuc_ice.cpp