Node-based smoothed radial point interpolation method for electromagnetic-thermal coupled analysis

The node-based smoothed radial point interpolation method for solving the transient responses of magneto-electro-elastic structures in thermal environment is proposed. Considering the coupling relations between the elasticity, magnetism, electricity and heat, the generalized displacement (displacement, electric potential and magnetic potential) is calculated using the modified Newmark method. G space theory and the weakened weak formulation are applied to derive the equations of node-based smoothed radial point interpolation method for the magneto-electro-thermo-elastic multi-physics coupling problems. We use triangular background meshes as they could be generated more easily for structures with complex geometry. In some cases, they could even be created automatically. Detailed numerical study has shown that node-based smoothed radial point interpolation method not only successfully overcomes the overly-stiff behavior in the FEM and provides more accurate results, but also works well with distorted meshes. Therefore, the node-based smoothed radial point interpolation method could be adopted to solve the magneto-electro-thermo-elastic multi-physics coupling problems in practical engineering.     

A model-coupling framework for nearshore waves,currents, sediment transport,and seabed morphology

This paper presents a framework for synchronously coupling wave, current, sediment transport, and seabed morphology for the accurate simulation of multi-physics coastal ocean processes. The governing equations, which represent models that are commonly adopted in practical simulations, are discretized using finite-difference methods. The resulting system is validated against analytical solutions. In order to test the performance of the proposed framework and the numerical methods, dam-break flow over a mobile-bed and evolution of a wave-driven sand dune are simulated. The interactions among waves, currents, and seabed morphology are illustrated.     

Fundamental elastic field in an infinite plane of two-dimensional piezoelectric quasicrystal subjected to multi-physics loads

By utilizing the extended Stroh formalism, the Green's function of infinite plane is obtained for the problem of two-dimensional decagonal quasicrystals with the piezoelectric effect subjected to multi-physics loads. By numerical computations, the piezoelectric effect of the two-dimensional decagonal quasicrystals is revealed; the changes of the stress and displacement fields with multi-physics loads are discussed. The variation laws of material constants in stress and displacement fields are investigated. The results show that the effect of the phason field on the generalized displacement is larger than that on the generalized stress; and the effects of material parameters are different in diverse field.     

Parallel multi-frontal solver for p adaptive finite element modeling of multi-physics computational problems

The paper presents a parallel direct solver for multi-physics problems. The solver is dedicated for solving problems resulting from adaptive finite element method computations. The concept of finite element is actually replaced by the concept of the node. The computational mesh consists of several nodes, related to element vertices, edges, faces and interiors. The ordering of unknowns in the solver is performed on the level of nodes. The concept of the node can be efficiently utilized in order to recognize unknowns that can be eliminated at a given node of the elimination tree. The solver is tested on the exemplary three-dimensional multi-physics problem involving the computations of the linear acoustics coupled with linear elasticity. The three-dimensional tetrahedral mesh generation and the solver algorithm are modeled by using graph grammar formalism. The execution time and the memory usage of the solver are compared with the MUMPS solver.     

An adaptive integration scheme for heat conduction in additive manufacturing

Solidification dynamics are important for determining final microstructure in additively manufactured parts. Recently, researchers have adopted semi-analytical approaches for predicting heat conduction effects at length and time scales not accessible to complex multi-physics numerical models. The present work focuses on improving a semi-analytical heat conduction model for additive manufacturing by designing an adaptive integration technique. The proposed scheme considers material properties, process conditions, and the inherent physical behavior of the transient heat conduction around both stationary and moving heat sources. Both the adaptive integration scheme and a technique for calculating only the points within the melt pools are described in detail. The full algorithm is then implemented and compared against a simple Riemann sum integration scheme for a variety of cases that highlight process and material variations relevant to additive manufacturing. The new scheme is shown to have significant improvements in computational efficiency, solution accuracy, and usability.     

Construction of a high order fluid-structure interaction solver

Accuracy is critical if we are to trust simulation predictions. In settings such as fluid-structure interaction, it is all the more important to obtain reliable results to understand, for example, the impact of pathologies on blood flows in the cardiovascular system. In this paper, we propose a computational strategy for simulating fluid structure interaction using high order methods in space and time.First, we present the mathematical and computational core framework, Life, underlying our multi-physics solvers. Life is a versatile library allowing for 1D, 2D and 3D partial differential solves using h/p type Galerkin methods. Then, we briefly describe the handling of high order geometry and the structure solver. Next we outline the high-order space-time approximation of the incompressible Navier-Stokes equations and comment on the algebraic system and the preconditioning strategy. Finally, we present the high-order Arbitrary Lagrangian Eulerian (ALE) framework in which we solve the fluid-structure interaction problem as well as some initial results.     

EWE: Toward electro-mechanical cardiac simulations with MOOSE

We present the software framework EWE, which is designed for coupled electromechanical simulations in computational cardiology. EWE is build on the multi-physics framework MOOSE. Numerical simulations of coupled problems on an idealized geometry for a left ventricle are shown. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An experimental study of interface relaxation methods for composite elliptic differential equations

The recently proposed simulation framework of interface relaxation for developing multi-domain multi-physics simulation engines is considered. An experimental study of the behavior of two representative interface relaxation methods is presented. Three linear and one non-linear elliptic two-dimensional PDE problems are considered and they are coupled with both cartesian and general decompositions. The characteristics and the effectiveness of the proposed collaborative PDE solving framework in general, and of the two interface relaxation methods in particular are shown.     

Thermo-Mechanical FE2 Simulation Scheme for Abaqus

A fully coupled multi-physics (thermo-mechanical) multi-length scale (macro-micro) finite element based simulation scheme for the commercial FE code Abaqus is introduced. The coupling of scales, initial boundary value problems on the macro- and micro-scale as well as the micro-scale testing procedure, which completes the consistent loop across the scales, are outlined. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

关于M/M/1动态排队系统解的半稳定性

本文讨论了两类 M/M/1 动态系统的数学模型 ,利用常微分方程所描述的 M/M/1 系统的结果证明了较复杂的偏微分方程所描述的 M/M/1 系统的一些性质 ,该方法简化了已有结果     

Necessary and sufficient conditions for the existence of symmetry groups for systems of differential-difference equations with partial derivatives

In the present paper, we obtain necessary and sufficient conditions for the existence of symmetry groups for systems of differential-difference equations with partial derivatives. Constructions of symmetry groups and recursion operators for various differential-difference equations are described. The corresponding symmetry generators that transform solutions of given equations into other solutions are obtained.     

Auxiliary equation method for solving nonlinear Wick-type partial differential equations

Using the solutions of an auxiliary differential equation, a direct algebraic method is described to construct several kinds of exact travelling wave solutions for some Wick-type nonlinear partial differential equations. By this method some physically important nonlinear equations are investigated and new exact travelling wave solutions are explicitly obtained. In addition, the links between Wick-type partial differential equations and variable coefficient partial differential equations are also clarified generally.     

High order nonlinear parabolic equations

The basic results and methods of the theory of high order nonlinear parabolic equations are described. In the first chapter boundary problems for quasilinear parabolic equations having divergent form are considered. In the second chapter nonlinear parabolic equations of general form are considered. Attention is mainly paid to methods of study of nonlinear parabolic problems. In particular, the methods of monotonicity and compactness, the method of a priori estimates, the functional-analytic method, etc. are described.Translated from Itogi Nauki i Tekhniki, Seriya Sovremennye Problemy Matematiki, Noveishie Dostizheniya, Vol. 37, pp. 89–166, 1990.     

Dimension reduction in fluid dynamics equations

A method for transforming the Euler and Navier-Stokes equations and a complete system of fluid dynamics equations in three dimensions to a closed system on any moving surface is proposed. As a result, for an arbitrary geometric configuration, the dimension of the equations is reduced by one, which makes them convenient for numerical simulation. The general principles of the method are described, and verifying examples are presented.     

On RF-pairs, Bäcklund transformations and linearization of nonlinear equations

Some classes of nonlinear equations of mathematical physics are described that admit order reduction through the use of a hydrodynamic-type transformation, where the unknown function is taken as a new independent variable and an appropriate partial derivative is taken as the new dependent variable. RF-pairs and associated Bäcklund transformations are constructed for evolution equations of general form. The results obtained are used for order reduction of hydrodynamic equations (Navier-Stokes and boundary layer) and constructing exact solutions to these equations. A generalized Calogero equation and a number of other new linearizable nonlinear differential equations of the second, third and forth orders are considered. Some integro-differential equations are analyzed.     

Structural theory of special functions

A block diagram is suggested for classifying differential equations whose solutions are special functions of mathematical physics. Three classes of these equations are identified: the hypergeometric, Heun, and Painlevé classes. The constituent types of equations are listed for each class. The confluence processes that transform one type into another are described. The interrelations between the equations belonging to different classes are indicated. For example, the Painlevé-class equations are equations of classical motion for Hamiltonians corresponding to Heun-class equations, and linearizing the Painlevé-class equations leads to hypergeometric-class equations. The “confluence principle” is stated, and an example of its application is given. Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 119, No. 1, pp. 3–19, April, 1999.     

The One-Shot Method in SU2

For PDE-constrained optimization tasks the one-shot approach is well-suited. The solution of the state equation is augmented with an adjoint solver to provide approximate reduced derivatives that can be used in an optimization iteration to change the design. In this paper different one-shot strategies for SU2 are implemented and tested for aerodynamic shape optimization in the open-source multi-physics package SU2 to extend the choice of optimization algorithms. These include the application of the one-shot approach when considering additional constraints apart from the state equation. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A general adjoint relation for functional differential and Volterra integral equations with application to control

A general adjoint relation is developed between solutions of linear functional differential equations and linear Volterra integral equations. Several useful representations for solutions of such equations arise as a consequence of the adjoint relationship. These representations are then used to obtain directly several results for controlling systems described by either linear functional differential equations or linear Volterra integral equations.This work was supported by the National Science Foundation under Grant No. GK-5798.