• 제목/요약/키워드: Adjoint Problem

검색결과 94건 처리시간 0.023초

SZEGÖ PROJECTIONS FOR HARDY SPACES IN QUATERNIONIC CLIFFORD ANALYSIS

  • He, Fuli;Huang, Song;Ku, Min
    • 대한수학회보
    • /
    • 제59권5호
    • /
    • pp.1215-1235
    • /
    • 2022
  • In this paper we study Szegö kernel projections for Hardy spaces in quaternionic Clifford analysis. At first we introduce the matrix Szegö projection operator for the Hardy space of quaternionic Hermitean monogenic functions by the characterization of the matrix Hilbert transform in the quaternionic Clifford analysis. Then we establish the Kerzman-Stein formula which closely connects the matrix Szegö projection operator with the Hardy projection operator onto the Hardy space, and we get the matrix Szegö projection operator in terms of the Hardy projection operator and its adjoint. At last, we construct the explicit matrix Szegö kernel function for the Hardy space on the sphere as an example, and get the solution to a Diriclet boundary value problem for matrix functions.

AN ASYMPTOTIC FINITE ELEMENT METHOD FOR SINGULARLY PERTURBED HIGHER ORDER ORDINARY DIFFERENTIAL EQUATIONS OF CONVECTION-DIFFUSION TYPE WITH DISCONTINUOUS SOURCE TERM

  • Babu, A. Ramesh;Ramanujam, N.
    • Journal of applied mathematics & informatics
    • /
    • 제26권5_6호
    • /
    • pp.1057-1069
    • /
    • 2008
  • We consider singularly perturbed Boundary Value Problems (BVPs) for third and fourth order Ordinary Differential Equations(ODEs) of convection-diffusion type with discontinuous source term and a small positive parameter multiplying the highest derivative. Because of the type of Boundary Conditions(BCs) imposed on these equations these problems can be transformed into weakly coupled systems. In this system, the first equation does not have the small parameter but the second contains it. In this paper a computational method named as 'An asymptotic finite element method' for solving these systems is presented. In this method we first find an zero order asymptotic approximation to the solution and then the system is decoupled by replacing the first component of the solution by this approximation in the second equation. Then the second equation is independently solved by a fitted mesh Finite Element Method (FEM). Numerical experiments support our theoritical results.

  • PDF

A novel sensitivity method to structural damage estimation in bridges with moving mass

  • Mirzaee, Akbar;Shayanfar, Mohsenali;Abbasnia, Reza
    • Structural Engineering and Mechanics
    • /
    • 제54권6호
    • /
    • pp.1217-1244
    • /
    • 2015
  • In this research a theoretical and numerical study on a bridge damage detection procedure is presented based on vibration measurements collected from a set of accelerometers. This method, referred to as "Adjoint Variable Method", is a sensitivity-based finite element model updating method. The approach relies on minimizing a penalty function, which usually consists of the errors between the measured quantities and the corresponding predictions attained from the model. Moving mass is an interactive model and includes inertia effects between the model and mass. This interactive model is a time varying system and the proposed method is capable of detecting damage in this variable system. Robustness of the proposed method is illustrated by correct 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 comparative study on common sensitivity and the proposed method confirms its efficiency and performance improvement in sensitivity-based damage detection methods. In addition various possible sources of error, including the effects of measurement noise and initial assumption error in stability of method are also discussed.

Analysis of alpha modes in multigroup diffusion

  • Sanchez, Richard;Tomatis, Daniele;Zmijarevic, Igor;Joo, Han Gyu
    • Nuclear Engineering and Technology
    • /
    • 제49권6호
    • /
    • pp.1259-1268
    • /
    • 2017
  • The alpha eigenvalue problem in multigroup neutron diffusion is studied with particular attention to the theoretical analysis of the model. Contrary to previous literature results, the existence of eigenvalue and eigenflux clustering is investigated here without the simplification of a unique fissile isotope or a single emission spectrum. A discussion about the negative decay constants of the neutron precursors concentrations as potential eigenvalues is provided. An in-hour equation is derived by a perturbation approach recurring to the steady state adjoint and direct eigenvalue problems of the effective multiplication factor and is used to suggest proper detection criteria of flux clustering. In spite of the prior work, the in-hour equation results give a necessary and sufficient condition for the existence of the eigenvalue-eigenvector pair. A simplified asymptotic analysis is used to predict bands of accumulation of eigenvalues close to the negative decay constants of the precursors concentrations. The resolution of the problem in one-dimensional heterogeneous problems shows numerical evidence of the predicted clustering occurrences and also confirms previous theoretical analysis and numerical results.

The methods of CADIS-NEE and CADIS-DXTRAN in NECP-MCX and their applications

  • Qingming He;Zhanpeng Huang;Liangzhi Cao;Hongchun Wu
    • Nuclear Engineering and Technology
    • /
    • 제56권7호
    • /
    • pp.2748-2755
    • /
    • 2024
  • This paper presents two new methods for variance reduction for shielding calculation in Monte Carlo radiation transport. One method is CADIS-NEE, which combines Consistent Adjoint Driven Importance Sampling (CADIS) and next-event estimator (NEE) methods to increase the calculation efficiency of tallies at points. The other is CADIS-deterministic transport (DXTRAN), which combines CADIS and DXTRAN to obtain higher performance than using CADIS and DXTRAN separately. The combination processes are derived and implemented in the hybrid Monte-Carlo-Deterministic particle-transport code NECP-MCX. Various problems are tested to demonstrate the effectiveness of the two methods. According to the results, the two combination methods have higher efficiency than using CADIS, NEE or DXTRAN separately. In a long-distance photon-transport problem, CADIS-NEE converges faster than NEE and the figure of merit (FOM) of CADIS-NEE is 75.6 times of NEE. In a labyrinthine problem, CADIS-DXTRAN's FOM surpasses that of DXTRAN and CADIS by a factor of 45.3 and 17.7, respectively. Therefore, it is advisable to employ these two novel methods selectively in appropriate scenarios to reduce variance.

구조물 최적화를 위한 선형화 기법 (Application of linearization method for large-scale structure optimizations)

  • 이희각
    • 전산구조공학
    • /
    • 제1권1호
    • /
    • pp.87-94
    • /
    • 1988
  • 반복 비선형 계획법의 하나인 선형화 기법을 절대수렴의 전제아래 합성 구조물의 최적 설계에 응용한다. 선형화 기법은 설계문제의 제약조건을 선형화된 등가 제약조건으로 변형시키며 active-set 정책을 구사한다. 결과, 매 설계단계에서 풀어야 할 상태 및 수반 방정식의 수를 줄임으로써 실질적인 계산의 절감을 기한다. 기둥으로 받쳐진 판-보 구조물은 최적화 기법의 능력을 시험키 위한 합성구조물의 좋은 예로서, 설계결과 선형화 기법은 만족할만한 수렴치로써 최적해를 산출함을 알 수 있고 나아가 이 방법은 모든 종류의 최적화 문제에 적용될 수 있을 것으로 보인다.

  • PDF

고 정밀 항공우주 유동해석 및 설계를 위한 공력계산 툴 (Essential Computational Tools for High-Fidelity Aerodynamic Simulation and Design)

  • 김종암
    • 유체기계공업학회:학술대회논문집
    • /
    • 유체기계공업학회 2006년 제4회 한국유체공학학술대회 논문집
    • /
    • pp.33-36
    • /
    • 2006
  • As the computing environment is rapidly improved, the interests of CFD are gradually focused on large-scale computation over complex geometry. Keeping pace with the trend, essential computational tools to obtain solutions of complex aerospace flow analysis and design problems are examined. An accurate and efficient flow analysis and design codes for large-scale aerospace problem are presented in this work. With regard to original numerical schemes for flow analysis, high-fidelity flux schemes such as RoeM, AUSMPW+ and higher order interpolation schemes such as MLP (Multi-dimensional Limiting Process) are presented. Concerning the grid representation method, a general-purpose basis code which can handle multi-block system and overset grid system simultaneously is constructed. In respect to design optimization, the importance of turbulent sensitivity is investigated. And design tools to predict highly turbulent flows and its sensitivity accurately by fully differentiating turbulent transport equations are presented. Especially, a new sensitivity analysis treatment and geometric representation method to resolve the basic flow characteristics are presented. Exploiting these tools, the capability of the proposed approach to handle complex aerospace simulation and design problems is tested by computing several flow analysis and design problems.

  • PDF

Goal-oriented multi-collision source algorithm for discrete ordinates transport calculation

  • Wang, Xinyu;Zhang, Bin;Chen, Yixue
    • Nuclear Engineering and Technology
    • /
    • 제54권7호
    • /
    • pp.2625-2634
    • /
    • 2022
  • Discretization errors are extremely challenging conundrums of discrete ordinates calculations for radiation transport problems with void regions. In previous work, we have presented a multi-collision source method (MCS) to overcome discretization errors, but the efficiency needs to be improved. This paper proposes a goal-oriented algorithm for the MCS method to adaptively determine the partitioning of the geometry and dynamically change the angular quadrature in remaining iterations. The importance factor based on the adjoint transport calculation obtains the response function to get a problem-dependent, goal-oriented spatial decomposition. The difference in the scalar fluxes from one high-order quadrature set to a lower one provides the error estimation as a driving force behind the dynamic quadrature. The goal-oriented algorithm allows optimizing by using ray-tracing technology or high-order quadrature sets in the first few iterations and arranging the integration order of the remaining iterations from high to low. The algorithm has been implemented in the 3D transport code ARES and was tested on the Kobayashi benchmarks. The numerical results show a reduction in computation time on these problems for the same desired level of accuracy as compared to the standard ARES code, and it has clear advantages over the traditional MCS method in solving radiation transport problems with reflective boundary conditions.

Second order of average current nodal expansion method for the neutron noise simulation

  • Poursalehi, N.;Abed, A.
    • Nuclear Engineering and Technology
    • /
    • 제53권5호
    • /
    • pp.1391-1402
    • /
    • 2021
  • The aim of this work is to prepare a neutron noise calculator based on the second order of average current nodal expansion method (ACNEM). Generally, nodal methods have the ability to fulfill the neutronic analysis with adequate precision using coarse meshes as large as a fuel assembly size. But, for the zeroth order of ACNEM, the accuracy of neutronic simulations may not be sufficient when coarse meshes are employed in the reactor core modeling. In this work, the capability of second order ACNEM is extended for solving the neutron diffusion equation in the frequency domain using coarse meshes. For this purpose, two problems are modeled and checked including a slab reactor and 2D BIBLIS PWR. For validating of results, a semi-analytical solution is utilized for 1D test case, and for 2D problem, the results of both forward and adjoint neutron noise calculations are exploited. Numerical results indicate that by increasing the order of method, the errors of frequency dependent coarse mesh solutions are considerably decreased in comparison to the reference. Accordingly, the accuracy of second order ACNEM can be acceptable for the neutron noise calculations by using coarse meshes in the nuclear reactor core.

무요소법을 이용한 균열진전 문제의 형상 최적설계 (Shape Design Optimization of Crack Propagation Problems Using Meshfree Methods)

  • 김재현;하승현;조선호
    • 한국전산구조공학회논문집
    • /
    • 제27권5호
    • /
    • pp.337-343
    • /
    • 2014
  • 본 논문에서는 재생 커널 기법을 사용하여 혼합모드 균열진전 문제에 대한 연속체 기반의 형상 설계민감도 해석을 수행하였다. 재생 커널 기법은 기존의 유한요소법과 달리 요소망을 재구성할 필요가 없어, 커널 함수의 연속성을 증가시켰을 때 높은 정밀도의 형상함수를 얻을 수 있다는 장점을 가지고 있다. 균열선단 주변에서 J-적분을 수행하기 위해 선형탄성 조건이 고려되었다. 변위장과 응력 확대 계수의 설계변수에 대한 감도해석을 위하여 물질도함수를 도입하였으며 직접 미분법보다 효율적인 애조인 방법을 사용하여 설계민감도를 유도하였다. 수치 예제들을 통해서 재생 커널 기법을 이용한 균열진전 해석결과의 타당성을 확인하였으며 애조인 방법을 이용한 형상 설계민감도 해석 결과를 유한차분법과 비교하여 매우 정확하고 효율적인 결과를 얻을 수 있음을 알 수 있었다. 이를 바탕으로 간단한 모델에 대하여 형상 최적설계를 수행하여 균열이 발생될 수 있는 구조물에 대해서 균열에 의한 피해를 최소화할 수 있도록 균열을 제어할 수 있는 최적의 형상을 도출하였다.