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.

Estimating material elasticity by spherical indentation load-relaxation tests on viscoelastic samples of finite thickness

A two-step viscoelastic spherical indentation method is proposed to compensate for 1) material relaxation and 2) sample thickness. In the first step, the indenter is moved at a constant speed and the reaction force is measured. In the second step, the indenter is held at a constant position and the relaxation response of the material is measured. Then the relaxation response is fit with a multi-exponential function which corresponds to a three-branch general Maxwell model. The relaxation modulus is derived by correcting the finite ramp time introduced in the first step. The proposed model takes into account the sample thickness, which is important for applications in which the sample thickness is less than ten times the indenter radius. The model is validated numerically by finite element simulations. Experiments are carried out on a 10% gelatin phantom and a chicken breast sample with the proposed method. The results for both the gelatin phantom and the chicken breast sample agree with the results obtained from a surface wave method. Both the finite element simulations and experimental results show improved elasticity estimations by incorporating the sample thickness into the model. The measured shear elasticities of the 10% gelatin sample are 6.79 and 6.93 kPa by the proposed finite indentation method at sample thickness of 40 and 20 mm, respectively. The elasticity of the same sample is estimated to be 6.53 kPa by the surface wave method. For the chicken breast sample, the shear elasticity is measured to be 4.51 and 5.17 kPa by the proposed indentation method at sample thickness of 40 and 20 mm, respectively. Its elasticity is measured by the surface wave method to be 4.14 kPa. © 2011 IEEE.

A FEniCS-based programming framework for modeling turbulent flow by the Reynolds-averaged Navier-Stokes equations

Finding an appropriate turbulence model for a given flow case usually calls for extensive experimentation with both models and numerical solution methods. This work presents the design and implementation of a flexible, programmable software framework for assisting with numerical experiments in computational turbulence. The framework targets Reynolds-averaged Navier-Stokes models, discretized by finite element methods. The novel implementation makes use of Python and the FEniCS package, the combination of which leads to compact and reusable code, where model- and solver-specific code resemble closely the mathematical formulation of equations and algorithms. The presented ideas and programming techniques are also applicable to other fields that involve systems of nonlinear partial differential equations. We demonstrate the framework in two applications and investigate the impact of various linearizations on the convergence properties of nonlinear solvers for a Reynolds-averaged Navier-Stokes model. © 2011 Elsevier Ltd.

Building flexible user interfaces for solving PDEs

FEniCS is a collection of software tools for the automated solution of differential equations by finite element methods. In this note, we describe how FEniCS can be used to solve a simple nonlinear model problem with varying levels of automation. At one extreme, FEniCS provides tools for the fully automated and adaptive solution of nonlinear partial differential equations. At the other extreme, FEniCS provides a range of tools that allow the computational scientist to experiment with novel solution algorithms. © 2010 American Institute of Physics.

Study of iron saturation in Brushless Doubly-Fed induction machines

This paper presents an analytical modelling approach for the Brushless Doubly-Fed Machine (BDFM) taking iron saturation into account. A generalised coupled-circuit model is developed which considers stator and rotor teeth saturation effects. A method of calculating the machine inductance parameters is presented which can be implemented in time-stepping simulations. The model has been implemented in MATLAB/Simulink and verified by Finite Element analysis and experimental tests. The tests are carried out on a 180 frame size BDFM. Flux search coils have been utilised to measure airgap and teeth flux densities. © 2010 IEEE.

Experimental and numerical investigation of fracture mechanics of bitumen beams

This work presents an experimental and numerical investigation to characterise the fracture properties of pure bitumen (the binder in asphalt paving materials). The paper is divided into two parts. The first part describes an experimental study of fracture characterisation parameters of pure bitumen as determined by three-point bend tests. The second part deals with modelling of fracture and failure of bitumen by Finite Element analysis. Fracture mechanics parameters, stress intensity factor, KIC, fracture energy, GIC, and J-integral, JIC, are used for evaluation of bitumen's fracture properties. The material constitutive model developed by Ossa et al. [40,41] which was implemented into a FE code by Costanzi [18] is combined with cohesive zone models (CZM) to simulate the fracture behaviour of pure bitumen. Experimental and numerical results are presented in the form of failure mechanism maps where ductile, brittle and brittle-ductile transition regimes of fracture behaviour are classified. The FE predictions of fracture behaviour match well with experimental results. © 2012 Elsevier Ltd.

Mechanics of capillary forming of aligned carbon nanotube assemblies.

Elastocapillary self-assembly is emerging as a versatile technique to manufacture three-dimensional (3D) microstructures and complex surface textures from arrangements of micro- and nanoscale filaments. Understanding the mechanics of capillary self-assembly is essential to engineering of properties such as shape-directed actuation, anisotropic wetting and adhesion, and mechanical energy transfer and dissipation. We study elastocapillary self-assembly (herein called "capillary forming") of carbon nanotube (CNT) microstructures, combining in situ optical imaging, micromechanical testing, and finite element modeling. By imaging, we identify sequential stages of liquid infiltration, evaporation, and solid shrinkage, whose kinetics relate to the size and shape of the CNT microstructure. We couple these observations with measurements of the orthotropic elastic moduli of CNT forests to understand how the dynamic of shrinkage of the vapor-liquid interface is coupled to the compression of the forest. We compare the kinetics of shrinkage to the rate of evporation from liquid droplets having the same size and geometry. Moreover, we show that the amount of shrinkage during evaporation is governed by the ability of the CNTs to slip against one another, which can be manipulated by the deposition of thin conformal coatings on the CNTs by atomic layer deposition (ALD). This insight is confirmed by finite element modeling of pairs of CNTs as corrugated beams in contact and highlights the coupled role of elasticity and friction in shrinkage and stability of nanoporous solids. Overall, this study shows that nanoscale porosity can be tailored via the filament density and adhesion at contact points, which is important to the development of lightweight multifunctional materials.