Finite element analysis of human joints

Our work focuses on the development of finite element models (FEMs) that describe the biomechanics of human joints. Finite element modeling is becoming a standard tool in industrial applications. In highly complex problems such as those found in biomechanics research, however, the full potential of FEMs is just beginning to be explored, due to the absence of precise, high resolution medical data and the difficulties encountered in converting these enormous datasets into a form that is usable in FEMs. With increasing computing speed and memory available, it is now feasible to address these challenges. We address the first by acquiring data with a high resolution C-ray CT scanner and the latter by developing semi-automated method for generating the volumetric meshes used in the FEM. Issues related to tomographic reconstruction, volume segmentation, the use of extracted surfaces to generate volumetric hexahedral meshes, and applications of the FEM are described

Displacements and rotations of a body moving about an arbitrary axis in a global reference frame

Measurement of human joint motion frequently involves the use of markers to describe joint motion in a global reference frame. Results may be quite arbitrary if the reference frame is not properly chosen with respect to the joint`s rotational axis(es). In nature joint axes can exist at any orientation and location relative to an arbitrarily chosen global reference frame. An arbitrary axis is any axis that is not coincident with a reference coordinate. Calculations are made of the errors that result when joint motion occurs about an arbitrary axis in a global reference frame

Automated volumetric grid generation for finite element modeling of human hand joints

We are developing techniques for finite element analysis of human joints. These techniques need to provide high quality results rapidly in order to be useful to a physician. The research presented here increases model quality and decreases user input time by automating the volumetric mesh generation step

The chromosphere: gateway to the corona, or the purgatory of solar physics?

I argue that one should attempt to understand the solar chromosphere not only for its own sake, but also if one is interested in the physics of: the corona; astrophysical dynamos; space weather; partially ionized plasmas; heliospheric UV radiation; the transition region. I outline curious observations which I personally find puzzling and deserving of attention.Comment: To appear in the proceedings of the 25th NSO Workshop "Chromospheric Structure and Dynamics. From Old Wisdom to New Insights", Memorie della Societa' Astronomica Italiana, Eds. Tritschler et a

The clustering of polarity reversals of the geomagnetic field

Often in nature the temporal distribution of inhomogeneous stochastic point processes can be modeled as a realization of renewal Poisson processes with a variable rate. Here we investigate one of the classical examples, namely the temporal distribution of polarity reversals of the geomagnetic field. In spite of the commonly used underlying hypothesis, we show that this process strongly departs from a Poisson statistics, the origin of this failure stemming from the presence of temporal clustering. We find that a Levy statistics is able to reproduce paleomagnetic data, thus suggesting the presence of long-range correlations in the underlying dynamo process.Comment: 4 pages, in press on PRL (31 march 2006?

Correlative geochemical study of crude oils from southeastern and southern parts of the Pannonian Basin

Several crude oils from the southeastern Yugoslav part of the Pannonian Basin are correlated for the first time on the basis of many bulk and specific parameters. The investigations involved 15 oils from oil-gas fields A and B and oil field C from the Banat depression and oil field D from South Bačka depression, and four oils from recently discovered oil shows E (the Kostolac depression). On the basis of structural group analysis and other bulk parameters such as API gravity, contents of asphaltenes and sulfur, as well as content of n-alkanes and pristane to phytane ratio, the examined oils are classified into four groups. However, the distributions and relative abundances of n-alkanes, isoprenoids (C19, C20), steranes, tricyclic terpanes and pentacyclic triterpanes, obtained by capillary GC and computerized GC-MS analysis, suggest classification of the examined oils into only two genetic types: A-D and E1-E4. The oils are of relatively high maturity. They are not biodegraded, except for both oils from field D and the oil A9. © 1987

SLIM Ultrahigh Resolution Ion Mobility Spectrometry Separations of Isotopologues and Isotopomers Reveal Mobility Shifts due to Mass Distribution Changes

We report on separations of ion isotopologues and isotopomers using ultrahigh-resolution traveling wave-based Structures for Lossless Ion Manipulations with serpentine ultralong path and extended routing ion mobility spectrometry coupled to mass spectrometry (SLIM SUPER IMS-MS). Mobility separations of ions from the naturally occurring ion isotopic envelopes (e.g., [M], [M+1], [M+2], ... ions) showed the first and second isotopic peaks (i.e., [M+1] and [M+2]) for various tetraalkylammonium ions could be resolved from their respective monoisotopic ion peak ([M]) after SLIM SUPER IMS with resolving powers of ∼400–600. Similar separations were obtained for other compounds (e.g., tetrapeptide ions). Greater separation was obtained using argon versus helium drift gas, as expected from the greater reduced mass contribution to ion mobility described by the Mason–Schamp relationship. To more directly explore the role of isotopic substitutions, we studied a mixture of specific isotopically substituted (15N, 13C, and 2H) protonated arginine isotopologues. While the separations in nitrogen were primarily due to their reduced mass differences, similar to the naturally occurring isotopologues, their separations in helium, where higher resolving powers could also be achieved, revealed distinct additional relative mobility shifts. These shifts appeared correlated, after correction for the reduced mass contribution, with changes in the ion center of mass due to the different locations of heavy atom substitutions. The origin of these apparent mass distribution-induced mobility shifts was then further explored using a mixture of Iodoacetyl Tandem Mass Tag (iodoTMT) isotopomers (i.e., each having the same exact mass, but with different isotopic substitution sites). Again, the observed mobility shifts appeared correlated with changes in the ion center of mass leading to multiple monoisotopic mobilities being observed for some isotopomers (up to a ∼0.04% difference in mobility). These mobility shifts thus appear to reflect details of the ion structure, derived from the changes due to ion rotation impacting collision frequency or momentum transfer, and highlight the potential for new approaches for ion structural characterization

Time-dependent probability density functions and information geometry in stochastic logistic and Gompertz models

A probabilistic description is essential for understanding growth processes in non-stationary states. In this paper, we compute time-dependent probability density functions (PDFs) in order to investigate stochastic logistic and Gompertz models, which are two of the most popular growth models. We consider different types of short-correlated multiplicative and additive noise sources and compare the time-dependent PDFs in the two models, elucidating the effects of the additive and multiplicative noises on the form of PDFs. We demonstrate an interesting transition from a unimodal to a bimodal PDF as the multiplicative noise increases for a fixed value of the additive noise. A much weaker (leaky) attractor in the Gompertz model leads to a significant (singular) growth of the population of a very small size. We point out the limitation of using stationary PDFs, mean value and variance in understanding statistical properties of the growth in non-stationary states, highlighting the importance of time-dependent PDFs. We further compare these two models from the perspective of information change that occurs during the growth process. Specifically, we define an infinitesimal distance at any time by comparing two PDFs at times infinitesimally apart and sum these distances in time. The total distance along the trajectory quantifies the total number of different states that the system undergoes in time, and is called the information length. We show that the time-evolution of the two models become more similar when measured in units of the information length and point out the merit of using the information length in unifying and understanding the dynamic evolution of different growth processes