Continuity of Lyapunov functions and of energy level for generalized gradient semigroup

The global attractor of a gradient-like semigroup has a Morse decomposition. Associated to this Morse decomposition there is a Lyapunov function (differentiable along solutions)-defined on the whole phase space- which proves relevant information on the structure of the attractor. In this paper we prove the continuity of these Lyapunov functions under perturbation. On the other hand, the attractor of a gradient-like semigroup also has an energy level decomposition which is again a Morse decomposition but with a total order between any two components. We claim that, from a dynamical point of view, this is the optimal decomposition of a global attractor; that is, if we start from the finest Morse decomposition, the energy level decomposition is the coarsest Morse decomposition that still produces a Lyapunov function which gives the same information about the structure of the attractor. We also establish sufficient conditions which ensure the stability of this kind of decomposition under perturbation. In particular, if connections between different isolated invariant sets inside the attractor remain under perturbation, we show the continuity of the energy level Morse decomposition. The class of Morse-Smale systems illustrates our results.

Invariant Solutions of Richards' Equation for Water Movement in Dissimilar Soils

Scaling methods allow a single solution to Richards' equation (RE) to suffice for numerous specific cases of water flow in unsaturated soils. During the past half-century, many such methods were developed for similar soils. In this paper, a new method is proposed for scaling RE for a wide range of dissimilar soils. Exponential-power (EP) functions are used to reduce the dependence of the scaled RE on the soil hydraulic properties. To evaluate the proposed method, the scaled RE was solved numerically considering two test cases: infiltration into relatively dry soils having initially uniform water content distributions, and gravity-dominant drainage occurring from initially wet soil profiles. Although the results for four texturally different soils ranging from sand to heavy clay (adopted from the UNSODA database) showed that the scaled solution were invariant for a wide range of flow conditions, slight deviations were observed when the soil profile was initially wet in the infiltration case or deeply wet in the drainage case. The invariance of the scaled RE makes it possible to generalize a single solution of RE to many dissimilar soils and conditions. Such a procedure reduces the numerical calculations and provides additional opportunities for solving the highly nonlinear RE for unsaturated water flow in soils.

A generalized finite element method with global-local enrichment functions for confined plasticity problems

The main feature of partition of unity methods such as the generalized or extended finite element method is their ability of utilizing a priori knowledge about the solution of a problem in the form of enrichment functions. However, analytical derivation of enrichment functions with good approximation properties is mostly limited to two-dimensional linear problems. This paper presents a procedure to numerically generate proper enrichment functions for three-dimensional problems with confined plasticity where plastic evolution is gradual. This procedure involves the solution of boundary value problems around local regions exhibiting nonlinear behavior and the enrichment of the global solution space with the local solutions through the partition of unity method framework. This approach can produce accurate nonlinear solutions with a reduced computational cost compared to standard finite element methods since computationally intensive nonlinear iterations can be performed on coarse global meshes after the creation of enrichment functions properly describing localized nonlinear behavior. Several three-dimensional nonlinear problems based on the rate-independent J (2) plasticity theory with isotropic hardening are solved using the proposed procedure to demonstrate its robustness, accuracy and computational efficiency.