• International Journal of Vascular Biomedical Engineering
• /
• v.2 no.1
• /
• pp.17-24
• /
• 2004

Tumor angiogenesis was simulated using a two-dimensional computational model. The equation that governed angiogenesis comprised a tumor angiogenesis factor (TAF) conservation equation in time and space, which was solved numerically using the Galerkin finite element method. The time derivative in the equation was approximated by a forward Euler scheme. A stochastic process model was used to simulate vessel formation and vessel elongation towards a paracrine site, i.e., tumor-secreted basic fibroblast growth factor (bFGF). In this study, we assumed a two-dimensional model that represented a thin (1.0 mm) slice of the tumor. The growth of the tumor over time was modeled according to the dynamic value of bFGF secreted within the tumor. The data used for the model were based on a previously reported model of a brain tumor in which four distinct stages (namely multicellular spherical, first detectable lesion, diagnosis, and death of the virtual patient) were modeled. In our study, computation was not continued beyond the 'diagnosis' time point to avoid the computational complexity of analyzing numerous vascular branches. The numerical solutions revealed that no bFGF remained within the region in which vessels developed, owing to the uptake of bFGF by endothelial cells. Consequently, a sharp, declining gradient of bFGF existed near the surface of the tumor. The vascular architecture developed numerous branches close to the tumor surface (the brush-border effect). Asymmetrical tumor growth was associated with a greater degree of branching at the tumor surface.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.1
• /
• pp.88-99
• /
• 1996

The explicit scheme for finite element analysis of sheet metal forming problems has been widely used for providing practical solutions since it improves the convergency problem, memory size and computational time especially for the case of complicated geometry and large element number. The explicit schemes in general use are based on the elastic-plastic modeling of material requiring large computataion time. In the present work, a basic formulation for rigid-plastic explicit finite element analysis of plain strain sheet metal forming problems has been proposed. The effect of some basic parameters involved in the dynamic analysis has been studied in detail. Thus, the effective ranges of parameters have been proposed for numerical simultion by the rigid-plastic explicit finite element method. A direct trial-and-error method is introduced to treat contact and friction. In computation, sheet material is assumed to possess normal anisotropy and rigid-plastic workhardening characteristics. In order to show the validity and effectiveness of the proposed explicit scheme, computations are carried out for cylindrical punch stretching and the computational results are compared with those by the implicit scheme as well as with a commercial code. The proposed rigid-plastic exlicit finite element method can be used as a robust and efficient computational method for analysis of sheet metal forming.

• Journal of Soil and Groundwater Environment
• /
• v.17 no.5
• /
• pp.56-67
• /
• 2012

In the present study, an enhanced subsurface prediction algorithm based on a non-parametric geostatistical model and a history matching technique through Gibbs sampler is developed and the iterative prediction improvement procedure is proposed. The developed model is applied to a simple two-dimensional synthetic case where domain is composed of three different hydrogeologic media with $500m{\times}40m$ scale. In the application, it is assumed that there are 4 independent pumping tests performed at different vertical interval and the history curves are acquired through numerical modeling. With two hypothetical borehole information and pumping test data, the proposed prediction model is applied iteratively and continuous improvements of the predictions with reduced uncertainties of the media distribution are observed. From the results and the qualitative/quantitative analysis, it is concluded that the proposed model is good for the subsurface prediction improvements where the history data is available as a supportive information. Once the proposed model be a matured technique, it is believed that the model can be applied to many groundwater, geothermal, gas and oil problems with conventional fluid flow simulators. However, the overall development is still in its preliminary step and further considerations needs to be incorporated to be a viable and practical prediction technique including multi-dimensional verifications, global optimization, etc. which have not been resolved in the present study.

• Geophysics and Geophysical Exploration
• /
• v.24 no.2
• /
• pp.35-44
• /
• 2021

As the execution speed of Python is slower than those of other programming languages (e.g., C, C++, and FORTRAN), Python is not considered to be efficient for writing numerical geodynamic code that requires numerous iterations. Recently, many computational techniques, such as the Just-In-Time (JIT) compiler, have been developed to enhance the calculation speed of Python. Here, we developed two-dimensional (2D) numerical geodynamic code that was optimized for the JIT compiler, based on Python. Our code simulates mantle convection by combining the Particle-In-Cell (PIC) scheme and the finite element method (FEM), which are both commonly used in geodynamic modeling. We benchmarked well-known mantle convection problems to evaluate the reliability of our code, which confirmed that the root mean square velocity and Nusselt number obtained from our numerical modeling were consistent with those of the mantle convection problems. The matrix assembly and PIC processes in our code, when run with the JIT compiler, successfully achieved a speed-up 30× and 258× faster than without the JIT compiler, respectively. Our Python-based FEM-PIC code shows the high potential of Python for geodynamic modeling cases that require complex computations.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.24 no.7
• /
• pp.939-944
• /
• 2018

Clean fire extinguishing agents refer to chemical that can replace Halon 1211 and Halon 1310 according to the Montreal Protocol fermented to protect the Earth's ozone layer. In Korea and abroad, system standardization and performance evaluation of clean fire extinguishing agents are being carried out. This paper proposes an optimal nozzle shape by modeling and numerical analysis of various nozzle shapes based on general clean fire extinguishing system. The ejection speed of the nozzle can be improved by studying three - dimensional modeling of the nozzle for two shapes, Type A and B. Flow analysis was performed on the two types of nozzles and the gas velocity and pressure distribution were measured with different nozzle diameters. It was confirmed that the jetting speed was changed at the nozzle outlet according to the number and diameter of the nozzle holes. The flow rate increased with increasing the pressure regardless of the nozzle hole diameter. Based on the results obtained from the experiment, the K-factor value was deduced. Finally, a nozzle with a 12-hole structure with a 5-mm nozzle hole was proposed.

• Journal of the Korean Geotechnical Society
• /
• v.39 no.5
• /
• pp.27-40
• /
• 2023

A pile is a structural element that is used to transfer external loads from superstructures and has been widely utilized in construction fields all over the world. The method of installing a pile into the ground should be selected based on geotechnical conditions, location, site status, environmental factors, and construction costs, among others. It can be divided into two types: direct hammering and preboring. The direct hammering method installs a pile into the bearing layer, such as rock, using a few types of hammer, generating a considerable amount of pile driving-induced vibration. The vibration from pile driving influences adjacent structures and the ground; therefore, quantitatively investigating the effects of vibration is inevitably required. In this study, two-dimensional dynamic numerical modeling and analysis are performed using the finite difference method to investigate the influence on the adjacent slope, including temporary supporting system. Time-dependent loading induced by pile driving is estimated and used in the numerical analysis. Consequently, large surface displacement is estimated due to surface waves and less wave deflection, and refraction at the surface. The total displacement decreases with the increase of the distance from the source. However, lateral displacement at the top of the slope shows a larger value than vertical displacement, and the overall displacement tends to be concentrated near the face of the slope.

• Geomechanics and Engineering
• /
• v.16 no.3
• /
• pp.295-308
• /
• 2018

The excavation of twin tunnels is a process that destabilizes the ground. The stability of the tunnel lining, the control of ground displacements around the tunnel resulting from each excavation and the interaction between them must be controlled. This paper provides a new approach for replacing the costly 3D analyses with the equivalent 2D analyses that closely reflects the in-situ measurements when excavating twin tunnels. The modeling was performed in two dimensions using the FLAC2D finite difference code. The three-dimensional effect of excavation is taken into account through the deconfinement rate ${\lambda}$ of the soil surrounding the excavation by applying the convergence-confinement method. A comparison between settlements derived by the proposed 2D analysis and the settlements measured in a real project in Algeria shows an acceptable agreement. Also, this paper reports the investigation into the changes in deformations on tunnel linings and surface settlements which may be expected if the twin tunnels of T4 El-Harouche Skikda were constructed with a tunneling machine. Special attention was paid to the influence of the excavation phase shift distance between the two mechanized tunnel faces. It is revealed that the ground movements and the lining deformations during tunnel excavation depend on the distance between the tunnels' axis and the excavation phase shift.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.3 no.1
• /
• pp.27-36
• /
• 2011

In the proposed paper numerical calculations are carried out using two versions of a three-dimensional, timedomain panel method developed by the group of Prof. P. Sclavounos at MIT, i.e. the linear code SWAN2, enabling optionally the use of the instantaneous non-linear Froude-Krylov and hydrostatic forces and the fully non-linear SWAN4. The analytical results are compared with experimental results for three hull forms with increasing geometrical complexity, the Series 60, a reefer vessel with stern bulb and a modern fast ROPAX hull form with hollow bottom in the stern region. The details of the geometrical modeling of the hull forms are discussed. In addition, since SWAN4 does not support transom sterns, only the two versions of SWAN2 were evaluated over experimental results for the parent hull form of the NTUA double-chine, wide-transom, high-speed monohull series. The effect of speed on the numerical predictions was investigated. It is concluded that both versions of SWAN2 the linear and the one with the non-linear Froude-Krylov and hydrostatic forces provide a more robust tool for prediction of the dynamic response of the vessels than the non-linear SWAN4 code. In general, their results are close to what was expected on the basis of experience. Furthermore, the use of the option of non-linear Froude-Krylov and hydrostatic forces is beneficial for the accuracy of the predictions. The content of the paper is based on the Diploma thesis of the second author, supervised by the first one and further refined by the third one.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.11a
• /
• pp.729-732
• /
• 2010

The present study has been motivated by the development of a reliable numerical methodology for simulation of kerosene/LOx coaxial swirl injectors. To deal with thermodynamic non-ideality and anomalies of transport properties pronounced at supercritical pressures, a set of subroutine libraries has been constructed based on the cubic equations of state, and applied to an existing flamelet analysis code. For computational efficiency, two-dimensional axisymmetric RANS formulation with swirl was adopted and validated successfully against an isothermal coaxial swirling jet. For the actual problem with high pressure combustion, however, numerical results show that the RANS models yield excessive production of turbulence probably due to high density gradient magnitude in the vicinity of mixing layer of swirling film flow, and imply strongly further improvement of the turbulence models.

• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1059-1069
• /
• 2023

Dam-break flow occurs when an elevated dam suddenly collapses, resulting in the catastrophic release of rapid and uncontrolled impounded water. This study compares laminar and turbulent closure models for simulating three-dimensional dam-break flows using OpenFOAM. The Reynolds-Averaged Navier-Stokes (RANS) model, specifically the k-ε model, is employed to capture turbulent dissipation. Two scenarios are evaluated based on a laboratory experiment and a modified multi-layered block obstacle scenario. Both models effectively represent dam-break flows, with the turbulent closure model reducing oscillations. However, excessive dissipation in turbulent models can underestimate water surface profiles. Improving numerical schemes and grid resolution enhances flow recreation, particularly near structures and during turbulence. Model stability is more significantly influenced by numerical schemes and grid refinement than the use of turbulence closure. The k-ε model's reliance on time-averaging processes poses challenges in representing dam-break profiles with pronounced discontinuities and unsteadiness. While simulating turbulence models requires extensive computational efforts, the performance improvement compared to laminar models is marginal. To achieve better representation, more advanced turbulence models like Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) are recommended, necessitating small spatial and time scales. This research provides insights into the applicability of different modeling approaches for simulating dam-break flows, emphasizing the importance of accurate representation near structures and during turbulence.