• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.1
• /
• pp.75-81
• /
• 2002

In the re-entry comtrol system, errors apt to induce because the time deriviative of drag acceleration is analytically estimated. Still more, the difficulty of estimation of th exact drag coefficient in hypersonic velocity and the non-reality of the scale height cause a steady-state drag errer. In the Space-Shuttle, a steady-state drag error is reduced by the addition of the integral term of drag acceleation error into the control system. This method, however, induces a difficulties in respect to the modern controller composition due to the multi-poles in a closed-loop system. Thus, this paper proposes the additional method of the disturbance observer. This reduces the steady-state drag error according to the following by the analytic calculation, and then creates the new drag acceleration time derivative using the estimated error. The performance of the re-entry control system is verified about 32 refernce trajectories.

• Nuclear Engineering and Technology
• /
• v.23 no.4
• /
• pp.387-400
• /
• 1991

An ocean compartment model simulating transport of nuclides by advection due to ocean circulation and intertaction with suspended sediments is developed, by which concentration breakthrough curves of nuclides can be calculated as a function of time. Dividing ocean into arbitrary number of characteristic compartments and performing a balance of mass of nuclides in each ocean compartment, the governing equation for the concentration in the ocean is obtained and a solution by the numerical integration is obtained. The integration method is specially useful for general stiff systems. For transfer coefficients describing advective transport between adjacent compartments by ocean circulation, the ocean turnover time is calculated by a two-dimensional numerical ocean model. To exemplify the compartment model, a reference case calculation for breakthrough curves of three nuclides in low-level radioactive wastes, Tc-99, Cs-137, and Pu-238 released from hypothetical repository under the seabed is carried out with five ocean compartments. Sensitivity analysis studies for some parameters to the concentration breakthrough curves are also made, which indicates that parameters such as ocean turnover time and ocean water volume of compartments have an important effect on the breakthrough curves.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.18 no.9
• /
• pp.1023-1029
• /
• 2007

In this paper, we propose the method to solve the BLT equations using Laguerre polynomials in time domain. The solution of BLT equations is obtained by recursive, differential and integral properties of Laguerre polynomials. The verification of the proposed method is tested by applying it to the two-wired transmission line with resistors and capacitors, which is illuminated by the electromagnetic plane wave pulse. And the result is compared with the corresponding transient responses obtained from inverse fast Fourier transform(IFFT) of the frequency domain solutions of BLT equations.

• Journal of the Korea Society for Simulation
• /
• v.25 no.3
• /
• pp.15-28
• /
• 2016

It is necessary to analyze underwater acoustics channel(UAC) modeling and simulation for underwater weapon system development and acquisition. In order to analyze UAC, there are underwater acoustics propagation numerical analysis models(Ray theory, Parabolic equation, Normal-mode, Wavenumber integration). However, If these models are used for multiple frequency signal analysis, they are inaccurate to calculate result of analysis effectiveness and restricted for signal processing and analysis. In this paper, to overcome this problem, we propose simple/multiple frequency signal analysis model of the Pseudospectral Time-Domain Method synthetic environment UAC applied to underwater environment noise model as like as realistic underwater environment. In order to confirm the validation of the model, we performed the 9 scenarios simulation(4 scenarios of single frequency signal, 4 scenarios of multiple frequency signal, 1 scenario of single/multiple frequency signal like submarine radiated noise) for validation and confirmed the validation of this model through the simulation model.

• Journal of the Korean Magnetics Society
• /
• v.24 no.1
• /
• pp.18-21
• /
• 2014

The microwave absorption ($P_{tot}$), which is the double integration value of ferromagnetic resonance signal, propositional to the saturation magnetization, and the increase of the $P_{tot}$ measured during the thermal reaction time expect the growth process of the nanoparticles. Therefore, in this work, we measured the $P_{tot}$ in order to obtain the growth time of iron oxide nanoparticles after thermal decomposition of $Fe(acac)_3$ precursor at aging temperature $T_a=273$, 300 and $324^{\circ}C$, respectively. The best condition for monodisperse nanoparticles was obtained at $T_a=300^{\circ}C$, which condition showed the most rapid increase of $P_{tot}$ with thermal reaction time. Finally, the rapid growth rate was necessary condition for the synthesis of iron-oxide monodisperse nanoparticles.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.6
• /
• pp.331-337
• /
• 2018

This paper presents dynamic algorithm of the MLS(moving least squares) difference method using first order differential Approximation. The governing equations are only discretized by the first order MLS derivative approximation. The system equation consists of an assembly of the approximate function, so the shape of system equation is similar to FEM(finite element method). The CDM(central difference method) is used for time integration of dynamic equilibrium equation. The natural frequency analyses of the MLS difference method and FEM are performed, and two analysis results are compared. Also, the accuracy of the proposed numerical method is verified by displaying the dynamic analysis results together with the results by the existing second order differential approximation. In the process of assembling the first order MLS derivative approximation, the oscillation error was suppressed and the stress distribution was interpreted as relatively uniform.

• The Journal of the Acoustical Society of Korea
• /
• v.31 no.1
• /
• pp.19-28
• /
• 2012

The time-domain solution from the Kirchhoff integral equation for an exterior problem is not unique at certain eigen-frequencies associated with the fictitious internal modes as happening in frequency-domain analysis. One of the solution methods is the CHIEF (Combined Helmholtz Integral Equation Formulation) approach, which is based on employing additional zero-pressure constraints at some interior points inside the body. Although this method has been widely used in frequency-domain boundary element method due to its simplicity, it was not used in time-domain analysis. In this work, the CHIEF approach is formulated appropriately for time-domain acoustic boundary element method by constraining the unknown surface pressure distribution at the current time, which was obtained by setting the pressure at the interior point to be zero considering the shortest retarded time between boundary nodes and interior point. Sound radiation of a pulsating sphere was used as a test example. By applying the CHIEF method, the low-order fictitious modes could be damped down satisfactorily, thus solving the non-uniqueness problem. However, it was observed that the instability due to high-order fictitious modes, which were beyond the effective frequency, was increased.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.14 no.3
• /
• pp.11-22
• /
• 2010

An analysis method is proposed for the transient linear or nonlinear analysis of dynamic interactions between a flexible dam body and reservoir impounding compressible water under earthquake loadings. The coupled dam-reservoir system consists of three substructures: (1) a dam body with linear or nonlinear behavior; (2) a semi-infinite fluid region with constant depth; and (3) an irregular fluid region between the dam body and far field. The dam body is modeled with linear and/or nonlinear finite elements. The far field is formulated as a displacement-based transmitting boundary in the frequency domain that can radiate energy into infinity. Then the transmitting boundary is transformed for the direct coupling in the time domain. The near field region is modeled as a compressible fluid contained between two substructures. The developed method is verified and applied to various earthquake response analyses of dam-reservoir systems. Also, the method is applied to a nonlinear analysis of a concrete gravity dam. The results show the location and severity of damage demonstrating the applicability to the seismic evaluation of existing and new dams.

• The Journal of the Acoustical Society of Korea
• /
• v.27 no.3
• /
• pp.140-148
• /
• 2008

In this study, a numerical method for a time-domain acoustic wave backscattering analysis is established based on a physical optics and a Fourier transform. The frequency responses of underwater targets are calculated based on physical optics derived from the Kirchhoff-Helmholtz integral equation by applying Kirchhoff approximation and the time-domain signals are simulated taking inverse fast Fourier transform to the obtained frequency responses. Particularly, the adaptive triangular beam method is introduced to calculate the areas impinged directly by acoustic incident wave and the virtual surface concept is adopted to consider the multiple reflection effect. The numerical analysis result for an acoustic plane wave field incident normally upon a square flat plate is coincident with the result by the analytic time-domain physical optics derived theoretically from a conventional physical optics. The numerical simulation result for a hemi-spherical end-capped cylinder model is compared with the measurement result, so that it is recognized that the presented method is valid when the specular reflection effect is predominant, but, for small targets, gives errors due to higher order scattering components. The numerical analysis of an idealized submarine shows that the established method is effectively applicable to large and complex-shaped underwater targets.

• Korean Journal of Optics and Photonics
• /
• v.16 no.3
• /
• pp.182-190
• /
• 2005

The fiber Raman amplifier(FRA) is a distinctly advantageous technology. Due to its wider, flexible gain bandwidth, and intrinsically lower noise characteristics, FRA has become an indispensable technology of today. Various FRA modeling methods, with different levels of convergence speed and accuracy, have been proposed in order to gain valuable insights for the FRA dynamics and optimum design before real implementation. Still, all these approaches share the common platform of coupled ordinary differential equations(ODE) for the Raman equation set that must be solved along the long length of fiber propagation axis. The ODE platform has classically set the bar for achievable convergence speed, resulting exhaustive calculation efforts. In this work, we propose an alternative, highly efficient framework for FRA analysis. In treating the Raman gain as the perturbation factor in an adiabatic process, we achieved implementation of the algorithm by deriving a recursive relation for the integrals of power inside fiber with the effective length and by constructing a matrix formalism for the solution of the given FRA problem. Finally, by adiabatically turning on the Raman process in the fiber as increasing the order of iterations, the FRA solution can be obtained along the iteration axis for the whole length of fiber rather than along the fiber propagation axis, enabling faster convergence speed, at the equivalent accuracy achievable with the methods based on coupled ODEs. Performance comparison in all co-, counter-, bi-directionally pumped multi-channel FRA shows more than 102 times faster with the convergence speed of the Average power method at the same level of accuracy(relative deviation < 0.03dB).