Boundary Element and Finite Element Coupling for Aeroacoustics Simulations

We consider the scattering of acoustic perturbations in a presence of a flow. We suppose that the space can be split into a zone where the flow is uniform and a zone where the flow is potential. In the first zone, we apply a Prandtl-Glauert transformation to recover the Helmholtz equation. The well-known setting of boundary element method for the Helmholtz equation is available. In the second zone, the flow quantities are space dependent, we have to consider a local resolution, namely the finite element method. Herein, we carry out the coupling of these two methods and present various applications and validation test cases. The source term is given through the decomposition of an incident acoustic field on a section of the computational domain's boundary.Comment: 25 page

Equivalent drawbead model in finite element simulations

In 3D simulations of the deep drawing process the drawbead geometries are seldom included. Therefore equivalent drawbeads are used. In order to investigate the drawbead behaviour a 2D plane strain finite element model was used. For verification of this model experiments were performed. The analyses showed that not only the restraining force should be applied but also the strain changes. The effects of the restraining force and the strain change were implemented in an equivalent drawbead. The effect of using the equivalent drawbead is demonstrated with a few example

Modelling drawbeads with finite elements and verification

Drawbeads are commonly used in deep drawing processes to control the flow of the blank during the forming operation. In finite element simulations of deep drawing the drawbead geometries are seldom included because of the small radii; because of these small radii a very large number of elements is required in 3-D simulations. To cope with this problem, a 2-D analysis of the drawbead has been performed and the calculated restraining force will be applied in the near future in 3-D simulations with an equivalent drawbead element. Modelling drawbeads by only applying an additional restraining force is not satisfactory. During the flow of the material through a drawbead, the strain distribution changes and the material usually becomes thinner. These effects must be incorporated in the equivalent drawbead element.\ud \ud For the modelling of the drawbead a 2-D plane strain finite element model was developed. Several simulations were carried out to investigate the behaviour of the drawbead. Various geometries were investigated, the friction coefficient was varied and also the frictionless case was taken into account.\ud \ud To verify the model an experimental set-up was built. An extensive set of drawbead geometries was used. The results are compared with the finite element calculations and the similarity is very satisfactory

Crystal plasticity finite element simulations of cast α-uranium

α-uranium, the stable phase of uranium up to 670◦C, has a base-centred orthorombic crystal structure. This crystal structure gives rise to elastic and thermal anisotropy, meaning α-uranium exhibits complex deformation and fracture behaviour. Understanding the relationship between the microstructure and mechanical properties is important to prevent fracture during manufacture and usage of components. The lattice of α-uranium corresponds to a distorted close-packed-hexagonal crystal structure and it exhibits twins of both the 1st and 2nd kind. Therefore, detailed examination of the behaviour of α-uranium can also contribute to the general understanding of the interaction between plasticity, twinning and fracture in hcp crystals. Plastic deformation in α-uranium can be accommodated by 4 slip systems and 3 twin systems, previously identified by McCabe et al. These deformation modes are implemented into a crystal plasticity finite element (CPFE) material model. A temperature dependent, dislocation density based law is implemented to describe the critical resolved shear stress on the different slip/twin systems. The strong anisotropic thermal expansion behaviour is taken into account to simulate the development of internal residual stresses following casting of the material. During cooling, the internal stresses in α-uranium are sufficient to induce plasticity. This effect is quantified using polycrystal simulations, in which first the temperature is decreased, then plastic relaxation takes place, followed by application of a mechanical load. The asymmetry between mechanical properties in tension and compression, due to the presence of twins, is investigated. The model is calibrated using stress strain curves and the lattice strain found from published neutron diffraction experiments carried out on textured samples at ISIS. The strength of the slip systems is found to be lower than in fine grained material, while the strength of the twin system is similar to single crystals. The CPFE method allows the heterogeneity of the strain between neighbouring grains and its influence on the evolution of the internal stress state to be investigated

Time Domain Simulations of EMRIs using Finite Element Methods

This is a brief report on time-domain numerical simulations of extreme-mass-ratio binaries based on finite element methods. We discuss a new technique for solving the perturbative equations describing a point-like object orbiting a non-rotating massive black hole and the prospects of using it for the evaluation of the gravitational self-force responsible of the inspiral of these binary systems. We also discuss the perspectives of transferring this technology to the more astrophysically relevant case of a central rotating massive black hole.Comment: 5 pages. Submitted to the proceedings of the 6th LISA symposiu

Time-resolved investigation of magnetization dynamics of arrays of non-ellipsoidal nanomagnets with a non-uniform ground state

We have performed time-resolved scanning Kerr microscopy (TRSKM) measurements upon arrays of square ferromagnetic nano-elements of different size and for a range of bias fields. The experimental results were compared to micromagnetic simulations of model arrays in order to understand the non-uniform precessional dynamics within the elements. In the experimental spectra two branches of excited modes were observed to co-exist above a particular bias field. Below the so-called crossover field, the higher frequency branch was observed to vanish. Micromagnetic simulations and Fourier imaging revealed that modes from the higher frequency branch had large amplitude at the center of the element where the effective field was parallel to the bias field and the static magnetization. Modes from the lower frequency branch had large amplitude near the edges of the element perpendicular to the bias field. The simulations revealed significant canting of the static magnetization and the effective field away from the direction of the bias field in the edge regions. For the smallest element sizes and/or at low bias field values the effective field was found to become anti-parallel to the static magnetization. The simulations revealed that the majority of the modes were de-localized with finite amplitude throughout the element, while the spatial character of a mode was found to be correlated with the spatial variation of the total effective field and the static magnetization state. The simulations also revealed that the frequencies of the edge modes are strongly affected by the spatial distribution of the static magnetization state both within an element and within its nearest neighbors

Poor-man's model of hollow-core anti-resonant fibers

We investigate various methods for extending the simple analytical capillary model to describe the dispersion and loss of anti-resonant hollow-core fibers without the need of detailed finite-element simulations across the desired wavelength range. This poor-man's model can with a single fitting parameter quite accurately mimic dispersion and loss resonances and anti-resonances from full finite-element simulations. Due to the analytical basis of the model it is easy to explore variations in core size and cladding wall thickness, and should therefore provide a valuable tool for numerical simulations of the ultrafast nonlinear dynamics of gas-filled hollow-core fibers.Comment: In preparatio

Finite-Element Simulations of Light Propagation through Circular Subwavelength Apertures

Light transmission through circular subwavelength apertures in metallic films with surrounding nanostructures is investigated numerically. Numerical results are obtained with a frequency-domain finite-element method. Convergence of the obtained observables to very low levels of numerical error is demonstrated. Very good agreement to experimental results from the literature is reached, and the utility of the method is demonstrated in the investigation of the influence of geometrical parameters on enhanced transmission through the apertures

A finite element method for the resolution of the Reduced Navier-Stokes/Prandtl equations

A finite element method to solve the bidimensional Reduced Navier-Stokes Prandtl (RNS/P) equations is described. These equations are an asymptotical simplification of the full Navier-Stokes equations, obtained when one dimension of the domain is of one order smaller than the others. These aretherefore of particular interest to describe flows in channels or pipes of small diameter. A low order finite element discretization, based on a piecewise constant approximation of the pressure, is proposed and analyzed. Numerical experiments which consist in fluid flow simulations within a constricted pipe are provided. Comparisons with Navier-Stokes simulations allow to evaluate the performance of prediction of the finite element method, and of the model itself