Geometric integration applied to moving mesh methods and degenerate Lagrangians

Moving mesh methods (also called r-adaptive methods) are space-adaptive strategies used for the numerical simulation of time-dependent partial differential equations. These methods keep the total number of mesh points fixed during the simulation, but redistribute them over time to follow the areas where a higher mesh point density is required. There are a very limited number of moving mesh methods designed for solving field-theoretic partial differential equations, and the numerical analysis of the resulting schemes is challenging. In this thesis we present two ways to construct r-adaptive variational and multisymplectic integrators for (1+1)-dimensional Lagrangian field theories. The first method uses a variational discretization of the physical equations and the mesh equations are then coupled in a way typical of the existing r-adaptive schemes. The second method treats the mesh points as pseudo-particles and incorporates their dynamics directly into the variational principle. A user-specified adaptation strategy is then enforced through Lagrange multipliers as a constraint on the dynamics of both the physical field and the mesh points. We discuss the advantages and limitations of our methods. The proposed methods are readily applicable to (weakly) non-degenerate field theories---numerical results for the Sine-Gordon equation are presented.

In an attempt to extend our approach to degenerate field theories, in the last part of this thesis we construct higher-order variational integrators for a class of degenerate systems described by Lagrangians that are linear in velocities. We analyze the geometry underlying such systems and develop the appropriate theory for variational integration. Our main observation is that the evolution takes place on the primary constraint and the 'Hamiltonian' equations of motion can be formulated as an index 1 differential-algebraic system. We then proceed to construct variational Runge-Kutta methods and analyze their properties. The general properties of Runge-Kutta methods depend on the 'velocity' part of the Lagrangian. If the 'velocity' part is also linear in the position coordinate, then we show that non-partitioned variational Runge-Kutta methods are equivalent to integration of the corresponding first-order Euler-Lagrange equations, which have the form of a Poisson system with a constant structure matrix, and the classical properties of the Runge-Kutta method are retained. If the 'velocity' part is nonlinear in the position coordinate, we observe a reduction of the order of convergence, which is typical of numerical integration of DAEs. We also apply our methods to several models and present the results of our numerical experiments.

Experimental Investigation of Thrust Fault Rupture Mechanics

Thrust fault earthquakes are investigated in the laboratory by generating dynamic shear ruptures along pre-existing frictional faults in rectangular plates. A considerable body of evidence suggests that dip-slip earthquakes exhibit enhanced ground motions in the acute hanging wall wedge as an outcome of broken symmetry between hanging and foot wall plates with respect to the earth surface. To understand the physical behavior of thrust fault earthquakes, particularly ground motions near the earth surface, ruptures are nucleated in analog laboratory experiments and guided up-dip towards the simulated earth surface. The transient slip event and emitted radiation mimic a natural thrust earthquake. High-speed photography and laser velocimeters capture the rupture evolution, outputting a full-field view of photo-elastic fringe contours proportional to maximum shearing stresses as well as continuous ground motion velocity records at discrete points on the specimen. Earth surface-normal measurements validate selective enhancement of hanging wall ground motions for both sub-Rayleigh and super-shear rupture speeds. The earth surface breaks upon rupture tip arrival to the fault trace, generating prominent Rayleigh surface waves. A rupture wave is sensed in the hanging wall but is, however, absent from the foot wall plate: a direct consequence of proximity from fault to seismometer. Signatures in earth surface-normal records attenuate with distance from the fault trace. Super-shear earthquakes feature greater amplitudes of ground shaking profiles, as expected from the increased tectonic pressures required to induce super-shear transition. Paired stations measure fault parallel and fault normal ground motions at various depths, which yield slip and opening rates through direct subtraction of like components. Peak fault slip and opening rates associated with the rupture tip increase with proximity to the fault trace, a result of selective ground motion amplification in the hanging wall. Fault opening rates indicate that the hanging and foot walls detach near the earth surface, a phenomenon promoted by a decrease in magnitude of far-field tectonic loads. Subsequent shutting of the fault sends an opening pulse back down-dip. In case of a sub-Rayleigh earthquake, feedback from the reflected S wave re-ruptures the locked fault at super-shear speeds, providing another mechanism of super-shear transition.

Edge function method applied to thin plates resting on Winkler foundation

The Edge Function method formerly developed by Quinlan⁽²⁵⁾ is applied to solve the problem of thin elastic plates resting on spring supported foundations subjected to lateral loads the method can be applied to plates of any convex polygonal shapes, however, since most plates are rectangular in shape, this specific class is investigated in this thesis. The method discussed can also be applied easily to other kinds of foundation models (e.g. springs connected to each other by a membrane) as long as the resulting differential equation is linear. In chapter VII, solution of a specific problem is compared with a known solution from literature. In chapter VIII, further comparisons are given. The problems of concentrated load on an edge and later on a corner of a plate as long as they are far away from other boundaries are also given in the chapter and generalized to other loading intensities and/or plates springs constants for Poisson's ratio equal to 0.2

Condições mecânicas do sistema respiratório em adultos com fibrose cística: estudo pela espirometria, pletismografia e técnica de oscilações forçadas

A fibrose cística (FC) é a doença autossômica recessiva mais comum na população branca que leva à redução na expectativa de vida. A doença pulmonar é a maior causa de morbidade e mortalidade. A relevância do presente estudo se dá diante de alguns fatores: aumento drástico da sobrevida média nos últimos 60 anos na FC, a fisiopatologia pulmonar não é bem compreendida, ausência de estudos reportados na literatura, até o momento, utilizando a técnica de oscilações forçadas (TOF) exclusivamente em adultos com FC. Assim sendo os objetivos deste estudo são: analisar as alterações da mecânica respiratória em adultos com FC através da espirometria, pletismografia e TOF; correlacionar os resultados da TOF aos espirométricos e pletismográficos e avaliar a sensibilidade e especificidade da TOF nestes indivíduos. É um estudo de corte transversal descritivo, no qual foram analisados dois grupos de indivíduos: controle (n=23) e FC (n=27). Os resultados foram expressos através média desvio-padrão. As técnicas funcionais respiratórias foram realizadas na seguinte sequência: TOF, espirometria, pletismografia. Na pletismografia foram avaliados os parâmetros: CPT (capacidade pulmonar total), CRF (capacidade residual funcional) e VR (volume residual), CRF/CPT e VR/CPT, resistência (Rva) e condutância específica das vias aéreas (SGva). Na espirometria: volume expiratório forçado no primeiro segundo (VEF1), capacidade vital forçada (CVF), fluxo expiratório entre 25% e 75% (FEF25%-75%) da CVF (FEF25%-75%) e razões VEF1/CVF (%) e FEF/CVF (%). Na TOF: propriedades resistivas do sistema respiratório- R0 (resistência no intercepto), Rm (resistência média) e S (inclinação da reta de resistência) e propriedades reativas: Cdin,sr (complacência dinâmica do sistema respiratório), Xm (reatância média), frequência de ressonância (fr); e o módulo da impedância em 4 Hz (׀Zrs4Hz׀). Na espirometria o distúrbio ventilatório obstrutivo (DVO) com CVF reduzida foi predominante, com marcante redução do FEF25%-75% no grupo FC (p<0,0001) em relação ao controle. Na pletismografia: destacou-se a elevação de VR, na presença de CPT normal e elevação da Rva e redução da SGva no grupo FC. Alterações da TOF ocorridas no grupo FC em relação ao controle: aumento de R0 e Rm (p<0,0001) e fr (p<0,0002), relacionados à obstrução das vias aéreas; redução de S (p<0,0006), Xm (p<0,0001) associadas à não-homogeneidade do sistema respiratório e Cdin,sr (p<0,0001), relacionada à redução da complacência pulmonar; aumento do módulo da impedância em 4 Hz (׀Zrs4Hz׀) (representando a carga mecânica total do sistema respiratório) resultante da interação das demais alterações da TOF citadas. Os parâmetros da TOF apresentaram correlações muito boas com a espirometria e moderadas com a pletismografia. Rm foi o único parâmetro que não se relacionou com nenhuma destas técnicas. A sensibilidade e especificidade da TOF em adultos com FC apresentaram valores elevados, sobretudo nos parâmetros reativos, em especial, Xm (85,2% e 73,9% respectivamente e área sob a curva de 0,86).

Advanced density functional theory methods for materials science

In this work we chiefly deal with two broad classes of problems in computational materials science, determining the doping mechanism in a semiconductor and developing an extreme condition equation of state. While solving certain aspects of these questions is well-trodden ground, both require extending the reach of existing methods to fully answer them. Here we choose to build upon the framework of density functional theory (DFT) which provides an efficient means to investigate a system from a quantum mechanics description.

Zinc Phosphide (Zn₃P₂) could be the basis for cheap and highly efficient solar cells. Its use in this regard is limited by the difficulty in n-type doping the material. In an effort to understand the mechanism behind this, the energetics and electronic structure of intrinsic point defects in zinc phosphide are studied using generalized Kohn-Sham theory and utilizing the Heyd, Scuseria, and Ernzerhof (HSE) hybrid functional for exchange and correlation. Novel 'perturbation extrapolation' is utilized to extend the use of the computationally expensive HSE functional to this large-scale defect system. According to calculations, the formation energy of charged phosphorus interstitial defects are very low in n-type Zn₃P₂ and act as 'electron sinks', nullifying the desired doping and lowering the fermi-level back towards the p-type regime. Going forward, this insight provides clues to fabricating useful zinc phosphide based devices. In addition, the methodology developed for this work can be applied to further doping studies in other systems.

Accurate determination of high pressure and temperature equations of state is fundamental in a variety of fields. However, it is often very difficult to cover a wide range of temperatures and pressures in an laboratory setting. Here we develop methods to determine a multi-phase equation of state for Ta through computation. The typical means of investigating thermodynamic properties is via ’classical’ molecular dynamics where the atomic motion is calculated from Newtonian mechanics with the electronic effects abstracted away into an interatomic potential function. For our purposes, a ’first principles’ approach such as DFT is useful as a classical potential is typically valid for only a portion of the phase diagram (i.e. whatever part it has been fit to). Furthermore, for extremes of temperature and pressure quantum effects become critical to accurately capture an equation of state and are very hard to capture in even complex model potentials. This requires extending the inherently zero temperature DFT to predict the finite temperature response of the system. Statistical modelling and thermodynamic integration is used to extend our results over all phases, as well as phase-coexistence regions which are at the limits of typical DFT validity. We deliver the most comprehensive and accurate equation of state that has been done for Ta. This work also lends insights that can be applied to further equation of state work in many other materials.

Suspension mechanics: I. Inertial and non-Newtonian migration of neutrally buoyant rigid spheres in two-dimensional unidirectional flows; II. The effect of viscoelasticity on the creeping motion of a train of neutrally buoyant Newtonian drops through a circular tube

The lateral migration of neutrally buoyant rigid spheres in two-dimensional unidirectional flows was studied theoretically. The cases of both inertia-induced migration in a Newtonian fluid and normal stress-induced migration in a second-order fluid were considered. Analytical results for the lateral velocities were obtained, and the equilibrium positions and trajectories of the spheres compared favorably with the experimental data available in the literature. The effective viscosity was obtained for a dilute suspension of spheres which were simultaneously undergoing inertia-induced migration and translational Brownian motion in a plane Poiseuille flow. The migration of spheres suspended in a second-order fluid inside a screw extruder was also considered.

The creeping motion of neutrally buoyant concentrically located Newtonian drops through a circular tube was studied experimentally for drops which have an undeformed radius comparable to that of the tube. Both a Newtonian and a viscoelastic suspending fluid were used in order to determine the influence of viscoelasticity. The extra pressure drop due to the presence of the suspended drops, the shape and velocity of the drops, and the streamlines of the flow were obtained for various viscosity ratios, total flow rates, and drop sizes. The results were compared with existing theoretical and experimental data.

High-resolution automatic water quality monitoring systems applied to catchment and reservoir monitoring

Automatic recording instruments provide the ideal means of recording the responses of rivers, lakes and reservoirs to short-term changes in the weather. As part of the project ‘Using Automatic Monitoring and Dynamic Modelling for the Active Management of Lakes and Reservoirs', a family of three automatic monitoring stations were designed by engineers at the Centre for Ecology and Hydrology in Windermere to monitor such responses. In this article, the authors describe this instrument network in some detail and present case studies that illustrate the value of high resolution automatic monitoring in both catchment and reservoir applications.

Comparison of marketing policies applied by the same company in different countries and analysis of the reasons for the differences

The focus of this project is going to be, as the title indicates, on the comparison of marketing policies applied by the same company in different countries and analysis of the reasons for the differences. In order to do that, I have selected the company Nestlé to analyze the marketing decisions it makes across national boundaries to market the brand of Kit Kat and keep it as a leader snack worldwide. After having analyzed the brand in all continents, I can say the execution of the strategy used by Nestlé with Kit Kat really matches the planning of the strategy which is to think globally and act locally. Nestlé uses global brand identity but, from the internal point of view, it uses local ingredients and gives autonomy to its local branches based in different countries to make pricing and distributions decisions and therefore satisfy different consumers’ needs and preferences in local environments where changes happen very rapidly. The “glocal” approach to marketing is an effective way for Nestlé Kit Kat to stay focused on the consumer in worldwide markets.