Matching Items (5)
Filtering by
- All Subjects: Turbulence
- Creators: Kim, Jeonglae
Description
The current work aims to understand the influence of particles on scalar transport in particle-laden turbulent jets using point-particle direct numerical simulations (DNS). Such turbulence phenomena are observed in many applications, such as aircraft and rocket engines (e.g., engines operating in dusty environments and when close to the surface) and geophysical flows (sediment-laden rivers discharging nutrients into the oceans), etc.This thesis looks at systematically understanding the fundamental interplay between (1) fluid turbulence, (2) inertial particles, and (3) scalar transport. This work considers a temporal jet of Reynolds number of 5000 filled with the point-particles and the influence of Stokes number (St). Three Stokes numbers, St = 1, 7.5, and 20, were considered for the current work. The simulations were solved using the NGA solver, which solves the Navier-Stokes, advection-diffusion, and particle transport equations.
The statistical analysis of the mean and turbulence quantities, along with the Reynolds stresses, are estimated for the fluid and particle phases throughout the domain. The observations do not show a significant influence of St in the mean flow evolution of fluid, scalar, and particle phases. The scalar mixture fraction variance and the turbulent kinetic energy (TKE) increase slightly for the St = 1 case, compared to the particle-free and higher St cases, indicating that an optimal St exists for which the scalar variation increases. The current preliminary study establishes that the scalar variance is influenced by particles under the optimal particle St. Directions for future studies based on the current observations are presented.
ContributorsPaturu, Venkata Sai Sushant (Author) / Pathikonda, Gokul (Thesis advisor) / Kasbaoui, Mohamed Houssem (Committee member) / Kim, Jeonglae (Committee member) / Prabhakaran, Prasanth (Committee member) / Arizona State University (Publisher)
Created2023
Description
This dissertation presents a volume filtering framework to solve particle-laden flows. Particle-laden flows are studied, employing the well-established Euler-Lagrange method, using the point-particle approximation. This approach requires the filter width to be much larger than the particle diameter. The method assumes that the particle is smaller than the Kolmogorov length scale. This thesis investigates how inertial particles at semi-dilute volume fractions modulate the flow characteristics for particles smaller than 1 in wall units, when dispersed within wall-bounded channel flows at friction Reynolds number of 180. The simulations are performed with 4 way coupling in order to account for high local concentration of particles, to capture mechanisms such as turbophoresis and preferential concentration. We show that drag attenuation or augmentation is determined by the particle inertia. As particle size is increased greater than 1 in wall units, the regime becomes finite-sized, requiring an interface-resolved description. To do this a novel Immersed Boundaries (IB) framework based on the concept of volume-filtering called the Volume-Filtered Immersed Boundary (VF-IB) method is presented. Transport equations are obtained by volume-filtering the Navier-Stokes equation and accounting for the stresses at the solid-fluid interface. Boundary conditions are transformed into bodyforces that appear as surface integrals on the right hand side of the filtered equation. The approach requires the filter width to be much smaller than the particle diameter in order to accurately resolve the interfacial dynamics. Several canonical tests are conducted for both stationary and moving immersed solids and report comparable results to the experimental and/or body-fitted simulations. Keep in mind, the VF-IB method reverts back to the Euler-Lagrange formulation if the filter width is significantly greater than the particle diameter. An artifact of volume-filtering is the emergence of unclosed terms we define as the sub-filter scale term. In order to characterize the contribution of this term on the solution, a more simpler case of a 2-D varying coefficient hyperbolic equation that has an exact solution is looked into. It is observed that the sub-filter scale term scales inversely with the square of the filter width. For fine interface resolution (i.e. small filter width), this value can be ignored with negligible effect to the accuracy of the numerical solution. However for coarse interface resolution (i.e. large filter width), including the sub-filter scale term significantly increases the accuracy of the numerical solution
ContributorsDave, Himanshu (Author) / Kasbaoui, Mohamed Houssem (Thesis advisor) / Herrmann, Marcus (Thesis advisor) / Dahm, Werner (Committee member) / Kim, Jeonglae (Committee member) / Lopez, Juan (Committee member) / Arizona State University (Publisher)
Created2024
Description
Advancements to a dual scale Large Eddy Simulation (LES) modeling approach for immiscible turbulent phase interfaces are presented. In the dual scale LES approach, a high resolution auxiliary grid, used to capture a fully resolved interface geometry realization, is linked to an LES grid that solves the filtered Navier-Stokes equations. Exact closure of the sub-filter interface terms is provided by explicitly filtering the fully resolved quantities from the auxiliary grid. Reconstructing a fully resolved velocity field to advance the phase interface requires modeling several sub-filter effects, including shear and accelerational instabilities and phase change. Two sub-filter models were developed to generate these sub-filter hydrodynamic instabilities: an Orr-Sommerfeld model and a Volume-of-Fluid (VoF) vortex sheet method. The Orr-Sommerfeld sub-filter model was found to be incompatible with the dual scale approach, since it is unable to generate interface rollup and a process to separate filtered and sub-filter scales could not be established. A novel VoF vortex sheet method was therefore proposed, since prior vortex methods have demonstrated interface rollup and following the LES methodology, the vortex sheet strength could be decomposed into its filtered and sub-filter components. In the development of the VoF vortex sheet method, it was tested with a variety of classical hydrodynamic instability problems, compared against prior work and linear theory, and verified using Direct Numerical Simulations (DNS). An LES consistent approach to coupling the VoF vortex sheet with the LES filtered equations is presented and compared against DNS. Finally, a sub-filter phase change model is proposed and assessed in the dual scale LES framework with an evaporating interface subjected to decaying homogeneous isotropic turbulence. Results are compared against DNS and the interplay between surface tension forces and evaporation are discussed.
ContributorsGoodrich, Austin Chase (Author) / Herrmann, Marcus (Thesis advisor) / Dahm, Werner (Committee member) / Kim, Jeonglae (Committee member) / Huang, Huei-Ping (Committee member) / Kostelich, Eric (Committee member) / Arizona State University (Publisher)
Created2023
Description
Autonomic closure is a new general methodology for subgrid closures in large eddy simulations that circumvents the need to specify fixed closure models and instead allows a fully- adaptive self-optimizing closure. The closure is autonomic in the sense that the simulation itself determines the optimal relation at each point and time between any subgrid term and the variables in the simulation, through the solution of a local system identification problem. It is based on highly generalized representations of subgrid terms having degrees of freedom that are determined dynamically at each point and time in the simulation. This can be regarded as a very high-dimensional generalization of the dynamic approach used with some traditional prescribed closure models, or as a type of “data-driven” turbulence closure in which machine- learning methods are used with internal training data obtained at a test-filter scale at each point and time in the simulation to discover the local closure representation.
In this study, a priori tests were performed to develop accurate and efficient implementations of autonomic closure based on particular generalized representations and parameters associated with the local system identification of the turbulence state. These included the relative number of training points and bounding box size, which impact computational cost and generalizability of coefficients in the representation from the test scale to the LES scale. The focus was on studying impacts of these factors on the resulting accuracy and efficiency of autonomic closure for the subgrid stress. Particular attention was paid to the associated subgrid production field, including its structural features in which large forward and backward energy transfer are concentrated.
More than five orders of magnitude reduction in computational cost of autonomic closure was achieved in this study with essentially no loss of accuracy, primarily by using efficient frame-invariant forms for generalized representations that greatly reduce the number of degrees of freedom. The recommended form is a 28-coefficient representation that provides subgrid stress and production fields that are far more accurate in terms of structure and statistics than are traditional prescribed closure models.
In this study, a priori tests were performed to develop accurate and efficient implementations of autonomic closure based on particular generalized representations and parameters associated with the local system identification of the turbulence state. These included the relative number of training points and bounding box size, which impact computational cost and generalizability of coefficients in the representation from the test scale to the LES scale. The focus was on studying impacts of these factors on the resulting accuracy and efficiency of autonomic closure for the subgrid stress. Particular attention was paid to the associated subgrid production field, including its structural features in which large forward and backward energy transfer are concentrated.
More than five orders of magnitude reduction in computational cost of autonomic closure was achieved in this study with essentially no loss of accuracy, primarily by using efficient frame-invariant forms for generalized representations that greatly reduce the number of degrees of freedom. The recommended form is a 28-coefficient representation that provides subgrid stress and production fields that are far more accurate in terms of structure and statistics than are traditional prescribed closure models.
ContributorsKshitij, Abhinav (Author) / Dahm, Werner J.A. (Thesis advisor) / Herrmann, Marcus (Committee member) / Hamlington, Peter E (Committee member) / Peet, Yulia (Committee member) / Kim, Jeonglae (Committee member) / Arizona State University (Publisher)
Created2019
Description
Autonomic closure is a recently-proposed subgrid closure methodology for large eddy simulation (LES) that replaces the prescribed subgrid models used in traditional LES closure with highly generalized representations of subgrid terms and solution of a local system identification problem that allows the simulation itself to determine the local relation between each subgrid term and the resolved variables at every point and time. The present study demonstrates, for the first time, practical LES based on fully dynamic implementation of autonomic closure for the subgrid stress and the subgrid scalar flux. It leverages the inherent computational efficiency of tensorally-correct generalized representations in terms of parametric quantities, and uses the fundamental representation theory of Smith (1971) to develop complete and minimal tensorally-correct representations for the subgrid stress and scalar flux. It then assesses the accuracy of these representations via a priori tests, and compares with the corresponding accuracy from nonparametric representations and from traditional prescribed subgrid models. It then assesses the computational stability of autonomic closure with these tensorally-correct parametric representations, via forward simulations with a high-order pseudo-spectral code, including the extent to which any added stabilization is needed to ensure computational stability, and compares with the added stabilization needed in traditional closure with prescribed subgrid models. Further, it conducts a posteriori tests based on forward simulations of turbulent conserved scalar mixing with the same pseudo-spectral code, in which velocity and scalar statistics from autonomic closure with these representations are compared with corresponding statistics from traditional closure using prescribed models, and with corresponding statistics of filtered fields from direct numerical simulation (DNS). These comparisons show substantially greater accuracy from autonomic closure than from traditional closure. This study demonstrates that fully dynamic autonomic closure is a practical approach for LES that requires accuracy even at the smallest resolved scales.
ContributorsStallcup, Eric Warren (Author) / Dahm, Werner J.A. (Thesis advisor) / Herrmann, Marcus (Committee member) / Calhoun, Ronald (Committee member) / Kim, Jeonglae (Committee member) / Kostelich, Eric J. (Committee member) / Arizona State University (Publisher)
Created2020