A second order control-volume finite-element least-squares strategy for simulating diffusion in strongly anisotropic media

An unstructured mesh �nite volume discretisation method for simulating di�usion in anisotropic media in two-dimensional space is discussed. This technique is considered as an extension of the fully implicit hybrid control-volume �nite-element method and it retains the local continuity of the ux at the control volume faces. A least squares function recon- struction technique together with a new ux decomposition strategy is used to obtain an accurate ux approximation at the control volume face, ensuring that the overall accuracy of the spatial discretisation maintains second order. This paper highlights that the new technique coincides with the traditional shape function technique when the correction term is neglected and that it signi�cantly increases the accuracy of the previous linear scheme on coarse meshes when applied to media that exhibit very strong to extreme anisotropy ratios. It is concluded that the method can be used on both regular and irregular meshes, and appears independent of the mesh quality.

Pulse shaping properties of volume holographic gratings in anisotropic media

Based on a modified coupled wave theory, the pulse shaping properties of volume holographic gratings (VHGs) in anisotropic media VHGs are studied systematically. Taking photorefractive LiNbO3 crystals as an example, the combined effect that the grating parameters, the dispersion and optical anisotropy of the crystal, the pulse width, and the polarization state of the input ultrashort pulsed beam (UPB) have on the pulse shaping properties are considered when the input UPB with arbitrary polarization state propagates through the VHG. Under the combined effect, the diffraction bandwidth, pulse profiles of the diffracted and transmitted pulsed beams, and the total diffraction efficiency are shown. The studies indicate that the properties of the shaping of the o and e components of the input UPB in the crystal are greatly different; this difference can be used for pulse shaping applications. (c) 2006 Optical Society of America.

The resolution of the thermodynamic paradox and the theory of guided wave propagation in anisotropic media

The resolution of the so-called thermodynamic paradox is presented in this paper. It is shown, in direct contradiction to the results of several previously published papers, that the cutoff modes (evanescent modes having complex propagation constants) can carry power in a waveguide containing ferrite. The errors in all previous “proofs” which purport to show that the cutoff modes cannot carry power are uncovered. The boundary value problem underlying the paradox is studied in detail; it is shown that, although the solution is somewhat complicated, there is nothing paradoxical about it.

The general problem of electromagnetic wave propagation through rectangular guides filled inhomogeneously in cross-section with transversely magnetized ferrite is also studied. Application of the standard waveguide techniques reduces the TM part to the well-known self-adjoint Sturm Liouville eigenvalue equation. The TE part, however, leads in general to a non-self-adjoint eigenvalue equation. This equation and the associated expansion problem are studied in detail. Expansion coefficients and actual fields are determined for a particular problem.

3D raytracing through homogeneous anisotropic media with smooth interfaces

O raio conectando dois pontos em um meio anisotrópico, homogêneo por partes e com variação lateral, é calculado utilizando-se técnicas de continuação em 3D. Se combinado com algoritmos para solução do problema de valor inicial, o método pode ser estendido para o cálculo de eventos qS₁ e qS₂. O algoritmo apresenta a mesma eficiência e robustez que implementações de técnicas de continuação em meios isotrópicos. Rotinas baseadas neste algoritmo têm várias aplicações de interesse. Primeiramente, na modelagem e inversão de parâmetros elásticos na presença de anisotropia. Em segundo lugar, as iterações de Newton-Raphson produzem atributos da frente de onda como vetor vagarosidade e a matrix hessiana do tempo de trânsito, quantidades que permitem determinar o espalhamento geométrico e aproximações de segunda ordem para o tempo de trânsito. Estes atributos permitem calcular as amplitudes ao longo do raio e investigar os efeitos da anisotropia no empilhamento CRS em modelos de velocidade simples.

Numerical Simulation Of Ultrasonic Wave Propagation In Elastically Anisotropic Media

The ultrasonic non-destructive testing of components may encounter considerable difficulties to interpret some inspections results mainly in anisotropic crystalline structures. A numerical method for the simulation of elastic wave propagation in homogeneous elastically anisotropic media, based on the general finite element approach, is used to help this interpretation. The successful modeling of elastic field associated with NDE is based on the generation of a realistic pulsed ultrasonic wave, which is launched from a piezoelectric transducer into the material under inspection. The values of elastic constants are great interest information that provide the application of equations analytical models, until small and medium complexity problems through programs of numerical analysis as finite elements and/or boundary elements. The aim of this work is the comparison between the results of numerical solution of an ultrasonic wave, which is obtained from transient excitation pulse that can be specified by either force or displacement variation across the aperture of the transducer, and the results obtained from a experiment that was realized in an aluminum block in the IEN Ultrasonic Laboratory. The wave propagation can be simulated using all the characteristics of the material used in the experiment evaluation associated to boundary conditions and from these results, the comparison can be made.

Design of anisotropic reflector with birefringent thin films

A novel design for dielectric anisotropic mirrors with birefringent thin films for normal incidence is presented. This mirror consists of a stack of quarter-wave biaxial layers. The biaxial anisotropic layers can be fabricated by oblique deposition. The reflectance is different for two linear polarizations of light incidence on the mirrors. As a numerical example, the design is carried out on glass with TiO2 and ZrO2. These thin films could be applied to anisotropic reflective devices for lasers.

Modal analysis of anisotropic diffused-channel waveguide by a scalar finite element method

A procedure to model optical diffused-channel waveguides is presented in this work. The dielectric waveguides present anisotropic refractive indexes which are calculated from the proton concentration. The proton concentration inside the channel is calculated by the anisotropic 2D-linear diffusion equation and converted to the refractive indexes using mathematical relations obtained from experimental data, the arbitrary refractive index profile is modeled by a. nodal expansion in the base functions. The TE and TM-like propagation properties (effective index) and the electromagnetic fields for well-annealed proton-exchanged (APE) LiNbO3 waveguides are computed by the finite element method.

Finite element analysis of anisotropic optical waveguide with arbitrary index profile

This work presents the application of a scalar finite element formulation for Ex (TE-like) modes in anisotropic planar and channel waveguides with diagonal permittivity tensor, diffused in both transversal directions. This extended formulation considers explicitly both the variations of the refractive index and their spatial derivates inside of each finite element. Dispersion curves for Ex modes in planar and channel waveguides are shown, and the results compared with solutions obtained by other formulations.

An anisotropic linear elastic boundary element formulation for masonry wall analysis

This work presents an application of a Boundary Element Method (BEM) formulation for anisotropic body analysis using isotropic fundamental solution. The anisotropy is considered by expressing a residual elastic tensor as the difference of the anisotropic and isotropic elastic tensors. Internal variables and cell discretization of the domain are considered. Masonry is a composite material consisting of bricks (masonry units), mortar and the bond between them and it is necessary to take account of anisotropy in this type of structure. The paper presents the formulation, the elastic tensor of the anisotropic medium properties and the algebraic procedure. Two examples are shown to validate the formulation and good agreement was obtained when comparing analytical and numerical results. Two further examples in which masonry walls were simulated, are used to demonstrate that the presented formulation shows close agreement between BE numerical results and different Finite Element (FE) models. © 2012 Elsevier Ltd.

Crack extension and kinking in laminates and bicrystals

Singular fields at the tip of an interface crack in anisotropic solids are reviewed with emphasis on establishing a framework to quantify fracture resistance under mixed mode conditions. The concepts of mode mixity and surface toughness are unified by using generalized interface traction components. The similarity between the anisotropic theory and existing isotropic theory is shown. Explicit formulae are given for misoriented orthotropic bimaterials with potential applications envisioned including composite laminates and semiconductor crystals. Competition between crack extension along the interface and kinking into the substrate is investigated using a boundary layer formulation. Several case studies reveal the role of anisotropy. An explicit complex variable representation for orthotropic materials and a solution to a dislocation interacting with a crack are presented in two self-contained Appendices.

Analysis of characteristic matrix for a uniaxially birefringent thin film

As the principal section consistent with the principal plane, electromagnetic propagation in a uniaxially birefingent thin film can be described with a concise 2 x 2 characteristic matrix, in which the refractive indices of the forward and backward propagating extraordinary rays are different and dependent. In this letter, based on Huygen's construction, the refractive indices and effective optical admittances of the forward and backward propagating extraordinary rays are discussed further, and the 2 x 2 characteristic matrix is simplified. Furthermore, the input optical admittance, reflectance and transmittance of assembly is presented, just as an isotropic thin film. The result can be extended to the general case of multilayer uniaxially birefringent thin films with their optic axes in the incident plane. Copyright (c) EPLA, 2007.

Propagation properties of beams generated by Gaussian mirror resonator in uniaxial crystals

Based on the paraxial vectorial theory of beams propagating in uniaxially anisotropic media, we have derived the analytical propagation equations of beams generated by Gaussian mirror resonator (GMR) in uniaxial crystals, and given the typical numerical example to illustrate our analytical results. Due to the anisotropy crystals, the ordinary and extraordinary beams originated by incident beams generated by GMR propagate with different diffraction lengths, thus the linear polarization state and axial symmetry of the incident beams generated by GMR do not remain during propagating in crystals. (c) 2006 Elsevier GmbH. All rights reserved.

Propagation and polarization properties of hollow Gaussian beams in uniaxial crystals

Based on the paraxial vectorial theory of beams propagating in uniaxially anisotropic media, we have derived the analytical propagation equations of hollow Gaussian beams (HGBs) in uniaxial crystals, and given the typical numerical example to illustrate our analytical results. Due to the anisotropy crystals, the ordinary and extraordinary beams originated by incident HGBs propagate with different diffraction lengths, thus the linear polarization state and axial symmetry of incident HGBs do not remain during propagating in crystals. (c) 2007 Published by Elsevier B.V.

地球介质非线性与各向异性的地震波响应及其成像效应

The practice of geophysical prospecting shows us the complex interior earth. The studies of the complexity play an important role and practical guide for the subsurface structure. At present, the complexity of the earth mainly means lateral and vertical homogeneity, anisotropy and non-linear quality. And the anisotropy and non-linear media studies become the frontier in seismology and exploration seismology. This paper summarizes the development of complexities and presents the forward and inverse in the non-linear and anisotropic media. Firstly, the paper introduces the theory of seismic wave propagation in the non-linear and anisotropic media, the theoretical basis for simulation and inversion research. Secondly, high quality numerical simulation method with little dispersion has been developed to investigate the influence of complexity including anisotropy and non-linear multi-component seismograms. Because most real data in seismology have a single component, we developed two aspects work on anisotropic multi-component imaging. One is prestack reflection migration. The result show that distorted images are obtained if data from anisotropic media are migrated using isotropic extrapolation. Moreover, image quality will be improved greatly after considering anisotropy in subsurface layers. The other one is the we take advantage of multi-component data to inversion of the anisotropic parameters jointly seimic reflection travel time and polarization information. Based on these research works, we get the following results: 1.Combing numerical simulation, systematical studies indicate that anisotropy and non-linear seismograms characters are significant to detect cracked belts in the earth and to understand deformation field and mechanism. 2.Based on anisotropic media models, we developed an efficient prestack migration method for subsurface structure and different observation methods seismic data, which improving the imaging quality with VSP, seismograms and real data. 3.Jointly seismic inversion combining seismic anisotropic reflection traveltimes and polarizations data show that the complete wrong inversion and the following explanation will be resulted by ignoring anisotropy.