Quantification of deformation processes in the Torlesse accretionary wedge, New Zealand

Zusammenfassung:In dieser Studie werden Deformationsprozesse im mesozoischen Torlesse Akkretionskeil (Neuseeland) quantifiziert, um Aufschluß über die Dynamik in Akkretionskeilen zu erhalten. Absolute und relative Verformungsmessungen zeigen sowohl im lokalen als auch regionalen Maßstab eine stark heterogene Deformation des Torlesse Keils. Die regionale Deformation wurde mit Hilfe einer Tensordurchschnittsberechnung, unter Benutzung einzelner lokaler Verformungsdaten, als uniaxiale Verkürzung entlang einer subvertikalen, maximalen Verkürzungsachse charakterisiert. Absolute Verformungsmessungen an niedriggradigen Metasandsteinen belegen darüber hinaus durchschnittliche Volumenverluste von ca. 20% SiO2. Volumenveränderungen in tieferkrustalen Aufschlüssen wurden mittels einer geochemischen Massenbilanzanalyse abgeschätzt. Chemische Zusammensetzungen höhergradiger Zonen weichen je nach Grad der Volumenverformung von der Protolitzusammensetzung ab und zeigen somit Verluste von 15% SiO2 an. Da Speicherorte für das gelöste Material nicht bekannt sind, muss angenommen werden, dass das Material aus dem Keil abtransportiert wurde. Die Verformungsergebnisse geben weiterhin Aufschluß über den Grad der Kopplung zwischen Akkretionskeil und subduzierter Platte. Die ermittelten Scherwerte in den Gesteinen liegen deutlich unter den zu erwartenden Scherwerten, die mittels eines einfachen Modells berechnet wurden, das sowohl verschiedene Konvergenzgeschwindigkeiten als auch Exhumierungsraten berücksichtigt. Dies belegt, dass der Torlesse Keil stark von der subduzierten pazifischen Platte entkoppelt war und die Deformation hauptsächlich durch den Fluß der Sedimente in und aus dem Keil bestimmt wurde.

Kinematic models of interseismic deformation from inversion of GPS and InSAR measurements to estimate fault parameters and coupling degree

We have used kinematic models in two Italian regions to reproduce surface interseismic velocities obtained from InSAR and GPS measurements. We have considered a Block modeling, BM, approach to evaluate which fault system is actively accommodating the occurring deformation in both considered areas. We have performed a study for the Umbria-Marche Apennines, obtaining that the tectonic extension observed by GPS measurements is explained by the active contribution of at least two fault systems, one of which is the Alto Tiberina fault, ATF. We have estimated also the interseismic coupling distribution for the ATF using a 3D surface and the result shows an interesting correlation between the microseismicity and the uncoupled fault portions. The second area analyzed concerns the Gargano promontory for which we have used jointly the available InSAR and GPS velocities. Firstly we have attached the two datasets to the same terrestrial reference frame and then using a simple dislocation approach, we have estimated the best fault parameters reproducing the available data, providing a solution corresponding to the Mattinata fault. Subsequently we have considered within a BM analysis both GPS and InSAR datasets in order to evaluate if the Mattinata fault may accommodate the deformation occurring in the central Adriatic due to the relative motion between the North-Adriatic and South-Adriatic plates. We obtain that the deformation occurring in that region should be accommodated by more that one fault system, that is however difficult to detect since the poor coverage of geodetic measurement offshore of the Gargano promontory. Finally we have performed also the estimate of the interseismic coupling distribution for the Mattinata fault, obtaining a shallow coupling pattern. Both of coupling distributions found using the BM approach have been tested by means of resolution checkerboard tests and they demonstrate that the coupling patterns depend on the geodetic data positions.

Shear band evolution in Zr/Hf-based bulk metallic glasses under static and dynamic indentations

Bulk metallic glasses (BMGs) exhibit superior mechanical properties as compared with other conventional materials and have been proposed for numerous engineering and technological applications. Zr/Hf-based BMGs or tungsten reinforced BMG composites are considered as a potential replacement for depleted uranium armor-piercing projectiles because of their ability to form localized shear bands during impact, which has been known to be the dominant plastic deformation mechanism in BMGs. However, in conventional tensile, compressive and bending tests, limited ductility has been observed because of fracture initiation immediately following the shear band formation. To fully investigate shear band characteristics, indentation tests that can confine the deformation in a limited region have been pursued. In this thesis, a detailed investigation of thermal stability and mechanical deformation behavior of Zr/Hf-based BMGs is conducted. First, systematic studies had been implemented to understand the influence of relative compositions of Zr and Hf on thermal stability and mechanical property evolution. Second, shear band evolution under indentations were investigated experimentally and theoretically. Three kinds of indentation studies were conducted on BMGs in the current study. (a) Nano-indentation to determine the mechanical properties as a function of Hf/Zr content. (b) Static Vickers indentation on bonded split specimens to investigate the shear band evolution characteristics beneath the indention. (c) Dynamic Vickers indentation on bonded split specimens to investigate the influence of strain rate. It was found in the present work that gradually replacing Zr by Hf remarkably increases the density and improves the mechanical properties. However, a slight decrease in glass forming ability with increasing Hf content has also been identified through thermodynamic analysis although all the materials in the current study were still found to be amorphous. Many indentation studies have revealed only a few shear bands surrounding the indent on the top surface of the specimen. This small number of shear bands cannot account for the large plastic deformation beneath the indentations. Therefore, a bonded interface technique has been used to observe the slip-steps due to shear band evolution. Vickers indentations were performed along the interface of the bonded split specimen at increasing loads. At small indentation loads, the plastic deformation was primarily accommodated by semi-circular primary shear bands surrounding the indentation. At higher loads, secondary and tertiary shear bands were formed inside this plastic zone. A modified expanding cavity model was then used to predict the plastic zone size characterized by the shear bands and to identify the stress components responsible for the evolution of the various types of shear bands. The applicability of various hardness—yield-strength ( H −σγ ) relationships currently available in the literature for bulk metallic glasses (BMGs) is also investigated. Experimental data generated on ZrHf-based BMGs in the current study and those available elsewhere on other BMG compositions were used to validate the models. A modified expanding-cavity model, employed in earlier work, was extended to propose a new H −σγ relationship. Unlike previous models, the proposed model takes into account not only the indenter geometry and the material properties, but also the pressure sensitivity index of the BMGs. The influence of various model parameters is systematically analyzed. It is shown that there is a good correlation between the model predictions and the experimental data for a wide range of BMG compositions. Under dynamic Vickers indentation, a decrease in indentation hardness at high loading rate was observed compared to static indentation hardness. It was observed that at equivalent loads, dynamic indentations produced more severe deformation features on the loading surface than static indentations. Different from static indentation, two sets of widely spaced semi-circular shear bands with two different curvatures were observed. The observed shear band pattern and the strain rate softening in indentation hardness were rationalized based on the variations in the normal stress on the slip plane, the strain rate of shear and the temperature rise associated with the indentation deformation. Finally, a coupled thermo-mechanical model is proposed that utilizes a momentum diffusion mechanism for the growth and evolution of the final spacing of shear bands. The influence of strain rate, confinement pressure and critical shear displacement on the shear band spacing, temperature rise within the shear band, and the associated variation in flow stress have been captured and analyzed. Consistent with the known pressure sensitive behavior of BMGs, the current model clearly captures the influence of the normal stress in the formation of shear bands. The normal stress not only reduces the time to reach critical shear displacement but also causes a significant temperature rise during the shear band formation. Based on this observation, the variation of shear band spacing in a typical dynamic indentation test has been rationalized. The temperature rise within a shear band can be in excess of 2000K at high strain rate and high confinement pressure conditions. The associated drop in viscosity and flow stress may explain the observed decrease in fracture strength and indentation hardness. The above investigations provide valuable insight into the deformation behavior of BMGs under static and dynamic loading conditions. The shear band patterns observed in the above indentation studies can be helpful to understand and model the deformation features under complex loading scenarios such as the interaction of a penetrator with armor. Future work encompasses (1) extending and modifying the coupled thermo-mechanical model to account for the temperature rise in quasistatic deformation; and (2) expanding this model to account for the microstructural variation-crystallization and free volume migration associated with the deformation.

Properties of granular analogue model materials: A community wide survey

We report the material properties of 26 granular analogue materials used in 14 analogue modelling laboratories. We determined physical characteristics such as bulk density, grain size distribution, and grain shape, and performed ring shear tests to determine friction angles and cohesion, and uniaxial compression tests to evaluate the compaction behaviour. Mean grain size of the materials varied between c. 100 and 400 μm. Analysis of grain shape factors shows that the four different classes of granular materials (14 quartz sands, 5 dyed quartz sands, 4 heavy mineral sands and 3 size fractions of glass beads) can be broadly divided into two groups consisting of 12 angular and 14 rounded materials. Grain shape has an influence on friction angles, with most angular materials having higher internal friction angles (between c. 35° and 40°) than rounded materials, whereas well-rounded glass beads have the lowest internal friction angles (between c. 25° and 30°). We interpret this as an effect of intergranular sliding versus rolling. Most angular materials have also higher basal friction angles (tested for a specific foil) than more rounded materials, suggesting that angular grains scratch and wear the foil. Most materials have an internal cohesion in the order of 20–100 Pa except for well-rounded glass beads, which show a trend towards a quasi-cohesionless (C < 20 Pa) Coulomb-type material. The uniaxial confined compression tests reveal that rounded grains generally show less compaction than angular grains. We interpret this to be related to the initial packing density after sifting, which is higher for rounded grains than for angular grains. Ring-shear test data show that angular grains undergo a longer strain-hardening phase than more rounded materials. This might explain why analogue models consisting of angular grains accommodate deformation in a more distributed manner prior to strain localisation than models consisting of rounded grains.

Contact angle effects on microdroplet deformation using CFD

In this paper we use computational fluid dynamics (CFD) to study the effect of contact angle on droplet shape as it moves through a contraction. A new non-dimensional number is proposed in order to predict situations where the deformed droplet will form a slug in the contraction and thus have the opportunity to interact with the channel wall. It is proposed that droplet flow into a contraction is a useful method to ensure that a droplet will wet a channel surface without a trapped lubrication film, and thus help ensure that a slug will remain attached to the wall downstream of the contraction. We demonstrate that when a droplet is larger than a contraction, capillary and Reynolds numbers, and fluid properties may not be sufficient to fully describe the droplet dynamics through a contraction. We show that, with everything else constant, droplet shape and breakup can be controlled simply by changing the wetting properties of the channel wall. CFD simulations with contact angles ranging from 30 degrees to 150 degrees show that lower contact angles can induce droplet breakup while higher contact angles can form slugs with contact angle dependent shape. Crown Copyright (c) 2005 Published by Elsevier Inc. All rights reserved.

Anisotropic viscous models of large-deformation Mohr-Coulomb failure

We have developed a way to represent Mohr-Coulomb failure within a mantle-convection fluid dynamics code. We use a viscous model of deformation with an orthotropic viscoplasticity (a different viscosity is used for pure shear to that used for simple shear) to define a prefered plane for slip to occur given the local stress field. The simple-shear viscosity and the deformation can then be iterated to ensure that the yield criterion is always satisfied. We again assume the Boussinesq approximation, neglecting any effect of dilatancy on the stress field. An additional criterion is required to ensure that deformation occurs along the plane aligned with maximum shear strain-rate rather than the perpendicular plane, which is formally equivalent in any symmetric formulation. We also allow for strain-weakening of the material. The material can remember both the accumulated failure history and the direction of failure. We have included this capacity in a Lagrangian-integration-point finite element code and show a number of examples of extension and compression of a crustal block with a Mohr-Coulomb failure criterion. The formulation itself is general and applies to 2- and 3-dimensional problems.

Shapes and Dynamics of Blood Cells in Poiseuille and Shear Flows

The dynamics, shape, deformation, and orientation of red blood cells in microcirculation affect the rheology, flow resistance and transport properties of whole blood. This leads to important correlations of cellular and continuum scales. Furthermore, the dynamics of RBCs subject to different flow conditions and vessel geometries is relevant for both fundamental research and biomedical applications (e.g drug delivery). In this thesis, the behaviour of RBCs is investigated for different flow conditions via computer simulations. We use a combination of two mesoscopic particle-based simulation techniques, dissipative particle dynamics and smoothed dissipative particle dynamics. We focus on the microcapillary scale of several μm. At this scale, blood cannot be considered at the continuum but has to be studied at the cellular level. The connection between cellular motion and overall blood rheology will be investigated. Red blood cells are modelled as viscoelastic objects interacting hydrodynamically with a viscous fluid environment. The properties of the membrane, such as resistance against bending or shearing, are set to correspond to experimental values. Furthermore, thermal fluctuations are considered via random forces. Analyses corresponding to light scattering measurements are performed in order to compare to experiments and suggest for which situations this method is suitable. Static light scattering by red blood cells characterises their shape and allows comparison to objects such as spheres or cylinders, whose scattering signals have analytical solutions, in contrast to those of red blood cells. Dynamic light scattering by red blood cells is studied concerning its suitability to detect and analyse motion, deformation and membrane fluctuations. Dynamic light scattering analysis is performed for both diffusing and flowing cells. We find that scattering signals depend on various cell properties, thus allowing to distinguish different cells. The scattering of diffusing cells allows to draw conclusions on their bending rigidity via the effective diffusion coefficient. The scattering of flowing cells allows to draw conclusions on the shear rate via the scattering amplitude correlation. In flow, a RBC shows different shapes and dynamic states, depending on conditions such as confinement, physiological/pathological state and cell age. Here, two essential flow conditions are studied: simple shear flow and tube flow. Simple shear flow as a basic flow condition is part of any more complex flow. The velocity profile is linear and shear stress is homogeneous. In simple shear flow, we find a sequence of different cell shapes by increasing the shear rate. With increasing shear rate, we find rolling cells with cup shapes, trilobe shapes and quadrulobe shapes. This agrees with recent experiments. Furthermore, the impact of the initial orientation on the dynamics is studied. To study crowding and collective effects, systems with higher haematocrit are set up. Tube flow is an idealised model for the flow through cylindric microvessels. Without cell, a parabolic flow profile prevails. A single red blood cell is placed into the tube and subject to a Poiseuille profile. In tube flow, we find different cell shapes and dynamics depending on confinement, shear rate and cell properties. For strong confinements and high shear rates, we find parachute-like shapes. Although not perfectly symmetric, they are adjusted to the flow profile and maintain a stationary shape and orientation. For weak confinements and low shear rates, we find tumbling slippers that rotate and moderately change their shape. For weak confinements and high shear rates, we find tank-treading slippers that oscillate in a limited range of inclination angles and strongly change their shape. For the lowest shear rates, we find cells performing a snaking motion. Due to cell properties and resultant deformations, all shapes differ from hitherto descriptions, such as steady tank-treading or symmetric parachutes. We introduce phase diagrams to identify flow regimes for the different shapes and dynamics. Changing cell properties, the regime borders in the phase diagrams change. In both flow types, both the viscosity contrast and the choice of stress-free shape are important. For in vitro experiments, the solvent viscosity has often been higher than the cytosol viscosity, leading to a different pattern of dynamics, such as steady tank-treading. The stress-free state of a RBC, which is the state at zero shear stress, is still controversial, and computer simulations enable direct comparisons of possible candidates in equivalent flow conditions.

Glulam beams reinforced with FRP externally-bonded: theoretical and experimental evaluation

The glued- laminated lumber (glulam) technique is an efficient process for the rational use of wood. Fiber-reinforced polymer (FRPs) associated with glulam beams provide significant improvements in strength and stiffness and alter the failure mode of these structural elements. In this context, this paper presents guidance for glulam beam production, an experimental analysis of glulam beams made of Pinus caribea var. hondurensis species without and with externally-bonded FRP and theoretical models to evaluate reinforced glulam beams (bending strength and stiffness). Concerning the bending strength of the beams, this paper aims only to analyze the limit state of ultimate strength in compression and tension. A specific disposal was used in order to avoid lateral buckling, once the tested beams have a higher ratio height-to-width. The results indicate the need of production control so as to guarantee a higher efficiency of the glulam beams. The FRP introduced in the tensile section of glulam beams resulted in improvements on their bending strength and stiffness due to the reinforcement thickness increase. During the beams testing, two failure stages were observed. The first was a tensile failure on the sheet positioned under the reinforcement layer, while the second occurred as a result of a preliminary compression yielding on the upper side of the lumber, followed by both a shear failure on the fiber-lumber interface and a tensile failure in wood. The model shows a good correlation between the experimental and estimated results.

Fretting fatigue in overhead conductors: Rig design and failure analysis of a Grosbeak aluminium cable steel reinforced conductor

The performance optimisation of overhead conductors depends on the systematic investigation of the fretting fatigue mechanisms in the conductor/clamping system. As a consequence, a fretting fatigue rig was designed and a limited range of fatigue tests was carried out at the middle high cycle fatigue regime in order to access an exploratory S-N curve for a Grosbeak conductor, which was mounted on a mono-articulated aluminium clamping system. Subsequent to these preliminary fatigue tests, the components of the conductor/clamping system, such as ACSR conductor, upper and lower clamps, bolt and nuts, were subjected to a failure analysis procedure in order to investigate the metallurgical free variables interfering on the fatigue test results, aiming at the optimisation of the testing reproducibility. The results indicated that the rupture of the planar fracture surfaces observed in the external At strands of the conductor tested under lower bending amplitude (0.9 mm) occurred by fatigue cracking (I mm deep), followed by shear overload. The V-type fracture surfaces observed in some At strands of the conductor tested under higher bending amplitude (1.3 mm) were also produced by fatigue cracking (approximately 400 mu m deep), followed by shear overload. Shear overload fracture (45 degrees fracture surface) was also observed on the remaining At wires of the conductor tested under higher bending amplitude (1.3 mm). Additionally, the upper and lower Al-cast clamps presented microstructure-sensitive cracking, which was folowed by particle detachment and formation of abrasive debris on the clamp/conductor tribo-interface, promoting even further the fretting mechanism. The detrimental formation of abrasive debris might be inhibited by the selection of a more suitable class of as-cast At alloy for the production of clamps. Finally, the bolt/nut system showed intense degradation of the carbon steel nut (fabricated in ferritic-pearlitic carbon steel, featuring machined threads with 190 HV), with intense plastic deformation and loss of material. Proper selection of both the bolt and nut materials and the finishing processing might prevent the loss in the clamping pressure during the fretting testing. It is important to control the specification of these components (clamps, bolt and nuts) prior to the start of large scale fretting fatigue testing of the overhead conductors in order to increase the reproducibility of this assessment. (c) 2008 Elsevier Ltd. All rights reserved.

Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

CHARACTERIZATION OF Kiss1 NEURONS USING TRANSGENIC MOUSE MODELS

Humans and mice with loss-of-function mutations of the genes encoding kisspeptins (Kiss1) or kisspeptin receptor (Kiss1r) are infertile due to hypogonadotropic hypogonadism. Within the hypothalamus, Kiss1 mRNA is expressed in the anteroventral periventricular nucleus (AVPV) and the arcuate nucleus (Arc). In order to better study the different populations of kisspeptin cells we generated Kiss1-Cre transgenic mice. We obtained one line with Cre activity specifically within Kiss1 neurons (line J2-4), as assessed by generating mice with Cre-dependent expression of green fluorescent protein or beta-galactosidase. Also, we demonstrated Kiss1 expression in the cerebral cortex and confirmed previous data showing Kiss1 mRNA in the medial nucleus of amygdala and anterodorsal preoptic nucleus. Kiss1 neurons were more concentrated towards the caudal levels of the Arc and higher leptin-responsivity was observed in the most caudal population of Arc Kiss1 neurons. No evidence for direct action of leptin in AVPV Kiss1 neurons was observed. Me lanocortin fibers innervated subsets of Kiss1 neurons of the preoptic area and Arc, and both populations expressed melanocortin receptors type 4 (MC4R). Specifically in the preoptic area, 18-28% of Kiss1 neurons expressed MC4R. In the Arc, 90% of Kiss1 neurons were glutamatergic, 50% of which also were GABAergic. In the AVPV, 20% of Kiss1 neurons were glutamatergic whereas 75% were GABAergic. The differences observed between the Kiss1 neurons in the preoptic area and the Arc likely represent neuronal evidence for their differential roles in metabolism and reproduction. (C) 2011 IBRO. Published by Elsevier Ltd. All rights reserved.

A sensitive vorticity gauge using rotated porphyroblasts, and its application to rocks adjacent to the Alpine Fault, New Zealand

Variable aspect ratio porphyroblasts deformed in non-coaxial flow. and internally containing rotated relicts of an external foliation, can be used to characterise plane strain flow regimes. The distribution obtained by plotting the orientation of the long axis of such grains, classified by aspect ratio, against the orientation of the internal foliation is potentially a sensitive gauge of both the bulk shear strain (as previously suggested) and kinematic vorticity number. We illustrate the method using rotated biotite porphyroblasts in the Alpine Schist: a sequence of mid-crustal rocks that have been ramped to the surface along the Alpine Fault. a major transpressional plate boundary. Results indicate that, at distances greater than or equal to similar to1 km from the fault, the rocks have undergone a combination of irrotational fattening and dextral-oblique, normal-sense shear, with a bulk shear strain of similar to0.6 and kinematic vorticity number of similar to0.2. The vorticity analysis is compatible with estimates of strongly oblate bulk strain of similar to 75% maximum shortening. Dextral-reverse transpressional flow characterises higher strain S-tectonite mylonite within similar to1 km of the Alpine Fault. These relationships provide insight into the kinematics of flow and distribution of strain in the hangingwall of the Alpine Fault and place constraints on numerical mechanical models for the exhumation of these mid-crustal rocks. (C) 2001 Elsevier Science Ltd. All rights reserved.

Stress correlation function evolution in lattice solid elasto-dynamic models of shear and fracture zones and earthquake prediction

It has been argued that power-law time-to-failure fits for cumulative Benioff strain and an evolution in size-frequency statistics in the lead-up to large earthquakes are evidence that the crust behaves as a Critical Point (CP) system. If so, intermediate-term earthquake prediction is possible. However, this hypothesis has not been proven. If the crust does behave as a CP system, stress correlation lengths should grow in the lead-up to large events through the action of small to moderate ruptures and drop sharply once a large event occurs. However this evolution in stress correlation lengths cannot be observed directly. Here we show, using the lattice solid model to describe discontinuous elasto-dynamic systems subjected to shear and compression, that it is for possible correlation lengths to exhibit CP-type evolution. In the case of a granular system subjected to shear, this evolution occurs in the lead-up to the largest event and is accompanied by an increasing rate of moderate-sized events and power-law acceleration of Benioff strain release. In the case of an intact sample system subjected to compression, the evolution occurs only after a mature fracture system has developed. The results support the existence of a physical mechanism for intermediate-term earthquake forecasting and suggest this mechanism is fault-system dependent. This offers an explanation of why accelerating Benioff strain release is not observed prior to all large earthquakes. The results prove the existence of an underlying evolution in discontinuous elasto-dynamic, systems which is capable of providing a basis for forecasting catastrophic failure and earthquakes.

Volatility forecasting with range models: An evaluation of new alternatives to the CARR model

The aim of this paper is to analyze the forecasting ability of the CARR model proposed by Chou (2005) using the S&P 500. We extend the data sample, allowing for the analysis of different stock market circumstances and propose the use of various range estimators in order to analyze their forecasting performance. Our results show that there are two range-based models that outperform the forecasting ability of the GARCH model. The Parkinson model is better for upward trends and volatilities which are higher and lower than the mean while the CARR model is better for downward trends and mean volatilities.

Legal and political regulation of higher education in Portugal: an evolutionary perspective on the current policy issues and the problem of order

The higher education system in Europe is currently under stress and the debates over its reform and future are gaining momentum. Now that, for most countries, we are in a time for change, in the overall society and the whole education system, the legal and political dimensions have gained prominence, which has not been followed by a more integrative approach of the problem of order, its reform and the issue of regulation, beyond the typical static and classical cost-benefit analyses. The two classical approaches for studying (and for designing the policy measures of) the problem of the reform of the higher education system - the cost-benefit analysis and the legal scholarship description - have to be integrated. This is the argument of our paper that the very integration of economic and legal approaches, what Warren Samuels called the legal-economic nexus, is meaningful and necessary, especially if we want to address the problem of order (as formulated by Joseph Spengler) and the overall regulation of the system. On the one hand, and without neglecting the interest and insights gained from the cost-benefit analysis, or other approaches of value for money assessment, we will focus our study on the legal, social and political aspects of the regulation of the higher education system and its reform in Portugal. On the other hand, the economic and financial problems have to be taken into account, but in a more inclusive way with regard to the indirect and other socio-economic costs not contemplated in traditional or standard assessments of policies for the tertiary education sector. In the first section of the paper, we will discuss the theoretical and conceptual underpinning of our analysis, focusing on the evolutionary approach, the role of critical institutions, the legal-economic nexus and the problem of order. All these elements are related to the institutional tradition, from Veblen and Commons to Spengler and Samuels. The second section states the problem of regulation in the higher education system and the issue of policy formulation for tackling the problem. The current situation is clearly one of crisis with the expansion of the cohorts of young students coming to an end and the recurrent scandals in private institutions. In the last decade, after a protracted period of extension or expansion of the system, i. e., the continuous growth of students, universities and other institutions are competing harder to gain students and have seen their financial situation at risk. It seems that we are entering a period of radical uncertainty, higher competition and a new configuration that is slowly building up is the growth in intensity, which means upgrading the quality of the higher learning and getting more involvement in vocational training and life-long learning. With this change, and along with other deep ones in the Portuguese society and economy, the current regulation has shown signs of maladjustment. The third section consists of our conclusions on the current issue of regulation and policy challenge. First, we underline the importance of an evolutionary approach to a process of change that is essentially dynamic. A special attention will be given to the issues related to an evolutionary construe of policy analysis and formulation. Second, the integration of law and economics, through the notion of legal economic nexus, allows us to better define the issues of regulation and the concrete problems that the universities are facing. One aspect is the instability of the political measures regarding the public administration and on which the higher education system depends financially, legally and institutionally, to say the least. A corollary is the lack of clear strategy in the policy reforms. Third, our research criticizes several studies, such as the one made by the OECD in late 2006 for the Ministry of Science, Technology and Higher Education, for being too static and neglecting fundamental aspects of regulation such as the logic of actors, groups and organizations who are major players in the system. Finally, simply changing the legal rules will not necessary per se change the behaviors that the authorities want to change. By this, we mean that it is not only remiss of the policy maker to ignore some of the critical issues of regulation, namely the continuous non-respect by academic management and administrative bodies of universities of the legal rules that were once promulgated. Changing the rules does not change the problem, especially without the necessary debates form the different relevant quarters that make up the higher education system. The issues of social interaction remain as intact. Our treatment of the matter will be organized in the following way. In the first section, the theoretical principles are developed in order to be able to study more adequately the higher education transformation with a modest evolutionary theory and a legal and economic nexus of the interactions of the system and the policy challenges. After describing, in the second section, the recent evolution and current working of the higher education in Portugal, we will analyze the legal framework and the current regulatory practices and problems in light of the theoretical framework adopted. We will end with some conclusions on the current problems of regulation and the policy measures that are discusses in recent years.