background image
1D ANSYS® based busbar model
Like NEWBUS [1] and BUSCAL [2], it is possible to write a
FORTRAN program to solve the coupled voltage and temperature
set of non-linear equations generated from the discretisation of a
3D busbar network geometry into a network of 1D line elements.
Nowadays, it is also possible to create an Excel spreadsheet similar
to the one presented in figure 1.
But building a versatile and user friendly in-house application takes
time. For efficiency reasons, it was decided to use 1D thermo-
electric line elements available in ANSYS
®
to build a 1D line
elements network model similar to NEWBUS and BUSCAL (see
figure 6).
Figure 6. Mesh of the 1D ANSYS® based busbar model.
The model contains 108 line elements and 120 nodes. Except for the
V shaped current collectors at the base of the negative risers, each
busbar has been represented by a 1D line element of identical
section and length. The external busbar surface area needs to be
calculated as an element property and again except for the V
shaped negative current collectors those surface areas are also
identical to the one automatically computed by ANSYS
®
in the 3D
busbar model. Of course, the aluminum busbar material properties
and the boundary conditions are identical to the one set in the 3D
model.
The maximum temperature reached in the positive busbar is 98.4
ºC and in the negative busbar is 186.0 ºC. The global busbar
network voltage drop is 230 mV with 235.6 kA or 47.1% passing
in the negative side and 264.4 kA or 52.9% passing in the positive
side.
By comparing the results between the 3D and the 1D model, it
appears obvious that the negative risers section is more resistive in
the 1D model version. Clearly the 1D line elements representation
of the V shaped current collectors, which is really a 2D feature, is
not accurately enough represented.
It is of course possible to improve the 1D model representation of
the V shaped negative current conductors (now that it is known that
it is required to do so following the comparison exercise). This only
requires the construction of a more complex 1D line elements
network (see figure 7).
In this improved 1D model version, the maximum temperature
reach in the positive busbar is 95.5 ºC and in the negative busbar is
190.9 ºC. The global busbar network voltage drop is 224 mV with
241.2 kA or 48.2% passing in the negative side and 258.8 kA or
51.8% passing in the positive side.
Figure 7. Improved mesh of the 1D ANSYS® based busbar model.
This is clearly a better representation of the real 3D busbar network
but still some discrepancies remain. This comparison exercise
highlights very well that it is not that easy to perfectly represent a
real 3D busbar geometry using only a 1D line elements
discretisation. But considering that little loss of accuracy, is the 1D
ANSYS
®
based model a more efficient tool to carry out busbar
sizing optimization studies compared to the 3D ANSYS
®
based
model?
It takes of course a lot more efforts to build a 3D model geometry
than a 1D model geometry. An even if it is possible to build a
parametric 3D model, there is a limit in the versatility of a single
model topology i.e. it will not be possible to analyze all the
possible busbar configuration one might want to analyze using a
single 3D model topology. So using the 3D ANSYS
®
based model
will translate into long delays spent in building different 3D busbar
configuration model topology often to use them only a few times!
This is clearly not a very practical approach.
Unfortunately, using an ANSYS
®
based 1D line elements model is
not a very practical solution either because ANSYS
®
is not offering
a convenient way to apply the global heat transfer boundary
conditions. Each busbar external surface area must be calculated
manually by the model developer and input as a material property.
This means that without creating a macro to streamline the model
setup task, it would be very cumbersome to use that approach to
carry out a busbar sizing optimization study. So neither the 3D
ANSYS
®
based model and the 1D ANSYS
®
based model would be
the efficient and versatile tool required.
Versatile 1D busbar model part of MHD-Valdis
MHD-Valdis is a fully non-linear MHD cell stability solver that
was initially developed as an in-house code for the former USSR
aluminum industry and later adapted for the in-house use at
Reynolds Metals. But it has now become available, quite
exceptionally, as a commercial software of certain unique
advantages. Because it is a specialized MHD cell stability solver, it
has been designed to perform this single task efficiently. Being a
commercial package, it is quite user friendly and comes with
customer training and support if required.
The data input for the network current modeling and the magneto-
hydrodynamic programs are unified. The MHD package generates
the very large set of Kirchhoff equations set automatically from the
relatively simple unified data input. The computed temperatures of
the electrically heated bars permit to accommodate to the resistivity
changes. Additional flexibility permits to simulate anode changes,