Método da propagação de feixe de ângulo largo para análise de guias de ondas ópticos não-lineares

Este trabalho propõe uma extensão do método de propagação de feixe (BPM - Beam Propagation Method) para a análise de guias de ondas ópticos e acopladores baseados em materiais não-lineares do tipo Kerr. Este método se destina à investigação de estruturas onde a utilização da equação escalar de Helmholtz (EEH) em seu limite paraxial não mais se aplica. Os métodos desenvolvidos para este fim são denominados na literatura como métodos de propagação de feixe de ângulo largo. O formalismo aqui desenvolvido é baseado na técnica das diferenças finitas e nos esquemas de Crank-Nicholson (CN) e Douglas generalizado (GD). Estes esquemas apresentam como característica o fato de apresentarem um erro de truncamento em relação ao passo de discretização transversal, Δx, proporcional a O(Δx2) para o primeiro e O(Δx4). A convergência do método em ambos esquemas é otimizada pela utilização de um algoritmo interativo para a correção do campo no meio não-linear. O formalismo de ângulo largo é obtido pela expansão da EEH para os esquemas CN e GD em termos de polinômios aproximantes de Padé de ordem (1,0) e (1,1) para CN e GD, e (2,2) e (3,3) para CN. Os aproximantes de ordem superior a (1,1) apresentam sérios problemas de estabilidade. Este problema é eliminado pela rotação dos aproximantes no plano complexo. Duas condições de contorno nos extremos da janela computacional são também investigadas: 1) (TBC - Transparent Boundary Condition) e 2) condição de contorno absorvente (TAB - Transparent Absorbing Boundary). Estas condições de contorno possuem a facilidade de evitar que reflexões indesejáveis sejam transmitidas para dentro da janela computacional. Um estudo comparativo da influência destas condições de contorno na solução de guias de ondas ópticos não-lineares é também abordada neste trabalho.

Avaliação de métodos numéricos de análise linear de estabilidade para perfis de aço formados a frio.

Para o projeto de estruturas com perfis de aço formados a frio, é fundamental a compreensão dos fenômenos da instabilidade local e global, uma vez que estes apresentam alta esbeltez e baixa rigidez à torção. A determinação do carregamento crítico e a identificação do modo de instabilidade contribuem para o entendimento do comportamento dessas estruturas. Este trabalho avalia três metodologias para a análise linear de estabilidade de perfis de aço formados a frio isolados, com o objetivo de determinar os carregamentos críticos elásticos de bifurcação e os modos de instabilidade associados. Estritamente, analisa-se perfis de seção U enrijecido e Z enrijecido isolados, de diversos comprimentos e diferentes condições de vinculação e carregamento. Determinam-se os carregamentos críticos elásticos de bifurcação e os modos de instabilidade globais e locais por meio de: (i) análise com o Método das Faixas Finitas (MFF), através do uso do programa computacional CUFSM; (ii) análise com elementos finitos de barra baseados na Teoria Generalizada de Vigas (MEF-GBT), via uso do programa GBTUL; e (iii) análise com elementos finitos de casca (MEF-cascas) por meio do uso do programa ABAQUS. Algumas restrições e ressalvas com relação ao uso do MFF são apresentadas, assim como limitações da Teoria Generalizada de Viga e precauções a serem tomadas nos modelos de cascas. Analisa-se também a influência do grau de discretização da seção transversal. No entanto, não é feita avaliação em relação aos procedimentos normativos e tampouco análises não lineares, considerando as imperfeições geométricas iniciais, tensões residuais e o comportamento elastoplástico do material.

Tropically-driven climate shifts in southwestern Europe during MIS 19, a low eccentricity interglacial

The relative roles of high- versus low-latitude forcing of millennial-scale climate variability are still not well understood. Here we present terrestrial–marine climate profiles from the southwestern Iberian margin, a region particularly affected by precession, that show millennial climate oscillations related to a nonlinear response to the Earth's precession cycle during Marine Isotope Stage (MIS) 19. MIS 19 has been considered the best analogue to our present interglacial from an astronomical point of view due to the reduced eccentricity centred at 785 ka. In our records, seven millennial-scale forest contractions punctuated MIS 19 superimposed to two orbitally-driven Mediterranean forest expansions. In contrast to our present interglacial, we evidence for the first time low latitude-driven 5000-yr cycles of drying and cooling in the western Mediterranean region, along with warmth in the subtropical gyre related to the fourth harmonic of precession. These cycles indicate repeated intensification of North Atlantic meridional moisture transport that along with decrease in boreal summer insolation triggered ice growth and may have contributed to the glacial inception, at ∼774 ka. The freshwater fluxes during MIS 19ab amplified the cooling events in the North Atlantic promoting further cooling and leading to MIS 18 glaciation. The discrepancy between the dominant cyclicity observed during MIS 1, 2500-yr, and that of MIS 19, 5000-yr, challenges the similar duration of the Holocene and MIS 19c interglacials under natural boundary conditions.

Sharp Sobolev inequality involving a critical nonlinearity on a boundary

We consider the solvability of the Neumann problem for the equation -Delta u + lambda u = 0, partial derivative u/partial derivative v = Q(x)vertical bar u vertical bar(q-2)u on partial derivative Omega, where Q is a positive and continuous coefficient on partial derivative Omega, lambda is a parameter and q = 2(N - 1)/(N - 2) is a critical Sobolev exponent for the trace embedding of H-1(Omega) into L-q(partial derivative Omega). We investigate the joint effect of the mean curvature of partial derivative Omega and the shape of the graph of Q on the existence of solutions. As a by product we establish a sharp Sobolev inequality for the trace embedding. In Section 6 we establish the existence of solutions when a parameter lambda interferes with the spectrum of -Delta with the Neumann boundary conditions. We apply a min-max principle based on the topological linking.

Random vibration of the functionally graded laminates in thermal environments

This work deals with the random free vibration of functionally graded laminates with general boundary conditions and subjected to a temperature change, taking into account the randomness in a number of independent input variables such as Young's modulus, Poisson's ratio and thermal expansion coefficient of each constituent material. Based on third-order shear deformation theory, the mixed-type formulation and a semi-analytical approach are employed to derive the standard eigenvalue problem in terms of deflection, mid-plane rotations and stress function. A mean-centered first-order perturbation technique is adopted to obtain the second-order statistics of vibration frequencies. A detailed parametric study is conducted, and extensive numerical results are presented in both tabular and graphical forms for laminated plates that contain functionally graded material which is made of aluminum and zirconia, showing the effects of scattering in thermo-clastic material constants, temperature change, edge support condition, side-to-thickness ratio, and plate aspect ratio on the stochastic characteristics of natural frequencies. (c) 2005 Elsevier B.V. All rights reserved.

T-Q relation and exact solution for the XYZ chain with general non-diagonal boundary terms

We propose that the Baxter's Q-operator for the quantum XYZ spin chain with open boundary conditions is given by the j -> infinity limit of the corresponding transfer matrix with spin-j (i.e., (2j + I)-dimensional) auxiliary space. The associated T-Q relation is derived from the fusion hierarchy of the model. We use this relation to determine the Bethe Ansatz solution of the eigenvalues of the fundamental transfer matrix. The solution yields the complete spectrum of the Hamiltonian. (c) 2006 Elsevier B.V. All rights reserved.

A new relaxation method for roll forming problems

Finite element analysis (FEA) of nonlinear problems in solid mechanics is a time consuming process, but it can deal rigorously with the problems of both geometric, contact and material nonlinearity that occur in roll forming. The simulation time limits the application of nonlinear FEA to these problems in industrial practice, so that most applications of nonlinear FEA are in theoretical studies and engineering consulting or troubleshooting. Instead, quick methods based on a global assumption of the deformed shape have been used by the roll-forming industry. These approaches are of limited accuracy. This paper proposes a new form-finding method - a relaxation method to solve the nonlinear problem of predicting the deformed shape due to plastic deformation in roll forming. This method involves applying a small perturbation to each discrete node in order to update the local displacement field, while minimizing plastic work. This is iteratively applied to update the positions of all nodes. As the method assumes a local displacement field, the strain and stress components at each node are calculated explicitly. Continued perturbation of nodes leads to optimisation of the displacement field. Another important feature of this paper is a new approach to consideration of strain history. For a stable and continuous process such as rolling and roll forming, the strain history of a point is represented spatially by the states at a row of nodes leading in the direction of rolling to the current one. Therefore the increment of the strain components and the work-increment of a point can be found without moving the object forward. Using this method we can find the solution for rolling or roll forming in just one step. This method is expected to be faster than commercial finite element packages by eliminating repeated solution of large sets of simultaneous equations and the need to update boundary conditions that represent the rolls.

The boundary lubricated friction and wear of low alloy steel

Pin on disc wear machines were used to study the boundary lubricated friction and wear of AISI 52100 steel sliding partners. Boundary conditions were obtained by using speed and load combinations which resulted in friction coefficients in excess of 0.1. Lubrication was achieved using zero, 15 and 1000 ppm concentrations of an organic dimeric acid additive in a hydrocarbon base stock. Experiments were performed for sliding speeds of 0.2, 0.35 and 0.5 m/s for a range of loads up to 220 N. Wear rate, frictional force and pin temperature were continually monitored throughout tests and where possible complementary methods of measurement were used to improve accuracy. A number of analytical techniques were used to examine wear surfaces, debris and lubricants, namely: Scanning Electron Microscopy (SEM), Auger Electron Spectroscopy (AES), Powder X-ray Diffraction (XRD), X-ray Photoelectron Spectroscopy (XPS), optical microscopy, Back scattered Electron Detection (BSED) and several metallographic techniques. Friction forces and wear rates were found to vary linearly with load for any given combination of speed and additive concentration. The additive itself was found to act as a surface oxidation inhibitor and as a lubricity enhancer, particularly in the case of the higher (1000 ppm) concentration. Wear was found to be due to a mild oxidational mechanism at low additive concentrations and a more severe metallic mechanism at higher concentrations with evidence of metallic delamination in the latter case. Scuffing loads were found to increase with increasing additive concentration and decrease with increasing speed as would be predicted by classical models of additive behaviour as an organo-metallic soap film. Heat flow considerations tended to suggest that surface temperature was not the overriding controlling factor in oxidational wear and a model is proposed which suggests oxygen concentration in the lubricant is the controlling factor in oxide growth and wear.

Effect of aspect ratio on thermal convection in open and closed systems

Internally heated fluids are found across the nuclear fuel cycle. In certain situations the motion of the fluid is driven by the decay heat (i.e. corium melt pools in severe accidents, the shutdown of liquid metal reactors, molten salt and the passive control of light water reactors) as well as normal operation (i.e. intermediate waste storage and generation IV reactor designs). This can in the long-term affect reactor vessel integrity or lead to localized hot spots and accumulation of solid wastes that may prompt local increases in activity. Two approaches to the modeling of internally heated convection are presented here. These are based on numerical analysis using codes developed in-house and simulations using widely available computational fluid dynamics solvers. Open and closed fluid layers at around the transition between conduction and convection of various aspect ratios are considered. We determine optimum domain aspect ratio (1:7:7 up to 1:24:24 for open systems and 5:5:1, 1:10:10 and 1:20:20 for closed systems), mesh resolutions and turbulence models required to accurately and efficiently capture the convection structures that evolve when perturbing the conductive state of the fluid layer. Note that the open and closed fluid layers we study here are bounded by a conducting surface over an insulating surface. Conclusions will be drawn on the influence of the periodic boundary conditions on the flow patterns observed. We have also examined the stability of the nonlinear solutions that we found with the aim of identifying the bifurcation sequence of these solutions en route to turbulence.

Determining the temperature from incomplete boundary data

An iterative procedure for determining temperature fields from Cauchy data given on a part of the boundary is presented. At each iteration step, a series of mixed well-posed boundary value problems are solved for the heat operator and its adjoint. A convergence proof of this method in a weighted L2-space is included, as well as a stopping criteria for the case of noisy data. Moreover, a solvability result in a weighted Sobolev space for a parabolic initial boundary value problem of second order with mixed boundary conditions is presented. Regularity of the solution is proved. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

On the Stabilization of the Wave Equation by the Boundary

* Partially supported by CNPq (Brazil)

Asymptotic Property of Eigenvalues and Eigenfunctions of the Laplace Operator in Domain with a Perturbed Boundary

2000 Mathematics Subject Classification: 35J05, 35C15, 44P05

Anti-Periodic Boundary Value Problem for Impulsive Fractional Integro Differential Equations

MSC 2010: 34A37, 34B15, 26A33, 34C25, 34K37

Asymptotic Expansion of Solution for Almost Regular and Weakly Perturbed Systems of Ordinary Differential Equations

Л. И. Каранджулов, Н. Д. Сиракова - В работата се прилага методът на Поанкаре за решаване на почти регулярни нелинейни гранични задачи при общи гранични условия. Предполага се, че диференциалната система съдържа сингулярна функция по отношение на малкия параметър. При определени условия се доказва асимптотичност на решението на поставената задача.

Crack Theory with Singularities at the Boundary

2002 Mathematics Subject Classification: 35S15, 35J70, 35J40, 38J40