• Title/Summary/Keyword: adjoint model

Search Result 53, Processing Time 0.026 seconds

Estimation of Optimal Weight in Tidal Modeling with the Adjoint Method (조석 모델링에서 adjoint 방법 적용시 적정 가중치 산정)

  • Lee, Jae-Hak;Park, Kyeong;Song, Yong-Sik
    • The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
    • /
    • v.5 no.3
    • /
    • pp.177-185
    • /
    • 2000
  • The adjoint method is a method of data assimilation to improve the model results by seeking for model parameters that minimize the cost function and satisfy the governing equations of a model simultaneously. An adjoint package was set up for the two-dimensional linear tidal model and was applied to an idealized domain for an optimal estimation of the open boundary conditions. The assimilating data were selected from the results of forward modeling. Attention is paid on the response of the adjoint package to weighting parameters, the importance of initial estimates of model parameters and the applicability of the adjoint package to the case with varying depth. A procedure to determine optimal weight is presented based on the relationships between weights and other factors.

  • PDF

AERODYNAMIC SENSITIVITY ANALYSIS FOR NAVIER-STOKES EQUATIONS

  • Kim, Hyoung-Jin;Kim, Chongam;Rho, Oh-Hyun;Lee, Ki Dong
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • v.3 no.2
    • /
    • pp.161-171
    • /
    • 1999
  • Aerodynamic sensitivity analysis codes are developed via the hand-differentiation using a direct differentiation method and an adjoint method respectively from discrete two-dimensional compressible Navier-Stokes equations. Unlike previous other researches, Baldwin-Lomax algebraic turbulence model is also differentiated by hand to obtain design sensitivities with respect to design variables of interest in turbulent flows. Discrete direct sensitivity equations and adjoint equations are efficiently solved by the same time integration scheme adopted in the flow solver routine. The required memory for the adjoint sensitivity code is greatly reduced at the cost of the computational time by allowing the large banded flux jacobian matrix unassembled. Direct sensitivity code results are found to be exactly coincident with sensitivity derivatives obtained by the finite difference. Adjoint code results of a turbulent flow case show slight deviations from the exact results due to the limitation of the algebraic turbulence model in implementing the adjoint formulation. However, current adjoint sensitivity code yields much more accurate sensitivity derivatives than the adjoint code with the turbulence eddy viscosity being kept constant, which is a usual assumption for the prior researches.

  • PDF

Development of Tools for calculating Forecast Sensitivities to the Initial Condition in the Korea Meteorological Administration (KMA) Unified Model (UM) (통합모델의 초기 자료에 대한 예측 민감도 산출 도구 개발)

  • Kim, Sung-Min;Kim, Hyun Mee;Joo, Sang-Won;Shin, Hyun-Cheol;Won, DukJin
    • Atmosphere
    • /
    • v.21 no.2
    • /
    • pp.163-172
    • /
    • 2011
  • Numerical forecasting depends on the initial condition error strongly because numerical model is a chaotic system. To calculate the sensitivity of some forecast aspects to the initial condition in the Korea Meteorological Administration (KMA) Unified Model (UM) which is originated from United Kingdom (UK) Meteorological Office (MO), an algorithm to calculate adjoint sensitivities is developed by modifying the adjoint perturbation forecast model in the KMA UM. Then the new algorithm is used to calculate adjoint sensitivity distributions for typhoon DIANMU (201004). Major initial adjoint sensitivities calculated for the 48 h forecast error are located horizontally in the rear right quadrant relative to the typhoon motion, which is related with the inflow regions of the environmental flow into the typhoon, similar to the sensitive structures in the previous studies. Because of the upward wave energy propagation, the major sensitivities at the initial time located in the low to mid- troposphere propagate upward to the upper troposphere where the maximum of the forecast error is located. The kinetic energy is dominant for both the initial adjoint sensitivity and forecast error of the typhoon DIANMU. The horizontal and vertical energy distributions of the adjoint sensitivity for the typhoon DIANMU are consistent with those for other typhoons using other models, indicating that the tools for calculating the adjoint sensitivity in the KMA UM is credible.

Variational Data Assimilation for Optimal Initial Conditions in Air Quality Modeling

  • Park, Seon-Ki
    • Journal of Korean Society for Atmospheric Environment
    • /
    • v.19 no.E2
    • /
    • pp.75-81
    • /
    • 2003
  • Variational data assimilation, which is recently introduced to the air quality modeling, is a promising tool for obtaining optimal estimates of initial conditions and other important parameters such as emission and deposition rates. In this paper. two advanced techniques for variational data assimilation, based on the adjoint and quasi-inverse methods, are tested for a simple air quality problem. The four-dimensional variational assimilation (4D-Var) requires to run an adjoint model to provide the gradient information in an iterative minimization process, whereas the inverse 3D-Var (I3D-Var) seeks for optimal initial conditions directly by running a quasi -inverse model. For a process with small dissipation, I3D-Vu outperforms 4D-Var in both computing time and accuracy. Hybrid application which combines I3D-Var and standard 4D-Var is also suggested for efficient data assimilation in air quality problems.

An Affordable Implementation of Kalman Filter by Eliminating the Explicit Temporal Evolution of the Background Error Covariance Matrix (칼만필터의 자료동화 활용을 위한 배경오차 공분산의 명시적 시간 진전 제거)

  • Lim, Gyu-Ho;Suh, Ae-Sook;Ha, Ji-Hyun
    • Atmosphere
    • /
    • v.23 no.1
    • /
    • pp.33-37
    • /
    • 2013
  • In meteorology, exploitation of Kalman filter as a data assimilation system is virtually impossible due to simultaneous requirements of adjoint model and large computer resource. The other substitute of utilizing ensemble Kalman filter is only affordable by compensating an enormous usage of computing resource. Furthermore, the latter employs ensemble integration sets for evolving the background error covariance matrix by compensating the dynamical feature of the temporal evolution of weather conditions. We propose a new implementation method that works without the adjoint model by utilizing the explicit expression of the background error covariance matrix in backward evolution. It will also break a barrier in the evolution of the covariance matrix. The method may be applied with a slight modification to the real time assimilation or the retrospective analysis.

A Study on the Optimal Position for the Secondary Neutron Source in Pressurized Water Reactors

  • Sun, Jungwon;Yahya, Mohd-Syukri;Kim, Yonghee
    • Nuclear Engineering and Technology
    • /
    • v.48 no.6
    • /
    • pp.1291-1302
    • /
    • 2016
  • This paper presents a new and efficient scheme to determine the optimal neutron source position in a model near-equilibrium pressurized water reactor, which is based on the OPR1000 Hanul Unit 3 Cycle 7 configuration. The proposed scheme particularly assigns importance of source positions according to the local adjoint flux distribution. In this research, detailed pin-by-pin reactor adjoint fluxes are determined by using the Monte Carlo KENO-VI code from solutions of the reactor homogeneous critical adjoint transport equations. The adjoint fluxes at each allowable source position are subsequently ranked to yield four candidate positions with the four highest adjoint fluxes. The study next simulates ex-core detector responses using the Monte Carlo MAVRIC code by assuming a neutron source is installed in one of the four candidate positions. The calculation is repeated for all positions. These detector responses are later converted into an inverse count rate ratio curve for each candidate source position. The study confirms that the optimal source position is the one with very high adjoint fluxes and detector responses, which is interestingly the original source position in the OPR1000 core, as it yields an inverse count rate ratio curve closest to the traditional 1/M line. The current work also clearly demonstrates that the proposed adjoint flux-based approach can be used to efficiently determine the optimal geometry for a neutron source and a detector in a modern pressurized water reactor core.

ADJOINT METHOD FOR CONTROLLED CAVITATION INVERSE NOZZLE DESIGN

  • Petropoulou, S.;Gavaises, M.;Theodorakakos, A.
    • International Journal of Automotive Technology
    • /
    • v.7 no.3
    • /
    • pp.283-288
    • /
    • 2006
  • A mathematical methodology is proposed for designing nozzle hole shapes producing controlled geometric cavitation. The proposed methodology uses an unstructured RANS flow solver, with the ability to compute sensitivity derivatives via an adjoint algorithm. The adjoint formulation for the N-S equations is presented while variation of the turbulence viscosity is not taken into account during the geometry modifications. The sensitivities are calculated in a mode independently of the shape parameterisation. The method is used to develop and evaluate conceptual shapes for nozzle hole cavitation reduction. The localized region at the hole inlet producing cavitation, is parameterised using its radius of curvature, while a cost function is formulated to eliminate the negative pressures present at this location. Sensitivity derivatives are used to assess the dependence of the localized region on the minimum pressure, and to drive the geometry to the targeted shape. The results show that the computer model can provide nozzle hole entry shapes that produce predefined flow characteristics, and thus can be used as an inverse design tool for nozzle hole cavitation control.

Simultaneous identification of damage in bridge under moving mass by Adjoint variable method

  • Mirzaee, Akbar;Abbasnia, Reza;Shayanfar, Mohsenali
    • Smart Structures and Systems
    • /
    • v.21 no.4
    • /
    • pp.449-467
    • /
    • 2018
  • In this paper, a theoretical and numerical study on bridge simultaneous damage detection procedure for identifying both the system parameters and input excitation mass, are presented. This method is called 'Adjoint Variable Method' which is an iterative gradient-based model updating method based on the dynamic response sensitivity. The main advantage of proposed method is inclusion of an analytical method to augment the accuracy and speed of the solution. Moving mass is a model which takes into account the inertia effects of the vehicle. This interaction model is a time varying system and proposed method is capable of detecting damage in this variable system. Robustness of proposed method is illustrated by correctly detection of the location and extension of predetermined single, multiple and random damages in all ranges of speed and mass ratio of moving vehicle. A comparison study of common sensitivity and proposed method confirms its efficiency and performance improvement in sensitivity-based damage detection methods. Various sources of errors including the effects of measurement noise and initial assumption error in stability of method are also discussed.

Robust Design Optimization of a Fighter Wing Using an Uncertainty Model Constructed by Neural Network (신경망으로 구축된 불확실성 모델을 이용한 전투기 날개의 강건 최적 설계)

  • Kim, Ju-Hyun;Kim, Byung-Kon;Jun, Sang-Ook;Jeon, Yong-Hee;Lee, Dong-Ho
    • Journal of the Korean Society for Aeronautical & Space Sciences
    • /
    • v.36 no.2
    • /
    • pp.99-104
    • /
    • 2008
  • This study performed robust design optimization of fighter wing planform, considering uncertainty based on neural network model. To construct uncertainty model, aerodynamic performance and their sensitivity were evaluated by 3-dimensional Euler equations and adjoint variable method at experimental points selected from central composite design. In addition, because a neural network model has the advantage of capturing non-linear characteristic, it was possible to predict sensitivity of the aerodynamic performance efficiently and accurately . From the results of robust design optimization, it could be confirmed that the robustness of the objective function and constraints were improved if the variation of uncertainty and sigma level were increased.