Perturbed steepest-descent technique in multiextremal problems

The steepest-descent technique dealing with the perturbed values of the objective function and its gradients and with nonexact line searches is considered. Attention is given to the case where the perturbations do not decrease on the algorithm trajectories; the aim is to investigate how perturbations at every iteration of the algorithm perturb its original attractor set.Based on the Liapunov direct method for attraction analysis of discrete-time processes, a sharp estimation of the attractor set generated by a perturbed steepest-descent technique with respect to the perturbation magnitudes is obtained. Some global optimization properties of finite-difference analogues of the gradient method are discovered. These properties are not inherent in methods which use exact gradients.The author is grateful to the referees for many useful suggestions.     

Propagating beam analysis of fiber wavelength filters with a depressed-index section

A wavelength filter consisting of single-mode and few-mode fibers is investigated numerically. A simple finite-difference beam-propagation method, in which a transparent boundary condition can be imposed, is developed for circularly symmetric waveguides. After confirming the validity of the numerical method by the mode-mismatch loss, we calculate the propagating field in the fiber wavelength filter, in which interference between LP₀₁, and LP₀₂ modes occurs. To improve the filtering operation, a depressed-index fiber is employed for the few-mode fiber. The effects of the radius and refractive index of the depressed section on the transmission power are revealed and discussed. Power is suppressed to less than 0.1% at 1.3 μm, while maintaining power transmission of more than 85% at 1.55 μm. It is also found that the filtering operation shifts to higher wavelengths as the input power is increased when we choose a self-focusing nonlinear material in the depressed section.     

Turbulence investigations during ergodic divertor operation with induced 2/1 tearing mode

Within the last year of TEXTOR operation a major part of the experiments were dedicated to the interaction of the ergodized plasma edge, induced by the dynamic ergodic divertor (DED), with the core plasma. Systematic experiments are performed in the 3/1 configuration of the DED. Depending on , the DED generates a locked 2/1 mode in the plasma at a critical current in the DED coils. Different diagnostics were used to measure the plasma response in the core and edge. The investigations in this paper are focused on the measurement of density fluctuation and the analysis of the coherent-mode structure outside the q = 2 surface using O-mode poloidal correlation reflectometry. A reduction of the quasi-coherent mode amplitude and frequency is observed, caused by the 2/1 mode together with a decrease of the density scale length. Moreover, the deduced poloidal rotation velocity at the q = 3 surface changes from the electron to the ion diamagnetic drift direction with the onset of the 2/1 island. Also a reduction of the density fluctuations is observed as well as a decrease in the poloidal correlation length. From Ar-injection a delayed inward propagation of the emission maximum is observed for different Ar ionization stages when the 2/1 mode is generated.Similar experiments with a slowly rotating DED show that the steepened density gradients are spatially localized and confirm the hypothesis that the 2/1 mode is responsible for the observations.Presented at the Workshop Electric Fields Structures and Relaxation in Edge Plasmas, Nice, France, October 26–27, 2004.     

Non-destructive testing of tubes using a time reverse numerical simulation (TRNS) method

A method for the detection of defects in cylindrical structures and the determination of their positions and orientations is presented in this paper. The scattered field, which is generated by the interaction of excited guided waves with a defect, is evaluated with an approach named time reverse numerical simulation method (TRNS). Since the excited waves and the scattered field propagate along the sample, the time-consuming scanning of the whole tube can be eliminated. The scattered displacement field is measured in three dimensions over time with a laser vibrometer at different locations distributed equally around the circumference at a fixed axial coordinate far away from the defect. Instead of analyzing the complicated time signals directly, they are played back in time. If the recorded displacement histories of the scattered field are reversed in time and played back in an identical structure, the waves travel back the same path and interfere to a maximum at their origin. The result is an amplitude increase at the position of the defect where the scattered field was generated. Instead of playing back the recorded time signals in an experiment, this step is replaced by a numerical simulation. Only this enables the visualization and detection of the amplitude increase. As long as the simulation is of high accuracy, the position of the maximum interference corresponds exactly to the location of the defect in the experiment, although no defect is implemented in the simulation.     

Optical Characterization of Iridescent Wings of <Emphasis Type="Italic">Morpho</Emphasis> Butterflies using a High Accuracy Nonstandard Finite-Difference Time-Domain Algorithm

In certain species of moths and butterflies iridescent colors arise from subwavelength diffractive surface corrugation of the wing-scales. The optical properties of such structures depend strongly on wavelength, incidence angle, and state of polarization of illuminating radiation, and the viewing angle. In this paper, we study the reflection spectra of the wings of the Morpho didius butterfly by simulating a double-layered model of a transverse cross-section comprised of the ground scale and the cover scale. Each layer contains a certain quasi-periodic arrangement of tree-like subwavelength microstructures. The simulation is done using a high accuracy nonstandard finite-difference time-domain (FDTD) method in two dimensions. We assume that the structure is made of a slightly lossy dielectric material. The wavelength dependence of the complex refractive index for the ground scale of Morpho didius is assumed to be similar to that of Morpho sulkowskyi. The complex refractive index in the latter case was obtained by comparing the computed reflection/transmission spectra with corresponding experimental measurements at normal incidence.     

Numerical Analysis of Electromagnetic Fields in Multiscale Model

Modeling technique for electromagnetic fields excited by antennas is an important topic in computational electromagnetics, which is concerned with the numerical solution of Maxwell's equations. In this paper, a novel hybrid technique that combines method of moments(MoM) with finite-difference time-domain(FDTD) method is presented to handle the problem. This approach employed Huygen's principle to realize the hybridization of the two classical numerical algorithms. For wideband electromagnetic data, the interpolation scheme is used in the MoM based on the dyadic Green's function. On the other hand, with the help of equivalence principle, the scattered electric and magnetic fields on the Huygen's surface calculated by MoM are taken as the sources for FDTD. Therefore, the electromagnetic fields in the environment can be obtained by employing finite-difference time-domain method. Finally, numerical results show the validity of the proposed technique by analyzing two canonical samples.     

Analytical perfect fluid perturbation for the matter dominated epoch

We present the analytical solution to the linear evolution equation of a one component Friedmann perturbation using an equation of state of the form p = (1/3)μσ²(t), where μ is the mass density and σ(t) is the root mean square (rms) velocity in the matter dominated epoch. It is assumed that this rms velocity depends only on the time coordinate and decreases as 1/a, a being the expansion factor of the Friedmann background. The evolution equations are written for scales below the horizon using the longitudinal gauge. The general solution, in the coordinate space, of the evolution equation for the scalar mode is obtained. In the case of spherical symmetry, this solution is expressed in terms of unidimensional integrals of the initial conditions: the initial values of the Newtonian potential and its first time derivative. This perfect fluid solution is a good approximation to the evolution of warm dark matter perturbations obtained by solving the Vlasov’s equation for collisionless particles.     

横向划痕对光场调制的近场模拟

熔石英亚表面划痕对入射激光的近场调制是导致光学元件低阈值损伤的主要因素之一. 用三维时域有限差分方法研究了连续横向划痕的近场分布, 对比了尖锐截面与光滑截面场调制的差异, 着重探讨了光场调制与划痕宽深比R的关系. 研究表明: 酸蚀后的光滑截面有助于减弱近场调制, 这类划痕的R>10.0时调制较弱且相互接近, R<5.0时调制显著增强. 当R取1---3时, 亚表面的调制达最大值, 最大电场幅值为入射波幅值的4.3倍. 当R取1.0---3.5时, 缺陷附近有80%以上取样点的最大电场幅值超过入射波幅值的2倍. 随着深度的增大, 强场区具有明显的"趋肤效应": 位于划痕正下方的强场区首先往左右两侧移动, 然后移向抛物口界面以及水平界面, 同时衍生出的多条增强线诱导整个亚表面层的光场增强.     

Control of the Electronic Ground State on an Electron‐Transfer Copper Site by Second‐Sphere Perturbations

The Cu_A center is a dinuclear copper site that serves as an optimized hub for long‐range electron transfer in heme–copper terminal oxidases. Its electronic structure can be described in terms of a σ_u* ground‐state wavefunction with an alternative, less populated ground state of π_u symmetry, which is thermally accessible. It is now shown that second‐sphere mutations in the Cu_A containing subunit of Thermus thermophilus ba₃ oxidase perturb the electronic structure, which leads to a substantial increase in the population of the π_u state, as shown by different spectroscopic methods. This perturbation does not affect the redox potential of the metal site, and despite an increase in the reorganization energy, it is not detrimental to the electron‐transfer kinetics. The mutations were achieved by replacing the loops that are involved in protein–protein interactions with cytochrome c, suggesting that transient protein binding could also elicit ground‐state switching in the oxidase, which enables alternative electron‐transfer pathways.