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

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.

2016 ◽  
Vol 788 ◽  
Author(s):  
Ehud Yariv ◽  
Itzchak Frankel

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.


2012 ◽  
Vol 226-228 ◽  
pp. 1829-1834 ◽  
Author(s):  
Jing Yuan Tang ◽  
Jian Ming Chen ◽  
Hong Bin Ma ◽  
Guang Yu Tang

The flow field characteristics in U-typed bend has been extensively studied for transit time ultrasonic flowmeters designing, but for the flowmeter with three-Z-shaped round pipe there is still lack of corresponding research. This paper presents a computational fluid dynamics (CFD) approach for modeling of the three-Z-shaped ultrasonic flowmeter and studying of internal fluid field characteristics based on Reynolds stress model (RSM). The fluid velocity profile in the three ultrasound path is obtained using CFD and secondary flow in cross section also is analyzed. The simulation results show that the internal flow fields in the flowmeter are not fully developed turbulence with asymmetric axial velocity distribution and dramatic changes along the flow direction, and there are obvious secondary cross flows on theirs cross-sections. The CFD simulations provide useful insights into the flow field associated with ultrasonic flowmeters design.


Author(s):  
Foukeea Qasim ◽  
Tian-Chuan Sun ◽  
S. Z. Abbas ◽  
W. A. Khan ◽  
M. Y. Malik

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.


2021 ◽  
Vol 63 ◽  
pp. 420-433
Author(s):  
Benjamin J. Maldon ◽  
Bishnu Lamichhane ◽  
Ngamta Thamwattana

Dye-sensitized solar cells consistently provide a cost-effective avenue for sources of renewable energy, primarily due to their unique utilization of nanoporous semiconductors. Through mathematical modelling, we are able to uncover insights into electron transport to optimize the operating efficiency of the dye-sensitized solar cells. In particular, fractional diffusion equations create a link between electron density and porosity of the nanoporous semiconductors. We numerically solve a fractional diffusion model using a finite-difference method and a finite-element method to discretize space and an implicit finite-difference method to discretize time. Finally, we calculate the accuracy of each method by evaluating the numerical errors under grid refinement. doi:10.1017/S1446181121000353


2019 ◽  
Vol 224 ◽  
pp. 02003
Author(s):  
Andrey Shobukhov

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.


Geophysics ◽  
2014 ◽  
Vol 79 (3) ◽  
pp. T125-T141 ◽  
Author(s):  
Josep de la Puente ◽  
Miguel Ferrer ◽  
Mauricio Hanzich ◽  
José E. Castillo ◽  
José M. Cela

Finite-difference methods for modeling seismic waves are known to be inaccurate when including a realistic topography, due to the large dispersion errors that appear in the modelled surface waves and the scattering introduced by the staircase approximation to the topography. As a consequence, alternatives to finite-difference methods have been proposed to circumvent these issues. We present a new numerical scheme for 3D elastic wave propagation in the presence of strong topography. This finite-difference scheme is based upon a staggered grid of the Lebedev type, or fully staggered grid (FSG). It uses a grid deformation strategy to make a regular Cartesian grid conform to a topographic surface. In addition, the scheme uses a mimetic approach to accurately solve the free-surface condition and hence allows for a less restrictive grid spacing criterion in the computations. The scheme can use high-order operators for the spatial derivatives and obtain low-dispersion results with as few as six points per minimum wavelength. A series of tests in 2D and 3D scenarios, in which our results are compared to analytical and numerical solutions obtained with other numerical approaches, validate the accuracy of our scheme. The resulting FSG mimetic scheme allows for accurate and efficient seismic wave modelling in the presence of very rough topographies with the advantage of using a structured staggered grid.


Author(s):  
Yufeng Xu ◽  
Om Agrawal

AbstractIn this paper, numerical solutions of Burgers equation defined by using a new Generalized Time-Fractional Derivative (GTFD) are discussed. The numerical scheme uses a finite difference method. The new GTFD is defined using a scale function and a weight function. Many existing fractional derivatives are the special cases of it. A linear recurrence relationship for the numerical solutions of the resulting system of linear equations is found via finite difference approach. Burgers equations with different fractional orders and coefficients are computed which show that this numerical method is simple and effective, and is capable of solving the Burgers equation accurately for a wide range of viscosity values. Furthermore, we study the influence of the scale and the weight functions on the diffusion process of Burgers equation. Numerical simulations illustrate that a scale function can stretch or contract the diffusion on the time domain, while a weight function can change the decay velocity of the diffusion process.


Author(s):  
Richard H. J. Willden ◽  
J. Michael R. Graham

The paper reports results from two strip theory CFD investigations of the Vortex-Induced Vibrations of model riser pipes. The first investigation is concerned with the vibrations of a vertical riser pipe that was subjected to a stepped current profile. An axial spatial resolution study was conducted to determine the number of simulation planes required to achieve tolerably converged numerical solutions. It was found that six to seven simulation planes are required per half-wavelength of pipe vibration in order to obtain convergence. The second investigation is concerned with the simultaneous in-plane and out-of-plane vibrations of a model Steel Catenary Riser that was subjected to a uniform current profile. The pipe’s simulated vibrations were found to agree very well with those determined experimentally. This result was achieved despite the questionable usage of simulation planes at high angles to the flow direction.


Sign in / Sign up

Export Citation Format

Share Document