diff --git a/contributing/contributing.html b/contributing/contributing.html index 59e88b2cc..4f1882361 100644 --- a/contributing/contributing.html +++ b/contributing/contributing.html @@ -1055,21 +1055,21 @@

Notebook Creation#

Both sets of equations of motion have common aerodynamic models.

First, the lift coefficient is computed using a model which assumes linearity in lift wrt the angle of attack.

-
-(21)#\[\begin{align} +
+(21)#\[\begin{align} C_L &= C_{L0} + \frac{\alpha}{\alpha_{max}} \left(C_{L-max} - C_{L0}\right) \end{align}\]

Next, the drag-due-to-lift factor is computed (Equations 12.48 and 12.61 in Raymer[@raymer2012aircraft]).

-
-(22)#\[\begin{align} +
+(22)#\[\begin{align} K_{nom} &= \frac{1}{ \pi AR e} \\ b &= \frac{span}{2} \\ K &= 33 K_{nom} \frac{ \left(\frac{h + h_w}{b} \right) ^{\frac{3}{2}}}{1.0 + 33 \left( \frac{h + h_w}{b}\right) ^{\frac{3}{2}}} @@ -956,8 +956,8 @@

The Aerodynamic Model

Finally, the lift and drag are computed after computing the dynamic pressure.

-
-(23)#\[\begin{align} +
+(23)#\[\begin{align} q &= 0.5 \rho v^2 \\ L &= q S C_L \\ D &= q S \left( C_{D0} + K C_{L}^2 \right) @@ -968,8 +968,8 @@

Stall Speed\(\frac{v}{v_{stall}}\)). This constraint is used to trigger the beginning of rotation and is used as a boundary constraint at the end of the initial climb. Stall speed is given by Equation 5.6 in Raymer[@raymer2012aircraft].

-
-(24)#\[\begin{align} +
+(24)#\[\begin{align} W &= m g \\ v_{stall} &= \sqrt{\frac{2 * W}{\rho S C_{L-max}}} \end{align}\]
@@ -977,8 +977,8 @@

Stall Speed

Runway Equations of Motion#

The runway equations of motion are used to integrate range and speed as the vehicle rolls along the runway.

-
-(25)#\[\begin{align} +
+(25)#\[\begin{align} F_r &= mg - L \cos \alpha - T \sin \alpha \\ \dot{v} &= \frac{T \cos \alpha - D - F_r \mu_r}{m} \\ \dot{r} &= v @@ -1005,8 +1005,8 @@

Runway Equations of Motion

Flight Equations of Motion#

The flight equations of motion include two additional state variables: the flight-path angle (\(\gamma\)) and altitude (\(h\)).

-
-(26)#\[\begin{align} +
+(26)#\[\begin{align} \dot{v} &= \frac{T}{m} \cos \alpha - \frac{D}{m} - g \sin \gamma \\ \dot{\gamma} &= \frac{T}{m v} \sin \alpha + \frac{L}{m v} - \frac{g \cos \gamma}{v} \\ \dot{h} &= v \sin \gamma \\ @@ -1531,7 +1531,7 @@

Building and running the problem -
<dymos.phase.phase.Phase at 0x7f0e800f6b50>
+
<dymos.phase.phase.Phase at 0x7f9769437390>
@@ -1727,8 +1727,8 @@

Building and running the problem
Coloring for 'traj.phases.v1_to_vr.rhs_all' (class BalancedFieldODEComp)
+
+Coloring for 'traj.phases.v1_to_vr.rhs_all' (class BalancedFieldODEComp)
 
 Jacobian shape: (132, 60)  (10.15% nonzero)
 FWD solves: 16   REV solves: 0
 Total colors vs. total size: 16 vs 60  (73.33% improvement)
 
 Sparsity computed using tolerance: 1e-25
-Time to compute sparsity:   0.0157 sec
-Time to compute coloring:   0.0098 sec
+Time to compute sparsity:   0.0156 sec
+Time to compute coloring:   0.0097 sec
 Memory to compute coloring:   0.0000 MB
 
 Coloring for 'traj.phases.rotate.rhs_all' (class BalancedFieldODEComp)
@@ -1761,8 +1760,8 @@ 
Building and running the problem
Full total jacobian for problem 'problem' was computed 3 times, taking 0.9825835510000047 seconds.
+
Full total jacobian for problem 'problem' was computed 3 times, taking 1.009388161000004 seconds.
 Total jacobian shape: (178, 166) 
 
 
@@ -1786,10 +1785,10 @@ 
Building and running the problem
/usr/share/miniconda/envs/test/lib/python3.11/site-packages/openmdao/core/total_jac.py:1652: DerivativesWarning:Constraints or objectives [('traj.phases.climb->path_constraint->gam', inds=[(0, 0)])] cannot be impacted by the design variables of the problem.
@@ -1810,11 +1809,11 @@ 
Building and running the problem
/usr/share/miniconda/envs/test/lib/python3.11/site-packages/openmdao/core/system.py:1701: DerivativesWarning:'traj.phases.rto.rhs_all' <class ExplicitFuncComp>: No partials found but coloring was requested.  Declaring ALL partials as dense (method='jax')
@@ -1270,8 +1270,8 @@ 
Building and running the problem
Full total jacobian for problem 'problem' was computed 3 times, taking 3.2737298330000044 seconds.
+
Full total jacobian for problem 'problem' was computed 3 times, taking 3.247299850000047 seconds.
 Total jacobian shape: (178, 166) 
 
 
@@ -1329,10 +1329,10 @@ 
Building and running the problem
/usr/share/miniconda/envs/test/lib/python3.11/site-packages/openmdao/core/total_jac.py:1652: DerivativesWarning:Constraints or objectives [('traj.phases.climb->path_constraint->gam', inds=[(0, 0)])] cannot be impacted by the design variables of the problem.
diff --git a/examples/brachistochrone/brachistochrone.html b/examples/brachistochrone/brachistochrone.html
index d56c2360c..16829554a 100644
--- a/examples/brachistochrone/brachistochrone.html
+++ b/examples/brachistochrone/brachistochrone.html
@@ -866,8 +866,8 @@ 
State variables
 
System dynamics#

 
From the free-body diagram above, the evolution of the state variables is given by the following ordinary differential equations (ODE).

-

-(13)#\[\begin{align}
+

+(13)#\[\begin{align}
     \frac{d x}{d t} &= v \sin(\theta) \\
     \frac{d y}{d t} &= -v \cos(\theta) \\
     \frac{d v}{d t} &= g \cos(\theta)
@@ -885,8 +885,8 @@ 
The initial and final conditions
-(14)#\[\begin{align}
+

+(14)#\[\begin{align}
     x_0 &= 0 \\
     y_0 &= 10 \\
     v_0 &= 0
@@ -894,8 +894,8 @@ 
The initial and final conditions
-(15)#\[\begin{align}
+

+(15)#\[\begin{align}
     x_f &= 10 \\
     y_f &= 5 \\
     v_f &= \mathrm{free}
@@ -1061,27 +1061,27 @@ 
Testing the ODE
Full total jacobian for problem 'problem2' was computed 3 times, taking 0.021816589000025033 seconds.
+
Full total jacobian for problem 'problem2' was computed 3 times, taking 0.021945925000011357 seconds.
 Total jacobian shape: (40, 50) 
 
 
@@ -1188,10 +1188,10 @@ 
Solving the problem with Legendre-Gauss-Lobatto collocation in Dymos

 

 
Optimization terminated successfully    (Exit mode 0)
@@ -1337,7 +1337,7 @@ 
Solving the problem with Legendre-Gauss-Lobatto collocation in Dymos

 

-
Full total jacobian for problem 'problem5' was computed 3 times, taking 0.014237258000036945 seconds.
+
Full total jacobian for problem 'problem5' was computed 3 times, taking 0.014270430999999917 seconds.
 Total jacobian shape: (12, 51) 
 
 No coloring was computed successfully.
diff --git a/examples/brachistochrone/brachistochrone_tandem_phases.html b/examples/brachistochrone/brachistochrone_tandem_phases.html
index 97aa0b528..112fa540c 100644
--- a/examples/brachistochrone/brachistochrone_tandem_phases.html
+++ b/examples/brachistochrone/brachistochrone_tandem_phases.html
@@ -994,7 +994,7 @@ 
Setup and run
Model viewer data has already been recorded for Driver.
 

 

-
Full total jacobian for problem 'problem' was computed 3 times, taking 0.1541596110000114 seconds.
+
Full total jacobian for problem 'problem' was computed 3 times, taking 0.15848857100002078 seconds.
 Total jacobian shape: (223, 287) 
 

 

@@ -1003,10 +1003,10 @@ 
Setup and run
/usr/share/miniconda/envs/test/lib/python3.11/site-packages/openmdao/visualization/opt_report/opt_report.py:625: UserWarning: Attempting to set identical low and high ylims makes transformation singular; automatically expanding.
@@ -1854,10 +1854,10 @@ 
Plots#
     
   
   
-    

+    

   
-    
     
   
   
-    

+    

   
-    
     
   
   
-    

+    

   
-    
     
   
   
-    

+    

   
-