Effects of Cilia Movement on Fluid Velocity: II Numerical Solutions Over a Fixed Domain

2020 ◽  
Vol 134 (2) ◽  
pp. 471-489
Author(s):  
Kanognudge Wuttanachamsri ◽  
Lynn Schreyer
Keyword(s):  
Numerical Solutions ◽  
Fluid Velocity ◽  
Fixed Domain

On steady compressible flows with compact vorticity; the compressible Hill's spherical vortex

1998 ◽  
Vol 374 ◽  
pp. 285-303 ◽  
Author(s):  
D. W. MOORE ◽  
D. I. PULLIN
Keyword(s):  
Finite Difference ◽  
Compressible Flows ◽  
Numerical Solutions ◽  
Flow Direction ◽  
Fluid Velocity ◽  
Major Axis ◽  
Sonic Point ◽  
Euler Flow ◽  
Compressible Euler ◽  
Spherical Vortex

We consider steady compressible Euler flow corresponding to the compressible analogue of the well-known incompressible Hill's spherical vortex (HSV). We first derive appropriate compressible Euler equations for steady homentropic flow and show how these may be used to define a continuation of the HSV to finite Mach number M∞=U∞/C∞, where U∞, C∞ are the fluid velocity and speed of sound at infinity respectively. This is referred to as the compressible Hill's spherical vortex (CHSV). It corresponds to axisymmetric compressible Euler flow in which, within a vortical bubble, the azimuthal vorticity divided by the product of the density and the distance to the axis remains constant along streamlines, with irrotational flow outside the bubble. The equations are first solved numerically using a fourth-order finite-difference method, and then using a Rayleigh–Janzen expansion in powers of M2∞ to order M4∞. When M∞>0, the vortical bubble is no longer spherical and its detailed shape must be determined by matching conditions consisting of continuity of the fluid velocity at the bubble boundary. For subsonic compressible flow the bubble boundary takes an approximately prolate spheroidal shape with major axis aligned along the flow direction. There is good agreement between the perturbation solution and Richardson extrapolation of the finite difference solutions for the bubble boundary shape up to M∞ equal to 0.5. The numerical solutions indicate that the flow first becomes locally sonic near or at the bubble centre when M∞≈0.598 and a singularity appears to form at the sonic point. We were unable to find shock-free steady CHSVs containing regions of locally supersonic flow and their existence for the present continuation of the HSV remains an open question.

Electrohydrodynamic rotation of drops at large electric Reynolds numbers

2016 ◽  
Vol 788 ◽  
Author(s):  
Ehud Yariv ◽  
Itzchak Frankel
Keyword(s):  
Electric Fields ◽  
Numerical Solutions ◽  
Fluid Velocity ◽  
Reynolds Numbers ◽  
Mirror Image ◽  
Liquid Drops ◽  
Weakly Conducting ◽  
Strong Electric Fields ◽  
Quincke Rotation ◽  
Asymmetric Solution

When subject to sufficiently strong electric fields, particles and drops suspended in a weakly conducting liquid exhibit spontaneous rotary motion. This so-called Quincke rotation is a fascinating example of nonlinear symmetry-breaking phenomena. To illuminate the rotation of liquid drops we here analyse the asymptotic limit of large electric Reynolds numbers, $\mathit{Re}\gg 1$, within the framework of a two-dimensional Taylor–Melcher electrohydrodynamic model. A non-trivial dominant balance in this singular limit results in both the fluid velocity and surface-charge density scaling as $\mathit{Re}^{-1/2}$. The flow is governed by a self-contained nonlinear boundary-value problem that does not admit a continuous fore–aft symmetric solution, thus necessitating drop rotation. Furthermore, thermodynamic arguments reveal that a fore–aft asymmetric solution exists only when charge relaxation within the suspending liquid is faster than that in the drop. The flow problem possesses both mirror-image (with respect to the direction of the external field) and flow-reversal symmetries; it is transformed into a universal one, independent of the ratios of electric conductivities and dielectric permittivities in the respective drop phase and suspending liquid phase. The rescaled angular velocity is found to depend weakly upon the viscosity ratio. The corresponding numerical solutions of the exact equations indeed collapse at large $\mathit{Re}$ upon the asymptotically calculated universal solution.

Numerical analysis of time-dependent stagnation point flow of Oldroyd-B fluid subject to modified Fourier’s law

2021 ◽  
pp. 2150187
Author(s):  
Foukeea Qasim ◽  
Tian-Chuan Sun ◽  
S. Z. Abbas ◽  
W. A. Khan ◽  
M. Y. Malik
Keyword(s):  
Differential Equation ◽  
Fluid Flow ◽  
Stagnation Point ◽  
Numerical Solutions ◽  
Fluid Velocity ◽  
Nonlinear Partial Differential Equation ◽  
Thermal Relaxation ◽  
Stagnation Point Flow ◽  
Time Dependent ◽  
Parameter Values

This paper aims to investigate the time-dependent stagnation point flow of an Oldroyd-B fluid subjected to the modified Fourier law. The flow into a vertically stretched cylinder at the stagnation point is discussed. The heat flux model of a non-Fourier is intended for the transfer of thermal energy in fluid flow. The study is carried out on the surface heating source, namely the surface temperature. The developed nonlinear partial differential equation for regulating fluid flow and heat transport is transformed via appropriate similarity variables into a nonlinear ordinary differential equation. The development and analysis of convergent series solutions were considered for velocity and temperature. Prandtl number numerical values are computed and investigated. This study’s findings are compared to the previous findings. By making use of the bvp4c Matlab method, numerical solutions are obtained. Besides, high buoyancy parameter values are found to increase the fluid velocity for the stimulating approach. By improving the thermal relaxation time parameter values, heat transfer in the fluid flow decreases. The temperature field effects are displayed graphically.

scholarly journals Mathematical Model for the Electrokinetic Instability of Electrolyte Flow

2019 ◽  
Vol 224 ◽  
pp. 02003
Author(s):  
Andrey Shobukhov
Keyword(s):  
Electric Potential ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Fluid Velocity ◽  
Steady State Solution ◽  
Stability Loss ◽  
Navier Stokes ◽  
Ion Distribution ◽  
Navier Stokes Equations ◽  
Dilute Aqueous Solution

We study a one-dimensional model of the dilute aqueous solution of KCl in the electric field. Our model is based on a set of Nernst-Planck-Poisson equations and includes the incompressible fluid velocity as a parameter. We demonstrate instability of the linear electric potential variation for the uniform ion distribution and compare analytical results with numerical solutions. The developed model successfully describes the stability loss of the steady state solution and demonstrates the emerging of spatially non-uniform distribution of the electric potential. However, this model should be generalized by accounting for the convective movement via the addition of the Navier-Stokes equations in order to substantially extend its application field.

Numerical simulation of two-dimensional compressible convection

1975 ◽  
Vol 70 (4) ◽  
pp. 689-703 ◽  
Author(s):  
Eric Graham
Keyword(s):  
Rayleigh Number ◽  
Numerical Solutions ◽  
Fluid Velocity ◽  
Critical Rayleigh Number ◽  
Finite Amplitude ◽  
Two Dimensions ◽  
Onset Of Convection ◽  
Local Sound ◽  
Constant Viscosity ◽  
Convection Problem

A procedure for obtaining numerical solutions to the equations describing thermal convection in a compressible fluid is outlined. The method is applied to the case of a perfect gas with constant viscosity and thermal conductivity. The fluid is considered to be confined in a rectangular region by fixed slippery boundaries and motions are restricted to two dimensions. The upper and lower boundaries are maintained at fixed temperatures and the side boundaries are thermally insulating. The resulting convection problem can be characterized by six dimension-less parameters. The onset of convection has been studied both by obtaining solutions to the nonlinear equations in the neighbourhood of the critical Rayleigh number Rc and by solving the linear stability problem. Solutions have been obtained for values of the Rayleigh number up to 100Rc and for pressure variations of a factor of 300 within the fluid. In some cases the fluid velocity is comparable to the local sound speed. The Nusselt number increases with decreasing Prandtl number for moderate values of the depth parameter. Steady finite amplitude solutions have been found in all the cases considered. As the horizontal dimension A of the rectangle is increased, the length of time needed to reach a steady state also increases. For large values of A the solution consists of a number of rolls. Even for small values of A, no solutions have been found where one roll is vertically above another.

Fluid Flow and Lateral Fluid Dispersion in Bounded Granular Beds

2002 ◽  
Author(s):  
Azita Soleymani ◽  
Eveliina Takasuo ◽  
Piroz Zamankhan ◽  
William Polashenski
Keyword(s):  
Reynolds Number ◽  
Porous Media ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Numerical Study ◽  
Fluid Velocity ◽  
Three Dimensional ◽  
Numerical Results ◽  
Transition To Turbulence ◽  
Wide Range

Results are presented from a numerical study examining the flow of a viscous, incompressible fluid through random packing of nonoverlapping spheres at moderate Reynolds numbers (based on pore permeability and interstitial fluid velocity), spanning a wide range of flow conditions for porous media. By using a laminar model including inertial terms and assuming rough walls, numerical solutions of the Navier-Stokes equations in three-dimensional porous packed beds resulted in dimensionless pressure drops in excellent agreement with those reported in a previous study (Fand et al., 1987). This observation suggests that no transition to turbulence could occur in the range of Reynolds number studied. For flows in the Forchheimer regime, numerical results are presented of the lateral dispersivity of solute continuously injected into a three-dimensional bounded granular bed at moderate Peclet numbers. Lateral fluid dispersion coefficients are calculated by comparing the concentration profiles obtained from numerical and analytical methods. Comparing the present numerical results with data available in the literature, no evidence has been found to support the speculations by others for a transition from laminar to turbulent regimes in porous media at a critical Reynolds number.

scholarly journals Finite Element Framework for Computational Fluid Dynamics in FEBio

2018 ◽  
Vol 140 (2) ◽  
Author(s):  
Gerard A. Ateshian ◽  
Jay J. Shim ◽  
Steve A. Maas ◽  
Jeffrey A. Weiss
Keyword(s):  
Fluid Dynamics ◽  
Computational Fluid Dynamics ◽  
Finite Element ◽  
Solid Mechanics ◽  
Numerical Solutions ◽  
Virtual Work ◽  
Biological Fluids ◽  
Fluid Velocity ◽  
Benchmark Problems ◽  
Future Extension

The mechanics of biological fluids is an important topic in biomechanics, often requiring the use of computational tools to analyze problems with realistic geometries and material properties. This study describes the formulation and implementation of a finite element framework for computational fluid dynamics (CFD) in FEBio, a free software designed to meet the computational needs of the biomechanics and biophysics communities. This formulation models nearly incompressible flow with a compressible isothermal formulation that uses a physically realistic value for the fluid bulk modulus. It employs fluid velocity and dilatation as essential variables: The virtual work integral enforces the balance of linear momentum and the kinematic constraint between fluid velocity and dilatation, while fluid density varies with dilatation as prescribed by the axiom of mass balance. Using this approach, equal-order interpolations may be used for both essential variables over each element, contrary to traditional mixed formulations that must explicitly satisfy the inf-sup condition. The formulation accommodates Newtonian and non-Newtonian viscous responses as well as inviscid fluids. The efficiency of numerical solutions is enhanced using Broyden's quasi-Newton method. The results of finite element simulations were verified using well-documented benchmark problems as well as comparisons with other free and commercial codes. These analyses demonstrated that the novel formulation introduced in FEBio could successfully reproduce the results of other codes. The analogy between this CFD formulation and standard finite element formulations for solid mechanics makes it suitable for future extension to fluid–structure interactions (FSIs).

Role of heatlines on thermal management during Rayleigh-Bénard heating within enclosures with concave/convex horizontal walls

2017 ◽  
Vol 27 (9) ◽  
pp. 2070-2104 ◽  
Author(s):  
Pratibha Biswal ◽  
Tanmay Basak
Keyword(s):  
Heat Flow ◽  
Flow Visualization ◽  
Numerical Solutions ◽  
Fluid Velocity ◽  
Content Type ◽  
Bénard Convection ◽  
Benard Convection ◽  
First Time ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

Purpose This study aims to carry out the analysis of Rayleigh-Bénard convection within enclosures with curved isothermal walls, with the special implication on the heat flow visualization via the heatline approach. Design/methodology/approach The Galerkin finite element method has been used to obtain the numerical solutions in terms of the streamlines (ψ ), heatlines (Π), isotherms (θ), local and average Nusselt number ( Nut¯) for various Rayleigh numbers (103 ≤ Ra ≥ 105), Prandtl numbers (Pr = 0.015 and 7.2) and wall curvatures (concavity/convexity). Findings The presence of the larger fluid velocity within the curved cavities resulted in the larger heat transfer rates and thermal mixing compared to the square cavity. Case 3 (high concavity) exhibits the largest Nut¯ at the low Ra for all Pr. At the high Ra, Nut¯ is the largest for Case 3 (high concavity) at Pr = 0.015, whereas at Pr = 7.2, Nut¯ is the largest for Case 1 (high concavity and convexity). Practical implications The results may be useful for the material processing applications. Originality/value The study of Rayleigh-Bénard convection in cavities with the curved isothermal walls is not carried out till date. The heatline approach is used for the heat flow visualization during Rayleigh-Benard convection within the curved walled enclosures for the first time. Also, the existence of the enhanced fluid and heat circulation cells within the curved walled cavities during Rayleigh-Benard heating is illustrated for the first time.

scholarly journals Accurate Solutions to Water Wave Scattering by Vertical Thin Porous Barriers

2015 ◽  
Vol 2015 ◽  
pp. 1-11 ◽  
Author(s):  
Ai-jun Li ◽  
Yong Liu ◽  
Hua-jun Li
Keyword(s):  
Water Wave ◽  
Wave Scattering ◽  
Numerical Solutions ◽  
Fluid Velocity ◽  
Wave Interaction ◽  
Wave Forces ◽  
Pressure Jump ◽  
Water Wave Scattering ◽  
Porous Barrier ◽  
Porous Barriers

The water wave scattering by vertical thin porous barriers is accurately solved in this study. Two typical structures of a surface-piercing barrier and a submerged bottom-standing barrier are considered. The solution procedure is based on the multi-term Galerkin method, in which the pressure jump across a porous barrier is expanded in a set of basis functions involving the Chebychev polynomials. Then, the square-root singularity of fluid velocity at the edge of the porous barrier is correctly modeled. The present solutions have the merits of very rapid convergence. Accurate results for both the reflection and the transmission coefficients and wave forces are presented. This study not only gives a promising procedure to tackle wave interaction with vertical thin porous barriers but also provides a reliable benchmark for complicated numerical solutions.

On the estimation of sound produced by complex fluid–structure interactions, with application to a vortex interacting with a shrouded rotor

1991 ◽  
Vol 433 (1889) ◽  
pp. 573-598 ◽  
Keyword(s):  
Numerical Solutions ◽  
Stokes Equations ◽  
Fluid Velocity ◽  
Leading Edge ◽  
Surface Forces ◽  
Fluid Structure Interactions ◽  
Fluid Structure ◽  
Frequency Sound ◽  
Complex Fluid ◽  
Shrouded Rotor

An approximate analytical method is described for determining the sound produced by a class of complex fluid-structure interactions in low Mach number flows. This can be used to model noise sources in practical systems, and to check the accuracy of predictions based on time accurate numerical solutions of the Navier-Stokes equations. The dominant acoustic sources are dipoles whose strengths are dependent on the unsteady surface forces, and are expressed in terms of fluid velocity and vorticity, and a set of harmonic functions determined by the shapes of the structural elements interacting with the flow. (The theory of surface forces for arbitrary motion of a rigid body in viscous, incompressible flow in the presence of a fixed system of boundaries is discussed in an appendix.) These elements can influence both the generation and propagation of sound, and are frequently sources of new vorticity shed into wakes. The procedure is illustrated by application to a model problem in which sound is generated by a vortex interacting with a shrouded rotor in a duct. High-frequency sound is generated when the vortex is draw n into the rotor disc and ‘chopped’ by the blades. Sound is also produced through indirect blade-vortex interactions which, in this case, occur as a result of unsteady blade loadings produced when the core of the vortex is close to the leading edge of the shroud. This is relatively low-frequency sound and is the only component of blade-vortex interaction noise when the vortex is convected through the gap between the shroud and wall of the duct.

Sign in / Sign up

Export Citation Format

Share Document