ME 322 MATLAB Tools
This page provides links to MATLAB mfiles for use with ME 322. The mfiles are organized into two Toolboxes.
The mfiles may be downloaded as a single zip archive or as individuals files from the following list.
| m-file name | Description | ||
|---|---|---|---|
| contractionLoss.m | computes the minor loss coefficient for a sudden contraction. | ||
| demoFlowRate.m | demonstrates how to solve a flow rate problem. | ||
| demoHeadLoss.m | demonstrates how to solve a head loss problem. | ||
| demoPumpCurve.m | demonstrates how to construct a polynomial curve fit to tabulated pump curve data. | ||
| demoSysBalance.m | demonstrates how to find the balance point for a pipe system and a pump, where the pump curve is defined as a polynomial. | ||
| demoSystemCurve.m | demonstrates how to generate a system curve. | ||
| entranceLength.m | plots the entrance length correlations, Le/D = fcn(Re), for laminar and turbulent flow. | ||
| makePumpCurve.m | uses least squares method to find a polynomial curve fit to tabulated pump curve data. | ||
| moody.m | evaluates the friction factor by solving the Colebrook equation (Moody Chart). | ||
| pipeFlowSolve.m | uses root-finding to solve a flow rate problem. | ||
| pipeHeadBal.m | evaluates the imbalance in the energy equation for a pipe run. The energy equation is put in the form f(Q) = 0 for use with root-finding algorithms. pipeHeadBal returns the value of f(Q) for a given Q. When the correct Q is supplied as input, pipeHeadBal will return f(Q)=0. pipeHeadBal is used by pipeFlowSolve. | ||
| pipeLoss.m | computes viscous and minor losses for a pipe run. |
The mfiles may be downloaded as a single zip archive or as individuals files from the following list. A short PDF file explaining how to use these mfiles is also available.
| m-file name | Description | ||
|---|---|---|---|
| aasmaResidual.m | computes residual of relationship between A/A* and Ma. Used by the built-in fzero root-finding routine in the solution for computing A/A* as a function of Ma. | ||
| isenAAs.m | computes area ratio A/A* as a function of Ma and k. | ||
| isenMaaas.m | computes Ma as function of area ratio A/A* and k. The appropriate branch --- subsonic or supersonic --- is chosen with an optional input parameter. | ||
| isenMapp0.m | computes Ma as function of pressure ratio p/p0 and k. | ||
| isenMarr0.m | computes Ma as function of density ratio r/r0 and k. | ||
| isenMaTT0.m | computes Ma as function of temperature ratio T/T0 and k. | ||
| isenpp0.m | computes pressure ratio p/p0 as function of Ma and k. | ||
| isenrr0.m | computes density ratio r/r0 as function of Ma and k. | ||
| isenTT0.m | computes temperature ratio T/T0 as function of Ma and k. | ||
| MYO_11_38.m | End of chapter exercise number 11.38 from Fundamentals of Fluid Mechanics by Munson, Young and Okiishi, 3rd edition, 1998, Wiley, p. 774. | ||
| testIsenProps.m | exercises all of the stagnation property routines by performing calculations in both directions, e.g. T/T0 = f(Ma,k) and Ma = g(T/T0,k). | ||
| white_E9_3.m | Exercise 9.3 from Fluid Mechanics by F.M. White, 4th edition, 1999, McGraw-Hill, pp. 582-583. |
ME 322 Home Page Gerry's Other Classes Gerry's Home Page