Slabsim: Simulador de dispositivos ópticos integrados

En los diseños y desarrollos de ingeniería, antes de comenzar la construcción e implementación de los objetivos de un proyecto, es necesario realizar una serie de análisis previos y simulaciones que corroboren las expectativas de la hipótesis inicial, con el fin de obtener una referencia empírica que satisfaga las condiciones de trabajo o funcionamiento de los objetivos de dicho proyecto. A menudo, los resultados que satisfacen las características deseadas se obtienen mediante la iteración de métodos de ensayo y error. Generalmente, éstos métodos utilizan el mismo procedimiento de análisis con la variación de una serie de parámetros que permiten adaptar una tecnología a la finalidad deseada. Hoy en día se dispone de computadoras potentes, así como algoritmos de resolución matemática que permiten resolver de forma veloz y eficiente diferentes tipos de problemas de cálculo. Resulta interesante el desarrollo de aplicaciones que permiten la resolución de éstos problemas de forma rápida y precisa en el análisis y síntesis de soluciones de ingeniería, especialmente cuando se tratan expresiones similares con variaciones de constantes, dado que se pueden desarrollar instrucciones de resolución con la capacidad de inserción de parámetros que definan el problema. Además, mediante la implementación de un código de acuerdo a la base teórica de una tecnología, se puede lograr un código válido para el estudio de cualquier problema relacionado con dicha tecnología. El desarrollo del presente proyecto pretende implementar la primera fase del simulador de dispositivos ópticos Slabsim, en cual se puede representar la distribución de la energía de una onda electromagnética en frecuencias ópticas guiada a través de una una guía dieléctrica plana, también conocida como slab. Este simulador esta constituido por una interfaz gráfica generada con el entorno de desarrollo de interfaces gráficas de usuario Matlab GUIDE, propiedad de Mathworks©, de forma que su manejo resulte sencillo e intuitivo para la ejecución de simulaciones con un bajo conocimiento de la base teórica de este tipo de estructuras por parte del usuario. De este modo se logra que el ingeniero requiera menor intervalo de tiempo para encontrar una solución que satisfaga los requisitos de un proyecto relacionado con las guías dieléctricas planas, e incluso utilizarlo para una amplia diversidad de objetivos basados en esta tecnología. Uno de los principales objetivos de este proyecto es la resolución de la base teórica de las guías slab a partir de métodos numéricos computacionales, cuyos procedimientos son extrapolables a otros problemas matemáticos y ofrecen al autor una contundente base conceptual de los mismos. Por este motivo, las resoluciones de las ecuaciones diferenciales y características que constituyen los problemas de este tipo de estructuras se realizan por estos medios de cálculo en el núcleo de la aplicación, dado que en algunos casos, no existe la alternativa de uso de expresiones analíticas útiles. ABSTRACT. The first step in engineering design and development is an analysis and simulation process which will successfully corroborate the initial hypothesis that was made and find solutions for a particular. In this way, it is possible to obtain empirical evidence which suitably substantiate the purposes of the project. Commonly, the characteristics to reach a particular target are found through iterative trial and error methods. These kinds of methods are based on the same theoretical analysis but with a variation of some parameters, with the objective to adapt the results for a particular aim. At present, powerful computers and mathematical algorithms are available to solve different kinds of calculation problems in a fast and efficient way. Computing application development is useful as it gives a high level of accurate results for engineering analysis and synthesis in short periods of time. This is more notable in cases where the mathematical expressions on a theoretical base are similar but with small variations of constant values. This is due to the ease of adaptation of the computer programming code into a parameter request system that defines a particular solution on each execution. Additionally, it is possible to code an application suitable to simulate any issue related to the studied technology. The aim of the present project consists of the construction of the first stage of an optoelectronics simulator named Slabsim. Slabism is capable of representing the energetic distribution of a light wave guided in the volume of a slab waveguide. The mentioned simulator is made through the graphic user interface development environment Matlab GUIDE, property of Mathworks©. It is designed for an easy and intuitive management by the user to execute simulations with a low knowledge of the technology theoretical bases. With this software it is possible to achieve several aims related to the slab waveguides by the user in low interval of time. One of the main purposes of this project is the mathematical solving of theoretical bases of slab structures through computing numerical analysis. This is due to the capability of adapting its criterion to other mathematical issues and provides a strong knowledge of its process. Based on these advantages, numerical solving methods are used in the core of the simulator to obtain differential and characteristic equations results that become represented on it.

Some aspects of lubrication in heavy regimes: thermal effects, stability and turbulence

In this work, a combination of numerical methods applied to thermohydrodynamic lubrication problems with cavitation is presented. It should be emphasized the difficulty of the nonlinear mathematical coupled model involving a free boundary problem, but also the simplicity of the algorithms employed to solve it. So, finite element discretizations for the hydrodynamic and thermal equations combined with upwind techniques for the convection terms and duality methods for nonlinear features are proposed. Additionally, a model describing the movement of the shaft is provided. Considering the shaft as a rigid body this model will consist of an ODE system relating acceleration of the center of gravity and external and pressure loads. The numerical experiments of mechanical stability try to clarify the position of the neutral stability curve. Finally, a rotating machine for ship propulsion involving both axial and radial bearings operating with nonconventional lubricants (seawater to avoid environmental pollution) is analyzed by using laminar and turbulent inertial flows.

Simulation and geometrical design of multi-section tapered semiconductor optical amplifiers at 1.57 µm

Fully integrated semiconductor master-oscillator power-amplifiers (MOPA) with a tapered power amplifier are attractive sources for applications requiring high brightness. The geometrical design of the tapered amplifier is crucial to achieve the required power and beam quality. In this work we investigate by numerical simulation the role of the geometrical design in the beam quality and in the maximum achievable power. The simulations were performed with a Quasi-3D model which solves the complete steady-state semiconductor and thermal equations combined with a beam propagation method. The results indicate that large devices with wide taper angles produce higher power with better beam quality than smaller area designs, but at expenses of a higher injection current and lower conversion efficiency.

Characterization of genes modulated during pheomelanogenesis using differential display

Molecular and biochemical mechanisms that modulate the production of eumelanin or pheomelanin pigments involve the opposing effects of two intercellular signaling molecules, α-melanocyte stimulating hormone (MSH) and agouti signal protein (ASP). ASP is an antagonist of MSH signaling through the melanocyte-specific MSH receptor, although its mechanism(s) of action is controversial. We previously have reported significant down-regulation of all known melanogenic genes during the eumelanin to pheomelanin switch in murine hair follicle melanocytes and in cultured melanocytes treated with recombinant ASP. To identify factors that might be involved in the switch to pheomelanogenesis, we screened ASP-treated melanocytes by using differential display and identified three up-regulated genes: a DNA replication control protein, a basic helix–loop–helix transcription factor, and a novel gene. We have simultaneously identified six down-regulated genes in ASP-treated melanocytes; two of those encode tyrosinase and TRP2, melanogenic genes known to be down-regulated during pheomelanogenesis, which provide good internal controls for this approach. These results suggest that there are complex mechanisms involved in the switch to pheomelanin production, and that these modulated genes might be involved in the pleiotropic changes seen in yellow mice, including the change in coat color.

Localization and Recycling of gp27 (hp24γ3): Complex Formation with Other p24 Family Members

We report here the characterization of gp27 (hp24γ3), a glycoprotein of the p24 family of small and abundant transmembrane proteins of the secretory pathway. Immunoelectron and confocal scanning microscopy show that at steady state, gp27 localizes to the cis side of the Golgi apparatus. In addition, some gp27 was detected in COPI- and COPII-coated structures throughout the cytoplasm. This indicated cycling that was confirmed in three ways. First, 15°C temperature treatment resulted in accumulation of gp27 in pre-Golgi structures colocalizing with anterograde cargo. Second, treatment with brefeldin A caused gp27 to relocate into peripheral structures positive for both KDEL receptor and COPII. Third, microinjection of a dominant negative mutant of Sar1p trapped gp27 in the endoplasmic reticulum (ER) by blocking ER export. Together, this shows that gp27 cycles extensively in the early secretory pathway. Immunoprecipitation and coexpression studies further revealed that a significant fraction of gp27 existed in a hetero-oligomeric complex. Three members of the p24 family, GMP25 (hp24α2), p24 (hp24β1), and p23 (hp24δ1), coprecipitated in what appeared to be stochiometric amounts. This heterocomplex was specific. Immunoprecipitation of p26 (hp24γ4) failed to coprecipitate GMP25, p24, or p23. Also, very little p26 was found coprecipitating with gp27. A functional requirement for complex formation was suggested at the level of ER export. Transiently expressed gp27 failed to leave the ER unless other p24 family proteins were coexpressed. Comparison of attached oligosaccharides showed that gp27 and GMP25 recycled differentially. Only a very minor portion of GMP25 displayed complex oligosaccharides. In contrast, all of gp27 showed modifications by medial and trans enzymes at steady state. We conclude from these data that a portion of gp27 exists as hetero-oligomeric complexes with GMP25, p24, and p23 and that these complexes are in dynamic equilibrium with individual p24 proteins to allow for differential recycling and distributions.

Isolation and characterization of a cell surface albumin-binding protein from vascular endothelial cells.

Albumin-binding proteins identified in vascular endothelial cells have been postulated to contribute to the transport of albumin via a process involving transcytosis. In the present study, we have purified and characterized a 57- to 60-kDa (gp60) putative albumin-binding protein from bovine pulmonary microvessel endothelial cells. The endothelial cell membranes were isolated from cultured cells by differential centrifugation and solubilized with sodium cholate and urea. The solubilized extract was concentrated after dialysis by ethanol precipitation and reextracted with Triton X-100, and the resulting extract was subjected to DEAE-cellulose column chromatography. Proteins eluted from this column were further separated using preparative sodium dodecyl sulfate/polyacrylamide gel electrophoresis and used for immunizing rabbits. Fluorescence-activated cell sorter analysis using the anti-gp60 antibodies demonstrated the expression of gp60 on the endothelial cell surface. Affinity-purified anti-gp60 antibodies inhibited approximately 90% of the specific binding of 125I-labeled albumin to bovine pulmonary microvessel endothelial cell surface. The anti-gp60 antibodies reacted with gp60 from bovine pulmonary artery, bovine pulmonary microvessel, human umbilical vein, and rat lung endothelial cell membranes. Bovine anti-gp60 antibodies also reacted with bovine secreted protein, acidic and rich in cysteine (SPARC). However, bovine SPARC NH2-terminal sequence (1-56 residues) antibodies did not react with gp60, indicating that the endothelial cell-surface-associated albumin-binding protein gp60 was different from the secreted albumin-binding protein SPARC. We conclude that the endothelial cell-surface-associated gp60 mediates the specific binding of native albumin to endothelial cells and thus may regulate the uptake of albumin and its transcytosis.

Discovery of adrenomedullin in rat ischemic cortex and evidence for its role in exacerbating focal brain ischemic damage.

Focal brain ischemia is the most common event leading to stroke in humans. To understand the molecular mechanisms associated with brain ischemia, we applied the technique of mRNA differential display and isolated a gene that encodes a recently discovered peptide, adrenomedullin (AM), which is a member of the calcitonin gene-related peptide (CGRP) family. Using the rat focal stroke model of middle cerebral artery occlusion (MCAO), we determined that AM mRNA expression was significantly increased in the ischemic cortex up to 17.4-fold at 3 h post-MCAO (P < 0.05) and 21.7-fold at 6 h post-MCAO (P < 0.05) and remained elevated for up to 15 days (9.6-fold increase; P < 0.05). Immunohistochemical studies localized AM to ischemic neuronal processes, and radioligand (125I-labeled CGRP) displacement revealed high-affinity (IC50 = 80.3 nmol) binding of AM to CGRP receptors in brain cortex. The cerebrovascular function of AM was studied using synthetic AM microinjected onto rat pial vessels using a cranial window or applied to canine basilar arteries in vitro. AM, applied abluminally, produced dose-dependent relaxation of preconstricted pial vessels (P < 0.05). Intracerebroventricular (but not systemic) AM administration at a high dose (8 nmol), prior to and after MCAO, increased the degree of focal ischemic injury (P < 0.05). The ischemia-induced expression of both AM mRNA and peptide in ischemic cortical neurons, the demonstration of the direct vasodilating effects of the peptide on cerebral vessels, and the ability of AM to exacerbate ischemic brain damage suggests that AM plays a significant role in focal ischemic brain injury.

Mutagenesis of palmitoylation sites in endothelial nitric oxide synthase identifies a novel motif for dual acylation and subcellular targeting.

The endothelial nitric oxide synthase (ec-NOS) plays a key role in the transduction of signals from the bloodstream to the underlying smooth muscle. ecNOS undergoes a complex series of covalent modifications, including myristoylation and palmitoylation, which appear to play a role in ecNOS membrane association. Mutagenesis of the myristoylation site, which prevents both myristoylation and palmitoylation, blocks ecNOS targeting to cell membranes. Further, as described for some G-protein alpha subunits, both membrane association and palmitoylation of ecNOS are dynamically regulated: in response to agonists, the enzyme undergoes partial redistribution to the cell cytosol concomitant with depalmitoylation. To clarify the role of palmitoylation in determining ecNOS subcellular localization, we have constructed palmitoylation-deficient mutants of ecNOS. Serine was substituted for cysteine at two potential palmitoylation sites (Cys-15 and Cys-26) by site-directed mutagenesis. Immunoprecipitation of ecNOS mutants following cDNA transfection and biosynthetic labeling with [3H]palmitate revealed that mutagenesis of either cysteine residue attenuated palmitoylation, whereas replacement of both residues completely eliminated palmitoylation. Analysis of N-terminal deletion mutations of ecNOS demonstrated that the region containing these two cysteine residues is both necessary and sufficient for enzyme palmitoylation. The cysteines thus identified as the palmitoylation sites for ecNOS are separated by an unusual (Gly-Leu)5 sequence and appear to define a sequence motif for dual acylation. We analyzed the subcellular distribution of ecNOS mutants by differential ultracentrifugation and found that mutagenesis of the ecNOS palmitoylation sites markedly reduced membrane association of the enzyme. These results document that ecNOS palmitoylation is an important determinant for the subcellular distribution of ecNOS and identify a new motif for the reversible palmitoylation of signaling proteins.

Restriction fragment length polymorphism-coupled domain-directed differential display: a highly efficient technique for expression analysis of multigene families.

In this paper, a reverse-transcriptase PCR-based protocol suitable for efficient expression analysis of multigene families is presented. The method combines restriction fragment length polymorphism (RFLP) technology with a gene family-specific version of mRNA differential display and hence is called "RFLP-coupled domain-directed differential display. "With this method, expression of all members of a multigene family at many different developmental stages, in diverse tissues and even in different organisms, can be displayed on one gel. Moreover, bands of interest, representing gene family members, are directly accessible to sequence analysis, without the need for subcloning. The method thus enables a detailed, high-resolution expression analysis of known gene family members as well as the identification and characterization of new ones. Here the technique was used to analyze differential expression of MADS-box genes in male and female inflorescences of maize (Zea mays ssp. mays). Six different MADS-box genes could be identified, being either specifically expressed in the female sex or preferentially expressed in male or female inflorescences, respectively. Other possible applications of the method are discussed.

Real time optimization in chemical process: evaluation of strategies, improvements and industrial application.

The increasing economic competition drives the industry to implement tools that improve their processes efficiencies. The process automation is one of these tools, and the Real Time Optimization (RTO) is an automation methodology that considers economic aspects to update the process control in accordance with market prices and disturbances. Basically, RTO uses a steady-state phenomenological model to predict the process behavior, and then, optimizes an economic objective function subject to this model. Although largely implemented in industry, there is not a general agreement about the benefits of implementing RTO due to some limitations discussed in the present work: structural plant/model mismatch, identifiability issues and low frequency of set points update. Some alternative RTO approaches have been proposed in literature to handle the problem of structural plant/model mismatch. However, there is not a sensible comparison evaluating the scope and limitations of these RTO approaches under different aspects. For this reason, the classical two-step method is compared to more recently derivative-based methods (Modifier Adaptation, Integrated System Optimization and Parameter estimation, and Sufficient Conditions of Feasibility and Optimality) using a Monte Carlo methodology. The results of this comparison show that the classical RTO method is consistent, providing a model flexible enough to represent the process topology, a parameter estimation method appropriate to handle measurement noise characteristics and a method to improve the sample information quality. At each iteration, the RTO methodology updates some key parameter of the model, where it is possible to observe identifiability issues caused by lack of measurements and measurement noise, resulting in bad prediction ability. Therefore, four different parameter estimation approaches (Rotational Discrimination, Automatic Selection and Parameter estimation, Reparametrization via Differential Geometry and classical nonlinear Least Square) are evaluated with respect to their prediction accuracy, robustness and speed. The results show that the Rotational Discrimination method is the most suitable to be implemented in a RTO framework, since it requires less a priori information, it is simple to be implemented and avoid the overfitting caused by the Least Square method. The third RTO drawback discussed in the present thesis is the low frequency of set points update, this problem increases the period in which the process operates at suboptimum conditions. An alternative to handle this problem is proposed in this thesis, by integrating the classic RTO and Self-Optimizing control (SOC) using a new Model Predictive Control strategy. The new approach demonstrates that it is possible to reduce the problem of low frequency of set points updates, improving the economic performance. Finally, the practical aspects of the RTO implementation are carried out in an industrial case study, a Vapor Recompression Distillation (VRD) process located in Paulínea refinery from Petrobras. The conclusions of this study suggest that the model parameters are successfully estimated by the Rotational Discrimination method; the RTO is able to improve the process profit in about 3%, equivalent to 2 million dollars per year; and the integration of SOC and RTO may be an interesting control alternative for the VRD process.

On the development of advanced techniques for mixed-elastohydrodynamic lubrification modelling of journal and sliding bearing systems.

The present thesis is focused on the development of a thorough mathematical modelling and computational solution framework aimed at the numerical simulation of journal and sliding bearing systems operating under a wide range of lubrication regimes (mixed, elastohydrodynamic and full film lubrication regimes) and working conditions (static, quasi-static and transient conditions). The fluid flow effects have been considered in terms of the Isothermal Generalized Equation of the Mechanics of the Viscous Thin Films (Reynolds equation), along with the massconserving p-Ø Elrod-Adams cavitation model that accordingly ensures the so-called JFO complementary boundary conditions for fluid film rupture. The variation of the lubricant rheological properties due to the viscous-pressure (Barus and Roelands equations), viscous-shear-thinning (Eyring and Carreau-Yasuda equations) and density-pressure (Dowson-Higginson equation) relationships have also been taken into account in the overall modelling. Generic models have been derived for the aforementioned bearing components in order to enable their applications in general multibody dynamic systems (MDS), and by including the effects of angular misalignments, superficial geometric defects (form/waviness deviations, EHL deformations, etc.) and axial motion. The bearing exibility (conformal EHL) has been incorporated by means of FEM model reduction (or condensation) techniques. The macroscopic in fluence of the mixedlubrication phenomena have been included into the modelling by the stochastic Patir and Cheng average ow model and the Greenwood-Williamson/Greenwood-Tripp formulations for rough contacts. Furthermore, a deterministic mixed-lubrication model with inter-asperity cavitation has also been proposed for full-scale simulations in the microscopic (roughness) level. According to the extensive mathematical modelling background established, three significant contributions have been accomplished. Firstly, a general numerical solution for the Reynolds lubrication equation with the mass-conserving p - Ø cavitation model has been developed based on the hybridtype Element-Based Finite Volume Method (EbFVM). This new solution scheme allows solving lubrication problems with complex geometries to be discretized by unstructured grids. The numerical method was validated in agreement with several example cases from the literature, and further used in numerical experiments to explore its exibility in coping with irregular meshes for reducing the number of nodes required in the solution of textured sliding bearings. Secondly, novel robust partitioned techniques, namely: Fixed Point Gauss-Seidel Method (PGMF), Point Gauss-Seidel Method with Aitken Acceleration (PGMA) and Interface Quasi-Newton Method with Inverse Jacobian from Least-Squares approximation (IQN-ILS), commonly adopted for solving uid-structure interaction problems have been introduced in the context of tribological simulations, particularly for the coupled calculation of dynamic conformal EHL contacts. The performance of such partitioned methods was evaluated according to simulations of dynamically loaded connecting-rod big-end bearings of both heavy-duty and high-speed engines. Finally, the proposed deterministic mixed-lubrication modelling was applied to investigate the in fluence of the cylinder liner wear after a 100h dynamometer engine test on the hydrodynamic pressure generation and friction of Twin-Land Oil Control Rings.

Análise da abrangência de modelo modificado para mancais curtos com deformação.

Este trabalho apresenta uma discussão sobre o estudo dos efeitos térmicos e elásticos decorrentes da pressão de sustentação presentes nos mancais. Para tanto, propõe-se um modelo matemático baseado nas equações para mancais curtos considerando a região de cavitação e utilizando o princípio da continuidade de massa. Com isto, deduzem-se as equações para o mancal a partir das equações de Reynolds e da energia, aplicando uma solução modificada para a solução de Ocvirk, sendo as equações resolvidas numericamente pelo Método das Diferenças Finitas. Somado o tratamento de mecânica dos fluidos, o trabalho discute dois modelos térmicos de previsão de temperatura média do fluido e sua influência no campo de pressão, apresentando gráficos representativos do campo de pressão e de temperatura, assim como as diferenças e implicações das diferenças. Para o cálculo de deformação da estrutura, utiliza-se um Modelo de Elementos Finitos para uma dada geometria, fazendo-se uma avaliação da variação do campo de pressão e o quanto essa diferença afeta as demais propriedades do fluido. Por fim, com o modelo completo, calcula-se o quanto esse modelamento para mancais curtos se aproxima de soluções para mancais finitos, com base em resultados da literatura, chegando a desvios quase oito vezes menores que os previstos pela literatura. Além disso, pode-se estabelecer a abrangência do modelo, ou seja, prever as condições em que suas propriedades são válidas e podem ser utilizadas para estudos iniciais.

No anticorrelation between cyclotron line energy and X-ray flux in 4U 0115+634

We report on an outburst of the high mass X-ray binary 4U 0115+634 with a pulse period of 3.6 s in 2008 March/April as observed with RXTE and INTEGRAL. During the outburst the neutron star’s luminosity varied by a factor of 10 in the 3–50 keV band. In agreement with earlier work we find evidence of five cyclotron resonance scattering features at ~10.7, 21.8, 35.5, 46.7, and 59.7 keV. Previous work had found an anticorrelation between the fundamental cyclotron line energy and the X-ray flux. We show that this apparent anticorrelation is probably due to the unphysical interplay of parameters of the cyclotron line with the continuum models used previously, e.g., the negative and positive exponent power law (NPEX). For this model, we show that cyclotron line modeling erroneously leads to describing part of the exponential cutoff and the continuum variability, and not the cyclotron lines. When the X-ray continuum is modeled with a simple exponentially cutoff power law modified by a Gaussian emission feature around 10 keV, the correlation between the line energy and the flux vanishes, and the line parameters remain virtually constant over the outburst. We therefore conclude that the previously reported anticorrelation is an artifact of the assumptions adopted in the modeling of the continuum.

Suspended matter on the meridional section in the Laptev Sea in August-September 1991

Suspended matter concentration along a meridian section from the Lena River delta to 78°N (~500 km) at ten stations from the surface to the bottom was studied with weight and optical (light attenuation index) techniques. At seven stations average residence time of suspended matter in surface waters was determined by the disequilibrium 234Th method. Average residence time of suspended matter in other depth intervals was calculated by regression dependence between the 234Th/238U ratio and suspended matter concentration. Differential and integral fluxes of suspended matter in the water column were also calculated. Nepheloid matter dominates in suspended matter composition in surface waters. Calculations indicate that, before being buried on the bottom, solid river run-off is resuspended 2.3 times (aver.). Redistribution of nepheloid suspended matter in the near-bottom layer results in formation of a strongly pronounced depocentre - an area of maximal discharge of solid river run-off within the Laptev Sea.

Introduction to the calculus,

"A revision of the author's A first course in the differential and integral calculus."--Pref.