We simulate diffusion in multimaterial systems with a cell-centered Eulerian mesh in two dimensions.
A system with immiscible fluids contains sharp interfaces. An Eulerian mesh is fixed in space and does not move with the material. Therefore, cells with an interface contain multiple fluids; these are known as mixed cells.
The treatment of mixed cells can vary in computational cost and accuracy. In some cases, the primary source of inaccuracy can be attributed to approximations made in modeling the mixed cells. This thesis focuses on the treatment of mixed cells based on the diffusion approximation of the transport equation. We introduce five subgrid, mixed-cell models. Two models have a single temperature for each cell, while the other three allow a separate temperature for each phase. The single-temperature models are implemented using the Support-Operators Method, which is derived herein. The first single-temperature model utilizes an effective tensor diffusivity that distinguishes diffusion tangent and normal to the interface. The second single-temperature model specifies a unique diffusivity in each corner of a mixed cell, which is effectively a mesh refinement of the mixed cell. The three multi-temperature models have increasingly accurate levels of approximation of the flux: (i) flux is calculated between cell-centers for each phase, (ii) flux is calculated between the centroid of each phase, and (iii) flux normal to an interface is calculated between centroids of each phase. The physical interpretations of these models are: (i) each phase occupies the entire cell, (ii) oblique flux is continuous, (iii) only normal flux is continuous. The standard approximation, using the harmonic mean of the diffusivities present in a mixed cell as an effective diffusivity, is also tested for comparison.
We also derive two time-dependent analytical solutions for diffusion in a two-phase system, in both one and two dimensions. With the standard model as a reference point, the accuracy of the new models is quantified, and the convergence rates of the error are determined between pairs of spatial resolutions for the two problems with analytical solutions. Simulations of multiphysics and multimaterial phenomenon may benefit from increased mixed-cell fidelity achieved in this dissertation.PHDApplied PhysicsUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/107150/1/leftynm_1.pd
