Mathematical simulation of thermomechanical processes in bodies of low electrical conductivity in an external quasistationary electromagnetic field

We give a mathematical formulation and computational scheme for the problem of determining the thermoelastic state of bodies of low electrical conductivity caused by an external quasistationary electromagnetic field. We study the characteristics of these processes in a layer and a hollow cylinder.Translated fromMatematicheskie Metody i Fiziko-Mekhanicheskie Polya, Issue 36, 1992, pp. 34–38.     

Phase field simulation on mechanically induced 180 degree switching in nanomagnets

The mechanically induced 180° switching of the magnetization in cylindrical nanomagnets is demonstrated through the switching dynamics simulation by a constraint-free phase field model which is fully magneto-mechanical coupled. By using the overrun behavior of the magnetization during the mechanically induced precessional switching, a deterministic 180° switching can be achieved. It is found that the residual torque resulting from the azimuthal-angle dependence of the magneto-elastic coupling and the overrun and precessional dynamics are responsible for the switching mechanism. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase field simulation of the intercalation-induced stresses in the hyperelastic solids via isogeometric analysis

In this work the Cahn-Hilliard-type diffusion in a hyperelastic solid will be studied and simulated using isogeometric analysis. The Cahn-Hilliard model is employed to account for the phase segregation phenomenon, in which a fourth-order partial differential operator exists. In order to deal with this high-order operator, isogeometric analysis is applied. Finally, simulation results of a disc is given. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A phase field model for fracture

The variational formulation of brittle fracture as formulated for example by Francfort and Marigo in [1], where the total energy is minimized with respect to any admissible crack set and displacement field, allows the identification of crack paths, branching of preexisting cracks and even crack initiation without additional criteria. For its numerical treatment a continuous approximation of the model in the sense of Γ-convergence has been presented by Bourdin in [2]. In the regularized Francfort–Marigo model cracks are represented by an additional field variable (secondary variable) s∈[0,1] which is 0 if the material is cracked and 1 if it is undamaged. In this work, we reinterpret the crack variable as a phase field order parameter and address cracking as a phase transition problem. The crack growth is governed by the evolution equation of the order parameter which resembles the Ginzburg–Landau equation. The numerical treatment is done by finite elements combined with an implicit Euler scheme for the time integration. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase field simulations in ferroelectric materials

The hindering of domain wall movement by defects in ferroelectric materials is closely connected to electric fatigue. A movable domain wall in a ferroelectric material in most cases is modelled as a singular surface which allows the use of configurational forces. In contrast, the present approach treats the polarization as an order parameter, extending the total energy by a phase separation energy and a domain wall energy. The polarization then no longer has a discontinuity at the domain wall but is a continuous vector field (phase field). As an example, a numerical simulation of domain evolution under stress free boundary conditions is presented. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Interpretation of parameters in phase field models for fracture

Phase field fracture models typically feature a length parameter, which controls the width of the diffuse transition zone between broken and undamaged material. In the limit case of a vanishing length parameter, these models converge to a sharp crack formulation. From this point of view, the length scale parameter is a purely auxiliary numerical quantity. However, the study of the stability of homogeneous solutions in a one dimensional setting permits a different interpretation. Since the length parameter is directly related to the critical stress at which the homogeneous solution becomes unstable and crack nucleation occurs, it can be related to the strength of the material. In this regard, the length parameter itself may be seen as a material parameter. These analytical findings are approved by finite element simulations. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Discharge simulation in a high-frequency quasistatic field

Translated from Metody Matematicheskogo Modelirovaniya i Vychislitel'noi Diagnostika, pp. 275–281, Izd. Moskovskogo Universiteta, Moscow, 1990.     

Phase field modeling of complex crack patterns in multi-physics problems

Recently developed continuum phase field models for brittle fracture show excellent modeling capability in situations with complex crack topologies including branching in the small and large strain applications. This work presents a generalization towards fully coupled multi-physics problems at large strains. A modular concept is outlined for the linking of the diffusive crack modeling with complex multi field material response, where the focus is put on the model problem of finite thermo-elasticity. This concerns a generalization of crack driving forces from the energetic definitions towards stress-based criteria, the constitutive modeling of degradation of non-mechanical fluxes on generated crack faces. Particular assumptions are made on the generation of convective heat exchanges approximating surface load integrals of the sharp crack approach by distinct volume integrals. The coupling effect is also shown in generation of cracks due to thermally induced stress states. We finally demonstrate the performance of the phase field formulation of fracture at large strains by means of representative numerical examples. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On a thermomechanical milling model

This paper deals with a new mathematical model to characterize the interaction between machine and workpiece in a milling process. The model consists of a harmonic oscillator equation for the dynamics of the cutter and a linear thermoelastic workpiece model. The coupling through the cutting force adds delay terms and further nonlinear effects. After a short derivation of the governing equations it is shown that the complete system admits a unique weak solution. A numerical solution strategy is outlined and complemented by numerical simulations of stable and unstable cutting conditions.     

A numerical approach for the simulation of ductile fracture with embedded discontinuities

In the present study, a computational approach for the numerical simulation of ductile fracture within the framework of the finite element method is proposed. In the developed macroscopic formulation, the inelastic behavior in the bulk of the material is described by the finite elasto‐plastic material model proposed in [4]. The failure process is modeled by introducing discontinuities when a special local fracture criterion is satisfied. The discontinuities are incorporated via special triangular finite elements with embedded interfaces following the line of [2]. Finally, the numerical procedure is evaluated for a twodimensional representative test problem. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analysis of a thermomechanical phase transition model

Subject of this work is a macroscopic thermomechanical model of phase transitions in steel. Effects like transformation strain and transformation plasticity induced by the phase transitions are considered and used to formulate a consistent thermomechanical model. The resulting system of state equations consists of a quasistatic momentum balance coupled with a nonlinear stress-strain relation, a nonlinear energy balance equation and a system of ODEs for the phase volume fractions. We prove the existence of a unique weak solution using fixed-point arguments. A key issue is a regularity analysis for the mechanical subsystem to obtain continuity of the stress tensor. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Simulation of Hertzian cone cracks using a phase field description for fracture

The present contribution focuses on fracture caused by indentation loading on the surface of a brittle solid. Its theoretical prediction is a challenging task due to the fact that crack nucleation is not geometrically induced, but is caused by the stress concentration in the contact near-field. The application of the phase field model requires constitutive assumptions to ensure a tension-compression asymmetric material response and prevent damage in compressed regions. This is achieved at the cost of giving up the variational concept of brittle fracture. We simulate the indentation of a cylindrical flat-ended punch on brittle materials like silicate glass. In order to reduce the numerical effort, we exploit axisymmetric conditions for the finite element formulation. After crack initiation stable propagation of a cone crack can be observed in good agreement with experiments. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Non-local modeling of thermomechanical localization in metals

The purpose of this work is to compare and contrast two different modeling approaches for high-speed, dynamic loading of metals. A so-called “local” model including hardening, strain-rate hardening as well as temperature and ductile damage softening behaviour is extended by gradient plasticity. In particular, both models are compared in the context of their application to the modeling and simulation of dynamic shear-banding in the alloy Inconel 718. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical simulation of the flow field of a diffused pneumatic silencer

A diffused pneumatic silencer had been widely used in the pneumatic fields due to its small dimensions and high level of performance in noise reduction. A numerical simulation of its interior and exterior flow field was important for studying the gas flow in the silencer and the flow structure outside the silencer, as well as for understanding the mechanism of the silencer’s noise reduction. A porous media model and the Darcy–Forchheimer principle were used as the basic theoretical models in this paper. The unified governing equations were used here to describe the compressible flow in and out of the silencer. A robust numerical scheme was used to discritize the equations, and the TDBC (Time-dependent boundary conditions) was used to treat the non-reflecting boundaries. The detailed structures of the inner and outer flow fields of the diffused pneumatic silencer were obtained. The simulation results displayed the characteristics of the flow in the silencer. The nature of the flow outside the silencer, comparable with the experimental data, was also obtained.     

A novel treatment of crack boundary conditions in phase field models of fracture

In the phase field approach for fracture an additional scalar field is introduced in order to describe the state of the material between intact and fully broken. So far, for the loading dependent degradation of stiffness (damage) either the volumetric-deviatoric split of strain [1, 2] or the spectral decomposition [3, 4] is used. In contrast to such an isotropic degradation of stiffness, the fully broken state represents a crack with a particular orientation. Both aforementioned approaches do not take the crack orientation into account. This may lead to the violation of the crack boundary conditions. In order to satisfy these conditions the phase field approach is modified here by taking the orientation of the crack into account. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase field simulations on the interaction of ferroelectric domains with point defects

One of the suspected micro-mechanical mechanisms causing electric fatigue in ferroelectric materials is the hindering and blocking of domain wall movement. These blocking or pinning phenomena are thought to be due to point defects which interact with domain walls and applied external loads. A phase field model employing the spontaneous polarization as an order parameter is used to simulate the inhomogeneous material behavior. The coupled field equations are solved using the Finite Element Method. The influence of a stationary point defect on a domain wall is shown in a numerical simulation. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dissipation function in thermomechanical phenomena of viscoelastic solids

Zusammenfassung Eine explizite Form der Dissipationsfunktion für thermomechanische Vorgänge in viskoelastischen Körpern wird abgeleitet. Diese Funktion hängt von beobachtbaren Grössen, wie z. B. Spannungs- und Dehnungsgeschwindigkeiten und ihren zeitlichen Ableitungen ab. Die Formulierung basiert auf der irreversiblen Thermodynamik deformierbarer Körper, die vonBiot entwickelt und vonZiegler undSchapery verallgemeinert wurde. Die vorgeschlagene Dissipationsfunktion ist sowohl zum experimentellen als auch zum theoretischen Gebrauch geeignet.

---

---

The results presented in this paper were obtained in the course of research conducted under Contract N00014-67-A-0109-0003, Task NR 064-496 with the Office of Naval Research, Washington, D.C.     

Unsaturated porous media flow with thermomechanical interaction

We propose a model for unsaturated poro‐plastic flow derived from the thermodynamic principles. For the isothermal case, the problem consists of a degenerate coupled system of two PDEs with two independent hysteresis operators describing hysteresis phenomena in both the solid and the pore fluids. Under natural hypotheses, we prove the existence of a global strong solution for this system. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical simulation of inviscid bubble dynamics in a centrally symmetric gravitational field

The motion of bubbles in a centrally symmetric gravitational field is numerically simulated using two-dimensional conservation laws (Euler equations). The dynamics of bubbles with various numbers of modes in the initial perturbation are studied. The numerical results reveal features that are substantially different from the plane case in a homogeneous gravitational field. Bubble perturbations nearly do not interact at the formation stage. The lowest modes are amplified in the course of the bubble evolution.