2016-11-22 13 views
0

FitzHugh–Nagumo 반응 확산 모델을 해결하기 위해 Alternating direction implicit method을 구현 중입니다. 나는 blog에 대한 파이썬 구현 예를 발견했다,하지만 난 방법에 오류가 있다고 생각 - 스텐실에 여기에 제시된 : 는 enter image description here 그것은 f를 반응에게 용어 을 곱 하프 타임 스텝 크기 없을 겁니다?파이썬에서 pde의 유한 차분 솔버를위한 대체 방향 암시 적 메서드

+1

http://math.stackexchange.com/에서 더 많은 답변을 얻으실 수 있습니다. 이는 수학 문제 일뿐입니다. – Daneel

+1

또는 http://scicomp.stackexchange.com에서 대규모 과학 컴퓨팅을위한 전문가가 있습니다. - 또한 분모'D/(2 * dp^2) '의 요소 2가 의심 스러우며 링크 된 위키 피 디아 페이지에서이를 확인한 것 같습니다. - 그리고 왜 'dt'가있는 것입니까? 다른 용어는 모두 완전 차액입니다. – LutzL

답변

1

차동 몫으로 차이 몫 교체, 하나는 원래 제기 작업이었다 방정식

U_t = D * (U_xx + U_yy) + f 

하지 않은 두 경우에

U_t = D/2 * U_xx + D/2 * U_yy + Δt*f 

을 가져옵니다. 이 f위한 U_pp, p=x,y1에서 U_t, D/(Δp^2)에 있다는

따라서 계수 1/(Δt/2)이어야한다.

수식은 차이 차수를 가진 수식과 Δt/2을 곱한 다음 단계가 혼합 된 것 같습니다. f의 요인 Δt/2해야한다는

그리고 그 다음 수식에 하나가 참으로 α_p=σ_p, p=x,y 같은 새로운 상수를 필요로하지 않습니다

는 당신은 맞다.