Fluid - Structure Interaction Modeling of Elastohydrodynamically Lubricated Line Contacts

2020 ◽  
pp. 1-39
Author(s):  
Kushagra Singh ◽  
Farshid Sadeghi ◽  
Thomas Russell ◽  
Steven J. Lorenz ◽  
Wyatt L. Peterson ◽  
...  

Abstract This paper presents a partitioned fluid-structure interaction (FSI) solver to model elastohydrodynamic lubrication (EHL) of line contacts. The FSI model was constructed using the multiphysics simulation software ANSYS wherein an iterative implicit coupling scheme is implemented to facilitate the interaction between fluid and solid components. The model employs a finite volume method (FVM) based computational fluid dynamics (CFD) solver to determine the lubricant flow behavior using the Navier-Stokes equations. Additionally, the finite element method (FEM) is utilized to model the structural response of the solid. Fluid cavitation, compressibility, non-Newtonian lubricant rheology, load balance algorithm and dynamic meshing were incorporated in the FSI model. The pressure and film thickness results obtained from the model are presented for a wide range of loads, speeds, slide to roll ratios (SRR), surface dent, material properties (elastic plastic), etc. The model presents a detailed understanding of EHL contacts by removing any assumptions relative to the Reynolds equation. It provides the (i) two-dimensional variation of pressure, velocity, viscosity etc. in the fluid, and (ii) stress, elastic/plastic strain in the solid, simultaneously. The FSI model is robust, easy to implement and computationally efficient. It provides an effective approach to solve sophisticated EHL problems. The FSI model was used to investigate the effects of surface dents, plasticity and material inclusions under heavily loaded lubricated line contacts as can be found in gears and rolling element bearings. The results from the model exhibit excellent corroboration with published results based on the Reynolds equation solvers.

Author(s):  
Lucia Sargentini ◽  
Benjamin Cariteau ◽  
Morena Angelucci

This paper is related to fluid-structure interaction analysis of sodium cooled fast reactors core (Na-FBR). Sudden liquid evacuation between assemblies could lead to overall core movements (flowering and compaction) causing variations of core reactivity. The comprehension of the structure behavior during the evacuation could improve the knowledge about some SCRAMs for negative reactivity occurred in PHÉNIX reactor and could contribute on the study of the dynamic behavior of a FBR core. An experimental facility (PISE-2c) is designed composed by a Poly-methyl methacrylate hexagonal rods (2D-plan similitude with PHÉNIX assembly) with a very thin gap between assemblies. Another experimental device (PISE-1a) is designed and composed by a single hexagonal rod for testing the dynamic characteristics. Different experiments are envisaged: free vibrations and oscillations during water injection. A phenomenological analysis is reported showing the flow behavior in the gap and the structure response. Also computational simulations are presented in this paper. An efficient numerical method is used to solve Navier-Stokes equations coupled with structure dynamic equation. The numerical method is verified by the comparison of analytic models and experiments.


Author(s):  
Jeremiah N. Mpagazehe ◽  
C. Fred Higgs

In many tribological applications, such as journal bearings and gears, a fluid film is used to accommodate velocity between moving surfaces. To model the behavior of this film and to predict its ability to carry load, the Reynolds equation is predominantly employed. As computational processing power continues to increase, computational fluid dynamics (CFD) is increasingly being employed to predict the fluid behavior in lubrication environments. Using CFD is advantageous in that it can provide a more general approximation to the Navier-Stokes equations than the Reynolds equation. Moreover, using CFD allows for the simulation of multiphase flows as could occur during bearing contamination and bearing exit conditions. Because the bearing surfaces move relative to each other as they obtain equilibrium with the fluid pressure, there is a need to incorporate the moving boundary into the CFD calculation, which is a non-trivial task. In this work, a fluid-structure interaction (FSI) technique is explored as an approach to model the dynamic coupling between the moving bearing surfaces and the lubricant. The benefits of using an FSI approach are discussed and the results of its implementation in a lubricated sliding contact model are presented.


2021 ◽  
Vol 2118 (1) ◽  
pp. 012003
Author(s):  
H E Calderón ◽  
L M Rada ◽  
J S De Plaza

Abstract This research focuses on improving the hydraulic behavior of a traditionally design bottom rack intake, from variations in roughness parameters, free height, and the inclusion of chamfers, establishing a contribution to the contrast between classical physics and the physics that takes over the partial resolution of the Navier-Stokes equations. To make possible the structure in OpenFOAM, it is necessary to use the geometric tool Salome-Meca, as well as a meshing tool (snappyHexMesh), and the InterFOAM solver in the processing stage. In the same way, through the turbulence model (K-E) local effects are evidenced in the Fluid-Structure interaction, as well as the identification of events and the development of the phenomenon of vorticity. The results show the improvement presented in some areas of the structure from the stabilization of the water flow through of the fluid-structure interaction change, the modification of the geometry and roughness, minimizing the presence of vertical vortices, cavitation, and surrounding areas. This allows us to conclude that traditional hydraulic do not consider the real physical flow behavior within the structure and neither the subsequent phenomena that develop, establishing as a starting point the need to rethink the design of the bottom rack intakes.


2016 ◽  
Vol 08 (08) ◽  
pp. 1650095 ◽  
Author(s):  
H. Devaraj ◽  
Kean C. Aw ◽  
E. Haemmerle ◽  
R. Sharma

3D printed hair-like micro-structures have been previously demonstrated in a novel micro-fluidic flow sensor aimed at sensing air flows down to rates of a few milliliters per second. However, there is a lack of in-depth understanding of the structural response of these ‘micro-hairs' under a fluid flow field. This paper demonstrates the use of lattice Boltzmann methods (LBM) to understand this structural response towards a better optimization of the micro-hair flow sensors designed to suit the end applications' needs. The LBM approach was chosen as an efficient alternative to simulate Navier–Stokes equations for modeling fluid flow around complex geometries primarily for improved accuracy and simplicity with lesser computational costs. As the spatial dimensions of the sensor's flow channel are much larger in comparison to the actual micro-hairs (the sensing element), a multidimensional approach of combining two-dimensional (D2Q9) and three-dimensional (D3Q19) lattice configurations were implemented for improved computational speeds and efficiency. The drag force on the micro-hairs was estimated using the momentum-exchange method in the D3Q19 configuration and this drag force is transferred to the structural analysis model which determines the micro-hair deformation using Euler–Bernoulli beam theory. The entirety of the LBM Fluid–Structure Interaction (FSI) model was implemented within MATLAB and the obtained results are compared against the numerical model implemented on a commercially available software package.


Author(s):  
M. Benaouicha ◽  
S. Guillou ◽  
A. Santa Cruz ◽  
H. Trigui

The study deals with a 3D Fluid-Structure Interaction (FSI) numerical model of a rectangular cantilevered flexible hydrofoil subjected to a turbulent fluid flow regime. The structural response and dynamic deformations are studied by analyzing the oscillations frequencies and amplitudes, under a hydrodynamics loads. The obtained numerical results are confronted with experimental ones, for validation. The numerical model is performed in the same geometric, physical and material conditions as the experimental set-up carried out in a hydrodynamic tunnel. A polyacetal (POM) flexible hydrofoil NACA0015 with an angle of attack of 8° is considered to be immersed in a fluid flow at a Reynold number of 3 × 105. The structure is initially at rest and then moved by the action of the fluid flow. The numerical model is based on a strong coupling procedure for solving the Fluid-Structure Interaction problem. The Arbitrary Lagrangian-Eulerian (ALE) formulation of the Navier-Stokes equations is used and an anisotropic diffusion equation is solved to compute the fluid mesh velocity and position at each time step. The finite volume method is used for the numerical resolution of the fluid dynamics equations. The structure deformations are described by the linear elasticity equation which is solved by the finite elements method. The Fluid-Structure coupled problem is solved by using the partitioned FSI implicit algorithm. A good agreement between numerical and experimental results for the hydrodynamics coefficients and hydrofoil deformations, maximum deflection and frequencies is obtained. The added mass and damping are analyzed and then the FSI effect on the dynamic deformations of the structure is highlighted.


Author(s):  
Fande Kong ◽  
Xiao-Chuan Cai

Fluid-structure interaction (FSI) problems are computationally very challenging. In this paper we consider the monolithic approach for solving the fully coupled FSI problem. Most existing techniques, such as multigrid methods, do not work well for the coupled system since the system consists of elliptic, parabolic and hyperbolic components all together. Other approaches based on direct solvers do not scale to large numbers of processors. In this paper, we introduce a multilevel unstructured mesh Schwarz preconditioned Newton–Krylov method for the implicitly discretized, fully coupled system of partial differential equations consisting of incompressible Navier–Stokes equations for the fluid flows and the linear elasticity equation for the structure. Several meshes are required to make the solution algorithm scalable. This includes a fine mesh to guarantee the solution accuracy, and a few isogeometric coarse meshes to speed up the convergence. Special attention is paid when constructing and partitioning the preconditioning meshes so that the communication cost is minimized when the number of processor cores is large. We show numerically that the proposed algorithm is highly scalable in terms of the number of iterations and the total compute time on a supercomputer with more than 10,000 processor cores for monolithically coupled three-dimensional FSI problems with hundreds of millions of unknowns.


Author(s):  
Yogesh Ramesh Patel

This paper provides a brief overview of the research in the field of Fluid-structure interaction in Wind Turbines. Fluid-Structure Interaction (FSI) is the interplay of some movable or deformable structure with an internal or surrounding fluid flow. Flow brought about vibrations of two airfoils used in wind turbine blades are investigated by using a strong coupled fluid shape interplay approach. The approach is based totally on a regularly occurring Computational Fluid Dynamics (CFD) code that solves the Navier-Stokes equations defined in Arbitrary Lagrangian-Eulerian (ALE) coordinates by way of a finite extent method. The need for the FSI in the wind Turbine system is studied and comprehensively presented.


2019 ◽  
Vol 142 (4) ◽  
Author(s):  
Srinivas Swaroop Kolla ◽  
Ram S. Mohan ◽  
Ovadia Shoham

Abstract The gas-liquid cylindrical cyclone (GLCC©, The University of Tulsa, 1994) is a simple, compact, and low-cost separator, which provides an economically attractive alternative to conventional gravity-based separators over a wide range of applications. The GLCC© inlet section design is a key parameter, which is crucial for its performance and proper operation. An in-depth evaluation of specific design modifications and their effect on safety and structural robustness are carried out in this study using finite element analysis (FEA). Fluid–structure interaction (FSI) analysis is also carried out using the results of computational fluid dynamics (CFD) aimed at investigating the effect of fluid flow on the inlet section structural integrity. The selected design modifications are based on feasibility of GLCC© manufacturing and assembly for field applications. Different case studies incorporating sustained GLCC© internal pressure, dead weight loading, forces generated because of slug flow and high temperatures are analyzed and presented in this paper. The concept of holes cut out in baffle has been effective with no stresses or deformation in the baffle area. FSI simulation of slug flow has proved that FEA direct loading case studies are far more conservative.


Fluids ◽  
2019 ◽  
Vol 4 (2) ◽  
pp. 94 ◽  
Author(s):  
Cornel Marius Murea

A monolithic semi-implicit method is presented for three-dimensional simulation of fluid–structure interaction problems. The updated Lagrangian framework is used for the structure modeled by linear elasticity equation and, for the fluid governed by the Navier–Stokes equations, we employ the Arbitrary Lagrangian Eulerian method. We use a global mesh for the fluid–structure domain where the fluid–structure interface is an interior boundary. The continuity of velocity at the interface is automatically satisfied by using globally continuous finite element for the velocity in the fluid–structure mesh. The method is fast because we solve only a linear system at each time step. Three-dimensional numerical tests are presented.


2008 ◽  
Vol 130 (3) ◽  
Author(s):  
C. G. Giannopapa ◽  
G. Papadakis

In the conventional approach for fluid-structure-interaction problems, the fluid and solid components are treated separately and information is exchanged across their interface. According to the conventional terminology, the current numerical methods can be grouped in two major categories: partitioned methods and monolithic methods. Both methods use separate sets of equations for fluid and solid that have different unknown variables. A unified solution method has been presented in the previous work of Giannopapa and Papadakis (2004, “A New Formulation for Solids Suitable for a Unified Solution Method for Fluid-Structure Interaction Problems,” ASME PVP 2004, San Diego, CA, July, PVP Vol. 491–1, pp. 111–117), which is different from these methods. The new approach treats both fluid and solid as a single continuum; thus, the whole computational domain is treated as one entity discretized on a single grid. Its behavior is described by a single set of equations, which are solved fully implicitly. In this paper, the elastodynamic equations are reformulated so that they contain the same unknowns as the Navier–Stokes equations, namely, velocities and pressure. Two time marching and one spatial discretization scheme, widely used for fluid equations, are applied for the solution of the reformulated equations for solids. Using linear stability analysis, the accuracy and dissipation characteristics of the resulting difference equations are examined. The aforementioned schemes are applied to a transient structural problem (beam bending) and the results compare favorably with available analytic solutions and are consistent with the conclusions of the stability analysis. A parametric investigation using different meshes, time steps, and beam dimensions is also presented. For all cases examined, the numerical solution was stable and robust and therefore is suitable for the next stage of application to full fluid-structure-interaction problems.


Sign in / Sign up

Export Citation Format

Share Document