<p>Computational homogenization of particle-matrix composites with material interfaces is considered, whereby linear transient diffusion driven by a chemical potential is used as a model problem. Due to the models linearity, it is beneficial to assume micro-stationarity in order to provide direct upscaling and avoid excessive computational cost. First order homogenization with a single macroscale chemical potential is taken as the most basic approach; however, the accuracy is negatively affected whenever relevant micro-transient effects can not be captured using the stationary sub-scale problem. To improve the accuracy while still upscaling from a stationary sub-scale problem, different formulations based on dual macroscale potentials, one for each phase, are proposed and investigated in this paper. As to the prolongation order within the particles and matrix phase, respectively, two types are considered: constant-linear and linear-linear. Most importantly, for the case of linear prolongation, different ways of defining the macroscale variables (acting as loading on the RVE-problem) in terms of suitable measures of the chemical potential can be envisioned: (1) averaging of <InlineEquation ID="IEq1"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="466_2025_2605_Article_IEq1.gif" Format="GIF" Height="14" Rendition="HTML" Resolution="72" Type="Linedraw" Width="25" /> </InlineMediaObject> <EquationSource Format="TEX">\(0\text {th}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>0</mn> <mtext>th</mtext> </mrow> </math></EquationSource> </InlineEquation> and <InlineEquation ID="IEq2"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="466_2025_2605_Article_IEq2.gif" Format="GIF" Height="14" Rendition="HTML" Resolution="72" Type="Linedraw" Width="22" /> </InlineMediaObject> <EquationSource Format="TEX">\(1\text {st}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>1</mn> <mtext>st</mtext> </mrow> </math></EquationSource> </InlineEquation> gradient of the potential, (2) averaging of <InlineEquation ID="IEq3"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="466_2025_2605_Article_IEq1.gif" Format="GIF" Height="14" Rendition="HTML" Resolution="72" Type="Linedraw" Width="25" /> </InlineMediaObject> <EquationSource Format="TEX">\(0\text {th}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>0</mn> <mtext>th</mtext> </mrow> </math></EquationSource> </InlineEquation> and <InlineEquation ID="IEq4"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="466_2025_2605_Article_IEq2.gif" Format="GIF" Height="14" Rendition="HTML" Resolution="72" Type="Linedraw" Width="22" /> </InlineMediaObject> <EquationSource Format="TEX">\(1\text {st}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>1</mn> <mtext>st</mtext> </mrow> </math></EquationSource> </InlineEquation> moment of the potential. The pros and cons of the different approaches were assessed in a numerical study and compared to a reference solution from Direct Numerical Simulation (DNS) for an example problem. It was concluded that the moment-based linear-linear method was the only one that could match the DNS solution for all considered material parameters. However, for sufficiently large interface resistance, leading to a more pronounced potential jump across the interfaces, the constant-linear prolongation gave comparable results.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Variationally consistent homogenization of diffusion in particle composites with material interfaces using dual macroscale chemical potentials

  • David René Rollin,
  • Fredrik Larsson,
  • Kenneth Runesson,
  • Ralf Jänicke

摘要

Computational homogenization of particle-matrix composites with material interfaces is considered, whereby linear transient diffusion driven by a chemical potential is used as a model problem. Due to the models linearity, it is beneficial to assume micro-stationarity in order to provide direct upscaling and avoid excessive computational cost. First order homogenization with a single macroscale chemical potential is taken as the most basic approach; however, the accuracy is negatively affected whenever relevant micro-transient effects can not be captured using the stationary sub-scale problem. To improve the accuracy while still upscaling from a stationary sub-scale problem, different formulations based on dual macroscale potentials, one for each phase, are proposed and investigated in this paper. As to the prolongation order within the particles and matrix phase, respectively, two types are considered: constant-linear and linear-linear. Most importantly, for the case of linear prolongation, different ways of defining the macroscale variables (acting as loading on the RVE-problem) in terms of suitable measures of the chemical potential can be envisioned: (1) averaging of \(0\text {th}\) 0 th and \(1\text {st}\) 1 st gradient of the potential, (2) averaging of \(0\text {th}\) 0 th and \(1\text {st}\) 1 st moment of the potential. The pros and cons of the different approaches were assessed in a numerical study and compared to a reference solution from Direct Numerical Simulation (DNS) for an example problem. It was concluded that the moment-based linear-linear method was the only one that could match the DNS solution for all considered material parameters. However, for sufficiently large interface resistance, leading to a more pronounced potential jump across the interfaces, the constant-linear prolongation gave comparable results.