User’s Guide to For2R: A Module of Fortran 95 Output Routines Compatible with the R Statistics Language

For2R is a collection of Fortran routines for saving complex data structures into a file that can be read in the R statistics environment with a single command.1 For2R provides both the means to transfer data structures significantly more complex than simple tables, and an archive mechanism to store data for future reference. We developed this software because we write and run computationally intensive numerical models in Fortran, C++, and AD Model Builder. We then analyse results with R. We desired to automate data transfer to speed diagnostics during working-group meetings. We thus developed the For2R interface to write an R data object (of type list) to a plain-text file. The master list can contain any number of matrices, values, dataframes, vectors or lists, all of which can be read into R with a single call to the dget function. This allows easy transfer of structured data from compiled models to R. Having the capacity to transfer model data, metadata, and results has sharply reduced the time spent on diagnostics, and at the same time, our diagnostic capabilities have improved tremendously. The simplicity of this interface and the capabilities of R have enabled us to automate graph and table creation for formal reports. Finally, the persistent storage in files makes it easier to treat model results in analyses or meta-analyses devised months—or even years—later. We offer For2R to others in the hope that they will find it useful. (PDF contains 31 pages)

Numerical investigation of crack growth in concrete subjected to compression by the generalized beam lattice model

The beam lattice-type models, such as the Euler-Bernoulli (or Timoshenko) beam lattice and the generalized beam (GB) lattice, have been proved very effective in simulating failure processes in concrete and rock due to its simplicity and easy implementation. However, these existing lattice models only take into account tensile failures, so it may be not applicable to simulation of failure behaviors under compressive states. The main aim in this paper is to incorporate Mohr-Coulomb failure criterion, which is widely used in many kinds of materials, into the GB lattice procedure. The improved GB lattice procedure has the capability of modeling both element failures and contact/separation of cracked elements. The numerical examples show its effectiveness in simulating compressive failures. Furthermore, the influences of lateral confinement, friction angle, stiffness of loading platen, inclusion of aggregates on failure processes are respectively analyzed in detail.

The elasto-damage theory of the components assembling model

The potential energy in materials is well approximated by pair functional which is composed of pair potentials and embedding energy. During calculating material potential energy, the orientational component and the volumetric component are derived respectively from pair potentials and embedding energy. The sum of energy of all these two kinds of components is the material potential. No matter how microstructures change, damage or fracture, at the most level, they are all the changing and breaking atomic bonds. As an abstract of atomic bonds, these components change their stiffness during damaging. Material constitutive equations have been formulated by means of assembling all components' response functions. This material model is called the component assembling model. Theoretical analysis and numerical computing indicate that the proposed model has the capacity of reproducing some results satisfactorily, with the advantages of great conceptual simplicity, physical explicitness, and intrinsic induced anisotropy, etc.

Seismic reflection experiments imaging the physical nature of crustal structures in southern California

Seismic reflection methods have been extensively used to probe the Earth's crust and suggest the nature of its formative processes. The analysis of multi-offset seismic reflection data extends the technique from a reconnaissance method to a powerful scientific tool that can be applied to test specific hypotheses. The treatment of reflections at multiple offsets becomes tractable if the assumptions of high-frequency rays are valid for the problem being considered. Their validity can be tested by applying the methods of analysis to full wave synthetics.

Three studies illustrate the application of these principles to investigations of the nature of the crust in southern California. A survey shot by the COCORP consortium in 1977 across the San Andreas fault near Parkfield revealed events in the record sections whose arrival time decreased with offset. The reflectors generating these events are imaged using a multi-offset three-dimensional Kirchhoff migration. Migrations of full wave acoustic synthetics having the same limitations in geometric coverage as the field survey demonstrate the utility of this back projection process for imaging. The migrated depth sections show the locations of the major physical boundaries of the San Andreas fault zone. The zone is bounded on the southwest by a near-vertical fault juxtaposing a Tertiary sedimentary section against uplifted crystalline rocks of the fault zone block. On the northeast, the fault zone is bounded by a fault dipping into the San Andreas, which includes slices of serpentinized ultramafics, intersecting it at 3 km depth. These interpretations can be made despite complications introduced by lateral heterogeneities.

In 1985 the Calcrust consortium designed a survey in the eastern Mojave desert to image structures in both the shallow and the deep crust. Preliminary field experiments showed that the major geophysical acquisition problem to be solved was the poor penetration of seismic energy through a low-velocity surface layer. Its effects could be mitigated through special acquisition and processing techniques. Data obtained from industry showed that quality data could be obtained from areas having a deeper, older sedimentary cover, causing a re-definition of the geologic objectives. Long offset stationary arrays were designed to provide reversed, wider angle coverage of the deep crust over parts of the survey. The preliminary field tests and constant monitoring of data quality and parameter adjustment allowed 108 km of excellent crustal data to be obtained.

This dataset, along with two others from the central and western Mojave, was used to constrain rock properties and the physical condition of the crust. The multi-offset analysis proceeded in two steps. First, an increase in reflection peak frequency with offset is indicative of a thinly layered reflector. The thickness and velocity contrast of the layering can be calculated from the spectral dispersion, to discriminate between structures resulting from broad scale or local effects. Second, the amplitude effects at different offsets of P-P scattering from weak elastic heterogeneities indicate whether the signs of the changes in density, rigidity, and Lame's parameter at the reflector agree or are opposed. The effects of reflection generation and propagation in a heterogeneous, anisotropic crust were contained by the design of the experiment and the simplicity of the observed amplitude and frequency trends. Multi-offset spectra and amplitude trend stacks of the three Mojave Desert datasets suggest that the most reflective structures in the middle crust are strong Poisson's ratio (σ) contrasts. Porous zones or the juxtaposition of units of mutually distant origin are indicated. Heterogeneities in σ increase towards the top of a basal crustal zone at ~22 km depth. The transition to the basal zone and to the mantle include increases in σ. The Moho itself includes ~400 m layering having a velocity higher than that of the uppermost mantle. The Moho maintains the same configuration across the Mojave despite 5 km of crustal thinning near the Colorado River. This indicates that Miocene extension there either thinned just the basal zone, or that the basal zone developed regionally after the extensional event.

Angular distributions of multiphoton detachment of H- in various infrared laser fields

Using a nonperturbative quantum scattering theory, the photoelectron angular distributions (PADs) from the multiphoton detachment of H- ions in strong, linearly polarized infrared laser fields are obtained to interpret recent experimental observations. In our theoretical treatment, the PADs in n-photon detachment are determined by the nth-order generalized phased Bessel (GPB) functions X-n(Z(f),eta). The advantage of using the GPB scenario to calculate PADs is its simplicity: a single special function (GPB) without any mixing coefficient can express PADs observed by recent experiments. Thus, the GPB scenario can be called a parameterless scenario.

Modeling and analysis of hysteretic structural behavior

For damaging response, the force-displacement relationship of a structure is highly nonlinear and history-dependent. For satisfactory analysis of such behavior, it is important to be able to characterize and to model the phenomenon of hysteresis accurately. A number of models have been proposed for response studies of hysteretic structures, some of which are examined in detail in this thesis. There are two popular classes of models used in the analysis of curvilinear hysteretic systems. The first is of the distributed element or assemblage type, which models the physical behavior of the system by using well-known building blocks. The second class of models is of the differential equation type, which is based on the introduction of an extra variable to describe the history dependence of the system.

Owing to their mathematical simplicity, the latter models have been used extensively for various applications in structural dynamics, most notably in the estimation of the response statistics of hysteretic systems subjected to stochastic excitation. But the fundamental characteristics of these models are still not clearly understood. A response analysis of systems using both the Distributed Element model and the differential equation model when subjected to a variety of quasi-static and dynamic loading conditions leads to the following conclusion: Caution must be exercised when employing the models belonging to the second class in structural response studies as they can produce misleading results.

The Massing's hypothesis, originally proposed for steady-state loading, can be extended to general transient loading as well, leading to considerable simplification in the analysis of the Distributed Element models. A simple, nonparametric identification technique is also outlined, by means of which an optimal model representation involving one additional state variable is determined for hysteretic systems.

Neural network design for switching network control

A neural network is a highly interconnected set of simple processors. The many connections allow information to travel rapidly through the network, and due to their simplicity, many processors in one network are feasible. Together these properties imply that we can build efficient massively parallel machines using neural networks. The primary problem is how do we specify the interconnections in a neural network. The various approaches developed so far such as outer product, learning algorithm, or energy function suffer from the following deficiencies: long training/ specification times; not guaranteed to work on all inputs; requires full connectivity.

Alternatively we discuss methods of using the topology and constraints of the problems themselves to design the topology and connections of the neural solution. We define several useful circuits-generalizations of the Winner-Take-All circuitthat allows us to incorporate constraints using feedback in a controlled manner. These circuits are proven to be stable, and to only converge on valid states. We use the Hopfield electronic model since this is close to an actual implementation. We also discuss methods for incorporating these circuits into larger systems, neural and nonneural. By exploiting regularities in our definition, we can construct efficient networks. To demonstrate the methods, we look to three problems from communications. We first discuss two applications to problems from circuit switching; finding routes in large multistage switches, and the call rearrangement problem. These show both, how we can use many neurons to build massively parallel machines, and how the Winner-Take-All circuits can simplify our designs.

Next we develop a solution to the contention arbitration problem of high-speed packet switches. We define a useful class of switching networks and then design a neural network to solve the contention arbitration problem for this class. Various aspects of the neural network/switch system are analyzed to measure the queueing performance of this method. Using the basic design, a feasible architecture for a large (1024-input) ATM packet switch is presented. Using the massive parallelism of neural networks, we can consider algorithms that were previously computationally unattainable. These now viable algorithms lead us to new perspectives on switch design.

Estudio y soluciones a la problemática de contaminación acústica en el entorno del aeropuerto de Loiu

[ES]El rápido aumento de las nuevas tecnologías, así como el aumento del uso tanto del transporte público como del privado ha conllevado un aumento en los niveles de ruido. Además, la proximidad de las industrias a las ciudades y los trabajos diarios a pie de calle no han mejorado esta situación, sino que la ha empeorado. Debido a este aumento de la intensidad sonora, las enfermedades que derivan de ello han aumentado. Por ello, las soluciones para disminuir los niveles de ruido que llegan a la población se han ido desarrollando cada vez más. En el caso de las viviendas cercanas a los aeropuertos, la solución adoptada comúnmente es la insonorización acústica del edificio por su relativa sencillez. En este trabajo, se estudia un caso concreto de un edificio residencial cercano al aeropuerto de Loiu, Bilbao y se resuelve la opción más sencilla posible de aislamiento acústico para la vivienda.

Injection Locked Clocking and Transmitter Equalization Techniques for Chip to Chip Interconnects

Semiconductor technology scaling has enabled drastic growth in the computational capacity of integrated circuits (ICs). This constant growth drives an increasing demand for high bandwidth communication between ICs. Electrical channel bandwidth has not been able to keep up with this demand, making I/O link design more challenging. Interconnects which employ optical channels have negligible frequency dependent loss and provide a potential solution to this I/O bandwidth problem. Apart from the type of channel, efficient high-speed communication also relies on generation and distribution of multi-phase, high-speed, and high-quality clock signals. In the multi-gigahertz frequency range, conventional clocking techniques have encountered several design challenges in terms of power consumption, skew and jitter. Injection-locking is a promising technique to address these design challenges for gigahertz clocking. However, its small locking range has been a major contributor in preventing its ubiquitous acceptance.

In the first part of this dissertation we describe a wideband injection locking scheme in an LC oscillator. Phase locked loop (PLL) and injection locking elements are combined symbiotically to achieve wide locking range while retaining the simplicity of the latter. This method does not require a phase frequency detector or a loop filter to achieve phase lock. A mathematical analysis of the system is presented and the expression for new locking range is derived. A locking range of 13.4 GHz–17.2 GHz (25%) and an average jitter tracking bandwidth of up to 400 MHz are measured in a high-Q LC oscillator. This architecture is used to generate quadrature phases from a single clock without any frequency division. It also provides high frequency jitter filtering while retaining the low frequency correlated jitter essential for forwarded clock receivers.

To improve the locking range of an injection locked ring oscillator; QLL (Quadrature locked loop) is introduced. The inherent dynamics of injection locked quadrature ring oscillator are used to improve its locking range from 5% (7-7.4GHz) to 90% (4-11GHz). The QLL is used to generate accurate clock phases for a four channel optical receiver using a forwarded clock at quarter-rate. The QLL drives an injection locked oscillator (ILO) at each channel without any repeaters for local quadrature clock generation. Each local ILO has deskew capability for phase alignment. The optical-receiver uses the inherent frequency to voltage conversion provided by the QLL to dynamically body bias its devices. A wide locking range of the QLL helps to achieve a reliable data-rate of 16-32Gb/s and adaptive body biasing aids in maintaining an ultra-low power consumption of 153pJ/bit.

From the optical receiver we move on to discussing a non-linear equalization technique for a vertical-cavity surface-emitting laser (VCSEL) based optical transmitter, to enable low-power, high-speed optical transmission. A non-linear time domain optical model of the VCSEL is built and evaluated for accuracy. The modelling shows that, while conventional FIR-based pre-emphasis works well for LTI electrical channels, it is not optimum for the non-linear optical frequency response of the VCSEL. Based on the simulations of the model an optimum equalization methodology is derived. The equalization technique is used to achieve a data-rate of 20Gb/s with power efficiency of 0.77pJ/bit.

Discovery of viruses lysing blue-green algae in the Dneprovsk reservoirs. [Translation of: Water blooms (ed. A. V. Topachevskii) pp.171-174. Kiev, Naukova Dumka, 1968]

Viruses, which are characterised by a relative simplicity of chemical composition, are involved with all the groups of the animal and plant world. The discovery of viruses of lower organisms has special interest. Along with the already known viruses lysing bacteria and actinomycetes, viruses have been discovered in recent years which lyse algae. During investigations of water from water-bloom patches and of mud taken from zones of massive accumulation of blue-green algae in the Dneprovsk reservoirs, the authors obtained viruses lysing algae. The revealing of viruses producing lysis of blue-green algae, which one could use in the control of water-blooms, has the greatest interest. With this aim, samples of water were collected from various zones of water-bloom patches in the Kremenchug, Dneprovsk and Kukhov reservoirs. For viruses lysing algae we propose the name 'algophages'. Along with the existence of viruses of algae of the phage type, one cannot deny the possibility of the existence of viruses of another type, multiplying in the cells of algae and causing their virus illnesses.

Resistência à flexão e microestrutura de vidro-cerâmicas feldspáticas processadas por injeção em moldes à alta temperatura

Enquanto os métodos tradicionais de processamento de restaurações cerâmicas tornaram-se notórios por sua complexidade, as técnicas mais modernas vêm privilegiando a simplicidade de execução e a automação. Dentre estas, destaca-se a injeção em moldes, que recentemente, foi associada a métodos CAD-CAM. Estudos anteriores demonstraram a viabilidade de utilização de um vidro feldspático de baixa expansão térmica, Alpha (Vita Zahnfabrik), para injeção, porém, faltam informações quanto às propriedades mecânicas e a microestrutura deste material quando submetido à injeção. Os objetivos deste estudo são: produzir pastilhas vidrocerâmicas para injeção a quente a partir de Alpha e da mistura deste vidro com partículas de alumina e zircônia; avaliar a resistência à flexão dos materiais processados, e compará-la a um material compatível existente no mercado (PM9 - Vita Zahnfabrik); estudar a estrutura microscópica dos materiais e correlacioná-la com suas propriedades mecânicas; identificar por meio de difração de raios X a formação de fases cristalinas durante as diferentes etapas de processamento. A injeção aumentou a resistência do vidro Alpha devido à redução na quantidade e tamanho dos defeitos internos, principalmente porosidades. Apesar de ter sido observada nucleação de cristais nos dois materiais, durante o processamento, não foi possível determinar de que forma este fenômeno afetou as propriedades mecânicas dos materiais. Não foi detectada alteração no padrão de distribuição das fases cristalinas observadas em microscópio eletrônico de varredura antes e depois da injeção. Não foi verificada diferença estatística significante entre a resistência à flexão de Alpha injetado e PM9. A adição de partículas de alumina e zirconia ao vidro Alpha provocou redução da resistência, devido à formação de aglomerados durante a confecção das pastilhas e a incapacidade da injeção em dispersá-los. Tais aglomerados funcionaram como concentradores de tensões, enfraquecendo o material.

On a method of determination of chlorophyll in cells of algae collected on a membrane filter. [Translation from: Trudy Instituta Biologii Vnutrennykh Vodnany 11(14), 198-202, 1966.]

The original method, proposed by Yentsch (1957), of determination of chlorophyll directly in the cells, attracts attention by its simplicity. In order to measure the content of chlorophyll by this method, a determined volume of suspension of algae is filtered through a membrane filter. The latter is dried a little, clarified by immersion oil, clamped between two glasses, and spectrophotometrized. Extinction is read off at , wavelengths equal to 670 millimicrons (around the maximum absorption of chlorophyll a in the cell) and 750 millimicrons (correction for non- specific absorption and dispersion of light by particles of the preparation). The method of Yentsch was employed by the authors for determination of chlorophyll-a in samples of phytoplankton. They conclude that in spite of the simplicity and convenience of determination the method must be applied sufficiently carefully. It is more suitable for analysis of cultures of algae, where, non-specific absorption of light is insignificant.

Full and Model-Reduced Structure-Preserving Simulation of Incompressible Fluids

This thesis outlines the construction of several types of structured integrators for incompressible fluids. We first present a vorticity integrator, which is the Hamiltonian counterpart of the existing Lagrangian-based fluid integrator. We next present a model-reduced variational Eulerian integrator for incompressible fluids, which combines the efficiency gains of dimension reduction, the qualitative robustness to coarse spatial and temporal resolutions of geometric integrators, and the simplicity of homogenized boundary conditions on regular grids to deal with arbitrarily-shaped domains with sub-grid accuracy.

Both these numerical methods involve approximating the Lie group of volume-preserving diffeomorphisms by a finite-dimensional Lie-group and then restricting the resulting variational principle by means of a non-holonomic constraint. Advantages and limitations of this discretization method will be outlined. It will be seen that these derivation techniques are unable to yield symplectic integrators, but that energy conservation is easily obtained, as is a discretized version of Kelvin's circulation theorem.

Finally, we outline the basis of a spectral discrete exterior calculus, which may be a useful element in producing structured numerical methods for fluids in the future.

Epidermis differentiation in Salmo irideus Gilb. development [Translation from: Archo ital.Anat.Embriol. 76, 19-30, 1971]

Recent histochemical and histophysiological research on the skin of fish has posed interesting problems both with regard to the presence of specialized types of cell and with regard to the function of this organ. The present paper intends to study the development and the functional differentiation of the trout's skin, both from the histomorphological and the histochemical point of view. The skin of Salmonids is devoid of specialized cells; the granular cells of Petromyzonti, the serous cells of Selacii and the clavate cells described in many other teleosts and it lacks keratin. As such it can be considered a good working model which can be used to show the eventual histomorphological and histochemical changes occurring both in the transformation from the endovular to the aqueous environment. Histomorphological observations were carried out demonstrate the great structural simplicity of the trout's skin at all stages. The article concludes that glycogen increases greatly when the epidermis thickens and therefore when it becomes necessary to guarantee resistance, amongst which is a support mechanism.

Wave-front analysis method of circular aperture sampling for collimation testing

A new type of wave-front analysis method for the collimation testing of laser beams is proposed. A concept of wave-front height is defined, and, on this basis, the wave-front analysis method of circular aperture sampling is introduced. The wave-front height of the tested noncollimated wave can be estimated from the distance between two identical fiducial diffraction planes of the sampled wave, and then the divergence is determined. The design is detailed, and the experiment is demonstrated. The principle and experiment results of the method are presented. Owing to the simplicity of the method and its low cost, it is a promising method for checking the collimation of a laser beam with a large divergence. © 2005 Optical Society of America.