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.

High-Pressure Behavior and Phase Stability of Al5BO9, a Mullite-Type Ceramic Material

Phase stability, elastic behavior, and pressure-induced structural evolution of synthetic boron-mullite Al5BO9 (a = 5.6780(7), b = 15.035(6), and c =7.698(3) Å, space group Cmc21, Z = 4) were investigated up to 25.6(1) GPa by in situ single-crystal synchrotron X-ray diffraction with a diamond anvil cell (DAC) under hydrostatic conditions. No evidence of phase transition was observed up to 21.7(1) GPa. At 25.6(1) GPa, the refined unit-cell parameters deviated significantly from the compressional trend, and the diffraction peaks appeared broader than at lower pressure. At 26.7(1) GPa, the diffraction pattern was not indexable, suggesting amorphization of the material or a phase transition to a high-pressure polymorph. Fitting the P–V data up to 21.7(1) GPa with a second-order Birch–Murnaghan Equation-of-State, we obtained a bulk modulus KT0 = 164(1) GPa. The axial compressibilities, here described as linearized bulk moduli, are as follows: KT0(a) = 244(9), KT0(b) = 120(4), and KT0(c) = 166(11) GPa (KT0(a):KT0(b):KT0(c) = 2.03:1:1.38). The structure refinements allowed a description of the main deformation mechanisms in response to the applied pressure. The stiffer crystallographic direction appears to be controlled by the infinite chains of edge-sharing octahedra running along [100], making the structure less compressible along the a-axis than along the b- and c-axis.

Dynamic high-resolution computer simulation of isotachophoretic enantiomer separation and zone stability

The development of electrophoretic computer models and their use for simulation of electrophoretic processes has increased significantly during the last few years. Recently, GENTRANS and SIMUL5 were extended with algorithms that describe chemical equilibria between solutes and a buffer additive in a fast 1:1 interaction process, an approach that enables simulation of the electrophoretic separation of enantiomers. For acidic cationic systems with sodium and H3 0(+) as leading and terminating components, respectively, acetic acid as counter component, charged weak bases as samples, and a neutral CD as chiral selector, the new codes were used to investigate the dynamics of isotachophoretic adjustment of enantiomers, enantiomer separation, boundaries between enantiomers and between an enantiomer and a buffer constituent of like charge, and zone stability. The impact of leader pH, selector concentration, free mobility of the weak base, mobilities of the formed complexes and complexation constants could thereby be elucidated. For selected examples with methadone enantiomers as analytes and (2-hydroxypropyl)-β-CD as selector, simulated zone patterns were found to compare well with those monitored experimentally in capillary setups with two conductivity detectors or an absorbance and a conductivity detector. Simulation represents an elegant way to provide insight into the formation of isotachophoretic boundaries and zone stability in presence of complexation equilibria in a hitherto inaccessible way.

Comparison of theoretical fixation stability of three devices employed in medial opening wedge high tibial osteotomy: a finite element analysis.

BACKGROUND Medial open wedge high tibial osteotomy is a well-established procedure for the treatment of unicompartmental osteoarthritis and symptomatic varus malalignment. We hypothesized that different fixation devices generate different fixation stability profiles for the various wedge sizes in a finite element (FE) analysis. METHODS Four types of fixation were compared: 1) first and 2) second generation Puddu plates, and 3) TomoFix plate with and 4) without bone graft. Cortical and cancellous bone was modelled and five different opening wedge sizes were studied for each model. Outcome measures included: 1) stresses in bone, 2) relative displacement of the proximal and distal tibial fragments, 3) stresses in the plates, 4) stresses on the upper and lower screw surfaces in the screw channels. RESULTS The highest load for all fixation types occurred in the plate axis. For the vast majority of the wedge sizes and fixation types the shear stress (von Mises stress) was dominating in the bone independent of fixation type. The relative displacements of the tibial fragments were low (in μm range). With an increasing wedge size this displacement tended to increase for both Puddu plates and the TomoFix plate with bone graft. For the TomoFix plate without bone graft a rather opposite trend was observed.For all fixation types the occurring stresses at the screw-bone contact areas pulled at the screws and exceeded the allowable threshold of 1.2 MPa for at least one screw surface. Of the six screw surfaces that were studied, the TomoFix plate with bone graft showed a stress excess of one out of twelve and without bone graft, five out of twelve. With the Puddu plates, an excess stress occurred in the majority of screw surfaces. CONCLUSIONS The different fixation devices generate different fixation stability profiles for different opening wedge sizes. Based on the computational simulations, none of the studied osteosynthesis fixation types warranted an intransigent full weight bearing per se. The highest fixation stability was observed for the TomoFix plates and the lowest for the first generation Puddu plate. These findings were revealed in theoretical models and need to be validated in controlled clinical settings.

Storage stability of SBS/sulfur modified bitumens at high temperature: influence of bitumen composition and structure

Polymer modified bitumens, PMBs, are usually prepared at high temperature and subsequently stored for a period of time, also at high temperature. The stability of PMBs, in these conditions, has a decisive influence in order to obtain the adequate performances for practical applications. In this article the attention is focused in the analysis of the factors that determine the stability of styrene–butadiene–styrene copolymer (SBS)/sulfur modified bitumens when the mixtures are maintained at high temperature. Bitumens from different crude oil sources were used to prepare SBS/sulfur modified bitumens. Changes in the values of viscosity, softening point, as well as in the morphology of PMB samples, stored at 160 °C, were related to the bitumen chemical composition and to the amount of asphaltene micelles present in the neat bitumen used in their preparation El trabajo se centra en el estudio de la influencia de la estructura /composición del betún sobre la compatibilidad del sistema betún/SBS. Cuatro betunes provenientes de dos crudos distintos se seleccionaron y sus mezclas se utilizaron para preparar betunes modificados con contenidos de SBS del 3% en peso

Global Stability Analysis of a Compressible Turbulent Flow around a High-Lift Configuration

The purpose of this work is to analyze a complex high lift configuration for which significant regions of separated flow are present. Current state of the art methods have some diffculty to predict the origin and the progression of this separated flow when increasing the angle of attack. The mechanisms responsible for the maximum lift limit on multi-element wing con?gurations are not clear; this stability analysis could help to understand the physics behind the phenomenon and to find a relation between the flow separation and the instability onset. The methodology presented herein consists in the computation of a steady base flow solution based on a finite volume discretization and a proposal of the solution for a generalized eigenvalue problem corresponding to the perturbed and linearized problem. The eigenvalue problem has been solved with the Arnoldi iterative method, one of the Krylov subspace projection methods. The described methodology was applied to the NACA0012 test case in subsonic and in transonic conditions and, finally, for the first time to the authors knowledge, on an industrial multi-component geometry, such as the A310 airfoil, in order to identify low frequency instabilities related to the separation. One important conclusion is that for all the analyzed geometries, one unstable mode related to flow separation appears for an angle of attack greater than the one correspondent to the maximum lift coe?cient condition. Finally, an adjoint study was carried out in order to evaluate the receptivity and the structural sensitivity of the geometries, giving an indication of the domain region that could be modified resulting in the biggest change of the flowfield.

A quantitative, high-throughput screen for protein stability

In proteomic research, it is often necessary to screen a large number of polypeptides for the presence of stable structure. Described here is a technique (referred to as SUPREX, stability of unpurified proteins from rates of H/D exchange) for measuring the stability of proteins in a rapid, high-throughput fashion. The method uses hydrogen exchange to estimate the stability of microgram quantities of unpurified protein extracts by using matrix-assisted laser desorption/ionization MS. The stabilities of maltose binding protein and monomeric λ repressor variants determined by SUPREX agree well with stability data obtained from conventional CD denaturation of purified protein. The method also can detect the change in stability caused by the binding of maltose to maltose binding protein. The results demonstrate the precision of the method over a wide range of stabilities.

Role of aromatic residues in stabilization of the [Fe4S4] cluster in high-potential iron proteins (HiPIPs): physical characterization and stability studies of Tyr-19 mutants of Chromatium vinosum HiPIP.

The functional role of residue Tyr-19 of Chromatium vinosum HiPIP has been evaluated by site-directed mutagenesis experiments. The stability of the [Fe4S4] cluster prosthetic center is sensitive to side-chain replacements. Polar residues result in significant instability, while nonpolar residues (especially with aromatic side chains) maintain cluster stability. Two-dimensional NMR data of native and mutant HiPIPs are consistent with a model where Tyr-19 serves to preserve the structural rigidity of the polypeptide backbone, thereby maintaining a hydrophobic barrier for exclusion of water from the cluster cavity. Solvent accessibility results in more facile oxidation of the cluster by atmospheric oxygen, with subsequent rapid hydrolysis of the [Fe4S4]3+ core.

The high-bias stability of monatomic chains

For the metals Au, Pt and Ir it is possible to form freely suspended monatomic chains between bulk electrodes. The atomic chains sustain very large current densities, but finally fail at high bias. We investigate the breaking mechanism, that involves current-induced heating of the atomic wires and electromigration forces. We find good agreement of the observations for Au based on models due to Todorov and co-workers. The high-bias breaking of atomic chains for Pt can also be described by the models, although here the parameters have not been obtained independently. In the limit of long chains the breaking voltage decreases inversely proportional to the length.

Measurement and stability of FTY720 in human whole blood by high-performance liquid chromatography-atmospheric pressure chemical ionization-tandem mass spectrometry

We report here a validated method for the quantification of a new immunosuppressant drug FTY720, using HPLC-tandem mass spectrometry. Whole blood samples (500 mu l) were subjected to liquid-liquid extraction, in the presence of an internal standard (Y-32919). Mass spectrometric detection was by selected reaction monitoring with an atmospheric pressure chemical ionization source in positive ionization mode (FTY720: m/z 308.3 -> 255.3). The assay was linear from 0.2 to 25 mu g/l (r(2) > 0.997, n = 5). The inter- and intra-day analytical recovery and imprecision for quality control samples (0.5, 7 and 15 mu g/l) were 95.8-103.2 and < 5.5%, respectively. At the lower limit of quantification (0.2 mu g/l) the interand intra-day analytical recovery was 99.0-102.8% with imprecision of < 7.6% (n = 5). The assay had a mean relative recovery of 100.5 +/- 5.8% (n = 15). Extracted samples were stable for 16 h. IFTY720 quality control samples were stable at room temperature for 16 h at 4 degrees C for at least 8 days and when taken through at least three freeze-thaw cycles. In conclusion, the method described displays analytical performance characteristics that are suitable for pharmacokinetic studies in humans. (c) 2006 Elsevier B.V. All rights reserved.