新手求教：双马赫问题编程

zhujiejie

问题描述：      如图所示，假设初始时刻，Mach=10的正激波向平板运动。激波与平板之间成60°夹角。 试计算流场的演化，画出t=0.2时刻流场的密度及压力分布，并进行流场分析。
主要疑惑：
1，Steger-Warming 分裂
对二维无量纲Euler方程组进行分裂，但方程形式负责，求解特征根、特征向量困难，如何解决？
2，激波运动：
这一问题是非定常的，可否用时间推进法？边界条件如何应用？
望赐教！

通流

（1）线性化以后无所谓简单复杂。
（2）时间推进法本来就是非定常的。边界条件按非定常线性化的方程给。就是根据特征值特征向量确定那些值从计算域中取，那些是外面传进来的。

zhujiejie

谢谢指教！

zhujiejie

你好！
        我对二维无量纲Euler方程组，通过MATLAB符号运算得到F和G对U的雅克比矩阵A.B，但用MATLAB符号运算对A、B求特征根时得到的结果太过冗长，明显错误。难以进一步分裂。
    请问如何解决Steger-Warming 分裂这个问题？是否有已经推导得到的结果？
    期待你的回复！谢谢!