skip to main content


Title: Prediction of Hydroplaning Potential Using Fully Coupled Finite Element-Computational Fluid Dynamics Tire Models
Abstract Hydroplaning is a phenomenon that occurs when a layer of water between the tire and pavement pushes the tire upward. The tire detaches from the pavement, preventing it from providing sufficient forces and moments for the vehicle to respond to driver control inputs such as breaking, accelerating, and steering. This work is mainly focused on the tire and its interaction with the pavement to address hydroplaning. Using a tire model that is validated based on results found in the literature, fluid–structure interaction (FSI) between the tire-water-road surfaces is investigated through two approaches. In the first approach, the coupled Eulerian–Lagrangian (CEL) formulation was used. The drawback associated with the CEL method is the laminar assumption and that the behavior of the fluid at length scales smaller than the smallest element size is not captured. To improve the simulation results, in the second approach, an FSI model incorporating finite element methods (FEMs) and the Navier–Stokes equations for a two-phase flow of water and air, and the shear stress transport k–ω turbulence model, was developed and validated, improving the prediction of real hydroplaning scenarios. With large computational and processing requirements, a grid dependence study was conducted for the tire simulations to minimize the mesh size yet retain numerical accuracy. The improved FSI model was applied to hydroplaning speed and cornering force scenarios.  more » « less
Award ID(s):
1650423
NSF-PAR ID:
10215338
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Journal of Fluids Engineering
Volume:
142
Issue:
10
ISSN:
0098-2202
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT Tire-pavement interaction noise (TPIN) is one of the main sources of exterior noise produced by vehicles traveling at greater than 50 kph. The dominant frequency content is typically within 500–1500 Hz. Structural tire vibrations are among the principal TPIN mechanisms. In this work, the structure of the tire is modeled and a new wave propagation solution to find its response is proposed. Multiple physical effects are accounted for in the formulation. In an effort to analyze the effects of curvature, a flat plate and a cylindrical shell model are presented. Orthotropic and nonuniform structural properties along the tire's transversal direction are included to account for differences between its sidewalls and belt. Finally, the effects of rotation and inflation pressure are also included in the formulation. Modeled frequency response functions are analyzed and validated. In addition, a new frequency-domain formulation is presented for the computation of input tread pattern contact forces. Finally, the rolling tire's normal surface velocity response is coupled with a boundary element model to demonstrate the radiated noise at the leading and trailing edge locations. These results are then compared with experimental data measured with an on-board sound intensity system. 
    more » « less
  2. ABSTRACT Being able to estimate tire/rubber friction is very important to tire engineers, materials developers, and pavement engineers. This is because of the need for estimating forces generated at the contact, optimizing tire and vehicle performance, and estimating tire wear. Efficient models for contact area and interfacial separation are key for accurate prediction of friction coefficient. Based on the contact mechanics and surface roughness, various models were developed that can predict real area of contact and penetration depth/interfacial separation. In the present work, we intend to compare the analytical contact mechanics models using experimental results and numerical analysis. Nano-indentation experiments are performed on the rubber compound to obtain penetration depth data. A finite element model of a rubber block in contact with a rough surface was developed and validated using the nano-indentation experimental data. Results for different operating conditions obtained from the developed finite element model are compared with analytical model results, and further model improvements are discussed. 
    more » « less
  3. Abstract

    A new general contact model is proposed for preventing inter‐leaflet penetration of bio‐prosthetic heart valves (BHV) at the end of the systole, which has the advantage of applying kinematic constraints directly and creating smooth free edges. At the end of each time step, the impenetrability constraints and momentum exchange between the impacting bodies are applied separately based on the coefficient of restitution. The contact method is implemented in a rotation‐free, large deformation, and thin shell finite‐element (FE) framework based on loop's subdivision surfaces. A nonlinear, anisotropic material model for a BHV is employed which uses Fung‐elastic constitutive laws for in‐plane and bending responses, respectively. The contact model is verified and validated against several benchmark problems. For a BHV‐specific validation, the computed strains on different regions of a BHV under constant pressure are compared with experimentally measured data. Finally, dynamic simulations of BHV under physiological pressure waveform are performed for symmetrical and asymmetrical fiber orientations incorporating the new contact model and compared with the penalty contact method. The proposed contact model provides the coaptation area of a functioning BHV during the closing phase for both of the fiber orientations. Our results show that fiber orientation affects the dynamic of leaflets during the opening and closing phases. A swirling motion for the BHV with asymmetrical fiber orientation is observed, similar to experimental data. To include the fluid effects, fluid–structure interaction (FSI) simulation of the BHV is performed and compared to the dynamic results.

     
    more » « less
  4. This study presents the first 3D two-way coupled fluid structure interaction (FSI) simulation of a hybrid anechoic wind tunnel (HAWT) test section with modeling all important effects, such as turbulence, Kevlar wall porosity and deflection, and reveals for the first time the complete 3D flow structure associated with a lifting model placed into a HAWT. The Kevlar deflections are captured using finite element analysis (FEA) with shell elements operated under a membrane condition. Three-dimensional RANS CFD simulations are used to resolve the flow field. Aerodynamic experimental results are available and are compared against the FSI results. Quantitatively, the pressure coefficients on the airfoil are in good agreement with experimental results. The lift coefficient was slightly underpredicted while the drag was overpredicted by the CFD simulations. The flow structure downstream of the airfoil showed good agreement with the experiments, particularly over the wind tunnel walls where the Kevlar windows interact with the flow field. A discrepancy between previous experimental observations and juncture flow-induced vortices at the ends of the airfoil is found to stem from the limited ability of turbulence models. The qualitative behavior of the flow, including airfoil pressures and cross-sectional flow structure is well captured in the CFD. From the structural side, the behavior of the Kevlar windows and the flow developing over them is closely related to the aerodynamic pressure field induced by the airfoil. The Kevlar displacement and the transpiration velocity across the material is dominated by flow blockage effects, generated aerodynamic lift, and the wake of the airfoil. The airfoil wake increases the Kevlar window displacement, which was previously not resolved by two-dimensional panel-method simulations. The static pressure distribution over the Kevlar windows is symmetrical about the tunnel mid-height, confirming a dominantly two-dimensional flow field. 
    more » « less
  5. As a step towards addressing a scarcity of references on this topic, we compared the Eulerian and Lagrangian Computational Fluid Dynamics (CFD) approaches for the solution of free-surface and Fluid–Solid Interaction (FSI) problems. The Eulerian approach uses the Finite Element Method (FEM) to spatially discretize the Navier–Stokes equations. The free surface is handled via the volume-of-fluid (VOF) and the level-set (LS) equations; an Immersed Boundary Method (IBM) in conjunction with the Nitsche’s technique were applied to resolve the fluid–solid coupling. For the Lagrangian approach, the smoothed particle hydrodynamics (SPH) method is the meshless discretization technique of choice; no additional equations are needed to handle free-surface or FSI coupling. We compared the two approaches for a flow around cylinder. The dam break test was used to gauge the performance for free-surface flows. Lastly, the two approaches were compared on two FSI problems—one with a floating rigid body dropped into the fluid and one with an elastic gate interacting with the flow. We conclude with a discussion of the robustness, ease of model setup, and versatility of the two approaches. The Eulerian and Lagrangian solvers used in this study are open-source and available in the public domain. 
    more » « less