第6次作业leveque书习题41abadeterminethematricesand.docVIP

第6次作业 Leveque书习题 4.1 （a）（b） （a）Determine the matricesandas defined in (4.45) for the acoustics equations (2.50) 解： （2.50） ，其中 系数矩阵 其特征值为 ， 对应的特征向量为 ， 定义右特征向量矩阵，同时有 １，当时，，则有，，故 ， ２，同样分析，得到时， ， ３，时 （ｂ）Determine the wavesandthat result from arbitrary dataandfor this system. 解： 故 Leveque书习题 4.2 （a）（b）（c） if we apply the upwind method（4.25）to the advection equationwith, and choose the time step so that, then the method reduces to . The initial data simply shifts one grid cell each time step and the exact solution is obtained, up to the accuracy of the initial data. ( If the datais the exact cell average of , then the numerical solution will be the exact cell average for every step. ) This is a nice property for a numerical method to have and is sometimes called the unit CFL condition. (a) Sketch figures analogous to Figure 4.5(a) for this case to illustrate the wave-propagation interpretation of this result. (b) Does the Lax-Friedrichs method (4.20) satisfy the unit CFL condition? Does the two-step Lax-Wendroff method of Section 4.7? (c) Show that the exact solution ( in the same sense as above ) is also obtained for the constant-coefficient acoustics equations (2.50) withif we choose the time step so thatand apply Godunov’s method. Determine the formulas forandthat result in this case, and show how they are related to the solution obtained from characteristic theory. 解： （a） （b） Lax-Friedrichs method （4.20） 将格式（4.20）应用到热传导方程，即为 当时上式变为 所以Lax-Friedrichs格式（4.20）满足unit CLF条件 two-step Lax-Wendroff method 当时 故 同样有 此时 故两步 Lax-Wendroff 方法满足unit CLF条件 （c） （2.50） 求系数矩阵的特征值和特征向量，根据右特征向量定义 Godunov方法一种是守恒型有限体积格式，如下 其中 这里需要在处近似或者精确求解Riemann问题（）得到。 对于线性常系数问题（2.50），Riemann问题和上述积分都可以精确计算： 其中的计算如下 另一方面，由于，则 取，可得到 下面证明上述即为特征理论得到的精确解。 如左图，为计算，需要知道沿特征线方向传播的特征分量，和沿特征线传播的t特征分量。 由于，这里计算如下 所以 同样，这里 故解是特征向量的线性组合 即证。 （d）Is it possible to obtain a similar exact result by a suitable choice of in the case where in acoustics ? 解：由4.1知系数矩阵右特征向量与取值无关，进而根据（c）计算过程可知时精确解仍然为 另一方面，考虑Godunov方法 1，当时，，，故 故无论如何