Optimal design of a 7 T highly homogeneous superconducting magnet for a Penning trap

A Penning trap system called Lanzhou Penning Trap (LPT) is now being developed for precise mass measurements at the Institute of Modern Physics (IMP). One of the key components is a 7 T actively shielded superconducting magnet with a clear warm bore of 156 mm. The required field homogeneity is 3 x 10(-7) over two 1 cubic centimeter volumes lying 220 mm apart along the magnet axis. We introduce a two-step method which combines linear programming and a nonlinear optimization algorithm for designing the multi-section superconducting magnet. This method is fast and flexible for handling arbitrary shaped homogeneous volumes and coils. With the help of this method an optimal design for the LPT superconducting magnet has been obtained.

基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

粘弹性介质中的地震波传播与波形反演研究

The real earth is far away from an ideal elastic ball. The movement of structures or fluid and scattering of thin-layer would inevitably affect seismic wave propagation, which is demonstrated mainly as energy nongeometrical attenuation. Today, most of theoretical researches and applications take the assumption that all media studied are fully elastic. Ignoring the viscoelastic property would, in some circumstances, lead to amplitude and phase distortion, which will indirectly affect extraction of traveltime and waveform we use in imaging and inversion. In order to investigate the response of seismic wave propagation and improve the imaging and inversion quality in complex media, we need not only consider into attenuation of the real media but also implement it by means of efficient numerical methods and imaging techniques. As for numerical modeling, most widely used methods, such as finite difference, finite element and pseudospectral algorithms, have difficulty in dealing with problem of simultaneously improving accuracy and efficiency in computation. To partially overcome this difficulty, this paper devises a matrix differentiator method and an optimal convolutional differentiator method based on staggered-grid Fourier pseudospectral differentiation, and a staggered-grid optimal Shannon singular kernel convolutional differentiator by function distribution theory, which then are used to study seismic wave propagation in viscoelastic media. Results through comparisons and accuracy analysis demonstrate that optimal convolutional differentiator methods can solve well the incompatibility between accuracy and efficiency, and are almost twice more accurate than the same-length finite difference. They can efficiently reduce dispersion and provide high-precision waveform data. On the basis of frequency-domain wavefield modeling, we discuss how to directly solve linear equations and point out that when compared to the time-domain methods, frequency-domain methods would be more convenient to handle the multi-source problem and be much easier to incorporate medium attenuation. We also prove the equivalence of the time- and frequency-domain methods by using numerical tests when assumptions with non-relaxation modulus and quality factor are made, and analyze the reason that causes waveform difference. In frequency-domain waveform inversion, experiments have been conducted with transmission, crosshole and reflection data. By using the relation between media scales and characteristic frequencies, we analyze the capacity of the frequency-domain sequential inversion method in anti-noising and dealing with non-uniqueness of nonlinear optimization. In crosshole experiments, we find the main sources of inversion error and figure out how incorrect quality factor would affect inverted results. When dealing with surface reflection data, several frequencies have been chosen with optimal frequency selection strategy, with which we use to carry out sequential and simultaneous inversions to verify how important low frequency data are to the inverted results and the functionality of simultaneous inversion in anti-noising. Finally, I come with some conclusions about the whole work I have done in this dissertation and discuss detailly the existing and would-be problems in it. I also point out the possible directions and theories we should go and deepen, which, to some extent, would provide a helpful reference to researchers who are interested in seismic wave propagation and imaging in complex media.

Design and optimization of highly nonlinear low-dispersion crystal fiber with high birefringence for four-wave mixing

We have proposed a novel type of photonic crystal fiber (PCF) with low dispersion and high nonlinearity for four-wave mixing. This type of fiber is composed of a solid silica core and a cladding with a squeezed hexagonal lattice elliptical airhole along the fiber length. Its dispersion and nonlinearity coefficient are investigated simultaneously by using the full vectorial finite element method. Numerical results show that the proposed highly nonlinear low-dispersion fiber has a total dispersion as low as +/- 2.5 ps nm(-1) km(-1) over an ultrabroad wavelength range from 1.43 to 1.8 mu m, and the corresponding nonlinearity coefficient and birefringence are about 150 W-1 km(-1) and 2.5 x 10(-3) at 1.55 mu m, respectively. The proposed PCF with low ultraflattened dispersion, high nonlinearity, and high birefringence can have important application in four-wave mixing. (C) 2010 Optical Society of America

Optimization Of Off-Null Ellipsometry For Air/Solid Interfaces

The optimization of off-null ellipsometry is described with emphasis on the improvement of sample thickness sensitivity. Optimal conditions are dependent on azimuth angle settings of the polarizer, compensator, and analyzer in a polarizer-compensator-sample-analyzer ellipsometer arrangement. Numerical simulation utilized offers an approach to present the dependence of the sensitivity on the azimuth angle settings, from which optimal settings corresponding to the best sensitivity are derived. For a series of samples of SiO2 layer (thickness in the range of 1.8-6.5 nm) on silicon substrate, the theory analysis proves that sensitivity at the optimal settings is increased 20 times compared to that at null settings used in most works, and the relationship between intensity and thickness is simplified as a linear type instead of the original nonlinear type, with the relative error reduced to similar to 1/100 at the optimal settings. Furthermore the discussion has been extended toward other factors affecting the sensitivity of the practical system, such as the linear dynamic range of the detector, the signal-to-noise ratio and the intensity from the light source, etc. Experimental results from the investigation Of SiO2 layer on silicon substrate are chosen to verify the optimization. (c) 2007 Optical Society of America.

Optimization of single attosecond x-ray pulses by genetic algorithm control of the chirp and initial phase of 5 fs laser pulses

We show that the peak intensity of single attosecond x-ray pulses is enhanced by 1 or 2 orders of magnitude, the pulse duration is greatly compressed, and the optimal propagation distance is shortened by genetic algorithm optimization of the chirp and initial phase of 5 fs laser pulses. However, as the laser intensity increases, more efficient nonadiabatic self-phase matching can lead to a dramatically enhanced harmonic yield, and the efficiency of optimization decreases in the enhancement and compression of the generated attosecond pulses. (c) 2006 Optical Society of America.

Optimization of high-order harmonic by genetic algorithm for the chirp and phase of few-cycle pulses

The brightness of a particular harmonic order is optimized for the chirp and initial phase of the laser pulse by genetic algorithm. The influences of the chirp and initial phase of the excitation pulse on the harmonic spectra are discussed in terms of the semi-classical model including the propagation effects. The results indicate that the harmonic intensity and cutoff have strong dependence on the chirp of the laser pulse, but slightly on its initial phase. The high-order harmonics can be enhanced by the optimal laser pulse and its cutoff can be tuned by optimization of the chirp and initial phase of the laser pulse.

Optimization and characterization of the guest-host polyetherketone polymer films for application in integrated optical devices

In order to optimize the loading of 3-(1, 1-dicyanothenyl)-1-phenyl-4, 5-dihydro-1H-pryazole (DCNP) in polyetherketone (PEK-c) guest-host polymer films, ten kinds of DCNP/PEK-c thin films, in which the weight per cent of DCNP changes from 5 to 50, were prepared. Their second-order nonlinear optical coefficients chi(33)((2)) at 1064 nm were measured by Using Maker fringe method after poling under the optimal poling condition. Their optical waveguide transmission losses were measured at 632.8 nm. Optimal weight per cent of the chromophore for the DCNP/PEK-c guest-host polymer system has been determined as about 20 for use in the integrated optical devices.

THEORETICAL TREATMENT OF A3B5 CHLORIDE VAPOR-PHASE EPITAXY - GROWTH, DOPING, OPTIMIZATION

A theoretical description of chloride vapour-phase epitaxy (CVPE) has been proposed which contains two-dimensional (2D) gas-dynamic equations for transport of reactive components and kinetic equations for surface growth processes connected by nonlinear adiabatic boundary conditions. No one of these stages is supposed to be the limiting one. Calculated variations of growth rate and impurity concentrations along the growing layer fit experimental data well.

Algorithm Studies on How to Obtain a Conditional Nonlinear Optimal Perturbation (CNOP)

The conditional nonlinear optimal perturbation (CNOP), which is a nonlinear generalization of the linear singular vector (LSV), is applied in important problems of atmospheric and oceanic sciences, including ENSO predictability, targeted observations, and ensemble forecast. In this study, we investigate the computational cost of obtaining the CNOP by several methods. Differences and similarities, in terms of the computational error and cost in obtaining the CNOP, are compared among the sequential quadratic programming (SQP) algorithm, the limited memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm, and the spectral projected gradients (SPG2) algorithm. A theoretical grassland ecosystem model and the classical Lorenz model are used as examples. Numerical results demonstrate that the computational error is acceptable with all three algorithms. The computational cost to obtain the CNOP is reduced by using the SQP algorithm. The experimental results also reveal that the L-BFGS algorithm is the most effective algorithm among the three optimization algorithms for obtaining the CNOP. The numerical results suggest a new approach and algorithm for obtaining the CNOP for a large-scale optimization problem.

Identification of nonlinear dynamic systems: Time-frequency filtering and skeleton curves

The nonlinear behavior varying with the instantaneous response was analyzed through the joint time-frequency analysis method for a class of S. D. O. F nonlinear system. A masking operator an definite regions is defined and two theorems are presented. Based on these, the nonlinear system is modeled with a special time-varying linear one, called the generalized skeleton linear system (GSLS). The frequency skeleton curve and the damping skeleton curve are defined to describe the main feature of the non-linearity as well. Moreover, an identification method is proposed through the skeleton curves and the time-frequency filtering technique.

FEM Simulations and Optimization about Residual Stresses in Coating Structures with Functionally Graded Materials Layer

An elasto-plastic finite element method is developed to predict the residual stresses of thermal spraying coatings with functionally graded material layer. In numerical simulations, temperature sensitivity of various material constants is included and mix

Nonlinear Analysis of Submarine Pipelines During Single Point Lifting

基于管道微单元体平衡建立了海管单点提升的非线性力学模型的控制微分方程组,使用变弧长的无量纲代换将动边界问题化为固定边界的两点边值问题,利用maple环境下编制的两点边值问题的打靶法程序得到了该问题在各个提升阶段的数值解答和在单点提升过程中管道的极限弯矩约为0.71q~{1/3}(EI)~{2/3}。

Application of the Conjugate Gradient Optimization Algorithm of Neural Network on Laser Machining

在应用激光技术加工复杂曲面时,通常以采样点集为插值点来建立曲面函数,然后实现曲面上任意坐标点的精确定位。人工神经网络的BP算法能实现函数插值,但计算精度偏低,往往达不到插值精确要求,造成较大的加工误差。提出人工神经网络的共轭梯度最优化插值新算法,并通过实例仿真,证明了这种曲面精确定位方法的可行性,从而为激光加工的三维精确定位提供了一种良好解决方案。这种方法已经应用在实际中。