1. |
- Berggren, Martin, et al.
(författare)
-
Weak material approximation of holes with traction-free boundaries
- 2012
-
Ingår i: SIAM Journal on Numerical Analysis. - : SIAM Publications Online. - 0036-1429 .- 1095-7170. ; 50:4, s. 1827-1848
-
Tidskriftsartikel (refereegranskat)abstract
- Consider the solution of a boundary-value problem for steady linear elasticity in which the computational domain contains one or several holes with traction-free boundaries. The presence of holes in the material can be approximated using a weak material; that is, the relative density of material rho is set to 0 < epsilon = rho << 1 in the hole region. The weak material approach is a standard technique in the so-called material distribution approach to topology optimization, in which the inhomogeneous relative density of material is designated as the design variable in order to optimize the spatial distribution of material. The use of a weak material ensures that the elasticity problem is uniquely solvable for each admissible value rho is an element of [epsilon, 1] of the design variable. A finite-element approximation of the boundary-value problem in which the weak material approximation is used in the hole regions can be viewed as a nonconforming but convergent approximation of a version of the original problem in which the solution is continuously and elastically extended into the holes. The error in this approximation can be bounded by two terms that depend on epsilon. One term scales linearly with epsilon with a constant that is independent of the mesh size parameter h but that depends on the surface traction required to fit elastic material in the deformed holes. The other term scales like epsilon(1/2) times the finite-element approximation error inside the hole. The condition number of the weak material stiffness matrix scales like epsilon(-1), but the use of a suitable left preconditioner yields a matrix with a condition number that is bounded independently of epsilon. Moreover, the preconditioned matrix admits the limit value epsilon -> 0, and the solution of corresponding system of equations yields in the limit a finite-element approximation of the continuously and elastically extended problem.
|
|
2. |
- Fotios, Kasolis, et al.
(författare)
-
Analysis of fictitious domain approximations of hard scatterers
- 2015
-
Ingår i: SIAM Journal on Numerical Analysis. - Philadelphia : Siam publications. - 0036-1429 .- 1095-7170. ; 53:5, s. 2347-2362
-
Tidskriftsartikel (refereegranskat)abstract
- Consider the Helmholtz equation del center dot alpha del p+k(2 alpha)p = 0 in a domain that contains a so-called hard scatterer. The scatterer is represented by the value alpha = epsilon, for 0 < epsilon << 1, whereas alpha = 1 whenever the scatterer is absent. This scatterer model is often used for the purpose of design optimization and constitutes a fictitious domain approximation of a body characterized by homogeneous Neumann conditions on its boundary. However, such an approximation results in spurious resonances inside the scatterer at certain frequencies and causes, after discretization, ill-conditioned system matrices. Here, we present a stabilization strategy that removes these resonances. Furthermore, we prove that, in the limit epsilon -> 0, the stabilized problem provides linearly convergent approximations of the solution to the problem with an exactly modeled scatterer. Numerical experiments indicate that a finite element approximation of the stabilized problem is free from internal resonances, and they also suggest that the convergence rate is indeed linear with respect to epsilon.
|
|
3. |
|
|
4. |
- Kasolis, Fotios, et al.
(författare)
-
Fixed-mesh curvature-parameterized shape optimization of an acoustic horn
- 2012
-
Ingår i: Structural and multidisciplinary optimization (Print). - : Springer. - 1615-147X .- 1615-1488. ; 46:5, s. 727-738
-
Tidskriftsartikel (refereegranskat)abstract
- We suggest a boundary shape optimization approach in which the optimization is carried out on the coefficients in a boundary parameterization based on a local, discrete curvature. A fixed mesh is used to numerically solve the governing equations, in which the geometry is represented through inhomogeneous coefficients, similarly as done in the material distribution approach to topology optimization. The method is applied to the optimization of an acoustic horn in two space dimensions. Numerical experiments show that this method can calculate the horn's transmission properties as accurately as a traditional, body-fitted approach. Moreover, the use of a fixed mesh allows the optimization to create shapes that would be difficult to handle with a traditional approach that uses deformations of a body-fitted mesh. The parameterization inherently promotes smooth designs without unduly restriction of the design flexibility. The optimized, smooth horns consistently show favorable transmission properties.
|
|
5. |
|
|
6. |
- Kasolis, Fotios, 1983-
(författare)
-
The material distribution method : analysis and acoustics applications
- 2014
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- For the purpose of numerically simulating continuum mechanical structures, different types of material may be represented by the extreme values {,1}, where 0<1, of a varying coefficient in the governing equations. The paramter is not allowed to vanish in order for the equations to be solvable, which means that the exact conditions are approximated. For example, for linear elasticity problems, presence of material is represented by the value = 1, while = provides an approximation of void, meaning that material-free regions are approximated with a weak material. For acoustics applications, the value = 1 corresponds to air and = to an approximation of sound-hard material using a dense fluid. Here we analyze the convergence properties of such material approximations as !0, and we employ this type of approximations to perform design optimization.In Paper I, we carry out boundary shape optimization of an acoustic horn. We suggest a shape parameterization based on a local, discrete curvature combined with a fixed mesh that does not conform to the generated shapes. The values of the coefficient , which enters in the governing equation, are obtained by projecting the generated shapes onto the underlying computational mesh. The optimized horns are smooth and exhibit good transmission properties. Due to the choice of parameterization, the smoothness of the designs is achieved without imposing severe restrictions on the design variables.In Paper II, we analyze the convergence properties of a linear elasticity problem in which void is approximated by a weak material. We show that the error introduced by the weak material approximation, after a finite element discretization, is bounded by terms that scale as and 1/2hs, where h is the mesh size and s depends on the order of the finite element basis functions. In addition, we show that the condition number of the system matrix scales inversely proportional to , and we also construct a left preconditioner that yields a system matrix with a condition number independent of .In Paper III, we observe that the standard sound-hard material approximation with = gives rise to ill-conditioned system matrices at certain wavenumbers due to resonances within the approximated sound-hard material. To cure this defect, we propose a stabilization scheme that makes the condition number of the system matrix independent of the wavenumber. In addition, we demonstrate that the stabilized formulation performs well in the context of design optimization of an acoustic waveguide transmission device.In Paper IV, we analyze the convergence properties of a wave propagation problem in which sound-hard material is approximated by a dense fluid. To avoid the occurrence of internal resonances, we generalize the stabilization scheme presented in Paper III. We show that the error between the solution obtained using the stabilized soundhard material approximation and the solution to the problem with exactly modeled sound-hard material is bounded proportionally to .
|
|