A parallel monolithic algorithm for the numerical simulation of large‐scale fluid structure interaction problems

A novel parallel monolithic algorithm has been developed for the numerical simulation of large‐scale fluid structure interaction problems. The governing incompressible Navier–Stokes equations for the fluid domain are discretized using the arbitrary Lagrangian–Eulerian formulation‐based side‐centered unstructured finite volume method. The deformation of the solid domain is governed by the constitutive laws for the nonlinear Saint Venant–Kirchhoff material, and the classical Galerkin finite element method is used to discretize the governing equations in a Lagrangian frame. A special attention is given to construct an algorithm with exact total fluid volume conservation while obeying both the global and the local discrete geometric conservation law. The resulting large‐scale algebraic nonlinear equations are multiplied with an upper triangular right preconditioner that results in a scaled discrete Laplacian instead of a zero block in the original system. Then, a one‐level restricted additive Schwarz preconditioner with a block‐incomplete factorization within each partitioned sub‐domains is utilized for the modified system. The accuracy and performance of the proposed algorithm are verified for the several benchmark problems including a pressure pulse in a flexible circular tube, a flag interacting with an incompressible viscous flow, and so on. John Wiley & Sons, Ltd.     

A modular,operator‐splitting scheme for fluid–structure interaction problems with thick structures

We present an operator‐splitting scheme for fluid–structure interaction (FSI) problems in hemodynamics, where the thickness of the structural wall is comparable to the radius of the cylindrical fluid domain. The equations of linear elasticity are used to model the structure, while the Navier–Stokes equations for an incompressible viscous fluid are used to model the fluid. The operator‐splitting scheme, based on the Lie splitting, separates the elastodynamics structure problem from a fluid problem in which structure inertia is included to achieve unconditional stability. We prove energy estimates associated with unconditional stability of this modular scheme for the full nonlinear FSI problem defined on a moving domain, without requiring any sub‐iterations within time steps. Two numerical examples are presented, showing excellent agreement with the results of monolithic schemes. First‐order convergence in time is shown numerically. Modularity, unconditional stability without temporal sub‐iterations, and simple implementation are the features that make this operator‐splitting scheme particularly appealing for multi‐physics problems involving FSI. Copyright © 2013 John Wiley & Sons, Ltd.     

Magnetohydrodynamic and heat transfer effects on the stagnation-point flow of an electrically conducting nanofluid past a porous vertical shrinking/stretching sheet in the presence of variable stream conditions

A steady two-dimensional magnetohydrodynamic stagnation-point flow of an electrically conducting fluid and heat transfer with thermal radiation of a nanofluid past a shrinking and stretching sheet is investigated numerically. The model used for the nanofluid incorporates the effects of the Brownian motion and thermophoresis. A similarity transformation is used to convert the governing nonlinear boundary-layer equations into coupled higher-order nonlinear ordinary differential equations. The result shows that the velocity, temperature, and concentration profiles are significantly influenced by the Brownian motion, heat radiation, and thermophoresis particle deposition.     

Reduced order modelling method via proper orthogonal decomposition (POD) for flow around an oscillating cylinder

This paper presents reduced order modelling (ROM) in fluid–structure interaction (FSI). The ROM via the proper orthogonal decomposition (POD) method has been chosen, due to its efficiency in the domain of fluid mechanics. POD-ROM is based on a low-order dynamical system obtained by projecting the nonlinear Navier–Stokes equations on a smaller number of POD modes. These POD modes are spatial and temporally independent. In FSI, the fluid and structure domains are moving, owing to which the POD method cannot be applied directly to reduce the equations of each domain. This article proposes to compute the POD modes for a global velocity field (fluid and solid), and then to construct a low-order dynamical system. The structure domain can be decomposed as a rigid domain, with a finite number of degrees of freedom. This low-order dynamical system is obtained by using a multiphase method similar to the fictitious domain method. This multiphase method extends the Navier–Stokes equations to the solid domain by using a penalisation method and a Lagrangian multiplier. By projecting these equations on the POD modes obtained for the global velocity field, a nonlinear low-order dynamical system is obtained and tested on a case of high Reynolds number.     

Numerical studies for flow and heat transfer of the Powell-Eyring fluid thin film over an unsteady stretching sheet with internal heat generation using the chebyshev finite difference method

An analysis is carried out to study the unsteady two-dimensional Powell-Eyring flow and heat transfer to a laminar liquid film from a horizontal stretching surface in the presence of internal heat generation. The flow of a thin fluid film and subsequent heat transfer from the stretching surface is investigated with the aid of a similarity transformation. The transformation enables to reduce the unsteady boundary layer equations to a system of nonlinear ordinary differential equations. A numerical solution of the resulting nonlinear differential equations is found by using an efficient Chebyshev finite difference method. A comparison of numerical results is made with the earlier published results for limiting cases. The effects of the governing parameters on the flow and thermal fields are thoroughly examined and discussed.     

Chemically reactive and radiative von Kármán swirling flow due to a rotating disk

A new mathematical model is presented to study the heat and mass transfer characteristics of magnetohydrodynamic (MHD) Maxwell fluid flow over a convectively heated stretchable rotating disk. To regulate the fluid temperature at the surface, a simple isothermal model of homogeneous-heterogeneous reactions is employed. The impact of nonlinear thermal radiative heat flux on thermal transport features is studied. The transformed nonlinear system of ordinary differential equations is solved numerically with an efficient method, namely, the Runge-Kutta-Felberg fourth-order and fifth-order (RKF45) integration scheme using the MAPLE software. Achieved results are validated with previous studies in an excellent way. Major outcomes reveal that the magnetic flux reduces the velocity components in the radial, angular, and axial directions, and enhances the fluid temperature. Also, the presence of radiative heat flux is to raise the temperature of fluid. Further, the strength of homogeneous–heterogeneous reactions is useful to diminish the concentration of reaction.     

Rotational flow of Oldroyd-B nanofluid subject to Cattaneo-Christov double diffusion theory

A nanofluid is composed of a base fluid component and nanoparticles, in which the nanoparticles are dispersed in the base fluid. The addition of nanoparticles into a base fluid can remarkably improve the thermal conductivity of the nanofluid, and such an increment of thermal conductivity can play an important role in improving the heat transfer rate of the base fluid. Further, the dynamics of non-Newtonian fluids along with nanoparticles is quite interesting with numerous industrial applications. The present predominately predictive modeling studies the flow of the viscoelastic Oldroyd-B fluid over a rotating disk in the presence of nanoparticles. A progressive amendment in the heat and concentration equations is made by exploiting the Cattaneo-Christov heat and mass flux expressions. The characteristic of the Lorentz force due to the magnetic field applied normal to the disk is studied. The Buongiorno model together with the Cattaneo-Christov theory is implemented in the Oldroyd-B nanofluid flow to investigate the heat and mass transport mechanism. This theory predicts the characteristics of the fluid thermal and solutal relaxation time on the boundary layer flow. The von K′arm′an similarity functions are utilized to convert the partial differential equations(PDEs) into ordinary differential equations(ODEs). A homotopic approach for obtaining the analytical solutions to the governing nonlinear problem is carried out. The graphical results are obtained for the velocity field, temperature, and concentration distributions. Comparisons are made for a limiting case between the numerical and analytical solutions, and the results are found in good agreement. The results reveal that the thermal and solutal relaxation time parameters diminish the temperature and concentration distributions, respectively. The axial flow decreases in the downward direction for higher values of the retardation time parameter. The impact of the thermophoresis parameter boosts the temperature distribution.     

Jeffrey fluid flow due to curved stretching surface with Cattaneo-Christov heat flux

The two-dimensional (2D) motion of the Jeffrey fluid by the curved stretching sheet coiled in a circle is investigated. The non-Fourier heat flux model is used for the heat transfer analysis. Feasible similarity variables are used to transform the highly nonlinear ordinary equations to partial differential equations (PDEs). The homotopy technique is used for the convergence of the velocity and temperature equations. The effects of the involved parameters on the physical properties of the fluid are described graphically. The results show that the curvature parameter is an increasing function of velocity and temperature, and the temperature is a decreasing function of the thermal relaxation time. Besides, the Deborah number has a reverse effect on the pressure and surface drag force.     

属性与温度相关材料的广义热弹性问题有限元分析

计及材料特性与温度的相关性，基于Lord和Shulman(L-S)广义热弹性理 论，建立了此类问题的有限元控制方程. 由于材料属性的温度相关性，温度控制方程具有非 线性，积分变换求解方法难以采用，因而将有限元方程直接在时间域求解. 利用所建立方法 研究了材料特性与温度相关、带有孔洞的无限大体在热冲击和机械冲击作用下的广义热弹性 问题. 分析表明，在时间域直接求解材料属性与温度相关的广义热弹性问题是可行的，所得 结果具有很高的精度，热的波动性得到充分的展现. 同时发现，热冲击载荷作用时，材料属 性与温度的相关性对结构的机械响应影响显著，对温度响应影响很小；机械载荷作用时，材 料参数与温度的相关性对所有响应影响都很小. 因此，研究热冲击载荷作用的机械响应时， 必须考虑材料属性的温度相关性，而研究温度响应时，无论何种冲击载荷，都可以不考虑材 料属性的温度相关性.     

Upshot of ohmically dissipated Darcy-Forchheimer slip flow of magnetohydrodynamic Sutterby fluid over radiating linearly stretched surface in view of Cash and Carp method

The present work concerns the momentum and heat transmission of the electro-magnetohydrodynamic(E-MHD) boundary layer Darcy-Forchheimer flow of a Sutterby fluid over a linear stretching sheet with slip. The nonlinear equations for the proposed model are analyzed numerically. Suitable techniques are used to transform the coupled nonlinear partial differential equations(PDEs) conforming to the forced balance law, energy, and concentration equations into a nonlinear coupled system of ordinary differential equations(ODEs). Numerical solutions of the transformed nonlinear system are obtained using a shooting method, improved by the Cash and Carp coefficients. The influence of important physical variables on the velocity, the temperature, the heat flux coefficient, and the skin-friction coefficient is verified and analyzed through graphs and tables. From the comprehensive analysis of the present work, it is concluded that by intensifying the magnitude of the Hartmann number, the momentum distribution decays,whereas the thermal profile of fluid increases. Furthermore, it is also shown that by augmenting the values of the momentum slip parameter, the velocity profile diminishes. It is found that the Sutterby fluid model shows shear thickening and shear thinning behaviors.The momentum profile shows that the magnitude of velocity for the shear thickening case is dominant as compared with the shear thinning case. It is also demonstrated that the Sutterby fluid model reduces to a Newtonian model by fixing the fluid parameter to zero.In view of the limiting case, it is established that the surface drag in the case of the Sutterby model shows a trifling pattern as compared with the classical case.     

Effects of temperature dependent fluid properties and variable Prandtl number on the transient convective flow due to a porous rotating disk

In this paper we have studied the effects of temperature dependent fluid properties such as density, viscosity and thermal conductivity and variable Prandtl number on unsteady convective heat transfer flow over a porous rotating disk. Using similarity transformations we reduce the governing nonlinear partial differential equations for flow and heat transfer into a system of ordinary differential equations which are then solved numerically by applying Nachtsheim–Swigert shooting iteration technique along with sixth-order Runge–Kutta integration scheme. Comparison with previously published work for steady case of the problem were performed and found to be in very good agreement. The obtained numerical results show that the rate of heat transfer in a fluid of constant properties is higher than in a fluid of variable properties. The results further show that consideration of Prandtl number as constant within the boundary layer for variable fluid properties lead unrealistic results. Therefore, modeling thermal boundary layers with temperature dependent fluid properties Prandtl number must treated as variable inside the boundary layer.     

Fractional four-step finite element method for analysis of thermally coupled fluid-solid interaction problems

An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin (SUPG) method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid are performed by the Galerkin method. The second-order semiimplicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are linearized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effiectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.     

A comparison of two formulations of the fin efficiency for straight fins

A formulation of the fin efficiency based on Newton’s law of cooling is compared with a formulation based on a ratio of heat transferred from the fin surface to the surrounding fluid to the heat conducted through the base.The first formulation requires that the solution of the nonlinear fin equations for constant heat transfer coefficient and constant thermal conductivity is known,whilst the second formulation of the fin efficiency requires only that a first integral of the model equation is known.This paper shows the first formulation of the fin efficiency contains approximation errors as only power series and approximate solutions to the nonlinear fin equations have been determined.The second formulation of the fin efficiency is exact when the first integrals can be determined.     

Effect of Viscous Dissipation on the Boundary Layer Flow and Heat Transfer Past a Permeable Stretching Surface Embedded in a Porous Medium with a Second-Order Slip Using Chebyshev Finite Difference Method

This article investigates a theoretical and numerical study for the effect of viscous dissipation on the steady flow with heat transfer of Newtonian fluid toward a permeable stretching surface embedded in a porous medium with a second-order slip and thermal slip. The governing nonlinear partial differential equations are converted into nonlinear ordinary differential equations (ODEs) using similarity variables. The resulting ODEs are successfully solved numerically with the help of Chebyshev finite difference method. Graphically results are shown for non-dimensional velocities and temperature. The effects of the porous parameter, the suction (injection) parameter, Eckert number, first- and second-order velocity slip parameter, the thermal slip parameter and the Prandtl number on the flow and temperature profiles are presented. Moreover, the local skin-friction and Nusselt numbers are presented. A comparison of numerical results is made with the earlier published results under limiting cases.     

On the viscous dissipation modeling of thermal fluid flow in a porous medium

The problem of viscous dissipation and thermal dispersion in saturated porous medium is numerically investigated for the case of non-Darcy flow regime. The fluid is induced to flow upward by natural convection as a result of a semi-infinite vertical wall that is immersed in the porous medium and is kept at constant higher temperature. The boundary layer approximations were used to simplify the set of the governing, nonlinear partial differential equations, which were then non-dimensionalized and solved using the finite elements method. The results for the details of the governing parameters are presented and investigated. It is found that the irreversible process of transforming the kinetic energy of the moving fluid to heat energy via the viscosity of the moving fluid (i.e., viscous dissipation) resulted in insignificant generation of heat for the range of parameters considered in this study. On the other hand, thermal dispersion has shown to disperse heat energy normal to the wall more effectively compared with the normal diffusion mechanism.     

Series solutions of annular axisymmetric stagnation flow and heat transfer on moving cylinder

The steady, laminar, incompressible flow and heat transfer of a viscous fluid between two circular cylinders for two different types of thermal boundary conditions are investigated. The governing Navier-Stokes and thermal equations of the flow are reduced to a nonlinear system of ordinary differential equations. The equations are solved analyt- ically using the homotopy analysis method （HAM）. Convergence of the HAM solutions is discussed in detail. These solutions are then compared with recently obtained numericM and perturbative solutions. Plots of the velocity and temperature profiles are provided for various values of the relevant parameters.     

Thermal radiation effect on a mixed convection flow and heat transfer of the Williamson fluid past an exponentially shrinking permeable sheet with a convective boundary condition

The thermal radiation effect on a steady mixed convective flow with heat transfer of a nonlinear (non-Newtonian) Williamson fluid past an exponentially shrinking porous sheet with a convective boundary condition is investigated numerically. In this study, both an assisting flow and an opposing flow are considered. The governing equations are converted into nonlinear ordinary differential equations by using a suitable transformation. A numerical solution of the problem is obtained by using the Matlab software package for different values of the governing parameters. The results show that dual nonsimilar solutions exist for the opposing flow, whereas the solution for the assisting flow is unique. It is also observed that the dual nonsimilar solutions exist only if a certain amount of mass suction is applied through the porous sheet, which depends on the Williamson parameter, convective parameter, and radiation parameter.     

Development of a hybrid model based on mesh and meshfree methods and its application to fluid–elastic structure interaction for free surface waves

In this paper, a hybrid scheme, Fluid–Fluid–Elastic Structure (FFES) model was developed in the time domain to address the wave breaking impact on the structure. The model is developed based on the partitioned approach with different governing equations that describe various regions of the model domain. The fluid–fluid model denotes that two different fluid models were used to describe fluid in the actual physical domain. The method is a physics-based approximation to reduce the computational time, i.e. in the far-field inviscid fluid (fully nonlinear potential flow theory model), and near to the structure, viscous fluid (Navier Stokes model) is used. The coupled model then interacts with the elastic structure (based on Euler–Bernoulli beam theory). The system of equations is strongly coupled both in space and time. The Fluid–Fluid coupling uses an implicit predictor–corrector scheme, and the fluid–structure coupling works based on an iterative scheme. This approach makes the method more robust and for future extension. Three different possibilities for introducing the coupling was identified and implemented. The model was validated against results from the analytical solution and literature. The method proposed is a reliable, robust, and efficient alternative for simulating fluid–structure interaction problems.     

Free Convection Boundary Layer Flow Past a Vertical Surface in a Porous Medium with Temperature-Dependent Properties

Numerical solutions for the free convection heat transfer in a viscous fluid at a permeable surface embedded in a saturated porous medium, in the presence of viscous dissipation with temperature-dependent variable fluid properties, are obtained. The governing equations for the problem are derived using the Darcy model and the Boussinesq approximation (with nonlinear density temperature variation in the buoyancy force term). The coupled non-linearities arising from the temperature-dependent density, viscosity, thermal conductivity, and viscous dissipation are included. The partial differential equations of the model are reduced to ordinary differential equations by a similarity transformation and the resulting coupled, nonlinear ordinary differential equations are solved numerically by a second order finite difference scheme for several sets of values of the parameters. Also, asymptotic results are obtained for large values of | f _w|. Moreover, the numerical results for the velocity, the temperature, and the wall-temperature gradient are presented through graphs and tables, and are discussed. It is observed that by increasing the fluid variable viscosity parameter, one could reduce the velocity and thermal boundary layer thickness. However, quite the opposite is true with the non-linear density temperature variation parameter.     

Viscoelastic MHD flow boundary layer over a stretching surface with viscous and ohmic dissipations

In this study the momentum and heat transfer characteristics in an incompressible electrically conducting viscoelastic boundary layer fluid flow over a linear stretching sheet are considered. Highly non-linear momentum and thermal boundary layer equations are reduced to set of nonlinear ordinary differential equations by appropriate transformation.