
Citation: Kevin Caffrey, Mari Chinn, Matthew Veal, Michael Kay. Biomass supply chain management in North Carolina (part 2): biomass feedstock logistical optimization[J]. AIMS Energy, 2016, 4(2): 280-299. doi: 10.3934/energy.2016.2.280
[1] | Bauyrzhan Rakhadilov, Lyaila Bayatanova, Sherzod Kurbanbekov, Ravil Sulyubayev, Nurdaulet Shektibayev, Nurbol Berdimuratov . Investigation on the effect of technological parameters of electrolyte-plasma cementation method on phase structure and mechanical properties of structural steel 20X. AIMS Materials Science, 2023, 10(5): 934-947. doi: 10.3934/matersci.2023050 |
[2] | Andrea Di Schino . Corrosion behavior of new generation super-ferritic stainless steels. AIMS Materials Science, 2019, 6(5): 646-656. doi: 10.3934/matersci.2019.5.646 |
[3] | Julien G. Mahy, Stéphanie D. Lambert, Jérémy Geens, Alain Daniel, David Wicky, Catherine Archambeau, Benoît Heinrichs . Large scale production of photocatalytic TiO2 coating for volatile organic compound (VOC) air remediation. AIMS Materials Science, 2018, 5(5): 945-956. doi: 10.3934/matersci.2018.5.945 |
[4] | Ojo Sunday Isaac Fayomi, Adedamola Sode, Itopa Godwin Akande, Abimbola Patricia Idowu Popoola, Oluranti Agboola . Improving the structural properties and corrosion behaviour of electroless deposited Ni-P-Zn coatings on mild steel for advanced processes. AIMS Materials Science, 2020, 7(4): 441-452. doi: 10.3934/matersci.2020.4.441 |
[5] | Bauyrzhan Rakhadilov, Rinat Kussainov, Aisulu Kalitova, Zarina Satbayeva, Aibek Shynarbek . The impact of technological parameters of electrolytic-plasma treatment on the changes in the mechano-tribological properties of steel 45. AIMS Materials Science, 2024, 11(4): 666-683. doi: 10.3934/matersci.2024034 |
[6] | Mochamad Asrofi, Mochammad Agus Vian Hidayatulloh, Gaguk Jatisukamto, Hary Sutjahjono, Rahma Rei Sakura . The effect of temperature and volume fraction of mahoni (Swietenia mahogani) wood charcoal on SS400 steel using pack carburizing method: Study of hardness and microstructure characteristics. AIMS Materials Science, 2020, 7(3): 354-363. doi: 10.3934/matersci.2020.3.354 |
[7] | ME Makhatha . Effect of titanium addition on sub-structural characteristics of low carbon copper bearing steel in hot rolling. AIMS Materials Science, 2022, 9(4): 604-616. doi: 10.3934/matersci.2022036 |
[8] | Xuan Luc Le, Nguyen Dang Phu, Nguyen Xuan Duong . Enhancement of ferroelectricity in perovskite BaTiO3 epitaxial thin films by sulfurization. AIMS Materials Science, 2024, 11(4): 802-814. doi: 10.3934/matersci.2024039 |
[9] | A. A. Ayoola, N. Auta-Joshua, B. M. Durodola, O. J. Omodara, E. A. Oyeniyi . Combating A36 mild steel corrosion in 1 M H2SO4 medium using watermelon seed oil inhibitor. AIMS Materials Science, 2021, 8(1): 130-143. doi: 10.3934/matersci.2021009 |
[10] | Yu Shiratsuchi, Yoshinori Kotani, Saori Yoshida, Yasunori Yoshikawa, Kentaro Toyoki, Atsushi Kobane, Ryoichi Nakatani, Tetsuya Nakamura . Direct observations of ferromagnetic and antiferromagnetic domains in Pt/Co/Cr2O3/Pt perpendicular exchange biased film. AIMS Materials Science, 2015, 2(4): 484-496. doi: 10.3934/matersci.2015.4.484 |
Oil and gas fields are very severe environments for metals, and high corrosion resistance and strength steels are generally used. Stress corrosion cracking (SCC), sulfide stress cracking (SSC) and hydrogen embrittlement (HE) are principal corrosion problems of steels used in oil and gas environments. The corrosion behavior of steels has been examined in simulated oil and gas environments [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16]. Kimura et al. reported that the modified 13 mass% Cr stainless steel showed excellent resistance to CO2 corrosion and resistance to SSC [2]. Sunaba et al. investigated the effects of chloride ion concentration on the passive and corrosion product films of modified martensitic stainless steels at high temperatures in CO2 environments [7]. They reported that molybdenum improved stability of the corrosion product films in high temperature CO2 and the presence of chloride ions. Hashizume et al. introduced a solution flow type micro-droplet cell and a glow discharge optical emission spectrometer (GDS) to investigate the electrochemical behavior of low carbon 13 mass% Cr welded joints and the effect of a post welded heat treatment (PWHT) in oil and gas environments [17]. They reported that the PWHT eliminated the Cr depleted layer formed during welding, and reduced the SCC susceptibility.
Hashizume et al. also investigated the effects of chemical composition and tempering of martensitic stainless steel on corrosion rates and pitting potentials in simulated CO2 gas wells [13]. They reported that the corrosion rate and the pitting potential in environments simulating CO2 gas wells were governed mainly by the amount of Cr in the steel matrix. The effect of acetate ions on the structure of surface oxide films and corrosion behavior of 13 mass% Cr stainless steel in 25% NaCl solutions were investigated by X-ray photoelectron spectroscopy (XPS) and electrochemical techniques [18]. It was found that the intensity ratios of both Fe0/(Fe3+ + Fe2+ + Fe0) and Cr0/(Cr3+ + Cr0) decrease with increasing concentration of CH3COONa, and after a minimum increased with further increasing in CH3COONa concentration. Fierro et al. applied XPS for the analysis of the corrosion film formed on 13 mass% Cr martensitic stainless steel in CO2-H2S-Cl− [19,20]. The presence of H2S resulted in more the formation of the less protective Fe sulfide-rich and Fe oxide-rich layers, while the presence of CO2 resulted in Cr oxide enrichment in the outermost layers of the film. Islam et al. investigated the corrosion layer on pipeline steel in CO2 containing environments (sweet environments) [21].
However, the effects of environmental parameters such as pH, the concentration of CH3COONa, and CO2 on corrosion behavior of stainless steels have not been fully elucidated. The purpose of this study is to investigate the effects of the concentrations of CH3COONa and CO2 on the corrosion behavior of 13 mass% Cr martensitic stainless steels in simulated oil and gas environments.
Plates of 13 mass% Cr martensitic stainless steel (12.4 mass% Cr, 0.2 mass% C) with 3 mm thickness were used as samples. The sample size for electrochemical and XPS analysis was 7 × 15 mm, and for Fourier transform infrared spectroscopy (FT-IR) analysis was 15 × 15 mm.
For electrochemical measurement, Cu wires were soldered to the samples and then embedded in epoxy resin. The exposed surface of the embedded samples was ground with SiC abrasive paper from #240 to #1200 grit size. For immersion corrosion tests, samples were embedded in epoxy resin. The exposed surfaces of the samples were ground with SiC abrasive paper from #240 to #2400 grit size. The samples were ultrasonically cleaned in ethanol and in highly purified water, before the tests.
Solutions of 5% NaCl with 0, 0.4, 4.0 g/L CH3COONa were used as test solutions. The pH of the solutions was adjusted to 3.5 by HCl. Before the experiments, N2 gas was bubbled through the solutions for 3.6 ks to reduce the concentration of dissolved O2 as low as possible, and then CO2 gas was bubbled as needed. The solution temperature was 353 K, and detailed procedures of the gas bubbling and pH adjustment are illustrated in Figure 1. During experiments, N2 or CO2 gas flowed over the solutions.
Anodic potentiodynamic polarization measurements were carried out in a three electrode cell. A Pt plate and Ag/AgCl sat. KCl was used as counter and reference electrodes respectively. The samples were immersed in the solutions for 21.6 ks and measured the immersion potential, Eim, and then the polarization measurements with a scan rate of 0.5 mV/s were started from (Eim - 10) mV until current density reached 1 mA/cm2. After the measurements, samples were rinsed with highly purified water and dried in a desiccator, and then the surface observations were performed with an optical microscope. Samples with crevice corrosion were discarded.
The samples were immersed in the solutions for 21.6 ks. After the immersion, samples were rinsed with highly purified water and dried in the desiccator. The sample surfaces were analyzed by XPS (JEOL Ltd., JPS-9200). For the XPS analysis, Mg Kα (1253.6 eV) was used as the X-ray source, and the analysis area was 1 × 1 mm. The XPS depth analysis was carried out by sputtering with Ar ions. The sputtering time was converted to depth by the sputtering rate of SiO2. The wide range of depth profiles of Fe, Cr, and O were obtained by GDS (Horiba-Jobin Yvon, JY 5000 RF). The analysis area was 4 mm in diameter. The FT-IR (JASCO, FT-IR 4700) was used to determine the organic substances on the samples. The IR spectra were measured by reflection absorption spectroscopy (RAS) with the resolution of 1 cm−1. The background spectra were obtained with gold evaporation mirror and the subtracted the obtained samples spectra. The surface observation was carried out with a scanning electron microscope, SEM (JEOL Ltd., JSL6510-LA).
Anodic polarization curves were measured after 21.6 ks immersion and are shown in Figure 2. All curves except for 0 g/L CH3COONa in the presence of CO2 have the plateau region in the current, which implies that the corrosion reaction reaches a temporary steady state. The currents of plateau regions are higher than that of general passive current [22]. This result suggests that the passive state may not be achieved in the experimental conditions of this study. The current suddenly increased as the potential moves to a positive direction, which may imply the initiation of pitting corrosion. Hereafter, the plateau region in the current is called the average current density of the flat region, iflat, and the potential at 1 mA/cm2 is called pitting potential, Epit.
The relationships between these two parameters, iflat and Epit, and the CH3COONa concentration are shown in Figure 3. In Figure 3a, the iflat decreases with increasing CH3COONa concentration independent of CO2 bubbling, implying that adding CH3COONa inhibits the corrosion reaction. In Figure 3b, Epit obtained in the solutions of the presence of CO2 increases with increasing CH3COONa concentration, however, the effect of CH3COONa addition is not significant in the solutions of an absence of CO2.
The surface optical images of the samples after the polarization tests are shown in Figure 4. Dark spots which may be related to pitting corrosion are observed, and the number of the spot decreased with increasing concentration of CH3COONa. This result coincides with Figure 3b.
The sample surfaces were examined by XPS after immersion in 4.0 g/L CH3COONa absence of CO2 for 21.6 ks. Figure 5 shows XPS narrow scan spectra of (a) Fe 2p3/2, (b) Cr 2p3/2, (c) O 1s, and (d) C 1s at different depths. In the narrow scan spectra of Fe 2p3/2 (Figure 5a), no peak is observed from 0 and 10 nm, while peaks around 707 eV are measured depth deeper than 10 nm. The peaks are attributed to Fe0+ [23,24,25,26,27], and the intensity of the peak increases with depth. In the narrow scan spectra of Cr 2p3/2 (Figure 5b), the peaks around 577 eV are observed depth from 0 nm to 70 nm. The peaks are attributed to Cr3+ [23,25,26,27], and the intensity decreases with depth. The peaks around 574 eV are appeared after sputtering. The peaks are attributed to Cr0+ [23,25,26,27], and the intensity of the peak increases with depth.
In Figure 5c, the peaks around 532 eV attributed to OH− [23,25,26,27] is observed on the surface (at 0 nm depth). As sputtering, the peaks around 530 eV appear. The peaks are attributed to O2− [23,25,26,27], and the intensity changed with depth. Form Fe, Cr and O narrow scan spectra, corrosion product layer formed on the steel surface has a double layer structure as an outer layer; Cr hydroxide, and an inner layer; Cr oxide. In Figure 5d, peaks are observed on the surface (at 0 nm depth), and no peak is observed after sputtering. The peaks could not be attributed, however, the result suggests that the CH3COO− in the solution may exist only on the surface.
To investigate the composition of sample surface after immersion, a composition ration of Fe and Cr was calculated from XPS narrow scan spectra. Figure 6 shows composition ratio of Fe/Cr as a function of depth. The ratio increases with increasing depth, however, no significant difference with solution composition is observed. These results indicate that the composition ratio of Fe/Cr ratio at the surface is independent of the chemical composition of the used solution and CO2 bubbling.
It is needed to investigate the depth profile of elements in the sample after the immersion tests. GDS was applied to measure the depth profiles of Fe, Cr and O. The depth profile of the sample immersed in 4.0 g/L CH3COONa solution absence of CO2 is shown in Figure 7a. The Fe intensity is low at the initial and it increases sharply as sputter time longer than 5 s. The Cr intensity shows a peak about 2 s and a minimum about 5 s. The O intensity shows peak around 1 s and becomes low after 5 s sputtering. These results suggest that the Cr rich layer exists until Cr intensity reached minimum. The sputter time at a minimum of Cr intensity decided from GDS spectra. Figure 7b shows the sputtering time at a minimum of Cr intensity as a function of CH3COONa concentration. The sputter time at minimum of Cr intensity decreases with increasing CH3COONa concentration. This result indicates that Cr rich layer becomes thinner with increasing CH3COONa concentration.
It was suggested that organic substances in the solution adsorbed on the surface from XPS analysis. FT-IR applied to analyze the organic substances on the surface, and obtained spectra are shown in Figure 8. The background spectrum was subtracted to reduce peaks related to CO2 and H2O in the environments. However, the residue of CO2 related peaks are still observed around 2400~2300 cm−1 [28,29]. The peaks around 1700~1300 cm−1 are attributed to CH3COO− or HCO3− [30,31], and the intensity of the peaks increases with increasing CH3COONa concentration. The intensity obtained in presence of CO2 shows lager value independent of CH3COONa concentration. These results indicate that CH3COO− or CO32− adsorbs on the surfaces.
Surface morphologies of the samples after 21.6 ks immersion are shown in Figure 9. In Figure 9a, b, uneven surface with dark lines are observed in 0 g/L CH3COONa. In Figure 9c, d, a small number of dark lines are observed in 0.4 g/L CH3COONa. However, their size and length are almost the same as observed in 0 g/L CH3COONa. Smooth surfaces without dark lines are observed in 4.0 g/L CH3COONa (Figure 9e, f). There are precipitates on the surfaces in all samples, and size and number are changed with CH3COONa concentration. Little difference of the surface morphology is observed in the presence and absence of CO2 in the solutions.
It is known that dissolved CO2 produced carbonic acid and H+.
$ {\rm{C}}{{\rm{O}}_{\rm{2}}}{\rm{ + }}{{\rm{H}}_{\rm{2}}}{\rm{O}} \to {{\rm{H}}_{\rm{2}}}{\rm{C}}{{\rm{O}}_{\rm{3}}} $ | (1) |
$ {{\rm{H}}_{\rm{2}}}{\rm{C}}{{\rm{O}}_3} \to {\rm{HC}}{{\rm{O}}_{\rm{3}}}^{\rm{ - }}{\rm{ + }}{{\rm{H}}^{\rm{ + }}} $ | (2) |
In this experiment, it is considered that the cathodic reaction is the hydrogen evolution reaction because the pH of the used solutions was 3.5. The reaction may change the pH of the solution during polarization measurement. If CO2 gas flowing during the experiments, Eqs 1 and 2 may occur continuously and the carbonic acid acts as the source of H+. Therefore, the pH may be maintained during the polarization tests. This is a possible reason for iflat of 0 g/L CH3COONa in presence of CO2 showed the largest value as shown in Figure 3a. The supply of H+ by Eq 2 in presence of CO2, make it possible to increase corrosiveness of the solution and may reduce Epit as shown in Figure 3b.
In this experiment, iflat decreased with increasing concentration of CH3COONa, while a small difference in surface composition was observed by surface analysis. The thickness of the Cr rich layer changed with the concentration of CH3COONa that was observed by GDS, and the adsorptions of CH3COO− or HCO3− on the surface were observed by FT-IR. These results may closely relate to corrosion process. Schematic representation of the corrosion process in these experimental conditions is shown in Figure 10. A highly protective layer is formed on the steel surface when the absence of CH3COONa and CO2 in the solution (Figure 10a). The corrosion is progressed by Cl−. The Cr rich layer is formed on the surface because of preferential Fe dissolution [32,33]. The thickness of the Cr rich layer may become thicker and then reaches a certain value as metal dissolution. In the case of CH3COONa and CO2 are present in the solutions, CH3COO− and HCO3− adsorb on the steel surface (Figure 10b). The absorbed layer may inhibit the dissolution of metals, therefore, the amount of dissolved metals decreases and relatively thin Cr rich layer is formed.
The influences of CH3COONa and CO2 on the corrosion behavior of 13 mass% Cr martensitic steels in simulated oil and gas environments were investigated with electrochemical and surface analysis techniques. Following conclusions may be drawn:
(1) There was the plateau region of the current density, iflat, in the anodic polarization curves, and the iflat decreased with increasing CH3COONa concentration. In the case of the solutions with CH3COONa, little effect of CO2 on the iflat was observed. A pitting potential, Epit shifted to the positive direction as increasing CH3COONa concentration, and showed the negative value by adding CO2.
(2) The Cr enriched layer formed on the surface of the sample after the immersion test, and the thickness of the layer became thinner as increasing CH3COONa concentration. The adsorbed CH3COO− and HCO3− on the sample surfaces prevented the corrosion.
An observation of SEM was conducted at the Laboratory of XPS analysis, Joint-use facilities, Hokkaido University, supported by Material Analysis and Structure Analysis Open Unit (MASAOU). XPS analysis was conducted at the Laboratory of XPS analysis, Joint-use facilities, Hokkaido University, supported by 'Nanotechnology Platform' Program of the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan.
All authors have no conflict of interest directly relevant to the contents of this article.
[1] | US EPA (2011) Biofuels and the Environment: First Triennial Report to Congress. United States Environmental Protection Agency, National Center for Environmental Assessment, Office of Research and Development. EPA/600/R-10/183F. |
[2] | US Cong (2007) Energy Independence and Security Act of 2007. 100th Congress, 1st session HR 6.4. |
[3] | US DOE (2011) US Billion Ton Update: Biomass Supply for a Bioenergy and Bioproducts Industry. Leads: RD Perlack and BJ Stokes. United States Department of Energy, Oak Ridge National Laboratory. ORNL/TM-2011/224. |
[4] | US DOE (2005) Biomass as Feedstock for a Bioenergy and Bioproducts Industry: The Technical Feasibility of a Billion Ton Annual Supply. Leads: RD Perlack, BJ Stokes, and DC Erbach. United States Department of Energy, Oak Ridge National Laboratory ORNL/TM-2005/66. |
[5] | BR & DB (2008) National Biofuels Action Plan October 2008. Biomass Research and Development Board, Biomass Research and Development Initiative. |
[6] |
Haldi J, David W (1967) Economies of Scale in Industrial Plants. J polit econ 75: 373-385. doi: 10.1086/259293
![]() |
[7] | US DOE (2002) Lignocellulosic Biomass to Ethanol Process Design and Economics Utilizing Co-Current Dilute Acid Prehydrolysis and Enzymatic Hydrolysis for Corn Stover. A Aden, M Ruth, K Ibsen, et al., United States Department of Energy, National Renewable Energy Laboratory, NREL/TP-510-32438. |
[8] |
Resop JP, Cundiff JS, Heatwole CD (2011) Spatial analysis of site satellite storage locations for herbaceous biomass in the piedmont of the southeast. Appl eng agric 27: 25-32. doi: 10.13031/2013.36221
![]() |
[9] | Judd JD, Sarin SC, Cundiff JS (2012) Design, Modeling and Analysis of a Feedstock Logistics System. Bioresource Technol 103: 209-218. |
[10] | Kumar A, Sokhansanj S (2007) Switchgrass (Panicum virgatum, L.) delivery to a biorefinery using integrated biomass supply analysis and logistics (IBSAL) model. Bioresource Technol 98: 1033-1044. |
[11] |
Worley JS, Cundiff JS (1991) System analysis of sweet sorghum harvest for ethanol production in the piedmont. Transactions of the ASAE 34: 539-547. doi: 10.13031/2013.31697
![]() |
[12] | US DOE (2009) Commodity-Scale Production of an Infrastructure-Compatible Bulk Solid from Herbaceous Lignocellulosic Biomass: Uniform-Format Bioenergy Feedstock Supply System Design Report Series. Hess JR, Kenney KL, Ovard LP, et al. United States Department of Energy, National Renewable Energy Laboratory. INL/EXT-09-17527. |
[13] |
Shastri Y, Hansen A, Rodriguez L, et al. (2011) Development and application of BioFeed model for optimization of herbaceous feedstock production. Biomass Bioenerg 35: 2961-2974. doi: 10.1016/j.biombioe.2011.03.035
![]() |
[14] | Shastri YN, Rodriguez LF, Hansen AC, et al. (2012) Impact of distributed storage and pre-processing on Miscanthus production and provision systems. Biofuel Bioprod Bior 6: 21-31. |
[15] |
Ravula PP, Grisso RD, Cundiff JS (2008) Cotton logistics as a model for a biomass transportation system. Biomass Bioenerg 32: 314-325. doi: 10.1016/j.biombioe.2007.10.016
![]() |
[16] | Cundiff JS, Grisso RD, Ravula PP (2004) Management for Biomass Delivery at a Conversion Plant. 2004 American Society of Agricultural Engineers and Canadian Society of Agricultural Engineers Annual International Meeting August 1-4, 2004 Ottawa, Ontario. |
[17] | Cundiff JS, Grisso RD (2008) Containerized Handling to Minimize Hauling Cost of Herbaceous Biomass. Biomass Bioenerg 32: 308-313. |
[18] | An H, Willhelm WE, Searcy SW (2011) A mathematical model to design a lignocellulosic biofuel supply chain system with a case study based on a region in Central Texas. Bioresource Technol 102: 7860-7870. |
[19] | Kay, MG (2014) Matlog: Logistics engineering Matlab toolbox. North Carolina State University, Fitts Department of Industrial and Systems Engineering. Available from: http://www.ise.ncsu.edu/kay/matlog/ |
[20] | Caffrey KR, Chinn MS, Veal MW (2015). Biomass supply chain management in North Carolina (part 1): predictive model for cropland conversion to biomass feedstocks. AIMS Energy 4: 256-279. |
[21] | Hughes SR, Gibbons WR, Moser BR, et al., (2013) Sustainable Multipurpose Biorefineries for Third-Generation Biofuels and Value-Added Co-Products. InTech. Available from: www.intechopen.com |
[22] | Saxe C (2004) Big Bale Storage Losses. University of Wisconsin Extension, November 2004. |
[23] | Lazarus WF (2009) Machinery Cost Estimates. University of Minnesota Extension. September 2014. |
[24] | Lazarus WF (2014) Machinery Cost Estimates. University of Minnesota Extension. June 2014. |
[25] | USDA (2015) Cash Rents by County. United States Department of Agriculture, National Agricultural Statistics Service. August 2014. Available at: www.nass.usda.gov |
[26] | Atkinson AD, Rich BA, Tungate KD, et al. (2006) North Carolina Canola Production. North Carolina State University, North Carolina Solar Center & College of Agricultural and Life Sciences. SJS/KEL-9/06-W07. |
[27] |
Xiao R, Cai Z, Zhang X (2012) A production optimization model of supply-driven chain quality uncertainty. J Syst Sci Systems Eng 21: 144-160. doi: 10.1007/s11518-011-5184-8
![]() |
[28] |
Benson GA, Green JT (2013) Corn Silage 2013. North Carolina State University, Department of Agricultural and Resource Economics. Available from: http://ag-econ.ncsu.edu/extension/budgets. |
[29] |
Hwang S, Epplin FM, Lee BH, et al. (2009) A probabilistic estimate of the frequency of mowing and baling days available in Oklahoma USA for the harvest of switchgrass for use in biorefineries. Biomass Bioenerg 33: 1037-1045. doi: 10.1016/j.biombioe.2009.03.003
![]() |
[30] | Crouse D (2003) Realistic yields and nitrogen application factors for North Carolina crops. http://nutrients.soil.ncsu.edu/yields/ North Carolina State University, North Carolina Department of Agriculture and Consumer Services, North Carolina Department of Environment and Natural Resources, Natural Resources Conservation Service. Raleigh NC |
[31] | USDA (2014) Cropland Data Layer 2013. Published crop-specific data layer [online]. United States Department of Agriculture, National Agricultural Statistics Service. Available from: http://nassgeodata.gmu.edu/CropScape. Accessed: July 2014. |
[32] | Bullen G, Weddington E (2012) Enterprise Budgets: Corn- Conventional Till-NC, Coastal Plain 2012. North Carolina State University, Department of Agricultural and Resource Economics. Available from: http://ag-econ.ncsu.edu/extension/budgets. |
[33] | Bullen G, Weddington E (2012). Enterprise Budgets: Corn- No Till-NC, Coastal Plain 2012. North Carolina State University, Department of Agricultural and Resource Economics. Available from: http://ag-econ.ncsu.edu/extension/budgets. |
[34] | NCSU (2014) Enterprise Budget: No-Till Grain Sorghum for the Coastal Plain Region of NC. North Carolina State University, Department of Agricultural and Resource Economics. Available from: http://ag-econ.ncsu.edu/extension/budgets. |
[35] | Bullen G, Dunphy J (2012). Enterprise Budgets: Soybeans- Full Season, Conventional Tillage 2012. North Carolina State University, Department of Agricultural and Resource Economics. Available from: http://ag-econ.ncsu.edu/extension/budgets. |
[36] | Bullen G, Dunphy J (2012) Enterprise Budgets: Soybeans- Full Season, No Till 2012. North Carolina State University, Department of Agricultural and Resource Economics. Available from: http://ag-econ.ncsu.edu/extension/budgets. |
[37] | Bullen G, Weddington E (2012) Enterprise Budgets: Wheat for Grain Conventional 2012. North Carolina State University, Department of Agricultural and Resource Economics. Available from: http://ag-econ.ncsu.edu/extension/budgets. |
[38] | Bullen G, Weddington E (2012) Enterprise Budgets: Wheat for Grain No-Till 2012. North Carolina State University, Department of Agricultural and Resource Economics. Available from: http://ag-econ.ncsu.edu/extension/budgets. |
[39] |
Green JT, Benson GA (2013) Bermuda Grass for Pasture. North Carolina State University, Department of Agricultural and Resource Economics. Available from: http://ag-econ.ncsu.edu/extension/budgets. |
[40] |
Green JT, Benson GA (2013) Bermuda Grass for Hay. North Carolina State University, Department of Agricultural and Resource Economics. Available from: http://ag-econ.ncsu.edu/extension/budgets. |
[41] |
Green JT, Benson GA (2013) Cool Season Perennial Grass for Pasture. North Carolina State University, Department of Agricultural and Resource Economics. Available from: http://ag-econ.ncsu.edu/extension/budgets. |
[42] |
Green JT, Benson GA (2013) Hay Harvest Cost, Large Round Baler. North Carolina State University, Department of Agricultural and Resource Economics. Available from: http://ag-econ.ncsu.edu/extension/budgets. |
[43] | Cooper L (1964) Heuristic Method for Location-Allocation Problems. SIAM Review 6: 36-53. |
[44] |
US BOC (2014) Topologically Integrated Geographic Encoding and Referencing 2010. United States Bureau of the Census. Available from: https://www.census.gov/geo/maps-data/data/tiger.html |
[45] |
US DOI (2014) National Land Cover Database 2011. United States Department of the Interior, Multi-Resolution Land Characteristics Consortium (MRLC). Available from: http://www.mrlc.gov/nlcd2011.php |
[46] | Love R, Morris J, Wesolowsky G (1988) Facilities Location: Models and Methods, New York: North-Holland. |
1. | C. Thinaharan, R. P. George, John Philip, In-Situ Detection of Early Corrosion of Ferritic Cr-Mo Steel in Aqueous Solutions of different Anions using Laser Raman Spectroscopy, 2021, 168, 0013-4651, 041508, 10.1149/1945-7111/abf5a9 | |
2. | Masaki Kato, Masatoshi Sakairi, Analysis of the Oxide Film Formed on Stainless Steel Containing Mo in Simulated Condensed Water of Oil and Gas Environments, 2023, 72, 0917-0480, 64, 10.3323/jcorr.72.64 |