scholarly journals Numerical Simulation of Micromixing of Particles and Fluids with Galloping Cylinder

Symmetry ◽  
2020 ◽  
Vol 12 (4) ◽  
pp. 580 ◽  
Author(s):  
Zahra Abdelmalek ◽  
Mohammad Yaghoub Abdollahzadeh Jamalabadi

Micromixers are significant segments inside miniaturized scale biomedical frameworks. Numerical investigation of the effects of galloping cylinder characteristics inside a microchannel Newtonian, incompressible fluid in nonstationary condition is performed. Governing equations of the system include the continuity equation, and Navier–Stokes equations are solved within a moving mesh domain. The symmetry of laminar entering the channel is broken by the self-sustained motion of the cylinder. A parameter study on the amplitude and frequency of passive moving cylinder on the mixing of tiny particles in the fluid is performed. The results show a significant increase to the index of mixing uses of the galloping body in biomedical frameworks in the course of micro-electromechanical systems (MEMS) devices.

1997 ◽  
Vol 52 (4) ◽  
pp. 358-368 ◽  
Author(s):  
Michio Nishida ◽  
Masashi Matsumotob

Abstract • This paper describes a computational study of the thermal and chemical nonequilibrium occuring in a rapidly expanding flow of high-temperature air transported as a free jet from an orifice into low-density stationary air. Translational, rotational, vibrational and electron temperatures are treated separately, and in particular the vibrational temperatures are individually treated; a multi-vibrational temperature model is adopted. The governing equations are axisymmetric Navier-Stokes equations coupled with species vibrational energy, electron energy and species mass conservation equations. These equations are numerically solved, using the second order upwind TVD scheme of the Harten-Yee type. The calculations were carried out for two different orifice temperatures and also two different orifice diameters to investigate the effects of such parameters on the structure of a nonequilibrium free jet.


1999 ◽  
Vol 387 ◽  
pp. 227-254 ◽  
Author(s):  
VALOD NOSHADI ◽  
WILHELM SCHNEIDER

Plane and axisymmetric (radial), horizontal laminar jet flows, produced by natural convection on a horizontal finite plate acting as a heat dipole, are considered at large distances from the plate. It is shown that physically acceptable self-similar solutions of the boundary-layer equations, which include buoyancy effects, exist in certain Prandtl-number regimes, i.e. 0.5<Pr[les ]1.470588 for plane, and Pr>1 for axisymmetric flow. In the plane flow case, the eigenvalues of the self-similar solutions are independent of the Prandtl number and can be determined from a momentum balance, whereas in the axisymmetric case the eigenvalues depend on the Prandtl number and are to be determined as part of the solution of the eigenvalue problem. For Prandtl numbers equal to, or smaller than, the lower limiting values of 0.5 and 1 for plane and axisymmetric flow, respectively, the far flow field is a non-buoyant jet, for which self-similar solutions of the boundary-layer equations are also provided. Furthermore it is shown that self-similar solutions of the full Navier–Stokes equations for axisymmetric flow, with the velocity varying as 1/r, exist for arbitrary values of the Prandtl number.Comparisons with finite-element solutions of the full Navier–Stokes equations show that the self-similar boundary-layer solutions are asymptotically approached as the plate Grashof number tends to infinity, whereas the self-similar solution to the full Navier–Stokes equations is applicable, for a given value of the Prandtl number, only to one particular, finite value of the Grashof number.In the Appendices second-order boundary-layer solutions are given, and uniformly valid composite expansions are constructed; asymptotic expansions for large values of the lateral coordinate are performed to study the decay of the self-similar boundary-layer flows; and the stability of the jets is investigated using transient numerical solutions of the Navier–Stokes equations.


Author(s):  
Wolfgang Höhn

During the design of the compressor and turbine stages of today’s aeroengines, aerodynamically induced vibrations become increasingly important since higher blade load and better efficiency are desired. In this paper the development of a method based on the unsteady, compressible Navier-Stokes equations in two dimensions is described in order to study the physics of flutter for unsteady viscous flow around cascaded vibrating blades at stall. The governing equations are solved by a finite difference technique in boundary fitted coordinates. The numerical scheme uses the Advection Upstream Splitting Method to discretize the convective terms and central differences discretizing the viscous terms of the fully non-linear Navier-Stokes equations on a moving H-type mesh. The unsteady governing equations are explicitly and implicitly marched in time in a time-accurate way using a four stage Runge-Kutta scheme on a parallel computer or an implicit scheme of the Beam-Warming type on a single processor. Turbulence is modelled using the Baldwin-Lomax turbulence model. The blade flutter phenomenon is simulated by imposing a harmonic motion on the blade, which consists of harmonic body translation in two directions and a rotation, allowing an interblade phase angle between neighboring blades. Non-reflecting boundary conditions are used for the unsteady analysis at inlet and outlet of the computational domain. The computations are performed on multiple blade passages in order to account for nonlinear effects. A subsonic massively stalled unsteady flow case in a compressor cascade is studied. The results, compared with experiments and the predictions of other researchers, show reasonable agreement for inviscid and viscous flow cases for the investigated flow situations with respect to the Steady and unsteady pressure distribution on the blade in separated flow areas as well as the aeroelastic damping. The results show the applicability of the scheme for stalled flow around cascaded blades. As expected the viscous and inviscid computations show different results in regions where viscous effects are important, i.e. in separated flow areas. In particular, different predictions for inviscid and viscous flow for the aerodynamic damping for the investigated flow cases are found.


It is shown that the boundary layer approximation to the flow of a viscous fluid past a flat plate of length l , generally valid near the plate when the Reynolds number Re is large, fails within a distance O( lRe -3/4 ) of the trailing edge. The appropriate governing equations in this neighbourhood are the full Navier- Stokes equations. On the basis of Imai (1966) these equations are linearized with respect to a uniform shear and are then completely solved by means of a Wiener-Hopf integral equation. The solution so obtained joins smoothly on to that of the boundary layer for a flat plate upstream of the trailing edge and for a wake downstream of the trailing edge. The contribution to the drag coefficient is found to be O ( Re -3/4 ) and the multiplicative constant is explicitly worked out for the linearized equations.


2019 ◽  
Vol 141 (12) ◽  
Author(s):  
Filipe S. Pereira ◽  
Luís Eça ◽  
Guilherme Vaz

The importance of the turbulence closure to the modeling accuracy of the partially-averaged Navier–Stokes equations (PANS) is investigated in prediction of the flow around a circular cylinder at Reynolds number of 3900. A series of PANS calculations at various degrees of physical resolution is conducted using three Reynolds-averaged Navier–Stokes equations (RANS)-based closures: the standard, shear-stress transport (SST), and turbulent/nonturbulent (TNT) k–ω models. The latter is proposed in this work. The results illustrate the dependence of PANS on the closure. At coarse physical resolutions, a narrower range of scales is resolved so that the influence of the closure on the simulations accuracy increases significantly. Among all closures, PANS–TNT achieves the lowest comparison errors. The reduced sensitivity of this closure to freestream turbulence quantities and the absence of auxiliary functions from its governing equations are certainly contributing to this result. It is demonstrated that the use of partial turbulence quantities in such auxiliary functions calibrated for total turbulent (RANS) quantities affects their behavior. On the other hand, the successive increase of physical resolution reduces the relevance of the closure, causing the convergence of the three models toward the same solution. This outcome is achieved once the physical resolution and closure guarantee the precise replication of the spatial development of the key coherent structures of the flow.


1989 ◽  
Vol 56 (1) ◽  
pp. 47-50 ◽  
Author(s):  
C. Y. Wang

Melting of a disk is facilitated by rotation. The problem is governed by a nondimensional parameter α which represents the relative importance of injection (melt) rate and rotation times viscosity. The nonlinear governing equations are solved by perturbations for small α and numerical integration for arbitrary α. Torque and heat transfer rates are found. The solution is one of the rare exact similarity solutions of the Navier-Stokes equations.


2005 ◽  
Author(s):  
Cenk Evren Ku¨krer ◽  
I˙lker Tarı

Numerical Simulations of liquid flow in a micro-channel between two horizontal plates are performed. The channel is infinite in streamwise and spanwise directions and its height is taken as 3.1×10−4 m which falls within the dimension ranges of micro-channels. The Navier-Stokes equations with the addition of Brinkman number (Br) to the energy equation are used as the governing equations and a spectral methods based approach is applied to obtain the required accuracy to handle liquid flow in the micro-channel. It is known for micro-channels that Br combines the effects of conduction and viscous dissipation in liquids and is also a way of comparing the importance of later relative to former. A laminar flow of a liquid in a micro-channel shows different characteristics compared to a similar flow in a macro-channel. To observe the differences, three different cases are run over each of a range of Reynolds numbers: one with no axial conduction assumption that correspond to a case similar to macro-channel flow, another case with axial conduction included in the energy equation to simulate one of the main differences and lastly a case with inclusion of Br number in the governing equations. The results are compared with each other to see the effects of axial conduction and Br inclusion. A qualitative comparison is made with the previous results in literature.


2020 ◽  
Vol 221 (2) ◽  
pp. 1264-1280
Author(s):  
Rene Gassmöller ◽  
Juliane Dannberg ◽  
Wolfgang Bangerth ◽  
Timo Heister ◽  
Robert Myhill

SUMMARY Mantle convection and long-term lithosphere dynamics in the Earth and other planets can be treated as the slow deformation of a highly viscous fluid, and as such can be described using the compressible Navier–Stokes equations. Since on Earth-sized planets the influence of compressibility is not a dominant effect, density deviations from a reference profile are at most on the order of a few percent and using the full governing equations poses numerical challenges, most modelling studies have simplified the governing equations. Common approximations assume a temporally constant, but depth-dependent reference profile for the density (the anelastic liquid approximation), or drop compressibility altogether and use a constant reference density (the Boussinesq approximation). In most previous studies of mantle convection and crustal dynamics, one can assume that the error introduced by these approximations was small compared to the errors that resulted from poorly constrained material behaviour and limited numerical accuracy. However, as model parametrizations have become more realistic, and model resolution has improved, this may no longer be the case and the error due to using simplified conservation equations might no longer be negligible: while such approximations may be reasonable for models of mantle plumes or slabs traversing the whole mantle, they may be unsatisfactory for layered materials experiencing phase transitions or materials undergoing significant heating or cooling. For example, at boundary layers or close to dynamically changing density gradients, the error arising from the use of the aforementioned compressibility approximations can be the dominant error source, and common approximations may fail to capture the physical behaviour of interest. In this paper, we discuss new formulations of the continuity equation that include dynamic density variations due to temperature, pressure and composition without using a reference profile for the density. We quantify the improvement in accuracy relative to existing formulations in a number of benchmark models and evaluate for which practical applications these effects are important. Finally, we consider numerical aspects of the new formulations. We implement and test these formulations in the freely available community software aspect, and use this code for our numerical experiments.


2021 ◽  
Vol 158 (A3) ◽  
Author(s):  
A Lavrov ◽  
C Guedes Soares

The laminar flow around heaving axisymmetric and three-dimensional cylinders with damping plates is numerically studied for various Keulegan-Carpenter numbers. The Navier-Stokes equations are solved using OpenFOAM, which is applied to the flow on a moving mesh. For processing of results the semi-empirical Morison equation is used. Calculations are conducted for one cylinder, one cylinder with one disk, one cylinder with two disks, and one cylinder with one pentagonal plate. The calculated values are compared against experimental data.


Sign in / Sign up

Export Citation Format

Share Document