A subdivision-based implementation of the hierarchical b-spline finite element method

A novel technique is presented to facilitate the implementation of hierarchical b-splines and their interfacing with conventional finite element implementations. The discrete interpretation of the two-scale relation, as common in subdivision schemes, is used to establish algebraic relations between the basis functions and their coefficients on different levels of the hierarchical b-spline basis. The subdivision projection technique introduced allows us first to compute all element matrices and vectors using a fixed number of same-level basis functions. Their subsequent multiplication with subdivision matrices projects them, during the assembly stage, to the correct levels of the hierarchical b-spline basis. The proposed technique is applied to convergence studies of linear and geometrically nonlinear problems in one, two and three space dimensions. © 2012 Elsevier B.V.

Finite element modelling of the sismic behaviour of water front structures

Water front structures have suffered significant damage in many of the recent earthquakes. These include gravity type quay walls, vertically composite walls, cantilever retaining walls, anchored bulkheads and similar structures. One of the primary causes for the poor performance of these classes of structures is the liquefaction of the foundation soil and in some instances liquefaction of the backfill soil. The liquefaction of the soil in-front of the quay wall tends to cause large lateral displacements and rotation of the wall. Often such gravity walls are placed on rubble mound deposited onto the sea bed.This paper presents finite element analyses of such a problem in which strength degradation of the foundation soil and the backfill material will be modelled using PZ mark III constitutive model. The performance of the wall in terms of its lateral displacement, vertical settlement and/or the rotation suffered by the wall will be presented. In addition, the contours of the horizontal and vertical effective stresses and the excess pore pressure ratio will be presented at different time instants together with hyrdraulic gradients. Immediately after the earthquake, the hydraulic gradients indicate migration of pore water into the region below the wall, suggesting further softening of the foundation soil below the wall.

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.

Linear quadratic game and non-cooperative predictive methods for potential application to modelling driver-AFS interactive steering control

This paper is concerned with the modelling of strategic interactions between the human driver and the vehicle active front steering (AFS) controller in a path-following task where the two controllers hold different target paths. The work is aimed at extending the use of mathematical models in representing driver steering behaviour in complicated driving situations. Two game theoretic approaches, namely linear quadratic game and non-cooperative model predictive control (non-cooperative MPC), are used for developing the driver-AFS interactive steering control model. For each approach, the open-loop Nash steering control solution is derived; the influences of the path-following weights, preview and control horizons, driver time delay and arm neuromuscular system (NMS) dynamics are investigated, and the CPU time consumed is recorded. It is found that the two approaches give identical time histories as well as control gains, while the non-cooperative MPC method uses much less CPU time. Specifically, it is observed that the introduction of weight on the integral of vehicle lateral displacement error helps to eliminate the steady-state path-following error; the increase in preview horizon and NMS natural frequency and the decline in time delay and NMS damping ratio improve the path-following accuracy. © 2013 Copyright Taylor and Francis Group, LLC.

Semi-implicit second order finite volume jet stream computations

Semi-implicit, second order temporal and spatial finite volume computations of the flow in a differentially heated rotating annulus are presented. For the regime considered, three cyclones and anticyclones separated by a relatively fast moving jet of fluid or "jet stream" are predicted. Two second order methods are compared with, first order spatial predictions, and experimental measurements. Velocity vector plots are used to illustrate the predicted flow structure. Computations made using second order central differences are shown to agree best with experimental measurements, and to be stable for integrations over long time periods (> 1000s). No periodic smoothing is required to prevent divergence.

Finite element analysis of compressor disc heat transfer using the transient conduction method with sparse boundary collocation points

Surface temperature measurements from two discs of a gas turbine compressor rig are used as boundary conditions for the transient conduction solution (inverse heat transfer analysis). The disc geometry is complex, and so the finite element method is used. There are often large radial temperature gradients on the discs, and the equations are therefore solved taking into account the dependence of thermal conductivity on temperature. The solution technique also makes use of a multigrid algorithm to reduce the solution time. This is particularly important since a large amount of data must be analyzed to obtain correlations of the heat transfer. The finite element grid is also used for a network analysis to calculate the radiant heat transfer in the cavity formed between the two compressor discs. The work discussed here proved particularly challenging as the disc temperatures were only measured at four different radial locations. Four methods of surface temperature interpolation are examined, together with their effect on the local heat fluxes. It is found that the choice of interpolation method depends on the available number of data points. Bessel interpolation gives the best results for four data points, whereas cubic splines are preferred when there are considerably more data points. The results from the analysis of the compressor rig data show that the heat transfer near the disc inner radius appears to be influenced by the central throughflow. However, for larger radii, the heat transfer from the discs and peripheral shroud is found to be consistent with that of a buoyancy-induced flow.

State estimation schemes for independent component coupled hidden Markov models

Conventional Hidden Markov models generally consist of a Markov chain observed through a linear map corrupted by additive noise. This general class of model has enjoyed a huge and diverse range of applications, for example, speech processing, biomedical signal processing and more recently quantitative finance. However, a lesser known extension of this general class of model is the so-called Factorial Hidden Markov Model (FHMM). FHMMs also have diverse applications, notably in machine learning, artificial intelligence and speech recognition [13, 17]. FHMMs extend the usual class of HMMs, by supposing the partially observed state process is a finite collection of distinct Markov chains, either statistically independent or dependent. There is also considerable current activity in applying collections of partially observed Markov chains to complex action recognition problems, see, for example, [6]. In this article we consider the Maximum Likelihood (ML) parameter estimation problem for FHMMs. Much of the extant literature concerning this problem presents parameter estimation schemes based on full data log-likelihood EM algorithms. This approach can be slow to converge and often imposes heavy demands on computer memory. The latter point is particularly relevant for the class of FHMMs where state space dimensions are relatively large. The contribution in this article is to develop new recursive formulae for a filter-based EM algorithm that can be implemented online. Our new formulae are equivalent ML estimators, however, these formulae are purely recursive and so, significantly reduce numerical complexity and memory requirements. A computer simulation is included to demonstrate the performance of our results. © Taylor & Francis Group, LLC.

Synchronization with partial state feedback on SO(n)

In this paper we consider the problem of constructing a distributed feedback law to achieve synchronization for a group of k agents whose states evolve on SO(n) and which exchange only partial state information along communication links. The partial state information is given by the action of the state on reference vectors in ℝn. We propose a gradient based control law which achieves exponential local convergence to a synchronization configuration under a rank condition on a generalized Laplacian matrix. Furthermore, we discuss the case of time-varying reference vectors and provide a convergence result for this case. The latter helps reach synchronization, requiring less communication links and weaker conditions on the instantaneous reference vectors. Our methods are illustrated on an attitude synchronization problem where agents exchange only their relative positions observed in the respective body frames. ©2009 IEEE.

Trading the stability of finite zeros for global stabilization of nonlinear cascade systems

This note analyzes the stabilizability properties of nonlinear cascades in which a nonminimum phase linear system is interconnected through its output to a Stable nonlinear system. It is shown that the instability of the zeros of the linear System can be traded with the stability of the nonlinear system up to a limit fixed by the growth properties of the cascade interconnection term. Below this limit, global stabilization is achieved by smooth static-state feedback. Beyond this limit, various examples illustrate that controllability of the cascade may be lost, making it impossible to achieve large regions of attractions.

Trading the stability of finite zeros for global stabilization of nonlinear cascade systems

This paper analyzes the stabilizability properties of nonlinear cascades in which a nonminimum phase linear system is interconnected through its output to a stable nonlinear system. It is shown that the instability of the zeros of the linear system can be traded with the stability of the nonlinear system up to a limit fixed by the growth properties of the cascade interconnection term. Below this limit, global stabilization is achieved by smooth static state feedback. Beyond this limit, various examples illustrate that controllability of the cascade may be lost, making it impossible to achieve large regions of attractions.

Efficient decoy-state quantum key distribution with quantified security

We analyse the finite-size security of the efficient Bennett-Brassard 1984 protocol implemented with decoy states and apply the results to a gigahertz-clocked quantum key distribution system. Despite the enhanced security level, the obtained secure key rates are the highest reported so far at all fibre distances.

Statistics anxiety, state anxiety during an examination, and academic achievement.

BACKGROUND: A large proportion of students identify statistics courses as the most anxiety-inducing courses in their curriculum. Many students feel impaired by feelings of state anxiety in the examination and therefore probably show lower achievements. AIMS: The study investigates how statistics anxiety, attitudes (e.g., interest, mathematical self-concept) and trait anxiety, as a general disposition to anxiety, influence experiences of anxiety as well as achievement in an examination. SAMPLE: Participants were 284 undergraduate psychology students, 225 females and 59 males. METHODS: Two weeks prior to the examination, participants completed a demographic questionnaire and measures of the STARS, the STAI, self-concept in mathematics, and interest in statistics. At the beginning of the statistics examination, students assessed their present state anxiety by the KUSTA scale. After 25 min, all examination participants gave another assessment of their anxiety at that moment. Students' examination scores were recorded. Structural equation modelling techniques were used to test relationships between the variables in a multivariate context. RESULTS: Statistics anxiety was the only variable related to state anxiety in the examination. Via state anxiety experienced before and during the examination, statistics anxiety had a negative influence on achievement. However, statistics anxiety also had a direct positive influence on achievement. This result may be explained by students' motivational goals in the specific educational setting. CONCLUSIONS: The results provide insight into the relationship between students' attitudes, dispositions, experiences of anxiety in the examination, and academic achievement, and give recommendations to instructors on how to support students prior to and in the examination.

A fixed-grid b-spline finite element technique for fluid-structure interaction

We present a fixed-grid finite element technique for fluid-structure interaction problems involving incompressible viscous flows and thin structures. The flow equations are discretised with isoparametric b-spline basis functions defined on a logically Cartesian grid. In addition, the previously proposed subdivision-stabilisation technique is used to ensure inf-sup stability. The beam equations are discretised with b-splines and the shell equations with subdivision basis functions, both leading to a rotation-free formulation. The interface conditions between the fluid and the structure are enforced with the Nitsche technique. The resulting coupled system of equations is solved with a Dirichlet-Robin partitioning scheme, and the fluid equations are solved with a pressure-correction method. Auxiliary techniques employed for improving numerical robustness include the level-set based implicit representation of the structure interface on the fluid grid, a cut-cell integration algorithm based on marching tetrahedra and the conservative data transfer between the fluid and structure discretisations. A number of verification and validation examples, primarily motivated by animal locomotion in air or water, demonstrate the robustness and efficiency of our approach. © 2013 John Wiley & Sons, Ltd.

An important pitfall of pseudo-static finite element analysis

Finite Element (FE) pseudo-static analysis can provide a good compromise between simplified methods of dynamic analysis and time domain analysis. The pseudo-static FE approach can accurately model the in situ, stresses prior to seismic loading (when it follows a static analysis simulating the construction sequence) is relatively simple and not as computationally expensive as the time domain approach. However this method should be used with caution as the results can be sensitive to the choice of the mesh dimensions. In this paper two simple examples of pseudo-static finite element analysis are examined parametrically, a homogeneous slope and a cantilever retaining wall, exploring the sensitivity of the pseudo-static analysis results on the adopted mesh size. The mesh dependence was found to be more pronounced for problems with high critical seismic coefficients values (e.g. gentle slopes or small walls), as in these cases a generalised layer failure mechanism is developed simultaneously with the slope or wall mechanism. In general the mesh width was found not to affect notably the predicted value of critical seismic coefficient but to have a major impact on the predicted movements. © 2012 Elsevier Ltd.

Heavy metals in plants and phytoremediation - A state-of-the-art report with special reference to literature published in Chinese journals

Goal, Scope and Background. In some cases, soil, water and food are heavily polluted by heavy metals in China. To use plants to remediate heavy metal pollution would be an effective technique in pollution control. The accumulation of heavy metals in plants and the role of plants in removing pollutants should be understood in order to implement phytoremediation, which makes use of plants to extract, transfer and stabilize heavy metals from soil and water. Methods. The information has been compiled from Chinese publications stemming mostly from the last decade, to show the research results on heavy metals in plants and the role of plants in controlling heavy metal pollution, and to provide a general outlook of phytoremediation in China. Related references from scientific journals and university journals are searched and summarized in sections concerning the accumulation of heavy metals in plants, plants for heavy metal purification and phytoremediation techniques. Results and Discussion. Plants can take up heavy metals by their roots, or even via their stems and leaves, and accumulate them in their organs. Plants take up elements selectively. Accumulation and distribution of heavy metals in the plant depends on the plant species, element species, chemical and bioavailiability, redox, pH, cation exchange capacity, dissolved oxygen, temperature and secretion of roots. Plants are employed in the decontamination of heavy metals from polluted water and have demonstrated high performances in treating mineral tailing water and industrial effluents. The purification capacity of heavy metals by plants are affected by several factors, such as the concentration of the heavy metals, species of elements, plant species, exposure duration, temperature and pH. Conclusions. Phytoremediation, which makes use of vegetation to remove, detoxify, or stabilize persistent pollutants, is a green and environmentally-friendly tool for cleaning polluted soil and water. The advantage of high biomass productive and easy disposal makes plants most useful to remediate heavy metals on site. Recommendations and Outlook. Based on knowledge of the heavy metal accumulation in plants, it is possible to select those species of crops and pasturage herbs, which accumulate fewer heavy metals, for food cultivation and fodder for animals; and to select those hyperaccumulation species for extracting heavy metals from soil and water. Studies on the mechanisms and application of hyperaccumulation are necessary in China for developing phytoremediation.