Recycling inflow method for simulations of spatially evolving turbulent boundary layers over rough surfaces

The technique by Lund et al. to generate turbulent inflow for simulations of developing boundary layers over smooth flat plates is extended to the case of surfaces with roughness elements. In the Lund et al. method, turbulent velocities on a sampling plane are rescaled and recycled back to the inlet as inflow boundary condition. To rescale mean and fluctuating velocities, appropriate length scales need be identified and for smooth surfaces, the viscous scale l_ν = ν/u_τ (where ν is the kinematic viscosity and u_τ is the friction velocity) is employed for the inner layer. Different from smooth surfaces, in rough wall boundary layers the length scale of the inner layer, i.e. the roughness sub-layer scale l_d, must be determined by the geometric details of the surface roughness elements and the flow around them. In the proposed approach, it is determined by diagnosing dispersive stresses that quantify the spatial inhomogeneity caused by the roughness elements in the flow. The scale l_d is used for rescaling in the inner layer, and the boundary layer thickness δ is used in the outer region. Both parts are then combined for recycling using a blending function. Unlike the blending function proposed by Lund et al. which transitions from the inner layer to the outer layer at approximately 0.2δ, here the location of blending is shifted upwards to enable simulations of very rough surfaces in which the roughness length may exceed the height of 0.2δ assumed in the traditional method. The extended rescaling–recycling method is tested in large eddy simulation of flow over surfaces with various types of roughness element shapes.     

Inflow conditions for spatial direct numerical simulation of turbulent boundary layers

The inflow conditions for spatial direct numerical simulation (SDNS) of turbulent boundary layers should reflect the characteristics of upstream turbulence, which is a puzzle. In this paper a new method is suggested, in which the flow field obtained by using temporal direct numerical simulation (TDNS) for fully developed turbulent flow (only flow field for a single moment is sufficient) can be used as the inflow of SDNS with a proper transformation. The calculation results confirm that this method is feasible and effective. It is also found that, under a proper time-space transformation, all statistics of the fully developed turbulence obtained by both temporal mode and spatial mode DNS are in excellent agreement with each other, not only qualitatively, but also quantitatively. The normal-wise distributions of mean flow profile, turbulent Mach number and the root mean square (RMS) of the fluctuations of various variables, as well as the Reynolds stresses of the fully developed turbulence obtained by using SDNS, bear similarity in nature. Supported by the National Natural Science Foundation of China (Grant No. 90205021), the China Postdoctoral Science Foundation (Grant No. 20060400707), and the Foundation for the Author of National Excellent Doctoral Dissertation of China (Grant No. 200328), and partially supported by Liu-Hui Center of Applied Mathematics, Nankai University and Tianjin University     

Non—equilibrium extrapolation method for velocity and pressure boundary conditions in the lattice Boltzmann method

In this paper, we propose a new approach to implementing boundary conditions in the lattice Boltzmann method (LBM). The basic idea is to decompose the distribution function at the boundary node into its equilibrium and non-equilibrium parts, and then to approximate the non-equilibrium part with a first-order extrapolation of the non-equilibrium part of the distribution at the neighbouring fluid node. Schemes for velocity and pressure boundary conditions are constructed based on this method. The resulting schemes are of second-order accuracy. Numerical tests show that the numerical solutions of the LBM together with the present boundary schemes are in excellent agreement with the analytical solutions. Second-order convergence is also verified from the results. It is also found that the numerical stability of the present schemes is much better than that of the original extrapolation schemes proposed by Chen et al. (1996 Phys. Fluids 8 2527).     

POD-based wall boundary conditions for the numerical simulation of turbulent channel flows

Simulation of turbulent wall-bounded flows requires a high spatial resolution in the wall region, which limits the range of Reynolds numbers which can be effectively reached. In previous work, we proposed proper orthogonal decomposition (POD) based wall boundary conditions to bypass the simulation of the inner wall region. Tests were carried out for direct numerical simulation at a low Reynolds number Re_τ = 180. The boundary condition is based on the POD spatial eigenfunctions which are determined a priori in the full channel. It consists of a three-component velocity field on the plane y⁺ = 50 which is reconstructed at each instant from a combination of selected eigenfunctions. The coefficients of the combination are estimated from the simulation in the reduced domain using the threshold-based reconstruction method described in Podvin et al. The study is now extended to large-eddy simulation at higher Reynolds numbers Re_τ = 295 and Re_τ = 590. Two versions of the reconstruction method are considered. In the first version, both the phases and the moduli of the coefficients are allowed to vary. In the second version, only the phases are adjusted. We find that the latter method is associated with improved statistics and is relatively robust with respect to the reconstruction threshold. However, it is sensitive to the details of the numerical simulation, unlike the former method, which is associated with less accurate statistics and is more dependent on the reconstruction threshold.     

Characteristic boundary conditions for direct simulations of turbulent counterflow flames

Improved Navier–Stokes characteristic boundary conditions (NSCBC) are formulated for the direct numerical simulations (DNS) of laminar and turbulent counterflow flame configurations with a compressible flow formulation. The new boundary scheme properly accounts for multi-dimensional flow effects and provides nonreflecting inflow and outflow conditions that maintain the mean imposed velocity and pressure, while substantially eliminating spurious acoustic wave reflections. Applications to various counterflow configurations demonstrate that the proposed boundary conditions yield accurate and robust solutions over a wide range of flow and scalar variables, allowing high fidelity in detailed numerical studies of turbulent counterflow flames.     

Derivation of Slip boundary conditions for the Navier-Stokes system from the Boltzmann equation

Rarefied gas flow behavior is usually described by the Boltzmann equation, the Navier-Stokes system being valid when the gas is less rarefied. Slip boundary conditions for the Navier-Stokes equations are derived in a rigorous and systematic way from the boundary condition at the kinetic level (Boltzmann equation). These slip conditions are explicitly written in terms of asymptotic behavior of some linear half-space problems. The validity of this analysis is established in the simple case of the Couette flow, for which it is proved that the right boundary conditions are obtained.     

Fast and accurate point-based method for time-harmonic maxwell problems involving thin layer materials

We present a high-order hybrid boundary-finite elements method well-suited for solving time-harmonic electromagnetic scattering problems. Actually, this method is specially devoted to perfect electric conductors coated with a thin layer material. On such class of problems this method is shown to be fast and accurate. The fast feature is due to the joint use of finite elements of anisotropic order fitting the layer thickness, and of a point-based boundary element method on the skin. The accuracy is ensured, first by a discretization scheme satisfying the H_curl–H_div conformity required by the integro-differential equation and, secondly, by an adaptive technique of integration based on the detection of some local potential trouble on the geometry such as sharp edges or high dilatation of the elements. This algorithm does not need further information from the user and does not deteriorate the computation time. Numerical examples confirm the efficiency of this approach.     

A box-type scheme for fractional sub-diffusion equation with Neumann boundary conditions

Combining order reduction approach and L1 discretization, a box-type scheme is presented for solving a class of fractional sub-diffusion equation with Neumann boundary conditions. A new inner product and corresponding norm with a Sobolev embedding inequality are introduced. A novel technique is applied in the proof of both stability and convergence. The global convergence order in maximum norm is O(τ^2−α + h²). The accuracy and efficiency of the scheme are checked by two numerical tests.     

Artificial boundary conditions for the numerical solution of the Euler equations by the discontinuous galerkin method

We present artificial boundary conditions for the numerical simulation of compressible flows using high-order accurate discretizations with the discontinuous Galerkin (DG) finite element method. The construction of the proposed boundary conditions is based on characteristic analysis and applied for boundaries with arbitrary shape and orientation. Numerical experiments demonstrate that the proposed boundary treatment enables to convect out of the computational domain complex flow features with little distortion. In addition, it is shown that small-amplitude acoustic disturbances could be convected out of the computational domain, with no significant deterioration of the overall accuracy of the method. Furthermore, it was found that application of the proposed boundary treatment for viscous flow over a cylinder yields superior performance compared to simple extrapolation methods.     

A Cartesian method for fitting the bathymetry and tracking the dynamic position of the shoreline in a three-dimensional, hydrodynamic model

This paper presents a Cartesian method for the simultaneous fitting of the bathymetry and shorelines in a three-dimensional, hydrodynamic model for free-surface flows. The model, named LESS3D (Lake & Estuarine Simulation System in Three Dimensions), solves flux-based finite difference equations in the Cartesian-coordinate system (x,y,z). It uses a bilinear bottom to fit the bottom topography and keeps track the dynamic position of the shoreline. The resulting computational cells are hybrid: interior cells are regular Cartesian grid cells with six rectangular faces, and boundary/bottom cells (at least one face is the water–solid interface) are unstructured cells whose faces are generally not rectangular. With the bilinear interpolation, the shape of a boundary/bottom cell can be determined at each time step. This allows the Cartesian coordinate model to accurately track the dynamic position of the shorelines. The method was tested with a laboratory experiment of a Tsunami runup case on a circular island. It was also tested for an estuary in Florida, USA. Both model applications demonstrated that the Cartesian method is quite robust. Because the present method does not require any coordinate transformation, it can be an attractive alternative to curvilinear grid model.     

A primitive-variable Riemann method for solution of the shallow water equations with wetting and drying

A Riemann flux that uses primitive variables rather than conserved variables is developed for the shallow water equations with nonuniform bathymetry. This primitive-variable flux is both conservative and well behaved at zero depth. The unstructured finite-volume discretization used is suitable for highly nonuniform grids that provide resolution of complex geometries and localized flow structures. A source-term discretization is derived for nonuniform bottom that balances the discrete flux integral both for still water and in dry regions. This primitive-variable formulation is uniformly valid in wet and dry regions with embedded wetting and drying fronts. A fully nonlinear implicit scheme and both nonlinear and time-linearized explicit schemes are developed for the time integration. The implicit scheme is solved by a parallel Newton-iterative algorithm with numerically computed flux Jacobians. A concise treatment of characteristic-variable boundary conditions with source terms is also given. Computed results obtained for the one-dimensional dam break on wet and dry beds and for normal-mode oscillations in a circular parabolic basin are in very close agreement with the analytical solutions. Other results for a forced breaking wave with friction interacting with a sloped bottom demonstrate a complex wave motion with wetting, drying and multiple interacting wave fronts. Finally, a highly nonuniform, coastline-conforming unstructured grid is used to demonstrate an unsteady simulation that models an artificial coastal flooding due to a forced wave entering the Gulf of Mexico.     

An efficient method for solving highly anisotropic elliptic equations

Solving elliptic PDEs in more than one dimension can be a computationally expensive task. For some applications characterized by a high degree of anisotropy in the coefficients of the elliptic operator, such that the term with the highest derivative in one direction is much larger than the terms in the remaining directions, the discretized elliptic operator often has a very large condition number – taking the solution even further out of reach using traditional methods. This paper will demonstrate a solution method for such ill-behaved problems. The high condition number of the D-dimensional discretized elliptic operator will be exploited to split the problem into a series of well-behaved one and (D − 1)-dimensional elliptic problems. This solution technique can be used alone on sufficiently coarse grids, or in conjunction with standard iterative methods, such as Conjugate Gradient, to substantially reduce the number of iterations needed to solve the problem to a specified accuracy. The solution is formulated analytically for a generic anisotropic problem using arbitrary coordinates, hopefully bringing this method into the scope of a wide variety of applications.     

A method for determining foil thicknesses in TEM by using convergent beam electron diffraction under weak beam conditions

In the present work, convergent beam electron diffraction was studied in zirconium (a material of intermediate atomic number) at 300 keV, under weak beam diffraction conditions. For a particular thickness, the details in an observed low order disc were matched to those calculated using the multibeam dynamical theory. This presents the possibility of determining foil thickness over a wide range, with an estimated experimental accuracy of ≈7% or less. In contrast to other convergent beam techniques, the present method, which uses weak beam conditions, can employ commonly-occurring low order reflections to extract thicknesses.

A simple equation based on the two beam approximation, is derived to determine foil thickness (to within ≈ 10%) without resorting to detailed image matching. This equation can be used for a rough estimate of foil thickness while carrying out TEM observations.     

Application of the HN method to the critical slab problem for reflecting boundary conditions

The recently developed H_N method is used to solve the critical slab problem for a slab which is surrounded by a reflector. In the special case for R=0 (the reflection coefficient) the problem reduces to the one under vacuum boundary conditions. It is shown that the method is concise and leads to fast converging numerical results. The presented numerical results are compared with the data available in literature.     

A high-order cut-cell method for numerical simulation of hypersonic boundary-layer instability with surface roughness

Laminar-turbulent transition of hypersonic boundary layers can be affected significantly by the existence of surface roughness. Currently many important mechanisms of roughness-induced transition are not well understood. In recent years, direct numerical simulation (DNS) has been extensively applied for investigating instability and transition mechanisms of hypersonic boundary layers. Most of the past DNS studies, however, have been based on body-fitted grids for smooth surfaces without roughness. Due to complex geometry of embedded roughness, the use of body-fitted grids can be very difficult for flow with arbitrary surface roughness. In this paper, we present a new high-order cut-cell method to overcome the natural complexities in grid generation around arbitrary surface of roughness. The new method combines a non-uniform-grid finite-difference method for discrete grid points near the solid boundary and a shock-fitting method for the treatment of the bow shock. The non-uniform-grid finite-difference formulas are expressed in a general explicit form so that they can be applied to different multi-dimensional problems without any modification. The computational accuracy of new algorithms of up to O(h⁴) are tested on several one- and two-dimensional elliptic equations in irregular domains. In addition, the new method is applied to the simulation of the receptivity process of a Mach 5.92 flow over a flat plate under the combined effect of an isolated surface roughness element and surface blow and suction. A good agreement is found between the unsteady flow results and those obtained by a Linear Stability Theory (LST).     

A simple and accurate method for measuring program/erase speed in a memory capacitor structure

With the merits of a simple process and a short fabrication period, the capacitor structure provides a convenient way to evaluate memory characteristics of charge trap memory devices. However, the slow minority carrier generation in a capacitor often makes an underestimation of the program/erase speed. In this paper, illumination around a memory capacitor is proposed to enhance the generation of minority carriers so that an accurate measurement of the program/erase speed can be achieved. From the dependence of the inversion capacitance on frequency, a time constant is extracted to quantitatively characterize the formation of the inversion layer. Experimental results show that under a high enough illumination, this time constant is greatly reduced and the measured minority carrier-related program/erase speed is in agreement with the reported value in a transistor structure.     

A micro/nanocasting method for replicating soft micro/nanosurface structures using a dynamic electrical field

We propose a new electric field-induced micro/nanocasting method to replicate soft patterns using micro/nanocasting techniques without pressure. The process uses an alternating current (AC) electrical field and rotation of one electrode, generating a dynamic electrical field that induces electrokinetic flow motion in a dielectric solution (polydimethylsilane, PDMS). We used a lotus leaf as a replication template and characterised the PDMS flow motion to observe the effects of various process parameters (e.g., electrical field strength, rotation speed of an electrode, and electrode shape). The unstable flow motion was significantly dependent on the processing parameters, especially the rotation speed of the electrode. Using the optimised processing conditions, the replication efficiency was about 88%. We believe that this method has potential for fabricating soft micro/nanosized structures.     

A new technique for the measurement of the product CO/CO2 ratio at the surface of char particles burning in a fluidized bed

A new experimental technique is proposed to measure the product CO/CO₂ ratio at the surface of spherical char particles during fluidized bed combustion. It is based on the measurement of the burning rate of a single char particle under low oxygen concentration conditions and on the use of an accurate prediction of the particle Sherwood number. This technique was applied to spherical char particles obtained from a bituminous coal which is characterized by a low attrition and fragmentation propensity. The product CO/CO₂ ratio was measured at a bed temperature of 850 °C and at a fluidization velocity of 0.3 m/s in a lab-scale apparatus operated with a bed of 0.5–0.6 mm sand. The char particle size was varied between 2 and 7 mm and the inlet oxygen concentration between 0.1% and 2.0%. Results showed that under the experimental conditions investigated carbon was mostly oxidized to CO₂ within the particle boundary layer, with a maximum fraction of carbon escaping as CO of 10–20% at the lowest oxygen concentrations and largest particle sizes.