Numerical simulation of some viscoelastic fluids

Numerical simulation of the Oldroyd-B type viscoelastic fluids is a very challenging problem. rnThe well-known High Weissenberg Number Problem" has haunted the mathematicians, computer scientists, and rnengineers for more than 40 years. rnWhen the Weissenberg number, which represents the ratio of elasticity to viscosity, rnexceeds some limits, simulations done by standard methods break down exponentially fast in time. rnHowever, some approaches, such as the logarithm transformation technique can significantly improve rnthe limits of the Weissenberg number until which the simulations stay stable. rnrnWe should point out that the global existence of weak solutions for the Oldroyd-B model is still open. rnLet us note that in the evolution equation of the elastic stress tensor the terms describing diffusive rneffects are typically neglected in the modelling due to their smallness. However, when keeping rnthese diffusive terms in the constitutive law the global existence of weak solutions in two-space dimension rncan been shown. rnrnThis main part of the thesis is devoted to the stability study of the Oldroyd-B viscoelastic model. rnFirstly, we show that the free energy of the diffusive Oldroyd-B model as well as its rnlogarithm transformation are dissipative in time. rnFurther, we have developed free energy dissipative schemes based on the characteristic finite element and finite difference framework. rnIn addition, the global linear stability analysis of the diffusive Oldroyd-B model has also be discussed. rnThe next part of the thesis deals with the error estimates of the combined finite element rnand finite volume discretization of a special Oldroyd-B model which covers the limiting rncase of Weissenberg number going to infinity. Theoretical results are confirmed by a series of numerical rnexperiments, which are presented in the thesis, too.

Geometrical optimization of a short-stem hip implant for the reduction of proximal stress shielding (ottimizzazione geometrica di uno stelo corto di protesi d?anca per la riduzione della schermatura del carico prossimale)

Nowadays the number of hip joints arthroplasty operations continues to increase because the elderly population is growing. Moreover, the global life expectancy is increasing and people adopt a more active way of life. For this reasons, the demand of implant revision operations is becoming more frequent. The operation procedure includes the surgical removal of the old implant and its substitution with a new one. Every time a new implant is inserted, it generates an alteration in the internal femur strain distribution, jeopardizing the remodeling process with the possibility of bone tissue loss. This is of major concern, particularly in the proximal Gruen zones, which are considered critical for implant stability and longevity. Today, different implant designs exist in the market; however there is not a clear understanding of which are the best implant design parameters to achieve mechanical optimal conditions. The aim of the study is to investigate the stress shielding effect generated by different implant design parameters on proximal femur, evaluating which ranges of those parameters lead to the most physiological conditions.

The contribution of the Finite Volume Point Dilution Method to the determination of solute mass flux in an alluvial aquifer.

Groundwater represents one of the most important resources of the world and it is essential to prevent its pollution and to consider remediation intervention in case of contamination. According to the scientific community the characterization and the management of the contaminated sites have to be performed in terms of contaminant fluxes and considering their spatial and temporal evolution. One of the most suitable approach to determine the spatial distribution of pollutant and to quantify contaminant fluxes in groundwater is using control panels. The determination of contaminant mass flux, requires measurement of contaminant concentration in the moving phase (water) and velocity/flux of the groundwater. In this Master Thesis a new solute flux mass measurement approach, based on an integrated control panel type methodology combined with the Finite Volume Point Dilution Method (FVPDM), for the monitoring of transient groundwater fluxes, is proposed. Moreover a new adsorption passive sampler, which allow to capture the variation of solute concentration with time, is designed. The present work contributes to the development of this approach on three key points. First, the ability of the FVPDM to monitor transient groundwater fluxes was verified during a step drawdown test at the experimental site of Hermalle Sous Argentau (Belgium). The results showed that this method can be used, with optimal results, to follow transient groundwater fluxes. Moreover, it resulted that performing FVPDM, in several piezometers, during a pumping test allows to determine the different flow rates and flow regimes that can occurs in the various parts of an aquifer. The second field test aiming to determine the representativity of a control panel for measuring mass flus in groundwater underlined that wrong evaluations of Darcy fluxes and discharge surfaces can determine an incorrect estimation of mass fluxes and that this technique has to be used with precaution. Thus, a detailed geological and hydrogeological characterization must be conducted, before applying this technique. Finally, the third outcome of this work concerned laboratory experiments. The test conducted on several type of adsorption material (Oasis HLB cartridge, TDS-ORGANOSORB 10 and TDS-ORGANOSORB 10-AA), in order to determine the optimum medium to dimension the passive sampler, highlighted the necessity to find a material with a reversible adsorption tendency to completely satisfy the request of the new passive sampling technique.

Prediction of dental implant torque with a fast and automatic finite element analysis: a pilot study

Despite its importance, implant removal torque can be assessed at present only after implantation. This paper presents a new technique to help clinicians preoperatively evaluate implant stability.

On the terminal homologation of physiologically active peptides as a means of increasing stability in human serum--neurotensin, opiorphin, B27-KK10 epitope, NPY

The terminal homologation by CH(2) insertion into the peptides mentioned in the title is described. This involves replacement of the N-terminal amino acid residue by a β(2) - and of the C-terminal amino acid residue by a β(3) -homo-amino acid moiety (β(2) hXaa and β(3) hXaa, resp.; Fig. 1). In this way, the structure of the peptide chain from the N-terminal to the C-terminal stereogenic center is identical, and the modified peptide is protected against cleavage by exopeptidases (Figs. 2 and 3). Neurotensin (NT; 1) and its C-terminal fragment NT(8-13) are ligands of the G-protein-coupled receptors (GPCR) NT1, NT2, NT3, and NT analogs are promising tools to be used in cancer diagnostics and therapy. The affinities of homologated NT analogs, 2b-2e, for NT1 and NT2 receptors were determined by using cell homogenates and tumor tissues (Table 1); in the latter experiments, the affinities for the NT1 receptor are more or less the same as those of NT (0.5-1.3 vs. 0.6 nM). At the same time, one of the homologated NT analogs, 2c, survives in human plasma for 7 days at 37° (Fig. 6). An NMR analysis of NT(8-13) (Tables 2 and 4, and Fig. 8) reveals that this N-terminal NT fragment folds to a turn in CD(3) OH. - In the case of the human analgesic opiorphin (3a), a pentapeptide, and of the HIV-derived B27-KK10 (4a), a decapeptide, terminal homologation (→3b and 4b, resp.) led to a 7- and 70-fold half-life increase in plasma (Fig. 9). With N-terminally homologated NPY, 5c, we were not able to determine serum stability; the peptide consisting of 36 amino acid residues is subject to cleavage by endopetidases. Three of the homologated compounds, 2b, 2c, and 5c, were shown to be agonists (Fig. 7 and 11). A comparison of terminal homologation with other stability-increasing terminal modifications of peptides is performed (Fig. 5), and possible applications of the neurotensin analogs, described herein, are discussed.

Comparison of implant stability by means of resonance frequency analysis for flapless and conventionally inserted implants

Background: Resonance frequency analysis (RFA) is a noninvasive technique for the quantitative assessment of implant stability. Information on the implant stability quotient (ISQ) of transmucosally inserted implants is limited. Purpose: The aim of this investigation was to compare the ISQ of conventionally inserted implants by raising a muco-periostal flap with implants inserted using a flapless procedure. Materials and Methods: Forty elderly patients with complete edentulous maxilla were consecutively admitted for treatment with implant-supported prostheses. A computer tomography was obtained for the computer-assisted implant planning. One hundred ten implants were placed conventionally in 23 patients (flap-group) and 85 implants in 17 patients by means of the flapless method (flapless-group) using a stereolithographic template. RFA measurements were performed after implant placement (baseline) and after a healing time of 12 weeks (reentry). Results: All implants exhibited clinically and radiographically successful osseointegration. Bone level did not change significantly neither for genders nor type of surgical protocol. Mean ISQ values of the flapless-group were significantly higher at baseline (p < .001) and at reentry (p < .001) compared with the flap-group. The ISQ values were significantly lower at reentry compared with baseline for the flap-group (p = .028) but not for the flapless-group. This group showed a moderate, but insignificant increase. RFA measurements of males resulted in ISQ values that were thoroughly higher as compared with females at both time-points in both groups. Correlation between RFA and bone level was not found. Conclusions: The flapless procedure showed favorable conditions with regard to implant stability and crestal bone level. Some changes of the ISQ values that represent primary (mechanical) and secondary (bone remodeling) implant stability were observed in slight favor of the flapless method and male patients. In properly planned and well-selected cases, the minimal invasive transmucosal technique using a drill-guide is a safe procedure.

Uncovering the dynamics of cardiac systems using stochastic pacing and frequency domain analyses

Alternans of cardiac action potential duration (APD) is a well-known arrhythmogenic mechanism which results from dynamical instabilities. The propensity to alternans is classically investigated by examining APD restitution and by deriving APD restitution slopes as predictive markers. However, experiments have shown that such markers are not always accurate for the prediction of alternans. Using a mathematical ventricular cell model known to exhibit unstable dynamics of both membrane potential and Ca2+ cycling, we demonstrate that an accurate marker can be obtained by pacing at cycle lengths (CLs) varying randomly around a basic CL (BCL) and by evaluating the transfer function between the time series of CLs and APDs using an autoregressive-moving-average (ARMA) model. The first pole of this transfer function corresponds to the eigenvalue (λalt) of the dominant eigenmode of the cardiac system, which predicts that alternans occurs when λalt≤−1. For different BCLs, control values of λalt were obtained using eigenmode analysis and compared to the first pole of the transfer function estimated using ARMA model fitting in simulations of random pacing protocols. In all versions of the cell model, this pole provided an accurate estimation of λalt. Furthermore, during slow ramp decreases of BCL or simulated drug application, this approach predicted the onset of alternans by extrapolating the time course of the estimated λalt. In conclusion, stochastic pacing and ARMA model identification represents a novel approach to predict alternans without making any assumptions about its ionic mechanisms. It should therefore be applicable experimentally for any type of myocardial cell.

Evaluation of the Design of Komendant's Cycloidal, Thin Shell Roofs At the Kimbell Art Museum Using Finite Element Analysis

The analysis of Komendant's design of the Kimbell Art Museum was carried out in order to determine the effectiveness of the ring beams, edge beams and prestressing in the shells of the roof system. Finite element analysis was not available to Komendant or other engineers of the time to aid them in the design and analysis. Thus, the use of this tool helped to form a new perspective on the Kimbell Art Museum and analyze the engineer's work. In order to carry out the finite element analysis of Kimbell Art Museum, ADINA finite element analysis software was utilized. Eight finite element models (FEM-1 through FEM-8) of increasing complexity were created. The results of the most realistic model, FEM-8, which included ring beams, edge beams and prestressing, were compared to Komendant's calculations. The maximum deflection at the crown of the mid-span surface of -0.1739 in. in FEM-8 was found to be larger than Komendant's deflection in the design documents before the loss in prestressing force (-0.152 in.) but smaller than his prediction after the loss in prestressing force (-0.3814 in.). Komendant predicted a larger longitudinal stress of -903 psi at the crown (vs. -797 psi in FEM-8) and 37 psi at the edge (vs. -347 psi in FEM-8). Considering the strength of concrete of 5000 psi, the difference in results is not significant. From the analysis it was determined that both FEM-5, which included prestressing and fixed rings, and FEM-8 can be successfully and effectively implemented in practice. Prestressing was used in both models and thus served as the main contribution to efficiency. FEM-5 showed that ring and edge beams can be avoided, however an architect might find them more aesthetically appropriate than rigid walls.

Phosphatase and tensin homolog regulates stability and activity of EphB1 receptor

Deregulation of receptor tyrosine kinases (RTKs) is linked to a broad range of cancers, stressing the necessity of studying their regulatory pathways. We and others demonstrated previously that c-Cbl is necessary for the lysosomal degradation of erythropoietin-producing hepatocellular B1 (EphB1) carcinoma and epidermal growth factor receptor (EGFR) RTKs. Moreover, the tumor suppressor phosphatase and tensin homolog (PTEN) was shown to modulate c-Cbl-dependent EGFR degradation. We therefore investigated the involvement of PTEN in EphB1 signaling and degradation. We used PTEN mutants, PTEN, and NHERF1 small interfering RNA in CHO-EphB1 and SW480 cells endogenously expressing EphB1 to delineate EphB1-PTEN interactions. PTEN was constitutively associated with c-Cbl, protecting it from degradation. EphB1 stimulation triggered ∼50% serine-threonine PTEN dephosphorylation and PTEN-Cbl complex disruption, a process requiring PTEN protein phosphatase activity. Both proteins independently translocated to EphB1, with PTEN in association with the scaffold protein NHERF1. Biologically, PTEN lipid phosphatase activity impairs EphB1-dependent cell adhesion and chemotaxis. This study demonstrates for the first time in mammalian cells that the Eph receptor and PTEN associate and influence their signaling. Moreover, it contributes to the emerging concept that PTEN regulates expression of RTKs through modulation of their degradation. Finally, it reveals a new role for PTEN protein phosphatase activity involved in this process.

Eluding oblivion with smart stochastic selection of synaptic updates

The variables involved in the equations that describe realistic synaptic dynamics always vary in a limited range. Their boundedness makes the synapses forgetful, not for the mere passage of time, but because new experiences overwrite old memories. The forgetting rate depends on how many synapses are modified by each new experience: many changes means fast learning and fast forgetting, whereas few changes means slow learning and long memory retention. Reducing the average number of modified synapses can extend the memory span at the price of a reduced amount of information stored when a new experience is memorized. Every trick which allows to slow down the learning process in a smart way can improve the memory performance. We review some of the tricks that allow to elude fast forgetting (oblivion). They are based on the stochastic selection of the synapses whose modifications are actually consolidated following each new experience. In practice only a randomly selected, small fraction of the synapses eligible for an update are actually modified. This allows to acquire the amount of information necessary to retrieve the memory without compromising the retention of old experiences. The fraction of modified synapses can be further reduced in a smart way by changing synapses only when it is really necessary, i.e. when the post-synaptic neuron does not respond as desired. Finally we show that such a stochastic selection emerges naturally from spike driven synaptic dynamics which read noisy pre and post-synaptic neural activities. These activities can actually be generated by a chaotic system.

The long-term mechanical integrity of non-reinforced PEEK-OPTIMA polymer for demanding spinal applications: experimental and finite-element analysis

Polyetheretherketone (PEEK) is a novel polymer with potential advantages for its use in demanding orthopaedic applications (e.g. intervertebral cages). However, the influence of a physiological environment on the mechanical stability of PEEK has not been reported. Furthermore, the suitability of the polymer for use in highly stressed spinal implants such as intervertebral cages has not been investigated. Therefore, a combined experimental and analytical study was performed to address these open questions. A quasi-static mechanical compression test was performed to compare the initial mechanical properties of PEEK-OPTIMA polymer in a dry, room-temperature and in an aqueous, 37 degrees C environment (n=10 per group). The creep behaviour of cylindrical PEEK polymer specimens (n=6) was measured in a simulated physiological environment at an applied stress level of 10 MPa for a loading duration of 2000 hours (12 weeks). To compare the biomechanical performance of different intervertebral cage types made from PEEK and titanium under complex loading conditions, a three-dimensional finite element model of a functional spinal unit was created. The elastic modulus of PEEK polymer specimens in a physiological environment was 1.8% lower than that of specimens tested at dry, room temperature conditions (P<0.001). The results from the creep test showed an average creep strain of less than 0.1% after 2000 hours of loading. The finite element analysis demonstrated high strain and stress concentrations at the bone/implant interface, emphasizing the importance of cage geometry for load distribution. The stress and strain maxima in the implants were well below the material strength limits of PEEK. In summary, the experimental results verified the mechanical stability of the PEEK-OPTIMA polymer in a simulated physiological environment, and over extended loading periods. Finite element analysis supported the use of PEEK-OPTIMA for load-bearing intervertebral implants.

Statistical finite element model for bone shape and biomechanical properties

We present a framework for statistical finite element analysis combining shape and material properties, and allowing performing statistical statements of biomechanical performance across a given population. In this paper, we focus on the design of orthopaedic implants that fit a maximum percentage of the target population, both in terms of geometry and biomechanical stability. CT scans of the bone under consideration are registered non-rigidly to obtain correspondences in position and intensity between them. A statistical model of shape and intensity (bone density) is computed by means of principal component analysis. Afterwards, finite element analysis (FEA) is performed to analyse the biomechanical performance of the bones. Realistic forces are applied on the bones and the resulting displacement and bone stress distribution are calculated. The mechanical behaviour of different PCA bone instances is compared.

Significance of primary stability for osseointegration of dental implants

AIM: To investigate the significance of the initial stability of dental implants for the establishment of osseointegration in an experimental capsule model for bone augmentation. MATERIAL AND METHODS: Sixteen male rats were used in the study. In each rat, muscle-periosteal flaps were elevated on the lateral aspect of the mandibular ramus on both sides, resulting in exposure of the bone surface. Small perforations were then produced in the ramus. A rigid, hemispherical Teflon capsule with a diameter of 6 mm and a height of 4 mm and with a hole in its middle portion, prepared to fit the circumference of an ITI HC titanium implant of 2.8 mm in diameter, was fixed to the ramus using 4 mini-screws. On one side of the jaw, the implant was placed through the hole in such a way that its apex did not make contact with the mandibular ramus (test). This placement of the implant did not ensure primary stability. On the other side of the jaw, a similar implant was placed through the hole of the capsule in such a way that contact was made between the implant and the surface of the ramus (control). This provided primary stability of the implant. After placement of the implants, the soft tissues were repositioned over the capsules and sutured. After 1, 3, 6 and 9 months, four animals were sacrificed and subjected to histometric analysis. RESULTS: The mean height of direct bone-to-implant contact of implants with primary stability was 38.8%, 52.9%, 64.6% and 81.3% of the implant length at 1, 3, 6 and 9 months, respectively. Of the bone adjacent to the implant surface, 28.1%, 28.9%, 52.6% and 69.6%, respectively, consisted of mineralized bone. At the test implants, no bone-to-implant contact was observed at any observation time or in any of these non-stabilized specimens. CONCLUSION: The findings of the present study indicate that primary implant stability is a prerequisite for successful osseointegration, and that implant instability results in fibrous encapsulation, thus confirming previously made clinical observations.

On the Accelerated Failure Time Model for Current Status and Interval Censored Data

This paper introduces a novel approach to making inference about the regression parameters in the accelerated failure time (AFT) model for current status and interval censored data. The estimator is constructed by inverting a Wald type test for testing a null proportional hazards model. A numerically efficient Markov chain Monte Carlo (MCMC) based resampling method is proposed to simultaneously obtain the point estimator and a consistent estimator of its variance-covariance matrix. We illustrate our approach with interval censored data sets from two clinical studies. Extensive numerical studies are conducted to evaluate the finite sample performance of the new estimators.

On Time Series Analysis of Public Health and Biomedical Data

A time series is a sequence of observations made over time. Examples in public health include daily ozone concentrations, weekly admissions to an emergency department or annual expenditures on health care in the United States. Time series models are used to describe the dependence of the response at each time on predictor variables including covariates and possibly previous values in the series. Time series methods are necessary to account for the correlation among repeated responses over time. This paper gives an overview of time series ideas and methods used in public health research.