DOI QR코드

DOI QR Code

Quantum annealing을 통한 hybrid composite의 두께 방향 열전도 특성 개선

Improving Through-thickness Thermal Conductivity Characteristic of Hybrid Composite with Quantum Annealing

  • Sung wook Cho (Department of Mechanical Engineering, Graduated School, Kongju National University) ;
  • Seong S. Cheon (Department of Mechanical Engineering, Graduated School, Kongju National University)
  • 투고 : 2024.06.11
  • 심사 : 2024.06.20
  • 발행 : 2024.06.30

초록

본 연구에서는 탄소섬유 강화 플라스틱(CFRP)에 얇은 구리 필름(Cu film)이 배치된 hybrid composite이 제안되었으며, 두께방향 열전도도가 최대가 될 수 있는 Cu film 배치조합을 도출하는데 양자 어닐링(Quantum Annealing)이 적용되었다. CFRP의 각 ply와 Cu film간의 상관관계 분석이 유한요소 해석을 통해 수행되었으며, 수행된 결과를 바탕으로 조합 최적화 문제가 정의되었다. 정의된 문제를 양자 어닐링에 임베딩하기 위한 공식화 과정이 진행되었으며 이를 통해 CFRP의 각 ply에 투입될 수 있는 Cu film 수량에 관한 목적함수와 제약조건이 수식으로 구현되었다. 공식화된 수식은 D-Wave 양자 어닐러에 임베딩되기 위해 Ocean SDK(software development kit)와 Leap을 통해 프로그래밍 되었으며, 양자 어닐링 과정을 통해 두께 방향 열전도도가 최대를 만족하는 최적의 Cu film 배치 조합이 도출되었다. 도출된 배치 조합은 투입될 수 있는 Cu film의 수량이 적어질수록 단순한 배치 형태를 나타내었으며, 수량이 많아질수록 세밀한 배치를 보였다. Cu film의 배치 수량에 따라 생성된 최적 조합들은 두께 방향으로의 고유 열전도 경로를 나타내었으며, Cu film의 횡방향 배치 자유도가 두께 방향 열전도도 결과에 민감하게 나타날 수 있음을 보였다.

This study proposes a hybrid composite where a thin copper film (Cu film) is embedded in carbon fiber reinforced plastic (CFRP), and quantum annealing is applied to derive the combination of Cu film placement that maximizes the through-thickness thermal conductivity. The correlation between each ply of CFRP and the Cu film is analyzed through finite element analysis, and based on the results, a combination optimization problem is formulated. A formalization process is conducted to embed the defined problem into quantum annealing, resulting in the formulation of objective functions and constraints regarding the quantity of Cu films that can be inserted into each ply of CFRP. The formulated equations are programmed using Ocean SDK (Software Development Kit) and Leap to be embedded into D-Wave quantum annealer. Through the quantum annealing process, the optimal arrangement of Cu films that satisfies the maximum through-thickness thermal conductivity is determined. The resulting arrangements exhibit simpler patterns as the quantity of insertable Cu films decreases, while more intricate arrangements are observed as the quantity increases. The optimal combinations generated according to the quantity of Cu film placement illustrate the inherent thermal conductivity pathways in the thickness direction, indicating that the transverse placement freedom of the Cu film can significantly affect the results of through-thickness thermal conductivity.

키워드

1. 서론

탄소섬유 강화 플라스틱(Carbon-Fiber Reinforcement Plastics, CFRP)은 높은 비강도와 비강성 그리고 우수한 내식성을 나타내는 복합재이다. 또한 CFRP는 내열성이 높아 고온에서도 안정적인 기계적 특성 유지가 가능하다[1]. 이러한 장점들로 인해 항공 우주 및 우주 응용 분야에서는 지속적으로 광범위하게 CFRP가 적용되고 있다[2-4]. 특히 기술이 발전함에 따라 CFRP로 구성된 고성능 장치들은 효과적인 방열 성능이 요구되어진다. 열전도도에 있어서 CFRP는 섬유가 지배적으로 더 나은 면내 방향(In-plane)의 열전도도를 나타내지만[14] 기지의 절연특성으로 인해[5-7] 두께 방향의 열전도도는 매우 낮은 값을 나타내게 된다. 따라서 기지와 섬유의 서로 다른 열적 특성으로 인해 CFRP는 고유의 열적 이방성을 나타내게 된다. 이러한 고르지 못한 열전도 특성은 효과적이지 못한 방열 성능으로 이어지게 될 수 있다[8]. 이러한 문제는 나노 필러 형태의 탄소 나노튜브[9,10], 그래핀[11], 금속 나노 입자[12] 등을 첨가한 hybrid composite 형태로 개선되어질 수 있지만 대다수의 접근 방식들은 실제 실험이 주를 이루기에 시간과 비용이 많이 소요되는 단점이 존재한다. 이러한 현실적 제약을 극복하고자 Kim 등[13]은 meta-heuristic algorithm 중 하나인 유전 알고리즘(Genetic Algorithm)을 이용하여 복합재의 열전도도 향상을 위한 최적 구조 설계를 수행한 바 있다. 최근에는 AI와 머신러닝의 학습을 통한 접근 방식[14]으로 앞선 단점들의 극복과 실패에 대한 리스크의 최소화를 구현할 수 있지만 실제 해결하고자 하는 문제의 대다수는 그 복잡성에 의한 여러 로컬 최소값과 안장점의 존재로 인해 기존 알고리즘을 통한 효율적인 최적값 수렴이 어렵다[15]. 특히 변수가 다양 해질수록 경우의 수는 기하급수적으로 증가하게 되는데, 해당 대규모 전역 영역에서의 최적값 도출에 있어 양자 컴퓨터는 앞서 제시된 제약을 극복하기 좋은 전략이 될 수 있다[15-17]. 현재 양자 컴퓨터는 양자 논리 게이트(Quantum logical gate)와 양자 어닐링(Quantum Annealing, QA)으로 분류되며 조합 최적화에 특화되어 있는 양자 어닐러가 가장 빠르게 상용화되어 있는 상태이다[18]. 이러한 양자 어닐러를 통한 QA를 수행하기 위해서는 먼저 조합 최적화에 해당하는 문제가 정의되어야 하며, 양자 어닐러에 임베딩하기 위해서는 정의된 문제의 목적함수를 Quadratic Unconstrained Binary Optimization (QUBO) 형태나 Ising model 형태 중 하나로 공식화되어야 한다[19,20]. 따라서 공식화 과정이 QA의 핵심이라 할 수 있다. Imanaka 등[21]은 방열판 복합재의 혼합비에 따른 높은 방열성, 부착성, 경량화를 만족하는 조합 최적화를 수행하였다. 이를 위해 열전도도, 열팽창계수 및 비중을 QUBO형태로 공식화하였다. 공식화된 목적 함수는 QA와 호환이 가능하지만 해당 연구에서는 Fujitsu Digital Annealer(DA)를 통해 최적화가 수행되었다. 조합 후보로 선정된 10개 재료에 대한 총 변수 1024개의 조합 최적화는 10초 이내에 빠른 최적값 도출을 나타내었으며 도출된 조합은 혼합물 법칙을 통해 해당 복합재료의 특성 검증이 이루어졌다.

본 연구에서는 Fig. 1에서와 같이 제한된 개수의 얇은 구리 필름(Cu film)이 배치된 CFRP 기반 hybrid composite(Cu film CFRP)이 제안되었으며 이를 통해 최적의 두께 방향 열전도도를 나타낼 수 있는 Cu film의 배치를 QA를 통해 효율적으로 도출하고자 한다. 먼저 열전도도의 도출을 위해 Laser Flash Analysis(LFA) 방식이 유한요소해석(Finite Element Analysis, FEA)으로 구축 및 적용되었으며 QA를 위한 문제 정의를 위해 Cu film 배치 조건과 투입될 각 ply의 적층 순서간 상관관계 분석 이루어졌다. 해당 상관관계 결과를 바탕으로 조합 최적화 문제 정의가 수행되었으며 Cu film CFRP의 과도한 무게 증가를 제한하기 위한 Cu film 배치 제약이 적용된 목적함수가 공식화되었다. 이렇게 공식화된 목적함수는 양자 컴퓨터에 임베딩됨으로써 QA를 통한 조합 최적화가 이루어졌으며 도출된 Cu film의 최적배치 조합은 FEA를 통해 pure CFRP 대비 향상된 두께 방향 열전도도의 비교 분석으로 진행되었다.

BHJRB9_2024_v37n3_170_2_f0001.png 이미지

Fig. 1. Schematic of Cu film placement in hybrid composite

2. FEA 구축과 Ply별 Cu film 배치에 따른 상관 관계 분석

2.1 두께 방향 열전도도 도출을 위한 FEA

Cu film CFRP의 두께 방향 열전도도 도출을 위한 방식으로 Transient method 중 하나인 LFA 방식이 적용되었다. 측정 과정은 ASTM-E1461[22]에 따라 Fig. 2에서와 같이 시편에 laser pulse를 가하여 반대편의 온도 변화(∆T) 검출이 이뤄지게 된다. 측정된 ∆T는 Fig. 3과 같은 시간에 따른 온도 변화 그래프로 나타나며 시편 두께에 따른 열확산도 계산을 거쳐 열전도도가 도출되는 방식이다. 이상적인 단열의 경우 열확산도는 식 (1)에 의해 계산된다[23].

\(\begin{align}\alpha=0.1388 \frac{d^{2}}{t_{1 / 2}}\end{align}\)       (1)

BHJRB9_2024_v37n3_170_3_f0001.png 이미지

Fig. 2. Laser flash analysis method

Fig. 3. ∆T -time graph derived from the bottom side

여기서 α는 열확산도, 0.1388은 heating zone과 detector의 최대 온도차의 50%(1/2ΔTmax)에 해당하는 상수, d는 시편의 두께, t1/2은 1/2ΔTmax까지의 도달 시간을 나타낸다.

도출된 열확산도는 최종 식 (2)를 통해 열전도도로 계산되어진다.

λ = α · ρ · Cp       (2)

여기서 λ는 열전도도, ρ는 시편의 밀도, Cp는 시편의 비열을 나타낸다.

이러한 LFA 방식은 시편 가열 시간이 짧고 시편의 크기가 다른 열전도도 측정법에 비해 상대적으로 작다는 장점이 존재한다[24]. 이는 해석 시간 단축에 용이하게 작용될 수 있다.

LFA를 구현한 FEA 모델은 LS-DYNA를 통해 구축되었으며 Fig. 4과 같다. 시편의 형상은 가로, 세로 10 mm의 정사각형이며 요소의 총 개수는 7,500개이다. CFRP의 각 ply 두께는 0.1 mm이고 Uni-direction(UD) type에 [0]3T의 조건으로 구현되었으며 Bottom side에서 Top side 방향으로 1,2,3순의 ply 번호가 할당되게 된다. CFRP의 열적 이방성은 In-plane에서의 열전도도인 섬유방향 k11과 섬유의 횡방향 k22 그리고 두께방향 k33으로 정의될 수 있으며, MAT-021_ORTHOTROPIC_THERMAL을 통해 구현되었다. Cu film의 형상은 폭 0.2 mm, 두께 33 μm이며 길이는 CFRP와 동일한 10 mm이다. Cu film은 구리로 구성되어 있으므로 T01-THERMAL_ISOTROPIC을 통해 열적 등방성이 구현되었다. Cu film의 경우 배치되는 개수가 사이 간격을 나타내는 Pitch를 통해 결정되며 배치는 Fig. 4에서와 같이 ply 두께의 1/2 위치에 적용되게 된다. 또한 Cu film은 CFRP의 섬유 방향과 일치하도록 배치되게 된다. 배치되는 경우는 Table 1에서와 같이 총 4가지 type에 한해서 배치되어질 수 있다. Laser pulse에 의한 heating zone은 Top side 면적의 70%에 10 ms동안 20 mW/mm2의 heat flux가 투입되며 반대편인 Bottom side에서 Fig. 3에서와 같은 시간에 따른 온도변화 그래프가 도출되게 된다. 그 외 측면은 모두 단열로 설정되었으며 Top side와 Bottom side의 경우 대류와 복사에 의한 열 전달은 고려되지 않는다.

BHJRB9_2024_v37n3_170_3_f0003.png 이미지

Fig. 4. Example of hybrid composite FE model

Table 1. Cu film placement spacing and corresponding number

BHJRB9_2024_v37n3_170_3_t0001.png 이미지

2.2 Thermal properties

CFRP는 섬유 부피분율에 따라 그 열전도도가 다르게 나타난다. 따라서 FE 모델에 적용될 CFRP의 열적 특성은 혼합물 법칙을 통해 섬유와 기지 각각의 고유 열적 특성이 적용되어 도출되었다. 섬유와 기지는 각각 AS4와 8552 resin의 열적 특성이 적용되었으며 그 특성은 Table 2와 같다[25]. CFRP의 전체 밀도(ρc)는 식 (3)을 통해 계산되게 된다.

ρc = ρfυf+ ρmυm       (3)

Table 2. Properties of fiber & matrix

여기서 ρf와 ρm은 각각 섬유와 기지의 밀도를 의미하고 υf와 υm은 섬유와 기지의 부피분율을 나타내며 υf= 0.514이고 υm = (1 − υf)에 의해 0.4286으로 적용되었다. Fig. 4에서 정의된 CFRP의 열전도도 k11, k22, k33은 D.J Radcliffe 등[26]의 이론적 수식인 식 (4)와 식 (5)를 통해 도출될 수 있으며 k11의 경우 식 (4)를 통해 계산되어진다.

k11 = kfυf + kmυm       (4)

여기서 kf와 km은 각각 섬유와 기지의 열전도도를 의미한다. k22와 k33의 경우 UD 형태의 CFRP에서 각 열전도도는 서로 동일하다고 가정할 수 있으며[27] 이를 식 (5)와 같이 나타낼 수 있다.

\(\begin{align}k_{22}=k_{33}=\frac{k_{11}}{\varepsilon}\end{align}\)       (5)

여기서 ε은 경험적 비례 계수이며 본 연구에서는 AS4/8552의 두께 방향 열전도도인 k33의 실제 측정 결과[25]에 근거하여 그 평균 오차율 0.132%를 만족하는 ε인 8.109가 적용되었다.

비열(cp)의 경우 식 (6)을 통해 나타낼 수 있게 된다[28].

cp = cfυf + cmυm       (6)

여기서 cf와 cm은 각각 섬유와 기지의 비열을 나타낸다. 따라서 위 혼합물 법칙 수식들을 통해 도출된 pure CFRP의 열적 특성은 Table 3과 같다.

Table 3. Properties of pure CFRP

BHJRB9_2024_v37n3_170_4_t0002.png 이미지

Cu film의 경우 열적 등방성을 만족하므로 Table 4와 같은 열적 특성이 적용되었다[29].

Table 4. Properties of Cu film

BHJRB9_2024_v37n3_170_4_t0003.png 이미지

2.3 Ply별 Cu film 배치에 따른 상관관계 분석

QA 적용을 위한 문제 정의와 공식화 수행을 위해 Cu film이 배치된 ply의 적층 순서에 따른 k33간 상관관계 분석이 선행되었다. 적용된 Cu film은 Table 1의 총 4가지 type이 각각 적용되었고 Cu film이 배치된 각각의 Cu film CFRP는 pure CFRP(k33=0.572 W/m℃)와의 비교를 통해 k33의 증가율로 나타내었다. Fig. 5(a)는 1에서 3번 ply에 Type-A가 각각 배치되었을 때의 결과를 나타내며 2번 ply에서 18.48%(0.659 W/m℃), 3번 ply에서 16.17% (0.665 W/m℃), 1번 ply에서 15.24% (0.659 W/m℃)로 도출되었다. Type-A는 Cu film의 투입량이 가장 많으며 배치 간격이 매우 조밀하다. 이는 열 확산에 있어 매우 유리하며 곧 열전도도의 증가로 나타나게 된다. 따라서 heating zone인 3번 ply와 1번 ply 사이에서 효율적으로 열전달이 가능한 2번 ply에서의 k33가 최대로 나타나게 된 것으로 사료되어진다. Fig. 5(b)는 Type-B에 해당하며 2번 ply에서 12.72% (0.645 W/m℃), 1번 ply에서 10.90% (0.635 W/m℃), 3번 ply에서 10.63% (0.633 W/m℃) 순으로 증가율이 도출되었다. Type-B 역시 Type-A와 같이 2번 ply에서의 k33 증가율이 가장 크게 나타났으며 오히려 1번 ply에서의 k33가 그 다음으로 크게 도출되었다. Fig. 5(c)의 경우 Type-C에 해당하며 3번 ply에서 8.79% (0.622 W/m℃), 2번 ply에서 8.64% (0.621 W/m\degc), 1번 ply에서 8.30% (0.619 W/m℃)의 증가 경향을 보였다. Fig. 5(d)의 Type-D 역시 3번 ply에서 6.37% (0.609 W/m℃), 2번 ply에서 6.19% (0.608 W/m℃), 1번 ply에서 6.14% (0.607 W/m℃)로 Type-C와 유사한 증가 추세를 나타내었다. Type-C와 Type-D의 경우 Cu film의 배치 간격이 Type-A와 B에 비해 크며 투입량 역시 적다. 이로 인해 ply의 층간 열확산이 느리게 작용한 결과로 보이며 Type-A와 B에 비해 열에 민감한 반응을 나타내지 못한 것으로 이해할 수 있다. 따라서 heating zone에 가까워질수록 k33가 점점 증가하는 추세를 나타낸 것으로 사료되어진다.

BHJRB9_2024_v37n3_170_4_f0001.png 이미지

Fig. 5. k33 according to Cu film type and stacking sequence: (a) Type-A, (b) Type-B, (c) Type-C, (d) Type-D

3. 문제 정의와 공식화를 통한 QA적용

3.1 조합 최적화 문제 정의

본 조합 최적화 문제는 Cu film CFRP의 ply수 확장([0]3T→[0]5T)을 통해 해결하고자 하는 조합 최적화 문제의 복잡성을 증가시켜 정의되었다. 총 5개의 ply에 개별적으로 적용되어질 수 있는 조합 후보들은 Cu film의 Type-A부터 Type-D와 아무것도 투입이 안된 상태인 None을 포함하여 총 5가지이고 Fig. 6와 같이 총 25가지의 조합 후보들로 구성되었으며 전체 가용 조합의 수는 55=3,125이다. ply 수가 확장됨에 따라 각 Type별 Cu film의 각기 다른 k33 증가율 표현은 Fig. 5의 결과를 바탕으로 각 ply 사이 k33를 선형보간 하여 유도되었으며, 조합 후보간 각 ply에서의 k33 우세함을 나타낼 수 있는 가중치(ki)로 적용되게 된다. 각 ply에서의 k33와 Cu film 배치 수량 간 데이터 일관성을 위한 정규화(Normalization)가 수행되었으며 전체 조합 후보들은 Table 5와 같이 나타낼 수 있다. 이러한 조합 후보들을 갖는 최적화 문제에는 전체 배치 가능한 Cu film의 무게제한을 두는 제약 조건이 부여되게 되며 이를 위해 Type별 Cu film의 배치 수량이 제약조건의 가중치(ωi)로 적용되게 된다. 또한 각 ply에서는 해당 ply 번호에 해당하는 Cu film이 중복으로 선정될 수 없도록 하는 제약조건이 추가로 적용되게 된다.

BHJRB9_2024_v37n3_170_5_f0001.png 이미지

Fig. 6. Conditions that can be placed on each ply

Table 5. Combination candidates that can be applied to each ply

BHJRB9_2024_v37n3_170_5_t0001.png 이미지

3.2 정의된 문제 공식화

QA는 목적함수로 시스템의 전체 에너지 함수(E(x))를 양자의 물리적 현상(superposition, entanglement, tunneling)을 이용해 최소값을 도출하는 과정이다. 이러한 에너지 함수는 Ising model과 QUBO로 표현될 수 있으며, Ising model의 경우 -1, 1인 spin 형태의 변수가 적용되고 QUBO의 경우 binary 형태인 0, 1의 변수가 적용된다. 위 문제 정의에서 각 ply에 할당될 Cu film 배치 조건은 배치가 안되었을 경우 0, 배치되었을 경우 1인 binary 형태로 표현될 수 있다. 따라서 binary 변수를 사용하는 QUBO 형태로 목적함수를 공식화할 수 있으며 식 (7)과 같다[30].

\(\begin{align}E(x)=\sum_{i}^{n} a_{i} x_{i}+\sum_{i<j}^{n} b_{i j} x_{i} x_{j}\end{align}\)       (7)

여기서 xi와 xj는 xi, xj∈{0,1}를 만족하는 binary 변수를 나타낸다. 𝑎i는 linear coefficient를 의미하며 앞서 선정한 가중치가 해당 항에 적용되게 된다. bij는 quadratic coefficient를 의미하는데, 본 연구에서는 Cu film의 배치 유무를 나타내는 binary 변수항만 적용되기에 xi만 사용된다. 이로 인해 QUBO형식이 선형계획법(linear programming)형식으로 정의될 수 있게 된다. 따라서 xi 항만이 사용되기에 bij항은 적용되지 않게 된다.

앞서 정의된 문제에 대한 목적함(Ek(k))는 아래 식 (8)과 같이 나타낼 수 있다.

\(\begin{align}E_{k}(x)=-\sum_{i=1}^{n} k_{i} x_{i}\end{align}\)       (8)

여기서 i는 i∈{1, 2, ..., 25}를 만족하며 Table 5의 Index No.를 의미한다. ki는 i에 해당하는 k33가 k로 적용되게 된다. 각 ply에 적용될 Cu film의 배치 조합은 최대가 되어야 하지만 QA에서는 에너지 함수가 최저인 바닥상태가 최적 솔루션으로 도출되므로 전체 수식에 -부호를 부여하여 최대값 도출이 이뤄지도록 적용되었다.

Cu film 무게 제한을 나타내는 제약조건은 식 (9)와 같이 나타낼 수 있다.

\(\begin{align}\sum_{i=1}^{n} \omega_{i} x_{i} \leq W\end{align}\)       (9)

여기서 ωi는 Table 5의 정규화 된 Cu film 투입 수량을 의미하며 W는 Cu film CFRP에 투입될 수 있는 전체 Cu film의 수량을 나타낸다.

각 ply에서 선정될 수 있는 Cu film 후보들은 총 5개를 만족해야 하며 각 ply에서의 binary 변수의 합이 1이 되는 것으로 표현될 수 있다. 따라서 해당 제약조건은 식 (10)과 같이 나타낼 수 있다.

\(\begin{align}\sum_{i=5 p-4}^{5 p} x_{i}=1\end{align}\)       (10)

여기서 p는 Table 5에서와 같이 각 ply의 번호를 의미하며 p∈{1, 2, 3, 4, 5}를 만족한다.

3.3 QA를 통한 조합 최적화 수행

정의된 조합 최적화 문제로부터 공식화된 수식은 D-Wave에서 제공하는 python 기반 SDK(software development kit)인 Ocean의 dimod API(Application Programming Interface)를 통한 프로그래밍과 D-Wave Leap 양자 클라우드 서비스를 통한 양자 컴퓨터로의 임베딩 과정이 Fig. 7의 순서로 수행되었다. 정의된 문제는 제약조건이 포함되어 있으므로 초기 Initialize에서 CQM(constrained quadratic model) 선언을 통해 D-wave에서 제공하는 CQM sampler가 적용되었다. 조합 후보의 정의는 dictionary 형태로 입력되게 되며 각각의 dictionary 요소들은 binary항으로 표현되게 된다. 따라서 식 (8)의 목적함수는 binary항으로 표현된 모든 조합후보와 ki간 곱의 합으로 프로그래밍 되며, 식 (9)의 경우 목적함수와 동일하게 표현되나 배치가능한 Cu film의 최대 설정 개수 이하를 만족하도록 제약조건이 적용되게 된다. 본 연구에서는 6가지 경우(Cu film의 최대 배치 제한 개수: 24, 32, 40, 48, 56, 64)에 대해서만 수행되었다. 식 (10)의 경우는 정해진 index No.안에서 모든 binary항의 합이 1이 되도록 작성되었다. LeapHybridCQMSampler의 선언을 통해 sampler를 정의하고 작성된 목적함수와 제약조건들은 양자컴퓨터에 임베딩되어 최적화가 수행되게 되며 각 경우마다 5번씩 반복 수행되었다.

BHJRB9_2024_v37n3_170_6_f0001.png 이미지

Fig. 7. Flow chart for embedding into quantum computer

4. 최적화 결과

Cu film의 최대 배치 제한 개수로 선정된 6가지 경우에 대한 QA 결과를 통해 도출된 최적의 Cu film 배치 순서는 Table 6과 같이 도출되었다. 각 Quantity of Cu films는 제약조건으로 지정된 가용 배치 수량을 의미하며 도출된 최적 배치 조합은 ply 번호 순으로 나열되었다. 각 최적 조합들은 LFA 방식의 FEA를 통해 k33가 계산되었으며, 밀도(ρ)는 각 최적 조합이 적용된 Cu film CFRP의 전체 밀도를 나타낸다. 도출된 각각의 최적 조합을 보면 배치될 수 있는 Cu film 수량이 많아질수록 배치 간격이 조밀한 Type 위주로 조합이 형성된 것을 볼 수 있으며, 반대로 수량이 적어질수록 조합이 단조롭게 구성되는 양상을 볼 수 있다. 이는 Cu film 수량 제한이 제약조건으로써 우세하게 작용한 결과로 사료되어진다.

Table 6. Results of QA optimization based on constraints

BHJRB9_2024_v37n3_170_7_t0001.png 이미지

Fig. 8은 Table 6의 결과를 바탕으로 Cu film 배치 수량에 따른 pure CFRP 대비 k33와 밀도의 증가율을 나타낸다. Cu film 수량이 증가함에 따라 Cu film CFRP의 밀도 역시 선형적으로 증가하는 것을 볼 수 있다. 하지만 k33증가율의 경우 40개 다음으로 48개(0.800 W/m℃)에서 감소되었다가 56개(0.817 W/m℃)에서 다시 증가하는 결과를 보였다. 해당 수량(40, 48, 56개)에서의 두께 방향 열전달 상태를 나타내기 위해 Fig. 9에서 가리키는 단면 영역에서의 온도 분포 결과가 Fig. 10과 같이 도출되었다. 해당 결과는 Cu film CFRP의 Top과 Bottom side에서의 최대 온도차(∆Tmax)에 도달하였을 때의 온도 분포 결과이다. 각 온도 분포 결과는 특정 조합으로 형성된 Cu film 배치 패턴의 영향을 나타내며 곧 두께 방향으로의 열전도 경로(Heat conductive pathway)가 형성되었음을 보여준다. 여기서 Fig. 10(b)에서의 Cu film 배치는 Fig. 10(a)와 Fig. 10(c)에 비해 상대적으로 Y축 방향으로의 배치 자유도가 낮음을 알 수 있다. Cu film은 CFRP 대비 빠른 열전도가 가능하기에 서로 인접하게 배치될 수록 연쇄적인 열전도 과정이 유도될 수 있지만 Fig. 10(c)와 같은 경우 열전도 경로가 국부적으로 형성되어 Fig. 10(a)와 Fig. 10(b)에 비해 상대적으로 감소한 k33 증가율이 도출된 것으로 사료되어진다.

BHJRB9_2024_v37n3_170_7_f0001.png 이미지

Fig. 8. Results of k33 and ρ increase rate

BHJRB9_2024_v37n3_170_7_f0002.png 이미지

Fig. 9. Heat transfer area to be observed

BHJRB9_2024_v37n3_170_7_f0003.png 이미지

Fig. 10. Heat conductive pathway in thickness direction according to combination: (a) 40 Cu films, (b) 48 Cu films, (c) 56 Cu films

5. 결론

본 연구에서는 CFRP(UD, [0]3T)의 k33 향상을 위해 제한된 개수의 얇은 Cu film이 배치된 Cu film CFRP가 제안되었으며, 한정된 배치 수량에 따른 최적의 Cu film 배치 조합 도출이 QA를 통해 수행되었다.

조합 최적화 문제 정의를 위해 CFRP의 ply 순서와 각 ply에 배치될 Cu film 개수간 상관관계 분석이 진행되었고 이러한 분석은 FEA를 통해 k33를 도출함으로써 수행되었다. 이를 통해 [0]5T로 확장된 CFRP에서의 k33가 최대가 되는 Cu film 배치에 대한 조합 최적화 문제가 정의되었다. QA 적용을 위해 정의된 문제의 목적함수와 제약조건에 대한 공식화가 수행되었으며, 양자 컴퓨터에 공식화된 문제를 임베딩함으로써 QA를 통한 6가지 경우의 Cu film 배치에 대한 최적 조합이 도출되었다. 주요 결론은 다음과 같다.

Cu film의 배치 조건과 적층 순서에 따른 k33의 증가량은 Cu film의 배치 간격이 좁을 수록 중간 ply에서 높게 나오며, 배치 간격이 넓을 수록 상부 ply에서 높게 도출되게 된다. 중간 ply의 Cu film 배치 간격이 좁은 경우 상부 ply에서 하부 ply로의 층간 열확산에 있어 효율적으로 작용될 수 있기에 k33가 높게 도출된 것으로 사료된다. 반면 Cu film 배치 간격이 넓은 경우에는 간격 사이에 열전도도가 낮은 CFRP로 인해 중간과 하부 ply에서의 열확산이 더디게 발생하게 된다. 따라서 heating zone의 열을 바로 전도할 수 있는 상부 ply에서의 k33가 높게 도출된 것으로 사료된다.

각 Cu film 최적 배치 조합은 Cu film 배치 수량이 많아질수록 배치간격이 조밀한 Type 위주로 형성되었으며, 수량이 적을수록 단조로운 구성의 조합이 형성되었다. 이는 Cu film의 수량 제한이 제약 조건으로 강하게 작용한 결과에 의한 것으로 사료된다.

Cu film의 조합을 통해 형성된 각 배치 패턴은 두께 방향으로의 고유 열전도 경로를 나타내었다. 이러한 열전도 경로의 형성은 Cu film의 횡방향 배치 자유도에 따라 다르게 나타나며 배치 자유도가 높을 수록 더욱 개선된 k33가 도출될 수 있을 것으로 판단된다.

향후 연구에서는 섬유의 배향 조건까지 고려한다면 더욱 다양한 조합에 의한 개선된 k33 결과를 도출할 수 있을 것으로 사료되며 이러한 연구 결과는 향후 hybrid composite의 열 관리 성능을 개선하는데 기여할 수 있을 것으로 기대된다.

후기

본 논문은 한국 정부(MSIT)의 재원으로 국가과학기술연구회(NST)의 지원을 받아 수행되었습니다(CRC23011-210). 또한 이 논문은 2023년도 정부(과학기술정보통신부)의 재원으로 한국연구재단 양자정보과학 연구개발생태계 조성 사업의 지원을 받아 수행된 연구임(2020M3H3A1110365).

참고문헌

  1. Ren, Y., Wang, H., Guan, Z., and Yang, K., "Evaluation of the Properties and Applications of FRP Bars and Anchors: A Review," Reviews on Advanced Materials Science, Vol. 62, No. 1, 2023.
  2. Pozegic, T.R., Anguita, J.V., Hamerton, I., Jayawardena, K.D.G.I., Chen, J.S., Stolojan, V., Ballocchi, P., Walsh, R., and Silva, S.R.P., "Multi-Functional Carbon Fibre Composites using Carbon Nanotubes as an Alternative to Polymer Sizing," Scientific Reports, Vol. 6, 2016, pp. 1-11. https://doi.org/10.1038/s41598-016-0001-8
  3. Maiti, S., Islam, M.R., Uddin, M.A., Afroj, S., Eichhorn, S.J., and Karim, N., "Sustainable fiber-reinforced composites: A Review," Advanced Sustainable Systems, Vol. 6, No. 11, 2022, 2200258.
  4. Borges, C., Chicharo, A., Araujo, A., Silva, J., and Santos, R.M., "Designing of Carbon Fiber-reinforced Polymer (CFRP) Composites for a Second-life in the Aeronautic Industry: Strategies Towards a More Sustainable Future," Frontiers in Materials, Vol. 10, 2023, pp. 1-14.
  5. Nasri, W., Driss, Z., Djebali, R., Lee, K.Y., Park, H.H., Bezazi, A., and Reis, P.N.B., "Thermal Study of Carbon-Fiber-Reinforced Polymer Composites Using Multiscale Modeling," Materials (Basel), Vol. 16, No. 22, 2023, 7233.
  6. Ouyang, Z., Rao, Q., and Peng, X., "Significantly Improving Thermal Conductivity of Carbon Fiber Polymer Composite by Weaving Highly Conductive Films," Composites Part A: Applied Science and Manufacturing, Vol. 163, 2022, 107183.
  7. Wei, J., Liao, M., Ma, A., Chen, Y., Duan, Z., Hou, X., Li, M., Jiang, N., ang Yu, J., "Enhanced Thermal Conductivity of Polydimethylsiloxane Composites with Carbon Fiber," Composites Communications, Vol. 17, 2020, pp. 141-146. https://doi.org/10.1016/j.coco.2019.12.004
  8. Sun, D., "Multi-gating Injection Molding to Enhance the Thermal Conductivity of Carbon Fiber/polysulfone Composite," Polymer Composites, Vol. 38, No. 1, 2017, pp. 185-191. https://doi.org/10.1002/pc.23574
  9. Chen, J., Han, J., and Xu, D., "Thermal and Electrical Properties of the Epoxy Nanocomposites Reinforced with Purified Carbon Nanotubes," Materials Letters, Vol. 246, 2019, pp. 20-23. https://doi.org/10.1016/j.matlet.2019.03.037
  10. Tiwari, M., Billing, B.K., Bedi, H.S., and Agnihotri, P.K., "Quantification of Carbon Nanotube Dispersion and Its Correlation with Mechanical and Thermal Properties of Epoxy Nanocomposites," Journal of Applied Polymer Science, Vol. 137, No. 29, 2020, 48879.
  11. Caradonna, A., "Electrical and Thermal Conductivity of Epoxy-carbon Filler Composites Processed by Calendaring," Materials (Basel), Vol. 12, No. 9, 2019, 1522.
  12. Wang, H., Li, L., Chen, Y., Li, M., Fu, H., Hou, X., Wu, X., Lin, C.T., Jiang, N., and Yu, J., "Efficient Thermal Transport Highway Construction within Epoxy Matrix via Hybrid Carbon Fibers and Alumina Particles," ACS Omega, Vol. 5, No. 2, 2020, pp. 1170-1177. https://doi.org/10.1021/acsomega.9b03465
  13. Kim, M., Sung, D.H., Park, Y.B., and Park, K., "Structural Optimization for Improvement of Thermal Conductivity of Woven Fabric Composites," Composites Research, Vol. 30, No. 1, 2017, pp. 26-34. https://doi.org/10.7234/composres.2017.30.1.026
  14. Lee, W., Kim, S., Sim, H.J., Lee, J.H., An, B.H., Kim, Y.J., Jeong, S.Y., and Shin, H., "Development of Homogenization Databased Transfer Learning Framework to Predict Effective Mechanical Properties and Thermal Conductivity of Foam Structures," Composites Research, Vol. 36, No. 3, 2023, pp. 205-210. https://doi.org/10.7234/COMPOSRES.2023.36.3.205
  15. Mw, J., "Next Generation Quantum Annealing System," 2019.
  16. Adachi, S.H., and Henderson, M.P., "Application of Quantum Annealing to Training of Deep Neural Networks," ArXiv Prepr. ArXiv 1510.06356 (2015).
  17. Boyda, E., Basu, S., Ganguly, S., Michaelis, A., Mukhopadhyay, S., and Nemani, R.R., "Deploying a Quantum Annealing Processor to Detect Tree Cover in Aerial Imagery of California," PLoS One, Vol. 12, No. 2, 2017, pp. 1-22. https://doi.org/10.1371/journal.pone.0172505
  18. Rajak, A., Suzuki, S., and Dutta, A., "Quantum Annealing: An Overview," Philosophical Transactions of the Royal Society A, Vol. 381, No. 2241, 2023.
  19. Hauke, P., Katzgraber, H.G., Lechner, W., Nishimori, H., and Oliver, W.D., "Perspectives of Quantum Annealing: Methods and Implementations," Reports on Progress in Physics, Vol. 83, No. 5, 2020, 054401.
  20. Albash, T., and Lidar, D.A., "Adiabatic Quantum Computation," Reviews of Modern Physics, Vol. 90, No. 1, 2018, 15002.
  21. Imanaka, Y., Anazawa, T., Kumasaka, F., and Jippo, H., "Optimization of the Composition in a Composite Material for Microelectronics Application Using the Ising Model," Scientific Reports,, Vol. 11, No. 1, 2021, pp. 1-7. https://doi.org/10.1038/s41598-021-81243-2
  22. ASTM E1461-13 Standard Test Method for Thermal Diffusivity by the Flash Method, (n.d.).
  23. Lian, T.W., Kondo, A., Akoshima, M., Abe, H., Ohmura, T., Tuan, W.H., and Naito, M., "Rapid Thermal Conductivity Measurement of Porous Thermal Insulation Material by Laser Flash Method," Advanced Powder Technology, Vol. 27, No. 3, 2016, pp. 882-885. https://doi.org/10.1016/j.apt.2016.01.008
  24. Palacios, A., Cong, L., Navarro, M.E., Ding, Y., and Barreneche, C., "Thermal Conductivity Measurement Techniques for Characterizing Thermal Energy Storage Materials - A Review," Renewable and Sustainable Energy Reviews, Vol. 108, 2019, pp. 32-52. https://doi.org/10.1016/j.rser.2019.03.020
  25. Johnston, A.A., "An Integrated Model of the Development of Process-induced Deformation in Autoclave Processing of Composite Structures," Univ. Br. Columbia (April), 1997, pp. 1-367.
  26. Radcliffe, D.J., and Rosenberg, H.M., "The Thermal Conductivity of Glass-fibre and Carbon-fibre/epoxy Composites from 2 to 80 K," Cryogenics (Guildf), Vol. 22, No. 5, 1982, pp. 245-249. https://doi.org/10.1016/0011-2275(82)90122-9
  27. Athanasopoulos, N., Koutsoukis, G., Vlachos, D., and Kostopoulos, V., "Temperature Uniformity Analysis and Development of Open Lightweight Composite Molds Using Carbon Fibers as Heating Elements," Composites Part B: Engineering, Vol. 50, 2013, pp. 279-289. https://doi.org/10.1016/j.compositesb.2013.02.038
  28. He, C., and Xu, J., "Finite Element Analysis of the Thermal Conductivity and the Specific Heat of Carbon Fiber Reinforced Plastic (CFRP) Composites," Proc. - Int. Conf. Artif. Intell. Electromechanical Autom. AIEA 2020 (2020) 771-774.
  29. "MIT Library. 6.777J/2.751J Material Properties Database; MIT Library: Cambridge, MA, USA," (n.d.).
  30. "D-Wave Ocean Documentation, https://docs.ocean.dwavesys.com/en/stable/concepts/index.html," (n.d.).