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 PDFInfo
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling 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
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.
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)
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)
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 |
-
2021
- 2021-10-22 CN CN202111230670.9A patent/CN113887069B/en active Active
Patent Citations (2)
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)
Title |
---|
宋馨等: "防冰表面的对流换热计算分析", 《北京航空航天大学学报》 * |
Cited By (1)
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 |