Non-reflective boundary conditions for non-uniform mean velocity inlet and outlet flows

The numerical solution of problems in unbounded physical space requires a truncation of the computational domain to a reasonable size. As a result, the conditions on the artificial boundaries are generally unknown. Assumptions like constant pressure or velocities are only valid in the far field and lead to spurious reflections if applied on the boundaries of the truncated domain. A number of attempts have been made over the past decades to design conditions that prevent such reflections. One approach is based on characteristics. The standard analysis assumes a spatially uniform mean flow field but this is often impractical. In the present paper we show how to extend the formulation to the more general case of a non-uniform mean velocity field. A number of test cases are provided and our results compare favourably with other boundary conditions. In principle the present approach can be extended to include non-uniformities in all variables.

Off-wall boundary conditions for turbulent flows obtained from buffer-layer minimal flow units

There is strong evidence that the transport processes in the buffer region of wall-bounded turbulence are common across various flow configurations, even in the embryonic turbulence in transition (Park et al., Phys. Fl. 24). We use this premise to develop off-wall boundary conditions for turbulent simulations. Boundary conditions are constructed from DNS databases using periodic minimal flow units and reduced order modeling. The DNS data was taken from a channel at Reτ=400 and a zero-pressure gradient transitional boundary layer (Sayadi et al., submitted to J. Fluid Mech.). Both types of boundary conditions were first tested on a DNS of the core of the channel flow with the aim of extending their application to LES and to spatially evolving flows.

Direct numerical simulation of turbulent flows over superhydrophobic surfaces with gas pockets using linearized boundary conditions

Superhydrophobic surfaces are shown to be effective for surface drag reduction under laminar regime by both experiments and simulations (see for example, Ou and Rothstein, Phys. Fluids 17:103606, 2005). However, such drag reduction for fully developed turbulent flow maintaining the Cassie-Baxter state remains an open problem due to high shear rates and flow unsteadiness of turbulent boundary layer. Our work aims to develop an understanding of mechanisms leading to interface breaking and loss of gas pockets due to interactions with turbulent boundary layers. We take advantage of direct numerical simulation of turbulence with slip and no-slip patterned boundary conditions mimicking the superhydrophobic surface. In addition, we capture the dynamics of gas-water interface, by deriving a proper linearized boundary condition taking into account the surface tension of the interface and kinematic matching of interface deformation and normal velocity conditions on the wall. We will show results from our simulations predicting the dynamical behavior of gas pocket interfaces over a wide range of dimensionless surface tensions.

Boundary conditions in physical model tests-the influence of deck pinning on the response of piled bridge abutments in laterally spreading soils

Physical models are widely used in the study of geotechnical earthquake engineering phenomena, and the comparison of modelling results to observations from field reconnaissance provides a transparent means of evaluating the design of our physical models. This paper compares centrifuge tests of pile groups in laterally spreading slopes with the response of piled bridge abutments in the 2011 Christchurch earthquake. We show that the model foundation's fixity conditions strongly affect the success with which the mechanism of response of the real abutments is replicated in the tests. © 2012 American Society of Civil Engineers.

The influence of boundary conditions on tip leakage flow

Most of the current understanding of tip leakage flows has been derived from detailed cascade experiments. However, the cascade model is inherently approximate since it is difficult to simulate the boundary conditions present in a real machine, particularly the secondary flows convecting from the upstream stator row and the relative motion of the casing and blade. This problem is further complicated when considering the high pressure turbine rotors of aero engines, where the high Mach numbers must also be matched in order to correctly model the aerodynamics and heat transfer. More realistic tests can be performed on high-speed turbines, but the experimental fidelity and resolution achievable in such set-ups is limited. In order to examine the differences between cascade models and real-engine behavior, the influence of boundary conditions on the tip leakage flow in an unshrouded high pressure turbine rotor is investigated using RANS calculations. This study examines the influence of the rotor inlet condition and relative casing motion. A baseline calculation with a simplified inlet condition and no relative endwall motion exhibits similar behavior to cascade studies. Only minor changes to the leakage flow are induced by introducing either a more realistic inlet condition or relative casing motion. However when both of these conditions are applied simultaneously the pattern of leakage flow is very different, with ingestion of flow over much of the early suction surface. The paper explores the physical processes driving this change and the impact on leakage losses and modeling requirements. Copyright © 2013 by ASME.

Scaling of velocity fluctuations in off-wall boundary conditions for turbulent flows

A model for off-wall boundary conditions for turbulent flow is investigated. The objective of such a model is to circumvent the need to resolve the buffer layer near the wall, by providing conditions in the logarithmic layer for the overlying flow. The model is based on the self-similarity of the flow at different heights in the logarithmic layer. It was first proposed by Mizuno and Jiménez (2013), imposing at the boundary plane a velocity field obtained on-the-fly from an overlying region. The key feature of the model was that the lengthscales of the field were rescaled to account for the self-similarity law. The model was successful at sustaining a turbulent logarithmic layer, but resulted in some disagreements in the flow statistics, compared to fully-resolved flows. These disagreements needed to be addressed for the model to be of practical application. In the present paper, a more refined, wavelength-dependent rescaling law is proposed, based on the wavelength-dependent dynamics in fully-resolved flows. Results for channel flow show that the new model eliminates the large artificial pressure fluctuations found in the previous one, and a better agreement is obtained in the bulk properties, the flow fluctuations, and their spectral distribution across the whole domain. © Published under licence by IOP Publishing Ltd.

Boundary conditions at the solid-liquid surface over the multiscale channel size from nanometer to micron

The boundary condition at the solid surface is one of the important problems for the microfluidics. In this paper we study the effects of the channel sizes on the boundary conditions (BC), using the hybrid computation scheme adjoining the molecular dynamics (MD) simulations and the continuum fluid mechanics. We could reproduce the three types of boundary conditions (slip, no-slip and locking) over the multiscale channel sizes. The slip lengths are found to be mainly dependent on the interfacial parameters with the fixed apparent shear rate. The channel size has little effects on the slip lengths if the size is above a critical value within a couple of tens of molecular diameters. We explore the liquid particle distributions nearest the solid walls and found that the slip boundary condition always corresponds to the uniform liquid particle distributions parallel to the solid walls, while the no-slip or locking boundary conditions correspond to the ordered liquid structures close to the solid walls. The slip, no-slip and locking interfacial parameters yield the positive, zero and negative slip lengths respectively. The three types of boundary conditions existing in "microscale" still occur in "macroscale". However, the slip lengths weakly dependent on the channel sizes yield the real shear rates and the slip velocity relative to the solid wall traveling speed approaching those with the no-slip boundary condition when the channel size is larger than thousands of liquid molecular diameters for all of the three types of interfacial parameters, leading to the quasi-no-slip boundary conditions.

A New Criterion Number for the Boundary Conditions at the Solid/Liquid Interface in Nanoscale

A simple, but important three-atom model was proposed at the solid/liquid interface, leading to a new criterion number, lambda, governing the boundary conditions (BCs) in nanoscale. The solid wall is considered as the face-centered-cubic (fcc) structure. The fluid is the liquid argon with the well-known LJ potential. Based on the concept, the two micro-systems have the same BCs if they have The same criterion number. The degree of the locking BCs is enhanced when lambda equals to 0.757. Such critical criterion number results in the substantial epitaxial ordering and one, two, or even three liquid layers are locked by the solid wall, depending on the coupling energy scale ratio of the solid and liquid atoms. With deviation from the critical criterion number, the flow approaches the slip BCs and there are little ordering structures within the liquid. Always at the same criterion number, the degree of the slip is decreased or the locking is enhanced with increasing the coupling energy scale ratio of the solid and liquid atoms. The above analysis is well confirmed by the molecular dynamics (MD) simulation. The slip length is well correlated in terms of the new criterion number. The future work is suggested to extend the present theory for other microstructures of the solid wall atoms and quasi-LJ potentials.

Boundary Element Method (BEM) Analysis for Galvanic Corrosion of Hot Dip Galvanized Steel Immersed in Seawater

A numerical analysis of galvanic corrosion of hot-dip galvanized steel immersed in seawater was presented. The analysis was based on the boundary element methods (BEMs) coupled with Newton-Raphson iterative technique to treat the nonlinear boundary conditions, which were determined by the experimental polarization curves. Results showed that galvanic current density concentrates on the boundary of steel substrate and zinc coating, and the sacrificial protection of zinc coating to steel substrate results in overprotection of steel cathode. Not only oxygen reduction but also hydrogen reduction could occur as cathode reactions, which probably led up to the adsorption and absorption of hydrogen atoms. Flat galvanized steel tensile sample shows a brittle behavior similar to hydrogen embrittlement according to the SSRT (show strain rate test) in seawater.

Monotone scheme and boundary conditions for finite volume simulation of magnetohydrodynamic internal flows at high Hartmann number

A monotone scheme for finite volume simulation of magnetohydrodynamic internal flows at high Hartmann number is presented. The numerical stability is analysed with respect to the electromagnetic force. Standard central finite differences applied to finite volumes can only be numerically stable if the vector products involved in this force are computed with a scheme using a fully staggered grid. The electromagnetic quantities (electric currents and electric potential) must be shifted by half the grid size from the mechanical ones (velocity and pressure). An integral treatment of the boundary layers is used in conjunction with boundary conditions for electrically conducting walls. The simulations are performed with inhomogeneous electrical conductivities of the walls and reach high Hartmann numbers in three-dimensional simulations, even though a non-adaptive grid is used.

Multi-zone TAP-reactors theory and application - IV. Ideal and non-ideal boundary conditions

We study the influence of non-ideal boundary and initial conditions (BIC) of a temporal analysis of products (TAP) reactor model on the data (observed exit flux) analysis. The general theory of multi-response state-defining experiments for a multi-zone TAP reactor is extended and applied to model several alternative boundary and initial conditions proposed in the literature. The method used is based on the Laplace transform and the transfer matrix formalism for multi-response experiments. Two non-idealities are studied: (1) the inlet pulse not being narrow enough (gas pulse not entering the reactor in Dirac delta function shape) and (2) the outlet non-ideality due to imperfect vacuum. The effect of these non-idealities is analyzed to the first and second order of approximation. The corresponding corrections were obtained and discussed in detail. It was found that they are negligible. Therefore, the model with ideal boundary conditions is proven to be completely adequate to the description and interpretation of transport-reaction data obtained with TAP-2 reactors.

Solving the problem of diffraction of an optical-band electromagnetic wave by metal nanostructured aperture arrays with the use of the method of impedance boundary conditions

The problem of diffraction of an optical wave by a 2D periodic metal aperture array with square, circular, and ring apertures is solved with allowance for the finite permittivity of a metal in the optical band. The correctness of the obtained results is verified through comparison with experimental data. It is shown that the transmission coefficient can be substantially greater than the corresponding value reached in the case of diffraction by a grating in a perfectly conducting screen.