Reducing Undesired Wave Reflection at Domain Boundaries in 3D Finite Volume-Based Flow Simulations via Forcing Zones

2020 ◽  
Vol 64 (01) ◽  
pp. 23-47
Author(s):  
Robinson Peric ◽  
Moustafa Abdel-Maksoud
Keyword(s):  
Finite Volume ◽  
Wave Reflection ◽  
Three Dimensional ◽  
Navier Stokes ◽  
Special Focus ◽  
Computational Domain ◽  
Wave Reflections ◽  
Flow Simulations ◽  
Large Eddy ◽  
Domain Boundaries

This article reviews different types of forcing zones (sponge layers, damping zones, relaxation zones, etc.) as used in finite volume-based flow simulations to reduce undesired wave reflections at domain boundaries, with special focus on the case of strongly reflecting bodies subjected to long-crested incidence waves. Limitations and possible sources of errors are discussed. A novel forcing-zone arrangement is presented and validated via three-dimensional (3D) flow simulations. Furthermore, a recently published theory for predicting the forcing-zone behavior was investigated with regard to its relevance for practical 3D hydrodynamics problems. It was found that the theory can be used to optimally tune the case-dependent parameters of the forcing zones before running the simulations. 1. Introduction Wave reflections at the boundaries of the computational domain can cause significant errors in flow simulations, and must therefore be reduced. In contrast to boundary element codes, where much progress in this respect has been made decades ago (see e.g., Clement 1996; Grilli &Horillo 1997), for finite volume-based flow solvers, there are many unresolved questions, especially:How to reliably reduce reflections and disturbances from the domain boundaries?How to predict the amount of undesired wave reflection before running the simulation? This work aims to provide further insight to these questions for flow simulations based on Navier-Stokes-type equations (Reynolds-averaged Navier-Stokes, Euler equations, Large Eddy Simulations, etc.), when using forcing zones to reduce undesired reflections. The term "forcing zones" is used here to describe approaches that gradually force the solution in the vicinity of the boundary towards some reference solution, as described in Section 2; some examples are absorbing layers, sponge layers, damping zones, relaxation zones, or the Euler overlay method (Mayer et al. 1998; Park et al. 1999; Chen et al. 2006; Choi &Yoon 2009; Jacobsen et al. 2012; Kimet al. 2012; Schmitt & Elsaesser 2015; Perić & Abdel-Maksoud 2016a; Vukčević et al. 2016).

Large Eddy Simulation of Cross Flow in Pipe Junction

2018 ◽  
Author(s):  
Jiajun Chen ◽  
Yue Sun ◽  
Hang Zhang ◽  
Dakui Feng ◽  
Zhiguo Zhang
Keyword(s):  
Large Eddy Simulation ◽  
Stokes Equations ◽  
Numerical Study ◽  
Three Dimensional ◽  
Cross Flow ◽  
Exciting Force ◽  
Navier Stokes ◽  
Pipe Network ◽  
Eddy Simulation ◽  
Large Eddy

Mixing in pipe junctions can play an important role in exciting force and distribution of flow in pipe network. This paper investigated the cross pipe junction and proposed an improved plan, Y-shaped pipe junction. The numerical study of a three-dimensional pipe junction was performed for calculation and improved understanding of flow feature in pipe. The filtered Navier–Stokes equations were used to perform the large-eddy simulation of the unsteady incompressible flow in pipe. From the analysis of these results, it clearly appears that the vortex strength and velocity non-uniformity of centerline, can be reduced by Y-shaped junction. The Y-shaped junction not only has better flow characteristic, but also reduces head loss and exciting force. The results of the three-dimensional improvement analysis of junction can be used in the design of pipe network for industry.

scholarly journals Development of a Three-Dimensional Geometry Optimization Method for Turbomachinery Applications

2004 ◽  
Vol 10 (5) ◽  
pp. 373-385
Author(s):  
Steffen Kämmerer ◽  
Jürgen F. Mayer ◽  
Heinz Stetter ◽  
Meinhard Paffrath ◽  
Utz Wever ◽  
...  
Keyword(s):  
Optimization Problem ◽  
Stokes Equations ◽  
Turbine Blades ◽  
Three Dimensional ◽  
Optimization Techniques ◽  
Navier Stokes ◽  
Physical Parameters ◽  
Flow Solver ◽  
Flow Simulations ◽  
Sensitivity Equations

This article describes the development of a method for optimization of the geometry of three-dimensional turbine blades within a stage configuration. The method is based on flow simulations and gradient-based optimization techniques. This approach uses the fully parameterized blade geometry as variables for the optimization problem. Physical parameters such as stagger angle, stacking line, and chord length are part of the model. Constraints guarantee the requirements for cooling, casting, and machining of the blades.The fluid physics of the turbomachine and hence the objective function of the optimization problem are calculated by means of a three-dimensional Navier-Stokes solver especially designed for turbomachinery applications. The gradients required for the optimization algorithm are computed by numerically solving the sensitivity equations. Therefore, the explicitly differentiated Navier-Stokes equations are incorporated into the numerical method of the flow solver, enabling the computation of the sensitivity equations with the same numerical scheme as used for the flow field solution.This article introduces the components of the fully automated optimization loop and their interactions. Furthermore, the sensitivity equation method is discussed and several aspects of the implementation into a flow solver are presented. Flow simulations and sensitivity calculations are presented for different test cases and parameters. The validation of the computed sensitivities is performed by means of finite differences.

scholarly journals PARALLEL COMPUTATIONS IN STUDIES OF DEPENDENCE OF GAS DYNAMIC PARAMETERS OF UPWARD SWIRLING FLOW OF GAS ON BLOWING VELOCITY

2016 ◽  
pp. 92-97
Author(s):  
R. E. Volkov ◽  
A. G. Obukhov
Keyword(s):  
Stokes Equations ◽  
Swirling Flow ◽  
Initial Conditions ◽  
Exact Analytical Solution ◽  
Three Dimensional ◽  
Navier Stokes ◽  
Heat Conducting ◽  
Computational Domain ◽  
Navier Stokes Equations ◽  
Gas Dynamic

The rectangular parallelepiped explicit difference schemes for the numerical solution of the complete built system of Navier-Stokes equations. These solutions describe the three-dimensional flow of a compressible viscous heat-conducting gas in a rising swirling flows, provided the forces of gravity and Coriolis. This assumes constancy of the coefficient of viscosity and thermal conductivity. The initial conditions are the features that are the exact analytical solution of the complete Navier-Stokes equations. Propose specific boundary conditions under which the upward flow of gas is modeled by blowing through the square hole in the upper surface of the computational domain. A variant of parallelization algorithm for calculating gas dynamic and energy characteristics. The results of calculations of gasdynamic parameters dependency on the speed of the vertical blowing by the time the flow of a steady state flow.

scholarly journals Modeling the 3-d flow around a cylinder using the SAS hybrid model

2014 ◽  
Vol 16 (5) ◽  
pp. 901-918 ◽  
Keyword(s):  
Turbulence Modeling ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Cross Flow ◽  
Numerical Code ◽  
Navier Stokes ◽  
Computational Domain ◽  
Mean Flow ◽  
Unstable Flow ◽  
Vortex Size

<div> <p>Three-dimensional calculations were performed to simulate the flow around a cylindrical vegetation element using the Scale Adaptive Simulation (SAS) model; commonly, this is the first step of the modeling of the flow through multiple vegetation elements. SAS solves the Reynolds Averaged Navier-Stokes equations in stable flow regions, while in regions with unstable flow it goes unsteady producing a resolved turbulent spectrum after reducing eddy viscosity according to the locally resolved vortex size represented by the von Karman length scale. A finite volume numerical code was used for the spatial discretisation of the rectangular computational domain with stream-wise, cross-flow and vertical dimensions equal to 30D, 11D and 1D, respectively, which was resolved with unstructured grids. Calculations were compared with experiments and Large Eddy Simulations (LES). Predicted overall flow parameters and mean flow velocities exhibited a very satisfactory agreement with experiments and LES, while the agreement of predicted turbulent stresses was satisfactory. Calculations showed that SAS is an efficient and relatively fast turbulence modeling approach, especially in relevant practical problems, in which the very high accuracy that can be achieved by LES at the expense of large computational times is not required.</p> </div> <p>&nbsp;</p>

scholarly journals Three-Dimensional Turbulent Vortex Shedding From a Surface-Mounted Square Cylinder: Predictions With Large-Eddy Simulations and URANS

2014 ◽  
Vol 136 (6) ◽  
Author(s):  
B. A. Younis ◽  
A. Abrishamchi
Keyword(s):  
Experimental Data ◽  
Turbulence Model ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Large Eddy Simulations ◽  
Navier Stokes ◽  
Reynolds Stresses ◽  
Mean Flow ◽  
Navier Stokes Equations ◽  
Large Eddy

The paper reports on the prediction of the turbulent flow field around a three-dimensional, surface mounted, square-sectioned cylinder at Reynolds numbers in the range 104–105. The effects of turbulence are accounted for in two different ways: by performing large-eddy simulations (LES) with a Smagorinsky model for the subgrid-scale motions and by solving the unsteady form of the Reynolds-averaged Navier–Stokes equations (URANS) together with a turbulence model to determine the resulting Reynolds stresses. The turbulence model used is a two-equation, eddy-viscosity closure that incorporates a term designed to account for the interactions between the organized mean-flow periodicity and the random turbulent motions. Comparisons with experimental data show that the two approaches yield results that are generally comparable and in good accord with the experimental data. The main conclusion of this work is that the URANS approach, which is considerably less demanding in terms of computer resources than LES, can reliably be used for the prediction of unsteady separated flows provided that the effects of organized mean-flow unsteadiness on the turbulence are properly accounted for in the turbulence model.

Dielectrophoretic Mixing With Novel Electrode Geometry

2009 ◽  
Author(s):  
G. Naga Siva Kumar ◽  
Sushanta K. Mitra ◽  
Subir Bhattacharjee
Keyword(s):  
Electric Field ◽  
Cell Activation ◽  
Colloidal Particles ◽  
Colloidal Suspension ◽  
Three Dimensional ◽  
Colloidal Suspensions ◽  
Navier Stokes ◽  
Mixing Efficiency ◽  
Computational Domain ◽  
Element Analysis

Electrokinetic mixing of analytes at micro-scale is important in several biochemical applications like cell activation, DNA hybridization, protein folding, immunoassays and enzyme reactions. This paper deals with the modeling and numerical simulation of micromixing of two different types of colloidal suspensions based on principle of dielectrophoresis (DEP). A mathematical model is developed based on Laplace, Navier-Stokes, and convection-diffusion-migration equations to calculate electric field, velocity, and concentration distributions, respectively. Mixing of two colloidal suspensions is simulated in a three-dimensional computational domain using finite element analysis considering dielectrophoretic, gravitational and convective (advective)–diffusive forces. Phase shifted AC signal is applied to the alternating electrodes for achieving the mixing of two different colloidal suspensions. The results indicate that the electric field and DEP forces are maximum at the edges of the electrodes and become minimum elsewhere. As compared to curved edges, straight edges of electrodes have lower electric field and DEP forces. The results also indicate that DEP force decays exponentially along the height of the channel. The effect of DEP forces on the concentration profile is studied. It is observed that, the concentration of colloidal particles at the electrodes edges is very less compared to elsewhere. Mixing of two colloidal suspensions due to diffusion is observed at the interface of the two suspensions. The improvement in mixing after applying the repulsive DEP forces on the colloidal suspension is observed. Most of the mixing takes place across the slant edges of the triangular electrodes. The effect of electrode pairs and the mixing length on degree of mixing efficiency are also observed.

Study of Afterburning in Solid Rocket Booster Jet: Towards a Model for Afterburning in Large Eddy Simulations

2014 ◽  
Author(s):  
Adèle Poubeau ◽  
Roberto Paoli ◽  
Daniel Cariolle
Keyword(s):  
Three Dimensional ◽  
Axial Direction ◽  
Large Eddy Simulations ◽  
Chemical Mechanism ◽  
Computational Domain ◽  
Time Step ◽  
Simple Method ◽  
Large Eddy ◽  
Solid Rocket ◽  
Comparison Of The Results

This paper focuses on two decisive steps towards Large Eddy Simulation of a solid rocket booster jet. First, three-dimensional Large Eddy Simulations of a non-reactive booster jet including the nozzle were obtained at flight conditions of 20 km of altitude. A particularly long computational domain (400 nozzle exit diameters in the jet axial direction) was simulated, thanks to an innovative local time-stepping method via coupling multi instances of a fluid solver. The dynamics of the jet is analysed and comparison of the results with previous knowledge validates the simulations and confirms that this computational setup can be applied for Large Eddy Simulations of a reactive booster jet. The second part of this paper details the implementation of a simple method to study the hot plume chemistry. Despite its limitations, it is accurate enough to observe the various steps of the chemical mechanism and assess the effect of uncertainties of the rate parameters on chlorine reactions. It was also used to reduce the set of chemical reactions into a short scheme involving a minimum of species and having a limited impact on the physical time step of the Large Eddy Simulations.

scholarly journals ICCM2015: A Comparison of RANS and LES Computational Methods in Analyzing Ventilation Flow Through a Room Fitted with a Two-Sided Windcatcher

2017 ◽  
Vol 14 (03) ◽  
pp. 1750021 ◽  
Author(s):  
A. Niktash ◽  
B. P. Huynh
Keyword(s):  
Natural Ventilation ◽  
Three Dimensional ◽  
Navier Stokes ◽  
Interior Space ◽  
Eddy Simulation ◽  
Large Eddy ◽  
Computational Fluid Dynamics Cfd ◽  
Flow Through ◽  
Cfd Method ◽  
Good Agreement

A windcatcher is a structure for providing natural ventilation using wind power; it is usually fitted on the roof of a building to exhaust the inside stale air to the outside and supplies the outside fresh air into the building interior space working by pressure difference between outside and inside of the building. In this paper, the behavior of free wind flow through a three-dimensional room fitted with a centered position two-canal bottom shape windcatcher model is investigated numerically, using a commercial computational fluid dynamics (CFD) software package and LES (Large Eddy Simulation) CFD method. The results have been compared with the obtained results for the same model but using RANS (Reynolds Averaged Navier–Stokes) CFD method. The model with its surrounded space has been considered in both method. It is found that the achieved results for the model from LES method are in good agreement with RANS method’s results for the same model.

Sign in / Sign up

Export Citation Format

Share Document