25 results
Search Results
2. Discrete Adjoint Optimization Method for Low-Boom Aircraft Design Using Equivalent Area Distribution.
- Author
-
Ma, Chuang, Huang, Jiangtao, Li, Daochun, Deng, Jun, Liu, Gang, Zhou, Lin, and Chen, Cheng
- Subjects
AERODYNAMIC load ,DRAG coefficient ,DRAG force ,AIRPLANE motors ,EQUATIONS - Abstract
This paper introduces a low-boom aircraft optimization design method guided by equivalent area distribution, which effectively improves the intuitiveness and refinement of inverse design. A gradient optimization method based on discrete adjoint equations is proposed to achieve the fast solution of the gradient information of target equivalent area distribution relative to design variables and to drive the aerodynamic shape update to the optimal solution. An optimization experiment is carried out based on a self-developed supersonic civil aircraft configuration with engines. The results show that the equivalent area distribution adjoint equation can accurately solve the gradient information. After optimization, the sonic boom level of the aircraft was reduced by 13.2 PLdB, and the drag coefficient was reduced by 60.75 counts. Moreover, the equivalent area distribution adjoint optimization method has outstanding advantages, such as high sensitivity and fast convergence speed, and can take both the low sonic boom and the low drag force of the aircraft into account, providing a powerful tool for the comprehensive optimization design of supersonic civil aircraft by considering sonic boom and aerodynamic force. [ABSTRACT FROM AUTHOR]
- Published
- 2024
- Full Text
- View/download PDF
3. Lagrange-Hamilton Approach in Optimization Problems with Isoperimetric-Type Constraints.
- Author
-
Treanţă, Savin
- Subjects
ADJOINT differential equations ,INTEGRAL functions ,EQUATIONS of state ,HAMILTON-Jacobi equations ,LAGRANGE multiplier - Abstract
This paper derives necessary optimality conditions for a certain class of optimal control problems without linearity or convexity assumptions. The optimal control problem has a general objective function of integral type and a finite number of isoperimetric type constraints. For proving the main result derived in this paper, the Lagrange function and the control Hamiltonian are introduced and an adjoint differential equation is stated. In addition, we formulate some examples where the derived necessary optimality conditions are applied. [ABSTRACT FROM AUTHOR]
- Published
- 2022
- Full Text
- View/download PDF
4. Determination of a Time-Varying Point Source in Cauchy Problems for the Convection–Diffusion Equation.
- Author
-
Georgiev, Slavi and Vulkov, Lubin
- Subjects
TRANSPORT equation ,CAUCHY problem ,INVERSE problems ,WATER pollution ,POROUS materials ,ADJOINT differential equations ,PROBLEM solving - Abstract
Featured Application: This work could be readily used in modeling contaminant transportation in homogeneous porous media with constant mean transport velocity, first-order decay and linear equilibrium sorption, e.g., the spread of a solute in water. In particular, the proposed method recovers the time-dependent pollution source strength from point measurements in a(n) (un)bounded domain. In this paper, we suggest a method for recovering the unknown time-dependent strength of a contaminant concentration source from measurements of the concentration inside an unbounded domain. This problem is formulated as a Cauchy parabolic inverse problem. For its efficient numerical processing, the problem is solved by reduction of the Cauchy problem to a Dirichet one on a bounded domain using the method of the fundamental (potential) solutions in combination with an adjoint equation technique. A numerical solution to this approach is explained. Next, by choosing the source strength in the form of a finite series of shape functions with unknown constant coefficients and using a linear-square method, the term concentration source is estimated. Computational simulations using model examples from water pollution are discussed. [ABSTRACT FROM AUTHOR]
- Published
- 2023
- Full Text
- View/download PDF
5. Non-zero-sum differential games of delayed backward doubly stochastic systems and their application.
- Author
-
Jie Xu, Rui Zhang, and Ruiqiang Lin
- Subjects
DIFFERENTIAL games ,STOCHASTIC systems ,STOCHASTIC differential equations ,NASH equilibrium ,DIFFERENTIAL equations ,DUALITY theory (Mathematics) ,MAXIMUM principles (Mathematics) - Abstract
In this paper, we investigate the non-zero-sum differential game problem based on the delayed backward doubly stochastic differential equation. We discuss the case where the controlled system contains time-delayed variables, and the time-delayed variables are different. We construct a new kind of adjoint equation consisting of a doubly stochastic differential equation and three simple differential equations. Under appropriate premises for such equations, we deduce a stochastic maximum principle as a necessary condition and a verification theorem for the Nash equilibrium point by using the duality method and convex variation technique. As an application, we apply our result to a linear delayed backward doubly stochastic non-zero-sum game problem to verify the effectiveness of our results. [ABSTRACT FROM AUTHOR]
- Published
- 2024
- Full Text
- View/download PDF
6. A Functional Approach to Interpreting the Role of the Adjoint Equation in Machine Learning.
- Author
-
Fekete, Imre, Molnár, András, and Simon, Péter L.
- Abstract
The connection between numerical methods for solving differential equations and machine learning has been revealed recently. Differential equations have been proposed as continuous analogues of deep neural networks, and then used in handling certain tasks, such as image recognition, where the training of a model includes learning the parameters of systems of ODEs from certain points along their trajectories. Treating this inverse problem of determining the parameters of a dynamical system that minimize the difference between data and trajectory by a gradient-based optimization method presents the solution of the adjoint equation as the continuous analogue of backpropagation that yields the appropriate gradients. The paper explores an abstract approach that can be used to construct a family of loss functions with the aim of fitting the solution of an initial value problem to a set of discrete or continuous measurements. It is shown, that an extension of the adjoint equation can be used to derive the gradient of the loss function as a continuous analogue of backpropagation in machine learning. Numerical evidence is presented that under reasonably controlled circumstances the gradients obtained this way can be used in a gradient descent to fit the solution of an initial value problem to a set of continuous noisy measurements, and a set of discrete noisy measurements that are recorded at uncertain times. [ABSTRACT FROM AUTHOR]
- Published
- 2024
- Full Text
- View/download PDF
7. Parallel Implementation of a Sensitivity Operator-Based Source Identification Algorithm for Distributed Memory Computers.
- Author
-
Penenko, Alexey and Rusin, Evgeny
- Subjects
DISTRIBUTED algorithms ,INVERSE problems ,OPERATOR equations ,AIR quality ,MESSAGE passing (Computer science) ,PARALLEL algorithms - Abstract
Large-scale inverse problems that require high-performance computing arise in various fields, including regional air quality studies. The paper focuses on parallel solutions of an emission source identification problem for a 2D advection–diffusion–reaction model where the sources are identified by heterogeneous measurement data. In the inverse modeling approach we use, a source identification problem is transformed to a quasi-linear operator equation with a sensitivity operator, which allows working in a unified way with heterogeneous measurement data and provides natural parallelization of numeric algorithms by concurrent calculation of the rows of a sensitivity operator matrix. The parallel version of the algorithm implemented with a message passing interface (MPI) has shown a 40× speedup on four Intel Xeon Gold 6248R nodes in an inverse modeling scenario for the Lake Baikal region. [ABSTRACT FROM AUTHOR]
- Published
- 2022
- Full Text
- View/download PDF
8. Near‐optimal control of a stochastic model for mountain pine beetles with pesticide application.
- Author
-
Xu, Jing, Yu, Zhensheng, Zhang, Tonghua, and Yuan, Sanling
- Subjects
MOUNTAIN pine beetle ,PESTICIDES ,DROUGHT management ,STOCHASTIC models - Abstract
In this paper, considering the effect of environmental factors such as climate change, temperature and drought on the growth of forests and beetles as well as the pine forests death and economic loss caused by outbreaks of mountain pine beetles, we first develop a stochastic mountain pine beetle model with pesticide application, and then consider its near‐optimal control (NOC) problem. The problem NOC is formulated by decreasing the mountain pine beetles while keeping the cost of pesticide application to a minimum. By Pontryagin stochastic maximum principle, we establish some sufficient and necessary conditions of the near‐optimality. Some numerical simulations are carried out to illustrate the theoretical results, which indicate that the pesticide application could effectively control the number of beetles and relieve the pressure from the attacking of beetles on the pines. To some extent, this control strategy could protect the pine forests and reduce economic losses caused by the mountain pine beetles. [ABSTRACT FROM AUTHOR]
- Published
- 2022
- Full Text
- View/download PDF
9. Variational Data Assimilation for the Sea Thermodynamics Model and Sensitivity of Marine Characteristics to Observation Errors
- Author
-
Shutyaev, V. P. and Parmuzin, E. I.
- Published
- 2023
- Full Text
- View/download PDF
10. A Mean-Field Optimal Control for Fully Coupled Forward-Backward Stochastic Control Systems with Lévy Processes
- Author
-
Huang, Zhen, Wang, Ying, and Wang, Xiangrong
- Published
- 2022
- Full Text
- View/download PDF
11. Near‐optimal control of a stochastic pine wilt disease model with prevention strategies.
- Author
-
Xu, Jing and Yuan, Sanling
- Subjects
- *
CONIFER wilt , *WILT diseases , *PINEWOOD nematode , *MEDICAL model , *HAMILTON'S principle function , *PREVENTIVE medicine , *INSECTICIDE application , *PINACEAE - Abstract
Pine wood nematode disease has become nowadays a major threat to forest ecosystems and caused huge economic and environmental losses to many countries and regions. Considering the influence of environmental factors such as temperature, drought, and humidity, as well as possible prevention strategies, including trunk injection, cutting infected pines, and insecticide application, in this paper, we develop a stochastic pine wilt disease model with prevention strategies, aiming to obtain the conditions for its near‐optimal control. The near‐optimality problem is constructed by decreasing the infected pine forests and beetle population while minimizing the costs of prevention strategies. We first verify that, under certain convexity conditions, the near‐minimum condition of the Hamiltonian function is sufficient for the near‐optimal control. Then, by utilizing the spike variational technique, we derive the necessary conditions of the near‐optimal control. Finally, we present an example to illustrate the validity of the main theoretical results, which indicates that reasonable prevention strategies could effectively reduce pine wilt infection, and therefore, to some extent, could protect the pine forests from pine wilt disease and reduce economic losses caused by this disease. [ABSTRACT FROM AUTHOR]
- Published
- 2023
- Full Text
- View/download PDF
12. Inverse Boundary Value Problem Solution for Deflected Beams Joined Together by Elastic Medium.
- Author
-
Alomari, Omar and Dubovski, Pavel B.
- Subjects
- *
BOUNDARY value problems , *EULER-Bernoulli beam theory , *ADJOINT differential equations - Abstract
In this paper, we extend the Euler-Bernoulli beam theory for bending boundary value problem into mechanically coupled system. We follow the inverse approach to find the exerted force on two beams separated by elastic material. The theory was utilized in two ways: in the first approach, we calculate the force exerted on the beams using known values for the stiffness constant and measured values for the beam deflections. In the second method, we calculate the stiffness constant using a single known force and measured deflections. These problems are typically illposed problems whose solution does not depend continuously on the boundary data. To minimize the variational functional, we develop an iterative algorithm based on the system of three equations: the direct, adjoint, and control equations. Then, we present numerical examples to obtain the solutions. [ABSTRACT FROM AUTHOR]
- Published
- 2023
- Full Text
- View/download PDF
13. The general maximum principle for discrete-time stochastic control problems.
- Author
-
Song, Yuanzhuo and Wu, Zhen
- Subjects
- *
NOISE control , *CONVEX domains , *DISCRETE-time systems - Abstract
In this paper, we discuss the problem of discrete-time stochastic control with multiplicative noise, where the control domains may not be convex. We propose a novel approach inspired by the classical spike variation in continuous-time cases, but adapted to the discrete-time setting by perturbing in a random scale instead of time scale. Our approach allows us to obtain the maximum principle recursively without the need for variation equations or adjoint equations, which are typically required in previous methods. Moreover, we provide two maximum principles, one for the general case and the other when feedback forms exist. Our approach covers the previous results when the control domains are convex. Finally, we demonstrate the effectiveness of our new results through an example. [ABSTRACT FROM AUTHOR]
- Published
- 2024
- Full Text
- View/download PDF
14. Discrete Adjoint Optimization Method for Low-Boom Aircraft Design Using Equivalent Area Distribution
- Author
-
Chuang Ma, Jiangtao Huang, Daochun Li, Jun Deng, Gang Liu, Lin Zhou, and Cheng Chen
- Subjects
supersonic civil aircraft ,sonic boom ,adjoint equation ,optimization ,Motor vehicles. Aeronautics. Astronautics ,TL1-4050 - Abstract
This paper introduces a low-boom aircraft optimization design method guided by equivalent area distribution, which effectively improves the intuitiveness and refinement of inverse design. A gradient optimization method based on discrete adjoint equations is proposed to achieve the fast solution of the gradient information of target equivalent area distribution relative to design variables and to drive the aerodynamic shape update to the optimal solution. An optimization experiment is carried out based on a self-developed supersonic civil aircraft configuration with engines. The results show that the equivalent area distribution adjoint equation can accurately solve the gradient information. After optimization, the sonic boom level of the aircraft was reduced by 13.2 PLdB, and the drag coefficient was reduced by 60.75 counts. Moreover, the equivalent area distribution adjoint optimization method has outstanding advantages, such as high sensitivity and fast convergence speed, and can take both the low sonic boom and the low drag force of the aircraft into account, providing a powerful tool for the comprehensive optimization design of supersonic civil aircraft by considering sonic boom and aerodynamic force.
- Published
- 2024
- Full Text
- View/download PDF
15. Stochastic near-optimal control for a system with Markovian switching and Lévy noise.
- Author
-
Kuang, Daipeng, Li, Jianli, Gao, Dongdong, and Luo, Danfeng
- Subjects
- *
ERROR functions , *HAMILTON'S principle function , *STOCHASTIC systems , *VARIATIONAL principles , *STOCHASTIC control theory , *INTEGRAL functions - Abstract
The near-optimal control conditions of a stochastic system are the main subject of this research work, which will extend some previous results. At first, the stochastic system and its adjoint equations are estimated based on certain mild assumptions on a convex control set. Additionally, utilizing various analysis techniques and Ekeland's variational principle, sufficient and necessary near-optimality conditions are obtained, independent of the second-order adjoint equation. Ultimately, an example of a biological, mathematical model is used to demonstrate the correctness of the theoretical analysis. • Based on some standard assumptions, the necessary as well as sufficient conditions of near optimality of stochastic system are given, which may be used to find the development of stochastic optimal control. • The systems can be seen as a generalization of [8–10, 21](e.g., the system degenerates to the model (1) in [8] when g = 0,n = 1.), and the results in [9,10,21] are included in our paper. • Especially, one of main results shows that any ϵ –optimal control nearly minimizes the Hamiltonian function in an integral form with an error order of "almost" ϵ 3 16 . • The another result shows that ϵ –minimum condition in terms of the Hamiltonian function in the integral form is sufficient for the near- optimality of order ϵ 1 2 . [ABSTRACT FROM AUTHOR]
- Published
- 2024
- Full Text
- View/download PDF
16. OPTIMAL CONTROL OF THE LANDAU--DE GENNES MODEL OF NEMATIC LIQUID CRYSTALS.
- Author
-
SUROWIEC, THOMAS M. and WALKER, SHAWN W.
- Subjects
NEMATIC liquid crystals ,OPTIMIZATION algorithms ,NUMERICAL analysis ,LIQUID crystals ,FINITE element method - Abstract
We present an analysis and numerical study of an optimal control problem for the Landau--de Gennes (LdG) model of nematic liquid crystals (LCs), which is a crucial component in modern technology. They exhibit long range orientational order in their nematic phase, which is represented by a tensor-valued (spatial) order parameter Q=Q(x). Equilibrium LC states correspond to Q functions that (locally) minimize an LdG energy functional. Thus, we consider an L2-gradient flow of the LdG energy that allows for finding local minimizers and leads to a semilinear parabolic PDE, for which we develop an optimal control framework. We then derive several a priori estimates for the forward problem, including continuity in space-time, that allow us to prove existence of optimal boundary and external "force" controls and to derive optimality conditions through the use of an adjoint equation. Next, we present a simple finite element scheme for the LdG model and a straightforward optimization algorithm. We illustrate optimization of LC states through numerical experiments in two and three dimensions that seek to place LC defects (where Q(x) = 0) in desired locations, which is desirable in applications. [ABSTRACT FROM AUTHOR]
- Published
- 2023
- Full Text
- View/download PDF
17. Backward-Eulerian Footprint Modelling Based on the Adjoint Equation for Atmospheric and Urban-Terrain Dispersion
- Author
-
Jia, Hongyuan and Kikumoto, Hideki
- Published
- 2023
- Full Text
- View/download PDF
18. Estimation of the point source parameters by the adjoint equation in the time-varying atmospheric environment with unknown turn-on time.
- Author
-
Zhu, Jianjie, Zhou, Xuanyi, Cong, Beihua, and Kikumoto, Hideki
- Subjects
FIX-point estimation ,COMPUTATIONAL fluid dynamics ,BAYESIAN field theory ,EMERGENCY management - Abstract
Bayesian inference coupled with computational fluid dynamics (CFD) has been widely used in the source term estimation (STE). At present, most scholars have studied the pollution source released in the time-invariant flow field and the turn-on time of the source is often regarded as a known parameter which requires no estimation. Nevertheless, it is of greater practical significance to estimate the source parameters by considering the time-varying characteristics of flow field. Besides, estimating the turn-on time is very important as it provides the critical information for post-disaster rescue and accident cause analysis. In this paper, the predicted concentrations are calculated by the adjoint equation in the time-varying flow field, and the turn-on time, taken as an unknown parameter, is estimated together with the source location and the release rate. Accurate estimated results are thereby obtained. Then, how the sensors' sampling time intervals and the measured data in different dispersion stages influence the estimated results is investigated. It is found that reducing the sampling time interval can decrease the uncertainties of the estimations. In addition, in order to estimate the turn-on time accurately, the measured information in the developing stage of dispersion is indispensable. However, only using the measured information in the stable stage of dispersion cannot predict the turn-on time. Finally, the earliest time when the source parameters can be estimated with accuracy is also explored. The results show that the method of the STE proposed can obtain accurate source information at the early stage of dispersion. • The source parameters are estimated in the time-varying flow field. • The turn-on time treated as an unknown parameter is estimated. • The predicted concentrations are obtained by adjoint equation in time-varying flow. • The impacts of key factors of measured data on the estimated results are studied. • The source parameters can be estimated accurately at early stage of dispersion. [ABSTRACT FROM AUTHOR]
- Published
- 2023
- Full Text
- View/download PDF
19. Locating leakage in pipelines based on the adjoint equation of inversion modeling
- Author
-
Chang Chang, Xiangli Li, Lin Duanmu, Hongwei Li, and Wenbin Zhou
- Subjects
Pipeline leakage ,Inversion modeling ,Sensitivity analysis ,Adjoint equation ,Science (General) ,Q1-390 ,Social sciences (General) ,H1-99 - Abstract
This paper presents an adjoint method for locating potential leakage in a single-phase fluid pipeline based on the analytic solution of inversion modeling. By studying the mechanism of pipeline leakage pressure, the adjoint equation based on the governing equation of transient flow is established in the single-liquid phase aspect using inverse adjoint theory and sensitivity analysis method. The inverse transient adjoint equation is primarily derived from the single linear fluid pipeline in the semi-infinite domain. The Laplace method is then used to obtain an analytical solution that determines the location of pipeline leakage. The experimental results indicate that the analytic solution can quickly and accurately judge the leakage location of the pipeline. Furthermore, it presents a new approach to engineering applications, such as gas-liquid two-phase flow complex pipe networks, etc.
- Published
- 2023
- Full Text
- View/download PDF
20. The maximum principle for stochastic control problem with Markov chain in progressive structure.
- Author
-
Chen, Tian, Song, Yuanzhuo, and Wu, Zhen
- Subjects
- *
MARKOV processes , *STOCHASTIC control theory , *STOCHASTIC systems , *BROWNIAN motion - Abstract
In this paper, we obtain the stochastic maximum principle for progressive optimal control of the model driven by Brownian motions and Markov chain. The maximum principle in progressive structure is essentially different from that in the classical case. Our maximum principle has two parts, including continuous part and jump part, which are used to characterize the characteristics of control at continuous time and jump time respectively. This shows that the progressive structure can better describe the optimal control of stochastic system with jumps. Moreover, we give an LQ example to show that the cost functional can indeed get a smaller value in progressive structure. • Necessary and sufficient maximum principle in progressive structure. • Two maximum principle includes continuous part and jump part. • Progressive structure is essentially different from predictable structure. • Progressive structure accurately describes the feature of control for system with jump. • Progressive structure can find better control than predictable structure. [ABSTRACT FROM AUTHOR]
- Published
- 2022
- Full Text
- View/download PDF
21. Determination of a Time-Varying Point Source in Cauchy Problems for the Convection–Diffusion Equation
- Author
-
Slavi Georgiev and Lubin Vulkov
- Subjects
parabolic Cauchy problem ,boundary conditions ,dispersion modeling ,point source ,inverse problem ,adjoint equation ,Technology ,Engineering (General). Civil engineering (General) ,TA1-2040 ,Biology (General) ,QH301-705.5 ,Physics ,QC1-999 ,Chemistry ,QD1-999 - Abstract
In this paper, we suggest a method for recovering the unknown time-dependent strength of a contaminant concentration source from measurements of the concentration inside an unbounded domain. This problem is formulated as a Cauchy parabolic inverse problem. For its efficient numerical processing, the problem is solved by reduction of the Cauchy problem to a Dirichet one on a bounded domain using the method of the fundamental (potential) solutions in combination with an adjoint equation technique. A numerical solution to this approach is explained. Next, by choosing the source strength in the form of a finite series of shape functions with unknown constant coefficients and using a linear-square method, the term concentration source is estimated. Computational simulations using model examples from water pollution are discussed.
- Published
- 2023
- Full Text
- View/download PDF
22. Parallel Implementation of a Sensitivity Operator-Based Source Identification Algorithm for Distributed Memory Computers
- Author
-
Alexey Penenko and Evgeny Rusin
- Subjects
air quality ,large-scale inverse problem ,source identification ,sensitivity operator ,adjoint equation ,Lake Baikal region ,Mathematics ,QA1-939 - Abstract
Large-scale inverse problems that require high-performance computing arise in various fields, including regional air quality studies. The paper focuses on parallel solutions of an emission source identification problem for a 2D advection–diffusion–reaction model where the sources are identified by heterogeneous measurement data. In the inverse modeling approach we use, a source identification problem is transformed to a quasi-linear operator equation with a sensitivity operator, which allows working in a unified way with heterogeneous measurement data and provides natural parallelization of numeric algorithms by concurrent calculation of the rows of a sensitivity operator matrix. The parallel version of the algorithm implemented with a message passing interface (MPI) has shown a 40× speedup on four Intel Xeon Gold 6248R nodes in an inverse modeling scenario for the Lake Baikal region.
- Published
- 2022
- Full Text
- View/download PDF
23. Reduced-order methods for dynamic problems in topology optimization: A comparative study
- Author
-
Quhao Li, Jakob Søndergaard Jensen, Niels Aage, and Ole Sigmund
- Subjects
Optimization problem ,Computer science ,Mechanical Engineering ,Topology optimization ,Computational Mechanics ,Stability (learning theory) ,General Physics and Astronomy ,Krylov subspace ,Computer Science Applications ,Dynamic problem ,Mechanics of Materials ,Adjoint equation ,Applied mathematics ,Subspace topology ,Curse of dimensionality - Abstract
The dynamics of engineering structures are of great importance for topology optimization problems in both academia and industry. However, for design problems where broadband frequency responses are required, the computational burden becomes enormous, especially for large-scale applications. To remedy this numerical bottleneck, using the Reduced-Order Methods (ROMs) is an efficient approach by recasting the original problem into a subspace with a much smaller dimensionality than the full model. In this paper, a systematic comparative study of some typical and potential ROMs for solving the broadband frequency response optimization problems is provided, including the Quasi-Static Ritz Vector (QSRV), the Pade expansion and the second-order Krylov subspace method . Furthermore, the effects of the orthonormalization processes are discussed. Two representative test problems, a vibration problem and a wave propagation problem, are solved, analyzed, and compared based on the ROMs’ accuracy, their stability in approximating the state and adjoint equations and the applicability to topology optimization problems. From the extensive numerical results, we find that the second-order Krylov subspace with moment-matching Gram–Schmidt orthonormalization (SOMMG) and the Second-Order Arnoldi method (SOAR) provides superior accuracy and stability. Moreover, the results verify that the basis vectors computed for the state equation cannot be reused for solving the adjoint equation, and hence, that new basis vectors should be constructed. Analysis of the computational cost for the 3D test problems shows an improvement in numerical performance in the order of 100–10000 for the ROMs compared to the full approach.
- Published
- 2021
24. Development and verification of a higher-order mathematical adjoint nodal diffusion solver
- Author
-
Kostadin Ivanov, René van Geemert, Muhammad Ramzy Altahhan, and Maria Avramova
- Subjects
symbols.namesake ,Quadratic equation ,Nuclear Energy and Engineering ,Adjoint equation ,Lagrange multiplier ,symbols ,Applied mathematics ,Node (circuits) ,Sensitivity (control systems) ,Solver ,System of linear equations ,Eigenvalues and eigenvectors ,Mathematics - Abstract
In this paper, we derive a mathematical formulation of the higher order adjoint NEM-M2B2 equations by preconditioning the nodal interface neutron currents equations of the forward equations system, and by using the Lagrangian Multipliers analysis method. In the NEM-M2B2 system of equations, the quadratic transverse leakage approximation is used to model the leakage of neutrons between each node in the system. The solution of the adjoint equation can be used to perform adjoint-based predictive sensitivity/perturbation analysis. As an example, we use the mathematical adjoint solution as sensitivity weighting for predicting the response of the IAEA-3D benchmark’s eigenvalue to a perturbation in the independent parameters of the system (i.e., cross-sections). We also derive perturbation equations associated with the particular NEM-M2B2 model we are using. These perturbation-equations are used in predicting the model eigenvalue change without resorting to recalculating the forward NEM-M2B2 system of equations again (labeled as exact calculations). They also enabled construction of a reactivity sensitivity map showing the importance of each calculation node of the benchmark depending on its spatial and spectral coordinates. Perturbations were imposed on both the absorption cross-sections (fast and thermal) and the scattering cross-section of the IAEA-3D benchmark problem. Several verification steps were taken to ensure that the developed mathematical adjoint solver is adequate for adjoint analysis (e.g., commutativity checks, and comparison against exact calculations).
- Published
- 2021
25. An efficient adjoint method for the aero-stealth shape optimization design
- Author
-
Ming Li, Jiaojiao Chen, Junqiang Bai, Feng Qu, and Xiaoyu Feng
- Subjects
Radar cross-section ,Automatic differentiation ,Computer science ,Adjoint equation ,Finite difference ,Aerospace Engineering ,Shape optimization ,Sensitivity (control systems) ,Residual ,Reduction (mathematics) ,Algorithm - Abstract
Sensitivity analysis is an important part of gradient-based aero-stealth shape optimization design. An efficient Radar Cross Section (RCS) gradient evaluation method based on the adjoint approach is proposed in this paper. The merits of this proposed method mainly include the following two aspects. First, the adjoint field will be obtained from a single solution of the scattering problem without solving the adjoint equation. Second, the residual derivatives computing code is generated through the adjoint mode of Automatic Differentiation (AD) technique and is optimized by changing the “two-sweeps” architecture into the “inner-loop two-sweeps” architecture. The time consumption and memory usage are tested and the results show that the AD adjoint mode with “inner-loop two-sweeps” architecture is high-efficiency while the memory consumption is friendly. In addition, the accuracy of the adjoint method is verified through the finite differences method. After that, a gradient-based aero-stealth optimization design method is developed and employed to study the shape requirements of aerodynamic and stealth performance. The cost function is composed of the drag coefficient and the RCS. The results show that several trade-off configurations would be obtained by selecting different weight factors and the pure stealth optimization design would lead to a screwy shape. Further studies illustrate that the RCS reduction comes from two aspects. On the one hand, the sharp leading edge could reduce the surface current magnitude and avoid the specular back-scattering. On the other hand, the undulations could change the phases and weaken the superimposed effect.
- Published
- 2021
Discovery Service for Jio Institute Digital Library
For full access to our library's resources, please sign in.