Modeling Frequency-Dependent Boundaries as Digital Impedance Filters in FDTD and K-DWM Room Acoustics Simulations

A new method for modeling-frequency-dependent boundaries in finite-difference time-domain (FDTD) and Kirchhoff variable digital waveguide mesh (K-DWM) room acoustics simulations is presented. The proposed approach allows the direct incorporation of a digital impedance filter (DIF) in the Multidimensional (2D or 3D) FDTD boundary model of a locally reacting surface. An explicit boundary update equation is obtained by carefully constructing a Suitable recursive formulation. The method is analyzed in terms of pressure wave reflectance for different wall impedance filters and angles of incidence. Results obtained from numerical experiments confirm the high accuracy of the proposed digital impedance filter boundary model, the reflectance of which matches locally reacting surface (LRS) theory closely. Furthermore a numerical boundary analysis (NBA) formula is provided as a technique for an analytic evaluation of the numerical reflectance of the proposed digital impedance filter boundary formulation.

Wideband and Isotropic Room Acoustics Simulation Using 2-D Interpolated FDTD Schemes

In this paper, a complete method for finite-difference time-domain modeling of rooms in 2-D using compact explicit schemes is presented. A family of interpolated schemes using a rectilinear, nonstaggered grid is reviewed, and the most accurate and isotropic schemes are identified. Frequency-dependent boundaries are modeled using a digital impedance filter formulation that is consistent with locally reacting surface theory. A structurally stable and efficient boundary formulation is constructed by carefully combining the boundary condition with the interpolated scheme. An analytic prediction formula for the effective numerical reflectance is given, and a stability proof provided. The results indicate that the identified accurate and isotropic schemes are also very accurate in terms of numerical boundary reflectance, and outperform directly related methods such as Yee's scheme and the standard digital waveguide mesh. In addition, one particular scheme-referred to here as the interpolated wideband scheme-is suggested as the best scheme for most applications.

A Phase Grating Approach to Modeling Surface Diffusion in FDTD Room Acoustics Simulations

In this paper, a method for modeling diffusive boundaries in finite difference time domain (FDTD) room acoustics simulations with the use of impedance filters is presented. The proposed technique is based on the concept of phase grating diffusers, and realized by designing boundary impedance filters from normal-incidence reflection filters with added delay. These added delays, that correspond to the diffuser well depths, are varied across the boundary surface, and implemented using Thiran allpass filters. The proposed method for simulating sound scattering is suitable for modeling high frequency diffusion caused by small variations in surface roughness and, more generally, diffusers characterized by narrow wells with infinitely thin separators. This concept is also applicable to other wave-based modeling techniques. The approach is validated by comparing numerical results for Schroeder diffusers to measured data. In addition, it is proposed that irregular surfaces are modeled by shaping them with Brownian noise, giving good control over the sound scattering properties of the simulated boundary through two parameters, namely the spectral density exponent and the maximum well depth.

Discretization and solution of the inhomogeneous Bogoliubov-de Gennes equations

In the presence of inhomogeneities, defects and currents, the equations describing a Bose-condensed ensemble of alkali atoms have to be solved numerically. By combining both linear and nonlinear equations within a Discrete Variable Representation framework, we describe a computational scheme for the solution of the coupled Bogoliubov-de Gennes (BdG) and nonlinear Schrodinger (NLS) equations for fields in a 3D spheroidal potential. We use the method to calculate the collective excitation spectrum and quasiparticle mode densities for excitations of a Bose condensed gas in a spheroidal trap. The method is compared against finite-difference and spectral methods, and we find the DVR computational scheme to be superior in accuracy and efficiency for the cases we consider. (C) 2004 Elsevier B.V. All rights reserved.

Finite Element Application to Analysis, Testing, and Manufacturing of Satellite Structures

Designing satellite structures poses an ongoing challenge as the interaction between analysis, experimental testing, and manufacturing phases is underdeveloped. Finite Element Analysis for Satellite Structures: Applications to Their Design, Manufacture and Testing explains the theoretical and practical knowledge needed to perform design of satellite structures. By layering detailed practical discussions with fully developed examples, Finite Element Analysis for Satellite Structures: Applications to Their Design, Manufacture and Testing provides the missing link between theory and implementation.
Computational examples cover all the major aspects of advanced analysis; including modal analysis, harmonic analysis, mechanical and thermal fatigue analysis using finite element method. Test cases are included to support explanations an a range of different manufacturing simulation techniques are described from riveting to shot peening to material cutting. Mechanical design of a satellites structures are covered in three steps: analysis step under design loads, experimental testing to verify design, and manufacturing.
Stress engineers, lecturers, researchers and students will find Finite Element Analysis for Satellite Structures: Applications to Their Design, Manufacture and Testing a key guide on with practical instruction on applying manufacturing simulations to improve their design and reduce project cost, how to prepare static and dynamic test specifications, and how to use finite element method to investigate in more details any component that may fail during testing.

Reliability Assessment of 3-D Space Frame Structures Applying Stochastic Finite Element Analysis

Throughout design development of satellite structure, stress engineer is usually challenged with randomness in applied loads and material properties. To overcome such problem, a risk-based design is applied which estimates satellite structure probability of failure under static and thermal loads. Determining probability of failure can help to update initially applied factors of safety that were used during structure preliminary design phase. These factors of safety are related to the satellite mission objective. Sensitivity-based analysis is to be implemented in the context of finite element analysis (probabilistic finite element method or stochastic finite element method (SFEM)) to determine the probability of failure for satellite structure or one of its components.

A Physically-Constrained Source Model For FDTD Acoustic Simulation

The Finite Difference Time Domain (FDTD) method is becoming increasingly popular for room acoustics simulation. Yet, the literature on grid excitation methods is relatively sparse, and source functions are traditionally implemented in a hard or additive form
using arbitrarily-shaped functions which do not necessarily obey the physical laws of sound generation. In this paper we formulate
a source function based on a small pulsating sphere model. A physically plausible method to inject a source signal into the grid
is derived from first principles, resulting in a source with a near-flat spectrum that does not scatter incoming waves. In the final
discrete-time formulation, the source signal is the result of passing a Gaussian pulse through a digital filter simulating the dynamics of the pulsating sphere, hence facilitating a physically correct means to design source functions that generate a prescribed sound field.

Tip-enhanced Raman microscopy: practicalities and limitations

The feasibility of apertureless scanning near-field Raman microscopy, exploiting the local enhancement in Raman scattering in the vicinity of a silver or gold tip, was investigated. Using the finite difference time domain method we calculated the enhancement of electric field strength, and hence Raman scattering, achieved through the resonant excitation of local modes in the tip. By modelling the frequency-dependent dielectric response of the metal tip we were able to highlight the resonant nature of the tip-enhancement and determine the excitation wavelength required for the strongest electric field enhancement, and hence Raman scattering intensity, which occurs for the excitation of modes localized at the tip apex. It is demonstrated that a peak Raman enhancement of 10(7)-fold should be achievable with <5 nm spatial resolution. We show that surface-enhanced Raman scattering from carbon contamination on a silver or gold tip can be significant. However, we find for a tip of radius of curvature 20 nm that the Raman enhancement should decay totally within 20 nm from the tip. Hence withdrawal of the tip by this distance should lead to the disappearance of the tip-enhanced signal, leaving only that from carbon contamination on the tip itself and the intrinsic signal from the sample. Copyright (C) 2003 John Wiley Sons, Ltd.

The development of a surface defect machining method for hard turning processes

Hard turning (HT) is a material removal process employing a combination of a single point cutting tool and high speeds to machine hard ferrous alloys which exhibit hardness values over 45 HRC. In this paper, a surface defect machining (SDM) method for HT is proposed which harnesses the combined advantages of porosity machining and pulsed laser pre-treatment processing. From previous experimental work, this was shown to provide better controllability of the process and improved quality of the machined surface. While the experiments showed promising results, a comprehensive understanding of this new technique could only be achieved through a rigorous, in depth theoretical analysis. Therefore, an assessment of the SDM technique was carried out using both finite element method (FEM) and molecular dynamics (MD) simulations.
FEM modelling was used to compare the conventional HT of AISI 4340 steel (52 HRC) using an Al2O3 insert with the proposed SDM method. The simulations showed very good agreement with the previously published experimental results. Compared to conventional HT, SDM provided favourable machining outcomes, such as reduced shear plane angle, reduced average cutting forces, improved surface roughness, lower residual stresses on the machined surface, reduced tool–chip interface contact length and increased chip flow velocity. Furthermore, a scientific explanation of the improved surface finish was revealed using a state-of-the-art MD simulation model which suggested that during SDM, a combination of both the cutting action and rough polishing action help improve the machined surface finish.

Energy conserving schemes for the simulation of musical instrument contact dynamics

Collisions are an innate part of the function of many musical instruments. Due to the nonlinear nature of contact forces, special care has to be taken in the construction of numerical schemes for simulation and sound synthesis. Finite difference schemes and other time-stepping algorithms used for musical instrument modelling purposes are normally arrived at by discretising a Newtonian description of the system. However because impact forces are non-analytic functions of the phase space variables, algorithm stability can rarely be established this way. This paper presents a systematic approach to deriving energy conserving schemes for frictionless impact modelling. The proposed numerical formulations follow from discretising Hamilton׳s equations of motion, generally leading to an implicit system of nonlinear equations that can be solved with Newton׳s method. The approach is first outlined for point mass collisions and then extended to distributed settings, such as vibrating strings and beams colliding with rigid obstacles. Stability and other relevant properties of the proposed approach are discussed and further demonstrated with simulation examples. The methodology is exemplified through a case study on tanpura string vibration, with the results confirming the main findings of previous studies on the role of the bridge in sound generation with this type of string instrument.

Impact of finite element idealisation on the prediction of welded fuselage stiffened panel buckling

Lap joints are widely used in the manufacture of stiffened panels and influence local panel sub-component stability, defining buckling unit dimensions and boundary conditions. Using the Finite Element method it is possible to model joints in great detail and predict panel buckling behaviour with accuracy. However, when modelling large panel structures such detailed analysis becomes computationally expensive. Moreover, the impact of local behaviour on global panel performance may reduce as the scale of the modelled structure increases. Thus this study presents coupled computational and experimental analysis, aimed at developing relationships between modelling fidelity and the size of the modelled structure, when the global static load to cause initial buckling is the required analysis output. Small, medium and large specimens representing welded lap-joined fuselage panel structure are examined. Two element types, shell and solid-shell, are employed to model each specimen, highlighting the impact of idealisation on the prediction of welded stiffened panel initial skin buckling.

Finite element modeling of debonding failures in FRP-strengthened RC beams: A dynamic approach

This paper is concerned with the finite element simulation of debonding failures in FRP-strengthened concrete beams. A key challenge for such simulations is that common solution techniques such as the Newton-Raphson method and the arc-length method often fail to converge. This paper examines the effectiveness of using a dynamic analysis approach in such FE simulations, in which debonding failure is treated as a dynamic problem and solved using an appropriate time integration method. Numerical results are presented to show that an appropriate dynamic approach effectively overcomes the convergence problem and provides accurate predictions of test results.

A Simplified Method for Ductility Calculation in RC Jacketed Columns

Reinforced concrete (RC) jacketing is a common method to retrofit existing columns with poor structural performance. It can be applied in two different ways: if the continuity of the jacket is ensured, the axial load of the column can be transferred to the jacket, which will be directly loaded; conversely, if no continuity is provided, the jacket induces only confinement action. In both cases the strength and ductility evaluation is rather complex, due to the different physical phenomena included, such as confinement, composite action core-jacket, preload, buckling of longitudinal bars.
Although different theoretical studies have been carried out to calculate the confinement effects, a practical approach to evaluate the flexural capacity and ductility is still missing. The calculation of these quantities is often related to the use of commercial computer programs, taking advantage of numerical methods such as fiber method or finite element method.
This paper presents a simplified approach to calculate the flexural strength and ductility of square RC jacketed sections subjected to axial load and bending moment. In particular the proposed approach is based on the calibration of the stress-block parameters including the confinement effect. Equilibrium equations are determined and buckling of longitudinal bars is modeled with a suitable stress-strain law. Moment-curvature curves are derived with simple calculations. Finally, comparisons are made with numerical analyses carried out with the code OpenSees and with experimental data available in the literature, showing good agreement.

Discrete-Time Conserved Quantities for Damped Oscillators

Numerical sound synthesis is often carried out using the finite difference time domain method. In order to analyse the stability of the derived models, energy methods can be used for both linear and nonlinear settings. For Hamiltonian systems the existence of a conserved numerical energy-like quantity can be used to guarantee the stability of the simulations. In this paper it is shown how to derive similar discrete conservation laws in cases where energy is dissipated due to friction or in the presence of an energy source due to an external force. A damped harmonic oscillator (for which an analytic solution is available) is used to present the proposed methodology. After showing how to arrive at a conserved quantity, the simulation of a nonlinear single reed shows an example of an application in the context of musical acoustics.

Manipulating plasmon modes in tip-enhanced Raman microscopy

The influence of polarization on the plasmon modes excited in tip-enhanced near-field optical microscopy have been investigated using the Finite Difference Time Domain Method. Analysis of the calculated results have laid particular emphasis on the ability to align local field enhancements with the orientation of molecules in order to optimize Raman signals, with particular relevance to recent experimental work on carbon nanotubes.