• Steel and Composite Structures
• /
• v.10 no.2
• /
• pp.129-149
• /
• 2010

The objective of this study is to formulate a general 3D material-structural analysis framework for the thermomechanical behavior of steel-concrete structures in a fire environment. The proposed analysis framework consists of three sequential modeling parts: fire dynamics simulation, heat transfer analysis, and a thermomechanical stress analysis of the structure. The first modeling part consists of applying the NIST (National Institute of Standards and Technology) Fire Dynamics Simulator (FDS) where coupled CFD (Computational Fluid Dynamics) with thermodynamics are combined to realistically model the fire progression within the steel-concrete structure. The goal is to generate the spatial-temporal (ST) solution variables (temperature, heat flux) on the surfaces of the structure. The FDS-ST solutions are generated in a discrete form. Continuous FDS-ST approximations are then developed to represent the temperature or heat-flux at any given time or point within the structure. An extensive numerical study is carried out to examine the best ST approximation functions that strike a balance between accuracy and simplicity. The second modeling part consists of a finite-element (FE) transient heat analysis of the structure using the continuous FDS-ST surface variables as prescribed thermal boundary conditions. The third modeling part is a thermomechanical FE structural analysis using both nonlinear material and geometry. The temperature history from the second modeling part is used at all nodal points. The ABAQUS (2003) FE code is used with external user subroutines for the second and third simulation parts in order to describe the specific heat temperature nonlinear dependency that drastically affects the transient thermal solution especially for concrete materials. User subroutines are also developed to apply the continuous FDS-ST surface nodal boundary conditions in the transient heat FE analysis. The proposed modeling framework is applied to predict the temperature and deflection of the well-documented third Cardington fire test.

• Proceedings of the KSR Conference
• /
• 2008.11b
• /
• pp.66-72
• /
• 2008

In this study, transient 3D numerical simulations were performed to analyze the characteristics of fire driven flow for smoke ventilation system operating conditions in the deeply underground subway station. The smoke flow patterns were compared and discussed under smoke fan operating mode and off mode in the platform. Soongsil Univ. station(line number 7)was chosen for simulation which was the one of the deepest underground subway stations in the each lines of Seoul. The geometry for model is 365m in length include railway, 23.5m for width, 47m for depth. Therefore 10,000,000 structured grids were used for fire simulation. The parallel computational method for fast calculation was employed to compute the heat and mass transfer eqn's with 6 CPUs(Intel 3.0GHz Dual CPU, 12Cores) of the linux clustering machine. The fire driven flow was simulated with using FDS code in which LES method was applied. The Heat release rate was 10MW and The Ultrafast model was applied for the growing model of the fire source.

• Nuclear Engineering and Technology
• /
• v.42 no.6
• /
• pp.620-635
• /
• 2010

In this paper we describe current activities on the project Multi-Scale Modeling and Analysis of convective boiling (MSMA), conducted jointly by the Paul Scherrer Institute (PSI) and the Swiss Nuclear Utilities (Swissnuclear). The long-term aim of the MSMA project is to formulate improved closure laws for Computational Fluid Dynamics (CFD) simulations for prediction of convective boiling and eventually of the Critical Heat Flux (CHF). As boiling is controlled by the competition of numerous phenomena at various length and time scales, a multi-scale approach is employed to tackle the problem at different scales. In the MSMA project, the scales on which we focus range from the CFD scale (macro-scale), bubble size scale (meso-scale), liquid micro-layer and triple interline scale (micro-scale), and molecular scale (nano-scale). The current focus of the project is on micro- and meso-scales modeling. The numerical framework comprises a highly efficient, parallel DNS solver, the PSI-BOIL code. The code has incorporated an Immersed Boundary Method (IBM) to tackle complex geometries. For simulation of meso-scales (bubbles), we use the Constrained Interpolation Profile method: Conservative Semi-Lagrangian $2^{nd}$ order (CIP-CSL2). The phase change is described either by applying conventional jump conditions at the interface, or by using the Phase Field (PF) approach. In this work, we present selected results for flows in complex geometry using the IBM, selected bubbly flow simulations using the CIP-CSL2 method and results for phase change using the PF approach. In the subsequent stage of the project, the importance of effects of nano-scale processes on the global boiling heat transfer will be evaluated. To validate the models, more experimental information will be needed in the future, so it is expected that the MSMA project will become the seed for a long-term, combined theoretical and experimental program.

• The Journal of Advanced Prosthodontics
• /
• v.11 no.3
• /
• pp.169-178
• /
• 2019

PURPOSE. While dental implants have displayed high success rates, poor mechanical fixation is a common complication, and their biomechanical response to occlusal loading remains poorly understood. This study aimed to develop and validate a computational model of a natural first premolar and a dental implant with matching crown morphology, and quantify their mechanical response to loading at the occlusal surface. MATERIALS AND METHODS. A finite-element model of the stomatognathic system comprising the mandible, first premolar and periodontal ligament (PDL) was developed based on a natural human tooth, and a model of a dental implant of identical occlusal geometry was also created. Occlusal loading was simulated using point forces applied at seven landmarks on each crown. Model predictions were validated using strain gauge measurements acquired during loading of matched physical models of the tooth and implant assemblies. RESULTS. For the natural tooth, the maximum vonMises stress (6.4 MPa) and maximal principal strains at the mandible ($1.8m{\varepsilon}$, $-1.7m{\varepsilon}$) were lower than those observed at the prosthetic tooth (12.5 MPa, $3.2m{\varepsilon}$, and $-4.4m{\varepsilon}$, respectively). As occlusal load was applied more bucally relative to the tooth central axis, stress and strain magnitudes increased. CONCLUSION. Occlusal loading of the natural tooth results in lower stress-strain magnitudes in the underlying alveolar bone than those associated with a dental implant of matched occlusal anatomy. The PDL may function to mitigate axial and bending stress intensities resulting from off-centered occlusal loads. The findings may be useful in dental implant design, restoration material selection, and surgical planning.

• Wind and Structures
• /
• v.31 no.6
• /
• pp.533-548
• /
• 2020

Despite the current technologic developments, failures in existent tensile fabric structures (TFS) subjected to wind do happen. However, design pressure coefficients are only obtained for large projects. Moreover, studies on TFSs with realistic supporting frames, comparing static and dynamic analyses and discussing the design implications, are lacking. In this study, fluid-Structure analyses of a TFS supported by masts and inclined cables, by subjecting it to different wind speeds, are carried out, to gain more understanding in the above-referred aspects. Wind-induced stresses in the fabric and axial forces in masts and cables are assessed for a hypar by using computational fluid dynamics. Comparisons are carried out versus an equivalent static analysis and also versus loadings deemed representative for design. The procedure includes the so-called form-finding, a finite element formulation for the TFS and the fluid formulation. The selected structure is deemed realistic, since the supporting frame is included and the shape and geometry of the TFS are not uncommon. It is found that by carrying out an equivalent static analysis with the determined pressure coefficients, differences of up to 24% for stresses in the fabric, 5.4% for the compressive force in the masts and 21% for the tensile force in the cables are found with respect to results of the dynamic analysis. If wind loads commonly considered for design are used, significant differences are also found, specially for the reactions at the supporting frame. The results in this study can be used as an aid by designers and researchers.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.24 no.2
• /
• pp.211-218
• /
• 2021

This paper deals with the computational work to characterize the formation and pinching of a plasma in an X-pinch configuration. A resistive magnetohydrodynamic model of a single fluid and two temperature is adopted assuming a hollow conical structure in the (r,z) domain. The model includes the thermodynamic parameter of tungsten from the corrected Thomas-Fermi EOS(equation of state), determining the average ionization charge, pressure, and internal energy. The transport coefficients, resistivity and thermal conductivity, are obtained by the corrected Lee & More model and a simple radiation loss rate by recombination process is considered in the simulation. The simulation demonstrated the formation of a core-corona plasma and intense compression process near the central region which agrees with the experimental observation in the X-pinch device at Seoul National University. In addition, it confirmed the increase in radiation loss rate with the density and temperature of the core plasma.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.13 no.1
• /
• pp.236-245
• /
• 2021

In this study, direct numerical and large eddy simulations of transitional flows around studs were conducted to investigate the effectiveness of turbulence stimulators at very low speeds for the minimum propulsion power condition of four knots. For simplicity, the studs were assumed to be installed on a flat plate, while the wake was observed up to 0.23 m downstream behind the second stud. For applicability to a model ship, we also studied the flow characteristics behind the first and second studs installed on a curved plate, which was designed to describe the geometry of a bulbous bow. A laminar-to-turbulent transition was observed in the wake at Re_D ≥ 921 (U_∞≥0.290 m/s), and the wall shear stress at Re_D = 1162 (U_∞ = 0.366 m/s) in the second wake was similar to that of the fully developed turbulent boundary layer after a laminar-to-turbulent transition in the first wake. At Re_D = 581 (U_∞ = 0.183 m/s), no turbulence was stimulated in the wake behind the first and second studs on the flat plate, while a cluster of vortical structures was observed in the first wake over the curved plate. However, a cluster of vortical structures was revealed to be generated by the reattachment process of the separated shear layer, which was disturbed by the first stud rather than directly initiated by the first stud. It was quite different from a typical process of transition, which was observed at relatively high Re_D that the spanwise scope of the turbulent vortical structures expanded gradually as it went downstream.

• Korean Chemical Engineering Research
• /
• v.52 no.6
• /
• pp.826-833
• /
• 2014

Driven by both environmental and economic reasons, the development of small to medium scale GTL(gas-to-liquid) process for offshore applications and for utilizing other stranded or associated gas has recently been studied increasingly. Microchannel GTL reactors have been prefrered over the conventional GTL reactors for such applications, due to its compactness, and additional advantages of small heat and mass transfer distance desired for high heat transfer performance and reactor conversion. In this work, multi-microchannel reactor was simulated by using commercial CFD code, ANSYS FLUENT, to study the geometric effect of the microchannels on the heat transfer phenomena. A heat generation curve was first calculated by modeling a Fischer-Tropsch reaction in a single-microchannel reactor model using Matlab-ASPEN integration platform. The calculated heat generation curve was implemented to the CFD model. Four design variables based on the microchannel geometry namely coolant channel width, coolant channel height, coolant channel to process channel distance, and coolant channel to coolant channel distance, were selected for calculating three dependent variables namely, heat flux, maximum temperature of coolant channel, and maximum temperature of process channel. The simulation results were visualized to understand the effects of the design variables on the dependent variables. Heat flux and maximum temperature of cooling channel and process channel were found to be increasing when coolant channel width and height were decreased. Coolant channel to process channel distance was found to have no effect on the heat transfer phenomena. Finally, total heat flux was found to be increasing and maximum coolant channel temperature to be decreasing when coolant channel to coolant channel distance was decreased. Using the qualitative trend revealed from the present study, an appropriate process channel and coolant channel geometry along with the distance between the adjacent channels can be recommended for a microchannel reactor that meet a desired reactor performance on heat transfer phenomena and hence reactor conversion of a Fischer-Tropsch microchannel reactor.

• Korean Journal of Ecology and Environment
• /
• v.42 no.1
• /
• pp.85-94
• /
• 2009

This study examines the operational effectiveness of a Convection Current Aeration System (CCAS) in reservoir. CCAS was run from June, 2008 when the thermocline begun forming in the reservoir. This paper reviews the influence of stratification, dissolved oxygen dynamics and temperature in the lake's natural state from June to October 2008. The survey was done on a week basis. Upwelling flow effects a radius of $7{\sim}10m$ at a surface directly and was irrelevant to the strength of thermocline. On the other hand, it was affected the number of working days, and strength of thermocline at vertical profiles of the reservoir. Longer CCAS run, the deeper was the vertical direct flow area. However it didn't break the thermocline during summer season of 2008. The operating efficiency of the CCAS in the reservoir depends on hydraulics and meteological conditions. Computational Fluid Dynamics (CFD) is a very useful tool for evaluating the operating efficiency of fluid dynamics. The geometry for CFD simulation consists of a cylindrical vessel 25 m radius and 40 m height. The CCAS is located in center of domain. The non-uniform tetrahedral meshes had a bulk of the geometry. The meshes ranged from the coarse to the very fine. This is attributed to the cold water flowing into the downcomer and rising, creating a horizontal flow to the top of the CCAS. The result of CFD demonstrate a closer agreement with surveyed data for temperature and flow velocity. Theoretical dispersion volume were calculated at 8m depth, 120 m diameter working for 30 days and 10 m depth, 130 m diameter working for 50 days.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.18 no.12
• /
• pp.1970-1982
• /
• 1993

This paper describes an efficient array algorithm for parallel computation of vector-radix two-dimensional (2-D) fast discrete cosine transform (VR-FCT), and its VLSI implementation. By mapping the 2-D VR-FCT onto a 2-D array of processing elements (PEs), the butterfly structure of the VR-FCT can be efficiently importanted with high concurrency and local communication geometry. The proposed array algorithm features architectural modularity, regularity and locality, so that it is very suitable for VLSI realization. Also, no transposition memory is required, which is invitable in the conventional row-column decomposition approach. It has the time complexity of O(N+Nnzp-log2N) for (N*N) 2-D DCT, where Nnzd is the number of non-zero digits in canonic-signed digit(CSD) code, By adopting the CSD arithmetic in circuit desine, the number of addition is reduced by about 30%, as compared to the 2`s complement arithmetic. The computational accuracy analysis for finite wordlength processing is presented. From simulation result, it is estimated that (8*8) 2-D DCT (with Nnzp=4) can be computed in about 0.88 sec at 50 MHz clock frequency, resulting in the throughput rate of about 72 Mega pixels per second.