 # Delivered to the U.S. Government in Accordance with Contract FA8650-10-C-7080.
 #
 # EXPORT CONTROLLED
 # ECCN#   EAR99  
 # Warning:  This document contains technical data which is controlled under the
 # Bureau of Industry & Security’s Export Administration Regulations and may not be
 # exported or re-exported without authorization by the Bureau of Industry & Security.
 # Diversion contrary to US law is prohibited.
 #
 # (c) 2011 United Technologies Corp.
 # Developed with the sponsorship of the Defense Advanced Research Projects Agency (DARPA).
 #
 # Permission is hereby granted, free of charge, to any person obtaining a copy of this data,
 # including any software or models in source or binary form, as well as any drawings, 
 # specifications, and documentation (collectively "the Data"), to deal in the Data without 
 # restriction, including without limitation the rights to use, copy, modify, merge, publish,
 # distribute, sublicense, and/or sell copies of the Data, and to permit persons to whom the Data
 # is furnished to do so, subject to the following conditions:
 #  
 # The above copyright notice and this permission notice shall be included in all copies or
 # substantial portions of the Data.
 #  
 # THE DATA IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING
 # BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
 # NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS, SPONSORS, DEVELOPERS, CONTRIBUTORS, OR
 # COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
 # OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE DATA OR THE
 # USE OR OTHER DEALINGS IN THE DATA.

######## Single_Front_Prop$:  			
param Np$;
param Transfer_Efficiency$;
var rp$>=1.05, <=1.3;
var bpr>=2.0, <=20;
var ma >= 0.01; 
var To$>=1;
var Po$>=1;
var cp$>=900,<=1400;
var gam$>=1.2,<=1.5;
var M$ >=0;
var cost$>=0;

var V_prop>=1;
var A_prop>=0.1;

	
subject to conPo$:Po$ = Pac;
subject to conTo$:To$ = Tac;

subject to conV_prop$: V_prop = sqrt(2*Pac*(1+(gam1-1)/2*Mac^2)^(gam1/(gam1-1))*(rp$-1)/(Pac/(R*Tac))+Vac^2);
subject to conA_prop$: A_prop = ma_core*(1+bpr)/((Pac/(R*Tac))*0.5*(Vac+V_prop));
subject to conma$: ma = ma_core*(1+bpr);


subject to comM$: M$ = 0.5*(Vac+V_prop)/(sqrt(gam1*R*Tac));
subject to conThrust$: Thto = ma_core*(1+bpr)*(V_prop-Vac)+ma_core*(V_core-Vac);
subject to conW:W2=0.5*ma_core*(1+bpr)*(V_prop^2-Vac^2)/(Transfer_Efficiency$*Np$);			 


subject to concp$:cp$ = cp1;
subject to congam$: gam$ = gam1;
subject to conchangeV_prop$: V_prop-Vac>=10;

subject to concost$: cost$=bpr*ma_core*Np$*700;


