초록
P1 비순응 요소를 이용하여 정상 비압축성 Navier-Stokes 유동의 위상최적화 문제를 푸는 방법을 제시한다. 본 연구는 Stokes 유동의 위상최적화 문제에 P1 비순응 요소를 적용하여 그 수치적 효용성을 보인바 있는 이전 연구에 대한 후속 연구이다. 비압축성 물질 해석에서 잠김현상이 발생하지 않으며 선형형상함수를 가지는 P1 비순응 요소의 장점이 관성항을 가지는 유체 문제의 해석과 설계에도 유효한 지를 파악하고자 한다. 일반적으로 사용되는 혼합정식화법과 비교하여 P1 비순응 요소의 사용은 벌칙 함수를 이용하여 연속 방정식을 따로 사용하지 않고 운동방정식에 부과할 수 있기 때문에 자유도의 개수를 감소시킬 수 있다. 벌칙 파라미터가 해의 정확도에 주는 영향과 적정 범위는 수치적으로 검토하도록 한다. 또한 보통의 사각 비순응 요소들이 요소면의 중앙에 절점을 가지고 고차의 형상함수를 지니는데 비하여, 본 연구에서 제시하는 P1 비순응 요소는 요소의 꼭지점에 절점을 가지고 {1, x, y}의 P1 형상함수로 구성됨으로써 수치적인 구현의 용이함이 일반 선형 사절점 요소와 동일하다. 제안한 방법의 효용성을 다양한 레이놀즈수에 따른 유동최적화 문제들을 살펴봄으로써 검증하도록 한다.
An alternative approach for topology optimization of steady incompressible Navier-Stokes flow problems is presented by using P1 nonconforming finite elements. This study is the extended research of the earlier application of P1 nonconforming elements to topology optimization of Stokes problems. The advantages of the P1 nonconforming elements for topology optimization of incompressible materials based on locking-free property and linear shape functions are investigated if they are also valid in fluid equations with the inertia term. Compared with a mixed finite element formulation, the number of degrees of freedom of P1 nonconforming elements is reduced by using the discrete divergence-free property; the continuity equation of incompressible flow can be imposed by using the penalty method into the momentum equation. The effect of penalty parameters on the solution accuracy and proper bounds will be investigated. While nodes of most quadrilateral nonconforming elements are located at the midpoints of element edges and higher order shape functions are used, the present P1 nonconforming elements have P1, {1, x, y}, shape functions and vertex-wisely defined degrees of freedom. So its implentation is as simple as in the standard bilinear conforming elements. The effectiveness of the proposed formulation is verified by showing examples with various Reynolds numbers.