The discrete Green's function paradigm for two-way coupled Euler–Lagrange simulation

We outline a methodology for the simulation of two-way coupled particle-laden flows. The drag force that couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. In addition to hydrodynamics, the method can be extended to other physics present in particle-laden flows such as heat transfer and electromagnetism. The discrete Green's functions for the Navier–Stokes equations are obtained at low particle Reynolds number in a two-plane channel geometry. We perform verification at different Reynolds numbers for a particle settling under gravity parallel to a plane wall, for different wall-normal separations. Compared with other point-particle schemes, the Stokesian discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations. To account for degradation in accuracy away from the wall at finite Reynolds number, we extend the present methodology to an Oseen-like discrete Green's function. The extended discrete Green's function method is found to be accurate within $6\,\%$ at all wall-normal separations for particle Reynolds numbers up to 24. The discrete Green's function approach is well suited to dilute systems with significant mass loading and this is highlighted by comparison against other Euler–Lagrange as well as particle-resolved simulations of gas–solid turbulent channel flow. Strong particle–turbulence coupling is observed in the form of turbulence modification and turbophoresis suppression, and these observations are placed in context of the different methods.

Conjugate Heat Transfer Analysis Using the Discrete Green's Function

Conjugate heat transfer problems generally require a coupled solution of the temperature fields in the fluid and solid domains. Implementing the boundary condition at the surface of the solid using a discrete Green's function (DGF) decouples the solutions. A DGF is determined first considering only the fluid domain with prescribed thermal boundary conditions, then the temperature distribution in the solid is calculated using standard numerical methods. The only compatibility requirement is that the DGF must be specified with the same discretization as the surface of the solid. The method is demonstrated for both steady-state and transient heating of a thin plate with laminar boundary layers flowing over both sides. The resulting set of linear algebraic equations for the steady-state problem or linear ordinary differential equations for the transient problem are easily solved using conventional scientific programming packages. The method converges with nearly second-order accuracy as the discretization resolution is increased.

The Discrete Green's Function for Convective Heat Transfer—Part 2: Semi-Analytical Estimates of Boundary Layer Discrete Green's Function

This is the second paper in a set that defines the discrete Green's function (DGF). This paper focuses first on the turbulent boundary layer and presents two different methods to estimate the DGF. The long-element formulation defines the DGF with just two simple algebraic equations, but it is not quantitatively accurate for short element lengths. A short element correction is derived, but must be recalculated for each selection of flow parameters and element lengths. A similarity solution is derived that allows accurate estimates of the DGF diagonal elements for laminar boundary layers and for turbulent boundary layers discretized with short element lengths. To illustrate other methods to derive DGFs in more complex flows, a low-resolution DGF for laminar stagnation line boundary layers is determined using the skin-friction formulation combined with similarity solutions for two different thermal boundary conditions. Stagnation line flow is shown to be highly sensitive to the thermal boundary condition, and this can be analyzed effectively using the DGF.

The Discrete Green's Function for Convective Heat Transfer—Part 1: Definition and Physical Understanding

The discrete Green's function (DGF) is a superposition-based descriptor of the relationship between the surface temperature and the convective heat transfer from a surface. The surface is discretized into a finite number of elements and the DGF matrix elements relate the heat transfer out of any element i to the temperature rise on every element j of the surface. For a given flow situation, the DGF is insensitive to the thermal boundary condition so it allows direct calculation of the heat transfer for any temperature distribution and noniterative solution of conjugate heat transfer problems. The diagonal elements of the matrix are determined solely by the local velocity field while the off-diagonals are determined by the spread of the thermal wake downstream of a heated element. An analytical DGF for the laminar flat-plate boundary layers is included as an example.

Wave scattering on lattice structures involving an array of cracks

Scattering of waves as a result of a vertical array of equally spaced cracks on a square lattice is studied. The convenience of Floquet periodicity reduces the study to that of scattering of a specific wave-mode from a single crack in a waveguide. The discrete Green’s function, for the waveguide, is used to obtain the semi-analytical solution for the scattering problem in the case of finite cracks whereas the limiting case of semi-infinite cracks is tackled by an application of the Wiener–Hopf technique. Reflectance and transmittance of such an array of cracks, in terms of incident wave parameters, is analysed. Potential applications include construction of tunable atomic-scale interfaces to control energy transmission at different frequencies.