High flow and pressure in peripheral reagion of Hele-Shaw cell

David_MicroFluidics · May 19, 2017, 10:07

Dear co-foamers,

I am using twoLiquidMixingFoam to simulate laminar, incompressible flow of two miscible liquids in a Hele-Shaw cell.
The cell is open on the edges (BC U: pressureinletOutletVelocity p: totalPressure) and has two connections close to the center with flow directions perpendicular to the orientation of the cell. The cell is filled with species A and species B is injected through on aperture and reaspirated at a higher flow rate through the second neighbouring aperture (BC U: fixedValue, p: fixedFluxPressure). Walls are walls (BC U: noSlip, p: fixedFluxPressure).

When simulating this, I nicely see the confinement of the injected species B and thigs in general look OK, but I get a max in flow and pressure in a region of my cell where actually nothing much happens and flow should be just radially directed inwards. Can anyone please give me a hint on what I might be doing wrong (the Mesh seems to be OK)? Thank you very much in advance!

Best regards,
David

HeleShawTrouble_David.jpg

MeshRect250_10.jpg

PS: one of the pioneering papers regarding this idea confining flows in such an arrangement, in case you are interested http://www.nature.com/nmat/journal/v.../nmat1435.html).

David_MicroFluidics · May 26, 2017, 11:04

Edit: solved the issue by going into fvSolution and switching from nOuterCorrectors = 0 to nOuterCorrectors = 50 (and implementing an abort condition when the residuals become smaller than a defined threshold, so that not all 50 iterations have to be made each cycle). If I understand correctly, with the previous setting I was running my PIMPLE solver in PISO mode, now it's true PIMPLE. Simulations are still running, but stuff looks promising. Fingers crossed

May 19, 2017, 10:07	High flow and pressure in peripheral reagion of Hele-Shaw cell	#1
David_MicroFluidics New Member David Taylor Join Date: Feb 2017 Location: Zürich Posts: 2 Rep Power: 0	Dear co-foamers, I am using twoLiquidMixingFoam to simulate laminar, incompressible flow of two miscible liquids in a Hele-Shaw cell. The cell is open on the edges (BC U: pressureinletOutletVelocity p: totalPressure) and has two connections close to the center with flow directions perpendicular to the orientation of the cell. The cell is filled with species A and species B is injected through on aperture and reaspirated at a higher flow rate through the second neighbouring aperture (BC U: fixedValue, p: fixedFluxPressure). Walls are walls (BC U: noSlip, p: fixedFluxPressure). When simulating this, I nicely see the confinement of the injected species B and thigs in general look OK, but I get a max in flow and pressure in a region of my cell where actually nothing much happens and flow should be just radially directed inwards. Can anyone please give me a hint on what I might be doing wrong (the Mesh seems to be OK)? Thank you very much in advance! Best regards, David HeleShawTrouble_David.jpg MeshRect250_10.jpg PS: one of the pioneering papers regarding this idea confining flows in such an arrangement, in case you are interested http://www.nature.com/nmat/journal/v.../nmat1435.html).

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
BC for Actuator Disk - Flow Vel from High Pressure to Low Pressure	vinguva	OpenFOAM Running, Solving & CFD	2	March 8, 2016 00:46
sonicFoam - pressure driven pipe: flow continuity violation and waveTransmissive BC	Endel	OpenFOAM Running, Solving & CFD	3	September 11, 2014 17:29
reference pressure and compressible flow	bingo10	CFX	0	September 11, 2013 08:32
lid-driven cavity in matlab using BiCGStab	Don456	Main CFD Forum	1	January 19, 2012 16:00
Hydrostatic pressure in 2-phase flow modeling (long)	DS & HB	Main CFD Forum	0	January 8, 2000 16:00

May 26, 2017, 11:04		#2
David_MicroFluidics New Member David Taylor Join Date: Feb 2017 Location: Zürich Posts: 2 Rep Power: 0	Edit: solved the issue by going into fvSolution and switching from nOuterCorrectors = 0 to nOuterCorrectors = 50 (and implementing an abort condition when the residuals become smaller than a defined threshold, so that not all 50 iterations have to be made each cycle). If I understand correctly, with the previous setting I was running my PIMPLE solver in PISO mode, now it's true PIMPLE. Simulations are still running, but stuff looks promising. Fingers crossed