hp-dGFEM for Second-Order Elliptic Problems in Plyhedra II: Exponential Convergence

The goal of this paper is to establish exponential convergence of $hp$-version interior penalty (IP) discontinuous Galerkin (dG) finite element methods for the numerical approximation of linear second-order elliptic boundary-value problems with homogeneous Dirichlet boundary conditions and piecewise analytic data in three-dimensional polyhedral domains. More precisely, we shall analyze the convergence of the $hp$-IP dG methods considered in [D. Schötzau, C. Schwab, T. P. Wihler, SIAM J. Numer. Anal., 51 (2013), pp. 1610--1633] based on axiparallel $\sigma$-geometric anisotropic meshes and $\bm{s}$-linear anisotropic polynomial degree distributions.

Energy harvesting through arterial wall deformation: A FEM approach to fluid–structure interactions and magneto-hydrodynamics

Engineers are confronted with the energy demand of active medical implants in patients with increasing life expectancy. Scavenging energy from the patient’s body is envisioned as an alternative to conventional power sources. Joining in this effort towards human-powered implants, we propose an innovative concept that combines the deformation of an artery resulting from the arterial pressure pulse with a transduction mechanism based on magneto-hydrodynamics. To overcome certain limitations of a preliminary analytical study on this topic, we demonstrate here a more accurate model of our generator by implementing a three-dimensional multiphysics finite element method (FEM) simulation combining solid mechanics, fluid mechanics, electric and magnetic fields as well as the corresponding couplings. This simulation is used to optimize the generator with respect to several design parameters. A first validation is obtained by comparing the results of the FEM simulation with those of the analytical approach adopted in our previous study. With an expected overall conversion efficiency of 20% and an average output power of 30 μW, our generator outperforms previous devices based on arterial wall deformation by more than two orders of magnitude. Most importantly, our generator provides sufficient power to supply a cardiac pacemaker.

Subharmonic Mixer Based on EBG Technology

We present the design of a submillimeter-wave mixer based on electromagnetic band gap (EBG) technology and using subharmonic local oscillator (LO) injection. The indicated device converts an incoming submilimeter wavelength signal into a 1-5 GHz intermediate frequency (IF) signal by mixing it with a subharmonic LO signal. The mixer consists of a dual-band receiver and two coplanar stripline (CPS) filters, collocated on top of a three-dimensional (3-D) EBG structure. A four-element array of the proposed receivers was designed, fabricated and tested. The configuration demonstrated reasonable performance: conversion loss below 8 dB and noise temperature below 3000 K. The presented concept can be used for higher frequencies, provided the availability of sufficiently powerful LO sources.

Patient-specific spinal stiffness in AIS: a preoperative and noninvasive method.

INTRODUCTION The clinical tests currently used to assess spinal biomechanics preoperatively are unable to assess true mechanical spinal stiffness. They rely on spinal displacement without considering the force required to deform a patient's spine. We propose a preoperative method for noninvasively quantifying the three-dimensional patient-specific stiffness of the spines of adolescent idiopathic scoliosis patients. METHODS The technique combines a novel clinical test with numerical optimization of a finite element model of the patient's spine. RESULTS A pilot study conducted on five patients showed that the model was able to provide accurate 3D reconstruction of the spine's midline and predict the spine's stiffness for each patient in flexion, bending, and rotation. Statistically significant variation of spinal stiffness was observed between the patients. CONCLUSION This result confirms that spinal biomechanics is patient-specific, which should be taken into consideration to individualize surgical treatment.

Structure of the histone mRNA hairpin required for cell cycle regulation of histone gene expression.

Expression of replication-dependent histone genes requires a conserved hairpin RNA element in the 3' untranslated regions of poly(A)-less histone mRNAs. The 3' hairpin element is recognized by the hairpin-binding protein or stem-loop-binding protein (HBP/SLBP). This protein-RNA interaction is important for the endonucleolytic cleavage generating the mature mRNA 3' end. The 3' hairpin and presumably HBP/SLBP are also required for nucleocytoplasmic transport, translation, and stability of histone mRNAs. RNA 3' processing and mRNA stability are both regulated during the cell cycle. Here, we have determined the three-dimensional structure of a 24-mer RNA comprising a mammalian histone RNA hairpin using heteronuclear multidimensional NMR spectroscopy. The hairpin adopts a novel UUUC tetraloop conformation that is stabilized by base stacking involving the first and third loop uridines and a closing U-A base pair, and by hydrogen bonding between the first and third uridines in the tetraloop. The HBP interaction of hairpin RNA variants was analyzed in band shift experiments. Particularly important interactions for HBP recognition are mediated by the closing U-A base pair and the first and third loop uridines, whose Watson-Crick functional groups are exposed towards the major groove of the RNA hairpin. The results obtained provide novel structural insight into the interaction of the histone 3' hairpin with HBP, and thus the regulation of histone mRNA metabolism.

Two-Dimensional Crystallisation of Membrane Proteins and Structural Assessment

Two-dimensional (2D) crystallisation of Membrane proteins reconstitutes them into their native environment, the lipid bilayer. Electron crystallography allows the structural analysis of these regular protein–lipid arrays up to atomic resolution. The crystal quality depends on the protein purity, ist stability and on the crystallisation conditions. The basics of 2D crystallisation and different recent advances are reviewed and electron crystallography approaches summarised. Progress in 2D crystallisation, sample preparation, image detectors and automation of the data acquisition and processing pipeline makes 2D electron crystallography particularly attractive for the structural analysis of membrane proteins that are too small for single-particle analyses and too unstable to form three-dimensional (3D) crystals.

Correlation between alveolar ventilation and electrical properties of lung parenchyma.

One key problem in modern medical imaging is linking measured data and actual physiological quantities. In this article we derive such a link between the electrical bioimpedance of lung parenchyma, which can be measured by electrical impedance tomography (EIT), and the magnitude of regional ventilation, a key to understanding lung mechanics and developing novel protective ventilation strategies. Two rat-derived three-dimensional alveolar microstructures obtained from synchrotron-based x-ray tomography are each exposed to a constant potential difference for different states of ventilation in a finite element simulation. While the alveolar wall volume remains constant during stretch, the enclosed air volume varies, similar to the lung volume during ventilation. The enclosed air, serving as insulator in the alveolar ensemble, determines the resulting current and accordingly local tissue bioimpedance. From this we can derive a relationship between lung tissue bioimpedance and regional alveolar ventilation. The derived relationship shows a linear dependence between air content and tissue impedance and matches clinical data determined from a ventilated patient at the bedside.

Modelos de elementos finitos con aplicación tridimensional de la analogía fluido-fibra, para predecir el comportamiento mecánico de piezas estructurales de madera bajo la presencia de los principales tipos de nudos

Los nudos son los defectos que más disminuyen la resistencia de piezas de madera en la escala estructural, al ocasionar no solo una discontinuidad material, sino también la desviación de las fibras que se encuentran a su alrededor. En la década de los 80 se introdujo la teoría de la analogía fluido-fibra, como un método que aproximaba adecuadamente todas estas desviaciones. Sin embargo en su aplicación tridimensional, nunca se consideraron las diferencias geométricas en el sentido perpendicular al eje longitudinal de las piezas estructurales, lo cual imposibilitaba la simulación numérica de algunos de los principales tipos de nudos, y disminuía la precisión obtenida en aquellos nudos en los que la modelización sí era viable. En este trabajo se propone un modelo programado en lenguaje paramétrico de un software de elementos finitos que, bajo una formulación en tres dimensiones más general, permitirá estudiar de forma automatizada el comportamiento estructural de la madera bajo la influencia de los principales tipos de nudos, a partir de la geometría visible de los mismos y la posición de la médula en la pieza, y el cual ha sido contrastado experimentalmente, simulando de forma muy precisa el comportamiento mecánico de vigas sometidas a ensayos de flexión a cuatro puntos. Knots are the defects that most reduce the strength of lumber at the structural level, by causing not only a material discontinuity but also the deviation of the fibers that surround them. In the 80's it was introduced the theory of the flow-grain analogy as a method to approximating adequately these deviations. However, in three-dimensional applications, geometrical differences in the direction perpendicular to the longitudinal axis of the structural specimens were never considered before, which prevented the numerical simulation of some of the main types of knots, and decreased the achieved precision in those kind of knots where modeling itself was possible. This paper purposes a parametric model programmed in a finite element software, in the way that with a more general three-dimensional formulation, an automated study of the structural behavior of timber under the influence of the main types of knots is allowed by only knowing the visible geometry of such defects, and the position of the pith. Furthermore that has been confirmed experimentally obtaining very accurately simulations of the mechanical behavior of beams under four points bending test.

Improved design of an overhead rail current conductor for application in underground lines

Overhead rail current collector systems for railway traction offer certain features, such as low installation height and reduced maintenance, which make them predominantly suitable for use in underground train infrastructures. Due to the increased demands of modern catenary systems and higher running speeds of new vehicles, a more capable design of the conductor rail is needed. A new overhead conductor rail has been developed and its design has been patented [13]. Modern simulation and modelling techniques were used in the development approach. The new conductor rail profile has a dynamic behaviour superior to that of the system currently in use. Its innovative design permits either an increase of catenary support spacing or a higher vehicle running speed. Both options ensure savings in installation or operating costs. The simulation model used to optimise the existing conductor rail profile included both a finite element model of the catenary and a three-dimensional multi-body system model of the pantograph. The contact force that appears between pantograph and catenary was obtained in simulation. A sensitivity analysis of the key parameters that influence in catenary dynamics was carried out, finally leading to the improved design.

Dynamic stiffness analysis of bridge abutments

Strong motion records obtained in instrumented short-span bridges show the importance of the abutments in the dynamic response of the structure. Existing models study the pier foundation influence but not the abutment performance. This work proposes two and three dimensional boundary element models in the frequency domain and studies the dimensionless dynamic stiffness of standard bridge abutments.

Analytic free-form lens design in 3D: coupling three ray sets using two lens surfaces

The two-dimensional analytic optics design method presented in a previous paper [Opt. Express 20, 5576–5585 (2012)] is extended in this work to the three-dimensional case, enabling the coupling of three ray sets with two free-form lens surfaces. Fermat’s principle is used to deduce additional sets of functional differential equations which make it possible to calculate the lens surfaces. Ray tracing simulations demonstrate the excellent imaging performance of the resulting free-form lenses described by more than 100 coefficients.

A methodology to measure the interface shear strength by means of the fiber push-in test

A methodology is presented to measure the fiber/matrix interface shear strength in composites. The strategy is based on performing a fiber push-in test at the central fiber of highly-packed fiber clusters with hexagonal symmetry which are often found in unidirectional composites with a high volume fraction of fibers. The mechanics of this test was analyzed in detail by means of three-dimensional finite element simulations. In particular, the influence of different parameters (interface shear strength, toughness and friction as well as fiber longitudinal elastic modulus and curing stresses) on the critical load at the onset of debonding was established. From the results of the numerical simulations, a simple relationship between the critical load and the interface shear strength is proposed. The methodology was validated in an unidirectional C/epoxy composite and the advantages and limitations of the proposed methodology are indicated.

Low-dimensional modeling of streaks in a wedge flow boundary layer

This paper is concerned with the low dimensional structure of optimal streaks in a wedge flow boundary layer, which have been recently shown to consist of a unique (up to a constant factor) three-dimensional streamwise evolving mode, known as the most unstable streaky mode. Optimal streaks exhibit a still unexplored/unexploited approximate self-similarity (not associated with the boundary layer self-similarity), namely the streamwise velocity re-scaled with their maximum remains almost independent of both the spanwise wavenumber and the streamwise coordinate; the remaining two velocity components instead do not satisfy this property. The approximate self-similar behavior is analyzed here and exploited to further simplify the description of optimal streaks. In particular, it is shown that streaks can be approximately described in terms of the streamwise evolution of the scalar amplitudes of just three one-dimensional modes, providing the wall normal profiles of the streamwise velocity and two combinations of the cross flow velocity components; the scalar amplitudes obey a singular system of three ordinary differential equations (involving only two degrees of freedom), which approximates well the streamwise evolution of the general streaks.

The Dynamics of High-Speed Railway Bridges: A Review of Design Issues and New Research for Lateral Dynamics

The dynamic effects of high-speed trains on viaducts are important issues for the design of the structures, as well as for the consideration of safe running conditions for the trains. In this work we start by reviewing the relevance of some basic design aspects. The significance of impact factor envelopes for moving loads is considered first. Resonance which may be achieved for high-speed trains requires dynamic analysis, for which some key aspects are discussed. The relevance of performing a longitudinal distribution of axle loads, the number of modes taken in analysis, and the consideration of vehicle-structure interaction are discussed with representative examples. The lateral dynamic effects of running trains on bridges is of importance for laterally compliant viaducts, such as some very tall structures erected in new high-speed lines. The relevance of this study is mainly for the safety of the traffic, considering both internal actions such as the hunting motion as well as external actions such as wind or earthquakes [1]. These studies require three-dimensional dynamic coupled vehicle-bridge models, and consideration of wheel to rail contact, a phenomenon which is complex and costly to model in detail. We describe here a fully nonlinear coupled model, described in absolute coordinates and incorporated into a commercial finite element framework [2]. The wheel-rail contact has been considered using a FastSim algorithm which provides a compromise between accuracy and computational cost, and captures the main nonlinear response of the contact interface. Two applications are presented, firstly to a vehicle subject to a strong wind gust traversing a bridge, showing the relevance of the nonlinear wheel-rail contact model as well as the dynamic interaction between bridge and vehicle. The second application is to a real HS viaduct with a long continuous deck and tall piers and high lateral compliance [3]. The results show the safety of the traffic as well as the importance of considering features such as track alignment irregularities.

A comparison between PML, infinite elements and an iterative BEM as mesh truncation methods for HP self-adaptive procedures in electromagnetics

Finite element hp-adaptivity is a technology that allows for very accurate numerical solutions. When applied to open region problems such as radar cross section prediction or antenna analysis, a mesh truncation method needs to be used. This paper compares the following mesh truncation methods in the context of hp-adaptive methods: Infinite Elements, Perfectly Matched Layers and an iterative boundary element based methodology. These methods have been selected because they are exact at the continuous level (a desirable feature required by the extreme accuracy delivered by the hp-adaptive strategy) and they are easy to integrate with the logic of hp-adaptivity. The comparison is mainly based on the number of degrees of freedom needed for each method to achieve a given level of accuracy. Computational times are also included. Two-dimensional examples are used, but the conclusions directly extrapolated to the three dimensional case.