CN111931370B - COMSOL-based ceramic insulator flash firing method - Google Patents

COMSOL-based ceramic insulator flash firing method Download PDF

Info

Publication number
CN111931370B
CN111931370B CN202010780521.9A CN202010780521A CN111931370B CN 111931370 B CN111931370 B CN 111931370B CN 202010780521 A CN202010780521 A CN 202010780521A CN 111931370 B CN111931370 B CN 111931370B
Authority
CN
China
Prior art keywords
sample
flash
temperature
simulation
current
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.)
Active
Application number
CN202010780521.9A
Other languages
Chinese (zh)
Other versions
CN111931370A (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.)
Tianjin University
Original Assignee
Tianjin 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 Tianjin University filed Critical Tianjin University
Priority to CN202010780521.9A priority Critical patent/CN111931370B/en
Publication of CN111931370A publication Critical patent/CN111931370A/en
Application granted granted Critical
Publication of CN111931370B publication Critical patent/CN111931370B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like

Abstract

The method utilizes simulation software COMSOL, takes alumina as a sample, establishes a flash combustion model, obtains the temperature distribution of the surface of the sample, and analyzes the influence of the electric field intensity and the current density on the density and the size of the sample. The method is beneficial to deepening understanding of flash firing, and has important significance for producing high-performance ceramic materials and high-performance ceramic insulators and improving the safety and stability of high-voltage transmission lines in China.

Description

COMSOL-based ceramic insulator flash firing method
Technical Field
The invention relates to the technical field of ceramic insulator preparation, in particular to a research on a flash firing process for preparing a ceramic insulator by utilizing COMSOL simulation software.
Background
The ceramic insulator has high hardness, high melting point, good chemical stability at high temperature, good electrical insulation, corrosion resistance and the like, and is widely applied to a power grid. However, the ceramic insulator has poor toughness and high brittleness, and is easy to break, and the ceramic insulator is easy to have surface flashover in the operation process according to the field operation. The properties of a ceramic insulator are closely related to its manufacturing process. The traditional sintering process not only needs high sintering temperature and long sintering time, consumes a large amount of energy, but also generates a large amount of greenhouse gas and is very unfavorable for the environment. And the density of the prepared sample is not high enough, which can affect the performance of the ceramic material. In order to solve the above problems, many new technologies have been proposed in recent years. Compared with other novel technologies, the flash sintering has the characteristics of low sintering temperature, high speed, short heat preservation time and the like, so that the flash sintering process is widely concerned.
At present, scholars at home and abroad have deeply studied flash burning process parameters and flash burning mechanisms and obtain a large number of stage achievements. However, students set up flash burning experimental models to analyze the flash burning process. Due to different built experiment platforms, the obtained experiment results do not reach the same. In addition, because the flash burning process is short in time, some parameters of the sample cannot be accurately measured, so that the research on the flash burning mechanism is controversial. Therefore, through establishing the simulation model of the flash combustion, people can conveniently analyze the influence of the experimental parameters, the measurement accuracy is improved, the research of the flash combustion mechanism is hopefully promoted, the actual experimental condition is predicted at first, the cost spent in the actual experiment is reduced, the judgment on the actual experimental result is facilitated, the error or the mistake of the actual experiment is found more quickly, and some dangerous actual operations can be avoided.
Disclosure of Invention
The invention aims to establish a flash combustion model taking alumina as a sample in COMSOL, introduce the steps of establishing a simulation model in detail, obtain the surface temperature of the sample in a stable stage and discuss the influence of electric field intensity and current density on the sample. The research result can be used for researching a flash combustion mechanism, analyzing the influence of factors such as electric field intensity, current density and the like on a flash combustion sample in a flash combustion experiment, opening up a new path for producing high-performance ceramic materials and high-performance ceramic insulators, and having important significance for improving the safety and stability of high-voltage transmission lines in China.
In order to achieve the purpose, the invention adopts the following technical scheme that a ceramic insulator flash firing method based on COMSOL comprises the following specific steps:
1. selection of physical fields
COMSOLMIC is a powerful multi-physical simulation software, and can simulate physical processes in multiple fields. The flash process involves electrical, thermal, etc. physical fields, and thus "joule heat and thermal expansion" physical fields are chosen in COMSOL. The "joule heating and thermal expansion" multi-physical interface involves a combination of thermal, electrical and structural multi-physical effects. The current distribution, the electric field distribution and the potential distribution in the dielectric medium can be calculated in a "current" interface. The "solid mechanics" interface can calculate strain, stress and displacement of the sample. The "solid phase heat transfer" interface can simulate conduction heat transfer, radiation heat transfer and current heat transfer in a model, and the temperature in the interface can cause thermal expansion as the heat load of the "solid mechanics" interface.
2. Establishment of geometric model
The flash-burned sample is a simple three-dimensional geometry, thus creating a cuboid to simulate the sample, which helps us to better understand the flash-burning process, as shown in fig. 1. Wherein the length, width and height of the cuboid are respectively 2.5 μm, 2.5 μm and 1 μm. The two planes along the X-axis direction were fixed, and the material of the sample was set to alumina and the surrounding atmosphere was set to air.
3. Application of boundary conditions
This step is the most critical step in the COMSOL software, and whether the boundary condition setting is reasonable directly determines whether the simulation can be successful. In the flash combustion experiment process, the variables applied externally comprise two variables of a direct current voltage source and a high temperature furnace, so the two conditions also need to be added in the simulation process.
In practical process, the conductivity of alumina will increase with the increase of temperature, so we need to find the conduction current option in the setting window of current conservation under the 'current' interface, and input the relation between the conductivity of alumina and temperature. Thus the current through the sample during the simulation increases with increasing temperature. Two planes along the X-axis are set, one as the input for voltage and the other to ground. In the simulation process, the voltage value can be changed, namely the voltage intensity is changed, and further the influence of the voltage intensity on the sample can be analyzed. The input voltage was set to 0.0375V and 0.075V, and the flash process was studied in both cases.
The initial temperature of the sample was set to 800K at the solid heat transfer node of the model developer window. The purpose of this is to provide an initial temperature of the flash, at which the entire geometric model is maintained. In order to simulate the cooling of the sample with convection air in the flash-burn test, a heat flux boundary condition with a heat transfer coefficient h of 5W/(m2 · K) was set.
4. Partitioning of a grid
The grid of the geometric model directly determines the solving mode of the model, and different grid dividing modes can cause different model solving time, different solving precision and different memory amount required by calculating the simulation model. The finer the grid division is, the higher the solving precision is, and the obtained simulation result is closer to an actual value. Therefore, for different geometric models and different problems, a proper mesh division mode needs to be selected. The simulation accuracy and the simulation time are comprehensively considered, the self-existing tetrahedral mesh is selected for division, and the mesh size setting column is selected to be more detailed. The mesh partition of the geometric model is shown in fig. 2.
The method uses COMSOL to simulate the flash firing process for preparing the ceramic insulator, has important significance for guaranteeing the stability and the reliability of the power transmission line, and has the following specific beneficial effects:
1. although flash firing has been applied and a range of ceramic materials produced, the understanding and study of the flash firing mechanism is limited by the uncertainty of the temperature measurement, and the exploration of the flash firing mechanism is greatly limited by the sharp change and difficulty of measurement of the temperature. The simulation obtains the surface temperature of the flash burning sample in the stable stage, is beneficial to the exploration of the flash burning mechanism, and provides support for the flash burning mechanism proposed by the scholars.
2. The influence of voltage intensity and current density on the density of the sample and the size of the sample is analyzed by utilizing the built flash burning model of the ceramic insulator, the relation between the voltage intensity and the current density is researched, the density of the sample can be really improved by flash burning, and the flash burning model has guiding significance for selecting the voltage intensity, the current density and the temperature in a flash burning experiment.
Drawings
Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention. Also, like reference numerals are used to refer to like parts throughout the drawings. In the drawings:
FIG. 1 is a geometric model diagram of a flash burn simulation;
FIG. 2 is a geometric model mesh partition diagram of a flash simulation;
FIG. 3 is a temperature profile of the surface of a flash-burned sample;
FIG. 4 is a graph of current density over time at different electric field strengths;
FIG. 5 is a graph of surface displacements of flash-burned samples at different voltage levels;
FIG. 6 is a graph showing the surface displacement of the flash sample at a current density of 1.74A/m 2.
Detailed Description
Exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
1. Flash-burned sample surface temperature analysis:
the surface temperature profile of the flash-burned sample at the stabilization stage when the electric field strength was 300V/cm by changing the DC voltage value is shown in FIG. 3. It can be seen that there is a temperature gradient across the surface of the sample and that the temperature at many locations on the surface of the sample is greater than 800K. Due to the fact that the surface temperature of the sample is increased, the internal void volume is reduced, the particle distance is shortened, the density of the sample is increased, and a stable crystal structure is formed more quickly. This is a good explanation for the short flash time and the rapid densification.
The increase in the temperature of the sample surface can be explained by the joule heating effect. The sample generates joule heat under the action of the current, and the joule heat generates a large amount of heat to raise the ambient temperature. Since the conductivity of alumina is directly proportional to temperature, an increase in temperature increases the conductivity of the sample. When the electric field strength is constant, the electrical conductivity increases, i.e., the resistance decreases, causing an increase in the thermal effect of the electrical energy and an increase in the temperature. This process is a positive feedback and causes the temperature of the sample to rise continuously. The temperature of the sample will increase endlessly if the current density is not limited.
2. Influence of electric field intensity:
when the applied DC voltage was 0.075V, i.e., the electric field intensity was 300V/cm, the change of the current density with time is shown in FIG. 4 (a). According to the simulation 2 result diagram, the flash burning can be divided into three stages: 0-0.25s, 0.25-0.5s, and later 0.5 s. In the first phase, the electric field strength remains constant and the current density and power increase slowly, this phase, called the inoculation phase, is a phase before flash-over occurs. In the second phase, the temperature of the sample is constantly raised by the effect of the joule heating effect, reaching the threshold of flash, so that the sample undergoes flash, rapidly shrinking and accompanied by electroluminescence. Therefore, the second stage is called the flash-burning generation stage, and the most obvious phenomenon in this stage is that the current density rapidly increases. When the current density reaches the critical value, the current density is kept unchanged, and the third stage is entered. The third stage is called the holding stage of flash firing, the current density is not changed any more, and the resistivity of the sample is kept unchanged until the sample is cooled to room temperature.
When the electric field strength was 150V/cm, the current density through the sample was plotted against time as shown in FIG. 4 (b). By comparing fig. 4(a) and fig. 4(b), we can see that the electric field strength has an influence on the incubation time, and the higher the electric field strength, the shorter the incubation time. When the value of the input dc voltage is changed, we find that the displacement maps of the sample surface are the same, as shown in fig. 5. From this result, it can be derived: the electric field intensity can not influence the compactness of the sample.
Thus, the electric field strength does not affect the variation in density of the sample, it only causes the variation in time of the flash inoculation stage.
3. Influence of current density:
the current density is an important factor because it is directly related to the temperature of the sample. The greater the current density, the higher the temperature of the sample. The current density was set to 1.74A/m 2 The displacement of the sample surface is shown in fig. 6. By comparing fig. 5 and fig. 6, we have found that differences in current density lead to differences in sample surface displacement. In other words, the greater the current density, the greater the density of the sample, but the inevitable increase in the size of the sample.
When the chemical components constituting the ceramic insulator are the same, the microstructure of the ceramic is a major factor affecting the performance of the ceramic. In order to improve the insulating property of the ceramic insulator, the performance of the insulator can be well improved by preparing the insulator through a flash firing process. In the flash combustion process, the inoculation time can be adjusted by controlling the voltage intensity, and the density and the size of the sample can be adjusted by controlling the current density.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the appended claims.

Claims (2)

1. A ceramic insulator flash-firing method based on COMSOL is characterized by comprising the following specific steps:
1) selection of physical fields
COMSOLMIC is a powerful multi-physical simulation software, and can simulate physical processes in multiple fields; the flash burning process involves physical fields such as electricity, heat and the like, so that the physical fields of 'joule heat and thermal expansion' are selected in COMSOL; the "joule heat and thermal expansion" multi-physical interface comprises a combination of thermal, electrical and structural multi-physical effects; wherein the current distribution, electric field distribution and electric potential distribution in the dielectric medium can be calculated in a 'current' interface; the "solid mechanics" interface can calculate strain, stress and displacement of the sample; the solid phase heat transfer interface can simulate conduction heat transfer, radiation heat transfer and current heat transfer in a model, and the temperature in the interface is used as the heat load of the solid mechanical interface to cause thermal expansion;
2) establishment of geometric model
The flash-burned sample is a simple three-dimensional geometry, so a cuboid is established to simulate the sample, which helps us to better understand the flash-burning process; wherein the length, the width and the height of the cuboid are respectively 2.5 micrometers, 2.5 micrometers and 1 micrometer; fixing two planes along the X-axis direction, setting the material of the sample to be alumina and setting the surrounding environment to be air;
3) application of boundary conditions
The step is the most critical step in COMSOL software, and whether the boundary condition setting is reasonable or not directly determines whether the simulation can be successful or not; in the flash combustion experiment process, the externally applied variables comprise a direct current voltage source and a high temperature furnace, so that the two conditions also need to be added in the simulation process; in the practical process, the conductivity of the alumina is increased along with the increase of the temperature, so that a conduction current option needs to be found in a setting window of current conservation under a current interface, and the relation between the conductivity of the alumina and the temperature is input; thus the current through the sample during the simulation increases with increasing temperature; setting two planes along the X-axis direction as an input end of voltage, and grounding the other plane; in the simulation process, the voltage value can be changed, namely the voltage intensity is changed, so that the influence of the voltage intensity on the sample can be analyzed; the input voltage was set to 0.0375V and 0.075V, the flash process was studied in both cases; setting the initial temperature of the sample to 800K at the solid heat transfer node of the model developer window; the purpose of this is to provide an initial temperature of the flash, at which the entire geometric model is maintained; in order to simulate the cooling of the sample and the convection air in the flash combustion experiment, a heat flux boundary condition with the heat transfer coefficient h being 5W/(m 2K) is set;
4) partitioning of a grid
The grid of the geometric model directly determines the solving mode of the model, and different grid dividing modes can cause different model solving time, different solving precision and different memory amount required by calculating the simulation model; the finer the grid division is, the higher the solving precision is, and the obtained simulation result is closer to an actual value; therefore, for different geometric models and different problems, a proper mesh division mode needs to be selected; the simulation accuracy and the simulation time are comprehensively considered, the self-existing tetrahedral mesh is selected for division, and the mesh size setting column is selected to be finer.
2. The COMSOL-based ceramic insulator flash firing method according to claim 1, characterized in that:
1) changing the direct current voltage value to obtain a surface temperature distribution diagram of the flash-burned sample at a stable stage; it can be found that the surface of the sample has a temperature gradient, and the temperature of a plurality of parts on the surface of the sample is more than 800K;
2) the electric field intensity does not influence the change of the density of the sample, and only causes the change of the time of a flash inoculation stage;
3) the greater the current density, the greater the density of the sample, but inevitably the larger the sample size.
CN202010780521.9A 2020-08-06 2020-08-06 COMSOL-based ceramic insulator flash firing method Active CN111931370B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010780521.9A CN111931370B (en) 2020-08-06 2020-08-06 COMSOL-based ceramic insulator flash firing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010780521.9A CN111931370B (en) 2020-08-06 2020-08-06 COMSOL-based ceramic insulator flash firing method

Publications (2)

Publication Number Publication Date
CN111931370A CN111931370A (en) 2020-11-13
CN111931370B true CN111931370B (en) 2022-08-02

Family

ID=73306809

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010780521.9A Active CN111931370B (en) 2020-08-06 2020-08-06 COMSOL-based ceramic insulator flash firing method

Country Status (1)

Country Link
CN (1) CN111931370B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109408901A (en) * 2018-09-29 2019-03-01 国网山西省电力公司太原供电公司 The modeling method of the three-dimensional simulation model of cable connector field distribution
CN109635442A (en) * 2018-12-13 2019-04-16 武汉理工大学 A kind of face inner mold thermo-electric device construction design method based on COMSOL software
CN110357633A (en) * 2019-07-11 2019-10-22 上海交通大学 A kind of method that room temperature quickly prepares titanium aluminium-carbon ceramic
CN110427637A (en) * 2019-06-06 2019-11-08 中国石油大学(华东) A kind of emulation mode for the direct current cables distribution of space charge considering temperature and electric-force gradient influence
CN110442907A (en) * 2019-07-02 2019-11-12 浙江中科电声研发中心 The Numerical Simulation Analysis method of piezoelectric type MEMS loudspeaker fundamental characteristics
CN110577399A (en) * 2019-07-12 2019-12-17 北京科技大学 Multi-field coupling flash sintering system based on induction heating
CN110750892A (en) * 2019-10-14 2020-02-04 河海大学 Multi-physical-field simulation method for lignite microwave drying
CN111079328A (en) * 2019-12-03 2020-04-28 上海理工大学 Motion simulation method for slow wire-moving electrospark wire-electrode cutting machining electroerosion product
CN111460713A (en) * 2020-03-31 2020-07-28 东北大学 Finite element analysis method based on temperature distribution of cladding material under electromagnetic induction heating condition

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109408901A (en) * 2018-09-29 2019-03-01 国网山西省电力公司太原供电公司 The modeling method of the three-dimensional simulation model of cable connector field distribution
CN109635442A (en) * 2018-12-13 2019-04-16 武汉理工大学 A kind of face inner mold thermo-electric device construction design method based on COMSOL software
CN110427637A (en) * 2019-06-06 2019-11-08 中国石油大学(华东) A kind of emulation mode for the direct current cables distribution of space charge considering temperature and electric-force gradient influence
CN110442907A (en) * 2019-07-02 2019-11-12 浙江中科电声研发中心 The Numerical Simulation Analysis method of piezoelectric type MEMS loudspeaker fundamental characteristics
CN110357633A (en) * 2019-07-11 2019-10-22 上海交通大学 A kind of method that room temperature quickly prepares titanium aluminium-carbon ceramic
CN110577399A (en) * 2019-07-12 2019-12-17 北京科技大学 Multi-field coupling flash sintering system based on induction heating
CN110750892A (en) * 2019-10-14 2020-02-04 河海大学 Multi-physical-field simulation method for lignite microwave drying
CN111079328A (en) * 2019-12-03 2020-04-28 上海理工大学 Motion simulation method for slow wire-moving electrospark wire-electrode cutting machining electroerosion product
CN111460713A (en) * 2020-03-31 2020-07-28 东北大学 Finite element analysis method based on temperature distribution of cladding material under electromagnetic induction heating condition

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"8YSZ陶瓷闪烧过程的温度场研究";刘光鑫;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20180915;B015-127 *
"Flash sintering of sodium niobate ceramics";Xinghua Su等;《Elsevier》;20190115;pp15-18 *

Also Published As

Publication number Publication date
CN111931370A (en) 2020-11-13

Similar Documents

Publication Publication Date Title
Oreshkin et al. Numerical calculation of the current specific action integral at the electrical explosion of wires
Lin et al. Thermal analysis on electromagnetic launcher under transient conditions
CN111931370B (en) COMSOL-based ceramic insulator flash firing method
Qiao et al. Heat transfer analysis of different conditions for SF $ _ {6} $/N $ _2 $ gas-insulated transmission lines
Liang et al. Pressure‐assisted flash sintering of ZnO ceramics
CN110045247B (en) Experimental method for accelerating thermal aging of transformer bushing
Ma et al. Heat dissipation performance of coaxial cable for electromagnetic launch
CN110083905B (en) Hot spot temperature analysis method for dynamic capacity increase of isolating switch
Zhao et al. Calculation on heating effect due to void discharge in polymers in cumulative breakdown process
Lee et al. Hydrodynamic Modeling for Discharge Analysis in a Dielectric Medium with the Finite Element Method under Lightning Impulse
Sun et al. Simulation of the temperature field distribution in medium-voltage vacuum interrupter and experimental verification
Bespalov et al. Determination of temperature parameters of IGBT, used in converter of AC electric drive
Zhang et al. Investigation on space charge behavior in water tree aged crosslinked polyethylene (XLPE) cable by experiment and simulation
Si et al. Experimental analysis on influence of temperature and bolted torque to DC contact resistance of terminal connector in UHV valve hall
Qi et al. Simulation of Electrical Tree Breakdown Process in Epoxy Resin with Phase Field Modeling
Muñoz-Hernández et al. Modeling of thermal runaway of carbonaceous materials: Graphite, biochar, and wood
Wang et al. A resistance model of fault current limiting coils under dc impact
CN213482354U (en) High-temperature ferroelectric performance test fixture
Huang et al. Effects of Flash Sintering Parameters on Insulation Performance of Ceramic Insulator
Zeng et al. Heat transfer characteristics of graphite welding pool in penetration brazing of Cu/Ni dissimilar metal wires
Haque Modeling supercritical fluids and fabricating electret films to address dielectric challenges in high-power-density systems
Kadim et al. Nanoparticle Effects on the Electrical and Thermal Performance of the Transformer Mineral Oil
Guan et al. Analysis of Abnormal Heating of Outdoor High Voltage Isolating Switch
Tong et al. Electromagnetic-thermal coupled analysis of a 10 kV silicone rubber cast dry-type transformer
Feng et al. Research on Temperature Field of All-film Pulsed Capacitors in Application of Repetitive Pulse Discharge

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