• Korean Journal of Optics and Photonics
• /
• v.25 no.4
• /
• pp.193-199
• /
• 2014

The net effect of the emitter orientation, Mie scatters, and pillow lenses on the outcoupling efficiency (OCE) of a bottom-emitting OLED having an internal photonic crystal layer was investigated by a combined optical simulation based on the finite-difference time-domain method (FDTD) and the ray-tracing technique. The simulation showed that when the emitter orientation was horizontal with respect to the OLED surface, the OCE could be increased by 54% when a photonic crystal layer was employed, while it could be improved by 86% under optimized conditions of Mie scatters and pillow lenses applied to the glass substrate. The peculiar intensity distribution of the OLED, caused by the periodic lattice structure of the photonic crystal layer, could be ameliorated by inserting Mie scatters into the glass substrate. This study suggests that conventional outcoupling structures combined with control of the emitter orientation could improve the OCE substantially.

• Journal of electromagnetic engineering and science
• /
• v.6 no.1
• /
• pp.47-52
• /
• 2006

The objective of this study is to evaluate the effects of size and permittivity on the specific absorption rate(SAR) values of rat brains during microwave exposure at mobile phone frequency bands. A finite difference time domain (FDTD) technique with perfect matching layer(PML) absorbing boundaries is used for this evaluation process. A color coded digital image of the Sprague Dawley(SD) rat based on magnetic resonance imaging(MRI) is used in FDTD calculation with appropriate permittivity values corresponding to different tissues for 3, 4, 7, and 10 week old rats. This study is comprised of three major parts. First, the rat model structure is scaled uniformly, i.e., the rat size is increased without change in permittivity. The simulated SAR values are compared with other experimental and numerical results. Second, the effect of permittivity on SAR values is examined by simulating the microwave exposure on rat brains with various permittivity values for a fixed rat size. Finally, the SAR distributions in depth, and the brain-averaged SAR and brain 1 voxel peak SAR values are computed during the microwave exposure on a rat model structure when both size and permittivity have varied corresponding to different ages ranging from 3 to 10 weeks. At 900 MHz, the simulation results show that the brain-averaged SAR values decreased by about 54 % for size variation from the 3 week to the 10 week-old rat model, while the SAR values decreased only by about 16 % for permittivity variation. It is found that the brain averaged SAR values decreased by about 63 % when the variations in size and permittivity are taken together. At 1,800 MHz, the brain-averaged SAR value is decreased by 200 % for size variation, 9.7 % for permittivity variation, and 207 % for both size and permittivity variations.

• The Journal of Engineering Geology
• /
• v.14 no.1
• /
• pp.21-28
• /
• 2004

Three-dimensional (3-D) resistivity inversion including a topographic effect can be considered theoretically to be the technique of acquiring the most accurate image in the interpretation of resistivity data, because it includes characteristic image that the actual subsurface structure is 3-D. In this study, a finite-element method was used as the numerical method in modeling, and the efficiency of Jacobian calculation has been maximized with sensitivity analysis for the destination block in inversion process. Also, during the iterative inversion, the resolution of inversion can be improved with the method of selecting the optimal value of Lagrange multiplier yielding minimum RMS(root mean square) error in the parabolic equation. In this paper, we present synthetic examples to compare the difference between the case which has the toprographic effect and the other case which has not the effect in the inversion process.

• Journal of Korea Water Resources Association
• /
• v.36 no.1
• /
• pp.129-142
• /
• 2003

In this study, the conduct of underground water flow system with 3 dimensions is interpreted in order to examine closely the actual condition regarding the flow of the underground water which is forecast from the tunnel segment and the interpretation result which selects the design and the construction technique of the tunnel segment was applied. Also, an obstacle to construct that relates with the underground water flow in construction duration in advance will be able to apply with information that is necessary in order to establish the countermeasure. The objective tunnel is the BEOPGI tunnel segment that is 2 parallel tunnels that are a one-way 2 lane and the parameters of the MODFLOW model executing the boring investigation and the permeability examination were presumed. The underground water flow of the excavation tunnel inside was interpreted by the MODFLOW model using the parameters which is presumed and two values which compared with calculated value and observed value are the same almost. Also, when the underground water discharge quantify that followed in tunnel excavation tries to compare, the underground water total discharge quantity from tunnel point of start until destination was presumed as 0.0269㎥/day/$m^2$.

• Transactions of the Korean hydrogen and new energy society
• /
• v.32 no.4
• /
• pp.228-235
• /
• 2021

Polymer electrolyte membrane fuel cells (PEMFC) are currently being used in various transport applications such as drones, unmanned aerial vehicles, and automobiles. The power required is different according to the type of use, purpose, and the conditions adjusted using a cell stack. The fuel cell stack is compressed to reduce the size and prevent fuel leakage. The unit cells that make up the cell stack are subjected to compression by clamping force, which makes geometrical changes in the porous media and it impacts on cell performance. In this study, finite elements method (FEM) and computational fluid dynamics (CFD) analysis for the deformed unit cell considering the effects of clamping force is performed. First, structural analysis using the FEM technique over the deformed gas diffusion layer (GDL) considering compression is carried out, and the resulting porosity changed in the GDL is calculated. The PEMFC model is then verified by a three-dimensional, two-phase fuel cell simulation applying the physical properties and geometry obtained before and after compression. The detailed simulation results showed different concentration distributions of fuel between the original and deformed geometry, resulting in the difference in the distribution of current density is represented at compressed GDL region with low oxygen concentration.

• Journal of the Korean Geotechnical Society
• /
• v.23 no.8
• /
• pp.117-127
• /
• 2007

Slope stability analysis is a geotechnical engineering problem characterized by many sources of uncertainty. Some of them are connected to the uncertainties of soil properties involved in the analysis. In this paper, a numerical procedure for integrating commercial finite difference method into probabilistic analysis of slope stability is presented. Since the limit-state function cannot be expressed in an explicit form, the ANN-based response surface method is adopted to approximate the limit-state function and the first-, second-order reliability method and the Monte Carlo simulation technique are used to calculate the probability of failure. Probabilistic stability assessments for a hypothetical two-layer slope and the Sugar Creek embankment were performed to verify the application potential to the slope stability problems. The examples show the successful implementation and the possibility of the extension of the proposed procedure to the variety of geotechnical engineering problems.

• Journal of the Korean Society of Groundwater Environment
• /
• v.6 no.1
• /
• pp.23-32
• /
• 1999

A new numerical method using solution-adaptive grids (SAG) is developed to solve the Richards' equation (RE) for unsaturated flow in porous media. Using a grid generation technique, the SAG method automatically redistributes a fixed number of grid points during the flow process, so that more grid points are clustered in regions of large solution gradients. The method uses the coordinate transformation technique to employ a new transformed RE, which is solved with the standard finite difference method. The movement of grid points is incorporated into the transformed RE, and therefore all computation is performed on fixed grid points of the transformed domain without using any interpolation techniques. Thus, numerical difficulties arising from the movement of the wetting front during the infiltration process have been substantially overcome by the new method. Numerical experiments for an one-dimensional infiltration problem are presented to compare the SAG method to the modified Picard method using a fixed grid. Results show that accuracy of a SAG solution using 41 nodes is comparable with the solution of the fixed grid method using 201 nodes, while it requires only 50% of the CPU time. The global mass balance and the convergence of SAG solutions are strongly affected by the time step size (Δt) and the weighting parameter (${\gamma}$) used for generating solution-adaptive grids. Thus, the method requires automated readjustment of Δt and ${\gamma}$ to yield mass-conservative and convergent solutions, although it may increase computational costs. The method can be effective especially for simulating unsaturated flow and other transport problems involving the propagation of a sharp-front.

• International Journal of Aeronautical and Space Sciences
• /
• v.8 no.1
• /
• pp.10-20
• /
• 2007

For spacecraft attitude control, reaction wheel (RW) steering laws with more than three wheels for three-axis attitude control can be derived by using a control allocation (CA) approach.1-2 The CA technique deals with a problem of distributing a given control demand to available sets of actuators.3-4 There are many references for CA with applications to aerospace systems. For spacecraft, the control torque command for three body-fixed reference frames can be constructed by a combination of multiple wheels, usually four-wheel pyramid sets. Multi-wheel configurations can be exploited to satisfy a body-axis control torque requirement while satisfying objectives such as minimum control energy.1-2 In general, the reaction wheel steering laws determine required torque command for each wheel in the form of matrix pseudo-inverse. In general, the attitude control command is generated in the form of a feedback control. The spacecraft body angular rate measured by gyros is used to estimate angular displacement also.⁵ Combination of the body angular rate and attitude parameters such as quaternion and MRPs(Modified Rodrigues Parameters) is typically used in synthesizing the control command which should be produced by RWs.¹ The attitude sensor signals are usually corrupted by noise; gyros tend to contain errors such as drift and random noise. The attitude determination system can estimate such errors, and provide best true signals for feedback control.⁶ Even if the attitude determination system, for instance, sophisticated algorithm such as the EKF(Extended Kalman Filter) algorithm⁶, can eliminate the errors efficiently, it is quite probable that the control command still contains noise sources. The noise and/or other high frequency components in the control command would cause the wheel speed to change in an undesirable manner. The closed-loop system, governed by the feedback control law, is also directly affected by the noise due to imperfect sensor characteristics. The noise components in the sensor signal should be mitigated so that the control command is isolated from the noise effect. This can be done by adding a filter to the sensor output or preventing rapid change in the control command. Dynamic control allocation(DCA), recently studied by Härkegård, is to distribute the control command in the sense of dynamics⁴: the allocation is made over a certain time interval, not a fixed time instant. The dynamic behavior of the control command is taken into account in the course of distributing the control command. Not only the control command requirement, but also variation of the control command over a sampling interval is included in the performance criterion to be optimized. The result is a control command in the form of a finite difference equation over the given time interval.⁴ It results in a filter dynamics by taking the previous control command into account for the synthesis of current control command. Stability of the proposed dynamic control allocation (CA) approach was proved to ensure the control command is bounded at the steady-state. In this study, we extended the results presented in Ref. 4 by adding a two-step dynamic CA term in deriving the control allocation law. Also, the strict equality constraint, between the virtual and actual control inputs, is relaxed in order to construct control command with a smooth profile. The proposed DCA technique is applied to a spacecraft attitude control problem. The sensor noise and/or irregular signals, which are existent in most of spacecraft attitude sensors, can be handled effectively by the proposed approach.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2013.02a
• /
• pp.334-334
• /
• 2013

In the past decade, the infrared detectors based on intersubband transition in quantum dots (QDs) have attracted much attention due to lower dark currents and increased lifetimes, which are in turn due a three-dimensional confinement and a reduction of scattering, respectively. In parallel, focal plane array development for infrared imaging has proceeded from the first to third generations (linear arrays, 2D arrays for staring systems, and large format with enhanced capabilities, respectively). For a step further towards the next generation of FPAs, it is envisioned that a two-dimensional metal hole array (2D-MHA) structures will improve the FPA structure by enhancing the coupling to photodetectors via local field engineering, and will enable wavelength filtering. In regard to the improved performance at certain wavelengths, it is worth pointing out the structural difference between previous 2D-MHA integrated front-illuminated single pixel devices and back-illuminated devices. Apart from the pixel linear dimension, it is a distinct difference that there is a metal cladding (composed of a number of metals for ohmic contact and the read-out integrated circuit hybridization) in the FPA between the heavily doped gallium arsenide used as the contact layer and the ROIC; on the contrary, the front-illuminated single pixel device consists of two heavily doped contact layers separated by the QD-absorber on a semi-infinite GaAs substrate. This paper is focused on analyzing the impact of a two dimensional metal hole array structure integrated to the back-illuminated quantum dots-in-a-well (DWELL) infrared photodetectors. The metal hole array consisting of subwavelength-circular holes penetrating gold layer (2DAu-CHA) provides the enhanced responsivity of DWELL infrared photodetector at certain wavelengths. The performance of 2D-Au-CHA is investigated by calculating the absorption of active layer in the DWELL structure using a finite integration technique. Simulation results show the enhanced electric fields (thereby increasing the absorption in the active layer) resulting from a surface plasmon, a guided mode, and Fabry-Perot resonances. Simulation method accomplished in this paper provides a generalized approach to optimize the design of any type of couplers integrated to infrared photodetectors.

• Applied Biological Chemistry
• /
• v.43 no.2
• /
• pp.116-123
• /
• 2000

A numerical model of three-dimensional soil water distribution for drip irrigation management under cropped conditions was developed using Richards equation in Cartesian coordinates. The model accounts for both seasonal and diurnal changes in evaporation and transpiration, and the growth of plant root and the shape of root zone. Solutions were numerically approximated using the Crank-Nicolson implicit finite difference technique on the block-centered grid system and the Gauss-Seidel elimination in tandem. The model was tested under several conditions to allow the flow rates and configurations of drip emitters vary. In general, simulation results agreed well with experimental results and were as follows. The velocity of soil-water flow decreased drastically with distance from the drip source, and the rate of expansion of the wetted zone decreased rapidly during irrigation. The wetting front of wetted zone from a surface drip emitter traveled farther in vertical direction than in horizontal direction. Under this experimental weather condition, water use efficiency of a drip-irrigated apple field was greatest for 4-drip-emitter system buried at 25 cm, resulting from 10% increase in transpiration but 20% reduction in soil evaporation compared to those for surface 1-drip emitter system. Soil moisture retention curve obtained using disk tension infiltrometer showed significant difference from the curve obtained with pressure plate extractor.