Binxin Yang¹, Jie Ouyang¹, Tao Jiang¹, Chuntai Liu²

CMES-Computer Modeling in Engineering & Sciences, Vol.63, No.3, pp. 191-222, 2010, DOI:10.3970/cmes.2010.063.191

Abstract A gas-solid-liquid three-phase model is proposed for fiber reinforced composites mold filling process. The fluid flow is described in Eulerian coordinate while the dynamics of fibers is described in Langrangian coordinate. The interaction of fluid flow and fibers are enclosed in the model. The influence of fluid flow on fibers is described by the resultant forces imposed on fibers and the influence of fibers on fluid flow is described by the momentum exchange source term in the model. A finite volume method coupled with a level set method for viscoelastic-Newtonian fluid flow is used to More >

Yue-Ting Zhou¹, Xing Li², De-Hao Yu³, Kang Yong Lee^1,4

CMES-Computer Modeling in Engineering & Sciences, Vol.63, No.2, pp. 163-190, 2010, DOI:10.3970/cmes.2010.063.163

Abstract In this paper, a coupled crack/contact model is established for the composite material with arbitrary periodic cracks indented by periodic punches. The contact of crack faces is considered. Frictional forces are modeled to arise between the punch foundation and the composite material boundary. Kolosov-Muskhelisvili complex potentials with Hilbert kernels are constructed, which satisfy the continuity conditions of stress and displacement along the interface identically. The considered problem is reduced to a system of singular integral equations of first and second kind with Hilbert kernels. Bounded functions are defined so that singular integral equations of Hilbert More >

Y. Inoue¹, R. Kobayashi¹, S. Ogata¹, T. Gotoh¹

CMES-Computer Modeling in Engineering & Sciences, Vol.63, No.2, pp. 137-162, 2010, DOI:10.3970/cmes.2010.063.137

Abstract Vibration of a single graphene and a pair of graphenes at micro meter scale induced by air flow is numerically simulated and examined by using a hybrid computational method starting from a microscopic level of description for the graphene. In order to bridge a huge gap in spatial and time scales in their motions, the carbon atoms of the graphene are represented by a small number of coarse grained particles, the fluid motion is described by the lattice Boltzmann equation and the momentum exchange at the boundary is treated by the time averaged immersed boundary… More >

A.N.M. Imroz Choudhury¹, Michael D. Steffen¹, James E. Guilkey², Steven G.Parker³

CMES-Computer Modeling in Engineering & Sciences, Vol.63, No.2, pp. 117-136, 2010, DOI:10.3970/cmes.2010.063.117

Abstract We present a physically based method for visualizing deformation in particle simulations, such as those describing structural mechanics simulations. The method uses the deformation gradient tensor to transform carefully chosen glyphs representing each particle. The visualization approximates how simulated objects responding to applied forces might look in reality, allowing for a better understanding of material deformation, an important indicator of, for example, material failure. It can also help highlight possible errors and numerical deficiencies in the simulation itself, suggesting how simulations might be changed to yield more accurate results. More >

J.H. Tian^1,2, X. Han², S.Y. Long², G.Y. Sun², Y. Cao¹, G.Q. Xie³

CMES-Computer Modeling in Engineering & Sciences, Vol.63, No.2, pp. 101-116, 2010, DOI:10.3970/cmes.2010.063.101

Abstract A transient heat conduction analysis of the functionally graded material (FGM) plates has been investigated based on the hybrid numerical method (HNM). HNM combines the layer element method with the method of Fourier transforms and proves to be efficient and reliable. The FGM plates are infinite large and the material properties vary continuously through thickness. The transient heat source acted on the FGM plates. The temperature distribution of the FGM plates is obtained in different time and different position. Some useful results for transient heat conduction are shown in figures. Applications of HNM to transient More >

Alex Ferguson¹, Matthew R. Smith², J.-S. Wu³

CMES-Computer Modeling in Engineering & Sciences, Vol.63, No.1, pp. 79-100, 2010, DOI:10.3970/cmes.2010.063.079

Abstract A novel approach for the use of multiple continuous uniform distributions for reconstruction of the Maxwell-Boltzmann equilibrium probability distribution function is used for the solution of one and two dimensional Euler equations. The Uniform distribution Equilibrium Flux Method (UEFM) is a kinetic-theory based flux solver which calculates true directional, volume to volume fluxes based on integration (over velocity space and physical space) of a sum of uniform probability distribution functions working to approximate the equilibrium distribution function. The resulting flux expressions contain only the Heaviside unit step function and do not require the evaluation of More >

Roland Martin¹, Carlos Couder-Castaneda¹

CMES-Computer Modeling in Engineering & Sciences, Vol.63, No.1, pp. 47-78, 2010, DOI:10.3970/cmes.2010.063.047

Abstract We develop an unsplit convolutional perfectly matched layer (CPML) technique to absorb efficiently compressible viscous flows and their related supersonic or subsonic regimes at the outer boundary of a distorted computational domain. More particularly subsonic outgoing flows or subsonic wall-boundary layers close to the PML are well absorbed, which is difficult to obtain without creating numerical instabilities over long time periods. This new PML (CPML) introduces the calculation of auxiliary memory variables at each time step and allows an unsplit formulation of the PML. Damping functions involving a high shift in the frequency domain allow… More >

B.W. Kim¹, H.G. Sung¹, S.Y. Hong¹, H.J. Jung²

CMES-Computer Modeling in Engineering & Sciences, Vol.63, No.1, pp. 29-46, 2010, DOI:10.3970/cmes.2010.063.029

Abstract This paper numerically investigates the behavior of sag and tension of inclined catenary structure considering elastic deformation. Equilibrium equation for computing elastic catenary is formulated by employing finite element method (FEM). Minimum potential energy principle and the Lagrange multiplier method are used in the formulation to derive equilibrium equation with constraint condition for catenary length. Since stiffness and loading forces of catenary are dependent on its own geometry, the equilibrium equation is nonlinear. Using the iterative scheme such as fixed point iteration or bisection, equilibrium position and tension are found. Based on the formulation, a More >

Ching Kong Chao^1,2, Chin Kun Chen¹, Fu Mo Chen³

CMES-Computer Modeling in Engineering & Sciences, Vol.63, No.1, pp. 1-28, 2010, DOI:10.3970/cmes.2010.063.001

Abstract A general analytical solution for a reinforced elliptical hole embedded in an infinite matrix subjected to a point heat source is provided in this paper. Based on the technique of conformal mapping and the method of analytical continuation in conjunction with the alternating technique, the general expressions of the temperature and stresses in the reinforcement layer and the matrix are derived explicitly in a series form. Some numerical results are provided to investigate the effects of the material combinations and geometric configurations on the interfacial stresses. The solution obtained can be treated as Green's functions More >

Sung-Eun Kim¹, Shin Hyung Rhee²

CMES-Computer Modeling in Engineering & Sciences, Vol.62, No.3, pp. 291-310, 2010, DOI:10.3970/cmes.2010.062.291

Abstract Turbulent flow past a finite wing has been computed to assess the fidelity of modern computational fluid dynamics in predicting tip vortex flows. The efficacy of a feature-adaptive local mesh refinement to resolve the steep gradients in the flow field near the tip vortex is demonstrated. The impact of turbulence modeling is evaluated using several popular eddy viscosity models and a Reynolds stress transport model. The results indicate that the combination of a computational mesh with an adequate resolution, high-order spatial discretization scheme along with the use of advanced turbulence models can predict tip vortex More >