Weakly coupled analysis of a blade in multiphase mixing vessel

Two or more physical systems frequently interact with each other, where the independent solution of one system is impossible without a simultaneous solution of the others. An obvious coupled system is that of a dynamic fluid-structure interaction. [8] In this paper a computational analysis of the fluid-structure interaction in a mixing vessel is presented. In mixing vessels the fluid can have a significant influence on the deformation of blades during mixing, depending on speed of mixing blades and fluid viscosity. For this purpose a computational weakly coupled analysis has been performed to determine the multiphase fluid influences on the mixing vessel structure. The multiphase fluid field in the mixing vessel was first analyzed with the computational fluid dynamics (CFD) code CFX. The results in the form of pressure were then applied to the blade model, which was the analysed with the structural code MSC.visualNastran forWindows, which is based on the finite element method (FEM). (© 2004 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Preconditioned iterative methods for coupled discretizations of fluid flow problems

Computational fluid dynamics, where simulations require largecomputation times, is one of the areas of application of highperformance computing. Schemes such as the SIMPLE (semi-implicitmethod for pressure-linked equations) algorithm are often usedto solve the discrete Navier-Stokes equations. Generally theseschemes take a short time per iteration but require a largenumber of iterations. For simple geometries (or coarser grids)the overall CPU time is small. However, for finer grids or morecomplex geometries the increase in the number of iterationsmay be a drawback and the decoupling of the differential equationsinvolved implies a slow convergence of rotationally dominatedproblems that can be very time consuming for realistic applications.So we analyze here another approach, DIRECTO, that solves theequations in a coupled way. With recent advances in hardwaretechnology and software design, it became possible to solvecoupled Navier-Stokes systems, which are more robust but implyincreasing computational requirements (both in terms of memoryand CPU time). Two approaches are described here (band blockLU factorization and preconditioned GMRES) for the linear solverrequired by the DIRECTO algorithm that solves the fluid flowequations as a coupled system. Comparisons of the effectivenessof incomplete factorization preconditioners applied to the GMRES(generalized minimum residual) method are shown. Some numericalresults are presented showing that it is possible to minimizeconsiderably the CPU time of the coupled approach so that itcan be faster than the decoupled one.     

The numerical simulation of a typical fluid-structure interaction problem

To determine the dynamic response of a structure under the influence of the fluid flow one must solve a coupled computational fluid dynamics (CFD) and computational structural dynamics (CSD) mathematical problem. This paper presents the comparison of two methods for the calculation of the fluid-structure interaction. The first one is of explicit-implicit type and uses a staggered time advancement of the fluid and structure problems. The second uses a fully implicit discretization in the physical time of the fluid-structure equations and an explicit advancement in the dual-time. The physical fluid-structure problem is accompanied by the equations of the mesh motion, which are written as for a pseudo-structural system with its own dynamics. Representative numerical results are presented for the two degrees of freedom tipical section in unsteady transonic flow. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

非一致性界面热流固耦合作用的整体求解

提出了非一致性界面热流固耦合作用整体求解的一种方法．热流体求解基于Boussinesq假设和不可压缩的Navier-Stokes方程．流体区域的运动采用任意Lagrange-Euler（ALE）方法．拟固体元方法实现流体区域的变形．使用几何非线性的热弹性动力学描述固体运动．为了保证界面处应力和传热的平衡,采用了基于Gauss积分点的数据交换方法,对热流固耦合最终形成的强非线性方程实现整体求解．数值实例分析表明该方法的健壮性和有效性．     

Numerical method of mixed finite volume-modified upwind fractional step difference for three-dimensional semiconductor device transient behavior problems

Transient behavior of three-dimensional semiconductor device with heat conduction is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an elliptic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two concentrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L² norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.     

2D simulation of fluid-structure interaction using finite element method

This paper deals with pressure-based finite element analysis of fluid–structure systems considering the coupled fluid and structural dynamics. The present method uses two-dimensional fluid elements and structural line elements for the numerical simulation of the problem. The equations of motion of the fluid, considered inviscid and compressible, are expressed in terms of the pressure variable alone. The solution of the coupled system is accomplished by solving the two systems separately with the interaction effects at the fluid–solid interface enforced by an iterative scheme. Non-divergent pressure and displacement are obtained simultaneously through iterations. The Galerkin weighted residual method-based FE formulation and the iterative solution procedure are explained in detail followed by some numerical examples. Numerical results are compared with the existing solutions to validate the code for sloshing with fluid–structure coupling.     

A Long Wave Approximation for Capillary-Gravity Waves and an Effect of the Bottom

The Korteweg–de Vries (KdV) equation is known as a model of long waves in an infinitely long canal over a flat bottom and approximates the 2-dimensional water wave problem, which is a free boundary problem for the incompressible Euler equation with the irrotational condition. In this article, we consider the validity of this approximation in the case of the presence of the surface tension. Moreover, we consider the case where the bottom is not flat and study an effect of the bottom to the long wave approximation. We derive a system of coupled KdV like equations and prove that the dynamics of the full problem can be described approximately by the solution of the coupled equations for a long time interval. We also prove that if the initial data and the bottom decay at infinity in a suitable sense, then the KdV equation takes the place of the coupled equations.     

矩形截面梁中波传播的精确分析

将傅立叶级数方法推广应用于矩形截面梁中波传播的精确分析．不仅试着直接从三维弹性动力学方程出发,导出了矩形截面梁中波传播的一般解析解,而且给出了弹性波在自由矩形截面梁中的传播特性．波传播精确模型的提出,为实现梁波的耦合控制奠定了坚实基础．     

Strong Interactions between Solitary Waves Belonging to Different Wave Modes

Strong interactions between weakly nonlinear long waves are studied. Strong interactions occur when the linear long wave phase speeds are nearly equal although the waves belong to different modes. Specifically we study this situation in the context of internal wave modes propagating in a density stratified fluid. The interaction is described by two coupled Korteweg-deVries equations, which possess both dispersive and nonlinear coupling terms. It is shown that the coupled equations possess an exact analytical solution involving the characteristic “sech²” profile of the Korteweg-deVries equation. It is also shown that when the coefficients satisfy some special conditions, the coupled equations possess an n-solition solution analogous to the Korteweg-deVries n-solition solution. In general though the coupled equations are found not to be amenable to solution by the inverse scattering transform technique, and thus a numerical method has been employed in order to find solutions. This method is described in detail in Appendix A. Several numerical solutions of the coupled equations are presented.     

A coupled far-field formulation for time-periodic numerical problems in fluid dynamics

Consider uniform flow past an oscillating body generating a time-periodic motion in an exterior domain, modelled by a numerical fluid dynamics solver in the near field around the body. A far-field formulation, based on the Oseen equations, is presented for coupling onto this domain thereby enabling the whole space to be modelled. In particular, examples for formulations by boundary elements and infinite elements are described.     

MIP-based instantaneous control of mixed-integer PDE-constrained gas transport problems

We study the transient optimization of gas transport networks including both discrete controls due to switching of controllable elements and nonlinear fluid dynamics described by the system of isothermal Euler equations, which are partial differential equations in time and 1-dimensional space. This combination leads to mixed-integer optimization problems subject to nonlinear hyperbolic partial differential equations on a graph. We propose an instantaneous control approach in which suitable Euler discretizations yield systems of ordinary differential equations on a graph. This networked system of ordinary differential equations is shown to be well-posed and affine-linear solutions of these systems are derived analytically. As a consequence, finite-dimensional mixed-integer linear optimization problems are obtained for every time step that can be solved to global optimality using general-purpose solvers. We illustrate our approach in practice by presenting numerical results on a realistic gas transport network.     

包含泥沙冲淤的浅水方程的混合有限元法(Ⅱ)——时间沿特征方向离散的全离散情形

进一步研究由水动力学方程、泥沙输运方程和河床变化方程组成的浅水方程混合有限元法,给出时间沿特征方向离散的一种全离散格式,并证明全离散的水流速度、床底高度、水体厚度、水中泥沙含量的混合有限元解的存在性和收敛性(误差估计)．     

Numerical simulation of fluid-structure interaction problems on hybrid meshes with algebraic multigrid methods

Fluid-structure interaction problems arise in many fields of application such as flows around elastic structures and blood flow in arteries. The method presented in this paper for solving such a problem is based on a reduction to an equation at the interface, involving the so-called Steklov-Poincaré operators. This interface equation is solved by a Newton iteration, for which directional derivatives involving shape derivatives with respect to the interface perturbation have to be evaluated appropriately. One step of the Newton iteration requires the solution of several decoupled linear sub-problems in the structure and the fluid domains. These sub-problems are spatially discretized by a finite element method on hybrid meshes. For the time discretization, implicit first-order methods are used for both sub-problems. The discretized equations are solved by algebraic multigrid methods.     

Computation of flow-induced motion of floating bodies

A computational procedure for the prediction of motion of rigid bodies floating in viscous fluids and subjected to currents and waves is presented. The procedure is based on a coupled iterative solution of equations of motion of a rigid body with up to six degrees of freedom and the Reynolds-averaged Navier–Stokes equations describing the two- or three-dimensional fluid flow. The fluid flow is predicted using a commercial CFD package which can use moving grids made of arbitrary polyhedral cells and allows sliding interfaces between fixed and moving grid blocks. The computation of body motion is coupled to the CFD code via user-coding interfaces. The method is used to compute the 2D motion of floating bodies subjected to large waves and the results are compared to available experimental data, showing favorable agreement.     

Numerical simulation of thermal problems coupled with magnetohydrodynamic effects in aluminium cell

A system of partial differential equations describing the thermal behavior of aluminium cell coupled with magnetohydrodynamic effects is numerically solved. The thermal model is considered as a two-phases Stefan problem which consists of a non-linear convection–diffusion heat equation with Joule effect as a source. The magnetohydrodynamic fields are governed by Navier–Stokes and by static Maxwell equations. A pseudo-evolutionary scheme (Chernoff) is used to obtain the stationary solution giving the temperature and the frozen layer profile for the simulation of the ledges in the cell. A numerical approximation using a finite element method is formulated to obtain the fluid velocity, electrical potential, magnetic induction and temperature. An iterative algorithm and 3-D numerical results are presented.     

Solutions of the Stokes,Oseen, and Navier-Stokes Equations in Three Dimensions

A constructive analytical procedure is described for generating a class of solutions pertaining to the fluid velocity field satisfying the three-dimensional Navier-Stokes equations. The arbitrary functions occurring in the solution are related to those that arise in the general solution of the Stokes and Oseen equations.     

Gradient stability of numerical algorithms in local nonequilibrium problems of critical dynamics

The critical dynamics of a spatially inhomogeneous system are analyzed with allowance for local nonequilibrium, which leads to a singular perturbation in the equations due to the appearance of a second time derivative. An extension is derived for the Eyre theorem, which holds for classical critical dynamics described by first-order equations in time and based on the local equilibrium hypothesis. It is shown that gradient-stable numerical algorithms can also be constructed for second-order equations in time by applying the decomposition of the free energy into expansive and contractive parts, which was suggested by Eyre for classical equations. These gradient-stable algorithms yield a monotonically nondecreasing free energy in simulations with an arbitrary time step. It is shown that the gradient stability conditions for the modified and classical equations of critical dynamics coincide in the case of a certain time approximation of the inertial dynamics relations introduced for describing local nonequilibrium. Model problems illustrating the extended Eyre theorem for critical dynamics problems are considered.     

Comparison of preconditioned Krylov subspace iteration methods for PDE-constrained optimization problems

The governing dynamics of fluid flow is stated as a system of partial differential equations referred to as the Navier-Stokes system. In industrial and scientific applications, fluid flow control becomes an optimization problem where the governing partial differential equations of the fluid flow are stated as constraints. When discretized, the optimal control of the Navier-Stokes equations leads to large sparse saddle point systems in two levels. In this paper, we consider distributed optimal control for the Stokes system and test the particular case when the arising linear system can be compressed after eliminating the control function. In that case, a system arises in a form which enables the application of an efficient block matrix preconditioner that previously has been applied to solve complex-valued systems in real arithmetic. Under certain conditions, the condition number of the so preconditioned matrix is bounded by 2. The numerical and computational efficiency of the method in terms of number of iterations and execution time is favorably compared with other published methods.     

On coupled problems in geomechanics

Geomechanical problems are generally based on the category of granular, cohesive-frictional materials with a fluid pore content. At the macroscopic scale of continuum mechanics, these materials can be successfully described on the basis of the well-founded Theory of Porous Media (TPM) [1]. The present contribution touches fundamental problems of coupled media by investigating the interacting behaviour of an elasto-viscoplastic porous solid skeleton, the soil, and two pore fluids, water and air. Furthermore, electro-chemical reactions are considered in order to include the swelling behaviour of active soil. In conclusion, this leads to a system of strongly coupled partial differential equations (PDE) that can be solved by use of the finite element method (FEM). In particular, the presentation includes fluid-flow situations in the fully or the partially saturated range, swelling phenomena of active clay [3] as well as localisation phenomena [2] as a result of fluid flow or heavy rainfall events. The computations are carried out by use of the single-processor FE tool PANDAS [4] and, in case of large 3-d problems, by coupling PANDAS with the multi-processor solver M++. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Galerkin and weighted Galerkin methods for a forward-backward heat equation

Summary. Galerkin and weighted Galerkin methods are proposed for the numerical solution of parabolic partial differential equations where the diffusion coefficient takes different signs. The approach is based on a simultaneous discretization of space and time variables by using continuous finite element methods. Under some simple assumptions, error estimates and some numerical results for both Galerkin and weighted Galerkin methods are presented. Comparisons with the previous methods show that new methods not only can be used to solve a wider class of equations but also require less regularity for the solution and need fewer computations. Received March 3, 1995