2.â«Î©
2.â«Î©
2.â«Î©
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
5 th week<br />
A. 一 維 問 題 之 理 解 : 概 念 和 步 驟<br />
(3 種 方 法 :<br />
○1 Trefftz method,<br />
○2 Fundamental solution method)<br />
加 權 殘 差 法 -weak formulation, 二 次 部 分 積 分 .<br />
⎛ d w ⎞<br />
⎝ dx ⎠<br />
dw 1<br />
dx 0<br />
du<br />
dw<br />
2<br />
1<br />
1<br />
○1<br />
∫ u − − + = 0<br />
0<br />
⎜<br />
2<br />
⎟dx<br />
∫ xwdx u w<br />
0<br />
2<br />
d w<br />
使 用 = 0 求 w. 代 入 積 分 式 , 使 用 Bc’s 求 解 , 求 出 邊 界 上 的 未 知 值 。<br />
2<br />
dx<br />
1<br />
0<br />
2<br />
d w<br />
○2 = −Δ( x − x i<br />
) < 負 號 為 配 合 積 分 式 符 號 ><br />
2<br />
dx<br />
求 出 w 表 示 式 ,w 稱 為 基 本 解 ,w 代 入 積 分 式 , 把 x 定 在 邊 界 上 , 求 解<br />
⎛ dw ⎞<br />
得 到 邊 界 上 的 未 知 函 數 值 ⎜u or ⎟ ( 此 時 , 邊 界 積 分 式 合 為 x 的 函 數 )<br />
⎝ dx ⎠ i<br />
du( xi<br />
)<br />
因 此 , 可 以 計 算 任 意 位 置 x<br />
i<br />
的 函 數 值 u<br />
i<br />
, 且 可 以 計 算 微 分 值 在 任 意 位 置<br />
d( xi<br />
)<br />
上 。<br />
B. 二 維 問 題 : 以 Elliptic 型 態 問 題 先 作 說 明<br />
使 用 G.E. 配 合 加 權 殘 差 法 , 二 次 部 分 積 分 得 到<br />
∫<br />
Ω<br />
2 ∂w<br />
∂u<br />
u ∇ wdΩ − ∫ u dΓ +<br />
Γ ∫ wdΓ = 0<br />
∂n<br />
Γ<br />
∂n<br />
* ∂w<br />
使 用 基 本 解 名 稱 w = u , = q<br />
∂n<br />
邊 界 積 分 式 寫 為<br />
令 ∇ 2 u<br />
*<br />
= −Δ<br />
i<br />
∫<br />
Ω<br />
∫<br />
Γ<br />
∗<br />
∂u , = q<br />
∂n<br />
2 *<br />
*<br />
*<br />
u ∇ u dΩ − uq dΓ + qu dΓ = 0<br />
代 入 積 分 式<br />
∫<br />
Γ<br />
−<br />
u i<br />
−<br />
∫<br />
Γ<br />
*<br />
uq dΓ +<br />
∫<br />
Γ<br />
*<br />
qu dΓ = 0<br />
∫<br />
*<br />
*<br />
or − u i − uq dΓ + qu dΓ = 0<br />
Γ<br />
∫<br />
後 續 型 式 改 為 ( )u=( )q<br />
2<br />
Note:1. 二 維 ∇ u = 0 的 基 本 解 ( 推 導 )<br />
Γ<br />
2 *<br />
2.∫ u∇<br />
u dΩ = −u<br />
Ω<br />
i<br />
( 積 分 出 來 的 型 式 必 須 相 同 , 會 有 singularity) 在 奇 異
點 積 分 上 需 要 討 論 , 對 ∇ 2 u = 0 而 言 , 其 基 本 解 可 表 示 為<br />
1<br />
(1) 3D Prob. u<br />
* 2 * i<br />
= ( ∇ u = −Δ<br />
)<br />
4πr<br />
以 圖 形 表 示<br />
r<br />
= x − x 1<br />
r 為<br />
1<br />
x 到 x 的 距 離<br />
r<br />
=<br />
x<br />
−<br />
x 1<br />
* 1 1<br />
(2) 2D Prob. u = ln<br />
2π<br />
r<br />
使 用 前 述 基 本 解 (3D,2D)<br />
∫<br />
Ω<br />
u∇<br />
2<br />
u<br />
*<br />
dΩ<br />
= −u<br />
i<br />
Ω 為 半 徑 ε 的 圓 (2D)<br />
球 (3D)<br />
最 後 令 ε → 0<br />
∫<br />
Ω<br />
u∇<br />
2<br />
u<br />
*<br />
dΩ<br />
∫<br />
=<br />
∫<br />
Ω ε<br />
其 中 上 式 ( )<br />
Ω − Ω ε<br />
( ) dΩ<br />
+ ( ) dΩ<br />
d Ω 項 為 0<br />
∫<br />
Ω − Ω ε<br />
2 *<br />
i<br />
積 分 結 果 ∫ u∇<br />
u dΩ<br />
= −u<br />
Ω<br />
在 邊 界 方 法 中 , Δ i 函 數 必 須 定 義 到 邊 界 上 , 引 進 邊 界 條 件 , 求 解<br />
邊 界 分 為 smooth 和 corner
i<br />
x 定 義 在 邊 界 上 時 ,∫<br />
Ω<br />
i<br />
2 * u<br />
u ∇ u dΩ = − (smooth 邊 界 )<br />
2<br />
(corner 邊 界 另 外 討 論 )<br />
邊 界 積 分 式 為<br />
u i<br />
∫<br />
*<br />
*<br />
− − q udΓ + u qdΓ = 0<br />
2<br />
Γ<br />
∫<br />
Γ<br />
i<br />
此 式 每 一 次 都 會 有 有 基 本 解 x 的 定 義 , 部 分 為 隱 含 表 示<br />
綜 合 整 理 :<br />
u i<br />
i<br />
*<br />
*<br />
x ∈ Γ , − − ∫q<br />
udΓ + ∫u<br />
qdΓ = 0<br />
2<br />
Γ<br />
i<br />
*<br />
*<br />
x ∈ Ω , − u i − ∫ q udΓ + ∫u<br />
qdΓ = 0<br />
Γ<br />
or − = ∫q<br />
udΓ − ∫<br />
u i *<br />
Γ<br />
Γ<br />
Γ<br />
* u qdΓ<br />
( 此 時 u 和 q 均 為 已 知 )<br />
Γ<br />
上 兩 式 均 含 一 般 的 任 意 形 態 邊 界 的 積 分 , 因 此 , 在 計 算 上 則 使 用 " 元 素 ",<br />
(element) 的 概 念 作 計 算 。 此 元 素 的 概 念 和 " 有 限 元 素 法 " 的 概 念 類 似 , 差 異 只 在<br />
於 基 本 解 的 正 確 使 用 。