• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.51-56
• /
• 2008

The waveform inversion for isotropic media has ever been studied since the 1980s, but there has been few studies for anisotropic media. We present a seismic waveform inversion algorithm for 2-D heterogeneous transversely isotropic structures. A cell-based finite difference algorithm for anisotropic media in time domain is adopted. The steepest descent during the non-linear iterative inversion approach is obtained by backpropagating residual errors using a reverse time migration technique. For scaling the gradient of a misfit function, we use the pseudo Hessian matrix which is assumed to neglect the zero-lag auto-correlation terms of impulse responses in the approximate Hessian matrix of the Gauss-Newton method. We demonstrate the use of these waveform inversion algorithm by applying them to a two layer model and the anisotropic Marmousi model data. With numerical examples, we show that it's difficult to converge to the true model when we assumed that anisotropic media are isotropic. Therefore, it is expected that our waveform inversion algorithm for anisotropic media is adequate to interpret real seismic exploration data.

• Journal of Astronomy and Space Sciences
• /
• v.20 no.1
• /
• pp.63-82
• /
• 2003

Using the Schwarzschild's linear programming technique, we obtained the general solutions of the collisionless Bolzmann equation describing the spherical galaxy in dynamical equilibrium. From this calculation we have confirmed the existence of isotropic spherical galaxies obeying a de Vaucouleurs'law which includes a dark halo. The flattening profile of the velocity dispersion curve seen in the elliptical galaxies can be explained as the increase of mass to light ratio in this dark matter. The space density distribution of this dark matter shows that the core radius of the dark matter is smaller than the effective radius of the galaxy.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.3B
• /
• pp.247-257
• /
• 2009

In the present paper, turbulent open-channel flows with alternate vegetated zones are numerically simulated using threedimensional model. The Reynolds-averaged Navier-Stokes Equations are solved with the ${\kappa}-{\varepsilon}$ model. The CFD code developed by Olsen(2004) is used for the present study. For model validation, the partly vegetated channel flows are simulated, and the computed depth-averaged mean velocity and Reynolds stress are compared with measured data in the literature. Comparisons reveal that the present model successfully predicts the mean flow and turbulent structures in vegetated open-channel. However, it is found that the ${\kappa}-{\varepsilon}$ model cannot accurately predict the momentum transfer at the interface between the vegetated zone and the non-vegetated zone. It is because the ${\kappa}-{\varepsilon}$ model is the isotropic turbulence model. Next, the open channel flows with alternate vegetated zones are simulated. The computed mean velocities are compared well with the previously reported measured data. Good agreement between the simulated results and the experimental data was found. Also, the turbulent flows are computed for different densities of vegetation. It is found that the vegetation curves the flow and the meandering flow pattern becomes more obvious with increasing vegetation density. When the vegetation density is 9.97%, the recirculation flows occur at the locations opposite to the vegetation zones. The impacts of vegetation on the flow velocity and the water surface elevation are also investigated.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.37 no.5A
• /
• pp.298-304
• /
• 2012

In this paper, we propose robust estimators for the frequency offset of orthogonal frequency division multiplexing in non-Gaussian noise environments. We first propose a maximum-likelihood (ML) estimator in non-Gaussian noise modeled as a complex isotropic Cauchy process, and then, we present a simpler suboptimal estimator based on the ML estimator. From numerical results, it is demonstrated that the proposed estimators not only outperform the conventional estimators, but also have a robustness in non-Gaussian noise environments.

• Journal of Dental Rehabilitation and Applied Science
• /
• v.25 no.3
• /
• pp.255-265
• /
• 2009

Recently many studies have been published on application of immediate loaded implants. However, the immediate loading protocol has not been well documented. The purpose of the present study was to evaluate the stress distribution between bone-implant interfaces and the effect of implant length in the anterior maxilla using 3 dimensional finite element analyses. The diameter 4.0 mm threaded type implants with different length(8.5 mm, 10.0 mm, 11.5 mm, 13.0 mm, 15.0 mm) were used in this study. The bone quality of anterior maxillary bone block was assumed to D3 bone. Bone-implant interfaces of immediately loaded implant were constructed using a contact element for simulating the non osseointegration status. For simplification of all the processing procedures, all of the material assumed to be homogenous, isotropic, and linearly elastic. The 178 N of static force was applied on the middle of the palatoincisal line angle of the abutment with $120^{\circ}$ angle to the long axis of abutment. Maximum von Mises stress were concentrated on the labial cortical bone of the implant neck area, especially at the cortical-cancellous bone interfaces. Compared the different length, highest peak stress value was observed at the 8.5 mm implants and the results indicated a tendency towards favorable stress distribution on the bone, when the length was increased. Presence of cortical bone was very important to immediate loading, and it appears that implants of a length more than 13 mm are preferable for immediate loading at the anterior maxilla.

• Geophysics and Geophysical Exploration
• /
• v.13 no.4
• /
• pp.307-314
• /
• 2010

Although conventional seismic data processing is based on the assumption that the media are isotropic, the subsurface is often anisotropy in shale formation or carbonate with cracks and fractures. This paper presents the anisotropic parameter and seismic modeling in transversely isotropic media with a vertical symmetry axis using seismic physical modeling. The experiment was successfully carried out with VTI media, laminated bakelite material, using contact transducer of p and s-wave transmission. The variation of velocities with angle of incidence was clearly shown in anisotropic material. Comparing these velocities with the calculated phase velocities, the (P) and (S)-wave velocity observed in anisotropic material was a very good agreement with the calculated values. Anisotropic parameter ${\varepsilon}$, ${\delta}$, ${\gamma}$ was estimated by using Lame's constant calculated from the observed velocity. For the purpose of testing (S)-wave polarization, a birefringence experiment was carried out. The higher velocity was associated with the polarization parallel to the fracture, and the lower velocity was associated with the polarization perpendicular to the fracture.

• Applied Chemistry for Engineering
• /
• v.9 no.6
• /
• pp.822-828
• /
• 1998

Binary system of water and polyoxyethylene dodecyl ether, $C_{12}H_{25}(OCH_2CH_2)nOH$, have been studied by $^1H$ NMR techniques. For n=5($C_{12}EO_5$) and n=8($C_{12}EO_8$), the self-diffusion coefficients of nonionic surfactants in the isotropic phase($L_1$) have been measured by using pulsed field gradient technique for a range of temperature and concentrations. In addition the line widths of the different proton signals have been monitored, and samples of some liquid crystalline characteristic were also studied. Dramatic Broadening of the methylene signals of the alkyl($C_{12}H_{25}$) chain is observed as the hexagonal liquid crystalline phase is approached in the $C_{12}EO_5-$water system, while only small broadening is observed in the $C_{12}EO_8-$water system. It was shown that there was a growth of $C_{12}EO_5$ micelles to rods with increasing concentrations, while the $C_{12}EO_8-$ micelles at low temperature remain small in the concentration range. The self-diffusion coefficients of the surfactants decrease rapidly with increasing concentration until a minimum is reached after which there is slow increase. The location of the minimum point occurs at lower concentrations the temperature is close to the cloud point, where the system separate into two isotropic phase. In the line width studies, broadening is found at a certain temperature interval when the concentration is increased in the $C_{12}EO_5$ system. The results indicate that the surfactant aggregates grow in size at the cloud point is approached. The aggregates seem to be flexible and probably not to be of a definite shape close to the cloud point. In the $C_{12}EO_8$ system, the micelles are much less affected by an increase in temperature and micellar growth can't be unambiguously established. The methylene signals of the ethylene oxide moieties consistantly show narrower $^1H$ signals, showing that in the aggregates they are less ordered than the chain methylenes. The various changes in aggregate size and shape are correlated with the stability ranges of the isotropic and liquid crystalline phases according to phase diagrams from the literature. Both aggregate size and phase structure are in qualitative agreement with concentration based on the effective shape of the molecules at different temperature and concentration.

• Geophysics and Geophysical Exploration
• /
• v.16 no.2
• /
• pp.71-78
• /
• 2013

The one-way wave equation migration is much more computationally efficient comparing with reverse time migration and it can provide better image than the migration algorithm based on the ray theory. We have developed the prestack depth migration module adopting (GS) propagator designed for vertical transverse isotropic media. Since GS propagator considers the higher-order term by expanding the Taylor series of the vertical slowness in the thin slab of the phase-screen propagator, the GS migration can offer more correct image for the complex subsurface with large lateral velocity variation or steep dip. To verify the validity of the developed GS migration module, we analyzed the accuracy with the order of the GS propagator for VTI media (GSVTI propagator) and confirmed that the accuracy of the wavefield propagation with the wide angles increases as the order of the GS propagator increases. Using the synthetic seismic data, we compared the migration results obtained from the isotropic GS migration module with the anisotropic GS migration module. The results show that the anisotropic GS migration provides better images and the improvement is more evident on steeply dipping structures and in a strongly anisotropic medium.

• Geophysics and Geophysical Exploration
• /
• v.18 no.2
• /
• pp.74-84
• /
• 2015

To yield physically meaningful images through elastic reverse-time migration, the wavefield separation which extracts P- and S-waves from reconstructed vector wavefields by using elastic wave equation is prerequisite. For expanding the application of the elastic reverse-time migration to anisotropic media, not only the anisotropic modelling algorithm but also the anisotropic wavefield separation is essential. The anisotropic wavefield separation which uses pseudo-derivative filters determined according to vertical velocities and anisotropic parameters of elastic media differs from the Helmholtz decomposition which is conventionally used for the isotropic wavefield separation. Since applying these pseudo-derivative filter consumes high computational costs, we have developed the efficient anisotropic wavefield separation algorithm which has capability of parallel computing by using GPUs (Graphic Processing Units). In addition, the highly efficient anisotropic elastic reverse-time migration algorithm using MPI (Message-Passing Interface) and incorporating the developed anisotropic wavefield separation algorithm with GPUs has been developed. To verify the efficiency and the validity of the developed anisotropic elastic reverse-time migration algorithm, a VTI elastic model based on Marmousi-II was built. A synthetic multicomponent seismic data set was created using this VTI elastic model. The computational speed of migration was dramatically enhanced by using GPUs and MPI and the accuracy of image was also improved because of the adoption of the anisotropic wavefield separation.

• Journal of radiological science and technology
• /
• v.32 no.3
• /
• pp.299-306
• /
• 2009

Purpose : To compare radiation dose for coronary CT angiography (CTA) obtained with 6 examination protocols such as a retrospectively ECG gated helical scan, a prospectively ECG gated sequential scan, low kVp technique, and cardiac dose modulation technique. Materials and Methods : Coronary CTA was performed by using 6 current clinical protocols to evaluate effective dose and organ dose in primary beam area with anthropomorphic female phantom and glass dosimetric system in 64 channel multi-detector CT. After acquiring topograms of frontal and lateral projection with 80 kVp and 10 mA, main coronary scan was done with 0.35 sec tube rotation time, 40 mm collimation ($0.625\;mm{\times}64\;ea$), small scan field of view (32 cm diameter), 105 mm scan length. Heart beat rate of phantom was maintained 60 bpm in ECG gating. In constant mAs technique 120 kVp, 600 mA was used, and 100 kVp for low kVp technique. In a retrospectively ECG gated helical CT technique 0.22 pitch was used, peak mA (600 mA) was adopted in range of $40{\sim}80%$ of R-R interval and 120mA(80% reduction) in others with cardiac dose modulation. And 210 mAs was used without cardiac dose modulation. In a prospectively ECG gated sequential CT technique data were acquired at 75% R-R interval (middle diastolic phase in cardiac cycle), and 120 msec additional padding of the tube-on time was used. For effective dose calculation region specific conversion factor of dose length product in thorax was used, which was recommended by EUR 16262. Results : The mean effective dose for conventional coronary CTA without cardiac dose modulation in a retrospectively ECG gated helical scan was 17.8 mSv, and mean organ dose of heart was 103.8 mGy. With low kVp and cardiac dose modulation the mean effective dose showed 54.5% reduction, and heart dose showed 52.3% reduction, compared with that of conventional coronary CTA. And at the sequential scan(SnapShot pulse mode) under prospective ECG gating the mean effective dose was 4.9 mSv, this represents an 72.5% reduction compared with that of conventional coronary CTA. And heart dose was 33.8 mGy, this represents 67.4% reduction. In the sequential scan technique under prospective ECG gating with low kVp the mean effective dose was 3.0 mSv, this represents an 83.2% reduction compared with that of conventional coronary CTA. And heart dose was 17.7 mGy, this represents an 82.9% reduction. Conclusion : In coronary CTA at retrospectively ECG gated helical scan, cardiac dose modulation technique using low kVp reduced dose to 50% above compared with the conventional helical scan. And the prospectively ECG gated sequential scan offers substantially reduced dose compared with the traditional retrospectively ECG gated helical scan.