# After i generate random numbers of x, how can i use this random x values to the rest of the code one by one. Because after that, I want to plot a friction factor-x number grapt with this different numbers

1 view (last 30 days)
isilay aydin on 18 May 2020
Edited: Stephen Cobeldick on 19 May 2020
x = 0.05 + (0.1-0.05).*rand(10,1);
Perimeter = (pi*HydraulicDiameter);
NumberOfCoolantChannels = 3;
Npass = 1; ...Number of Passes
% Ac = (pi*(HydraulicDiameter^2)/(4)); ...Single channel flow cross section
Ac = 1.571e-5;
Atotal = NumberOfCoolantChannels*Ac; ... Total Flow Cross Section
RibHeighHydraulicDiameterRatio = x;
FrictionFactor2 = 1.447*(1.14-2*log10(x))^-2; %transverse ribs

Geoff Hayes on 18 May 2020
isilay - since x is an array then you can use element-wise power to calculate the friction factor. Just change your line of code to
FrictionFactor2 = 1.447*(1.14-2*log10(x)).^-2; %transverse ribs
where a '.' is placed in front of the exponent operator.

isilay aydin on 19 May 2020
this is my final version still idk if it is true
x = 0.05 + (0.1-0.05).*rand(10,1);
Perimeter = (pi*HydraulicDiameter);
NumberOfCoolantChannels = 3;
Npass = 1; ...Number of Passes
% Ac = (pi*(HydraulicDiameter^2)/(4)); ...Single channel flow cross section
Ac = 1.571e-5;
Atotal = NumberOfCoolantChannels*Ac; ... Total Flow Cross Section
RibHeighHydraulicDiameterRatio = x;
FrictionFactor2 = 1.447*(1.14-2*log10(x)).^-2; %transverse ribs
OverallPressureDrop = 8300;
ChannelCoolantMassFlow = CoolantDensity*Ac*BulkVelocity;
Re = (ChannelCoolantMassFlow/Ac)*(HydraulicDiameter/CoolantViscosity);
Nu = 0.15*Re.^0.7; ...based on the Hydraulic Diameter
InternalHeatTransferCoefficient = Nu*(CoolantConductivity/HydraulicDiameter);
ConvectionCoolingEfficiency = 1 - exp((-TechnologyFactor)/(MassFlowFunction));
OverallEffectinevess = (MassFlowFunction*ConvectionCoolingEfficiency)/(1 + MassFlowFunction*ConvectionCoolingEfficiency);
MetalTemperature = (ExternalGasTemperature) - OverallEffectinevess*(ExternalGasTemperature-CoolantInletTemperature);
Geoff Hayes on 19 May 2020
isilay - I think that you need to step through each line of code (using the MATLAB debugger) and verify that all variables are as you expect and that, given their dimensions, the code is performing the correct operation.
Stephen Cobeldick on 19 May 2020
"...i thought that i should put dot every variables which affect the X and Friction factor"
But those are not the only non-scalar variables in your code. These ones are too (or probably should be):
BulkVelocity
ChannelCoolantMassFlow
PercentageCoolantBleed
Re
Nu
InternalHeatTransferCoefficient
TechnologyFactor
MassFlowFunction
etc. etc.
Why? Lets have a look at the definition of BulkVelocity: