Finite volume distance field solution applied to medial axis transform

Accurate and efficient computation of the nearest wall distance d (or level set) is important for many areas of computational science/engineering. Differential equation-based distance/ level set algorithms, such as the hyperbolic-natured Eikonal equation, have demonstrated valuable computational efficiency. Here, in the context, as an 'auxiliary' equation to the main flow equations, the Eikonal equation is solved efficiently with two different finite volume approaches (the cell vertex and cell-centered). Application of the distance solution is studied for various geometries. Moreover, a procedure using the differential field to obtain the medial axis transform (MAT) for different geometries is presented. The latter provides a skeleton representation of geometric models that has many useful analysis properties. As an alternative approach to the pure geometric methods (e.g. the Voronoi approach), the current d-MAT procedure bypasses many difficulties that are usually encountered by pure geometric methods, especially in three dimensional space. It is also shown that the d-MAT approach provides the potential to sculpt/control the MAT form for specialized solution purposes. Copyright © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Theoretical simulation studies of pulsed field magnetisation of (RE)BCO bulk superconductors

The magnetisation of bulk high temperature superconductors (HTS), such as RE-Ba-Cu-O [(RE)BCO, where RE is a rare earth element or Y], by a practical technique is essential for their application in high field, permanent magnet-like devices. Research to-date into the pulsed field magnetisation (PFM) of these materials, however, has been limited generally to experimental techniques, with relatively little progress in the development of theoretical models. This is because not only is a multi-physics approach needed to take account of the heating of the samples but also the high electric fields generated are well above the regime in which there are reliable experimental results. This paper describes a framework of theoretical simulation using the finite element method (FEM) that is applicable to both single- and multi-pulse magnetisation processes of (RE)BCO bulk superconductors. The model incorporates the heat equation and provides a convenient way of determining the distribution of trapped field, current density and temperature change within a bulk superconductor at each stage of the magnetisation process. An example of the single-pulse magnetisation of a (RE)BCO bulk is described. Potentially, the model may serve as a cost-effective tool for the optimisation of the bulk geometry and the magnetisation profile in multi-pulse magnetisation processes. © 2010 IOP Publishing Ltd.

The direct field boundary impedance of two-dimensional periodic structures with application to high frequency vibration prediction.

Large sections of many types of engineering construction can be considered to constitute a two-dimensional periodic structure, with examples ranging from an orthogonally stiffened shell to a honeycomb sandwich panel. In this paper, a method is presented for computing the boundary (or edge) impedance of a semi-infinite two-dimensional periodic structure, a quantity which is referred to as the direct field boundary impedance matrix. This terminology arises from the fact that none of the waves generated at the boundary (the direct field) are reflected back to the boundary in a semi-infinite system. The direct field impedance matrix can be used to calculate elastic wave transmission coefficients, and also to calculate the coupling loss factors (CLFs), which are required by the statistical energy analysis (SEA) approach to predicting high frequency vibration levels in built-up systems. The calculation of the relevant CLFs enables a two-dimensional periodic region of a structure to be modeled very efficiently as a single subsystem within SEA, and also within related methods, such as a recently developed hybrid approach, which couples the finite element method with SEA. The analysis is illustrated by various numerical examples involving stiffened plate structures.

Evolution of localized blobs of swirling or buoyant fluid with and without an ambient magnetic field.

We investigate the evolution of localized blobs of swirling or buoyant fluid in an infinite, inviscid, electrically conducting fluid. We consider the three cases of a strong imposed magnetic field, a weak imposed magnetic field, and no magnetic field. For a swirling blob in the absence of a magnetic field, we find, in line with others, that the blob bursts radially outward under the action of the centrifugal force, forming a thin annular vortex sheet. A simple model of this process predicts that the vortex sheet thins exponentially fast and that it moves radially outward with constant velocity. These predictions are verified by high-resolution numerical simulations. When an intense magnetic field is applied, this phenomenon is suppressed, with the energy and angular momentum of the blob now diffusing axially along the magnetic field lines, converting the blob into a columnar structure. For modest or weak magnetic fields, there are elements of both types of behavior, with the radial bursting dominating over axial diffusion for weak fields. However, even when the magnetic field is very weak, the flow structure is quite distinct to that of the nonmagnetic case. In particular, a small but finite magnetic field places a lower bound on the thickness of the annular vortex sheet and produces an annulus of counter-rotating fluid that surrounds the vortex core. The behavior of the buoyant blob is similar. In the absence of a magnetic field, it rapidly develops the mushroomlike shape of a thermal, with a thin vortex sheet at the top and sides of the mushroom. Again, a simple model of this process predicts that the vortex sheet at the top of the thermal thins exponentially fast and rises with constant velocity. These predictions are consistent with earlier numerical simulations. Curiously, however, it is shown that the net vertical momentum associated with the blob increases linearly in time, despite the fact that the vertical velocity at the front of the thermal is constant. As with the swirling blob, an imposed magnetic field inhibits the formation of a vortex sheet. A strong magnetic field completely suppresses the phenomenon, replacing it with an axial diffusion of momentum, while a weak magnetic field allows the sheet to form, but places a lower bound on its thickness. The magnetic field does not, however, change the net vertical momentum of the blob, which always increases linearly with time.

Immersed Boundary Method for generalised finite volume and finite difference Navier-Stokes solvers

In Immersed Boundary Methods (IBM) the effect of complex geometries is introduced through the forces added in the Navier-Stokes solver at the grid points in the vicinity of the immersed boundaries. Most of the methods in the literature have been used with Cartesian grids. Moreover many of the methods developed in the literature do not satisfy some basic conservation properties (the conservation of torque, for instance) on non-uniform meshes. In this paper we will follow the RKPM method originated by Liu et al. [1] to build locally regularized functions that verify a number of integral conditions. These local approximants will be used both for interpolating the velocity field and for spreading the singular force field in the framework of a pressure correction scheme for the incompressible Navier-Stokes equations. We will also demonstrate the robustness and effectiveness of the scheme through various examples. Copyright © 2010 by ASME.

Energy stable and momentum conserving hybrid finite element method for the incompressible Navier-Stokes equations

A hybrid method for the incompressible Navier-Stokes equations is presented. The method inherits the attractive stabilizing mechanism of upwinded discontinuous Galerkin methods when momentum advection becomes significant, equal-order interpolations can be used for the velocity and pressure fields, and mass can be conserved locally. Using continuous Lagrange multiplier spaces to enforce flux continuity across cell facets, the number of global degrees of freedom is the same as for a continuous Galerkin method on the same mesh. Different from our earlier investigations on the approach for the Navier-Stokes equations, the pressure field in this work is discontinuous across cell boundaries. It is shown that this leads to very good local mass conservation and, for an appropriate choice of finite element spaces, momentum conservation. Also, a new form of the momentum transport terms for the method is constructed such that global energy stability is guaranteed, even in the absence of a pointwise solenoidal velocity field. Mass conservation, momentum conservation, and global energy stability are proved for the time-continuous case and for a fully discrete scheme. The presented analysis results are supported by a range of numerical simulations. © 2012 Society for Industrial and Applied Mathematics.

3D modeling of high-T c superconductors by finite element software

A three-dimensional (3D) numerical model is proposed to solve the electromagnetic problems involving transport current and background field of a high-T c superconducting (HTS) system. The model is characterized by the E-J power law and H-formulation, and is successfully implemented using finite element software. We first discuss the model in detail, including the mesh methods, boundary conditions and computing time. To validate the 3D model, we calculate the ac loss and trapped field solution for a bulk material and compare the results with the previously verified 2D solutions and an analytical solution. We then apply our model to test some typical problems such as superconducting bulk array and twisted conductors, which cannot be tackled by the 2D models. The new 3D model could be a powerful tool for researchers and engineers to investigate problems with a greater level of complicity.

A finite element model of magnetization of superconducting bulks using a solid-state flux pump

Superconductors have a bright future; they are able to carry very high current densities, switch rapidly in electronic circuits, detect extremely small perturbations in magnetic fields, and sustain very high magnetic fields. Of most interest to large-scale electrical engineering applications are the ability to carry large currents and to provide large magnetic fields. There are many projects that use the first property, and these have concentrated on power generation, transmission, and utilization; however, there are relatively few, which are currently exploiting the ability to sustain high magnetic fields. The main reason for this is that high field wound magnets can and have been made from both BSCCO and YBCO, but currently, their cost is much higher than the alternative provided by low-Tc materials such as Nb3Sn and NbTi. An alternative form of the material is the bulk form, which can be magnetized to high fields. This paper explains the mechanism, which allows superconductors to be magnetized without the need for high field magnets to perform magnetization. A finite-element model is presented, which is based on the E-J current law. Results from this model show how magnetization of the superconductor builds up cycle upon cycle when a traveling magnetic wave is induced above the superconductor. © 2011 IEEE.

Modification of the trapped field in bulk high-temperature superconductors as a result of the drilling of a pattern of artificial columnar holes

The trapped magnetic field is examined in bulk high-temperature superconductors that are artificially drilled along their c-axis. The influence of the hole pattern on the magnetization is studied and compared by means of numerical models and Hall probe mapping techniques. To this aim, we consider two bulk YBCO samples with a rectangular cross-section that are drilled each by six holes arranged either on a rectangular lattice (sample I) or on a centered rectangular lattice (sample II). For the numerical analysis, three different models are considered for calculating the trapped flux: (i), a two-dimensional (2D) Bean model neglecting demagnetizing effects and flux creep, (ii), a 2D finite-element model neglecting demagnetizing effects but incorporating magnetic relaxation in the form of an E-J power law, and, (iii), a 3D finite element analysis that takes into account both the finite height of the sample and flux creep effects. For the experimental analysis, the trapped magnetic flux density is measured above the sample surface by Hall probe mapping performed before and after the drilling process. The maximum trapped flux density in the drilled samples is found to be smaller than that in the plain samples. The smallest magnetization drop is found for sample II, with the centered rectangular lattice. This result is confirmed by the numerical models. In each sample, the relative drops that are calculated independently with the three different models are in good agreement. As observed experimentally, the magnetization drop calculated in the sample II is the smallest one and its relative value is comparable to the measured one. By contrast, the measured magnetization drop in sample (1) is much larger than that predicted by the simulations, most likely because of a change of the microstructure during the drilling process.

Study of thermal effects in bulk RE-BCO superconductors submitted to a variable magnetic field

When bulk RE-BCO superconductors are used as permanent magnets in engineering applications, they are likely to experience transient variations of the applied magnetic field. The resulting vortex motion may cause a significant temperature increase. As a consequence the initial trapped flux is reduced. In the present work, we first focus on the cause of a temperature increase. The temperature distribution within a superconducting finite cylinder subjected to an alternating magnetic field is theoretically predicted. Results are compared to experimental data obtained by two temperature sensors attached to a bulk YBCO pellet. Second, we consider curative methods for reducing the effect of heat flux on the temperature increase. Hall-probe mappings on YBCO samples maintained out of the thermal equilibrium are performed for two different morphologies : a plain single domain and a single domain with a regularly spaced hole array. The drilled single-domain displays a trapped induction which is weakly affected by the local heating while displaying a high trapped field. © 2006 IOP Publishing Ltd.

Analysis of fields in an air-cored superconducting synchronous motor with an HTS racetrack field winding

High temperature superconducting (HTS) synchronous motors can offer significant weight and size reductions, as well as improved efficiency, over conventional copper-wound machines due to the higher current density of high temperature superconducting (HTS) materials. In order to optimise the design parameters and performance of such a machine, this paper proposes a basic physical model of an air-cored HTS synchronous motor with a copper armature winding and HTS field winding. An analytical method for the field analysis in the synchronous motor is then presented, followed by a numerical finite element analysis (FEA) model to verify the analytical solution. The model is utilised to study the influence of the geometry of the HTS coils on the magnetic field at the armature winding, and geometrical parameter optimisation is carried out using this theoretical model to obtain a more sinusoidal magnetic field at the armature, which has a major influence on the performance of the motor.

Field programmable gate array based predictive control system for spacecraft rendezvous in elliptical orbits

Copyright © 2014 John Wiley & Sons, Ltd. Copyright © 2014 John Wiley & Sons, Ltd. Summary A field programmable gate array (FPGA) based model predictive controller for two phases of spacecraft rendezvous is presented. Linear time-varying prediction models are used to accommodate elliptical orbits, and a variable prediction horizon is used to facilitate finite time completion of the longer range manoeuvres, whilst a fixed and receding prediction horizon is used for fine-grained tracking at close range. The resulting constrained optimisation problems are solved using a primal-dual interior point algorithm. The majority of the computational demand is in solving a system of simultaneous linear equations at each iteration of this algorithm. To accelerate these operations, a custom circuit is implemented, using a combination of Mathworks HDL Coder and Xilinx System Generator for DSP, and used as a peripheral to a MicroBlaze soft-core processor on the FPGA, on which the remainder of the system is implemented. Certain logic that can be hard-coded for fixed sized problems is implemented to be configurable online, in order to accommodate the varying problem sizes associated with the variable prediction horizon. The system is demonstrated in closed-loop by linking the FPGA with a simulation of the spacecraft dynamics running in Simulink on a PC, using Ethernet. Timing comparisons indicate that the custom implementation is substantially faster than pure embedded software-based interior point methods running on the same MicroBlaze and could be competitive with a pure custom hardware implementation.

Modelling and comparison of trapped fields in (RE)BCO bulk superconductors for activation using pulsed field magnetization

The ability to generate a permanent, stable magnetic field unsupported by an electromotive force is fundamental to a variety of engineering applications. Bulk high temperature superconducting (HTS) materials can trap magnetic fields of magnitude over ten times higher than the maximum field produced by conventional magnets, which is limited practically to rather less than 2 T. In this paper, two large c-axis oriented, single-grain YBCO and GdBCO bulk superconductors are magnetized by the pulsed field magnetization (PFM) technique at temperatures of 40 and 65 K and the characteristics of the resulting trapped field profile are investigated with a view of magnetizing such samples as trapped field magnets (TFMs) in situ inside a trapped flux-type superconducting electric machine. A comparison is made between the temperatures at which the pulsed magnetic field is applied and the results have strong implications for the optimum operating temperature for TFMs in trapped flux-type superconducting electric machines. The effects of inhomogeneities, which occur during the growth process of single-grain bulk superconductors, on the trapped field and maximum temperature rise in the sample are modelled numerically using a 3D finite-element model based on the H-formulation and implemented in Comsol Multiphysics 4.3a. The results agree qualitatively with the observed experimental results, in that inhomogeneities act to distort the trapped field profile and reduce the magnitude of the trapped field due to localized heating within the sample and preferential movement and pinning of flux lines around the growth section regions (GSRs) and growth sector boundaries (GSBs), respectively. The modelling framework will allow further investigation of various inhomogeneities that arise during the processing of (RE)BCO bulk superconductors, including inhomogeneous Jc distributions and the presence of current-limiting grain boundaries and cracks, and it can be used to assist optimization of processing and PFM techniques for practical bulk superconductor applications. © 2014 IOP Publishing Ltd.