<?xml version="1.0" encoding="utf-8" standalone="yes"?><rss version="2.0" xmlns:atom="http://www.w3.org/2005/Atom"><channel><title>CFD | Joseph Mouallem</title><link>https://josephmouallem.github.io/tag/cfd/</link><atom:link href="https://josephmouallem.github.io/tag/cfd/index.xml" rel="self" type="application/rss+xml"/><description>CFD</description><generator>Wowchemy (https://wowchemy.com)</generator><language>en-us</language><lastBuildDate>Fri, 01 Oct 2021 00:00:00 +0000</lastBuildDate><image><url>https://josephmouallem.github.io/media/icon_hu08dff4d70575caa8b25c1fc7498ce3a4_155912_512x512_fill_lanczos_center_3.png</url><title>CFD</title><link>https://josephmouallem.github.io/tag/cfd/</link></image><item><title>Targeted Particle Delivery via Vortex Ring Reconnection</title><link>https://josephmouallem.github.io/research/vortex-ring-delivery/</link><pubDate>Fri, 01 Oct 2021 00:00:00 +0000</pubDate><guid>https://josephmouallem.github.io/research/vortex-ring-delivery/</guid><description>&lt;h2 id="motivation">Motivation&lt;/h2>
&lt;p>Precise particle targeting is important in manufacturing, propulsion, and medical applications. This work shows how vortex ring dynamics can coherently transport particles to specific wall locations, with design parameters (ring size, Stokes number) controlling delivery accuracy.&lt;/p>
&lt;h2 id="concept">Concept&lt;/h2>
&lt;p>A conceptual model for &lt;strong>targeted particle delivery&lt;/strong> is proposed using controlled
vortex ring reconnection. Particles entrained in the core of a vortex ring are
efficiently transported as the ring advects by self-induction.&lt;/p>
&lt;figure id="figure-left-initial-particle-seeding-within-the-cores-of-the-two-vortex-rings-right-the-reconnection-sequence--first-reconnection-second-reconnection-and-pinch-off--that-redirects-the-particles-toward-the-sidewalls">
&lt;div class="d-flex justify-content-center">
&lt;div class="w-100" >&lt;img alt="Left: initial particle seeding within the cores of the two vortex rings. Right: the reconnection sequence — first reconnection, second reconnection, and pinch off — that redirects the particles toward the sidewalls." srcset="
/research/vortex-ring-delivery/particle-delivery_huaa8523547833b1a2ad240d91967b2ebc_139263_17b3d66b54d88b30a09d305782e132e0.webp 400w,
/research/vortex-ring-delivery/particle-delivery_huaa8523547833b1a2ad240d91967b2ebc_139263_d905aa71e7eaae07affc96b01b66c937.webp 760w,
/research/vortex-ring-delivery/particle-delivery_huaa8523547833b1a2ad240d91967b2ebc_139263_1200x1200_fit_q100_h2_lanczos_3.webp 1200w"
src="https://josephmouallem.github.io/research/vortex-ring-delivery/particle-delivery_huaa8523547833b1a2ad240d91967b2ebc_139263_17b3d66b54d88b30a09d305782e132e0.webp"
width="720"
height="262"
loading="lazy" data-zoomable />&lt;/div>
&lt;/div>&lt;figcaption>
Left: initial particle seeding within the cores of the two vortex rings. Right: the reconnection sequence — first reconnection, second reconnection, and pinch off — that redirects the particles toward the sidewalls.
&lt;/figcaption>&lt;/figure>
&lt;h2 id="mechanism">Mechanism&lt;/h2>
&lt;p>A pair of these particle-transporting vortex rings traveling in the streamwise
direction along parallel trajectories will mutually interact, resulting in a pair
of &lt;strong>vortex reconnection events&lt;/strong>. The reconnection causes a topological change to
the ring, accompanied by a rapid repulsion in the plane perpendicular to the
direction of travel. This effectively transports the particles toward the desired
location on the sidewalls of a ducted flow.&lt;/p>
&lt;figure class="video-figure">
&lt;video
class="video-figure-media"
autoplay loop muted playsinline controls
preload="metadata"
poster="/research/vortex-ring-delivery/vortex-ring-reconnection_poster.jpg">
&lt;source src="https://josephmouallem.github.io/research/vortex-ring-delivery/vortex-ring-reconnection.mp4" type="video/mp4">
Your browser does not support the video tag.
&lt;/video>
&lt;figcaption>Simulation of the two particle-laden vortex rings advecting, interacting, and reconnecting inside the duct.&lt;/figcaption>
&lt;/figure>
&lt;p>In addition to proposing this conceptual model, we identify the dominant physics
of the process and the design considerations required to achieve targeted
delivery.&lt;/p>
&lt;h2 id="particle-inertia-sets-the-delivery">Particle inertia sets the delivery&lt;/h2>
&lt;p>Whether the particles actually follow the rings is governed by the Stokes number.
Low-inertia particles stay locked to the ring cores and are carried coherently
through both reconnections, while high-inertia particles decouple from the
vortex and disperse before reaching the wall.&lt;/p>
&lt;figure id="figure-mean-streamwise-particle-position-left-and-mean-wall-normal-particle-position-right-for-st--01-1-and-10-shaded-bands-show-one-standard-deviation-and-the-dashed-lines-mark-the-two-reconnection-events-at-st--01-and-st--1-the-particles-track-the-rings-and-are-delivered-together-at-st--10-they-lag-and-spread-widely">
&lt;div class="d-flex justify-content-center">
&lt;div class="w-100" >&lt;img alt="Mean streamwise particle position (left) and mean wall-normal particle position (right) for St = 0.1, 1 and 10. Shaded bands show one standard deviation and the dashed lines mark the two reconnection events. At St = 0.1 and St = 1 the particles track the rings and are delivered together; at St = 10 they lag and spread widely." srcset="
/research/vortex-ring-delivery/stokes-number_hue17089a22b4a7f39de24dea8dfa8dedc_159826_5530cf3cd6bf21947fa8be16e44b2706.webp 400w,
/research/vortex-ring-delivery/stokes-number_hue17089a22b4a7f39de24dea8dfa8dedc_159826_3b284a51dc9c2157862bdf82d373f301.webp 760w,
/research/vortex-ring-delivery/stokes-number_hue17089a22b4a7f39de24dea8dfa8dedc_159826_1200x1200_fit_q100_h2_lanczos_3.webp 1200w"
src="https://josephmouallem.github.io/research/vortex-ring-delivery/stokes-number_hue17089a22b4a7f39de24dea8dfa8dedc_159826_5530cf3cd6bf21947fa8be16e44b2706.webp"
width="760"
height="284"
loading="lazy" data-zoomable />&lt;/div>
&lt;/div>&lt;figcaption>
Mean streamwise particle position (left) and mean wall-normal particle position (right) for St = 0.1, 1 and 10. Shaded bands show one standard deviation and the dashed lines mark the two reconnection events. At St = 0.1 and St = 1 the particles track the rings and are delivered together; at St = 10 they lag and spread widely.
&lt;/figcaption>&lt;/figure>
&lt;h2 id="design-considerations">Design considerations&lt;/h2>
&lt;p>The ring radius, relative to the channel width, controls how strongly the rings
repel after reconnection and therefore where on the sidewall the particles
arrive. This makes ring size the primary design parameter for aiming the
delivery.&lt;/p>
&lt;figure id="figure-spreading-angle-θ-after-reconnection-for-three-ring-sizes-with-r-set-by-the-channel-width-w-r--015w-blue-r--0125w-green-and-r--01w-red-larger-rings-separate-more-aggressively-moving-the-impact-point-further-upstream">
&lt;div class="d-flex justify-content-center">
&lt;div class="w-100" >&lt;img alt="Spreading angle θ after reconnection for three ring sizes, with R set by the channel width W: R = 0.15W (blue), R = 0.125W (green) and R = 0.1W (red). Larger rings separate more aggressively, moving the impact point further upstream." srcset="
/research/vortex-ring-delivery/ring-size_hu22716c85b7a01cbc9234f09b70e82e54_206918_0e568b846196ede59c5c4dcc30500bd2.webp 400w,
/research/vortex-ring-delivery/ring-size_hu22716c85b7a01cbc9234f09b70e82e54_206918_9a633b30c6207a713fbdbb9b09efdc8c.webp 760w,
/research/vortex-ring-delivery/ring-size_hu22716c85b7a01cbc9234f09b70e82e54_206918_1200x1200_fit_q100_h2_lanczos_3.webp 1200w"
src="https://josephmouallem.github.io/research/vortex-ring-delivery/ring-size_hu22716c85b7a01cbc9234f09b70e82e54_206918_0e568b846196ede59c5c4dcc30500bd2.webp"
width="760"
height="485"
loading="lazy" data-zoomable />&lt;/div>
&lt;/div>&lt;figcaption>
Spreading angle θ after reconnection for three ring sizes, with R set by the channel width W: R = 0.15W (blue), R = 0.125W (green) and R = 0.1W (red). Larger rings separate more aggressively, moving the impact point further upstream.
&lt;/figcaption>&lt;/figure>
&lt;h2 id="reference">Reference&lt;/h2>
&lt;p>Mouallem, J., Daryan, H., Wawryk, J., Pan, Z., and Hickey, J.-P.: &lt;em>Targeted
particle delivery via vortex ring reconnection&lt;/em>, &lt;strong>Physics of Fluids&lt;/strong>, 33(10),
2021.
&lt;a href="https://doi.org/10.1063/5.0066443" target="_blank" rel="noopener">https://doi.org/10.1063/5.0066443&lt;/a>&lt;/p></description></item><item><title>Macro-Scale Effects on Sub-Grid Closures in Gas-Solid Riser Flows</title><link>https://josephmouallem.github.io/research/multiphase-riser-flows/</link><pubDate>Fri, 01 Jun 2018 00:00:00 +0000</pubDate><guid>https://josephmouallem.github.io/research/multiphase-riser-flows/</guid><description>&lt;h2 id="motivation">Motivation&lt;/h2>
&lt;p>Sub-grid closures for two-fluid models are traditionally derived at a single scale, but this work demonstrates that flow topology at the system scale strongly affects closure accuracy. Accounting for these effects is essential for predictive industrial riser flow simulations.&lt;/p>
&lt;h2 id="background">Background&lt;/h2>
&lt;p>Filtered two-fluid formulations of gas-solid fluidized flows require closure
models to deal with sub-grid filtered parameters. These closures are derived by
filtering the results of meso-scale highly resolved simulations (HRS) with
two-fluid modeling, and then applying them on the coarse large-scale simulation
(LSS) grid.&lt;/p>
&lt;figure id="figure-closure-strategy-a-highly-resolved-simulation-with-the-two-fluid-model-and-microscopic-closures-is-filtered-to-provide-sub-grid-closures-for-the-coarse-large-scale-filtered-two-fluid-model">
&lt;div class="d-flex justify-content-center">
&lt;div class="w-100" >&lt;img alt="Closure strategy: a highly resolved simulation with the two-fluid model and microscopic closures is filtered to provide sub-grid closures for the coarse, large-scale filtered two-fluid model." srcset="
/research/multiphase-riser-flows/filtered-parameters_huf04d603bef027e576909f98379d5017e_97768_e7c97ecc5a354493e5edd076b6658052.webp 400w,
/research/multiphase-riser-flows/filtered-parameters_huf04d603bef027e576909f98379d5017e_97768_c6e40025e761b84a7546c3f8d7a63ad3.webp 760w,
/research/multiphase-riser-flows/filtered-parameters_huf04d603bef027e576909f98379d5017e_97768_1200x1200_fit_q100_h2_lanczos_3.webp 1200w"
src="https://josephmouallem.github.io/research/multiphase-riser-flows/filtered-parameters_huf04d603bef027e576909f98379d5017e_97768_e7c97ecc5a354493e5edd076b6658052.webp"
width="760"
height="391"
loading="lazy" data-zoomable />&lt;/div>
&lt;/div>&lt;figcaption>
Closure strategy: a highly resolved simulation with the two-fluid model and microscopic closures is filtered to provide sub-grid closures for the coarse, large-scale filtered two-fluid model.
&lt;/figcaption>&lt;/figure>
&lt;h2 id="the-question">The question&lt;/h2>
&lt;p>Trusting in scale separation, the correlation of filtered parameters has
traditionally been performed against meso-scale filtered data only, disregarding
any macro-scale effects. &lt;strong>In this work, the correctness of that practice is
tested — and it fails.&lt;/strong>&lt;/p>
&lt;h2 id="approach">Approach&lt;/h2>
&lt;p>Two macro-scale parameters associated with flow topology are considered for their
effects on the relevant filtered parameters: the average solid volume fraction and
the average gas Reynolds number. Highly resolved simulations are filtered while
holding each of these macro-scale parameters constant at various levels. The
interest is directed toward the dilute conditions typical of riser flows.&lt;/p>
&lt;figure id="figure-instantaneous-solid-volume-fraction-for-increasing-solids-loading-left-to-right-and-increasing-gas-reynolds-number-top-to-bottom-the-cluster-structure--and-therefore-the-sub-grid-closure--depends-strongly-on-both-macro-scale-parameters">
&lt;div class="d-flex justify-content-center">
&lt;div class="w-100" >&lt;img alt="Instantaneous solid volume fraction for increasing solids loading (left to right) and increasing gas Reynolds number (top to bottom). The cluster structure — and therefore the sub-grid closure — depends strongly on both macro-scale parameters." srcset="
/research/multiphase-riser-flows/riser-flow_hu4db40483cf06ee0181befb9bdcc0261d_242512_76b7384a04a339bf73435483702f7850.webp 400w,
/research/multiphase-riser-flows/riser-flow_hu4db40483cf06ee0181befb9bdcc0261d_242512_01cdd1b61e40386ad3edc6e5381e535a.webp 760w,
/research/multiphase-riser-flows/riser-flow_hu4db40483cf06ee0181befb9bdcc0261d_242512_1200x1200_fit_q100_h2_lanczos_3.webp 1200w"
src="https://josephmouallem.github.io/research/multiphase-riser-flows/riser-flow_hu4db40483cf06ee0181befb9bdcc0261d_242512_76b7384a04a339bf73435483702f7850.webp"
width="490"
height="760"
loading="lazy" data-zoomable />&lt;/div>
&lt;/div>&lt;figcaption>
Instantaneous solid volume fraction for increasing solids loading (left to right) and increasing gas Reynolds number (top to bottom). The cluster structure — and therefore the sub-grid closure — depends strongly on both macro-scale parameters.
&lt;/figcaption>&lt;/figure>
&lt;h2 id="the-drag-closure-is-not-scale-separated">The drag closure is not scale separated&lt;/h2>
&lt;p>The drag coefficient correction H, the single most influential sub-grid term, is
the clearest evidence of the failure of scale separation. Correlating H against
the meso-scale filtered variables alone leaves a systematic spread that is set
entirely by the macro-scale state of the flow.&lt;/p>
&lt;figure id="figure-drag-coefficient-correction-h-against-the-filtered-solid-volume-fraction-a-effect-of-the-domain-average-gas-reynolds-number-at-fixed-solids-loading-b-effect-of-the-domain-average-solid-volume-fraction-at-fixed-reynolds-number-curves-that-should-collapse-if-scale-separation-held-instead-fan-out-by-a-factor-of-two-or-more">
&lt;div class="d-flex justify-content-center">
&lt;div class="w-100" >&lt;img alt="Drag coefficient correction H against the filtered solid volume fraction. (a) Effect of the domain average gas Reynolds number at fixed solids loading; (b) effect of the domain average solid volume fraction at fixed Reynolds number. Curves that should collapse if scale separation held instead fan out by a factor of two or more." srcset="
/research/multiphase-riser-flows/drag-correction_huc01044e791deef848fd9405cab75ea4d_134539_3385d5cddfa722f869717c5ce958575e.webp 400w,
/research/multiphase-riser-flows/drag-correction_huc01044e791deef848fd9405cab75ea4d_134539_cf999186d89f7972b3dfc784723b59f0.webp 760w,
/research/multiphase-riser-flows/drag-correction_huc01044e791deef848fd9405cab75ea4d_134539_1200x1200_fit_q100_h2_lanczos_3.webp 1200w"
src="https://josephmouallem.github.io/research/multiphase-riser-flows/drag-correction_huc01044e791deef848fd9405cab75ea4d_134539_3385d5cddfa722f869717c5ce958575e.webp"
width="760"
height="273"
loading="lazy" data-zoomable />&lt;/div>
&lt;/div>&lt;figcaption>
Drag coefficient correction H against the filtered solid volume fraction. (a) Effect of the domain average gas Reynolds number at fixed solids loading; (b) effect of the domain average solid volume fraction at fixed Reynolds number. Curves that should collapse if scale separation held instead fan out by a factor of two or more.
&lt;/figcaption>&lt;/figure>
&lt;h2 id="the-same-holds-for-the-stress-closures">The same holds for the stress closures&lt;/h2>
&lt;p>The effect is not limited to drag. The filtered solid pressure, which closes the
solid-phase momentum equation, shifts by roughly an order of magnitude across the
range of macro-scale conditions at an otherwise identical filtered state.&lt;/p>
&lt;figure id="figure-dimensionless-filtered-solid-pressure-against-the-filtered-solid-volume-fraction-a-varying-the-domain-average-gas-reynolds-number-b-varying-the-domain-average-solid-volume-fraction-the-vertical-spread-at-fixed-filtered-state-is-the-macro-scale-signature-that-traditional-closures-ignore">
&lt;div class="d-flex justify-content-center">
&lt;div class="w-100" >&lt;img alt="Dimensionless filtered solid pressure against the filtered solid volume fraction. (a) Varying the domain average gas Reynolds number; (b) varying the domain average solid volume fraction. The vertical spread at fixed filtered state is the macro-scale signature that traditional closures ignore." srcset="
/research/multiphase-riser-flows/filtered-solid-pressure_hu35ff66501e13845ed4c67b4ec7342ad0_126256_654a40ad182d6f9a2843c886c7f955e1.webp 400w,
/research/multiphase-riser-flows/filtered-solid-pressure_hu35ff66501e13845ed4c67b4ec7342ad0_126256_66a8bb822df64b7f8e0788e19249cd8c.webp 760w,
/research/multiphase-riser-flows/filtered-solid-pressure_hu35ff66501e13845ed4c67b4ec7342ad0_126256_1200x1200_fit_q100_h2_lanczos_3.webp 1200w"
src="https://josephmouallem.github.io/research/multiphase-riser-flows/filtered-solid-pressure_hu35ff66501e13845ed4c67b4ec7342ad0_126256_654a40ad182d6f9a2843c886c7f955e1.webp"
width="760"
height="266"
loading="lazy" data-zoomable />&lt;/div>
&lt;/div>&lt;figcaption>
Dimensionless filtered solid pressure against the filtered solid volume fraction. (a) Varying the domain average gas Reynolds number; (b) varying the domain average solid volume fraction. The vertical spread at fixed filtered state is the macro-scale signature that traditional closures ignore.
&lt;/figcaption>&lt;/figure>
&lt;h2 id="conclusion">Conclusion&lt;/h2>
&lt;p>Results show that &lt;strong>both&lt;/strong> macro-scale parameters should be accounted for in
sub-grid correlations if higher accuracy is to be achieved.&lt;/p>
&lt;h2 id="reference">Reference&lt;/h2>
&lt;p>Mouallem, J., Chavez-Cussy, N., Niaki, S. R. A., Milioli, C. C., and Milioli,
F. E.: &lt;em>On the effects of the flow macro-scale over meso-scale filtered parameters
in gas-solid riser flows&lt;/em>, &lt;strong>Chemical Engineering Science&lt;/strong>, 182, 200-211, 2018.
&lt;a href="https://doi.org/10.1016/j.ces.2018.02.039" target="_blank" rel="noopener">https://doi.org/10.1016/j.ces.2018.02.039&lt;/a>&lt;/p></description></item><item><title>Targeted particle delivery via vortex ring reconnection</title><link>https://josephmouallem.github.io/publication/vortex-ring-particle-delivery/</link><pubDate>Fri, 01 Oct 2021 00:00:00 +0000</pubDate><guid>https://josephmouallem.github.io/publication/vortex-ring-particle-delivery/</guid><description/></item><item><title>Macro-scale effects over filtered and residual stresses in gas-solid riser flows</title><link>https://josephmouallem.github.io/publication/filtered-residual-stresses-riser/</link><pubDate>Sat, 23 Feb 2019 00:00:00 +0000</pubDate><guid>https://josephmouallem.github.io/publication/filtered-residual-stresses-riser/</guid><description/></item><item><title>On the effects of the flow macro-scale over meso-scale filtered parameters in gas-solid riser flows</title><link>https://josephmouallem.github.io/publication/macro-scale-riser-flows/</link><pubDate>Fri, 01 Jun 2018 00:00:00 +0000</pubDate><guid>https://josephmouallem.github.io/publication/macro-scale-riser-flows/</guid><description/></item></channel></rss>