• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.46 no.3
• /
• pp.210-218
• /
• 2018

In order to generate meshes automatically for finite element analysis of complex structures such as aircraft, a large number of triangular elements are typically created. However, triangular elements are less accurate than rectangular elements, so it is difficult to obtain a reliable solution. This problem can be improved through the meshfree method using the back cell integration. However, this method also causes some problems such as over-use of the integration points and inefficiency of the integral domain. In order to improve these problems, a method of performing integration by setting the integral area based on a node basis has been proposed, but in the case of incompressible material problems, the numerical accuracy deteriorates due to the vibration phenomenon of the solution. Therefore, in this paper, the modified meshfree method is proposed which sets the integral domain as an element domain instead of the nodal domain, and the proposed method improves the numerical instability caused by the conventional meshfree method without decreasing the accuracy regardles of the shape of integral domain. The effectiveness of the modified meshfree method is verified by using 2-D examples.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.25-33
• /
• 2001

Loop-loop electromagnetic (EM) survey in frequency domain has been carried out in order to provide basic solution to geotechnical applications. Source and receiver configuration may be horizontal co-planar (HCP) and/or vertical co-planar (VCP). Three quadrature components of mutual impedance ratio for each configuration are used to construct the subsurface image. For the purpose of obtaining the model response and validating the reasonable performance of the inversion, we obtained each responses of two-layered and three-layered earth models and two-dimensional (2-D) isolated anomalous body. The response of 2-D isolated anomalous body has been calculated using extended Born approximation for the solution of 2.5-D integral equation describing EM scattering problem. As a result of the least-squares inversion with variable Lagrangian multiplier, we could construct more resolvable image from HCP data than VCP data. Furthermore, joint inversion of HCP and VCP data made better stability and resolution of the inversion. Resistivity values, however, did not exactly match the true ones. Loop-loop EM field data was obtained with EM34-3XL system manufactured by Geonics Ltd. (Canada). Electrical resistivity survey was conducted on the same line for the comparison in advance. Since the constructed image from loop-loop EM data by 2-D inversion algorithm showed almost similar resistivity distribution to that from electrical resistivity one, we expect the developed 2.5-D loop-loop EM inversion program can be applied for the reconnaissance site survey.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.236-249
• /
• 2002

We have extended the three-dimensional (3-D) resistivity imaging algorithm to cover the 3-D resistivity tomography problem, where resistivity data are acquired using electrodes installed in several boreholes as well as at the earth surface. The imaging algorithm consists of the 3-D finite element forward modeling and least-squares inversion scheme, where the ACB (Active Constraint Balancing) is adopted to enhance the resolving power of the inversion. Sensitivity analysis with numerical verifications shows that 3-D resistivity tomography is a very appealing method and can be used to get 3-D attitude of subsurface structures with very high-resolution. Moreover, we could accurately handle the topography effect, which could cause artifacts in the resistivity tomography. In the application of 3-D resistivity tomography to the real field data set acquired at the quarry mine, we could derive a very reasonable and accurate image of the subsurface.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.4
• /
• pp.249-256
• /
• 2019

This paper presents an optimization framework of electrical impedance tomography for characterizing electrical conductivity profiles of concrete structures in two dimensions. The framework utilizes a partial-differential-equation(PDE)-constrained optimization approach that can obtain the spatial distribution of electrical conductivity using measured electrical potentials from several electrodes located on the boundary of the concrete domain. The forward problem is formulated based on a complete electrode model(CEM) for the electrical potential of a medium due to current input. The CEM consists of a Laplace equation for electrical potential and boundary conditions to represent the current inputs to the electrodes on the surface. To validate the forward solution, electrical potential calculated by the finite element method is compared with that obtained using TCAD software. The PDE-constrained optimization approach seeks the optimal values of electrical conductivity on the domain of investigation while minimizing the Lagrangian function. The Lagrangian consists of least-squares objective functional and regularization terms augmented by the weak imposition of the governing equation and boundary conditions via Lagrange multipliers. Enforcing the stationarity of the Lagrangian leads to the Karush-Kuhn-Tucker condition to obtain an optimal solution for electrical conductivity within the target medium. Numerical inversion results are reported showing the reconstruction of the electrical conductivity profile of a concrete specimen in two dimensions.

• Nuclear Medicine and Molecular Imaging
• /
• v.41 no.4
• /
• pp.317-325
• /
• 2007

Purpose: Biological parameters can be quantified using dynamic PET data with compartment modeling and Nonlinear Least Square (NLS) estimation. However, the generation of parametric images using the NLS is not appropriate because of the initial value problem and excessive computation time. In irreversible model, Patlak graphical analysis (PGA) has been commonly used as an alternative to the NLS method. In PGA, however, the start time ($t^*$, time where linear phase starts) has to be determined. In this study, we suggest a new Multiple Linear Analysis for irreversible radiotracer (MLAIR) to estimate fluoride bone influx rate (Ki). Methods: $[^{18}F]Fluoride$ dynamic PET scans was acquired for 60 min in three normal mini-pigs. The plasma input curve was derived using blood sampling from the femoral artery. Tissue time-activity curves were measured by drawing region of interests (ROls) on the femur head, vertebra, and muscle. Parametric images of Ki were generated using MLAIR and PGA methods. Result: In ROI analysis, estimated Ki values using MLAIR and PGA method was slightly higher than those of NLS, but the results of MLAIR and PGA were equivalent. Patlak slopes (Ki) were changed with different $t^*$ in low uptake region. Compared with PGA, the quality of parametric image was considerably improved using new method. Conclusion: The results showed that the MLAIR was efficient and robust method for the generation of Ki parametric image from $[^{18}F]Fluoride$ PET. It will be also a good alternative to PGA for the radiotracers with irreversible three compartment model.

• Economic and Environmental Geology
• /
• v.27 no.6
• /
• pp.557-561
• /
• 1994

Q estimates are made by direct measurements of energy loss per cycle from primary P and S waves, as a function of frequency. Assuming that intrinsic Q is frequency independent and scattering Q is frequency dependent over the frequencies of interest, the relative contributions of each, to a total observed Q, may be estimated. Test examples are produced by computing viscoelastic synthetic seismograms using a pseudo spectral solution with inclusion of relaxation mechanisms (for intrinsic Q) and a fractal distribution of scatterers (for scattering Q). The composite theory implies that when the total Q for S-waves is smaller than that for P-waves (the usual situation), intrinsic Q is dominating; when it is larger, scattering Q is dominating. In the inverse problem, performed by a global least squares search, intrinsic $Q_p$ and $Q_s$ estimates are reliable and unique when their absolute values are sufficiently low that their effects are measurable in the data. Large $Q_p$ and $Q_s$ have no measurable effect and hence are not resolvable. Standard deviation of velocity $({\sigma})$ and scatterer size (A) are less unique as they exhibit a tradeoff as predicted by Blair's equation. For the P-waves, intrinsic and scattering contributions are of approximately the same importance, for S-waves, the intrinsic contributions dominate.

• Analytical Science and Technology
• /
• v.31 no.4
• /
• pp.161-170
• /
• 2018

The epidemic of disorders associated with synthetic stimulants, such as methamphetamine (MA) and amphetamine (AP), is a health, social, legal, and financial problem. Owing to the high potential of their abuse and addiction, reliable analytical methods are required to detect and identify MA, AP, and their metabolites in biological samples. Thus, a dilute-and-shoot liquid chromatography-tandem mass spectrophotometry (LC-MS/MS) was developed for simultaneous determination of MA, 4-hydroxymethamphetamine (4HMA), AP, and 4-hydroxyamphetamine (4HA) in urine. Urine sample ($100{\mu}L$) was mixed with $50{\mu}L$ of mobile phase consisting of 0.4 % formic acid and methanol and $50{\mu}L$ of working internal-standard solution. Aliquots of $8{\mu}L$ diluted urine was injected into the LC-MS/MS system. For all analytes, chromatographic separation was performed using a C18 reversed-phase column with gradient elution and a total run time of 5 min. The identification and quantification were performed by multiple reaction monitoring (MRM). Linear least-squares regression was conducted to generate a calibration curve, with $1/x^2$ as the weighting factor. The linear ranges were 2.0-200, 1.0-800, and 10-2500 ng/mL for 4HA and 4HMA, AP, and MA, respectively. The inter- and intraday precisions were within 6.6 %, whereas the inter- and intraday accuracies ranged from -14.9 to 11.3 %. The low limits of quantification were 2.0 ng/mL (4HA and 4HMA), 1.0 ng/mL (AP), and 10 ng/mL (MA). The proposed method exhibited satisfactory selectivity, dilution integrity, matrix effect, and stability, which are required for validation. Moreover, the purification efficiency of high-speed centrifugation was clearly higher than 6-15 % for QC samples (n=5), which was higher than that of the membrane-filtration method. The applicability of the proposed method was tested by forensic analysis of urine samples from drug abusers.