• Nuclear Engineering and Technology
• /
• v.41 no.7
• /
• pp.885-892
• /
• 2009

This manuscript will discuss a numerical method where the six equations of two-phase flow, the solid heat conduction equations, and the two equations that describe neutron diffusion and precursor concentration are solved together in a tightly coupled, nonlinear fashion for a simplified model of a nuclear reactor core. This approach has two important advantages. The first advantage is a higher level of accuracy. Because the equations are solved together in a single nonlinear system, the solution is more accurate than the traditional "operator split" approach where the two-phase flow equations are solved first, the heat conduction is solved second and the neutron diffusion is solved third, limiting the temporal accuracy to $1^{st}$ order because the nonlinear coupling between the physics is handled explicitly. The second advantage of the method described in this manuscript is that the time step control in the fully implicit system can be based on the timescale of the solution rather than a stability-based time step restriction like the material Courant limit required of operator-split methods. In this work, a pilot code was used which employs this tightly coupled, fully implicit method to simulate a reactor core. Results are presented from a simulated control rod movement which show $2^{nd}$ order accuracy in time. Also described in this paper is a simulated rod ejection demonstrating how the fastest timescale of the problem can change between the state variables of neutronics, conduction and two-phase flow during the course of a transient.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.97-106
• /
• 1989

This paper is on the modeling of two-dimensional groundwater flow, which is the first step of the development of Dynamic System Model for groundwater flow and pollutant transport in subsurface porous media. The particular features of the model are its versatility and flexibility to deal with as many real-world problems as possible. Points as well as distributed sources/sinks are included to represent recharges/pumping and rainfall infiltrations. All sources/sinks can be transient or steady state. Prescribed hydraulic head on the Dirichlet boundaries and fluxes on Neumann or Cauchy boundaries can be time-dependent or constant. Sources/sinks strength over each element and node, hydraulic head at each Dirichlet boundary node and flux at each boundary segment can vary independently of each other. Either completely confined or completely unconfined aquifers, or partially confined and partially unconfined aquifers can be dealt with effectively. Discretization of a compound region with very irregular curved boundaries is made easy by including both quadrilateral and triangular elements in the formulation. Large-field problems can be solved efficiently by including a pointwise iterative solution strategy as an optional alternative to the direct elimination solution methed for the matrix equation approximating the partial differential equation of groundwater flow. The model also includes transient flow through confining leaky aquifers lying above and/or below the aquifer of interest. The model is verified against three simple cases to which analytical solutions are available. The groundwater flow model shall be combined with the model of pollutant transport in subsurface porous media. Then the combined model, with the applications of the Eigenvalue technique and the Dynamic system theory, shall be improved to the Dynamic System Model which can simulate the real groundwater flow and the pollutant transport accurately and effectively for the analyses and predictions.

• Journal of the Korean Solar Energy Society
• /
• v.35 no.1
• /
• pp.35-43
• /
• 2015

A numerical simulation on the wake flow of a wind turbine which is a scaled version of a multi-megawatt wind turbine has been performed. Two different inlet conditions of averaged wind speed including one below and one above the rated wind speed were used in the simulation. Steady-state pitch angles of the blade associated with the two averaged wind speeds were imposed for the simulation. The steady state analysis based on the Reynolds averaged Navier-Stokes equations with the method of frame motion were used for the simulation to find the torque of the rotor and the wake field behind the wind turbine. The simulation results were compared with the results obtained from the wind tunnel testing. From comparisons, it was found that the simulation results on the turbine power are pretty close to the experimental values. Also, the wake results were relatively close to the experimental results but there existed some discrepancy in the shape of velocity deficit. The reason for the discrepancy is considered due to the steady state solution with the frame motion method used in the simulation. However, the method is considered useful for solutions with much reduced calculation time and reasonably good accuracy compared to the transient analysis.

• Geomechanics and Engineering
• /
• v.3 no.1
• /
• pp.1-16
• /
• 2011

Piles passing through sloping liquefiable deposits are prone to lateral loading if these deposits liquefy and flow during earthquakes. These lateral loads caused by the relative soil-pile movement will induce bending in the piles and may result in failure of the piles or excessive pile-head displacement. Whilst the weak nature of the flowing liquefied soil would suggest that only small loads would be exerted on the piles, it is known from case histories that piles do fail owing to the influence of laterally spreading soils. It will be shown, based on dynamic centrifuge test data, that dilatant behaviour of soil close to the pile is the major cause of these considerable transient lateral loads which are transferred to the pile. This paper reports the results of geotechnical centrifuge tests in which models of gently sloping liquefiable sand with pile foundations passing through them were subjected to earthquake excitation. The soil close to the pile was instrumented with pore-pressure transducers and contact stress cells in order to monitor the interaction between soil and pile and to track the soil stress state both upslope and downslope of the pile. The presence of instrumentation measuring pore-pressure and lateral stress close to the pile in the research described in this paper gives the opportunity to better study the soil stress state close to the pile and to compare the loads measured as being applied to the piles by the laterally spreading soils with those suggested by the JRA design code. This test data shows that lateral stresses much greater than one might expect from calculations based on the residual strength of liquefied soil may be applied to piles in flowing liquefied slopes owing to the dilative behaviour of the liquefied soil. It is shown at least for the particular geometry studied that the current JRA design code can be un-conservative by a factor of three for these dilation-affected transient lateral loads.

• Journal of ILASS-Korea
• /
• v.25 no.4
• /
• pp.178-187
• /
• 2020

The current emission regulations, US Tier-4 and EU Stage-V, are only able to satisfy the regulations when all currently mass-produced emission reduction technologies such as EGR, DOC, DPF, and SCR are applied. Therefore, in this study, for the application of the Urea-SCR system to non-road diesel engines, the database was established by measuring the NO, NO2 concentration and calculating the NO2/NOx ratio based on the catalyst temperature and exhaust mass flow rate. Also, based on the measured NO2/NOx ratio data, a mathematical model was proposed to predict the NO2/NOx ratio at SCR catalyst, and the suitability of the model was verified through steady-state and transient mode. As a result of comparing the NO2/NOx ratio measured at the DOC outlet under the steady-state condition to two model values separately, the R2 was 0.9811 for the 3D map model and 0.9303 for the mathematical model. And in the case of the NO2/NOx ratio measured at the DPF outlet, the R2 was 0.9797 for the 3D map model and 0.935 for the mathematical model. It was confirmed that the R2 with the model value of the 3D Map of the mathematical model in the transient mode is 0.957, which shows high reliability.

• International Journal of Fluid Machinery and Systems
• /
• v.1 no.1
• /
• pp.1-9
• /
• 2008

In the first part of the paper, Computational Fluid Dynamics analysis of the combusting flow within a high-swirl lean premixed gas turbine combustor and over the $1^{st}$ row nozzle guide vanes is presented. In this analysis, the focus of the investigation is the fluid dynamics at the combustor/turbine interface and its impact on the turbine. The predictions show the existence of a highly-rotating vortex core in the combustor, which is in strong interaction with the turbine nozzle guide vanes. This has been observed to be in agreement with the temperature indicated by thermal paint observations. The results suggest that swirling flow vortex core transition phenomena play a very important role in gas turbine combustors with modern lean-premixed dry low emissions technology. As the predictability of vortex core transition phenomena has not yet been investigated sufficiently, a fundamental validation study has been initiated, with the aim of validating the predictive capability of currently-available modelling procedures for turbulent swirling flows near the sub/supercritical vortex core transition. In the second part of the paper, results are presented which analyse such transitional turbulent swirling flows in two different laboratory water test rigs. It has been observed that turbulent swirling flows of interest are dominated by low-frequency transient motion of coherent structures, which cannot be adequately simulated within the framework of steady-state RANS turbulence modelling approaches. It has been found that useful results can be obtained only by modelling strategies which resolve the three-dimensional, transient motion of coherent structures, and do not assume a scalar turbulent viscosity at all scales. These models include RSM based URANS procedures as well as LES and DES approaches.

• Korea-Australia Rheology Journal
• /
• v.17 no.3
• /
• pp.125-130
• /
• 2005

Coalescence process of binary immiscible fluid mixtures having bimodal size distributions, prepared by mixing two pre-sheared samples at different shear rates, ${\gamma}_{pre1}\;and\;{\gamma}_{pre2}$, under shear flow at a final shear rate, ${\gamma}_f$, are examined by transient shear stress measurements and microscopic observations in comparison with the results for simply pre-sheared samples having narrow size distributions (unimodal distribution samples). Component fluids are a silicone oil (PDMS) and a hydrocarbon-formaldehyde resin (Genelite) and their viscosities are 14.1 and 21.0 $pa{\cdot}sec$ at room temperature $(ca.\;20^{\circ}C)$, respectively. The weight ratio of PDMS: Genelite was 7:3. Three cases, $({\gamma}_{pre1}=7.2sec^{-1},\;{\gamma}_{pre2}=12.0sec^{-1}\;and\;{\gamma}_f=2.4sec^{-1}),\;({\gamma}_{pre1}=0.8sec^{-1},\;{\gamma}_{pre2}=4.0sec^{-1}\;and\;{\gamma}_f=2.4sec^{-1}),\;and\;({\gamma}_{pre1}=7.2sec^{-1},\;{\gamma}_{pre2}=12.0^sec^{-1}\;and\;{\gamma}_f=7.2sec^{-1})$ the first case, transient shear stress did not show any significant difference but domains larger than the initial state are observed at short times. In the latter cases, there exist undershoot of shear stress, reflecting existence of deformed large domains, which is confirmed by the direct observation. It is concluded that coalescence between large and small domains more frequently occur than coalescence between the domains with similar size in the bimodal distribution samples.

• Nuclear Engineering and Technology
• /
• v.55 no.5
• /
• pp.1651-1664
• /
• 2023

It is very difficult to capture the multi-dimensional phenomena such as asymmetric flow and temperature distributions with the one-dimensional (1D) model, obviously, due to its inherent limitation. In order to overcome such a limitation of the 1D representation, many state-of-the-art system codes have equipped a three-dimensional (3D) component for multi-dimensional analysis capability. In this study, a standard multi-dimensional analysis model of APR1400 (Advanced Power Reactor 1400) has been developed using TRACE (TRAC/RELAP Advanced Computational Engine). The entire reactor pressure vessel (RPV) of APR1400 has been modeled using a single 3D component. The fuels in the reactor core have been described with detailed and coarse representations, respectively, to figure out the impact of the fuel description. Using both 3D RPV models, a comparative analysis has been performed postulating a double-ended guillotine break at a cold leg. Based on the results of comparative analysis, it is revealed that both models show no significant difference in general plant behavior and the model with coarse fuel model could be used for faster transient analysis without reactor kinetics coupling. The analysis indicates that the asymmetric temperature and flow distributions are captured during the transient, and such nonuniform distributions contribute to asymmetric quenching behaviors during blowdown and reflood phases. Such asymmetries are directly connected to the figure of merits in the LBLOCA analysis. Therefore, it is recommended to employ a multi-dimensional RPV model with a detailed fuel description for a realistic safety analysis with the consideration of the spatial configuration of the reactor core.

• Wind and Structures
• /
• v.16 no.1
• /
• pp.1-23
• /
• 2013

The wind-induced transient response of internal pressure following the creation of a sudden dominant opening during the occurrence of high external pressure, in low-rise residential and industrial buildings was numerically investigated. The values of the ill-defined parameters namely the flow contraction coefficient, loss coefficient and the effective slug length were calibrated by matching the analytical response with the computational fluid dynamics predictions. The effect of a sudden i.e., "instantaneously created" windward opening in the Texas Technical University (TTU) test building envelope was studied for two different envelope flexibility-leakage combinations namely: (1) a quasi-statically flexible and non-porous envelope and (2) a quasi-statically flexible and porous envelope. The responses forced by creating the openings at different time leads/lags with respect to the occurrence of the peak external pressure showed that for cases where the openings are created in close temporal proximity to the peak pressure, the transient overshoot values of internal pressure could be higher than the peak values of internal pressure in the pre-sequent or subsequent resonant response. In addition, the influence of time taken for opening creation on the level of overshoot was also investigated for the TTU building for the two different envelope characteristics. Non-dimensional overshoot factors are presented for a variety of cavity volume-opening area combinations for (1) buildings with rigid/quasi-statically flexible non-porous envelope, and (2) buildings with rigid/quasi-statically flexible and porous envelope (representing most low rise residential and industrial buildings). While the factors appear slightly on the high side due to conservative assumptions made in the analysis, a careful consideration regarding the implication of the timing and magnitude of such overshoots during strong gusts, in relation to the steady state internal pressure response in cyclonic regions, is warranted.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.6
• /
• pp.734-743
• /
• 1999

The local Nusselt numbers have been measured for a round turbulent jet impinging on the concave surface with and without rib. Liquid crystal/transient method was used to determine the Nusselt number distributions along the surface. The temperature on the surface was measured using liquid crystal and a digital color image processing system. The experiments were made for the jet Reynolds number (Re) 23,000, the dimensionless nozzle-to-surface distance (L/d) from 4 to 10, the dimensionless surface curvature (d/D) 0.056, and the rib type (height ($d_1$) 0.2 cm, pitch (p) from 1.2 to 3.2 cm). It was founded that only when $L/d{\geq}6$, the average Nusselt numbers on the concave surface with rib are higher than those without rib, mainly due to an increase in the turbulent intensity caused by the effect of rib attached to the wall surface. It was realized that the rib attached to the concave surface may no longer enhance the heat transfer rate or even lowers it depending on the rib type and flow conditions. In addition, the results by the steady-state method using the gold-film Intrex were in good agreement with those by the transient shroud method.