skip to main content

Explore Research Products in the NSF-PAR

Title: A Dual Scale Approach to Modeling Sub-Filter Velocities due to Shear-Induced Instabilities
A method to compute sub-filter velocities due to shear induced instabilities on a liquid-gas interface for use in a dual scale LES-DNS model is presented. The method reconstructs the sub-filter velocity field as the sum of a prescribed base velocity profile and a perturbation velocity field determined by the Orr-Sommerfeld equations. The base velocity profile is approximated as an error function appropriately scaled with flow parameters, and the perturbation velocity field is computed by solving the Orr-Sommerfeld equations with appropriate boundary and interface conditions. The perturbation velocities of the Orr-Sommerfeld equations are expanded into Chebyshev polynomials to create a linear eigenvalue problem as outlined by Schmid and Henningson (2001). Finally the eigenvalue problem is solved using a standard linear algebra package and used to evaluate the perturbation velocities. The Chebyshev method is tested under a variety of flow parameters and initial interface disturbances. Results are presented and compared against prior literature and asymptotic solutions.  more » « less
Award ID(s):
1803657
NSF-PAR ID:
10281846
Author(s) / Creator(s):
;
Date Published:
Journal Name:
ILASS-Americas 31st Annual Conference on Liquid Atomization and Spray Systems
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ( , ILASS-Americas 31st Annual Conference on Liquid Atomization and Spray Systems)
    null (Ed.)
    A method to predict sub-filter shear-induced velocities on a liquid-gas phase interface for use in a dual scale LES model is presented. The method reconstructs the sub-filter velocity field in the vicinity of the interface by introducing a vortex sheet at the interface. The vortex sheet is transported by an unsplit geometric volume and surface area advection scheme with a Piecewise Linear Interface Construction (PLIC) representation of the material interface. At each step and desired location the shear-induced velocities can be calculated by integrating the vortex sheet and other relevant quantities over the liquid-gas interface with the sub-grid velocity reconstruction limited to a small number of cells near the phase interface. The vortex sheet method is tested and compared against prior literature. 
    more » « less
  2. ; ( , ICLASS 2021, 15th Triennial International Conference on Liquid Atomization and Spray Systems)
    null (Ed.)
    A method to predict sub-filter shear-induced velocities on a liquid-gas phase interface for use in a dual scale LES model is presented. The method reconstructs the sub-filter velocity field in the vicinity of the interface by introducing a vortex sheet at the interface. The vortex sheet is transported by an unsplit geometric volume and surface area advection scheme with a Piece- wise Linear Interface Construction (PLIC) representation of the material interface. At each step and desired location the shear-induced velocities can be calculated by integrating the vortex sheet and other relevant quantities over the liquid-gas surface with the sub-grid velocity recon- struction limited to a small number of cells near the phase interface. The vortex sheet method is tested and compared against prior literature. 
    more » « less
  3. ; ( , Proceedings of the 32nd Annual Conference on Liquid Atomization and Spray Systems)
    A method to predict sub-filter shear-induced velocities on a liquid-gas phase interface for use in a dual scale LES model is presented and compared against prior work on Vortex Sheet methods. The method reconstructs the sub-filter velocity field in the vicinity of the interface by employing a vortex sheet at the interface location. The vortex sheet is transported by an unsplit geometric volume and surface area advection scheme with a Piecewise Linear Interface Construction (PLIC) representation of the material interface. At each step, the vorticity field is constructed by evaluating a volume integral of the vortex sheet and a numerical spreading parameter near the liquid-gas interface. A Poisson equation can then be constructed and solved for the vector potential; the self-induced velocities due to the vortex sheet are subsequently evaluated from the vector potential. The described vortex sheet method is tested and compared against prior literature. 
    more » « less
  4. ; ; ; ; ; ; ( , Earth Surface Processes and Landforms)
    Abstract

    Decimeter‐scale early‐stage aeolian bedforms represent topographic features that differ notably from their mature dune counterparts, with nascent forms exhibiting more gently sloping lee sides and a reverse asymmetry in their flow‐parallel bed profile compared to mature dunes. Flow associated with the development of these “protodunes,” wherein the crest gradually shifts downstream towards its mature state, was investigated by studying the perturbation of the turbulent boundary layer over a succession of representative bedforms. Rigid, three‐dimensional models were studied in a refractive‐index‐matched experimental flume that enabled near‐surface quantification of mean velocities and Reynolds stresses using particle‐image velocimetry in wall‐normal and wall‐parallel measurement planes. Data indicate strong, topographically induced flow perturbations over the protodunes, to a similar relative degree to that found over mature dunes, despite their low‐angled slopes. The shape of the crest is found to be an important factor in the development of flow perturbations, and only in the case with the flattest crest was maximal speed‐up of flow, and reduction in turbulent stresses, found to occur upstream of the crest. Investigation of the log‐linearity of the boundary layer profile over the stoss sides showed that, although the profile is strongly perturbed, a log‐linear region exists, but is shifted vertically. A streamwise trend in friction velocity is thus present, showing a behavior similar to the trends in mean velocity. Analysis of the growth of the internal boundary layer on the dune stoss sides, beginning at the toe region, reveals a similar development for all dune shapes, despite clear differences in mean velocity and turbulent stress perturbations in their toe regions. The data presented herein provide the first documentation of flow over morphologies broadly characteristic of subtle, low‐angle, aeolian protodunes, and indicate key areas where further study is required to yield a more complete quantitative understanding of flow–form–transport couplings that govern their morphodynamics.

     
    more » « less
  5. ; ; ( , IMA Journal of Numerical Analysis)
    Abstract We introduce and analyze a partially augmented fully mixed formulation and a mixed finite element method for the coupled problem arising in the interaction between a free fluid and a poroelastic medium. The flows in the free fluid and poroelastic regions are governed by the Navier–Stokes and Biot equations, respectively, and the transmission conditions are given by mass conservation, balance of fluid force, conservation of momentum and the Beavers–Joseph–Saffman condition. We apply dual-mixed formulations in both domains, where the symmetry of the Navier–Stokes and poroelastic stress tensors is imposed in an ultra-weak and weak sense. In turn, since the transmission conditions are essential in the fully mixed formulation, they are imposed weakly by introducing the traces of the structure velocity and the poroelastic medium pressure on the interface as the associated Lagrange multipliers. Furthermore, since the fluid convective term requires the velocity to live in a smaller space than usual, we augment the variational formulation with suitable Galerkin-type terms. Existence and uniqueness of a solution are established for the continuous weak formulation, as well as a semidiscrete continuous-in-time formulation with nonmatching grids, together with the corresponding stability bounds and error analysis with rates of convergence. Several numerical experiments are presented to verify the theoretical results and illustrate the performance of the method for applications to arterial flow and flow through a filter. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email