Heat Chap05-043

合集下载

Training Process ESEC 2009SSI Level 3 中文版

Training Process ESEC 2009SSI Level 3 中文版

Presentation Name
company confidential
Page 19
Machine working principle 机器工作原理
特点 Single or dual wire mode Fully motorized in X, Y, Z Solder wire diameter range
Welcome to the World of
Soft Solder Process Training 2009SSI
Softsolder Technology
Time Planning
09:00AM – 09:30AM 09:30AM – 10:00AM 10:30AM – 10:45AM 10:45AM – 11:15AM 11:15PM – 12:00PM 13:30PM – 14:30PM 14:45PM – 15:00PM 15:00PM – 17:00PM
Page 4
<Name Presentaion>
Presentation Name
company confidential
Page 4
Machine Overview
Page 5
<Name Presentaion>
Presentation Name
company confidential
Page 5
Flow of Material材料流程
Leadframe 框架
Solder wire 焊锡丝
Wafer晶元
Die Bonded Leadframes Loaded in magazine 粘有芯片的框架装载 入料盒中
Page 8
<Name Presentaion>

【听力训练】School-leavers’Party

【听力训练】School-leavers’Party

Listening activity—School-leavers’ Party
Listen to the conversation again and fill in the blanks with the words in the box.
balloons calendar fork knife paint spoon
The day for the school-leavers’ party is an important date in the school 1. _c_a_le_n_d_a_r_. Tony is going to bring some 2. __b_a_ll_o_o_n_s_ and 3. _p__a_in_t_ some pictures for the party. Everyone is going to bring a traditional dish that can be eaten with their fingers. Soup is no good because it is not finger food and people need a 4. _s_p_o_o_n_ for it. Anything that needs a 5. __k_n_i_f_e_ and 6. __fo_r_k__ is not finger food either.
Listening activity—School-leavers’ion and answer the following questions.
1. What was Daming chosen to do at the school-leavers party? Play the dance music.

ASTM A 743-2013 一般用耐腐蚀铬铁及镍铬铁合金铸件

ASTM A 743-2013 一般用耐腐蚀铬铁及镍铬铁合金铸件
4.1.2 Grade, 4.1.3 Heat treatment, 4.1.4 Options in the specification, 4.1.5 Whether castings are to be produced using the investment casting process, and
4. Ordering Information
4.1 Orders for material to this specification should include the following, as required, to describe the material adequately:
4.1.1 Description of the casting by pattern number or drawing,
2 For referenced ASTM standards, visit the ASTM website, , or contact ASTM Customer Service at service@. For Annual Book of ASTM Standards volume information, refer to the standard’s Document Summary page on the ASTM website.
1.2 The values stated in either inch-pound units or SI units are to be regarded separately as standard. Within the text, the SI units are shown in brackets. The values stated in each system are not exact equivalents; therefore, each system must be used independently of the other. Combining values from the two systems may result in nonconformance with the specification. Inch-pound units are applicable for material ordered to Specification A743 and SI units for material ordered to Specification A743M.

PANTONE国际色卡C面效果电子版

PANTONE国际色卡C面效果电子版

矿产资源开发利用方案编写内容要求及审查大纲
矿产资源开发利用方案编写内容要求及《矿产资源开发利用方案》审查大纲一、概述
㈠矿区位置、隶属关系和企业性质。

如为改扩建矿山, 应说明矿山现状、
特点及存在的主要问题。

㈡编制依据
(1简述项目前期工作进展情况及与有关方面对项目的意向性协议情况。

(2 列出开发利用方案编制所依据的主要基础性资料的名称。

如经储量管理部门认定的矿区地质勘探报告、选矿试验报告、加工利用试验报告、工程地质初评资料、矿区水文资料和供水资料等。

对改、扩建矿山应有生产实际资料, 如矿山总平面现状图、矿床开拓系统图、采场现状图和主要采选设备清单等。

二、矿产品需求现状和预测
㈠该矿产在国内需求情况和市场供应情况
1、矿产品现状及加工利用趋向。

2、国内近、远期的需求量及主要销向预测。

㈡产品价格分析
1、国内矿产品价格现状。

2、矿产品价格稳定性及变化趋势。

三、矿产资源概况
㈠矿区总体概况
1、矿区总体规划情况。

2、矿区矿产资源概况。

3、该设计与矿区总体开发的关系。

㈡该设计项目的资源概况
1、矿床地质及构造特征。

2、矿床开采技术条件及水文地质条件。

肿瘤细胞株查询

肿瘤细胞株查询

ATCC® No.
CCL-2 CCL-2.1 CCL-2.2 CRL-1958 CRL-7850* CRL-13002† CRL-13003† CRL-13011† CCL-249 CRL-7213* CRL-10632† HTB-40 CRL-1677 CRL-7911* HTB-45 CRL-5866 CRL-5868 CRL-7380* HTB-57 CRL-2869 CRL-2871 CRL-2868 CRL-7222* CRL-7226* CRL-7227* CRL-7245* CRL-7253* CRL-7477* CRL-7480* CRL-2116 CRL-1666 CRL-1743 CRL-1747 CRL-6602* HTB-75 CRL-1687 CRL-1997
Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma Adenocarcinoma, bronchogenic Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, colorectal Adenocarcinoma, ductal; cystic fibrosis Adenocarcinoma, ductal Adenocarcinoma, gastric Adenocarcinoma, ileocecal colorectal Adenocarcinoma, large cell, non-small cell lung cancer Adenocarcinoma, malignant, clear cell carcinoma Adenocarcinoma, malignant, endometrioid carcinoma

DNVGL-RP-0034_2015-02

DNVGL-RP-0034_2015-02

RECOMMENDED PRACTICEThe electronic pdf version of this document found through is the officially binding version.The documents are available free of charge in PDF format.DNVGL-RP-0034Edition February 2015Steel forgings for subsea applications© DNV GL ASAny comments may be sent by e-mail to rules@This service document has been prepared based on available knowledge, technology and/or information at the time of issuance of this document. The use of this document by others than DNV GL is at the user's sole risk. DNV GL does not accept any liability or responsibility for loss or damages resulting from any use of FOREWORDDNV GL recommended practices contain sound engineering practice and guidance.C h a n g e s – c u r r e n tCHANGES – CURRENTGeneralThis is a new document.Acknowledgements:This recommended practice was developed by a Joint Industry Project (JIP). The work was performed by DNV GL and discussed in regular project meetings and workshops with individuals from the participating companies. They are hereby acknowledged for their valuable and constructive input. In case consensus has not been achievable, DNV GL has sought to provide acceptable compromise.Sponsors of the JIP included the following organisations:Further organisations have participated in the review process. DNV GL is grateful for the valuable co-operations and discussions with individuals in these organisations.Aker Solutions Brück CelsaChevron Det Norske Dril-Quip Ellwood Group Eni ExxonMobilFMC FrisaGEJapan Steel Works LundinOneSubsea Petrobras Ringmill Scana SubseaShellStatoilTotalC o n t e n t sCONTENTSCHANGES – CURRENT .................................................................................................. 3Sec.1General ......................................................................................................... 61.1Introduction...........................................................................................61.2Scope and application ............................................................................61.3Steel forging classes ..............................................................................61.4Purchase order information....................................................................71.5Normative and informative references ...................................................71.6Definitions..............................................................................................91.7Verbal forms...........................................................................................91.8Abbreviations.........................................................................................9Sec.2Quality assurance and quality control (11)2.1Quality management system................................................................112.2Manufacturing procedure specification.................................................112.3Inspection and test plan.......................................................................122.4Marking and traceability.......................................................................122.5Certification and documentation ..........................................................12Sec.3Technical provisions (14)3.1Manufacturing practices.......................................................................143.1.1General ......................................................................................143.1.2Melting, refining, and casting.........................................................143.1.3Forging.......................................................................................143.1.4Heat treatment............................................................................153.2Chemical composition...........................................................................153.2.1Heat analysis...............................................................................153.2.2Product analysis...........................................................................153.3Mechanical testing ...............................................................................163.3.1Test laboratories..........................................................................163.3.2Test material...............................................................................163.3.3Test sampling..............................................................................163.3.4Test methods ..............................................................................173.3.5Mechanical properties...................................................................173.3.6Hardness testing..........................................................................173.3.7Re-testing...................................................................................183.3.8Testing after simulated post weld heat treatment ............................183.4Metallographic examination .................................................................183.5Non-destructive testing........................................................................183.5.1General ......................................................................................183.5.2Visual testing ..............................................................................193.5.3Magnetic particle testing...............................................................193.5.4Ultrasonic testing.........................................................................193.6Dimensional inspection ........................................................................203.7Repair ..................................................................................................20Sec.4Manufacturing procedure qualification (21)4.1General ................................................................................................214.2Qualification testing.............................................................................214.2.1Chemical composition. (21)C o n t e n t s4.2.2Tensile and Charpy V-notch testing ................................................214.2.3Hardness testing..........................................................................224.2.4Testing after simulated post weld heat treatment.............................224.2.5Metallographic examination...........................................................224.2.6Fracture toughness testing............................................................234.2.7Non-destructive testing.................................................................234.2.8Dimensional inspection .. (23)4.3Validity (23)SECTION 1 GENERAL1.1 IntroductionThis recommended practice (RP) contains criteria, technical requirements and guidance on qualification, manufacture and testing of carbon and low alloy steel forgings for subsea applications.The RP has been written for general world-wide application. Governmental regulations may include requirements in excess of the provisions given by this RP.The objectives of this RP are to:a)provide an internationally acceptable and harmonised standard for carbon and low alloy steel forgingsb)provide for reduced lead time, enhanced stock keeping and interchangeabilityc)provide consistent quality to increase reliability/integrity of subsea equipmentd)simplify the risk assessment processe)serve as a contractual reference document between manufacturers and purchasersf)serve as a guideline for designers, suppliers, purchasers and regulatorsg)comply with and complement existing industry codes for subsea equipment.The RP is divided into four main sections:—Sec.1 General: Contains introduction, scope and application, information to be supplied by purchaser, normative and informative references, definitions and abbreviations.—Sec.2 Quality assurance and quality control: Contains requirements for quality management system, manufacturing procedure specification, inspection and test plan, marking and traceability,documentation and certification.—Sec.3 Technical provisions: Contains requirements for manufacture, testing and inspection of production parts.—Sec.4 Manufacturing procedure qualification: Contains requirements for qualification of manufacturing processes, methods, procedures and validity.1.2 Scope and applicationThe steel forgings covered herein are intended for components in subsea equipment. Typical applications include:a)subsea wellhead and tree equipment as per ISO 13628-4 or API 17Db)completion/workover riser systems as per ISO 13628-7 or API 17Gc)subsea structures and manifolds as per ISO 13628-15 or API 17Pd)flexible pipe as per ISO 13628-11 or API 17Be)drill-through equipment as per ISO 13533 or API 16A.The material grades covered are carbon steels, micro-alloyed steels, low alloy steels and modified grades thereof as per material groups 1 and 2 in API 20B.This RP is not intended to inhibit a vendor from offering, or the purchaser from accepting, alternative materials or manufacturing processes. This can be particularly applicable where there is innovative or developing technology. Where an alternative is offered, it is the responsibility of the vendor to identify any variations from this RP and provide details to the purchaser.1.3 Steel forging classesThis RP establishes requirements for three steel forging classes (SFC) designated SFC 1, SFC 2, and SFC 3. These SFC designations define different levels of forged product technical, quality and qualification requirements.All classes are intended for equipment that shall meet product specification level (PSL) 3 and 3Grequirements in the referenced industry codes, e.g. ISO 10423 or API 6A.The three forging classes reflect increasing criticality as defined by the end user and, hence, increasing requirements from SFC 1 to SFC 3: —SFC 1 is intended for less critical components, e.g. components that are not subjected to continuous exposure to flowing hydrocarbons or components with simple shapes.—SFC 2 is intended for pressure containing and/or load bearing components that are of significant enough size and complexity to warrant additional mechanical testing and surveillance.—SFC 3 is intended for fatigue sensitive pressure containing and/or load bearing components.1.4 Purchase order informationThe purchaser shall provide at least the following information in the order:a)that the forging(s) shall be made according to this RPb)the steel forging class (SFC);c)the steel designation (name or number);d)the quantity of forgings required;e)the drawing number(s) containing the dimensions, tolerances and surface finish;f)the position and thickness of the critical section(s) in the part;g)the minimum design temperature (MDT)/Charpy V-notch (CVN) test temperature;h)the specified minimum yield strength (SMYS);i)the marking requirements for the forging(s);j)the type of certification document;k)whether the forging(s) shall be ISO 15156-2 or NACE MR0175 compliant;l)whether the forging(s) shall be subjected to simulated post weld heat treatment (SPWHT);m)whether forging plan/sketch shall be submitted for review or approval prior to production;n)whether test sample drawing shall be submitted for review or approval prior to production;o)whether ultrasonic testing (UT) procedure and scan plan shall be submitted for review or approval prior to production;p)whether magnetic particle testing (MT) procedure and test coverage description shall be submitted for review or approval prior to production;q)whether manufacturing procedure specification (MPS) shall be submitted for review or approval prior to production;r)whether inspection and test plan (ITP) shall be submitted for mark-up, review or approval prior to production;s)whether manufacturing procedure qualification (MPQ) report shall be submitted for review or approval prior to production;t)whether any additional requirements shall apply.1.5 Normative and informative referencesThe codes and standards in Table 1-1 and Table 1-2 include provisions and guidance which, through reference in this text, constitute provisions and guidance of this RP. The latest edition applies unless dated references are given.Other recognised codes and standards may be used provided it can be demonstrated that these meet or exceed the requirements of the referenced codes and standards.Any deviations, exceptions and modifications to the codes and standards shall be documented and agreed between the manufacturer and purchaser.Table 1-1 Normative referencesAPI RP 6HT Heat Treatment and Testing of Carbon and Low Alloy Steel Large Cross Section andCritical Section ComponentsAPI Spec 20B Open Die Shaped Forgings for Use in the Petroleum and Natural Gas IndustryACCP-CP-1American Society for Nondestructive Testing Central Certification ProgramASME BPVC-V ASME Boiler and Pressure Vessel Code (BPVC), Section V: Nondestructive Examination ASTM A370Standard Test Methods and Definitions for Mechanical Testing of Steel ProductsASTM A388Standard Practice for Ultrasonic Examination of Steel ForgingsASTM A604Standard Practice for Macroetch Testing of Consumable Electrode Remelted Steel Barsand BilletsASTM A694Standard Specification for Carbon and Alloy Steel Forgings for Pipe Flanges, Fittings,Valves, and Parts for High-Pressure Transmission ServiceASTM A707Standard Specification for Forged Carbon and Alloy Steel Flanges for Low-TemperatureServiceASTM A751Standard Test Methods, Practices, and Terminology for Chemical Analysis of SteelProductsASTM A788Standard Specification for Steel Forgings, General RequirementsASTM E45Standard Test Methods for Determining the Inclusion Content of SteelASTM E110Standard Test Method for Indentation Hardness of MetallicMaterials by Portable Hardness TestersASTM E112Standard Test Methods for Determining Average Grain SizeASTM E381Standard Method of Macroetch Testing Steel Bars, Billets, Blooms, and ForgingsASTM E709Standard Guide for Magnetic Particle TestingASTM E1820Standard Test Methods for Measurement of Fracture ToughnessEN 10204Metallic products - Types of inspection documentsISO 643Steels - Micrographic determination of the apparent grain sizeISO 9001Quality Management Systems – RequirementsISO 9712Non-destructive testing - Qualification and certification of NDT personnelISO 10423/API Spec 6A Petroleum and natural gas industries - Drilling and production equipment - Wellhead andchristmas tree equipment/Specification for Wellhead and Christmas Tree equipment ISO 10474Steel and steel products - Inspection documentsISO 12135Metallic materials - Unified method of test for the determination of quasistatic fracturetoughnessISO 15156-2/NACE MR0175Petroleum and natural gas industries - Materials for use in H2S-containing environmentsin oil and gas production - Part 2: Cracking-resistant carbon and low-alloy steels, and theuse of cast irons/Petroleum and natural gas industries - Materials for use in H2S-containing environments in oil and gas productionISO/IEC 17020Conformity assessment - General criteria for the operation of various types of bodiesperforming inspectionISO/IEC 17025General requirements for the competence of testing and calibration laboratoriesTable 1-2 Informative referencesAPI RP 17B Recommended Practice for Flexible PipeAPI Spec 16A Specification for Drill Through EquipmentISO 13533Petroleum and natural gas industries - Drilling and production equipment - Drill-throughequipmentISO 13628-1/API RP 17A Petroleum and natural gas industries - Design and operation of subseaproduction systems - Part 1: General requirements and recommendations/Design andOperation of Subsea Production Systems-General Requirements and Recommendations ISO 13628-4/API Spec 17D Petroleum and natural gas industries - Design and operation of subsea productionsystems - Part 4: Subsea wellhead and tree equipment/Design and Operation of SubseaProduction Systems-Subsea Wellhead and Tree EquipmentISO 13628-7/API RP 17G Petroleum and natural gas industries - Design and operation of subsea productionsystems - Part 7: Completion/workover riser systems/Recommended Practice forCompletion/Workover RisersISO 13628-11Petroleum and natural gas industries - Design and operation of subsea productionsystems - Part 11: Flexible pipe systems for subsea and marine applicationsISO 13628-15/API RP 17P Petroleum and natural gas industries - Design and operation of subsea productionsystems - Part 15: Subsea structures and manifolds/Design and Operation of SubseaProduction Systems - Subsea Structures and Manifolds1.6 DefinitionsTable 1-3 Definitionscritical section(s)section(s) of the forging in which mechanical properties have to meet the specified minimum requirements and are deemed critical to the design and safe operation of the component.1.7 Verbal formsTable 1-4 Verbal formsshall verbal form used to indicate requirements strictly to be followed in order to conform to the document. should verbal form used to indicate that among several possibilities one is recommended as particularly suitable, without mentioning or excluding others, or that a certain course of action is preferred but not necessarilyrequired.may verbal form used to indicate a course of action permissible within the limits of the document.1.8 AbbreviationsTable 1-5 AbbreviationsCE (IIW)carbon equivalent (international institute of welding)CS carbon steelCTOD crack tip opening displacementCVN Charpy V-notchDAC distance amplitude correctionDGS distance gain-sizeEAF electric arc furnaceESR electro slag re-meltingFBH flat bottom holeITP inspection and test planLAS low alloy steelLR ladle refiningMDT minimum design temperatureMPQ manufacturing procedure qualificationMPS manufacturing procedure specificationMT magnetic particle testingTable 1-5 Abbreviations (Continued)MTR material test reportNDT non-destructive testingPWHT post weld heat treatmentQMS quality management systemRP recommended practiceSFC steel forging classSMYS specified minimum yield strength SPWHT simulated post weld heat treatment UT ultrasonic testingVAR vacuum arc re-meltingVD vacuum degassingVT visual testingSECTION 2 QUALITY ASSURANCE AND QUALITY CONTROL2.1 Quality management systemThe forging manufacturers and steel manufacturers shall have a certified quality management system (QMS) conforming to ISO 9001.2.2 Manufacturing procedure specificationAll production shall be based on an manufacturing procedure specification (MPS) established by the forging manufacturer. The MPS shall describe how the specified properties will be achieved and verified. The MPS shall address all factors that affect the quality and reliability of production including all subcontractors applied. Every principal production step from starting material to shipment of finished product(s) shall be addressed. References to the detailed procedures used for the execution of all steps shall be included.As a minimum, the MPS shall include the information in Table 2-1.Table 2-1 Content of MPS1) General a) A descriptive title and a unique identification number with revision controlb)Description of product(s) and size range(s)c)Reference to applicable standards and specificationsd)Reference to MPQ(s)2) Starting materials a)Steel manufacturer, steel grade, melt practice, refining and casting methodb)Specification for chemical composition and carbon equivalent (CE)c)Ingot typed)Methods and practices for ingot discard3) Forging a)Forging method, e.g. open die, closed die or ring rollingb)Forging press/equipment capacity, as applicablec)Hot work temperature range and method of temperature monitoring during forgingd)Description of basic forging steps, i.e. sequence of upsetting, drawing, etce)Forging reduction calculation method for each step and minimum overall forging reductionratiof)Sketch of forging(s) including prolongation, if used, in as-forged condition4) Heat treatment a)Heat treat condition of supply, e.g. quenched and temperedb)Sketch of forging(s) including prolongation, if used, at time of heat treatment, if anydifferent from as-forged geometry. Finished geometry, if known, shall be given with dashedlinesc)Maximum thickness of forging(s) at time of heat treatmentd)Description of furnace loading practice with typical sketch(es) showing maximum loadingweight, location and minimum spacing of parts in the furnace, location of prolongation orsacrificial part and location of contact thermocouples/heat sinkse)Description of heat treatment cycles, temperatures and timesf)Minimum quench tank size/volumeg)Quenching medium and type of agitationh)Quenching medium start and finish temperature and maximum transfer time to quenchi)Maximum surface metal temperature at removal from quench tank including how and whentemperature is measured5) Mechanical testing, metallographic examination, and product analysis a)Specified tensile, Charpy V-notch (CVN) and surface hardness testingb)Sketch showing sampling position and specimen orientation in prolongation or sacrificialforging, as applicablec)Sketch with the locations for surface hardness testingd)Grain size determinatione)Inclusion rating for SFC 3f)Chemical composition determined as product analysis for SFC 32.3 Inspection and test planThe forging manufacturer shall establish an inspection and test plan (ITP) for SFC 2 and SFC 3. The ITP shall have a reference to the relevant MPS and shall list the sequence of activities contained in the MPS. Manufacturers may use their own ITP format, but it shall as a minimum include:a)all principal production, inspection and testing activities b)location of activityc)associated procedure or specification including acceptance criteria governing the activity d)verifying document to be used for recording inspection and test results e)forging manufacturer’s intervention activities f)columns for intervention by purchaser and 3rd party.2.4 Marking and traceabilityThe forgings shall be marked to ensure full traceability to the heat or re-melt ingot (as applicable), heat treatment lot and the certificate representing the forging. Additional marking shall be as specified in the purchase order.Each forging shall be marked with a low stress marking method on a position as stated in the MPS.2.5 Certification and documentationCertification of production forgings shall be as specified by the purchaser. Any of the following certification documents are applicable:— A manufacturer’s test report (MTR) giving the results of all specified tests and inspections.—An inspection document type 3.1 or 3.2 according to ISO 10474 or EN 10204.The certificates shall be supplied by the manufacturer to the purchaser and shall, as a minimum, give the following particulars:a)purchaser’s name, order number and part number b)forging manufacturer’s name and order numberc)description of forging(s) including quantity and drawing number(s)d)reference to this RP e)steel designation and SFC f)reference to the applicable MPS g)steelmaker’s nameh)steelmaking process including secondary refining i)forging method and forging reduction ratio j)name of heat treat subcontractor, if applicablek)heat number(s) and heat treat lot number(s)l)Heat treatment temperatures, soaking times, quenching medium and transfer time from furnace to quench tankm)heat analysis and, where applicable, product analysis6) Non-destructive testing, as applicablea)Visual testingb)Magnetic particle testing, including sketch or description of test coverage c)Ultrasonic testing, including scan plan7) Dimensionalcontrol and marking a)Dimensional controlb)Method, extent, and position of marking 8) Final certificationa)Type of certificateb)Associated records and documentationTable 2-1 Content of MPS (Continued)n)prolongation or sacrificial part dimensionso)results of tensile, CVN, surface hardness, metallographic and any other testing requiredp)NDT report(s)q)marking of forging(s).The certificates shall be accompanied by documentation and records as follows:a)heat treatment chart(s) and furnace loading sketch(es)b)test plan/sketch showing locations for tensile, CVN, surface hardness, metallographic and any othertesting requiredc)NDT procedure(s)d)dimensional inspection report(s)e)for SFC 2 and SFC 3, copy of certificate(s) from steelmaker.The manufacturer shall maintain documentation and records of relevant manufacturing procedure qualification (MPQ).SECTION 3 TECHNICAL PROVISIONS3.1 Manufacturing practices3.1.1 GeneralAll manufacturing shall be based on the MPS, the associated MPQ (for SFC 2 and SFC 3), and the requirements of this RP.3.1.2 Melting, refining, and casting3.1.2.1 The steel shall be melted using the electric arc furnace (EAF) followed by secondary refining such as ladle refining (LR) and vacuum degassing (VD). Secondary re-melt processes such as electro slag re-melting (ESR) and vacuum arc re-melting (VAR) may also be used.3.1.2.2 The steel shall be fully killed and made to a fine grain practice. See also [3.4].3.1.2.3 The steel for SFC 3 shall be treated for inclusion shape control. When Ca treatment is used, Ca shall not exceed 0.005%. See also [3.4].3.1.2.4 The steel shall be ingot cast with bottom pouring, ingot cast with top pouring in vacuum, or continuous cast. Adequate top and bottom ingot discards shall be made to ensure freedom from piping and harmful segregations in the finished forgings. Surface and skin defects, which may be detrimental during the subsequent working and forming operations, shall be removed.3.1.2.5 Repair by welding on ingots, blooms or billets shall not be permitted.3.1.2.6 Melting, refining and casting practices shall be as stated in the MPS.3.1.3 Forging3.1.3.1 Forgings shall be made by any of the following methods: Open die, closed die or ring rolling. 3.1.3.2 The material shall be hot worked, and shall be forged as close as practical to the finished shape and size.3.1.3.3 The hot work temperature shall be monitored during the forging process by pyrometer or equivalent equipment.3.1.3.4 The overall forging (total hot work) reduction ratio shall be minimum4.0:1 for all classes, with the following considerations:a)The initial free upsetting operations of the as cast ingot shall not be considered as part of the overallforging reduction ratio.b)Upsetting following cogging or drawing may be considered as part of the overall forging reduction ratio.If upsetting following cogging or drawing is to be considered, the overall forging reduction ratio shall be minimum 6.0:1.c)For ring rolling or mandrel forging, punching or piercing of holes shall not be considered as part of theoverall forging reduction ratio.d)The overall reduction ratio shall be sufficient to produce a wrought structure throughout the entire part.3.1.3.5 The hot work reduction ratio for a single hot work operation and the total hot work reduction ratio shall be calculated according to API 20B.3.1.3.6 For ring rolling or mandrel forging, the initial cross-sectional area shall be as-punched or pierced wall thickness times as-punched or pierced height. The final cross-sectional area shall be final wall thickness times final height.3.1.3.7 For closed die forging, the reduction ratio shall meet the requirement of this RP prior to the closed die operation.3.1.3.8 Repair by welding on forgings shall not be permitted.。

numerical analysis of heat and mass transfer in the cappillary structure of a loop heat pipe

Numerical analysis of heat and mass transfer in the capillarystructure of a loop heat pipeTarik Kaya *,John GoldakCarleton University,Department of Mechanical and Aerospace Engineering,1125Colonel By Drive,Ottawa,Ont.,Canada K1S 5B6Received 24February 2005;received in revised form 20December 2005Available online 31March 2006AbstractThe heat and mass transfer in the capillary porous structure of a loop heat pipe (LHP)is numerically studied and the LHP boiling limit is investigated.The mass,momentum and energy equations are solved numerically using the finite element method for an evapo-rator cross section.When a separate vapor region is formed inside the capillary structure,the shape of the free boundary is calculated by satisfying the mass and energy balance conditions at the interface.The superheat limits in the capillary structure are estimated by using the cluster nucleation theory.An explanation is provided for the robustness of LHPs to the boiling limit.Ó2006Elsevier Ltd.All rights reserved.Keywords:Two-phase heat transfer;Boiling in porous media;Boiling limit;Loop heat pipes;Capillary pumped loops1.IntroductionTwo-phase capillary pumped heat transfer devices are becoming standard tools to meet the increasingly demand-ing thermal control problems of high-end electronics.Among these devices,loop heat pipes (LHPs)are particu-larly interesting because of several advantages in terms of robust operation,high heat transport capability,operabil-ity against gravity,flexible transport lines and fast diode action.As shown in Fig.1,a typical LHP consists of an evaporator,a reservoir (usually called a compensation chamber),vapor and liquid transport lines and a con-denser.The cross section of a typical evaporator is also shown in Fig.1.The evaporator consists of a liquid-pas-sage core,a capillary porous wick,vapor-evacuation grooves and an outer casing.In many LHPs,a secondarywick between the reservoir and the evaporator is also used to ensure that liquid remains available to the main wick at all times.Heat is applied to the outer casing of the evapo-rator,leading to the evaporation of the liquid inside the wick.The resulting vapor is collected in the vapor grooves and pushed through the vapor transport line towards the condenser.The meniscus formed at the surface or inside the capillary structure naturally adjusts itself to establish a capillary head that matches the total pressure drop in the LHP.The subcooled liquid from the condenser returns to the evaporator core through the reservoir,completing the cycle.Detailed descriptions of the main characteristics and working principles of the LHPs can be found in Maidanik et al.[1]and Ku [2].In this present work,the heat and mass transfer inside the evaporator of an LHP is considered.The formulation of the problem is similar to a previous work performed by Demidov and Yatsenko [3],where the capillary struc-ture contains a vapor region under the fin separated from the liquid region by a free boundary as shown in Fig.2.Demidov and Yatsenko [3]have developed a numerical procedure and studied the growth of the vapor region under increasing heat loads.They also present a qualitative0017-9310/$-see front matter Ó2006Elsevier Ltd.All rights reserved.doi:10.1016/j.ijheatmasstransfer.2006.01.028*Corresponding author.E-mail addresses:tkaya@mae.carleton.ca (T.Kaya),jgoldak@mrco2.carleton.ca (J.Goldak)./locate/ijhmtanalysis of the additional evaporation from the meniscus formed in thefin–wick corner when the vapor region is small without exceeding thefin surface.They report that the evaporation from this meniscus could be much higher than that from the surface of the wick and designs facilitating the formation of the meniscus would be desir-able.Figus et al.[4]have also presented a numerical solu-tion for the problem posed by Demidov and Yatsenko[3] using to a certain extent similar boundary conditions and a different method of solution.First,the solutions are obtained for a single pore-size distribution by using the Darcy model.Then,the solution method is extended to a wick with a varying pore-size distribution by using a two-dimensional pore network model.An important conclusion of this work is that the pore network model results are nearly identical to those of the Darcy model for an ordered single pore-size distribution.On the basis of this study,we consider a capillary structure with an ordered pore distri-bution possessing a characteristics single pore size.A sim-ilar problem has also been studied analytically by Cao and Faghri[5].Unlike[3,4],a completely liquid-saturated wick is considered.Therefore,the interface is located at the sur-face of the wick.They indicate that the boiling limit inside the wick largely depends on the highest temperature under thefin.This statement needs further investigation espe-cially when a vapor region under thefin is present.In a later study,Cao and Faghri[6]have extended their work to a three-dimensional geometry,where a two-dimensional liquid in the wick and three-dimensional vaporflow in the grooves separated by aflat interface at the wick surface is considered.A qualitative discussion of the boiling limit in a capillary structure is provided.They also compare the results of the two-dimensional model without the vapor flow in the grooves and three-dimensional model and con-clude that reasonably accurate results can be obtained by a two-dimensional model especially when the vapor velo-cities are small for certain workingfluids such as Freon-11 and ammonia.Based on these results,in our work,we con-sider a two-dimensional geometry to simplify the formula-tion of the problem.All these referenced works assume a steady-state process.Dynamic phenomena and specifically start-up is also extensively studied[7,8].The superheat at the start-up and temperature overshoots is still not well understood.In this work,the transient regimes and start-up are not investigated.One of the goals of the present study is a detailed inves-tigation of the boiling limit in a capillary structure.There-fore,the completely liquid-saturated and vapor–liquid wick cases are both studied.The boiling limit in a porous struc-ture is calculated by using the method developed by Mish-kinis and Ochterbeck[9]based on the cluster nucleation theory of Kwak and Panton[10].Our primary interest in this study is LHPs.In comparison,the previously refer-enced works focus primarily on capillary pumped loops (CPLs),a closely related two-phase heat transfer device to an LHP.Unlike in a CPL,the proximity of the reservoir to the evaporator in an LHP ensures that the wick is con-tinuously supplied with liquid.However,there is no signif-icant difference in the mathematical modeling of both devices especially because only a cross section of the evap-orator is studied.The main difference here is that LHPs easily tolerate the use of metallic wicks with very small pore sizes,with a typical effective pore radius of1l m,resulting in larger available capillary pressure heads.Nomenclaturec p specific heat at constant pressure[J kgÀ1KÀ1] h c convection heat transfer coefficient[W mÀ2KÀ1] h i interfacial heat transfer coefficient[W mÀ2KÀ1] h fg latent heat of evaporation[J kgÀ1]J nc critical nucleation rate[nuclei mÀ3sÀ1]k thermal conductivity[W mÀ1KÀ1]K permeability[m2]L length[m]p pressure[Pa]D p pressure drop across wick[Pa]Pe Peclet numberQ b heat load for boiling limit[W]q in applied heatflux[W mÀ2]Q in applied heat load[W]r radius[m]r p pore radius[m]Re Reynolds numbert thickness[m]T temperature[K]u velocity vector[m sÀ1]Greek symbolsh angle[degrees]l viscosity[Pa s]q density[kg mÀ3]u porosityr liquid–vapor surface tension[N mÀ1] Subscriptsc casingeffeffectiveg groovein inletint interfacel liquidmax maximumn normal componentsat saturationv vaporw wick3212T.Kaya,J.Goldak/International Journal of Heat and Mass Transfer49(2006)3211–32202.Mathematical formulationA schematic of the computational model for the wick segment studied is shown in Fig.3.Because of the symme-try,a segment of the evaporator cross section is considered,which is between the centerlines of the fin and adjacent vapor groove.The numerical solutions for thisgeometryFig.2.Schematic of evaporation inside theevaporator.Fig.1.Schematic of a typical LHP and cross section of the evaporator.T.Kaya,J.Goldak /International Journal of Heat and Mass Transfer 49(2006)3211–32203213are obtained for two separate configurations.At low heat loads,the wick is entirely saturated by the liquid.At higher heat loads,the wick contains two regions divided by an interface as shown in Fig.3:an all-vapor region in the vicinity of thefin and a liquid region in the remaining part of the wick.Heat is applied on the exterior walls of the cas-ing and it is transferred through thefin and wick to the vapor–liquid interface.This leads to the evaporation of the liquid at the interface and thus theflow of the vapor into the grooves.For the vapor–liquid wick,the vapor formed inside the wick is pushed towards the grooves through a small region at the wick–groove border.In both of the cases,as a result of the pressure difference across the wick,the liquid from the core replaces the outflowing vapor.Under a given heat load,the system reaches the steady state and the operation is maintained as long as the heat load is applied.The mathematical model adopted in this work is based on the following assumptions:the process is steady state; the capillary structure is homogenous and isotropic;radia-tive and gravitational effects are negligible;thefluid is Newtonian and has constant properties at each phase; and there is local thermal equilibrium between the porous structure and the workingfluid.Many of these assump-tions are similar to those made in Demidov and Yatsenko [3]and Figus et al.[4].In addition,we also take into account convective terms in the energy(advection–diffu-sion)equation.The validity of the Darcy equation for the problem studied is also discussed.Under these assump-tions,the governing equations for vapor and liquid phases (continuity,Darcy and energy)are as follows:rÁu¼0ð1Þu¼ÀKlr pð2Þq c p rðu TÞ¼k eff r2Tð3ÞIt should be noted that the Darcy solverfirst calculates the pressure from the Laplace equation for pressure($2p=0), which is obtained by combining Eqs.(1)and(2).The vapor flow in the groove region is not solved to simplify the prob-lem.The boundary conditions for the liquid-saturated wick are described as follows:At r=r ip¼pcore;T¼T satð4ÞAt r=r o and h A6h6h Cu n¼Àk effq l h fgo To n;k effo To n¼h iðTÀT vÞð5ÞAt r=r o and h C6h6h Do p o n ¼0;k co To n¼k effo To nð6ÞAt r=r g and h A6h6h CÀk c o To n¼h cðTÀT vÞð7ÞAt r=r ck co To n¼q inð8ÞAt h=h A and r i6r6r o and r g6r6r co po h¼0;o To h¼0ð9ÞAt h=h C and r o6r6r gÀk co To n¼h cðTÀT vÞð10ÞAt h=h D and r i6r6r co po h¼0;o To h¼0ð11ÞIn the equations above,(o/o n)represents the differentialoperator along the normal vector to a boundary.Theboundary conditions for the wick with the separate vaporand liquid regions are identical to the above equations ex-cept along the wick–groove boundary and for the vapor–liquid interface inside the wick.The following equationssummarize these additional boundary conditions for thevapor–liquid wick:At r=r o and h A6h6h Bu n¼Àk effq l h fgo To n;k effo To n¼h iðTÀT vÞð12ÞAt r=r o and h B6h6h Cp¼pv;o To n¼0ð13ÞThe interface is assumed to have zero thickness.Sharpdiscontinuities of the material properties are maintainedacross the interface.The interfacial conditions are writtenas follows:The mass continuity conditionðu nÞvq v¼ðu nÞlq lð14ÞThe energy conservation conditionðk effÞvo T vo nÀðk effÞlo T lo n¼ðu nÞvqvh fgð15ÞFor the interface temperature condition,we assumelocal thermal equilibrium at the interface inside the wick:T int¼T v¼T lð16ÞHere,we assume that the interface temperature T int is givenby the vapor temperature.This condition is used to locatethe vapor–liquid interface as explained in the followingsection.For the interface at the wick–groove border,a convectiveboundary condition is used,Eqs.(5)and(12).A temper-ature boundary condition ignoring the interfacial resistanceis also possible.The interfacial heat transfer coefficient iscalculated by using the relation given in Carey[11]basedon the equation suggested by Silver and Simpson[12].The heat transfer coefficient h c between the cover plateand the vaporflow is calculated by using a correlation sug-3214T.Kaya,J.Goldak/International Journal of Heat and Mass Transfer49(2006)3211–3220gested by Sleicher and Rouse[13]for fully developedflows in round ducts.It is extremely difficult to experimentally determine the heat transfer coefficient h c and a three-dimen-sional model is necessary to solve the vaporflow in the grooves.A convective boundary condition is more realistic since the use of temperature boundary condition implies h c!1.The convective boundary condition here with a reasonable heat transfer coefficient also allows some heat flux through the groove rather than assuming the entire heat load is transferred to the wick through thefin.3.Numerical procedureThe governing equations and associated boundary con-ditions described previously are solved by using the Galer-kinfinite element method.The computational domain under consideration is discretized with isoparametric and quadratic triangular elements.The numerical solution sequence for the all-liquid wick is straightforward.As the entire process is driven by the liquid evaporation at the vapor–liquid front,the energy equation isfirst solved.The numerical solution sequence is as follows:1.Initialize the problem by solving the energy equationassuming zero velocity inside the wick.2.Calculate the normal component of the outflow velocityat the interface between the wick and groove from the results of the energy equation,which is then used as an outflow boundary condition for the Darcy solver.3.Solve the Darcy equation to obtain the liquid velocityfield inside the wick.4.Solve the energy equation on the entire domain with theDarcy velocities.5.Return to step2until all equations and boundary con-ditions are satisfied to a desired level of accuracy.At high heat loads,when a separate vapor region devel-ops in the wick,the numerical procedure is more compli-cated since the location of the interface is also an unknown of the problem.Therefore,a more involved iter-ative scheme is necessary.The numerical solution proce-dure is summarized as follows:1.Initialize the problem by solving the Laplace equationfor temperature($2T=0)on the entire domain for a liquid-saturated wick.2.Choose an arbitrary temperature isoline close to thefinas the initial guess for the location of the vapor–liquid interface.3.Solve the energy equation for two separate domains:casing-vapor region and liquid region.Calculate the normal conductive heatflux at the vapor–liquid interface.4.Solve the Darcy equation separately in the vapor andliquid regions to calculate the vapor and liquid velocities inside the wick.5.Solve the energy equation with the Darcy velocities onthe entire domain by imposing the energy conservation boundary condition at the interface.6.Check if the temperature condition at the interface issatisfied.If it is not satisfied,the interface shape needs to be modified.7.Return to step3until all equations and boundary con-ditions are satisfied to obtain a preset level of accuracy.After each interface update at step6,the solution domain needs to be remeshed.As the transient terms are not maintained in the governing equations,the numerical procedure presented is not a moving boundary technique and only the converged solutions have a physical meaning. For each solution,the static pressure drop across the inter-face is calculated to make sure that the difference in pres-sures is less than the maximum available capillary pressure in the wick(P vÀP l62r/r p),where the normal viscous stress discontinuity and inertial forces are neglected.Thus,the momentum jump condition across the interface is satisfied as long as the maximum capillary pressure is not exceeded.The accommodation coefficient for all the calculations is assumed to be0.1,leading to a typical value of h i=3.32·106W mÀ2KÀ1.To test the influence of this parameter,the results are also obtained with the accommodation coeffi-cients of0.01and1.Since the resulting interfacial heat transfer coefficients are sufficiently large,the change in the maximum temperature is negligibly small,on the order of less than0.01%.A typical value for the convection heat transfer coefficient h c is100W mÀ2KÀ1.The change of h c from100to50results in an increase of less than3%in the cover plate maximum temperature.However,the over-all change in the wick temperatures is negligibly small. 4.Results and discussionNumerical calculations are performed for the evaporator section with an outer diameter of25.4·10À3m as shown in Fig.3.The porous wick inside the evaporator has an outer diameter of21.9·10À3and a thickness of7.24·10À3m. The wick permeability and porosity are K=4·10À14m2 and u=60%,respectively.The workingfluid is ammonia. The LHP saturation temperature and pressure difference on both sides of the wick are calculated by using a one-dimensional mathematical model.The model is based on the steady-state energy conservation equations and the pressure drop calculations along thefluid path inside the LHP.The details of this mathematical model are presented in[14].Fig.4represents the calculated saturation temper-ature and pressure drop values across the wick as a function of the applied power.The pressure drops and heat transfer coefficients in the two-phase regions of the LHP are calculated by using the interfacial shear model of Chen [15].Incompressible fully developedfluidflow relations are used to calculate the pressure drop for the single phase regions.T.Kaya,J.Goldak/International Journal of Heat and Mass Transfer49(2006)3211–32203215Fig.5represents the temperature field and liquid velo-city vectors when the wick is completely saturated by liquid at Q in =100W.The solution is obtained by solving the mass conservation,Darcy and energy equations.At this heat load,by using the one-dimensional mathematical model,it is calculated that T sat =7.81°C and D p =247Pa.As the vapor flows along the grooves,it becomes super-heated due to the heating from the wall.Without solving the vapor flow in the grooves using a three-dimensional model,it is not possible to calculate the vapor temperature in the grooves.Accurate experimental measurements are also difficult although a range for the vapor superheat can be deduced based on the wall-temperature measure-ments.In our calculations,the vapor in the grooves is assumed to be superheated by 3°C.A similar approach is also used in Figus et al.[4].Thus,T v =10.81°C and other related parameters for the calculations are as follows:P v =569784.8Pa,P core =569537.8Pa,q in =1254W m À2,k c =k w =14.5W m À1K À1,k eff=6.073W m À1K À1,and h i =2.733·106W m À2K À1.The thermal properties of ammonia are calculated at the saturation temperature for a given applied power using the relations in [16].In the numerical calculations,the thermal properties are assumed constant for a given saturation temperature.It can be seen from Fig.5that,the working fluid evap-orates at the wick interface under the applied heat load.The liquid flows from the evaporator core into the wick and turns toward the interface under the fin.The heat flux along the fin–wick interface is not constant and varies around 2000W m À2.In comparison,in the previously ref-erenced works [3–5]with an exception in [6],an estimated constant heat flux is directly applied at the fin–wick surface and the temperature drop across the casing is ignored due to the low thermal resistance.Applying the heat at the cas-ing allows the calculation of the temperature distribution at the casing surface.At low heat loads,the liquid velocity is relatively small as well as the corresponding Peclet number (Pe =q in L w c pl /h fg k eff).For example,at Q in =100W,Pe is on the order of 10À2.Therefore,the contribution from the convective terms could be neglected.Therefore,in the earlier solutions [3–5],the Laplace equation for the temper-ature is solved instead of the full energy equation.With this assumption,the Darcy and energy equations are also decoupled,which significantly simplifies the solution algo-rithm.However,at higher heat loads,the convective terms need to be taken account as is done in [6].In our study,we keep the convective terms in the governing equations and solve together the mass conservation,Darcy and energy equations as a coupled problem.The determination of the effective thermal conductivity of the wick k effis not trivial as it depends in a complex man-ner on the geometry of the porous medium.The solution on Fig.5is obtained by assuming that there is no heat transfer between the solid porous matrix and fluid (heat transfer in parallel).This is a well-known correlation obtained by the weighted arithmetic mean of k l and k w (k eff=u k l +(1Àu )k w ),where u is the wick porosity.A number of relations for the prediction of k effis proposed in the literature.To inves-tigate the effect of k effon the results,the same problem is solved for the all-liquid wick case by using six different correlations in addition to the weighted arithmetic mean.These are weighted harmonic (heat transfer in series)and geometric means of the thermal conductivities of k w and k l ,and other relations developed by Maxwell [17],Krupiczka [18],Zehner and Schlunder [19],and Alexander [20].Fig.6represents the results obtained by using the dif-ferent k effvalues at an arbitrarily chosen location of h =80°.The change of slope indicates the wick and fin interface.The temperature profiles directly depend on k eff.The series and parallel arrangements represent the highest and lowest con-ductivities,respectively.The other relations are intermedi-ate between these two.One specific difficulty is that the correlations produce significantly different values when the thermal conductivities of the porous medium and fluid are greatly different from each other as previously studied in Nield [21].As an example,the ratio of the thermal con-ductivities for the liquid and vapor regions of the wick at T sat =7.81°C are k l /k w =0.0362and k v /k w =0.0017,respectively.There is therefore further difficulty when both of the phases are present inside the wick.A given relation for k effwill not have the same accuracy for the liquidandFig.5.Velocity vectors and temperature field at Q in =100W.3216T.Kaya,J.Goldak /International Journal of Heat and Mass Transfer 49(2006)3211–3220vapor regions.The effective thermal conductivity k effobtained by using different relations are given in Table1. The results vary significantly.There is clearly a need for experimental data for an accurate determination of k eff.In the lack of experimental data,we use the parallel arrange-ment for the rest of the numerical calculations.This is also used in several previous works[3–6].As shown in Fig.6,the different values of k efflead to the qualitatively similar tem-perature profiles.The largest temperature difference between the parallel and series solutions was within0.5K. The difference in temperature is small because of the low Peclet numbers.At the lower limit,as Pe!0,the energy equation reduces to the Laplace equation for temperature and the influence of k effon the temperature distribution is primarily through theflux boundary conditions.Note that the temperature at the core is imposed as a boundary con-dition and it has the same value for all cases.The adequacy of using the Darcy’s law for describing theflow inside the wick is also considered.For example, Cao and Faghri[6]use an expression from analogy with Navier–Stokes equation for theflow inside the porous medium,which takes into account the convective(uÆ$)u/u and viscous transport l/u($2u)terms in addition to the Darcy’s law.Especially at high heatflux rates,the non-Darcyflow behavior could be important.Beck[22]has showed that the inclusion of the convection term in the Darcy equation may lead to an under or over specified sys-tem of equations.Similar conclusions have also been reported in[23].For these reasons,we do not take into account the convective terms.The maximum Reynolds number based on the effective pore diameter of the wick of2.4l m is on the order of10À2,which occurs in the vapor region near thefin edge.Therefore,the quadratic inertia terms are negligible in both the vapor and liquid regions.A comparison of the results obtained from Darcy and Brinkman equations for the all-liquid wick case showed that contribution from the Brinkman terms can also be safely neglected.As a result,the non-Darcyflow effects could be ignored without penalty.At sufficiently high heatflux values,it is expected that the nucleation will start at the microscopic cavities at the fin–wick interface.The boiling can initiate at small super-heat values as a result of trapped gas in these cavities. The vapor bubbles formed at thefin–wick interface unite and lead to the formation of a vapor–liquid interface inside the wick as originally suggested by Demidov and Yatsenko [3].With increasing heatflux,the vapor–liquid interface recedes further into the wick because of the increased evap-oration and insufficient supply of the returning subcooled liquid.Thus,the vapor zone under thefin continues to grow in size and starts connecting with the vapor grooves. For a given heat load,there exists a steady-state solution for which the heat transferred to the wick from thefin sur-face is balanced by the convective heat output to the vapor–groove interface where the evaporation takes place. As the applied heat load is increased,the vapor region under thefin grows.For sufficiently large applied heat load,no converged solution is possible unless the removal of vapor from the interface inside the wick is allowed from the wick–groove interface.For the transition from the all-liquid wick to the vapor–liquid wick,a boiling incipient superheat value is assumed. It is difficult to predict the incipient superheat,which depends on several parameters in a complex manner.In our calculations,when the liquid temperature under the fin is4°C higher than T sat,it is assumed that a vapor region will form under thefin.Then,a new solution is obtained by using the numerical procedure outlined for the vapor–liquid wick.These results provide a reference base for the boiling analysis of the LHP using nuclear clus-ter theory,which will be addressed later in the paper.Fig.7 represents the results obtained at a heat load of Q in=300W.The LHP saturation temperature and pres-sure drop is T sat=11.03°C and D p=2181Pa,respec-tively.It should be noted that the one-dimensional model does not take into account the presence of a vapor region inside the wick.The change on the wick effective thermal conductivity in the presence of vapor zone needs to be esti-mated to improve the calculations of the boundary condi-tions from the one-dimensional model.An iterative procedure between the one-and two-dimensional models could be more representative.However,this would be com-putationally intensive and no significant change in the overall results is expected.Other required numerical valuesTable1The effective thermal conductivity values for the liquid and vapor regions computed from different correlationsRelation k l(W mÀ1KÀ1)k v(W mÀ1KÀ1) Harmonic mean(series arrangement)0.8490.040Alexandre[20] 1.0540.093Zehner and Schlunder[19] 1.6840.148Krupiczka[18] 1.7560.152Geometric mean 1.9670.309Maxwell[17] 4.845 4.450Arithmetic mean (parallel arrangement)6.073 5.774T.Kaya,J.Goldak/International Journal of Heat and Mass Transfer49(2006)3211–32203217。

NADCAP热处理工作组审核手册

NADCAP热处理工作组审核手册发布日期:2013年10月28日目录1. 范围 (3)2. 定义 (3)3. 依据Nadcap审核准则实施审核 (6)4. 质量体系................................................................................. .. (21)5. 热处理审核清单 (23)6. 附录A:铝合金标准 (32)7. 附录B: 主要客户标准矩阵图 (38)8. 附录C: 主要客户对保温开始时间与保温时间的定义 (42)9. 附录D: 审核员公告 (43)9.1 主题 (43)9.1.1 高温测量 (43)9.1.2 工序控制测试 (51)9.1.3 主要客户的关注内容 (53)9.1.3.1审核员公告11-001 HT–罗罗公司 (53)9.1.4 作业审核 (59)9.1.5 审核管理 (60)9.1.6 其他内容 (64)10. 附录E:文件修订历史 (65)1. 范围This Audit Handbook has been prepared to assist the Heat Treating Auditors and Suppliers as follows:此审核手册在以下方面给予热处理审核员和供应商相应指导:1.Where necessary, provide clarification on the intent and rationale of the Heat TreatingTask Group as it pertains to specific questions contained in the current revision of AC7102 and its corresponding slash sheets, as well as heat treating specific issues with Materials Testing Laboratories (MTL) checklists commonly used during heat treating audits.需要时,说明热处理工作组的目的和理论依据。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

Two-Dimensional Steady Heat Conduction5-43C For a medium in which the finite difference formulation of a general interior node is given in itssimplest form as T T T T T gl k left top right bottom node node +++-+=402 :(a ) Heat transfer is steady, (b ) heat transfer is two-dimensional, (c ) there is heat generation in the medium, (d ) the nodal spacing is constant, and (e ) the thermal conductivity of the medium is constant.5-44C For a medium in which the finite difference formulation of a general interior node is given in its simplest form as T T T T T node left top right bottom =+++()/4:(a ) Heat transfer is steady, (b ) heat transfer is two-dimensional, (c ) there is no heat generation in the medium, (d ) the nodal spacing is constant, and (e ) the thermal conductivity of the medium is constant.5-45C A region that cannot be filled with simple volume elements such as strips for a plane wall, and rectangular elements for two-dimensional conduction is said to have irregular boundaries . A practical way of dealing with such geometries in the finite difference method is to replace the elements bordering the irregular geometry by a series of simple volume elements.temperatures and the rate of heat loss from the bottom surface through a 1-m long section are to be determined.Assumptions 1 Heat transfer through the body is given to be steady and two-dimensional. 2 Heat is generated uniformly in the body. 3 Radiation heat transfer is negligible. Properties The thermal conductivity is given to be k = 45 W/m ⋅°C.Analysis The nodal spacing is given to be ∆x =∆x =l =0.05 m, and the general finite difference form of an interior node for steady two-dimensional heat conduction is expressed asT T T T T gl kleft top right bottom node node +++-+=402whereC 5.93CW/m 214)m 05.0)( W/m 108(236202node ︒=︒⋅⨯==k l g k l g The finite difference equations for boundary nodes are obtained by applying an energy balance on the volume elements and taking the direction of all heat transfers to be towards the node under consideration: 04-200240290260 :(interior) 3 Node 04-290325290350 :(interior) 2 Node 02)(32522902402 :)convection ( 1 Node 203202201111=++++=++++=+-+-+-+-∞kl gT k l gT kl g T T hl l T l k l T kl l T l k whereC 20 , W/m 108 C,. W/m 50 C, W/m .45362︒=⨯=︒=︒=∞T gh k Substituting,T 1 = 280.9°C , T 2 = 397.1°C , T 3 = 330.8°C ,(b ) The rate of heat loss from the bottom surface through a 1-m long section isW1808=︒+++⨯︒⋅=-+-+-+-=-==∞∞∞∞∞∑∑C 20)/2]-(32520)-(280.920)-(24020)/2-m)[(200 1 m C)(0.05 W/m 50()325)(2/()()240()200)(2/()(21surface, element,T l h T T hl T hl T l h T T hAQQmm m mmh, T ∞ConvectionAssumptions 1 Heat transfer through the body is given to be steady and two-dimensional. 2 There is no heat generation in the body.Properties The thermal conductivity is given to be k = 45 W/m ⋅°C.Analysis The nodal spacing is given to be ∆x =∆x =l =0.01 m, and the general finite difference form of an interior node for steady two-dimensional heat conduction for the case of no heat generation is expressed as4/)( 04bottom right top left node 2node nodebottom right top left T T T T T kl g T T T T T +++=→=+-+++There is symmetry about the horizontal, vertical, and diagonal lines passing through the midpoint, and thus we need to consider only 1/8th of the region. Then,86429731T T T T T T T T ======Therefore, there are there are only 3 unknown nodaltemperatures, 531 and ,,T T T , and thus we need only 3 equations to determine them uniquely. Also, we can replace the symmetry lines by insulation and utilize the mirror-image concept when writing the finite difference equations for the interior nodes.225152214/4 :(interior) 3 Node 4/)2200( :(interior) 2 Node 4/)2180180( :(interior) 1 Node T T T T T T T T ==++=++= Solving the equations above simultaneously givesC190C 185︒=====︒====865429731T T T T T T T T TDiscussion Note that taking advantage of symmetry simplified the problem greatly.180200180 150 180 200 180 150150 180 200 180 150180 200 180Assumptions 1 Heat transfer through the body is given to be steady and two-dimensional. 2 There is no heat generation in the body.Properties The thermal conductivity is given to be k = 20 W/m ⋅°C.Analysis The nodal spacing is given to be ∆x =∆x =l =0.02 m, and the general finite difference form of an interior node for steady two-dimensional heat conduction for the case of no heat generation is expressed as4/)( 04bottom right top left node 2node nodebottom right top left T T T T T kl g T T T T T +++=→=+-+++(a ) There is symmetry about the insulated surfaces as well as about the diagonal line. Therefore 23T T =, and 421 and ,,T T T are the only 3 unknown nodal temperatures. Thus we need only 3 equations to determine them uniquely. Also, we can replace the symmetry lines by insulation and utilize the mirror-image concept when writing the finite difference equations for the interior nodes.4/)22( :(interior) 4 Node 4/)2200( :(interior) 2 Node 4/)180180( :(interior) 1 Node 324142321T T T T T T T T T +=++=+++=Also,23T T =Solving the equations above simultaneously givesC185C 190︒=︒===1432T T T T(b ) There is symmetry about the insulated surface as well as the diagonal line. Replacing the symmetry linesby insulation, and utilizing the mirror-image concept, the finite difference equations for the interior nodes can be written as4/)21402( :(interior) 4 Node 4/)12140( :(interior) 3 Node 4/)120120( :(interior) 2 Node 4/)120120( :(interior) 1 Node 324243142321T T T T T T T T T T T T T ++==++=+++=+++=Solving the equations above simultaneously givesC128.6C 122.9︒==︒==4321T T T TDiscussion Note that taking advantage of symmetry simplified the problem greatly.180 200 150 180 200100120140 100120 1405-49 Starting with an energy balance on a volume element, the steady two-dimensional finite difference equation for a general interior node in rectangular coordinates for T (x , y ) for the case of variable thermal conductivity and uniform heat generation is to be obtained.Analysis We consider a volume element of size ∆∆x y ⨯⨯1 centered about a general interior node (m , n ) ina region in which heat is generated at a constant rate of gand the thermal conductivity k is variable (see Fig. 5-24 in the text). Assuming the direction of heat conduction to be towards the node under consideration at all surfaces, the energy balance on the volume element can be expressed asQ Q Q Q G E tcond, left cond, top cond, right cond, bottom element element++++==∆∆0 for the steady case. Again assuming the temperatures between the adjacent nodes to vary linearly and noting that the heat transfer area is ∆y ⨯1 in the x direction and ∆x ⨯1 in the y direction, the energy balance relation above becomes0)1()1(+ )1()1(+)1(0,1,,,,1,,1,,,,1,=⨯∆⨯∆+∆-⨯∆∆-⨯∆+∆-⨯∆∆-⨯∆-++-y x gyT T x k xT T y k y T T x k xT T y k nm n m n m nm n m n m nm n m n m nm n m n mDividing each term by ∆∆x y ⨯⨯1 and simplifying gives022,021,,1,2,1,,1=+∆+-+∆+-+-+-nm n m n m n m nm n m n m k gy T T T x T T T For a square mesh with ∆x = ∆y = l, and the relation above simplifies to04,20,1,1,,1,1=+-+++--+-nm nm n m n m n m n m k l g T T T T T It can also be expressed in the following easy-to-remember form:04node20nodebottom right top left =+-+++k l g T T T T T5-50 A long solid body is subjected to steady two-dimensional heat transfer. The unknown nodal temperatures and the rate of heat loss from the top surface are to be determined.Assumptions 1 Heat transfer through the body is given to be steady and two-dimensional. 2 Heat is generated uniformly in the body.Properties The thermal conductivity is given to be k = 180 W/m ⋅°C.Analysis (a ) The nodal spacing is given to be ∆x =∆x =l =0.1 m, and the general finite difference form of an interior node equation for steady two-dimensional heat conduction for the case of constant heat generation is expressed as042node nodebottom right top left =+-+++kl g T T T T TThere is symmetry about a vertical line passing through the middle of the region, and thus we need to consider only half of the region. Then,4321 and T T T T ==Therefore, there are there are only 2 unknown nodal temperatures, T 1 and T 3, and thus we need only 2 equations to determine them uniquely. Also, we can replace the symmetry lines by insulation and utilize the mirror-image concept when writing the finite difference equations for the interior nodes.4200150 :(interior) 3 Node 04120100 :(interior) 1 Node 23412132=+-+++=+-+++kl gT T T kl gT T TNoting that 4321 and T T T T == and s ubstituting,C W/m 180m))(0.1 W/m 10(33500CW/m 180m))(0.1 W/m 10(32202373123713=︒⋅+-+=︒⋅+-+T T T TThe solution of the above system isC439.0C 411.5︒==︒==4321T T T T(b ) The total rate of heat transfer from the top surface topQ can be determined from an energy balance on a volume element at the top surface whose height is l /2, length 0.3 m, and depth 1 m:depth)(per C 100)-m)(411.5 1(C )100120(2m 1)C W/m 180(2m )2/1.03.0)( W/m 10(010012100120212)2/13.0(337top10top m Q l T kl l l k l g Q W 265,750=⎪⎭⎫ ⎝⎛︒+︒-︒⋅-⨯-==⎪⎪⎭⎫ ⎝⎛-⨯+-⨯+⨯⨯+100 100 100 100 120150 200 200 200 200120 1505-51"!PROBLEM 5-51""GIVEN"k=180 "[W/m-C], parameter to be varied"g_dot=1E7 "[W/m^3], parameter to be varied"DELTAx=0.10 "[m]"DELTAy=0.10 "[m]"d=1 "[m], depth""Temperatures at the selected nodes are also given in the figure""ANALYSIS""(a)"l=DELTAxT_1=T_2 "due to symmetry"T_3=T_4 "due to symmetry""Using the finite difference method, the two equations for the two unknown temperatures are determined to be"120+120+T_2+T_3-4*T_1+(g_dot*l^2)/k=0150+200+T_1+T_4-4*T_3+(g_dot*l^2)/k=0"(b)""The rate of heat loss from the top surface can be determined from an energy balance on a volume element whose height is l/2, length 3*l, and depth d=1 m"-Q_dot_top+g_dot*(3*l*d*l/2)+2*(k*(l*d)/2*(120-100)/l+k*l*d*(T_1-100)/l)=05-52 A long solid body is subjected to steady two-dimensional heat transfer. The unknown nodal temperatures are to be determined.Assumptions 1 Heat transfer through the body is given to be steady and two-dimensional. 2 There is no heat generation in the body.Properties The thermal conductivity is given to be k = 20 W/m ⋅°C.Analysis The nodal spacing is given to be ∆x =∆x =l =0.01m, and the general finite difference form of an interior node for steady two-dimensional heat conduction for the case of no heat generation is expressed as04 04node bottom right top left 2node nodebottom right top left =-+++→=+-+++T T T T T kl g T T T T T(a ) There is symmetry about a vertical line passing through the nodes 1 and 3. Therefore, 23T T =, 46T T =, and 5421 and ,,,T T T T are the only 4 unknown nodal temperatures, and thus we need only 4 equations to determine them uniquely. Also, we can replace the symmetry lines by insulation and utilize the mirror-image concept when writing the finite difference equations for the interior nodes.44 :(interior) 5 Node4250250 :(interior) 4 Node 4200 :(interior) 2 Node 042150150 :(interior) 1 Node 52452245112=--+++=-+++=-++T T T T T T T T T T T C200C 225C 200C 175︒=︒==︒==︒=564321T T T T T T(b ) There is symmetry about a vertical line passing through the middle. Therefore, 23T T =and 14T T =. Replacing the symmetry lines by insulation and utilizing the mirror-image concept, the finite difference equations for the interior nodes 1 and 2 are determined to be 4200100100 :(interior) 2 Node 042200100 :(interior) 1 Node 4112-+++=-++T T T T Solving the 2 equations abovesimultaneously givesC136C 143︒==︒==3241T T T TDiscussion Note that taking advantage of symmetry simplified the problem greatly.100300100 100 100 100200 200 200 200 2005-53 Heat conduction through a long L-shaped solid bar with specified boundary conditions is considered. The unknown nodal temperatures are to be determined with the finite difference method. √Assumptions 1 Heat transfer through the body is given to be steady and two-dimensional. 2 Thermal conductivity is constant. 3 Heat generation is uniform. Properties The thermal conductivity is given to be k = 45 W/m ⋅°C. Analysis (a ) The nodal spacing is given to be ∆x =∆x =l =0.015 m, and the general finite difference form of an interior node for steady two-dimensional heat conduction for the case of constant heat generation is expressed as0420nodebottom right top left =+-+++kl gT T T T T We observe that all nodes are boundary nodes except node 5 that is an interior node. Therefore, we will have to rely on energy balances to obtain the finite difference equations. Using energy balances, the finite difference equations for each of the 8 nodes are obtained as follows:Node 1: 0422)(222014121=+-+-+-+∞l g l T T l k l T T l k T T l h l q L Node 2: 0222)(202523212=+-+-+-+-∞l g l T T kl l T T l k l T T l kT T hl Node 3: 0422)(2036323=+-+-+-∞l g l T T l k l T T l kT T hl Node 4: 02120222045441=+-+-+-+l g l T T kl l T l k l T T l k l qL Node 5: 0412*******=+-+++kl gT T T T Node 6: 04321202)(2067665636=+-+-+-+-+-∞l g l T T l k l T kl l T T kl l T T l kT T hl Node 7: 0212022)(20778767=+-+-+-+-∞l g l T kl l T T l k l T T l kT T hl Node 8: 0412022)(2208878=+-+-+-∞l g l T l k l T T l kT T l h where , W/m 8000 , W/m 1052360=⨯=L q g l = 0.015 m, k = 45 W/m ⋅︒C, h = 55 W/m 2⋅︒C, and T ∞ =30︒C.This system of 8 equations with 8 unknowns is the finite difference formulation of the problem.(b ) The 8 nodal temperatures under steady conditions are determined by solving the 8 equations above simultaneously with an equation solver to beT 1 =163.6︒C, T 2 =160.5︒C, T 3 =156.4︒C, T 4 =154.0︒C, T 5 =151.0︒C, T 6 =144.4︒C, T 7 =134.5︒C, T 8 =132.6︒CDiscussion The accuracy of the solution can be improved by using more nodal points.5-54E A long solid bar is subjected to steady two-dimensional heat transfer. The unknown nodal temperatures and the rate of heat loss from the bar through a 1-ft long section are to be determined. Assumptions 1 Heat transfer through the body is given to be steady and two-dimensional. 2 Heat is generated uniformly in the body. 3 The heat transfer coefficient also includes the radiation effects. Properties The thermal conductivity is given to be k = 16 Btu/h.ft ⋅°C.Analysis The nodal spacing is given to be ∆x =∆x =l =0.2 ft, and the general finite difference form of an interior node for steady two-dimensional heat conduction is expressed asT T T T T gl kleft top rightbottom node node +++-+=402(a ) There is symmetry about the vertical, horizontal, and diagonal lines passing through the center. Therefore, 9731T T T T ===and8642T T T T ===, and 521 and ,,T T T are the only 3 unknown nodal temperatures, and thus we need only 3 equations to determine them uniquely. Also, we can replace the symmetry lines by insulation and utilize the mirror-image concept for the interior nodes.The finite difference equations for boundary nodes are obtained by applying an energy balance onthe volume elements and taking the direction of all heat transfers to be towards the node under consideration:44 :(interior) 5 Node 02)( 22 :)convection ( 2 Node 04)(2222 :)convection ( 1 Node 2052202252120112=+-=+-+-+-=+-+-∞∞kl gT T l g T T hl l T T kl l T T l k l g T T lh l T T l kwhere 350f t Btu/h 1019.0⋅⨯=g, l = 0.2 ft, k = 16 Btu/h.ft ⋅︒F, h =7.9 Btu/h.ft 2⋅︒F, and T ∞ =70︒F. The 3 nodal temperatures under steady conditions are determined by solving the 3 equations above simultaneously with an equation solver to be9731T T T T ====304.85︒F,8642T T T T ====316.16︒F, T 5 =328.04︒F(b ) The rate of heat loss from the bar through a 1-ft long section is determined from an energy balance on one-eight section of the bar, and multiplying the result by 8:[]flength)ft per (F 70]2-16.316ft)[304.85 ft)(1 F)(0.2/2ft Btu/h 9.7(8ft) 1()228ft) 1()(2)(28822121conv section,eight one Btu/h 3040=︒⨯+︒⋅⋅=-+⨯=⎥⎦⎤⎢⎣⎡-+-⨯=⨯=∞∞∞-T T T l h T T l h T T l h Q QDiscussion Under steady conditions, the rate of heat loss from the bar is equal to the rate of heat generationwithin the bar per unit length, and is determined to belength)f t per (Btu/h 3040f t) 1 f t 0.4 f t 4.0)(Btu/h.f t 1019.0(350gen =⨯⨯⨯===V gE Qwhich confirms the results obtained by the finite difference method.h , T ∞h , T ∞ h , T ∞ h , T ∞5-55 Heat transfer through a square chimney is considered. The nodal temperatures and the rate of heat loss per unit length are to be determined with the finite difference method.Assumptions 1 Heat transfer is given to be steady and two-dimensional since the height of the chimney is large relative to its cross-section, and thus heat conduction through the chimney in the axial direction is negligible. It is tempting to simplify the problem further by considering heat transfer in each wall to be one dimensional which would be the case if the walls were thin and thus the corner effects were negligible. This assumption cannot be justified in this case since the walls are very thick and the corner sections constitute a considerable portion of the chimney structure. 2 There is no heat generation in the chimney. 3 Thermal conductivity is constant.Properties The thermal conductivity and emissivity are given to be k = 1.4 W/m ⋅°C and ε = 0.9. Analysis (a ) The most striking aspect of this problem is the apparent symmetry about the horizontal and vertical lines passing through the midpoint of the chimney. Therefore, we need to consider only one-fourth of the geometry in the solution whose nodal network consists of 10 equally spaced nodes. No heat can cross a symmetry line, and thus symmetry lines can be treated as insulated surfaces and thus “mirrors” in the finite -difference formulation. Considering a unit depth and using the energy balance approach for the boundary nodes (again assuming all heat transfer to be into the volume element for convenience), the finite difference formulation is obtained to beNode 1: 0])273([222)(24141512100=+-+-+-+-T T l l T T l k l T T l k T T l h surr εσ Node 2: 0])273([22)(424262321200=+-+-+-+-+-T T l l T T kl l T T l k l T T l k T T l h surr εσNode 3: 0])273([22)(434373432300=+-+-+-+-+-T T l lT T kl l T T l k l T T l k T T l h surr εσNode 4: 0])273([22)(4444843400=+-εσ+-+-+-T T l lT T l k l T T l k T T l h surrNode 5: 022)(251565=-+-+-lT T l k l T T l k T T l h i iNode 6: 022)(6267656=-+-+-+-lT T kl l T T l k l T T l k T T l h i iNode 7: 022)(787379767=-+-+-+-+-lT T kl l T T kl l T T l k l T T l k T T l h i iNode 8: 0])273([22)(4848781084800=+-+-+-+-+-T T l lT T kl l T T l k l T T l k T T l h surr εσNode 9: 022)(2910979=-+-+-lT T l k l T T l k T T l h i iNode 10: 0])273([222)(241041091081000=+-+-+-+-T T l l T T l k l T T l k T T l h surr εσ1 2 3 4 T surrh o , T owhere l = 0.1 m, k = 1.4 W/m ⋅︒C, h i = 75 W/m 2⋅︒C, T i =280︒C, h o = 18 W/m 2⋅︒C, T 0 =15︒C, T surr =250 K, ε = 0.9, and σ = 5.67⨯10-8 W/m 2.K 4. This system of 10 equations with 10 unknowns constitutes the finite difference formulation of the problem.(b ) The 10 nodal temperatures under steady conditions are determined by solving the 10 equations above simultaneously with an equation solver to beT 1 =94.5︒C, T 2 =93.0︒C, T 3 =82.1︒C, T 4 =36.1︒C, T 5 =250.6︒C, T 6 =249.2︒C, T 7 =229.7︒C, T 8 =82.3︒C, T 9 =261.5︒C, T 10 =94.6︒C(c ) The rate of heat loss through a 1-m long section of the chimney is determined fromW3153=︒+++⨯︒⋅=-+-+-+-=-===∑∑∑C 261.5)/2]-(280229.7)-(280249.2)-(280250.6)/2-m)[(280 1 m C)(0.1 W/m 75(4)])(2/()()())(2/([4)(44429765surface,surfaceinner element,chimneyof fourth -one T T l h T T l h T T l h T T l h T T Ah QQQi i i i i i i i mm i m iDiscussion The rate of heat transfer can also be determined by calculating the heat loss from the outersurface by convection and radiation.5-56 Heat transfer through a square chimney is considered. The nodal temperatures and the rate of heat loss per unit length are to be determined with the finite difference method.Assumptions 1 Heat transfer is given to be steady and two-dimensional since the height of the chimney is large relative to its cross-section, and thus heat conduction through the chimney in the axial direction is negligible. It is tempting to simplify the problem further by considering heat transfer in each wall to be one dimensional which would be the case if the walls were thin and thus the corner effects were negligible. This assumption cannot be justified in this case since the walls are very thick and the corner sections constitute a considerable portion of the chimney structure. 2 There is no heat generation in the chimney. 3 Thermal conductivity is constant. 4 Radiation heat transfer is negligible. Properties The thermal conductivity of chimney is given to be k = 1.4 W/m ⋅°C. Analysis (a ) The most striking aspect of this problem is the apparent symmetry about the horizontal and vertical lines passing through the midpoint of the chimney. Therefore, we need to consider only one-fourth of the geometry in the solution whose nodal network consists of 10 equally spaced nodes. No heat can cross a symmetry line, and thus symmetry lines can be treated as insulated surfaces and thus “mirrors” in the finite -difference formulation. Considering a unit depth and using the energy balance approach for the boundary nodes (again assuming all heat transfer to be into the volume element for convenience), the finite difference formulation is obtained to beNode 1: 022)(21512100=-+-+-lT T l k l T T l k T T l h Node 2: 022)(262321200=-+-+-+-l T T kl l T T l k l T T l k T T l hNode 3: 022)(373432300=-+-+-+-l T T kl l T T l k l T T l k T T l hNode 4: 022)(4843400=-+-+-lT T l k l T T l k T T l hNode 5: 022)(251565=-+-+-lT T l k l T T l k T T l h i iNode 6: 022)(6267656=-+-+-+-lT T kl l T T l k l T T l k T T l h i iNode 7: 022)(787379767=-+-+-+-+-l T T kl l T T kl l T T l k l T T l k T T l h i iNode 8: 022)(8781084800=-+-+-+-lT T kl l T T l k l T T l k T T l hNode 9: 022)(2910979=-+-+-lT T l k l T T l k T T l h i iNode 10: 022)(21091081000=-+-+-lTT l k l T T l k T T l h1 2 3 4h o , T owhere l = 0.1 m, k = 1.4 W/m ⋅︒C, h i = 75 W/m 2⋅︒C, T i =280︒C, h o = 18 W/m 2⋅︒C, T 0 =15︒C, and σ = 5.67⨯10-8 W/m 2.K 4. This system of 10 equations with 10 unknowns constitutes the finite difference formulation of the problem.(b ) The 10 nodal temperatures under steady conditions are determined by solving the 10 equations above simultaneously with an equation solver to beT 1 =118.8︒C, T 2 =116.7︒C, T 3 =103.4︒C, T 4 =53.7︒C, T 5 =254.4︒C, T 6 =253.0︒C, T 7 =235.2︒C, T 8 =103.5︒C, T 9 =263.7︒C, T 10 =117.6︒C(c ) The rate of heat loss through a 1-m long section of the chimney is determined fromW2783=︒+++⨯︒⋅=-+-+-+-=-===∑∑∑C 263.7)/2]-(280235.2)-(280253.0)-(280254.4)/2-m)[(280 1 m C)(0.1 W/m 75(4)])(2/()()())(2/([4)(44429765surface,surfaceinner element,chimneyof fourth -one T T l h T T l h T T l h T T l h T T Ah QQQi i i i i i i i mm i m iDiscussion The rate of heat transfer can also be determined by calculating the heat loss from the outersurface by convection.5-57"!PROBLEM 5-57""GIVEN"k=1.4 "[W/m-C]"A_flow=0.20*0.40 "[m^2]"t=0.10 "[m]"T_i=280 "[C], parameter to ve varied"h_i=75 "[W/m^2-C]"T_o=15 "[C]"h_o=18 "[W/m^2-C]"epsilon=0.9 "parameter to ve varied"T_sky=250 "[K]"DELTAx=0.10 "[m]"DELTAy=0.10 "[m]"d=1 "[m], unit depth is considered"sigma=5.67E-8 "[W/m^2-K^4], Stefan-Boltzmann constant""ANALYSIS""(b)"l=DELTAx"We consider only one-fourth of the geometry whose nodal network consists of 10 nodes. Using the finite difference method, 10 equations for 10 unknown temperatures are determined to be"h_o*l/2*(T_o-T_1)+k*l/2*(T_2-T_1)/l+k*l/2*(T_5-T_1)/l+epsilon*sigma*l/2*(T_sky^4-(T_1+273)^4)=0 "Node 1"h_o*l*(T_o-T_2)+k*l/2*(T_1-T_2)/l+k*l/2*(T_3-T_2)/l+k*l*(T_6-T_2)/l+epsilon*sigma*l*(T_sky^4-(T_2+273)^4)=0 "Node 2"h_o*l*(T_o-T_3)+k*l/2*(T_2-T_3)/l+k*l/2*(T_4-T_3)/l+k*l*(T_7-T_3)/l+epsilon*sigma*l*(T_sky^4-(T_3+273)^4)=0 "Node 3"h_o*l*(T_o-T_4)+k*l/2*(T_3-T_4)/l+k*l/2*(T_8-T_4)/l+epsilon*sigma*l*(T_sky^4-(T_4+273)^4)=0 "Node 4"h_i*l/2*(T_i-T_5)+k*l/2*(T_6-T_5)/l+k*l/2*(T_1-T_5)/l=0 "Node 5"h_i*l*(T_i-T_6)+k*l/2*(T_5-T_6)/l+k*l/2*(T_7-T_6)/l+k*l*(T_2-T_6)/l=0 "Node 6"h_i*l*(T_i-T_7)+k*l/2*(T_6-T_7)/l+k*l/2*(T_9-T_7)/l+k*l*(T_3-T_7)/l+k*l*(T_8-T_7)/l=0 "Node 7" h_o*l*(T_o-T_8)+k*l/2*(T_4-T_8)/l+k*l/2*(T_10-T_8)/l+k*l*(T_7-T_8)/l+epsilon*sigma*l*(T_sky^4-(T_8+273)^4)=0 "Node 8"h_i*l*(T_i-T_9)+k*l/2*(T_7-T_9)/l+k*l/2*(T_10-T_9)/l=0 "Node 9"h_o*l/2*(T_o-T_10)+k*l/2*(T_8-T_10)/l+k*l/2*(T_9-T_10)/l+epsilon*sigma*l/2*(T_sky^4-(T_10+273)^4)=0 "Node 10""Right top corner is considered. The locations of nodes are as follows:""Node 1: Middle of top surfaceNode 2: At the right side of node 1Node 3: At the right side of node 2Node 4: Corner nodeNode 5: The node below node 1, at the middle of inner top surfaceNode 6: The node below node 2Node 7: The node below node 3, at the inner cornerNode 8: The node below node 4Node 9: The node below node 7,at the middle of inner right surfaceNode 10: The node below node 8, at the middle of outer right surface"T_corner=T_4T_inner_middle=T_9"(c)""The rate of heat loss through a unit depth d=1 m of the chimney is"Q_dot=4*(h_i*l/2*d*(T_i-T_5)+h_i*l*d*(T_i-T_6)+h_i*l*d*(T_i-T_7)+h_i*l/2*d*(T_i-T_9))5-58 The exposed surface of a long concrete damn of triangular cross-section is subjected to solar heat flux and convection and radiation heat transfer. The vertical section of the damn is subjected to convection with water. The temperatures at the top, middle, and bottom of the exposed surface of the damn are to be determined.Assumptions 1 Heat transfer through the damn is given to be steady and two-dimensional. 2 There is no heat generation within the damn. 3 Heat transfer through the base is negligible. 4 Thermal properties and heat transfer coefficients are constant.Properties The thermal conductivity and solar absorptivity are given to be k = 0.6 W/m ⋅°C and αs = 0.7. Analysis The nodal spacing is given to be ∆x =∆x =l =1 m, and all nodes are boundary nodes. Node 5 on the insulated boundary can be treated as an interior node for which 04node bottom right top left =-+++T T T T T . Using the energy balance approach and taking the direction of all heat transfer to be towards the node, the finite difference equations for the nodes are obtained to be as follows: Node 1: []0)(45sin 2/2)(2100121=-++-+-T T h q l l T T l k T T l h s s i iα Node 2: 022)(2324211=-+-+-+-lT T kl l T T l k l T T l k T T l h i iNode 3: []0)(45sin 3003532=-++-+-T T h q ll T T kl l T T kls s α Node 4: 022)(245424=-+-+-lT T l k l T T l k T T l h i iNode 5: 0425634=-++T T T T Node 6: []0)(45sin 2/260065=-++-T T h q l l T T l ks s α where l = 1 m, k = 0.6 W/m ⋅︒C, h i =150 W/m 2⋅︒C, T i =15︒C, h o = 30 W/m 2⋅︒C, T 0 =25︒C, αs = 0.7, and2 W/m 800=s q. The system of 6 equations with 6 unknowns constitutes the finite difference formulation of the problem. The 6 nodal temperatures under steady conditions are determined by solving the 6 equations above simultaneously with an equation solver to beT 1 = T top =21.3︒C, T 2 =15.1︒C, T 3 = T middle =43.2︒C, T 4 =15.1︒C, T 5 =36.3︒C, T 6 = T bottom =43.6︒C Discussion Note that the highest temperature occurs at a location furthest away from the water, as expected.InsulatedWater h i , T i1。

相关文档
最新文档