5.0

5.0 | 1 rating Rate this file 35 Downloads (last 30 days) File Size: 4.96 KB File ID: #28135
image thumbnail

STDATMO: Standard Atmosphere Function

by Sky Sartorius

 

05 Jul 2010 (Updated 07 Jul 2011)

Gas properties at altitudes from N-dim input. Non-standard atmospheres capable. Dual unit support.

| Watch this File

File Information
Description

A standard atmosphere function based on the 1976 Standard Atmosphere returns density, speed of sound, temperature, pressure, and viscosity for a given altitude input up to 86km.

This function is designed to be useful for those designing and analyzing aircraft and has the following features that hitherto have not been combined in a single standard atmosphere function:

•Input may be a scalar, vector, matrix, or n-dimensional array of altitudes. The function is vectorized and fast for computing atmospheric conditions at a large number of points simultaneously.

•Temperature offset option for non-standard atmospheres, e.g analyzing hot day aircraft performance.

•Input and output may independently be either SI or imperial units.

•Returns everything needed to easily determine important parameters such as dynamic pressure, Mach number, Reynolds number, stagnation temperature, etc.

•Option for geometric instead of geopotential altitude input, with other choice as optional output.

Examples included in help:

•Create tables of atmospheric properties up to 30000ft for a cold (-15degC), standard, and hot (+15degC) day with columns
[h(ft) Z(ft) rho(slug/ft3) sigma a(ft/s) T(R) P(psf) µ(slug/ft-s) nu(ft2/s)]
using 3-dimensional array inputs.

•Find the SI dynamic pressure, Mach number, Reynolds number, and stagnation temperature of an aircraft flying at flight level FL500 (50000ft) with speed 250m/s and characteristic length of 2m.

•Find atmospheric properties at every 100m of geometric height for an off-standard atmosphere with temperature offset varying +/- 25 degC sinusoidally with a period of 4km.
  

References: ESDU 77022; www.pdas.com/atmos.html

Acknowledgements

The author wishes to acknowledge the following in the creation of this submission:
Standard Atmosphere Calculator, Complete 1976 Standard Atmosphere, ISA chart

MATLAB release MATLAB 7.10 (2010a)
Tags for This File  
Everyone's Tags
Tags I've Applied
Add New Tags Please login to tag files.
Comments and Ratings (2)
05 Oct 2010 Sander Aerts

very clean code with a good help. Haven't yet been able to check the correctness, but it seems OK. The code seems to be implemented way better than the "Complete 1976 Standard Atmosphere" implementation.
Great job

09 May 2011 Sky Sartorius

Right now STDATMO expects the temperature offset to be in degrees Kelvin/Celsius regardless of which length unit is used for altitude input. Keep in mind that dT is an offset, so when converting between Celsius and Fahrenheit, use only the scaling factor (offset C/F = K/R = 5/9).

It was originally thought that non-standard atmospheres would rarely be based on temperature offsets in degrees Fahrenheit, but I would greatly appreciate user feedback: Should the temperature offset for altitude input in meters be in degrees C/K and in degrees F/R for altitude input in feet?

Please login to add a comment or rating.
Updates
07 Jul 2010

Update to references and improved input handling

12 Jul 2010

Changed handling of units to be more intuitive and added feature to allow different units for input and output. Earlier version syntax works as before.
Two examples added to help.

26 Jul 2010

Now substantially faster than earlier versions and ATMOSISA. Cleaned up help and setup for Example 1.

09 Aug 2010

Final output, Hgeop, changed to ZorH, and is geometric altitude for geopotential input, and geopotential for geometric input. Code made a bit more succinct.

27 Aug 2010

Added a section such that troposphere-only simple function calls are a bit faster

27 Aug 2010

A bit quicker for troposphere only problems; fixed small units bug.

23 Dec 2010

Minor changes. Tested in R2010a. Added another example to help.

13 May 2011

fixed error that occurred when using 5 output arguments

07 Jul 2011

Imperial unit temperature offset now in F/R instead of C/K

Tag Activity for this File
Tag Applied By Date/Time
aerospace Sky Sartorius 06 Jul 2010 12:24:07
atmosphere Sky Sartorius 06 Jul 2010 12:24:07
density Sky Sartorius 06 Jul 2010 12:24:08
1976 standard atmosphere Sky Sartorius 06 Jul 2010 12:24:08
standard atmosphere Sky Sartorius 06 Jul 2010 12:24:08
aeronautics Sky Sartorius 06 Jul 2010 12:24:08
international standard atmosphere Sky Sartorius 06 Jul 2010 12:24:08
pressure Sky Sartorius 06 Jul 2010 12:24:08
altitude Sky Sartorius 06 Jul 2010 12:24:08
temperature Sky Sartorius 06 Jul 2010 12:24:08
aircraft design Sky Sartorius 06 Jul 2010 12:24:08
standard atmosphere table Sky Sartorius 06 Jul 2010 12:24:08
isa Sky Sartorius 06 Jul 2010 13:27:47
standard Sky Sartorius 08 Jul 2010 07:47:39
viscosity Sky Sartorius 10 Jul 2010 10:56:51
speed of sound Sky Sartorius 10 Jul 2010 10:56:51
geopotential Sky Sartorius 10 Jul 2010 10:56:51
speed of so Sky Sartorius 09 Aug 2010 10:38:55
atmosphere calculator Sky Sartorius 10 Aug 2010 12:17:58
atmosnonstd Sky Sartorius 26 Aug 2010 06:56:53
atmosisa Sky Sartorius 26 Aug 2010 06:56:53
atmoscira Sky Sartorius 26 Aug 2010 06:56:53
atmoscoesa Sky Sartorius 26 Aug 2010 06:56:53
atmo Sky Sartorius 26 Aug 2010 06:56:53
standard atmosphere table Richard 23 Sep 2011 09:36:41

Contact us at files@mathworks.com