점탄성 고체 내 유체 확산 모델

초록

본 논문은 큰 변형을 겪는 점탄성 고체와 뉴턴 유체의 혼합계에 대한 확산 현상을 열역학적 엔트로피 생산률 최대화 원리를 이용해 모델링한다. 헬름홀츠 자유에너지와 소산률을 구체적으로 가정하고, 고체와 유체의 부분 응력, 상호작용력, 고체의 자연구성 변화 방정식을 도출한다. 혼합 부피 보존 조건을 제약으로 두어 실험 데이터와의 일치를 검증했으며, 고온 폴리이미드에 대한 용매 확산 및 외력에 의한 팽창 현상을 성공적으로 설명한다.

상세 분석

이 연구는 고체와 유체가 동시에 존재하는 복합 시스템을 고전적인 혼합물 이론에 기반하여 재구성하고, 열역학적 비가역 과정의 핵심인 엔트로피 생산률을 최대화하는 원칙을 적용함으로써 새로운 구성 관계를 제시한다. 저자는 먼저 두 성분의 질량 보존식과 운동량 보존식을 각각 정의하고, 전체 혼합물의 자유에너지 밀도를 헬름홀츠 자유에너지 형태로 가정한다. 여기서 중요한 점은 고체의 자연구성(Reference configuration)을 시간에 따라 진화시키는 변형 텐서를 도입함으로써, 점탄성 특성을 비선형적으로 기술한다는 것이다. 자유에너지 함수는 변형 이력과 용매 농도에 대한 의존성을 포함하도록 설계되었으며, 이를 통해 고체의 부분 응력 텐서와 유체의 점성 응력 텐서를 각각 도출한다.

소산률(Rate of Dissipation)은 두 부분으로 나뉘는데, 하나는 고체 내부의 점탄성 손실, 다른 하나는 유체 흐름에 의한 점성 손실이다. 저자는 이 두 소산항을 각각 제곱형식으로 가정하고, 라그랑주 승수를 이용해 엔트로피 생산률 최대화 조건을 적용한다. 결과적으로 얻어지는 구성 관계는 다음과 같다. 첫째, 고체의 부분 응력은 자유에너지의 변형 이력에 대한 편미분과 점탄성 소산항에 대한 비례항을 포함한다. 둘째, 유체의 부분 응력은 전형적인 뉴턴 유체식, 즉 점성 계수와 속도 구배의 곱으로 표현된다. 셋째, 두 성분 사이의 상호작용력은 상대 속도와 농도 구배에 의존하는 항으로 나타나며, 이는 다공성 매질에서의 Darcy 흐름과 유사한 형태를 가진다. 넷째, 고체의 자연구성 진화 방정식은 소산률을 최소화하는 방향으로, 즉 변형 이력과 용매 흡수에 의해 구동되는 비선형 미분 방정식으로 도출된다.

또한, 혼합물 전체 부피가 두 성분의 자연 상태 부피 합과 동일하다는 가정은 체적 보존 제약을 제공한다. 이를 라그랑주 승수 형태로 도입함으로써, 압력과 화학 포텐셜이 자연스럽게 연계된다. 저자는 이 제약을 통해 얻은 압력식이 고체와 유체의 부분 압력의 가중 평균임을 확인하고, 실험적 부피 변화와 일치함을 보였다.

모델 검증 단계에서는 고온 폴리이미드(Polyimide) 시료에 다양한 유기용매를 침투시킨 실험 데이터를 사용하였다. 실험에서 측정된 용매 흡수율, 팽창 변형, 그리고 응력-변형 곡선이 제시된 이론적 예측과 높은 상관관계를 보였으며, 특히 용매 종류에 따른 확산 계수와 점탄성 파라미터의 민감도가 모델에 의해 정량적으로 설명되었다. 마지막으로 외부 하중을 가했을 때 고체의 팽창 거동을 시뮬레이션한 결과, 응력 이완과 용매 재분배가 복합적으로 작용하여 비선형 팽창 특성을 나타냄을 확인하였다. 이러한 결과는 항공기 엔진 부품 등 고온 환경에서 용매 침투와 구조적 변형이 동시에 일어나는 상황을 예측하는 데 유용한 도구가 될 수 있다.