Many-body effect in electrorheological responses

The many-body effect in the kinetic responses of ER fluids is studied by a molecular-dynamic simulation method. The mutual polarization effects of the particles are considered by self-consistently calculating the dipole strength on each particle according to the external field and the dipole field due to all the other particles in the fluids. The many-body effect is found to increase with the enhancement of the particle concentration and the permittivity ratio between the solvent and the particles. The calculated response times are shorter than that predicted with the 'point-dipole' model and agree very well with experimental results. The many-body effect enhances the shear stresses of the fluids by several times. But they are not proportional to the many-body correction factor lambda as expected. This is due to the fact that larger interaction forces between the particles lead to coarsening of the fibers formed in the suspensions. The results show that the many-body and multipolar interaction between the particles must be treated comprehensively in the simulations in order to get more reliable results.

Cool roofs for passive cooling: performance in different climates and for different insulation levels in Italy

Roofs are severely hit by solar radiation in summer; hence the use of cool materials on the finishing layer provides a significant reduction in the heat flow entering the building, with sensible attenuation in the building cooling load. In this paper, a case study is presented, based on the dynamic simulation of an existing office building in Catania (southern Italy). Here, a part of the roof has been recently treated with a commercial cool paint, with the aim of improving thermal comfort in summer. Hence, the simulations represent a preliminary study that will allow assessing the expected effectiveness of the intervention. More in detail, the results of the simulations will be discussed in terms of both thermal comfort and energy savings, through the evaluation of parameters such as the roof surface temperature, the operative temperature and the cooling load for both conditions, i.e. with and without the cool paint. The paper also discusses the potential increase in the energy needs for winter heating, and looks at the overall annual balance in terms of primary energy; this is made by considering different climatic conditions and envelope characteristics. These aspects are usually not well highlighted in the current scientific literature.

Structural and thermodynamic analysis of thrombin:suramin interaction in solution and crystal phases

Suramin is a hexasulfonated naphthylurea which has been recently characterized as a non-competitive inhibitor of human alpha-thrombin activity over fibrinogen, although its binding site and mode of interaction with the enzyme remain elusive. Here, we determined two X-ray structure of the thrombin: suramin complex, refined at 2.4 angstrom resolution. While a single thrombin: suramin complex was found in the asymmetric unit cell of the crystal, some of the crystallographic contacts with symmetrically related molecules are mediated by both the enzyme and the ligand. Molecular dynamics simulations with the 1:1 complex demonstrate a large rearrangement of suramin in the complex, but with the protein scaffold and the more extensive protein-ligand regions keep unchanged. Small-angle X-ray scattering measurements at high micromolar concentration demonstrate a suramin-induced dimerization of the enzyme. These data indicating a dissimilar binding mode in the monomeric and oligomeric states, with a monomeric, 1:1 complex to be more likely to exist at the thrombin physiological, nanomolar concentration range. Collectively, close understanding on the structural basis for interaction is given which might establish a basis for design of suramin analogues targeting thrombin. Crown Copyright (C) 2009 Published by Elsevier B.V. All rights reserved.

Constructing a Structural Macroeconomic Model Using Multiple Techniques

This project constructs a structural model of the United States Economy. This task is tackled in two separate ways: first econometric methods and then using a neural network, both with a structure that mimics the structure of the U.S. economy. The structural model tracks the performance of U.S. GDP rather well in a dynamic simulation, with an average error of just over 1 percent. The neural network performed well, but suffered from some theoretical, as well as some implementation issues.

Simulação computacional de medidas estereológicas em estruturas de metal duro (WC-Co)

This work focuses on the creation and applications of a dynamic simulation software in order to study the hard metal structure (WC-Co). The technological ground used to increase the GPU hardware capacity was Geforce 9600 GT along with the PhysX chip created to make games more realistic. The software simulates the three-dimensional carbide structure to the shape of a cubic box where tungsten carbide (WC) are modeled as triangular prisms and truncated triangular prisms. The program was proven effective regarding checking testes, ranging from calculations of parameter measures such as the capacity to increase the number of particles simulated dynamically. It was possible to make an investigation of both the mean parameters and distributions stereological parameters used to characterize the carbide structure through cutting plans. Grounded on the cutting plans concerning the analyzed structures, we have investigated the linear intercepts, the intercepts to the area, and the perimeter section of the intercepted grains as well as the binder phase to the structure by calculating the mean value and distribution of the free path. As literature shows almost consensually that the distribution of the linear intercepts is lognormal, this suggests that the grain distribution is also lognormal. Thus, a routine was developed regarding the program which made possible a more detailed research on this issue. We have observed that it is possible, under certain values for the parameters which define the shape and size of the Prismatic grain to find out the distribution to the linear intercepts that approach the lognormal shape. Regarding a number of developed simulations, we have observed that the distribution curves of the linear and area intercepts as well as the perimeter section are consistent with studies on static computer simulation to these parameters.

Desenvolvimento de Modelos no Programa DigSILENT PowerFactoryTM para o Estudo de Estabilidade Transitória em Sistemas Elétricos de Potência com Aerogerador

This work presents a description of models development at DigSILENT PowerFactoryT M program for the transient stability study in power systems with wind turbine. The main goal is to make available means to use a dynamic simulation program in power systems, widely published, and utilize it as a tool that helps in programs results evaluations used for this intent. The process of simulations and analyses results starts after the models setting description phase. The results obtained by the DigSILENT PowerFactoryT M and ATP, program chosen to the validation also international recognized, are compared during this phase. The main tools and guide lines of PowerFactoryT M program use are presented here, directing these elements to the solution of the approached problem. For the simulation it is used a real system which it will be connected a wind farm. Two different technologies of wind turbines were implemented: doublyfed induction generator with frequency converter, connecting the rotor to the stator and to the grid, and synchronous wind generator with frequency converter, interconnecting the generator to the grid. Besides presenting the basic conceptions of dynamic simulation, it is described the implemented control strategies and models of turbine and converters. The stability of the wind turbine interconnected to grid is analyzed in many operational conditions, resultant of diverse kinds of disturbances

Desenvolvimento de modelos no programa digSILENT powerfactoryTM para o estudo de estabilidade transitória em sistemas elétricos de potência com aerogerador

This work presents a description of models development at DigSILENT PowerFactoryTM program for the transient stability study in power systems with wind turbine. The main goal is to make available means to use a dynamic simulation program in power systems, widely published, and utilize it as a tool that helps in programs results evaluations used for this intent. The process of simulations and analyses results starts after the models setting description phase. The results obtained by the DigSILENT PowerFactoryTM and ATP, program chosen to the validation also international recognized, are compared during this phase. The main tools and guide lines of PowerFactoryTM program use are presented here, directing these elements to the solution of the approached problem. For the simulation it is used a real system which it will be connected a wind farm. Two different technologies of wind turbines were implemented: doubly-fed induction generator with frequency converter, connecting the rotor to the stator and to the grid, and synchronous wind generator with frequency converter, interconnecting the generator to the grid. Besides presenting the basic conceptions of dynamic simulation, it is described the implemented control strategies and models of turbine and converters. The stability of the wind turbine interconnected to grid is analyzed in many operational conditions, resultant of diverse kinds of disturbances

Simulação computacional da interação fluido-estrutura em bombas de cavidades progressivas

The pumping through progressing cavities system has been more and more employed in the petroleum industry. This occurs because of its capacity of elevation of highly viscous oils or fluids with great concentration of sand or other solid particles. A Progressing Cavity Pump (PCP) consists, basically, of a rotor - a metallic device similar to an eccentric screw, and a stator - a steel tube internally covered by a double helix, which may be rigid or deformable/elastomeric. In general, it is submitted to a combination of well pressure with the pressure generated by the pumping process itself. In elastomeric PCPs, this combined effort compresses the stator and generates, or enlarges, the clearance existing between the rotor and the stator, thus reducing the closing effect between their cavities. Such opening of the sealing region produces what is known as fluid slip or slippage, reducing the efficiency of the PCP pumping system. Therefore, this research aims to develop a transient three-dimensional computational model that, based on single-lobe PCP kinematics, is able to simulate the fluid-structure interaction that occurs in the interior of metallic and elastomeric PCPs. The main goal is to evaluate the dynamic characteristics of PCP s efficiency based on detailed and instantaneous information of velocity, pressure and deformation fields in their interior. To reach these goals (development and use of the model), it was also necessary the development of a methodology for generation of dynamic, mobile and deformable, computational meshes representing fluid and structural regions of a PCP. This additional intermediary step has been characterized as the biggest challenge for the elaboration and running of the computational model due to the complex kinematic and critical geometry of this type of pump (different helix angles between rotor and stator as well as large length scale aspect ratios). The processes of dynamic generation of meshes and of simultaneous evaluation of the deformations suffered by the elastomer are fulfilled through subroutines written in Fortan 90 language that dynamically interact with the CFX/ANSYS fluid dynamic software. Since a structural elastic linear model is employed to evaluate elastomer deformations, it is not necessary to use any CAE package for structural analysis. However, an initial proposal for dynamic simulation using hyperelastic models through ANSYS software is also presented in this research. Validation of the results produced with the present methodology (mesh generation, flow simulation in metallic PCPs and simulation of fluid-structure interaction in elastomeric PCPs) is obtained through comparison with experimental results reported by the literature. It is expected that the development and application of such a computational model may provide better details of the dynamics of the flow within metallic and elastomeric PCPs, so that better control systems may be implemented in the artificial elevation area by PCP

An approach of optimal sensitivity applied in the tertiary loop of the automatic generation control

This paper proposes an approach of optimal sensitivity applied in the tertiary loop of the automatic generation control. The approach is based on the theorem of non-linear perturbation. From an optimal operation point obtained by an optimal power flow a new optimal operation point is directly determined after a perturbation, i.e., without the necessity of an iterative process. This new optimal operation point satisfies the constraints of the problem for small perturbation in the loads. The participation factors and the voltage set point of the automatic voltage regulators (AVR) of the generators are determined by the technique of optimal sensitivity, considering the effects of the active power losses minimization and the network constraints. The participation factors and voltage set point of the generators are supplied directly to a computational program of dynamic simulation of the automatic generation control, named by power sensitivity mode. Test results are presented to show the good performance of this approach. (C) 2008 Elsevier B.V. All rights reserved.

A novel set of ss-N-biaryl ether sulfonamide hydroxamates as potential MMPs inhibitors: Molecular dynamics simulations and molecular properties evaluation

Matrix metalloproteinases (MMPs) constitute a family of zinc-dependent proteases involved in the extracellular matrix degradation. MMP-2 and MMP9 are overexpressed in several human cancer types, including melanoma, thus the development of new compounds to inhibit MMPs' activity is desirable. Molecular dynamic simulation and molecular properties calculations were performed on a set of novel beta-N-biaryl ether sulfonamide-based hydroxamates, reported as MMP-2 and MMP-9 inhibitors, for providing data to develop an exploratory analysis. Thermodynamic, electronic, and steric descriptors have significantly discriminated highly active from moderately and less active inhibitors of MMP-2 whereas apparent partition coefficient at pH 1.5 was also significant for the MMP-9 data set. Compound 47 was considered an outlier in all analysis, indicating the presence of a bulky substituent group in R3 is crucial to this set of inhibitors for the establishment of molecular interactions with the S1 subsite of both enzymes, but there is a limit. (C) 2012 Wiley Periodicals, Inc.

Metodi numerici per lo studio di centri di lavoro ad alta velocità

A Machining Centre is nowadays a complex mechanical, electronic, electrical system that needs integrated design capabilities which very often require a high time-consuming effort. Numerical techniques for designing and dimensioning the machine structure and components usually requires different knowledge according to the system that have to be designed. This Ph. D Thesis is related about the efforts of the Authors to develop a system that allows to perform the complete project of a new machine optimized in its dynamic behaviour. An integration of the different systems developed, each of which respond to specific necessities of designer, is here presented. In particular a dynamic analysis system, based on a lumped mass approach, that rapidly allows to setup the drives of the machine and an Integrated Dynamic Simulation System, based on a FEM approach, that permit a dynamic optimization, are shown. A multilevel Data Base, and an operator interface module provide to complete the designing platform. The proposed approach represents a significant step toward the virtual machining for the prediction of the quality of the worked surface.

Modellazione multibody di veicoli ferroviari: sviluppo ed implementazione di modelli innovativi per l'analisi del contatto ruota - rotaia

The wheel - rail contact analysis plays a fundamental role in the multibody modeling of railway vehicles. A good contact model must provide an accurate description of the global contact phenomena (contact forces and torques, number and position of the contact points) and of the local contact phenomena (position and shape of the contact patch, stresses and displacements). The model has also to assure high numerical efficiency (in order to be implemented directly online within multibody models) and a good compatibility with commercial multibody software (Simpack Rail, Adams Rail). The wheel - rail contact problem has been discussed by several authors and many models can be found in the literature. The contact models can be subdivided into two different categories: the global models and the local (or differential) models. Currently, as regards the global models, the main approaches to the problem are the so - called rigid contact formulation and the semi – elastic contact description. The rigid approach considers the wheel and the rail as rigid bodies. The contact is imposed by means of constraint equations and the contact points are detected during the dynamic simulation by solving the nonlinear algebraic differential equations associated to the constrained multibody system. Indentation between the bodies is not permitted and the normal contact forces are calculated through the Lagrange multipliers. Finally the Hertz’s and the Kalker’s theories allow to evaluate the shape of the contact patch and the tangential forces respectively. Also the semi - elastic approach considers the wheel and the rail as rigid bodies. However in this case no kinematic constraints are imposed and the indentation between the bodies is permitted. The contact points are detected by means of approximated procedures (based on look - up tables and simplifying hypotheses on the problem geometry). The normal contact forces are calculated as a function of the indentation while, as in the rigid approach, the Hertz’s and the Kalker’s theories allow to evaluate the shape of the contact patch and the tangential forces. Both the described multibody approaches are computationally very efficient but their generality and accuracy turn out to be often insufficient because the physical hypotheses behind these theories are too restrictive and, in many circumstances, unverified. In order to obtain a complete description of the contact phenomena, local (or differential) contact models are needed. In other words wheel and rail have to be considered elastic bodies governed by the Navier’s equations and the contact has to be described by suitable analytical contact conditions. The contact between elastic bodies has been widely studied in literature both in the general case and in the rolling case. Many procedures based on variational inequalities, FEM techniques and convex optimization have been developed. This kind of approach assures high generality and accuracy but still needs very large computational costs and memory consumption. Due to the high computational load and memory consumption, referring to the current state of the art, the integration between multibody and differential modeling is almost absent in literature especially in the railway field. However this integration is very important because only the differential modeling allows an accurate analysis of the contact problem (in terms of contact forces and torques, position and shape of the contact patch, stresses and displacements) while the multibody modeling is the standard in the study of the railway dynamics. In this thesis some innovative wheel – rail contact models developed during the Ph. D. activity will be described. Concerning the global models, two new models belonging to the semi – elastic approach will be presented; the models satisfy the following specifics: 1) the models have to be 3D and to consider all the six relative degrees of freedom between wheel and rail 2) the models have to consider generic railway tracks and generic wheel and rail profiles 3) the models have to assure a general and accurate handling of the multiple contact without simplifying hypotheses on the problem geometry; in particular the models have to evaluate the number and the position of the contact points and, for each point, the contact forces and torques 4) the models have to be implementable directly online within the multibody models without look - up tables 5) the models have to assure computation times comparable with those of commercial multibody software (Simpack Rail, Adams Rail) and compatible with RT and HIL applications 6) the models have to be compatible with commercial multibody software (Simpack Rail, Adams Rail). The most innovative aspect of the new global contact models regards the detection of the contact points. In particular both the models aim to reduce the algebraic problem dimension by means of suitable analytical techniques. This kind of reduction allows to obtain an high numerical efficiency that makes possible the online implementation of the new procedure and the achievement of performance comparable with those of commercial multibody software. At the same time the analytical approach assures high accuracy and generality. Concerning the local (or differential) contact models, one new model satisfying the following specifics will be presented: 1) the model has to be 3D and to consider all the six relative degrees of freedom between wheel and rail 2) the model has to consider generic railway tracks and generic wheel and rail profiles 3) the model has to assure a general and accurate handling of the multiple contact without simplifying hypotheses on the problem geometry; in particular the model has to able to calculate both the global contact variables (contact forces and torques) and the local contact variables (position and shape of the contact patch, stresses and displacements) 4) the model has to be implementable directly online within the multibody models 5) the model has to assure high numerical efficiency and a reduced memory consumption in order to achieve a good integration between multibody and differential modeling (the base for the local contact models) 6) the model has to be compatible with commercial multibody software (Simpack Rail, Adams Rail). In this case the most innovative aspects of the new local contact model regard the contact modeling (by means of suitable analytical conditions) and the implementation of the numerical algorithms needed to solve the discrete problem arising from the discretization of the original continuum problem. Moreover, during the development of the local model, the achievement of a good compromise between accuracy and efficiency turned out to be very important to obtain a good integration between multibody and differential modeling. At this point the contact models has been inserted within a 3D multibody model of a railway vehicle to obtain a complete model of the wagon. The railway vehicle chosen as benchmark is the Manchester Wagon the physical and geometrical characteristics of which are easily available in the literature. The model of the whole railway vehicle (multibody model and contact model) has been implemented in the Matlab/Simulink environment. The multibody model has been implemented in SimMechanics, a Matlab toolbox specifically designed for multibody dynamics, while, as regards the contact models, the CS – functions have been used; this particular Matlab architecture allows to efficiently connect the Matlab/Simulink and the C/C++ environment. The 3D multibody model of the same vehicle (this time equipped with a standard contact model based on the semi - elastic approach) has been then implemented also in Simpack Rail, a commercial multibody software for railway vehicles widely tested and validated. Finally numerical simulations of the vehicle dynamics have been carried out on many different railway tracks with the aim of evaluating the performances of the whole model. The comparison between the results obtained by the Matlab/ Simulink model and those obtained by the Simpack Rail model has allowed an accurate and reliable validation of the new contact models. In conclusion to this brief introduction to my Ph. D. thesis, we would like to thank Trenitalia and the Regione Toscana for the support provided during all the Ph. D. activity. Moreover we would also like to thank the INTEC GmbH, the society the develops the software Simpack Rail, with which we are currently working together to develop innovative toolboxes specifically designed for the wheel rail contact analysis.

Simulazione dinamica di un veicolo dotato di powertrain ibrido endotermico-elettrico

Nel panorama motoristico ed automobilistico moderno lo sviluppo di motori a combustione interna e veicoli è fortemente influenzato da diverse esigenze che spesso sono in contrasto le une con le altre. Infatti gli obiettivi di economicità e riduzione dei costi riguardanti la produzione e la commercializzazione dei prodotti sono in contrasto con gli sforzi che devono essere operati dalle case produttrici per soddisfare le sempre più stringenti normative riguardanti le emissioni inquinanti ed i consumi di carburante dei veicoli. Fra le numerose soluzioni presenti i veicoli ibridi rappresentano una alternativa che allo stato attuale è già presente sul mercato in varie forme, a seconda della tipologie di energie accoppiate. In letteratura è possibile trovare numerosi studi che trattano l’ottimizzazione dei componenti o delle strategie di controllo di queste tipologie di veicoli: in moltissimi casi l’obiettivo è quello di minimizzare consumi ed emissioni inquinanti. Normalmente non viene posta particolare attenzione agli effetti che l’aggiunta delle macchine elettriche e dei componenti necessari per il funzionamento delle stesse hanno sulla dinamica del veicolo. Il presente lavoro di tesi è incentrato su questi aspetti: si è considerata la tipologia di veicoli ibridi termici-elettrici di tipo parallelo andando ad analizzare come cambiasse il comportamento dinamico del veicolo in funzione del tipo di installazione considerato per la parte elettrica del powertrain. In primo luogo è stato quindi necessario costruire ed implementare un modello dinamico di veicolo che permettesse di applicare coppie alle quattro ruote in maniera indipendente per considerare diverse tipologie di powertrain. In seguito si sono analizzate le differenze di comportamento dinamico fra il veicolo considerato e l’equivalente versione ibrida e i possibili utilizzi delle macchine elettriche per correggere eventuali deterioramenti o cambiamenti indesiderati nelle prestazioni del veicolo.

OPTIMIZING A SHORT-END ELECTROPHORETICALLY MEDIATED MICRO-ANALYSIS (EMMA) ASSAY FOR CREATININE

In this work electrophoretically mediated micro-analysis (EMMA) is used in conjunction with short end injection to improve the in-capillary Jaffé assay for creatinine. Key advances over prior work include (i) using simulation to ensure intimate overlap of reagent plugs, (ii) using OH- to drive the reaction, (iii) using short-end injection to minimize analysis time and in-line product degradation. The potential-driven overlapping time with the EMMA approach, as well as the borate buffer background electrolyte (BGE) concentration and pH are optimized with the short end approach. The best conditions for short-end analyses would not have been predicted by the prior long end work, owing to a complex interplay of separation time and product degradation rates. Raw peak areas and flow-adjusted peak areas for the Jaffé reaction product (at 505 nm) are used to assess the sensitivity of the short-end EMMA approach. Optimal overlap conditions depend heavily on local conductivity differences within the reagent zone(s), as these differences cause dramatic voltage field differences, which effect reagent overlap dynamics. Simul 5.0, a dynamic simulation program for capillary electrophoresis (CE) systems, is used to understand the ionic boundaries and profiles that give rise to the experimentally obtained data for EMMA analysis. Overall, fast migration of hydroxide ions from the picrate zone makes difficult reagent overlap. In addition, the challenges associated with the simultaneous overlapping of three reagent zones are considered, and experimental results validate the predictions made by the simulation. With one set of “optimized” conditions including OH- (253 mM) as the third reagent zone the response was linear with creatinine concentration (R2 = 0.998) and reproducible over the clinically relevant range (0.08 to 0.1 mM) of standard creatinine concentrations. An LOD (S/N = 3) of 0.02 mM and LOQ (S/N=10) of 0.08 mM were determined. A significant improvement (43%) in assay sensitivity was obtained compared to prior work that considered only two reagents in the overlap.

Proteins in silico-modeling and sampling

Proteins are linear chain molecules made out of amino acids. Only when they fold to their native states, they become functional. This dissertation aims to model the solvent (environment) effect and to develop & implement enhanced sampling methods that enable a reliable study of the protein folding problem in silico. We have developed an enhanced solvation model based on the solution to the Poisson-Boltzmann equation in order to describe the solvent effect. Following the quantum mechanical Polarizable Continuum Model (PCM), we decomposed net solvation free energy into three physical terms– Polarization, Dispersion and Cavitation. All the terms were implemented, analyzed and parametrized individually to obtain a high level of accuracy. In order to describe the thermodynamics of proteins, their conformational space needs to be sampled thoroughly. Simulations of proteins are hampered by slow relaxation due to their rugged free-energy landscape, with the barriers between minima being higher than the thermal energy at physiological temperatures. In order to overcome this problem a number of approaches have been proposed of which replica exchange method (REM) is the most popular. In this dissertation we describe a new variant of canonical replica exchange method in the context of molecular dynamic simulation. The advantage of this new method is the easily tunable high acceptance rate for the replica exchange. We call our method Microcanonical Replica Exchange Molecular Dynamic (MREMD). We have described the theoretical frame work, comment on its actual implementation, and its application to Trp-cage mini-protein in implicit solvent. We have been able to correctly predict the folding thermodynamics of this protein using our approach.