[go: up one dir, main page]

CN118211517B - Jet flow pneumatic/infrared radiation coupling gradient calculation method based on accompanying method - Google Patents

Jet flow pneumatic/infrared radiation coupling gradient calculation method based on accompanying method Download PDF

Info

Publication number
CN118211517B
CN118211517B CN202410614660.2A CN202410614660A CN118211517B CN 118211517 B CN118211517 B CN 118211517B CN 202410614660 A CN202410614660 A CN 202410614660A CN 118211517 B CN118211517 B CN 118211517B
Authority
CN
China
Prior art keywords
grid
grids
infrared radiation
equation
accompanying
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
CN202410614660.2A
Other languages
Chinese (zh)
Other versions
CN118211517A (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.)
Institute of Aerospace Technology of China Aerodynamics Research and Development Center
Original Assignee
Institute of Aerospace Technology of China Aerodynamics Research and Development Center
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 Institute of Aerospace Technology of China Aerodynamics Research and Development Center filed Critical Institute of Aerospace Technology of China Aerodynamics Research and Development Center
Priority to CN202410614660.2A priority Critical patent/CN118211517B/en
Publication of CN118211517A publication Critical patent/CN118211517A/en
Application granted granted Critical
Publication of CN118211517B publication Critical patent/CN118211517B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)

Abstract

The invention discloses a jet flow pneumatic/infrared radiation coupling gradient calculation method based on an accompanying method, which comprises the following steps: based on an integral form jet radiation transfer equation, a Lagrangian multiplier is introduced to obtain an accompanying equation of the infrared radiation problem; and solving the adjoint equation of the infrared radiation problem to obtain adjoint variables, and obtaining the gradient based on the adjoint method according to the adjoint variables. The invention can realize the rapid solving of the infrared radiation gradient.

Description

Jet flow pneumatic/infrared radiation coupling gradient calculation method based on accompanying method
Technical Field
The invention relates to the technical field of jet infrared radiation optimization design of aircrafts, in particular to a jet pneumatic/infrared radiation coupling gradient calculation method based on an accompanying method.
Background
The tail jet is one of the important radiation sources of the infrared radiation of the exhaust system of the aircraft, the size of the tail jet can be several times of the size of the aircraft, and the infrared radiation can be detected in a wide angle range. The spray pipe appearance design can increase the mixing between the high-temperature jet flow and the environment atmosphere by increasing the contact area between the outlet jet flow and the external environment, and reduce the temperature of the tail jet flow, thereby achieving the purpose of inhibiting the infrared radiation illuminance. The data show that the infrared radiation illuminance of the tail jet can be reduced by 60% by adopting the blending technology. At present, the research on the spray tube appearance optimization design method aiming at the infrared radiation of the tail jet is less, the published research at present mainly adopts a trial-and-error method, the influence of the geometric parameters on the infrared radiation illuminance is discussed through manual change of the geometric parameters, and the automation of the design process is difficult to realize, so that the development of the design method suitable for the infrared radiation optimization of the tail jet is needed.
The accompanying method is an important method for precisely designing the high-dimensional design problem based on the gradient by deducing the accompanying equation of the control equation, solving the control equation once and solving the accompanying equation once to accurately and efficiently acquire the gradient of the objective function. How to apply the accompanying method to the spray pipe appearance design to realize the jet flow infrared radiation illuminance inhibition is the key for improving the pneumatic/infrared meter efficiency of the aircraft spray pipe.
Disclosure of Invention
In view of the above, the invention provides a jet flow pneumatic/infrared radiation coupling gradient calculation method based on an accompanying method, which solves the pneumatic/infrared radiation coupling accompanying equation to realize efficient and high-precision solving of the infrared radiation gradient.
The invention discloses a jet flow pneumatic/infrared radiation coupling gradient calculation method based on an accompanying method, which comprises the following steps:
Based on an integral form jet radiation transfer equation, a Lagrangian multiplier is introduced to obtain an accompanying equation of the infrared radiation problem; and solving the adjoint equation of the infrared radiation problem to obtain adjoint variables, and obtaining the gradient based on the adjoint method according to the adjoint variables.
Further, solving an integral form of the jet radiation transfer equation, ignoring wall contribution and atmospheric attenuation, the integral form of the jet radiation transfer equation being as follows:
Wherein, Is the target total irradiance; n is the number of boundary grids except the wall surface; Is the included angle between the normal vector of the infinitesimal plane m and the connecting line of the infinitesimal plane m and the detection point; The solid angle is formed by the infinitesimal plane m to the detection point;
when solving the jet radiation transfer equation in the integral form, rays are emitted from the detection points to boundary grids except the wall surface, the grids through which the light paths from the detection points to the boundary grids pass are calculated through the light path tracking algorithm of the structural grids, the light paths emitted from the detection points to the infinitesimal plane m are recorded as light paths m, and the target radiation illuminance of the light paths needs to be calculated Wherein n is the number of grids through which the optical path m passes,Is the blackbody spectrum radiation brightness; For the gas transmission of the ith cell through which the optical path passes, Is the gas transmissivity of the j-th unit along the light path.
Further, the objective function of the infrared radiation optimization is to minimize the radiation illuminance, and the optimization objective is to:
Wherein, As a function of the minimum value of the function,To be constrained, Q is a flow field variable, D is a design variable,Meaning that Q is a function of D and R is an N-S equation.
Further, lagrangian multipliers are introducedA new objective function L is constructed, since r=0, the value of L is equal toThe values of (2) are the same:
Then L derives the design variable D as:
Wherein, To derive a symbol;
Order the The previous coefficients are zero, deducing the concomitant equation for the infrared radiation problem:
Wherein, The partial derivative of the flow field residual error with respect to the flow field variable is calculated by an N-S equation accompanying solver.
Further, solving the accompanying equation can obtain the accompanying variableGradient calculations based on the accompanying method are:
the key to optimization based on concomitant jet infrared radiation is to derive The expression of (2) givesAnd then leading into a flow field adjoint solver, and solving adjoint equations of the infrared radiation problem to obtain adjoint variables.
Further, the method comprises the steps of,The development form of (2) is as follows:
regarding all grid points along the path m, for the path start grid i=1, there is:
when i=2, there are:
and so on to get Variation of the flow field variable in any grid point k along the path:
Further, calculate In (a) and (b)The process of (1) is as follows:
the calculation formula of the radiation brightness of the mth light path is as follows:
Construction calculation The first term i=1 has:
the second term i=2 is:
The third term i=3 is:
and so on:
I.e. Is given by the recurrence formula of (2)
Further according toObtaining; According toObtaining
Calculation of
Calculation ofWhen the optical path is followed, each grid point i is calculated and storedObtainingAfter that, only need toThe first term of (1), sum1, constructs a recursive formula without processing sum2, for sum1 when i=1: When i=2 has; And so on i=k:
further, the algorithm for tracking the optical path through the structural grid calculates a grid through which the optical path from the detection point to the boundary grid passes, including:
according to the surface division strategy of the structural grid, the grid through which the light path from the detection point to the boundary grid passes is calculated by adopting the following steps:
Step 101: decomposing each blocking grid according to the surface division strategy of the structural grid, and judging whether the characteristic lines intersect all triangles obtained according to the surface division strategy of the structural grid and the positions of the intersection points by using an algorithm;
step 102: traversing grids of a plane where the intersection points are located in the partitioned grids with the intersection points, and finding grids through which the characteristic lines pass;
Step 103: decomposing the grid penetrated by the characteristic line according to the surface division strategy of the structural grid, and solving by using a ray triangle intersection algorithm to obtain the intersection point of the characteristic line and the grid; judging the next grid to which the characteristic line goes according to the logic coordinate relation between the intersection points and the structural grids;
step 104: repeating step 103 until the next grid exceeds the boundary of the partitioned grid; and after all the block grids are calculated, sequencing the obtained intersection points according to time.
Further, the surface division strategy of the structural grid is:
Dividing 6 faces of the structural grid into 12 triangles according to the logical coordinates of the structural grid, wherein three vertexes of each triangle are vertexes of the structural grid, and numbering the triangles according to the logical coordinates of the vertexes of the triangle; when the surface division strategy of the structural grid is adopted, the numbers of the sharing surfaces of the adjacent grids are fixed in the logic coordinate systems i, j and k.
Due to the adoption of the technical scheme, the invention has the following advantages: the method can realize rapid solution of jet infrared radiation, and can realize rapid solution of infrared radiation gradient by solving a pneumatic/infrared radiation coupling variation equation, thereby providing support for pneumatic/infrared coupling optimization based on the accompanying method.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for the description of the embodiments will be briefly described below, and it is apparent that the drawings in the following description are only some embodiments described in the embodiments of the present invention, and other drawings may be obtained according to these drawings for those skilled in the art.
FIG. 1 is a schematic view of an optical path through a structural grid in accordance with an embodiment of the present invention;
Fig. 2 is a schematic diagram of vertex numbering and triangle division of a structural mesh according to an embodiment of the present invention.
Detailed Description
The present invention will be further described with reference to the accompanying drawings and examples, wherein the examples are shown only in a partial, but not in all embodiments of the invention. All other embodiments obtained by those skilled in the art are intended to fall within the scope of the embodiments of the present invention.
Referring to fig. 1, the present invention provides an embodiment of a jet pneumatic/infrared radiation coupling gradient calculation method based on an accompanying method, which includes:
Based on an integral form jet radiation transfer equation, a Lagrangian multiplier is introduced to obtain an accompanying equation of the infrared radiation problem; and solving the adjoint equation of the infrared radiation problem to obtain adjoint variables, and obtaining the gradient based on the adjoint method according to the adjoint variables.
The specific implementation steps of this embodiment are as follows:
solving an integral form of the jet radiation transfer equation, ignoring wall contribution and atmospheric attenuation, the integral form of the jet radiation transfer equation being as follows:
Wherein, Is the target total irradiance; n is the number of boundary grids except the wall surface; Is the included angle between the normal vector of the infinitesimal plane m and the connecting line of the infinitesimal plane m and the detection point; the solid angle of the infinitesimal plane m to the detection point. When solving the jet radiation transfer equation, a light path tracking method is adopted to send out rays from the detection points to the boundary grids except the wall surface, and the grids through which the light paths from the detection points to the boundary grids pass are calculated through a light path tracking algorithm. For the optical path (denoted as optical path m) emitted from the detection point to the infinitesimal m, it is necessary to calculate the optical path Where n is the number of grids traversed by the optical path m,Is the blackbody spectrum radiation brightness; is the gas transmittance of the ith cell through which the optical path passes.
The objective function of the infrared radiation optimization is to minimize the radiation illuminance, and the optimization objective is that:
Wherein, As a function of the minimum value of the function,To be constrained, Q is a flow field variable, D is a design variable,Representing Q as a function of D, R as an N-S equation, incorporating Lagrangian multipliersA new objective function L is constructed, since r=0, the value of L is equal toThe values of (2) are the same:
Then L derives the design variable D as:
Wherein, To derive a symbol;
Order the The previous coefficients are zero, deducing the concomitant equation for the infrared radiation problem:
Wherein, For the partial derivative of the flow field residual error with respect to the flow field variable, the partial derivative is calculated by an N-S equation accompanying solver, and an accompanying variable can be obtained by solving an accompanying equation. Gradient calculation based on the accompanying method is:
the key to optimization based on concomitant jet infrared radiation is to derive The expression of (2) givesAnd then leading into a flow field adjoint solver, and solving adjoint equations of the infrared radiation problem to obtain adjoint variables.
By the variation method of jet infrared radiation illuminance to along Cheng Liuchang variable, deducingThe process of the expression of (2) is:
The calculation formula of the target total radiation illuminance is to sum all light paths, and the target radiation illuminance of the light path m is as follows:
where n is the number of grids traversed by the optical path m (as in figure 1), Is the blackbody spectrum radiation brightness; For the gas transmission of the ith cell through which the optical path passes, Is the gas transmissivity of the j-th unit along the light path.
The development form of (2) is as follows:
regarding all grid points along the path m, for the path start grid i=1, there is:
when i=2, there are:
and so on to get Variation of the flow field variable in any grid point k along the path:
Referring to fig. 2, taking an example that the structural grid has 6 faces, dividing the 6 faces of the structural grid into 12 triangles according to logical coordinates, wherein the triangle 1 consists of points 4, 8 and 5; triangle number 2 consists of points 4, 5, 1; triangle No. 3 is composed of points 5, 6, 2; triangle No. 4 consists of points 5, 2, 1; triangle number 5 consists of points 4,1, 2; triangle 6 is composed of points 4, 2 and 3; triangle 7 is composed of points 3, 7, 6; triangle 8 is composed of points 3, 6 and 2; triangle 9 consists of points 3, 7 and 8; triangle number 10 is composed of points 3, 8, 4; the triangle 11 is composed of points 8, 7 and 6; triangle No. 12 consists of points 8, 6, 5. When the meshing strategy is adopted, the shared surface numbers of adjacent grids in the i, j and k directions are fixed.
Based on the surface division strategy, the optical path tracking method mainly comprises the following 6 steps:
(1) Decomposing each block grid (block) in the above manner and using The algorithm judges whether the characteristic line intersects 12 triangles or not and the position of the intersection point;
(2) Traversing grids of a plane where the intersection points are located in the partitioned grids with the intersection points, and finding grids (meshes) through which the characteristic lines pass;
(3) Decomposing the grid according to the grid decomposition strategy, and solving by using a ray triangle intersection algorithm to obtain a second intersection point of the characteristic line and the grid;
(4) Judging the next grid to which the characteristic line goes according to the logic coordinate relation between the intersection points and the structural grids;
(5) Repeating the steps 3 and 4; until the next grid exceeds the grid boundary of the block;
(6) And after all the block grids are calculated, sequencing the obtained intersection points according to time.
Recursive formula-based radiance variation solving method and calculationIn (a) and (b)The process of (1) is as follows:
the calculation formula of the radiation brightness of the mth light path is as follows:
Construction calculation The first term i=1 has:
the second term i=2 is:
The third term i=3 is:
and so on:
I.e. Is given by the recurrence formula of (2)
The expression of (2) is:
Calculation of When the optical path is followed, each grid point i is calculated and stored. ObtainingAfter that, only need toThe first term (sum 1) of (2) constructs a recurrence formula without special handling of sum2. For sum1, there is when i=1: When i=2 has; And so on i=k:
finally, it should be noted that: the above embodiments are only for illustrating the technical aspects of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the above embodiments, it should be understood by those of ordinary skill in the art that: modifications and equivalents may be made to the specific embodiments of the invention without departing from the spirit and scope of the invention, which is intended to be covered by the claims.

Claims (7)

1. A method for computing a jet aerodynamic/infrared radiation coupling gradient based on a concomitant method, comprising:
Based on an integral form jet radiation transfer equation, a Lagrangian multiplier is introduced to obtain an accompanying equation of the infrared radiation problem; solving an accompanying equation of the infrared radiation problem to obtain an accompanying variable, and obtaining gradient based on an accompanying method according to the accompanying variable;
solving an integral form of the jet radiation transfer equation, ignoring wall contribution and atmospheric attenuation, the integral form of the jet radiation transfer equation being as follows:
Wherein, Is the target total irradiance; n is the number of boundary grids except the wall surface; Is the included angle between the normal vector of the infinitesimal plane m and the connecting line of the infinitesimal plane m and the detection point; The solid angle is formed by the infinitesimal plane m to the detection point;
when solving the jet radiation transfer equation in the integral form, rays are emitted from the detection points to boundary grids except the wall surface, the grids through which the light paths from the detection points to the boundary grids pass are calculated through the light path tracking algorithm of the structural grids, the light paths emitted from the detection points to the infinitesimal plane m are recorded as light paths m, and the target radiation illuminance of the light paths needs to be calculated Wherein n is the number of grids through which the optical path m passes,Is the blackbody spectrum radiation brightness; For the gas transmission of the ith cell through which the optical path passes, The gas transmittance of the jth unit through which the light path passes;
The objective function of the infrared radiation optimization is to minimize the radiation illuminance, and the optimization objective is that:
Wherein, As a function of the minimum value of the function,To be constrained, Q is a flow field variable, D is a design variable,Representing Q as a function of D, R as an N-S equation;
Introducing Lagrangian multipliers A new objective function L is constructed, since r=0, the value of L is equal toThe values of (2) are the same:
Then L derives the design variable D as:
Wherein, To derive a symbol;
Order the The previous coefficients are zero, deducing the concomitant equation for the infrared radiation problem:
Wherein, The partial derivative of the flow field residual error with respect to the flow field variable is calculated by an N-S equation accompanying solver.
2. The method of claim 1, wherein solving the accompanying equation yields an accompanying variableGradient calculations based on the accompanying method are:
the key to optimization based on concomitant jet infrared radiation is to derive The expression of (2) givesAnd then leading into a flow field adjoint solver, and solving adjoint equations of the infrared radiation problem to obtain adjoint variables.
3. The method of claim 1, wherein the step of determining the position of the substrate comprises,The development form of (2) is as follows:
regarding all grid points along the path m, for the path start grid i=1, there is:
when i=2, there are:
and so on to get Variation of the flow field variable in any grid point k along the path:
4. a method according to claim 3, characterized in that the calculation is performed In (a) and (b)The process of (1) is as follows:
the calculation formula of the radiation brightness of the mth light path is as follows:
Construction calculation The first term i=1 has:
the second term i=2 is:
The third term i=3 is:
and so on:
I.e. Is given by the recurrence formula of (2)
5. The method according to claim 4, characterized in thatObtaining; According toObtaining
Calculation of
Calculation ofWhen the optical path is followed, each grid point i is calculated and storedObtainingAfter that, only need toThe first term of (1), sum1, constructs a recursive formula without processing sum2, for sum1 when i=1: When i=2 has; And so on i=k:
6. the method of claim 1, wherein the computing a grid through which the optical path of the probe point to the boundary grid passes through the optical path tracking algorithm of the structural grid comprises:
according to the surface division strategy of the structural grid, the grid through which the light path from the detection point to the boundary grid passes is calculated by adopting the following steps:
Step 101: decomposing each blocking grid according to the surface division strategy of the structural grid, and judging whether the characteristic lines intersect all triangles obtained according to the surface division strategy of the structural grid and the positions of the intersection points by using an algorithm;
step 102: traversing grids of a plane where the intersection points are located in the partitioned grids with the intersection points, and finding grids through which the characteristic lines pass;
Step 103: decomposing the grid penetrated by the characteristic line according to the surface division strategy of the structural grid, and solving by using a ray triangle intersection algorithm to obtain the intersection point of the characteristic line and the grid; judging the next grid to which the characteristic line goes according to the logic coordinate relation between the intersection points and the structural grids;
step 104: repeating step 103 until the next grid exceeds the boundary of the partitioned grid; and after all the block grids are calculated, sequencing the obtained intersection points according to time.
7. The method of claim 6, wherein the surface partitioning strategy of the structural grid is:
Dividing 6 faces of the structural grid into 12 triangles according to the logical coordinates of the structural grid, wherein three vertexes of each triangle are vertexes of the structural grid, and numbering the triangles according to the logical coordinates of the vertexes of the triangle; when the surface division strategy of the structural grid is adopted, the numbers of the sharing surfaces of the adjacent grids are fixed in the logic coordinate systems i, j and k.
CN202410614660.2A 2024-05-17 2024-05-17 Jet flow pneumatic/infrared radiation coupling gradient calculation method based on accompanying method Active CN118211517B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410614660.2A CN118211517B (en) 2024-05-17 2024-05-17 Jet flow pneumatic/infrared radiation coupling gradient calculation method based on accompanying method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410614660.2A CN118211517B (en) 2024-05-17 2024-05-17 Jet flow pneumatic/infrared radiation coupling gradient calculation method based on accompanying method

Publications (2)

Publication Number Publication Date
CN118211517A CN118211517A (en) 2024-06-18
CN118211517B true CN118211517B (en) 2024-07-23

Family

ID=91448776

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410614660.2A Active CN118211517B (en) 2024-05-17 2024-05-17 Jet flow pneumatic/infrared radiation coupling gradient calculation method based on accompanying method

Country Status (1)

Country Link
CN (1) CN118211517B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116956472A (en) * 2023-09-20 2023-10-27 中国空气动力研究与发展中心空天技术研究所 RCS surface sensitivity calculation method for MLFMA concomitant solution
CN117669413A (en) * 2023-11-28 2024-03-08 西安电子科技大学 Offshore foam-containing and ocean wave infrared sea surface simulation method based on physics

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2850647B1 (en) * 2012-05-18 2020-02-26 Micromass UK Limited Orthogonal acceleration coaxial cylinder time of flight mass analyser
US9495516B2 (en) * 2012-11-21 2016-11-15 The Trustees Of Columbia University In The City Of New York Systems, methods, and devices for image reconstruction using combined PDE-constrained and simplified spherical harmonics algorithm
CN105468846B (en) * 2015-11-24 2019-01-11 北京宇航系统工程研究所 A method of rocket bottom hot-fluid is determined using RADIATION ANGLE COEFFICIENT
CN116956483B (en) * 2023-07-10 2024-07-02 北京航空航天大学 Turbine optimization method and system based on S2 flow surface accompanying equation

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116956472A (en) * 2023-09-20 2023-10-27 中国空气动力研究与发展中心空天技术研究所 RCS surface sensitivity calculation method for MLFMA concomitant solution
CN117669413A (en) * 2023-11-28 2024-03-08 西安电子科技大学 Offshore foam-containing and ocean wave infrared sea surface simulation method based on physics

Also Published As

Publication number Publication date
CN118211517A (en) 2024-06-18

Similar Documents

Publication Publication Date Title
CN108681617B (en) An optimal design method for the layout of spacecraft multi-satellite sensors
CN116738891B (en) LU-SGS improvement method for enhancing simulation stability of aircraft flow field
CN116305892B (en) A method, system, device and medium for predicting radiation heat transfer of non-ash gas based on Monte Carlo method
CN118211517B (en) Jet flow pneumatic/infrared radiation coupling gradient calculation method based on accompanying method
CN108470079B (en) Simulation method for radiation safety evaluation of nuclear operation of space launching field
CN113408688A (en) Unknown environment-oriented multi-radioactive source online searching method
CN115480032B (en) Point source emission intensity prediction method based on ground remote sensing measurement
Liao et al. Analysis of flow field aero-optical effects on the imaging by near-earth space all-time short-wave infrared star sensors
CN115861414A (en) Attitude calculation method and system for fixed star observation of remote sensing satellite
CN114329777B (en) A numerical simulation method for ship thermal characteristics
CN107344630B (en) A Method for Determining On-Star Configuration Layout of Multi-field Star Sensors
CN114820912B (en) Simulation method of space-based optical imaging of complex targets based on multi-level distributed beam tracking
CN114398768A (en) Method, system and equipment for determination of OH ultraviolet spectral radiation characteristics in non-equilibrium flow
CN111780745B (en) Short arc ellipse fitting optimization method for deep space exploration optical navigation
Mohammadnejad et al. Baffle simulation and design with the aim of use in space cameras
CN118067571A (en) Solar field angle value inhibition method and system for high-rise atmospheric density in X-ray sun-shading detection
CN117848506A (en) Satellite infrared radiation analysis method considering engine flow field
Li et al. Simulation analysis on the characteristics of aerosol particles to inhibit the infrared radiation of exhaust plumes
CN117289203A (en) Aerial target TDOA iterative positioning method based on satellite radiation source
CN113849953B (en) Micro focal spot device design optimization method for space X-ray communication
CN112926237B (en) Space target key feature identification method based on photometric signals
Baker Mesh generation for the computation of flowfields over complex aerodynamic shapes
Gao et al. An improved reverse Monte Carlo method for the investigation of aerodynamic and infrared radiation characteristics of a flying wing UAV
CN119558013B (en) A method for error correction of turbine blade radiation temperature measurement based on reflected radiation model
CN119046597A (en) Method for improving ray travel tracking success rate

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