skip to main content

Title: Laminar and Turbulent Behavior Captured by A 3-D Kinetic-Based Discrete Dynamic System
We have derived a 3-D kinetic-based discrete dynamic system (DDS) from the lattice Boltzmann equation (LBE) for incompressible flows through a Galerkin procedure. Expressed by a poor-man lattice Boltzmann equation (PMLBE), it involves five bifurcation parameters including relaxation time from the LBE, splitting factor of large and sub-grid motion scales, and wavevector components from the Fourier space. Numerical experiments have shown that the DDS can capture laminar behaviors of periodic, subharmonic, n-period, and quasi-periodic and turbulent behaviors of noisy periodic with harmonic, noisy subharmonic, noisy quasi-periodic, and broadband power spectra. In this work, we investigated the effects of bifurcation parameters on the capturing of the laminar and turbulent flows in terms of the convergence of time series and the pattern of power spectra. We have found that the 2nd order and 3rd order PMLBEs are both able to capture laminar and turbulent flow behaviors but the 2nd order DDS performs better with lower computation cost and more flow behaviors captured. With the specified ranges of the bifurcation parameters, we have identified two optimal bifurcation parameter sets for laminar and turbulent behaviors. Beyond this work, we are exploring the regime maps for a deeper understanding of the contributions of the bifurcation parameters to the capturing of laminar and turbulent behaviors. Surrogate models (to replace the PMLBE) are being developed using deep learning techniques to overcome the overwhelming computation cost for the regime maps. Meanwhile, the DDS is being employed in the large eddy simulation of turbulent pulsatile flows to provide dynamic sub-grid scale information.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ;
Brehm, Christoph; Pandya, Shishir
Date Published:
Journal Name:
Eleventh International Conference on Computational Fluid Dynamics (ICCFD11)
Page Range / eLocation ID:
1 - 9
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Image-based computational fluid dynamics (CFD) has become a new capability for determining wall stresses of pulsatile flows. However, a computational platform that directly connects image information to pulsatile wall stresses is lacking. Prevailing methods rely on manual crafting of a hodgepodge of multidisciplinary software packages, which is usually laborious and error-prone. We present a new computational platform, to compute wall stresses in image-based pulsatile flows using the volumetric lattice Boltzmann method (VLBM). The novelty includes: (1) a unique image processing to extract flow domain and local wall normality, (2) a seamless connection between image extraction and VLBM, (3) an en-route calculation of strain-rate tensor, and (4) GPU acceleration (not included here). We first generalize the streaming operation in the VLBM and then conduct application studies to demonstrate its reliability and applicability. A benchmark study is for laminar and turbulent pulsatile flows in an image-based pipe (Reynolds number: 10 to 5000). The computed pulsatile velocity and shear stress are in good agreements with Womersley's analytical solutions for laminar pulsatile flows and concurrent laboratory measurements for turbulent pulsatile flows. An application study is to quantify the pulsatile hemodynamics in image-based human vertebral and carotid arteries including velocity vector, pressure, and wall-shear stress. The computed velocity vector fields are in reasonably well agreement with MRA (magnetic resonance angiography) measured ones. This computational platform is good for image-based CFD with medical applications and pore-scale porous media flows in various natural and engineering systems. 
    more » « less
  2. Abstract In this article, the recently discovered phenomenon of delayed Hopf bifurcations (DHB) in reaction–diffusion partial differential equations (PDEs) is analysed in the cubic Complex Ginzburg–Landau equation, as an equation in its own right, with a slowly varying parameter. We begin by using the classical asymptotic methods of stationary phase and steepest descents on the linearized PDE to show that solutions, which have approached the attracting quasi-steady state (QSS) before the Hopf bifurcation remain near that state for long times after the instantaneous Hopf bifurcation and the QSS has become repelling. In the complex time plane, the phase function of the linearized PDE has a saddle point, and the Stokes and anti-Stokes lines are central to the asymptotics. The non-linear terms are treated by applying an iterative method to the mild form of the PDE given by perturbations about the linear particular solution. This tracks the closeness of solutions near the attracting and repelling QSS in the full, non-linear PDE. Next, we show that beyond a key Stokes line through the saddle there is a curve in the space-time plane along which the particular solution of the linear PDE ceases to be exponentially small, causing the solution of the non-linear PDE to diverge from the repelling QSS and exhibit large-amplitude oscillations. This curve is called the space–time buffer curve. The homogeneous solution also stops being exponentially small in a spatially dependent manner, as determined also by the initial data and time. Hence, a competition arises between these two solutions, as to which one ceases to be exponentially small first, and this competition governs spatial dependence of the DHB. We find four different cases of DHB, depending on the outcomes of the competition, and we quantify to leading order how these depend on the main system parameters, including the Hopf frequency, initial time, initial data, source terms, and diffusivity. Examples are presented for each case, with source terms that are a uni-modal function, a smooth step function, a spatially periodic function and an algebraically growing function. Also, rich spatio-temporal dynamics are observed in the post-DHB oscillations. Finally, it is shown that large-amplitude source terms can be designed so that solutions spend substantially longer times near the repelling QSS, and hence, region-specific control over the delayed onset of oscillations can be achieved. 
    more » « less
  3. Adibi, Ali ; Lin, Shawn-Yu ; Scherer, Axel (Ed.)
    Periodic optical lattices consisting of isolated-particle arrays in vacuum are treated with rigorous electromagnetics. These structures possess a wealth of interesting properties including perfect reflection across small or large spectral bandwidths depending on the choice of materials and design parameters. Pertinent spectral expressions have been observed theoretically and experimentally via one-dimensional (1D) and two-dimensional (2D) structures commonly known as resonant gratings, metamaterials, and metasurfaces. The physical cause of perfect reflection and related properties is guided-mode resonance mediated by lateral Bloch modes excited by evanescent diffraction orders in the subwavelength regime. Here, we review recent results on differentiation of local Mie resonance and guided-mode lattice resonance in causing resonant reflection by periodic particle assemblies. We treat a classic 2D periodic array consisting of dielectric spheres. To disable Mie resonance, we apply antireflection (AR) coatings to the spheres. Reflectance maps for coated and uncoated spheres demonstrate that perfect reflection persists in both cases. We find that the Mie scattering efficiency of an AR-coated sphere is greatly diminished. Additionally, in a 1D cylindrical rod-type lattice, we investigate and compare local field profiles in periodic assemblies and in the constituent isolated particles. In general, the lattice and particle resonance wavelengths differ. When the lateral leaky-mode field profiles approach the isolated-particle Mie field profiles, the resonance locus tends towards the Mie resonance wavelength. This correspondence is referred to as Mie modal memory. These fundamentals may help distinguish Mie effects and leaky-mode lattice effects in generating the observed spectra in this class of optical devices while elucidating the basic resonance properties across the entire spectral domain. 
    more » « less
  4. Direct numerical simulation (DNS) of turbulent flows is computationally expensive and cannot be applied to flows with large Reynolds numbers. Low-resolution large eddy simulation (LES) is a popular alternative, but it is unable to capture all of the scales of turbulent transport accurately. Reconstructing DNS from low-resolution LES is critical for large-scale simulation in many scientific and engineering disciplines, but it poses many challenges to existing super-resolution methods due to the complexity of turbulent flows and computational cost of generating frequent LES data. We propose a physics-guided neural network for reconstructing frequent DNS from sparse LES data by enhancing its spatial resolution and temporal frequency. Our proposed method consists of a partial differential equation (PDE)-based recurrent unit for capturing underlying temporal processes and a physics-guided super-resolution model that incorporates additional physical constraints. We demonstrate the effectiveness of both components in reconstructing the Taylor-Green Vortex using sparse LES data. Moreover, we show that the proposed recurrent unit can preserve the physical characteristics of turbulent flows by leveraging the physical relationships in the Navier-Stokes equation. 
    more » « less
  5. We present numerical simulations of laminar and turbulent channel flow of an elastoviscoplastic fluid. The non-Newtonian flow is simulated by solving the full incompressible Navier–Stokes equations coupled with the evolution equation for the elastoviscoplastic stress tensor. The laminar simulations are carried out for a wide range of Reynolds numbers, Bingham numbers and ratios of the fluid and total viscosity, while the turbulent flow simulations are performed at a fixed bulk Reynolds number equal to 2800 and weak elasticity. We show that in the laminar flow regime the friction factor increases monotonically with the Bingham number (yield stress) and decreases with the viscosity ratio, while in the turbulent regime the friction factor is almost independent of the viscosity ratio and decreases with the Bingham number, until the flow eventually returns to a fully laminar condition for large enough yield stresses. Three main regimes are found in the turbulent case, depending on the Bingham number: for low values, the friction Reynolds number and the turbulent flow statistics only slightly differ from those of a Newtonian fluid; for intermediate values of the Bingham number, the fluctuations increase and the inertial equilibrium range is lost. Finally, for higher values the flow completely laminarizes. These different behaviours are associated with a progressive increases of the volume where the fluid is not yielded, growing from the centreline towards the walls as the Bingham number increases. The unyielded region interacts with the near-wall structures, forming preferentially above the high-speed streaks. In particular, the near-wall streaks and the associated quasi-streamwise vortices are strongly enhanced in an highly elastoviscoplastic fluid and the flow becomes more correlated in the streamwise direction. 
    more » « less