Numerical investigations on the strain-adaptive bone remodelling in the periprosthetic femur: Influence of the boundary conditions

Background: There are several numerical investigations on bone remodelling after total hip arthroplasty (THA) on the basis of the finite element analysis (FEA). For such computations certain boundary conditions have to be defined. The authors chose a maximum of three static load situations, usually taken from the gait cycle because this is the most frequent dynamic activity of a patient after THA. Materials and methods: The numerical study presented here investigates whether it is useful to consider only one static load situation of the gait cycle in the FE calculation of the bone remodelling. For this purpose, 5 different loading cases were examined in order to determine their influence on the change in the physiological load distribution within the femur and on the resulting strain-adaptive bone remodelling. First, four different static loading cases at 25%, 45%, 65% and 85% of the gait cycle, respectively, and then the whole gait cycle in a loading regime were examined in order to regard all the different loadings of the cycle in the simulation. Results: The computed evolution of the apparent bone density (ABD) and the calculated mass losses in the periprosthetic femur show that the simulation results are highly dependent on the chosen boundary conditions. Conclusion: These numerical investigations prove that a static load situation is insufficient for representing the whole gait cycle. This causes severe deviations in the FE calculation of the bone remodelling. However, accompanying clinical examinations are necessary to calibrate the bone adaptation law and thus to validate the FE calculations.

Distribution of Light in the Human Retina under Natural Viewing Conditions

Age-related macular degeneration (AMD) is the leading cause of blindness inAmerica. The fact that AMD wreaks most of the damage in the center of the retina raises the question of whether light, integrated over long periods, is more concentrated in the macula. A method, based on eye-tracking, was developed to measure the distribution of light in the retina under natural viewing conditions. The hypothesis was that integrated over time, retinal illumination peaked in the macula. Additionally a possible relationship between age and retinal illumination was investigated. The eye tracker superimposed the subject’s gaze position on a video recorded by a scene camera. Five informed subjects were employed in feasibility tests, and 58 naïve subjects participated in 5 phases. In phase 1 the subjects viewed a gray-scale image. In phase 2, they observed a sequence of photographic images. In phase 3 they viewed a video. In phase 4, they worked on a computer; in phase 5, the subjects walked around freely. The informed subjects were instructed to gaze at bright objects in the field of view and then at dark objects. Naïve subjects were allowed to gaze freely for all phases. Using the subject’s gaze coordinates, and the video provided by the scene camera, the cumulative light distribution on the retina was calculated for ~15° around the fovea. As expected for control subjects, cumulative retinal light distributions peaked and dipped in the fovea when they gazed at bright or dark objects respectively. The light distribution maps obtained from the naïve subjects presented a tendency to peak in the macula for phases 1, 2, and 3, a consistent tendency in phase 4 and a variable tendency in phase 5. The feasibility of using an eye-tracker system to measure the distribution of light in the retina was demonstrated, thus helping to understand the role played by light exposure in the etiology of AMD. Results showed that a tendency for light to peak in the macula is a characteristic of some individuals and of certain tasks. In these situations, risk of AMD could be increased. No significant difference was observed based on age.

Measure-valued mass evolution problems with flux boundary conditions and solution-dependent velocities

In this paper we prove well-posedness for a measure-valued continuity equation with solution-dependent velocity and flux boundary conditions, posed on a bounded one-dimensional domain. We generalize the results of an earlier paper [J. Differential Equations, 259 (2015), pp. 10681097] to settings where the dynamics are driven by interactions. In a forward-Euler-like approach, we construct a time-discretized version of the original problem and employ those results as a building block within each subinterval. A limit solution is obtained as the mesh size of the time discretization goes to zero. Moreover, the limit is independent of the specific way of partitioning the time interval [0, T]. This paper is partially based on results presented in Chapter 5 of [Evolution Equations for Systems Governed by Social Interactions, Ph.D. thesis, Eindhoven University of Technology, 2015], while a number of issues that were still open there are now resolved.

Automatic finite element formulation and assembly of hyperelastic higher order structural models

The formulation of higher order structural models and their discretization using the finite element method is difficult owing to their complexity, especially in the presence of non-linearities. In this work a new algorithm for automating the formulation and assembly of hyperelastic higher-order structural finite elements is developed. A hierarchic series of kinematic models is proposed for modeling structures with special geometries and the algorithm is formulated to automate the study of this class of higher order structural models. The algorithm developed in this work sidesteps the need for an explicit derivation of the governing equations for the individual kinematic modes. Using a novel procedure involving a nodal degree-of-freedom based automatic assembly algorithm, automatic differentiation and higher dimensional quadrature, the relevant finite element matrices are directly computed from the variational statement of elasticity and the higher order kinematic model. Another significant feature of the proposed algorithm is that natural boundary conditions are implicitly handled for arbitrary higher order kinematic models. The validity algorithm is illustrated with examples involving linear elasticity and hyperelasticity. (C) 2013 Elsevier Inc. All rights reserved.

Cálculo das variações fraccional generalizado

Nesta tese de doutoramento apresentamos um cálculo das variações fraccional generalizado. Consideramos problemas variacionais com derivadas e integrais fraccionais generalizados e estudamo-los usando métodos directos e indirectos. Em particular, obtemos condições necessárias de optimalidade de Euler-Lagrange para o problema fundamental e isoperimétrico, condições de transversalidade e teoremas de Noether. Demonstramos a existência de soluções, num espaço de funções apropriado, sob condições do tipo de Tonelli. Terminamos mostrando a existência de valores próprios, e correspondentes funções próprias ortogonais, para problemas de Sturm- Liouville.

Tropically-driven climate shifts in southwestern Europe during MIS 19, a low eccentricity interglacial

The relative roles of high- versus low-latitude forcing of millennial-scale climate variability are still not well understood. Here we present terrestrial–marine climate profiles from the southwestern Iberian margin, a region particularly affected by precession, that show millennial climate oscillations related to a nonlinear response to the Earth's precession cycle during Marine Isotope Stage (MIS) 19. MIS 19 has been considered the best analogue to our present interglacial from an astronomical point of view due to the reduced eccentricity centred at 785 ka. In our records, seven millennial-scale forest contractions punctuated MIS 19 superimposed to two orbitally-driven Mediterranean forest expansions. In contrast to our present interglacial, we evidence for the first time low latitude-driven 5000-yr cycles of drying and cooling in the western Mediterranean region, along with warmth in the subtropical gyre related to the fourth harmonic of precession. These cycles indicate repeated intensification of North Atlantic meridional moisture transport that along with decrease in boreal summer insolation triggered ice growth and may have contributed to the glacial inception, at ∼774 ka. The freshwater fluxes during MIS 19ab amplified the cooling events in the North Atlantic promoting further cooling and leading to MIS 18 glaciation. The discrepancy between the dominant cyclicity observed during MIS 1, 2500-yr, and that of MIS 19, 5000-yr, challenges the similar duration of the Holocene and MIS 19c interglacials under natural boundary conditions.

Higher-order variational problems of Herglotz type

We obtain a generalized Euler–Lagrange differential equation and transversality optimality conditions for Herglotz-type higher-order variational problems. Illustrative examples of the new results are given.

Natural convection boundary-layer adjacent to an inclined flat plate subject to sudden and ramp heating

The natural convection thermal boundary layer adjacent to an inclined flat plate subject to sudden heating and a temperature boundary condition which follows a ramp function up until a specified time and then remains constant is investigated. The development of the flow from start-up to a steady-state has been described based on scaling analyses and verified by numerical simulations. Different flow regimes based on the Rayleigh number are discussed with numerical results for both boundary conditions. For ramp heating, the boundary layer flow depends on the comparison of the time at which the ramp heating is completed and the time at which the boundary layer completes its growth. If the ramp time is long compared with the steady state time, the layer reaches a quasi steady mode in which the growth of the layer is governed solely by the thermal balance between convection and conduction. On the other hand, if the ramp is completed before the layer becomes steady; the subsequent growth is governed by the balance between buoyancy and inertia, as for the case of instantaneous heating.

A review of natural convection and heat transfer in attic-shaped space

This paper presents a review of studies on natural convection heat transfer in the triangular enclosure namely, in attic-shaped space. Much research activity has been devoted to this topic over the last three decades with a view to providing thermal comfort to the occupants in attic-shaped buildings and to minimising the energy costs associated with heating and air-conditioning. Two basic thermal boundary conditions of attic are considered to represent hot and cold climates or day and night time. This paper also reports on a significant number of studies which have been performed recently on other topics related to the attic space, for example, attics subject to localized heating and attics filled with porous media.

Natural frequencies of circumferentially truncated sector plates with simply supported straight edges

The classical Rayleigh-Ritz method in conjunction with suitable co-ordinate transformations is found to be effective for accurate estimation of natural frequencies of circumferentially truncated circular sector plates with simply supported straight edges. Numerical results are obtained for all the nine combinations of clamped, simply supported and free boundary conditions at the circular edges and presented in the form of graphs. The analysis confirms an earlier observation that the plate behaves like a long rectangular strip as the width of the plate in the radial direction becomes small.

Natural frequencies of circumferentially truncated sector plates with simply supported straight edges

The classical Rayleigh-Ritz method in conjunction with suitable co-ordinate transformations is found to be effective for accurate estimation of natural frequencies of circumferentially truncated circular sector plates with simply supported straight edges. Numerical results are obtained for all the nine combinations of clamped, simply supported and free boundary conditions at the circular edges and presented in the form of graphs. The analysis confirms an earlier observation that the plate behaves like a long rectangular strip as the width of the plate in the radial direction becomes small.

Analysis of continuum using boundary compatibility conditions of integrated force method

Static and vibration problems of an indeterminate continuum are traditionally analyzed by the stiffness method. The force method is more or less non-existent for such problems. This situation is primarily due to the incomplete state of development of the compatibility conditions which are essential for the analysis of indeterminate structures by the flexibility method. The understanding of the Compatibility Conditions (CC) has been substantially augmented. Based on the understanding of CC, a novel formulation termed the Integrated Force Method (IFM) has been established. In this paper IFM has been extended for the static and vibration analyses of a continuum. The IFM analysis is illustrated taking three examples: 1. (1) rectangular plate in flexure 2. (2) analysis of a cantilevered dam 3. (3) free vibration analysis of a beam. From the examples solved it is observed that the force response of an indeterminate continuum with mixed boundary conditions can be generated by IFM without any reference to displacements in the field or on the boundary. Displacements if required can be calculated by back substitution.

Boundary and internal conditions for adjoint fluid-flow problems

The ever-increasing robustness and reliability of flow-simulation methods have consolidated CFD as a major tool in virtually all branches of fluid mechanics. Traditionally, those methods have played a crucial role in the analysis of flow physics. In more recent years, though, the subject has broadened considerably, with the development of optimization and inverse design applications. Since then, the search for efficient ways to evaluate flow-sensitivity gradients has received the attention of numerous researchers. In this scenario, the adjoint method has emerged as, quite possibly, the most powerful tool for the job, which heightens the need for a clear understanding of its conceptual basis. Yet, some of its underlying aspects are still subject to debate in the literature, despite all the research that has been carried out on the method. Such is the case with the adjoint boundary and internal conditions, in particular. The present work aims to shed more light on that topic, with emphasis on the need for an internal shock condition. By following the path of previous authors, the quasi-1D Euler problem is used as a vehicle to explore those concepts. The results clearly indicate that the behavior of the adjoint solution through a shock wave ultimately depends upon the nature of the objective functional.

Dynamic stationary response of reinforced plates by the boundary element method

A direct version of the boundary element method (BEM) is developed to model the stationary dynamic response of reinforced plate structures, such as reinforced panels in buildings, automobiles, and airplanes. The dynamic stationary fundamental solutions of thin plates and plane stress state are used to transform the governing partial differential equations into boundary integral equations (BIEs). Two sets of uncoupled BIEs are formulated, respectively, for the in-plane state ( membrane) and for the out-of-plane state ( bending). These uncoupled systems are joined to formamacro-element, in which membrane and bending effects are present. The association of these macro-elements is able to simulate thin-walled structures, including reinforced plate structures. In the present formulation, the BIE is discretized by continuous and/or discontinuous linear elements. Four displacement integral equations are written for every boundary node. Modal data, that is, natural frequencies and the corresponding mode shapes of reinforced plates, are obtained from information contained in the frequency response functions (FRFs). A specific example is presented to illustrate the versatility of the proposed methodology. Different configurations of the reinforcements are used to simulate simply supported and clamped boundary conditions for the plate structures. The procedure is validated by comparison with results determined by the finite element method (FEM).

CH and C2 radicals characterization in natural gas turbulent diffusion flames

This paper reports the construction of an axisymmetric nonpremixed piloted jet burner, with well-defined initial and boundary conditions, known as the Delft burner, to assess turbulence-chemistry interaction in non-premixed turbulent flames. Detailed experimental information is described, involving hot-wire anemometry, thin-wire thermocouples and chemiluminescence visualization measurements. Radial profile of the axial mean velocity indicates excellent agreement between flow patterns developed within Delft installation and the one described herein. Chemiluminescence emissions from CH and C2 free-radicals were acquired with a CCD camera. Tomography reconstruction analysis was utilised to compare radical emissions and temperature spatial distributions. There was a strong dependence between temperature and CH/C 2 emissions. This is an indication that these radicals can be used in flame front studies.