CN114861134B - Step length determination method for calculating water drop motion track and storage medium - Google Patents

Step length determination method for calculating water drop motion track and storage medium Download PDF

Info

Publication number
CN114861134B
CN114861134B CN202210797044.6A CN202210797044A CN114861134B CN 114861134 B CN114861134 B CN 114861134B CN 202210797044 A CN202210797044 A CN 202210797044A CN 114861134 B CN114861134 B CN 114861134B
Authority
CN
China
Prior art keywords
time step
current time
value
water drop
error
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
CN202210797044.6A
Other languages
Chinese (zh)
Other versions
CN114861134A (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.)
Chengdu Jiangde Technology Co ltd
Original Assignee
Sichuan 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 Sichuan University filed Critical Sichuan University
Priority to CN202210797044.6A priority Critical patent/CN114861134B/en
Publication of CN114861134A publication Critical patent/CN114861134A/en
Application granted granted Critical
Publication of CN114861134B publication Critical patent/CN114861134B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • 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)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Algebra (AREA)
  • Evolutionary Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Measuring Volume Flow (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

The embodiment of the application discloses a step length determining method for calculating a water drop motion track and a storage medium, and relates to the technical field of water drop impact characteristic research. The method comprises the following steps: acquiring a single step error under the current time step according to the last time step of the water drop and the motion parameter variation of the current time step, wherein the motion parameter variation comprises position variation and/or speed variation; if the single step error under the current time step is smaller than the minimum value of the preset interval, increasing the value of the current time step as the value of the next time step; and if the single step error under the current time step is within a preset interval, taking the value of the current time step as the value of the next time step. And if the single step error under the current time step is larger than the maximum value of the preset interval, reducing the value of the current time step as the value of the next time step. Therefore, the problem that the motion trail of the water drop cannot meet the requirements of high precision and small calculated amount at the same time can be solved.

Description

Step length determination method for calculating water drop motion track and storage medium
Technical Field
The present disclosure relates to the field of research technologies for water drop impact characteristics, and more particularly, to a method for determining a step length for calculating a water drop motion trajectory and a storage medium.
Background
When the plane passes through the cloud layer, supercooled water drops impact the plane body, phase change may occur and icing may be caused. Icing can change the appearance and the streaming flow field of the airplane, destroy the pneumatic performance, reduce the maneuverability and the stability, threaten the flight safety and cause air crash accidents in severe cases. The water drop collection rate can quantitatively represent water drops impacting on an object and is an important parameter for researching the impact characteristics of the water drops. Therefore, the water drop collection rate can be used for researching the icing characteristic of the airplane and can also be used for designing an anti-icing system of the airplane.
At present, a numerical simulation method can be adopted to study the water drop collection rate, wherein the lagrangian method is one of the common methods for studying the water drop collection rate. And obtaining the motion trail of the water drop released by the far field plane by using a Lagrange method, thereby obtaining the condition that the water drop impacts the object plane. The water drop collection rate can then be obtained by using an area ratio method or a particle statistical method.
In the prior art, when a motion trajectory of a water drop released from a far-field plane is obtained, the same time step is usually adopted, if the accuracy of the obtained motion trajectory of the water drop is high enough to obtain the accuracy of the obtained water drop collection rate, the time step needs to be very small, and when the motion trajectory of the water drop is obtained, the number of calculated steps is very large, which results in a large calculation amount; if the calculation amount for obtaining the movement trajectory of the water droplet is to be reduced, the time step needs to be increased, but this will result in insufficient accuracy of the obtained movement trajectory of the water droplet.
Therefore, the prior art has the problem that high precision and small calculation amount cannot be simultaneously met when the motion trail of the water drop is obtained.
Disclosure of Invention
The application provides a step length determining method for calculating a water drop motion track, which is characterized by obtaining a single-step error under a current time step according to a motion parameter variation of the last time step of a water drop and a motion parameter variation of the current time step, wherein the motion parameter variation comprises a position variation and/or a speed variation, and the single-step error comprises a position variation error and/or a speed variation error; and if the single step error under the current time step is smaller than the minimum value of the preset interval, increasing the value of the current time step as the value of the next time step. And if the single step error under the current time step is larger than the maximum value of the preset interval, reducing the value of the current time step as the value of the next time step. And if the single step error under the current time step is within a preset interval, taking the value of the current time step as the value of the next time step. According to the method provided by the application, the value of the next time step is determined according to the error in each time step, when the error is small, namely the calculation precision of the water drop motion track can be ensured, the time step is increased, the times of calculating the water drop motion track are reduced, and therefore the calculation amount is reduced under the condition of ensuring the precision; when the error is large, the time step length is reduced, and the calculation precision of the movement track of the water drop is improved; and when the error is proper, keeping the original time step to continue the calculation. Therefore, the problem that in the prior art, when the motion trail of the water drop is obtained, high precision and small calculated amount cannot be met simultaneously can be effectively solved.
In a first aspect, an embodiment of the present application provides a step size determining method for calculating a movement trajectory of a water droplet, where the method includes: s110, obtaining a single step error under the current time step according to the motion parameter variation of the last time step of the water drop and the motion parameter variation of the current time step, wherein the motion parameter variation comprises position variation and/or speed variation, and the single step error comprises position variation error
Figure 196548DEST_PATH_IMAGE001
And/or speed variance error
Figure 595168DEST_PATH_IMAGE002
(ii) a S120, if the single step error under the current time step is smaller than the minimum value of the preset interval, increasing the value of the current time step as the value of the next time step, wherein the preset interval is
Figure 179733DEST_PATH_IMAGE003
Figure 35825DEST_PATH_IMAGE004
For the purpose of the upper limit of the single-step error,
Figure 844381DEST_PATH_IMAGE005
is an error control factor, and
Figure 402401DEST_PATH_IMAGE006
(ii) a And S130, if the single step error under the current time step is within a preset interval, taking the value of the current time step as the value of the next time step.
In a second aspect, the present application further provides a computer-readable storage medium, in which a program code is stored, where the program code can be called by a processor to execute the above method.
In summary, the present application has at least the following technical effects:
1. the value of the next time step is determined according to the error in each time step, when the error is small, namely the calculation precision of the water drop motion track can be guaranteed, the time step is increased, the number of times of calculating the water drop motion track is reduced, and therefore the calculation amount is reduced under the condition of guaranteeing the precision.
2. This application is through the error in according to every time step, confirms the value of next time step, and when the error is great, reduces the time step, improves the computational accuracy of water droplet movement track.
3. The value of the next time step is determined according to the error in each time step, when the error is appropriate, the original time step is kept to continue to be calculated, and meanwhile, high precision and small calculated amount are met.
Therefore, the scheme provided by the application can effectively solve the problem that the prior art cannot meet the requirements of high precision and small calculated amount at the same time when the motion trail of the water drop is obtained.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a schematic flow chart illustrating a step size determining method for calculating a water drop motion trajectory, provided by embodiment 1 of the present application;
fig. 2 shows a schematic diagram of a far field plane releasing water droplets provided in example 1 of the present application;
fig. 3 is a schematic diagram illustrating a last time step and a next time step provided in embodiment 1 of the present application;
FIG. 4 is a schematic diagram showing the current time step provided in embodiment 1 of the present application;
fig. 5 shows a storage unit for storing or carrying program codes for implementing the step size determining method for calculating a water drop motion trajectory according to the embodiment of the present application, provided in embodiment 2 of the present application.
Detailed Description
In order to make the technical solutions of the present application better understood, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
Currently, the lagrangian method is one of the common methods for studying the water droplet collection rate. And obtaining the motion trail of the water drop released by the far field plane by using a Lagrange method, thereby obtaining the condition that the water drop impacts the object plane. And then the water drop collection rate can be obtained by adopting an area ratio method or a particle statistical method. When the motion trail of the water drop released from the far-field plane is obtained, the same time step is usually adopted, if the precision of the obtained motion trail of the water drop is high enough to ensure that the precision of the obtained water drop collection rate is high enough, the time step needs to be very small, and when the motion trail of the water drop is obtained, the number of calculated steps is very large, so that the calculation amount is very large; if the calculation amount for obtaining the movement trajectory of the water droplet is to be reduced, the time step needs to be increased, but this will result in insufficient accuracy of the obtained movement trajectory of the water droplet.
Therefore, in order to solve the above-mentioned drawback, an embodiment of the present application provides a step size determining method for calculating a movement trajectory of a water droplet, where the method includes: obtaining a single step error under the current time step according to the motion parameter variation of the last time step of the water drop and the motion parameter variation of the current time step, wherein the motion parameter variation comprises a position variation and/or a speed variation, and the single step error comprises a position variation error and/or a speed variation error; if the single step error under the current time step is smaller than the minimum value of the preset interval, increasing the value of the current time step as the value of the next time step; and if the single step error under the current time step is within a preset interval, taking the value of the current time step as the value of the next time step. When the error is small, namely the calculation precision of the water drop motion track can be ensured, the time step length is increased, the frequency of calculating the water drop motion track is reduced, and when the error is appropriate, the original time step length is kept to continue calculation, so that the calculation amount is reduced under the condition of ensuring the precision.
The following describes a step length determination method for calculating a water drop motion trajectory according to the present application.
Example 1
Referring to fig. 1, fig. 1 is a schematic flowchart of a method for determining a step length for calculating a motion trajectory of a water droplet according to embodiment 1 of the present application. In this embodiment, the step length determining method for calculating the movement trajectory of the water drop may include the following steps:
step S110: obtaining a single step error under the current time step according to the motion parameter variation of the last time step and the motion parameter variation of the current time step of the water drop, wherein the motion parameter variation comprises a position variation and/or a speed variation, and the single step error comprises a position variation error
Figure 39925DEST_PATH_IMAGE001
And/or speed variation error
Figure 62107DEST_PATH_IMAGE002
In the embodiment of the present application, far field flat is shown in FIG. 2The Lagrange method judges whether each water drop can impact the object plane or not by obtaining the motion trail of each water drop released by the far field plane. The movement locus of the water drop can be obtained according to the stress condition of the water drop. Specifically, depending on the gravity, buoyancy and resistance to which the water droplets are subjected, it is possible to obtain:
Figure 244827DEST_PATH_IMAGE007
wherein, in the process,
Figure 572034DEST_PATH_IMAGE008
is the density of the air flow field and,
Figure 232823DEST_PATH_IMAGE009
is the density of the water droplets and,
Figure 578353DEST_PATH_IMAGE010
is the volume of the water droplets,
Figure 197554DEST_PATH_IMAGE011
is the displacement of the water droplet or droplets,
Figure 973574DEST_PATH_IMAGE012
is the velocity of the air flow field and,
Figure 438054DEST_PATH_IMAGE013
is the speed of the water droplets,
Figure 638091DEST_PATH_IMAGE014
is the area of the water drop facing the wind,
Figure 428192DEST_PATH_IMAGE015
is the acceleration of the force of gravity,
Figure 198833DEST_PATH_IMAGE016
is a water drop resistance coefficient, and
Figure 732583DEST_PATH_IMAGE017
Figure 583864DEST_PATH_IMAGE018
is the relative reynolds number for the gas,
Figure 528555DEST_PATH_IMAGE019
and is made of
Figure 301339DEST_PATH_IMAGE020
Figure 373200DEST_PATH_IMAGE021
Is the equivalent diameter of the water drops,
Figure 282251DEST_PATH_IMAGE022
is the air viscosity coefficient. From the above formula, one can obtain:
Figure 164887DEST_PATH_IMAGE023
wherein, in the step (A),
Figure 159388DEST_PATH_IMAGE024
is the acceleration of the water droplet. According to the acceleration of the water drop, the position of the water drop in each time step can be obtained until the water drop moves from a far-field plane to an object plane or moves from the far-field plane to the periphery of the object plane, and then the movement track of the water drop can be obtained.
As an alternative, the variation of the motion parameter includes a variation of a position or a variation of a speed, and the single step error is a variation of a position error
Figure 34940DEST_PATH_IMAGE025
Or error in variation of speed
Figure 64076DEST_PATH_IMAGE026
As another alternative, the motion parameter variation includes a position variation and a speed variation, and the single-step error is a position variation error
Figure 350570DEST_PATH_IMAGE025
And velocityError of variation
Figure 97946DEST_PATH_IMAGE026
As still another alternative, the motion parameter variation includes a position variation and a speed variation, and the single-step error is a position variation error
Figure 511610DEST_PATH_IMAGE025
And error in speed variation
Figure 129673DEST_PATH_IMAGE026
Maximum value of (2).
In the embodiment of the present application,
Figure 88533DEST_PATH_IMAGE025
may be a relative error of the amount of position change, i.e.
Figure 57626DEST_PATH_IMAGE027
Figure 274980DEST_PATH_IMAGE025
It may also be the absolute error of the amount of change in position, i.e.
Figure 13129DEST_PATH_IMAGE028
Wherein, in the step (A),
Figure 907005DEST_PATH_IMAGE029
is the amount of position change of the current time step,
Figure 160132DEST_PATH_IMAGE030
the position variation of the last time step is used as the position variation of the last time step.
In the embodiment of the present application,
Figure 915598DEST_PATH_IMAGE026
may be a relative error in the amount of speed change, i.e.
Figure 524565DEST_PATH_IMAGE031
Figure 340074DEST_PATH_IMAGE026
It may also be the absolute error of the speed variation, i.e.
Figure 18180DEST_PATH_IMAGE032
Wherein, in the process,
Figure 577338DEST_PATH_IMAGE033
is the speed variation of the current time step,
Figure 279626DEST_PATH_IMAGE034
the speed variation of the last time step is adopted.
In an exemplary embodiment, if the last time step is followed by the next time step, step S110 may include sub-step S111 and sub-step S112.
Substep S111: and obtaining the acceleration of the last time step of the water drop and the acceleration in the current time step after the last time step of the water drop according to the stress condition of the water drop, wherein the current time step is equal to the current time step.
Substep S112: and obtaining the motion parameter variation of the previous time step according to the acceleration of the previous time step on the water drop and the value of the previous time step, and obtaining the motion parameter variation of the current time step according to the acceleration in the current time step after the previous time step on the water drop and the value of the current time step.
In the embodiment of the present application, as shown in fig. 3, the last time step has a value of
Figure 266037DEST_PATH_IMAGE035
The next time step has a value of
Figure 697018DEST_PATH_IMAGE036
. The last time step is followed by the next time step. In this embodiment, the current time stepThe length is not the time step actually used in the process of acquiring the motion track of the water drop, but is an assumed value used for determining the next time step, and the value of the current time step is equal to the value of the previous time step.
In an exemplary embodiment, if the last time step is followed by the current time step and the current time step is followed by the next time step, step S110 may further include substeps S113 and S114.
Substep S113: and obtaining the acceleration of the last time step of the water drop and the acceleration of the current time step according to the stress condition of the water drop, wherein the current time step is determined by the single step error of the last time step.
Substep S114: and obtaining the motion parameter variation of the previous time step according to the acceleration of the previous time step of the water drop and the value of the previous time step, and obtaining the motion parameter variation of the current time step according to the acceleration of the current time step of the water drop and the value of the current time step.
In the embodiment of the present application, as shown in fig. 4, the last time step has a value of
Figure 794287DEST_PATH_IMAGE035
The current time step has a value of
Figure 112267DEST_PATH_IMAGE037
The next time step has a value of
Figure 269579DEST_PATH_IMAGE036
. The last time step is followed by the current time step, and the next time step is followed by the current time step. In this embodiment, the current time step is a time step actually used in the process of obtaining the movement trajectory of the water drop, and the current time step is determined by a single-step error in the previous time step, if the single-step error in the previous time step is within a preset interval, the current time step is equal to the previous time step, and if the single-step error in the previous time step is not within the preset interval, the current time step is not within the preset intervalThe time step may not be equal to the previous time step, but the value of the previous time step is decreased or increased according to the single step error as the value of the current time step.
Step S120: if the single step error under the current time step is smaller than the minimum value of the preset interval, increasing the value of the current time step as the value of the next time step, wherein the preset interval is
Figure 922277DEST_PATH_IMAGE003
Figure 823237DEST_PATH_IMAGE004
For the purpose of the upper limit of the single-step error,
Figure 759838DEST_PATH_IMAGE005
is an error control factor, and
Figure 88051DEST_PATH_IMAGE006
in an exemplary embodiment, the value of the upper limit of the single step error and the value of the error control factor may be set empirically, for example, the value of the upper limit of the single step error may be
Figure 24783DEST_PATH_IMAGE038
The value of the error control factor may be 1.3.
As an alternative, the current time step value may be increased empirically.
As another alternative, the value of the current time step may be increased by using a binary method.
For example: if the last time step is later than the next time step, setting
Figure 214587DEST_PATH_IMAGE039
Figure 756427DEST_PATH_IMAGE040
Figure 255541DEST_PATH_IMAGE041
. The last time step has a value of
Figure 882832DEST_PATH_IMAGE035
Let the current time step have a value of
Figure 374862DEST_PATH_IMAGE042
And at this time
Figure 36787DEST_PATH_IMAGE043
If the single-step error in the current time step is smaller than the minimum value of the preset interval, the implementation manner may be: order to
Figure 706803DEST_PATH_IMAGE044
Figure 555811DEST_PATH_IMAGE045
Then give an order again
Figure 352996DEST_PATH_IMAGE046
Increase the current time step to
Figure 869428DEST_PATH_IMAGE047
And is used as the value of the next time step; the implementation may also be: order to
Figure 444766DEST_PATH_IMAGE048
Figure 781070DEST_PATH_IMAGE049
Then give an order again
Figure 614902DEST_PATH_IMAGE050
Increasing the value of the current time step to
Figure 985841DEST_PATH_IMAGE051
And at this time
Figure 794397DEST_PATH_IMAGE052
Repeating the steps until the single step error under the increased current time step is in a preset interval, and increasing the current time step to the current time step
Figure 634308DEST_PATH_IMAGE051
As the value of the next time step.
Another example is: if the last time step is followed by the current time step and the next time step is followed by the current time step, setting
Figure 756985DEST_PATH_IMAGE053
Figure 248009DEST_PATH_IMAGE054
Figure 430729DEST_PATH_IMAGE055
. Current step of time and that time
Figure 8470DEST_PATH_IMAGE051
If the single step error in the current time step is smaller than the minimum value of the preset interval, the embodiment may be: order to
Figure 669258DEST_PATH_IMAGE056
Figure 14789DEST_PATH_IMAGE057
Then give an order again
Figure 368410DEST_PATH_IMAGE058
Increase the current time step to
Figure 917334DEST_PATH_IMAGE059
And is used as the value of the next time step; the implementation may also be: order to
Figure 647393DEST_PATH_IMAGE060
Figure 847430DEST_PATH_IMAGE061
Then give an order again
Figure 371952DEST_PATH_IMAGE062
Increasing the value of the current time step to
Figure 906708DEST_PATH_IMAGE051
And at this time
Figure 440457DEST_PATH_IMAGE063
Repeating the above steps until the single step error of the increased current time step is within a preset interval, and increasing the current time step to the current time step
Figure 495001DEST_PATH_IMAGE051
As the value of the next time step.
The value of the next time step is determined according to the error in each time step, when the error is small, namely the calculation precision of the water drop motion track can be guaranteed, the time step is increased, the number of times of calculating the water drop motion track is reduced, and therefore the calculation amount is reduced under the condition of guaranteeing the precision.
Step S130: and if the single step error under the current time step is within a preset interval, taking the value of the current time step as the value of the next time step.
The value of the next time step is determined according to the error in each time step, when the error is appropriate, the original time step is kept to continue to be calculated, and meanwhile, high precision and small calculated amount are met.
Step S140: and if the single step error under the current time step is larger than the maximum value of the preset interval, reducing the value of the current time step as the value of the next time step.
As an alternative, the current time step value may be reduced by empirically reducing the current time step value.
As another alternative, the current time step value may be decreased by a binary method. The step of decreasing the value of the current time step by the bisection method may refer to the step of increasing the value of the current time step by the bisection method, which is not described in detail herein.
This application is through the error in according to every time step, confirms the value of next time step, and when the error is great, reduces the time step, improves the computational accuracy of water droplet movement track.
As an alternative embodiment, step S110, step S120, step S130 and step S140 are performed once every n time steps have elapsed, wherein,
Figure 190424DEST_PATH_IMAGE064
in the present embodiment, the first time step during the movement of the water droplet may be determined empirically.
In the embodiment of the present application, the first time step in the movement process of the water drop can also be determined by using a dichotomy.
For example, set
Figure 713941DEST_PATH_IMAGE065
Figure 785802DEST_PATH_IMAGE066
Figure 694852DEST_PATH_IMAGE067
Figure 826756DEST_PATH_IMAGE068
Can be any time value set according to experience, and the value of the first time step is
Figure 867262DEST_PATH_IMAGE069
And at this time
Figure 742814DEST_PATH_IMAGE070
Let us orderThe value of two time steps is also
Figure 506371DEST_PATH_IMAGE069
And obtaining the motion parameter variation of the first time step according to the value of the first time step and the acceleration of the water drop in the first time step, obtaining the motion parameter variation of the second time step according to the value of the second time step and the acceleration of the water drop in the second time step, and obtaining the single-step error under the first time step according to the motion parameter variation of the first time step and the motion parameter variation of the second time step.
If the single step error is less than the minimum value of the preset interval, the order is given
Figure 559909DEST_PATH_IMAGE071
Figure 41706DEST_PATH_IMAGE072
Then give a new order
Figure 455369DEST_PATH_IMAGE073
Repeating the steps until the single step error under the increased first time step is within a preset interval, and adding the increased single step error at the moment
Figure 339012DEST_PATH_IMAGE069
As the value of the first time step.
If the single step error is within a preset interval, the error will be detected
Figure 61986DEST_PATH_IMAGE069
As the value of the first time step, and this time
Figure 31079DEST_PATH_IMAGE074
If the single step error is larger than the minimum value of the preset interval, order
Figure 982855DEST_PATH_IMAGE075
Figure 721003DEST_PATH_IMAGE076
Then give an order again
Figure 116344DEST_PATH_IMAGE077
Repeating the steps until the single step error under the reduced first time step is within a preset interval, and reducing the reduced single step error at the moment
Figure 572733DEST_PATH_IMAGE069
As the value of the first time step.
The method for determining the step length for calculating the movement track of the water drops is suitable for calculating the water drop collection rate on the airplane and also suitable for calculating the water drop collection rate on other equipment such as a wind driven generator.
Example 2
Referring to fig. 5, fig. 5 is a block diagram illustrating a structure of a computer-readable storage medium according to embodiment 2 of the present application. The computer-readable storage medium 600 has stored therein program code that can be called by a processor to execute the method described in the above-described method embodiments.
The computer-readable storage medium 600 may be an electronic memory such as a flash memory, an EEPROM (electrically erasable programmable read only memory), an EPROM (erasable programmable read only memory), a hard disk, or a ROM. Alternatively, the computer-readable storage medium 600 includes a non-volatile computer-readable storage medium. The computer readable storage medium 600 has storage space for program code 610 for performing any of the method steps of the method described above. The program code can be read from or written to one or more computer program products. The program code 610 may be compressed, for example, in a suitable form.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solutions of the present application, and not to limit the same; although the present application has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; such modifications and substitutions do not necessarily depart from the spirit and scope of the corresponding technical solutions in the embodiments of the present application.

Claims (6)

1. A step length determination method for calculating a movement locus of a water drop is characterized by comprising the following steps:
s110, obtaining a single step error under the current time step according to the motion parameter variation of the last time step of the water drop and the motion parameter variation of the current time step, wherein the motion parameter variation comprises position variation and/or speed variation, and the single step error comprises position variation error
Figure 585798DEST_PATH_IMAGE001
And/or speed variance error
Figure 991503DEST_PATH_IMAGE002
S120, if the single step error under the current time step is smaller than the minimum value of a preset interval, increasing the value of the current time step as the value of the next time step, wherein the preset interval is
Figure 541433DEST_PATH_IMAGE003
Figure 485118DEST_PATH_IMAGE004
For the purpose of the upper limit of the single-step error,
Figure 293543DEST_PATH_IMAGE005
is an error control factor, and
Figure 6284DEST_PATH_IMAGE006
s130, if the single step error under the current time step is within the preset interval, taking the value of the current time step as the value of the next time step;
s140, if the single step error under the current time step is larger than the maximum value of the preset interval, reducing the value of the current time step as the value of the next time step.
2. The method for determining step length for calculating trajectory of water droplet according to claim 1, wherein if the previous time step is followed by the next time step, step S110 further comprises:
obtaining the acceleration of the last time step of the water drop and the acceleration in the current time step after the last time step of the water drop according to the stress condition of the water drop, wherein the value of the current time step is equal to the value of the last time step;
and obtaining the motion parameter variation of the last time step according to the acceleration of the last time step of the water drop and the value of the last time step, and obtaining the motion parameter variation of the current time step according to the acceleration of the water drop in the current time step after the last time step and the value of the current time step.
3. The method as claimed in claim 1, wherein if the previous time step is followed by the current time step, and the current time step is followed by the next time step, step S110 further includes:
obtaining the acceleration of the last time step of the water drop and the acceleration of the current time step according to the stress condition of the water drop, wherein the current time step is determined by the single step error under the last time step;
and obtaining the motion parameter variation of the previous time step according to the acceleration of the previous time step of the water drop and the value of the previous time step, and obtaining the motion parameter variation of the current time step according to the acceleration of the current time step of the water drop and the value of the current time step.
4. The method as claimed in any one of claims 1 to 3, wherein if the motion parameter variation includes a position variation and a speed variation, the single step error is a position variation error
Figure 992695DEST_PATH_IMAGE001
And error in speed variation
Figure 892518DEST_PATH_IMAGE002
Maximum value of (2).
5. The step size determining method for calculating a trajectory of a water droplet of claim 1, wherein the steps S110, S120, S130 and S140 are performed every n time steps, wherein,
Figure 740519DEST_PATH_IMAGE007
6. a computer-readable storage medium, characterized in that a program code is stored in the computer-readable storage medium, which program code can be called by a processor to execute the method of any of claims 1-5.
CN202210797044.6A 2022-07-08 2022-07-08 Step length determination method for calculating water drop motion track and storage medium Active CN114861134B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210797044.6A CN114861134B (en) 2022-07-08 2022-07-08 Step length determination method for calculating water drop motion track and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210797044.6A CN114861134B (en) 2022-07-08 2022-07-08 Step length determination method for calculating water drop motion track and storage medium

Publications (2)

Publication Number Publication Date
CN114861134A CN114861134A (en) 2022-08-05
CN114861134B true CN114861134B (en) 2022-09-06

Family

ID=82626909

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210797044.6A Active CN114861134B (en) 2022-07-08 2022-07-08 Step length determination method for calculating water drop motion track and storage medium

Country Status (1)

Country Link
CN (1) CN114861134B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101178747A (en) * 2007-12-18 2008-05-14 东北大学 Method for forecasting transient state temperature field with S type step length changing method in the process of plate belt hot rolling
CN103227623A (en) * 2013-03-29 2013-07-31 北京邮电大学 Step value-variable LMS (Least Mean Square) self-adaptation filtering algorithm and filter
CN107272597A (en) * 2017-07-14 2017-10-20 福建工程学院 A kind of nurbs curve interpolation based on advance and retreat method quickly pre-reads processing method
CN108052783A (en) * 2018-01-29 2018-05-18 济南大学 A kind of unsaturated soil dynamic numerical calculation method based on adaptive step
CN108733622A (en) * 2018-05-31 2018-11-02 南京信息工程大学 Ionized layer TEC single station forecasting method based on auto-correlation unequal interval iterative method
CN111396269A (en) * 2020-06-08 2020-07-10 中国空气动力研究与发展中心低速空气动力研究所 Multi-time-step unsteady icing calculation method and system and storage medium

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140257770A1 (en) * 2011-11-30 2014-09-11 Ming Lu Numerical simulation method for the flight-icing of helicopter rotary-wings
US20140257771A1 (en) * 2011-11-30 2014-09-11 Ming Lu Numerical simulation method for aircrasft flight-icing

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101178747A (en) * 2007-12-18 2008-05-14 东北大学 Method for forecasting transient state temperature field with S type step length changing method in the process of plate belt hot rolling
CN103227623A (en) * 2013-03-29 2013-07-31 北京邮电大学 Step value-variable LMS (Least Mean Square) self-adaptation filtering algorithm and filter
CN107272597A (en) * 2017-07-14 2017-10-20 福建工程学院 A kind of nurbs curve interpolation based on advance and retreat method quickly pre-reads processing method
CN108052783A (en) * 2018-01-29 2018-05-18 济南大学 A kind of unsaturated soil dynamic numerical calculation method based on adaptive step
CN108733622A (en) * 2018-05-31 2018-11-02 南京信息工程大学 Ionized layer TEC single station forecasting method based on auto-correlation unequal interval iterative method
CN111396269A (en) * 2020-06-08 2020-07-10 中国空气动力研究与发展中心低速空气动力研究所 Multi-time-step unsteady icing calculation method and system and storage medium

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Numerical simulation for three-dimensional rotor icing in forward flight;Zhengzhi Wang 等;《Advances in Mechanical Engineering》;20180423;1-12 *
二元翼型水滴收集率计算研究;周 峰 等;《民用飞机设计与研究》;20081231;15-18,33 *
翼型积冰的数值模拟;易 贤 等;《空气动力学学报》;20021231;第20卷(第4期);428-433 *

Also Published As

Publication number Publication date
CN114861134A (en) 2022-08-05

Similar Documents

Publication Publication Date Title
CN101059821B (en) Simulation method and simulator
CN107451319B (en) Modeling method of space debris environment long-term evolution model
CN114861134B (en) Step length determination method for calculating water drop motion track and storage medium
US8452577B2 (en) Golf ball trajectory simulation method
CN114118537A (en) Combined prediction method for carbon emission of airspace flight
CN115292656A (en) Aircraft ice accretion prediction method and device based on fuzzy logic
CN115983160A (en) Lagrange method-based water drop trajectory simulation calculation method
CN103921957B (en) The energy management method of jumping up that a kind of lunar exploration airship great-jump-forward reenters
WO2019186151A1 (en) Methods and apparatus for simulating liquid collection on aerodynamic components
Rudnick-Cohen et al. Modeling unmanned aerial system (UAS) risks via Monte Carlo simulation
EP3074713B1 (en) System integration
CN116562192B (en) Method, device, equipment and storage medium for predicting icing ice shape of airplane
CN103760769A (en) Small unmanned aerial vehicle control object modeling method based on test data
CN111259337B (en) Heavy debris real-time drop point forecasting method based on statistics
CN115081121B (en) Icing simulation method considering stream phenomenon and storage medium
CN114556437A (en) Method of generating a grid of parts, method of using a grid of parts, computer program and computer readable medium
KR102489288B1 (en) Meteorological prediction apparatus and method capable of adjusting the parameters of cloud microphysics, and recorded media recording program to perform the same
WO2021213541A1 (en) Method for assessing personal injury risk of unmanned aerial vehicle crash in out-of-control or power-loss fault state
CN115099048A (en) Messinger model-based icing numerical simulation method, system, equipment and medium
Pueyo et al. An Eulerian Approach with Mesh Adaptation for Highly Accurate 3D Droplet Dynamics Simulations
CN108280043A (en) A kind of method and system of fast prediction flight path
Sóbester et al. Notes on meteorological balloon mission planning
CN109003470B (en) Method and device for monitoring and alarming track consistency
KR101418490B1 (en) Analysis method of droplet interaction based on eulerian for impingement prediction of supercooled large droplets
CN104657595A (en) Coefficient calibrating method for single particle drag model

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
TR01 Transfer of patent right

Effective date of registration: 20231011

Address after: No. 888, Section 2, Chenglong Avenue, Chengdu Economic and Technological Development Zone, 610000, Sichuan

Patentee after: CHENGDU JIANGDE TECHNOLOGY Co.,Ltd.

Address before: 610000, No. 24, south section of Ring Road, Sichuan, Chengdu

Patentee before: SICHUAN University

TR01 Transfer of patent right