Multiscale Modeling and Simulation of Stepped Crystal Surfaces

A primary goal of this dissertation is to understand the links between mathematical models that describe crystal surfaces at three fundamental length scales: The scale of individual atoms, the scale of collections of atoms forming crystal defects, and macroscopic scale. Characterizing connections between different classes of models is a critical task for gaining insight into the physics they describe, a long-standing objective in applied analysis, and also highly relevant in engineering applications. The key concept I use in each problem addressed in this thesis is coarse graining, which is a strategy for connecting fine representations or models with coarser representations. Often this idea is invoked to reduce a large discrete system to an appropriate continuum description, e.g. individual particles are represented by a continuous density. While there is no general theory of coarse graining, one closely related mathematical approach is asymptotic analysis, i.e. the description of limiting behavior as some parameter becomes very large or very small. In the case of crystalline solids, it is natural to consider cases where the number of particles is large or where the lattice spacing is small. Limits such as these often make explicit the nature of links between models capturing different scales, and, once established, provide a means of improving our understanding, or the models themselves. Finding appropriate variables whose limits illustrate the important connections between models is no easy task, however. This is one area where computer simulation is extremely helpful, as it allows us to see the results of complex dynamics and gather clues regarding the roles of different physical quantities. On the other hand, connections between models enable the development of novel multiscale computational schemes, so understanding can assist computation and vice versa. Some of these ideas are demonstrated in this thesis. The important outcomes of this thesis include: (1) a systematic derivation of the step-flow model of Burton, Cabrera, and Frank, with corrections, from an atomistic solid-on-solid-type models in 1+1 dimensions; (2) the inclusion of an atomistically motivated transport mechanism in an island dynamics model allowing for a more detailed account of mound evolution; and (3) the development of a hybrid discrete-continuum scheme for simulating the relaxation of a faceted crystal mound. Central to all of these modeling and simulation efforts is the presence of steps composed of individual layers of atoms on vicinal crystal surfaces. Consequently, a recurring theme in this research is the observation that mesoscale defects play a crucial role in crystal morphological evolution.

Offshore Wind Energy System with DC Transmission Discrete Mass: Modeling and Simulation

This paper presents an integrated model for an offshore wind energy system taking into consideration a contribution for the marine wave and wind speed with perturbations influences on the power quality of current injected into the electric grid. The paper deals with the simulation of one floating offshore wind turbine equipped with a PMSG and a two-level converter connected to an onshore electric grid. The use of discrete mass modeling is accessed in order to reveal by computing the THD on how the perturbations of the captured energy are attenuated at the electric grid injection point. Two torque actions are considered for the three-mass modeling, the aerodynamic on the flexible part and on the rigid part of the blades. Also, a torque due to the influence of marine waves in deep water is considered. PI fractional-order control supports the control strategy. A comparison between the drive train models is presented.

Offshore Wind Energy Conversion System Connected to the Electric Grid: Modeling and Simulation

This paper is on modeling and simulation for an offshore wind system equipped with a semi-submersible floating platform, a wind turbine, a permanent magnet synchronous generator, a multiple point clamped four level or five level full-power converter, a submarine cable and a second order filter. The drive train is modeled by three mass model considering the resistant stiffness torque, structure and tower in deep water due to the moving surface elevation. The system control uses PWM by space vector modulation associated with sliding mode and proportional integral controllers. The electric energy is injected into the electric grid either by an alternated current link or by a direct current link. The model is intend to be a useful tool for unveil the behavior and performance of the offshore wind system, especially for the multiple point clamped full-power converter, under normal operation or under malfunctions.

Estudo do processo de drenagem gravitacional do óleo assistido com injeção de vapor e solvente

Como os recursos de hidrocarbonetos convencionais estão se esgotando, a crescente demanda mundial por energia impulsiona a indústria do petróleo para desenvolver mais reservatórios não convencionais. Os recursos mundiais de betume e óleo pesado são estimados em 5,6 trilhões de barris, dos quais 80% estão localizados na Venezuela, Canadá e EUA. Um dos métodos para explorar estes hidrocarbonetos é o processo de drenagem gravitacional assistido com injeção de vapor e solvente (ES-SAGD Expanding Solvent Steam Assisted Gravity Drainage). Neste processo são utilizados dois poços horizontais paralelos e situados verticalmente um acima do outro, um produtor na base do reservatório e um injetor de vapor e solvente no topo do reservatório. Este processo é composto por um método térmico (injeção de vapor) e um método miscível (injeção de solvente) com a finalidade de causar a redução das tensões interfaciais e da viscosidade do óleo ou betume. O objetivo deste estudo é analisar a sensibilidade de alguns parâmetros operacionais, tais como: tipo de solvente injetado, qualidade do vapor, distância vertical entre os poços, porcentagem de solvente injetado e vazão de injeção de vapor sobre o fator de recuperação para 5, 10 e 15 anos. Os estudos foram realizados através de simulações concretizadas no módulo STARS (Steam Thermal, and Advanced Processes Reservoir Simulator) do programa da CMG (Computer Modelling Group), versão 2010.10, onde as interações entre os parâmetros operacionais, estudados em um modelo homogêneo com características de reservatórios semelhantes aos encontrados no Nordeste Brasileiro, foram observadas. Os resultados obtidos neste estudo mostraram que os melhores fatores de recuperação ocorreram para níveis máximos do percentual de solvente injetado e da distância vertical entre os poços. Observou-se também que o processo será rentável dependendo do tipo e do valor do solvente injetado

Estudo do processo de drenagem gravitacional do óleo assistido com injeção de vapor e solvente

Como os recursos de hidrocarbonetos convencionais estão se esgotando, a crescente demanda mundial por energia impulsiona a indústria do petróleo para desenvolver mais reservatórios não convencionais. Os recursos mundiais de betume e óleo pesado são estimados em 5,6 trilhões de barris, dos quais 80% estão localizados na Venezuela, Canadá e EUA. Um dos métodos para explorar estes hidrocarbonetos é o processo de drenagem gravitacional assistido com injeção de vapor e solvente (ES-SAGD Expanding Solvent Steam Assisted Gravity Drainage). Neste processo são utilizados dois poços horizontais paralelos e situados verticalmente um acima do outro, um produtor na base do reservatório e um injetor de vapor e solvente no topo do reservatório. Este processo é composto por um método térmico (injeção de vapor) e um método miscível (injeção de solvente) com a finalidade de causar a redução das tensões interfaciais e da viscosidade do óleo ou betume. O objetivo deste estudo é analisar a sensibilidade de alguns parâmetros operacionais, tais como: tipo de solvente injetado, qualidade do vapor, distância vertical entre os poços, porcentagem de solvente injetado e vazão de injeção de vapor sobre o fator de recuperação para 5, 10 e 15 anos. Os estudos foram realizados através de simulações concretizadas no módulo STARS (Steam Thermal, and Advanced Processes Reservoir Simulator) do programa da CMG (Computer Modelling Group), versão 2010.10, onde as interações entre os parâmetros operacionais, estudados em um modelo homogêneo com características de reservatórios semelhantes aos encontrados no Nordeste Brasileiro, foram observadas. Os resultados obtidos neste estudo mostraram que os melhores fatores de recuperação ocorreram para níveis máximos do percentual de solvente injetado e da distância vertical entre os poços. Observou-se também que o processo será rentável dependendo do tipo e do valor do solvente injetado

Liquid-liquid slug flow: Hydrodynamics and pressure drop

In this paper, the hydrodynamics and the pressure drop of liquid-liquid slug flow in round microcapillaries are presented. Two liquid-liquid flow systems are considered, viz. water-toluene and ethylene glycol/water-toluene. The slug lengths of the alternating continuous and dispersed phases were measured as a function of the slug velocity (0.03-0.5 m/s), the organic-to-aqueous flow ratio (0.1-4.0), and the microcapillary internal diameter (248 and 498 mu m). The pressure drop is modeled as the sum of two contributions: the frictional and the interface pressure drop. Two models are presented, viz, the stagnant film model and the moving film model. Both models account for the presence of a thin liquid film between the dispersed phase slug and the capillary wall. It is found that the film velocity is of negligible influence on the pressure drop. Therefore, the stagnant film model is adequate to accurately predict the liquid-liquid slug flow pressure drop. The influence of inertia and the consequent change of the slug cap curvature are accounted for by modifying Bretherton's curvature parameter in the interface pressure drop equation. The stagnant film model is in good agreement with experimental data with a mean relative error of less than 7%.

Experimental assessment of pressure drop in air stripping columns

One important step in the design of air stripping operations for the removal of VOC is the choice of operating conditions, which are based in the phase ratio. This parameter sets on directly the stripping factor and the efficiency of the operation. Its value has an upper limit determined by the flooding regime, which is previewed using empirical correlations, namely the one developed by Eckert. This type of approach is not suitable for the development of algorithms. Using a pilot scale column and a convenient solution, the pressure drop was determined in different operating conditions and the experimental values were compared with the estimations. This particular research will be incorporated in a global model for simulating the dynamics of air stripping using a multi variable distributed parameter system.

Evaluation of the accuracy of an impact load cell using lumped parameter modeling and analysis

The objective of the current study is to evaluate the fidelity of load cell reading during impact testing in a drop-weight impactor using lumped parameter modeling. For the most common configuration of a moving impactor-load cell system in which dynamic load is transferred from the impactor head to the load cell, a quantitative assessment is made of the possible discrepancy that can result in load cell response. A 3-DOF (degrees-of-freedom) LPM (lumped parameter model) is considered to represent a given impact testing set-up. In this model, a test specimen in the form of a steel hat section similar to front rails of cars is represented by a nonlinear spring while the load cell is assumed to behave in a linear manner due to its high stiffness. Assuming a given load-displacement response obtained in an actual test as the true behavior of the specimen, the numerical solution of the governing differential equations following an implicit time integration scheme is shown to yield an excellent reproduction of the mechanical behavior of the specimen thereby confirming the accuracy of the numerical approach. The spring representing the load cell, however,predicts a response that qualitatively matches the assumed load-displacement response of the test specimen with a perceptibly lower magnitude of load.

Open Circuit Interaction of Formic Acid with Oxidized Pt Surfaces: Experiments, Modeling, and Simulations

We report in this work the study of the interaction between formic acid and an oxidized platinum surface under open circuit conditions. The investigation was carried out with the aid of in situ infrared spectroscopy, and results analyzed in terms of a mathematical model and numerical simulations. It has been found that during the first seconds of the interaction a small amount of CO(2) is produced and absolutely no adsorbed CO was observed. A sudden drop in potential then follows, which is accompanied by a steep increase first of CO(2) production and then by adsorbed CO. The steep transient was rationalized in terms of an autocatalytic production of free platinum sites which enhances the overall rate of reaction. Modeling and simulation showed nearly quantitative agreement with the experimental observations and provided further insight into some experimentally inaccessible variables such as surface free sites. Finally, based on the understanding provided from the combined experimental and theoretical approach, we discuss the general aspects influencing the open circuit transient.

Modeling and simulating a mathematical tool for multi-robot pattern transformation

The work reported in this paper is motivated by the need to investigate general methods for pattern transformation. A formal definition for pattern transformation is provided and four special cases namely, elementary and geometric transformation based on repositioning all and some agents in the pattern are introduced. The need for a mathematical tool and simulations for visualizing the behavior of a transformation method is highlighted. A mathematical method based on the Moebius transformation is proposed. The transformation method involves discretization of events for planning paths of individual robots in a pattern. Simulations on a particle physics simulator are used to validate the feasibility of the proposed method.

Modeling and simulating a mathematical tool for multi-robot pattern transformation

The work reported in this paper is motivated by the need to investigate general methods for pattern transformation. A formal definition for pattern transformation is provided and four special cases namely, elementary and geometric transformation based on repositioning all and some agents in the pattern are introduced. The need for a mathematical tool and simulations for visualizing the behavior of a transformation method is highlighted. A mathematical method based on the Moebius transformation is proposed. The transformation method involves discretization of events for planning paths of individual robots in a pattern. Simulations on a particle physics simulator are used to validate the feasibility of the proposed method.

Feature-oriented regional modeling and simulations (FORMS) for the western South Atlantic: Southeastern Brazil region

The multi-scale synoptic circulation system in the southeastern Brazil (SEBRA) region is presented using a feature-oriented approach. Prevalent synoptic circulation structures, or ""features,"" are identified from previous observational studies. These features include the southward-flowing Brazil Current (BC), the eddies off Cabo Sao Tome (CST - 22 degrees S) and off Cabo Frio (CF - 23 degrees S), and the upwelling region off CF and CST. Their synoptic water-mass (T-S) structures are characterized and parameterized to develop temperature-salinity (T-S) feature models. Following [Gangopadhyay, A., Robinson, A.R., Haley, PJ., Leslie, W.J., Lozano, C.j., Bisagni, J., Yu, Z., 2003. Feature-oriented regional modeling and simulation (forms) in the gulf of maine and georges bank. Cont. Shelf Res. 23 (3-4), 317-353] methodology, a synoptic initialization scheme for feature-oriented regional modeling and simulation (FORMS) of the circulation in this region is then developed. First, the temperature and salinity feature-model profiles are placed on a regional circulation template and objectively analyzed with available background climatology in the deep region. These initialization fields are then used for dynamical simulations via the Princeton Ocean Model (POM). A few first applications of this methodology are presented in this paper. These include the BC meandering, the BC-eddy interaction and the meander-eddy-upwelling system (MEUS) simulations. Preliminary validation results include realistic wave-growth and eddy formation and sustained upwelling. Our future plan includes the application of these feature models with satellite, in-situ data and advanced data-assimilation schemes for nowcasting and forecasting the SEBRA region. (c) 2008 Elsevier Ltd. All rights reserved.

The Boundary Element Method in elastodynamics

As is well known B.E.M. is obtained as a mixture of the integral representation formula of classical elasticity and the discretization philosophy of the finite element method (F.E.M.). The paper presents the application of B.E.M. to elastodynamic problems. Both the transient and steady state solutions are presented as well as some techniques to simplify problems with a free-stress boundary.

Formal modeling and analysis techniques for high level Petri nets

Petri Nets are a formal, graphical and executable modeling technique for the specification and analysis of concurrent and distributed systems and have been widely applied in computer science and many other engineering disciplines. Low level Petri nets are simple and useful for modeling control flows but not powerful enough to define data and system functionality. High level Petri nets (HLPNs) have been developed to support data and functionality definitions, such as using complex structured data as tokens and algebraic expressions as transition formulas. Compared to low level Petri nets, HLPNs result in compact system models that are easier to be understood. Therefore, HLPNs are more useful in modeling complex systems. ^ There are two issues in using HLPNs—modeling and analysis. Modeling concerns the abstracting and representing the systems under consideration using HLPNs, and analysis deals with effective ways study the behaviors and properties of the resulting HLPN models. In this dissertation, several modeling and analysis techniques for HLPNs are studied, which are integrated into a framework that is supported by a tool. ^ For modeling, this framework integrates two formal languages: a type of HLPNs called Predicate Transition Net (PrT Net) is used to model a system's behavior and a first-order linear time temporal logic (FOLTL) to specify the system's properties. The main contribution of this dissertation with regard to modeling is to develop a software tool to support the formal modeling capabilities in this framework. ^ For analysis, this framework combines three complementary techniques, simulation, explicit state model checking and bounded model checking (BMC). Simulation is a straightforward and speedy method, but only covers some execution paths in a HLPN model. Explicit state model checking covers all the execution paths but suffers from the state explosion problem. BMC is a tradeoff as it provides a certain level of coverage while more efficient than explicit state model checking. The main contribution of this dissertation with regard to analysis is adapting BMC to analyze HLPN models and integrating the three complementary analysis techniques in a software tool to support the formal analysis capabilities in this framework. ^ The SAMTools developed for this framework in this dissertation integrates three tools: PIPE+ for HLPNs behavioral modeling and simulation, SAMAT for hierarchical structural modeling and property specification, and PIPE+Verifier for behavioral verification.^

Formal Modeling and Analysis Techniques for High Level Petri Nets

Petri Nets are a formal, graphical and executable modeling technique for the specification and analysis of concurrent and distributed systems and have been widely applied in computer science and many other engineering disciplines. Low level Petri nets are simple and useful for modeling control flows but not powerful enough to define data and system functionality. High level Petri nets (HLPNs) have been developed to support data and functionality definitions, such as using complex structured data as tokens and algebraic expressions as transition formulas. Compared to low level Petri nets, HLPNs result in compact system models that are easier to be understood. Therefore, HLPNs are more useful in modeling complex systems. There are two issues in using HLPNs - modeling and analysis. Modeling concerns the abstracting and representing the systems under consideration using HLPNs, and analysis deals with effective ways study the behaviors and properties of the resulting HLPN models. In this dissertation, several modeling and analysis techniques for HLPNs are studied, which are integrated into a framework that is supported by a tool. For modeling, this framework integrates two formal languages: a type of HLPNs called Predicate Transition Net (PrT Net) is used to model a system's behavior and a first-order linear time temporal logic (FOLTL) to specify the system's properties. The main contribution of this dissertation with regard to modeling is to develop a software tool to support the formal modeling capabilities in this framework. For analysis, this framework combines three complementary techniques, simulation, explicit state model checking and bounded model checking (BMC). Simulation is a straightforward and speedy method, but only covers some execution paths in a HLPN model. Explicit state model checking covers all the execution paths but suffers from the state explosion problem. BMC is a tradeoff as it provides a certain level of coverage while more efficient than explicit state model checking. The main contribution of this dissertation with regard to analysis is adapting BMC to analyze HLPN models and integrating the three complementary analysis techniques in a software tool to support the formal analysis capabilities in this framework. The SAMTools developed for this framework in this dissertation integrates three tools: PIPE+ for HLPNs behavioral modeling and simulation, SAMAT for hierarchical structural modeling and property specification, and PIPE+Verifier for behavioral verification.