A.B. Tatomir^1,2, A.Szymkiewicz³, H. Class¹, R. Helmig¹

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.2, pp. 81-112, 2011, DOI:10.3970/cmes.2011.077.081

Abstract We present a two phase flow conceptual model, the corresponding simulator (2pMINC) and a workflow for large-scale fractured reservoirs, based on a continuum fracture approach which uses the multiple interacting continua (MINC) method complemented with an improved upscaling technique. The complex transient behavior of the flow processes in fractured porous media is captured by subgridding the coarse blocks in nested volume elements which have effective properties calculated from the detailed representation of the fracture system. In this way, we keep a physically based approach, preserve the accuracy of the model, avoid the common use of empirically derived transfer functions and… More >

W. A. El-Askary¹

CMES-Computer Modeling in Engineering & Sciences, Vol.74, No.3&4, pp. 203-232, 2011, DOI:10.3970/cmes.2011.074.203

Abstract Large eddy simulation (LES) is a viable and powerful tool to analyze unsteady three- dimensional turbulent flows. In this paper, the method of LES is used to compute a plane turbulent supersonic boundary layer subjected to different pressure gradients. The pressure gradients are generated by allowing the flow to pass in the vicinity of an expansion-compression ramp (inclined backward-facing step with leeward-face angle of 25 degrees) for an upstream Mach number of 2.9. The inflow boundary condition is the main problem for all turbulent wall-bounded flows. An approach to solve this problem is to extract instantaneous velocities, temperature and density… More >

Xiaowei Zeng¹

CMES-Computer Modeling in Engineering & Sciences, Vol.74, No.3&4, pp. 183-202, 2011, DOI:10.3970/cmes.2011.074.183

Abstract This paper presents an atomistic field theory and its application in modeling and simulation of nano/micro materials. Atomistic formulation and finite element implementation of the atomistic field theory is briefly introduced. Numerical simulations based on the field theory are performed to investigate the material behaviors of bcc iron at coarse-grained scale and we have obtained the mechanical strength and elastic modulus, which are in good agreement with results by first principles calculations. Also the nanoscale deformation and failure mechanism are revealed in bcc iron nanorods under simple tension. It is interesting to observe that under tensile loading, iron has gone… More >

Shiqin He¹, Pengfei Li¹, Feng Shang²

CMES-Computer Modeling in Engineering & Sciences, Vol.73, No.4, pp. 387-394, 2011, DOI:10.3970/cmes.2011.073.387

Abstract A three-dimensional finite element (FE) and analytical approach for the simulation of the shear properties of steel-concrete composite beams are presented in this paper. To simulate the interfacial behavior between steel girders and concrete slabs, we apply an interface slip model in the simulation. This model has been used in analyzing the flexural properties of composite beams. Both simply supported beam and continuous composite beam experiments reported in literature are simulated. The load deflection and slip rule between steel girders and concrete slabs, as well as the crack pattern and contour at the ultimate load, are analyzed. The results obtained… More >

T. Nakamura¹, R. Kobyashi¹, S. Ogata¹

CMES-Computer Modeling in Engineering & Sciences, Vol.73, No.4, pp. 357-386, 2011, DOI:10.3970/cmes.2011.073.357

Abstract The coarse-grained particle (CGP) method has been proposed to coarse-grain a crystalline system of atoms to meso-scale. In the method, virtual particles are distributed in the system, and the inter-particle interaction is calculated through the constrained statistical ensemble average of the atomic Hamiltonian at a given temperature. For simplicity, however, the harmonic approximation has been used for the inter-atomic interaction and hence anharmonicity at finite temperatures has been ignored. We improve the former CGP method to incorporate the anharmonicity of atomic system at finite temperatures into the inter-particle interaction. Also the divide-and-conquer strategy is applied to calculate the inter-particle interaction… More >

A. Gakwaya¹, H. Sharifi², M. Guillot¹, M. Souli³, F. Erchiqui⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.73, No.3, pp. 209-266, 2011, DOI:10.3970/cmes.2011.073.209

Abstract A mechanical computer aided design and engineering system can be used to reduce the design-to-manufacture cycle time in metal forming process. Such a system could be built upon a solid modeling geometry engine and an efficient finite element (FE) solver. The maintenance of a high-quality mesh throughout the analysis is an essential feature of an efficient finite element simulation of large strain metal forming problems. In this paper, a mesh adaptation technique employing the Arbitrary Lagrangian-Eulerian formulation (ALE) is applied to some industrial metal forming problems. An ACIS boundary representation of the solid model is employed. This type of representation… More >

Chun Yang^1,2, Dalin Tang², Satya Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.72, No.1, pp. 53-78, 2011, DOI:10.3970/cmes.2011.072.053

Abstract Previously, we introduced a computational procedure based on three-dimensional meshless generalized finite difference (MGFD) method and serial magnetic resonance imaging (MRI) data to quantify patient-specific carotid atherosclerotic plaque growth functions and simulate plaque progression. Structure-only models were used in our previous report. In this paper, fluid-stricture interaction (FSI) was added to improve on prediction accuracy. One participating patient was scanned three times (T1, T2, and T3, at intervals of about 18 months) to obtain plaque progression data. Blood flow was assumed to laminar, Newtonian, viscous and incompressible. The Navier-Stokes equations with arbitrary Lagrangian-Eulerian (ALE) formulation were used as the governing… More >

B. Ren^1,2, J. Qian¹, X. Zeng¹, A. K. Jha³, S. Xiao⁴, S. Li^1,5

CMES-Computer Modeling in Engineering & Sciences, Vol.71, No.3, pp. 253-278, 2011, DOI:10.3970/cmes.2011.071.253

Abstract Ductile failure is a complex multi-scale phenomenon evolved from the micro-voids to macro-crack. There are three main failure mechanisms behinds a ductile failure: adiabatic shear band (ASB), spall fracture, and crack. Since this type of thermo-mechanical phenomena involves large deformation and large scale plastic yielding, a meshfree method has intrinsic advantages in solving this kind of problems over the conventional finite element method. In this paper, the numerical methodologies including multi-physics approach for ASB, parametric visibility condition for crack propagation, and multi-scale approach to determine spall strength in simulating ductile failure have been reviewed. A thermo-mechanical coupling algorithm is proposed… More >

S. Abbasbandy^1,2, V. Sladek³, A. Shirzadi¹, J. Sladek³

CMES-Computer Modeling in Engineering & Sciences, Vol.71, No.1, pp. 15-38, 2011, DOI:10.3970/cmes.2011.071.015

Abstract This paper deals with the development of a new method for solution of initial-boundary value problems governed by a couple of nonlinear diffusion equations occurring in the theory of self-organization in non-equilibrium systems. The time dependence is treated by linear interpolation using the finite difference method and the semi-discrete partial differential equations are considered in a weak sense by using the local integral equation method with approximating 2-d spatial variations of the field variables by the Moving Least Squares. The evaluation techniques are discussed and the applicability of the presented method is demonstrated on two illustrative examples with exact solutions… More >

Yuta Kitamura¹, Noriyuki Miyazaki¹, Takahito Kumazaki², Naoto Nagakura³, Yasuhiro Hashimoto³, Isao Masada³

CMES-Computer Modeling in Engineering & Sciences, Vol.68, No.2, pp. 151-166, 2010, DOI:10.3970/cmes.2010.068.151

Abstract CaF₂ single crystal is used as high performance optical elements. We developed an analysis system for simulating birefringence of CaF₂ single crystal used as a chamber window of a gas laser light source. The analysis system consists of a stress analysis and a birefringence analysis. In the stress analysis, the finite element method was applied to obtain the mechanical stress caused by a window holder and gas pressure. In the birefringence analysis, the photo-elastic effect gives the change of refractive indices, from which the optical path difference and the fast axis are calculated by using the average stress method. The… More >