CN110705029A - Flow field prediction method of oscillating flapping wing energy acquisition system based on transfer learning - Google Patents
Flow field prediction method of oscillating flapping wing energy acquisition system based on transfer learning Download PDFInfo
- Publication number
- CN110705029A CN110705029A CN201910838392.1A CN201910838392A CN110705029A CN 110705029 A CN110705029 A CN 110705029A CN 201910838392 A CN201910838392 A CN 201910838392A CN 110705029 A CN110705029 A CN 110705029A
- Authority
- CN
- China
- Prior art keywords
- domain
- information
- fluid
- flow field
- data
- 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
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/044—Recurrent networks, e.g. Hopfield networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for predicting a flow field of an oscillating flapping wing energy acquisition system based on transfer learning. The method comprises the steps of establishing a generation and discrimination network of generation type antagonistic neural network transfer learning by acquiring and preprocessing known fluid domain information, realizing fluid domain information prediction under different parameter values through training, and realizing continuous prediction of the fluid domain information along with time by taking the predicted fluid domain information as the known fluid domain information; the flow field is restored by the predicted fluid domain information, the flow field prediction is realized, and fluid dynamics parameters and performance parameters at the prediction time point can be calculated. The method can accurately predict the flow field of the oscillating flapping wing energy acquisition system and obtain fluid dynamics parameters, performance parameters and the like while greatly reducing time and economic cost, and has important significance for the research and development of the oscillating flapping wing energy acquisition system.
Description
Technical Field
The invention belongs to the technical field of fluid dynamics, and particularly relates to a flow field prediction method of an oscillating flapping wing energy acquisition system based on transfer learning.
Background
With the rapid development of various social industries, the world energy situation is increasingly severe, and the exploration and application of renewable energy sources gradually become the development theme of various countries all over the world. The oscillating flapping wing energy acquisition system has a huge development prospect due to the fact that the oscillating flapping wing energy acquisition system is simple in structure and can effectively utilize low-quality renewable energy sources such as wind energy and ocean energy, and the oscillating flapping wing energy acquisition system becomes a research hotspot in recent years.
In the process of researching, designing and developing the oscillating flapping wing energy collecting system, two basic methods of numerical simulation and experimental research are mainly relied on at the present stage. The actual working condition of the oscillating flapping wing energy collecting system is determined by a complex external environment, so that the working conditions are numerous, and a large amount of numerical simulation or experimental research is bound to be carried out in the process of researching, designing and developing the oscillating flapping wing energy collecting system. For the numerical simulation method, because the calculation of the oscillating flapping wing energy collecting system belongs to unsteady flow field calculation, a large amount of calculation by using the traditional commercial CFD software brings huge time cost. For the experimental research method, because the actual working conditions of the oscillating flapping wing energy collecting system are complicated, different experimental conditions are needed for the simulation of different working conditions, and the experimental research method brings huge economic cost.
In recent years, the field of artificial intelligence rises rapidly, and simultaneously along with the rapid development of the field of computers, the computing capability is greatly improved, so that new solutions are brought to many traditional fields, and a new revolution of revolution is brought to many traditional industries. For the field of fluid dynamics, the application of artificial intelligence is increasingly wide, and great progress is made from single-parameter prediction to multi-parameter prediction of flow performance and from steady-state flow field prediction to unsteady-state flow field prediction of a flow field state.
Disclosure of Invention
The invention aims to provide a flow field prediction method of an oscillating flapping wing energy acquisition system based on transfer learning, which can greatly reduce time cost and economic cost while ensuring the accuracy of a predicted flow field.
The invention is realized by adopting the following technical scheme:
a method for predicting a flow field of an oscillating flapping wing energy acquisition system based on transfer learning is disclosed, wherein unknown fluid domain information of a next time point is predicted from known fluid domain information of a plurality of continuous time points after time dispersion; specifically, by acquiring node position information and node flow field parameter information of a fluid domain after spatial dispersion, performing matrixing processing and standardization processing, distributing a source domain and a target domain according to a set proportion, further extracting input data sets and corresponding theoretical output data sets of the source domain and the target domain, then respectively establishing a generation network and a discrimination network of generation type antagonistic neural network transfer learning, after training, realizing fluid domain information prediction under different parameter values, and taking unknown fluid domain information obtained by prediction as new known fluid domain information, namely realizing continuous prediction of the fluid domain information along with time under different parameter values; the whole flow field is restored through node position information contained in the fluid domain information obtained through prediction and flow field parameter information on the nodes, flow field prediction is achieved, fluid dynamics parameters and performance parameters of the oscillating flapping wing energy collecting system at the prediction time point are calculated, and therefore researchers can analyze the oscillating flapping wing energy collecting system conveniently.
The invention is further improved in that the method specifically comprises the following steps:
1) obtaining model data of oscillating flapping wing energy collecting system
Respectively collecting known fluid domain information of a plurality of continuous time-dispersed time points of the model aiming at different parameter values, wherein the known fluid domain information comprises node position information and node flow field parameter information of the fluid domain after space dispersion; carrying out space discretization on the model to obtain grid information Hj,qPerforming unsteady state numerical simulation calculation, calculating an application layer flow model, spatially adopting a second-order difference format, temporally adopting a first-order difference format, solving by using a dynamic grid technology, and obtaining node position information at each time point according to a solving resultAnd flow field parameter information on the nodeJ is 1, 2.. the J, J is the number of grid units, Q is 1, 2.. the Q, Q is the number of nodes of each grid unit, R is 1, 2.. the R, R is the total number of different selected parameter values, T is 1, 2.. the T, T is the number of time points after time dispersion, I is 1, 2.. the I, I is the total number of nodes, N is 1, 2.. the N, N is the number of coordinate dimensions, M is 1, 2.. the M, M is the number of flow field parameters to be considered;
2) model data preprocessing of oscillating flapping wing energy collecting system
Regular matrixing is carried out on the model data of the oscillating flapping wing energy acquisition system, and then node position information is rearrangedAnd flow field parameter information on the nodeCorresponding elements in the node are transformed to obtain the position information of the nodes after matrixingAnd flow field parameter information on the nodeWherein a is 1,2, a, B is 1,2, a, B, C is 1,2, C, A, B, C are the number of nodes in three mutually orthogonal directions, respectively, and a × B × C is satisfied;
for the node position information after matrixingAnd flow field parameter information on the nodePerforming standardization treatment to obtainNode position information G after normalization processingr,t,a,b,c,nAnd flow field parameter information F on the noder,t,a,b,c,m
3) Model data splitting source domain and target domain of oscillating flapping wing energy acquisition system
The normalized node position information G corresponding to different values selected by the parametersr,t,a,b,c,nAnd flow field parameter information F on the noder,t,a,b,c,mSplitting at random or according to the requirement according to the set proportion, using one part as the source domain and the node position information of the source domainAnd flow field parameter information on the nodeAnother part is used as target domain, target domain node position informationAnd flow field parameter information on the nodeRS, RS are different total numbers of values selected as parameters of a source domain, RT is 1,2,., RT are different total numbers of values selected as parameters of a target domain, and RS + RT is satisfied;
4) extracting respective input data set and theoretical output data set in source domain and target domain
The number N of continuous time points required according to actual input dataTBCarrying out further data extraction on the source domain and the target domain to obtain respective input data set and theoretical output data set;
for the source domain, the node position information of the source domain is firstlyAnd flow field parameter information on the nodeSplicing to obtain the total data information of the source domainWhere K1, 2,., K N + M, a source domain Input dataset Input is obtainedSDAnd a corresponding source domain theoretical output data set LableSD;
For the target domain, firstly, the node position information of the target domainAnd flow field parameter information on the nodeSplicing to obtain the total data information of the target domainWhere K is 1,2,., K is N + M, the target domain Input dataset Input is obtainedTDAnd a corresponding target domain theoretical output data set LableTD;
5) Fluid domain information prediction GAN migration learning network for constructing oscillating flapping wing energy acquisition system
The fluid domain information prediction GAN migration learning network of the oscillating flapping wing energy acquisition system comprises a generation network and a judgment network, wherein the generation network is called G Net for short, and the judgment network is called D Net for short;
firstly, a convolution residual error network is adopted to establish a characteristic extraction part of G Net, and Input data Input from a source domain is realizedSDOr target field Input data InputTDExtracting corresponding source domain fluid domain Feature information FeatureSDOr target domain fluid domain Feature information FeatureTDThen, a deconvolution residual error network is adopted to establish a predicted fluid domain information part of G Net, and the Feature information Feature of the source domain fluid domain is realizedSDOr target domain fluid domain Feature information FeatureTDGenerating corresponding source domain actual Output data OutputSDOr target domain actual Output data OutputTD(ii) a Wherein each source domain fluid domain characteristic informationES, each source domain actually outputs data asES, characteristic information of each target domain fluid domainET 1,2, ET, each target field actually outputting data as ET 1, 2.. ET, G Net total 21 layers of network; secondly, a D Net network is established by adopting a convolution residual error network to judge whether the fluid domain information is true or false, and the actual output data of each source domain is realizedAnd each source domain theoretical output data corresponding to the source domain theoretical output dataTo obtain the corresponding discrimination resultIs 0 and1, the data sets of the discrimination result are respectively expressed as D _ OutputSDAnd D _ LableSDSimultaneously undertake the actual output of data from each target domainObtaining a discrimination resultThe corresponding data set of the discrimination result is D _ OutputTDThe two-classification loss function part is used for calculating the network training of G Net, and D Net has 19 layers of networks in total;
with the quantity SSDA batch of source domain input data and quantity STDA batch of target domain input data for calculating maximum mean difference loss;
g Net in GAN network adopts two-norm loss function L2Loss of class, two-class loss function LbceGLoss function of loss and maximum mean difference LmmdA weighted average of _ loss as a total loss function;
d Net in GAN network adopts two-classification loss of source domain actual output dataTwo-classification loss L of source domain theoretical output databceLbLoss of class two for the _lossand actual output data of the target domainThe sum as a function of total loss;
6) fluid domain information prediction GAN migration learning network for training oscillating flapping wing energy acquisition system
Firstly, the number of the oscillating flapping wing energy acquisition systems is SSDThe fluid domain information of a batch of source domain input data is predicted through G Net in a GAN transfer learning network, and the obtained quantity is SSDA batch of source domain fluid domain characteristic information and a batch of source domain actual output data, and then the number of the oscillating flapping wing energy acquisition systems is STDThe fluid domain information prediction is carried out on a batch of target domain input data for calculating the maximum mean difference loss through G Net in a GAN migration learning network, and the obtained quantities are STDA batch of target domain fluid domain characteristic information and a batch of target domain actual output data;
then the obtained quantity is SSDA batch of source domain actual output data, the corresponding quantity is SSDA batch of source domain theoretical output data and the quantity is STDThe actual output data of a batch of target domains are respectively distinguished by DNet in the GAN migration learning network, and the D Net is subjected to parameter updating once according to the obtained distinguishing result, so that the D Net can better distinguish the prediction of the G NetThe fluid domain information prediction result and the fluid domain information theoretical result are obtained;
then the obtained quantity is SSDA batch of source domain actual output data and the number STDThe actual output data of a batch of target domains are respectively judged through D Net which is subjected to parameter updating once in the GAN migration learning network, and according to the obtained judgment result and the obtained quantity are SSDA batch of source domain fluid domain characteristic information and a batch of source domain actual output data, the quantity is STDIs known to correspond to the actual output of the source domain, and has a quantity SSDThe G Net is subjected to parameter updating once by a batch of source domain theoretical output data, so that the fluid domain information prediction result of the G Net is closer to the fluid domain information theoretical result, and the judgment of the fluid domain information prediction result and the fluid domain information theoretical result of the D Net is confused as much as possible;
through the continuous and repeated confrontation of G Net and D Net, the judgment of the fluid domain information prediction result and the fluid domain information theoretical result by the D Net is more accurate, and the fluid domain information prediction result of the G Net is forced to be closer to the fluid domain information theoretical result;
7) flow field prediction result processing and analysis of oscillating flapping wing energy acquisition system
The trained fluid domain information prediction GAN migration learning network of the oscillating flapping wing energy acquisition system is used for predicting fluid domain information under other parameter values, and the predicted unknown fluid domain information at the next time point is used as new known fluid domain information, so that the fluid domain information prediction at the next time point can be continuously carried out, and the continuous prediction of the fluid domain information along with time is realized under different parameter values; the whole flow field can be restored through node position information and flow field parameter information on the nodes contained in the fluid domain information obtained through prediction; furthermore, fluid dynamics parameters and performance parameters such as torque, lift coefficient and resistance coefficient of the oscillating flapping wing energy collecting system at a prediction time point are calculated according to a flow field prediction result so as to provide convenience for researchers to analyze the oscillating flapping wing energy collecting system.
The invention is further improved in that in step 2), the position information of the nodes after matrixing is processedAnd flow field parameter information on the nodeThe normalization process was performed as follows:
for the node position information after matrixingThe standardization processing method comprises the following steps:
wherein, MaxGnRepresents the maximum value of the nth coordinate dimension of the node in all data, wherein N is 1,2nRepresenting the minimum value of the nth coordinate dimension of the node in all data of the nth coordinate dimension;
to the parameter information of the flow field on the nodes after matrixingThe standardization processing method comprises the following steps:
wherein, MaxFmRepresents the maximum value of the mth parameter of the node in all data, wherein M is 1,2mRepresents the minimum value of the mth parameter of the node in all data thereof.
In a further development of the invention, in step 4), the source domain Input dataset Input is obtainedSDAnd a corresponding source domain theoretical output data set LableSDThe method comprises the following steps:
the source domain input dataset is:
wherein each source domain input data isES, which is the total amount of source domain data and satisfies
ES=RS×(T-NTB)
The corresponding source domain theoretical output data set is:
wherein each source domain theoretical output data isES, which is the total amount of source domain data and satisfies ES RS × (T-N)TB)。
The invention further improves that in the step 4), a target domain Input data set Input is obtainedTDAnd a corresponding target domain theoretical output data set LableTDThe method comprises the following steps:
the target domain input dataset is:
wherein each target domain inputs data as ET 1,2, ET is the total amount of target domain data, and there are multiple time period data missing, i.e. ET RT × (T-N) may not be satisfiedTB);
The corresponding target domain theoretical output data set is:
wherein each target domain theoretical output data is ET 1, 2.. ET, ET is the total amount of target domain data, and there may be a plurality of time segment data deletions corresponding to the target domain input data set, i.e., ET RT × (T-N) may not be satisfiedTB)。
The invention further improves that in the step 5), G Net in the GAN network adopts a two-norm loss function L2Loss of class, two-class loss function LbceGLoss function of loss and maximum mean difference LmmdThe weighted average of _ loss is used as the total loss function, and the specific loss function formula is as follows:
two-norm loss function L2The formula of _ loss is:
binary loss function LbceGThe formula of _ loss is:
maximum mean difference loss function LmmdThe formula of _ loss is:
whereinRepresenting the Regenerated Kernel Hilbert Space (RKHS), phi (-) represents the mapping function of the original feature space to RKHS, further developed to obtain:
for ease of calculation, a simplified equation is introduced:
wherein delta is a certain value parameter;
the G Net total loss function is as follows:
LGNet_loss=L2_loss+WbceG×LbceG_loss+Wmmd×Lmmd_loss
wherein, WbceGIs a binary loss function LbceGWeight of _ loss, WmmdIs a maximum mean difference loss function LmmdWeight of _ loss, according to experience WbceG∈[0.04,0.2],Wmmd∈[0.04,0.2]。
The invention further improves that in the step 5), D Net in the GAN network adopts the two-classification loss of the actual output data of the source domainTwo-classification loss L of source domain theoretical output databceLbLoss of class two for the _lossand actual output data of the target domainThe sum is taken as a total loss function, and the specific loss function formula is as follows:
binary loss of source domain actual output dataThe formula is as follows:
two-classification loss L of source domain theoretical output databceLbThe formula of _ loss is:
the D Net total loss function is as follows:
the invention further improves the method that in the step 6), during the training of the GAN transfer learning network, the optimizer selection and learning rate setting of G Net and the optimizer selection and learning rate setting of D Net are as follows:
in the training process, an optimizer of G Net selects an Adam optimization algorithm, the initial learning rate is set to be 0.002, after 200 times of data cycle training, the learning rate is set to be 0.0002, and then the learning rate is reduced by 90% after each 100 times of data cycle training;
the optimizer of D Net selects an Adam optimization algorithm, the initial learning rate is set to be 0.002, after 200 times of data cycle training, the learning rate is set to be 0.0002, and then the learning rate is reduced by 90% after every 100 times of data cycle training.
The invention has the following beneficial technical effects:
the invention provides a flow field prediction method of an oscillating flapping wing energy acquisition system based on transfer learning, which is characterized in that unknown fluid domain information of the next time point is predicted by transfer learning of a generative antagonistic neural network (GAN) from known fluid domain information of a plurality of time points after time dispersion, including node position information of a fluid domain after space dispersion and flow field parameter information on a node, and the unknown fluid domain information obtained by prediction is used as new known fluid domain information, so that the continuous prediction of the fluid domain information along with time can be realized under different parameter values. The whole flow field can be restored through node position information contained in the fluid domain information obtained through prediction and flow field parameter information on the nodes, flow field prediction is achieved, and fluid dynamics parameters, performance parameters and the like of the oscillating flapping wing energy collecting system at the prediction time point can be calculated. The method can greatly reduce time cost and economic cost, can efficiently and accurately predict the flow field of the oscillating flapping wing energy acquisition system, further can obtain fluid dynamics parameters, performance parameters and the like of the oscillating flapping wing energy acquisition system, and has important significance for research, design and development of the oscillating flapping wing energy acquisition system.
The method has the advantages that on one hand, compared with the traditional commercial CFD software, the time cost can be greatly saved, for similar flow fields or similar flow conditions, the method can directly predict the corresponding flow field change conditions, and can further calculate the fluid dynamics parameters, the performance parameters and the like of the oscillating flapping wing energy collecting system under the corresponding conditions. On the other hand, compared with experimental research, the method can greatly save economic cost and does not need expensive experimental instruments. In addition, compared with other artificial intelligence technologies, the method has the advantages that direct conversion from the fluid domain information at the known time point to the fluid domain information at the next unknown time point can be realized, so that the data structure is clear, and the further processing and analysis of the data are very convenient.
Drawings
FIG. 1 is a flow chart of the present invention for flow field prediction for an oscillating flapping wing energy harvesting system based on transfer learning;
FIG. 2 is a schematic diagram of a physical model of the oscillating flapping wing energy harvesting system of the present invention;
FIG. 3 is a spatial discrete model and a partial enlarged view of a grid of the oscillating flapping wing energy harvesting system of the present invention;
FIG. 4 is an illustration of a GAN migration learning flow field prediction network of the oscillating flapping wing energy harvesting system of the present invention;
FIG. 5 is an illustration of the oscillating flapping wing energy harvesting system generation network (G Net) architecture of the present invention;
FIG. 6 is an illustration of the oscillating flapping wing energy harvesting system discrimination network (D Net) architecture of the present invention;
FIG. 7 is a schematic diagram of loss functions and influence ranges of G Net and D Net of the oscillating flapping wing energy collecting system of the present invention;
FIG. 8 is a graph comparing the theoretical results and the predicted results of the pressure field of the oscillating flapping wing energy acquisition system of the present invention; wherein FIG. 8(a) is theoretical fluency and FIG. 8(b) is predicted fluency;
FIG. 9 is a comparison graph of the theoretical result and the predicted result of the X-direction velocity field of the oscillating flapping wing energy collecting system of the present invention; wherein FIG. 9(a) is theoretical fluency and FIG. 9(b) is predicted fluency;
FIG. 10 is a comparison graph of the theoretical results and the predicted results of the velocity field in the Y direction of the oscillating flapping wing energy collecting system of the present invention; where FIG. 10(a) is theoretical fluency and FIG. 10(b) is predicted fluency.
Detailed Description
The invention is described in more detail below with reference to a specific example according to the summary of the invention. The following is an application of the present invention, but not limited to this, and the implementer may change the parameters according to the specific problem and the actual application.
The invention provides a method for predicting a flow field of an oscillating flapping wing energy acquisition system based on transfer learning. The method comprises the steps of acquiring node position information and flow field parameter information on nodes after a fluid domain is spatially dispersed, performing matrixing processing and standardization processing, distributing a source domain and a target domain according to a set proportion, further extracting input data sets and corresponding theoretical output data sets of the source domain and the target domain, then respectively establishing a generation network (G Net) and a discrimination network (D Net) of generative antagonistic neural network (GAN) transfer learning, training, realizing fluid domain information prediction under different parameter values, and taking predicted unknown fluid domain information as new known fluid domain information, namely realizing continuous prediction of the fluid domain information along with time under different parameter values; the whole flow field is restored through node position information contained in the fluid domain information obtained through prediction and flow field parameter information on the nodes, flow field prediction is achieved, fluid dynamics parameters and performance parameters of the oscillating flapping wing energy collecting system at the prediction time point are calculated, and therefore researchers can analyze the oscillating flapping wing energy collecting system conveniently.
According to the flow chart of the flow field prediction of the oscillation flapping wing energy acquisition system based on the transfer learning in the figure 1, the flow field prediction of different Re values is carried out on the oscillation flapping wing energy acquisition system of a specific wing section, the physical model schematic diagram of the system is shown in the figure 2, and the system comprises the wing section 1, an inner area 2 which does rigid motion along with the wing section, an inner area interface 3, an outer area 4 which adopts a moving grid technology, an inlet boundary 5, an outlet boundary 6 and a symmetrical plane 7. The spatially discrete model and the grid partial magnification of the system are shown in fig. 3.
Firstly, obtaining model data of an oscillating flapping wing energy acquisition system;
and respectively collecting known fluid domain information of a plurality of continuous time-dispersed time points of the model aiming at different Re values, wherein the known fluid domain information comprises node position information and node flow field parameter information of the fluid domain after space dispersion. The node position information comprises an x coordinate and a y coordinate of a node, and the flow field parameter information comprises a speed u of a fluid on the node in the x direction, a speed v of the fluid on the node in the y direction and a pressure value p. Carrying out space discretization on the model to obtain grid information Hj,qPerforming unsteady state numerical simulation calculation, calculating an application layer flow model, spatially adopting a second-order difference format, temporally adopting a first-order difference format, solving by using a dynamic grid technology, and obtaining node position information at each time point according to a solving resultAnd flow field parameter information on the nodeWhere J is 1, 2., J is the number of mesh units, Q is 1, 2., Q is the number of nodes per mesh unit, here a two-dimensional four-node structured mesh, Q is 1,2,3,4, which respectively represents the number of four nodes per mesh unit, R is 1, 2., R is the total number of different selected Re values, T is 1, 2., T is the number of time points after time dispersion, I is 1, 2., I is the total number of nodes, N is 1, 2., N is the number of coordinate dimensions, here a two-dimensional problem, N is 1,2, which respectively represents the x and y coordinates of the nodes, M is 1, 2., M is the number of flow parameters to be considered, and where the x-direction of the fluid on the nodes is considered The three parameters, i.e., the velocity v in the y direction and the pressure value p, m is 1,2, and 3, and represents the velocity u in the x direction, the velocity v in the y direction, and the pressure value p, respectively, of the fluid at the node.
Secondly, preprocessing the model data of the oscillating flapping wing energy collecting system;
regular matrixing is carried out on the model data of the oscillating flapping wing energy acquisition system, and then node position information is rearrangedAnd flow field parameter information on the nodeCorresponding elements in the node are transformed to obtain the position information of the nodes after matrixingAnd flow field parameter information on the nodeWhere a 1,2,., a, B1, 2,., B, C1, 2,., C, A, B, C are the number of nodes in three mutually orthogonal directions, respectively, which satisfies a × B × C ═ I, in the problem here, thenA is the number of circumferential nodes along the airfoil surface, B is the number of nodes along the airfoil radial, and here is a two-dimensional plane problem, then C is 1.
For the node position information after matrixingAnd flow field parameter information on the nodeCarrying out standardization processing to obtain node position information G after standardization processingr,t,a,b,c,nAnd flow field parameter information F on the noder,t,a,b,c,m。
For the node position information after matrixingThe standardization processing method comprises the following steps:
wherein, MaxGnDenotes the maximum value of the nth coordinate dimension of a node among all its data, N is 1,2nRespectively representing the maximum value of the x-coordinate and the y-coordinate of the node, MinG, respectively, in all its datanDenotes the minimum value of the nth coordinate dimension of the node among all its data, N is 1,2nRespectively, the x-coordinate and the y-coordinate of the node each represent the minimum of all of its data.
To the parameter information of the flow field on the nodes after matrixingThe standardization processing method comprises the following steps:
wherein, MaxFmThe M-th parameter of a node is represented as the maximum value of all data, M is 1,2,.. the M is the number of flow field parameters to be considered, three parameters of the speed u of the fluid on the node in the x direction, the speed v of the fluid on the node in the y direction and the pressure value p are considered, and M is 1,2,3, MaxFmThe maximum value, MinF, of each of the three parameters, i.e., the velocity u in the x-direction, the velocity v in the y-direction, and the pressure value p, of the fluid at the node pointmThe M-th parameter of a node is represented as the minimum value of all data, M is 1,2,.. the M is the number of flow field parameters to be considered, three parameters of the speed u of the fluid on the node in the x direction, the speed v of the fluid on the node in the y direction and the pressure value p are considered, and M is 1,2,3, MinFmThe three parameters, i.e. the velocity u in the x-direction of the fluid at the node, the velocity v in the y-direction and the pressure value p, respectively, represent the minimum of all their data.
Thirdly, splitting a source domain and a target domain of the model data of the oscillating flapping wing energy collecting system;
the normalized node position information G corresponding to different values selected by the parametersr,t,a,b,c,nAnd flow field parameter information F on the noder,t,a,b,c,mSplitting at random or according to requirement according to a certain proportion, using one part as source domain and source domain node position informationAnd flow field parameter information on the nodeAnother part is used as target domain, target domain node position informationAnd flow field parameter information on the nodeWherein RS is the total number of different values selected as the parameters of the source domain, RT is 1,2, and RT is the total number of different values selected as the parameters of the target domain, and RS + RT is satisfied. Here, a 1:1 relationship is selected to select the source domain and the target domain, and RS ═ RT ═ R/2.
Fourthly, respectively extracting respective input data sets and theoretical output data sets in the source domain and the target domain;
the number N of continuous time points required according to actual input dataTBFor example, where N is selectedTBAnd 2, further extracting data of the source domain and the target domain to obtain respective input data sets and theoretical output data sets.
For the source domain, the node position information of the source domain is firstlyAnd flow field parameter information on the nodeSplicing to obtain the total data information of the source domainWhere K is 1,2, the., K is N + M, where N is 2, where three parameters of the velocity u in the x direction, the velocity v in the y direction, and the pressure value p of the fluid on the node are considered, and M is 3, where K is 5, K is 1,2,3,4,5, which respectively represent the x coordinate and the y coordinate of the node and the velocity u in the x direction, the velocity v in the y direction, and the pressure value p of the fluid on the node.
The source domain input dataset is:
wherein each source domain input data isES, which is the total amount of source domain data and satisfies ES RS × (T-N)TB)。
The corresponding source domain theoretical output data set is:
wherein each source domain theoretical output data isES, which is the total amount of source domain data and satisfies ES RS × (T-N)TB)。
For the target domain, firstly, the node position information of the target domainAnd flow field parameter information on the nodeSplicing to obtain the total data information of the target domainWhere K is 1,2, the., K is N + M, where N is 2, where three parameters of the velocity u in the x direction, the velocity v in the y direction, and the pressure value p of the fluid on the node are considered, and M is 3, where K is 5, K is 1,2,3,4,5, which respectively represent the x coordinate and the y coordinate of the node and the velocity u in the x direction, the velocity v in the y direction, and the pressure value p of the fluid on the node.
The target domain input dataset is:
wherein each target domain inputs data as ET 1, 2.. ET, ET is the total amount of target domain data, and there may be a plurality of time period data missing, that is, ET RT × (T-N) may not be satisfiedTB)。
The corresponding target domain theoretical output data set is:
wherein each target domain theoretical output data is ET 1, 2.. ET, ET is the total amount of target domain data, and there may be a plurality of time segment data deletions corresponding to the target domain input data set, i.e., ET RT × (T-N) may not be satisfiedTB)。
Fifthly, constructing a fluid domain information prediction GAN migration learning network of the oscillating flapping wing energy acquisition system;
the fluid domain information prediction GAN transfer learning network of the oscillating flapping wing energy acquisition system comprises a generation network (GNet) and a discrimination network (D Net), and the specific GAN transfer learning network structure is shown in FIG. 4.
Firstly, a convolution residual error network is adopted to establish a characteristic extraction part of G Net, and Input data Input from a source domain is realizedSDOr target field Input data InputTDExtracting corresponding source domain fluid domain Feature information FeatureSDOr target domain fluid domain Feature information FeatureTDThen, a deconvolution residual error network is adopted to establish a predicted fluid domain information part of G Net, and the Feature information Feature of the source domain fluid domain is realizedSDOr target domain fluid domain Feature information FeatureTDGenerating corresponding source domain actual Output data OutputSDOr target domain actual Output data OutputTD. Wherein each source domain fluid domain characteristic informationES, each source domain actually outputs data asES, characteristic information of each target domain fluid domainET 1,2, ET, each target field actually outputting data asThe ET 1,2, the ET and the G Net have 21 layers of networks, and a specific G Net framework is shown in fig. 5. Secondly, a convolution residual error network is adopted to establish a DNet network to judge the truth of the fluid domain information, and the actual output data of each source domain is realizedAnd each source domain theoretical output data corresponding to the source domain theoretical output dataTo obtain the corresponding discrimination resultIs 0 and1, the data sets of the discrimination result are respectively expressed as D _ OutputSDAnd D _ LableSDSimultaneously undertake the actual output of data from each target domainObtaining a discrimination resultThe corresponding data set of the discrimination result is D _ OutputTDA two-classification loss function part for calculating network training of G Net,the D Net has 19 layers of network, and the specific D Net framework is shown in FIG. 6.
With the quantity SSDA batch of source domain input data and quantity STDIs used to calculate the loss of the maximum mean difference, where S is takenSD=16,STD=16。
G Net in GAN network adopts two-norm loss function L2Loss of class, two-class loss function LbceGLoss function of loss and maximum mean difference LmmdA weighted average of loss as a function of total loss.
Two-norm loss function L2The formula of _ loss is:
binary loss function LbceGThe formula of _ loss is:
maximum mean difference loss function LmmdThe formula of _ loss is:
whereinRepresenting the Regenerated Kernel Hilbert Space (RKHS), phi (-) represents the mapping function of the original feature space to RKHS, further developed to obtain:
for ease of calculation, a simplified equation is introduced:
where δ is a parameter of value.
The G Net total loss function is as follows:
LGNet_loss=L2_loss+WbceG×LbceG_loss+Wmmd×Lmmd_loss
wherein, WbceGIs a binary loss function LbceGWeight of _ loss, WmmdIs a maximum mean difference loss function LmmdWeight of _ loss, according to experience WbceG∈[0.04,0.2],Wmmd∈[0.04,0.2]At which W is takenbceG=0.1,WmmdThe loss functions and their range of influence are shown in fig. 7, which is 0.1.
D Net in GAN network adopts two-classification loss of source domain actual output dataTwo-classification loss L of source domain theoretical output databceLbLoss of class two for the _lossand actual output data of the target domainThe sum as a function of total loss.
two-classification loss L of source domain theoretical output databceLbThe formula of _ loss is:
the D Net total loss function is as follows:
the loss functions and their range of influence are shown in figure 7.
Sixthly, training a fluid domain information prediction GAN migration learning network of the oscillating flapping wing energy acquisition system;
firstly, the number of the oscillating flapping wing energy acquisition systems is SSDThe fluid domain information of a batch of source domain input data is predicted through G Net in a GAN transfer learning network, and the obtained quantity is SSDA batch of source domain fluid domain characteristic information and a batch of source domain actual output data, and then the number of the oscillating flapping wing energy acquisition systems is STDThe fluid domain information prediction is carried out on a batch of target domain input data for calculating the maximum mean difference loss through G Net in a GAN migration learning network, and the obtained quantities are STDA batch of target domain fluid domain characteristic information and a batch of target domain actual output data.
Then the obtained quantity is SSDA batch of source domain actual output data, the corresponding quantity is SSDA batch of source domain theoretical output data and the quantity is STDAnd (3) distinguishing the batch of target domain actual output data by DNet in the GAN migration learning network, and updating parameters of D Net once according to the obtained distinguishing result, so that the D Net can better distinguish the fluid domain information prediction result predicted by G Net and the fluid domain information theoretical result.
Then the obtained quantity is SSDA batch of source domain actual output data and the number STDThe actual output data of a batch of target domains are respectively judged through D Net which is subjected to parameter updating once in the GAN migration learning network, and according to the obtained judgment result and the obtained quantity are SSDA batch of source domain fluid domain characteristic information and a batch of source domain actual output data, the quantity is STDIs known to correspond to the actual output of the source domain, and has a quantity SSDThe G Net is subjected to parameter updating once by a batch of source domain theoretical output data, so that the fluid domain information prediction result of the G Net is closer to the fluid domain information theoretical result, and the judgment of the fluid domain information prediction result and the fluid domain information theoretical result of the D Net is confused as much as possible.
Through the continuous and repeated confrontation of G Net and D Net, the judgment of the fluid domain information prediction result and the fluid domain information theoretical result by the D Net is more accurate, and the fluid domain information prediction result of the G Net is forced to be closer to the fluid domain information theoretical result.
In the training process, the optimizer of G Net selects an Adam optimization algorithm, the initial learning rate is set to be 0.002, after 200 times of data loop training, the learning rate is set to be 0.0002, and then the learning rate is reduced by 90% after every 100 times of data loop training. The optimizer of DNet selects Adam optimization algorithm, initial learning rate is set to be 0.002, after 200 times of training, learning rate is set to be 0.0002, and then learning rate is reduced by 90% after every 100 times of training.
Seventhly, processing and analyzing flow field prediction results of the oscillating flapping wing energy acquisition system
The trained fluid domain information prediction GAN migration learning network of the oscillating flapping wing energy acquisition system is used for predicting fluid domain information under other Re, and the predicted unknown fluid domain information at the next time point is used as new known fluid domain information, so that the fluid domain information prediction at the next time point can be continuously carried out, and the continuous prediction of the fluid domain information along with time is realized under different parameter values. The whole flow field can be restored through node position information and flow field parameter information on the nodes contained in the fluid domain information obtained through prediction, the flow field prediction results are shown in fig. 8-10, theoretical flow field conditions and predicted flow field results of a pressure field, an X-direction speed field and a Y-direction speed field at a certain time point are respectively given, the predicted flow field results are well matched with the theoretical flow field conditions, and errors are within a reasonable range. Furthermore, fluid dynamics parameters and performance parameters such as torque, lift coefficient and resistance coefficient of the oscillating flapping wing energy collecting system at a prediction time point can be calculated through a flow field prediction result, so that researchers can analyze the oscillating flapping wing energy collecting system.
Claims (8)
1. A method for predicting a flow field of an oscillating flapping wing energy acquisition system based on transfer learning is characterized in that unknown fluid domain information of a next time point is predicted from known fluid domain information of a plurality of continuous time points after time dispersion; specifically, by acquiring node position information and node flow field parameter information of a fluid domain after spatial dispersion, performing matrixing processing and standardization processing, distributing a source domain and a target domain according to a set proportion, further extracting input data sets and corresponding theoretical output data sets of the source domain and the target domain, then respectively establishing a generation network and a discrimination network of generation type antagonistic neural network transfer learning, after training, realizing fluid domain information prediction under different parameter values, and taking unknown fluid domain information obtained by prediction as new known fluid domain information, namely realizing continuous prediction of the fluid domain information along with time under different parameter values; the whole flow field is restored through node position information contained in the fluid domain information obtained through prediction and flow field parameter information on the nodes, flow field prediction is achieved, fluid dynamics parameters and performance parameters of the oscillating flapping wing energy collecting system at the prediction time point are calculated, and therefore researchers can analyze the oscillating flapping wing energy collecting system conveniently.
2. The method for predicting the flow field of the oscillating flapping wing energy collecting system based on the transfer learning of claim 1, comprising the following steps:
1) obtaining model data of oscillating flapping wing energy collecting system
Respectively collecting known fluid domain information of a plurality of continuous time-dispersed time points of the model aiming at different parameter values, wherein the known fluid domain information comprises node position information and node flow field parameter information of the fluid domain after space dispersion;carrying out space discretization on the model to obtain grid information Hj,qPerforming unsteady state numerical simulation calculation, calculating an application layer flow model, spatially adopting a second-order difference format, temporally adopting a first-order difference format, solving by using a dynamic grid technology, and obtaining node position information at each time point according to a solving resultAnd flow field parameter information on the nodeJ is 1, 2.. the J, J is the number of grid units, Q is 1, 2.. the Q, Q is the number of nodes of each grid unit, R is 1, 2.. the R, R is the total number of different selected parameter values, T is 1, 2.. the T, T is the number of time points after time dispersion, I is 1, 2.. the I, I is the total number of nodes, N is 1, 2.. the N, N is the number of coordinate dimensions, M is 1, 2.. the M, M is the number of flow field parameters to be considered;
2) model data preprocessing of oscillating flapping wing energy collecting system
Regular matrixing is carried out on the model data of the oscillating flapping wing energy acquisition system, and then node position information is rearrangedAnd flow field parameter information on the nodeCorresponding elements in the node are transformed to obtain the position information of the nodes after matrixingAnd flow field parameter information on the nodeWherein a 1,2, a, B1, 2, B, C1, 2, C, A, B, C are the number of nodes in three mutually orthogonal directions, respectively, and a × B × C ═ I;
For the node position information after matrixingAnd flow field parameter information on the nodeCarrying out standardization processing to obtain node position information G after standardization processingr,t,a,b,c,nAnd flow field parameter information F on the noder,t,a,b,c,m;
3) Model data splitting source domain and target domain of oscillating flapping wing energy acquisition system
The normalized node position information G corresponding to different values selected by the parametersr,t,a,b,c,nAnd flow field parameter information F on the noder,t,a,b,c,mSplitting at random or according to the requirement according to the set proportion, using one part as the source domain and the node position information of the source domainAnd flow field parameter information on the nodeAnother part is used as target domain, target domain node position informationAnd flow field parameter information on the nodeRS, RS are different total numbers of values selected as parameters of a source domain, RT is 1,2,., RT are different total numbers of values selected as parameters of a target domain, and RS + RT is satisfied;
4) extracting respective input data set and theoretical output data set in source domain and target domain
The number N of continuous time points required according to actual input dataTBTo the sourceFurther data extraction is carried out on the domain and the target domain, and respective input data sets and theoretical output data sets are obtained;
for the source domain, the node position information of the source domain is firstlyAnd flow field parameter information on the nodeSplicing to obtain the total data information of the source domainWhere K1, 2,., K N + M, a source domain Input dataset Input is obtainedSDAnd a corresponding source domain theoretical output data set LableSD;
For the target domain, firstly, the node position information of the target domainAnd flow field parameter information on the nodeSplicing to obtain the total data information of the target domainWhere K is 1,2,., K is N + M, the target domain Input dataset Input is obtainedTDAnd a corresponding target domain theoretical output data set LableTD;
5) Fluid domain information prediction GAN migration learning network for constructing oscillating flapping wing energy acquisition system
The fluid domain information prediction GAN migration learning network of the oscillating flapping wing energy acquisition system comprises a generation network and a judgment network, wherein the generation network is called G Net for short, and the judgment network is called D Net for short;
firstly, a convolution residual error network is adopted to establish a characteristic extraction part of G Net, and Input data Input from a source domain is realizedSDOr target domainInput data InputTDExtracting corresponding source domain fluid domain Feature information FeatureSDOr target domain fluid domain Feature information FeatureTDThen, a deconvolution residual error network is adopted to establish a predicted fluid domain information part of G Net, and the Feature information Feature of the source domain fluid domain is realizedSDOr target domain fluid domain Feature information FeatureTDGenerating corresponding source domain actual Output data OutputSDOr target domain actual Output data OutputTD(ii) a Wherein each source domain fluid domain characteristic informationES, each source domain actually outputs data asES, characteristic information of each target domain fluid domainET 1,2, ET, each target field actually outputting data asET 1, 2.. ET, G Net total 21 layers of network; secondly, a D Net network is established by adopting a convolution residual error network to judge whether the fluid domain information is true or false, and the actual output data of each source domain is realizedAnd each source domain theoretical output data corresponding to the source domain theoretical output dataTo obtain the corresponding discrimination resultIs 0 and1, the data sets of the discrimination result are respectively expressed as D _ OutputSDAnd D _ LableSDSimultaneously undertake the actual output of data from each target domainObtaining a discrimination resultThe corresponding data set of the discrimination result is D _ OutputTDThe two-classification loss function part is used for calculating the network training of G Net, and D Net has 19 layers of networks in total;
with the quantity SSDA batch of source domain input data and quantity STDA batch of target domain input data for calculating maximum mean difference loss;
g Net in GAN network adopts two-norm loss function L2Loss of class, two-class loss function LbceGLoss function of loss and maximum mean difference LmmdA weighted average of _ loss as a total loss function;
d Net in GAN network adopts two-classification loss of source domain actual output dataTwo-classification loss L of source domain theoretical output databceLbLoss of class two for the _lossand actual output data of the target domainThe sum as a function of total loss;
6) fluid domain information prediction GAN migration learning network for training oscillating flapping wing energy acquisition system
Firstly, the number of the oscillating flapping wing energy acquisition systems is SSDThe fluid domain information of a batch of source domain input data is predicted through G Net in a GAN transfer learning network, and the obtained quantity is SSDA batch of source domain fluid domain characteristic information and a batch of source domain actual output data, and then the number of the oscillating flapping wing energy acquisition systems is STDThe fluid domain information prediction is carried out on a batch of target domain input data for calculating the maximum mean difference loss through G Net in a GAN migration learning network, and the obtained quantities are STDA batch of target domain fluid domain characteristic information and a batch of target domain actual output data;
then the obtained quantity is SSDA batch of source domain actual output data, the corresponding quantity is SSDA batch of source domain theoretical output data and the quantity is STDThe batch of actual output data of the target domain are respectively judged through D Net in the GAN migration learning network, and the D Net is subjected to parameter updating once according to the obtained judgment result, so that the D Net can better distinguish a fluid domain information prediction result predicted by GNet and a fluid domain information theoretical result;
then the obtained quantity is SSDA batch of source domain actual output data and the number STDThe actual output data of a batch of target domains are respectively judged through D Net which is subjected to parameter updating once in the GAN migration learning network, and according to the obtained judgment result and the obtained quantity are SSDA batch of source domain fluid domain characteristic information and a batch of source domain actual output data, the quantity is STDIs known to correspond to the actual output of the source domain, and has a quantity SSDThe G Net is subjected to parameter updating once by a batch of source domain theoretical output data, so that the fluid domain information prediction result of the G Net is closer to the fluid domain information theoretical result, and the judgment of the fluid domain information prediction result and the fluid domain information theoretical result of the D Net is confused as much as possible;
through the continuous and repeated confrontation of G Net and D Net, the judgment of the fluid domain information prediction result and the fluid domain information theoretical result by the D Net is more accurate, and the fluid domain information prediction result of the G Net is forced to be closer to the fluid domain information theoretical result;
7) flow field prediction result processing and analysis of oscillating flapping wing energy acquisition system
The trained fluid domain information prediction GAN migration learning network of the oscillating flapping wing energy acquisition system is used for predicting fluid domain information under other parameter values, and the predicted unknown fluid domain information at the next time point is used as new known fluid domain information, so that the fluid domain information prediction at the next time point can be continuously carried out, and the continuous prediction of the fluid domain information along with time is realized under different parameter values; the whole flow field can be restored through node position information and flow field parameter information on the nodes contained in the fluid domain information obtained through prediction; furthermore, fluid dynamics parameters and performance parameters such as torque, lift coefficient and resistance coefficient of the oscillating flapping wing energy collecting system at a prediction time point are calculated according to a flow field prediction result so as to provide convenience for researchers to analyze the oscillating flapping wing energy collecting system.
3. The method for predicting the flow field of the oscillating flapping wing energy collecting system based on the transfer learning of claim 2, wherein in the step 2), the position information of the nodes after the matrixing is carried outAnd flow field parameter information on the nodeThe normalization process was performed as follows:
for the node position information after matrixingThe standardization processing method comprises the following steps:
wherein, MaxGnRepresents the maximum value of the nth coordinate dimension of the node in all data, wherein N is 1,2nRepresenting the minimum value of the nth coordinate dimension of the node in all data of the nth coordinate dimension;
to the parameter information of the flow field on the nodes after matrixingThe standardization processing method comprises the following steps:
wherein, MaxFmRepresents the maximum value of the mth parameter of the node in all data, wherein M is 1,2mRepresents the minimum value of the mth parameter of the node in all data thereof.
4. The method for predicting the flow field of the oscillatory flapping wing energy collecting system based on the transfer learning of claim 3, wherein in the step 4), the source domain Input data set Input is obtainedSDAnd a corresponding source domain theoretical output data set LableSDThe method comprises the following steps:
the source domain input dataset is:
wherein each source domain input data isES, which is the total amount of source domain data and satisfies
ES=RS×(T-NTB)
The corresponding source domain theoretical output data set is:
5. The method for predicting the flow field of the oscillatory flapping wing energy collecting system based on the transfer learning of claim 4, wherein in the step 4), the target domain Input data set Input is obtainedTDAnd a corresponding target domain theoretical output data set LableTDThe method comprises the following steps:
the target domain input dataset is:
wherein each target domain inputs data asET 1,2, ET is the total amount of target domain data, and there are multiple time period data missing, i.e. ET RT × (T-N) may not be satisfiedTB);
The corresponding target domain theoretical output data set is:
6. The method for predicting the flow field of the oscillatory flapping wing energy collecting system based on the transfer learning of claim 5, wherein in the step 5), G Net in a GAN network adopts a two-norm loss function L2Loss of class, two-class loss function LbceGLoss function of loss and maximum mean difference LmmdThe weighted average of _ loss is used as the total loss function, and the specific loss function formula is as follows:
two-norm loss function L2The formula of _ loss is:
binary loss function LbceGThe formula of _ loss is:
maximum mean difference loss function LmmdThe formula of _ loss is:
whereinRepresenting the Regenerated Kernel Hilbert Space (RKHS), phi (-) represents the mapping function of the original feature space to RKHS, further developed to obtain:
for ease of calculation, a simplified equation is introduced:
the G Net total loss function is as follows:
LGNet_loss=L2_loss+WbceG×LbceG_loss+Wmmd×Lmmd_loss
wherein, WbceGIs a binary loss function LbceGWeight of _ loss, WmmdIs a maximum mean difference loss function LmmdWeight of _ loss, according to experience WbceG∈[0.04,0.2],Wmmd∈[0.04,0.2]。
7. The method for predicting the flow field of the oscillatory flapping wing energy collecting system based on the transfer learning of claim 6, wherein in the step 5), D Net in the GAN network adopts the binary classification loss of the actual output data of the source domainTwo-classification loss L of source domain theoretical output databceLbLoss of class two for the _lossand actual output data of the target domainThe sum is taken as a total loss function, and the specific loss function formula is as follows:
two-classification loss L of source domain theoretical output databceLbThe formula of _ loss is:
the D Net total loss function is as follows:
8. the method for predicting the flow field of the oscillatory flapping wing energy collecting system based on the transfer learning of claim 7, wherein in the step 6), in the process of training the GAN transfer learning network, the method for selecting the optimizer and setting the learning rate of G Net and the method for selecting the optimizer and setting the learning rate of D Net are as follows:
in the training process, an optimizer of G Net selects an Adam optimization algorithm, the initial learning rate is set to be 0.002, after 200 times of data cycle training, the learning rate is set to be 0.0002, and then the learning rate is reduced by 90% after each 100 times of data cycle training;
the optimizer of D Net selects an Adam optimization algorithm, the initial learning rate is set to be 0.002, after 200 times of data cycle training, the learning rate is set to be 0.0002, and then the learning rate is reduced by 90% after every 100 times of data cycle training.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910838392.1A CN110705029B (en) | 2019-09-05 | 2019-09-05 | Flow field prediction method of oscillating flapping wing energy acquisition system based on transfer learning |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910838392.1A CN110705029B (en) | 2019-09-05 | 2019-09-05 | Flow field prediction method of oscillating flapping wing energy acquisition system based on transfer learning |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110705029A true CN110705029A (en) | 2020-01-17 |
CN110705029B CN110705029B (en) | 2021-09-07 |
Family
ID=69194721
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910838392.1A Active CN110705029B (en) | 2019-09-05 | 2019-09-05 | Flow field prediction method of oscillating flapping wing energy acquisition system based on transfer learning |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110705029B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111444614A (en) * | 2020-03-26 | 2020-07-24 | 西安交通大学 | Flow field reconstruction method based on graph convolution |
CN111581784A (en) * | 2020-04-17 | 2020-08-25 | 浙江大学 | Flapping wing motion parameter optimization method based on data-driven self-adaptive quasi-steady-state model |
CN111724487A (en) * | 2020-06-19 | 2020-09-29 | 广东浪潮大数据研究有限公司 | Flow field data visualization method, device, equipment and storage medium |
CN111859746A (en) * | 2020-07-10 | 2020-10-30 | 西安交通大学 | Method for predicting variable working condition performance of turbomachinery based on flow field reconstruction |
CN112784508A (en) * | 2021-02-12 | 2021-05-11 | 西北工业大学 | Deep learning-based airfoil flow field rapid prediction method |
CN114996658A (en) * | 2022-07-20 | 2022-09-02 | 中国空气动力研究与发展中心计算空气动力研究所 | Projection-based hypersonic aircraft aerodynamic heat prediction method |
CN115758911A (en) * | 2022-12-07 | 2023-03-07 | 中国石油大学(华东) | Fusion point cloud residual error network and flow field and pressure field prediction method considering slippage |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102135760A (en) * | 2010-12-16 | 2011-07-27 | 天津工业大学 | Neural network energy coordinated controller for microgrid |
CN106251401A (en) * | 2016-07-06 | 2016-12-21 | 西安交通大学 | A kind of generation system and method for hexahedron 20 node unit grid |
CN107039975A (en) * | 2017-05-27 | 2017-08-11 | 上海电气分布式能源科技有限公司 | A kind of distributed energy resource system energy management method |
CN108563906A (en) * | 2018-05-02 | 2018-09-21 | 北京航空航天大学 | A kind of short fiber reinforced composite macro property prediction technique based on deep learning |
CN108932554A (en) * | 2017-05-26 | 2018-12-04 | 西安交通大学 | The method for optimizing configuration and device of a kind of wind power plant flow field measuring point |
CN109060001A (en) * | 2018-05-29 | 2018-12-21 | 浙江工业大学 | A kind of multiple operating modes process soft-measuring modeling method based on feature transfer learning |
CN109102126A (en) * | 2018-08-30 | 2018-12-28 | 燕山大学 | One kind being based on depth migration learning theory line loss per unit prediction model |
CN109918752A (en) * | 2019-02-26 | 2019-06-21 | 华南理工大学 | Mechanical failure diagnostic method, equipment and medium based on migration convolutional neural networks |
US20190243933A1 (en) * | 2018-02-07 | 2019-08-08 | Incucomm, Inc. | System and method that characterizes an object employing virtual representations thereof |
-
2019
- 2019-09-05 CN CN201910838392.1A patent/CN110705029B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102135760A (en) * | 2010-12-16 | 2011-07-27 | 天津工业大学 | Neural network energy coordinated controller for microgrid |
CN106251401A (en) * | 2016-07-06 | 2016-12-21 | 西安交通大学 | A kind of generation system and method for hexahedron 20 node unit grid |
CN108932554A (en) * | 2017-05-26 | 2018-12-04 | 西安交通大学 | The method for optimizing configuration and device of a kind of wind power plant flow field measuring point |
CN107039975A (en) * | 2017-05-27 | 2017-08-11 | 上海电气分布式能源科技有限公司 | A kind of distributed energy resource system energy management method |
US20190243933A1 (en) * | 2018-02-07 | 2019-08-08 | Incucomm, Inc. | System and method that characterizes an object employing virtual representations thereof |
CN108563906A (en) * | 2018-05-02 | 2018-09-21 | 北京航空航天大学 | A kind of short fiber reinforced composite macro property prediction technique based on deep learning |
CN109060001A (en) * | 2018-05-29 | 2018-12-21 | 浙江工业大学 | A kind of multiple operating modes process soft-measuring modeling method based on feature transfer learning |
CN109102126A (en) * | 2018-08-30 | 2018-12-28 | 燕山大学 | One kind being based on depth migration learning theory line loss per unit prediction model |
CN109918752A (en) * | 2019-02-26 | 2019-06-21 | 华南理工大学 | Mechanical failure diagnostic method, equipment and medium based on migration convolutional neural networks |
Non-Patent Citations (3)
Title |
---|
A.S. KEÇELI;A. KAYA: "Violent activity detection with transfer learning method", 《ELECTRONICS LETTERS 》 * |
熊晔颖: "基于深度学习的稀疏流场处理方法的研究与实现", 《中国优秀硕士学位论文全文数据库(基础科学辑)》 * |
谢永慧 等: "振荡扑翼流场能量采集研究进展", 《中国电机工程学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111444614A (en) * | 2020-03-26 | 2020-07-24 | 西安交通大学 | Flow field reconstruction method based on graph convolution |
CN111581784A (en) * | 2020-04-17 | 2020-08-25 | 浙江大学 | Flapping wing motion parameter optimization method based on data-driven self-adaptive quasi-steady-state model |
CN111581784B (en) * | 2020-04-17 | 2021-12-21 | 浙江大学 | Flapping wing motion parameter optimization method based on data-driven self-adaptive quasi-steady-state model |
CN111724487A (en) * | 2020-06-19 | 2020-09-29 | 广东浪潮大数据研究有限公司 | Flow field data visualization method, device, equipment and storage medium |
WO2021253666A1 (en) * | 2020-06-19 | 2021-12-23 | 广东浪潮智慧计算技术有限公司 | Flow field data visualization method, apparatus and device, and storage medium |
CN111724487B (en) * | 2020-06-19 | 2023-05-16 | 广东浪潮大数据研究有限公司 | Flow field data visualization method, device, equipment and storage medium |
CN111859746A (en) * | 2020-07-10 | 2020-10-30 | 西安交通大学 | Method for predicting variable working condition performance of turbomachinery based on flow field reconstruction |
CN112784508A (en) * | 2021-02-12 | 2021-05-11 | 西北工业大学 | Deep learning-based airfoil flow field rapid prediction method |
CN114996658A (en) * | 2022-07-20 | 2022-09-02 | 中国空气动力研究与发展中心计算空气动力研究所 | Projection-based hypersonic aircraft aerodynamic heat prediction method |
CN115758911A (en) * | 2022-12-07 | 2023-03-07 | 中国石油大学(华东) | Fusion point cloud residual error network and flow field and pressure field prediction method considering slippage |
Also Published As
Publication number | Publication date |
---|---|
CN110705029B (en) | 2021-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110705029B (en) | Flow field prediction method of oscillating flapping wing energy acquisition system based on transfer learning | |
Wang et al. | Hybrid physics-based and data-driven models for smart manufacturing: Modelling, simulation, and explainability | |
CN110288048B (en) | Submarine pipeline risk assessment method of SVM directed acyclic graph | |
CN104408518B (en) | Based on the neural network learning optimization method of particle swarm optimization algorithm | |
CN112541572B (en) | Residual oil distribution prediction method based on convolutional encoder-decoder network | |
CN112784508A (en) | Deep learning-based airfoil flow field rapid prediction method | |
CN110348059B (en) | Channel internal flow field reconstruction method based on structured grid | |
CN112084727A (en) | Transition prediction method based on neural network | |
CN103455612B (en) | Based on two-stage policy non-overlapped with overlapping network community detection method | |
CN113283186B (en) | Universal grid self-adaption method for CFD | |
CN112016253A (en) | High-fidelity chaotic polynomial correction method suitable for CFD uncertainty quantification | |
CN113221475A (en) | Grid self-adaption method for high-precision flow field analysis | |
CN116629352A (en) | Hundred million-level parameter optimizing platform | |
CN117854643B (en) | MEMS membrane simulation method and system based on graphic neural network | |
Vadrot et al. | Survey of machine-learning wall models for large-eddy simulation | |
Durasov et al. | Debosh: Deep bayesian shape optimization | |
CN116894180B (en) | Product manufacturing quality prediction method based on different composition attention network | |
CN114240687A (en) | Energy hosting efficiency analysis method suitable for comprehensive energy system | |
CN113609763A (en) | Uncertainty-based satellite component layout temperature field prediction method | |
CN109116300A (en) | A kind of limit learning position method based on non-abundant finger print information | |
Choubineh et al. | An innovative application of deep learning in multiscale modeling of subsurface fluid flow: Reconstructing the basis functions of the mixed GMsFEM | |
CN116415482A (en) | Wing flow field analysis method based on graph neural network MeshGraphNet | |
CN115438575A (en) | Analysis method for high-precision airfoil flow field prediction | |
CN114723877A (en) | Countermeasure sample generation method for three-dimensional sparse convolution network | |
CN110889396A (en) | Energy internet disturbance classification method and device, electronic equipment and storage medium |
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 |