skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on May 1, 2026

Title: Influence of magnetic field penetration in the hydromagnetic flow past a sphere
Abstract In this paper, viscous incompressible hydromagnetic flow around a sphere has been investigated by considering the penetration of the magnetic field inside it. Earlier researchers have found it difficult because it not only adds an extra equation to the system of governing equations but also needs a proper matching of the components of the magnetic fields at the interface. However, using a higher‐order compact finite difference scheme, we have successfully solved the governing highly nonlinear and coupled system of partial differential equations and have obtained converged solutions throughout the domain of the parameter space. In this novel numerical investigation, we have calculated the magnetic field throughout the whole domain, that is, both inside the sphere and within the fluid, with a suitable matching at the interface–a feature that has allowed us to capture the actual interactions occurring between the fluid flow and the magnetic field and unfurl several new characteristics of scientific and technological value. In fact, we have found that the magnetic field penetrating inside the sphere can effectively cause the critical value of Reynolds number to increase and can help to suppress flow separation more effectively than otherwise.  more » « less
Award ID(s):
2307563
PAR ID:
10644154
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
ZAMM
Date Published:
Journal Name:
ZAMM - Journal of Applied Mathematics and Mechanics / Zeitschrift für Angewandte Mathematik und Mechanik
Volume:
105
Issue:
5
ISSN:
0044-2267
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    The development of the governing equations for fluid flow in a surface-following coordinate system is essential to investigate the fluid flow near an interface deformed by propagating waves. In this paper, the governing equations of fluid flow, including conservation of mass, momentum and energy balance, are derived in an orthogonal curvilinear coordinate system relevant to surface water waves. All equations are further decomposed to extract mean, wave-induced and turbulent components. The complete transformed equations include explicit extra geometric terms. For example, turbulent stress and production terms include the effects of coordinate curvature on the structure of fluid flow. Furthermore, the governing equations of motion were further simplified by considering the flow over periodic quasi-linear surface waves wherein the wavelength of the disturbance is large compared to the wave amplitude. The quasi-linear analysis is employed to express the boundary layer equations in the orthogonal wave-following curvilinear coordinates with the corresponding decomposed equations for the mean, wave and turbulent fields. Finally, the vorticity equations are also derived in the orthogonal curvilinear coordinates in order to express the corresponding velocity–vorticity formulations. The equations developed in this paper proved to be useful in the analysis and interpretation of experimental data of fluid flow over wind-generated surface waves. Experimental results are presented in a companion paper. 
    more » « less
  2. Abstract We consider the finite element approximation of a coupled fluid‐structure interaction (FSI) system, which comprises a three‐dimensional (3D) Stokes flow and a two‐dimensional (2D) fourth‐order Euler–Bernoulli or Kirchhoff plate. The interaction of these parabolic and hyperbolic partial differential equations (PDE) occurs at the boundary interface which is assumed to be fixed. The vertical displacement of the plate dynamics evolves on the flat portion of the boundary where the coupling conditions are implemented via the matching velocities of the plate and fluid flow, as well as the Dirichlet boundary trace of the pressure. This pressure term also acts as a coupling agent, since it appears as a forcing term on the flat, elastic plate domain. Our main focus in this work is to generate some numerical results concerning the approximate solutions to the FSI model. For this, we propose a numerical algorithm that sequentially solves the fluid and plate subsystems through an effective decoupling approach. Numerical results of test problems are presented to illustrate the performance of the proposed method. 
    more » « less
  3. The coupled equations governing whistler waves propagating along a duct with cylindrical cross section are derived and then solved numerically. These equations are expressed in terms of magnetic and current flux functions and show that it is possible to have a solution where the waves are finite in the duct and decay exponentially outside the duct. This solution has the property of having zero radial Poynting flux everywhere, so, as required for whistler waves to bounce back and forth losslessly between magnetically conjugate terrestrial hemispheres, no wave power leaks from the duct. The coupled equations are solved numerically for a tangible realistic situation by dividing the radial domain into an inner and an outer region, where the interface between these regions is at a mode conversion location, where fast and slow modes inside the duct merge and effectively reflect. The result of this effective reflection is that there are fast and slow standing waves in the duct. In the region external to the duct, the wave solutions are also a form of standing waves, but with a strong exponential decay and a radial wavelength that is intermediate between that of the fast and slow waves in the duct. The numerical solution is shown to be in good quantitative agreement with estimates made from analytic models. Detailed examination of the solutions in the vicinity of the mode conversion location shows that the classic plane wave assumption fails to describe the true nature of the modes. 
    more » « less
  4. In this paper we study the self-induced low-Reynolds-number flow generated by a cylinder immersed in a stratified fluid. In the low Péclet limit, where the Péclet number is the ratio of the radius of the cylinder and the Phillips length scale, the flow is captured by a set of linear equations obtained by linearising the governing equations with respect to the prescribed far field conditions. We specifically focus on the low Péclet regime and develop a Green's function approach to solve the linearised equations governing the flow over the cylinder. We cross check our analytical solution against numerical solution of the nonlinear equations to obtain the range of the Péclet numbers for which the linear solution is valid. We then take advantage of the analytical solution to find explicit far-field decay rates of the flow. Our detailed analysis points out that the streamfunction and the velocity field decays algebraically in the far field. Intriguingly, this algebraic decay of the flow is much slower when compared with the exponential decay of the flow generated by a slow moving cylinder in the homogeneous Stokes regime, in the absence of stratification. Consequently, the flow generated by a cylinder in the stratified Stokes regime will have a larger domain of influence when compared with the flow generated by a cylinder in the homogeneous Stokes regime. 
    more » « less
  5. In this paper we report on the numerical analysis of convection patterns - due to changes in humidity - of air inside a two-dimensional square enclosure. The enclosure, with dimensions of 10 cm by 10 cm, and atmospheric air as the working fluid, is placed in a horizontal position with the gravitational force acting directly downward on it. Thermally, the system and its boundaries are at a constant temperature of 20°C, whereas the humidity varies with position inside the cavity. At the top and bottom walls, the relative humidity is set at 0 and 1, respectively, while the vertical walls are considered as impermeable. The mathematical model is based on two-dimensional versions of the conservation equations for mass, momentum and moisture, in Cartesian coordinates, under laminar flow and steady-state conditions. The governing equations were discretized and solved over the computational domain with the Finite Element method for different system conditions. The results, given in terms of velocity, density and humidity fields, show that convection patterns form as a result of buoyancy forces generated by humidity gradients, just as they do in thermal convection. Further comparison to the thermal convection process of dry air on the same system, show that the two are closely related. 
    more » « less