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. 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
  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. Abstract We consider the interaction between a poroelastic structure, described using the Biot model in primal form, and a free-flowing fluid, modelled with the time-dependent incompressible Stokes equations. We propose a diffuse interface model in which a phase field function is used to write each integral in the weak formulation of the coupled problem on the entire domain containing both the Stokes and Biot regions. The phase field function continuously transitions from one to zero over a diffuse region of width $$\mathcal{O}(\varepsilon)$$ around the interface; this allows the weak forms to be integrated uniformly across the domain, and obviates tracking the subdomains or the interface between them. We prove convergence in weighted norms of a finite element discretization of the diffuse interface model to the continuous diffuse model; here the weight is a power of the distance to the diffuse interface. We, in turn, prove convergence of the continuous diffuse model to the standard, sharp interface, model. Numerical examples verify the proven error estimates, and illustrate application of the method to fluid flow through a complex network, describing blood circulation in the circle of Willis. 
    more » « less