Computation of head–disk interface gap micro flowfields using DSMC and continuum–atomistic hybrid methods

This paper discusses computational modeling of micro flow in the head–disk interface (HDI) gap using the direct simulation Monte Carlo (DSMC) method. Modeling considerations are discussed in detail both for a stand‐alone DSMC computation and for the case of a hybrid continuum–atomistic simulation that couples the Navier–Stokes (NS) equation to a DSMC solver. The impact of the number of particles and number of cells on the accuracy of a DSMC simulation of the HDI gap is investigated both for two‐ and three‐dimensional configurations. An appropriate implicit boundary treatment method for modeling inflow and outflow boundaries is used in this work for a three‐dimensional DSMC micro flow simulation. As the flow outside the slider is in the continuum regime, a hybrid continuum–atomistic method based on the Schwarz alternating method is used to couple the DSMC model in the slider bearing region to the flow outside the slider modeled by NS equation. Schwarz coupling is done in two dimensions by taking overlap regions along two directions and the Chapman–Enskog distribution is employed for imposing the boundary condition from the continuum region to the DSMC region. Converged hybrid flow solutions are obtained in about five iterations and the hybrid DSMC–NS solutions show good agreement with the exact solutions in the entire domain considered. An investigation on the impact of the size of the overlap region on the convergence behavior of the Schwarz method indicates that the hybrid coupling by the Schwarz method is weakly dependent on the size of the overlap region. However, the use of a finite overlap region will facilitate the exchange of boundary conditions as the hybrid solution has been found to diverge in the absence of an overlap region for coupling the two models. Copyright © 2009 John Wiley & Sons, Ltd.     

A novel efficient hybrid DSMC–dynamic collision limiter algorithm for multiscale transitional flows

A new 2D parallel multispecies polyatomic particle–based hybrid flow solver is developed by coupling the Direct Simulation Monte Carlo (DSMC) method with a novel Dynamic Collision Limiter (DCL) approach to solve multiscale transitional flows. The hybrid DSMC‐DCL solver can solve nonequilibrium multiscale flows with length scales ranging from continuum to rarefied. The DCL method, developed in this work, dynamically assigns different number of collisions in cells, which is based on the local value of K‐S parameter such that the number of collisions per time step is limited in near‐equilibrium flow regions. Present hybrid solver uses the Kolmogorov‐Smirnov statistical test as the continuum breakdown parameter, based on which, the solution domain is decomposed into near‐equilibrium and nonequilibrium flow regions. Direct Simulation Monte Carlo is used where nonequilibrium flow regions are encountered, while the DCL method is used where flow regions are found to be in near‐equilibrium state. In this work, we have studied hypersonic flow of nitrogen over a blunt body with an aerospike and supersonic flow of argon through a micronozzle. The results obtained by the hybrid DSMC‐DCL solver are compared and shown to agree well with the experimental data and with those obtained from DSMC, with significant savings in the computational cost.     

Empirical slip and viscosity model performance for microscale gas flow

For the simple geometries of Couette and Poiseuille flows, the velocity profile maintains a similar shape from continuum to free molecular flow. Therefore, modifications to the fluid viscosity and slip boundary conditions can improve the continuum based Navier–Stokes solution in the non‐continuum non‐equilibrium regime. In this investigation, the optimal modifications are found by a linear least‐squares fit of the Navier–Stokes solution to the non‐equilibrium solution obtained using the direct simulation Monte Carlo (DSMC) method. Models are then constructed for the Knudsen number dependence of the viscosity correction and the slip model from a database of DSMC solutions for Couette and Poiseuille flows of argon and nitrogen gas, with Knudsen numbers ranging from 0.01 to 10. Finally, the accuracy of the models is measured for non‐equilibrium cases both in and outside the DSMC database. Flows outside the database include: combined Couette and Poiseuille flow, partial wall accommodation, helium gas, and non‐zero convective acceleration. The models reproduce the velocity profiles in the DSMC database within an L₂ error norm of 3% for Couette flows and 7% for Poiseuille flows. However, the errors in the model predictions outside the database are up to five times larger. Copyright © 2005 John Wiley & Sons, Ltd.     

DSMC／EPSM混合算法研究

将EPSM算法与DSMC方法结合，构造了可模拟含近连续流区及过渡流区的DSMC/EPSM混合算法。运用混合算法模拟了马赫数等于5时超音速竖板绕流及马赫数等于4时超音速平板绕流，并将结果与DSMC算法的结果进行比较，证明了DSMC／EPSM混合算法的有效性，同时将EPSM算法与DSMC算法的效率进行了比较。     

Investigation on a coupled Navier–Stokes–Direct Simulation Monte Carlo method for the simulation of plume flowfield of a conical nozzle

To assess the plume effects of space thrusters, the accurate plume flowfield is indispensable. The plume flow of thrusters involves both continuum and rarefied flow regimes. Coupled Navier–Stokes–Direct Simulation Monte Carlo (NS–DSMC) method is a major approach to the simulation of continuum‐rarefied flows. An axisymmetric coupled NS–DSMC solver, possessing adaptive‐interface and two‐way coupling features, is investigated in this paper for the simulation of the nozzle and plume flows of thrusters. The state‐based coupling scheme is adopted, and the gradient local Knudsen number is used to indicate the breakdown of continuum solver. The nitrogen flows in a conical nozzle and its plume are chosen as the reference case to test the coupled solver. The threshold value of the continuum breakdown parameter is studied based on both theoretical kinetic velocity sampling and coupled numerical tests. Succeeding comparisons between coupled and full DSMC results demonstrate their conformities, meanwhile, the former saves 58.8% computational time. The pitot pressure evaluated from the coupled simulation result is compared with the experimental data proposed in literatures, revealing that the coupled method makes precise predictions on the experimental pitot pressure. Copyright © 2014 John Wiley & Sons, Ltd.     

Multi‐scale time integration for transient conjugate heat transfer

The quasi‐steady assumption is commonly adopted in existing transient fluid–solid‐coupled convection–conduction (conjugate) heat transfer simulations, which may cause non‐negligible errors in certain cases of practical interest. In the present work, we adopt a new multi‐scale framework for the fluid domain formulated in a triple‐timing form. The slow‐varying temporal gradient corresponding to the time scales in the solid domain has been effectively included in the fluid equations as a source term, whilst short‐scale unsteadiness of the fluid domain is captured by a local time integration at a given ‘frozen’ large scale time instant. For concept proof, validation and demonstration purposes, the proposed methodology has been implemented in a loosely coupled procedure in conjunction with a hybrid interfacing treatment for coupling efficiency and accuracy. The present results indicate that a much enhanced applicability can be achieved with relatively small modifications of existing transient conjugate heat transfer methods at little extra cost. Copyright © 2016 John Wiley & Sons, Ltd.     

A quasi‐dynamic procedure for coupled thermal simulations

This paper outlines the development and adaptation of a coupling strategy for transient temperature analysis in a solid via a conjugate heat transfer method. This study proposes a quasi‐dynamic coupling procedure to bridge the temporal disparities between the fluid and the solid. In this approach, dynamic thermal modeling in the solid is coupled with a sequence of steady states in the fluid. This quasi‐dynamic algorithm has been applied to the problem of convective heat transfer over, and transient conduction heat transfer within, a flat plate using the severe thermal conditions of a solid propellant rocket. Two different coupled thermal computations have been performed. In the first one—referred to as the reference computation—the coupling period is equal to the smallest solid time constant. In the second one, a very large coupling period is used. The results show that the procedure can predict accurate transient temperature fields at a reasonable computational cost. The simulation CPU time is approximately reduced by up to 90%, while maintaining a very good accuracy. All the details of the numerical test case are given in the paper. This application illustrates the capabilities and the overall efficiency of this coupled approach in a solid transient problem using long term simulations of time dependent flows. Copyright © 2013 John Wiley & Sons, Ltd.     

A gas‐kinetic scheme for the simulation of turbulent flows on unstructured meshes

This paper presents a numerical method for simulating turbulent flows via coupling the Boltzmann BGK equation with Spalart–Allmaras one equation turbulence model. Both the Boltzmann BGK equation and the turbulence model equation are carried out using the finite volume method on unstructured meshes, which is different from previous works on structured grid. The application of the gas‐kinetic scheme is extended to the simulation of turbulent flows with arbitrary geometries. The adaptive mesh refinement technique is also adopted to reduce the computational cost and improve the efficiency of meshes. To organize the unstructured mesh data structure efficiently, a non‐manifold hybrid mesh data structure is extended for polygonal cells. Numerical experiments are performed on incompressible flow over a smooth flat plate and compressible turbulent flows around a NACA 0012 airfoil using unstructured hybrid meshes. These numerical results are found to be in good agreement with experimental data and/or other numerical solutions, demonstrating the applicability of the proposed method to simulate both subsonic and transonic turbulent flows. Copyright © 2016 John Wiley & Sons, Ltd.     

Multi‐objective optimization of a rotating cooling channel with staggered pin‐fins for heat transfer augmentation

A rotating channel with staggered pin‐fins is formulated numerically and optimized for performance (heat transfer/required pumping power) using a Kriging meta‐model and hybrid multi‐objective evolutionary algorithm. Two design variables related to cooling channel height, pin diameter, and spacing between the pins are selected for optimization, and two‐objective functions related to the heat transfer and friction loss are employed. A design of experiment is performed, and 20 designs are generated by Latin hypercube sampling. The objective function values are evaluated using a Reynolds‐averaged Navier–Stokes solver, and a Kriging model is constructed to obtain a Pareto‐optimal front through a multi‐objective evolutionary algorithm. Rotation in a cooling channel with staggered pin‐fins induces Coriolis force that causes a heat transfer discrepancy between the trailing (pressure) and leading (suction) surfaces, with a higher Nusselt number on the trailing surface. The tradeoff between the two competing objective functions is determined, and the distribution of the Pareto‐optimal solutions in the design space is discussed through k‐means clustering. In the optimal designs, with a decrease in spacing between the pins, heat transfer is enhanced whereas friction loss is increased. Copyright © 2011 John Wiley & Sons, Ltd.     

An exact non‐linear Navier–Stokes compressible‐flow solution for CFD code verification

An exact similarity solution of the compressible‐flow Navier–Stokes equations is presented, which embeds supersonic, transonic, and subsonic regions. Describing the viscous and heat‐conducting high‐gradient flow in a shock wave, the solution accommodates non‐linear temperature‐dependent viscosity as well as heat‐conduction coefficients and provides the variation of all the flow variables and their derivatives. Also presented are methods to obtain time‐dependent and/or multi‐dimensional solutions as well as verification benchmarks of increasing severity. Comparisons between the developed analytical solution and CFD solutions of the Navier–Stokes equations, with determination of convergence rates and orders of accuracy of these solutions, illustrate the utility of the developed exact solution for verification purposes. Copyright © 2012 John Wiley & Sons, Ltd.     

Local maximum principle satisfying high‐order non‐oscillatory schemes

The main contribution of this work is to classify the solution region including data extrema for which high‐order non‐oscillatory approximation can be achieved. It is performed in the framework of local maximum principle (LMP) and non‐conservative formulation. The representative uniformly second‐order accurate schemes are converted in to their non‐conservative form using the ratio of consecutive gradients. Using the local maximum principle, these non‐conservative schemes are analyzed for their non‐linear LMP/total variation diminishing stability bounds which classify the solution region where high‐order accuracy can be achieved. Based on the bounds, second‐order accurate hybrid numerical schemes are constructed using a shock detector. The presented numerical results show that these hybrid schemes preserve high accuracy at non‐sonic extrema without exhibiting any induced local oscillations or clipping error. Copyright © 2015 John Wiley & Sons, Ltd.     

Investigation of the sensitivity of turbulent closures and coupling of hybrid RANS‐LES models for predicting flow fields with separation and reattachment

Hybrid models have found widespread applications for simulation of wall‐bounded flows at high Reynolds numbers. Typically, these models employ Reynolds‐averaged Navier–Stokes (RANS) and large eddy simulation (LES) in the near‐body and off‐body regions, respectively. A number of coupling strategies between the RANS and LES regions have been proposed, tested, and applied in the literature with varying degree of success. Linear eddy‐viscosity models (LEVM) are often used for the closure of turbulent stress tensor in RANS and LES regions. LEVM incorrectly predicts the anisotropy of Reynolds normal stress at the RANS‐LES interface region. To overcome this issue, use of non‐linear eddy‐viscosity models (NLEVM) have started receiving attention. In this study, a generic non‐linear blended modeling framework for performing hybrid simulations is proposed. Flow over the periodic hills is used as the test case for model evaluation. This case is chosen due to complex flow physics with simplified geometry. Analysis of the simulations suggests that the non‐linear hybrid models show a better performance than linear hybrid models. It is also observed that the non‐linear closures are less sensitive to the RANS‐LES coupling and grid resolution. Copyright © 2016 John Wiley & Sons, Ltd.     

气体化学反应流动的DSMC／EPSM混合算法研究

发展了平衡粒子模拟方法(EPSM)，建立了与高温气体化学反应动力学理论相匹配的：EPSM耦合模型，并通过混合参数进行流区的自动识别，将：EPSM方法与蒙特卡罗直接模拟方法(OSMC)结合，构造了可模拟化学反应流动的DSMC／EPSM混合算法。应用该算法对汲及化学反应的二维高超音速竖板绕流流场进行模拟，将结果与DSMC方法的结果进行比较，验证了新算法对求解化学反应流动的可行性。将混合算法的计算效率与DSMC方法的计算效率进行比较，发现混合算法能够大大提高计算效率。     

A spectral‐element discontinuous Galerkin thermal lattice Boltzmann method for conjugate heat transfer applications

We present a spectral‐element discontinuous Galerkin thermal lattice Boltzmann method for fluid–solid conjugate heat transfer applications. Using the discrete Boltzmann equation, we propose a numerical scheme for conjugate heat transfer applications on unstructured, non‐uniform grids. We employ a double‐distribution thermal lattice Boltzmann model to resolve flows with variable Prandtl (Pr) number. Based upon its finite element heritage, the spectral‐element discontinuous Galerkin discretization provides an effective means to model and investigate thermal transport in applications with complex geometries. Our solutions are represented by the tensor product basis of the one‐dimensional Legendre–Lagrange interpolation polynomials. A high‐order discretization is employed on body‐conforming hexahedral elements with Gauss–Lobatto–Legendre quadrature nodes. Thermal and hydrodynamic bounce‐back boundary conditions are imposed via the numerical flux formulation that arises because of the discontinuous Galerkin approach. As a result, our scheme does not require tedious extrapolation at the boundaries, which may cause loss of mass conservation. We compare solutions of the proposed scheme with an analytical solution for a solid–solid conjugate heat transfer problem in a 2D annulus and illustrate the capture of temperature continuities across interfaces for conductivity ratio γ > 1. We also investigate the effect of Reynolds (Re) and Grashof (Gr) number on the conjugate heat transfer between a heat‐generating solid and a surrounding fluid. Steady‐state results are presented for Re = 5?40 and Gr = 10⁵?10⁶. In each case, we discuss the effect of Re and Gr on the heat flux (i.e. Nusselt number Nu) at the fluid–solid interface. Our results are validated against previous studies that employ finite‐difference and continuous spectral‐element methods to solve the Navier–Stokes equations. Copyright © 2016 John Wiley & Sons, Ltd.     

Enhanced multi‐objective optimization of a dimpled channel through evolutionary algorithms and multiple surrogate methods

Using a multi‐objective evolutionary algorithm (MOEA) and enhanced surrogate approximations, the present study demonstrates the numerical analysis and optimization of staggered‐dimple channels. Two surrogates, the response surface approximation (RSA) model and the Kriging (KRG) model, are applied in light of the surrogate fidelity of the approximate analysis. An enhanced Pareto‐optimal front is obtained by performing local resampling of the Pareto‐optimal front, which provides relatively more accurate Pareto‐optimal solutions in the design space for each surrogate model. Three dimensionless design variables are selected, which are related to geometric parameters, namely, the channel height, dimple print diameter, dimple spacing, and dimple depth. Two objective functions are selected that are related to the heat transfer and pressure loss, respectively. The objective‐function values are numerically evaluated through Reynolds‐averaged Navier–Stokes analysis at the design points that are selected through the Latin hypercube sampling method. Using these numerical simulations two surrogates, viz, the RSA and Kriging models, are constructed for each objective function and a hybrid MOEA is applied to obtain the Pareto‐optimal front. For the particular implementation of surrogate models, it is observed that Pareto‐optimal predictions of the RSA model are better than those of the KRG model, whereas the KRG model predicts equally well at the off‐Pareto‐region (region away from the Pareto‐optimal solutions), which is not the case with the RSA model. The local resampling of the Pareto‐optimal front increases the fidelity of the approximate solutions near the Pareto‐optimal region. The ratios of the channel height to the dimple print diameter and of the dimple print diameter to the dimple pitch are found to be more sensitive along the Pareto‐optimal front than the ratio of the dimple depth to the print diameter. The decrease of the ratio of the channel height to the dimple diameter and the increase of the ratio of the dimple print diameter to the pitch lead to greater heat transfer at the expense of the pressure loss, whereas the ratio of the dimple depth to the print diameter is rather insensitive to Pareto‐optimal solutions. Pareto‐optimal solutions at higher values of the Nusselt number are associated with higher values of the pressure loss due to the increased recirculation, mixing of fluid and vorticity generation. Copyright © 2010 John Wiley & Sons, Ltd.     

A simple high‐resolution advection scheme

A simple, robust, mass‐conserving numerical scheme for solving the linear advection equation is described. The scheme can estimate peak solution values accurately even in regions where spatial gradients are high. Such situations present a severe challenge to classical numerical algorithms. Attention is restricted to the case of pure advection in one and two dimensions since this is where past numerical problems have arisen. The authors' scheme is of the Godunov type and is second‐order in space and time. The required cell interface fluxes are obtained by MUSCL interpolation and the exact solution of a degenerate Riemann problem. Second‐order accuracy in time is achieved via a Runge–Kutta predictor–corrector sequence. The scheme is explicit and expressed in finite volume form for ease of implementation on a boundary‐conforming grid. Benchmark test problems in one and two dimensions are used to illustrate the high‐spatial accuracy of the method and its applicability to non‐uniform grids. Copyright © 2007 John Wiley & Sons, Ltd.     

An implicit two‐dimensional non‐hydrostatic model for free‐surface flows

An implicit finite volume model in sigma coordinate system is developed to simulate two‐dimensional (2D) vertical free surface flows, deploying a non‐hydrostatic pressure distribution. The algorithm is based on a projection method which solves the complete 2D Navier–Stokes equations in two steps. First the pressure term in the momentum equations is excluded and the resultant advection–diffusion equations are solved. In the second step the continuity and the momentum equation with only the pressure terms are solved to give a block tri‐diagonal system of equation with pressure as the unknown. This system can be solved by a direct matrix solver without iteration. A new implicit treatment of non‐hydrostatic pressure, similar to the lower layers is applied to the top layer which makes the model free of any hydrostatic pressure assumption all through the water column. This treatment enables the model to evaluate both free surface elevation and wave celerity more accurately. A series of numerical tests including free‐surface flows with significant vertical accelerations and nonlinear behaviour in shoaling zone are performed. Comparison between numerical results, analytical solutions and experimental data demonstrates a satisfactory performance. Copyright © 2006 John Wiley & Sons, Ltd.     

A two‐dimensional vertical non‐hydrostatic σ model with an implicit method for free‐surface flows

An implicit finite difference model in the σ co‐ordinate system is developed for non‐hydrostatic, two‐dimensional vertical plane free‐surface flows. To accurately simulate interaction of free‐surface flows with uneven bottoms, the unsteady Navier–Stokes equations and the free‐surface boundary condition are solved simultaneously in a regular transformed σ domain using a fully implicit method in two steps. First, the vertical velocity and pressure are expressed as functions of horizontal velocity. Second, substituting these relationship into the horizontal momentum equation provides a block tri‐diagonal matrix system with the unknown of horizontal velocity, which can be solved by a direct matrix solver without iteration. A new treatment of non‐hydrostatic pressure condition at the top‐layer cell is developed and found to be important for resolving the phase of wave propagation. Additional terms introduced by the σ co‐ordinate transformation are discretized appropriately in order to obtain accurate and stable numerical results. The developed model has been validated by several tests involving free‐surface flows with strong vertical accelerations and non‐linear waves interacting with uneven bottoms. Comparisons among numerical results, analytical solutions and experimental data show the capability of the model to simulate free‐surface flow problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Linear and non‐linear simulations of feedback control in plane Poiseuille flow

This paper examines the performance of optimal linear quadratic state and output feedback controllers in stabilizing two‐dimensional perturbations in a plane Poiseuille flow. The synthesis of the controllers is based on a linearized model of the flow using a new set of interpolating polynomials in the wall‐normal direction, which automatically satisfy the homogeneous Dirichlet and Neumann boundary conditions at the walls and eliminate spurious eigenvalues. The controllers are implemented into a non‐linear Navier–Stokes solver, which is modified to compute the evolution of the flow perturbations. Two cases are examined, one with small initial disturbances that do not violate the linearity assumptions and the other with much larger disturbances that trigger the non‐linear convection terms. For the smallest disturbances, the solver accurately reproduced the results of the linear simulations of open‐ and closed‐loop systems. The simulations for the larger disturbances without control showed a rapid initial growth but the flow soon reached a saturated state in agreement with previous findings in the literature. The large initial growth is a consequence of the non‐normal nature of the system dynamics. The state feedback and output feedback controllers were able to reduce significantly the perturbation energy. For the larger disturbances, the energy calculated from the state variables is well below the energy evaluated by direct integration of the velocity field. This is probably due to the non‐linear terms transferring energy to harmonics of the considered wavenumber, which are not sensed by the linear controller. Copyright © 2008 John Wiley & Sons, Ltd.     

Fully non‐linear free‐surface simulations by a 3D viscous numerical wave tank

A finite difference scheme using a modified marker‐and‐cell (MAC) method is applied to investigate the characteristics of non‐linear wave motions and their interactions with a stationary three‐dimensional body inside a numerical wave tank (NWT). The Navier–Stokes (NS) equation is solved for two fluid layers, and the boundary values are updated at each time step by a finite difference time marching scheme in the frame of a rectangular co‐ordinate system. The viscous stresses and surface tension are neglected in the dynamic free‐surface condition, and the fully non‐linear kinematic free‐surface condition is satisfied by the density function method developed for two fluid layers. The incident waves are generated from the inflow boundary by prescribing a velocity profile resembling flexible flap wavemaker motions, and the outgoing waves are numerically dissipated inside an artificial damping zone located at the end of the tank. The present NS–MAC NWT simulations for a vertical truncated circular cylinder inside a rectangular wave tank are compared with the experimental results of Mercier and Niedzwecki, an independently developed potential‐based fully non‐linear NWT, and the second‐order diffraction computation. Copyright © 1999 John Wiley & Sons, Ltd.