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.

基于信息重建技术的地震信号处理方法研究

The dissertation addressed the problems of signals reconstruction and data restoration in seismic data processing, which takes the representation methods of signal as the main clue, and take the seismic information reconstruction (signals separation and trace interpolation) as the core. On the natural bases signal representation, I present the ICA fundamentals, algorithms and its original applications to nature earth quake signals separation and survey seismic signals separation. On determinative bases signal representation, the paper proposed seismic dada reconstruction least square inversion regularization methods, sparseness constraints, pre-conditioned conjugate gradient methods, and their applications to seismic de-convolution, Radon transformation, et. al. The core contents are about de-alias uneven seismic data reconstruction algorithm and its application to seismic interpolation. Although the dissertation discussed two cases of signal representation, they can be integrated into one frame, because they both deal with the signals or information restoration, the former reconstructing original signals from mixed signals, the later reconstructing whole data from sparse or irregular data. The goal of them is same to provide pre-processing methods and post-processing method for seismic pre-stack depth migration. ICA can separate the original signals from mixed signals by them, or abstract the basic structure from analyzed data. I surveyed the fundamental, algorithms and applications of ICA. Compared with KL transformation, I proposed the independent components transformation concept (ICT). On basis of the ne-entropy measurement of independence, I implemented the FastICA and improved it by covariance matrix. By analyzing the characteristics of the seismic signals, I introduced ICA into seismic signal processing firstly in Geophysical community, and implemented the noise separation from seismic signal. Synthetic and real data examples show the usability of ICA to seismic signal processing and initial effects are achieved. The application of ICA to separation quake conversion wave from multiple in sedimentary area is made, which demonstrates good effects, so more reasonable interpretation of underground un-continuity is got. The results show the perspective of application of ICA to Geophysical signal processing. By virtue of the relationship between ICA and Blind Deconvolution , I surveyed the seismic blind deconvolution, and discussed the perspective of applying ICA to seismic blind deconvolution with two possible solutions. The relationship of PC A, ICA and wavelet transform is claimed. It is proved that reconstruction of wavelet prototype functions is Lie group representation. By the way, over-sampled wavelet transform is proposed to enhance the seismic data resolution, which is validated by numerical examples. The key of pre-stack depth migration is the regularization of pre-stack seismic data. As a main procedure, seismic interpolation and missing data reconstruction are necessary. Firstly, I review the seismic imaging methods in order to argue the critical effect of regularization. By review of the seismic interpolation algorithms, I acclaim that de-alias uneven data reconstruction is still a challenge. The fundamental of seismic reconstruction is discussed firstly. Then sparseness constraint on least square inversion and preconditioned conjugate gradient solver are studied and implemented. Choosing constraint item with Cauchy distribution, I programmed PCG algorithm and implement sparse seismic deconvolution, high resolution Radon Transformation by PCG, which is prepared for seismic data reconstruction. About seismic interpolation, dealias even data interpolation and uneven data reconstruction are very good respectively, however they can not be combined each other. In this paper, a novel Fourier transform based method and a algorithm have been proposed, which could reconstruct both uneven and alias seismic data. I formulated band-limited data reconstruction as minimum norm least squares inversion problem where an adaptive DFT-weighted norm regularization term is used. The inverse problem is solved by pre-conditional conjugate gradient method, which makes the solutions stable and convergent quickly. Based on the assumption that seismic data are consisted of finite linear events, from sampling theorem, alias events can be attenuated via LS weight predicted linearly from low frequency. Three application issues are discussed on even gap trace interpolation, uneven gap filling, high frequency trace reconstruction from low frequency data trace constrained by few high frequency traces. Both synthetic and real data numerical examples show the proposed method is valid, efficient and applicable. The research is valuable to seismic data regularization and cross well seismic. To meet 3D shot profile depth migration request for data, schemes must be taken to make the data even and fitting the velocity dataset. The methods of this paper are used to interpolate and extrapolate the shot gathers instead of simply embedding zero traces. So, the aperture of migration is enlarged and the migration effect is improved. The results show the effectiveness and the practicability.

Spectroscopic properties of Ce : YAP scintillation crystal grown by temperature gradient technique

Spectroscopic properties of Ce-doped yttrium orthoaluminate (Ce:YAlO3 or Ce:YAP) crystals grown by temperature gradient technique (TGT) were investigated, and the effects of the growth conditions on the properties were analyzed.. Methods of optical absorption (OA), photoluminescence (PL), photoluminescence decay (PLD), X-ray excited luminescence (XL) and cathodeluminescence (CL) were used in these investigations. The results showed that the absorption band peak at 202, 394 and 532 nm originated from F and F+ color center induced by the weak reducing growth atmosphere, green emission band near 500 ran derived from Ce3+ -Ce3+ pairs and band at 650 similar to 850 run from some unintentional impurity in crystals.

INVESTIGATION OF C-60 SINGLE-CRYSTAL BY X-RAY-METHODS

C-60 Single crystals grown by a single-temperature-gradient technique were characterized by synchrotron radiation white beam x-ray topography and x-ray double crystal diffraction with Cu K-alpha 1 radiation on conventional x-ray source. The results show that the crystal is rather well crystallized, The x-ray topographies give an evidence of dendritic growth mechanism of C-60 Single crystal, and x-ray double crystal diffraction rocking curve shows that there are mosaic structural defects in the sample. A phase transition st 249+/-1.5% K from a simple cubic to a face centered cubic structure is confirmed by in situ observation of synchrotron radiation white beam x-ray topography with the temperature varing from 230 to 295 K.

Soil microbial community analysis using denaturing gradient gel electrophoresis

The most biological diversity on this planet is probably harbored in soils. Understanding the diversity and function of the microbiological component of soil poses great challenges that are being overcome by the application of molecular biological approaches. This review covers one of many approaches being used: separation of polymerase chain reaction (PCR) amplicons using denaturing gradient gel electrophoresis (DGGE). Extraction of nucleic acids directly from soils allows the examination of a community without the limitation posed by cultivation. Polymerase chain reaction provides a means to increase the numbers of a target for its detection on gels. Using the rRNA genes as a target for PCR provides phylogenetic information on populations comprising communities. Fingerprints produced by this method have allowed spatial and temporal comparisons of soil communities within and between locations or among treatments. Numerous samples can be compared because of the rapid high throughput nature of this method. Scientists now have the means to begin addressing complex ecological questions about the spatial, temporal, and nutritional interactions faced by microbes in the soil environment.

Enrichment of the hydrogen-producing microbial community from marine intertidal sludge by different pretreatment methods

To determine the effects of pretreatment on hydrogen production and the hydrogen-producing microbial community, we treated the sludge from the intertidal zone of a bathing beach in Tianjin with four different pretreatment methods, including acid treatment, heat-shock, base treatment as well as freezing and thawing. The results showed that acid pretreatment significantly promoted the hydrogen production by sludge and provided the highest efficiency of hydrogen production among the four methods. The efficiency of the hydrogen production of the acid-pretreated sludge was 0.86 +/- 0.07 mol H-2/mol glucose (mean +/- S.E.), whereas that of the sludge treated with heat-shock, freezing and thawing, base method and control was 0.41 +/- 0.03 mol H-2/mol glucose, 0.17 +/- 0.01 mol H-2/mol glucose, 0.11 +/- 0.01 mol H-2/mol glucose and 0.20 +/- 0.04 mol H-2/mol glucose, respectively. The result of denaturing gradient gel electrophoresis (DGGE) showed that pretreatment methods altered the composition of the microbial community that accounts for hydrogen production. Acid and heat pretreatments were favorable to enrich the dominant hydrogen-producing bacterium, i.e. Clostridium sp., Enterococcus sp. and Bacillus sp., However, besides hydrogen-producing bacteria, much non-hydrogen-producing Lactobacillus sp. was also found in the sludge pretreated with base, freezing and thawing methods. Therefore, based on our results, we concluded that, among the four pretreatment methods using acid, heat-shock, base or freezing and thawing, acid pretreatment was the most effective method for promoting hydrogen production of microbial community. (C) 2009 Professor T. Nejat Veziroglu. Published by Elsevier Ltd. All rights reserved.

Evaluation of the Effectiveness of Representative Methods for Determining Young's Modulus and Hardness from Instrumented Indentation Data

The effectiveness of Oliver & Pharr's (O&P's) method, Cheng & Cheng's (C&C's) method, and a new method developed by our group for estimating Young's modulus and hardness based on instrumented indentation was evaluated for the case of yield stress to reduced Young's modulus ratio (sigma(y)/E-r) >= 4.55 x 10(-4) and hardening coefficient (n) <= 0.45. Dimensional theorem and finite element simulations were applied to produce reference results for this purpose. Both O&P's and C&C's methods overestimated the Young's modulus under some conditions, whereas the error can be controlled within +/- 16% if the formulation was modified with appropriate correction functions. Similar modification was not introduced to our method for determining Young's modulus, while the maximum error of results was around +/- 13%. The errors of hardness values obtained from all the three methods could be even larger and were irreducible with any correction scheme. It is therefore suggested that when hardness values of different materials are concerned, relative comparison of the data obtained from a single standard measurement technique would be more practically useful. It is noted that the ranges of error derived from the analysis could be different if different ranges of material parameters sigma(y)/E-r and n are considered.

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

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

Strain gradient theory with couple stress for crystalline solids

A new phenomenological strain gradient theory for crystalline solid is proposed. It fits within the framework of general couple stress theory and involves a single material length scale Ics. In the present theory three rotational degrees of freedom omega (i) are introduced, which denote part of the material angular displacement theta (i) and are induced accompanying the plastic deformation. omega (i) has no direct dependence upon u(i) while theta = (1 /2) curl u. The strain energy density omega is assumed to consist of two parts: one is a function of the strain tensor epsilon (ij) and the curvature tensor chi (ij), where chi (ij) = omega (i,j); the other is a function of the relative rotation tensor alpha (ij). alpha (ij) = e(ijk) (omega (k) - theta (k)) plays the role of elastic rotation reason The anti-symmetric part of Cauchy stress tau (ij) is only the function of alpha (ij) and alpha (ij) has no effect on the symmetric part of Cauchy stress sigma (ij) and the couple stress m(ij). A minimum potential principle is developed for the strain gradient deformation theory. In the limit of vanishing l(cs), it reduces to the conventional counterparts: J(2) deformation theory. Equilibrium equations, constitutive relations and boundary conditions are given in detail. For simplicity, the elastic relation between the anti-symmetric part of Cauchy stress tau (ij), and alpha (ij) is established and only one elastic constant exists between the two tensors. Combining the same hardening law as that used in previously by other groups, the present theory is used to investigate two typical examples, i.e., thin metallic wire torsion and ultra-thin metallic beam bend, the analytical results agree well with the experiment results. While considering the, stretching gradient, a new hardening law is presented and used to analyze the two typical problems. The flow theory version of the present theory is also given.

Effect of Microtopography, Slope Length and Gradient, and Vegetative Cover on Overland Flow Through Simulation

Overland flow on a hillslope is significantly influenced by its microtopography, slope length and gradient, and vegetative cover. A 1D kinematic wave model in conjunction with a revised form of the Green-Ampt infiltration equation was employed to evaluate the effect of these surface conditions. The effect of these conditions was treated through the resistance parameter in the kinematic wave model. The resistance in this paper was considered to be made up of grain resistance, form resistance, and wave resistance. It was found that irregular slopes with microtopography eroded more easily than did regular slopes. The effect of the slope gradient on flow velocity and flow shear stress could be negative or positive. With increasing slope gradient, the flow velocity and shear stress first increased to a peak value, then decreased again, suggesting that there exists a critical slope gradient for flow velocity and shear stress. The vegetative cover was found to protect soil from erosion primarily by enhancing erosion-resisting capacity rather than by decreasing the eroding capability of overland flow.

The Flow Theory of Mechanism-Based Strain Gradient Plasticity

The flow theory of mechanism-based strain gradient (MSG) plasticity is established in this paper following the same multiscale, hierarchical framework for the deformation theory of MSG plasticity in order to connect with the Taylor model in dislocation mechanics. We have used the flow theory of MSG plasticity to study micro-indentation hardness experiments. The difference between deformation and flow theories is vanishingly small, and both agree well with experimental hardness data. We have also used the flow theory of MSG plasticity to investigate stress fields around a stationary mode-I crack tip as well as around a steady state, quasi-statically growing crack tip. At a distance to crack tip much larger than dislocation spacings such that continuum plasticity still applies, the stress level around a stationary crack tip in MSG plasticity is significantly higher than that in classical plasticity. The same conclusion is also established for a steady state, quasi-statically growing crack tip, though only the flow theory can be used because of unloading during crack propagation. This significant stress increase due to strain gradient effect provides a means to explain the experimentally observed cleavage fracture in ductile materials [J. Mater. Res. 9 (1994) 1734, Scripta Metall. Mater. 31 (1994) 1037; Interface Sci. 3(1996) 169].

Stress Wave Propagation in a Gradient Elastic Medium

The gradient elastic constitutive equation incorporating the second gradient of the strains is used to determine the monochromatic elastic plane wave propagation in a gradient infinite medium and thin rod. The equation of motion, together with the internal material length, has been derived. Various dispersion relations have been determined. We present explicit expressions for the relationship between various wave speeds, wavenumber and internal material length.

A new hardening law for strain gradient plasticity

A new hardening law of the strain gradient theory is proposed in this paper, which retains the essential structure of the incremental version of conventional J(2) deformation theory and obeys thermodynamic restrictions. The key feature of the new proposal is that the term of strain gradient plasticity is represented as an internal variable to increase the tangent modulus. This feature which is in contrast to several proposed theories, allows the problem of incremental equilibrium equations to be stated without higher-order stress, higher-order strain rates or extra boundary conditions. The general idea is presented and compared with the theory given by Fleck and Hutchinson (Adv. in Appl. Mech. (1997) 295). The new hardening law is demonstrated by two experimental tests i.e. thin wire torsion and ultra-thin beam bending tests. The present theoretical results agree well with the experiment results.