直径任意分布球填充的数值模拟

提出球填充数值算法的新分类方法.改进原有的松弛算法,使其能够模拟直径任意分布的球填充问题,采用可变循环周期使不同球数情形下的填充率基本保持不变.算例数据表明,该算法的填充率和配位数均高于原算法.由于采用背景网格搜索和双向链表组数据结构,使得邻接球搜索效率有相当大的提高,算法的时间复杂度为O（N）（N为球数）.在一台AMD Athlon 3200＋PC上,对于10000个等径球的随机密排列,只需217s,填充率即可达到0.64.

Dilemas institucionais e cenários políticos : análise do discurso da alteração do sobrestamento da pauta na Câmara dos Deputados

Analisa a Questão de Ordem (QO) 411/09, respondida por Michel Temer, então Presidente da Câmara dos Deputados, que alterou o entendimento a respeito do sobrestamento da pauta na Câmara. A opção pela análise de uma questão de ordem justifica-se pelo entendimento de que se trata de um dos instrumentos regimentais mais legítimos e definidores do processo legislativo, pois tem o potencial de orientar os caminhos de uma proposição e criar jurisprudência legislativa, podendo até mesmo definir os rumos da atividade parlamentar, como é o caso da QO 411/09. O estudo pretende, a partir da abordagem teórica e metodológica da vertente francesa da análise do discurso, desvelar as relações de poder no plano das formações discursivas e suas respectivas formações ideológicas no interior do Legislativo vis-à-vis o Executivo, com ênfase no poder de agenda e no equilíbrio entre os Poderes da República. Para tanto, propõe-se analisar pronunciamentos proferidos em Plenário por ocasião da supracitada Questão de Ordem, tendo como referência o discurso político como cena teatral indicado por Charaudeau, de forma a tentar fornecer novos elementos interpretativos aos dilemas institucionais e às escolhas políticas dos atores envolvidos nesse debate. Conclui-se que as estratégias discursivas empregadas pelos atores políticos evidenciaram o processo de significação, circulação e entrecruzamento dos saberes a respeito da relação institucional entre os Poderes. A análise do discurso contribuiu para a compreensão do contexto social da produção discursiva da QO e da alteração do sobrestamento da pauta, além de reiterar a ideia de que a formação discursiva está intimamente relacionada à prática institucional das relações de poder dentro e fora da Câmara dos Deputados. Os embates dos pronunciamentos evidenciam esta premissa. Além disso, a análise da QO foi capaz de revelar a força política e retórica dos pronunciamentos em Plenário no sentido de modificar e, sobretudo, legitimar o entendimento sobre o trancamento da pauta. Ao utilizar a análise do discurso, tendo como pano de fundo o processo legislativo, direito constitucional e a ciência política, pretendeu-se dar novo enfoque à temática do sobrestamento da pauta e, dessa forma, tentar fornecer novos elementos interpretativos aos dilemas institucionais e às escolhas políticas dos atores envolvidos nesse contexto.

"A casa dos ímpios se desfará, mas a tenda dos retos florescerá" : a participação da Frente Parlamentar Evangélica no legislativo brasileiro

Analisa, a partir do trabalho de campo empreendido nas atividades da Frente Parlamentar Evangélica do Congresso Nacional, como estes parlamentares participam de tempos e espaços do legislativo. Deste modo, trago a descrição etnográfica os discursos destes deputados religiosos que denotam, mormente, uma "guerra santa" entre os "ungidos" e o Inimigo que ocupa o mundo da política e o espaço público brasileiro. Assim, suspeito que presenciamos contemporaneamente uma maior dessectarização de grupos evangélicos na política quando estes agentes procuram não apenas se tornar a Igreja do amanhã (concorrendo com a Igreja Católica pela posição de religião civil). Mas, sobretudo, tornar, através da participação na política, o mundo dos ímpios o mundo dos justos. Ora, considero, por isso mesmo, as formas particulares de aquisição de sutilezas e de competências necessárias a fim de profissionalizar o político e as relações políticas nas quais estes políticos se engajam. Deste modo, os evangélicos também têm que aprender tais habilidades especialmente por que estão no processo de construção de uma herança histórica na política. Por fim, analiso como os evangélicos, enquanto políticos profissionais, agenciam uma forma de representação e participação políticas que considera tanto o exercício cotidiano da fé quanto uma espécie de investimento nas atividades do século, reinterpretando, assim, o histórico processo de secularização do Estado.

Participação popular na produção legislativa : o caso da legislação urbanística no Brasil pós-Constituinte

Na década de 1980 a questão urbana no Brasil ganhou visibilidade a partir de uma gama de ocupações de imóveis ociosos que ocorreram em diversas cidades importantes do país. Às mobilizações pela redemocratização somou-se a luta pela reforma urbana. Com o advento do Processo Constituinte, ganhou importância a intervenção dos movimentos populares na esfera institucional, cuja principal ferramenta foi a apresentação da Emenda Popular da Reforma Urbana, convertida numa espécie de manifesto-programa pelos defensores da causa. Com a Constituição já promulgada, o movimento pela reforma urbana celebrou o fato de, pela primeira vez na História Constitucional brasileira, a questão urbana ter sido contemplada. O Capítulo de Política Urbana necessitou, contudo, de regulamentação para ter efeitos práticos, o que veio a acontecer com a aprovação do Estatuto da Cidade em 2001. Fazer um balanço da luta pela reforma urbana no Brasil e a da influência do movimento popular no delineamento da legislação urbanística desde o Processo Constituinte é o objetivo primeiro deste trabalho. A intenção de fundo é refletir acerca do modelo de democracia brasileira, tendo por pressuposto a ideia de que a participação popular modelou um regime democrático que avança em relação ao clássico modelo da democracia representativa vigente no mundo ocidental.

Slope Stability Analysis Using Anisotropic and Nonlinear Failure Criteria

The shear strength of soils or rocks developed in a landslide usually exhibits anisotropic and nonlinear behavior. The process of sedimentation and subsequent consolidation can cause anisotropy of sedimentary soils or rocks, for instance. Nonlinearity of failure envelope could be attributed to "interlocking" or "dilatancy" of the material, which is generally dependent upon the stress level. An analytical method considering both anisotropy and nonlinearity of the failure envelops of soil and rocks is presented in the paper. The nonlinearfailure envelopes can be determined from routine triaxial tests. A spreadsheet program, which uses the Janbu's Generalized Procedure of Slice and incorporates anisotropic, illustrates the implementation of the approach and nonlinearfailure envelops. In the analysis, an equivalent Mohr-Coulomb linear failure criterion is obtained by drawing a tangent to the nonlinear envelope of an anisotropic soil at an appropriate stress level. An illustrative example is presented to show the feasibility and numerical efficiency of the method.

Microstructure Evolution of Laser Pulse processed Ductile Iron Surface

Pulsed laser beam was used to modify surface processing for ductile iron. The microstructures of processed specimen were observed using optical microscope (OM). Nanoindentation and micro-hardness of microstructures were measured from surface to inner of sample. The experimental results show that, modification zone is consisted of light melted zone, phase transformation hardening area and transient area. The light melt area is made up of coarse dendrite crystalline with a thickness less than 20um, phase transformation hardening area mainly of laminal or acicular martensite, retained austenite and graphite, i.e. M+A prime+ G. The cow-eye microstructure around graphite sphere always is formed in phase transformation hardening area zone, which consisting of a variety structure with the distance from the surface. So, it maybe as a obvious sign distinguishing modification zone border. Finally, the microstructures evolution of laser pulse processed ductile iron was analyzed coupling with beam energy distribution in space and laser pulse heating procession characteristics. The analysis shows that energy distribution of laser pulse has an important effect on microstructure during laser pulse modified ductile iron. Multi-scale and interlace arrangement are the important features for laser pulse modified ductile iron. Of microstructure.

Mechanics of adhesive contact on a power-law graded elastic half-space

We consider adhesive contact between a rigid sphere of radius R and a graded elastic half-space with Young's modulus varying with depth according to a power law E = E-0(z/c(0))(k) (0 < k < 1) while Poisson's ratio v remaining a constant. Closed-form analytical solutions are established for the critical force, the critical radius of contact area and the critical interfacial stress at pull-off. We highlight that the pull-off force has a simple solution of P-cr= -(k+3)pi R Delta gamma/2 where Delta gamma is the work of adhesion and make further discussions with respect to three interesting limits: the classical JKR solution when k = 0, the Gibson solid when k --> 1 and v = 0.5, and the strength limit in which the interfacial stress reaches the theoretical strength of adhesion at pull-off. (C) 2009 Elsevier Ltd. All rights reserved.

An extension of the two dimensional JKR theory to the case with a large contact width

In the Hertz and JKR theories, parabolic assumptions for the rounded profiles of the sphere or cylinder are adopted under the condition that the contact radius (width) should be very small compared to the radius of the sphere or cylinder. However, a large contact radius (width) is often found in experiments even under a zero external loading. We aim at extending the plane strain JKR theory to the case with a large contact width. The relation between the external loading and the contact width is given. Solutions for the Hertz, JKR and rounded-profile cases are compared and analyzed. It is found that when the ratio of a/R is approximately larger than about 0.4, the parabolic assumptions in the Hertz and JKR theories are no longer valid and the exact rounded profile function should be used.

Numerical study on heavy rigid particle motion in a plane wake flow by spectral element method

Experimental particle dispersion patterns in a plane wake flow at a high Reynolds number have been predicted numerically by discrete vortex method (Phys. Fluids A 1992; 4:2244-2251; Int. J. Multiphase Flow 2000; 26:1583-1607). To address the particle motion at a moderate Reynolds number, spectral element method is employed to provide an instantaneous wake flow field for particle dynamics equations, which are solved to make a detail classification of the patterns in relation to the Stokes and Froude numbers. It is found that particle motion features only depend on the Stokes number at a high Froude number and depend on both numbers at a low Froude number. A ratio of the Stokes number to squared Froude number is introduced and threshold values of this parameter are evaluated that delineate the different regions of particle behavior. The parameter describes approximately the gravitational settling velocity divided by the characteristic velocity of wake flow. In order to present effects of particle density but preserve rigid sphere, hollow sphere particle dynamics in the plane wake flow is investigated. The evolution of hollow particle motion patterns for the increase of equivalent particle density corresponds to that of solid particle motion patterns for the decrease of particle size. Although the thresholds change a little, the parameter can still make a good qualitative classification of particle motion patterns as the inner diameter changes.

Movement of a Spherical Brownian Particle Between Infinite Parallel Plates: Hindered Dispersion and Sedimentation

The dispersion of an isolated, spherical, Brownian particle immersed in a Newtonian fluid between infinite parallel plates is investigated. Expressions are developed for both a 'molecular' contribution to dispersion, which arises from random thermal fluctuations, and a 'convective' contribution, arising when a shear flow is applied between the plates. These expressions are evaluated numerically for all sizes of the particle relative to the bounding plates, and the method of matched asymptotic expansions is used to develop analytical expressions for the dispersion coefficients as a function of particle size to plate spacing ratio for small values of this parameter.

It is shown that both the molecular and convective dispersion coefficients decrease as the size of the particle relative to the bounding plates increase. When the particle is small compared to the plate spacing, the coefficients decrease roughly proportional to the particle size to plate spacing ratio. When the particle closely fills the space between the plates, the molecular dispersion coefficient approaches zero slowly as an inverse logarithmic function of the particle size to plate spacing ratio, and the convective dispersion coefficent approaches zero approximately proportional to the width of the gap between the edges of the sphere and the bounding plates.

Invariants, Lie-Bäcklund operators and Bäcklund transformations

This thesis is mainly concerned with the application of groups of transformations to differential equations and in particular with the connection between the group structure of a given equation and the existence of exact solutions and conservation laws. In this respect the Lie-Bäcklund groups of tangent transformations, particular cases of which are the Lie tangent and the Lie point groups, are extensively used.

In Chapter I we first review the classical results of Lie, Bäcklund and Bianchi as well as the more recent ones due mainly to Ovsjannikov. We then concentrate on the Lie-Bäcklund groups (or more precisely on the corresponding Lie-Bäcklund operators), as introduced by Ibragimov and Anderson, and prove some lemmas about them which are useful for the following chapters. Finally we introduce the concept of a conditionally admissible operator (as opposed to an admissible one) and show how this can be used to generate exact solutions.

In Chapter II we establish the group nature of all separable solutions and conserved quantities in classical mechanics by analyzing the group structure of the Hamilton-Jacobi equation. It is shown that consideration of only Lie point groups is insufficient. For this purpose a special type of Lie-Bäcklund groups, those equivalent to Lie tangent groups, is used. It is also shown how these generalized groups induce Lie point groups on Hamilton's equations. The generalization of the above results to any first order equation, where the dependent variable does not appear explicitly, is obvious. In the second part of this chapter we investigate admissible operators (or equivalently constants of motion) of the Hamilton-Jacobi equation with polynornial dependence on the momenta. The form of the most general constant of motion linear, quadratic and cubic in the momenta is explicitly found. Emphasis is given to the quadratic case, where the particular case of a fixed (say zero) energy state is also considered; it is shown that in the latter case additional symmetries may appear. Finally, some potentials of physical interest admitting higher symmetries are considered. These include potentials due to two centers and limiting cases thereof. The most general two-center potential admitting a quadratic constant of motion is obtained, as well as the corresponding invariant. Also some new cubic invariants are found.

In Chapter III we first establish the group nature of all separable solutions of any linear, homogeneous equation. We then concentrate on the Schrodinger equation and look for an algorithm which generates a quantum invariant from a classical one. The problem of an isomorphism between functions in classical observables and quantum observables is studied concretely and constructively. For functions at most quadratic in the momenta an isomorphism is possible which agrees with Weyl' s transform and which takes invariants into invariants. It is not possible to extend the isomorphism indefinitely. The requirement that an invariant goes into an invariant may necessitate variants of Weyl' s transform. This is illustrated for the case of cubic invariants. Finally, the case of a specific value of energy is considered; in this case Weyl's transform does not yield an isomorphism even for the quadratic case. However, for this case a correspondence mapping a classical invariant to a quantum orie is explicitly found.

Chapters IV and V are concerned with the general group structure of evolution equations. In Chapter IV we establish a one to one correspondence between admissible Lie-Bäcklund operators of evolution equations (derivable from a variational principle) and conservation laws of these equations. This correspondence takes the form of a simple algorithm.

In Chapter V we first establish the group nature of all Bäcklund transformations (BT) by proving that any solution generated by a BT is invariant under the action of some conditionally admissible operator. We then use an algorithm based on invariance criteria to rederive many known BT and to derive some new ones. Finally, we propose a generalization of BT which, among other advantages, clarifies the connection between the wave-train solution and a BT in the sense that, a BT may be thought of as a variation of parameters of some. special case of the wave-train solution (usually the solitary wave one). Some open problems are indicated.

Most of the material of Chapters II and III is contained in [I], [II], [III] and [IV] and the first part of Chapter V in [V].

Two new integral transforms and their applications

This thesis is in two parts. In Part I the independent variable θ in the trigonometric form of Legendre's equation is extended to the range ( -∞, ∞). The associated spectral representation is an infinite integral transform whose kernel is the analytic continuation of the associated Legendre function of the second kind into the complex θ-plane. This new transform is applied to the problems of waves on a spherical shell, heat flow on a spherical shell, and the gravitational potential of a sphere. In each case the resulting alternative representation of the solution is more suited to direct physical interpretation than the standard forms.

In Part II separation of variables is applied to the initial-value problem of the propagation of acoustic waves in an underwater sound channel. The Epstein symmetric profile is taken to describe the variation of sound with depth. The spectral representation associated with the separated depth equation is found to contain an integral and a series. A point source is assumed to be located in the channel. The nature of the disturbance at a point in the vicinity of the channel far removed from the source is investigated.

On the settling speed of dilute arrays of spheres

A method is developed to calculate the settling speed of dilute arrays of spheres for the three cases of: I, a random array of freely moving particles; II, a random array of rigidly held particles; and III, a cubic array of particles. The basic idea of the technique is to give a formal representation for the solution and then manipulate this representation in a straightforward manner to obtain the result. For infinite arrays of spheres, our results agree with the results previously found by other authors, and the analysis here appears to be simpler. This method is able to obtain more terms in the answer than was possible by Saffman's unified treatment for point particles. Some results for arbitrary two sphere distributions are presented, and an analysis of the wall effect for particles settling in a tube is given. It is expected that the method presented here can be generalized to solve other types of problems.

Compressible flows at small Reynolds numbers

The problem of the slow viscous flow of a gas past a sphere is considered. The fluid cannot be treated incompressible in the limit when the Reynolds number Re, and the Mach number M, tend to zero in such a way that Re ~ o(M^2 ). In this case, the lowest order approximation to the steady Navier-Stokes equations of motion leads to a paradox discovered by Lagerstrom and Chester. This paradox is resolved within the framework of continuum mechanics using the classical slip condition and an iteration scheme that takes into account certain terms in the full Navier-Stokes equations that drop out in the approximation used by the above authors. It is found however that the drag predicted by the theory does not agree with R. A. Millikan's classic experiments on sphere drag.

The whole question of the applicability of the Navier-Stokes theory when the Knudsen number M/Re is not small is examined. A new slip condition is proposed. The idea that the Navier-Stokes equations coupled with this condition may adequately describe small Reynolds number flows when the Knudsen number is not too large is looked at in some detail. First, a general discussion of asymptotic solutions of the equations for all such flows is given. The theory is then applied to several concrete problems of fluid motion. The deductions from this theory appear to interpret and summarize the results of Millikan over a much wider range of Knudsen numbers (almost up to the free molecular or kinetic limit) than hitherto Believed possible by a purely continuum theory. Further experimental tests are suggested and certain interesting applications to the theory of dilute suspensions in gases are noted. Some of the questions raised in the main body of the work are explored further in the appendices.

Response of multiple simultaneously recorded macaque area LIP neurons in a memory saccade task

Cells in the lateral intraparietal cortex (LIP) of rhesus macaques respond vigorously and in spatially-tuned fashion to briefly memorized visual stimuli. Responses to stimulus presentation, memory maintenance, and task completion are seen, in varying combination from neuron to neuron. To help elucidate this functional segmentation a new system for simultaneous recording from multiple neighboring neurons was developed. The two parts of this dissertation discuss the technical achievements and scientific discoveries, respectively.

Technology. Simultanous recordings from multiple neighboring neurons were made with four-wire bundle electrodes, or tetrodes, which were adapted to the awake behaving primate preparation. Signals from these electrodes were partitionable into a background process with a 1/f-like spectrum and foreground spiking activity spanning 300-6000 Hz. Continuous voltage recordings were sorted into spike trains using a state-of-the-art clustering algorithm, producing a mean of 3 cells per site. The algorithm classified 96% of spikes correctly when tetrode recordings were confirmed with simultaneous intracellular signals. Recording locations were verified with a new technique that creates electrolytic lesions visible in magnetic resonance imaging, eliminating the need for histological processing. In anticipation of future multi-tetrode work, the chronic chamber microdrive, a device for long-term tetrode delivery, was developed.

Science. Simultaneously recorded neighboring LIP neurons were found to have similar preferred targets in the memory saccade paradigm, but dissimilar peristimulus time histograms, PSTH). A majority of neighboring cell pairs had a difference in preferred directions of under 45° while the trial time of maximal response showed a broader distribution, suggesting homogeneity of tuning with het erogeneity of function. A continuum of response characteristics was present, rather than a set of specific response types; however, a mapping experiment suggests this may be because a given cell's PSTH changes shape as well as amplitude through the response field. Spike train autocovariance was tuned over target and changed through trial epoch, suggesting different mechanisms during memory versus background periods. Mean frequency-domain spike-to-spike coherence was concentrated below 50 Hz with a significant maximum of 0.08; mean time-domain coherence had a narrow peak in the range ±10 ms with a significant maximum of 0.03. Time-domain coherence was found to be untuned for short lags (10 ms), but significantly tuned at larger lags (50 ms).