Direct finite element method

Elliptic Equations Elhptic equations can be solved with both finite difference and finite element methods. One-dimensional elhptic problems are two-point boundary value problems. Two- and three-dimensional elliptic problems are often solved with iterative methods when the finite difference method is used and direct methods when the finite element method is used. So there are two aspects to consider howthe equations are discretized to form sets of algebraic equations and howthe algebraic equations are then solved. 

These apparent restrictions in size and length of simulation time of the fully quantum-mechanical methods or molecular-dynamics methods with continuous degrees of freedom in real space are the basic reason why the direct simulation of lattice models of the Ising type or of solid-on-solid type is still the most popular technique to simulate crystal growth processes. Consequently, a substantial part of this article will deal with scientific problems on those time and length scales which are simultaneously accessible by the experimental STM methods on one hand and by Monte Carlo lattice simulations on the other hand. Even these methods, however, are too microscopic to incorporate the boundary conditions from the laboratory set-up into the models in a reahstic way. Therefore one uses phenomenological models of the phase-field or sharp-interface type, and finally even finite-element methods, to treat the diffusion transport and hydrodynamic convections which control a reahstic crystal growth process from the melt on an industrial scale. 

The beauty of finite-element modelling is that it is very flexible. The system of interest may be continuous, as in a fluid, or it may comprise separate, discrete components, such as the pieces of metal in this example. The basic principle of finite-element modelling, to simulate the operation of a system by deriving equations only on a local scale, mimics the physical reality by which interactions within most systems are the result of a large number of localised interactions between adjacent elements. These interactions are often bi-directional, in that the behaviour of each element is also affected by the system of which it forms a part. The finite-element method is particularly powerful because with the appropriate choice of elements it is easy to accurately model complex interactions in very large systems because the physical behaviour of each element has a simple mathematical description. 

Direct Energy Minimization Versus the Finite Element Methods. 144... 

An alternative approach to the finite element approach is one, introduced as a concept by Courant as early as 1943 [197], in which the total energy functional, implicit in the finite element method, is directly minimized with respect to all nodal positions. The approach is conjugate to the finite element method and merely differs in its procedural approach. It parallels, however, methods often used in atomistic modeling schemes where the potential energy functional of a system (e. g., given by the force field ) is minimized with respect to the position of all (or at least many) atoms of the system. A simple example of this emerging technique is given below. 

In the finite element method, Petrov-Galerkin methods are used to minimize the unphysical oscillations. The Petrov-Galerkin method essentially adds a small amount of diffusion in the flow direction to smooth the unphysical oscillations. The amount of diffusion is usually proportional to Ax so that it becomes negligible as the mesh size is reduced. The value of the Petrov-Galerkin method lies in being able to obtain a smooth solution when the mesh size is large, so that the computation is feasible. This is not so crucial in one-dimensional problems, but it is essential in two- and three-dimensional problems and purely hyperbolic problems. 

The first group of methods then manipulate a very small subset of vector elements Vi at a time, and a direct method continually updates the affected elements rj. Such methods are collectively known as relaxation methods, and they are primarily used in situations where, for each elements >< to be changed, the set of affected rjt and the matrix elements Ay, are immediately known. This applies in particular to difference approximations and also to the so-called Finite Element Method for obtaining tabular descriptions of the wave function, i.e. a list of values of the wave function at a set of electron positions. (Far some reason, such a description is commonly referred to as a numerical solution to the Schrodinger equation). Relaxation methods have also been applied to the Cl problem in the past, but due to their slow convergence they have been replaced by analytical methods. Even for numerical problems, the relaxation methods are slowly yielding to analytical methods. 

Continuum wavefunctions can also be generated by solving the partial differential equation (3.3) directly without first transforming it into a set of ordinary differential equations. One possible scheme is the finite elements method (Askar and Rabitz 1984 Jaquet 1987). Another method, which has been applied for the calculation of multi-dimensional scattering wavefunctions, is the 5-matrix version of the Kohn variational principle (Zhang and Miller 1990). 

With the development of high-speed personal computers, it is very convenient to use numerical techniques to solve heat transfer problems. The finite-difference method and the finite-element method are two popular and useful methods. The finite-element method is not as direct, conceptually, as the finite-difference method. It has some advantages over the finite-difference method in solving heat transfer problems, especially for problems with complex geometries. 

Furthermore, the stress and local preferred direction may vary from point to point in the continuum. Thus Eqn. (7) implies that the stress field and unit vector field, i.e. (Ti, xk) and d xk), must be specified to define . Current state of the art in designing structural components uses finite element methods to characterize the stress field. The reader will see shortly that this is conveniently utilized in analyzing component reliability. 

In order to solve the first principles model, finite difference method or finite element method can be used but the number of states increases exponentially when these methods are used to solve the problem. Lee et u/.[8] used the model reduction technique to reslove the size problem. However, the information on the concentration distribution is scarce and the physical meaning of the reduced state is hard to be interpreted. Therefore, we intend to construct the input/output data mapping. Because the conventional linear identification method cannot be applied to a hybrid SMB process, we construct the artificial continuous input/output mapping by keeping the discrete inputs such as the switching time constant. The averaged concentrations of rich component in raffinate and extract are selected as the output variables while the flow rate ratios in sections 2 and 3 are selected as the input variables. Since these output variables are directly correlated with the product purities, the control of product purities is also accomplished. 

Neither the finite element method nor a discretised integration will "catch" eddy or other motion below a certain scale determined by the choice of mesh. Small-scale motion in many cases may be better described as random, in which case the transport of the quantity A is called diffusion. Diffusion can be described by Pick s law, assuming that the flux density /, i.e., the number of "particles" (here a small parcel of a gas), passing a unit square in a given direction, is proportional to the negative gradient of particle concentration n (Bockris and Despic, 2004),... 

The pairwise Brownian dynamics method has several advantages over numerical methods based on Smoluchowski s [9] approach (e.g., finite element method), and we discuss these here. The primary advantage of the method is the ease of mathematical formulation even for cases involving complex reaction site geometries, hydrodynamic interactions, charge effects, anisotropic diffusion and flow fields. Furthermore the method obviates the need to solve complex diffusion equations to obtain the concentration field from which the rate constant is calculated in the Smoluchowski method. In contrast, the rate constant is obtained directly in the pairwise Brownian dynamics method. The effective rate constants for different reaction conditions may be obtained from a single simulation this is not possible using the finite element method. 

The mathematical model of a MAT reactor, considering a 12 lump model, has been discretized using a finite element method in the direction of gas flow. The resulting system of differential-algebraic equations (DAEs) has been solved by an appropriate computer code (DASSL). 

One method to solve partial differential equations using the numerical schemes developed for solving time dependent ordinary differential methods is the method of lines. In this method, the spatial derivatives at time t are replaced by discrete approximations such as finite differences or finite element methods such as collocation or Galerkin. The reason for this approach is the advanced stage of development of schemes to solve ordinary differential equations. The resulting numerical schemes are frequently similar to those developed directly for partial differential equations. 

As a complementary approach a mean-field method is used in combination with the finite element method to investigate the FGM. To compare the predictions with the periodic unit cell simulations, the FGM part is divided into nine sublayers. Each of them consists of two bi-quadratic 8-node plane elements over the thickness, and only one element is used in the horizontal direction. The parts of pure alumina and pure nickel are modeled by three and 12 elements, respectively. In each sublayer the volume fractions of the phases are constant. In addition, the center sublayer can be split into a metal-matrix and a ceramic-matrix half sublayer. The properties of the particular material on the meso-structural level within the finite element calculation are described via a con-... 

---