Search for: All records

Motivated by applications to fluid flows with conjugate heat transfer and electrokinetic effects, we propose a direct forcing immersed boundary method for simulating general, discontinuous, Dirichlet and Robin conditions at the interface between two materials. In comparison to existing methods, our approach uses smaller stencils and accommodates complex geometries with sharp corners. The method is built on the concept of a “forcing pair,” defined as two grid points that are adjacent to each other, but on opposite sides of an interface. For 2D problems this approach can simultaneously enforce discontinuous Dirichlet and Robin conditions using a six-point stencil at one of the forcing points, and a 12-point stencil at the other. In comparison, prior work requires up to 14-point stencils at both points. We also propose two methods of accommodating surfaces with sharp corners. The first locally reduces stencils in sharp corners. The second uses the signed distance function to globally smooth all corners on a surface. The smoothing is defined to recover the actual corners as the grid is refined. We verify second-order spatial accuracy of our proposed methods by comparing to manufactured solutions to the Poisson equation with challenging dis- continuous fields across immersed surfaces. Next, to explore the performance of our method for simulating fluid flows with conjugate heat transport, we couple our method to the incompressible Navier–Stokes and continuity equations using a finite-volume projection method. We verify the spatial-temporal accuracy of the solver using manufactured solutions and an analytical solution for circular Couette flow with conjugate heat transfer. Finally, to demonstrate that our method can model moving surfaces, we simulate fluid flow and conjugate heat transport between a stationary cylinder and a rotating ellipse or square. 

We present a novel method for simulating unsteady, variable density, fluid flows in membrane desalination systems. By assuming the density varies only with concentration and temperature, the scheme decouples the solution of the governing equations into two sequential blocks. The first solves the governing equations for the temperature and concentration fields, which are used to compute all thermophysical properties. The second block solves the conservation of mass and momentum equations for the velocity and pressure. We show that this is computationally more efficient than schemes that iterate over the full coupled equations in one block. We verify that the method achieves second-order spatial–temporal accuracy, and we use the method to investigate buoyancy-driven convection in a desalination process called vacuum membrane distillation. Specifically, we show that with gravity properly oriented, variations in temperature and concentration can trigger a double-diffusive instability that enhances mixing and improves water recovery. We also show that the instability can be strengthened by providing external heating. 

Membrane distillation (MD) is a thermally-driven desalination process that can treat hypersaline brines. Considerable MD literature has focused on mitigating temperature and concentration polarization. This literature largely neglects that temperature and concentration polarization increase the feed density near the membrane. With gravity properly oriented, this increase in density could trigger buoyancy-driven convection and increase permeate production. Convection could also be strengthened by heating the feed channel wall opposite the membrane. To investigate that possibility, we perform a series of experiments using a plate-and-frame direct contact MD system with an active membrane area of 300 cm2 and a feed channel wall heated using a resistive heater. The experiments measure the average transmembrane permeate flux for two gravitational orientations, feed Reynolds numbers between 128 and 1128, and wall heat fluxes up to 12 kW/m2. The results confirm that with gravity properly oriented, wall-heating can trigger buoyancy-driven convection for a wide range of feed Reynolds numbers, and increase permeate production between roughly 20 and 130 %. We estimate, however, that at high Reynolds numbers (𝑅𝑒 > 800), more than 70 % of the wall heat is carried out of the MD system by the feed flow, without contributing to permeate production. This suggests the need for longer membranes and heat recovery steps in any future practical implementation.