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., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
83 **********************************************************************/
92 if (Simtime >= iceTime)
97 // slowly increase icing severity over period of transientTime
98 if (Simtime < (iceTime + transientTime))
100 slope = eta_ice_final / transientTime;
101 eta_ice = slope * (Simtime - iceTime);
105 eta_ice = eta_ice_final;
112 double uiuc_ice_filter( double Ca_clean, double kCa )
116 //cout << "Ice Model Engaged" << endl;
118 Ca_iced = Ca_clean * (1 + kCa * eta_ice);