CN109297304B - A kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis - Google Patents

A kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis Download PDF

Info

Publication number
CN109297304B
CN109297304B CN201811093706.1A CN201811093706A CN109297304B CN 109297304 B CN109297304 B CN 109297304B CN 201811093706 A CN201811093706 A CN 201811093706A CN 109297304 B CN109297304 B CN 109297304B
Authority
CN
China
Prior art keywords
heating
alkali metal
temperature
gas chamber
finite element
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201811093706.1A
Other languages
Chinese (zh)
Other versions
CN109297304A (en
Inventor
丁铭
张吉
陆吉玺
姚涵
蔡佳书
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201811093706.1A priority Critical patent/CN109297304B/en
Publication of CN109297304A publication Critical patent/CN109297304A/en
Application granted granted Critical
Publication of CN109297304B publication Critical patent/CN109297304B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F27FURNACES; KILNS; OVENS; RETORTS
    • F27BFURNACES, KILNS, OVENS, OR RETORTS IN GENERAL; OPEN SINTERING OR LIKE APPARATUS
    • F27B17/00Furnaces of a kind not covered by any preceding group
    • F27B17/0016Chamber type furnaces
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F27FURNACES; KILNS; OVENS; RETORTS
    • F27DDETAILS OR ACCESSORIES OF FURNACES, KILNS, OVENS, OR RETORTS, IN SO FAR AS THEY ARE OF KINDS OCCURRING IN MORE THAN ONE KIND OF FURNACE
    • F27D11/00Arrangement of elements for electric heating in or on furnaces
    • F27D11/02Ohmic resistance heating
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F27FURNACES; KILNS; OVENS; RETORTS
    • F27DDETAILS OR ACCESSORIES OF FURNACES, KILNS, OVENS, OR RETORTS, IN SO FAR AS THEY ARE OF KINDS OCCURRING IN MORE THAN ONE KIND OF FURNACE
    • F27D19/00Arrangements of controlling devices
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F27FURNACES; KILNS; OVENS; RETORTS
    • F27DDETAILS OR ACCESSORIES OF FURNACES, KILNS, OVENS, OR RETORTS, IN SO FAR AS THEY ARE OF KINDS OCCURRING IN MORE THAN ONE KIND OF FURNACE
    • F27D19/00Arrangements of controlling devices
    • F27D2019/0003Monitoring the temperature or a characteristic of the charge and using it as a controlling value
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F27FURNACES; KILNS; OVENS; RETORTS
    • F27DDETAILS OR ACCESSORIES OF FURNACES, KILNS, OVENS, OR RETORTS, IN SO FAR AS THEY ARE OF KINDS OCCURRING IN MORE THAN ONE KIND OF FURNACE
    • F27D19/00Arrangements of controlling devices
    • F27D2019/0028Regulation
    • F27D2019/0034Regulation through control of a heating quantity such as fuel, oxidant or intensity of current
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F27FURNACES; KILNS; OVENS; RETORTS
    • F27DDETAILS OR ACCESSORIES OF FURNACES, KILNS, OVENS, OR RETORTS, IN SO FAR AS THEY ARE OF KINDS OCCURRING IN MORE THAN ONE KIND OF FURNACE
    • F27D19/00Arrangements of controlling devices
    • F27D2019/0096Arrangements of controlling devices involving simulation means, e.g. of the treating or charging step

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

The present invention relates to a kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis, according to the asymmetrical feature of alkali metal gas chamber Roasting oven structure, the optimum heating temperature for being pasted on oven outer surface heating film is calculated using the method for finite element analysis, make alkali metal plenum interior thermo parameters method more evenly using the heating condition, to obtain atomicity Density Distribution more evenly.This method can greatly improve the uniformity of alkali metal plenum interior thermo parameters method in the case of atmospheric chamber multiposition measuring, to improve the uniformity of detection signal.Compared with thermal current heating means, this method can greatly improve heating efficiency;Compared with traditional Electric heating, the present invention can get the heating effect of thermo parameters method more evenly.The heating means can for atmospheric chamber, multiposition measuring atom magnetometer biomagnetic measurement field application reference is provided.

Description

A kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis
Technical field
The present invention relates to a kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis, can be used for atom magnetic strength The heating of alkali metal gas chamber is counted, to obtain alkali metal gas chamber thermo parameters method more evenly.
Background technique
Alkali metal gas chamber is the important Sensitive Apparatus of atom magnetometer.Due to the atom number density point of alkali metal plenum interior Cloth is directly related with temperature, and in the case of atmospheric chamber, the measurement of multidigit point, the uniformity of atomicity Density Distribution be will have a direct impact on The quality of measuring signal is influenced to the uniformity of multidigit point magnetic field sensing signal.Thus it is guaranteed that under the premise of heating efficiency, The uniformity of alkali metal gas chamber thermo parameters method is of great significance in raising atmospheric chamber, the measurement of multidigit point.
Alkali metal gas chamber is heated by the way of air-flow heating earliest, can effectively avoid using the heating method Heating source and alkali metal gas chamber directly contact, it can be achieved that preferably without magnetic heating effect, and can get uniform alkali metal gas Room temperature field distribution, disadvantage are that heating efficiency is low, and it also requires the external oven of a large volume generates thermal current not Conducive to the miniaturization of atom magnetometer;Current atom magnetometer generally using the electrically heated mode of no magnetic come to alkali metal gas chamber into Row heating the advantage is that heating efficiency height, its shortcoming is that the heating method of currently used multi-panel isothermal heating can be led The uneven of alkali metal plenum interior heating temperature field distribution is caused, will have a direct impact on atmospheric chamber, multidigit point magnetic-field measurement letter Number uniformity to influence the quality of transducing signal.
Summary of the invention
The technology of the present invention solves the problems, such as: in order to overcome the above-mentioned deficiencies of the prior art, providing a kind of based on finite element analysis Alkali metal gas chamber non-uniform heating methods, alkali metal plenum interior temperature in the case of atmospheric chamber multiposition measuring can be greatly improved The uniformity of field distribution, to improve the uniformity of detection signal.
The problem to be solved in the present invention is: brought under multi-disc heating film isothermal heating condition for the heating of no magnetoelectricity Heat the non-uniform problem of plenum interior thermo parameters method, provide it is a kind of under the premise of guaranteeing alkali metal gas chamber high heating efficiency simultaneously Obtain the heating means of the heating temperature field distribution inside it more evenly.
The technical scheme adopted by the invention is that: a kind of alkali metal gas chamber non-uniform heating side based on finite element analysis Method is calculated respectively by the method for introducing finite element analysis according to the asymmetric feature of alkali metal gas chamber heating structure itself The optimal heating temperature of heating film, and its optimum results is applied in practical heating process, to realize to alkali metal gas chamber Optimization heating.The specific implementation steps are as follows:
(1) physical model of alkali metal gas chamber and its heating structure is constructed according to the heating structure of alkali metal gas chamber first; Then constructed physical model is imported into finite element analysis software, according to alkali metal gas chamber and its heating structure material The thermodynamic parameter of finite element analysis is arranged in performance parameter, to complete as composed by alkali metal gas chamber and its heating structure The building of the finite element analysis model of heating system;
(2) based on finite element analysis model constructed by step (1), first fixed base metal air chamber outside heating structure is symmetrical The heating temperature T of part heating surface1, change the heating temperature T in the asymmetric face of heating structure outside alkali metal gas chamber2, carry out limited Meta analysis, and calculate in different T2At a temperature of alkali metal plenum interior Steady-State Thermal Field obtained coefficient of variation CV value; With (T2-T1) it is abscissa, with (T2-T1) corresponding to CV value be ordinate carry out described point and carry out curve fitting, according to fitting tie Fruit solves CV and obtains (T corresponding when minimum2-T1) temperature value, that is, obtain in T1At a temperature of optimum heating temperature T2Value;
(3) T obtained according to step (2)1At a temperature of optimum heating temperature T2It is accurately controlled using feedback control principle viscous It is attached to the heating temperature of heating film on the alkali metal gas chamber heating structure plane of symmetry and asymmetric face, to obtain in alkali metal gas chamber The most uniform heating effect of portion's thermo parameters method.
The step (1) physical model passes through SolidWorks, ProE, Autodesk Inventor software building; Heating system finite element analysis model is constructed by the completion of ANSYS, COMSOL, NASTRAN software.
In the step (2), the calculation formula of alkali metal plenum interior temperature field coefficient of variation CV is as follows:
Wherein f is the standard deviation of plenum interior thermo parameters method, TaverFor the mean temperature of plenum interior.
The calculation formula of the f is as follows:
Wherein T indicates to be located at the steady temperature value at coordinate points (x, y, z), and V indicates the volume of alkali metal gas chamber.
The mean temperature T of plenum interioraverCalculation formula it is as follows:
Taver=∫ ∫ ∫ Tdxdydz/V
Temperature used by the step (3) controls closed loop feedback system by temperature transducer, PID control circuit and power amplifier Composition;The heating temperature that heating film is monitored by temperature transducer, is controlled according to temperature monitoring result by PID control electronic feedback The output power of power amplifier is to realize the accurate control to heating film heating temperature.
The advantages of present invention is compared with existing heating technique is:
(1) present invention is calculated according to the asymmetrical feature of alkali metal gas chamber Roasting oven structure using the method for finite element analysis It is pasted on the optimum heating temperature of oven outer surface heating film out, makes alkali metal plenum interior temperature field point using the heating condition Cloth more evenly, to obtain atomicity Density Distribution more evenly.
(2) present invention can greatly improve the equal of alkali metal plenum interior thermo parameters method in the case of atmospheric chamber multiposition measuring Even property, to improve the uniformity of detection signal.Compared with thermal current heating means, this method can greatly improve heating efficiency; Compared with traditional Electric heating, this method can get the heating effect of thermo parameters method more evenly.The heating means can be big Gas chamber, multiposition measuring atom magnetometer biomagnetic measurement field application provide reference.
(3) under the premise of not influencing no magnetoelectricity heating heating efficiency, using the method for finite element analysis, optimization heating temperature Degree condition is, it can be achieved that alkali metal gas chamber heating temperature field uniformity is substantially improved.
(4) promotion of Uniformity of Temperature Field, so that the uniformity of atom magnetometer atmospheric chamber, more site primer signals obtains It greatly improves, is of great significance to its application in biomagnetic measurement field.
Detailed description of the invention
Fig. 1 is the flow chart that heating means of the present invention are realized;
Fig. 2 is the heating structure schematic diagram that heating means of the present invention are realized.
Specific embodiment
The following describes the present invention in detail with reference to the accompanying drawings and embodiments.
As shown in Figure 1, specific implementation step of the invention is as follows:
(1) Fig. 2 is a kind of typical alkali metal gas chamber heating structure, and similar four side structures are symmetrical structure, bottom Face and surface structure it is different be unsymmetric structure, first according to alkali metal gas chamber 4 and its outside heating ceramic 3 practical ruler The physical model of very little its analysis of building;Constructed physical model is imported into finite element analysis software, according to alkali metal gas chamber 4 Institute's inflatable body type and the material of heating ceramic 3, determine its respective thermal coefficient and convection coefficient, foundation in material, gas chamber The thermodynamic parameter of these parameter setting finite element analyses, to complete to add as composed by alkali metal gas chamber and its heating structure The building of the finite element analysis model of hot systems.
(2) according to finite element analysis model constructed in (1), the heating temperature T of four side heat films 2 is first fixed1 (such as: 180 DEG C) changes the heating temperature T of bottom surface and top heating film 12, to the certain heating temperature difference (T2-T1) in range (such as :- 40~40 DEG C) multiple temperature spots carry out finite element analysis, to gas indoor temperature field obtained under each heating temperature condition Steady-state distribution result seek coefficient of variation CV, with (T2-T1) it is abscissa, it is vertical sit with CV value corresponding under the heating condition Mark described point simultaneously carries out curve fitting, and solves corresponding (T when CV obtains minimum according to fitting result2-T1) temperature value, It obtains in T1At a temperature of optimum heating temperature T2Value.The calculation formula of temperature field coefficient of variation CV is as follows:
Wherein f is the standard deviation of plenum interior thermo parameters method, TaverFor the mean temperature of plenum interior.
The calculation formula of f is as follows:
Wherein T indicates to be located at the steady temperature value at coordinate points (x, y, z), and V indicates the volume of alkali metal gas chamber.
The mean temperature T of plenum interioraverCalculation formula it is as follows:
Taver=∫ ∫ ∫ Tdxdydz/V
(3) T obtained according to step (2)1At a temperature of optimum heating temperature T2, accurately controlled using feedback control principle It is pasted onto the heating temperature of heating film on the alkali metal gas chamber heating structure plane of symmetry and asymmetric face, used temperature control is closed Ring feedback system is by temperature transducer, PID control circuit and power amplifier composition;The heating temperature of heating film is monitored by temperature transducer Degree heats temperature to heating film to realize by the output power of PID control electronic feedback control power amplifier according to temperature monitoring result The accurate control of degree, to obtain the most uniform heating effect of alkali metal plenum interior thermo parameters method.
The content that description in the present invention is not described in detail belongs to the prior art well known to professional and technical personnel in the field. Although the illustrative specific embodiment of the present invention is described above, in order to which the technical staff of this technology neck understands this Invention, it should be apparent that coming the present invention is not limited to the range of specific embodiment to those skilled in the art It says, as long as various change is in the spirit and scope of the present invention that the attached claims limit and determine, these variations are aobvious And be clear to, all are using the innovation and creation of present inventive concept in the column of protection.

Claims (6)

1. a kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis, it is characterised in that: the following steps are included:
(1) physical model of alkali metal gas chamber and its heating structure is constructed according to the heating structure of alkali metal gas chamber first;Then Constructed physical model is imported into finite element analysis software, according to the performance of alkali metal gas chamber and its heating structure material The thermodynamic parameter of finite element analysis is arranged in parameter, to complete to heat as composed by alkali metal gas chamber and its heating structure The building of the finite element analysis model of system;
(2) based on finite element analysis model constructed by step (1), first fixed base metal air chamber outside heating structure symmetric part The heating temperature T of heating surface1, change the heating temperature T in the asymmetric face of heating structure outside alkali metal gas chamber2, carry out finite element fraction Analysis, and calculate in different T2At a temperature of alkali metal plenum interior Steady-State Thermal Field obtained coefficient of variation CV value;With (T2- T1) it is abscissa, with (T2-T1) corresponding to CV value be ordinate carry out described point and carry out curve fitting, solved according to fitting result CV obtains (T corresponding when minimum out2-T1) temperature value, that is, obtain in T1At a temperature of optimum heating temperature T2's Value;
(3) T obtained according to step (2)1At a temperature of optimum heating temperature T2It is accurately controlled and is pasted onto using feedback control principle The heating temperature of heating film on the alkali metal gas chamber heating structure plane of symmetry and asymmetric face, to obtain alkali metal plenum interior temperature Spend the most uniform heating effect of field distribution.
2. a kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis according to claim 1, feature Be: the step (1) physical model passes through SolidWorks, ProE or Autodesk Inventor software building;Add Hot systems finite element analysis model is constructed by the completion of ANSYS, COMSOL or NASTRAN software.
3. a kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis according to claim 1, feature Be: in the step (2), the calculation formula of alkali metal plenum interior temperature field coefficient of variation CV is as follows:
Wherein f is the standard deviation of plenum interior thermo parameters method, TaverFor the mean temperature of plenum interior.
4. a kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis according to claim 3, feature Be: the calculation formula of the f is as follows:
Wherein T indicates to be located at the steady temperature value at coordinate points (x, y, z), and V indicates the volume of alkali metal gas chamber.
5. a kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis according to claim 3, feature It is: the mean temperature T of plenum interioraverCalculation formula it is as follows:
Taver=∫ ∫ ∫ Tdxdydz/V.
6. a kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis according to claim 1, feature Be: temperature used by the step (3) controls closed loop feedback system by temperature transducer, PID control circuit and power amplifier group At;The heating temperature that heating film is monitored by temperature transducer controls function by PID control electronic feedback according to temperature monitoring result The output power put is to realize the accurate control to heating film heating temperature.
CN201811093706.1A 2018-09-19 2018-09-19 A kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis Expired - Fee Related CN109297304B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811093706.1A CN109297304B (en) 2018-09-19 2018-09-19 A kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811093706.1A CN109297304B (en) 2018-09-19 2018-09-19 A kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis

Publications (2)

Publication Number Publication Date
CN109297304A CN109297304A (en) 2019-02-01
CN109297304B true CN109297304B (en) 2019-08-30

Family

ID=65163310

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811093706.1A Expired - Fee Related CN109297304B (en) 2018-09-19 2018-09-19 A kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis

Country Status (1)

Country Link
CN (1) CN109297304B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110859328B (en) * 2019-10-29 2024-03-26 贵州中烟工业有限责任公司 Design method of heat-conducting pipe of non-burning cigarette heater

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105004884A (en) * 2015-07-03 2015-10-28 北京航空航天大学 SiC-based micro-optics high-temperature accelerometer and design method
CN105403322A (en) * 2015-12-11 2016-03-16 东南大学 Apparatus and method for measuring temperature distribution in alkali metal air chamber of atom magnetometer
JP2016079470A (en) * 2014-10-17 2016-05-16 日立Geニュークリア・エナジー株式会社 Defect progress suppression method
CN106019180A (en) * 2016-07-18 2016-10-12 北京航空航天大学 Alkali metal atomic magnetometer air chamber electric heating magnetic field measurement method
CN106949985A (en) * 2017-05-15 2017-07-14 北京航空航天大学 A kind of precision measurement method of the alkali metal plenum interior temperature based on mixing optical pumping
CN108445428A (en) * 2018-04-11 2018-08-24 北京航空航天大学 A kind of SERF atom magnetometers electronic polarizability measurement method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016079470A (en) * 2014-10-17 2016-05-16 日立Geニュークリア・エナジー株式会社 Defect progress suppression method
CN105004884A (en) * 2015-07-03 2015-10-28 北京航空航天大学 SiC-based micro-optics high-temperature accelerometer and design method
CN105403322A (en) * 2015-12-11 2016-03-16 东南大学 Apparatus and method for measuring temperature distribution in alkali metal air chamber of atom magnetometer
CN106019180A (en) * 2016-07-18 2016-10-12 北京航空航天大学 Alkali metal atomic magnetometer air chamber electric heating magnetic field measurement method
CN106949985A (en) * 2017-05-15 2017-07-14 北京航空航天大学 A kind of precision measurement method of the alkali metal plenum interior temperature based on mixing optical pumping
CN108445428A (en) * 2018-04-11 2018-08-24 北京航空航天大学 A kind of SERF atom magnetometers electronic polarizability measurement method

Also Published As

Publication number Publication date
CN109297304A (en) 2019-02-01

Similar Documents

Publication Publication Date Title
CN109297304B (en) A kind of alkali metal gas chamber non-uniform heating methods based on finite element analysis
CN203117139U (en) Hot plate instrument for testing hot performance of fabric
CN104215660B (en) A kind of method and system simultaneously can surveying solid material heat conductivity and thermal diffusivity
CN108844991B (en) Device and method for measuring heat transfer coefficient of concrete plate
CN201697888U (en) Thermal inertia index detection device for single material space enclosing structures
CN106017739A (en) Temperature sensor calibration device
CN112270109B (en) Method for simulating heating performance of graphite rod in high-temperature carbonization furnace
CN108519406A (en) A kind of computational methods of the conductor thermal resistance and thermal-convection resistance of aerial earth wire axial direction
CN108489855A (en) A kind of sound-absorbing material flow resistance measuring instrument of temperature-controllable
CN208297402U (en) Expansion coefficients of metal wire experiment instrument
CN106248253A (en) The temperature sensor calibrator (-ter) unit of uniform temperature fields
CN113252724B (en) Method for detecting heat preservation performance of external wall
CN113077692B (en) Air transverse single tube forced convection heat transfer teaching experiment device and use method
CN109269682A (en) A kind of caliberating device and scaling method of heat flow transducer
CN108663187A (en) Based on hand automatic double mode formula power match reflux wind-tunnel heating system and application method
CN103983363A (en) Optimal refrigerating plate for low-temperature infrared target source
CN208076054U (en) A kind of temperature measurement box of test COB light source glue surface temperature
CN205039016U (en) Enameled wire production is with line temperature control constant temperature equipment
CN208171483U (en) A kind of U-shaped positive and negative temperature difference exciting bank
CN109425633A (en) A kind of heat shrinkage film shrinking percentage test method and test device
CN206038581U (en) Cavity glass thermal insulation performance detection device
CN209093393U (en) Novel humidifying device energy efficiency temperature experimental box
CN214794555U (en) Detachable wall body thermal resistance detection device
CN204086184U (en) A kind of system simultaneously can surveying solid material heat conductivity and thermal diffusivity
CN208653673U (en) Short branch electric thermo-couple temperature verifies furnace

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190830

Termination date: 20210919