Part 1: Theoretic background information - Review of Hooke’s law for linear elastic materials - The strain energy density of linear elastic materials - Hyperelastic material - Material laws for hyperelastic materials - About selecting the material model and performing tests - Implementation of hyperelastic material laws in Mechanica - Defining hyperelastic material parameters in Mechanica - Test set-ups and specimen shapes of the supported material tests - The uniaxial compression test - Stress and strain definitions in the Mechanica LDA analysis Part 2: Application examples - A test specimen subjected to uniaxial loading - A volumetric compression test - A planar test - Influence of the material law Appendix - PTC Simulation Services Introduction - Dictionary Technical English-German / Teil 1: Theoretische Hintergrundinformation - Das Hookesche Gesetz für linear-elastische Werkstoffe - Die Dehnungsenergiedichte für linear-elastische Materialien - Hyperelastisches Material - Materialgesetze für Hyperelastizität - Auswählen des Materialgesetzes und Testdurchführung - Implementierung der hyperelastischen Materialgesetze in Mechanica - Definieren der hyperelastischen Materialparameter in Mechanica - Testaufbauten und Prüfkörper der unterstützten Materialtests - Der einachsige Druckversuch - Spannungs- und Dehnungsdefinition in der Mechanica-Analyse mit großen Verformungen Teil 2: Anwendungsbeispiele - Ein einachsig beanspruchter Prüfkörper - Ein volumetrischer Drucktest - Ein planarer Test - Einfluss des Materialgesetzes Anhang: - Kurzvorstellung der PTC Simulationsdienstleistungen - Wörterbuch technisches Englisch-Deutsch

We present a variational method for unfolding of the cortex based on a user-chosen point of view as an alternative to more traditional global flattening methods, which incur more distortion around the region of interest. Our approach involves three novel contributions. The first is an energy function and its corresponding gradient flow to measure the average visibility of a region of interest of a surface from a given viewpoint. The second is an additional energy function and flow designed to preserve the 3D topology of the evolving surface. This latter contribution receives significant focus in this thesis as it is crucial to obtain the desired unfolding effect derived from the first energy functional and flow. Without it, the resulting topology changes render the unconstrained evolution uninteresting for the purpose of cortical visualization, exploration, and inspection. The third is a method that dramatically improves the computational speed of the 3D topology-preservation approach by creating a tree structure of the triangulated surface and using a recursion technique.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Este trabalho apresenta uma família de novos métodos numéricos euleriano-lagrangeanos localmente conservativos para leis de conservação hiperbólicas escalares. Estes métodos não utilizam soluções analíticas de problemas de Riemann e são bastante precisos na captura de saltos nas soluções. Estes métodos foram introduzidos, implementados computacionalmente e testados para leis de conservação em uma e duas dimensões espaciais. Foram consideradas as equações de Burgers e Buckley-Leverett. Nossos experimentos numéricos indicaram que os métodos são pouco difusivos e que as soluções não apresentam oscilações espúrias.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Este trabalho apresenta uma família de novos métodos numéricos euleriano-lagrangeanos localmente conservativos para leis de conservação hiperbólicas escalares. Estes métodos não utilizam soluções analíticas de problemas de Riemann e são bastante precisos na captura de saltos nas soluções. Estes métodos foram introduzidos, implementados computacionalmente e testados para leis de conservação em uma e duas dimensões espaciais. Foram consideradas as equações de Burgers e Buckley-Leverett. Nossos experimentos numéricos indicaram que os métodos são pouco difusivos e que as soluções não apresentam oscilações espúrias.

Lie-detection has been an evergreen and evolving subject. Polygraph techniques have been the most popular and successful technique till date. The main drawback of the polygraph is that good results cannot be attained without maintaining a physical contact, of the subject under test. In general, this physical contact would induce extra consciousness in the subject. Also, any sort of arousal in the subject triggers false positives while performing the traditional polygraph based tests. With all these drawbacks in the polygraph, also, due to rapid developments in the fields of computer vision and artificial intelligence, with newer and faster algorithms, have compelled mankind to search and adapt to contemporary methods in lie-detection. Observing the facial expressions of emotions in a person without any physical contact and implementing these techniques using artificial intelligence is one such method. The concept of magnifying a micro expression and trying to decipher them is rather premature at this stage but would evolve in future. Magnification using EVM technique has been proposed recently and it is rather new to extract these micro expressions from magnified EVM based on HOG features. Till date, HOG features have been used in conjunction with SVM, and generally for person/pedestrian detection. A newer, simpler and contemporary method of applying EVM with HOG features and Back-propagation Neural Network jointly has been introduced and proposed to extract and decipher the micro-expressions on the face. Micro-expressions go unnoticed due to its involuntary nature, but EVM is used to magnify them and makes them noticeable. Emotions behind the micro-expressions are extracted and recognized using the HOG features \& Back-Propagation Neural Network. One of the important aspects that has to be dealt with human beings is a biased mind. Since, an investigator is also a human and, he too, has to deal with his own assumptions and emotions, a Neural Network is used to give the investigator an unbiased start in identifying the true emotions behind every micro-expression. On the whole, this proposed system is not a lie-detector, but helps in detecting the emotions of the subject under test. By further investigation, a lie can be detected. / This thesis uses a magnification technique to magnify the subtle, faint and spontaneous facial muscle movements or more precisely, micro-expressions. This magnification would help a system in classifying them and estimating the emotion behind them. This technique additionally magnifies the color changes, which could be used to extract the pulse without a physical contact with the subject. The results are presented in a GUI. / Gautam: +46(0)739528573, +91-9701534064 Tushal: +46(0)723219833, +91-9000242241 Venu: +46(0)734780266, +91-9298653191 Sai: +91-9989410111

Problems governed by partial differential equations (PDEs) in deformable domains, t Rd; d = 2; 3; are of fundamental importance in science and engineering. They are of particular relevance in the design of many engineering systems e.g., aircrafts and bridges as well as to the analysis of several biological phenomena e.g., blood ow in arteries. However, developing numerical scheme for such problems is still very challenging even when the deformation of the boundary of domain is prescribed a priori. Possibility of excessive mesh distortion is one of the major challenge when solving such problems with numerical methods using boundary tted meshes. The arbitrary Lagrangian- Eulerian (ALE) approach is a way to overcome this difficulty. Numerical simulations of convection-dominated problems have for long been the subject to many researchers. Galerkin formulations, which yield the best approximations for differential equations with high diffusivity, tend to induce spurious oscillations in the numerical solution of convection dominated equations. Though such spurious oscillations can be avoided by adaptive meshing, which is computationally very expensive on ne grids. Alternatively, stabilization methods can be used to suppress the spurious oscillations. In this work, the considered equation is designed within the framework of ALE formulation. In the first part, Streamline Upwind Petrov-Galerkin (SUPG) finite element method with conservative ALE formulation is proposed. Further, the first order backward Euler and the second order Crank-Nicolson methods are used for the temporal discretization. It is shown that the stability of the semi-discrete (continuous in time) ALE-SUPG equation is independent of the mesh velocity, whereas the stability of the fully discrete problem is unconditionally stable for implicit Euler method and is only conditionally stable for Crank-Nicolson time discretization. Numerical results are presented to support the stability estimates and to show the influence of the SUPG stabilization parameter in a time-dependent domain. In the second part of this work, SUPG stabilization method with non-conservative ALE formulation is proposed. The implicit Euler, Crank-Nicolson and backward difference methods are used for the temporal discretization. At the discrete level in time, the ALE map influences the stability of the corresponding discrete scheme with different time discretizations, and it leads to schemes where conservative and non-conservative formulations are no longer equivalent. The stability of the fully discrete scheme, irrespective of the temporal discretization, is only conditionally stable. It is observed from numerical results that the Crank-Nicolson scheme induces high oscillations in the numerical solution compare to the implicit Euler and the backward difference time discretiza-tions. Moreover, the backward difference scheme is more sensitive to the stabilization parameter k than the other time discretizations. Further, the difference between the solutions obtained with the conservative and non-conservative ALE forms is significant when the deformation of domain is large, whereas it is negligible in domains with small deformation. Finally, the local projection stabilization (LPS) and the higher order dG time stepping scheme are studied for convection dominated problems. The analysis is based on the quadrature formula for approximating the integrals in time. We considered the exact integration in time, which is impractical to implement and the Radau quadrature in time, which can be used in practice. The stability and error estimates are shown for the mathematical basis of considered numerical scheme with both time integration methods. The numerical analysis reveals that the proposed stabilized scheme with exact integration in time is unconditionally stable, whereas Radau quadrature in time is conditionally stable with time-step restriction depending on the ALE map. The theoretical estimates are illustrated with appropriate numerical examples with distinct features. The second order dG(1) time discretization is unconditionally stable while Crank-Nicolson gives the conditional stable estimates only. The convergence order for dG(1) is two which supports the error estimate.

Les écoulements diphasiques jouent un rôle prépondérant dans les moteurs-fusées à ergols liquides cryogéniques, équipant par exemple les lanceurs de la famille Ariane. L'étude expérimentale de tels engins propulsifs étant complexe et onéreuse, disposer d'outils numériques à même de simuler fidèlement leur fonctionnement se révèle être un objectif aussi important qu'ambitieux. La difficulté majeure réside dans le caractère fortement multi-échelles du problème, si bien qu’aucune approche numérique existante n'est capable à elle seule de décrire parfaitement l'ensemble des échelles liquides. Partant de ce constat, les travaux présentés dans cette thèse visent à mettre en place une stratégie de couplage entre des modèles bien adaptés aux différentes topologies d'écoulement diphasique, et ce dans le cadre de la plateforme logicielle multi-physique CEDRE développée par l'ONERA. La démarche adoptée consiste précisément à coupler un modèle à interface diffuse de type ``4 équations'' pour les zones à phases séparées, et un modèle cinétique eulérien pour la phase dispersée, rendant ainsi possible la description de l’atomisation primaire. Par ailleurs, les conditions sévères qui règnent dans les moteurs cryotechniques, où de forts gradients de température, vitesse et densité sont rencontrés, mettent à l'épreuve la robustesse des méthodes numériques. Une nouvelle méthode MUSCL multipente pour maillages non structurés généraux a ainsi été développée, permettant d’améliorer la robustesse et la précision des schémas de discrétisation spatiale. L’ensemble de la stratégie de couplage est finalement appliquée à la simulation du banc Mascotte de l'ONERA pour la combustion cryotechnique. / Two-phase flows play a significant role for the proper functioning of cryogenic liquid-propellant rocketengines, such as those that equip the launchers of the Ariane family. Since the experimental investigationof such propulsion devices is complex and expensive, developing numerical tools able to accuratelysimulate their functioning, is a crucial but nonetheless ambitious objective. The major difficulty is due tothe multiscale nature of the problem, as a result of which there is currently no numerical approach ableto perfectly describe all the liquid scales on its own. Based on this observation the work presented in thisthesis aims at setting up a coupling strategy between models well-adapted to each two-phase flowtopology, in the framework of the ONERA’s multiphysics CEDRE software. The approach adoptedprecisely consists in coupling a 4-equation diffuse interface model for the separated phases and aeulerian kinetic model for the dispersed phase, thus making it possible to describe primary atomization.Besides, the harsh conditions within cryogenic rocket engines, where large temperature, velocity anddensity gradients are encountered, severely challenge the robustness of numerical methods. A newmultislope MUSCL method for general unstructured meshes is thus developed in order to improve therobustness and accuracy of space discretization schemes. The whole coupling strategy is finally appliedto the numerical simulation of the ONERA’s Mascotte test bench for cryogenic combustion research.

In power plants using large utility coal-fired boilers for generation of electricity the coal is pulverised in coal mills and then it has to be pneumatically transported and distributed to a larger number of burners (e.g. 30-40) circumferentially arranged in several rows around the burning chamber of the boiler. Besides the large pipework flow splitting devices are necessary for distribution of an equal amount of pulverised fuel (PF) to each of the burners. So called trifurcators (without inner fittings or guiding vanes) and ''riffle'' type bifurcators are commonly used to split the gas-coal particle flow into two or three pipes/channels with an equal amount of PF mass flow rate in each outflow cross section of the flow splitting device. These PF flow splitting devices are subject of a number of problems. First of all an uneven distribution of PF over the burners of a large utility boiler leads to operational and maintenance problems, increased level of unburned carbon and higher rates of NOX emissions. Maldistribution of fuel between burners caused by non uniform concentration of the PF (particle roping) in pipe and channel bends prior to flow splitting devices leads to uncontrolled differences in the fuel to air ratio between burners. This results in localised regions in the furnace which are fuel rich, where insufficient air causes incomplete combustion of the fuel. Other regions in the furnace become fuel lean, forming high local concentrations of NOX due to the high local concentrations of O2. Otherwise PF maldistribution can impact on power plant maintenance in terms of uneven wear on PF pipework, flow splitters as well as the effects on boiler panels (PF deposition, corrosion, slagging). In order to address these problems in establishing uniform PF distribution over the outlet cross sections of flow splitting devices in the pipework of coal-fired power plants the present paper deals with numerical prediction and analysis of the complex gas and coal particle (PF) flow through trifurcators and ''riffle'' type bifurcators. The numerical investigation is based on a 3-dimensional Eulerian- Lagrangian approach (MISTRAL/PartFlow-3D) developed by Frank et al. The numerical method is capable to predict isothermal, incompressible, steady gas- particle flows in 3-dimensional, geometrically complex flow geometries using boundary fitted, block-structured, numerical grids. Due to the very high numerical effort of the investigated gas-particle flows the numerical approach has been developed with special emphasis on efficient parallel computing on clusters of workstations or other high performance computing architectures. Besides the aerodynamically interaction between the carrier fluid phase and the PF particles the gas-particle flow is mainly influenced by particle-wall interactions with the outer wall boundaries and the inner fittings and guiding vanes of the investigated flow splitting devices. In order to allow accurate quantitative prediction of the motion of the disperse phase the numerical model requires detailed information about the particle-wall collision process. In commonly used physical models of the particle-wall interaction this is the knowledge or experimental prediction of the restitution coefficients (dynamic friction coefficient, coefficient of restitution) for the used combination of particle and wall material, e.g. PF particles on steel. In the present investigation these parameters of the particle-wall interaction model have been obtained from special experiments in two test facilities. Basic experiments to clarify the details of the particle-wall interaction process were made in a test facility with a spherical disk accelerator. This test facility furthermore provides the opportunity to investigate the bouncing process under normal pressure as well as under vacuum conditions, thus excluding aerodynamically influences on the motion of small particles in the near vicinity of solid wall surfaces (especially under small angles of attack). In this experiments spherical glass beads were used as particle material. In a second test facility we have investigated the real impact of non-spherical pulverised fuel particles on a steel/ceramic target. In this experiments PF particles were accelerated by an injector using inert gas like e.g. CO2 or N2 as the carrier phase in order to avoid dust explosion hazards. The obtained data for the particle-wall collision models were compared to those obtained for glass spheres, where bouncing models are proofed to be valid. Furthermore the second test facility was used to obtain particle erosion rates for PF particles on steel targets as a function of impact angles and velocities. The results of experimental investigations has been incorporated into the numerical model. Hereafter the numerical approach MISTRAL/PartFlow-3D has been applied to the PF flow through a ''riffle'' type bifurcator. Using ICEM/CFD-Hexa as grid generator a numerical mesh with approximately 4 million grid cells has been designed for approximation of the complex geometry of the flow splitting device with all its interior fittings and guiding vanes. Based on a predicted gas flow field a large number of PF particles are tracked throughout the flow geometry of the flow-splitter. Besides mean quantities of the particle flow field like e.g. local particle concentrations, mean particle velocities, distribution of mean particle diameter, etc. it is now possible to obtain information about particle erosion on riffle plates and guiding vanes of the flow splitting device. Furthermore the influence of different roping patterns in front of the flow splitter on the uniformness of PF mass flow rate splitting after the bifurcator has been investigated numerically. Results show the efficient operation of the investigated bifurcator in absence of particle roping, this means under conditions of an uniform PF particle concentration distribution in the inflow cross section of the bifurcator. If particle roping occurs and particle concentration differs over the pipe cross section in front of the bifurcator the equal PF particle mass flow rate splitting can be strongly deteriorated in dependence on the location and intensity of the particle rope or particle concentration irregularities. The presented results show the importance of further development of efficient rope splitting devices for applications in coal-fired power plants. Numerical analysis can be used as an efficient tool for their investigation and further optimisation under various operating and flow conditions.

The transport or the separation of solid particles or droplets suspended in a fluid flow is a common task in mechanical and process engineering. To improve machinery and physical processes (e.g. for coal combustion, reduction of NO_x and soot) an optimization of complex phenomena by simulation applying the fundamental conservation equations is required. Fluid-particle flows are characterized by the ratio of density of the two phases gamma=rho_P/rho_F, by the Stokes number St=tau_P/tau_F and by the loading in terms of void and mass fraction. Those numbers (Stokes number, gamma) define the flow regime and which relevant forces are acting on the particle. Dependent on the geometrical configuration the particle-wall interaction might have a heavy impact on the mean flow structure. The occurrence of particle-particle collisions becomes also more and more important with the increase of the local void fraction of the particulate phase. With increase of the particle loading the interaction with the fluid phase can not been neglected and 2-way or even 4-way coupling between the continous and disperse phases has to be taken into account. For dilute to moderate dense particle flows the Euler-Lagrange method is capable to resolve the main flow mechanism. An accurate computation needs unfortunately a high number of numerical particles (1,...,10^7) to get the reliable statistics for the underlying modelling correlations. Due to the fact that a Lagrangian algorithm cannot be vectorized for complex meshes the only way to finish those simulations in a reasonable time is the parallization applying the message passing paradigma. Frank et al. describes the basic ideas for a parallel Eulererian-Lagrangian solver, which uses multigrid for acceleration of the flow equations. The performance figures are quite good, though only steady problems are tackled. The presented paper is aimed to the numerical prediction of time-dependend fluid-particle flows using the simultanous particle tracking approach based on the Eulerian-Lagrangian and the particle-source-in-cell (PSI-Cell) approach. It is shown in the paper that for the unsteady flow prediction efficiency and load balancing of the parallel numerical simulation is an even more pronounced problem in comparison with the steady flow calculations, because the time steps for the time integration along one particle trajectory are very small per one time step of fluid flow integration and so the floating point workload on a single processor node is usualy rather low. Much time is spent for communication and waiting time of the processors, because for cold flow particle convection not very extensive calculations are necessary. One remedy might be a highspeed switch like Myrinet or Dolphin PCI/SCI (500 MByte/s), which could balance the relative high floating point performance of INTEL PIII processors and the weak capacity of the Fast-Ethernet communication network (100 Mbit/s) of the Chemnitz Linux Cluster (CLIC) used for the presented calculations. Corresponding to the discussed examples calculation times and parallel performance will be presented. Another point is the communication of many small packages, which should be summed up to bigger messages, because each message requires a startup time independently of its size. Summarising the potential of such a parallel algorithm, it will be shown that a Beowulf-type cluster computer is a highly competitve alternative to the classical main frame computer for the investigated Eulerian-Lagrangian simultanous particle tracking approach.

Part 1: Theoretic background information - Review of Hooke’s law for linear elastic materials - The strain energy density of linear elastic materials - Hyperelastic material - Material laws for hyperelastic materials - About selecting the material model and performing tests - Implementation of hyperelastic material laws in Mechanica - Defining hyperelastic material parameters in Mechanica - Test set-ups and specimen shapes of the supported material tests - The uniaxial compression test - Stress and strain definitions in the Mechanica LDA analysis Part 2: Application examples - A test specimen subjected to uniaxial loading - A volumetric compression test - A planar test - Influence of the material law Appendix - PTC Simulation Services Introduction - Dictionary Technical English-German / Teil 1: Theoretische Hintergrundinformation - Das Hookesche Gesetz für linear-elastische Werkstoffe - Die Dehnungsenergiedichte für linear-elastische Materialien - Hyperelastisches Material - Materialgesetze für Hyperelastizität - Auswählen des Materialgesetzes und Testdurchführung - Implementierung der hyperelastischen Materialgesetze in Mechanica - Definieren der hyperelastischen Materialparameter in Mechanica - Testaufbauten und Prüfkörper der unterstützten Materialtests - Der einachsige Druckversuch - Spannungs- und Dehnungsdefinition in der Mechanica-Analyse mit großen Verformungen Teil 2: Anwendungsbeispiele - Ein einachsig beanspruchter Prüfkörper - Ein volumetrischer Drucktest - Ein planarer Test - Einfluss des Materialgesetzes Anhang: - Kurzvorstellung der PTC Simulationsdienstleistungen - Wörterbuch technisches Englisch-Deutsch