基于广义正交多项式迭积微分算子的地震波场数值模拟研究

Seismic Numerical Modeling is one of bases of the Exploratory Seismology and Academic Seismology, also is a research field in great demand. Essence of seismic numerical modeling is to assume that structure and parameters of the underground media model are known, simulate the wave-field and calculate the numerical seismic record that should be observed. Seismic numerical modeling is not only a means to know the seismic wave-field in complex inhomogeneous media, but also a test to the application effect by all kinds of methods. There are many seismic numerical modeling methods, each method has its own merits and drawbacks. During the forward modeling, the computation precision and the efficiency are two pivotal questions to evaluate the validity and superiority of the method. The target of my dissertation is to find a new method to possibly improve the computation precision and efficiency, and apply the new forward method to modeling the wave-field in the complex inhomogeneous media. Convolutional Forsyte polynomial differentiator (CFPD) approach developed in this dissertation is robust and efficient, it shares some of the advantages of the high precision of generalized orthogonal polynomial and the high speed of the short operator finite-difference. By adjusting the operator length and optimizing the operator coefficient, the method can involve whole and local information of the wave-field. One of main tasks of the dissertation is to develop a creative, generalized and high precision method. The author introduce convolutional Forsyte polynomial differentiator to calculate the spatial derivative of seismic wave equation, and apply the time staggered grid finite-difference which can better meet the high precision of the convolutional differentiator to substitute the conventional finite-difference to calculate the time derivative of seismic wave equation, then creating a new forward method to modeling the wave-field in complex inhomogeneous media. Comparing with Fourier pseudo-spectral method, Chebyshev pseudo-spectral method, staggered- grid finite difference method and finite element method, convolutional Forsyte polynomial differentiator (CFPD) method has many advantages: 1. Comparing with Fourier pseudo-spectral method. Fourier pseudo-spectral method (FPS) is a local operator, its results have Gibbs effects when the media parameters change, then arose great errors. Therefore, Fourier pseudo-spectral method can not deal with special complex and random heterogeneous media. But convolutional Forsyte polynomial differentiator method can cover global and local information. So for complex inhomogeneous media, CFPD is more efficient. 2. Comparing with staggered-grid high-order finite-difference method, CFPD takes less dots than FD at single wave length, and the number does not increase with the widening of the studying area. 3. Comparing with Chebyshev pseudo-spectral method (CPS). The calculation region of Chebyshev pseudo-spectral method is fixed in , under the condition of unchangeable precision, the augmentation of calculation is unacceptable. Thus Chebyshev pseudo-spectral method is inapplicable to large area. CFPD method is more applicable to large area. 4. Comparing with finite element method (FE), CFPD can use lager grids. The other task of this dissertation is to study 2.5 dimension (2.5D) seismic wave-field. The author reviews the development and present situation of 2.5D problem, expatiates the essentiality of studying the 2.5D problem, apply CFPD method to simulate the seismic wave-field in 2.5D inhomogeneous media. The results indicate that 2.5D numerical modeling is efficient to simulate one of the sections of 3D media, 2.5D calculation is much less time-consuming than 3D calculation, and the wave dispersion of 2.5D modeling is obviously less than that of 3D modeling. Question on applying time staggered-grid convolutional differentiator based on CFPD to modeling 2.5D complex inhomogeneous media was not studied by any geophysicists before, it is a fire-new creation absolutely. The theory and practices prove that the new method can efficiently model the seismic wave-field in complex media. Proposing and developing this new method can provide more choices to study the seismic wave-field modeling, seismic wave migration, seismic inversion, and seismic wave imaging.

Multistep Methods for Integrating the Solar System

High order multistep methods, run at constant stepsize, are very effective for integrating the Newtonian solar system for extended periods of time. I have studied the stability and error growth of these methods when applied to harmonic oscillators and two-body systems like the Sun-Jupiter pair. I have also tried to design better multistep integrators than the traditional Stormer and Cowell methods, and I have found a few interesting ones.

Correlated electron transport in molecular electronics

Theoretical and experimental values to date for the resistances of single molecules commonly disagree by orders of magnitude. By reformulating the transport problem using boundary conditions suitable for correlated many-electron systems, we approach electron transport across molecules from a new standpoint. Application of our correlated formalism to benzene-dithiol gives current-voltage characteristics close to experimental observations. The method can solve the open system quantum many-body problem accurately, treats spin exactly, and is valid beyond the linear response regime.

Superheavy nuclei in relativistic effective Lagrangian model

Isotopic and isotonic chains of superheavy nuclei are analyzed to search for spherical double shell closures beyond Z=82 and N=126 within the new effective field theory model of Furnstahl, Serot, and Tang for the relativistic nuclear many-body problem. We take into account several indicators to identify the occurrence of possible shell closures, such as two-nucleon separation energies, two-nucleon shell gaps, average pairing gaps, and the shell correction energy. The effective Lagrangian model predicts N=172 and Z=120 and N=258 and Z=120 as spherical doubly magic superheavy nuclei, whereas N=184 and Z=114 show some magic character depending on the parameter set. The magicity of a particular neutron (proton) number in the analyzed mass region is found to depend on the number of protons (neutrons) present in the nucleus.

La forma del afecto : el orden de los modos a cada atributo implica el conocimiento de la unión de los atributos de una sola substancia

Se trata de una totalidad que sólo es concebible cuando se han abordado todos los puntos de vista que la abarcan. Así, todo aquel que lleve a cabo una lectura comprometida de esta obra, deberá seguir una petición que hace su autor

Weak constraints in four-dimensional variational data assimilation

The formulation of four-dimensional variational data assimilation allows the incorporation of constraints into the cost function which need only be weakly satisfied. In this paper we investigate the value of imposing conservation properties as weak constraints. Using the example of the two-body problem of celestial mechanics we compare weak constraints based on conservation laws with a constraint on the background state.We show how the imposition of conservation-based weak constraints changes the nature of the gradient equation. Assimilation experiments demonstrate how this can add extra information to the assimilation process, even when the underlying numerical model is conserving.

Perspex machine VII: The universal perspex machine

The perspex machine arose from the unification of projective geometry with the Turing machine. It uses a total arithmetic, called transreal arithmetic, that contains real arithmetic and allows division by zero. Transreal arithmetic is redefined here. The new arithmetic has both a positive and a negative infinity which lie at the extremes of the number line, and a number nullity that lies off the number line. We prove that nullity, 0/0, is a number. Hence a number may have one of four signs: negative, zero, positive, or nullity. It is, therefore, impossible to encode the sign of a number in one bit, as floating-, point arithmetic attempts to do, resulting in the difficulty of having both positive and negative zeros and NaNs. Transrational arithmetic is consistent with Cantor arithmetic. In an extension to real arithmetic, the product of zero, an infinity, or nullity with its reciprocal is nullity, not unity. This avoids the usual contradictions that follow from allowing division by zero. Transreal arithmetic has a fixed algebraic structure and does not admit options as IEEE, floating-point arithmetic does. Most significantly, nullity has a simple semantics that is related to zero. Zero means "no value" and nullity means "no information." We argue that nullity is as useful to a manufactured computer as zero is to a human computer. The perspex machine is intended to offer one solution to the mind-body problem by showing how the computable aspects of mind and. perhaps, the whole of mind relates to the geometrical aspects of body and, perhaps, the whole of body. We review some of Turing's writings and show that he held the view that his machine has spatial properties. In particular, that it has the property of being a 7D lattice of compact spaces. Thus, we read Turing as believing that his machine relates computation to geometrical bodies. We simplify the perspex machine by substituting an augmented Euclidean geometry for projective geometry. This leads to a general-linear perspex-machine which is very much easier to pro-ram than the original perspex-machine. We then show how to map the whole of perspex space into a unit cube. This allows us to construct a fractal of perspex machines with the cardinality of a real-numbered line or space. This fractal is the universal perspex machine. It can solve, in unit time, the halting problem for itself and for all perspex machines instantiated in real-numbered space, including all Turing machines. We cite an experiment that has been proposed to test the physical reality of the perspex machine's model of time, but we make no claim that the physical universe works this way or that it has the cardinality of the perspex machine. We leave it that the perspex machine provides an upper bound on the computational properties of physical things, including manufactured computers and biological organisms, that have a cardinality no greater than the real-number line.

Modelling the secular evolution of migrating planet pairs

The subject of this paper is the secular behaviour of a pair of planets evolving under dissipative forces. In particular, we investigate the case when dissipative forces affect the planetary semimajor axes and the planets move inwards/outwards the central star, in a process known as planet migration. To perform this investigation, we introduce fundamental concepts of conservative and dissipative dynamics of the three-body problem. Based on these concepts, we develop a qualitative model of the secular evolution of the migrating planetary pair. Our approach is based on the analysis of the energy and the orbital angular momentum exchange between the two-planet system and an external medium; thus no specific kind of dissipative forces is invoked. We show that, under the assumption that dissipation is weak and slow, the evolutionary routes of the migrating planets are traced by the Mode I and Mode II stationary solutions of the conservative secular problem. The ultimate convergence and the evolution of the system along one of these secular modes of motion are determined uniquely by the condition that the dissipation rate is sufficiently smaller than the proper secular frequency of the system. We show that it is possible to reassemble the starting configurations and the migration history of the systems on the basis of their final states and consequently to constrain the parameters of the physical processes involved.

External Morphology of the Immatures of Polybia paulista (Hymenoptera: Vespidae)

The immatures of Polybia paulista Ihering were described using light and scanning electron microscopy and the results are compared with previous descriptions within the same or related wasps. This study is based on 2 whole nests collected in the municipality of Rio Claro, São Paulo, in Brazil. We have detected the existence of 5 larval instars. The main morphological alterations over development occur in the relative size of structures, yet certain structures appear with subsequent instars and become more evident later in development: increasing density in the number of body spines and papillae; the appearance of body setae in fifth-instar larvae; opening of spiracles upon second-instar larvae; 2 body shapes in fifth-instar larvae; the appearance of a lateral tooth on the mandibles of fourth instar; presence of spines on the maxillae of fifth-instar larvae; altered shape of galea and palps upon third-instar larvae from a cluster of sensilla to a conical elevation; and the appearance of spines on postmentum upon fourth-instar larvae. This way, the present study presents a detailed description of the immatures of P. paulista, and we hope the presented information can be useful to morphological, taxonomic, and phylogenetic studies.

Leucograma e teores plasmáticos de proteínas de fase aguda de eqüinos portadores de abdômen agudo e submetidos à laparotomia

Foram examinados 20 eqüinos adultos portadores de abdômen agudo e submetidos à laparotomia. Dez recuperaram-se sem intercorrência pós-operatória (G1) e 10 foram a óbito sete a 10 dias após a cirurgia, com sinais de choque séptico (G2). Avaliaram-se temperatura retal, freqüências cardíaca e respiratória, tempo de preenchimento capilar e teores plasmáticos das proteínas de fase aguda - fibrinogênio, ceruloplasmina, proteína C-reativa, antitripsina, haptoglobina e glicoproteína ácida -, antes e até sete dias após a laparotomia. As leucometrias às 72h e no sétimo dia pós-operatório dos eqüinos que foram a óbito foram, respectivamente, 34,6% e 57,1%, mais altas que a dos animais curados. Os maiores valores de proteína de fase aguda ocorreram no sétimo dia após a cirurgia; os percentuais de elevação de fibrinogênio, antitripsina, glicoproteina ácida, proteína C-reativa, ceruloplasmina e haptoglobina de eqüinos do G2 em relação ao G1 foram 46,8%, 67,9%, 91,9%, 112,2%, 126,9% e 186,2%, respectivamente.

Energy Expenditure, Lipid Profile, Oxidative Stress, and Cardiac Energy Metabolism After Growth Hormone Treatment in Obese Young Rats

Obesity is rampant in modern society and growth hormone (GH) could be useful as adjunct therapy to reduce the obesity-induced cardiovascular damage. To investigate GH effects on obesity, initially 32 male Wistar rats were divided into two groups (n = 16): control (C) was fed standard-chow and water and hyper-caloric (H) was fed hypercaloric chow and 30% sucrose in its drinking water. After 45 days, both C and H groups were divided into two subgroups (n = 8): C + PL was fed standard-chow, water and received saline subcutaneously; C + GH was fed standard-chow, water, and received 2 mg/kg/day GH subcutaneously; H + PL was fed hypercaloric diet, 30% sucrose in its drinking water, and received saline subcutaneously; and H + GH was fed hypercaloric diet, 30% sucrose in its drinking water, and received GH subcutaneously. After 75 days of total experimental period, H + PL rats were considered obese, having higher body weight, body mass index, Lee-index, and atherogenic index (AI) compared to C + PL. Obesity was accompanied by enhanced myocardial lipid hydroperoxide (LH) and lactate dehydrogenase (LDH), as well of depressed energy expenditure (RMR) and oxygen consumption(VO(2))/body weight. H + GH rats had higher fasting RMR, as well as lower AI and myocardial LH than H + PL. Comparing C + GH with C + PL, despite no effects on morphometric parameters, lipid profile, myocardial LH, and LDH activity, GH enhanced fed RMR and myocardial pyruvate dehydrogenase. In conclusion, the present study brought new insights into the GH effects on obesity related cardiovascular damage demonstrating, for the first time, that GH regulated cardiac metabolic pathways, enhanced energy expenditure and improved the lipid profile in obesity condition. Growth hormone in standard fed condition also offered promising therapeutic value enhancing pyruvate-dehydrogenase activity and glucose oxidation in cardiac tissue, thus optimizing myocardial energy metabolism.

Trapping neutral fermions with kink-like potentials

The intrinsically relativistic problem of neutral fermions subject to kink-like potentials (similar to tanh gamma x) is investigated and the exact bound-state solutions are found. Apart from the lonely hump solutions for E = +/- mc(2), the problem is mapped into the exactly solvable Sturm-Liouville problem with a modified Poschl-Teller potential. An apparent paradox concerning the uncertainty principle is solved by resorting to the concepts of effective mass and effective Compton wavelength. (c) 2005 Elsevier B.V. All rights reserved.

Relativistic confinement of neutral fermions with a trigonometric tangent potential

The problem of neutral fermions subject to a pseudoscalar potential is investigated. Apart from the solutions for E = +/- mc(2), the problem is mapped into the Sturm-Liouville equation. The case of a singular trigonometric tangent potential (similar to tan gamma x) is exactly solved and the complete set of solutions is discussed in some detail. It is revealed that this intrinsically relativistic and true confining potential is able to localize fermions into a region of space arbitrarily small without the menace of particle-antiparticle production.

On the evection resonance and its connection to the stability of outer satellites

In terms of stability around the primary, it is widely known that the semimajor axis of the retrograde satellites is much larger than the corresponding semimajor axis of the prograde satellites. Usually this conclusion is obtained numerically, since precise analytical derivation is far from being easy, especially, in the case of two or more disturbers. Following the seminal idea that what is unstable in the restricted three-body problem is also unstable in the general N-body problem, we present a simplified model which allows us to derive interesting resonant configurations. These configurations are responsible for cumulative perturbations which can give birth to strong instability that may cause the ejection of the satellite. Then we obtain, analytically, approximate bounds of the stability of prograde and retrograde satellites. Although we recover quite well previous results of other authors, we comment very briefly some weakness of these bounds. Copyright (c) 2008 Tadashi Yokoyama et al.

Exploring S-type orbits in the Pluto-Charon binary system

This work generates, through a sample of numerical simulations of the restricted three-body problem, diagrams of semimajor axis and eccentricity which defines stable and unstable zones for particles in S-type orbits around Pluto and Charon. Since we consider initial conditions with 0 <= e <= 0.99, we found several new stable regions. We also identified the nature of each one of these newly found stable regions. They are all associated to families of periodic orbits derived from the planar circular restricted three-body problem. We have shown that a possible eccentricity of the Pluto-Charon system slightly reduces, but does not destroy, any of the stable regions.