For more information, contact Joel W. Sills, Langley Research Center, joel.w.sills@nasa.gov
NASA’s Chemical Equilibrium with Applications (CEA) code is a foundational tool for propulsion system analysis. It provides equilibrium chemistry, rocket performance, shock, and detonation calculations used across NASA and the broader aerospace community. NASA Engineering and Safety Center (NESC) Activity TI-22-01730 modernized the legacy CEA2 Fortran code into CEA v3, a Fortran 2008, object-oriented software package with expanded interface support, updated thermochemical data, improved maintainability, and substantially improved workflow integration. The modernized code preserves backward compatibility with legacy CEA input workflows while enabling direct use from modern analysis environments, including Python, C, MATLAB, and automated design studies.
Why Complex Modes Can Appear at Low Damping
Classical modal analysis often treats lightly damped structures as nearly proportional systems with real mode shapes. That interpretation can become misleading when two modes are closely spaced. In modal coordinates, nonproportional damping introduces an off-diagonal coupling term ε. The relevant question is therefore not simply whether damping is small, but whether the coupling is large relative to the baseline half-frequency separation δ=(ω2 -ω1 )/2.
Reduced Perturbation Framework
For a closely spaced pair, define the central eigenvalue λ0 ≈-ζω +iω and write the coupled eigenvalues as λ1,2=λ0 ±s. Under lightdamping and close-mode asymptotic assumptions, the local quartic eigenvalue problem reduces to s²=ε²/4-δ². This single relation determines the direction of eigenvalue separation in the complex plane, as summarized in Figure. 1. For 0<ε 2δ, s = ±β and the split is predominantly horizontal (decay rate). As ε increases further, β increases and the two decay rates separate more strongly while their frequencies remain approximately equal.
Engineering Interpretation and Recommended Practice
Complexity in an identified mode shape should not automatically be interpreted as evidence of heavy damping, nonlinearity, or poor test data. For a close modal pair, examine the complex poles together with mode-shape phase. Frequency convergence accompanied by near-quadrature modal components is consistent with the imaginary-s regime; nearly equal frequencies accompanied by decay-rate separation is consistent with the real-s regime. The framework also provides a practical testanalysis bridge: s=(λ1 -λ2 )/2 can be obtained from identified poles, while δ can be estimated from the uncoupled or pretest analytical modal spacing. These quantities can then be used to estimate the effective interaction strength and determine proximity to the transition boundary.
References
1. Majed, A., Henkel, E. E., and Sills, J. W., Emergence of Complex Modes in Lightly Damped Structural Systems: A Reduced Perturbation Framework, NASA/TM-20260007553, August 2026.