基于Berreman矩阵的双折射薄膜光谱响应特性计算

采用Berreman特征矩阵方法，通过数值计算研究了双折射薄膜的反射、透射等光谱响应特性。依据电磁场理论的电场分量、磁场分量的界面连续条件，推导了光波在各向异性双轴薄膜中的Berreman转移矩阵，用以分析含有各向异性介质层的复杂薄膜系统的光学性质。这些矩阵递推关系包含了界面处的多点反射，适用于一般的各向异性的多层膜系统，包括入射媒质或基底为各向异性的情况。在文中给出了各向同性入射媒质双轴各向异性膜层一各向同性基底薄膜系统的计算结果，验证了该计算方法的可行性，以此作为进一步研究各向异性薄膜和相关光学薄膜器

光学薄膜中界面层和表面吸附层对相位延迟的影响

多层介质反射镜在非正入射的时候，两个不同的偏振态之间会产生不同的相移．根据空气与膜层、膜层之间的实际情况，建立了界面层和表面吸附层模型，并运用它分析相位延迟产生误差的原因．通过优化设计，入射角为54°，在1285～1345nm之间p，s波获得了270±1°的相移，同时也使反射率在99．5％以上．用离子束溅射技术制备相位延迟膜，用分光光度计测试了光谱特性和用椭偏仪测试了相位特性，在相应波段获得了262．4±1．8°的相移，同时也使反射率在99．6％以上．误差的主要来源是离子源工作特性会产生不均匀的过渡层和最

Mapeamento das regiões radiantes em placas retangulares vibrantes via método dos elementos de contorno com velocidade média e intensidade útil

Em engenharia, a modelagem computacional desempenha um papel importante na concepção de produtos e no desenvolvimento de técnicas de atenuação de ruído. Nesse contexto, esta tese investiga a intensidade acústica gerada pela radiação sonora de superfícies vibrantes. De modo específico, a pesquisa enfoca a identificação das regiões de uma fonte sonora que contribuem efetivamente para potência sonora radiada para o campo afastado, quando a frequência de excitação ocorre abaixo da frequência crítica de coincidência. São descritas as fundamentações teóricas de duas diferentes abordagens. A primeira delas, denominada intensidade supersônica (analítica) é calculada via transformadas de Fourier para fontes sonoras com geometrias separáveis. A segunda, denominada intensidade útil (numérica) é calculada através do método dos elementos de contorno clássico para fontes com geometrias arbitrárias. Em ambas, a identificação das regiões é feita pela filtragem das ondas não propagantes (evanescentes). O trabalho está centrado em duas propostas, a saber. A primeira delas, é a apresentação implementação e análise de uma nova técnica numérica para o cálculo da grandeza intensidade útil. Essa técnica constitui uma variante do método dos elementos de contorno (MEC), tendo como base o fato de as aproximações para as variáveis acústicas pressão e velocidade normal serem tomadas como constantes em cada elemento. E também no modo peculiar de obter a velocidade constante através da média de um certo número de velocidades interiores a cada elemento. Por esse motivo, a técnica recebe o nome de método de elemento de contorno com velocidade média (AVBEMAverage Velocity Boundary Element Method). A segunda, é a obtenção da solução forma fechada do campo de velocidade normal para placas retangulares com oito diferentes combinações de condições contorno clássicas. Então, a intensidade supersônica é estimada e comparada à intensidade acústica. Nos ensaios numéricos, a comparação da intensidade útil obtida via MEC clássico e via AVBEM é mostrada para ilustrar a eficiência computacional da técnica aqui proposta, que traz como benefício adicional o fato de poder ser utilizada uma malha menos refinada para as simulações e, consequentemente, economia significativa de recursos computacionais.

Condições de contorno tipo albedo para cálculos globais de reatores nucleares com o modelo de dois grupos de energia na formulação de transporte SN em geometria bidimensional cartesiana

Os eventos de fissão nuclear, resultados da interação dos nêutrons com os núcleos dos átomos do meio hospedeiro multiplicativo, não estão presentes em algumas regiões dos reatores nucleares, e.g., moderador, refletor, e meios estruturais. Nesses domínios espaciais não há geração de potência nuclear térmica e, além disso, comprometem a eficiência computacional dos cálculos globais de reatores nucleares. Propomos nesta tese uma estratégia visando a aumentar a eficiência computacional dessas simulações eliminando os cálculos numéricos explícitos no interior das regiões não-multiplicativas (baffle e refletor) em torno do núcleo ativo. Apresentamos algumas modelagens e discutimos a eficiência da aplicação dessas condições de contorno aproximadas tipo albedo para uma e duas regiões nãomultiplicativas, na formulação de ordenadas discretas (SN) para problemas de autovalor a dois grupos de energia em geometria bidimensional cartesiana. A denominação Albedo, palavra de origem latina para alvura, foi originalmente definida como a fração da luz incidente que é refletida difusamente por uma superfície. Esta denominação latina permaneceu como o termo científico usual em astronomia e, nesta tese, este conceito é estendido para reflexão de nêutrons. Estas condições de contorno tipo albedo SN não-convencional substituem aproximadamente as regiões de baffle e refletor no em torno do núcleo ativo do reator, desprezando os termos de fuga transversal no interior dessas regiões. Se o problema, em particular, não possui termos de fuga transversal, i.e., trata-se de um problema unidimensional, então as condições de contorno albedo, como propostas nesta tese, são exatas. Por eficiência computacional entende-se a análise da precisão dos resultados numéricos em comparação com o tempo de execução computacional de cada simulação de um dado problema-modelo. Resultados numéricos considerando dois problemas-modelo com de simetria são considerados para ilustrar esta análise de eficiência.

A discrete dislocation analysis of fatigue crack growth

Cyclic loading of a plane strain mode I crack under small scale yielding is analyzed using discrete dislocation dynamics. The dislocations are all of edge character, and are modeled as line singularities in an elastic solid. At each stage of loading, superposition is used to represent the solution in terms of solutions for edge dislocations in a half-space and a non-singular complementary solution that enforces the boundary conditions, which is obtained from a linear elastic, finite element solution. The lattice resistance to dislocation motion, dislocation nucleation, dislocation interaction with obstacles and dislocation annihilation are incorporated into the formulation through a set of constitutive rules. An irreversible relation between the opening traction and the displacement jump across a cohesive surface ahead of the initial crack tip is also specified, which permits crack growth to emerge naturally. It is found that crack growth can occur under cyclic loading conditions even when the peak stress intensity factor is smaller than the stress intensity required for crack growth under monotonic loading conditions; however below a certain threshold value of ΔKI no crack growth was seen.

On the use of combined compact scheme for numerical solution of the two-layer oceanic shallow water equations

Many types of oceanic physical phenomena have a wide range in both space and time. In general, simplified models, such as shallow water model, are used to describe these oceanic motions. The shallow water equations are widely applied in various oceanic and atmospheric extents. By using the two-layer shallow water equations, the stratification effects can be considered too. In this research, the sixth-order combined compact method is investigated and numerically implemented as a high-order method to solve the two-layer shallow water equations. The second-order centered, fourth-order compact and sixth-order super compact finite difference methods are also used to spatial differencing of the equations. The first part of the present work is devoted to accuracy assessment of the sixth-order super compact finite difference method (SCFDM) and the sixth-order combined compact finite difference method (CCFDM) for spatial differencing of the linearized two-layer shallow water equations on the Arakawa's A-E and Randall's Z numerical grids. Two general discrete dispersion relations on different numerical grids, for inertia-gravity and Rossby waves, are derived. These general relations can be used for evaluation of the performance of any desired numerical scheme. For both inertia-gravity and Rossby waves, minimum error generally occurs on Z grid using either the sixth-order SCFDM or CCFDM methods. For the Randall's Z grid, the sixth-order CCFDM exhibits a substantial improvement , for the frequency of the barotropic and baroclinic modes of the linear inertia-gravity waves of the two layer shallow water model, over the sixth-order SCFDM. For the Rossby waves, the sixth-order SCFDM shows improvement, for the barotropic and baroclinic modes, over the sixth-order CCFDM method except on Arakawa's C grid. In the second part of the present work, the sixth-order CCFDM method is used to solve the one-layer and two-layer shallow water equations in their nonlinear form. In one-layer model with periodic boundaries, the performance of the methods for mass conservation is compared. The results show high accuracy of the sixth-order CCFDM method to simulate a complex flow field. Furthermore, to evaluate the performance of the method in a non-periodic domain the sixth-order CCFDM is applied to spatial differencing of vorticity-divergence-mass representation of one-layer shallow water equations to solve a wind-driven current problem with no-slip boundary conditions. The results show good agreement with published works. Finally, the performance of different schemes for spatial differencing of two-layer shallow water equations on Z grid with periodic boundaries is investigated. Results illustrate the high accuracy of combined compact method.

Breakdown voltage for superjunction power devices with charge imbalance: An analytical model valid for both punch through and non punch through devices

An analytical model for the electric field and the breakdown voltage (BV) of an unbalanced superjunction (SJ) device is presented in this paper. The analytical technique uses a superposition approach treating the asymmetric charge in the pillars as an excess charge component superimposed on a balanced charge component. The proposed double-exponentialmodel is able to accurately predict the electric field and the BV for unbalanced SJ devices in both punch through and non punch through conditions. The model is also reasonably accurate at extremely high levels of charge imbalance when the devices behave similarly to a PiN diode or to a high-conductance layer. The analytical model is compared against numerical simulations of charge unbalanced SJ devices and against experimental results. © 2009 IEEE.

Finite element analysis of compressor disc heat transfer using the transient conduction method with sparse boundary collocation points

Surface temperature measurements from two discs of a gas turbine compressor rig are used as boundary conditions for the transient conduction solution (inverse heat transfer analysis). The disc geometry is complex, and so the finite element method is used. There are often large radial temperature gradients on the discs, and the equations are therefore solved taking into account the dependence of thermal conductivity on temperature. The solution technique also makes use of a multigrid algorithm to reduce the solution time. This is particularly important since a large amount of data must be analyzed to obtain correlations of the heat transfer. The finite element grid is also used for a network analysis to calculate the radiant heat transfer in the cavity formed between the two compressor discs. The work discussed here proved particularly challenging as the disc temperatures were only measured at four different radial locations. Four methods of surface temperature interpolation are examined, together with their effect on the local heat fluxes. It is found that the choice of interpolation method depends on the available number of data points. Bessel interpolation gives the best results for four data points, whereas cubic splines are preferred when there are considerably more data points. The results from the analysis of the compressor rig data show that the heat transfer near the disc inner radius appears to be influenced by the central throughflow. However, for larger radii, the heat transfer from the discs and peripheral shroud is found to be consistent with that of a buoyancy-induced flow.

An Arbitrary Lagrangian Eulerian Method for Three-Phase Flows with Triple Junction Points

We present a moving mesh method suitable for solving two-dimensional and axisymmetric three-liquid flows with triple junction points. This method employs a body-fitted unstructured mesh where the interfaces between liquids are lines of the mesh system, and the triple junction points (if exist) are mesh nodes. To enhance the accuracy and the efficiency of the method, the mesh is constantly adapted to the evolution of the interfaces by refining and coarsening the mesh locally; dynamic boundary conditions on interfaces, in particular the triple points, are therefore incorporated naturally and accurately in a Finite- Element formulation. In order to allow pressure discontinuity across interfaces, double-values of pressure are necessary for interface nodes and triple-values of pressure on triple junction points. The resulting non-linear system of mass and momentum conservation is then solved by an Uzawa method, with the zero resultant condition on triple points reinforced at each time step. The method is used to investigate the rising of a liquid drop with an attached bubble in a lighter liquid.

Numerical simulation of post-flame oxidation of hydrocarbons in spark ignition engines

About 50-90 percent of the hydrocarbons that escape combustion during flame passage in spark-ignition engine operation are oxidized in the cylinder before leaving the system. The process involves the transport of unreacted fuel from cold walls towards the hotter burned gas regions and subsequent reaction. In order to understand controlling factors in the process, a transient one-dimensional reactive-diffusive model has been formulated for simulating the oxidation processes taking place in the reactive layer between hot burned gases and cold unreacted air/fuel mixture, with initial and boundary conditions provided by the emergence of hydrocarbons from the piston top land crevice. Energy and species conservation equations are solved for the entire process, using a detailed chemical kinetic mechanism for propane. Simulation results show that the post-flame oxidation process takes place within a reactive layer where intermediate hydrocarbon products are formed at temperatures above 1100-1200 K, followed by a carbon monoxide conversion region closer to the hot burned gases. Model results show that most of hydrocarbons leaving the crevice are completely oxidized inside the cylinder. The largest contribution of remaining hydrocarbons are those leaving the crevice at temperatures below 1400 K. The largest fraction of non-fuel (intermediate) hydrocarbons results from hydrocarbons leaving the crevice when core temperatures are around 1400 K Copyright © 1997 Society of Automotive Engineers, Inc.

Simulations of thermal upgrading methods for oil shale reservoirs

This paper analyzes reaction and thermal front development in porous reservoirs with reacting flows, such as those encountered in shale oil extraction. A set of dimensionless parameters and a 3D code are developed in order to investigate the important physical and chemical variables of such reservoirs when heated by in situ methods. This contribution builds on a 1D model developed for the precursor study to this work. Theory necessary for this study is presented, namely shale decomposition chemical mechanisms, governing equations for multiphase flow in porous media and necessary closure models. Plotting the ratio of the thermal wave speed to the fluid speed allows one to infer that the reaction wave front ends where this ratio is at a minimum. The reaction front follows the thermal front closely, thus allowing assumptions to be made about the extent of decomposition solely by looking at thermal wave progression. Furthermore, this sensitivity analysis showed that a certain minimum permeability is required in order to ensure the formation of a traveling thermal wave. It was found that by studying the non-dimensional governing parameters of the system one can ascribe characteristic values for these parameters for given initial and boundary conditions. This allows one to roughly predict the performance of a particular method on a particular reservoir given approximate values for initial and boundary conditions. Channelling and flow blockage due to carbon residue buildup impeded each method's performance. Blockage was found to be a result of imbalanced heating. Copyright 2012, Society of Petroleum Engineers.

Climb-enabled discrete dislocation plasticity

A small strain two-dimensional discrete dislocation plasticity framework coupled to vacancy diffusion is developed wherein the motion of edge dislocations is by a combination of glide and climb. The dislocations are modelled as line defects in a linear elastic medium and the mechanical boundary value problem is solved by the superposition of the infinite medium elastic fields of the dislocations and a complimentary non-singular solution that enforces the boundary conditions. Similarly, the climbing dislocations are modelled as line sources/sinks of vacancies and the vacancy diffusion boundary value problem is also solved by a superposition of the fields of the line sources/sinks in an infinite medium and a complementary non-singular solution that enforces the boundary conditions. The vacancy concentration field along with the stress field provides the climb rate of the dislocations. Other short-range interactions of the dislocations are incorporated via a set of constitutive rules. We first employ this formulation to investigate the climb of a single edge dislocation in an infinite medium and illustrate the existence of diffusion-limited and sink-limited climb regimes. Next, results are presented for the pure bending and uniaxial tension of single crystals oriented for single slip. These calculations show that plasticity size effects are reduced when dislocation climb is permitted. Finally, we contrast predictions of this coupled framework with an ad hoc model in which dislocation climb is modelled by a drag-type relation based on a quasi steady-state solution. © 2013 Elsevier Ltd. All rights reserved.

Non-equilibirum condensation effects on the flow field and the performance of a low pressure steam turbine

The influence of non-equilibrium condensation on the flow field and performance of a three stage low pressure model steam turbine is examined using modern three dimensional CFD techniques. An equilibrium steam model and a non-equilibrium steam model, which accounts for both subcooling and condensation effects, are used, and have been verified by comparison with test data in an earlier publication [1]. The differences in the calculated flow field and turbine performance with these models show that the latent heat released during condensation influences both the thermodynamic and the aerodynamic performance of the turbine, leading to a change in inlet flow angles of about 5°. The calculated three dimensional flowfield is used to investigate the magnitude and distribution of the additional thermo-dynamic wetness loss arising from steam condensation under non-equilibrium flow conditions. Three simple methods are described to calculate this, and all show that this amounts to around 6.5% of the total losses at the design condition. At other load conditions the wetness losses change in magnitude and axial distribution in the turbine. © 2010 by ASME.

耗散粒子动力学处理复杂固体壁面的一种有效方法

耗散粒子动力学(dissipative particle dynamics,DPD)作为一种介观尺度拉格朗日型粒子方法,已经成功地应用于微纳米流动和生化科技的研究中. 复杂固体壁面的处理和壁面边界条件的实施一直是DPD方法发展及应用的一个障碍. 提出了处理复杂固体壁面的一种新的方法. 复杂固体区域通过冻结随机分布并且达到平衡状态的DPD粒子代表;所冻结的DPD粒子位于临近流动区域的一个截距内;在靠近固体壁面的流动区域中设置流动反弹层,当流动DPD粒子进入此流动层后反弹回流动区域. 应用这种固体壁面处理方法.

Boundary Element Method (BEM) Analysis for Galvanic Corrosion of Hot Dip Galvanized Steel Immersed in Seawater

A numerical analysis of galvanic corrosion of hot-dip galvanized steel immersed in seawater was presented. The analysis was based on the boundary element methods (BEMs) coupled with Newton-Raphson iterative technique to treat the nonlinear boundary conditions, which were determined by the experimental polarization curves. Results showed that galvanic current density concentrates on the boundary of steel substrate and zinc coating, and the sacrificial protection of zinc coating to steel substrate results in overprotection of steel cathode. Not only oxygen reduction but also hydrogen reduction could occur as cathode reactions, which probably led up to the adsorption and absorption of hydrogen atoms. Flat galvanized steel tensile sample shows a brittle behavior similar to hydrogen embrittlement according to the SSRT (show strain rate test) in seawater.