There are two main approaches to solve the problem of finding closed geodesics on a Riemannian manifold M. The variational approach views a closed geodesic as a closed curve which happens to be a geodesic and it looks for critical points of the energy functional, while the dynamical systems approach views a closed geodesic as a geodesic which happens to close up and looks for periodic orbits of the geodesic ow on the unit tangent bundle.

A better, more optimized processing pipeline for functional connectivity (FC) data will likely accelerate practical advances within the field of neuroimaging. When using correlation-based measures of FC, researchers have recently employed a few data-driven methods to maximize its predictive power. In this study, we apply a few of these post-processing methods in both task, twin, and subject identification problems. First, we employ PCA reconstruction of the original dataset, which has been successfully used to maximize subject-level identifiability. We show there is dataset-dependent optimal PCA reconstruction for task and twin identification. Next, we analyze FCs in their native geometry using tangent space projection with various mean covariance reference matrices. We demonstrate that the tangent projection of the original FCs can drastically increase subject and twin identification rates. For example, the identification rate of 106 MZ twin pairs increased from 0.487 of the original FCs to 0.943 after tangent projection with the logarithmic Euclidean reference matrix. We also use Schaefer’s variable parcellation sizes to show that increasing parcellation granularity in general increases twin and subject identification rates. Finally, we show that our custom convolutional neural network classifier achieves an average task identification rate of 0.986, surpassing state-of-the-art results. These post-processing methods are promising for future research in functional connectome predictive modeling and, if optimized further, can likely be extended into clinical applications.

No description available.

In this thesis, formulas for the local tangent stiffness matrix of a plane frame member are derived by differentiating the member resistance vector in the displaced position. This approach facilitates an analysis using only one element per member. The formulas are checked by finite difference. The derivation leads to the familiar elastic and geometric stiffness matrices used by other authors plus an additional higher order geometric stiffness matrix. Contributions of each of the three sub-matrices to the tangent stiffness matrix are studied on both the member and structure levels through two numerical examples. These same examples are analyzed three different ways for comparison. First, the examples are analyzed using the method presented in this thesis. Second, they are analyzed with the finite element modeling software ABAQUS/CAE using only one element per member. Third, they are analyzed with ABAQUS using 200 elements per member. Comparisons are made assuming the ABAQUS analysis which uses 200 elements per member is the most accurate. The element presented in this thesis performs much better than the ABAQUS analysis which uses one element per member, with maximum errors of 1.0% and 40.8% respectively, for a cantilever column example. The maximum error for the two story frame example using the ABAQUS analysis with one element per member is 42.8%, while the results from the analysis using the element presented in this thesis are within 1.5%. Using the element presented in this thesis with only one element per member gives good and computationally efficient results for second-order analysis.

The primary purpose of this research is to investigate the effectiveness of a scaled roller rig for accurately assessing the contact mechanics and dynamics between a profiled steel wheel and rail, as is commonly used in rail vehicles. The established creep models of Kalker and Johnson and Vermeulen are used to establish correction factors, scaling factors, and transformation factors that allow us to relate the results from a scaled rig to those of a tangent track. �Correction factors, which are defined as the ratios of a given quantity (such as creep coefficient) between a roller rig and a track, are derived and used to relate the results between a full-size rig and a full-size track. Scaling factors are derived to relate the same quantities between roller rigs of different scales. Finally, transformation factors are derived by combining scaling factors with correction factors in order to relate the results from a scaled roller rig to a full-size tangent track. Close-end formulae for creep force correction, scaling, and transformation factors are provided in the thesis, along with their full derivation and an explanation of their limitations; these formulae can be used to calculate the correction factors for any wheel-rail geometry and scaling. For Kalker's theory, it is shown that the correction factor for creep coefficients is strictly a function of wheel and rail geometry, primarily the wheel and roller diameter ratio. For Johnson and Vermeulen's theory, the effects of creepage, scale, and load on the creep force correction factor are demonstrated. �It is shown that INRETS' scaling strategy causes the normalized creep curve to be identical for both a full-size and a scaled roller rig. �It is also shown that the creep force correction factors for Johnson and Vermeulen's model increase linearly with creepage, starting with the values predicted by Kalker's theory. �Therefore, Kalker's theory provides a conservative estimate for creep force correction factors. �A case study is presented to demonstrate the creep curves, as well as the correction and transformation factors, for a typical wheel-rail configuration. �Additionally, two studies by other authors that calculate the correction factor for Kalker's creep coefficients for specific wheel-rail geometries are reviewed and show full agreement with the results that are predicted by the formulae derived in this study. �Based on a review of existing and past roller rigs, as well as the findings of this thesis, a number of recommendations are given for the design of a roller rig for the purpose of assessing the wheel-rail contact mechanics. �A scaling strategy (INRETS') is suggested, and equations for power consumption of a roller rig are derived. Recommendations for sensors and actuators necessary for such a rig are also given. Special attention is given to the resolution and accuracy of velocity sensors, which are required to properly measure and plot the creep curves. / Master of Science

La modélisation des équilibres thermodynamiques entre phases est essentielle pour le génie des procédés et le génie pétrolier. L’analyse de la stabilité des phases est un problème de la plus haute importance parmi les calculs d’équilibre des phases. Le calcul de stabilité décide si un système se présente dans un état monophasique ou multiphasique ; si le système se sépare en deux ou plusieurs phases, les résultats du calcul de stabilité fournissent une initialisation de qualité pour les calculs de flash (Michelsen, 1982b), et permettent la validation des résultats des calculs de flash multiphasique. Le problème de la stabilité des phases est résolu par une minimisation sans contraintes de la fonction distance au plan tangent à la surface de l’énergie libre de Gibbs (« tangent plane distance », ou TPD). Une phase est considérée comme étant thermodynamiquement stable si la fonction TPD est non- négative pour tous les points stationnaires, tandis qu’une valeur négative indique une phase thermodynamiquement instable. La surface TPD dans l’espace compositionnel est non- convexe et peut être hautement non linéaire, ce qui fait que les calculs de stabilité peuvent être extrêmement difficiles pour certaines conditions, notamment aux voisinages des singularités. On distingue deux types de singularités : (i) au lieu de la limite du test de stabilité (stability test limit locus, ou STLL), et ii) à la spinodale (la limite intrinsèque de la stabilité thermodynamique). Du point de vue géométrique, la surface TPD présente un point selle, correspondant à une solution non triviale (à la STLL) ou triviale (à la spinodale). Dans le voisinage de ces singularités, le nombre d’itérations de toute méthode de minimisation augmente dramatiquement et la divergence peut survenir. Cet inconvénient est bien plus sévère pour la STLL que pour la spinodale. Le présent mémoire est structuré sur trois grandes lignes : (i) après la présentation du critère du plan tangent à la surface de l’énergie libre de Gibbs, plusieurs solutions itératives (gradient et méthodes d’accélération de la convergence, méthodes de second ordre de Newton et méthodes quasi- Newton), du problème de la stabilité des phases sont présentées et analysées, surtout du point de vue de leurs comportement près des singularités; (ii) Suivant l’analyse des valeurs propres, du conditionnement de la matrice Hessienne et de l’échelle du problème, ainsi que la représentation de la surface de la fonction TPD, la résolution du calcul de la stabilité des phases par la minimisation des fonctions coût modifiées est adoptée. Ces fonctions « coût » sont choisies de telle sorte que tout point stationnaire (y compris les points selle) de la fonction TPD soit converti en minimum global; la Hessienne à la STLL est dans ce cas positif définie, et non indéfinie, ce qui mène a une amélioration des propriétés de convergence, comme montré par plusieurs exemples pour des mélanges représentatifs, synthétiques et naturels. Finalement, (iii) les calculs de stabilité sont menés par une méthode d’optimisation globale, dite de Tunneling. La méthode de Tunneling consiste à détruire (en plaçant un pôle) les minima déjà trouvés par une méthode de minimisation locale, et a tunneliser pour trouver un point situé dans une autre vallée de la surface de la fonction coût qui contient un minimum 9 à une valeur plus petite de la fonction coût; le processus continue jusqu'à ce que les critères du minimum global soient remplis. Plusieurs exemples soigneusement choisis montrent la robustesse et l’efficacité de la méthode de Tunneling pour la minimisation de la fonction TPD, ainsi que pour la minimisation des fonctions coût modifiées. / The thermodynamic phase equilibrium modelling is an essential issue for petroleum and process engineering. Phase stability analysis is a highly important problem among phase equilibrium calculations. The stability computation establishes whether a given mixture is in one or several phases. If a mixture splits into two or more phases, the stability calculations provide valuables initialisation sets for the flash calculations, and allow the validation of multiphase flash calculations. The phase stability problem is solved as an unconstrained minimisation of the tangent plan distance (TPD) function to the Gibbs free energy surface. A phase is thermodynamically stable if the TPD function is non-negative at all its stationary points, while a negative value indicates an unstable case. The TPD surface is non-convex and may be highly non-linear in the compositional space; for this reason, phase stability calculation may be extremely difficult for certain conditions, mainly within the vicinity of singularities. One can distinguish two types of singularities: (i) the stability test limit locus (STLL), and (ii) the intrinsic limit of stability (spinodal). Geometrically, the TPD surface exhibits a saddle point, corresponding to a non-trivial (at the STLL) or trivial solution (at the spinodal). In the immediate vicinity of these singularities, the number of iterations of all minimisation methods increases dramatically, and divergence could occur. This inconvenient is more severe for the STLL than for the spinodal. The work presented herein is structured as follow: (i) after the introduction to the concept of tangent plan distance to the Gibbs free energy surface, several iterative methods (gradient, acceleration methods, second-order Newton and quasi-Newton) are presented, and their behaviour analysed, especially near singularities. (ii) following the analysis of Hessian matrix eigenvalues and conditioning, of problem scaling, as well as of the TPD surface representation, the solution of phase stability computation using modified objective functions is adopted. The latter are chosen in such a manner that any stationary point of the TPD function becomes a global minimum of the modified function; at the STLL, the Hessian matrix is no more indefinite, but positive definite. This leads to a better scheme of convergence as will be shown in various examples for synthetic and naturally occurring mixtures. Finally, (iii) the so-called Tunneling global optimization method is used for the stability analysis. This method consists in destroying the minima already found (by placing poles), and to tunnel to another valley of the modified objective function to find a new minimum with a smaller value of the objective function. The process is resumed when criteria for the global minimum are fulfilled. Several carefully chosen examples demonstrate the robustness and the efficiency of the Tunneling method to minimize the TPD function, as well as the modified objective functions.

Computed tomography is a very useful technic which allow non-invasive diagnosis in many applications for example is used with success in industry and medicine. However, manual analysis of the interesting structures can be tedious and extremely time consuming, or even impossible due its complexity. Therefore in this thesis we study and develop discrete geometry and topology algorithms suitable for use in many practical applications, especially, in the problem of automatic quantitative analysis of the human airway trees based on computed tomography images. In the first part, we define basic notions used in discrete topology and geometry then we showed that several class of discrete methods like skeletonisation algorithms, medial axes, tunnels closing algorithms and tangent estimators, are widely used in several different practical application. The second part consist of a proposition and theory of a new methods for solving particular problems. We introduced two new medial axis filtering method. The hierarchical scale medial axis which is based on previously proposed scale axis transform, however, is free of drawbacks introduced in the previously proposed method and the discrete adaptive medial axis where the filtering parameter is dynamically adapted to the local size of the object. In this part we also introduced an efficient and parameter less new tangent estimators along three-dimensional discrete curves, called 3D maximal segment tangent direction. Finally, we showed that discrete geometry and topology algorithms can be useful in the problem of quantitative analysis of the human airway trees based on computed tomography images. According to proposed in the literature design of such system we applied discrete topology and geometry algorithms to solve particular problems at each step of the quantitative analysis process. First, we propose a robust method for segmenting airway tree from CT datasets. The method is based on the tunnel closing algorithm and is used as a tool to repair, damaged by acquisition errors, CT images. We also proposed an algorithm for creation of an artificial model of the bronchial tree and we used such model to validate algorithms presented in this work. Then, we compare the quality of different algorithms using set of experiments conducted on computer phantoms and real CT dataset. We show that recently proposed methods which works in cubical complex framework, together with methods introduced in this work can overcome problems reported in the literature and can be a good basis for the further implementation of the system for automatic quantification of bronchial tree properties

Estudaremos as singularidades genéricas de superfécies regradas em R3. O objetivo do trabalho é mostrar que as singularidades genéricas que ocorrem no conjunto das superfícies regradas são as mesmas que ocorrem no conjunto das aplicações diferenciáveis de R2 em R3, enquanto que as singularidades genéricas das superfícies desenvolvíveis, que formam um subconjunto das superfícies regradas, são mais degeneradas. / We study generic singularities of ruled surfaces in R3. In this work we show that generic singularities appearing in the set of ruled surfaces are the same that occur in the set of map germs from R2 to R3, while the generic singularities of developable surfaces are more degenerate.

Neste trabalho estudamos as aplicações do método do Gradiente Espectral Projetado (SPG) em Meteorologia nos campos de previsibilidade, estabilidade e sensibilidade. Inicialmente revisamos os Vetores Singulares Lineares (LSVs) e em seguida apresentamos a teoria das Perturbações Condicionais Não-Lineares Ótimas (CNOPs). Enquanto os métodos clássicos estão baseados no Modelo Tangente Linear, as CNOPs são uma formulação do mesmo problema baseado em Programação Não-Linear. As CNOPs são descritas na literatura como responsáveis por melhorias em relação aos métodos anteriores. Finalmente analisamos três exemplos de aplicação do método à problemas de previsibilidade, estabilidade e sensibilidade. / A revision about applications of Spectral Projected Gradient (SPG) in meteorology is done in the fields of predictability, stability and sensitivity. Initially we review about Linear Singular Vectos (LSVs) and we present the Conditional Nonlinear Optimal perturbations (CNOPs). While the classic methods are based on the Tangent Linear Model, CNOPs are another formulation of the problem based on Nonlinear Programming. CNOPs are described in bibliography as responsible by better results than older methods. Finally we analyze three applications in predictability, stability and sensibility.

The bulk modulus of pure hydraulic oil and its dependency on pressure and temperature has been studied extensively over the past years. A comprehensive review of some of the more common definitions of fluid bulk modulus is conducted and comments on some of the confusion over definitions and different methods of measuring the fluid bulk modulus are presented in this thesis. In practice, it is known that there is always some form of air present in hydraulic systems which substantially decreases the oil bulk modulus. The term effective bulk modulus is used to account for the effect of air and/or the compliance of transmission lines. A summary from the literature of the effective bulk modulus models for a mixture of hydraulic oil and air is presented. Based on the reviews, these models are divided into two groups: “compression only” models and “compression and dissolve” models. A comparison of various “compression only” models, where only the volumetric compression of air is considered, shows that the models do not match each other at the same operating conditions. The reason for this difference is explained and after applying some modifications to the models, a theoretical model of the “compression only” model is suggested. The “compression and dissolve” models, obtained from the literature review, include the effects of the volumetric compression of air and the volumetric reduction of air due to the dissolving of air into the oil. It is found that the existing “compression and dissolve” models have a discontinuity at some critical pressure and as a result do not match the experimental results very well. The reason for the discontinuity is discussed and a new “compression and dissolve” model is proposed by introducing some new parameters to the theoretical model. A new critical pressure (PC) definition is presented based on the saturation limit of oil. In the new definition, the air stops dissolving into the oil after this critical pressure is reached and any remaining air will be only compressed afterwards. An experimental procedure is successfully designed and fabricated to verify the new proposed models and to reproduce the operating conditions that underlie the model assumptions. The pressure range is 0 to 6.9 MPa and the temperature is kept constant at °C. Air is added to the oil in different forms and the amount of air varies from about 1 to 5%. Experiments are conducted in three different phases: baseline (without adding air to the oil), lumped air (air added as a pocket of air to the top of the oil column) and distributed air (air is distributed in the oil in the form of small air bubbles). The effect of different forms and amounts of air and various volume change rates are investigated experimentally and it is shown that the value of PC is strongly affected by the volume change rate, the form, and the amount of air. It is also shown that the new model can represent the experimental data with great accuracy. The new proposed “compression and dissolve” model can be considered as a general model of the effective bulk modulus of a mixture of oil and air where it is applicable to any form of a mixture of hydraulic oil and air. However, it is required to identify model parameters using experimental measurements. A method of identifying the model parameters is introduced and the modeling errors are evaluated. An attempt is also made to verify independently the value of some of the parameters. The new proposed model can be used in analyzing pressure variations and improving the accuracy of the simulations in low pressure hydraulic systems. The new method of modeling the air dissolving into the oil can be also used to improve the modeling of cavitation phenomena in hydraulic systems.