Temperature dependence of conformational properties of short polyelectrolytes from simulations at a single temperature

We use a series expansion method introduced recently by Rickman and Phillpot (Phys. Rev. Lett. 1991, 66, 349) to study the temperature dependent conformational properties of short ionized polyelectrolyte chains in ionic solutions by conducting simulations at a single temperature. The charged beads located at the sites of a cubic lattice interact through screened Coulombic interactions. It is shown that this method provides results that correlate with other Monte Carlo simulations, performed over a range of temperatures, where conformational transitions induced by thermal and screening effects occur. It is also shown that the method can be used successfully when the potential is weakly dependent on temperature. © 1994 American Chemical Society.

Finite-element discretization of 3D energy-transport equations for semiconductors

In this thesis a mathematical model was derived that describes the charge and energy transport in semiconductor devices like transistors. Moreover, numerical simulations of these physical processes are performed. In order to accomplish this, methods of theoretical physics, functional analysis, numerical mathematics and computer programming are applied. After an introduction to the status quo of semiconductor device simulation methods and a brief review of historical facts up to now, the attention is shifted to the construction of a model, which serves as the basis of the subsequent derivations in the thesis. Thereby the starting point is an important equation of the theory of dilute gases. From this equation the model equations are derived and specified by means of a series expansion method. This is done in a multi-stage derivation process, which is mainly taken from a scientific paper and which does not constitute the focus of this thesis. In the following phase we specify the mathematical setting and make precise the model assumptions. Thereby we make use of methods of functional analysis. Since the equations we deal with are coupled, we are concerned with a nonstandard problem. In contrary, the theory of scalar elliptic equations is established meanwhile. Subsequently, we are preoccupied with the numerical discretization of the equations. A special finite-element method is used for the discretization. This special approach has to be done in order to make the numerical results appropriate for practical application. By a series of transformations from the discrete model we derive a system of algebraic equations that are eligible for numerical evaluation. Using self-made computer programs we solve the equations to get approximate solutions. These programs are based on new and specialized iteration procedures that are developed and thoroughly tested within the frame of this research work. Due to their importance and their novel status, they are explained and demonstrated in detail. We compare these new iterations with a standard method that is complemented by a feature to fit in the current context. A further innovation is the computation of solutions in three-dimensional domains, which are still rare. Special attention is paid to applicability of the 3D simulation tools. The programs are designed to have justifiable working complexity. The simulation results of some models of contemporary semiconductor devices are shown and detailed comments on the results are given. Eventually, we make a prospect on future development and enhancements of the models and of the algorithms that we used.

Baryons in the chiral regime

Quantum Chromodynamics (QCD) is the theory of strong interactions, one of the four fundamental forces in our Universe. It describes the interaction of gluons and quarks which build up hadrons like protons and neutrons. Most of the visible matter in our universe is made of protons and neutrons. Hence, we are interested in their fundamental properties like their masses, their distribution of charge and their shape. \\rnThe only known theoretical, non-perturbative and {\it ab initio} method to investigate hadron properties at low energies is lattice Quantum Chromodynamics (lattice QCD). However, up-to-date simulations (especially for baryonic quantities) do not achieve the accuracy of experiments. In fact, current simulations do not even reproduce the experimental values for the form factors. The question arises wether these deviations can be explained by systematic effects in lattice QCD simulations.rnrnThis thesis is about the computation of nucleon form factors and other hadronic quantities from lattice QCD. So called Wilson fermions are used and the u- and d-quarks are treated fully dynamically. The simulations were performed using gauge ensembles with a range of lattice spacings, volumes and pion masses.\\rnFirst of all, the lattice spacing was set to be able to make contact between the lattice results and their experimental complement and to be able to perform a continuum extrapolation. The light quark mass has been computed and found to be $m_{ud}^{\overline{\text{MS}}}(2\text{ GeV}) = 3.03(17)(38)\text{ MeV}$. This value is in good agreement with values from experiments and other lattice determinations.\\rnElectro-magnetic and axial form factors of the nucleon have been calculated. From these form factors the nucleon radii and the coupling constants were computed. The different ensembles enabled us to investigate systematically the dependence of these quantities on the volume, the lattice spacing and the pion mass.\newpage Finally we perform a continuum extrapolation and chiral extrapolations to the physical point.\\rnIn addition, we investigated so called excited state contributions to these observables. A technique was used, the summation method, which reduces these effects significantly and a much better agreement with experimental data was achieved. On the lattice, the Dirac radius and the axial charge are usually found to be much smaller than the experimental values. However, due to the carefully investigation of all the afore-mentioned systematic effects we get $\langle r_1^2\rangle_{u-d}=0.627(54)\text{ fm}^2$ and $g_A=1.218(92)$, which is in agreement with the experimental values within the errors.rnrnThe first three chapters introduce the theoretical background of form factors of the nucleon and lattice QCD in general. In chapter four the lattice spacing is determined. The computation of nucleon form factors is described in chapter five where systematic effects are investigated. All results are presented in chapter six. The thesis ends with a summary of the results and identifies options to complement and extend the calculations presented. rn

Composition of hydrothermal manganese oxide deposits from Galapagos mounds, DSDP Leg 70, Hole 509B

During DSDP Leg 70, a 1.60 m thick manganese oxide layer was sampled in hole 509B. This deposit is formed of alternating layers of hard plates of pure todorokite, about 2 mm thick, and of a more powdery material deeply impregnated with manganese oxide, about 3 mm thick. A SEM study of the plates and the associated powder shows that the powdery material is a transformation of a pre-existing sediment, while the plates are a direct precipitation from a hydrothermal solution. The uranium series disequilibrium method was used to determine the ages of the plates. They are found to be in good chronological sequence and in accordance with the sedimentation rate of the area (4.9 cm/10^3 years) which implies that they have been formed at the sediment-seawater interface during a pulsed injection of hydrothermal solution. The powder presents systematically an "older age" which is explained by a slowing down of the injection while the normal sediment settles; the older age is due to the 230Th excess of the sediment.

Linear Approach to the Orbiting Spacecraft Thermal Problem

A linear method is developed for solving the nonlinear differential equations of a lumped-parameter thermal model of a spacecraft moving in a closed orbit. This method, based on perturbation theory, is compared with heuristic linearizations of the same equations. The essential feature of the linear approach is that it provides a decomposition in thermal modes, like the decomposition of mechanical vibrations in normal modes. The stationary periodic solution of the linear equations can be alternately expressed as an explicit integral or as a Fourier series. This method is applied to a minimal thermal model of a satellite with ten isothermal parts (nodes), and the method is compared with direct numerical integration of the nonlinear equations. The computational complexity of this method is briefly studied for general thermal models of orbiting spacecraft, and it is concluded that it is certainly useful for reduced models and conceptual design but it can also be more efficient than the direct integration of the equations for large models. The results of the Fourier series computations for the ten-node satellite model show that the periodic solution at the second perturbative order is sufficiently accurate.

Characterizing chaos in a type of fractional duffing's equation

We characterize the chaos in a fractional Duffing’s equation computing the Lyapunov exponents and the dimension of the strange attractor in the effective phase space of the system. We develop a specific analytical method to estimate all Lyapunov exponents and check the results with the fiduciary orbit technique and a time series estimation method.

Numerical resolution of Emden's equation using Adomian polynomials

Purpose: In this paper the authors aim to show the advantages of using the decomposition method introduced by Adomian to solve Emden's equation, a classical non‐linear equation that appears in the study of the thermal behaviour of a spherical cloud and of the gravitational potential of a polytropic fluid at hydrostatic equilibrium. Design/methodology/approach: In their work, the authors first review Emden's equation and its possible solutions using the Frobenius and power series methods; then, Adomian polynomials are introduced. Afterwards, Emden's equation is solved using Adomian's decomposition method and, finally, they conclude with a comparison of the solution given by Adomian's method with the solution obtained by the other methods, for certain cases where the exact solution is known. Findings: Solving Emden's equation for n in the interval [0, 5] is very interesting for several scientific applications, such as astronomy. However, the exact solution is known only for n=0, n=1 and n=5. The experiments show that Adomian's method achieves an approximate solution which overlaps with the exact solution when n=0, and that coincides with the Taylor expansion of the exact solutions for n=1 and n=5. As a result, the authors obtained quite satisfactory results from their proposal. Originality/value: The main classical methods for obtaining approximate solutions of Emden's equation have serious computational drawbacks. The authors make a new, efficient numerical implementation for solving this equation, constructing iteratively the Adomian polynomials, which leads to a solution of Emden's equation that extends the range of variation of parameter n compared to the solutions given by both the Frobenius and the power series methods.

An FDTD model for calculation of gradient-induced eddy currents in MRI system

In most magnetic resonance imaging (MRI) systems, pulsed magnetic gradient fields induce eddy currents in the conducting structures of the superconducting magnet. The eddy currents induced in structures within the cryostat are particularly problematic as they are characterized by long time constants by virtue of the low resistivity of the conductors. This paper presents a three-dimensional (3-D) finite-difference time-domain (FDTD) scheme in cylindrical coordinates for eddy-current calculation in conductors. This model is intended to be part of a complete FDTD model of an MRI system including all RF and low-frequency field generating units and electrical models of the patient. The singularity apparent in the governing equations is removed by using a series expansion method and the conductor-air boundary condition is handled using a variant of the surface impedance concept. The numerical difficulty due to the asymmetry of Maxwell equations for low-frequency eddy-current problems is circumvented by taking advantage of the known penetration behavior of the eddy-current fields. A perfectly matched layer absorbing boundary condition in 3-D cylindrical coordinates is also incorporated. The numerical method has been verified against analytical solutions for simple cases. Finally, the algorithm is illustrated by modeling a pulsed field gradient coil system within an MRI magnet system. The results demonstrate that the proposed FDTD scheme can be used to calculate large-scale eddy-current problems in materials with high conductivity at low frequencies.

Time delay estimation with hidden Markov models

Most traditional methods for extracting the relationships between two time series are based on cross-correlation. In a non-linear non-stationary environment, these techniques are not sufficient. We show in this paper how to use hidden Markov models to identify the lag (or delay) between different variables for such data. Adopting an information-theoretic approach, we develop a procedure for training HMMs to maximise the mutual information (MMI) between delayed time series. The method is used to model the oil drilling process. We show that cross-correlation gives no information and that the MMI approach outperforms maximum likelihood.

O sistema de informação do tribunal colectivo dos géneros alimentícios: Estudo orgânico-funcional, produção e tramitação documental

Este trabalho é um estudo sobre o sistema de informação Tribunal Colectivo dos Géneros Alimentícios, como nasceu e como se desenvolveu, através do estudo orgânico-funcional, análise da produção e tramitação documental, desde a sua criação (1930) até à sua extinção (1976), através da aplicação do método de investigação quadripolar. Este método vai desenrolar-se mediante três actividades, cujo exercício se integra numa abordagem sistémica: a) Estudo da natureza da instituição a partir da análise dos seus fins, dos seus órgãos, das suas funções e da sua prática administrativa, tendo em conta tanto a evolução destes elementos ao longo do tempo como as relações que se estabelecem entre eles, através da análise de leis orgânicas e outros diplomas legais. b) Identificação e delimitação das séries documentais a partir da análise dos modos de produção documental, as características externas e internas destes tipos de documentos e as operações de conservação, descrição e organização arquivística. c) Apresentação do quadro de classificação: estrutura hierárquica e lógica que reflicta as funções deste sistema de informação. Pretende-se dar resposta às seguintes questões: 1- como classificar arquivisticamente a documentação acumulada no Arquivo Tribunal Colectivo dos Géneros Alimentícios? 2- Classificar arquivisticamente os processos de infracções/crimes do Tribunal de Recurso da Intendência-Geral da Segurança Pública e os processos de infracções/crimes até 1936 no sistema Tribunal Colectivo dos Géneros Alimentícios. 3- A produção e o circuito documental do Tribunal Colectivo dos Géneros Alimentícios são semelhantes à produção documental dos tribunais de jurisdição ordinária? 4- Quais a(s) especificidade(s) a nível processual e de tipologias documentais. ABSTRACT: This thesis is a study on the information system of the Tribunal Colectivo dos Géneros Alimentícios: how it started and how it was developed, and we do so proceding with the study of the organizational system, of the channeling of the proceedings and documents, since its creation (1930) until its ending (1976), using the quadripolar research method. This method will develop into three different steps, following a systemic approach and using the Quadripolar Method: a) A study of the nature of the institution supported by the analysis of its purposes, its organs, its functions and its administrative practices, bearing in mind the evolution of these elements throughout time and the connections that are established between them, with the support of some organic laws and other legal diplomas. b) Identification and caracterization of the document series; the method that it will be used will include an analysis of the process of document production, and also of both the internal and external characteristics of that type of documents and the archival conservation, description and organization. c) Demonstration of the classification board: hierarchic structure and the logic that is behind the institution of this documental base. This work intends to answer the following questions: 1. How to classify archival documentation accumulated in the Tribunal Colectivo dos Géneros Alimentícios? 2. How to classify the offense/crimes processes at the Tribunal de Recurso at Intendência-Geral da Segurança Pública and the offense/crimes processes until 1936 at Tribunal Colectivo dos Géneros Alimentícios? 3. Are the circuit and the document production of the Tribunal Colectivo dos Géneros Alimentícios similar to the production of documents in regular jurisdiction courts? 4. What are the specific level procedural and documentary types?

Axial buckling load of partially encased columns under fire

Partially encased columns have significant fire resistant. However, it is not possible to assess the fire resistance of such members simply by considering the temperature of the steel. The presence of concrete increases the mass and thermal inertia of the member and the variation of temperature within the cross section, in both the steel and concrete components. The annex G of EN1994-1-2 allows to calculate the load carrying capacity of partially encased columns, for a specific fire rating time, considering the balanced summation method. New formulas will be used to calculate the plastic resistance to axial compression and the effective flexural stiffness. These two parameters are used to calculate the buckling resistance. The finite element method is used to compare the results of the elastic critical load for different fire ratings of 30 and 60 minutes. The buckling resistance is also calculated by the finite element method, using an incremental and iterative procedure. This buckling resistance is also compared with the simple calculation method, evaluating the design buckling curve that best fits the results.

Biorthogonal series method for Oseen type flows

A biorthogonal series method is developed to solve Oseen type flow problems. The theory leads to a new set of eigenfunctions for a specific class of linear non-selfadjoint operators containing the biharmonic one. These eigenfunctions differ from those given earlier in the literature for the biharmonic operator. The method is applied to the problem of thermocapillary flow in a cylindrical liquid bridge of finite length with axial through flow. Flow and temperature distributions are obtained at leading order of an expansion for small surface tension Reynolds number and Prandtl number. Another related problem considered is that of cylindrical cavity flow. Solutions for both cases are presented in terms of biorthogonal series. The effect of axial through flow on velocity and temperature fields is discussed by numerical evaluation of the truncated analytical series. The presence of axial through flow not only convectively shifts the vortices induced by surface forces in the direction of the through flow, but also moves their centers toward the outer cylindrical boundary. This process can lead to significantly asymmetric flow structures.

Block coherence: a method for measuring the interdependence between two blocks of neurobiological time series

Multisensor recordings are becoming commonplace. When studying functional connectivity between different brain areas using such recordings, one defines regions of interest, and each region of interest is often characterized by a set (block) of time series. Presently, for two such regions, the interdependence is typically computed by estimating the ordinary coherence for each pair of individual time series and then summing or averaging the results over all such pairs of channels (one from block 1 and other from block 2). The aim of this paper is to generalize the concept of coherence so that it can be computed for two blocks of non-overlapping time series. This quantity, called block coherence, is first shown mathematically to have properties similar to that of ordinary coherence, and then applied to analyze local field potential recordings from a monkey performing a visuomotor task. It is found that an increase in block coherence between the channels from V4 region and the channels from prefrontal region in beta band leads to a decrease in response time.

An ultra-short-term wind power prediction method using “offline classification and optimization, online model matching” based on time series features

The applicability of ultra-short-term wind power prediction (USTWPP) models is reviewed. The USTWPP method proposed extracts featrues from historical data of wind power time series (WPTS), and classifies every short WPTS into one of several different subsets well defined by stationary patterns. All the WPTS that cannot match any one of the stationary patterns are sorted into the subset of nonstationary pattern. Every above WPTS subset needs a USTWPP model specially optimized for it offline. For on-line application, the pattern of the last short WPTS is recognized, then the corresponding prediction model is called for USTWPP. The validity of the proposed method is verified by simulations.