A linearised hp–finite element framework for acousto‐magneto‐mechanical coupling in axisymmetric MRI scanners

We propose a new computational framework for the treatment of acousto‐magneto‐mechanical coupling that arises in low‐frequency electro‐magneto‐mechanical systems such as magnetic resonance imaging scanners. Our transient Newton–Raphson strategy involves the solution of a monolithic system obtained from the linearisation of the coupled system of equations. Moreover, this framework, in the case of excitation from static and harmonic current sources, allows us to propose a simple linearised system and rigorously motivate a single‐step strategy for understanding the response of systems under different frequencies of excitation. Motivated by the need to solve industrial problems rapidly, we restrict ourselves to solving problems consisting of axisymmetric geometries and current sources. Our treatment also discusses in detail the computational requirements for the solution of these coupled problems on unbounded domains and the accurate discretisation of the fields using hp–finite elements. We include a set of academic and industrially relevant examples to benchmark and illustrate our approach. Copyright © 2017 The Authors. International Journal for Numerical Methods in Engineering Published by John Wiley & Sons, Ltd.     

Cyclic steady states of nonlinear electro‐mechanical devices excited at resonance

We present an efficient numerical method to solve for cyclic steady states of nonlinear electro‐mechanical devices excited at resonance. Many electro‐mechanical systems are designed to operate at resonance, where the ramp‐up simulation to steady state is computationally very expensive – especially when low damping is present. The proposed method relies on a Newton–Krylov shooting scheme for the direct calculation of the cyclic steady state, as opposed to a naïve transient time‐stepping from zero initial conditions. We use a recently developed high‐order Eulerian–Lagrangian finite element method in combination with an energy‐preserving dynamic contact algorithm in order to solve the coupled electro‐mechanical boundary value problem. The nonlinear coupled equations are evolved by means of an operator split of the mechanical and electrical problem with an explicit as well as implicit approach. The presented benchmark examples include the first three fundamental modes of a vibrating nanotube, as well as a micro‐electro‐mechanical disk resonator in dynamic steady contact. For the examples discussed, we observe power law computational speed‐ups of the form S  = 0.6·ξ  ^? 0.8, where ξ is the linear damping ratio of the corresponding resonance frequency. Copyright © 2016 John Wiley & Sons, Ltd.     

Non‐weak/strong solutions in gas dynamics: a C11 p‐version STLSFEF in Eulerian frame of reference using ρ, u,p primitive variables

This paper presents a space–time least squares finite element formulation of one‐dimensional transient Navier–Stokes equations (governing differential equations: GDE) for compressible flow in Eulerian frame of reference using ρ, u, p as primitive variables with C¹¹ type p‐version hierarchical interpolations in space and time. Time marching procedure is utilized to compute time evolutions for all values of time. For high speed gas dynamics the C¹¹ type interpolations in space and time possess the same orders of continuity in space and time as the GDE. It is demonstrated that with this approach accurate numerical solutions of Navier–Stokes equations are possible without any assumptions or approximations. In the approach presented here SUPG, SUPG/DC, SUPG/DC/LS operators are neither used nor needed. Time accurate numerical simulations show resolution of shock structure (i.e. shock speed, shock relations and shock width) to be in excellent agreement with the analytical solutions. The role of diffusion i.e. viscosity (physical or artificial) and thermal conductivity on shock structure is demonstrated. Riemann shock tube is used as a model problem. True time evolutions are reported beginning with the first time step until steady shock conditions are achieved. In this approach, when the computed error functionals become zero (computationally), the computed non‐weak solutions have characteristics as those of the strong solutions of the gas dynamics equations. Copyright © 2001 John Wiley & Sons, Ltd.     

A combined FIC‐TDG finite element approach for the numerical solution of coupled advection–diffusion–reaction equations with application to a bioregulatory model for bone fracture healing

Numerical schemes for the approximative solution of advection–diffusion–reaction equations are often flawed because of spurious oscillations, caused by steep gradients or dominant advection or reaction. In addition, for strong coupled nonlinear processes, which may be described by a set of hyperbolic PDEs, established time stepping schemes lack either accuracy or stability to provide a reliable solution. In this contribution, an advanced numerical scheme for this class of problems is suggested by combining sophisticated stabilization techniques, namely the finite calculus (FIC‐FEM) scheme introduced by Oñate et al. with time‐discontinuous Galerkin (TDG) methods. Whereas the former one provides a stabilization technique for the numerical treatment of steep gradients for advection‐dominated problems, the latter ensures reliable solutions with regard to the temporal evolution. A brief theoretical outline on the superior behavior of both approaches will be presented and underlined with related computational tests. The performance of the suggested FIC‐TDG finite element approach will be discussed exemplarily on a bioregulatory model for bone fracture healing proposed by Geris et al., which consists of at least 12 coupled hyperbolic evolution equations. Copyright © 2012 John Wiley & Sons, Ltd.     

A rapid simulation of nano‐particle transport in a two‐dimensional human airway using POD/Galerkin reduced‐order models

An approach is proposed for the rapid prediction of nano‐particle transport and deposition in the human airway, which requires the solution of both the Navier–Stokes and advection–diffusion equations and for which computational efficiency is a challenge. The proposed method builds low‐order models that are representative of the fully coupled equations by means of the Galerkin projection and proper orthogonal decomposition technique. The obtained reduced‐order models (ROMs) are a set of ordinary differential equations for the temporal coefficients of the basis functions. The numerical results indicate that the ROMs are highly efficient for the computation (the speedup factor is approximately 3 × 10³) and have reasonable accuracy compared with the full model (relative error of ≈7 × 10^?3). Using ROMs, the deposition of particles is studied for 1≤d_n≤100 nm, where d_n is the diameter of a nano‐particle. The effectiveness of this approach is promising for applications of health risk assessment. Copyright © 2015 John Wiley & Sons, Ltd.     

An accurate and efficient three-dimensional high-order finite element methodology for the simulation of magneto-mechanical coupling in MRI scanners

Transient magnetic fields are generated by the gradient coils in an magnetic resonance imaging (MRI) scanner and induce eddy currents in their conducting components, which lead to vibrations, imaging artefacts, noise, and the dissipation of heat. Heat dissipation can boil off the helium used to cool the super conducting magnets and, if left unchecked, will lead to a magnet quench. Understanding the mechanisms involved in the generation of these vibrations, and the heat being deposited in the cryostat, are key for a successful MRI scanner design. This requires the solution of a coupled physics magneto-mechanical problem, which will be addressed in this work. A novel computational methodology is proposed for the accurate simulation of the magneto-mechanical problem using a Lagrangian approach, which, with a particular choice of linearisation, leads to a staggered scheme. This is discretised by high-order finite elements leading to accurate solutions. We demonstrate the success of our scheme by applying it to realistic MRI scanner configurations.     

Efficient uncertainty quantification of dynamical systems with local nonlinearities and uncertainties

Idealized modeling of most engineering structures yields linear mathematical models, i.e., linear ordinary or partial differential equations. However, features like nonlinear dampers and/or springs can render nonlinear an otherwise linear model. Often, the connectivity of these nonlinear elements is confined to only a few degrees-of-freedom (DOFs) of the structure. In such cases, treating the entire structure as nonlinear results in very computationally expensive solutions. Moreover, if system parameters are uncertain, their stochastic nature can render the analysis even more computationally costly. This paper presents an approach for computing the response of such systems in a very efficient manner. The proposed solution procedure first segregates the DOFs appearing in the nonlinear and/or stochastic terms from those DOFs that involve only linear deterministic operations. Second, the responses of nonlinear/stochastic terms are determined using a non-standard form of a nonlinear Volterra integral equation (NVIE). Finally, the responses of the remaining DOFs are computed through a convolution approach using the fast Fourier transform to further increase the computational efficiency. Three examples are presented to demonstrate the efficacy and accuracy of the proposed method. It is shown that, even for moderately sized systems (∼1000 DOFs), the proposed method is about three orders of magnitude faster than a conventional Monte Carlo sampling method (i.e., solving the system of ODEs repeatedly).     

Variational‐based modeling of micro‐electro‐elasticity with electric field‐driven and stress‐driven domain evolutions

Recently, increasing interest in so‐called functional or smart materials with electromechanical coupling has been shown such as ferroelectric piezoceramics. These materials are characterized by microstructural properties, which can be changed by external stress and electric field stimuli, and hence find use as the active components in sensors and actuators. The electromechanical coupling effects result from the existence and rearrangement of microstructural domains with uniformly oriented electric polarization. The understanding and efficient simulation of these highly nonlinear and dissipative mechanisms, which occur on the microscale of ferroelectric piezoceramics, are a key challenge of the current research. This paper does not offer a substantially new physical model of these phenomena but a new mathematical modeling approach based on a rigorous exploitation of rate‐type variational principles. This provides a new insight in the structure of the coupled problem, where the governing field equations appear as the Euler equations of a variational statement. We outline a variational‐based micro‐electro‐elastic model for the microstructural evolution of both electrically and mechanically driven electric domains in ferroelectric ceramics, which also incorporates the surrounding free space. To this end, we extend recently developed multifield incremental variational principles of electromechanics from local to gradient‐extended dissipative response and specialize it by a Ginzburg–Landau‐type phase field model, where the thickness of the domain walls enters the formulation as a length scale. This serves as a natural starting point for a canonical compact, symmetric finite element implementation, considering the mechanical displacement, the microscopic polarization, and the electric potential induced by the polarization as the primary fields. The latter is defined on both the solid domain and a surrounding free space. Numerical simulations treat domain wall motions for electric field‐driven and stress‐driven loading processes, including the expansion of the electric potential into the free space. Copyright © 2012 John Wiley & Sons, Ltd.     

Theory‐Guided Synthesis of a Metastable Lead‐Free Piezoelectric Polymorph

Many technologically critical materials are metastable under ambient conditions, yet the understanding of how to rationally design and guide the synthesis of these materials is limited. This work presents an integrated approach that targets a metastable lead‐free piezoelectric polymorph of SrHfO₃. First‐principles calculations predict that the previous experimentally unrealized, metastable P4mm phase of SrHfO₃ should exhibit a direct piezoelectric response (d₃₃) of 36.9 pC N^?1 (compared to d₃₃ = 0 for the ground state). Combining computationally optimized substrate selection and synthesis conditions lead to the epitaxial stabilization of the polar P4mm phase of SrHfO₃ on SrTiO₃. The films are structurally consistent with the theory predictions. A ferroelectric‐induced large signal effective converse piezoelectric response of 5.2 pm V^?1 for a 35 nm film is observed, indicating the ability to predict and target multifunctionality. This illustrates a coupled theory‐experimental approach to the discovery and realization of new multifunctional polymorphs.     

A two‐dimensional ordinary,state‐based peridynamic model for linearly elastic solids

Peridynamics is a non‐local mechanics theory that uses integral equations to include discontinuities directly in the constitutive equations. A three‐dimensional, state‐based peridynamics model has been developed previously for linearly elastic solids with a customizable Poisson's ratio. For plane stress and plane strain conditions, however, a two‐dimensional model is more efficient computationally. Here, such a two‐dimensional state‐based peridynamics model is presented. For verification, a 2D rectangular plate with a round hole in the middle is simulated under constant tensile stress. Dynamic relaxation and energy minimization methods are used to find the steady‐state solution. The model shows m‐convergence and δ‐convergence behaviors when m increases and δ decreases. Simulation results show a close quantitative matching of the displacement and stress obtained from the 2D peridynamics and a finite element model used for comparison. Copyright © 2014 John Wiley & Sons, Ltd.     

Symmetric coupling of multi‐zone curved Galerkin boundary elements with finite elements in elasticity

The coupling of Finite Element Method (FEM) with a Boundary Element Method (BEM) is a desirable result that exploits the advantages of each. This paper examines the efficient symmetric coupling of a Symmetric Galerkin Multi‐zone Curved Boundary Element Analysis method with a Finite Element Method for 2‐D elastic problems. Existing collocation based multi‐zone boundary element methods are not symmetric. Thus, when they are coupled with FEM, it is very difficult to achieve symmetry, increasing the computational work to solve the problem. This paper uses a fully Symmetric curved Multi‐zone Galerkin Boundary Element Approach that is coupled to an FEM in a completely symmetric fashion. The symmetry is achieved by symmetrically converting the boundary zones into equivalent ‘macro finite elements’, that are symmetric, so that symmetry in the coupling is retained. This computationally efficient and fast approach can be used to solve a wide range of problems, although only 2‐D elastic problems are shown. Three elasticity problems, including one from the FEM‐BEM literature that explore the efficacy of the approach are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

A coupled hp‐finite element scheme for the solution of two‐dimensional electrostrictive materials

As part of the ongoing research within the field of computational analysis for the coupled electro‐magneto‐mechanical response of smart materials, the problem of linearised electrostriction is revisited and analysed for the first time using the framework of hp‐finite elements. The governing equations modelling the physics of the dielectric are suitably modified by introducing a new total Cauchy stress tensor (A. Dorfmann and R.W. Ogden. Nonlinear electroelasticity. Acta Mechanica, 174:167–183, 2005), which includes the electrostrictive effect and a staggered partitioned scheme for the numerical solution of the coupling phenomena. With the purpose of benchmarking numerical results, the problem of an infinite electrostrictive plate with a circular/elliptical dielectric insert is revisited. The presented analytical solution is based on the theoretical framework for two‐dimensional electrostriction proposed by Knops (R.J. Knops. Two‐dimensional electrostriction. Quarterly Journal of Mechanics and Applied Mathematics, 16:377–388, 1963) and uses classical techniques of complex variable analysis. Our presentation, to the best of our knowledge, provides the first correct closed form expression for the solution to the infinite electrostrictive plate with a circular/elliptical dielectric insert, correcting the errors made in previous presentations of this problem. We use this analytical solution to assess the accuracy, efficiency and robustness of the hp‐formulation in the case of nearly incompressible electrostrictive materials. Copyright © 2012 John Wiley & Sons, Ltd.     

A posteriori finite element bounds to linear functional outputs of the three‐dimensional Navier–Stokes equations

An implicit a posteriori finite element error estimation method is presented to inexpensively calculate lower and upper bounds for a linear functional output of the numerical solutions to the three‐dimensional Navier–Stokes (N–S) equations. The novelty of this research is to utilize an augmented Lagrangian based on a coarse mesh linearization of the N–S equations and the finite element tearing and interconnecting (FETI) procedure. The latter approach extends the a posteriori bound method to the three‐dimensional Crouzeix–Raviart space for N–S problems. The computational advantage of the bound procedure is that a single coupled non‐symmetric large problem can be decomposed into several uncoupled symmetric small problems. A simple model problem, which is selected here to illustrate the procedure, is to find the lower and upper bounds of the average velocity of a pressure driven, incompressible, steady Newtonian fluid flow moving at low Reynolds numbers through an endless square channel which has an array of rectangular obstacles. Numerical results show that the bounds for this output are rigorous, i.e. always in the asymptotic certainty regime, that they are sharp and that the required computational resources decrease significantly. Parallel implementation on a Beowulf cluster is also reported. Copyright © 2004 John Wiley & Sons, Ltd.     

Variational multiscale enrichment for modeling coupled mechano‐diffusion problems

In this paper, a new multiscale–multiphysics computational methodology is devised for the analysis of coupled diffusion–deformation problems. The proposed methodology is based on the variational multiscale principles. The basic premise of the approach is accurate fine‐scale representation at a small subdomain where it is known a priori that important physical phenomena are likely to occur. The response within the remainder of the problem domain is idealized on the basis of coarse‐scale representation. We apply this idea to evaluate a coupled mechano‐diffusion problem that idealizes the response of titanium structures subjected to a thermo–chemo–mechanical environment. The proposed methodology is used to devise a multiscale model in which the transport of oxygen into titanium is modeled as a diffusion process, whereas the mechanical response is idealized using concentration‐dependent elasticity equations. A coupled solution strategy based on operator split is formulated to evaluate the coupled multiphysics and multiscale problem. Numerical experiments are conducted to assess the accuracy and computational performance of the proposed methodology. Numerical simulations indicate that the variational multiscale enrichment has reasonable accuracy and is computationally efficient in modeling the coupled mechano‐diffusion response. Copyright © 2011 John Wiley & Sons, Ltd.     

Fourier‐based numerical approximation of the Weertman equation for moving dislocations

This work addresses the numerical approximation of solutions to a dimensionless form of the Weertman equation, which models a steadily moving dislocation and is an important extension (with advection term) of the celebrated Peierls‐Nabarro equation for a static dislocation. It belongs to the class of nonlinear reaction‐advection‐diffusion integro‐differential equations with Cauchy‐type kernel, thus involving an integration over an unbounded domain. In the Weertman problem, the unknowns are the shape of the core of the dislocation and the dislocation velocity. The proposed numerical method rests on a time‐dependent formulation that admits the Weertman equation as its long‐time limit. Key features are (1) time iterations are conducted by means of a new, robust, and inexpensive Preconditioned Collocation Scheme in the Fourier domain, which allows for explicit time evolution but amounts to implicit time integration, thus allowing for large time steps; (2) as the integration over the unbounded domain induces a solution with slowly decaying tails of important influence on the overall dislocation shape, the action of the operators at play is evaluated with exact asymptotic estimates of the tails, combined with discrete Fourier transform operations on a finite computational box of size L; (3) a specific device is developed to compute the moving solution in a comoving frame, to minimize the effects of the finite‐box approximation. Applications illustrate the efficiency of the approach for different types of nonlinearities, with systematic assessment of numerical errors. Converged numerical results are found insensitive to the time step, and scaling laws for the combined dependence of the numerical error with respect to L and to the spatial step size are obtained. The method proves fast and accurate and could be applied to a wide variety of equations with moving fronts as solutions; notably, Weertman‐type equations with the Cauchy‐type kernel replaced by a fractional Laplacian.     

POD–ISAT: An efficient POD‐based surrogate approach with adaptive tabulation and fidelity regions for parametrized steady‐state PDE discrete solutions

A combination of proper orthogonal decomposition (POD) analysis and in situ adaptive tabulation (ISAT) is proposed for the representation of parameter‐dependent solutions of coupled partial differential equation problems. POD is used for the low‐order representation of the spatial fields and ISAT for the local representation of the solution in the design parameter space. The accuracy of the method is easily controlled by free threshold parameters that can be adjusted according to user needs. The method is tested on a coupled fluid‐thermal problem: the design of a simplified aircraft air control system. It is successfully compared with the standard POD; although the POD is inaccurate in certain areas of the design parameters space, the POD–ISAT method achieves accuracy thanks to trust regions based on residuals of the fluid‐thermal problem. The presented POD–ISAT approach provides flexibility, robustness and tunable accuracy to represent solutions of parametrized partial differential equations.Copyright © 2013 John Wiley & Sons, Ltd.     

Stochastic response determination of nonlinear structural systems with singular diffusion matrices: A Wiener path integral variational formulation with constraints

The Wiener path integral (WPI) approximate semi-analytical technique for determining the joint response probability density function (PDF) of stochastically excited nonlinear oscillators is generalized herein to account for systems with singular diffusion matrices. Indicative examples include (but are not limited to) systems with only some of their degrees-of-freedom excited, hysteresis modeling via additional auxiliary state equations, and energy harvesters with coupled electro-mechanical equations. In general, the governing equations of motion of the above systems can be cast as a set of underdetermined stochastic differential equations coupled with a set of deterministic ordinary differential equations. The latter, which can be of arbitrary form, are construed herein as constraints on the motion of the system driven by the stochastic excitation. Next, employing a semi-classical approximation treatment for the WPI yields a deterministic constrained variational problem to be solved numerically for determining the most probable path; and thus, for evaluating the system joint response PDF in a computationally efficient manner. This is done in conjunction with a Rayleigh-Ritz approach coupled with appropriate optimization algorithms. Several numerical examples pertaining to both linear and nonlinear constraint equations are considered, including various multi-degree-of-freedom systems, a linear oscillator under earthquake excitation and a nonlinear oscillator exhibiting hysteresis following the Bouc–Wen formalism. Comparisons with relevant Monte Carlo simulation data demonstrate a relatively high degree of accuracy.     

Efficient thermo‐mechanical model for solidification processes

A new, computationally efficient algorithm has been implemented to solve for thermal stresses, strains, and displacements in realistic solidification processes which involve highly nonlinear constitutive relations. A general form of the transient heat equation including latent‐heat from phase transformations such as solidification and other temperature‐dependent properties is solved numerically for the temperature field history. The resulting thermal stresses are solved by integrating the highly nonlinear thermo‐elastic‐viscoplastic constitutive equations using a two‐level method. First, an estimate of the stress and inelastic strain is obtained at each local integration point by implicit integration followed by a bounded Newton–Raphson (NR) iteration of the constitutive law. Then, the global finite element equations describing the boundary value problem are solved using full NR iteration. The procedure has been implemented into the commercial package Abaqus (Abaqus Standard Users Manuals, v6.4, Abaqus Inc., 2004) using a user‐defined subroutine (UMAT) to integrate the constitutive equations at the local level. Two special treatments for treating the liquid/mushy zone with a fixed grid approach are presented and compared. The model is validated both with a semi‐analytical solution from Weiner and Boley (J. Mech. Phys. Solids 1963; 11 :145–154) as well as with an in‐house finite element code CON2D (Metal. Mater. Trans. B 2004; 35B (6):1151–1172; Continuous Casting Consortium Website. http://ccc.me.uiuc.edu [30 October 2005]; Ph.D. Thesis, University of Illinois, 1993; Proceedings of the 76th Steelmaking Conference, ISS, vol. 76, 1993) specialized in thermo‐mechanical modelling of continuous casting. Both finite element codes are then applied to simulate temperature and stress development of a slice through the solidifying steel shell in a continuous casting mold under realistic operating conditions including a stress state of generalized plane strain and with actual temperature‐dependent properties. Other local integration methods as well as the explicit initial strain method used in CON2D for solving this problem are also briefly reviewed and compared. Copyright © 2006 John Wiley & Sons, Ltd.     

Dual‐grid‐based tree/cotree decomposition of higher‐order interpolatory H(∇∧, Ω) basis

This work extends the zeroth‐order tree/cotree (TC) decomposition method into higher order (HO) interpolatory elements and develops the constraints operator required for the elimination of spurious solutions for general HO spectral basis. Earlier methods explicitly enforce the divergence condition that requires a mixed finite element (FE) formulation with both H¹ and H(?∧) expansions and involves repeated solutions of the Poisson equation. A recent approach, which avoids the mixed formulation and the Poisson problem, uses TC decomposition of edge DoF over the primal graph and construction of integration and gradient matrices. The approach is easily applied to HO hierarchical elements but becomes quite complex for HO spectral elements. In the presence of internal DoF, it is difficult to utilize the primal graph for an explicit decomposition of the spectral DoF. In contrast, this work utilizes the dual grid, resulting in an explicit decomposition of DoF and construction of constraint equations from a fixed element matrix. Thus, mixed formulation and the Poisson problems are avoided while eliminating the need for evaluation of integration and gradient matrices. The proposed constraints matrix is element‐geometry independent and possesses an explicit sparsity formulation reducing the need for dynamic memory allocation. Numerical examples are included for verification. Copyright © 2010 John Wiley & Sons, Ltd.     

Thermo‐hydro‐mechanical modeling of fluid geological storage by Godunov‐mixed methods

The efficient solution to the coupled system of PDEs governing the mass and the energy balance in deformable porous media requires advanced numerical algorithms. A combination of mixed/Galerkin finite elements and finite volumes along with a staggered method are employed. Fluid flow and heat transfer are addressed iteratively by a fully coupled approach and the medium deformation by an explicitly coupled scheme, at each time step. Such formulation allows for stable numerical solutions, element‐wise conservative velocity fields and accurate prediction of sharp temperature convective fronts. The proposed model is experimented with in realistic applications of a deep aquifer fluid injection. Copyright © 2012 John Wiley & Sons, Ltd.