 Award ID(s):
 1854434
 NSFPAR ID:
 10249390
 Date Published:
 Journal Name:
 SIAM journal on scientific computing
 Volume:
 43
 Issue:
 1
 ISSN:
 10648275
 Page Range / eLocation ID:
 A636A662
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

It is expensive to compute residual diffusivity in chaotic incompressible flows by solving advectiondiffusion equation due to the formation of sharp internal layers in the advection dominated regime. Proper orthogonal decomposition (POD) is a classical method to construct a small number of adaptive orthogonal basis vectors for low cost computation based on snapshots of fully resolved solutions at a particular molecular diffusivity D0* . The quality of POD basis deteriorates if it is applied to D0<< D0* . To improve POD, we adapt a superresolution generative adversarial deep neural network (SRGAN) to train a nonlinear mapping based on snapshot data at two values of D0* . The mapping models the sharpening effect on internal layers as D0 becomes smaller. We show through numerical experiments that after applying such a mapping to snapshots, the prediction accuracy of residual diffusivity improves considerably that of the standard POD.more » « less

This paper develops a treetopological local mesh refinement (TLMR) method on Cartesian grids for the simulation of bioinspired flow with multiple moving objects. The TLMR nests refinement mesh blocks of structured grids to the target regions and arrange the blocks in a tree topology. The method solves the timedependent incompressible flow using a fractionalstep method and discretizes the NavierStokes equation using a finitedifference formulation with an immersed boundary method to resolve the complex boundaries. When iteratively solving the discretized equations across the coarse and fine TLMR blocks, for better accuracy and faster convergence, the momentum equation is solved on all blocks simultaneously, while the Poisson equation is solved recursively from the coarsest block to the finest ones. When the refined blocks of the same block are connected, the parallel Schwarz method is used to iteratively solve both the momentum and Poisson equations. Convergence studies show that the algorithm is secondorder accurate in space for both velocity and pressure, and the developed mesh refinement technique is benchmarked and demonstrated by several canonical flow problems. The TLMR enables a fast solution to an incompressible flow problem with complex boundaries or multiple moving objects. Various bioinspired flows of multiple moving objects show that the solver can save over 80% computational time, proportional to the grid reduction when refinement is applied.more » « less

Gequations are popular level set Hamilton–Jacobi nonlinear partial differential equations (PDEs) of first or second order arising in turbulent combustion. Characterizing the effective burning velocity (also known as the turbulent burning velocity) is a fundamental problem there. We review relevant studies of the Gequation models with a focus on both the existence of effective burning velocity (homogenization), and its dependence on physical and geometric parameters (flow intensity and curvature effect) through representative examples. The corresponding physical background is also presented to provide motivations for mathematical problems of interest.
The
lack of coercivity of Hamiltonian is a hallmark of Gequations. When either the curvature of the level set or the strain effect of fluid flows is accounted for, the Hamiltonian becomeshighly nonconvex and nonlinear . In the absence of coercivity and convexity, the PDE (Eulerian) approach suffers from insufficient compactness to establish averaging (homogenization). We review and illustrate a suite of Lagrangian tools, most notably minmax (maxmin) game representations of curvature and strain Gequations, working in tandem with analysis of streamline structures of fluid flows and PDEs. We discuss open problems for future development in this emerging area of dynamic game analysis for averaging noncoercive, nonconvex, and nonlinear PDEs such as geometric (curvaturedependent) PDEs with advection. 
This paper employs an input–output based approach to analyse the convective velocities and transport of fluctuations in turbulent channel flows. The convective velocity for a fluctuating quantity associated with streamwise–spanwise wavelength pairs at each wallnormal location is obtained through the maximization of the power spectral density associated with the linearized Navier–Stokes equations with a turbulent mean profile and deltacorrelated Gaussian forcing. We first demonstrate that the mean convective velocities computed in this manner agree well with those reported previously in the literature. We then exploit the analytical framework to probe the underlying mechanisms contributing to the local convective velocity at different wallnormal locations by isolating the contributions of each streamwise–spanwise wavelength pair (flow scale). The resulting analysis suggests that the behaviour of the convective velocity in the nearwall region is influenced by largescale structures further away from the wall. These structures resemble Townsend’s attached eddies in the crossplane, yet show incomplete similarity in the streamwise direction. We then investigate the role of each linear term in the momentum equation to isolate the contribution of the pressure, mean shear, and viscous effects to the deviation of the convective velocity from the mean at each flow scale. Our analysis highlights the role of the viscous effects, particularly in regards to large channel spanning structures whose influence extends to the nearwall region. The results of this work suggest the promise of an input–output approach for analysing convective velocity across a range of flow scales using only the mean velocity profile.more » « less

Inspired by the recent realization of a twodimensional (2D) chiral fluid as an active monolayer droplet moving atop a 3D Stokesian fluid, we formulate mathematically its freeboundary dynamics. The surface droplet is described as a general 2D linear, incompressible and isotropic fluid, having a viscous shear stress, an active chiral driving stress and a Hall stress allowed by the lack of timereversal symmetry. The droplet interacts with itself through its driven internal mechanics and by driving flows in the underlying 3D Stokes phase. We pose the dynamics as the solution to a singular integral–differential equation, over the droplet surface, using the mapping from surface stress to surface velocity for the 3D Stokes equations. Specializing to the case of axisymmetric droplets, exact representations for the chiral surface flow are given in terms of solutions to a singular integral equation, solved using both analytical and numerical techniques. For a discshaped monolayer, we additionally employ a semianalytical solution that hinges on an orthogonal basis of Bessel functions and allows for efficient computation of the monolayer velocity field, which ranges from a nearly solidbody rotation to a unidirectional edge current, depending on the subphase depth and the Saffman–Delbrück length. Except in the nearwall limit, these solutions have divergent surface shear stresses at droplet boundaries, a signature of systems with codimensionone domains embedded in a 3D medium. We further investigate the effect of a Hall viscosity, which couples radial and transverse surface velocity components, on the dynamics of a closing cavity. Hall stresses are seen to drive inward radial motion, even in the absence of edge tension.more » « less