Image-based vs. mesh-based statistical appearance models of the human femur: Implications for finite element simulations.

Statistical appearance models have recently been introduced in bone mechanics to investigate bone geometry and mechanical properties in population studies. The establishment of accurate anatomical correspondences is a critical aspect for the construction of reliable models. Depending on the representation of a bone as an image or a mesh, correspondences are detected using image registration or mesh morphing. The objective of this study was to compare image-based and mesh-based statistical appearance models of the femur for finite element (FE) simulations. To this aim, (i) we compared correspondence detection methods on bone surface and in bone volume; (ii) we created an image-based and a mesh-based statistical appearance models from 130 images, which we validated using compactness, representation and generalization, and we analyzed the FE results on 50 recreated bones vs. original bones; (iii) we created 1000 new instances, and we compared the quality of the FE meshes. Results showed that the image-based approach was more accurate in volume correspondence detection and quality of FE meshes, whereas the mesh-based approach was more accurate for surface correspondence detection and model compactness. Based on our results, we recommend the use of image-based statistical appearance models for FE simulations of the femur.

Nonlinear quasi-static finite element simulations predict in vitro strength of human proximal femora assessed in a dynamic sideways fall setup

Osteoporotic proximal femur fractures are caused by low energy trauma, typically when falling on the hip from standing height. Finite element simulations, widely used to predict the fracture load of femora in fall, usually include neither mass-related inertial effects, nor the viscous part of bone's material behavior. The aim of this study was to elucidate if quasi-static non-linear homogenized finite element analyses can predict in vitro mechanical properties of proximal femora assessed in dynamic drop tower experiments. The case-specific numerical models of thirteen femora predicted the strength (R2=0.84, SEE=540 N, 16.2%), stiffness (R2=0.82, SEE=233 N/mm, 18.0%) and fracture energy (R2=0.72, SEE=3.85 J, 39.6%); and provided fair qualitative matches with the fracture patterns. The influence of material anisotropy was negligible for all predictions. These results suggest that quasi-static homogenized finite element analysis may be used to predict mechanical properties of proximal femora in the dynamic sideways fall situation.

Extraordinary transmission through a single coaxial aperture in a thin metal film

We investigate experimentally the transmission properties of single sub-wavelength coaxial apertures in thin metal films (t = 110 nm). Enhanced transmission through a single sub-wavelength coaxial aperture illuminated with a strongly focused radially polarized light beam is reported. In our experiments we achieved up to four times enhanced transmission through a single coaxial aperture as compared to a (hollow) circular aperture with the same outer diameter.We attribute this enhancement of transmission to the excitation of a TEM-mode for illumination with radially polarized light inside the single coaxial aperture. A strong polarization contrast is observed between the transmission for radially and azimuthally polarized illumination. Furthermore, the observed transmission through a single coaxial aperture can be strongly reduced if surface plasmons are excited. The experimental results are in good agreement with finite difference time domain (FDTD) simulations.

Patient-specific finite-element simulation of the human cornea: A clinical validation study on cataract surgery

The planning of refractive surgical interventions is a challenging task. Numerical modeling has been proposed as a solution to support surgical intervention and predict the visual acuity, but validation on patient specific intervention is missing. The purpose of this study was to validate the numerical predictions of the post-operative corneal topography induced by the incisions required for cataract surgery. The corneal topography of 13 patients was assessed preoperatively and postoperatively (1-day and 30-day follow-up) with a Pentacam tomography device. The preoperatively acquired geometric corneal topography – anterior, posterior and pachymetry data – was used to build patient-specific finite element models. For each patient, the effects of the cataract incisions were simulated numerically and the resulting corneal surfaces were compared to the clinical postoperative measurements at one day and at 30-days follow up. Results showed that the model was able to reproduce experimental measurements with an error on the surgically induced sphere of 0.38D one day postoperatively and 0.19D 30 days postoperatively. The standard deviation of the surgically induced cylinder was 0.54D at the first postoperative day and 0.38D 30 days postoperatively. The prediction errors in surface elevation and curvature were below the topography measurement device accuracy of ±5μm and ±0.25D after the 30-day follow-up. The results showed that finite element simulations of corneal biomechanics are able to predict post cataract surgery within topography measurement device accuracy. We can conclude that the numerical simulation can become a valuable tool to plan corneal incisions in cataract surgery and other ophthalmosurgical procedures in order to optimize patients' refractive outcome and visual function.

On the initiation of boudinage and folding instabilities in layered elasto-visco-plastic solids - Insights from natural microstructures and numerical simulations

The present understanding of the initiation of boudinage and folding structures is based on viscosity contrasts and stress exponents, considering an intrinsically unstable state of the layer. The criterion of localization is believed to be prescribed by geometry-material interactions, which are often encountered in natural structures. An alternative localization phenomenon has been established for ductile materials, in which instability emerges for critical material parameters and loading rates from homogeneous conditions. In this thesis, conditions are sought under which this type of instability prevails and whether localization in geological materials necessarily requires a trigger by geometric imperfections. The relevance of critical deformation conditions, material parameters and the spatial configuration of instabilities are discussed in a geological context. In order to analyze boudinage geometries, a numerical eigenmode analysis is introduced. This method allows determining natural frequencies and wavelengths of a structure and inducing perturbations on these frequencies. In the subsequent coupled thermo-mechanical simulations, using a grain size evolution and end-member flow laws, localization emerges when material softening through grain size sensitive viscous creep sets in. Pinch-and-swell structures evolve along slip lines through a positive feedback between the matrix response and material bifurcations inside the layer, independent from the mesh-discretization length scale. Since boudinage and folding are considered to express the same general instability, both structures should arise independently of the sign of the loading conditions and for identical material parameters. To this end, the link between material to energy instabilities is approached by means of bifurcation analyses of the field equations and finite element simulations of the coupled system of equations. Boudinage and folding structures develop at the same critical energy threshold, where dissipative work by temperature-sensitive creep overcomes the diffusive capacity of the layer. This finding provides basis for a unified theory for strain localization in layered ductile materials. The numerical simulations are compared to natural pinch-and-swell microstructures, tracing the adaption of grain sizes, textures and creep mechanisms in calcite veins. The switch from dislocation to diffusion creep relates to strain-rate weakening, which is induced by dissipated heat from grain size reduction, and marks the onset of continuous necking. The time-dependent sequence uncovers multiple steady states at different time intervals. Microstructurally and mechanically stable conditions are finally expressed in the pinch-and-swell end members. The major outcome of this study is that boudinage and folding can be described as the same coupled energy-mechanical bifurcation, or as one critical energy attractor. This finding allows the derivation of critical deformation conditions and fundamental material parameters directly from localized structures in the field.

Optical properties of a nanomatch-like plasmonic structure

The optical properties of a match-like plasmonic nanostructure are numerically investigated using full-wave finite-difference time-domain analysis in conjunction with dispersive material models. This work is mainly motivated by the developed technique enabling reproducible fabrication of nanomatch structures as well as the growing applications that utilize the localized field enhancement in plasmonic nanostructures. Our research revealed that due to the pronounced field enhancement and larger resonance tunabilities, some nanomatch topologies show potentials for various applications in the field of, e.g., sensing as well as a novel scheme for highly reproducible tips in scanning near field optical microscopy, among others. Despite the additional degrees of freedom that are offered by the composite nature of the proposed nanomatch topology, the paper also reflects on a fundamental complication intrinsic to the material interfaces especially in the nanoscale: stoichiometric mixing. We conclude that the specificity in material modeling will become a significant issue in future research on functionalized composite nanostructures.

Local strain distribution in real three-dimensional alveolar geometries

Mechanical ventilation is not only a life saving treatment but can also cause negative side effects. One of the main complications is inflammation caused by overstretching of the alveolar tissue. Previously, studies investigated either global strains or looked into which states lead to inflammatory reactions in cell cultures. However, the connection between the global deformation, of a tissue strip or the whole organ, and the strains reaching the single cells lining the alveolar walls is unknown and respective studies are still missing. The main reason for this is most likely the complex, sponge-like alveolar geometry, whose three-dimensional details have been unknown until recently. Utilizing synchrotron-based X-ray tomographic microscopy, we were able to generate real and detailed three-dimensional alveolar geometries on which we have performed finite-element simulations. This allowed us to determine, for the first time, a three-dimensional strain state within the alveolar wall. Briefly, precision-cut lung slices, prepared from isolated rat lungs, were scanned and segmented to provide a three-dimensional geometry. This was then discretized using newly developed tetrahedral elements. The main conclusions of this study are that the local strain in the alveolar wall can reach a multiple of the value of the global strain, for our simulations up to four times as high and that thin structures obviously cause hotspots that are especially at risk of overstretching.

The influence of yield surface shape and damage in the depth-dependent response of bone tissue to nanoindentation using spherical and Berkovich indenters

Prevention and treatment of osteoporosis rely on understanding of the micromechanical behaviour of bone and its influence on fracture toughness and cell-mediated adaptation processes. Postyield properties may be assessed by nonlinear finite element simulations of nanoindentation using elastoplastic and damage models. This computational study aims at determining the influence of yield surface shape and damage on the depth-dependent response of bone to nanoindentation using spherical and conical tips. Yield surface shape and damage were shown to have a major impact on the indentation curves. Their influence on indentation modulus, hardness, their ratio as well as the elastic-to-total work ratio is well described by multilinear regressions for both tip shapes. For conical tips, indentation depth was not statistically significant (p<0.0001). For spherical tips, damage was not a significant parameter (p<0.0001). The gained knowledge can be used for developing an inverse method for identification of postelastic properties of bone from nanoindentation.

Analysis of fractal electrodes for efficient neural stimulation

Planar electrodes are increasingly used in therapeutic neural stimulation techniques such as functional electrical stimulation, epidural spinal cord stimulation (ESCS), and cortical stimulation. Recently, optimized electrode geometries have been shown to increase the efficiency of neural stimulation by increasing the variation of current density on the electrode surface. In the present work, a new family of modified fractal electrode geometries is developed to enhance the efficiency of neural stimulation. It is shown that a promising approach in increasing the neural activation function is to increase the "edginess" of the electrode surface, a concept that is explained and quantified by fractal mathematics. Rigorous finite element simulations were performed to compute electric potential produced by proposed modified fractal geometries. The activation of 256 model axons positioned around the electrodes was then quantified, showing that modified fractal geometries required a 22% less input power while maintaining the same level of neural activation. Preliminary in vivo experiments investigating muscle evoked potentials due to median nerve stimulation showed encouraging results, supporting the feasibility of increasing neural stimulation efficiency using modified fractal geometries.

An adjustable focusing system for a 2 MeV H- ion beam line based on permanent magnet quadrupoles

A compact adjustable focusing system for a 2 MeV H- RFQ Linac is designed, constructed and tested based on four permanent magnet quadrupoles (PMQ). A PMQ model is realised using finite element simulations, providing an integrated field gradient of 2.35 T with a maximal field gradient of 57 T/m. A prototype is constructed and the magnetic field is measured, demonstrating good agreement with the simulation. Particle track simulations provide initial values for the quadrupole positions. Accordingly, four PMQs are constructed and assembled on the beam line, their positions are then tuned to obtain a minimal beam spot size of (1.2 x 2.2) mm^2 on target. This paper describes an adjustable PMQ beam line for an external ion beam. The novel compact design based on commercially available NdFeB magnets allows high flexibility for ion beam applications.

Development of a balanced experimental–computational approach to understanding the mechanics of proximal femur fractures

The majority of people who sustain hip fractures after a fall to the side would not have been identified using current screening techniques such as areal bone mineral density. Identifying them, however, is essential so that appropriate pharmacological or lifestyle interventions can be implemented. A protocol, demonstrated on a single specimen, is introduced, comprising the following components; in vitro biofidelic drop tower testing of a proximal femur; high-speed image analysis through digital image correlation; detailed accounting of the energy present during the drop tower test; organ level finite element simulations of the drop tower test; micro level finite element simulations of critical volumes of interest in the trabecular bone. Fracture in the femoral specimen initiated in the superior part of the neck. Measured fracture load was 3760 N, compared to 4871 N predicted based on the finite element analysis. Digital image correlation showed compressive surface strains as high as 7.1% prior to fracture. Voxel level results were consistent with high-speed video data and helped identify hidden local structural weaknesses. We found using a drop tower test protocol that a femoral neck fracture can be created with a fall velocity and energy representative of a sideways fall from standing. Additionally, we found that the nested explicit finite element method used allowed us to identify local structural weaknesses associated with femur fracture initiation.

A numerical model for inert gas transport in the lung based on fractal airway morphology

Patients suffering from cystic fibrosis (CF) show thick secretions, mucus plugging and bronchiectasis in bronchial and alveolar ducts. This results in substantial structural changes of the airway morphology and heterogeneous ventilation. Disease progression and treatment effects are monitored by so-called gas washout tests, where the change in concentration of an inert gas is measured over a single or multiple breaths. The result of the tests based on the profile of the measured concentration is a marker for the severity of the ventilation inhomogeneity strongly affected by the airway morphology. However, it is hard to localize underlying obstructions to specific parts of the airways, especially if occurring in the lung periphery. In order to support the analysis of lung function tests (e.g. multi-breath washout), we developed a numerical model of the entire airway tree, coupling a lumped parameter model for the lung ventilation with a 4th-order accurate finite difference model of a 1D advection-diffusion equation for the transport of an inert gas. The boundary conditions for the flow problem comprise the pressure and flow profile at the mouth, which is typically known from clinical washout tests. The natural asymmetry of the lung morphology is approximated by a generic, fractal, asymmetric branching scheme which we applied for the conducting airways. A conducting airway ends when its dimension falls below a predefined limit. A model acinus is then connected to each terminal airway. The morphology of an acinus unit comprises a network of expandable cells. A regional, linear constitutive law describes the pressure-volume relation between the pleural gap and the acinus. The cyclic expansion (breathing) of each acinus unit depends on the resistance of the feeding airway and on the flow resistance and stiffness of the cells themselves. Special care was taken in the development of a conservative numerical scheme for the gas transport across bifurcations, handling spatially and temporally varying advective and diffusive fluxes over a wide range of scales. Implicit time integration was applied to account for the numerical stiffness resulting from the discretized transport equation. Local or regional modification of the airway dimension, resistance or tissue stiffness are introduced to mimic pathological airway restrictions typical for CF. This leads to a more heterogeneous ventilation of the model lung. As a result the concentration in some distal parts of the lung model remains increased for a longer duration. The inert gas concentration at the mouth towards the end of the expirations is composed of gas from regions with very different washout efficiency. This results in a steeper slope of the corresponding part of the washout profile.