CN113887069A - Method for solving surface temperature of anti-icing component through cross-dimension coupling - Google Patents

Method for solving surface temperature of anti-icing component through cross-dimension coupling Download PDF

Info

Publication number
CN113887069A
CN113887069A CN202111230670.9A CN202111230670A CN113887069A CN 113887069 A CN113887069 A CN 113887069A CN 202111230670 A CN202111230670 A CN 202111230670A CN 113887069 A CN113887069 A CN 113887069A
Authority
CN
China
Prior art keywords
dimensional model
temperature
surface temperature
low
value
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.)
Granted
Application number
CN202111230670.9A
Other languages
Chinese (zh)
Other versions
CN113887069B (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.)
AECC Shenyang Engine Research Institute
Original Assignee
AECC Shenyang Engine Research Institute
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 AECC Shenyang Engine Research Institute filed Critical AECC Shenyang Engine Research Institute
Priority to CN202111230670.9A priority Critical patent/CN113887069B/en
Publication of CN113887069A publication Critical patent/CN113887069A/en
Application granted granted Critical
Publication of CN113887069B publication Critical patent/CN113887069B/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
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

The application relates to the technical field of aircraft anti-icing components, in particular to a method for solving the surface temperature of an anti-icing component through cross-dimensional coupling, which comprises the steps of obtaining the temperature value of each node on the surface of a high-dimensional model initial temperature field of the anti-icing component; acquiring the surface temperature of the anti-icing part of the low-dimensional model; obtaining a correction value of the surface temperature of the low-dimensional model by adopting a low-dimensional model surface temperature correction formula; determining a heat exchange quantity Q based on the water phase change of the low-dimensional model surface; loading the surface heat exchange quantity Q based on the low-dimensional model on the high-dimensional model to obtain the temperature value of each node on the surface of the high-dimensional model in the iteration step; and comparing the temperature difference value of the latest round and the previous round of the surface node of the high-dimensional model, and if the difference value is smaller than the convergence judgment standard, taking the temperature value of the latest round of the surface as the final calculation result. By adopting a cross-dimension coupling mode, the surface temperature of the anti-icing component can be efficiently and accurately solved.

Description

Method for solving surface temperature of anti-icing component through cross-dimension coupling
Technical Field
The application belongs to the technical field of aircraft anti-icing, and particularly relates to a method for solving surface temperature of an anti-icing component through cross-dimension coupling.
Background
When the surface temperature of the anti-icing component is solved, in consideration of the limitations of analysis cost and technical means, a low-dimensional analysis method and a high-dimensional analysis method are combined to establish an analysis method for solving the surface temperature of the anti-icing component through cross-dimensional coupling. When the surface temperature of the anti-icing component with water phase change heat exchange is solved, iterative calculation is often needed, and low-dimensional system-level parameters are loaded on a high-dimensional temperature analysis model; and then, processing a calculation result on the high-dimensional temperature analysis model, and transmitting the result to a low-dimensional system-level analysis mathematical model to realize the coupling calculation between the two models. The geometric anisotropy of the high-dimensional model cannot be considered by the low-dimensional model, and the high-dimensional model needs dimension reduction and transformation to meet the requirement of the low-dimensional system-level analysis mathematical model. How to solve the problem of dimension transformation when solving the surface temperature of the anti-icing component is the technical problem to be solved by the patent.
In order to solve the problem of solving the surface temperature of the anti-icing component through cross-dimensional coupling, it is necessary to ensure that the parameters transmitted to the low-dimensional model by the high-dimensional model temperature field result do not have geometric dependence characteristics, and the parameters transmitted to the high-dimensional model by the low-dimensional model are only associated with the parameters transmitted from the high-dimensional model and do not have geometric dependence characteristics.
The commercial software and the professional module are completely based on the same dimensionality when solving the problem of solving the surface temperature of the anti-icing component, the process of processing high-dimensionality calculation analysis is complex, and the development difficulty is extremely high. Therefore, how to utilize cross-dimension coupling to reduce the development difficulty and effectively solve the surface temperature of the anti-icing component is a problem to be solved.
Disclosure of Invention
The application aims to provide a method for solving the surface temperature of an anti-icing component through cross-dimensional coupling, so as to solve the problem that the difficulty in calculating and analyzing the surface temperature of the anti-icing component through a high-dimensional meter is high in the prior art.
The technical scheme of the application is as follows: a method for solving the surface temperature of an anti-icing component through cross-dimension coupling comprises the steps of obtaining the initial temperature of a high-dimension model of the anti-icing componentTemperature value T of each node on surface of degree fieldW1(ii) a Obtaining surface temperature T of anti-icing component of low-dimensional modelW2(ii) a Obtaining a correction value of the surface temperature of the low-dimensional model by adopting a low-dimensional model surface temperature correction formula; determining a heat exchange quantity Q based on the water phase change of the low-dimensional model surface; loading the surface heat exchange quantity Q based on the low-dimensional model on the high-dimensional model to complete the calculation of the surface temperature of the high-dimensional model, and obtaining the temperature value T 'of each node of the surface of the high-dimensional model in the iteration step'W1(ii) a Comparing the temperature difference value of the latest round and the previous round of the surface node of the high-dimensional model, and if the difference value is smaller than the convergence judgment standard, taking the temperature value of the surface of the latest round as the final calculation result; otherwise, repeating iterative operation of the surface temperature value of the high-dimensional model until the temperature difference value of each node on the surface of the high-dimensional model after two iterations is within the set threshold range.
Preferably, the correction formula of the surface temperature of the low-dimensional model is T'W2=a*TW1+b*TW2Wherein a + b is 1.
Preferably, the phase-change heat exchange quantity Q based on the low-dimensional model surface temperature water is Q ═ f (T'W2)。
Preferably, the surface temperature of the low-dimensional model anti-icing part is acquired in an energy balance mode.
As a specific implementation manner, a system for solving the surface temperature of an anti-icing component through cross-dimensional coupling comprises a first surface temperature acquisition module, a second surface temperature acquisition module and a third surface temperature acquisition module, wherein the first surface temperature acquisition module is used for acquiring temperature values of nodes on the surface of a high-dimensional model; the second surface temperature acquisition module is used for acquiring a surface temperature value of the low-dimensional model; the surface temperature correction module is used for correcting the surface temperature value of the low-dimensional model; the heat exchange quantity determining module is used for obtaining the heat exchange quantity Q of the low-dimensional model; the surface temperature iteration module is used for loading the heat exchange quantity Q on the high-dimensional model, completing the temperature calculation of each node on the surface of the high-dimensional model and obtaining the temperature value of each node on the surface of the high-dimensional model in the iteration step; the convergence judging module is used for comparing the temperature difference value of the latest round and the previous round of the surface nodes of the high-dimensional model, and if the difference value is smaller than the convergence judging standard, the temperature value of the latest round of the surface is the final calculation result; otherwise, repeating iterative operation of the surface temperature value of the high-dimensional model until the temperature difference value of each node on the surface of the high-dimensional model after two iterations is within the set threshold range.
The application provides a method for solving surface temperature of anti-icing component through cross-dimension coupling, surface temperature is used as parameter transmission between different dimension models, the surface temperature of each node of the surface of a high dimension model and the surface temperature of a low dimension model are obtained firstly, heat exchange quantity Q is obtained through a surface temperature correction value of the low dimension model and is loaded to the high dimension model, the surface temperature of each node of the surface of the high dimension model after iteration is obtained, then the surface temperature of the surface of the low dimension model is corrected again through the surface temperature of each node of the surface of the high dimension model after the iteration, repeated iteration is carried out, training convergence is achieved, and the surface temperature of the anti-icing component can be obtained quickly and accurately.
Drawings
In order to more clearly illustrate the technical solutions provided by the present application, the following briefly introduces the accompanying drawings. It is to be expressly understood that the drawings described below are only illustrative of some embodiments of the invention.
FIG. 1 is a schematic view of the overall flow structure of the present application;
fig. 2 is a schematic diagram of the overall system structure of the present application.
1. A first surface temperature acquisition module; 2. a second surface temperature acquisition module; 3. a temperature correction module; 4. a heat exchange amount determination module; 5. a surface temperature iteration module; 6. and a convergence judging module.
Detailed Description
In order to make the implementation objects, technical solutions and advantages of the present application clearer, the technical solutions in the embodiments of the present application will be described in more detail below with reference to the drawings in the embodiments of the present application.
A method for solving a surface temperature of an anti-icing component across a dimensional coupling, as shown in fig. 1, comprises:
step S100, obtaining the temperature value T of each node on the surface of the initial temperature field of the high-dimensional model of the anti-icing componentW1
Step S200, obtaining the surface temperature T of the anti-icing part of the low-dimensional modelW2
Step S300, obtaining a correction value of the surface temperature of the low-dimensional model by adopting a low-dimensional model surface temperature correction formula;
step S400, determining heat exchange quantity Q based on water phase change of the surface of the low-dimensional model;
step S500, loading the surface heat exchange quantity Q based on the low-dimensional model on the high-dimensional model to complete the calculation of the surface temperature of the high-dimensional model, and obtaining the temperature value T 'of each node of the surface of the high-dimensional model in the iteration step'W1
Step S600, comparing the temperature difference value of the latest round and the previous round of the surface node of the high-dimensional model, and if the difference value is smaller than the convergence judgment standard, namely abs (T)W1-T'W1) If the temperature value is less than dt, the temperature value of the surface of the latest wheel is the final calculation result; otherwise, repeating iterative operation of the surface temperature value of the high-dimensional model until the temperature difference value of each node on the surface of the high-dimensional model after two iterations is within the set threshold range.
Because the corresponding relation exists between the temperature of each node of the high-dimensional model and the surface temperature of the anti-icing component of the low-dimensional model, and the surface temperature value of the anti-icing component can be obtained more conveniently and accurately, the surface temperature of the anti-icing component of the low-dimensional model is calculated firstly, and then the heat exchange quantity is loaded on the high-dimensional model through a correction formula, so that the temperature value of each node of the high-dimensional model can be effectively obtained, the temperature value of each node of the high-dimensional model is iterated through the correction formula, and the temperature value of each node of the high-dimensional model can be accurately obtained after the training is converged.
Through carrying out the coupling solution anti-icing part surface temperature of striding the dimension based on surface temperature, avoided the influence that different dimension geometric parameters brought ingeniously, the parameter transmission of low dimension and high dimension all characterizes through the surface temperature parameter, and the heat transfer volume that surface temperature decided does not rely on the influence of geometric parameters again for the heat flow item parameter of low dimension surface temperature characterization is unanimous with high dimension surface temperature characterization heat transfer volume, thereby realizes striding the dimension and solves anti-icing part surface temperature. The method is adopted for solving, namely, the efficiency of solving the surface temperature in a low-dimensional mode is achieved, the precision of solving the surface temperature in a high-dimensional mode is achieved, the adaptability is good, convenience and flexibility are achieved, and consumed resources and time are less.
Preferably, the surface temperature correction formula of the low-dimensional model is as follows: t'W2=a*TW1+b*TW2Wherein a + b is 1. When iteration is difficult to converge, the assignment of a can be increased, and the assignment of b can be reduced, so that the solution is converged; when iteration can be converged quickly, assignment of b can be increased, and assignment of a can be reduced, so that the solving result is more accurate.
Preferably, the phase-change heat exchange quantity Q based on the low-dimensional model surface temperature water is Q ═ f (T'W2) Specifically, the corrected surface temperature of the low-dimensional model is calculated by a phase-change heat exchange formula.
Preferably, the surface temperature of the anti-icing part of the low-dimensional model is obtained in an energy balance mode, and the calculation is accurate.
Preferably, the low-dimensional model is a two-dimensional calculation model, the high-dimensional model is a three-dimensional calculation model, and particularly, an Ansys solving control program can be adopted to solve the surface temperature of the anti-icing component.
As a specific embodiment, a system for solving a surface temperature of an anti-icing component through cross-dimensional coupling is shown in fig. 2, and includes a first surface temperature obtaining module 1, a second surface temperature obtaining module 2, a surface temperature correcting module 3, a heat exchange amount determining module 4, a surface temperature iteration module 5, and a convergence determining module 6.
The first surface temperature acquisition module 1 is used for acquiring the temperature of each node on the surface of the high-dimensional model; the second surface temperature obtaining module 2 is used for obtaining a surface temperature value of the low-dimensional model; the heat exchange quantity determining module 4 is used for obtaining the heat exchange quantity Q of the low-dimensional model; the surface temperature iteration module 5 is used for loading the heat exchange quantity Q on the high-dimensional model, completing the temperature calculation of each node on the surface of the high-dimensional model, and obtaining the temperature value of each node on the surface of the high-dimensional model in the iteration step; the convergence judging module 6 is used for comparing the temperature difference values of the latest round and the previous round of the surface nodes of the high-dimensional model, and if the difference value is smaller than the convergence judging standard, the temperature value of the surface of the latest round is the final calculation result; otherwise, repeating iterative operation of the surface temperature value of the high-dimensional model until the temperature difference value of each node on the surface of the high-dimensional model after two iterations is within the set threshold range.
The cross-dimension solving is carried out on the basis of the surface temperature, the influence of geometric parameters can be avoided, the heat flow item parameter represented by the surface temperature of the low-temperature surface and the heat exchange quantity represented by the surface temperature of the high-dimension surface are kept constant, and the cross-dimension coupling solving of the surface temperature of the anti-icing component is realized.
The above description is only for the specific embodiments of the present application, but the scope of the present application 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 application should be covered within the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.

Claims (5)

1. A method for solving the surface temperature of an anti-icing component through cross-dimensional coupling is characterized by comprising the following steps of: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
obtaining surface node temperature values T of the initial temperature field of the high-dimensional model of the anti-icing componentW1
Obtaining surface temperature T of anti-icing component of low-dimensional modelW2
Obtaining a correction value of the surface temperature of the low-dimensional model by adopting a low-dimensional model surface temperature correction formula;
determining a heat exchange quantity Q based on the water phase change of the low-dimensional model surface;
loading the surface heat exchange quantity Q based on the low-dimensional model on the high-dimensional model to complete the calculation of the surface temperature of the high-dimensional model, and obtaining the temperature value T 'of each node of the surface of the high-dimensional model in the iteration step'W1
Comparing the temperature difference value of the latest round and the previous round of the surface node of the high-dimensional model, and if the difference value is smaller than the convergence judgment standard, taking the temperature value of the surface of the latest round as the final calculation result; otherwise, repeating iterative operation of the surface temperature value of the high-dimensional model until the temperature difference value of each node on the surface of the high-dimensional model after two iterations is within the set threshold range.
2. The method of resolving ice protection component surface temperature across a dimensional coupling of claim 1, wherein: the correction formula of the surface temperature of the low-dimensional model is T'W2=a*TW1+b*TW2Wherein a + b is 1.
3. The method of resolving ice protection component surface temperature across a dimensional coupling of claim 1, wherein: the phase-change heat exchange quantity Q based on the low-dimensional model surface temperature water is Q ═ f (T'W2)。
4. The method of resolving ice protection component surface temperature across a dimensional coupling of claim 1, wherein: method for acquiring surface temperature T of anti-icing component of low-dimensional model in energy balance modeW2
5. A system for solving a surface temperature of an anti-icing component across a dimensional coupling, comprising: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
the first surface temperature acquisition module (1) is used for acquiring temperature values of all nodes on the surface of the high-dimensional model;
the second surface temperature acquisition module (2) is used for acquiring the surface temperature value of the low-dimensional model;
the surface temperature correction module (3) is used for correcting the surface temperature value of the low-dimensional model;
the heat exchange quantity determining module (4) is used for obtaining the heat exchange quantity Q of the low-dimensional model;
the surface temperature iteration module (5) is used for loading the heat exchange quantity Q on the high-dimensional model, completing the temperature calculation of each node on the surface of the high-dimensional model and obtaining the temperature value of each node on the surface of the high-dimensional model in the iteration step;
the convergence judging module (6) is used for comparing the temperature difference value of the latest round and the previous round of the surface nodes of the high-dimensional model, and if the difference value is smaller than the convergence judging standard, the temperature value of the surface of the latest round is the final calculation result; otherwise, repeating iterative operation of the surface temperature value of the high-dimensional model until the temperature difference value of each node on the surface of the high-dimensional model after two iterations is within the set threshold range.
CN202111230670.9A 2021-10-22 2021-10-22 Method for solving surface temperature of anti-icing component through cross-dimension coupling Active CN113887069B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111230670.9A CN113887069B (en) 2021-10-22 2021-10-22 Method for solving surface temperature of anti-icing component through cross-dimension coupling

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111230670.9A CN113887069B (en) 2021-10-22 2021-10-22 Method for solving surface temperature of anti-icing component through cross-dimension coupling

Publications (2)

Publication Number Publication Date
CN113887069A true CN113887069A (en) 2022-01-04
CN113887069B CN113887069B (en) 2022-08-19

Family

ID=79004194

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111230670.9A Active CN113887069B (en) 2021-10-22 2021-10-22 Method for solving surface temperature of anti-icing component through cross-dimension coupling

Country Status (1)

Country Link
CN (1) CN113887069B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115817822A (en) * 2023-02-09 2023-03-21 中国空气动力研究与发展中心低速空气动力研究所 Heat load distribution design method of electric heating anti-icing system

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108090246A (en) * 2017-11-14 2018-05-29 中国航发沈阳发动机研究所 A kind of anti-ice components Calculation Method of Temperature Field
CN112989727A (en) * 2021-05-10 2021-06-18 中国空气动力研究与发展中心低速空气动力研究所 Wall surface temperature simulation method of anti-icing system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108090246A (en) * 2017-11-14 2018-05-29 中国航发沈阳发动机研究所 A kind of anti-ice components Calculation Method of Temperature Field
CN112989727A (en) * 2021-05-10 2021-06-18 中国空气动力研究与发展中心低速空气动力研究所 Wall surface temperature simulation method of anti-icing system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
宋馨等: "防冰表面的对流换热计算分析", 《北京航空航天大学学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115817822A (en) * 2023-02-09 2023-03-21 中国空气动力研究与发展中心低速空气动力研究所 Heat load distribution design method of electric heating anti-icing system

Also Published As

Publication number Publication date
CN113887069B (en) 2022-08-19

Similar Documents

Publication Publication Date Title
CN113887069B (en) Method for solving surface temperature of anti-icing component through cross-dimension coupling
CN110543618A (en) roundness uncertainty evaluation method based on probability density function estimation
CN112380613B (en) One-dimensional three-dimensional joint simulation method for automobile engine cooling system
CN110225454B (en) Confidence transfer distributed type volume Kalman filtering cooperative positioning method
CN103237320B (en) Wireless sensor network quantizes based on mixing the method for tracking target that Kalman is merged
JP2009277226A (en) Method for real time calculation of process model and simulator therefor
CN106679745A (en) Temperature compensation method and system and water flow detection method and system for ultrasonic water meter
CN110032706A (en) A kind of the two stages method for parameter estimation and system of low order time lag system
CN111428335B (en) Joint simulation method and device for battery module
CN112697215B (en) Kalman filtering parameter debugging method for ultrasonic water meter data filtering
CN110567421B (en) Dynamic roundness uncertainty evaluation method based on Bayesian principle
CN110779544B (en) Double-task deep matching method for self-positioning and target positioning of multiple unmanned aerial vehicles
CN113158589A (en) Simulation model calibration method and device of battery management system
CN103092814A (en) Universal least square adjustment method used for measuring data
CN109598059B (en) Thermal protection system optimization design method and design system based on proxy model
CN108763611B (en) Wing structure random eigenvalue analysis method based on probability density evolution
Hendry et al. Optimal galaxy distance estimators
CN111259340A (en) Saturated load prediction method based on logistic regression
CN112329291B (en) Thermosetting coupling field calculation method of reflector antenna
CN116127638A (en) Method for generating three-dimensional solid grids of investment casting shell based on curvature method
CN117634278A (en) Design optimization method and design optimization device
WO2015070613A1 (en) Wireless positioning method and device
CN110866356A (en) Filter coupling matrix decoupling transformation method based on hybrid optimization algorithm
Theunissen et al. Autonomous spatially adaptive sampling in experiments based on curvature, statistical error and sample spacing with applications in LDA measurements
CN113361072A (en) Multi-disciplinary collaborative optimization method based on BLISS

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