Comparison of analytical and numerical solution for optimal vortex diffuser design

The contribution deals with optimization of wing tip devices, so called vortex diffusers. A comparison is given between an analytical approach for obtaining the optimal circulation loading and the results of a numerical investigation using a lifting line method. The purpose of most wing tip devices is to reduce the induced drag of the main wing by converting vortex energy into thrust. In order to achieve an optimal design, a variational formulation originally proposed by Betz and Prandtl for air screws is applied to the circulation distribution of the diffuser blades. In extension to the inviscid formulation, a viscous correction is applied in order to account for frictional forces. In an effort to validate the analytical results, a comparison is given with numerical solutions from a lifting line method. The loading of the diffuser blades is parametrized and optimized with respect to resulting thrust by use of a quasi-Newton gradient method. Comparison shows that, knowing the velocity distribution in the near wake of the wing, considerable decrease of induced drag may be achieved making use of vortex diffusers. Although actual circulation loading may differ between the analytical and numerical estimation, resulting thrust agrees within a few percent. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Comparison of numerical homogenization techniques in BT-CFO multiferroic particle composites

This work is restricted to linear material behavior, i.e. the structure is considered to be in a perfectly poled state. Different numerical homogenization methods are investigated and used to calculate effective properties of a 0-2 composite modelled in reprensentative volume elements. Bariumtitanate (BT) and cobaltferrite (CFO) are employed in the Finite Element model, where the roles of matrix and inclusion are mutable in principle. Mixed magnetoelectromechanical boundary conditions based on different homogenization theories are applied to the model. The calculated macroscopic behaviors described by the different approaches are compared and presented in the paper. The special focus is on the prediction of coefficients of magnetoelectric coupling with respect to an optimization of the structural arrangement of the composite. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the numerical analysis of analytical nodal methods

This article presents the basic numerical analysis of the analytical nodal methods, which were originally developed in the late 1970s in relation with static and dynamic nuclear reactor calculations but are actually applicable to the numerical solution of partial differential equations (PDEs) in general, over fairly regular domains. The basic idea consists in “transverse integrating” the original PDE over all the variables minus one, leading to sets of 1D equations which are then solved in an “analytical” way, using fundamentals as well as particular solutions of the corresponding 1D operators. After examining the existing analytical methods in a critical way, we propose a more satisfactory extended analytical formalism. Superconvergence results finally lead us to useful conclusions with respect to the choice of a particular scheme.     

Comparison of several numerical optimization methods

Several numerical methods for solving the nonlinear two-point boundary value problem associated with an optimum spacecraft trajectory are considered. A comparative evaluation of the methods is made to determine the relative merits of each method. Particular attention is given to such characteristics as the simplicity of formulation and implementation, the convergence sensitivity, the computing time required, and the computer storage requirements. The methods considered are the perturbation method, the quasilinearization method, and the gradient method. The numerical comparison is made by considering a two-dimensional, low-thrust, minimum-time, Earth-Mars trajectory.The authors are greatly indebted to Mr. Robert D. Witty, Lockheed Electronics Company, for providing the excellent programming support.     

Multiscale methods for elliptic homogenization problems

In this article we study two families of multiscale methods for numerically solving elliptic homogenization problems. The recently developed multiscale finite element method [Hou and Wu, J Comp Phys 134 (1997), 169–189] captures the effect of microscales on macroscales through modification of finite element basis functions. Here we reformulate this method that captures the same effect through modification of bilinear forms in the finite element formulation. This new formulation is a general approach that can handle a large variety of differential problems and numerical methods. It can be easily extended to nonlinear problems and mixed finite element methods, for example. The latter extension is carried out in this article. The recently introduced heterogeneous multiscale method [Engquist and Engquist, Comm Math Sci 1 (2003), 87–132] is designed for efficient numerical solution of problems with multiscales and multiphysics. In the second part of this article, we study this method in mixed form (we call it the mixed heterogeneous multiscale method). We present a detailed analysis for stability and convergence of this new method. Estimates are obtained for the error between the homogenized and numerical multiscale solutions. Strategies for retrieving the microstructural information from the numerical solution are provided and analyzed. Relationship between the multiscale finite element and heterogeneous multiscale methods is discussed. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Tutorial on analytical methods as a complement to numerical computing

The importance of using analytical methods is taught through the discussion of an example, where the analytical treatment of a partial differential equation provides not only a suitable time scale and an asymptotic solution, but also information important for the accuracy of the numerical solution.Dedicated to Peter Naur on the occasion of his 60th birthday     

Sensitivity analysis and Taylor expansions in numerical homogenization problems

Summary. Two-scale numerical homogenization problems are addressed, with particular application to the modified compressible Reynolds equation with periodic roughness. It is shown how to calculate sensitivities of the homogenized coefficients that come out from local problems. This allows for significant reduction of the computational cost by two means: The construction of accurate Taylor expansions, and the implementation of rapidly convergent nonlinear algorithms (such as Newton's) instead of fixed-point-like ones. Numerical tests are reported showing the quantitative accuracy of low-order Taylor expansions in practical cases, independently of the shape and smoothness of the roughness function. Received December 8, 1997 / Published online January 27, 2000     

Multiscale method based on discontinuous Galerkin methods for homogenization problems

We propose a multiscale method for elliptic problems with highly oscillating coefficients based on a coupling of macro and micro methods in the framework of the heterogeneous multiscale method. The macro method, defined on a macroscopic triangulation, aims at recovering the effective (homogenized) solution of an unknown macro model. The unspecified data of this model are computed by micro methods on sampling domains during the macro assembly process. In this Note, we show how to construct such a coupling with a discontinuous macro finite element space. We show that the flux information needed in this formulation in order to impose weak interelement continuity can be recovered from the known micro calculations on the sampling domains. A fully discrete analysis is presented. To cite this article: A. Abdulle, C. R. Acad. Sci. Paris, Ser. I 346 (2008).     

Order of convergence for polynomial and analytical nodal methods

In [6] we analyzed the direct analytical nodal methods (ANM) of indexl and show that the corresponding mathematical methods are equivalent to the physical ones when the components of the matrices are calculated by generalized Radau reduced integration. In this article we extend the theorem 8 of [7] to the polynomial nodal methods (PNM) (exact calculation of moments) which are thus the order ofO(h ^l+3?δ _l0. We also show that the analytical nodal methods are only the order ofO(h ^l+2). Forl = 0 our numerical results confirm our theoretical results.     

Simplified analytical expressions for numerical differentiation via cycle index

Based on the Lagrange interpolation to the function f[_x0,⋅]$f[x_0,\cdot ]$ for arbitrarily chosen _x0$x_0$ and logarithmic differentiation, we give a simple approach to analytical expressions for numerical differentiation using cycle index. A detailed analysis for the remainder is also included.     

Mean and full field homogenization of artificial long fiber reinforced thermoset polymers

Based on two artificial microstructures representing a long fiber reinforced thermoset material, the effective linear elastic material properties are calculated by both a mean and a full field homogenization method. Concerning the mean field method, the effective elastic material properties are approximated using the homogenization scheme by Mori and Tanaka, formulated explicitly in terms of orientation averages. This allows to use orienation tensors of 2nd and 4th order describing the orientation information on the micro level. The full field method is based on the fast Fourier transformation (FFT), for which the effective material properties are determined by volume averaging. The comparison between both methods show good agreements, the deviations are in the range between 2% and 12%. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Novel analytical and numerical techniques for fractional temporal SEIR measles model

In this paper, a fractional temporal SEIR measles model is considered. The model consists of four coupled time fractional ordinary differential equations. The time-fractional derivative is defined in the Caputo sense. Firstly, we solve this model by solving an approximate model that linearizes the four time fractional ordinary differential equations (TFODE) at each time step. Secondly, we derive an analytical solution of the single TFODE. Then, we can obtain analytical solutions of the four coupled TFODE at each time step, respectively. Thirdly, a computationally effective fractional Predictor-Corrector method (FPCM) is proposed for simulating the single TFODE. And the error analysis for the fractional predictor-corrector method is also given. It can be shown that the fractional model provides an interesting technique to describe measles spreading dynamics. We conclude that the analytical and Predictor-Corrector schemes derived are easy to implement and can be extended to other fractional models. Fourthly, for demonstrating the accuracy of analytical solution for fractional decoupled measles model, we applied GMMP Scheme (Gorenflo-Mainardi-Moretti-Paradisi) to the original fractional equations. The comparison of the numerical simulations indicates that the solution of the decoupled and linearized system is close enough to the solution of the original system. And it also indicates that the linearizing technique is correct and effective.     

Mean field homogenization and experimental investigation of short and long fiber reinforced polymers

The effective elastic material properties of short fiber reinforced polypropylen are determined by means of the self-consistent (SC) method and the interaction direct derivative (IDD) method. In order to account for thermoelastic effective material properties, a Hashin-Shtrikman (HS) based two-step homogenization method with variable reference stiffness is used. The influence of the reference stiffness, dependent on a scalar parameter is investigated. Information on the microstructure are derived by computed tomography scans (µCT) and considered within the homogenization schemes. Thermomechanical properties of a long fiber reinforced polymer (LFRP) and a short fiber reinforced polymer (SFRP) are obtained by means of dynamic mechanical analysis (DMA). Simulation results for SFRP are compared to experimental results. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

不可压缩流动的并行数值方法

不可压缩流动的数值模拟是计算流体力学的重要组成部分. 基于有限元离散方法, 本文设计了不可压缩Navier-Stokes (N-S)方程支配流的若干并行数值算法. 这些并行算法可归为两大类: 一类是基于两重网格离散方法, 首先在粗网格上求解非线性的N-S方程, 然后在细网格的子区域上并行求解线性化的残差方程, 以校正粗网格的解; 另一类是基于新型完全重叠型区域分解技巧, 每台处理器用一局部加密的全局多尺度网格计算所负责子区域的局部有限元解. 这些并行算法实现简单, 通信需求少, 具有良好的并行性能, 能获得与标准有限元方法相同收敛阶的有限元解. 理论分析和数值试验验证了并行算法的高效性     

Multistep numerical methods for functional-differential-algebraic equations

This paper deals with constructing multistep numerical methods for differential delay equations under additional algebraic constraints. Theorems on the orders of convergence of these methods are proved both for functional-differential equations with algebraic constraints and for singular functional-differential equations.     

Hermite spectral and pseudospectral methods for numerical differentiation

A numerical differentiation problem for a given function with noisy data is discussed in this paper. A mollification method based on spanned by Hermite functions is proposed and the mollification parameter is chosen by a discrepancy principle. The convergence estimates of the derivatives are obtained. To get a practical approach, we also derive corresponding results for pseudospectral (Hermite-Gauss interpolation) approximations. Numerical examples are given to show the efficiency of the method.