A large-time-stepping scheme for balance equations

A well-balanced, large-time-stepping method for conservation laws with source terms is presented. The numerical method is based on a local reformulation of the balance law as a conservation law with a discontinuous flux function, and the approximate solution of this equation by a front tracking method. This yields an unconditionally stable method which is particularly well suited to calculate stationary states. The viability of this approach is demonstrated by several numerical examples. KHK has been supported in part by an Outstanding Young Investigators Award from the Research Council of Norway.     

Implementation of energy-conservation policies: the case of electric heating conversion in Denmark

This article analyses which kinds of public planning, regulation and initiatives are suitable for the implementation of energy-conservation policies. The policies are important elements of the Danish plans for implementing CO₂-reduction targets. The plans are characterised by a move away from energy systems, based on large power stations and large gas and oil extraction systems, to far more decentralised systems based on energy efficiency at the individual places of consumption and renewable energy systems adapted to local conditions. This change demands not only technical modifications, but also large, organisational changes, which will often entail the establishment of completely new organisations. The nature and complexity of the needed technological changes call for public regulation instruments of the same nature, i.e. numerous, differentiated and multi-purpose. In the case of electric-heating conversion, such public regulation so far has managed to achieve CO₂ reductions together with the creation of 150 jobs. In the near future, more than 1000 permanent jobs can be created without negative consequences for the balance-of-payment. The same kind of job creation is possible in most other European countries.     

Creep in phosphorus alloyed copper during power-law breakdown

During the first phase of storage, creep will take place in the copper canisters in the KBS-3 package for nuclear waste. The temperatures are below 100 °C, and the creep is well inside the power-law breakdown regime. Creep models for this situation have been developed. The analysed material is pure copper with about 50 ppm phosphorus. Constitutive equations for creep and other plastic deformation have been set up based on a generalised Norton expression and Kocks-Mecking’s model for the back stress. A model for the minimum creep rate based on fundamental principles for climb and glide has been derived. This model gives the correct order of magnitude for the creep rate in the temperature range from 400 to 20 °C without the use of fitted parameters. The creep exponent varies from 5 to 105 in this interval. The constitutive equations have also been formulated for multiaxial stress states.     

The effect of phosphorus on creep in copper

Pure copper with an addition of about 50 ppm phosphorus is the planned material for the outer part of the waste package for spent nuclear fuel in Sweden. Phosphorus is added to improve the creep ductility but it also strongly increases the creep strength. In the present paper the influence of phosphorus on the strength properties of copper is analysed. Using the Labusch-Nabarro model it is demonstrated that 50 ppm has a negligible influence on the yield strength in accordance with observations. For slow moving dislocations, the interaction energy between the P-atoms and the dislocations gives rise to an agglomeration and a locking. The computed break away stresses are in agreement with the difference in creep stress of copper with and without P-additions.     

Electrokinetic remediation of copper mine tailings: Implementing bipolar electrodes

Important process parameters to optimize in electrokinetic soil remediation are those influencing remediation time and power consumption since these directly affect the cost of a remediation action. This work shows how the electrokinetic remediation (EKR) process could be improved by implementing bipolar electrodes in the porous material. The bipolar electrodes in EKR meant two improvements: (1) a shorter migration pathway for the contaminant, and (2) an increased electrical conductivity in the remediation system. All together the remediation proceeded faster with lower electrical resistance than in similar experiments but without the bipolar electrodes.The new electrokinetic remediation design was tested on copper mine tailings with different applied electric fields, remediation times and pre-treatment. The results showed that the copper removal was increased from 8% (applying 20 V for 8 days in sulphuric acidified tailings) without bipolar electrodes to 42% when bipolar electrodes were implemented.Furthermore, the results showed that in this system sulphuric acid addition prior to remediation was better than citric acid addition. In addition, applying a too strong electric field (even with bipolar electrodes) could cause a severe polarization (e.g. a high electrical resistance) in the remediation system.     

Electroless deposition and silicidation of Ni contacts into p-type Porous Silicon

Porous Silicon (PS) has attracted much attention since the discovery of its photo luminescent behavior. It has also been used for various other applications such as electroluminescent light emitting-diodes (LEDs), photodetectors and solar cells. For such devices, it is important to make good metallic Ohmic contacts to the PS in order to maximize the efficiency. In order to produce buried contacts, barrier layers, Schottky devices, etc. in PS, it is advantageous to deposit metal that covers not only the surface of the porous layer, but also the inside walls and the bottom of the pores. In this work experiments were performed to examine the morphology and properties of electroless deposition of Nickel into p-type PS and subsequent formation of Nickel silicide after heat treatment. Circular PS samples of 6 mm diameter were produced by anodizing p-type Silicon wafers for 15 min and were subsequently plated with Ni using three different plating baths. The pores are on average 20 μm deep and 4 μm wide. Two samples of each type were heat treated in an nitrogen atmosphere for one hour at 400 and 600°C respectively to produce Nickel silicide. Reference samples were made by means of electron beam evaporation of Ni. SEM micrographs show that the best pore coverage was achieved using the Ni plating bath containing hypophosphite. I–V characterization shows that different rectifying and Ohmic contacts can be formed between electroless deposited Ni and PS depending on the conditions of the heat treatment. XRD and EDX characterizations show that both the NiSi and Ni₂Si phases exist in the sample at the same time.      

Experimental investigation of an industrial scale black liquor gasifier. 1. The effect of reactor operation parameters on product gas composition

A novel technology to mitigate the climate changes and improve energy security is Pressurized Entrained flow High Temperature Black Liquor Gasification (PEHT-BLG) in combination with an efficient fuel synthesis using the resulting syngas. In order to optimise the technology for use in a pulp and paper mill based biorefinery, it is of great importance to understand how the operational parameters of the gasifier affect the product gas composition. The present paper is based on experiments where gas samples were withdrawn from the hot part of a 3 MW entrained flow pressurized black liquor gasifier of semi industrial scale using a high temperature gas sampling system. Specifically, the influence of process conditions on product gas composition (CO₂, CO, H₂, CH₄, H₂S, and COS) were examined by systematically varying the operational parameters: system pressure, oxygen to black liquor equivalence ratio, black liquor flow rate to pressure ratio and black liquor pre-heat temperature. Due to the harsh environment inside the gasification reactor, gas sampling is a challenging task. However, for the purpose of the current study, a specially designed high temperature gas sampling system was successfully developed and used. The results, obtained from two separate experimental campaigns, show that all of the investigated operational parameters have a significant influence on the product gas composition and present valuable information about to the process characteristics.     

A numerical model for fluid injection induced seismicity at Soultz-sous-Forêts

During fluid injection experiments at the geothermal site of Soultz-sous-Forêts (France), more than 114,000 induced seismic events with magnitudes between ?2.0 and +2.9 were detected by a local downhole monitoring network. Of these, 35,039 events are sufficiently constrained to be located. Hypocenters align along a sub-vertical, planar structure with the apparent width being dominated by data scattering indicating that seismic activity predominantly occurs along a (pre-existing) larger scale fault structure. For this scenario, we present a numerical model to simulate hydraulic overpressures and induced seismicity during hydraulic injection. The numerical model is based on the physical processes of fluid pressure and stress diffusion with triggering of the induced seismicity being controlled by Coulomb friction.Even in its simplest form of a fault zone without any structural heterogeneity, the numerical model reproduces typical observations at Soultz-sous-Forêts, such as number and magnitude of induced events, hypocenter locations (including the Kaiser effect), occurrence of post-injection seismicity, and the largest magnitude event occurring several days after shut-in.