Search

Incremental models for magnetic vector hysteresis have been developed in previous works in accordance with basic principles of thermodynamics. In this paper, we present an equivalent representation of the associated hysteresis operator in terms of a co-energy functional which is useful for magnetic field computations based on a scalar potential. Using convex duality, we further define the corresponding energy functional and the associated inverse hysteresis operator which is required for computations based on the vector potential. The equivalence of the two representations with the energy-based hysteresis models proposed in earlier works is demonstrated and numerical results for some typical test problems are presented obtained by finite element simulation of corresponding scalar and vector potential formulations.

Fixed-point or Newton-methods are typically employed for the numerical solution of nonlinear systems arising from discretization of nonlinear magnetic field problems. We here discuss an alternative strategy which uses local Quasi-Newton updates to construct appropriate linearizations of the material behavior during the nonlinear iteration. The resulting scheme shows similar fast convergence as the Newton-method but, like the fixed-point methods, does not require derivative information of the underlying material law. As a consequence, the method can be used for the efficient solution of models with hysteresis which involve nonsmooth material behavior. The implementation of the proposed scheme can be realized in standard finite-element codes in parallel to the fixed-point and the Newton method. A full convergence analysis of all three methods is established proving global mesh-independent convergence. The theoretical results and the performance of the nonlinear iterative schemes are evaluated by computational tests for a typical benchmark problem.

Porous acoustic absorbers have excellent properties in the low-frequency range when positioned in room edges, therefore they are a common method for reducing low-frequency reverberation. However, standard room acoustic simulation methods such as ray tracing and mirror sources are invalid for low frequencies in general which is a consequence of using geometrical methods, yielding a lack of simulation tools for these so-called edge absorbers. In this article, a validated finite element simulation model is presented, which is able to predict the effect of an edge absorber on the acoustic field. With this model, the interaction mechanisms between room and absorber can be studied by high-resolved acoustic field visualizations in both room and absorber. The finite element model is validated against transfer function data computed from impulse response measurements in a reverberation chamber in style of ISO 354. The absorber made of Basotect is modeled using the Johnson-Champoux-Allard-Lafarge model, which is fitted to impedance tube measurements using the four-microphone transfer matrix method. It is shown that the finite element simulation model is able to predict the influence of different edge absorber configurations on the measured transfer functions to a high degree of accuracy. The evaluated third-octave band error exhibits deviations of 3.25dB to 4.11dB computed from third-octave band averaged spectra., Comment: 20 pages, 16 figures, 3 tables, Preprint submitted to Acta Acustica

This article deals with large-eddy simulations of 3D incompressible laryngeal flow followed by acoustic simulations of human phonation of five cardinal english vowels /u, i, \textipa{A}, o, {\ae}/. The flow and aeroacoustic simulations were performed in OpenFOAM and in-house code openCFS, respectively. Given the large variety of scales in the flow and acoustics, the simulation is separated into two steps: (1) computing the flow in the larynx using the finite volume method on a fine 2.2M grid followed by (2) computing the sound sources separately and wave propagation to the radiation zone around the mouth using the finite element method on a coarse 33k acoustic grid. The numerical results showed that the anisotropic minimum dissipation model, which is not well known since it is not available in common CFD software, predicted stronger sound pressure levels at higher harmonics and especially at first two formants than the wall-adapting local eddy-viscosity model. We implemented the model as a new open library in OpenFOAM and deployed the model on turbulent flow in the larynx with positive impact on the quality of simulated vowels. Numerical simulations are in very good agreement with positions of formants from measurements., Comment: 15 pages, 20 figures, 20 audio files in att

This work compares two Nitsche-type approaches to treat non-conforming triangulations for a high-order discontinuous Galerkin (DG) solver for the acoustic conservation equations. The first approach (point-to-point interpolation) uses inexact integration with quadrature points prescribed by a primary element. The second approach uses exact integration (mortaring) by choosing quadratures depending on the intersection between non-conforming elements. In literature, some excellent properties regarding performance and ease of implementation are reported for point-to-point interpolation. However, we show that this approach can not safely be used for DG discretizations of the acoustic conservation equations since, in our setting, it yields spurious oscillations that lead to instabilities. This work presents a test case in that we can observe the instabilities and shows that exact integration is required to maintain a stable method. Additionally, we provide a detailed analysis of the method with exact integration. We show optimal spatial convergence rates globally and in each mesh region separately. The method is constructed such that it can natively treat overlaps between elements. Finally, we highlight the benefits of non-conforming discretizations in acoustic computations by a numerical test case with different fluids.

In low Mach number aeroacoustics, the well known disparity of scales makes it possible to apply efficient hybrid simulation models using different meshes for flow and acoustics, which leads to a powerful computational procedure. Our study applies the hybrid workflow to the computationally efficient perturbed convective wave equation with only one scalar unknown, the acoustic velocity potential. The workflow of this aeroacoustic approach is based on three steps: 1. perform unsteady incompressible flow computations on a sub-domain; 2. compute the acoustic sources; 3. simulate the acoustic field using a mesh specifically suited for computational aeroacoustics. In general, hybrid aeroacoustic methods seek for robust and conservative mesh-to-mesh transformation of the aeroacoustic sources while high computational efficiency is ensured. In this paper, we investigate the accuracy of a cell-centroid based conservative interpolation scheme compared to the more accurate cut-volume cell approach and their application to the computation of rotating systems, namely an axial fan. The capability and robustness of the cut-volume cell interpolation in a hybrid workflow on different meshes are investigated by a grid convergence study. The results of the acoustic simulation are in good agreement with measurements thus demonstrating the applicability of the conservative cut-volume cell interpolation to rotating systems.

In modern transport systems, passengers' comfort is greatly influenced by flow-induced noise. In this study we investigate a generic deep cavity with an overhanging lip, mimicking a door gap in a vehicle, that is overflowed by air at two different free stream velocities, 26.8m/s and 50m/s. The turbulent boundary layer and the acoustic waves interact with the cavity's geometry and form a strong feedback mechanism. In the present work, we focus on the details of the compressible turbulent flow structures and their variations concerning the velocity, the boundary layer as well as the domain dimensionality for a later acoustic simulation within a hybrid aeroacoustic workflow. Furthermore, we verify the feasibility of reducing the acoustic computational domain from 3D to 2D for this application by conducting a coherence study of acoustically active flow structures in the spanwise direction. The role of the three-dimensional Taylor-G\"ortler vortices from the recirculation regarding the vortex formation and the vortex-edge interaction was also evaluated. Remarkably for the lower approaching velocity (26.8m/s), we found a special vortex-edge interaction, namely an alternating sequence of complete clipping and a subsequent partial escape. Lastly, we assigned previous unknown peaks in the pressure spectrum to their corresponding mechanisms., Comment: 16 pages, 12 Figures

Hybrid aeroacoustic methods seek for computational efficiency and robust noise prediction. Using already existing aeroacoustic wave equations, we propose a general hybrid aeroacoustic method, based on compressible source data. The main differences to current state of the art aeroacoustic analogies are that an additional decomposition is used to compute the aeroacoustic source terms and their application is extended to the source formulation, based on compressible flow data. By applying the Helmholtz-Hodge decomposition on arbitrary domains, we extract the incompressible projection (non-radiating base flow) of a compressible flow simulation. This method maintains the favorable properties of the hybrid aeroacoustic method while incorporating compressible effects on the base flow. The capabilities are illustrated for the aeroacoustic benchmark case, "cavity with a lip", involving acoustic feedback. The investigation is based on the equation of vortex sound, to incorporate convective effects during wave propagation., Comment: 12 pages, 11 figures