• Structural Engineering and Mechanics
• /
• 제27권6호
• /
• pp.683-696
• /
• 2007

This article is concerned with the presentation of a time-domain BEM approach applied to the solution of the scalar wave equation for 2D problems. The basic idea is quite simple: the basic variables of the problem at time $t_n$ (potential and flux) are computed with the results related to the potential and to its time derivative at time $t_{n-1}$ playing the role of "initial conditions". This time-marching scheme needs the computation of the potential and its time derivative at all boundary nodes and internal points, as well as the entire discretization of the domain. The convolution integrals of the standard time-domain BEM formulation, however, are not computed; the matrices assembled, only at the initial time interval, are those related to the potential, flux and to the potential time derivative. Two examples are presented and discussed at the end of the article, in order to verify the accuracy and potentialities of the proposed formulation.

• Structural Engineering and Mechanics
• /
• 제23권3호
• /
• pp.263-278
• /
• 2006

This work presents a time-truncation scheme, based on the Lagrange interpolation polynomial, for the solution of the two-dimensional scalar wave problem by the time-domain boundary element method. The aim is to reduce the number of stored matrices, due to the convolution integral performed from the initial time to the current time, and to keep a compromise between computational economy and efficiency and the numerical accuracy. In order to verify the accuracy of the proposed formulation, three examples are presented and discussed at the end of the article.

• Current Optics and Photonics
• /
• 제5권5호
• /
• pp.491-499
• /
• 2021

A scalar Fourier modal method for the numerical analysis of the scalar wave equation in inhomogeneous space with an arbitrary permittivity profile, is proposed as a novel theoretical embodiment of Fourier optics. The modeling of devices and systems using conventional Fourier optics is based on the thin-element approximation, but this approach becomes less accurate with high numerical aperture or thick optical elements. The proposed scalar Fourier modal method describes the wave optical characteristics of optical structures in terms of the generalized transmittance function, which can readily overcome a current limitation of Fourier optics.

• 대한수학회논문집
• /
• 제16권3호
• /
• pp.509-518
• /
• 2001

The purpose of this paper is to propose a numerical algorithm for the problem of coefficient identification of the scalar wave equation based on a local observation on the boundary: Determine the unknown coefficient function with the knowledge of simultaneous Dirichlet and Neumann boundary values on a part of boundary. To find the unknown coefficient function, the unknown Neumann boundary value is also identified. We recast our inverse problem to variational problem. The gradient method is applied to find the minimizing functions. We confirm the effectiveness of our algorithm by numerical experiments.

• 전자공학회논문지A
• /
• 제32A권8호
• /
• pp.140-152
• /
• 1995

Derived was the scalar wave equation of optical fibers. Based on the derived equation, the dispersion characteristics of arbitrarily profiled fibers were analyzed. We applied the 1-D FEM employing quadratic interpolation fucntions to solve the scalar wave equation. To find the optimum index distribution of a fiber that has the ultra-low total dispersion, we analyzed QC fibers as objects. Adding 2$_{nd}$ and 3$_{rd}$ clads to DC fiber, we investigated the change of dispersion characteristics. We found the QC fiber parameters for which the dispersion was ultra-low flattened, less than 0.5 ps/km.nm for ${\lambda}=1.4~1.6{\mu}m$, and the dispersion value was as low as 0.20 ps/km.nm at ${\lambda}=1.55{\mu}m$.

• Structural Engineering and Mechanics
• /
• 제51권4호
• /
• pp.685-705
• /
• 2014

This paper presents a new formulation for forward scalar wave simulations in semi-infinite media. Perfectly-Matched-Layers (PMLs) are used as a wave absorbing boundary layer to surround a finite computational domain truncated from the semi-infinite domain. In this work, a hybrid formulation was developed for the simulation of scalar wave motion in two-dimensional PML-truncated domains. In this formulation, displacements and stresses are considered as unknowns in the PML domain, while only displacements are considered to be unknowns in the interior domain. This formulation reduces computational cost compared to fully-mixed formulations. To obtain governing wave equations in the PML region, complex coordinate stretching transformation was introduced to equilibrium, constitutive, and compatibility equations in the frequency domain. Then, equations were converted back to the time-domain using the inverse Fourier transform. The resulting equations are mixed (contain both displacements and stresses), and are coupled with the displacement-only equation in the regular domain. The Newmark method was used for the time integration of the semi-discrete equations.

• 대한수학회지
• /
• 제58권5호
• /
• pp.1131-1145
• /
• 2021

In this article, we study the Klein-Gordon-Maxwell equations arising from a semilocal gauge field model. This model describes the interaction of two complex scalar fields and one gauge field, and generalizes the classical Klein-Gordon equation coupled with the Maxwell electrodynamics. We prove that there exist infinitely many standing wave solutions for p ∈ (2, 6) which are radially symmetric. Here, p comes from the exponent of the potential of scalar fields. We also prove the nonexistence of nontrivial solutions for the critical case p = 6.

• 한국광학회지
• /
• 제6권2호
• /
• pp.156-164
• /
• 1995

광섬유의 스칼라 파동방정식을 토대로 임의의 굴절율 분포를 갖는 광섬유의 분산특성을 해석하였다. 스칼라 파동방정식을 해석하기 위해서 2차 보간함수를 사용한 1차원 유한요소법을 적용하였다. DC형 광섬유를 대상으로 하였으며, 총분산을 최소화할 수 있는 굴절율분포를 조사하였다. $\lambda=1.3, 1.55{\mu}m$두 파장에서 동시에 분산이 거의 0이 되는 설계치를 찾을 수 있었다. 그리고 $\lambda=1.4~1.7{\mu}m$에 걸쳐서 1.0ps/km.nm 이하의 평탄분산특성을 갖고 $\lambda=1.55{\mu}m$에서 0.65ps/km.nm의 분선을 갖는 설계치도 찾을 수 있었다.

• 자원환경지질
• /
• 제29권2호
• /
• pp.183-192
• /
• 1996

Much computing time and large computer memory are needed to solve the wave equation in a large complex subsurface layer using finite difference method. The time and memory can be reduced by decreasing the number of grid per minimun wave length. However, decrease of grid may cause numerical dispersion and poor accuracy. In this study, we present 49 points weighted average method which save the computing time and memory and improve the accuracy. This method applies a new weighted average to the coordinate determined by transforming the coordinate of conventional 5 points finite difference stars to $0^{\circ}$ and $45^{\circ}$, 25 points finite differenc stars to $0^{\circ}$, $26.56^{\circ}$, $45^{\circ}$, $63.44^{\circ}$ and 49 finite difference stars to $0^{\circ}$, $18.43^{\circ}$, $33.69^{\circ}$, $45^{\circ}$, $56.30^{\circ}$, $71.56^{\circ}$. By this method, the grid points per minimum wave length can be reduced to 2.5, the computing time to $(2.5/13)^3$, and the required core memory to $(2.5/13)^4$ computing with the conventional method.

• 한국전산구조공학회논문집
• /
• 제34권2호
• /
• pp.113-120
• /
• 2021

무한 매질에서의 파전파 현상은 공학과 자연과학의 여러 분야에서 다양한 물리적 현상을 서술하는데 활용되고 있고, 이 문제에 대한 해를 얻기 위하여 해석적 방법 또는 수치적 방법이 개발되어 사용되고 있다. 이 문제에 대한 정확한 해를 얻기 위해서는 무한 영역으로의 에너지 방사를 정확히 고려해야 하고, 이를 위해 다양한 수치적 또는 역학적 모형 또는 경계조건이 개발되었다. 이 연구에서는 층상 waveguide에서의 scalar wave 또는 SH파 전파 문제에 적용할 수 있는 새로운 경계조건을 제안하고자 한다. 이를 위해 waveguide의 수직방향으로 유한요소 이산화를 적용하여 얻은 SH파의 지배방정식을 변형하여 waveguide의 무한 영역의 영향을 나타내는 경계조건을 유도한다. 층상 waveguide에서의 SH파에 대한 고유모드의 직교성을 이용하여, 새로운 경계조건은 기존의 root-finding absorbing boundary condition와 동등함을 보이고, 이로부터 새로운 경계조건의 차수가 증가할수록 정확성이 증가하고, 또한 이산화된 수준에서도 안정함을 유도할 수 있다. 제안된 경계조건을 층상 waveguide에서의 파전파 문제에 적용하여 그 정확성과 안정성을 검증한다.