• Tunnel and Underground Space
• /
• v.23 no.6
• /
• pp.493-505
• /
• 2013

This numerical study investigates seismicity and fault slip induced by fluid injection in deep geothermal reservoir with pre-existing fractures and fault. Particle Flow Code 2D is used with additionally implemented hydro-mechanical coupled fluid flow algorithm and acoustic emission moment tensor inversion algorithm. The output of the model includes spatio-temporal evolution of induced seismicity (hypocenter locations and magnitudes) and fault deformation (failure and slip) in relation to fluid pressure distribution. The model is applied to a case of fluid injection with constant rates changing in three steps using different fluid characters, i.e. the viscosity, and different injection locations. In fractured reservoir, spatio-temporal distribution of the induced seismicity differs significantly depending on the viscosity of the fracturing fluid. In a fractured reservoir, injection of low viscosity fluid results in larger volume of induced seismicity cloud as the fluid can migrate easily to the reservoir and cause large number and magnitude of induced seismicity in the post-shut-in period. In a faulted reservoir, fault deformation (co-seismic failure and aseismic slip) can occur by a small perturbation of fracturing fluid (<0.1 MPa) can be induced when the injection location is set close to the fault. The presented numerical model technique can practically be used in geothermal industry to predict the induced seismicity pattern and magnitude distribution resulting from hydraulic stimulation of geothermal reservoirs prior to actual injection operation.

• Journal of the Korean earth science society
• /
• v.28 no.7
• /
• pp.857-870
• /
• 2007

The structure and evolution of a thunderstorm outflow in two dimensions with no environmental wind are investigated using a cloud-resolving model with explicit liquid-ice phase microphysical processes (ARPS: Advanced Regional Prediction System). The turbulence structure of the outflow is explicitly resolved with a high-resolution grid size of 50m. The simulated single-cell storm and its associated Kelvin-Helmholtz (KH) billows are found to have the lift stages of development maturity, and decay. The secondary pulsation and splitting of convective cells resulted from interactions between cloud dynamics and microphysics are observed. The cooled downdrafts caused by the evaporation of rain and hail in the relatively dry lower atmosphere result in thunderstorm cold-air outflow. The outflow head propagates with almost constant speed. The KH billows formed by the KH instability cause turbulence mixing from the top of the outflow and control the structure of the outflow. Ihe KH billows are initiated at the outflow head, and pow and decay as moving rearward relative to the gust front. The numerical simulation results of the ratio of the horizontal wavelength of the fastest growing perturbation to the critical shear-layer depth and the ratio of the horizontal wavelength of the billow to its maximum amplitude are matched well with the results of other studies.

• Journal of the Korean earth science society
• /
• v.22 no.6
• /
• pp.512-527
• /
• 2001

Separated flows passed complex geometries are modeled by discrete vortex techniques. The flows are assumed to be rotational and inviscid, and a new techlnique is described to determine the stream functions for linear shear profiles. The geometries considered are the snow cornice and the backward-facing step, whose edges allow for the separation of the flow and reattachment downstream of the recirculation regions. A point vortex has been added to the flows in order to constrain the separation points to be located at the edges, while the conformal mappings have been modified in order to smooth the sharp edges and to let the separation points free to oscillate around the points of maximum curvature. Unsteadiness is imposed to the flow by perturbing the vortex location, either by displacing the vortex from the equilibrium, or by imposing a random perturbation with zero mean to the vortex in equilibrium. The trajectories of passive scalars continuously released upwind of the separation point and trapped by the recirculating bubble are numerically integrated, and concentration time series are calculated at fixed locations downwind of the reattachment points. This model proves to be capable of reproducing the trapping and intermittent release of scalars, in agreement with the simulation of the flow passed a snow cornice performed by a discrete multi-vortex model, as well as with direct numerical simulations of the flow passed a backward-facing step. The results of simulation indicate that for flows undergoing separation and reattachment the unsteadiness of the recirculating bubble is the main mechanism responsible for the intense large-scale concentration fluctuations downstream.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.9
• /
• pp.859-866
• /
• 2010

Nonlinear dynamics of pulsating instability-diffusional-thermal instability with Lewis numbers sufficiently higher than unity-in counterflow diffusion flames, is numerically investigated by imposing a Damkohler number perturbation. The flame evolution exhibits three types of nonlinear behaviors, namely, decaying pulsating behavior, diverging behavior (which leads to extinction), and stable limit-cycle behavior. The stable limit-cycle behavior is observed in counterflow diffusion flames, but not in diffusion flames with a stagnant mixing layer. The critical value of the perturbed Damkohler number, which indicates the region where the three different flame behaviors can be observed, is obtained. A stable simple limit cycle, in which two supercritical Hopf bifurcations exist, is found in a narrow range of Damkohler numbers. As the flame temperature is increased, the stable simple limit cycle disappears and an unstable limit cycle corresponding to subcritical Hopf bifurcation appears. The period-doubling bifurcation is found to occur in a certain range of Damkohler numbers and temperatures, which leads to extend the lower boundary of supercritical Hopf bifurcation.

• Korean journal of food and cookery science
• /
• v.31 no.2
• /
• pp.175-184
• /
• 2015

The purpose of this study was to determine the optimal mixing ratio of Agaricus blazei Murill powder and butter in the preparation of cookies. The experimental design utilized herein was based on central composite design for response surface methodology, which included 10 experimental points, including 2 replicates for Agaricus blazei Murill and butter. The physical, mechanical, and sensory properties of the test were measured, and these values were applied to the mathematical models. A canonical form and perturbation plot showed the influence of each ingredient on the final mixed product. The spread ratio increased significantly with an increase in Agaricus blazei Murill powder and butter (p<0.05). The response surface methodology was applied to evaluate the effect of Agaricus blazei Murill powder and butter on cookie moisture and color (L, a) (p<0.001). Sensory evaluation showed significant values for color (p<0.05), flavor (p<0.05), texture (p<0.05) and overall quality (p<0.01) in the predicted model. The optimum formulation by numerical and graphical methods was calculated as follows: Agaricus blazei Murill powder 3.63 g, butter 55.37 g.

• Korean journal of food and cookery science
• /
• v.26 no.6
• /
• pp.707-716
• /
• 2010

Germinated brown rice is a useful material that contains natural nutrients and biologically active substances. This study was conducted to develop a optimal composite recipe for functional germinated brown rice cookie with added spinach (spinacia oleracea L) powder and with high preference in all age groups. The experimental design utilized herein was based on the Central composite design methodology of response surface, which included 16 experimental points, including 2 replicates for spinach, sugar, and butter. The physical, mechanical, and sensory properties of the test were measured, and these values were applied to the mathematical models. A canonical form and perturbation plot showed the influence of each ingredient on the mixed final product. The results of spread ratio decreased significantly with increases in spinach powder, sugar, and butter(p<0.01). L value shows that lightness decreased significantly with an increase in sugar(p<0.05), but the a and b values were not significant. Sensory evaluation found significant values for color(p<0.05), flavor(p<0.05), texture(p<0.05) and overall quality(p<0.05) in the predicted model. The optimum formulation by the numerical and graphical methods was calculated as follows: spinach powder 3 g, sugar 35.5 g butter 54.4 g.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.16 no.9
• /
• pp.6221-6226
• /
• 2015

In this paper, waveguide analysis was interpreted as an optical waveguide using the feedback perturbation method simple and easy to apply compared to the analysis method, while the other almost identical to the numerical calculation method. In addition, efficient coupling with an optical transmission line of optical fiber and the waveguide form the thin film of different functional elements is required in order to achieve the full optical communication system. However, problems arise, such as the light field (field) and the decrease of the access efficiency due to inconsistency in the distribution characteristics of the connection surface by the difference in size of the cross section thereof when connecting the optical fiber and the waveguide directly to the combination of a thin film. Therefore propose a new type of connector structure to increase the efficiency of the connection when connecting the optical fiber waveguide and the thin film was analyzed by applying a coupled mode theory, the binding efficiency of the modified contactor. And by diffusing Ti on the $LiNbO_3$ substrate and a wide range of applications in the manufacture of integrated optical material made of a current low-loss Ti: $LiNbO_3$ optical waveguide and making the Y-branch waveguide, and the properties were confirmed.

• Journal of Korean Society of Steel Construction
• /
• v.20 no.1
• /
• pp.21-31
• /
• 2008

response surface method (RSM) is widely used to evaluate th e extremely smal probability of ocurence or toanalyze the reliability of very complicated structures. Althoug h Monte-Carlo Simulation (MCS) technique can evaluate any system, the procesing time of MCS dependson the reciprocal num ber of the probability of failure. The stochastic finite element method could solve thislimitation. However, it is limit ed to the specific program, in which the mean and coeficient o f random variables are programed by a perturbation or by a weigh ted integral method. Therefore, it is not aplicable when erequisite programing. In a few number of stage analyses, RSM can construct a regresion model from the response of the c omplicated structural system, thus, saving time and efort significantly. However, the acuracy of RSM depends on the dist ance of the axial points and on the linearity of the limit stat e functions. To improve the convergence in exact solution regardl es of the linearity limit of state functions, an improved adaptive response surface method is developed. The analyzed res ults have ben verified using linear and quadratic forms of response surface functions in two examples. As a result, the be st combination of the improved RSM techniques is determined and programed in a numerical code. The developed linear adapti ve weighted response surface method (LAW-RSM) shows the closest converged reliability indices, compared with quadratic form or non-adaptive or non-weighted RSMs.