Numerical modelling of fluid-structure interaction with application in hemodynamics

The thesis deals with numerical algorithms for fluid-structure interaction problems with application in blood flow modelling. It starts with a short introduction on the mathematical description of incompressible viscous flow with non-Newtonian viscosity and a moving linear viscoelastic structure. The mathematical model consists of the generalized Navier-Stokes equation used for the description of fluid flow and the generalized string model for structure movement. The arbitrary Lagrangian-Eulerian approach is used in order to take into account moving computational domain. A part of the thesis is devoted to the discussion on the non-Newtonian behaviour of shear-thinning fluids, which is in our case blood, and derivation of two non-Newtonian models frequently used in the blood flow modelling. Further we give a brief overview on recent fluid-structure interaction schemes with discussion about the difficulties arising in numerical modelling of blood flow. Our main contribution lies in numerical and experimental study of a new loosely-coupled partitioned scheme called the kinematic splitting fluid-structure interaction algorithm. We present stability analysis for a coupled problem of non-Newtonian shear-dependent fluids in moving domains with viscoelastic boundaries. Here, we assume both, the nonlinearity in convective as well is diffusive term. We analyse the convergence of proposed numerical scheme for a simplified fluid model of the Oseen type. Moreover, we present series of experiments including numerical error analysis, comparison of hemodynamic parameters for the Newtonian and non-Newtonian fluids and comparison of several physiologically relevant computational geometries in terms of wall displacement and wall shear stress. Numerical analysis and extensive experimental study for several standard geometries confirm reliability and accuracy of the proposed kinematic splitting scheme in order to approximate fluid-structure interaction problems.

Study of the effect of minor compounds on the oxidative stability and physical properties of bulk oil, oil-in-water emulsion, and food emulsions

Minor components are of particular interest due to their antioxidant and biological properties. Various classes of lipophilic minor components (plant sterols (PS) and α-tocopherol) were selected as they are widely used in the food industry. A Fast GC-MS method for PS analysis in functional dairy products was set up. The analytical performance and significant reduction of the analysis time and consumables, demonstrated that Fast GC-MS could be suitable for the PS analysis in functional dairy products. Due to their chemical structure, PS can undergo oxidation, which could be greatly impacted by matrix nature/composition and thermal treatments. The oxidative stability of PS during microwave heating was evaluated. Two different model systems (PS alone and in combination) were heated up to 30 min at 1000 W. PS degraded faster when they were alone than in presence of TAG. The extent of PS degradation depends on both heating time and the surrounding medium, which can impact the quality and safety of the food product destined to microwave heating/cooking. Many minor lipid components are included in emulsion systems and can affect the rate of lipid oxidation. The oxidative stability of oil-in-water (O/W) emulsions containing PS esters, ω-3 FA and phenolic compounds, were evaluated after a 14-day storage at room temperature. Due to their surface active character, PS could be particularly prone to oxidation when they are incorporated in emulsions, as they are more exposed to water-soluble prooxidants. Finally, some minor lipophilic components may increase oxidative stability of food systems due to their antioxidant activity. á-tocopherol partitioning and antioxidant activity was determined in the presence of excess SDS in stripped soybean O/W emulsions. Results showed that surfactant micelles could play a key role as an antioxidant carrier, by potentially increasing the accessibility of hydrophobic antioxidant to the interface.

Stability analysis of linear thin shells

Shell structure is widely used in engineering area. The purpose of this dissertation is to show the behavior of a thin shell under external load, especially for long cylindrical shell under compressive load, I analyzed not only for linear elastic problem and also for buckling problem, and by using finite element analysis it shows that the imperfection of a cylinder could affect the critical load which means the buckling capability of this cylinder. For linear elastic problem, I compared the theoretical results with the results got from Straus7 and Abaqus, and the results are really close. For the buckling problem I did the same: compared the theoretical and Abaqus results, the error is less than 1%, but in reality, it’s not possible to reach the theoretical buckling capability due to the imperfection of the cylinder, so I put different imperfection for the cylinder in Abaqus, and found out that with the increasing of the percentage of imperfection, the buckling capability decreases, for example 10% imperfection could decrease 40% of the buckling capability, and the outcome meet the buckling behavior in reality.

Assessment of the population structure and temporal changes in spatial dynamics and genetic characteristics of the Atlantic bluefin tuna under a fishery independent framework.

As a large and long-lived species with high economic value, restricted spawning areas and short spawning periods, the Atlantic bluefin tuna (BFT; Thunnus thynnus) is particularly susceptible to over-exploitation. Although BFT have been targeted by fisheries in the Mediterranean Sea for thousands of years, it has only been in these last decades that the exploitation rate has reached far beyond sustainable levels. An understanding of the population structure, spatial dynamics, exploitation rates and the environmental variables that affect BFT is crucial for the conservation of the species. The aims of this PhD project were 1) to assess the accuracy of larval identification methods, 2) determine the genetic structure of modern BFT populations, 3) assess the self-recruitment rate in the Gulf of Mexico and Mediterranean spawning areas, 4) estimate the immigration rate of BFT to feeding aggregations from the various spawning areas, and 5) develop tools capable of investigating the temporal stability of population structuring in the Mediterranean Sea. Several weaknesses in modern morphology-based taxonomy including demographic decline of expert taxonomists, flawed identification keys, reluctance of the taxonomic community to embrace advances in digital communications and a general scarcity of modern user-friendly materials are reviewed. Barcoding of scombrid larvae revealed important differences in the accuracy of the taxonomic identifications carried out by different ichthyoplanktologists following morphology-based methods. Using a Genotyping-by-Sequencing a panel of 95 SNPs was developed and used to characterize the population structuring of BFT and composition of adult feeding aggregations. Using novel molecular techniques, DNA was extracted from bluefin tuna vertebrae excavated from late iron age, ancient roman settlements Byzantine-era Constantinople and a 20th century collection. A second panel of 96 SNPs was developed to genotype historical and modern samples in order to elucidate changes in population structuring and allele frequencies of loci associated with selective traits.

Numerical analysis of the relation between interactions and structure in a molecular fluid

Coarse graining is a popular technique used in physics to speed up the computer simulation of molecular fluids. An essential part of this technique is a method that solves the inverse problem of determining the interaction potential or its parameters from the given structural data. Due to discrepancies between model and reality, the potential is not unique, such that stability of such method and its convergence to a meaningful solution are issues.rnrnIn this work, we investigate empirically whether coarse graining can be improved by applying the theory of inverse problems from applied mathematics. In particular, we use the singular value analysis to reveal the weak interaction parameters, that have a negligible influence on the structure of the fluid and which cause non-uniqueness of the solution. Further, we apply a regularizing Levenberg-Marquardt method, which is stable against the mentioned discrepancies. Then, we compare it to the existing physical methods - the Iterative Boltzmann Inversion and the Inverse Monte Carlo method, which are fast and well adapted to the problem, but sometimes have convergence problems.rnrnFrom analysis of the Iterative Boltzmann Inversion, we elaborate a meaningful approximation of the structure and use it to derive a modification of the Levenberg-Marquardt method. We engage the latter for reconstruction of the interaction parameters from experimental data for liquid argon and nitrogen. We show that the modified method is stable, convergent and fast. Further, the singular value analysis of the structure and its approximation allows to determine the crucial interaction parameters, that is, to simplify the modeling of interactions. Therefore, our results build a rigorous bridge between the inverse problem from physics and the powerful solution tools from mathematics. rn

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.

Factor Structure of the Neck Disability Index

BACKGROUND CONTEXT: The Neck Disability Index frequently is used to measure outcomes of the neck. The statistical rigor of the Neck Disability Index has been assessed with conflicting outcomes. To date, Confirmatory Factor Analysis of the Neck Disability Index has not been reported for a suitably large population study. Because the Neck Disability Index is not a condition-specific measure of neck function, initial Confirmatory Factor Analysis should consider problematic neck patients as a homogenous group. PURPOSE: We sought to analyze the factor structure of the Neck Disability Index through Confirmatory Factor Analysis in a symptomatic, homogeneous, neck population, with respect to pooled populations and gender subgroups. STUDY DESIGN: This was a secondary analysis of pooled data. PATIENT SAMPLE: A total of 1,278 symptomatic neck patients (67.5% female, median age 41 years), 803 nonspecific and 475 with whiplash-associated disorder. OUTCOME MEASURES: The Neck Disability Index was used to measure outcomes. METHODS: We analyzed pooled baseline data from six independent studies of patients with neck problems who completed Neck Disability Index questionnaires at baseline. The Confirmatory Factor Analysis was considered in three scenarios: the full sample and separate sexes. Models were compared empirically for best fit. RESULTS: Two-factor models have good psychometric properties across both the pooled and sex subgroups. However, according to these analyses, the one-factor solution is preferable from both a statistical perspective and parsimony. The two-factor model was close to significant for the male subgroup (p<.07) where questions separated into constructs of mental function (pain, reading headaches and concentration) and physical function (personal care, lifting, work, driving, sleep, and recreation). CONCLUSIONS: The Neck Disability Index demonstrated a one-factor structure when analyzed by Confirmatory Factor Analysis in a pooled, homogenous sample of neck problem patients. However, a two-factor model did approach significance for male subjects where questions separated into constructs of mental and physical function. Further investigations in different conditions, subgroup and sex-specific populations are warranted.

Structure and thermodynamics of H3O+(H2O)8 clusters: A combined molecular dynamics and quantum mechanics approach

We have studied the structure and stability of (H3O+)(H2O)8 clusters using a combination of molecular dynamics sampling and high-level ab initio calculations. 20 distinct oxygen frameworks are found within 2 kcal/mol of the electronic or standard Gibbs free energy minimum. The impact of quantum zero-point vibrational corrections on the relative stability of these isomers is quite significant. The box-like isomers are favored in terms of electronic energy, but with the inclusion of zero-point vibrational corrections and entropic effects tree-like isomers are favored at higher temperatures. Under conditions from 0 to 298.15 K, the global minimum is predicted to be a tree-like structure with one dangling singly coordinated water molecule. Above 298.15 K, higher entropy tree-like isomers with two or more singly coordinated water molecules are favored. These assignments are generally consistent with experimental IR spectra of (H3O+)(H2O)8 obtained at 150 K.

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.

Total synthesis of (-)-zampanolide and structure-activity relationship studies on (-)-dactylolide derivatives

A new total synthesis of the marine macrolide (-)-zampanolide (1) and the structurally and stereochemically related non-natural levorotatory enantiomer of (+)-dactylolide (2), that is, ent-2, has been developed. The synthesis features a high-yielding, selective intramolecular Horner-Wadsworth-Emmons (HWE) reaction to close the 20-membered macrolactone ring of 1 and ent-2. The β-keto phosphonate/aldehyde precursor for the ring-closure reaction was obtained by esterification of a ω-diethylphosphono carboxylic acid fragment and a secondary alcohol fragment incorporating the THP ring that is embedded in the macrocyclic core structure of 1 and ent-2. THP ring formation was accomplished through a segment coupling Prins-type cyclization. Employing the same overall strategy, 13-desmethylene-ent-2 as well as the monocyclic desTHP derivatives of 1 and ent-2 were prepared. Synthetic 1 inhibited human cancer cell growth in vitro with nM IC(50) values, while ent-2, which lacks the diene-containing hemiaminal-linked side chain of 1, is 25- to 260-fold less active. 13-Desmethylene-ent-2 as well as the reduced versions of ent-2 and 13-desmethylene-ent-2 all showed similar cellular activity as ent-2 itself. The same activity level was attained by the monocyclic desTHP derivative of 1. Oxidation of the aldehyde functionality of ent-2 gave a carboxylic acid that was converted into the corresponding N-hexyl amide. The latter showed only μM antiproliferative activity, thus being several hundred-fold less potent than 1.

Structure and properties of graphene nano disks (GND) with and without edge-dopants

Graphene is one of the most important materials. In this research, the structures and properties of graphene nano disks (GND) with a concentric shape were investigated by Density Functional Theory (DFT) calculations, in which the most effective DFT methods - B3lyp and Pw91pw91 were employed. It was found that there are two types of edges - Zigzag and Armchair in concentric graphene nano disks (GND). The bond length between armchair-edge carbons is much shorter than that between zigzag-edge carbons. For C24 GND that consists of 24 carbon atoms, only armchair edge with 12 atoms is formed. For a GND larger than the C24 GND, both armchair and zigzag edges co-exist. Furthermore, when the number of carbon atoms in armchair-edge are always 12, the number of zigzag-edge atoms increases with increasing the size of a GND. In addition, the stability of a GND is enhanced with increasing its size, because the ratio of edge-atoms to non-edge-atoms decreases. The size effect of a graphene nano disk on its HOMO-LUMO energy gap was evaluated. C6 and C24 GNDs possess HOMO-LUMO gaps of 1.7 and 2.1eV, respectively, indicating that they are semi-conductors. In contrast, C54 and C96 GNDs are organic metals, because their HOMO-LUMO gaps are as low as 0.3 eV. The effect of doping foreign atoms to the edges of GNDs on their structures, stabilities, and HOMO-LUMO energy gaps were also examined. When foreign atoms are attached to the edge of a GND, the original unsaturated carbon atoms become saturated. As a result, both of the C-C bonds lengths and the stability of a GND increase. Furthermore, the doping effect on the HOMO-LUMO energy gap is dependent on the type of doped atoms. The doping H, F, or OH into the edge of a GND increases its HOMO-LUMO energy gap. In contrast, a Li-doped GND has a lower HOMO-LUMO energy gap than that without doping. Therefore, Li-doping can increase the electrical conductance of a GND, whereas H, F, or OH-doping decreases its conductance.

Reconstruction of the palatal aponeurosis with autogenous fascia lata in secondary radical intravelar veloplasty: a new method

Velopharyngeal insufficiency in cleft patients with muscular insufficiency detected by nasendoscopy is commonly treated by secondary radical intravelar veloplasty, in which the palatal muscles are reoriented and positioned backwards. The dead space between the retro-displaced musculature and the posterior borders of the palatal bone remains problematic. Postoperatively, the surgically achieved lengthening of the soft palate often diminishes due to scar tissue formation in the dead space, leading to reattachment of the reoriented muscles to the palatal bone and to decreased mobility of the soft palate. To avoid this, the dead space should be restored by a structure imitating the function of the missing palatal aponeurosis. The entire dead space was covered using a double layer of autogenous fascia lata harvested from the lateral thigh, which should allow sufficient and permanent sliding of the retro-positioned musculature. A clinical case of a 9-year-old boy who underwent the operation is reported. Postoperatively, marked functional improvements were observable in speech assessment, nasendoscopy and nasometry. The case reported here suggests that the restoration of the dead space may be beneficial for effective secondary palatal repair. Fascia lata seems to be a suitable graft for this purpose.

Effect of high order interpolation in the stability and efficiency of the time-integration process in vorticity-velocity CFD algorithms

The numerical solution of the incompressible Navier-Stokes Equations offers an effective alternative to the experimental analysis of Fluid-Structure interaction i.e. dynamical coupling between a fluid and a solid which otherwise is very complex, time consuming and very expensive. To have a method which can accurately model these types of mechanical systems by numerical solutions becomes a great option, since these advantages are even more obvious when considering huge structures like bridges, high rise buildings, or even wind turbine blades with diameters as large as 200 meters. The modeling of such processes, however, involves complex multiphysics problems along with complex geometries. This thesis focuses on a novel vorticity-velocity formulation called the KLE to solve the incompressible Navier-stokes equations for such FSI problems. This scheme allows for the implementation of robust adaptive ODE time integration schemes and thus allows us to tackle the various multiphysics problems as separate modules. The current algorithm for KLE employs a structured or unstructured mesh for spatial discretization and it allows the use of a self-adaptive or fixed time step ODE solver while dealing with unsteady problems. This research deals with the analysis of the effects of the Courant-Friedrichs-Lewy (CFL) condition for KLE when applied to unsteady Stoke’s problem. The objective is to conduct a numerical analysis for stability and, hence, for convergence. Our results confirmthat the time step ∆t is constrained by the CFL-like condition ∆t ≤ const. hα, where h denotes the variable that represents spatial discretization.

Influence of Mechanical and Thermal Boundary Conditions on Stabilizing/Destabilizing Mechanisms in Evaporating Liquid Films

Liquid films, evaporating or non-evaporating, are ubiquitous in nature and technology. The dynamics of evaporating liquid films is a study applicable in several industries such as water recovery, heat exchangers, crystal growth, drug design etc. The theory describing the dynamics of liquid films crosses several fields such as engineering, mathematics, material science, biophysics and volcanology to name a few. Interfacial instabilities typically manifest by the undulation of an interface from a presumed flat state or by the onset of a secondary flow state from a primary quiescent state or both. To study the instabilities affecting liquid films, an evaporating/non-evaporating Newtonian liquid film is subject to a perturbation. Numerical analysis is conducted on configurations of such liquid films being heated on solid surfaces in order to examine the various stabilizing and destabilizing mechanisms that can cause the formation of different convective structures. These convective structures have implications towards heat transfer that occurs via this process. Certain aspects of this research topic have not received attention, as will be obvious from the literature review. Static, horizontal liquid films on solid surfaces are examined for their resistance to long wave type instabilities via linear stability analysis, method of normal modes and finite difference methods. The spatiotemporal evolution equation, available in literature, describing the time evolution of a liquid film heated on a solid surface, is utilized to analyze various stabilizing/destabilizing mechanisms affecting evaporating and non-evaporating liquid films. The impact of these mechanisms on the film stability and structure for both buoyant and non-buoyant films will be examined by the variation of mechanical and thermal boundary conditions. Films evaporating in zero gravity are studied using the evolution equation. It is found that films that are stable to long wave type instabilities in terrestrial gravity are prone to destabilization via long wave instabilities in zero gravity.

A reinvestigation into the reaction of NH4OAc with acetyl derivatives of Bay-lis-Hillman adducts: Formation of tertiary and secondary allyl amines in-stead of primary allyl amines

A reinvestigation into the reaction between ammonium acetate and the acetyl derivatives of Baylis-Hillman adducts has led us to conclude that the products obtained are tertiary and secondary allyl amines and not the primary allyl amines. The unambiguous assignment of the structure of products using chemical and spectroscopic methods is described