Mixed finite-element methods for elliptic convection-dominated problems arising in semiconductor physics

In this work we develop and analyze an adaptive numerical scheme for simulating a class of macroscopic semiconductor models. At first the numerical modelling of semiconductors is reviewed in order to classify the Energy-Transport models for semiconductors that are later simulated in 2D. In this class of models the flow of charged particles, that are negatively charged electrons and so-called holes, which are quasi-particles of positive charge, as well as their energy distributions are described by a coupled system of nonlinear partial differential equations. A considerable difficulty in simulating these convection-dominated equations is posed by the nonlinear coupling as well as due to the fact that the local phenomena such as "hot electron effects" are only partially assessable through the given data. The primary variables that are used in the simulations are the particle density and the particle energy density. The user of these simulations is mostly interested in the current flow through parts of the domain boundary - the contacts. The numerical method considered here utilizes mixed finite-elements as trial functions for the discrete solution. The continuous discretization of the normal fluxes is the most important property of this discretization from the users perspective. It will be proven that under certain assumptions on the triangulation the particle density remains positive in the iterative solution algorithm. Connected to this result an a priori error estimate for the discrete solution of linear convection-diffusion equations is derived. The local charge transport phenomena will be resolved by an adaptive algorithm, which is based on a posteriori error estimators. At that stage a comparison of different estimations is performed. Additionally a method to effectively estimate the error in local quantities derived from the solution, so-called "functional outputs", is developed by transferring the dual weighted residual method to mixed finite elements. For a model problem we present how this method can deliver promising results even when standard error estimator fail completely to reduce the error in an iterative mesh refinement process.

Anwendung der geoelektrischen 3D-Tomographie für die Analyse thermisch induzierter Strömungen im Labor

Sowohl in der Natur als auch in der Industrie existieren thermisch induzierte Strömungen. Von Interesse für diese Forschungsarbeit sind dabei die Konvektionen im Erdmantel sowie in den Glasschmelzwannen. Der dort stattfindende Materialtransport resultiert aus Unterschieden in der Dichte, der Temperatur und der chemischen Konzentration innerhalb des konvektierenden Materials. Um das Verständnis für die ablaufenden Prozesse zu verbessern, werden von zahlreichen Forschergruppen numerische Modellierungen durchgeführt. Die Verifikation der dafür verwendeten Algorithmen erfolgt meist über die Analyse von Laborexperimenten. Im Vordergrund dieser Forschungsarbeit steht die Entwicklung einer Methode zur Bestimmung der dreidimensionalen Temperaturverteilung für die Untersuchung von thermisch induzierten Strömungen in einem Versuchsbecken. Eine direkte Temperaturmessung im Inneren des Versuchsmaterials bzw. der Glasschmelze beeinflusst allerdings das Strömungsverhalten. Deshalb wird die geodynamisch störungsfrei arbeitende Impedanztomographie verwendet. Die Grundlage dieser Methode bildet der erweiterte Arrhenius-Zusammenhang zwischen Temperatur und spezifischer elektrischer Leitfähigkeit. Während der Laborexperimente wird ein zähflüssiges Polyethylenglykol-Wasser-Gemisch in einem Becken von unten her erhitzt. Die auf diese Weise generierten Strömungen stellen unter Berücksichtigung der Skalierung ein Analogon sowohl zu dem Erdmantel als auch zu den Schmelzwannen dar. Über mehrere Elektroden, die an den Beckenwänden installiert sind, erfolgen die geoelektrischen Messungen. Nach der sich anschließenden dreidimensionalen Inversion der elektrischen Widerstände liegt das Modell mit der Verteilung der spezifischen elektrischen Leitfähigkeit im Inneren des Versuchsbeckens vor. Diese wird mittels der erweiterten Arrhenius-Formel in eine Temperaturverteilung umgerechnet. Zum Nachweis der Eignung dieser Methode für die nichtinvasive Bestimmung der dreidimensionalen Temperaturverteilung wurden mittels mehrerer Thermoelemente an den Beckenwänden zusätzlich direkte Temperaturmessungen durchgeführt und die Werte miteinander verglichen. Im Wesentlichen sind die Innentemperaturen gut rekonstruierbar, wobei die erreichte Messgenauigkeit von der räumlichen und zeitlichen Auflösung der Gleichstromgeoelektrik abhängt.

Speleogenesis and secondary cave minerals in quartz-sandstone and quartzite environment

The main objective of this research is to improve the comprehension of the processes controlling the formation of caves and karst-like morphologies in quartz-rich lithologies (more than 90% quartz), like quartz-sandstones and metamorphic quartzites. In the scientific community the processes actually most retained to be responsible of these formations are explained in the “Arenisation Theory”. This implies a slow but pervasive dissolution of the quartz grain/mineral boundaries increasing the general porosity until the rock becomes incohesive and can be easily eroded by running waters. The loose sands produced by the weathering processes are then evacuated to the surface through processes of piping due to the infiltration of waters from the fracture network or the bedding planes. To deal with these problems we adopted a multidisciplinary approach through the exploration and the study of several cave systems in different tepuis. The first step was to build a theoretical model of the arenisation process, considering the most recent knowledge about the dissolution kinetics of quartz, the intergranular/grain boundaries diffusion processes, the primary diffusion porosity, in the simplified conditions of an open fracture crossed by a continuous flow of undersatured water. The results of the model were then compared with the world’s widest dataset (more than 150 analyses) of water geochemistry collected till now on the tepui, in superficial and cave settings. All these studies allowed verifying the importance and the effectiveness of the arenisation process that is confirmed to be the main process responsible of the primary formation of these caves and of the karst-like superficial morphologies. The numerical modelling and the field observations allowed evaluating a possible age of the cave systems around 20-30 million of years.

Tetti verdi: Analisi sperimentale e Modellazione numerica

I tetti verdi rappresentano, sempre più frequentemente, una tecnologia idonea alla mitigazione alle problematiche connesse all’ urbanizzazione, tuttavia la conoscenza delle prestazioni dei GR estensivi in clima sub-Mediterraneo è ancora limitata. La presente ricerca è supportata da 15 mesi di analisi sperimentali su due GR situati presso la Scuola di Ingegneria di Bologna. Inizialmente vengono comparate, tra loro e rispetto a una superficie di riferimento (RR), le prestazioni idrologiche ed energetiche dei due GR, caratterizzati da vegetazione a Sedum (SR) e a erbe native perenni (NR). Entrambi riducono i volumi defluiti e le temperature superficiali. Il NR si dimostra migliore del SR sia in campo idrologico che termico, la fisiologia della vegetazione del NR determina l'apertura diurna degli stomi e conseguentemente una maggiore evapotraspirazione (ET). Successivamente si sono studiate la variazioni giornaliere di umidità nel substrato del SR riscontrando che la loro ampiezza è influenzata dalla temperatura, dall’umidità iniziale e dalla fase vegetativa. Queste sono state simulate mediante un modello idrologico basato sull'equazione di bilancio idrico e su due modelli convenzionali per la stima della ET potenziale combinati con una funzione di estrazione dell’ umidità dal suolo. Sono stati proposti dei coefficienti di correzione, ottenuti per calibrazione, per considerare le differenze tra la coltura di riferimento e le colture nei GR durante le fasi di crescita. Infine, con l’ausilio di un modello implementato in SWMM 5.1. 007 utilizzando il modulo Low Impact Development (LID) durante simulazioni in continuo (12 mesi) si sono valutate le prestazioni in termini di ritenzione dei plot SR e RR. Il modello, calibrato e validato, mostra di essere in grado di riprodurre in modo soddisfacente i volumi defluiti dai due plot. Il modello, a seguito di una dettagliata calibrazione, potrebbe supportare Ingegneri e Amministrazioni nella valutazioni dei vantaggi derivanti dall'utilizzo dei GR.

Methoden zur Analyse, Interpretation und Modellierung des elektrischen Widerstandes von archäologischen Steinstrukturen

In der Archäologie werden elektrische Widerstandsmessungen routinemäßig zur Prospektion von Fundstellen eingesetzt. Die Methode ist kostengünstig, leicht anwendbar und liefert in den meisten Fällen zuverlässige und leicht zu interpretierende Ergebnisse. Dennoch kann die Methode die archäologischen Strukturen in manchen Fällen nur teilweise oder gar nicht abbilden, wenn die bodenphysikalischen und bodenchemischen Eigenschaften des Bodens und der archäologischen Strukturen dies nicht zulassen. Der spezifische elektrische Widerstand wird durch Parameter wie Wassergehalt, Bodenstruktur, Bodenskelett, Bodentextur, Salinität und Bodentemperatur beeinflusst. Manche dieser Parameter, wie z.B. der Wassergehalt und die Bodentemperatur, unterliegen einer saisonalen Veränderung. Die vorliegende Arbeit untersucht den spezifischen elektrischen Widerstand von archäologischen Steinstrukturen und evaluiert die Möglichkeit, auf Grundlage von Geländemessungen und Laboranalysen archäologische Strukturen und Böden als numerische Modelle darzustellen. Dazu wurde eine Kombination von verschiedenen bodenkundlichen, geoarchäologischen und geophysikalischen Methoden verwendet. Um archäologische Strukturen und Bodenprofile als numerische Widerstandsmodelle darstellen zu können, werden Informationen zur Geometrie der Strukturen und ihren elektrischen Widerstandswerten benötigt. Dabei ist die Qualität der Hintergrundinformationen entscheidend für die Genauigkeit des Widerstandsmodells. Die Geometrie der Widerstandsmodelle basiert auf den Ergebnissen von Rammkernsondierungen und archäologische Ausgrabungen. Die an der Ausbildung des elektrischen Widerstands beteiligten Parameter wurden durch die Analyse von Bodenproben gemessen und ermöglichen durch Pedotransfer-Funktion, wie die Rhoades-Formel, die Abschätzung des spezifischen elektrischen Widerstandes des Feinbodens. Um den Einfluss des Bodenskeletts auf den spezifischen elektrischen Widerstand von Bodenprofilen und archäologischen Strukturen zu berechnen, kamen die Perkolationstheorie und die Effective Medium Theory zum Einsatz. Die Genauigkeit und eventuelle Limitierungen der Methoden wurden im Labor durch experimentelle Widerstandsmessungen an ungestörten Bodenproben und synthetischen Materialien überprüft. Die saisonale Veränderung des Wassergehalts im Boden wurde durch numerische Modelle mit der Software HYDRUS simuliert. Die hydraulischen Modelle wurden auf Grundlage der ermittelten bodenkundlichen und archäologischen Stratigraphie erstellt und verwenden die Daten von lokalen Wetterstationen als Eingangsparameter. Durch die Kombination der HYDRUS-Ergebnisse mit den Pedotransfer-Funktionen konnte der Einfluss dieser saisonalen Veränderung auf die Prospektionsergebnisse von elektrischen Widerstandsmethoden berechnet werden. Die Ergebnisse der Modellierungsprozesse wurden mit den Geländemessungen verglichen. Die beste Übereinstimmung zwischen Modellergebnissen und den Prospektionsergebnissen konnte für die Fallstudie bei Katzenbach festgestellt werden. Bei dieser wurden die Modelle auf Grundlage von archäologischen Grabungsergebnissen und detaillierten bodenkundlichen Analysen erstellt. Weitere Fallstudien zeigen, dass elektrische Widerstandsmodelle eingesetzt werden können, um den Einfluss von ungünstigen Prospektionsbedingungen auf die Ergebnisse der elektrischen Widerstandsmessungen abzuschätzen. Diese Informationen unterstützen die Planung und Anwendung der Methoden im Gelände und ermöglichen eine effektivere Interpretation der Prospektionsergebnisse. Die präsentierten Modellierungsansätze benötigen eine weitere Verifizierung durch den Vergleich der Modellierungsergebnisse mit detailliertem geophysikalischem Gelände-Monitoring von archäologischen Fundstellen. Zusätzlich könnten elektrische Widerstandsmessungen an künstlichen Mauerstrukturen unter kontrollierten Bedingungen zur Überprüfung der Modellierungsprozesse genutzt werden.

Performance analysis of a miniature turbine generator for intracorporeal energy harvesting

Replacement intervals of implantable medical devices are commonly dictated by battery life. Therefore, intracorporeal energy harvesting has the potential to reduce the number of surgical interventions by extending the life cycle of active devices. Given the accumulated experience with intravascular devices such as stents, heart valves, and cardiac assist devices, the idea to harvest a small fraction of the hydraulic energy available in the cardiovascular circulation is revisited. The aim of this article is to explore the technical feasibility of harvesting 1 mW electric power using a miniature hydrodynamic turbine powered by about 1% of the cardiac output flow in a peripheral artery. To this end, numerical modelling of the fluid mechanics and experimental verification of the overall performance of a 1:1 scale friction turbine are performed in vitro. The numerical flow model is validated for a range of turbine configurations and flow conditions (up to 250 mL/min) in terms of hydromechanic efficiency; up to 15% could be achieved with the nonoptimized configurations of the study. Although this article does not entail the clinical feasibility of intravascular turbines in terms of hemocompatibility and impact on the circulatory system, the numerical model does provide first estimates of the mechanical shear forces relevant to blood trauma and platelet activation. It is concluded that the time-integrated shear stress exposure is significantly lower than in cardiac assist devices due to lower flow velocities and predominantly laminar flow.

Snow shielding factors for cosmogenic nuclide dating inferred from long-term neutron detector monitoring

The depth-dependent attenuation of the secondary cosmic-ray particle flux due to snow cover and its effects on production rates of cosmogenic nuclides constitutes a potential source of uncertainty for studies conducted in regions characterized by frequent seasonal snow burial. Recent experimental and numerical modelling studies have yielded new constraints on the effect of hydrogen-rich media on the production rates of cosmogenic nuclides by low- and high-energy neutrons (<10(-3) MeV and >10(2) MeV, respectively). Here we present long-term neutron-detector monitoring data from a natural setting that we use to quantify the effect of snow cover on the attenuation of fast neutrons (0.1-10 MeV), which are responsible for the production of Ne-21 from Mg and Cl-36 from K. We use data measured between July 2001 and May 2008 at seven stations located throughout the Ecrins-Pelvoux massif (French Western Alps) and its surroundings, at elevations ranging from 200 to 2500 m a.s.l. From the cosmic-ray fluxes recorded during summer, when snow is absent, we infer an apparent attenuation length of 148 g cm(-2) in the atmosphere at a latitude of similar to 45 degrees N and for altitudes ranging from similar to 200 to 2500 m a.s.l. Using snow water-equivalent (SWE) values obtained through snow-coring campaigns that overlap in time the neutron monitoring for five stations, we show that fast neutrons are much more strongly attenuated in snow than predicted by a conventional mass-shielding formulation and the attenuation length estimated in the atmosphere. We suggest that such strong attenuation results from boundary effects at the atmosphere/snow interface induced by the high efficiency of water as a neutron moderator. Finally, we propose an empirical model that allows calculating snow-shielding correction factors as a function of SWE for studies using Ne-21 and Cl-36 analyses in Mg- and K-rich minerals, respectively. This empirical model is of interest for studies with a focus on cosmic-ray exposure dating, particularly if the target rocks are made up of mafic to ultramafic units where seasonal snow-cover is a common phenomenon.

Reconstructed glaciation limits in Eurasia (digitized from map Fig.1)

Based on field investigations in northern Russia and interpretation of offshore seismic data, we have made a preliminary reconstruction of the maximum ice-sheet extent in the Barents and Kara Sea region during the Early/Middle Weichselian and the Late Weichselian. Our investigations indicate that the Barents and Kara ice sheets attained their maximum Weichselian positions in northern Russia prior to 50 000 yr BP, whereas the northeastern flank of the Scandinavian Ice Sheet advanced to a maximum position shortly after 17000 calendar years ago. During the Late Weichselian (25 000-10000 yr BP), much of the Russian Arctic remained ice-free. According to our reconstruction, the extent of the ice sheets in the Barents and Kara Sea region during the Late Weichselian glacial maximum was less than half that of the maximum model which, up to now, has been widely used as a boundary condition for testing and refining General Circulation Models (GCMs). Preliminary numerical-modelling experiments predict Late Weichselian ice sheets which are larger than the ice extent implied for the Kara Sea region from dated geological evidence, suggesting very low precipitation.

Age and chemical composition of apatites from five DSDP and ODP Sites in the Pacific, Atlantic and Indian Oceans

Since studies on deep-sea cores were carried out in the early 1990s it has been known that ambient temperature may have a marked affect on apatite fission track annealing. Due to sluggish annealing kinetics, this effect cannot be quantified by laboratory annealing experiments. The unknown amount of low-temperature annealing remains one of the main uncertainties for extracting thermal histories from fission track data, particularly for samples which experienced slow cooling in shallow crustal levels. To further elucidate these uncertainties, we studied volcanogenic sediments from five deep-sea drill cores, that were exposed to maximum temperatures between ~10° and 70°C over geological time scales of ~15-120 Ma. Mean track lengths (MTL) and etch pit diameters (Dpar) of all samples were measured, and the chemical composition of each grain analyzed for age and track length measurements was determined by electron microprobe analysis. Thermal histories of the sampled sites were independently reconstructed, based on vitrinite reflectance measurements and/or 1D numerical modelling. These reconstructions were used to test the most widely used annealing models for their ability to predict low-temperature annealing. Our results show that long-term exposure to temperatures below the temperature range of the nominal apatite fission track partial annealing zone results in track shortening ranging between 4 and 11%. Both chlorine content and Dpar values explain the downhole annealing patterns equally well. Low chlorine apatite from one drill core revealed a systematic relation between Si-content and Dpar value. The question whether Si-substitution in apatite has direct and systematic effects on annealing properties however, cannot be addressed by our data. For samples, which remained at temperatures <30°C, and which are low in chlorine, the Laslett et al. [Laslett G., Green P., Duddy I. and Gleadow A. (1987) Thermal annealing of fission tracks in apatite. Chem. Geol. 65, 1-13] annealing model predicts MTL up to 0.6 µm longer than those actually measured, whereas for apatites with intermediate to high chlorine content, which experienced temperatures >30°C, the predictions of the Laslett et al. (1987) model agree with the measured MTL data within error levels. With few exceptions, predictions by the Ketcham et al. [Ketcham R., Donelick R. and Carlson W. (1999) Variability of apatite fission-track annealing kinetics. III: Extrapolation to geological time scales. Am. Mineral. 84/9, 1235-1255] annealing model are consistent with the measured data for samples which remained at temperatures below ~30°C. For samples which experienced maximum temperatures between ~30 and 70°C, and which are rich in chlorine, the Ketcham et al. (1999) model overestimates track annealing.

Fly-through movie for the reconstructed folds in the onset region of the Petermann Glacier, North Greenland

The increasing catalogue of high-quality ice-penetrating radar data provides a unique insight in the internal layering architecture of the Greenland ice sheet. The stratigraphy, an indicator of past deformation, highlights irregularities in ice flow and reveals large perturbations without obvious links to bedrock shape. In this work, to establish a new conceptual model for the formation process, we analysed the radar data at the onset of the Petermann Glacier, North Greenland, and created a three-dimensional model of several distinct stratigraphic layers. We demonstrate that the dominant structures are cylindrical folds sub-parallel to the ice flow. By numerical modelling, we show that these folds can be formed by lateral compression of mechanically anisotropic ice, while a general viscosity contrast between layers would not lead to folding for the same boundary conditions. We conclude that the folds primarily form by converging flow as the mechanically anisotropic ice is channelled towards the glacier.

(Table 1) Moisture content of buried snowpack deposit samples from Pearse Valley, McMurdo Dry Valleys

Buried snowpack deposits are found within the McMurdo Dry Valleys of Antarctica, which offers the opportunity to study these layered structures of sand and ice within a polar desert environment. Four discrete buried snowpacks are studied within Pearse Valley, Antarctica, through in situ observations, sample analyses, O-H isotope measurements and numerical modelling of snowpack stability and evolution. The buried snowpack deposits evolve throughout the year and undergo deposition, melt, refreeze, and sublimation. We demonstrate how the deposition and subsequent burial of snow can preserve the snowpacks in the Dry Valleys. The modelled lifetimes of the buried snowpacks are dependent upon subsurface stratigraphy but are typically less than one year if the lag thickness is less than c. 7 cm and snow thickness is less than c. 10 cm, indicating that some of the Antarctic buried snowpacks form annually. Buried snowpacks in the Antarctic polar desert may serve as analogues for similar deposits on Mars and may be applicable to observations of the north polar erg, buried ice at the Mars Phoenix landing site, and observations of buried ice throughout the martian Arctic. Numerical modelling suggests that seasonal snows and subsequent burial are not required to preserve the snow and ice on Mars.

Using ground-penetrating radar, topography and classification of vegetation to model the sediment and active layer thickness in a periglacial lake catchment, Western Greenland, link to shapefiles

The geometries of a catchment constitute the basis for distributed physically based numerical modeling of different geoscientific disciplines. In this paper results from ground-penetrating radar (GPR) measurements, in terms of a 3D model of total sediment thickness and active layer thickness in a periglacial catchment in western Greenland, is presented. Using the topography, thickness and distribution of sediments is calculated. Vegetation classification and GPR measurements are used to scale active layer thickness from local measurements to catchment scale models. Annual maximum active layer thickness varies from 0.3 m in wetlands to 2.0 m in barren areas and areas of exposed bedrock. Maximum sediment thickness is estimated to be 12.3 m in the major valleys of the catchment. A method to correlate surface vegetation with active layer thickness is also presented. By using relatively simple methods, such as probing and vegetation classification, it is possible to upscale local point measurements to catchment scale models, in areas where the upper subsurface is relatively homogenous. The resulting spatial model of active layer thickness can be used in combination with the sediment model as a geometrical input to further studies of subsurface mass-transport and hydrological flow paths in the periglacial catchment through numerical modelling.

(Supplement 3) Radar sounding observations including bedrock and surface elevation and ice thickness in the ASB region, Antarctica

The first Cenozoic ice sheets initiated in Antarctica from the Gamburtsev Subglacial Mountains and other highlands as a result of rapid global cooling ~34 million years ago. In the subsequent 20 million years, at a time of declining atmospheric carbon dioxide concentrations and an evolving Antarctic circumpolar current, sedimentary sequence interpretation and numerical modelling suggest that cyclical periods of ice-sheet expansion to the continental margin, followed by retreat to the subglacial highlands, occurred up to thirty times. These fluctuations were paced by orbital changes and were a major influence on global sea levels. Ice-sheet models show that the nature of such oscillations is critically dependent on the pattern and extent of Antarctic topographic lowlands. Here we show that the basal topography of the Aurora Subglacial Basin of East Antarctica, at present overlain by 2-4.5 km of ice, is characterized by a series of well-defined topographic channels within a mountain block landscape. The identification of this fjord landscape, based on new data from ice-penetrating radar, provides an improved under¬standing of the topography of the Aurora Subglacial Basin and its surroundings, and reveals a complex surface sculpted by a succession of ice-sheet configurations substantially different from today's. At different stages during its fluctuations, the edge of the East Antarctic Ice Sheet lay pinned along the margins of the Aurora Subglacial Basin, the upland boundaries of which are currently above sea level and the deepest parts of which are more than 1 km below sea level. Although the timing of the channel incision remains uncertain, our results suggest that the fjord landscape was carved by at least two ice- flow regimes of different scales and directions, each of which would have over-deepened existing topographic depressions, reversing valley floor slopes.

Dynamic Effects On Culverts for High Speed Trains

Culverts are very common in recent railway lines. Wild life corridors and drainage conducts often fall in this category of partially buried structures. Their dynamic behavior has received far less attention than other structures such as bridges but its large number makes that study an interesting challenge from the point of view of safety and savings. In this paper a complete study of a culvert, including on-site measurements as well as numerical modelling, will be presented. The structure belongs to the high speed railway line linking Segovia and Valladolid, in Spain. The line was opened to traffic in 2004. Its dimensions (3x3m) are the most frequent along the line. Other factors such as reduced overburden (0.6m) and an almost right angle with the track axis make it an interesting example to extract generalized conclusions. On site measurements have been performed in the structure recording the dynamic response at selected points of the structure during the passage of high speed trains at speeds ranging between 200 and 300km/h. The measurements by themselves provide a good insight into the main features of the dynamic behaviour of the structure. A 3D finite element model of the structure, representing its key features was also studied as it allows further understanding of the dynamic response to the train loads . In the paper the discrepancies between predicted and measured vibration levels will be analyzed and some advices on numerical modelling will be proposed

Study of bond stress-slip relationship and radial dilation in prestressed concrete

This paper presents two test procedures for evaluating the bond stress–slip and the slip–radial dilation relationships when the prestressing force is transmitted by releasing the steel (wire or strand) in precast prestressed elements. The bond stress–slip relationship is obtained with short length specimens, to guarantee uniform bond stress, for three depths of the wire indentation (shallow, medium and deep). An analytical model for bond stress–slip relationship is proposed and compared with the experimental results. The model is also compared with the experimental results of other researchers. Since numerical models for studying bond-splitting problems in prestressed concrete require experimental data about dilatancy angle (radial dilation), a test procedure is proposed to evaluate these parameters. The obtained values of the radial dilation are compared with the prior estimated by numerical modelling and good agreement is reached