11,713 research outputs found
Efficient solution of 3D electromagnetic eddy-current problems within the finite volume framework of OpenFOAM
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
Verification of Unstructured Grid Adaptation Components
Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. Laminar flow over a delta wing and turbulent flow over an ONERA M6 wing are verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and a moment. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. The wing adaptation cases in the current study document a clear improvement to existing grid adaptation procedures. The stage is set for the infusion of verified grid adaptation into production fluid flow simulations
A non-hybrid method for the PDF equations of turbulent flows on unstructured grids
In probability density function (PDF) methods of turbulent flows, the joint
PDF of several flow variables is computed by numerically integrating a system
of stochastic differential equations for Lagrangian particles. A set of
parallel algorithms is proposed to provide an efficient solution of the PDF
transport equation, modeling the joint PDF of turbulent velocity, frequency and
concentration of a passive scalar in geometrically complex configurations. An
unstructured Eulerian grid is employed to extract Eulerian statistics, to solve
for quantities represented at fixed locations of the domain (e.g. the mean
pressure) and to track particles. All three aspects regarding the grid make use
of the finite element method (FEM) employing the simplest linear FEM shape
functions. To model the small-scale mixing of the transported scalar, the
interaction by exchange with the conditional mean model is adopted. An adaptive
algorithm that computes the velocity-conditioned scalar mean is proposed that
homogenizes the statistical error over the sample space with no assumption on
the shape of the underlying velocity PDF. Compared to other hybrid
particle-in-cell approaches for the PDF equations, the current methodology is
consistent without the need for consistency conditions. The algorithm is tested
by computing the dispersion of passive scalars released from concentrated
sources in two different turbulent flows: the fully developed turbulent channel
flow and a street canyon (or cavity) flow. Algorithmic details on estimating
conditional and unconditional statistics, particle tracking and particle-number
control are presented in detail. Relevant aspects of performance and
parallelism on cache-based shared memory machines are discussed.Comment: Accepted in Journal of Computational Physics, Feb. 20, 200
- …