A surface oriented solid formulation based on a hybrid Galerkin-collocation method

The contribution is concerned with a numerical method to analyze the mechanical behavior of 3D solids. The method employs directly the geometry defined by the boundary representation modeling technique, which is frequently used in CAD to define solids. It combines the benefits of the isogeometric analysis methodology with the scaled boundary finite element method. In the present approach, only the boundary surfaces of the solid are discretized. No tensor-product structure of three-dimensional objects is exploited to parametrize the physical domain. The weak form is applied only on the boundary surfaces. The governing partial differential equations of elasticity are transformed to an ordinary differential equation (ODE) of Euler type. The isogeometric Galerkin approach is employed to approximate the displacement response at the boundary surfaces. It exploits the two-dimensional NURBS objects to parametrize the boundary surfaces. To solve the Euler type ODE, the NURBS based collocation approach is applied. The accuracy of the method is validated against the analytical solutions. The presented method is able to analyze solids, which are bounded by an arbitrary number of surfaces. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Application and modification of the POD method and the POD-DEIM for model reduction in porous-media simulations

Researchers with a continuum-mechanical background typically use a multi-phasic and multi-component modelling approach for materials with a saturated porous microstructure. Therefore, the mechanical behaviour is considered in a continuum-mechanical manner and solved using the finite-element method (FEM). The developed models need to be complex enough to capture the relevant properties of the considered materials, what often results in expensive simulations with a very large number of degrees of freedom (DOF). The aim of the present contribution is to reduce the computing time of these simulations through model-reduction methods, while the accuracy of the solution needs to be maintained. Therefore, the method of proper-orthogonal decomposition (POD) for linear problems and the discrete-empirical-interpolation method (DEIM) in combination with the POD method (POD-DEIM) for nonlinear problems are investigated. Using the POD method, a given data set is approximated with a low-dimensional subspace. To generate this data set, the vector of unknowns of the FE simulation is stored in a pre-computation in the full (unreduced) system in each time-step (so-called “snapshots” of the system). Dealing with porous-media problems, the primary variables are the solid displacement, the pore pressure and, depending on the particular problem, other primary variables. Following this, the primary variables have entries with very huge differences in their absolute values. As a result, non-negligible rounding errors may occur when applying the POD method. To overcome this problems, modifications of the classical POD method need to be performed for such problems. The present contribution discusses this issue and presents results for the reduced simulations of porous media. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Parallel solution of volume-coupled multi-field problems using an Abaqus-PANDAS software interface

Numerical simulations have proven to be a powerful tool in several engineering disciplines, such as mechanical, civil and biomechanical engineering, and are thus widely used. However, the reliability of the simulations strongly relies on the governing material model. These models are usually developed in academic or industrial research projects and are implemented into dedicated software packages to proof their concepts. A transfer of these models from the research into a production-related environment is often time consuming and prone to failures, and therefore a costly task. The present work introduces a general interface between the research code PANDAS, which is a dedicated multi-field finite-element solver based on a monolithic solution strategy, and the commercial finite-element package Abaqus. The coupling is based on the user-defined element subroutine (UEL) of Abaqus. This procedure, on the one hand, allows for a straight-forward embedding of the PANDAS material models into Abaqus. On the other hand, it provides, in comparison to the native UEL subroutine of Abaqus, a user-friendly programming environment for user-defined material models with an extended number of degrees of freedom. Furthermore, the coupling also supports the parallel-analysis capabilities for large-scale problems on high-performance computing clusters. The Abaqus-PANDAS linkage can be applied to various coupled multi-field problems. However, the present contribution addresses, in particular, volume-coupled multi-field problems as they arise when proceeding from the Theory of Porous Media (TPM) as a modelling framework. For instance, it can be used to model partially or fully saturated soils, or chemically or electro-chemically driven swelling phenomena as they appear, for example, within hydrogels. Additionally, discontinuities, such as cracks, can be described for instance via phase-field models or by the extended finite-element method (XFEM). (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A method for C(alpha) direct-detection in protonless NMR

Attempts are made to efficiently decouple (13)C nuclei without significant loss of coherence during the application of the decoupling package. Such attempts are based on the S(3)E spin-state selection method. A newly developed double S(3)E (DS(3)E) is particularly efficient for C(alpha) detection for proteins as large as 480 kDa.     

Locating saddle points of any index on potential energy surfaces by the generalized gentlest ascent dynamics

The system of ordinary differential equations for the method of the gentlest ascent dynamics (GAD) is tested to determine the saddle points of the potential energy surface of some molecules. The method has been proposed earlier [E and Zhou in Nonlinearity 24:1831 (2011)]. We additionally use the metric of curvilinear internal coordinates. By a number of examples, we explain the possibilities of a GAD curve; it can find the transition state of interest by a gentlest ascent, directly or indirectly, or not. A GAD curve can be a model of a reaction path, if it does not contain a turning point for the energy. We further discuss generalized GAD formulas for the search of saddle points of a higher index. We calculate diverse examples.     

A Strategy for the Preparation of Thioantimonates Based on the Concept of Weak Acids and Corresponding Strong Bases

By following a new synthetic approach, which is based on the in situ formation of a basic medium by the reaction between the strong base Sb(V)S₄^3? and the weak acid H₂O, it was possible to prepare three layered thioantimonate(III) compounds of composition [TM(2,2′‐bipyridine)₃][Sb₆S₁₀] (TM=Ni, Fe) and [Ni(4,4′‐dimethyl‐2,2′‐bipyridine)₃][Sb₆S₁₀] under hydrothermal conditions featuring two different thioantimonate(III) network topologies. The antimony source, Na₃SbS₄ ? 9 H₂O, undergoes several decomposition reactions and produces the Sb^IIIS₃ species, which condenses to generate the layered anion. The application of transition‐metal complexes avoids crystallization of dense phases. The reactions are very fast compared to conventional hydrothermal/solvothermal syntheses and are much less sensitive to changes of the reaction parameters.