Abstract:We discuss how Dirichlet boundary conditions can be directly imposed for the Moulinec–Suquet discretization on the boundary of rectangular domains in iterative schemes based on the fast Fourier transform (FFT) and computational homogenization problems in mechanics. Classically, computational homogenization methods based on the fast Fourier transform work with periodic boundary conditions. There are applications, however, when Dirichlet (or Neumann) boundary conditions are required. For thermal homogenization problems, it is straightforward to impose such boundary conditions by using discrete sine (and cosine) transforms instead of the FFT. This approach, however, is not readily extended to mechanical problems due to the appearance of mixed derivatives in the Lamé operator of elasticity. Thus, Dirichlet boundary conditions are typically imposed either by using Lagrange multipliers or a "buffer zone" with a high stiffness. Both strategies lead to formulations which do not share the computational advantages of the original FFT-based schemes. The work at hand introduces a technique for imposing Dirichlet boundary conditions directly without the need for indefinite systems. We use a formulation on the deformation gradient—also at small strains—and employ the Green's operator associated to the vector Laplacian. Then, we develop the Moulinec–Suquet discretization for Dirichlet boundary conditions—requiring carefully selected weights at boundary points—and discuss the seamless integration into existing FFT-based computational homogenization codes based on dedicated discrete sine/cosine transforms. The article culminates with a series of well-chosen numerical examples demonstrating the capabilities of the introduced technology.

A fast numerical method for the conductivity of heterogeneous media with Dirichlet boundary conditions based on discrete sine–cosine transforms

A Discrete Sine‐Cosine Transforms Galerkin Method for the Conductivity of Heterogeneous Materials With Mixed Dirichlet/Neumann Boundary Conditions

Discontinuous Galerkin Method Using Laguerre Polynomials for Solving a Time-Domain Electric Field Integral Equation

A Fast Time Domain PMCHW Equations for Predicting Transient Responses of Dielectric Objects

FFT-based simulations of heterogeneous conducting materials with combined non-uniform Neumann, periodic and Dirichlet boundary conditions

Finite Analytic Numerical Method for Three‐dimensional Quasi‐laplace Equation with Conductivity in Tensor Form

Fast integral equation methods for linear and semilinear heat equations in moving domains

Rapidly Convergent Expansion Method For Calculating The Effective Conductivity Of Three-Dimensional Lattices Of Symmetric Inclusions

Solution of 1D Poisson Equation with Neumann-Dirichlet and Dirichlet-Neumann Boundary Conditions, Using the Finite Difference Method

General Boundary Element Method for Poisson Equation with Spatially Varying Conductivity

Finite element approximation of data-driven problems in conductivity

Imposing different boundary conditions for thermal computational homogenization problems with FFT‐ and tensor‐train‐based Green's operator methods

Numerical simulation for non-linear thermal wave

A Fast Integral Equation Formulation for Superconducting Structures

Imposing Dirichlet boundary conditions directly for FFT-based computational micromechanics

Fourier-based schemes with modified Green operator for computing the electrical response of heterogeneous media with accurate local fields

A Fast Fourier-Galerkin Method Solving Boundary Integral Equations for the Helmholtz Equation with Exponential Convergence

A Multiresolution Method for Distributed Conductivity Estimation of Maxwell Equation

Numerical validation of an explicit P1 finite-element scheme for Maxwell's equations in a polygon with variable permittivity away from its boundary

Discontinuous Galerkin Method with a Novel Physics-Informed Flux for Elastic Wave Simulations in Heterogeneous Media