ω循环型预条件矩阵及其应用

In this paper, we apply the preconditioned conjugate gradient method to the solution of positive-definite Toeplitz systems, especially we introduce a new kind of co-circulant preconditioners Pn[ca] by the use of embedding method. We have also discussed the properties of these new preconditioners and proved that many of former preconditioners can be considered as some special cases of Pn[co\. Because of the introduction of co-circulant preconditioners pn[a>], we can greatly overcome the singularity caused by circulant preconditioners. We have discussed the oo-circulant series and functions. We compare the ordinary circularity with the co-circularity, showing that the latter one can be considered as the extended form of the former one; correspondingly, many methods and theorems of the ordinary circularity can be extended. Furthermore, we present the co-circulant decompositional method. By the use of this method, we can divide any co-circulant signal into a summation of many sub-signals; especially among those sub-signals, there are many subseries of which their period is just equal to 1, which are actually the frequency elements of the original co-circulant signal. In this way, we can establish the relationship between the signal and its frequency elements, that is, the frequency elements hi the frequency domain are actually signals with the period of 1 in the spatial domain. We have also proved that the co-circulant has already existed in the traditional Fourier theory. By the use of different criteria for constructing preconditioners, we can get many different preconditioned systems. From the preconditioned systems PN[

Toeplitz矩阵的循环延拓及其分析与应用

In this paper, we have presented the combined preconditioner which is derived from k =±-1~(1/2) circulant extensions of the real symmetric positive-definite Toeplitz matrices, proved it with great efficiency and stability and shown that it is easy to make error analysis and to remove the boundary effect with the combined preconditioner. This paper has also presented the methods for the direct and inverse computation of the real Toeplitz sets of equations and discussed many problems correspondingly, especially replaced the Toeplitz matrices with the combined preconditoners for analysis. The paper has also discussed the spectral analysis and boundary effect. Finally, as an application in geophysics, the paper makes some discussion about the squared root of a real matrix which comes from the Laplace algorithm.

Prevalent positive epistasis in Escherichia coli and Saccharomyces cerevisiae metabolic networks

Epistasis refers to the interaction between genes. Although high-throughput epistasis data from model organisms are being generated and used to construct genetic networks(1-3), the extent to which genetic epistasis reflects biologically meaningful interactions remains unclear(4-6). We have addressed this question through in silico mapping of positive and negative epistatic interactions amongst biochemical reactions within the metabolic networks of Escherichia coli and Saccharomyces cerevisiae using flux balance analysis. We found that negative epistasis occurs mainly between nonessential reactions with overlapping functions, whereas positive epistasis usually involves essential reactions, is highly abundant and, unexpectedly, often occurs between reactions without overlapping functions. We offer mechanistic explanations of these findings and experimentally validate them for 61 S. cerevisiae gene pairs.

Development of a positive-negative selection procedure for gene targeting in fish cells

The gene targeting technique is a powerful tool for analyzing functions of cloned genes and for generating transgenic animals with site-directed integration of foreign genes. In order to develop this technique in fish, positive-negative selection (PNS) and homologous recombination vectors were constructed, and their expression was examined in fish cells. A vector (pNK) for PNS consists of the neomycin resistance gene (neo) as a positive selectable marker gene and the herpes simplex virus (HSV) thymidine kinase (tk) gene as a negative selectable marker gene. Positive selection with geneticin (G418) of epithelioma papulosum of carp (EPC) cells transfected with linearized pNK vector yielded 350 colonies, while double selection of transfected EPC cells with G418 and gancyclovir (Gc) resulted in nearly complete cell death, demonstrating that the PNS procedure is effective in fish cells. Homologous recombination vectors consist of the Xiphophorus melanoma receptor kinase (X mrk(Y)) gene as homologous sequence in addition to the neo and tk genes. Conditions for homologous recombination vector transfection and drug selection were established. After verification of the feasibility of expression of homologous recombination vectors in EPC cells, the first gene targeting experiments were attempted in the Xiphophorus melanoma cell line, PSM. Positive-negative selection of the targeting vector-transfectants led to a low enrichment in this particular cell line. The reasons for the low enrichment in PSM cells were discussed. (C) 2002 Elsevier Science B.V. All rights reserved.

Forms and functions of inorganic carbon in the Jiaozhou Bay sediments

Inorganic carbon forms and their influencing factors, mutual transformation and contribution to carbon cycling in the Jiaozhou Bay sediments were discussed. The results show that inorganic carbon in sediments could be divided into five forms: NaCl form, NH3 center dot H2O form, NaOH form, NH2OH center dot HCl form and HCl form. Thereinto, NH2OH center dot HCl form and HCl form account for more than 70% of total inorganic carbon. There was close relationship among every form of inorganic carbon and their correlativity was clearly different with different sedimentary environment except the similar strong positive correlation among NH2OH center dot HCl form, HCl form and total inorganic carbon in all regions of the Jiaozhou Bay. All forms of inorganic carbon were influenced by organic carbon, pH, Eh, Es, nitrogen and phosphorus in sediments, but their influence had different characteristics in different regions. Every farm of inorganic carbon transformed into each other continuously during early diagenesis of sediments and the common phenomenon was that NaCl form, NH3 center dot H2O form, NaOH form and NH2OH center dot HCl form might transform into steady HCl form. NaCl form, NH3 center dot H2O form, NaOH form and NH2OH center dot HCl form could participate in carbon recycle and they are potential carbon source; HCl form may be buried for a long time in sediments, and it may be one of the final resting places of atmospheric CO2. Inorganic carbon which entered into sediments was about 4.98 x 10(10) g in the Jiaozhou Bay every year, in which about 1.47x10(10) g of inorganic carbon might be buried for a long time and about 3.51. x 10(10) g of inorganic carbon might return into seawater and take part in carbon recycling.

Hierarchy of structure functions for passive scalars advected by turbulent flows

A hierarchical model is proposed for the joint moments of the passive scalar dissipation and the velocity dissipation in fluid turbulence. This model predicts that the joint probability density function (PDF) of the dissipations is a bivariate log-Poisson. An analytical calculation of the scaling exponents of structure functions of the passive scalar is carried out for this hierarchical model, showing a good agreement with the results of direct numerical simulations and experiments.

Calculations of longitudinal and transverse velocity structure functions using a vortex model of isotropic turbulence

The longitudinal structure function (LSF) and the transverse structure function (TSF) in isotropic turbulence are calculated using a vortex model. The vortex model is composed of the Rankine and Burgers vortices which have the exponential distributions in the vortex Reynolds number and vortex radii. This model exhibits a power law in the inertial range and satisfies the minimal condition of isotropy that the second-order exponent of the LSF in the inertial range is equal to that of the TSF. Also observed are differences between longitudinal and transverse structure functions caused by intermittency. These differences are related to their scaling differences which have been previously observed in experiments and numerical simulations.

Scaling Functions In Conical Indentation Of Elastic-Plastic Solids

The finite element method was used to simulate the conical indentation of elastic-plastic solids with work hardening. The ratio of the initial yield strength to the Young's modulus Y/E ranged from 0 to 0.02. Based on the calculation results, two sets of scaling functions for non-dimensional hardness H/K and indenter penetration h are presented in the paper, which have closed simple mathematical form and can be used easily for engineering application. Using the present scaling functions, indentation hardness and indentation loading curves can be easily obtained for a given set of material properties. Meanwhile one can use these scaling functions to obtain material parameters by an instrumented indentation load-displacement curve for loading and unloading if Young's modulus E and Poisson's ratio nu are known.

Anisotropic, gradient and metal-like mechanical behaviors of teeth and their implications on tooth functions

The anisotropy and gradient of the elastic modulus and the hardness of teeth were investigated by means of instrumented indentation method. Such properties are attributed to the unique microstructures of teeth based on scanning electron microscopic analysis. By comparing the relationship between the ratio of hardness to the reduced elastic modulus and the ratio of elastic unloading work to the total work of teeth in course of indentation to those of other materials, we found that the material behaviors of teeth display metal-like characteristics rather than ceramics as considered traditionally. These material behaviors and relevant functions are discussed briefly.

Wave propagation and the frequency domain Green's functions in viscoelastic Biot/squirt (BISQ) media

In this paper, we examine the characteristics of elastic wave propagation in viscoelastic porous media, which contain simultaneously both the Biot-flow and the squirt-flow mechanisms (BISQ). The frequency-domain Green's functions for viscoelastic BISQ media are then derived based on the classic potential function methods. Our numerical results show that S-waves are only affected by viscoelasticity, but not by squirt-flows. However, the phase velocity and attenuation of fast P-waves are seriously influenced by both viscoelasticity and squirt-flows; and there exist two peaks in the attenuation-frequency variations of fast P-waves. In the low-frequency range, the squirt-flow characteristic length, not viscoelasticity, affects the phase velocity of slow P-waves, whereas it is opposite in the high-frequency range. As to the contribution of potential functions of two types of compressional waves to the Green's function, the squirt-flow length has a small effect, and the effects of viscoelastic parameter are mainly in the higher frequency range. Crown Copyright (C) 2006 Published by Elsevier Ltd. All rights reserved.

The positive floating potential of high-temperature particles in plasma

It has been predicted that the floating potential of particles in plasma may become positive when the particle surface temperature is high enough, but, to our knowledge, no positive floating potential has been obtained yet. In the present paper the floating potential theory of high-temperature particles in plasma is developed to cover the positive potential range for the first time, and a general approximate analytical formula for the positive floating potential with a thin plasma sheath and subsonic plasma flow is derived from the new model recently proposed by the authors. The results show that when the floating potential is positive, the net flux of charge incident on the particle approaches a constant similar to the 'electron saturation' phenomena in the case of the electric probes.

Traces of nonsmooth functions on polyhedral domains

introduced in this paper are the definitions of the traces for a class of nonsmooth functions on polyhedral domains. By analyzing their properties we get the structures of these traces.

Complex stress functions and plastic zone sizes for notch cracks subjected to various loading conditions

In this paper, the conformal mapping method is used to solve the plane problem of an infinite plate containing a central lip-shaped notch subjected to biaxial loading at a remote boundary or a surface uniform pressure on the notch. The stress intensity factors KI and KII are obtained by the derived complex stress functions. The simple analytical expressions can be applied to the situation of cracks originating from a circular or an elliptical notch. The plastic zone sizes for such notch cracks are subsequently evaluated in light of the Dugdale strip yield concept. The results are consistent with available numerical data.

Non-linear theory of a positive-column in a magnetic-field

A nonlinear theory of an intermediate pressure discharge column in a magnetic field is presented. Motion of the neutral gas is considered. The continuity and momentum transfer equations for charged particles and neutral particles are solved by numerical methods. The main result obtained is that the rotating velocities of ionic gas and neutral gas are approximately equal. Bohm's criterion and potential inversion in the presence of neutral gas motion are also discussed.

Energy Functions And Basic Equations For Ferroelectrics

Energy functions (or characteristic functions) and basic equations for ferroelectrics in use today are given by those for ordinary dielectrics in the physical and mechanical communications. Based on these basic equations and energy functions, the finite element computation of the nonlinear behavior of the ferroelectrics has been carried out by several research groups. However, it is difficult to process the finite element computation further after domain switching, and the computation results are remarkably deviating from the experimental results. For the crack problem, the iterative solution of the finite element calculation could not converge and the solutions for fields near the crack tip oscillate. In order to finish the calculation smoothly, the finite element formulation should be modified to neglect the equivalent nodal load produced by spontaneous polarization gradient. Meanwhile, certain energy functions for ferroelectrics in use today are not compatible with the constitutive equations of ferroelectrics and need to be modified. This paper proposes a set of new formulae of the energy functions for ferroelectrics. With regard to the new formulae of the energy functions, the new basic equations for ferroelectrics are derived and can reasonably explain the question in the current finite element analysis for ferroelectrics.