scholarly journals Architecture of Open Source Program for Numerical Modeling of Flows on Mountain Slopes

2020 ◽  
Vol 32 (6) ◽  
pp. 183-200
Author(s):  
Daria Igorevna Romanova

In this paper, we compare two approaches to describe the dynamics of flows on mountain slopes using the depth-averaged equations of continuum mechanics and using the complete, not depth-averaged, equations of continuum mechanics for three-dimensional modeling. Using these two approaches, a simulation of an experimental slush flow in the tank and the interaction of the flow with dam barrier protection was carried out. Numerical solutions are compared with experimental data. Also, both approaches are applied to the calculation of an avalanche in the 22nd avalanche cite of Mount Yukspor (Khibiny). Avalanche run-out distance and the shape of the avalanche deposits are compared with field data obtained from the measurement of a real avalanche. In the course of a numerical experiment, distributions of such quantities as flow velocity, depth, density, molecular and turbulent viscosity, values of the density of turbulent kinetic energy, dissipation of turbulent kinetic energy, and shear stress at the bottom of the flow were obtained. Using the obtained data a mathematical model is developed to describe the entrainment of the underlying material by the flow during slope erosion and the deposition of the flow material on the slope. To implement the obtained mathematical model, the architecture of the multiphaseEulerChangeFoam solver was developed, which implements a three-phase multi-velocity model with phase exchange between the material of the underlying surface and the material of the flow. The classic solver multiphaseEulerFoam from the OpenFOAM package is taken as a basis for the developed solver.

2017 ◽  
Vol 10 (12) ◽  
pp. 4511-4523 ◽  
Author(s):  
Tarandeep S. Kalra ◽  
Alfredo Aretxabaleta ◽  
Pranay Seshadri ◽  
Neil K. Ganju ◽  
Alexis Beudin

Abstract. Coastal hydrodynamics can be greatly affected by the presence of submerged aquatic vegetation. The effect of vegetation has been incorporated into the Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST) modeling system. The vegetation implementation includes the plant-induced three-dimensional drag, in-canopy wave-induced streaming, and the production of turbulent kinetic energy by the presence of vegetation. In this study, we evaluate the sensitivity of the flow and wave dynamics to vegetation parameters using Sobol' indices and a least squares polynomial approach referred to as the Effective Quadratures method. This method reduces the number of simulations needed for evaluating Sobol' indices and provides a robust, practical, and efficient approach for the parameter sensitivity analysis. The evaluation of Sobol' indices shows that kinetic energy, turbulent kinetic energy, and water level changes are affected by plant stem density, height, and, to a lesser degree, diameter. Wave dissipation is mostly dependent on the variation in plant stem density. Performing sensitivity analyses for the vegetation module in COAWST provides guidance to optimize efforts and reduce exploration of parameter space for future observational and modeling work.


1994 ◽  
Vol 98 (978) ◽  
pp. 311-318
Author(s):  
C.P. Yeung ◽  
L.C. Squire

SummaryThe three-dimensional vortex/boundary layer interaction of a type which may occur on a high-lift aerofoil has been studied. The experimental configuration simulates the trailing vortex system generated by two differentially-deflected slats which interact with an otherwise two-dimensional boundary layer developed on the wing surface under a nominal zero pressure gradient. The mean and turbulent flowfields are measured by a triple hot-wire system. The measurements show that the trailing vortex system includes the vortex sheets shed from the slats and the single vortex formed at the discontinuity between them. The single vortex moves sideways and interacts with the boundary layer as it develops downstream. During the interaction with the boundary layer, the low momentum, high turbulent-kinetic energy flow carrying negative longitudinal vorticity is entrained from the boundary layer and rolled into the vortex at the line of lateral convergence on the test surface. Likewise, at the line of lateral divergence, the high momentum, low turbulent kinetic energy flow carried by the vortex impinges on the boundary layer, suppressing the turbulent energy level and the growth of the boundary layer.


1986 ◽  
Vol 108 (4) ◽  
pp. 806-813 ◽  
Author(s):  
H. Ozoe ◽  
A. Mouri ◽  
M. Hiramitsu ◽  
S. W. Churchill ◽  
N. Lior

This paper presents a model and numerical results for turbulent natural convection in a cubical enclosure heated from below, cooled on a portion of one vertical side wall and insulated on all other surfaces. Three-dimensional balances were derived for material, energy, and the three components of momentum, as well as for the turbulent kinetic energy k and the rate of dissipation of turbulent kinetic energy ε. The constants used in the model were the same as those used by Fraikin et al. for two-dimensional convection in a channel. Illustrative transient calculations were carried out for Ra = 106 and 107 and Pr = 0.7. Both the dominant component of the vector potential and the Nusselt number were found to converge to a steady state. Isothermal lines and velocity vectors for vertical cross sections normal to the cooled wall indicated three-dimensional effects near the side walls. A top view of the velocity vectors revealed a downward spiral flow near the side walls along the cooled vertical wall. A weak spiral flow was also found along the side walls near the wall opposing the partially cooled one. The highest values of the eddy diffusivity were 2.6 and 5.8 times the molecular kinematic viscosity for Ra = 106 and 107, respectively. A coaxial double spiral movement, similar to that previously reported for laminar natural convection, was found for the time-averaged flow field. This computing scheme is expected to be applicable to other thermal boundary conditions.


2014 ◽  
Vol 760 ◽  
pp. 304-312 ◽  
Author(s):  
Farid Karimpour ◽  
Subhas K. Venayagamoorthy

AbstractIn this study, we revisit the consequence of assuming equilibrium between the rates of production ($P$) and dissipation $({\it\epsilon})$ of the turbulent kinetic energy $(k)$ in the highly anisotropic and inhomogeneous near-wall region. Analytical and dimensional arguments are made to determine the relevant scales inherent in the turbulent viscosity (${\it\nu}_{t}$) formulation of the standard $k{-}{\it\epsilon}$ model, which is one of the most widely used turbulence closure schemes. This turbulent viscosity formulation is developed by assuming equilibrium and use of the turbulent kinetic energy $(k)$ to infer the relevant velocity scale. We show that such turbulent viscosity formulations are not suitable for modelling near-wall turbulence. Furthermore, we use the turbulent viscosity $({\it\nu}_{t})$ formulation suggested by Durbin (Theor. Comput. Fluid Dyn., vol. 3, 1991, pp. 1–13) to highlight the appropriate scales that correctly capture the characteristic scales and behaviour of $P/{\it\epsilon}$ in the near-wall region. We also show that the anisotropic Reynolds stress ($\overline{u^{\prime }v^{\prime }}$) is correlated with the wall-normal, isotropic Reynolds stress ($\overline{v^{\prime 2}}$) as $-\overline{u^{\prime }v^{\prime }}=c_{{\it\mu}}^{\prime }(ST_{L})(\overline{v^{\prime 2}})$, where $S$ is the mean shear rate, $T_{L}=k/{\it\epsilon}$ is the turbulence (decay) time scale and $c_{{\it\mu}}^{\prime }$ is a universal constant. ‘A priori’ tests are performed to assess the validity of the propositions using the direct numerical simulation (DNS) data of unstratified channel flow of Hoyas & Jiménez (Phys. Fluids, vol. 18, 2006, 011702). The comparisons with the data are excellent and confirm our findings.


2013 ◽  
Vol 8 (2) ◽  
pp. 119-127 ◽  
Author(s):  
D. K. Chibwe ◽  
G. Akdogan ◽  
P. Taskinen

Abstract Typical current operation of lateral-blown Peirce-Smith converters (PSCs) has the common phenomenon of splashing and slopping due to air injection. The splashing and wave motion in these converters cause metal losses and potential production lost time due to intermittent cleaning of the converter mouth and thus reduced process throughput. Understanding of the effect of combined top and lateral blowing could possibly lead to alternative technology advancement for increased process efficiency. In this study, computational fluid dynamics (CFD) simulations of conventional common practice (lateral blowing) and combined (top and lateral blowing) in a PSC were carried out, and results of flow variables (bath velocity, turbulence kinetic energy, etc.) were compared. The two-dimensional (2-D) and three-dimensional (3-D) simulations of the three-phase system (air–matte–slag) were executed utilizing a commercial CFD numerical software code, ANSYS FLUENT 14.0. These simulations were performed employing the volume of fluid and realizable turbulence models to account for multiphase and turbulent nature of the flow, respectively. Upon completion of the simulations, the results of the models were analysed and compared by means of density contour plots, velocity vector plots, turbulent kinetic energy vector plots, average turbulent kinetic energy, turbulent intensity contour plots and average matte bulk velocity. It was found that both blowing configuration and slag layer thickness have significant effects on mixing propagation, wave formation and splashing in the PSC as the results showed wave formation and splashing significantly being reduced by employing combined top- and lateral-blowing configurations.


2010 ◽  
Vol 160-162 ◽  
pp. 1680-1684
Author(s):  
Xing Jun Gao ◽  
Qing Liu ◽  
Ping Zou ◽  
Jian Song ◽  
Ping Li

The fundamental principle of the twist drill conical grinding method was introduced. The mathematical model of the twist drill was established. Mathematical model to establish drill bit is the geometric design, manufacture, cutting analysis and modeling on the basis of the drilling process. According to the twist drill grinding principle, using Pro/E the three-dimensional modeling of the twist drill was completed, and the feature of the conical grinding method was analyzed.


Author(s):  
Hao Lin ◽  
Brian D. Storey ◽  
Juan G. Santiago

We have developed a generalized electrokinetic model suitable for the study of microchannel flows with conductivity gradients and shallow channel depths. An asymptotic analysis was performed with channel depth-to-width ratio as the smallness parameter, and three dimensional transport equations are reduced to a set to depth-averaged equations governing flow dynamics in the streamwise-spanwise plane of a shallow channel. The momentum equation uses a Darcy-Brinkman-Forchheimer type formulation, and the convective-diffusive transport of the conductivity field in the depth direction manifests itself as a dispersion effect on in-plane motion. Accuracy of the model was assessed by comparing the numerical results with direct numerical simulations. These depth-averaged equations provide the accuracy of three-dimensional modeling with a convenient quasi-two-dimensional equation set applicable to a fairly wide class of microfluidic devices.


2000 ◽  
Vol 423 ◽  
pp. 175-203 ◽  
Author(s):  
CHANDRASEKHAR KANNEPALLI ◽  
UGO PIOMELLI

A three-dimensional shear-driven turbulent boundary layer over a flat plate generated by moving a section of the wall in the transverse direction is studied using large-eddy simulations. The configuration is analogous to shear-driven boundary layer experiments on spinning cylinders, except for the absence of curvature effects. The data presented include the time-averaged mean flow, the Reynolds stresses and their budgets, and instantaneous flow visualizations. The near-wall behaviour of the flow, which was not accessible to previous experimental studies, is investigated in detail. The transverse mean velocity profile develops like a Stokes layer, only weakly coupled to the streamwise flow, and is self-similar when scaled with the transverse wall velocity, Ws. The axial skin friction and the turbulent kinetic energy, K, are significantly reduced after the imposition of the transverse shear, due to the disruption of the streaky structures and of the outer-layer vortical structures. The turbulent kinetic energy budget reveals that the decrease in production is responsible for the reduction of K. The flow then adjusts to the perturbation, reaching a quasi-equilibrium three-dimensional collateral state. Following the cessation of the transverse motion, similar phenomena take place again. The flow eventually relaxes back to a two-dimensional equilibrium boundary layer.


Sign in / Sign up

Export Citation Format

Share Document