Microwave heating in multiphase systems: evaluation of series solutions

The 1D electric field and heat-conduction equations are solved for a slab where the dielectric properties vary spatially in the sample. Series solutions to the electric field are obtained for systems where the spatial variation in the dielectric properties can be expressed as polynomials. The series solution is used to obtain electric-field distributions for a binary oil-water system where the dielectric properties are assumed to vary linearly within the sample. Using the finite-element method temperature distributions are computed in a three-phase oil, water and rock system where the dielectric properties vary due to the changing oil saturation in the rock. Temperature distributions predicted using a linear variation in the dielectric properties are compared with those obtained using the exact nonlinear variation.

Role of length scales on microwave thawing dynamics in 2D cylinders

Microwave (MW) thawing of 2D frozen cylinders exposed to uniform plane waves from one face, is modeled using the effective heat capacity formulation with the MW power obtained from the electric field equations. Computations are illustrated for tylose (23% methyl cellulose gel) which melts over a range of temperatures giving rise to a mushy zone. Within the mushy region the dielectric properties are functions of the liquid volume fraction. The resulting coupled, time dependent non-linear equations are solved using the Galerkin finite element method with a fixed mesh. Our method efficiently captures the multiple connected thawed domains that arise due to the penetration of MWs in the sample. For a cylinder of diameter D, the two length scales that control the thawing dynamics are D/D-p and D/lambda(m), where D-p and lambda(m) are the penetration depth and wavelength of radiation in the sample respectively. For D/D-p, D/lambda(m) much less than 1 power absorption is uniform and thawing occurs almost simultaneously across the sample (Regime I). For D/D-p much greater than 1 thawing is seen to occur from the incident face, since the power decays exponentially into the sample (Regime III). At intermediate values, 0.2 < D/D-p, D/lambda(m) < 2.0 (Regime II) thawing occurs from the unexposed face at smaller diameters, from both faces at intermediate diameters and from the exposed and central regions at larger diameters. Average power absorption during thawing indicates a monotonic rise in Regime I and a monotonic decrease in Regime III. Local maxima in the average power observed for samples in Regime II are due to internal resonances within the sample. Thawing time increases monotonically with sample diameter and temperature gradients in the sample generally increase from Regime I to Regime III. (C) 2002 Elsevier Science Ltd. All rights reserved.

Mapping of large scale 158 mu m [CII] line emission: Orion A

We present the first results of an observational programme undertaken to map the fine structure line emission of singly ionized carbon ([ CII] 157 : 7409 mum) over extended regions using a Fabry Perot spectrometer newly installed at the focal plane of a 100 cm balloon- borne far- infrared telescope. This new combination of instruments has a velocity resolution of similar to 200 km s(-1) and an angular resolution of 1.'5. During the first flight, an area of 30' x 15' in Orion A was mapped. These observations extend over a larger area than previous observations, the map is fully sampled and the spectral scanning method used enables reliable estimation of the continuum emission at frequencies adjacent to the [ CII] line. The total [ CII] line luminosity, calculated by considering up to 20% of the maximum line intensity is 0.04% of the luminosity of the far- infrared continuum. We have compared the [ CII] intensity distribution with the velocity- integrated intensity distributions of (CO)-C-13(1- 0), CI(1- 0) and CO( 3- 2) from the literature. Comparison of the [ CII], [ CI] and the radio continuum intensity distributions indicates that the largescale [ CII] emission originates mainly from the neutral gas, except at the position of M 43, where no [ CI] emission corresponding to the [ CII] emission is seen. Substantial part of the [ CII] emission from here originates from the ionized gas. The observed line intensities and ratios have been analyzed using the PDR models by Kaufman et al. ( 1999) to derive the incident UV flux and volume density at a few selected positions. The models reproduce the observations reasonably well at most positions excepting the [ CII] peak ( which coincides with the position of theta(1) Ori C). Possible reason for the failure could be the simplifying assumption of a homogeneous plane parallel slab in place of a more complicated geometry.

Effects of phase inhomogeneity and boundary conditions on the dynamic response of SMA wire actuators

This paper reports the simulation results from the dynamic analysis of a Shape Memory Alloy (SMA) actuator. The emphasis is on understanding the dynamic behavior under various loading rates and boundary conditions, resulting in complex scenarios such as thermal and stress gradients. Also, due to the polycrystalline nature of SMA wires, presence of microstructural inhomogeneity is inevitable. Probing the effect of inhomogeneity on the dynamic behavior can facilitate the prediction of life and characteristics of SMA wire actuator under varieties of boundary and loading conditions. To study the effect of these factors, an initial boundary value problem of SMA wire is formulated. This is subsequently solved using finite element method. The dynamic response of the SMA wire actuator is analyzed under mechanical loading and results are reported. Effect of loading rate, micro-structural inhomogeneity and thermal boundary conditions on the dynamic response of SMA wire actuator is investigated and the simulation results are reported.

Modified lattice model for mode-I fracture analysis of notched plain concrete beam using probabilistic approach

A modified lattice model using finite element method has been developed to study the mode-I fracture analysis of heterogeneous materials like concrete. In this model, the truss members always join at points where aggregates are located which are modeled as plane stress triangular elements. The truss members are given the properties of cement mortar matrix randomly, so as to represent the randomness of strength in concrete. It is widely accepted that the fracture of concrete structures should not be based on strength criterion alone, but should be coupled with energy criterion. Here, by incorporating the strain softening through a parameter ‘α’, the energy concept is introduced. The softening branch of load-displacement curves was successfully obtained. From the sensitivity study, it was observed that the maximum load of a beam is most sensitive to the tensile strength of mortar. It is seen that by varying the values of properties of mortar according to a normal random distribution, better results can be obtained for load-displacement diagram.

Implementation of a phase array diffuse optical tomographic imager

Diffuse optical tomography (DOT) using near-infrared (NIR) light is a promising tool for noninvasive imaging of deep tissue. This technique is capable of quantitative reconstructions of absorption coefficient inhomogeneities of tissue. The motivation for reconstructing the optical property variation is that it, and, in particular, the absorption coefficient variation, can be used to diagnose different metabolic and disease states of tissue. In DOT, like any other medical imaging modality, the aim is to produce a reconstruction with good spatial resolution and accuracy from noisy measurements. We study the performance of a phase array system for detection of optical inhomogeneities in tissue. The light transport through a tissue is diffusive in nature and can be modeled using diffusion equation if the optical parameters of the inhomogeneity are close to the optical properties of the background. The amplitude cancellation method that uses dual out-of-phase sources (phase array) can detect and locate small objects in turbid medium. The inverse problem is solved using model based iterative image reconstruction. Diffusion equation is solved using finite element method for providing the forward model for photon transport. The solution of the forward problem is used for computing the Jacobian and the simultaneous equation is solved using conjugate gradient search. The simulation studies have been carried out and the results show that a phase array system can resolve inhomogeneities with sizes of 5 mm when the absorption coefficient of the inhomogeneity is twice that of the background tissue. To validate this result, a prototype model for performing a dual-source system has been developed. Experiments are carried out by inserting an inhomogeneity of high optical absorption coefficient in an otherwise homogeneous phantom while keeping the scattering coefficient same. The high frequency (100 MHz) modulated dual out-of-phase laser source light is propagated through the phantom. The interference of these sources creates an amplitude null and a phase shift of 180° along a plane between the two sources with a homogeneous object. A solid resin phantom with inhomogeneities simulating the tumor is used in our experiment. The amplitude and phase changes are found to be disturbed by the presence of the inhomogeneity in the object. The experimental data (amplitude and the phase measured at the detector) are used for reconstruction. The results show that the method is able to detect multiple inhomogeneities with sizes of 4 mm. The localization error for a 5 mm inhomogeneity is found to be approximately 1 mm.

Diffuse optical tomographic imager using a single light source

Near-infrared diffuse optical tomography (DOT) technique has the capability of providing good quantitative reconstruction of tissue absorption and scattering properties with additional inputs such as input and output modulation depths and correction for the photon leakage. We have calculated the two-dimensional (2D) input modulation depth from three-dimensional (3D) diffusion to model the 2D diffusion of photons. The photon leakage when light traverses from phantom to the fiber tip is estimated using a solid angle model. The experiments are carried for single (5 and 6 mm) as well as multiple inhomogeneities (6 and 8 mm) with higher absorption coefficient in a homogeneous phantom. Diffusion equation for photon transport is solved using finite element method and Jacobian is modeled for reconstructing the optical parameters. We study the development and performance of DOT system using modulated single light source and multiple detectors. The dual source methods are reported to have better reconstruction capabilities to resolve and localize single as well as multiple inhomogeneities because of its superior noise rejection capability. However, an experimental setup with dual sources is much more difficult to implement because of adjustment of two out of phase identical light probes symmetrically on either side of the detector during scanning time. Our work shows that with a relatively simpler system with a single source, the results are better in terms of resolution and localization. The experiments are carried out with 5 and 6 mm inhomogeneities separately and 6 and 8 mm inhomogeneities both together with absorption coefficient almost three times as that of the background. The results show that our experimental single source system with additional inputs such as 2D input/output modulation depth and air fiber interface correction is capable of detecting 5 and 6 mm inhomogeneities separately and can identify the size difference of multiple inhomogeneities such as 6 and 8 mm. The localization error is zero. The recovered absorption coefficient is 93% of inhomogeneity that we have embedded in experimental phantom.

Gap model of one-way cyclic lateral load on vertical files in soft clay

For the analysis and design of pile foundation used for coastal structures the prediction of cyclic response, which is influenced by the nonlinear behavior, gap (pile soil separation) and degradation (reduction in strength) of soil becomes necessary. To study the effect of the above parameters a nonlinear cyclic load analysis program using finite element method is developed, incorporating the proposed gap and degradation model and adopting an incremental-iterative procedure. The pile is idealized using beam elements and the soil by number of elastoplastic sub-element springs at each node. The effect of gap and degradation on the load-deflection behavior. elasto-plastic sub-element and resistance of the soil at ground-line have been clearly depicted in this paper.

Structure and deformation correlation of closed-cell aluminium foam subject to uniaxial compression

We report the results of an experimental and numerical study conducted on a closed-cell aluminium foam that was subjected to uniaxial compression with lateral constraint. X-ray computed tomography was utilized to gain access into the three-dimensional (3-D) structure of the foam and some aspects of the deformation mechanisms. A series of advanced 3-D image analyses are conducted on the 3-D images aimed at characterizing the strain localization regions. We identify the morphological/geometrical features that are responsible for the collapse of the cells and the strain localization. A novel mathematical approach based on a Minkowski tensor analysis along with the mean intercept length technique were utilized to search for signatures of anisotropy across the foam sample and its evolution as a function of loading. Our results show that regions with higher degrees of anisotropy in the undeformed foam have a tendency to initiate the onset of cell collapse. Furthermore, we show that strain hardening occurs predominantly in regions with large cells and high anisotropy. We combine the finite element method with the tomographic images to simulate the mechanical response of the foam. We predict further deformation in regions where the foam is already deformed. Crown Copyright (C) 2012 Published by Elsevier Ltd. on behalf of Acta Materialia Inc. All rights reserved.

Ultrasonic Guided Wave Characterization and Damage Detection in Foam-core Sandwich Panel using PWAS and LDV

Lamb wave type guided wave propagation in foam core sandwich structures and detectability of damages using spectral analysis method are reported in this paper. An experimental study supported by theoretical evaluation of the guided wave characteristics is presented here that shows the applicability of Lamb wave type guided ultrasonic wave for detection of damage in foam core sandwich structures. Sandwich beam specimens were fabricated with 10 mm thick foam core and 0.3 mm thick aluminum face sheets. Thin piezoelectric patch actuators and sensors are used to excite and sense guided wave. Group velocity dispersion curves and frequency response of sensed signal are obtained experimentally. The nature of damping present in the sandwich panel is monitored by measuring the sensor signal amplitude at various different distances measured from the center of the linear phased array. Delaminations of increasing width are created and detected experimentally by pitch-catch interrogation with guided waves and wavelet transform of the sensed signal. Signal amplitudes are analyzed for various different sizes of damages to differentiate the damage size/severity. A sandwich panel is also fabricated with a planer dimension of 600 mm x 400 mm. Release film delamination is introduced during fabrication. Non-contact Laser Doppler Vibrometer (LDV) is used to scan the panel while exciting with a surface bonded piezoelectric actuator. Presence of damage is confirmed by the reflected wave fringe pattern obtained from the LDV scan. With this approach it is possible to locate and monitor the damages by tracking the wave packets scattered from the damages.

Transition due to base roughness in a dense granular flow down an inclined plane

Particle simulations based on the discrete element method are used to examine the effect of base roughness on the granular flow down an inclined plane. The base is composed of a random configuration of fixed particles, and the base roughness is decreased by decreasing the ratio of diameters of the base and moving particles. A discontinuous transition from a disordered to an ordered flow state is observed when the ratio of diameters of base and moving particles is decreased below a critical value. The ordered flowing state consists of hexagonally close packed layers of particles sliding over each other. The ordered state is denser (higher volume fraction) and has a lower coordination number than the disordered state, and there are discontinuous changes in both the volume fraction and the coordination number at transition. The Bagnold law, which states that the stress is proportional to the square of the strain rate, is valid in both states. However, the Bagnold coefficients in the ordered flowing state are lower, by more than two orders of magnitude, in comparison to those of the disordered state. The critical ratio of base and moving particle diameters is independent of the angle of inclination, and varies very little when the height of the flowing layer is doubled from about 35 to about 70 particle diameters. While flow in the disordered state ceases when the angle of inclination decreases below 20 degrees, there is flow in the ordered state at lower angles of inclination upto 14 degrees. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4710543]

Effect of isothermal aging and electromigration on the microstructural evolution of solder interconnections during thermomechanical loading

The effect of different pre-aging treatments on the microstructural evolution of lead-free solder and growth of interfacial intermetallic compound layers under thermal cycling has been investigated in this work. The results show that there are distinct differences in the microstructural changes between samples with no pretreatment, samples that have experienced thermal annealing at 125A degrees C for 750 h before thermal cycling, and those that have had direct current (DC) stressing for 750 h as pretreatment. The microstructural evolution of the solder matrix is rationalized by utilizing the science of microstructures and analysis of the influence of electron flow on the precipitation phenomena. The finite-element method is utilized to understand the loading conditions imposed on the solder interconnections during cyclic stressing. The growth of intermetallic reaction layers is further analyzed by utilizing quantitative thermodynamic calculations coupled with kinetic analysis. The latter is based on the changes in the intrinsic diffusion fluxes of elements induced by current flow and alloying elements present in the system. With this concurrent approach the differences seen in thermal cycling behavior between the different pre-aging treatments can be explained.

Size-dependent free flexural vibration behavior of functionally graded nanoplates

In this paper, size dependent linear free flexural vibration behavior of functionally graded (FG) nanoplates are investigated using the iso-geometric based finite element method. The field variables are approximated by non-uniform rational B-splines. The nonlocal constitutive relation is based on Eringen's differential form of nonlocal elasticity theory. The material properties are assumed to vary only in the thickness direction and the effective properties for the FG plate are computed using Mori-Tanaka homogenization scheme. The accuracy of the present formulation is demonstrated considering the problems for which solutions are available. A detailed numerical study is carried out to examine the effect of material gradient index, the characteristic internal length, the plate thickness, the plate aspect ratio and the boundary conditions on the global response of the FG nanoplate. From the detailed numerical study it is seen that the fundamental frequency decreases with increasing gradient index and characteristic internal length. (c) 2012 Elsevier B.V. All rights reserved.

An adaptive singular ES-FEM for mechanics problems with singular field of arbitrary order

This paper presents a singular edge-based smoothed finite element method (sES-FEM) for mechanics problems with singular stress fields of arbitrary order. The sES-FEM uses a basic mesh of three-noded linear triangular (T3) elements and a special layer of five-noded singular triangular elements (sT5) connected to the singular-point of the stress field. The sT5 element has an additional node on each of the two edges connected to the singular-point. It allows us to represent simple and efficient enrichment with desired terms for the displacement field near the singular-point with the satisfaction of partition-of-unity property. The stiffness matrix of the discretized system is then obtained using the assumed displacement values (not the derivatives) over smoothing domains associated with the edges of elements. An adaptive procedure for the sES-FEM is proposed to enhance the quality of the solution with minimized number of nodes. Several numerical examples are provided to validate the reliability of the present sES-FEM method. (C) 2012 Elsevier B.V. All rights reserved.

Application of Artificial Neural Networks and Particle Swarm Optimization for the Management of Groundwater Resources

Ground management problems are typically solved by the simulation-optimization approach where complex numerical models are used to simulate the groundwater flow and/or contamination transport. These numerical models take a lot of time to solve the management problems and hence become computationally expensive. In this study, Artificial Neural Network (ANN) and Particle Swarm Optimization (PSO) models were developed and coupled for the management of groundwater of Dore river basin in France. The Analytic Element Method (AEM) based flow model was developed and used to generate the dataset for the training and testing of the ANN model. This developed ANN-PSO model was applied to minimize the pumping cost of the wells, including cost of the pipe line. The discharge and location of the pumping wells were taken as the decision variable and the ANN-PSO model was applied to find out the optimal location of the wells. The results of the ANN-PSO model are found similar to the results obtained by AEM-PSO model. The results show that the ANN model can reduce the computational burden significantly as it is able to analyze different scenarios, and the ANN-PSO model is capable of identifying the optimal location of wells efficiently.