• Interaction and multiscale mechanics
• /
• v.1 no.1
• /
• pp.83-101
• /
• 2008

A derivation is given of two-scale models that are able to describe deformation and flow in a fluid-saturated and progressively fracturing porous medium. From the micromechanics of the flow in the cavity, identities are derived that couple the local momentum and the mass balances to the governing equations for a fluid-saturated porous medium, which are assumed to hold on the macroscopic scale. By exploiting the partition-of-unity property of the finite element shape functions, the position and direction of the fractures are independent from the underlying discretization. The finite element equations are derived for this two-scale approach and integrated over time. The resulting discrete equations are nonlinear due to the cohesive crack model and the nonlinearity of the coupling terms. A consistent linearization is given for use within a Newton-Raphson iterative procedure. Finally, examples are given to show the versatility and the efficiency of the approach.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.3
• /
• pp.22-29
• /
• 2004

Deformation of the porous media has influence on performance of a proton exchange membrane fuel cell (PEMFC). The stress distributions and deformation of the porous media are major factors for safe and efficient operation in the PEMFC. In this paper, nonlinear contact analysis of air plate and porous media is performed under a working condition to predict the performance characteristics of the air plates. Two kinds of models are suggested for this study. The first porous media model has nonlinear material properties. The second model has nonlinear material properties with contact condition between porous media and air plate. The numerical analysis results of the two models are somewhat different. It is shown that the nonlinear contact analysis is required for the design study of the PEMFC.

• Proceedings of the SAREK Conference
• /
• 2009.06a
• /
• pp.1159-1165
• /
• 2009

In the present study, the flow and heat transfer characteristics of a large plate heat exchanger are investigated numerically. The flow passages are very complicated due to the grooved corrugation patterns of the plate surface so that the detailed mesh and the large amount of the computation time have to be required in the numerical simulation for the conjugate heat transfer analysis. In order to accomplish the efficient and fast analysis of the heat transfer characteristics in the plate heat exchanger, a semimicroscopic method using the porous media model has been investigated numerically. The results showed that the characteristics of the heat transfer and pressure drop, which are respectively presented with Colburn j-factor and Fanning f-factor, are in a good agreement between the detailed mesh and the porous media model. The results of the present study could be applicable to the numerical analysis of entire flow passages in the large plate heat exchanger using porous media treatment.

• Nuclear Engineering and Technology
• /
• v.56 no.5
• /
• pp.1747-1753
• /
• 2024

The accurate determination of formation density and the physical properties of rocks is the most critical logging tasks which can be obtained using gamma-ray transport and detection tools. Though the simulation works published so far have considerably improved the knowledge of the parameters that govern the responses of the detectors in these tools, recent studies have found considerable differences between the results of using a conventional model of a homogeneous mixture of formation and fluid and an inhomogeneous fractured medium. It has increased concerns about the importance of the complexity of the model used for the medium in simulation works. In the present study, we have suggested two various models for the flow of the fluid in porous media and fractured rock to be used for logging purposes. For a typical gamma-gamma logging tool containing a ¹³⁷Cs source and two NaI detectors, simulated by using the MCNPX code, a simplified porous (SP) model in which the formation is filled with elongated rectangular cubes loaded with either mineral material or oil was investigated. In this model, the oil directly reaches the top of the medium and the connection between the pores is not guaranteed. In the other model, the medium is a large 3-D matrix of 1 cm³ randomly filled cubes. The designed algorithm to fill the matrix sites is so that this realistic random (RR) model provides the continuum growth of oil flow in various disordered directions and, therefore, fulfills the concerns about modeling the rock textures consist of extremely complex pore structures. For an arbitrary set of oil concentrations and various formation materials, the response of the detectors in the logging tool has been considered as a criterion to assess the effect of modeling for the distribution of pores in the formation on simulation studies. The results show that defining a RR model for describing heterogeneities of a porous medium does not effectively improve the prediction of the responses of logging tools. Taking into account the computational cost of the particle transport in the complex geometries in the Monte Carlo method, the SP model can be satisfactory for gamma-gamma logging purposes.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.2
• /
• pp.95-103
• /
• 2007

An analytical closed-form solution for wave propagation velocity and damping in saturated porous media is presented in this paper The fully coupled field model with compressible solid Brains and pore water were used to derive this solution. An engineering approach for the analysis of fully saturated porous media was adopted and closed-form solutions for one dimensional wave propagation in a homogeneous domain were derived. The solution is highly versatile in that it considers compression of the solid grains, compression of the pore water, deformation of the porous skeleton, and spatial damping and can be used to compute wavespeeds of first and second kind and damping coefficients in various geologic materials. This solution provides a means of analyzing the influence of material property variations on wavespeed and attenuation. In Part 2 of this work the theoretical solution is incorporated into the numerical code and the code is used in a parametric study on wave propagation velocity and damping.

• Coupled systems mechanics
• /
• v.4 no.1
• /
• pp.37-65
• /
• 2015

Modeling and simulation of mechanical response of infrastructure object, solids and structures, relies on the use of computational models to foretell the state of a physical system under conditions for which such computational model has not been validated. Verification and Validation (V&V) procedures are the primary means of assessing accuracy, building confidence and credibility in modeling and computational simulations of behavior of those infrastructure objects. Validation is the process of determining a degree to which a model is an accurate representation of the real world from the perspective of the intended uses of the model. It is mainly a physics issue and provides evidence that the correct model is solved (Oberkampf et al. 2002). Our primary interest is in modeling and simulating behavior of porous particulate media that is fully saturated with pore fluid, including cyclic mobility and liquefaction. Fully saturated soils undergoing dynamic shaking fall in this category. Verification modeling and simulation of fully saturated porous soils is addressed in more detail by (Tasiopoulou et al. 2014), and in this paper we address validation. A set of centrifuge experiments is used for this purpose. Discussion is provided assessing the effects of scaling laws on centrifuge experiments and their influence on the validation. Available validation test are reviewed in view of first and second order phenomena and their importance to validation. For example, dynamics behavior of the system, following the dynamic time, and dissipation of the pore fluid pressures, following diffusion time, are not happening in the same time scale and those discrepancies are discussed. Laboratory tests, performed on soil that is used in centrifuge experiments, were used to calibrate material models that are then used in a validation process. Number of physical and numerical examples are used for validation and to illustrate presented discussion. In particular, it is shown that for the most part, numerical prediction of behavior, using laboratory test data to calibrate soil material model, prior to centrifuge experiments, can be validated using scaled tests. There are, of course, discrepancies, sources of which are analyzed and discussed.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.4
• /
• pp.361-368
• /
• 2010

In this paper, the parallel algorithm using MPI(Message-Passing Interface) library is introduced in order to improve numerical efficiency for the staggered method for nearly incompressible and impermeable porous media which was introduced by Park and Tak(2010). The porous media theory and the staggered method are also briefly introduced in this paper. Moreover, we account for MPI library for blocking, non-blocking, and collective communication, and propose combined the staggered method with the blocking and nonblocking MPI library. And then, we present how to allocate CPUs on the staggered method and the MPI library, which is related with the numerical efficiency in order to solve unknown variables on nearly incompressible and impermeable porous media. Finally, the results comparing serial solution with parallel solution are verified by 2 dimensional saturated porous model according to the number of FEM meshes.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.6
• /
• pp.552-559
• /
• 2009

The effect of porosity on the crack propagation is studied by using the cohesive zone model. Standard mode I fracture test were done by using compact tension specimens with various porosities. Load-load line displacement curves and ${\delta}_5$-crack resistance curves for various porosities were obtained from experiments. The cohesive zone model proposed by Xu and Needleman was employed to describe the crack propagation in porous media, and the Gurson model is used for constitutive relation of porous materials. These models were implemented into user subroutines of a finite element program ABAQUS. The fracture mode changes from ductile fracture to brittle fracture as the porosity increases. Numerical calculations agree well with experimental results.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.20-25
• /
• 2008

The finite element analysis for porous media is severe job because constituents have different physical peoperties, and element's continuity and stability should be considered. Thus, we propose the new mixed finite element method in order to overcome the problems. In this method, multi time step, remeshing step, and sub iteration step are introduced. The multi time step and remeshing step make it possible to satisfy a stability and an accuracy during sub iteration in which global time is determined. Finally, the proposed method is compared with the ABAQUS(2007) software and exact solution(Schiffman 1967) through two dimensional consolidation model.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.129-134
• /
• 2009

In the present study, the characteristic of flow distribution in the channel of a plate heat exchanger is investigated numerically. In order to accomplish the efficient and fast analyses of the flow characteristics in the channel, a semi-microscopic analysis has been performed using a porous media model. For semi-microscopic analysis using porous media, the flow resistance coefficients are obtained through the result of pressure drop in the experimental data. The results showed that the variation of mass flow rate, geometry and chevron angle strongly depend on the flow distribution in the channel. Particularly, the chevron angle is most important factor for uniform flow distribution.