The Modified LS-STAG Method Application for Planar Viscoelastic Flow Computation in a 4:1 Contraction Channel

Author(s):  
I.K. Marchevsky ◽  
V.V. Puzikova

In this study we present the modification of the LS-STAG immersed boundary cut-cell method. This modification is designed for viscoelastic fluids. Linear and quasilinear viscoelastic fluid models of a rate type are considered. The obtained numerical method is implemented in the LS-STAG software package developed by the author. This software is created for viscous incompressible flows simulation both by the LS-STAG method and by it developed modifications. Besides of this, the software package is designed to compute extra-stresses for viscoelastic Maxwell, Jeffreys, upper-convected Maxwell, Maxwell-A, Oldroyd-B, Oldroyd-A, Johnson --- Segalman fluids on the LS-STAG mesh. The construction of convective derivatives discrete analogues is described for Oldroyd, Cotter --- Rivlin, Jaumann --- Zaremba --- Noll derivatives. The centers of base LS-STAG mesh cells are the locations for shear non-Newtonian stresses computation. The corners of these cells are the positions for normal non-Newtonian stresses computation. The first order predictor--corrector scheme is the basis for time-stepping numerical algorithm. Benchmark solutions for the planar flow of Oldroyd-B fluid in a 4:1 contraction channel are presented. A critical value of Weissenberg number is defined. Computational results are in good agreement with the data known in the literature

Author(s):  
Lin Sun ◽  
Sanjay R. Mathur ◽  
Jayathi Y. Murthy

A numerical method is developed for solving the 3D, unsteady, incompressible flows with immersed moving solids of arbitrary geometrical complexity. A co-located (non-staggered) finite volume method is employed to solve the Navier-Stokes governing equations for flow region using arbitrary convex polyhedral meshes. The solid region is represented by a set of material points with known position and velocity. Faces in the flow region located in the immediate vicinity of the solid body are marked as immersed boundary (IB) faces. At every instant in time, the influence of the body on the flow is accounted for by reconstructing implicitly the velocity the IB faces from a stencil of fluid cells and solid material points. Specific numerical issues related to the non-staggered formulation are addressed, including the specification of face mass fluxes, and corrections to the continuity equation to ensure overall mass balance. Incorporation of this immersed boundary technique within the framework of the SIMPLE algorithm is described. Canonical test cases of laminar flow around stationary and moving spheres and cylinders are used to verify the implementation. Mesh convergence tests are carried out. The simulation results are shown to agree well with experiments for the case of micro-cantilevers vibrating in a viscous fluid.


Entropy ◽  
2019 ◽  
Vol 21 (12) ◽  
pp. 1219 ◽  
Author(s):  
Mark Dostalík ◽  
Vít Průša ◽  
Karel Tůma

Using a Lyapunov type functional constructed on the basis of thermodynamical arguments, we investigate the finite amplitude stability of internal steady flows of viscoelastic fluids described by the Giesekus model. Using the functional, we derive bounds on the Reynolds and the Weissenberg number that guarantee the unconditional asymptotic stability of the corresponding steady internal flow, wherein the distance between the steady flow field and the perturbed flow field is measured with the help of the Bures–Wasserstein distance between positive definite matrices. The application of the theoretical results is documented in the finite amplitude stability analysis of Taylor–Couette flow.


2020 ◽  
Vol 8 (10) ◽  
pp. 744
Author(s):  
Yu-Hsien Lin ◽  
Xian-Chen Li

A computational fluid dynamics (CFD)-based simulation using a finite volume code for a full-appendage DARPA (Defense Advanced Research Projects Agency) SUBOFF model was investigated with a sliding mesh model in a multi-zone fluid domain. Unsteady Reynolds Averaged Navier–Stokes (URANS) equations were coupled with a Menter’s shear stress transport (SST) k-ω turbulence closure based on the Boussinesq approximation. In order to simulate unsteady motions and capture unsteady interactions, the sliding mesh model was employed to simulate flows in the fluid domain that contains multiple moving zones. The pressure-based solver, semi-implicit method for the pressure linked equations-consistent (SIMPLEC) algorithm was employed for incompressible flows based on the predictor-corrector approach in a segregated manner. After the grid independence test, the numerical simulation was validated by comparison with the published experimental data and other numerical results. In this study, the capability of the CFD simulation with the sliding mesh model was well demonstrated to conduct the straight-line towing tests by analyzing hydrodynamic characteristics, viz. resistance, vorticity, frictional coefficients, and pressure coefficients.


2017 ◽  
Vol 817 ◽  
pp. 388-405 ◽  
Author(s):  
Qiang Yang ◽  
Lisa Fauci

We study the dynamics and transport of an elastic fibre in a polymeric cellular flow. The macroscopic fibre is much larger than the infinitesimal immersed polymer coils distributed in the surrounding viscoelastic fluid. Here we consider low-Reynolds-number flow using the Navier–Stokes/Fene-P equations in a two-dimensional, doubly periodic domain. The macroscopic fibre supports both tensile and bending forces, and is fully coupled to the viscoelastic fluid using an immersed boundary framework. We examine the effects of fibre flexibility and polymeric relaxation times on fibre buckling and transport as well as the evolution of polymer stress. Non-dimensional control parameters include the Reynolds number, the Weissenberg number, and the elasto-viscous number of the macroscopic fibre. We find that large polymer stresses occur in the fluid near the ends of the fibre when it is compressed. In addition, we find that viscoelasticity hinders a fibre’s ability to traverse multiple cells in the domain.


Author(s):  
Chloé Mimeau ◽  
Iraj Mortazavi ◽  
Georges-Henri Cottet

In this work, a hybrid particle-penalization technique is proposed to achieve accurate and efficient computations of 3D incompressible flows past bluff bodies. This immersed boundary approach indeed maintains the efficiency and the robustness of vortex methods and allows to easily model complex media, like solid-fluid-porous ones, without prescribing any boundary condition. In this paper, the method is applied to implement porous coatings on a hemisphere in order to passively control the flow dynamics.


Author(s):  
Byeong Rog Shin ◽  
Satoru Yamamoto ◽  
Xin Yuan

A preconditioned numerical method for gas-liquid two-phase flows is applied to solve cavitating flow. The present method employs a finite-difference method of dual time-stepping integration procedure and Roe’s flux difference splitting approximation with MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. The present density based numerical method permits simple treatment of the whole gas-liquid two-phase flow field including wave propagation, large density changes and incompressible flows characteristics at low Mach number. By this method, two-dimensional internal flows through a backward-facing step duct, a venturi tube and decelerating cascades are computed. Comparisons of predicted results with experiments are provided and discussed.


Sign in / Sign up

Export Citation Format

Share Document