Esse trabalho teve como objetivo contribuir para o debate sobre a importância das políticas de incentivo à inovação no Brasil. Os resultados esperados do uso que as empresas fizeram dos diferentes tipos de instrumentos sobre os gastos em pesquisa e desenvolvimento (P&D) foram avaliados pelo método de diferenças em diferenças. O método permitiu obter as diferenças de gastos entre empresas beneficiárias de instrumentos e as não-beneficiárias em três períodos consecutivos: 2005 em relação à 2003; 2008 em relação à 2005 e de 2011 em relação à 2008. Ao fazer isso, foi possível identificar se tais diferenças foram positivas e significativas, podendo ser atribuídas às influências dos instrumentos. Os instrumentos utilizados foram: incentivos fiscais, Lei de Informática, financiamentos em parcerias, financiamentos sem parcerias e subvenção. E a utilização dos mesmos pelas empresas teve maior relevância no âmbito de diversos programas de apoio à inovação vigentes no país a partir da retomada das políticas industriais e tecnológicas, nos anos 2000. O estudo concluiu que os efeitos positivos e significativos são limitados à determinados grupos tecnológicos e à poucos instrumentos, em geral, de caráter fiscal. Além disso, esses efeitos positivos surgem em apenas um período, sendo que para cada grupo tecnológico foram efetuadas estimativas para três períodos. Também não houve evidências de que os instrumentos financeiros exerçam efeitos significativos sobre as decisões de gastos em pesquisa e desenvolvimento, apesar da maior ênfase dada aos mesmos no período estudado. Os resultados sugerem fraca influência dos mecanismos de apoio à P&D no Brasil sobre o aumento dos gastos privados, apesar dos avanços recentes. / This study aimed to contribute to the debate about the importance of encouraging innovation policies in Brazil. Understanding the role played by public support instruments to research and development activities in relation to private business spending has become fundamental to achieving this end. The expected results of the use that companies made of different types of instruments on spending were evaluated by the method of difference in differences. This allows you to compare, over a certain period of time, the differences in expenses between the companies that made use of any type of instrument and the companies that did not. Thus, the method yielded differences in spending between the beneficiary and nonbeneficiary firms of those instruments of three consecutive periods: 2005 compared to 2003; 2008 compared to 2005 and 2011 compared to 2008. In doing so, it was possible to identify whether these differences were positive and significant and can be attributed to the influence of the instruments. The instruments used were: tax incentives, Information Technology Law, financing partnerships, financing and grants. And their use by companies had most relevance within various programs to support innovation in force in the country from the resumption of industry and technology policy in the 2000s. The study concluded that the positive and significant effects are limited to certain technological groups and the few instruments in general tax character. In addition, these positive effects arise in one period, and for each technology group estimates were made for three periods. There was also no evidence that the financial instruments carry significant effects on spending decisions on research and development, despite the greater emphasis given to them during the study period. The results suggest weak influence of mechanisms to support R&D in Brazil on the rise in private spending, despite progress in recent years.

L'apparition de CUDA en 2007 a rendu les GPU hautement programmables permettant ainsi aux applications scientifiques et techniques de profiter de leur capacité de calcul élevée. Des solutions ultra-rapides pour la résolution des transferts de chaleur par rayonnement et par conduction sur GPU sont présentées dans ce travail. Tout d'abord, la méthode MACZM pour le calcul des facteurs de transferts radiatifs directs en 3D et en milieu semi-transparent est représentée et validée. Ensuite, une implémentation efficace de la méthode à la base d'algorithmes de géométrie discrète et d'une parallélisation optimisée sur GPU dans CUDA atteignant 300 à 600 fois d'accélération, est présentée. Ceci est suivi par la formulation du NRPA, une version non-récursive de l'algorithme des revêtements pour le calcul des facteurs d'échange radiatifs totaux. La complexité du NRPA est inférieure à celle du PA et sont exécution sur GPU est jusqu'à 750 fois plus rapide que l'exécution du PA sur CPU. D'autre part, une implémentation efficace de la LOD sur GPU est présentée, consistant d'une alternance optimisée des solveurs et schémas de parallélisation et achevant une accélération GPU de 75 à 250 fois. Finalement, toutes les méthodes sont appliquées ensemble pour la résolution des transferts de chaleur en 3D dans un four de réchauffage sidérurgique de brames d'acier. Dans ce but, MACZM est appliquée avec un maillage multi-grille et le NRPA est appliqué au four en le découpant en zones, permettant d'avoir un temps de calcul très rapide une précision élevée. Ceci rend les méthodes utilisées de très grande importance pour la conception de stratégies de contrôle efficaces et précises. / The release of CUDA by NVIDIA in 2007 has tremendously increased GPU programmability, thus allowing scientific and engineering applications to take advantage of the high GPU compute capability. In this work, we present ultra-fast solutions for radiation and diffusion heat transfer on the GPU. First, the Multiple Absorption Coefficient Zonal Method (MACZM) for computing direct radiative exchange factors in 3D semi-transparent media is reviewed and validated. Then, an efficient implementation for MACZM is presented, based on discrete geometry algorithms, and an optimized GPU CUDA parallelization. The CUDA implementation achieves 300 to 600 times speed-up. The Non-recursive Plating Algorithm (NRPA), a non-recursive version of the plating algorithm for computing total exchange factors is then formulated. Due to low-complexity matrix multiplication algorithms, the NRPA has lower complexity than the PA does and it runs up to 750 times faster on the GPU by comparison to the CPU PA. On the other hand, an efficient GPU implementation for the Locally One Dimensional (LOD) finite difference split method for solving heat diffusion is presented, based on an optimiwed alternation between parallelization schemes and equation solvers, achieving accelerations from 75 to 250 times. Finally, all the methods are applied together for solving 3D heat transfer in a steel reheating furnace. A multi-grid approach is applied for MACZM and a zone-by zone computation for the NRPA. As a result, high precision and very fast computation time are achieved, making the methods of high interest for building precise and efficient control units.

Les travaux de cette thèse concerne l’étude d’une méthode élément finis d’ordre spatial élevé que l’on peut assimilé à une extension du schéma de Yee. On parle alors de méthode différences finies d’ordre élevé. Après avoir donné, dans un premier chapitre, un historique non exhaustif des principales méthodes utilisées pour résoudre les équations de Maxwell dans le cadre de problèmes de CEM et montré l’ intérêt de disposer d’un solveur de type "différences finies d’ ordre élevé", nous présentons dans un deuxième chapitre le principe de la méthode. Nous donnons pour cela les caractéristiques du schéma spatial et temporel en précisant les conditions de stabilité de la méthode. En outre, dans une étude purement numérique, nous étudions la convergence du schéma. On se focalise ensuite sur la possibilité d’utiliser des ordres spatiaux variable par cellules dans chaque direction de l’espace. Des comparaisons avec le schéma de Yee et un schéma de Galerkin Discontinu particulier sont ensuite effectuées pour montrer les gains en coûts calcul et mémoire et donc l’intérêt de notre approche par rapport aux deux autres. Dans un troisième chapitre, nous nous intéressons à l’étude de modèles physiques indispensable au traitement d’un problème de CEM. Pour cela, nous nous focalisons particulièrement sur un modèle de fil mince oblique, des modèles de matériaux volumiques et minces et enfin sur la prise en compte de sol parfaitement métallique dans une agression de type onde plane. Chaque modèle est détaillé et validé par comparaison avec des solutions analytiques ou résultant de la littérature, sur des exemples canoniques. Le quatrième chapitre est dédié à une technique d’hybridation entre notre méthode et une approche Galerkin Discontinu en vue de traiter des géométries possédant des courbures. Nous donnons pour cela une stratégie d’hybridation basée sur l’échange de flux qui garantie au niveau continue la conservation d’une énergie. Nous présentons ensuite quelques exemples montrant la validité de notre approche dans une stratégie multi-domaines/multi-méthodes que nous précisons. Enfin le dernier chapitre de cette thèse concerne l’exploitation de notre méthode sur des cas industriels en comparaisons avec d’autres méthodes ou des résultats expérimentaux. / This thesis is about the study of a high spatial finite element method whichcan be assimilated at an extension of the Yee schema. In the next, this method is also called high order finite difference method. In the first chapter, we give a non exhaustive recall of the major methods used to treat EMC problems and we show the necessity to have this kind of schema to simulate efficiently some EMC configurations. In the second chapter, the principle of the numerical method is presented and a stability condition is given. A numerical study analysis of the schema convergence is also done. Next, we show the interest to have the possibility to use local spatial order by cell in each direction of the computational domain. Some canonic examples are given to show the advantages interms of CPU time and memory storage of the method by comparison with Yee’s scheme and DG approach. In the third chapter, we define and validate on several examples,some physical models as thin wire, materials and perfectly metallic ground in presence of a plane wave, to have the possibility to treat EMC problems. The fourth chapter is about a hybridization strategy between our high order FDTD method and a DG schema.We focalize our study on a hybrid method which provides an energy conservation of the continuous problem. A numerical example is given to validate the method. Finally, in the last chapter, we present some simulations on industrial problems to show the possibility of the method to treat realistic EMC problems.

Les viroles en acier microallié 18MND5, destinées aux générateurs de vapeur, présentent une composition hétérogène à plusieurs échelles. Un écart au procédé de fabrication ou une teneur en hydrogène excessive, peuvent conduire à la formation des Décohésions Dues à l’Hydrogène. Ces DDH résultent de la désorption de l’hydrogène interne lors du refroidissement jusqu’à température ambiante. La pression interne n’étant pas mesurables expérimentalement, une modélisation du phénomène est requise. Afin de préciser les mécanismes sous-jacents, il est proposé un scénario de formation de ces défauts s’appuyant conjointement sur une expertise et la modélisation des processus de diffusion-désorption-propagation. Les observations ont révélé une corrélation entre les DDH, les zones ségrégées et les amas de MnS (sites préférentiels d’initiation). Un modèle de diffusion dans un milieu hétérogène a été proposé afin d’évaluer la pression interne associée. La pression maximale excède ainsi 8600 bar en considérant une loi d’Abel-Noble optimisée du gaz réel. Le couplage de ce modèle avec la mécanique de la rupture a permis de quantifier l’évolution des paramètres relatifs à la propagation (pression interne, taille finale, vitesse, …). Un scénario de formation des DDH industriel a ainsi pu être formulé sur la base d’une étude paramétrique. Bien que les simulations préliminaires corroborent le retour d’expérience, le modèle raffiné et la prise en compte du gonflement de la DDH semblent sous-estimer la cinétique. Le caractère multi-fissuré des amas de MnS (homogénéisation des propriétés mécaniques) associé à un critère de rupture à l’échelle locale permettrait d’ajuster ce modèle. / Heat generators are manufactured from ingots of 18MND5 (A508cl3) low alloy steel and present composition heterogeneities at different scales. Under specific conditions (non-respect of guidelines or high initial content of H), Hydrogen Induced Cracks (HIC) may result from diffusion-desorption of internal hydrogen during cooling down to room temperature. Since neither hydrogen redistribution nor its internal pressure within cavities could be measured by experimental techniques, quantitative investigation is based on the modelling of related physical phenomena. A scenario of HIC formation, based on industrial feedback and modelling, has been proposed. A correlation between these defects, segregated areas and clusters of MnS (preferred initiation sites) has been revealed by expertise of HIC. A model of diffusion in heterogeneous alloys has then been proposed to assess the maximal pressure of H2 in such HIC. Simulation has shown that internal pressures above 860MPa are reached by considering an optimized Abel-Noble real gas behavior. The previous model has then been coupled to a failure mechanics procedure to characterize and quantify the crack growth parameters. Based on a parametric study, a scenario of HIC formation during the cooling has been proposed regarding process. Although results from preliminary simulations matched with feedback, the refined model based on the pressure induced elastic deformation of HIC has been developed but provided an underestimated kinetic of crack growth. Consequently, the multi-cracked nature of MnS clusters (homogenization of mechanical properties) and the updated local failure criterion appear to be a viable path to adjust predictions.

Neste trabalho são apresentados os resultados do desenvolvimento e teste de esquemas upwind de alta resolução para o controle da difusão numérica em leis de conservação gerais e problemas em dinâmica dos fluidos. Em particular, são derivados dois novos esquemas: o ALUS (Adaptive Linear Upwind Scheme) e o TOPUS (Third-Order Polynomial Upwind Scheme). Esses esquemas são testados no transporte de escalares, em equações 1D tipo convecção-difusão, em sistemas hiperbólicos 1D, nas equações de Euler 2D da dinâmica dos gases e nas equações de Navier-Stokes incompressíveis 2D/3D. Os esquemas são então associados a uma modelagem algébrica não linear para a simulação de problemas de escoamentos incompressíveis turbulentos 2D com/sem superfícies livres / In this work, results of the development and testing of high-resolution upwind schemes for controlling of the numerical diffusion for general conservation laws and fluid dynamics problems are presented. In particular, two new high-resolution upwind schemes are derived, namely, the ALUS (Adaptive Linear Upwind Scheme) and the TOPUS (Third-Order Polynomial Upwind Scheme). These schemes are tested in scalar transport, 1D convection-diffusion equations, 1D hyperbolic systems, 2D Euler equations of the gas dynamics, and in 2D/3D incompressible Navier-Stokes equations. The schemes are then combined with a nonlinear Reynolds stress algebraic equation model for the simulation of 2D incompressible turbulent flows with/without free surfaces

Les écoulements en eau peu profonde se rencontrent dans de nombreuses situations d’intérêts : écoulements de rivières et dans les lacs, mais aussi dans les mers et océans (courants de marée, tsunami, etc.). Ils sont modélisés par un système d’équations aux dérivées partielles, où les inconnues sont la vitesse de l’écoulement et la hauteur d’eau. On peut supposer que la composante verticale de la vitesse est petite devant les composantes horizontales et que ces dernières sont indépendantes de la profondeur. Le modèle est alors donné par les équations de shallow water (SWEs). Cette thèse se concentre sur la conception d’une nouvelle technique d’interaction de plusieurs grilles imbriquées pour modèle en eau peu profonde en utilisant des méthodes numériques. La première partie de cette thèse comprend, La dérivation complète de ces équations à partir des équations de Navier- Stokes est expliquée. Etudier le développement et l’évaluation des méthodes numériques en utilisant des méthodes de différences finies et plusieurs exemples numériques sont appliqués utilisant la condition initiale du niveau gaussien pour 2DSWEs. Dans la deuxième partie de la thèse, nous sommes intéressés à proposer une nouvelle technique d’interaction de plusieurs grilles imbriquées pour résoudre les modèles océaniques en utilisant quatre choix des opérateurs de restriction avec des résultats de haute précision. Notre travail s’est concentré sur la résolution numérique de SWE par grilles imbriquées. A chaque niveau de résolution, nous avons utilisé une méthode classique de différences finies sur une grille C d’Arakawa, avec un schéma de leapfrog complété par un filtre d’Asselin. Afin de pouvoir affiner les calculs dans les régions perturbées et de les alléger dans les zones calmes, nous avons considéré plusieurs niveaux de résolution en utilisant des grilles imbriquées. Ceci permet d’augmenter considérablement le rapport performance de la méthode, à condition de régler efficacement les interactions (spatiales et temporelles) entre les grilles. Dans la troisième partie de cette thèse, plusieurs exemples numéériques sont testés pour 2DSWE avec imbriqués 3:1 et 5:1. Finalement, la quatrième partie de ce travail, certaines applications de grilles imbriquées pour le modèle tsunami sont présentées. / Most flows in the rivers, seas, and ocean are shallow water flow in which the horizontal length andvelocity scales are much larger than the vertical ones. The mathematical formulation of these flows, so called shallow water equations (SWEs). These equations are a system of hyperbolic partial differentialequations and they are effective for many physical phenomena in the oceans, coastal regions, riversand canals. This thesis focuses on the design of a new two-way interaction technique for multiple nested grids 2DSWEs using the numerical methods. The first part of this thesis includes, proposing several ways to develop the derivation of shallow water model. The complete derivation of this system from Navier-Stokes equations is explained. Studying the development and evaluation of numerical methods by suggesting new spatial and temporal discretization techniques in a standard C-grid using an explicit finite difference method in space and leapfrog with Robert-Asselin filter in time which are effective for modeling in oceanic and atmospheric flows. Several numerical examples for this model using Gaussian level initial condition are implemented in order to validate the efficiency of the proposed method. In the second part of our work, we are interested to propose a new two-way interaction technique for multiple nested grids to solve ocean models using four choices of higher restriction operators (update schemes) for the free surface elevation and velocities with high accuracy results. Our work focused on the numerical resolution of SWEs by nested grids. At each level of resolution, we used explicit finite differences methods on Arakawa C-grid. In order to be able to refine the calculations in troubled regions and move them into quiet areas, we have considered several levels of resolution using nested grids. This makes it possible to considerably increase the performance ratio of the method, provided that the interactions (spatial and temporal) between the grids are effectively controlled. In the third part of this thesis, several numerical examples are tested to show and verify twoway interaction technique for multiple nested grids of shallow water models can works efficiently over different periods of time with nesting 3:1 and 5:1 at multiple levels. Some examples for multiple nested grids of the tsunami model with nesting 5:1 using moving boundary conditions are tested in the fourth part of this work.

Cílem diplomové práce je vytvoření 2D numerického modelu pohybujícího se řezu s proměnnými okrajovými podmínkami skutečné geometrie plynulého odlévání a chlazení předlitku v prostředí MATLAB. Model se zabývá vysoce nelineárními termofyzikálními podmínkami oceli během tuhnutí a chlazení. V práci je simulovaná nejen nelinearita termofyzikálních podmínek, ale také nelinearita při fázové změně. Fázová změna je modelovaná pomocí metody entalpie, metody zdánlivé kapacity a metody teplotního zotavení. Všechny výsledky práce jsou porovnány z více hledisek, jako např. z hlediska přesnosti, rychlosti výpočtu, nebo vhodnosti časového diskretizačního kroku pro nelineární problémy, a paralelizace.

In this thesis, important theories in financial mathematics will be explained and derived. These theories will later be used to value financial derivatives. An analytical formula for valuing European call and put option will be derived and European call options will be valued under the Black-Scholes partial differential equation using three different finite difference methods. The Crank-Nicholson method will then be used to value American call options and solve their corresponding free boundary value problem. The optimal exercise boundary can then be plotted from the solution of the free boundary value problem. The algorithm for valuing American call options will then be further developed to solve the stock loan problem. This will be achieved by exploiting a link that exists between American call options and stock loans. The Crank-Nicholson method will be used to value stock loans and their corresponding free boundary value problem. The optimal exit boundary can then be plotted from the solution of the free boundary value problem. The results that are obtained from the numerical calculations will finally be used to discuss how different parameters affect the valuation of American call options and the valuation of stock loans. In the end of the thesis, conclusions about the effect of the different parameters on the optimal prices will be presented. / I det här kandidatexamensarbetet kommer fundamentala teorier inom finansiell matematik förklaras och härledas. Dessa teorier kommer lägga grunden för värderingen av finansiella derivat i detta arbete. En analytisk formel för att värdera europeiska köp- och säljoptioner kommer att härledas. Dessutom kommer europeiska köpoptioner att värderas numeriskt med tre olika finita differensmetoder. Den finita differensmetoden Crank-Nicholson kommer sedan användas för att värdera amerikanska köpoptioner och lösa det fria gränsvärdesproblemet (free boundary value problem). Den optimala omvandlingsgränsen (Optimal Exercise Boundary) kan därefter härledas från det fria gränsvärdesproblemet. Algoritmen för att värdera amerikanska köpoptioner utökas därefter till att värdera lån med aktier som säkerhet. Detta kan åstadkommas genom att utnyttja ett samband mellan amerikanska köpoptioner med lån där aktier används som säkerhet. Den finita differensmetoden Crank-Nicholson kommer dessutom att användas för att värdera lån med aktier som säkerhet. Den optimala avyttringsgränsen (Optimal Exit Boundary) kan därefter härledas från det fria gränsvärdesproblemet. Resultaten från de numeriska beräkningarna kommer slutligen att användas för att diskutera hur olika parametrar påverkar värderingen av amerikanska köpoptioner, samt värdering av lån med aktier som säkerhet. Avslutningsvis kommer slutsatser om effekterna av dessa parametrar att presenteras.

The aim of this thesis is to calculate the stiffness and damping coefficients for radial hydrodynamic bearings. Cylindrical and lemon hydrodynamic bearings are considered. The solution to this problem mainly depends on the hydrodynamic pressure in the bearing. The numerical solution of the Reynolds equation is used to calculate the pressure. The effect of variable viscosity and density of the lubricant due to temperature changes is considered. The static equilibrium position of the journal centre is also solved. The stiffness and damping coefficients are determined using small amplitude journal motions about the equilibrium position. Three methods for determining these coefficients are presented. The outcome of this thesis is an algorithm for the calculation of stiffness and damping coefficients for cylindrical and lemon bearings. Results for lemon bearings are presented and comparison with the commercial software DynRot BR is made. The benefit of this thesis is the creation of an algorithm for the calculation of journal centre equilibrium position, a new way of incorporating the temperature changes in the viscosity and the density of the lubricant, and the modification of a method for calculating stiffness and damping coefficients based on experimental analogy.

Die vorliegende Arbeit beschäftigt sich mit Festkörperkontinuumsmechanik und Metall- bzw. Kristallplastizität auf verschiedenen Längenskalen. Diesbezüglich besteht die Arbeit aus drei größeren Teilen. Im ersten Teil werden Verformungsvorgänge mit expliziter FEM (Finite-Elemente-Methode) und einem makroskopischen phänomenologischen Modell der Viskoplastizität simuliert. Hierbei wird sich auf das Gleichkanalwinkelpressen (ECAP) eines Metallbarrens und die Stauchung einer sogenannten Crashbox konzentriert. In beiden Fällen gelingt es, die im Experiment bereits beobachtete Lokalisierung der Verformung korrekt wiederzugeben. Da bei den Simulationen die konkrete Mikrostruktur des Materials vernachlässigt wird, werden diese Lokalisierungsphänomene als verformungsinduziert angesehen. Der zweite Teil beschäftigt sich mit der Erweiterung des viskoplastischen Modells, sodass mikroskopische Vorgänge der Gitterdefektstruktur des Materials berücksichtigt werden können. Dazu wird ein Modell des dynamischen Verhaltens von Versetzungspopulationen entwickelt und an das makroskopische viskoplastische Modell gekoppelt. Auf diese Weise können Aspekte der sogenannten Kornfeinung – einem komplexen Strukturbildungsprozess von Versetzungen und anderen Gitterdefekten – erfasst werden. Allerdings kann die für die makroskopischen Eigenschaften entscheidende Bildung von Subkorngrenzen auf diese Weise nicht abgebildet werden. Um dies zu erreichen, wird im dritten Teil der Arbeit eine mesoskopische Theorie der Kristallplastizität mit kontinuierlich verteilten Versetzungen verwendet und weiterentwickelt. Hierbei werden die für eine Subkornbildung wesentlichen Freiheitsgrade hinzugenommen, die Anzahl phänomenologischer Ansätze und zugehöriger Materialparameter aber so klein wie möglich gehalten. Mit dieser Kontinuumsversetzungstheorie (KVT) gelingt es, die Bildung von Subkorngrenzen bei großen plastischen Verformungen eines Kristallits zu verfolgen. Bei den impliziten FEM-Simulationen wird ebenfalls eine Lokalisierung beobachtet, allerdings in Bezug auf die Aktivität der Versetzungen in verschiedenen Gleitebenen. Dementsprechend wird dieses Lokalisierungsphänomen als versetzungsinduziert angesehen. Der Beitrag der vorliegenden Arbeit liegt zum einen in der Aufarbeitung und Gegenüberstellung unterschiedlicher methodischer Herangehensweisen zur Modellierung verformungs- und versetzungsinduzierter Lokalisierungsphänomene. Zum anderen wird eine Analyse und Vereinheitlichung der geometrisch linearen KVT nach Berdichevsky & Le vorgenommen. Wie sich dabei zeigt, verhindern inhärente kinematische Einschränkungen der Theorie die Simulation einer Subkornbildung. Aus diesem Grund wird die konsistente geometrisch nichtlineare KVT von Gurtin aufgegriffen und erweitert. Mit einem daraus abgeleiteten elastisch und plastisch anisotropen Modell der Einkristallviskoplastizität wird der Nachweis erbracht, dass die Subkornbildung damit simuliert werden kann. Darüber hinaus wird eine Aufbereitung und Synthese von Algorithmen zur numerischen Lösung der zugehörigen Feldgleichungen mittels der Methode der finiten Differenzen und der finiten Elemente geliefert. Zudem werden beide Näherungsverfahren in Bezug auf Vor- und Nachteile sowie thermodynamische Konsistenz bei der Anwendung auf Mehrfeldprobleme miteinander verglichen. / The present thesis deals with solid continuum mechanics applied to metal and crystal plasticity on different length scales. In this respect, the work consists of three larger parts. In the first part, deformation processes are simulated with explicit FEM (Finite Element Method) and a macroscopic phenomenological model of viscoplasticity. Here the focus is on the Equal-Channel Angular Pressing (ECAP) of a metal billet and the compression of a so-called crash box. In both cases it is possible to correctly reproduce the localization of the deformation as already observed in the experiment. Since the concrete microstructure of the material is neglected in the simulations, these localization phenomena are regarded as deformation-induced. The second part deals with the extension of the viscoplastic model so that microscopic processes of the lattice defect structure of the material can be considered. A model of the dynamic behavior of dislocation populations is developed and coupled to the macroscopic viscoplastic model. In this way, aspects of the so-called grain refinement – a complex structure formation process of dislocations and other lattice defects – can be captured. However, the formation of subgrain boundaries, which is decisive for the macroscopic properties, cannot be predicted in this way. To achieve this, a mesoscopic theory of crystal plasticity with continuously distributed dislocations is used and further developed in the third part of the thesis. Here, the degrees of freedom essential for subgrain formation are added, while the number of phenomenological approaches and associated material parameters are kept as small as possible. With this continuum dislocation theory it is possible to follow the formation of subgrain boundaries during large plastic deformations of a crystallite. In the implicit FEM simulations, localization is also observed, but with respect to the dislocation activity in different slip planes. Accordingly, this localization phenomenon is considered dislocation-induced. The contribution of the present work lies on the one hand in the review and comparison of different methodical approaches to the modeling of deformation- and dislocation-induced localization phenomena. On the other hand, an analysis and unification of the geometrically linear continuum dislocation theory according to Berdichevsky & Le is carried out. As it turns out, inherent kinematic limitations of the theory prevent the simulation of subgrain formation. For this reason the consistent geometrically non-linear continuum dislocation theory from Gurtin is adopted and extended. With the derived model of elastically and plastically anisotropic single crystal viscoplasticity it is proven that subgrain formation can be simulated. Moreover, a preparation and synthesis of algorithms for the numerical solution of the associated field equations using the method of finite differences and finite elements is provided. In addition, both approximation methods are compared in terms of advantages and disadvantages as well as thermodynamic consistency when applied to multi-field problems.