5,579 research outputs found

    High performance computing of explicit schemes for electrofusion jointing process based on message-passing paradigm

    Get PDF
    The research focused on heterogeneous cluster workstations comprising of a number of CPUs in single and shared architecture platform. The problem statements under consideration involved one dimensional parabolic equations. The thermal process of electrofusion jointing was also discussed. Numerical schemes of explicit type such as AGE, Brian, and Charlies Methods were employed. The parallelization of these methods were based on the domain decomposition technique. Some parallel performance measurement for these methods were also addressed. Temperature profile of the one dimensional radial model of the electrofusion process were also given

    A new hybrid implicit-explicit FDTD method for local subgridding in multiscale 2-D TE scattering problems

    Get PDF
    The conventional finite-difference time-domain (FDTD) method with staggered Yee scheme does not easily allow including thin material layers, especially so if these layers are highly conductive. This paper proposes a novel subgridding technique for 2-D problems, based on a hybrid implicit-explicit scheme, which efficiently copes with this problem. In the subgrid, the new method collocates field components such that the thin layer boundaries are defined unambiguously. Moreover, aspect ratios of more than a million do not impair the stability of the method and allow for very accurate predictions of the skin effect. The new method retains the Courant limit of the coarse Yee grid and is easily incorporated into existing FDTD codes. A number of illustrative examples, including scattering by a metal grating, demonstrate the accuracy and stability of the new method

    ULTRA-SHARP nonoscillatory convection schemes for high-speed steady multidimensional flow

    Get PDF
    For convection-dominated flows, classical second-order methods are notoriously oscillatory and often unstable. For this reason, many computational fluid dynamicists have adopted various forms of (inherently stable) first-order upwinding over the past few decades. Although it is now well known that first-order convection schemes suffer from serious inaccuracies attributable to artificial viscosity or numerical diffusion under high convection conditions, these methods continue to enjoy widespread popularity for numerical heat transfer calculations, apparently due to a perceived lack of viable high accuracy alternatives. But alternatives are available. For example, nonoscillatory methods used in gasdynamics, including currently popular TVD schemes, can be easily adapted to multidimensional incompressible flow and convective transport. This, in itself, would be a major advance for numerical convective heat transfer, for example. But, as is shown, second-order TVD schemes form only a small, overly restrictive, subclass of a much more universal, and extremely simple, nonoscillatory flux-limiting strategy which can be applied to convection schemes of arbitrarily high order accuracy, while requiring only a simple tridiagonal ADI line-solver, as used in the majority of general purpose iterative codes for incompressible flow and numerical heat transfer. The new universal limiter and associated solution procedures form the so-called ULTRA-SHARP alternative for high resolution nonoscillatory multidimensional steady state high speed convective modelling

    Efficient solution of 3D electromagnetic eddy-current problems within the finite volume framework of OpenFOAM

    Full text link
    Eddy-current problems occur in a wide range of industrial and metallurgical applications where conducting material is processed inductively. Motivated by realising coupled multi-physics simulations, we present a new method for the solution of such problems in the finite volume framework of foam-extend, an extended version of the very popular OpenFOAM software. The numerical procedure involves a semi-coupled multi-mesh approach to solve Maxwell's equations for non-magnetic materials by means of the Coulomb gauged magnetic vector potential and the electric scalar potential. The concept is further extended on the basis of the impressed and reduced magnetic vector potential and its usage in accordance with Biot-Savart's law to achieve a very efficient overall modelling even for complex three-dimensional geometries. Moreover, we present a special discretisation scheme to account for possible discontinuities in the electrical conductivity. To complement our numerical method, an extensive validation is completing the paper, which provides insight into the behaviour and the potential of our approach.Comment: 47 pages, improved figures, updated references, fixed typos, reverse phase shift, consistent use of inner produc

    Extended Smoothed Boundary Method for Solving Partial Differential Equations with General Boundary Conditions on Complex Boundaries

    Full text link
    In this article, we describe an approach for solving partial differential equations with general boundary conditions imposed on arbitrarily shaped boundaries. A continuous function, the domain parameter, is used to modify the original differential equations such that the equations are solved in the region where a domain parameter takes a specified value while boundary conditions are imposed on the region where the value of the domain parameter varies smoothly across a short distance. The mathematical derivations are straightforward and generically applicable to a wide variety of partial differential equations. To demonstrate the general applicability of the approach, we provide four examples herein: (1) the diffusion equation with both Neumann and Dirichlet boundary conditions; (2) the diffusion equation with both surface diffusion and reaction; (3) the mechanical equilibrium equation; and (4) the equation for phase transformation with the presence of additional boundaries. The solutions for several of these cases are validated against corresponding analytical and semi-analytical solutions. The potential of the approach is demonstrated with five applications: surface-reaction-diffusion kinetics with a complex geometry, Kirkendall-effect-induced deformation, thermal stress in a complex geometry, phase transformations affected by substrate surfaces, and a self-propelled droplet.Comment: This document is the revised version of arXiv:0912.1288v

    Numerical Analysis of Transient Teflon Ablation with a Domain Decomposition Finite Volume Implicit Method on Unstructured Grids

    Get PDF
    This work investigates numerically the process of Teflon ablation using a finite-volume discretization, implicit time integration and a domain decomposition method in three-dimensions. The interest in Teflon stems from its use in Pulsed Plasma Thrusters and in thermal protection systems for reentry vehicles. The ablation of Teflon is a complex process that involves phase transition, a receding external boundary where the heat flux is applied, an interface between a crystalline and amorphous (gel) phase and a depolymerization reaction which happens on and beneath the ablating surface. The mathematical model used in this work is based on a two-phase model that accounts for the amorphous and crystalline phases as well as the depolymerization of Teflon in the form of an Arrhenius reaction equation. The model accounts also for temperature-dependent material properties, for unsteady heat inputs and boundary conditions in 3D. The model is implemented in 3D domains of arbitrary geometry with a finite volume discretization on unstructured grids. The numerical solution of the transient reaction-diffusion equation coupled with the Arrhenius-based ablation model advances in time using implicit Crank-Nicolson scheme. For each time step the implicit time advancing is decomposed into multiple sub-problems by a domain decomposition method. Each of the sub-problems is solved in parallel by Newton-Krylov non-linear solver. After each implicit time-advancing step, the rate of ablation and the fraction of depolymerized material are updated explicitly with the Arrhenius-based ablation model. After the computation, the surface of ablation front and the melting surface are recovered from the scalar field of fraction of depolymerized material and the fraction of melted material by post-processing. The code is verified against analytical solutions for the heat diffusion problem and the Stefan problem. The code is validated against experimental data of Teflon ablation. The verification and validation demonstrates the ability of the numerical method in simulating three dimensional ablation of Teflon

    Modelling of impulse loading in high-temperature superconductors. Assessment of accuracy and performance of computational techniques.

    No full text
    Purpose – The aim of this paper is to access performance of existing computational techniques to model strongly non-linear field diffusion problems. Design/methodology/approach – Multidimensional application of a finite volume front-fixing method to various front-type problems with moving boundaries and non-linear material properties is discussed. Advantages and implementation problems of the technique are highlighted by comparing the front-fixing method with computations using fixed grids. Particular attention is focused on conservation properties of the algorithm and accurate solutions close to the moving boundaries. The algorithm is tested using analytical solutions of diffusion problems with cylindrical symmetry with both spatial and temporal accuracy analysed. Findings – Several advantages are identified in using a front-fixing method for modelling of impulse phenomena in high-temperature superconductors (HTS), namely high accuracy can be obtained with a small number of grid points, and standard numerical methods for convection problems with diffusion can be utilised. Approximately, first order of spatial accuracy is found for all methods (stationary or mobile grids) for 2D problems with impulse events. Nevertheless, errors resulting from a front-fixing technique are much smaller in comparison with fixed grids. Fractional steps method is proved to be an effective algorithm for solving the equations obtained. A symmetrisation procedure has to be introduced to eliminate a directional bias for a standard asymmetric split in diffusion processes. Originality/value – This paper for the first time compares in detail advantages and implementation complications of a front-fixing method when applied to the front-type field diffusion problems common to HTS. Particular attention is paid to accurate solutions in the region close to the moving front where rapid changes in material properties are responsible for large computational errors. Keywords - Modelling, Numerical analysis, Diffusion, High temperatures, Superconductors Paper type - Research pape

    An overview of the proper generalized decomposition with applications in computational rheology

    Get PDF
    We review the foundations and applications of the proper generalized decomposition (PGD), a powerful model reduction technique that computes a priori by means of successive enrichment a separated representation of the unknown field. The computational complexity of the PGD scales linearly with the dimension of the space wherein the model is defined, which is in marked contrast with the exponential scaling of standard grid-based methods. First introduced in the context of computational rheology by Ammar et al. [3] and [4], the PGD has since been further developed and applied in a variety of applications ranging from the solution of the Schrödinger equation of quantum mechanics to the analysis of laminate composites. In this paper, we illustrate the use of the PGD in four problem categories related to computational rheology: (i) the direct solution of the Fokker-Planck equation for complex fluids in configuration spaces of high dimension, (ii) the development of very efficient non-incremental algorithms for transient problems, (iii) the fully three-dimensional solution of problems defined in degenerate plate or shell-like domains often encountered in polymer processing or composites manufacturing, and finally (iv) the solution of multidimensional parametric models obtained by introducing various sources of problem variability as additional coordinates
    corecore