• The Journal of Engineering Geology
• /
• v.19 no.2
• /
• pp.205-215
• /
• 2009

The conventional numerical models to analyze flow in subsurface porous media under the transient state usually generate numerical oscillation and unstability due to local flux domain for critical cases such as infiltration into initially dry soil during rainfall period. In this case, it is required refined mesh and small time step, but it decrease efficiency of computation. In this study, numerical unstability in discontinuity domain is removed by applying particle tracking algorithm to simulate unsteady subsurface flow with inflow boundary condition. Finally the hybrid LE FEM improving numerical stability is proposed. The hypothetical domains with unsteady uniform and nonuniform flow field were used to demonstrated algorithm verification. In comparison with analytic solution, we obtained reasonable results and conducted simulation of hypothetical 3-D recharge/pumping area. The proposed algorithm can simulate saturated/unsaturated porous media with more practical problems and will greatly contribute to accuracy and stability of numerical computation.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2001.05a
• /
• pp.201-204
• /
• 2001

A Confocal Laser Scanning Microscope (CLSM) is applied to determine three-dimensional fiber orientation states in injection-molded short fiber composites. Since the CLSM optically sections the composites, more than two planes either on or below the surface of composites can be obtained. Therefore, three dimensional fiber orientation states are determined without destruction. To predict the orientation states, velocity and temperature fields are calculated by using a hybrid FEM/FDM method. The change of orientation state during packing stage is also considered by employing a compressible Hele-Shaw model. The predicted orientation states show good agreement with measured ones. However, some differences are found at the end of cavity. They may result from other effects, which are not considered in the numerical analysis.