• Journal of computational fluids engineering
• /
• v.14 no.2
• /
• pp.46-51
• /
• 2009

The slip effect from the molecular interaction between fluid particles and solid surface atoms plays a key role in microscale fluid transport and heat transfer since the relative importance of surface forces increases as the size of the system decreases to the microscale. There exist two models to describe the slip effect: the Maxwell slip model in which the slip correction is made on the basis of the degree of shear stress near the wall surface and the Langmuir slip model based on a theory of adsorption of gases on solids. In this study, as the first step towards developing a general purpose numerical code of the compressible Navier-Stokes equations for computational simulations of microscale fluid flow and heat transfer, two slip models are implemented into a finite element numerical code of a simplified equation. In addition, a pressure-driven gas flow in a microchannel is investigated by the numerical code in order to validate numerical results.

• Economic and Environmental Geology
• /
• v.54 no.2
• /
• pp.247-257
• /
• 2021

The study area has been contaminated with explosive materials and heavy metals for several decades. For the design of the pump and treat remediation method, groundwater flow before and during groundwater pumping in a contaminated aquifer system was simulated, calibrated, and predicted using a generalized multidimensional hydrological numerical model. A three-dimensional geologic formation model representing the geology, hydrogeology, and topography of the aquifer system was established. A steady-state numerical simulation with model calibration was performed to obtain initial steady-state spatial distributions of groundwater flow and groundwater table in the aquifer system before groundwater pumping, and its results were illustrated and analyzed. A series of transient-state numerical simulations were then performed during groundwater pumping with the four different pumping rates at a potential location of the pumping well. Its results are illustrated and analyzed to provide primary reference data for the pump and treat remediation method. The results of both steady-state and transient-state numerical simulations show that the spatial distribution and properties of the geologic media and the topography have significant effects on the groundwater flow and thus depression zone.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.2
• /
• pp.426-437
• /
• 1993

The flowfield of Weis-Fogh type ship's Propulsion is visualized by numerical simulations using the discrete vortex method and by the hydrogen bubble technique. The simulations are performed by assuming that the separations occur at the trailing edge of the wing. The streak lines and time lines are calculated by introducing the tracers at adequate intervals. They agree well with experimental results. The flowfield is unsteady and complex, but the properties of the flow are clarified by numerical and experimental visualization.

• Advances in nano research
• /
• v.13 no.6
• /
• pp.599-617
• /
• 2022

In the present research, for the first time, the vibrational as well as buckling characteristics of a three-layered curved nanobeam including a core made of functionally graded (FG) material and two layers of smart material-piezo-magneto-electric-resting on a Winkler Pasternak elastic foundation are examined. The displacement field for the nanobeam is chosen via Timoshenko beam theory. Also, the size dependency is taken into account by using nonlocal strain gradient theory, aka NSGT. Then, by employing Hamilton's principle, energy procedure, the governing equations together with the boundary conditions are achieved. The solution procedure is a numerical solution called generalized differential quadrature method, or GDQM. The accuracy and reliability of the formulation alongside solution method is examined by using other published articles. Lastly, the parameter which can alter and affect the buckling or vocational behavior of the curved nanobeam is investigated in details.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.5
• /
• pp.400-410
• /
• 2019

In this study, the numerical code for the 3D nonlinear dynamic analysis of an SLWR (Steel Lazy Wave Riser) was developed using the lumped mass line model in a FORTRAN environment. Because the lumped mass line model is an explicit method, there is no matrix operation. Thus, the numerical algorithm is simple and fast. In the lumped mass line model, the equations of motion for the riser were derived by applying the various forces acting on each node of the line. The applied forces at the node of the riser consisted of the tension, shear force due to the bending moment, gravitational force, buoyancy force, riser/ground contact force, and hydrodynamic force based on the Morison equation. Time integration was carried out using a Runge-Kutta fourth-order method, which is known to be stable and accurate. To validate the accuracy of the developed numerical code, simulations using the commercial software OrcaFlex were carried out simultaneously and compared with the results of the developed numerical code. To understand the nonlinear dynamic characteristics of an SLWR, dynamic simulations of SLWRs excited at the hang-off point and of SLWRs in regular waves were carried out. From the results of these dynamic simulations, the displacements at the maximum bending moments at important points of the design, like the hang-off point, sagging point, hogging points, and touch-down point, were observed and analyzed.

• Geomechanics and Engineering
• /
• v.23 no.3
• /
• pp.189-206
• /
• 2020

The accurate modeling of boundary conditions is important in simulations of the discrete element method (DEM) and can affect the numerical results significantly. In conventional triaxial compression (CTC) tests, the specimens are wrapped by flexible membranes allowing to deform freely. To accurately model the boundary conditions of CTC, new flexible boundary algorithms for 2D and 3D DEM simulations are proposed. The new algorithms are computationally efficient and easy to implement. Moreover, both horizontal and vertical component of confining pressure are considered in the 2D and 3D algorithms, which can ensure that the directions of confining pressure are always perpendicular to the specimen surfaces. Furthermore, the boundaries are continuous and closed in the new algorithms, which can prevent the escape of particles from the specimens. The effectiveness of the proposed algorithms is validated by biaxial and triaxial simulations of granular materials. The results show that the algorithms allow the boundaries to deform non-uniformly on the premise of maintaining high control accuracy of confining pressure. Meanwhile, the influences of different lateral boundary conditions on the numerical results are discussed. It is indicated that the flexible boundary is more appropriate for the models with large strain or significant localization than rigid boundary.

• Structural Engineering and Mechanics
• /
• v.61 no.6
• /
• pp.693-700
• /
• 2017

Concrete cover cracking due to the corrosion of steel reinforcing bars is one of the main causes of deterioration in Reinforced Concrete (RC) structures. The oxidation level of the bars causes varying levels of expansion. The rebar expansions could lead to through-thickness cracking of the concrete cover, where depending on the cracking characteristics, the service life of the structures would be affected. In this paper, the effect of geometrical and material parameters, i.e., concrete cover thickness, reinforcing bar diameter, and concrete tensile strength, on the required pressure for concrete cover cracking due to corrosion has been investigated through detailed numerical simulations. ABAQUS finite element software is employed as a modeling platform where the concrete cracking is simulated by means of eXtended Finite Element Method (XFEM). The accuracy of the numerical simulations is verified by comparing the numerical results with experimental data obtained from the literature. Using a previously proposed empirical equation and the numerical model, the time from corrosion initiation to the cover cracking is predicted and then compared to the respective experimental data. Finally, a parametric study is undertaken to determine the optimum ratio of the rebar diameter to the reinforcing bars spacing in order to avoid concrete cover delamination.

• 한국전산유체공학회:학술대회논문집
• /
• 1996.05a
• /
• pp.85-92
• /
• 1996

Numerical simulations are made for the three-dimensional flow around a wing in ground effect craft haying the complex geometry. A numerical tool is developed for the primary design of hull and wing shape of practical Wing-In-Ground effect(WIG) stop. The finite-difference method is utilized to descretize the governing equations and pressure field is obtained by using Marker-And-Cell(MAC) method. The air and water flows are simultaneously simulated in the time-marching solution procedure for the Navier-Stokes equation. The porosity technique and the density function are devised for the implementation of the three-dimensional body-boundary and the free-surface conditions, respectively. In this paper, a craft is modeled simply by three blocks containing a wing mounted on a main body horizontally, with the endplate. The numerical calculations of a WIG advancing in a calm water are performed and the WIG-generated wave profiles are also obtained. In the final paper, details of the numerical methods employed for the present study and calculated results are discussed.

• Transactions of Materials Processing
• /
• v.14 no.3 s.75
• /
• pp.224-232
• /
• 2005

The high-pressure die-casting is one of the most effective methods to produce a large amount of products in short cycle time. This process, however, has a problem that the gas porosity defect appears easily. The generation of gas porosity is known mainly due to the air entrapment during the injection stage. Most of numerical simulations for the molten metal flow pattern observations have done in the treating of one phase fluid flow but the gas-liquid interface is essentially multi- phase phenomenon. In this paper, the two-phase fluid flow numerical simulation methods have been adapted to predict the gas porosity generations in the molten metal. The accuracy and the usefulness of the new simulation module have been emphasized and verified through some comparison experiments.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.1
• /
• pp.1-10
• /
• 2010

In the present paper, a numerical method based on the constraint interpolation profile (CIP) method is applied for simulating two-dimensional violent sloshing problems. The free surface boundary value problem is considered as a multiphase problem which includes water and air. A stationary Cartesian grid system is adopted, and an interface capturing method is used to trace the shape of free surface profile. The CIP combined unified procedure (CCUP) scheme is applied for flow solver, and the tangent of hyperbola for interface capturing (THINC) scheme is used for interface capturing. Numerical simulations have been carried out for partially-filled 2D tanks under forced sway and roll motions at various filling depths and frequencies. The computational results are compared with experiments and/or the other numerical results to validate the present numerical method.