Development of a Reproducible Powder Characterization Method Using a Powder Rheometer

Similar documents
Measuring the flow properties of powders. FT4 Powder Rheometer. freemantechnology

Application of Ring Shear Testing to Optimize Pharmaceutical Formulation and Process Development of Solid Dosage Forms

PTG-NIR Powder Characterisation System (not yet released for selling to end users)

Formulation of Low Dose Medicines - Theory and Practice

Technical brochure DuraLac H TABLETING DIRECT COMPRESSION ANHYDROUS LACTOSE

Electrostatics and cohesion: Cause or effect?

QbD QUANTITATIVE MEASUREMENTS OF CQAS IN SOLID DOSAGE FORM UNIT OPERATIONS

Title: Plan for Approach to the Understanding and Predicting Excipient Properties and Functionality

COMPARISON OF SOME PHYSICAL PARAMETERS OF WHOLE AND SCORED LISINOPRIL AND LISINOPRIL/ HYDROCHLORTHIAZIDE TABLETS

THE INTRODUCTION OF PROCESS ANALYTICAL TECHNOLOGY, USING NEAR INFRARED ANALYSIS, TO A PHARMACEUTICAL BLENDING PROCESS

A Digital Design Approach to Prediction of Powder Flowability

Thermal Analysis Excellence

A calibration method for a new type of rheometer

Performance evaluation of different model mixers by numerical simulation

Interfacial Shear Rheology of Films Formed by Coffee

Evaluating Pigment Dispersion Quality through Dynamic Oscillation Analysis John Meadows, Surface Specialties - UCB

Particle Characterization of Pharmaceutical Products by Dynamic Image Analysis

International Journal of Advanced Chemical Science and Applications (IJACSA)

Investigating Shear-Induced Particle Migration in Fresh Cement Mortars

Method validation for laser diffraction measurements

Mixing of powders. Wim Oostra

STANDARD PROCEDURE: NWSP R2 (15) Polyacrylate Superabsorbent Powders Determination of the Particle Size Distribution by Sieve Fractionation

A STUDY OF COMPRESSION PROCESS AND PROPERTIES OF TABLETS WITH MICROCRYSTALLINE CELLULOSE AND COLLOIDAL SILICON DIOXIDE

PLGA Conformational Behaviour Dependence of the Network Characterized by Rheological Study

A MECHANISTIC STUDY OF FLOW AND COMPRESSION BEHAVIOR OF PHARMACEUTICAL POWDERS

3 Flow properties of bulk solids

Prof. Jiri Zegzulka, Dr. Jan Necas Lucie Jezerska, Jan Vyletelek, Jiri Rozbroj Daniel Gelnar, Martin Žídek et al

Pharmaceutical Polymers for Tablets and Capsules

Setting Attainable and Practical Particle Size Specifications

Chapter 4: Verification of compendial methods

next generation Welcome to the Dispersion options Smarter particle sizing Your sample, expertly processed and prepared Particle size

BIO & PHARMA ANALYTICAL TECHNIQUES. Chapter 5 Particle Size Analysis

TEMPERATURE AND HUMIDITY DEPENDENCE ON STABILITY OF TORQUE MEASURING DEVICES Tassanai Sanponpute and Nittaya Arksonnarong

Analysis of Shear Rate inside a Concrete Truck Mixer

Research Article. Tao Feng, 1 Feng Wang, 1 Rodolfo Pinal, 1 Carl Wassgren, 2 and M. Teresa Carvajal 1,3

Particle Size and Shape Analysis in PHARMACEUTICAL INDUSTRY.

Influence of forced material in roller compactor parameters I.

SIMULATION IN MAGNETIC FIELD ENHANCED CENTRIFUGATION

ActiPix SDI300 Surface Dissolution Imaging System

Determination and Assessment of the Rheological Properties of Pastes for Screen Printing Ceramics

Powder Technology 200 (2010) Contents lists available at ScienceDirect. Powder Technology. journal homepage:

APPLICATION OF COMPACTION EQUATIONS FOR POWDERED PHARMACEUTICAL MATERIALS

Impact factor: 3.958/ICV: 4.10 ISSN:

Powder Technology 194 (2009) Contents lists available at ScienceDirect. Powder Technology. journal homepage:

The mechanism of plateau formation in the fractional quantum Hall effect

ASEAN GUIDELINES FOR VALIDATION OF ANALYTICAL PROCEDURES

Pharmaceutics I صيدالنيات 1. Unit 6

LAMMPS Simulation of a Microgravity Shear Cell 299r Progress Report Taiyo Wilson. Units/Parameters:

Inves4ga4on of Mechanical Property Variability in Lactose Products

AEROSIL Colloidal Silicon Dioxide for Pharmaceuticals

Shear tests on limestone powder CRM-116 with Ring Shear Testers RST-XS and RST-01.pc

Extended Material Characterization MCR

VŠB - Technical University of Ostrava, Faculty of Mechanical Engineering, Institute of Transportation,

Quality by Design and Analytical Methods

Available online at ScienceDirect. Procedia Engineering 136 (2016 )

Shear force and bending moment of beams 2.1 Beams 2.2 Classification of beams 1. Cantilever Beam Built-in encastre' Cantilever

Measuring the flow properties of small powder samples using an avalanche tester

Presenting author: S. Joseph Antony, Co-authors: Saeed Albaraki, Babatunde Arowosola, Neha Agarwal University of Leeds. Dr S Joseph Antony

Numerical Simulation on Concrete Pouring Process of Self-Compacting Concrete-Filled Steel Tube

Torsion/Axial Illustration: 1 (3/30/00)

APPENDIX 1 DESCRIPTION OF HOT WIRE ANEMOMETER

X-Ray Net Weight Control of Pharmaceutical Products

Mixing Process of Binary Polymer Particles in Different Type of Mixers

Determination of Design Space for Oral Pharmaceutical Drugs

Use of Near Infrared Spectroscopy for in- and off-line performance determination of continuous and batch powder mixers: opportunities & challenges

Rheology vs. slump and washout. First encounter with rheology

Table of Contents. Preface...xvii. Part 1. Level

Combined Stresses and Mohr s Circle. General Case of Combined Stresses. General Case of Combined Stresses con t. Two-dimensional stress condition

Beams. Beams are structural members that offer resistance to bending due to applied load

EA-10/14. EA Guidelines on the Calibration of Static Torque Measuring Devices. Publication Reference PURPOSE

Development and Validation of UV Spectrophotometric Estimation of Diclofenac Sodium Bulk and Tablet Dosage form using Area under Curve Method

Terra-Mechanical Simulation Using Distinct Element Method

Biotransport: Principles

On Relationship between PVT and Rheological Measurements of Polymer Melts

CDER Risk Assessment to Evaluate Potential Risks from the Use of Nanomaterials in Drug Products

ESTIMATION OF BINGHAM RHEOLOGICAL PARAMETERS OF SCC FROM SLUMP FLOW MEASUREMENT

Effect of Temperature and Pressure on Rheological Measurements of Cement Slurries

D : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.

DVS Intrinsic Compact and Economical Dynamic Vapor Sorption System

Prediction Quality Attributes For Mixtures From Single Component Data

International Journal of PharmTech Research CODEN (USA): IJPRIF, ISSN: , ISSN(Online): Vol.9, No.7, pp , 2016

Pharmaceutical compounding I Colloidal and Surface-Chemical Aspects of Dosage Forms Dr. rer. nat. Rebaz H. Ali

Haleh Azari, Ph.D. AASHTO Materials Reference Laboratory (AMRL) AASHTO Subcommittee on Materials Meeting August 2007

Chapter 1 General Introduction Instructor: Dr. Mürüde Çelikağ Office : CE Building Room CE230 and GE241

CONTRIBUTION OF THE COARSE AGGREGATE FRACTION TO RHEOLOGY

Impact factor: /ICV: PROCESS VALIDATION OF PARACETAMOL SUSPENSION Sisal Shah*, Dinesh G. Desai, A. K.

International Atomic Energy Agency. Department of Nuclear Sciences and Applications. IAEA Environment Laboratories

Torque and Rotational Equilibrium

4. SHAFTS. A shaft is an element used to transmit power and torque, and it can support

Landslide FE Stability Analysis

Dynamic Analysis for Needle Roller Bearings Under Planetary Motion

Enquiry. Demonstration of Uniformity of Dosage Units using Large Sample Sizes. Proposal for a new general chapter in the European Pharmacopoeia

Key Analytical Issues: Sample Preparation, Interferences and Variability. Tim Shelbourn, Eli Lilly and Company

Flow simulation of fiber reinforced self compacting concrete using Lattice Boltzmann method

Module 5: Theories of Failure

SHEAR STRENGTH OF SOIL

STUDY OF THE APPLICABILTY OF CONTENT UNIFORMITY AND DISSOLUTION VARIATION TEST ON ROPINIROLE HYDROCHLORIDE TABLETS

In-situ measurements of the complex acoustic impedance of materials in vehicle interiors

CHAPTER 8 SUMMARY, CONCLUSIONS AND RECOMMENDATIONS

Transcription:

ANNUAL TRANSACTIONS OF THE NORDIC RHEOLOGY SOCIETY, VOL. 20, 2012 Development of a Reproducible Powder Characterization Method Using a Powder Rheometer Søren Vinter Søgaard 1, Morten Allesø 2, Joergen Garnaes 3, Stefania Baldursdottír 1 and Jukka Rantanen 1 1 Department of Pharmacy, Faculty of Health and Medical Sciences, University of Copenhagen, Universitetsparken 2, DK-2100 Copenhagen, Denmark 2 Product Development and Life Cycle, H. Lundbeck A/S, Ottiliavej 9, DK-2500 Valby, Denmark 3 Danish Fundamental Metrology A/S, Matematiktorvet 307, DK-2500 Kgs. Lyngby, Denmark ABSTRACT In this study, a powder rheometer was used to measure flow characteristics of two pharmaceutical model powders. Precise measurements were obtained for one of the model powders whereas the results were less precise for the second powder. In conclusion, further work is needed to increase the mechanistic understanding of powder rheological measurements. INTRODUCTION Modern pharmaceuticals are primarily manufactured by empirical approaches and batch-wise unit operations. However, regulatory authorities have started to require that future manufacturing schemes should be developed, monitored and controlled based on a mechanistic understanding of how material attributes and process parameters affect the performance of the final product 1,2. Yet, the mechanistic understanding of pharmaceutical powder processing is inadequate and generally relies on a few methods that are not very discriminative, e.g. flow through a funnel or Carr s Compressibility Index as defined in several Pharmacopoeia tests. New tools for studying powders are therefore needed in order to increase the current understanding of underlying principles affecting powder behavior. The Freeman Technology 4 (FT4) bench-top powder rheometer, previously described elsewhere in the literature 3,4, might be such a tool due to the fact that it is designed to characterize powders under various conditions in ways that resemble large-scale production environments. In addition, the powder rheometer is able to measure several parameters related to process performance of powders. These methods include rheological, permeability, compressibility and torsional shear tests. Within the recent years, several studies have been published using a FT4 powder rheometer to investigate and characterize different powder properties 5-7. Though the FT4 rheometer does not provide direct inprocess measurements of powder properties, it is indeed a promising tool for increasing the mechanistic understanding of powder behavior in lab-scale. The knowledge obtained in the laboratory can then be transferred and applied to large-scale processes. In 2011, Dumarey et al. showed how methods provided by the FT4 powder rheometer can be of great value to understand how raw material attributes 239

affect a roll compaction process and thus the final tablet quality 6. Nonetheless, before such studies can be performed it is crucial to establish reproducible characterization methods, which is not a trivial task. Plenty of physical properties of powders, e.g. particle size and shape distributions, as well as environmental conditions in the laboratory, e.g. electrostatic charge, influence the measurements 5,7. These variables should be kept in mind when developing a powder characterization method. A commonly used method in the previous mentioned studies is the stability and variable flow rate method (SVFR method), which is a powder rheological method. The objective of this study is to develop an experimental methodology based on the SVFR method, and subsequently evaluate it in order to obtain a reproducible characterization method. Furthermore, the relative standard deviation of the methodology is estimated using 16 replicates. MATERIALS AND METHODS Materials The microcrystalline cellulose (MCC) was purchased as Avicel PH-200 from FMC BioPolymer, Philidelphia, USA. The anhydrous lactose (AHL) was a SuperTab 21AN quality provided from DMV-Fonterra Excipients, Nörten-Hardenburg, Germany. Sample preparation An identical sample preparation procedure was carried out for both MCC and AHL: 896 ml of bulk sample was poured into an open container and placed in an in-house built humidity control chamber with a fixed relative humidity and temperature at 50±5% and 21±1 C, respecttively. The bulk sample was left under these conditions for two days in order to make the water absorption on the powder material reach the equilibrium. After this the bulk sample was sieved through a 0.71 mm mesh to break down any agglomerates that may have arisen during the acclimatization. The bulk sample was then divided into 32 identical samples of 28 ml using a spinning riffler with eight divisions (Retsch PT100, Retsch GmbH, Haan, Germany). This was performed by dividing the bulk sample into 8 portions of 112 ml. Every portion was then further divided into 8 sub-portions of 14 ml. The sub-portions were added together two by two taking sub-portions placed opposite one another at the spinning riffler with the aim of obtaining a representative sample division. Finally, the samples with a volume of 28 ml were placed in the humidity control chamber under the previous mentioned conditions for two days prior to the measurements. 16 of the 32 samples were randomly selected and measured. The remaining samples were applied for other purposes not included in this study. Method The measurements were carried out in the humidity control chamber with a fixed relative humidity (RH) and temperature at 50±5% and 21±1 C, respectively, in order to reduce the measurement error related to static electricity 7,8. The samples were measured using a modified programme based on the SVFR method, which is a standard method provided by the FT4 powder rheometer (Freeman Technology Ltd., Tewkesbury, UK) The default SVFR method comprises a stability method containing seven test cycles and a variable flow rate method containing four test cycles (Fig. 1). However, in order for the powders to reach a stable state, which is recommended by the manufacturer of the rheometer, the stability method was extended to contain 12 test cycles 9. In that way, the modified SVFR method ended up consisting of 16 test cycles. The modified SVFR method was applied to measure the flow characteristics of each of the 16 replicates of AHL and MCC, respectively. 240

Figure 1. Structure of the default stability and variable flow rate method sequence. In the upper line of the sequence, the powder bed is conditioned followed by a splitting procedure to remove excess powder. The middle line illustrates the stability method consisting of seven pairs of conditioning and test cycles. This part was extended to 12 conditioning and test cycles in the modified SVFR method. The blade tip speed is -100 mm/s during the tests in the middle line. The lower line of the sequence illustrates the structure of the variable flow rate method. The numbers in the brackets are the blade tip speed during the test cycles. C: conditioning cycle, T: test cycle, Split: splitting of vessel to provide a precise volume of powder 12. The foundation for the rheological measurements is a geometrical complex stainless steel blade with a diameter of 23.5 mm, which rotates on a vertical axis while measuring the resistance of the powder to flow dynamically 8. During the stability method, the rotating blade moves into and through a powder bed contained in a 25 ml cylindrical vessel made from borosilicate (height: 51 mm, diameter: 25 mm). The blade operates with a tip speed of -100 mm/s (anti-clockwise) and helix angle of -5 10. During the movement of the blade, the torque (τ) and axial force (normal stress, σ) required to move the blade though the powder is recorded and used to calculate the total energy input required 5,11. Similarly, the force and torque can also be measured when the blade moves upward clockwise. The principle of the measurement for the variable flow rate method is the same as the downward movement of the blade, though in this part of the method the tip speed of the blade is varied (Fig. 1). Prior to all tests, a conditioning cycle is performed to improve the uniformity of the powder bed. This conditioning step consists of the blade moving downward through the powder with a helix angle of +5 after which the blade moves upward with a helix angle of -5. In that way, the powder is intended to be sliced and lifted and thereby prepared for the next measurement 5,13. The following equations define the output parameters from the SVFR method. M split and V split denote the mass and volume of the powder, respectively, after any excess powder has been removed. E test x denotes the energy measured during downward testing of test no. x, while E condition x,up denote the energy recorded from the upward clockwise movement of the blade during conditioning test x 12. Conditioned bulk density, CBD (g/ml): (1) Basic flowability energy, BFE (mj): (2) Specific energy, SE (mj/g):,, (3) Stability index, SI: (4) Flow rate index, FRI: (5) 241

RESULTS Fig. 2A and B show the measurements of the 16 individual replicates for AHL and MCC, respectively, which were obtained during the SVFR method. The measurements are expressed in total energy (mj). prior to the first conditioning and test cycle for AHL and MCC. However, based on the bulk densities alone it cannot be assessed if the powders were homogeneous, i.e. the powders were not segregated. Table 1. Flow characteristics of the powder samples obtained by the stability and variable flow rate method. Microcrystalline cellulose CBD, g/ml 0.375 ± 0.002 (0.5%) BFE (mj) 174 ± 23 (13.2%) SE (mj/g) 5.25 ± 0.65 (12.4%) SI 1.11 ± 0.10 (9.3%) FRI 1.31 ± 0.09 (6.5%) Average ± SD (RSD), n=16. Anhydrous lactose 0.662 ± 0.004 (0.6%) 266 ± 9 (3.6%) 5.11 ± 0.15 (2.9%) 1.17 ± 0.27 (23.8%) 1.25 ± 0.03 (2.3%) Figure 2. Measurements obtained from the stability and variable flow rate method. Test no. 1-12 constitute the stability method, while test no. 13-16 constitute the flow rate method. A: 16 replicates of anhydrous lactose. B: 16 replicates of microcrystalline cellulose. Table 1 contains the flow characteristics and the bulk densities of MCC and AHL obtained by the SVFR method. The results are shown as averages with standard deviations (SD) and relative standard deviations (RSD). DISCUSSION The bulk densities shown in Table 1 prove that a very uniform packing state of both powder types is obtained after the splitting procedure. This means that the packing state of the replicates is very similar The BFE and SE obtained for AHL seems to be within a reasonable precision for powders, i.e. less than 5% RSD. In contrast, the same parameters obtained for MCC have much larger RSDs. The observation is also illustrated in Fig. 2 where the total energy obtained from the initial test numbers of both powders differs between the replicates. Nevertheless, the AHL replicates tend to stabilize at the same level of total energy whereas the MCC replicates reach different levels of total energy in the end of the stability method. This influences the BFE and SE (a measure of flowability and cohesivity, respectively), because these parameters depend on the measurements in the end of the stability method 14,15 (Eq. 2 and 3). The same is the case for SI, which is a measure of how much the total energy for the sample is changing during the stability method (Eq. 4). Because all of the AHL replicates tend to end at the same level of total energy for the last test numbers of the stability method independent of the initial total energy, the SI of AHL has 242

a high RSD (Fig. 2A). On the other hand, the individual measurements of the MCC replicates elucidate that the individual replicates tend to start and end at the same level of total energy thereby obtaining SIs which are both closer to one and more precise than AHL. However, what the SI does not show is that the total energy obtained from the MCC replicates is spread over a large range of total energy (Fig. 2B). Looking at the FRI, which according to the manufacturer of the powder rheometer is a measure of the powder s sensitivity to the flow rate, again shows that AHL obtains a more precise result than MCC. In general, the individual replicates of AHL and MCC shows that the replicates of MCC are wider distributed than AHL. Yet, the total energy in both cases increases when the flow rate, i.e. tip speed, decreases (Fig. 2). The evaluation of the flow characteristics in the previous part shows that the parameters obtained for AHL overall are more precise compared to the parameters obtained for MCC. However, the SI results show the opposite although the SIs for AHL and MCC are fairly similar and in both cases above one. Since, test no. 1-12 are performed in the exact same way, this indicates that the measured total energy for the powders are increasing as a function of repeated testing. This could be caused by several phenomena. In the support documents of the powder rheometer, the manufacturer suggest that de-aeration, agglomeration, segregation as well as moisture uptake and electrostatic charge can cause this observation 9. Since, the moisture level of the powders was fixed two days prior to and also during the measurements this cause can most likely be neglected. Additionally, electrostatic charges as a factor is also not plausible as the amount of electrostatic charge decreases with increaseing relative humidity and is therefore barely present at 50 % RH 7. The observation is therefore probably caused by either deaeration, agglomeration or segregation. In the support documents, it is further mentioned that the BFE (or in general the measured total energy) is dependent on many physical powder properties, e.g. particle size/distribution, shape, density and elasticity 8. However, how these physical properties influence the BFE quantitatively is not stated. Furthermore, it is mentioned that a low BFE in some cases represents a powder with good flow properties while a high BFE represents a powder with poor flow. Yet, it is also stated that the opposite might be the case 8. The manufacturer claims that the reason is that the flow pattern is dependent on the powder s particle size. The shear zone in front of the blade is therefore not of constant size 8. This complicates the comparison of BFEs obtained from different powders with different particle sizes. Therefore, it does not seem feasible to establish a direct correlation between BFE and flowability when comparing a large range of different powders. The same seems to be the case for the other parameter measured in this study. In order for powder rheology to become a useful technique for powder characterrization, the theoretical foundation of the technique needs to be increased. Especially, the geometry of the shear zone and the stresses within this zone needs to be elucidated in order to obtain a quantitative statement of flowability and consolidation during testing 16. Based on the discussion above, methods which are well comprehended and based on bulk solids mechanics, e.g. shear testers, still seems to be a better solution than powder rheology for quantifying and predicting and in-process flowability of powders 16. 243

CONCLUSION In this study, two powders, anhydrous lactose (AHL) and microcrystalline cellulose (MCC), have been characterized using a powder rheometer. While the flow characteristics obtained for AHL in general were precise (RSD < 5%), the opposite was the case for MCC. The underlying principles of the powder mechanics, which are causing the variation of the rheological measurements, are not sufficiently understood. These principles need to be elucidated before powder rheology will be a useful characterization method for estimation of powder flowability. ACKNOWLEDGMENTS H. Lundbeck A/S, Danish Fundamental Metrology A/S and Drug Research Academy, University of Copenhagen are acknowledged for their collaboration and financial support of this work. REFERENCES 1. International conference on harmonisation of technical requirements for registration of pharmaceuticals for human use (ICH) 2005. Pharmaceutical Development Q8(R2)., Geneva. 2. U.S.Food and Drug Administration (FDA) 2004. Guidance for Industry. PAT - A Framework for Innovative Pharmaceutical Development, Manufacturing and Quality Assurance. 3. Freeman RE, Cooke JR, Schneider LCR (2009), "Measuring shear properties and normal stresses generated within a rotational shear cell for consolidated and non-consolidated powders", Powder Technology, 190, 65-69. 4. Freeman R. (2007), "Measuring the flow properties of consolidated, conditioned and aerated powders A comparative study using a powder rheometer and a rotational shear cell", Powder Technology, 174, 25-33. 5. Bharadwaj R, Ketterhagen WR, Hancock BC (2010), "Discrete element simulation study of a Freeman powder rheometer", Chemical Engineering Science, 65, 5747-5756. 6. Dumarey M, Wikström H, Fransson M, Sparén A, Tajarobi P, Josefson M, Trygg J (2011), "Combining experimental design and orthogonal projections to latent structures to study the influence of microcrystalline cellulose properties on roll compaction", International Journal of Pharmaceutics, 416, 110-119. 7. Léonard G, Abatzoglou N (2011), "Lubrication of pharmaceutical powder/wall interfaces and electrostatic effects", Powder Technology, 208, 54-62. 8. Freeman Technology (2008), The Basic Flowability Energy W7030. FT4 Support Documents, Castlemorton Common, Worcestershire (UK): Freeman Technology Ltd. p 1-8. 9. Freeman Technology (2007) Stability Method W7011. FT4 Support Documents, (UK): Freeman Technology Ltd. p 1-8. 10. Freeman Technology (2008) Getting Started W7001. FT4 Support Documents, (UK): Freeman Technology Ltd. p 1-6. 11. Freeman Technology (2008), Dynamic Energy Calculation INF21. FT4 Support Documents, Castlemorton Common, Worcestershire (UK): Freeman Technology Ltd. p 1-2. 12. Freeman Technology (2008) Density W7007. FT4 Support Documents, (UK): Freeman Technology Ltd. p 1-2. 13. Freeman Technology (2006) The Conditioning Process W7006. FT4 Support 244

Documents, Castlemorton Common, Worcestershire (UK): Freeman Technology Ltd. p 1. 14. Freeman Technology (2007), Stability & Variable Flow Rate Method W7013. FT4 Support Documents, Castlemorton Common, Worcestershire (UK): Freeman Technology Ltd. p 1-4. 15. Freeman Technology (2008), Specific Energy W7031. FT4 Support Documents, (UK): Freeman Technology Ltd. p 1-2. 16. Schulze D (2007), Powders and Bulk Solids - Behavior, Characterization, Storage and Flow, 1 ed., Berlin Heidelberg New York: Springer. 245