Erythrocyte ghost cell-alkaline phosphatase: construction and characterization of a vesicular system for use in biomineralization studies

Alkaline phosphatase is required for the mineralization of bone and cartilage. This enzyme is localized in the matrix vesicle, which plays a role key in calcifying cartilage. In this paper we standardize a method to construction a resealed ghost cell-alkaline phosphatase system to mimic matrix vesicles and examine the kinetic behavior of the incorporated enzyme. Polidocanol-solubilized alkaline phosphatase, free of detergent, was incorporated into resealed ghost cells. This process was time-dependent and practically 50% of the enzyme was incorporated into the vesicles in 40 h of incubation, at 25 degreesC. Alkaline phosphatase-ghost cell systems were relatively homogeneous with diameters of about 300 nm and were more stable when stored at -20 degreesC.Alkaline phosphatase was completely released from the resealed ghost cell-system using only phospholipase C. These experiments confirm that the interaction between alkaline phosphatase and the lipid bilayer of resealed ghost cell is exclusively via glycosylphosphatidylinositol (GPI) anchor of the enzyme.An important point shown is that an enzyme bound to resealed ghost cell does not lose the ability to hydrolyze ATP, pyrophosphate and p-nitrophenyl phosphate (PNPP), but the presence of a ghost membrane, as a support of the enzyme, affects its kinetic properties. Moreover, calcium ions stimulate and phosphate ions inhibit the PNPPase activity of alkaline phosphatase present in resealed ghost cells. (C) 2002 Elsevier B.V. B.V. All rights reserved.

Multiphase material modelling by multiscale particle-in-cell method

A particle-based method for multiscale modeling of multiphase materials such as Dual Phase (DP) and Transformation Induced Plasticity (TRIP) steels has been developed. The multiscale Particle-In-Cell (PIC) method benefits from the many advantages of the FEM and mesh-free methods, and to bridge the micro and macro scales through homogenization. The conventional mesh-based modeling methods fail to give reasonable and accurate predictions for materials with complex microstructures. Alternatively in the multiscale PIC method, the Lagrangian particles moving in an Eulerian grid represent the material deformation at both the micro and macro scales. The uniaxial tension test of two phase and three-phase materials was simulated and compared with FE based simulations. The predictions using multiscale PIC method showed that accuracy of field variables could be improved by up to 7%. This can lead to more accurate forming and springback predictions for materials with important multiphase microstructural effects.

Numerical and analytical vibro-acoustic modelling of porous materials: the Cell Method for poroelasticity and the case of inclusions

Porous materials are widely used in many fields of industrial applications, to achieve the requirements of noise reduction, that nowadays derive from strict regulations. The modeling of porous materials is still a problematic issue. Numerical simulations are often problematic in case of real complex geometries, especially in terms of computational times and convergence. At the same time, analytical models, even if partly limited by restrictive simplificative hypotheses, represent a powerful instrument to capture quickly the physics of the problem and general trends. In this context, a recently developed numerical method, called the Cell Method, is described, is presented in the case of the Biot's theory and applied for representative cases. The peculiarity of the Cell Method is that it allows for a direct algebraic and geometrical discretization of the field equations, without any reduction to a weak integral form. Then, the second part of the thesis presents the case of interaction between two poroelastic materials under the context of double porosity. The idea of using periodically repeated inclusions of a second porous material into a layer composed by an original material is described. In particular, the problem is addressed considering the efficiency of the analytical method. A analytical procedure for the simulation of heterogeneous layers based is described and validated considering both conditions of absorption and transmission; a comparison with the available numerical methods is performed. ---------------- I materiali porosi sono ampiamente utilizzati per diverse applicazioni industriali, al fine di raggiungere gli obiettivi di riduzione del rumore, che sono resi impegnativi da norme al giorno d'oggi sempre più stringenti. La modellazione dei materiali porori per applicazioni vibro-acustiche rapprensenta un aspetto di una certa complessità. Le simulazioni numeriche sono spesso problematiche quando siano coinvolte geometrie di pezzi reali, in particolare riguardo i tempi computazionali e la convergenza. Allo stesso tempo, i modelli analitici, anche se parzialmente limitati a causa di ipotesi semplificative che ne restringono l'ambito di utilizzo, rappresentano uno strumento molto utile per comprendere rapidamente la fisica del problema e individuare tendenze generali. In questo contesto, un metodo numerico recentemente sviluppato, il Metodo delle Celle, viene descritto, implementato nel caso della teoria di Biot per la poroelasticità e applicato a casi rappresentativi. La peculiarità del Metodo delle Celle consiste nella discretizzazione diretta algebrica e geometrica delle equazioni di campo, senza alcuna riduzione a forme integrali deboli. Successivamente, nella seconda parte della tesi viene presentato il caso delle interazioni tra due materiali poroelastici a contatto, nel contesto dei materiali a doppia porosità. Viene descritta l'idea di utilizzare inclusioni periodicamente ripetute di un secondo materiale poroso all'interno di un layer a sua volta poroso. In particolare, il problema è studiando il metodo analitico e la sua efficienza. Una procedura analitica per il calcolo di strati eterogenei di materiale viene descritta e validata considerando sia condizioni di assorbimento, sia di trasmissione; viene effettuata una comparazione con i metodi numerici a disposizione.

基于多矩的有限体积浸入边界法研究

浸入边界法（Immersed Boundary Method）是计算流体力学中求解具有复杂、移动边界流动问题的一类有效途径，该方法在笛卡尔坐标系上离散求解流体控制方程，并通过在控制方程中添加相应源相来代表浸入边界。尽管浸入边界法借助其简单、高效的显著特点在计算流体力学应用中显示出极强的生命力，特别是针对复杂的实际流动及动边界流动问题有着无可比拟的优势，但仍有许多问题需要进一步的研究。 本论文基于浸入边界方法及多矩VSIAM3（Volume/Surface Integrated Average Multi-Moment Method）格式提出了一种不可压缩流体求解数值格式。不可压N-S方程使用VSIAM3格式进行法进行离散，引入浸入边界法处理复杂、移动流动边界条件，使用虚拟网格方法计算动量方程修正项，同时还考虑了对连续方程的修正。VSIAM3格式是一种基于多矩的有限体积法，在方程的离散中总是使用两种或两种以上的矩，如：VIA（Volume Integrated Average）和SIA（Surface Integrated Average）。而不同的矩在求解过程中依据不同形式的控制方程使用不同的离散方法进行更新。VSIAM3格式更多的局地自由度及同时使用交错网格和同位网格的特点使浸入边界法的实施更加便利、高效。研究中，浸入边界法不仅应用于处理动力边界条件，同样可以处理热动力边界条件。 研究中对大量经典算例进行了数值实验，包括一维线性初始问题、方腔流问题、二维绕静止及振荡圆柱流动、三维绕球流动及热对流问题等。数值结果同实验值及其它计算结果保持一致，该算法可准确、高效处理具有复杂、移动边界及存在热对流的不可压流动问题，为实际应用打下了基础。

Introduction of high-content screening method for studying drug-induced mitochondrial dysfunction in hepatocytes

Drug induced liver injury is one of the frequent reasons for the drug removal from the market. During the recent years there has been a pressure to develop more cost efficient, faster and easier ways to investigate drug-induced toxicity in order to recognize hepatotoxic drugs in the earlier phases of drug development. High Content Screening (HCS) instrument is an automated microscope equipped with image analysis software. It makes the image analysis faster and decreases the risk for an error caused by a person by analyzing the images always in the same way. Because the amount of drug and time needed in the analysis are smaller and multiple parameters can be analyzed from the same cells, the method should be more sensitive, effective and cheaper than the conventional assays in cytotoxicity testing. Liver cells are rich in mitochondria and many drugs target their toxicity to hepatocyte mitochondria. Mitochondria produce the majority of the ATP in the cell through oxidative phosphorylation. They maintain biochemical homeostasis in the cell and participate in cell death. Mitochondria is divided into two compartments by inner and outer mitochondrial membranes. The oxidative phosphorylation happens in the inner mitochondrial membrane. A part of the respiratory chain, a protein called cytochrome c, activates caspase cascades when released. This leads to apoptosis. The aim of this study was to implement, optimize and compare mitochondrial toxicity HCS assays in live cells and fixed cells in two cellular models: human HepG2 hepatoma cell line and rat primary hepatocytes. Three different hepato- and mitochondriatoxic drugs (staurosporine, rotenone and tolcapone) were used. Cells were treated with the drugs, incubated with the fluorescent probes and then the images were analyzed using Cellomics ArrayScan VTI reader. Finally the results obtained after optimizing methods were compared to each other and to the results of the conventional cytotoxicity assays, ATP and LDH measurements. After optimization the live cell method and rat primary hepatocytes were selected to be used in the experiments. Staurosporine was the most toxic of the three drugs and caused most damage to the cells most quickly. Rotenone was not that toxic, but the results were more reproducible and thus it would serve as a good positive control in the screening. Tolcapone was the least toxic. So far the conventional analysis of cytotoxicity worked better than the HCS methods. More optimization needs to be done to get the HCS method more sensitive. This was not possible in this study due to time limit.

利用Ghost Fluid方法模拟激波与柱形界面相互作用

利用Ghost Fluid方法(后面简称Ghost方法)和γ-model方法，在同样的时空离散精度条件下，对激波与柱形界面相互作用的二维可压缩流场进行了直接模拟，并与实验结果相比较．从模拟结果看，在短时间内， Ghost方法和γ-model方法模拟的结果与实验结果基本相同，两种方法均正确地模拟出界面的位置、激波的强度和速度．但随着时间的发展，具有较大数值耗散的γ-model方法的计算结果与实验差别越来越大；而数值耗散较小的Ghost方法能较为正确地模拟界面的运动．

Comparison of terrain following and cut cell grids using a non-hydrostatic model

Terrain following coordinates are widely used in operational models but the cut cell method has been proposed as an alternative that can more accurately represent atmospheric dynamics over steep orography. Because the type of grid is usually chosen during model implementation, it becomes necessary to use different models to compare the accuracy of different grids. In contrast, here a C-grid finite volume model enables a like-for-like comparison of terrain following and cut cell grids. A series of standard two-dimensional tests using idealised terrain are performed: tracer advection in a prescribed horizontal velocity field, a test starting from resting initial conditions, and orographically induced gravity waves described by nonhydrostatic dynamics. In addition, three new tests are formulated: a more challenging resting atmosphere case, and two new advection tests having a velocity field that is everywhere tangential to the terrain following coordinate surfaces. These new tests present a challenge on cut cell grids. The results of the advection tests demonstrate that accuracy depends primarily upon alignment of the flow with the grid rather than grid orthogonality. A resting atmosphere is well-maintained on all grids. In the gravity waves test, results on all grids are in good agreement with existing results from the literature, although terrain following velocity fields lead to errors on cut cell grids. Due to semi-implicit timestepping and an upwind-biased, explicit advection scheme, there are no timestep restrictions associated with small cut cells. We do not find the significant advantages of cut cells or smoothed coordinates that other authors find.

Monte carlo simulation for molecular gas dynamics

The dynamics of low-density flows is governed by the Boltzmann equation of the kinetic theory of gases. This is a nonlinear integro-differential equation and, in general, numerical methods must be used to obtain its solution. The present paper, after a brief review of Direct Simulation Monte Carlo (DSMC) methods due to Bird, and Belotserkovskii and Yanitskii, studies the details of theDSMC method of Deshpande for mono as well as multicomponent gases. The present method is a statistical particle-in-cell method and is based upon the Kac-Prigogine master equation which reduces to the Boltzmann equation under the hypothesis of molecular chaos. The proposed Markoff model simulating the collisions uses a Poisson distribution for the number of collisions allowed in cells into which the physical space is divided. The model is then extended to a binary mixture of gases and it is shown that it is necessary to perform the collisions in a certain sequence to obtain unbiased simulation.

Mutual diffusion coefficients of some binary liquid systems: benzene-n-alkyl alcohol

The mutual diffusion coefficients for binary liquid systems of benzene-n-alkyl alcohol at various compositions have been determined by the diaphragm cell method at 28-degrees-C. The alcohols used were the members of n-paraffinic alcohols ranging from C1 to C8. The maximum possible experimental error is 14%. The data were fitted with a generalized correlation, giving the deviation from the experimental data to within 2.75%, on average.

Simulação numérica do escoamento de água em áreas alagáveis da floresta amazônica com a utilização da estrutura de dados Autonomous Leaves Graph

A Amazônia exibe uma variedade de cenários que se complementam. Parte desse ecossistema sofre anualmente severas alterações em seu ciclo hidrológico, fazendo com que vastos trechos de floresta sejam inundados. Esse fenômeno, entretanto, é extremamente importante para a manutenção de ciclos naturais. Neste contexto, compreender a dinâmica das áreas alagáveis amazônicas é importante para antecipar o efeito de ações não sustentáveis. Sob esta motivação, este trabalho estuda um modelo de escoamento em áreas alagáveis amazônicas, baseado nas equações de Navier-Stokes, além de ferramentas que possam ser aplicadas ao modelo, favorecendo uma nova abordagem do problema. Para a discretização das equações é utilizado o Método dos Volumes Finitos, sendo o Método do Gradiente Conjugado a técnica escolhida para resolver os sistemas lineares associados. Como técnica de resolução numérica das equações, empregou-se o Método Marker and Cell, procedimento explícito para solução das equações de Navier-Stokes. Por fim, as técnicas são aplicadas a simulações preliminares utilizando a estrutura de dados Autonomous Leaves Graph, que tem recursos adaptativos para manipulação da malha que representa o domínio do problema

Computational fluid-structure interaction of DGB parachutes in compressible fluid flow

This paper describes large-scale simulations of compressible flows over a supersonic disk-gap-band parachute system. An adaptive mesh refinement method is used to resolve the coupled fluid-structure model. The fluid model employs large-eddy simulation to describe the turbulent wakes appearing upstream and downstream of the parachute canopy and the structural model employed a thin-shell finite element solver that allows large canopy deformations by using subdivision finite elements. The fluid-structure interaction is described by a variant of the Ghost-Fluid method. The simulation was carried out at Mach number 1.96 where strong nonlinear coupling between the system of bow shocks, turbulent wake and canopy is observed. It was found that the canopy oscillations were characterized by a breathing type motion due to the strong interaction of the turbulent wake and bow shock upstream of the flexible canopy. Copyright © 2010 by ASME.

Stability of numerical schemes on staggered grids

This paper considers the stability of explicit, implicit and Crank-Nicolson schemes for the one-dimensional heat equation on a staggered grid. Furthemore, we consider the cases when both explicit and implicit approximations of the boundary conditions arc employed. Why we choose to do this is clearly motivated and arises front solving fluid flow equations with free surfaces when the Reynolds number can be very small. in at least parts of the spatial domain. A comprehensive stability analysis is supplied: a novel result is the precise stability restriction on the Crank-Nicolson method when the boundary conditions are approximated explicitly, that is, at t =n delta t rather than t = (n + 1)delta t. The two-dimensional Navier-Stokes equations were then solved by a marker and cell approach for two simple problems that had analytic solutions. It was found that the stability results provided in this paper were qualitatively very similar. thereby providing insight as to why a Crank-Nicolson approximation of the momentum equations is only conditionally, stable. Copyright (C) 2008 John Wiley & Sons, Ltd.

Numerical simulation of drop impact and jet buckling problems using the eXtended Pom-Pom model

This work presents numerical simulations of two fluid flow problems involving moving free surfaces: the impacting drop and fluid jet buckling. The viscoelastic model used in these simulations is the eXtended Pom-Pom (XPP) model. To validate the code, numerical predictions of the drop impact problem for Newtonian and Oldroyd-B fluids are presented and compared with other methods. In particular, a benchmark on numerical simulations for a XPP drop impacting on a rigid plate is performed for a wide range of the relevant parameters. Finally, to provide an additional application of free surface flows of XPP fluids, the viscous jet buckling problem is simulated and discussed. (C) 2011 Elsevier B.V. All rights reserved.