• Journal of the Korean Institute of Gas
• /
• v.20 no.1
• /
• pp.13-22
• /
• 2016

This paper presents development of an appropriate procedure and flow chart to analyze shale gas production data obtained from a multi-fractured horizontal well according to flow characteristics in order to calculate an estimated ultimate recovery. Also, the technical considerations were proposed when a rate transient analysis was performed with field production data occurred to only $1^{st}$ transient flow. If production data show the $1^{st}$ transient flow from log-log and square root time plot analysis, production forecasting must be performed by applying different method as before and after of the end of $1^{st}$ linear flow. It is estimated by an area of stimulated reservoir volume which can be calculated from analysis results of micro-seismic data. If there are no bottomhole pressure data or micro-seismic data, an empirical decline curve method can be used to forecast production performance. If production period is relatively short, an accuracy of production data analysis could be improved by analyzing except the early production data, if it is necessary, after evaluating appropriation with near well data. Also, because over- or under-estimation for stimulated reservoir volume could take place according to analysis method or analyzer's own mind, it is necessary to recalculate it with fracture modeling, reservoir simulation and rate transient analysis, if it is necessary, after adequacy evaluation for fracture stage, injection volume of fracture fluid and productivity of producers.

• Journal of Environmental Science International
• /
• v.8 no.4
• /
• pp.423-429
• /
• 1999

This study was performed to investigate the fluid flow characteristics, the temperature distribution, humidity and PMV(Predicted Mean Vote of thermal sensation) distributions of the 3-dimensional room with side wall exhaust. The finite volume method and turbulence k-$\varepsilon$ models with the SIMPLE computational algorithm are used. As the results of the three dimensional simulations, the region of exceeding Y=1.5m was high temperature and humidity. The inlet velocity and temperature were influenced to the floor strongly, and the room PMV was about -1.0 except the inlets.

• 유체기계공업학회:학술대회논문집
• /
• 2000.12a
• /
• pp.169-175
• /
• 2000

A design method for transonic turbine blades is developed based on Navier-Stokes equations. The present computing process is done on the four separate steps, 1.e., determination of the blade profile, generation of the computational grids, cascade flow simulation and analysis of the computed results in the sense of the aerodynamic performance. The blade shapes are designed using the cubic polynomials under the control of the design parameters. Numerical methods for the flow equations are based on Van-Leer's FVS with an upwind TVD scheme on the finite volume. Applications are made to the VKI transonic rotor blades. Computed results are analyzed with respect to the aerodynamic performance and are compared with the experimental data.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.9-13
• /
• 2006

In the present study, a numerical analysis of transient mixer flow is performed considering free surface formation. The flow patterns and free surface shape in a mixers formed by flat paddle and pitched paddle impellers are predicted. In a flat paddle mixer, two flow circulation regions are formed due to strong radial flow, whereas one large circulation is formed in a pitched paddle mixer due to axial downward flow. These differences affect the free surface evolution and shape. It is seen from the results that a flat paddle mixer gives deeper free surface at center region than a pitched paddle mixer. The free surface of 8-blades-flat-paddle mixer is also simulated to compare with the available experimental and simulation results. The present computational results agree reasonably well with the experimental data.

• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.271-278
• /
• 2020

The investigation into a full-scale 27 m high, by 6 m wide, thermosyphon loop. The simulation model is based on a one-dimensional axially-symmetrical control volume approach, where the loop is divided into a series of discreet control volumes. The three conservation equations, namely, mass, momentum and energy, were applied to these control volumes and solved with an explicit numerical method. The flow is assumed to be quasi-static, implying that the mass-flow rate changes over time. However, at any instant in time the mass-flow rate is constant around the loop. The boussinesq approximation was invoked, and a reasonable correlation between the experimental and theoretical results was obtained. Experimental results are presented and the flow regimes of the working fluid inside the loop identified. The results indicate that a series of such thermosyphon loops can be used as a cavity cooling system and that the one-dimensional theoretical model can predict the internal temperature and mass-flow rate of the thermosyphon loop.

• 유체기계공업학회:학술대회논문집
• /
• 2004.12a
• /
• pp.456-461
• /
• 2004

In this study, the characteristics of the cavitation behavior around hydrofoil are investigated with the commercial CFD code CFX-5. Calculations are performed for NACA 64108 hydrofoil by solving the time-averaged turbulent Navier-Stokes equations and discretized by finite volume method. We classify the different cavitating regimes in which appear at the hydrofoil according to the values of cavitation number, surface roughness and angle of attack. Special attention is paid to the following topics: cavity length, void fraction and lift breakdown.

• Journal of Welding and Joining
• /
• v.15 no.5
• /
• pp.84-91
• /
• 1997

The metal transfer phenomenon of the pulsed-GMAW is simulated by formulating the electromagnetic force incorporated with the Volume of Fluid algorithm. The free surface profiles, pressure and velocity distributions within the drop are computed numerically. Axial velocity and acceleration generated during peak current period are found to have a significant effect on drop detachment. Therefore, the accelerated inertia force becomes one of important factors affecting metal transfer in the pulsed-GMAW. When the pulse current parameters are selected properly, the molten drop is detached just after current pulse, and the operating range of the pulsing frequency increases with higher peak current and duty cycle. Calculated operating ranges show reasonably good agreements with the available experimental data.

• Journal of Magnetics
• /
• v.20 no.2
• /
• pp.110-113
• /
• 2015

Two kinds of magnetic fluids with different volume fractions are symmetrically filled into the W0.9 photonic crystal waveguide structure. The 2D plane-wave expansion method is used to investigate the slow light properties numerically. The constant group index criterion is employed to evaluate the slow light performance. The wavelength bandwidth ${\Delta}{\lambda}$ centering at ${\lambda}_0=1550nm$ varies from 32.4 to 44.2 nm when the magnetic field factor ${\alpha}_{\parallel}$ changes from 0 to 1. And the corresponding normalized delay bandwidth product can be tuned from 0.221 to 0.258. For comparison and optimization, two infiltration cases are investigated and the more advantageous infiltration scheme is found.

• 유체기계공업학회:학술대회논문집
• /
• 1999.12a
• /
• pp.179-186
• /
• 1999

Three-dimensional flow analysis is implemented to investigate the flow through transonic axial-flow compressor rotor(NASA R67), and to evaluate the performances of k-$\epsilon$ and Baldwin-Lomax turbulence models. A finite volume method is used for spatial discretization. And, the equations are solved implicitly in time with the use of approximate factorization. Upwind difference scheme is used for inviscid terms, but viscous terms are centrally differenced. The flux-difference-splitting of Roe is used to obtain fluxes at the cell faces. Numerical analysis is performed near peak efficiency and near stall. And, the results are compared with the experimental data for NASA R67 rotor. Blade-to-Blade Mach number distributions are compared to confirm the accuracy of the code. From the results, we conclude that k-$\epsilon$ model is better for the calculation of flow rate and efficiency than Baldwin-Lomax model. But, the predictions for Mach number and shock structure are almost same.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.379-392
• /
• 2002

Numerical flow computations around an aeroelastic 3D square cylinder immersed in the turbulent boundary layer are shown. Present computational code can be characterized by three numerical aspects which are 1) the method of artificial compressibility is adopted for the incompressible flow computations, 2) the domain decomposition technique is used to get better grid point distributions, and 3) to achieve the conservation law both in time and space when the flow is computed a with moving and transformed grid, the time derivatives of metrics are evaluated using the time-and-space volume. To provide time-dependant inflow boundary conditions satisfying prescribed time-averaged velocity profiles, a convenient way for generating inflow turbulence is proposed. The square cylinder is modeled as a 4-lumped-mass system and it vibrates with two-degree of freedom of heaving motion. Those blocks which surround the cylinder are deformed according to the cylinder's motion. Vigorous oscillations occur as the vortex shedding frequency approaches cylinder's natural frequencies.