Topology optimisation-based design of duct cross-sections for fully developed magnetohydrodynamic flows

In this communication, duct cross-sections for fully developed inductionless magnetohydrodynamics (MHD) flows are designed via topology optimisation (TO). The objective function to be maximised is the volumetric flow rate that can be achieved along the duct under a prescribed pressure gradient, subj...

Descripción completa

Detalles Bibliográficos
Autores: Álvarez Hostos, Juan Carlos|||0000-0002-4636-4948, Urgorri, Fernando Roca, Principe, Ricardo Javier|||0000-0002-1478-2651
Tipo de recurso: artículo
Fecha de publicación:2026
País:España
Institución:Universitat Politècnica de Catalunya (UPC)
Repositorio:UPCommons. Portal del coneixement obert de la UPC
Idioma:inglés
OAI Identifier:oai:upcommons.upc.edu:2117/456315
Acceso en línea:https://hdl.handle.net/2117/456315
https://dx.doi.org/10.1016/j.cma.2026.118787
Access Level:acceso abierto
Palabra clave:Magnetohydrodynamics
Topology optimisation
Inductionless
Carman-Kozeny
Developed flow
Àrees temàtiques de la UPC::Enginyeria mecànica::Mecànica de fluids
Descripción
Sumario:In this communication, duct cross-sections for fully developed inductionless magnetohydrodynamics (MHD) flows are designed via topology optimisation (TO). The objective function to be maximised is the volumetric flow rate that can be achieved along the duct under a prescribed pressure gradient, subject to an area constraint imposed on the cross-sectional geometry. The solid-liquid transition within the topology is modelled through a Darcy term, where the permeability is defined using the Carman-Kozeny equation as a function of an apparent liquid fraction, effectively enforcing zero velocity in solid regions. The well-known solid isotropic material with penalisation (SIMP) technique is used to define the apparent liquid fraction and to interpolate between the electrical conductivities of solid and fluid phases. A finite element framework is used for the multiphysics analysis, with a single design variable (artificial density) assigned to each finite element. Details of the adjoint-based sensitivity analysis within this multiphysics context are provided, and the influence of both the orientation and dominance of the magnetic field on the resulting topologies is also examined. The results demonstrate that the proposed optimisation strategy effectively yields solutions that maximise the flow rate, and that such topologies differ significantly from those expected in standard fluid dynamics problems, with strong dependence on magnetic field characteristics. Furthermore, the combined use of SIMP interpolation and the Carman-Kozeny-based Darcy term ensures stable convergence of the optimisation process within the inductionless MHD framework, even when accounting for the interpolation of electrical conductivity.