• Title/Summary/Keyword: Adjoint

Search Result 355, Processing Time 0.022 seconds

Singular Representation and Finite Element Methods

  • 김석찬
    • Proceedings of the Korean Society of Computational and Applied Mathematics Conference
    • /
    • 2003.09a
    • /
    • pp.9-9
    • /
    • 2003
  • Let $\Omega$ be a bounded, open, and polygonal domain in $R^2$ with re-entrant corners. We consider the following Partial Differential Equations: $$(I-\nabla\nabla\cdot+\nabla^{\bot}\nabla\times)u\;=\;f\;in\;\Omega$$, $$n\cdotu\;0\;0\;on\;{\Gamma}_{N}$$, $${\nabla}{\times}u\;=\;0\;on\;{\Gamma}_{N}$$, $$\tau{\cdot}u\;=\;0\;on\;{\Gamma}_{D}$$, $$\nabla{\cdot}u\;=\;0\;on\;{\Gamma}_{D}$$ where the symbol $\nabla\cdot$ and $\nabla$ stand for the divergence and gradient operators, respectively; $f{\in}L^2(\Omega)^2$ is a given vector function, $\partial\Omega=\Gamma_{D}\cup\Gamma_{N}$ is the partition of the boundary of $\Omega$; nis the outward unit vector normal to the boundary and $\tau$represents the unit vector tangent to the boundary oriented counterclockwise. For simplicity, assume that both $\Gamma_{D}$ and $\Gamma_{N}$ are nonempty. Denote the curl operator in $R^2$ by $$\nabla\times\;=\;(-{\partial}_2,{\partial}_1$$ and its formal adjoint by $${\nabla}^{\bot}\;=\;({-{\partial}_1}^{{\partial}_2}$$ Consider a weak formulation(WF): Find $u\;\in\;V$ such that $$a(u,v):=(u,v)+(\nabla{\cdot}u,\nabla{\cdot}v)+(\nabla{\times}u,\nabla{\times}V)=(f,v),\;A\;v{\in}V$$. (2) We assume there is only one singular corner. There are many methods to deal with the domain singularities. We introduce them shortly and we suggest a new Finite Element Methods by using Singular representation for the solution.

  • PDF

Sensitivity analysis for optimal design of piezoelectric structures (압전지능구조물의 최적설계를 위한 민감도 해석)

  • 김재환
    • Journal of KSNVE
    • /
    • v.8 no.2
    • /
    • pp.267-273
    • /
    • 1998
  • This study aims at performing sensitivity analysis of piezoelectric smart structure for minimizing radiated noise from the structure, The structure consists of a flat plate on which disk shaped piezoelectric actuator is mounted, and finite element modeling is used for the structure. The finite element modeling uses a combination of three dimensional piezoelectric, flat shell and transition elements so thus it can take into account the coupling effects of the piezoelectric device precisely and it can also reduce the degrees of freedom of the finite element model. Electric potential on the piezoelectric actuator is taken as a design variable and total radiated power of the structure is chosen as an objective function. The objective function can be represented as Rayleigh's integral equation and is a function of normal displacements of the structure. For the convenience of computation, all degrees of freedom of the finite element equation is condensed out except the normal displacements of the structure. To perform the design sensitivity analysis, the derivative of the objective function with respect to the normal displacements is found, and the derivative of the norma displacements with respect to the design variable is calculated from the finite element equation by using so called the adjoint variable method. The analysis results are compared with those of the finite difference method, and shows a good agreement. This sensitivity analysis is faster and more accurate than the finite difference method. Once the sensitivity analysis program is used for gradient-based optimizations, one could achieve a better convergence rate than non-derivative methods for optimal design of piezoelectric smart structures.

  • PDF

CONTINUITY OF HOMOMORPHISMS BETWEEN BANACH ALGEBRAS

  • Cho, Tae-Geun
    • Bulletin of the Korean Mathematical Society
    • /
    • v.20 no.2
    • /
    • pp.71-74
    • /
    • 1983
  • The problems of the continuity of homomorphisms between Banach algebras have been studied widely for the last two decades to obtain various fruitful results, yet it is far from characterizing the calss of Banach algebras for which each homomorphism from a member of the class into a Banach algebra is conitnuous. For commutative Banach algebras A and B a simple proof shows that every homomorphism .theta. from A into B is continuous provided that B is semi-simple, however, with a non semi-simple Banach algebra B examples of discontinuous homomorphisms from C(K) into B have been constructed by Dales [6] and Esterle [7]. For non commutative Banach algebras the problems of automatic continuity of homomorphisms seem to be much more difficult. Many positive results and open questions related to this subject may be found in [1], [3], [5] and [8], in particular most recent development can be found in the Lecture Note which contains [1]. It is well-known that a$^{*}$-isomorphism from a $C^{*}$-algebra into another $C^{*}$-algebra is an isometry, and an isomorphism of a Banach algebra into a $C^{*}$-algebra with self-adjoint range is continuous. But a$^{*}$-isomorphism from a $C^{*}$-algebra into an involutive Banach algebra is norm increasing [9], and one can not expect each of such isomorphisms to be continuous. In this note we discuss an isomorphism from a commutative $C^{*}$-algebra into a commutative Banach algebra with dense range via separating space. It is shown that such an isomorphism .theta. : A.rarw.B is conitnuous and maps A onto B is B is semi-simple, discontinuous if B is not semi-simple.

  • PDF

ON OPERATORS SATISFYING Tm(T|T|2kT)1/(k+1)Tm ≥ Tm|T|2Tm

  • Rashid, Mohammad H.M.
    • Communications of the Korean Mathematical Society
    • /
    • v.32 no.3
    • /
    • pp.661-676
    • /
    • 2017
  • Let T be a bounded linear operator acting on a complex Hilbert space ${\mathfrak{H}}$. In this paper we introduce the class, denoted ${\mathcal{Q}}(A(k),m)$, of operators satisfying $T^{m{\ast}}(T^{\ast}{\mid}T{\mid}^{2k}T)^{1/(k+1)}T^m{\geq}T^{{\ast}m}{\mid}T{\mid}^2T^m$, where m is a positive integer and k is a positive real number and we prove basic structural properties of these operators. Using these results, we prove that if P is the Riesz idempotent for isolated point ${\lambda}$ of the spectrum of $T{\in}{\mathcal{Q}}(A(k),m)$, then P is self-adjoint, and we give a necessary and sufficient condition for $T{\otimes}S$ to be in ${\mathcal{Q}}(A(k),m)$ when T and S are both non-zero operators. Moreover, we characterize the quasinilpotent part $H_0(T-{\lambda})$ of class A(k) operator.

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

  • Kim, Chong-Am
    • 유체기계공업학회:학술대회논문집
    • /
    • 2006.08a
    • /
    • 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

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
    • /
    • v.26 no.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

An Overview of Optimization of Structures Subjected to Transient Loads (동하중을 받는 구조물의 최적화에 관한 연구동향)

  • Park Gyung-Jin;Kang Byung-Soo
    • Transactions of the Korean Society of Mechanical Engineers A
    • /
    • v.29 no.3 s.234
    • /
    • pp.369-386
    • /
    • 2005
  • Various aspects of structural optimization techniques under transient loads are extensively reviewed. The main themes of the paper are treatment of time dependent constraints, calculation of design sensitivity, and approximation. Each subject is reviewed with the corresponding papers that have been published since 1970s. The treatment of time dependent constraints in both the direct method and the transformation method is discussed. Two ways of calculating design sensitivity of a structure under transient loads are discussed - direct differentiation method and adjoint variable method. The approximation concept mainly focuses on re- sponse surface method in crashworthiness and local approximation with the intermediate variable Especially, as an approximated optimization technique, Equivalent Static Load method which takes advantage of the well-established static response optimization technique is introduced. And as an application area of dynamic response optimization technique, the structural optimization in flexible multibody dynamic systems is re- viewed in the viewpoint of the above three themes

Design Sensitivity Analysis for the Vibration Characteristics of Vehicle Structure (수송체 구조물의 진동특성에 관한 설계민감도 해석)

  • 이재환
    • Computational Structural Engineering
    • /
    • v.7 no.1
    • /
    • pp.91-98
    • /
    • 1994
  • Design sensitivity analysis method for the vibration of vehicle structure is developed using adjoint variable method. A variational approach with complex response method is used to derive sensitivity expression. To evaluate sensitivity, FEM analysis of ship deck and vehicle structure are performed using MSC/NASTRAN installed in the super computer CRAY2S, and sensitivity computation is performed by PC. The accuracy of sensitivity is verified by the results of finite difference method. When compared to structural analysis time on CRAY2S, sensitivity computation is remarkably economical. The sensitivity of vehicle frame can be used to reduce the vibration responses such as displacement and acceleration of vehicle.

  • PDF

Level Set Based Shape Optimization of Linear Structures Using Topological Derivatives (Topological Derivative를 이용한 선형 구조물의 레벨셋 기반 형상 최적 설계)

  • Ha Seung-Hyun;Kim Min-Geun;Cho Seon-Ho
    • Proceedings of the Computational Structural Engineering Institute Conference
    • /
    • 2006.04a
    • /
    • pp.299-306
    • /
    • 2006
  • Using a level set method and topological derivatives, a topological shape optimization method that is independent of an initial design is developed for linearly elastic structures. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set function, which facilitates to handle complicated topological shape changes. The 'Hamilton-Jacobi (H-J)' equation and computationally robust numerical technique of 'up-wind scheme' lead the initial implicit boundary to an optimal one according to the normal velocity field while minimizing the objective function of compliance and satisfying the constraint of allowable volume. Based on the asymptotic regularization concept, the topological derivative is considered as the limit of shape derivative as the radius of hole approaches to zero. The required velocity field to update the H -J equation is determined from the descent direction of Lagrangian derived from optimality conditions. It turns out that the initial holes is not required to get the optimal result since the developed method can create holes whenever and wherever necessary using indicators obtained from the topological derivatives. It is demonstrated that the proper choice of control parameters for nucleation is crucial for efficient optimization process.

  • PDF

An Inverse Analysis of Two-Dimensional Heat Conduction Problem Using Regular and Modified Conjugate Gradient Method (표준공액구배법과 수정공액구배법을 이용한 2차원 열전도 문제의 역해석)

  • Choi, Eui-Rak;Kim, Woo-Seung
    • Transactions of the Korean Society of Mechanical Engineers B
    • /
    • v.22 no.12
    • /
    • pp.1715-1725
    • /
    • 1998
  • A two-dimensional transient inverse heat conduction problem involving the estimation of the unknown location, ($X^*$, $Y^*$), and timewise varying unknown strength, $G({\tau})$, of a line heat source embedded inside a rectangular bar with insulated boundaries has been solved simultaneously. The regular conjugate gradient method, RCGM and the modified conjugate gradient method, MCGM with adjoint equation, are used alternately to estimate the unknown strength $G({\tau})$ of the source term, while the parameter estimation approach is used to estimate the unknown location ($X^*$, $Y^*$) of the line heat source. The alternate use of the regular and the modified conjugate gradient methods alleviates the convergence difficulties encountered at the initial and final times (i.e ${\tau}=0$ and ${\tau}={\tau}_f$), hence stabilizes the computation and fastens the convergence of the solution. In order to examine the effectiveness of this approach under severe test conditions, the unknown strength $G({\tau})$ is chosen in the form of rectangular, triangular and sinusoidal functions.