Finite difference time domain modeling of phase grating diffusion

In this paper, a method for modeling diffusion caused by non-smooth boundary surfaces in simulations of room acoustics using finite difference time domain (FDTD) technique is investigated. The proposed approach adopts the well-known theory of phase grating diffusers to efficiently model sound scattering from rough surfaces. The variation of diffuser well-depths is attained by nesting allpass filters within the reflection filters from which the digital impedance filters used in the boundary implementation are obtained. The presented technique is appropriate for modeling diffusion at high frequencies caused by small surface roughness and generally diffusers that have narrow wells and infinitely thin separators. The diffusion coefficient was measured with numerical experiments for a range of fractional Brownian diffusers.

Modal Representation of the Resonant Body within a Finite Difference Framework for Simulation of String Instruments

This paper investigates numerical simulation of a string coupled
transversely to a resonant body. Starting from a complete nite
difference formulation, a second model is derived in which the
body is represented in modal form. The main advantage of this hybrid form is that the body model is scalable, i.e. the computational
complexity can be adjusted to the available processing power. Numerical results are calculated and discussed for simplied models
in the form of string-string coupling and string-plate coupling.

How can cells sense the elasticity of a substrate? An analysis using a cell tensegrity model

A eukaryotic cell attaches and spreads on substrates, whether it is the extracellular matrix naturally produced by the cell itself, or artificial materials, such as tissue-engineered scaffolds. Attachment and spreading require the cell to apply forces in the nN range to the substrate via adhesion sites, and these forces are balanced by the elastic response of the substrate. This mechanical interaction is one determinant of cell morphology and, ultimately, cell phenotype. In this paper we use a finite element model of a cell, with a tensegrity structure to model the cytoskeleton of actin filaments and microtubules, to explore the way cells sense the stiffness of the substrate and thereby adapt to it. To support the computational results, an analytical 1D model is developed for comparison. We find that (i) the tensegrity hypothesis of the cytoskeleton is sufficient to explain the matrix-elasticity sensing, (ii) cell sensitivity is not constant but has a bell-shaped distribution over the physiological matrix-elasticity range, and (iii) the position of the sensitivity peak over the matrix-elasticity range depends on the cytoskeletal structure and in particular on the F-actin organisation. Our model suggests that F-actin reorganisation observed in mesenchymal stem cells (MSCs) in response to change of matrix elasticity is a structural-remodelling process that shifts the sensitivity peak towards the new value of matrix elasticity. This finding discloses a potential regulatory role of scaffold stiffness for cell differentiation.

Cold-formed steel sections with web openings subjected to web crippling under two-flange loading conditions - part I: Tests and finite element analysis

The results of 82 web crippling tests are presented, with 20 tests conducted on channel sections without web openings and 62 tests conducted on channel sections with web openings. The tests consider both end-two-flange and interior-two-flange loading conditions. In the case of the tests with web openings, the hole was located directly under the concentrated load. The concentrated load was applied through bearing plates; the effect of different bearing lengths is investigated. In addition, the cases of both flanges fastened and unfastened to the support is considered. A non-linear elasto-plastic finite element model is described, and the results compared against the laboratory test results; a good agreement was obtained in terms of both strength and failure modes.

Optimal diet choice for large herbivores: an extended contingency model

1. A more general contingency model of optimal diet choice is developed, allowing for simultaneous searching and handling, which extends the theory to include grazing and browsing by large herbivores.</p><p>2. Foraging resolves into three modes: purely encounter-limited, purely handling-limited and mixed-process, in which either a handling-limited prey type is added to an encounter-limited diet, or the diet becomes handling-limited as it expands.</p><p>3. The purely encounter-limited diet is, in general, broader than that predicted by the conventional contingency model,</p><p>4. As the degree of simultaneity of searching and handling increases, the optimal diet expands to the point where it is handling-limited, at which point all inferior prey types are rejected,</p><p>5. Inclusion of a less profitable prey species is not necessarily independent of its encounter rate and the zero-one rule does not necessarily hold: some of the less profitable prey may be included in the optimal diet. This gives an optimal foraging explanation for herbivores' mixed diets.</p><p>6. Rules are shown for calculating the boundary between encounter-limited and handling-limited diets and for predicting the proportion of inferior prey to be included in a two-species diet,</p><p>7. The digestive rate model is modified to include simultaneous searching and handling, showing that the more they overlap, the more the predicted diet-breadth is likely to be reduced.</p>

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.

A progressive failure model for composite laminates subjected to low velocity impact damage

This paper presents a 3-D failure model for predicting the dynamic material response of composite laminates under impact loading. The formulation is based on the Continuum Damage Mechanics (CDM) approach and enables the control of the energy dissipation associated with each failure mode regardless of mesh refinement and fracture plane orientation. Internal thermodynamically irreversible damage variables were defined in order to quantify damage concentration associated with each possible failure mode and predict the gradual stiffness reduction during the impact damage process. The material model has been implemented into LS-DYNA explicit finite element code within solid elements and it has proven to be capable of reproducing experimental results with good accuracy in terms of static/dynamic responses, absorbed energy and extent of damage.

A simple non-linear global-local finite element methodology for composite structures

A simple non-linear global-local finite element methodology is presented. A global coarse model, using 2-D shell elements, is solved non-linearly and the displacements and rotations around a region of interest are applied, as displacement boundary conditions, to a refined local 3-D model using Kirchhoff plate assumptions. The global elements' shape functions are used to interpolate between nodes. The local model is then solved non-linearly with an incremental scheme independent of that used for the global model.

Global quantum correlations in finite-size spin chains

We perform an extensive study of the properties of global quantum correlations in finite-size one-dimensional quantum spin models at finite temperature. By adopting a recently proposed measure for global quantum correlations (Rulli and Sarandy 2011 Phys. Rev. A 84 042109), called global discord, we show that critical points can be neatly detected even for many-body systems that are not in their ground state. We consider the transverse Ising model, the cluster-Ising model where three-body couplings compete with an Ising-like interaction, and the nearest-neighbor XX Hamiltonian in transverse magnetic field. These models embody our canonical examples showing the sensitivity of global quantum discord close to criticality. For the Ising model, we find a universal scaling of global discord with the critical exponents pertaining to the Ising universality class.

Finite Element Analysis of Drilling of Carbon Fibre Reinforced Composites

Despite the increased applications of the composite materials in aerospace due to their exceptional physical and mechanical properties, the machining of composites remains a challenge. Fibre reinforced laminated composites are prone to different damages during machining process such as delamination, fibre pull-out, microcracks, thermal damages. Optimization of the drilling process parameters can reduces the probability of these damages. In the current research, a 3D finite element (FE) model is developed of the process of drilling in the carbon fibre reinforced composite (CFC). The FE model is used to investigate the effects of cutting speed and feed rate on thrust force, torque and delamination in the drilling of carbon fiber reinforced laminated composite. A mesoscale FE model taking into account of the different oriented plies and interfaces has been proposed to predict different damage modes in the plies and delamination. For validation purposes, experimental drilling tests have been performed and compared to the results of the finite element analysis. Using Matlab a digital image analysis code has been developed to assess the delamination factor produced in CFC as a result of drilling. © Springer Science+Business Media B.V. 2011.

Physical and numerical constraints in source modeling for finite difference simulation of room acoustics

In finite difference time domain simulation of room acoustics, source functions are subject to various constraints. These depend on the way sources are injected into the grid and on the chosen parameters of the numerical scheme being used. This paper addresses the issue of selecting and designing sources for finite difference simulation, by first reviewing associated aims and constraints, and evaluating existing source models against these criteria. The process of exciting a model is generalized by introducing a system of three cascaded filters, respectively, characterizing the driving pulse, the source mechanics, and the injection of the resulting source function into the grid. It is shown that hard, soft, and transparent sources can be seen as special cases within this unified approach. Starting from the mechanics of a small pulsating sphere, a parametric source model is formulated by specifying suitable filters. This physically constrained source model is numerically consistent, does not scatter incoming waves, and is free from zero- and low-frequency artifacts. Simulation results are employed for comparison with existing source formulations in terms of meeting the spectral and temporal requirements on the outward propagating wave.

Criticality, factorization, and long-range correlations in the anisotropic XY model

We study the long-range quantum correlations in the anisotropic XY model. By first examining the thermodynamic limit, we show that employing the quantum discord as a figure of merit allows one to capture the main features of the model at zero temperature. Furthermore, by considering suitably large site separations we find that these correlations obey a simple scaling behavior for finite temperatures, allowing for efficient estimation of the critical point. We also address ground-state factorization of this model by explicitly considering finite-size systems, showing its relation to the energy spectrum and explaining the persistence of the phenomenon at finite temperatures. Finally, we compute the fidelity between finite and infinite systems in order to show that remarkably small system sizes can closely approximate the thermodynamic limit.

Effect of the crown design and interface lute parameters on the stress-state of a machined crown-tooth system: A finite element analysis

The effect of preparation design and the physical properties of the interface lute on the restored machined ceramic crown-tooth complex are poorly understood. The aim of this work was to determine, by means of three-dimensional finite element analysis (3D FEA) the effect of the tooth preparation design and the elastic modulus of the cement on the stress state of the cemented machined ceramic crown-tooth complex. The three-dimensional structure of human premolar teeth, restored with adhesively cemented machined ceramic crowns, was digitized with a micro-CT scanner. An accurate, high resolution, digital replica model of a restored tooth was created. Two preparation designs, with different occlusal morphologies, were modeled with cements of 3 different elastic moduli. Interactive medical image processing software (mimics and professional CAD modeling software) was used to create sophisticated digital models that included the supporting structures; periodontal ligament and alveolar bone. The generated models were imported into an FEA software program (hypermesh version 10.0, Altair Engineering Inc.) with all degrees of freedom constrained at the outer surface of the supporting cortical bone of the crown-tooth complex. Five different elastic moduli values were given to the adhesive cement interface 1.8 GPa, 4 GPa, 8 GPa, 18.3 GPa and 40 GPa; the four lower values are representative of currently used cementing lutes and 40 GPa is set as an extreme high value. The stress distribution under simulated applied loads was determined. The preparation design demonstrated an effect on the stress state of the restored tooth system. The cement elastic modulus affected the stress state in the cement and dentin structures but not in the crown, the pulp, the periodontal ligament or the cancellous and cortical bone. The results of this study suggest that both the choice of the preparation design and the cement elastic modulus can affect the stress state within the restored crown-tooth complex.

Prediction of failure and fracture mechanisms of polymeric composites using finite element analysis. Part 1:Particulate filled composites

Micro-mechanical analysis of polymeric composites provides a powerful means for the quantitative assessment of their bulk behavior. In this paper we describe a robust finite element model (FEM) for the micro-structural modeling of the behavior of particulate filled polymer composites under external loads. The developed model is applied to simulate stress distribution in polymer composites containing particulate fillers. Quantitative information about the magnitude and location of maximum stress concentrations obtained from these simulations is used to predict the dominant failure and crack growth mechanisms in these composites. The model predictions are compared with the available experimental data and also with the values found using other methods reported in the literature. These comparisons show the range of the validity of the developed model and its predictive potential.

Finite element analysis of interfacial friction and slip in composites with fully unbonded filler

A finite element model is developed to predict the stress-strain behaviour of particulate composites with fully unbonded filler particles. This condition can occur because of the lack of adhesion property of the filler surface. Whilst part of the filler particle is separated from the matrix, another section of filler keeps in contact with the matrix because of the lateral compressive displacement of the matrix. The slip boundary condition is imposed on the section of the interface that remains closed. The states of stress and displacement fields are obtained. The location of any further deformation through crazing or shear band formation is identified. A completely unbonded inclusion with partial slip at a section of the interface reduces the concentration of the stress at the interface significantly. Whereas this might lead to slightly higher strength, it decreases the load transfer efficiency and stiffness of this type of composite.