ANALYTICAL MODELLING OF DRY-JET WET SPINNING

Similar documents
A DELAYED FRACTIONAL MODEL FOR COCOON HEAT-PROOF PROPERTY

High strength high modulus Fibres

NUMERICAL APPROACH TO INTER-FIBER FLOW IN NON-WOVENS WITH SUPER ABSORBENT FIBERS

DYEING BEHAVIOR PREDICTION OF COTTON FABRICS IN SUPERCRITICAL CO 2

Electrospinning of high-molecule PEO solution

MECHANISM OF NANOFIBER CRIMP

DEVELOPMENT OF AN EXPERIMENTAL AUTOMATIC WET- SPINNING MACHINE

NUMERICAL SIMULATION STUDY OF A STABLE JET SHAPE VARIATION IN ELECTROSPINNING. Donghua University, Shanghai , P. R. China

School of Fashion, Zhongyuan University of Technology, Zhengzhou , China

ISO 2076 INTERNATIONAL STANDARD. Textiles Man-made fibres Generic names. Textiles Fibres chimiques Noms génériques. Fourth edition

AN INTEGRAL TRANSFORM APPLIED TO SOLVE THE STEADY HEAT TRANSFER PROBLEM IN THE HALF-PLANE

EFFECT OF CONCENTRATION AND SALT ADDITIVE ON TAYLOR CONE STRUCTURE. Baturalp YALCINKAYA, Fatma YENER, Funda Cengiz-Çallıoğlu, Oldrich JIRSAK

Contents. Foreword by Darrell H. Reneker

NANOCOMPOSITE CATION EXCHANGE MEMBRANE WITH FOULING RESISTANCE AND ENHANCED SALINITY GRADIENT POWER GENERATION FOR REVERSE ELECTRODIALYSIS

Production, characterization and end products of dissolving pulp

Qualitative analysis of aramide polymers by FT-IR spectroscopy

THERMAL PROTECTION OF ELECTRONIC DEVICES WITH THE NYLON6/66-PEG NANOFIBER MEMBRANES

MODELING ON THE BREAKUP OF VISCO-ELASTIC LIQUID FOR EFFERVESCENT ATOMIZATION

M98-P2 (formerly C98-P1) Non-Newtonian Fluid Flow through Fabrics Matthew W. Dunn Philadelphia University

M98-D01 1. A Fundamental Investigation of the Formation and Properties of Electrospun Fibers

The Effect of PVAc Solution Viscosity on Diameter of PVAc Nanofibres Prepared by Technology of Electrospinning

(12) Patent Application Publication (10) Pub. No.: US 2005/ A1

A THEORETICAL ANALYSIS AND CFD SIMULATION ON THE CERAMIC MONOLITH HEAT EXCHANGER

Polymers in Modified Asphalt Robert Q. Kluttz KRATON Polymers

AN ANALYTICAL THERMAL MODEL FOR THREE-DIMENSIONAL INTEGRATED CIRCUITS WITH INTEGRATED MICRO-CHANNEL COOLING

The Preparation and Characterization of Poly(m-phenyleneisophthalamide)

Lecture No. (1) Introduction of Polymers

Polymer Reaction Engineering

MATERIALS SCIENCE POLYMERS

Three-dimensional Multi-recognition Flexible Wearable

Characterization of cellulose with NMR

FLUID MECHANICS PROF. DR. METİN GÜNER COMPILER

Fabrication of SiO 2, Al 2 O 3, and TiO 2 Microcapsules with Hollow Core and Mesoporous Shell Structure

Synthesis and Characterization of Innovative Multilayer, Multi Metal Oxide Thin Films by Modified Silar Deposition Method

EXPERIMENTAL STUDY ON MOISTURE TRANSFER THROUGH FIREFIGHTERS' PROTECTIVE FABRICS IN RADIANT HEAT EXPOSURES

Synthesis of Ultra-long Hollow Chalcogenide Nanofibers

Harmonic Distortion Modeling of Fully-BalancedThird-order Butterworth Gm-C Filter by a Block Diagram Approach

Original Articles. 1. Introduction. You Zhou 1,, Liwei Xue 1, Kai Yi 1, Li Zhang 1, Seung Kon Ryu 2 and Ri Guang Jin 1

Morphology-Selective Synthesis of Cu(NO3)2 2.5H2O. Micro/Nanostructures Achieved by Rational Manipulation

THE EFFECT OF TWO PHASE (AIR-WATER) FLOW CHARACTERISTICS ON MOMENTUM FLUX DUE TO FLOW TURNING ELEMENTS AT ATMOSPHERIC CONDITIONS

A PREDICTIVE MODEL OF THERMAL CONDUCTIVITY OF PLAIN WOVEN FABRICS

SYNTHESIS OF AN AZO DYE revisited (1 or 2 credits)

VALIDATION OF ACCURACY AND STABILITY OF NUMERICAL SIMULATION FOR 2-D HEAT TRANSFER SYSTEM BY AN ENTROPY PRODUCTION APPROACH

CPGAN # 006. The Basics of Filament Stretching Rheometry

Please do not adjust margins. Graphene oxide based moisture-responsive biomimetic film actuators with nacrelike layered structures

Electrospun Fibers in Catalysis

Prototype ring spinning tester with superconducting magnetic bearing system for high productivity

Supplementary information

Ultrasonic studies of aqueous solutions of poly diallyl dimethyl ammonium chloride

DOWNLOAD OR READ : MELT RHEOLOGY AND ITS APPLICATIONS IN THE PLASTICS INDUSTRY ENGINEERING MATERIALS AND PROCESSES PDF EBOOK EPUB MOBI

Hydrothermally Activated Graphene Fiber Fabrics for Textile. Electrodes of Supercapacitors

INVESTIGATION OF THE PROCESSING PARAMETERS OF A 3D WOVEN REINFORCEMENT

Mengying Li.

DESIGN AND PERFORMANCE OF THE CONVERGING-DIVERGING VORTEX FLOWMETER

Chapter 5 MASS AND ENERGY ANALYSIS OF CONTROL VOLUMES

SORPTION PROCESS USING POLYAMIDE NANOFIBRES TO REMOVE DYE FROM SIMULATED WASTEWATER. Jakub WIENER, Sihle NTAKA, P. S. NGCOBO, Roman KNÍŽEK

Experimental Measurements of sink terms for unsaturated flow during liquid composite molding

Supporting Information

The design and construction of 3D rose petal-shape MoS 2. hierarchical nanostructures with structure-sensitive. properties

Highly Stretchable and Transparent Thermistor Based on Self-Healing Double. Network Hydrogel

Contribution of the solid phase polymerization to the molecular weight distribution in acrylonitrile precipitation copolymerization

IN-LINE ULTRASONIC MONITORING OF POLYMER BLENDING IN A TWIN-SCREW EXTRUDER

A SPATIAL STRUCTURAL DERIVATIVE MODEL FOR ULTRASLOW DIFFUSION

Maximal Thermo-geometric Parameter in a Nonlinear Heat Conduction Equation

Supporting Information for Polybenzimidazolium Salts: A New Class of. Anion-Conducting Polymer

Functionalized lignin nanofibres: An alternative platform to valueadded lignin-based materials

Last update: 15 Oct 2007

SYNTHESIS OF AN AZO DYE revisited (1 or 2 credits)

Basic Fluid Mechanics

Fiber spinning and draw resonance: theoretical background

Introduction to Photolithography

Controllable Synthesis of Functional Polyaniline Nanotubes Via A Complex Template Ying WANG, Donghao SUN a and Yanfeng GUO

INVESTIGATION OF HYDRODYNAMIC PROCESSES IN GEOTHERMAL PLANT

IDENTIFICATION OF DEFLUIDIZATION REGION IN A GAS-SOLID FLUIDIZED BED USING A METHOD BASED ON PRESSURE FLUCTUATION MEASUREMENTS

Tuning the Shell Number of Multi-Shelled Metal Oxide. Hollow Fibers for Optimized Lithium Ion Storage

F r a U n h o F E r I n s T I T U T E F o r C h E M I C a l T E C h n o l o g y I C T EnErgEtic MatErials ParticlE technology

FLAME RETARDANCY AND ULTRAVIOLET RESISTANCE OF SILK FABRIC COATED BY GRAPHENE OXIDE

Supporting Information

Mechanical thermal expansion correction design for an ultrasonic flow meter

Supporting Information

AN ADSORPTION ISOTHERM MODEL FOR ADSORPTION PERFORMANCE OF SILVER-LOADED ACTIVATED CARBON

Semiconductor Optoelectronics Prof. M. R. Shenoy Department of Physics Indian Institute of Technology, Delhi. Lecture - 13 Band gap Engineering

THREE-DIMENSIONAL HAUSDORFF DERIVATIVE DIFFUSION MODEL FOR ISOTROPIC/ANISOTROPIC FRACTAL POROUS MEDIA

ENAS 606 : Polymer Physics

Electronic Supplementary Information. Enhanced Photocatalytic/photoelectrocatalytic Activities

Supplementary Material (ESI) for Chemical Communications This journal is (c) The Royal Society of Chemistry 2009

Assessing Technical Feasibility of Supercritical Extraction Processes Utilizing Laboratory Equipment

Experiment 2 - Using Physical Properties to Identify an Unknown Liquid

ELECTROBLEACHING OF COTTON FABRICS IN SODIUM CHLORIDE SOLUTION F. Vitero 1, P. Monllor 1, M. Bonet 1, E. Morallon 2, C. Quijada 1

Influence of the electrospinning parameters on the morphology of composite nanofibers

Lecture 5: Macromolecules, polymers and DNA

Supplementary Materials for

ICSE Board Class IX Chemistry Paper 3 Solution

Scheme 1: Reaction scheme for the synthesis of p(an-co-mma) copolymer

Test Study on Strength and Permeability Properties of Lime-Fly Ash Loess under Freeze-Thaw Cycles

Electrospinning of PVB Solved in Methanol and Isopropanol

By Tae-Seok Lee, Peter Witting, Harper International. Abstract

Particle Characterization of Pharmaceutical Products by Dynamic Image Analysis

SWNT and MWNT from a Polymeric Electrospun Nanofiber Precursor

The morphology of PVDF/1Gra and PVDF/1IL/1Gra was investigated by field emission scanning

Transcription:

THERMA SCIENCE, Year 017, Vol. 1, No. 4, pp. 1807-181 1807 ANAYTICA MODEING OF DRY-JET WET SPINNING by Hong-Yan IU a,b,c*, Zhi-M I d, Yan-Ju YAO b, c and Frank K. KO a National Engeerg aboratory for Modern Silk, College of Textile and Clothg Engeerg, Soochow University, Suzhou, Cha b School of Fashion Technology, Zhongyuan University of Technology, Zhengzhou, Cha c Department of Materials Engeerg, University of British Columbia, Vancouver, B. C., Canada d Rieter Textile Instrument Co., Shanghai, Cha Origal scientific paper https://doi.org/10.98/tsci16011007 This paper troduces an analytical method for the analysis and design of a dry-jet wet spng system. The 1-D mass conservation equation is used, and velocity distribution is assumed to derive a simple relationship among various spng parameters. The effect of spneret mass flow rate, solution density, spneret structure cludg velocity and air-gap length, and drawg velocity on the dry-jet wet spng was simulated usg the proposed analytical model. Theoretical prediction of fiber diameter is obtaed, which depends upon spng conditions, solution properties, and spneret structure. The theoretical results were verified by comparg experimental data with the numerical solution. It was found obviously that the theoretical prediction has comparable accuracy as that by numerical computation. The analytical model can be useful for prelimary design of a spng process for fabrication of fibers with controllable diameter by adjustg parameters spng conditions. Key words: analytical solution, dry-jet wet spng, fiber diameter Introduction Wet spng is a widely used technology polymer science and textile engeerg [1-7]. Dry-jet wet spng technology, a modification of wet spng, combes both the advantages of melt spng and wet spng, which can be used to fabricate high-performance fibers [8-1]. For the dry-jet wet spng process, the process is usually considered to be composed of two separate parts: the elongational flow the air-gap region, and the double diffusion the coagulation bath. The dynamics the coagulation bath can be regarded as the dynamics of wet spng, whereas the elongational flow the air-gap region is usually referred to as a model of the dynamics of melt spng or dry spng. Some of the benefits of dry-jet wet spng are: (a) high speed of spng, (b) high concentration of dope, (c) high degrees of jet-stretch ratios, and (d) control of coagulation ketics by monitorg coagulation bath parameters. Among these benefits, (a) to (c) are derived because of the use of dry-jet and the air-gap, while (d) is derived from the use of wet coagulation. * Correspondg authors, e-mail: phdliuhongyan@yahoo.com

1808 THERMA SCIENCE, Year 017, Vol. 1, No. 4, pp. 1807-181 In this method the polymer is dissolved an appropriate solvent to make the fiber solution. This solution is then extruded under heat and pressure to an air-gap before it enters a coagulation bath. The produced fiber is then washed and dried before it is heat treated and drawn. This is an alternative method to wet spng and is required as spng directg to the bath, for some fibers, creates micro-voids that negatively affect the fiber properties, this is due to the solvent beg drawn out of the liquid too quickly. An ert atmosphere may be required to prevent oxidization some polymers, if so fibers are extruded to a nitrogen atmosphere. This method is often required for high performance fibers with a liquid crystal structure. Due to their structural properties their melt temperature is either the same as, or dangerously close to their decomposition temperature, therefore they must be dissolved an appropriate solvent and extruded this manner. Several research works give valuable sight to dry-jet wet spng. Hauru et al. [1] explored spng stability and discussed the effects of extrusion velocity, draw ratio, spneret aspect ratio, and bath temperature on mechanical properties and orientation of cellulose filaments from Ionic liquid solution usg dry-jet wet spng method. Qian et al. [13] studied the mechanism and characteristics of dry-jet wet spng of acrylic fibers. Tan et al. [14] vestigated the spnability the dry-jet wet spng of polyacrylonitrile (PAN) precursor fiber. Yang et al. [15] produced ultra-fe fibers usg dry-jet wet spng. The effects of the polymer composition, coagulation bath temperature, and draw ratio on the cross-sectional morphology, structure, and tensile properties were reported by Bajaj et al. [16]. From the previous discussion, it is seen that a mathematical model for the spng process is useful for understandg the relationship between the spng conditions and fiber properties [17, 18]. Recently Xia et al. [, 4] suggested a full 3-D model to predict numerically the fiber diameter along the spng le, and important formation was revealed. The numerical results can be used for optimal design of the spng process, however, numerical simulation can not pick out an explicit relationship among spng conditions, spneret structure, and fiber properties. Therefore, an analytical model is very much needed for practical applications. We adopt the mass conservation equation the dry-jet wet spng process and obta a simple analytical equation for the model. The analytical prediction is compared with the experimental data and numerical results by Xia et al. [] to show the validity of our method. Figure 1. Schematic illustration of the dry-jet wet spng process One-dimensional model A 1-D model is widely used the morphology control fiber fabrication process cludg the bubble electrospng [19] and the bubbfil spng [0]. As shown fig. 1 the dry-jet wet spng process is described [4]. The spng solution enters the spneret. A polymer filament is extruded from a spneret to an air-gap region at a mass flow rate, Q. The filament is then elongated by a take-up device with a take-up speed, u out. The filament cools as it passes through the air-gap region with a specific length,, and orientation begs to occur at a certa take-up force. The gap between the

THERMA SCIENCE, Year 017, Vol. 1, No. 4, pp. 1807-181 1809 spneret and the coagulation bath surface, called air-gap, varies with the type of polymer and technology beg used. In acrylic spng, this gap may be as small as a few millimeters, while lyocell spng, it may be up to several centimeters. The jet, as illustrated fig. 1 durg the spng process can be approximately considered as a steady 1-D compressible flow. The mass conservation equation is shown [19, 1]: 1 π 4 D ρ u = Q (1) where D is the jet diameter, ρ the fluid density, u the jet velocity, and Q the flow rate. The jet velocity can be solved by the momentum equation and the energy equation. It is a complex solution process, but it can be solved numerically []. A simple assumption is made this paper to have a quick look to the relationship between spng conditions and fiber diameter. We assume that the jet velocity changes learly along the spng le: uout u uz ( ) = u + z () where u is the velocity of the spng le at the spneret exit, u out the drawg velocity as illustrated fig. 1, and the distance between the spneret exit and the roller. From eq. (1) we can obta: D Q = (3) u 4 Substitute eq. () to eq. (3) we acquire: D ( z) = u After a simple calculation, we have: u + out u z (4) Dz ( ) = u u + Equation (5) is an analytical model for predictg the fiber diameter without any complex numerical simulation or computation. This equation gives an explicit relationship among the spng condition (the flow rate, Q), solution property (density, ρ), spneret structure, u, and drawg velocity, u out. We use the experimental data [] to verify the prediction of eq. (5) by assumg that the density of the spng le remas constant durg spng process. The values of parameters volved eq. (5) are: Q = 0.136 g/m =.06 10 6 kg/s, ρ = 1.056.635 kg/m 3, u out = 0. m/s, and 18 cm. The diameter of the spneret orifice is D 0 = 0. mm, accordg to eq. (1), we have: out u z (5) u 3 D0 ρ 6 4.06 10 = = = 0.06 m/s π π(0. 10 ) 1056.635 (6)

1810 THERMA SCIENCE, Year 017, Vol. 1, No. 4, pp. 1807-181 Submittg all the previous parameters to eq. (5), we can predict the fal fiber diameter: D(18 cm) = = uout u u + z 6 4.06 10 3 = = 0.106 10 m π1056.635 0. (7) which is very closed to the experimental observation, 100 μm [], fig.. Results discussion By substitutg eq. (6) to eq. (5), we can re-write eq. (5) the form: Figure. Comparison of the theoretical prediction with experimental data and numerical solution by Xia et al. [] 1 Dz ( ) = = = D0 uout u 1 u + z uout πd πd 4 0 ρ + z Q 1 πd0 ρ 1+ 0 ρuout where D 0 is the diameter of the spneret orifice. Equation (8) reveals that the fiber size depends learly upon the size of spneret orifice. When D 0 is the range of hundreds of micrometers, nanofibers can be prepared as discussed by iu et al. [], where the spneret nozzle has a diameter of about 78 μm. Substitute z = to eq. (8), the fal fiber diameter can be written the form: z (8) Q D ( ) = (9) 1 π ρuout 4 This equation shows that a higher drawg velocity or a lower flow rate results a smaller fiber diameter. Similar phenomena were also observed the bubble electrospng [19], the bubbfil spng [0], and electrospng [, 3] Conclusions This paper suggests a simple analytical model for predictg the fiber diameter the dry-jet wet spng process. Neither complex numerical simulation nor complex derivation is volved, only simple mass conservation, eq. (1), is volved. The prediction equation, eq. (5), reveals the relationship among let conditions (flow rate), solution property (density), spneret structure (put velocity), spng length, and drawg velocity. Although the analytical model only provides an approximated solution, the calculated results of prediction has almost the same accuracy as numerical one as illustrated fig., this simple analytical model pro-

THERMA SCIENCE, Year 017, Vol. 1, No. 4, pp. 1807-181 1811 vides a helpful guide for prelimary design of a broad range of spng processes for fabrication of fibers with controllable diameter by adjustg parameters volved eq. (5). Acknowledgment The work is supported by Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), National Natural Science Foundation of Cha under Grant No. 113705, Jiangsu Planned Projects for Postdoctoral Research Funds under Grant No. 1401076B, Cha Postdoctoral Science Foundation under Grant Nos. 015M571806 and 016T90495, Cha National Textile and Apparel Council Project under Grant No. 015011, Key Scientific Research Projects of Henan Provce under Grant No. 16A540001, and program by Cha Scholarship Council. The support provided by the NSERC Discovery Grant is also appreciated. References [1] Hauru,. K., et al., Dry Jet-Wet Spng of Strong Cellulose Filaments from Ionic iquid Solution, Cellulose, 1 (014), 6, pp. 4471-4481 [] Xia, X., et al., Dynamic Modelg of Dry-Jet Wet Spng of Cellulose/[Bmim] Cl Solution: Complete Deformation the Air-Gap Region, Cellulose, (015), 3, pp. 1963-1976 [3] Zhang, W., et al., Effect of Epichlorohydr on the Wet Spng of Carrageenan Fibers under Optimal Parameter Conditions, Carbohydrate Polymers, 150 (016), Oct., pp. 3-40 [4] Xia, X., et al., Numerical Investigation of Spneret Geometric Effect on Spng Dynamics of Dry Jet Wet Spng of Cellulose/[Bmim] Cl Solution, Journal of Applied Polymer Science, 133, (016), 38, pp. 4396-4397 [5] Park, S. K., Farris, R. J., Dry-Jet Wet Spng of Aromatic Polyamic Acid Fiber Usg Chemical Imidization, Polymer, 4 (001), 6, pp. 10087-10093 [6] Xia, X., et al., Rheological Behaviors of Cellulose/[Bmim] Cl Solutions Varied with the Dissolvg Process, Journal of Polymer Research, 1 (014), 7, pp. 1-7 [7] Xia, X., et al., Simulation on Contraction Flow of Concentrated Cellulose/1-Butyl-3-Methylimidazolium Chloride Solution through Spneret Orifice, Materials Research Innovations, 18 (014), Supp., pp. S-874-S87-878 [8] Chae, D. W., et al., Physical Properties of yocell Fibers Spun from Isotropic Cellulose Dope Nmmo Monohydrate, Textile Research Journal, 7 (00), 4, pp. 335-340 [9] Choe, E. W., Kim, S. N., Synthesis, Spng, and Fiber Mechanical-Properties of Poly(Para- Phenylenebenzobisoxazole, Macromolecules, 14 (1981), 4, pp. 90-94 [10] Fk, H. P., et al., Structure Formation of Regenerated Cellulose Materials from Nmmo-Solutions, Progress Polymer Science, 6 (001), 9, pp. 1473-154 [11] Kim, D. B., et al., Dry Jet-Wet Spng of Cellulose/N-Methylmorphole N-Oxide Hydrate Solutions and Physical Properties of yocell Fibers, Textile Research Journal, 75 (005), 4, pp. 331-341 [1] Wan, S. X., et al., Acrylic Fibers Processg with Ionic iquid as Solvent, Polymers for Advanced Technologies, 0 (009), 11, pp. 857-86 [13] Qian, B., et al., The Mechanism and Characterstics of Dry Jet Wet Spng of Acrylic Fibers, Advances Polymer Technology, 6 (1986), 4, pp. 509-59 [14] Tan,., et al., Investigatg the Spnability the Dry Jet Wet Spng of Pan Precursor Fiber, Journal of Applied Polymer Science, 110 (008), 4, pp. 1997-000 [15] Yang, W., et al., Poly (M Phenylene Isophthalamide) Ultrafe Fibers from an Ionic iquid Solution by Dry Jet Wet Electrospng, Journal of Macromolecular Science, Part B, 45 (006), 4, pp. 573-579 [16] Bajaj, P., et al., Structure Development Durg Dry-Jet-Wet Spng of Acrylonitrile/Vyl Acids and Acrylonitrile/Methyl Acrylate Copolymers, Journal of Applied Polymer Science, 86 (00), 3, pp. 773-787 [17] Gupta, B., et al., Preparation of Poly (actic Acid) Fiber by Dry Jet Wet Spng. I. Influence of Draw Ratio on Fiber Properties, Journal of Applied Polymer Science, 100 (006),, pp. 139-146 [18] Gupta, B., et al., Preparation of Poly (actic Acid) Fiber by Dry-Jet-Wet Spng. Ii. Effect of Process Parameters on Fiber Properties, Journal of Applied Polymer Science, 101 (006), 6, pp. 3774-3780

181 THERMA SCIENCE, Year 017, Vol. 1, No. 4, pp. 1807-181 [19] iu, H. Y.,Wang, P., A Short Remark on Wan Model for Electrospng and Bubble Electrospng and Its Development, International Journal of Nonlear Sciences and Numerical Simulation, 16 (015), 1, pp. 1- [0] He, C.-H., et al., Bubbfil Spng for Fabrication of Pva Nanofibers, Thermal Science, 19 (015),, pp. 743-746 [1] Dou, H., et al., A Mathematical Model for the Blown Bubble-Spng and Stab-Proof of Nanofibrous Yarn, Thermal Science, 0 (016), 3, pp. 813-817 [] iu, Z., et al., Active Generation of Multiple Jets for Producg Nanofibers with High Quality and High Throughput, Materials & Design, 94 (016), Mar., pp. 496-501 [3] iu, Z., et al., Tunable Surface Morphology of Electrospun Pmma Fiber Usg Bary Solvent, Applied Surface Science, 364 (016), Feb., pp. 516-51 Paper submitted: January 10, 016 Paper revised: May 5, 016 Paper accepted: October 1, 016 017 Society of Thermal Engeers of Serbia. Published by the Vča Institute of Nuclear Sciences, Belgrade, Serbia. This is an open access article distributed under the CC BY-NC-ND 4.0 terms and conditions.