We present a Heterogeneous Multiscale Method for the Landau-Lifshitz equation with a highly oscillatory diffusion coefficient, a simple model for a ferromagnetic composite. A finite element macro scheme is combined with a finite difference micro model to approximate the effective equation corresponding to the original problem. This makes it possible to obtain effective solutions to problems with rapid material variations on a small scale, described by ε≪1, which would be too expensive to resolve in a conventional simulation.

We present a framework for adaptive finite element computation of turbulent flow and fluid structure interaction, with focus on general algorithms that allow for complex geometry and deforming domains. We give basic models and finite element discretization methods, adaptive algorithms and strategies for efficient parallel implementation. To illustrate the capabilities of the computational framework, we show a number of application examples from aerodynamics, aero-acoustics, biomedicine and geophysics. The computational tools are free to download open source as Unicorn, and as a high performance branch of the finite element problem solving environment DOLFIN, both part of the FEniCS project.

In this paper we present a computational study of the stability of time dependent dual problems for compressible flow at high Reynolds numbers in two dimensions. The dual problem measures the sensitivity of an output functional with respect to numerical errors and is a key part of goal oriented a posteriori error estimation. Our investigation shows that the dual problem associated with the computation of the drag force for the compressible Euler/Navier-Stokes equations, which are approximated numerically using different temporal discretization and stabilization techniques, is unstable and exhibits blow-up for several Mach regimes considered in this paper.

This chapter provides a description of the technology of Unicorn focusing on simple, efficient and general algorithms and software for the Unified Continuum (UC) concept and the adaptive General Galerkin (G2) discretization as a unified approach to continuum mechanics. We describe how Unicorn fits into the FEniCS framework, how it interfaces to other FEniCS components, what interfaces and functionality Unicorn provides itself and how the implementation is designed. We also present some examples in fluid–structure interaction and adaptivity computed with Unicorn. 

The massive computational cost for resolving all turbulent scales makes a direct numerical simulation of the underlying Navier-Stokes equations impossible in most engineering applications. We present recent advances in parallel adaptive finite element methodology that enable us to efficiently compute time resolved approximations for complex geometries with error control. In this paper we present a LES simulation of turbulent flow past a full car model, where we adaptively refine the unstructured mesh to minimize the error in drag prediction. The simulation was partly carried out on the new Cray XE6 at PDC/KTH where the solver shows near optimal strong and weak scaling for the entire adaptive process.