diff --git a/@BearImp/calcFilm.m b/@BearImp/calcFilm.m index 799296c8c09dd24c4ba093f44315afc0c4cb38b3..cdb25db26685a8c86e85a37ceb4d09ab4ec1c2a1 100644 --- a/@BearImp/calcFilm.m +++ b/@BearImp/calcFilm.m @@ -70,14 +70,18 @@ H.G = H.alpha_p.*E_red_tmp; H.W = Q_temp./(E_red_tmp .* R_x_temp.^2); switch method - case 'Hamrock/Dowson' + case {'Hamrock/Dowson','Dowson_withoutThermCorr'} k_tmp = [G.k_i;G.k_a]; H.H_0 = 2.69 * H.G.^0.53 .* H.U.^0.67 .* H.W.^-0.067 .* (1-0.61*exp(-0.73.*k_tmp)); H.H_min = 3.63 * H.G.^0.49 .* H.U.^0.68 .* H.W.^-0.073 .* (1- exp(-0.68.*k_tmp)); H.h_0raw = H.H_0 .*R_x_temp; % without correction factors H.h_min = H.H_min.*R_x_temp; - H.L = T.eta_0 .* T.alpha_etaT .* H.u.^2 ./ (4*S.lambda); - H.C_korr = 3.94 ./ (3.94 + H.L.^0.62); + if strcmp(method,'Hamrock/Dowson') + H.L = T.eta_0 .* T.alpha_etaT .* H.u.^2 ./ (4*S.lambda); + H.C_korr = 3.94 ./ (3.94 + H.L.^0.62); + else + H.C_korr = 1; + end H.h_0 = H.h_0raw .*H.C_korr; H.h_minth = H.h_min.*H.C_korr; diff --git a/possibleMethods.m b/possibleMethods.m index 0e9c5caf6e76d8d68cae5b3df6ab0438cc41a793..e96d4efc7c916f03ccbc9e8ab72066fcf25e95f5 100644 --- a/possibleMethods.m +++ b/possibleMethods.m @@ -20,7 +20,7 @@ classdef possibleMethods s = {'static','dynamic'}; end function s = H - s = {'Hamrock/Dowson','Moes'}; + s = {'Hamrock/Dowson','Hamrock/Dowson_withoutThermCorr','Moes'}; end function s = C s.unloadedRE = { 'neglect','stateOfTheArt','Leander_Parallel','Leander_Radial','LeanderSteffen','TobiasSteffen_Kugelfläche','TobiasSteffen_Laufbahnfläche','semianalytisch3D'};