1 /**********************************************************************
5 ----------------------------------------------------------------------
7 DESCRIPTION: reads in clean coefficient and icing severity
8 parameters and returns iced coefficient
10 ----------------------------------------------------------------------
14 ----------------------------------------------------------------------
18 ----------------------------------------------------------------------
20 HISTORY: 02/22/2000 initial release
21 04/25/2000 (JS) added uiuc_ice_eta function
22 (removed from uiuc_coefficients)
24 ----------------------------------------------------------------------
26 AUTHOR(S): Jeff Scott <jscott@mail.com>
28 ----------------------------------------------------------------------
32 ----------------------------------------------------------------------
37 -final icing severity (eta_ice_final)
40 -clean aero coefficient
41 -icing parameter for that coefficient (kC)
42 -icing severity (eta_ice)
44 ----------------------------------------------------------------------
46 OUTPUTS: uiuc_ice_eta:
47 -icing severity (eta_ice)
50 -iced aero coefficient
52 ----------------------------------------------------------------------
54 CALLED BY: uiuc_coefficients
62 ----------------------------------------------------------------------
66 ----------------------------------------------------------------------
68 COPYRIGHT: (C) 2000 by Michael Selig
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.
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.
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., 59 Temple Place - Suite 330, Boston, MA 02111-1307,
82 USA or view http://www.gnu.org/copyleft/gpl.html.
84 **********************************************************************/
93 if (Simtime >= iceTime)
98 // slowly increase icing severity over period of transientTime
99 if (Simtime < (iceTime + transientTime))
101 slope = eta_ice_final / transientTime;
102 eta_ice = slope * (Simtime - iceTime);
106 eta_ice = eta_ice_final;
113 double uiuc_ice_filter( double Ca_clean, double kCa )
117 //cout << "Ice Model Engaged" << endl;
119 Ca_iced = Ca_clean * (1 + kCa * eta_ice);