Consistent inlet and outlet boundary conditions for particle methods

In this paper, simple and consistent open boundary conditions are presented for the numerical simulation of viscous incompressible laminar flows. The present approach is based on an arbitrary Lagrangian-Eulerian particle method using upwind interpolation. Three kinds of inlet/outlet boundary conditions are proposed for particle methods, a pressure specified inlet/outlet condition, a velocity profile specified inlet/outlet condition, and a fully developed flow outlet condition. These inlet/outlet conditions are realized by using boundary particles and modification to the physical value such as velocity. Poiseuille flows, flows over a backward-facing step, and flows in a T-shape branch are calculated. The results are compared with those of mesh-based methods such as the finite volume method. The method presented herein exhibits accuracy and numerical stability.     

Cross influence of discharge and circulation on head loss of conduit of pump system with low head

The relationship between the head loss and the discharge and circulation of the conduit of a pump system with low head is an important problem with an obvious influence on the improvement of its hydraulic performance. The velocity circulation from the pump guide vane makes the relationship more complicated, which has to be understood comprehensively. The results indicate that, under the condition of zero circulation, the head loss of the inlet and outlet conduits is in proportion to the square of discharge. Under the condition that the Reynolds number is satisfied with the resistant square area, the conduit loss is in proportion to the square of discharge for the similar working points with different speeds in a certain rotational speed range, indicating that the pump system efficiency is constant. The outlet conduit loss of design discharge for a pump system with low head depends on the velocity circulation from the guide vane exit, and the relationship between the loss and the circulation is an open curve with an upward direction, meaning that there is an optimal circulation for the loss. Under the condition of various working points for a pump system with low head, the head loss of the outlet conduit is under the cross influence of both the discharge and the circulation. As a result, the relationship between the head loss and the discharge is almost linear, and the mechanism needs to be further studied.     

Characteristic‐based inlet and outlet boundary conditions for incompressible flows

Determining boundary conditions (BCs) for incompressible flows is such a delicate matter that affects the accuracy of the results. In this research, a new characteristic‐based BC for incompressible Navier‐Stokes equations is introduced. Discretization of equations has been done via finite volume. Additionally, artificial compressibility correction has been employed to deal with equations. Ordinary extrapolation from inner cells of a domain was used as a traditional way to estimate pressure and velocities on solid wall and inlet/outlet boundaries. Here, this method was substituted by the newly proposed BCs based on the characteristics of artificial compressibility equations. To follow this purpose, a computer code has been developed to carry out series of numerical tests for a flow over a backward‐facing step and was applied to a wide range of Reynolds numbers and grid combinations. Calculation of convective and viscous fluxes was done using Jameson's averaging scheme. Employing the characteristic‐based method for determining BCs has shown an improved convergence rate and reduced calculation time comparing with those of traditional ones. Furthermore, with the reduction of domain and computational cells, a similar accuracy was achieved for the results in comparison with the ones obtained from the traditional extrapolation method, and these results were in good agreement with the ones in the literature.     

改进的SPH边界处理方法与土体大变形模拟

沙土滑坡往往会造成重大的人身财产损失,研究这类土体大变形问题对防灾工程具有指导意义。光滑粒子流体动力学SPH（Smoothed Particle Hydrodynamics）方法是一种拉格朗日型无网格粒子法,十分适用于模拟大变形问题。在SPH方法中,合适的边界处理方法一直是个难点,传统的边界虚粒子法或排斥力法较难模拟复杂边界。本文引入了一种能处理任意形状边界的方法——统一半解析壁面边界条件处理方法USAW（unified semi-analytical wall boundary conditions）,通过在控制方程中引入修正因子并保留边界面积分项来弥补边界缺失。为了更准确模拟问题域边界,提出无质量边界粒子的新概念。利用该方法成功模拟了土体滑坡算例,验证了方法的可靠性,并避免了边界零粒子层问题。通过数值模拟,分析了内摩擦角和黏聚力等土体物性参数对滑坡过程的影响。最后,应用该方法研究了滑坡冲击楔形体时的压力响应。     

Three-dimensional flow of Oldroyd-B fluid over surface with convective boundary conditions

The present study addresses the three-dimensional flow of an Oldroyd-B fluid over a stretching surface with convective boundary conditions. The problem formulation is presented using the conservation laws of mass, momentum, and energy. The solutions to the dimensionless problems are computed. The convergence of series solutions by the homotopy analysis method (HAM) is discussed graphically and numerically. The graphs are plotted for various parameters of the temperature profile. The series solutions are verified by providing a comparison in a limiting case. The numerical values of the local Nusselt number are analyzed.     

A boundary condition for arbitrary shaped inlets in lattice–Boltzmann simulations

We introduce a mass‐flux‐based inlet boundary condition for the lattice‐Boltzmann method. The proposed boundary condition requires minimal amount of boundary data, it produces a steady‐state velocity field which is accurate close to the inlet even for arbitrary inlet geometries, and yet it is simple to implement. We demonstrate its capability for both simple and complex inlet geometries by numerical experiments. For simple inlet geometries, we show that the boundary condition provides very accurate inlet velocities when Re?1. Even with moderate Reynolds number, the inlet velocities are accurate for practical purposes. Furthermore, the potential of our boundary condition to produce inlet velocities which convincingly adapt to complex inlet geometries is highlighted with two specific examples. Copyright © 2009 John Wiley & Sons, Ltd.     

Characteristics of the critical heat flux for downward flow in a vertical tube at low flow rate and low pressure conditions

The characteristics of the critical heat flux (CHF) for downward flow were studied experimentally with an Inconel 600 circular tube test section in a water test loop at low-flow rate (0 200 kg/m²s) and low-pressure (0.1 0.7 MPa) conditions. The attention was given to the effects of upstream conditions—upper plenum and inlet throttling. Two totally different kinds of CHF behaviors were observed. It seems appropriate to interpret them as flooding-type CHF and dryout in annular flow. The CHF in downward flow may vary from extremely unstable flow CHF as low as near the flooding CHF value to stable flow CHF as high as that of upflow, depending on the upstream conditions of the test section. The CHF correlation by Mishima and that by Weber were proposed for the presentation of the lower and upper limits of the CHF for downward flow in a vertical tube at low-flow rate and low-pressure conditions.     

Flow of fluids with pressure dependent viscosities in an orthogonal rheometer subject to slip boundary conditions

We consider slow steady flows of Navier–Stokes-like fluids with pressure dependent viscosities between rotating infinite parallel plates with Navier slip boundary conditions. We derive exact solutions which correspond to flows in orthogonal and torsional rheometers, and investigate the effect of the slip coefficient and the material parameters on the solutions. We find that even when inertial effects are ignored, vorticity boundary layers develop at the upper boundary due to the pressure dependence of the viscosity. These boundary layers diminish and eventually disappear with increased slippage.     

On the implementation of symmetric and antisymmetric periodic boundary conditions for incompressible flow

In this paper we consider symmetric and antisymmetric periodic boundary conditions for flows governed by the incompressible Navier-Stokes equations. Classical periodic boundary conditions are studied as well as symmetric and antisymmetric periodic boundary conditions in which there is a pressure difference between inlet and outlet. The implementation of this type of boundary conditions in a finite element code using the penalty function formulation is treated and also the implementation in a finite volume code based on pressure correction. The methods are demonstrated by computation of a flow through a staggered tube bundle.     

Solutions for a class of Hamiltonian systems on time scales with non-local boundary conditions

In this work, the solvability of a class of second-order Hamiltonian systems on time scales is generalized to non-local boundary conditions. The measurements obtained by non-local conditions are more accurate than those given by local conditions in some problems. Compared with the known results, this work establishes the variational structure in an appropriate Sobolev’s space. Then, by applying the mountain pass theorem and symmetric mountain pass theorem, the existence and multiplicity of the s...     

Singular perturbation of general boundary value problem for nonlinear differential equation system

I.IntroductionAllphysicalsystemsarenonlineartosomeextent.Actually,Lineal.systemisimaginarymodelwherenonlinearfactorisomittedinnonlinearsystem.Insolvingtheautocontl'ol,nonlinearoscillationtheory,theboundarystagnationproblenloffluidIncchanicsandsollleproblemsofsemi-conducttheoryandquantummechanicsetc'.weOnlyncedtosolvethefollowingproblem,whichisnonlineardifferentialequationsystemwithtilesll,allparanletel'inhighestorderderivativeandnonlinearboundaryconditions.whereE>0isasmallparameter,teR,x,fi…     

Mathematically consistent boundary conditions and turbulence matching at block interfaces for computational aeroacoustics

The method of characteristics is used to implement the various boundary conditions (e.g. wall and interface) in a high-order computational aeroacoustic (CAA) code developed by the first author. Most characteristic methods do not satisfy Pfaff's condition (which needs to be satisfied for any mathematical relation to be valid). A mathematically consistent and valid method is used in this work to derive the characteristic boundary conditions. Also, a robust and efficient approach for the matching of turbulence quantities at multi-block interfaces is proposed. Various numerical simulation cases were run to validate the concepts. The computed results show that the proposed method is accurate, robust and is in excellent agreement with experimental data. The results also indicate that the matching of turbulence quantities is essential for accurate turbulent flow calculations.     

Improvement of boundary conditions for non-planar boundaries represented by polygons with an initial particle arrangement technique

The boundary conditions represented by polygons in moving particle semi-implicit (MPS) method (Koshizuka and Oka, Nuclear Science and Engineering, 1996 Koshizuka, S., and Y. Oka. 1996. “Moving-Particle Semi-Implicit Method for Fragmentation of Incompressible Fluid.” Nuclear Science and Engineering 123: 421–434.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) have been widely used in the industry simulations since it can simply simulate complex geometry with high efficiency. However, the inaccurate particle number density near non-planar wall boundaries dramatically affects the accuracy of simulations. In this paper, we propose an initial boundary particle arrangement technique coupled with the wall weight function method (Zhang et al. Transaction of JSCES, 2015 Zhang, T.G., S. Koshizuka, K. Shibata, K. Murotani, and E. Ishii. 2015. “Improved Wall Weight Function With Polygon Boundary in Moving Particle Semi-Implicit Method.” Transaction of JSCES. No. 20150012. [Google Scholar]) to improve the particle number density near slopes and curved surfaces with boundary conditions represented by polygons in three dimensions. Two uniform grids are utilized in the proposed technique. The grid points in the first uniform grid are used to construct boundary particles, and the second uniform grid stores the same information as in the work by Zhang et al. The wall weight functions of the grid points in the second uniform grid are calculated by newly constructed boundary particles. The wall weight functions of the fluid particles are interpolated from the values stored on the grid points in the second uniform grid. Because boundary particles are located on the polygons, complex geometries can be accurately represented. The proposed method can dramatically improve the particle number density and maintain the high efficiency. The performance of the previously proposed wall weight function (Zhang et al.) with the boundary particle arrangement technique is verified in comparison with the wall weight function without boundary particle arrangement by investigating two example geometries. The simulations of a water tank with a wedge and a complex geometry show the general applicability of the boundary particle arrangement technique to complex geometries and demonstrate its improvement of the wall weight function near the slopes and curved surfaces.     

High-order boundary conditions for the problems of Laplace equation in infinite region and their application

The high-order boundary conditions for the problems of Laplace equation in infiniteregion have been developed.The improvement in accuracy for numerical solution isachieved by imposing the high-order boundary conditions on the exterior boundary of areduced finite region in which the numerical method is used.So both the computing effortsand the required storage in computer are reduced.The numerical examples show that thelst-order boundary condition approaches to the exact boundary condition and it is clearlysuperior to the traditional boundary condition and the2nd-order boundary condition.     

Adomian-modified decomposition method for large-amplitude vibration analysis of stepped beams with elastic boundary conditions

The main objective of this paper is to apply an Adomian modified decomposition method for solving large amplitude vibration analysis of stepped beams with various general and elastic boundary conditions. Damaged or imperfect supports of beams can be modeled by using elastic boundary conditions composing of translational and rotational springs. For the beams subjected to dynamic severe loading, it is important to include the nonlinear term of axial stretching force developed by the large vibration amplitude in the governing equation for more accurate design. By using the method, the convergence studies for linear and nonlinear vibration analyses of stepped beams are shown for determining an appropriate number of terms in the solutions. The accuracy of the present results is validated numerically by comparing with some available results in the literature. New results of nonlinear frequency ratios of stepped beams with different boundary conditions are presented and discussed in detail. Aspects of step ratio, step location, boundary conditions, vibration amplitudes, etc., which have significant impact on linear and nonlinear frequencies of such beams are taken under investigation.     

DNS of turbulent channel flow with conjugate heat transfer: Effect of thermal boundary conditions on the second moments and budgets

Budgets of turbulent heat fluxes and temperature variance obtained from the Direct Numerical Simulation of an incompressible periodic channel flow with a Reynolds number of 150 (based on friction velocity) and a Prandtl number of 0.71 are presented and analysed for four cases: locally imposed temperature at the wall (constant Dirichlet), locally imposed heat flux (constant Neumann), heat exchange coefficient (Robin) and 3D conjugate heat transfer. The dissipation rate associated with the temperature variance is strongly impacted by the thermal boundary condition. For non-conjugate cases, a straightforward analytical analysis establishes the connection between the boundary condition, the temperature variance and the wall-normal part of the thermal dissipation rate at the wall. For the conjugate case, the two-point correlations of the thermal field in the solid domain confirms the existence of very large scale thermal structures.     

On the stability and dissipation of wall boundary conditions for compressible flows

Characteristic formulations for boundary conditions have demonstrated their effectiveness to handle inlets and outlets, especially to avoid acoustic wave reflections. At walls, however, most authors use simple Dirichlet or Neumann boundary conditions, where the normal velocity (or pressure gradient) is set to zero. This paper demonstrates that there are significant differences between characteristic and Dirichlet methods at a wall and that simulations are more stable when using walls modelled with a characteristic wave decomposition. The derivation of characteristic methods yields an additional boundary term in the continuity equation, which explains their increased stability. This term also allows to handle the two acoustic waves going towards and away from the wall in a consistent manner. Those observations are confirmed by stability matrix analysis and one‐ and two‐dimensional simulations of acoustic modes in cavities. Copyright © 2009 John Wiley & Sons, Ltd.     

Flow and heat transfer of nanofluid past stretching/shrinking sheet with partial slip boundary conditions

The boundary layer flow of a nanofluid past a stretching/shrinking sheet with hydrodynamic and thermal slip boundary conditions is studied. Numerical solutions to the governing equations are obtained using a shooting method. The results are found for the skin friction coefficient, the local Nusselt number, and the local Sherwood number as well as the velocity, temperature, and concentration profiles for some values of the velocity slip parameter, thermal slip parameter, stretching/shrinking parameter, thermophoresis parameter, and Brownian motion parameter. The results show that the local Nusselt number, which represents the heat transfer rate, is lower for higher values of thermal slip parameter, thermophoresis parameter, and Brownian motion parameter.     

High‐order boundary conditions for linearized shallow water equations with stratification,dispersion and advection

The two‐dimensional linearized shallow water equations are considered in unbounded domains with density stratification. Wave dispersion and advection effects are also taken into account. The infinite domain is truncated via a rectangular artificial boundary ??, and a high‐order open boundary condition (OBC) is imposed on ??. Then the problem is solved numerically in the finite domain bounded by ??. A recently developed boundary scheme is employed, which is based on a reformulation of the sequence of OBCs originally proposed by Higdon. The OBCs can easily be used up to any desired order. They are incorporated here in a finite difference scheme. Numerical examples are used to demonstrate the performance and advantages of the computational method, with an emphasis is on the effect of stratification. Published in 2004 by John Wiley & Sons, Ltd.