CN114021497A - Compressible turbulent flow fluid topology optimization method based on automatic differentiation - Google Patents
Compressible turbulent flow fluid topology optimization method based on automatic differentiation Download PDFInfo
- Publication number
- CN114021497A CN114021497A CN202111294234.8A CN202111294234A CN114021497A CN 114021497 A CN114021497 A CN 114021497A CN 202111294234 A CN202111294234 A CN 202111294234A CN 114021497 A CN114021497 A CN 114021497A
- Authority
- CN
- China
- Prior art keywords
- fluid
- equation
- optimization
- compressible
- adjoint
- 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.)
- Pending
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 93
- 239000012530 fluid Substances 0.000 title claims abstract description 84
- 238000000034 method Methods 0.000 title claims abstract description 72
- 230000004069 differentiation Effects 0.000 title claims abstract description 32
- 238000013461 design Methods 0.000 claims abstract description 27
- 230000035945 sensitivity Effects 0.000 claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims abstract description 17
- VVQNEPGJFQJSBK-UHFFFAOYSA-N Methyl methacrylate Chemical compound COC(=O)C(C)=C VVQNEPGJFQJSBK-UHFFFAOYSA-N 0.000 claims abstract description 11
- 238000005516 engineering process Methods 0.000 claims abstract description 7
- 238000013178 mathematical model Methods 0.000 claims abstract description 7
- 239000000463 material Substances 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 230000014509 gene expression Effects 0.000 claims description 3
- 230000009191 jumping Effects 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 6
- 238000009795 derivation Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 2
- 238000010206 sensitivity analysis Methods 0.000 description 2
- BQCADISMDOOEFD-UHFFFAOYSA-N Silver Chemical compound [Ag] BQCADISMDOOEFD-UHFFFAOYSA-N 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 229910052709 silver Inorganic materials 0.000 description 1
- 239000004332 silver Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- -1 which is neglected Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Data Mining & Analysis (AREA)
- Fluid Mechanics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Computing Systems (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
A compressible turbulent flow fluid topology optimization method based on automatic differentiation relates to topology optimization. Establishing a geometric model corresponding to the topological optimization of the flow channel; obtaining basic parameters of a topology optimization object, and constructing a CFD topology optimization model; constructing an objective function according to the basic parameters to form a nonlinear programming problem; deriving a sensitivity equation based on a companion method; solving a fluid control equation to obtain a flow field result and a value of an objective function, and outputting a new flow field result; constructing a Jacobian matrix and a gradient vector by utilizing an automatic differentiation technology; constructing a matrix form adjoint equation to obtain an adjoint multiplier; substituting the new flow field result and the adjoint multiplier into a sensitivity equation to solve the sensitivity; establishing a mathematical model by using an MMA (methyl methacrylate) numerical optimization method and combining a nonlinear programming problem; and carrying out optimization solution on the flow channel, updating the design variables to obtain an optimal solution, and outputting an optimal two-dimensional topological configuration. The problems that the manual derivation of compressible fluid is not enough in equation flexibility, complex in process, prone to error and the like are solved.
Description
Technical Field
The invention belongs to the technical field of fluid topology optimization in computational fluid mechanics, and particularly relates to a compressible turbulent fluid topology optimization method based on automatic differentiation.
Background
Compressibility is one of the inherent properties of all fluids, usually in water and other liquids, which is neglected, and gases like air, where compressibility is one of the most important considerations at high velocities and pressures. High velocity compressible fluid flow is widely found in the aerospace field, and fluid compressibility must be considered in aerodynamic designs such as engine flowpath parts, aircraft profiles, airfoil profiles, turbomachinery, and the like. Meanwhile, with the development of computer technology, optimization methods for compressible streams have emerged in succession, Cummings et al ([1] Cummings R M, Yang H T, Oh Y H supersonics, turbo flow calculation and drive optimization for axiomatic afterbodies [ J ]. Computers & flows, 1995,24(4):487 507) proposed size optimization for compressible Fluids, and burgreand bay ([2] burgreg w.three-dimensional affinity processing optimization using size optimization.
The fluid topology optimization method is different from size optimization and shape optimization, can realize the change of the topological configuration of an optimized object, and has great advantages in the initial concept design stage. Borrval and Petersson ([3] Borrval T, Petersson J. topology optimization of Fluids in stocks flow [ J ]. International Journal for Numerical Methods in Fluids,2003,41(1):77-107) have for the first time proposed a fluid topology optimization method and applied to the topology optimization problem of fluid channels. Since then, the fluid topology optimization method is widely applied to various flows. While the current fluid topology optimization is mainly focused on the field of incompressible flow. Recently, LFN S et al ("4 LFN S a, Okubo C M, silver e. poling Optimization of subsonic compressible flows [ J ] Structural and Multidisciplinary Optimization,2021(6)) achieved subsonic compressible fluid topology Optimization, but no description was given of the specific process of auto-differentiating to achieve dispersion, and no consideration was given to the effects of turbulence. In many scientific research and engineering fields, the working environments of wing airfoils, turbine blades, internal flow channels of engines and the like are high-speed compressible turbulent flow. Therefore, it is highly desirable to develop a compressible fluid topology optimization method that takes into account the turbulence model.
The solution of sensitivity in fluid topology optimization usually adopts a continuous or discrete adjoint method, and A. Carnarius et al. ([5] Carnarius A, Thiele F, Oezkaya E, et al. Adjoint errors For Optimal Flow Control [ C ]//5th Flow Control reference.2010) introduces the advantages and disadvantages of the two methods. For complex CFD solvers that employ separation and iterative solution methods, manually deriving discrete adjoint equations is an important and error-prone task. The automatic differentiation method is a technology for accurately solving the function gradient by using a computer program, and has the characteristics of strong flexibility, accurate calculation and the like. The automatic differentiation method is described in detail in R.general. ([6] R.general, T.Kaminski, strategies for additive code construction, ACM Trans.Math.Softw.24(4) (1998) 437-474). The Cetin B. Dilgen et al. ([7] Dilgen CB, Dilgen S B, Fuhrman D R, et al. topology optimization of turboflow [ J ]. Computer Methods in Applied Mechanics and Engineering,2018,331(APR.1):363-393) proposed the optimization of the incompressible turbulence topology based on the automatic differentiation method. The invention discloses a compressible fluid accompanying equation solving method by utilizing an automatic differentiation method, and provides a compressible fluid topological optimization method based on automatic differentiation and considering a turbulence model.
Disclosure of Invention
The invention aims to provide a compressible turbulent fluid topology optimization method based on automatic differentiation, which can realize fluid topology optimization aiming at compressible fluid, considers the influence of a turbulent model on the optimization, constructs an adjoint equation through an automatic differentiation method, efficiently and accurately displays and generates a Jacobian matrix by utilizing the characteristics of the automatic differentiation method, and avoids the problems of insufficient flexibility, complicated process, easy error and the like of manually deriving the adjoint equation of the compressible fluid.
The invention comprises the following steps:
step 1: establishing a geometric model corresponding to the topological optimization of the flow channel according to the actual requirement of the optimized object, and determining the positions of an optimized design region omega and the inlet and outlet of the flow channel;
step 2: obtaining basic parameters of a topological optimization object, and constructing a CFD topological optimization model based on a compressible fluid control equation;
and step 3: converting optimization objective and fluid control equation into weak form according to basic parametersConstructing a target function to form a nonlinear programming problem;
and 4, step 4: deriving a sensitivity equation based on a companion method;
and 5: setting appropriate design variables and initial values of intermediate variables, solving a fluid control equation to obtain a flow field result and a value of an objective function, and outputting a new flow field result;
step 6: in order to solve the adjoint equation, a Jacobian matrix and a gradient vector are constructed by utilizing an automatic differentiation technology;
and 7: constructing a matrix form adjoint equation by the Jacobian matrix and the gradient vector obtained in the step 6 to obtain an adjoint multiplier;
and 8: substituting the new flow field result obtained in the step 5 and the adjoint multiplier obtained in the step 7 into the sensitivity equation in the step 4, and solving the sensitivity;
and step 9: using an MMA (methyl methacrylate) numerical optimization method, and building a mathematical model by combining a known topological optimization nonlinear programming problem;
step 10: performing optimization solution on the flow channel by using an MMA algorithm, updating design variables, obtaining an optimal solution if iteration is converged, and outputting an optimal two-dimensional topological configuration; otherwise, jumping to step 5, and setting the flow field result as a new initial value.
In step 2, the obtaining of the basic parameters of the topology optimization object mainly includes: optimizing targets, fluid properties, inflow conditions, boundary conditions, heat conduction parameters; the specific steps of constructing the CFD topology optimization model based on the compressible fluid control equation may be:
(1) the material density is designated as a design variable gamma based on a variable density topological optimization method, a source term f is added in a fluid control equation to represent the resistance brought by the porous medium, and the function of the resistance and the design variable is as follows:
f=-au
α(γ)=αmin+(αmax-αmin)q(1-γ)/(q+γ)
wherein alpha is a permeability coefficient, and the larger the alpha value is, the larger the resistance caused by the porous medium to the fluid is; u is the fluid velocity field, q is a penalty parameter, and γ is the material density.
(2) Constructing a CFD topological optimization model based on a compressible fluid control equation:
Min:J(u,p,T,γ);
0≤γ≤1
wherein E is energy density, ρ is mass density, u is fluid velocity field, p is fluid pressure field, T is fluid temperature field, k is thermal conductivity coefficient,is the viscosity tensor, μ is the kinetic viscosity, μTSolving the resulting turbulence viscosity, V, for the turbulence modelθThe volume fraction upper limit.
In step 3, the specific steps of constructing the objective function and constructing the nonlinear programming problem may be: converting optimization objective and fluid control equation into weak form according to basic parametersAnd construct an objective functionConverting volume fraction limits to constraintsWherein, a design variable γ, i.e. a function independent variable x, then constitutes the following nonlinear programming problem:
Minimize f0(x)
Subject to fi(x)≤0,i=1,…,m
in step 4, the sensitivity equation is derived based on the adjoint method as follows:
wherein, U is an intermediate variable comprising U, p and T;
in step 6, the specific method of the jacobian matrix and the gradient vector may be:
(1) assembling a jacobian matrix with a highly sparse structure by an automatic differential method:the expansion form is as follows:
wherein: m is the number of intermediate variables, and n is the dimension of the design variable.
(2) Assembling a gradient vector L by an automatic differential method, wherein the value of the gradient vector L is the partial derivative of the target function to the intermediate variable:
in step 7, the matrix form is accompanied by the equation:
and solving the equation to obtain a value of the adjoint multiplier, and outputting the adjoint multiplier lambda.
In step 8, the solution sensitivity is as follows;
in step 9, the key parameters of the mathematical model include flow field results u, p, T, accompanying multiplier λ, sensitivityAnd first-order and second-order differential expressions of the objective function and the constraint condition, and the like.
In step 10, it is determined whether the iteration has converged, and the maximum relative design variable change value of each grid cell before and after updating is compared with a set iteration termination condition as a convergence criterion.
The method aims at the compressible fluid to carry out fluid topology optimization, considers a turbulence model in a fluid control equation, and solves an adjoint equation in sensitivity analysis by adopting an automatic differentiation method. The compressible Navier-Stokes equation is coupled with a penalty equation based on a variable density material; and carrying out fluid topology optimization according to set constraint conditions such as a target function, volume limitation and the like, wherein the gradient in the optimization algorithm realizes discrete accompaniment through automatic differentiation, and then solving is carried out. The method has the advantages that the automatic differentiation method is utilized to realize and simplify the sensitivity solution of compressible fluid topology optimization under the consideration of the turbulence model, the burden of software development and maintainers is reduced, and the maintainability and flexibility of the fluid topology optimization application program are improved. The compressible turbulent flow fluid topology optimization method based on automatic differentiation provides a solving means for fluid topology optimization application in the engineering fields of aerospace and the like. The invention is suitable for the fluid topology optimization of compressible laminar flow and turbulent flow.
Drawings
FIG. 1 is a flowchart of an embodiment of a compressible fluid topology optimization method based on automatic differential solution with sensitivity.
Fig. 2 is a schematic view of a topology optimization design of a compressible turbulent two-dimensional straight-through flow channel provided by an embodiment of the present invention.
FIG. 3 is a graph of the convergence of the objective function and the volume constraint during the optimization process according to an embodiment of the present invention. Wherein, (a) is an objective function convergence diagram; (b) is a volume fraction convergence diagram.
Fig. 4 is a diagram of an optimal flow channel topology after solving according to an embodiment of the present invention.
Detailed Description
In the following description, for purposes of explanation and not limitation, specific details are set forth in order to provide a thorough understanding of the present invention.
The compressible fluid topological optimization method based on the automatic differential solution with the sensitivity is mainly constructed in a flow chart shown in figure 1 by taking the compressible turbulence optimization in a simple two-dimensional straight channel as an embodiment. The main method of the embodiment is as follows:
1. a square design field omega is constructed as required by the straight channel embodiment, fig. 2. The left side is an inlet, the right side is an outlet, and the width of the inlet is twice of that of the outlet;
2. setting basic parameters of a topology optimization object, mainly comprising: and (3) specifying the material density as gamma of a design variable, an optimization target, fluid properties, inflow conditions, boundary conditions and the like based on the variable density topological optimization method. Wherein the optimization objective is flow loss, the fluid property is air, and the boundary conditions are no-slip and thermal insulation boundaries:
designing variables: material density gamma
Intermediate variables: u, p, T
an objective function:wherein μ is dynamic viscosity, μTSolving the obtained turbulence viscosity for the turbulence model;
constraint conditions are as follows:wherein, VθIs an artificial limit value for the volume fraction.
In this example, the reynolds number Re is 2.6 × 106The inlet velocity is parabolic and umax167m/s, Mach number Ma 0.5, volume constraint upper limit Vθ=0.8。
3. The constraint conditions of the fluid topology optimization are divided into equality constraint conditions and inequality constraint conditions, wherein the equality constraint conditions comprise a fluid control equation after a source term f ═ au is added:
the inequality constraints include volume fraction upper limits:in summary, the optimization objective and the fluid control equation are transformed into weak formAnd construct an objective functionConverting volume fraction limits to constraintsWherein, a design variable γ, i.e. a function independent variable x, then constitutes the following nonlinear programming problem:
Minimize f0(x)
Subject to fi(x)≤0,i=1,…,m
4. deriving a sensitivity equation based on the adjoint method:
wherein, U is an intermediate variable comprising U, p and T;
5. setting appropriate design variables and initial values of intermediate variables, solving a fluid control equation to obtain a flow field result and a value of an objective function, and outputting fluid velocity fields u, p and J (u, p, T and gamma);
a) assembling elegance with highly sparse structure by auto-differential methodComparable matrices:its unfolded shape
The formula is as follows:
wherein: m is the number of intermediate variables, and n is the dimension of the design variable.
b) Assembling a gradient vector L by an automatic differential method, wherein the value of the gradient vector L is the partial derivative of the target function to the intermediate variable:
7. and (3) forming a matrix form adjoint equation by the Jacobian matrix K obtained in the step 6 and the gradient vector L:
solving an equation to obtain a value of an adjoint multiplier, and outputting lambda;
8. substituting the flow field result u obtained in the step 4 and the adjoint multiplier lambda obtained in the step 7 into the sensitivity equation in the step 5, and solving the sensitivity;
9. and (3) constructing a mathematical model by using an MMA (methyl methacrylate) numerical optimization method and combining a known topological optimization nonlinear programming problem. The key parameters of the mathematical model comprise flow field results obtained by solving the steps: u, p, T, adjoint multiplier: λ, sensitivity:and first-order and second-order differential expressions of the objective function and the constraint condition, and the like.
10. And (3) performing optimization solution on the straight-through flow channel through an MMA algorithm, continuously updating design variables, comparing the maximum relative design variable change value of each grid unit before and after updating with a set iteration termination condition as a convergence basis, and if the iteration is not converged, jumping to the step 5 and setting a flow field result as a new initial value. And stopping calculation until iterative convergence to obtain an optimal solution, and outputting the two-dimensional straight-through flow channel topological configuration with the lowest energy loss. In the optimization process, penalty parameters q and alphamaxDetermines the quality of the optimized result, namely q and alphamaxWhen the selection is not proper, the problem of unclear fluid-solid boundary can occur. FIG. 3 shows an example of the present invention where q is 0.1 and α ismaxThe objective function and the volume constraint convergence map under 10000 parameter settings. As can be seen from the figure, the change amplitude of the first 10 steps in the optimization process is large, and gradually converges from the 10 th step to the 60 th step. Wherein, the objective function convergence diagram (a) shows the topological configuration under different iteration steps. Fig. 4 is the result of the topology configuration after the final optimization convergence.
The invention provides a compressible fluid topology optimization method based on automatic differentiation, which is characterized in that an automatic differentiation technology is utilized to solve an adjoint equation in sensitivity analysis. The method not only realizes the compressible fluid topology optimization considering the turbulence influence, but also has the characteristic of realizing the discrete adjoint by automatic differentiation, so that repeated derivation is not needed according to target functions and boundary conditions of different research objects, and the method has obvious superiority compared with a fussy manual derivation adjoint equation method.
It should be understood that the above examples are only for clarity of illustration and are not intended to limit the embodiments. This need not be, nor should it be exhaustive of all examples. Any modification and improvement within the spirit of the invention are considered to be within the scope of the invention.
Claims (10)
1. An automatic differential based compressible turbulent fluid topology optimization method is characterized by comprising the following steps:
step 1: establishing a geometric model corresponding to the topological optimization of the flow channel according to the actual requirement of the optimized object, and determining the positions of an optimized design region omega and the inlet and outlet of the flow channel;
step 2: obtaining basic parameters of a topological optimization object, and constructing a CFD topological optimization model based on a compressible fluid control equation;
and step 3: converting optimization objective and fluid control equation into weak form according to basic parametersConstructing a target function to form a nonlinear programming problem;
and 4, step 4: deriving a sensitivity equation based on a companion method;
and 5: setting appropriate design variables and initial values of intermediate variables, solving a fluid control equation to obtain a flow field result and a value of an objective function, and outputting a new flow field result;
step 6: in order to solve the adjoint equation, a Jacobian matrix and a gradient vector are constructed by utilizing an automatic differentiation technology;
and 7: constructing a matrix form adjoint equation by the Jacobian matrix and the gradient vector obtained in the step 6 to obtain an adjoint multiplier;
and 8: substituting the new flow field result obtained in the step 5 and the adjoint multiplier obtained in the step 7 into the sensitivity equation in the step 4, and solving the sensitivity;
and step 9: using an MMA (methyl methacrylate) numerical optimization method, and building a mathematical model by combining a known topological optimization nonlinear programming problem;
step 10: performing optimization solution on the flow channel by using an MMA algorithm, updating design variables, obtaining an optimal solution if iteration is converged, and outputting an optimal two-dimensional topological configuration; otherwise, jumping to step 5, and setting the flow field result as a new initial value.
2. The compressible turbulent fluid topological optimization method based on automatic differentiation according to claim 1, wherein in step 2, the obtaining of basic parameters of a topological optimization object comprises: optimization objectives, fluid properties, inflow conditions, boundary conditions, heat transfer parameters.
3. The compressible turbulent fluid topology optimization method based on automatic differentiation according to claim 1, wherein in step 2, the concrete steps of constructing the CFD topology optimization model based on the compressible fluid control equation are as follows:
(1) the material density is designated as a design variable gamma based on a variable density topological optimization method, a source term f is added in a fluid control equation to represent the resistance brought by the porous medium, and the function of the resistance and the design variable is as follows:
f=-au
α(γ)=αmin+(αmax-αmin)q(1-γ)/(q+γ)
wherein alpha is a permeability coefficient, and the larger the alpha value is, the larger the resistance caused by the porous medium to the fluid is; u is a fluid velocity field, q is a penalty parameter, and gamma is a material density;
(2) constructing a CFD topological optimization model based on a compressible fluid control equation:
Min:J(u,p,T,γ);
0≤γ≤1
wherein E is the energy density and ρ isMass density, u is the fluid velocity field, p is the fluid pressure field, T is the fluid temperature field, k is the thermal conductivity,is the viscosity tensor, μ is the kinetic viscosity, μTSolving the resulting turbulence viscosity, V, for the turbulence modelθThe volume fraction upper limit.
4. The compressible turbulent fluid topology optimization method based on automatic differentiation according to claim 1, wherein in step 3, the specific steps of constructing the objective function and constructing the nonlinear programming problem are as follows: converting optimization objective and fluid control equation into weak form according to basic parametersAnd construct an objective functionConverting volume fraction limits to constraintsWhere x represents the design variable γ, then constitutes the following nonlinear programming problem:
Minimize f0(x)
Subject to fi(x)≤0,i=1,…,m。
5. the method for optimizing the topology of the compressible turbulent fluid based on automatic differentiation according to claim 1, wherein in step 4, the sensitivity equation is derived based on the adjoint method as follows:
wherein, U is an intermediate variable comprising U, p and T;
6. the compressible turbulent fluid topology optimization method based on automatic differentiation according to claim 1, wherein in step 6, the specific method for constructing the jacobian matrix and the gradient vector by using the automatic differentiation technology is as follows:
(1) assembling a jacobian matrix with a highly sparse structure by an automatic differential method:the expansion form is as follows:
wherein: m is the number of intermediate variables, and n is the dimension of a design variable;
(2) assembling a gradient vector L by an automatic differential method, wherein the value of the gradient vector L is the partial derivative of the target function to the intermediate variable:
7. the method for optimizing the topology of the compressible turbulent fluid based on automatic differentiation as claimed in claim 1, wherein in step 7, the accompanying equation in the form of a matrix is as follows:
and solving the equation to obtain a value of the adjoint multiplier, and outputting the adjoint multiplier lambda.
9. the method for optimizing the topology of the compressible turbulent fluid based on automatic differentiation according to claim 1, wherein in step 9, the key parameters of the mathematical model comprise flow field results u, p, T, sensitivity with multiplier λAnd first-order and second-order differential expressions of the objective function and the constraint condition.
10. The method for optimizing the topology of the compressible turbulent fluid based on automatic differentiation according to claim 1, wherein in step 10, the iteration is converged, and the maximum relative design variable change value of each grid cell before and after the updating is compared with the set iteration termination condition as a convergence criterion.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111294234.8A CN114021497A (en) | 2021-11-03 | 2021-11-03 | Compressible turbulent flow fluid topology optimization method based on automatic differentiation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111294234.8A CN114021497A (en) | 2021-11-03 | 2021-11-03 | Compressible turbulent flow fluid topology optimization method based on automatic differentiation |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114021497A true CN114021497A (en) | 2022-02-08 |
Family
ID=80060283
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111294234.8A Pending CN114021497A (en) | 2021-11-03 | 2021-11-03 | Compressible turbulent flow fluid topology optimization method based on automatic differentiation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114021497A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114896747A (en) * | 2022-05-30 | 2022-08-12 | 四川启睿克科技有限公司 | Sensitivity calculation-based micro-channel structure optimization design method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112231900A (en) * | 2020-10-09 | 2021-01-15 | 中国科学院力学研究所 | Variable physical property topology optimization method |
AU2020103808A4 (en) * | 2020-01-17 | 2021-02-11 | Beijing University Of Technology | A design method of the fail-safe topology optimization of continuum structures with the frequency and displacement constraints |
CN113094944A (en) * | 2021-03-22 | 2021-07-09 | 华中科技大学 | Micro-channel radiator and microscopic scale topology optimization design method thereof |
-
2021
- 2021-11-03 CN CN202111294234.8A patent/CN114021497A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2020103808A4 (en) * | 2020-01-17 | 2021-02-11 | Beijing University Of Technology | A design method of the fail-safe topology optimization of continuum structures with the frequency and displacement constraints |
CN112231900A (en) * | 2020-10-09 | 2021-01-15 | 中国科学院力学研究所 | Variable physical property topology optimization method |
CN113094944A (en) * | 2021-03-22 | 2021-07-09 | 华中科技大学 | Micro-channel radiator and microscopic scale topology optimization design method thereof |
Non-Patent Citations (1)
Title |
---|
CHENYUE XIE: "Artificial neural network mixed model for large eddy simulation of compressible isotropic turbulence", PHYSICS OF FLUIDS, 20 August 2019 (2019-08-20) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114896747A (en) * | 2022-05-30 | 2022-08-12 | 四川启睿克科技有限公司 | Sensitivity calculation-based micro-channel structure optimization design method |
CN114896747B (en) * | 2022-05-30 | 2024-05-14 | 四川启睿克科技有限公司 | Microchannel structure optimization design method based on sensitivity calculation |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Giles et al. | Two-dimensional transonic aerodynamic design method | |
Nadarajah et al. | Survey of shape parameterization techniques and its effect on three-dimensional aerodynamic shape optimization | |
Ramanujam et al. | Improving airfoil drag prediction | |
CN107944172B (en) | Engine intake manifold design method based on fluid topology optimization | |
Yildirim et al. | Aeropropulsive design optimization of a boundary layer ingestion system | |
CN111709096A (en) | Design method of special-shaped fin structure for strengthening natural convection heat transfer | |
CN116305523B (en) | Laminar flow airfoil optimization method based on discrete accompaniment | |
Nili-Ahmadabadi et al. | Development and validation of a hybrid aerodynamic design method for curved diffusers using genetic algorithm and ball-spine inverse design method | |
CN105844025B (en) | A kind of non-probability thermographic curve reliability design approach for hypersonic rudder face | |
CN114021497A (en) | Compressible turbulent flow fluid topology optimization method based on automatic differentiation | |
Ruh et al. | Robust modeling and optimal design of rotors using blade element momentum theory | |
Yildirim et al. | Performance analysis of optimized STARC-ABL designs across the entire mission profile | |
Giugno et al. | Adjoint-based optimization of a modern jet-engine fan blade | |
Barrett et al. | Airfoil shape design and optimization using multifidelity analysis and embedded inverse design | |
Divo et al. | Glenn-ht/bem conjugate heat transfer solver for large-scale turbomachinery models | |
CN116029003A (en) | Binary air inlet channel design method based on fluid topology optimization | |
Wernick et al. | Rotor Blade Design Optimization for Boundary Layer Ingesting Inlet Fan | |
CN111413867B (en) | Rapid modeling and optimization method for equivalent controller of hydraulic control mechanism on Flowmaster platform | |
Ling et al. | Flow field fusion simulation method based on model features and its application in CRDM | |
Wang et al. | Performance Prediction of an Optimized Centrifugal Pump with High Efficiency. | |
Broatch et al. | An automatized methodology to generate velocity distortion panels for wind tunnel testing | |
Leng et al. | Variable-fidelity surrogate model based on transfer learning and its application in multidisciplinary design optimization of aircraft | |
Wang et al. | Aerodynamic Optimization Design of Airfoil Shape Using Entropy Generation as an Objective | |
Zhang et al. | Mixed Fidelity Aerodynamics and Aero-Structural Optimization for Wings | |
Lucio | Design of a high subsonic nozzle for a transonic linear cascade operating at near stall conditions |
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 |