問題の収束履歴を監視しています。記録出力は以下の通りです。
Time = 24
Calculate volume forces from actuator disk
Total thrust = -8.46832
Total torque = 1.03471
ADisk volume = 0.0632799
smoothSolver: Solving for Ux, Initial residual = 0.000447755, Final residual = 2.68745e-05, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.0107909, Final residual = 0.000812227, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.0103399, Final residual = 0.000786661, No Iterations 2
GAMG: Solving for p, Initial residual = 0.123954, Final residual = 0.00958268, No Iterations 6
time step continuity errors : sum local = 7.42808e-05, global = -4.25546e-05, cumulative = 0.000413527
smoothSolver: Solving for epsilon, Initial residual = 0.00197379, Final residual = 0.000172248, No Iterations 1
smoothSolver: Solving for k, Initial residual = 0.000510499, Final residual = 2.78594e-05, No Iterations 2
ExecutionTime = 124.63 s ClockTime = 125 s
Time = 25
Calculate volume forces from actuator disk
Total thrust = -8.49093
Total torque = 1.03723
ADisk volume = 0.0632799
smoothSolver: Solving for Ux, Initial residual = 0.000409002, Final residual = 2.59552e-05, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.0103191, Final residual = 0.00077024, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.00985658, Final residual = 0.000742227, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0390756, Final residual = 0.00247253, No Iterations 7
time step continuity errors : sum local = 5.39785e-05, global = 3.40394e-05, cumulative = 0.000447566
smoothSolver: Solving for epsilon, Initial residual = 0.00182397, Final residual = 0.000157739, No Iterations 1
smoothSolver: Solving for k, Initial residual = 0.000465916, Final residual = 2.75864e-05, No Iterations 2
ExecutionTime = 129.45 s ClockTime = 130 s
Time = 26
Calculate volume forces from actuator disk
Total thrust = -8.51463
Total torque = 1.03953
ADisk volume = 0.0632799
私がしたいのは、特定の時間に値をコピーすることです。たとえば、この場合はTime=26
値をコピーします。
Total thrust = -8.51463
Total torque = 1.03953
次の形式を使用してください。
1.03953 -8.51463.
シェルスクリプトを使用してこれを行うのに役立つ人はいますか?
答え1
シェルスクリプトを要求しましたが、次のことがawk
できるようにしたいと思います。
find_thrust_torque.awk:
/^Total thrust =/ {thrust = $4}
/^Total torque =/ {torque = $4}
/^Time =/ {if (found) exit; if ($3 == time) found=1}
END {print torque " " thrust}
テストを受けてください:
$ awk -v time=25 -f find_thrust_torque.awk file1
1.03723 -8.49093
$ awk -v time=26 -f find_thrust_torque.awk file1
1.03953 -8.51463
答え2
ここでは、スロット範囲と次のスロット(または最後のスロットの場合はeof)を選択します。これは、現在のスロットにトルク/推力値がない場合に、次のスロットへのデータの取得を防ぐために行われます。したがって、古いデータが得られ、エラー報告は発生しません。このH
コマンドは、パターンスペースを予約済みスペースに追加します。g
予約済み領域が検索され、パターン空間に配置されます。
tslot=26; # input the time you want the thrust/torque data for
sed -ne '
/^Time = '"${tslot}"'$/,/^Time =/!d
/^Total thrust =/{H;d;}
/^Total torque =/{H;d;}
g;/\n.*\n/!d
/thrust.*torque/s/\(\n.*\)\(\n.*\)/\2\1/
s/\n[^=]*=//gp;q
' logfile
結果
1.03953 -8.51463