1 The performance of the developed Java finite element code is compared to that of the C finite element code on the solution of three-dimensional elasticity problems using Intel Pentium 4 computer. Untuned Java code is approximately two times slower then analogous C code. It is shown that code tuning with the use of blocking technique can provide Java/C performance ratio 90% for the LDU solutio...

In the present work, we study various numerical aspects of higher-order finite-element discretizations of the non-linear saddle-point formulation of orbital-free density-functional theory. We first investigate the robustness of viable solution schemes by analyzing the solvability conditions of the discrete problem. We find that a staggered solution procedure where the potential fields are compu...

We consider a space-time finite element discretization of a time-dependent nonlinear hyperbolic conservation law in one space dimension (Burgers' equation). The finite element method is higher-order accurate and is a Petrov-Galerkin method based on the so-called streamline diffusion modification of the test functions giving added stability. We first prove that if a sequence of finite element so...

in this paper, a finite element model is developed for the fully hydro-mechanical analysis of hydraulic fracturing in partially saturated porous media. the model is derived from the framework of generalized biot theory. the fracture propagation is governed by a cohesive fracture model. the flow within the fracture zone is modeled by the lubrication equation. the displacement of solid phase, and...

We consider the finite element solution of the vector Laplace equation on a domain in two dimensions. For various choices of boundary conditions, it is known that a mixed finite element method, in which the rotation of the solution is introduced as a second unknown, is advantageous, and appropriate choices of mixed finite element spaces lead to a stable, optimally convergent discretization. How...

where ∇ and ∇· are the gradient and divergence operators, Ω ⊂ R is an open bounded cubic domain with boundary Γ, n indicates the outward unit normal vector along Γ, A−1 = (αij)3×3 is a full positive definite matrix uniformly in Ω. Mixed finite element methods [1] should be employed to discretize the system (1.1). The main content of this paper is to present an analysis for the extrapolation of ...

Displacement finite element models of various beam theories have been developed using traditional finite element interpolations (i.e., Hermite cubic or equi-spaced Lagrange functions). Various finite element models of beams differ from each other in the choice of the interpolation functions used for the transverse deflection w, total rotation φ and/or shear strain γxz, or in the integral form u...

The paper is concerned with the discontinuous Galerkin finite element method for the numerical solution of nonlinear conservation laws and nonlinear convection-diffusion problems with emphasis on applications to the simulation of compressible flows. We discuss two versions of this method: (a) Finite volume discontinuous Galerkin method, which is a generalization of the combined finite volume—fi...

One of the important objectives in teaching finite element method at introductory level is to bring students into the comprehension of finite element procedures. This study presents a strategy of teaching structural line elements involving an open source computer-aided learning tool FreeMat integrated with another open source CALFEM finite element toolbox. FreeMat, which is a programming based ...