Unsteady free convection flow under the influence of a magnetic field

The unsteady free convection boundary layer at the stagnation point of a two-dimensional body and an axisymmetric body with prescribed surface heat flux or temperature has been studied. The magnetic field is applied parallel to the surface and the effect of induced magnetic field has been considered. It is found that for certain powerlaw distribution of surface heat flux or temperature and magnetic field with time, the governing boundary layer equations admit a self-similar solution locally. The resulting nonlinear ordinary differential equations have been solved using a finite element method and a shooting method with Newton's corrections for missing initial conditions. The results show that the skin friction and heat transfer coefficients, and x-component of the induced magnetic field on the surface increase with the applied magnetic field. In general, the skin friction, heat transfer and x-component of the induced magnetic field for axisymmetric case are more than those of the two-dimensional case. Also they change more when the surface heat flux or temperature decreases with time than when it increases with time. The skin friction, heat transfer and x-component of the induced magnetic field are significantly affected by the magnetic Prandtl number and they increase as the magnetic Prandtl number decreases. The skin friction and x-component of the magnetic field increase with the dissipation parameter, but heat transfer decreases.

Plates on Elastic Foundation

Analysis of rectangular plates resting on a Winkler-type, one-parameter foundation is studied. The finite element method is applied and a 12-degree-of-freedom, nonconforming rectangular plate element is adopted. Based on shape functions of the plate element, an energy approach is used to derive a closed-form, 12-by-12, consistent foundation stiffness matrix for a rectangular plate on an elastic subgrade. A commonly used method of modeling structural elements on an elastic foundation is the application of discrete springs at the element nodes. The model developed in this paper is compared with the discrete spring model and the convergence of both models is discussed. The convergence of the models is compared with the well-known classical solution of plates on elastic foundation developed in the 1950s. Both models show good convergence to the classical solution. The continuous subgrade response model converges in a manner more consistent with the flexibility of the plate element.

Stability of laminated composite plates with cut-outs

Critical buckling loads of laminated fibre-reinforced plastic square panels have been obtained using the finite element method. Various boundary conditions, lay-up details, fibre orientations, cut-out sizes are considered. A 36 degrees of freedom triangular element, based on the classical lamination theory (CLT) has been used for the analysis. The performance of this element is validated by comparing results with some of those available in literature. New results have been given for several cases of boundary conditions for [0°/ ± 45°/90°]s laminates. The effect of fibre-orientation in the ply on the buckling loads has been investigated by considering [±?]6s laminates.

Stiffened composite cylindrical panels�Optimum lay-up for buckling by ranking

Buckling of discretely stiffened composite cylindrical panels made of repeated sublaminate construction is studied using a finite element method. In repeated sublaminate construction, a full laminate is obtained by repeating a basic sublaminate, which has a smaller number of plies. This paper deals with the determination of the optimum lay-up for buckling by ranking of such stiffened (longitudinal and hoop) composite cylindrical panels. For this purpose we use the particularized form of a four-noded, 48 degrees of freedom doubly curved quadrilateral thin shell finite element together with a fully compatible two-noded, 16 degrees of freedom composite stiffener element. The computer program developed has been used, after extensive checking for correctness, to obtain an optimum orientation scheme of the plies in the sublaminate so as to achieve maximum buckling load for a specified thickness of typical stiffened composite cylindrical panels.

Role of surface tension driven convection in pulsed - laser melting and solidification

Surface melting by a stationary, pulsed laser has been modelled by the finite element method. The role of the surface tension driven convection is investigated in detail. Numerical results are presented for a triangular laser pulse of durations 10, 50 and 200 ms. Though the magnitude of the velocity is high due to the surface tension forces, the present results indicate that a finite time is required for convection to affect the temperature distribution within the melt pool. The effect of convection is very significant for pulse durations longer than 10 ms.

Synthesis of path generating flexible-link mechanisms

Flexible-link mechanisms are those linkage mechanisms (or structures) which are capable of motion by virtue of elastic deformation of one or more;links. In such mechanisms a single flexible link; can replace several rigid links and joints resulting in fewer links, fewer pin joints, reduced overall weight and reduced mechanical error. In spite of such clear advantages, contributions toward flexible-link mechanisms remain very scarce. The area of flexible-link mechanisms offers much scope for further exploration. This paper attempts to show the potential of flexible-link mechanisms in accomplishing a kinematic task like path generation. Synthesis of a four-bar mechanism with a flexible rocker for circular and straight line path generation is carried out. Displacement analysis of the structure is carried out using finite element method (FEM) and synthesis is formulated and solved as an optimization problem. Several numerical examples are presented for illustration. Based on the results obtained with these examples, the flexible-link mechanism considered shows good promise for-path generation.

Unsteady flow and heat transfer on a semi-infinite flat plate with an aligned magnetic field

The unsteady laminar boundary layer flow of an electrically conducting fluid past a semi-infinite flat plate with an aligned magnetic field has been studied when at time t > 0 the plate is impulsively moved with a constant velocity which is in the same or opposite direction to that of free stream velocity. The effect of the induced magnetic field has been included in the analysis. The non-linear partial differential equations have been solved numerically using an implicit finite-difference method. The effect of the impulsive motion of the surface is found to be more pronounced on the skin friction but its effect on the x-component of the induced magnetic field and heat transfer is small. Velocity defect occurs near the surface when the plate is impulsively moved in the same direction as that of the free stream velocity. The surface shear stress, x-component of the induced magnetic field on the surface and the surface heat transfer decrease with an increasing magnetic field, but they increase with the reciprocal of the magnetic Prandtl number. However, the effect of the reciprocal of the magnetic Prandtl number is more pronounced on the x-component of the induced magnetic field. (C) 1999 Elsevier Science Ltd. All rights reserved.

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.

Three-dimensional computational modeling of momentum, heat, and mass transfer in a laser surface alloying process

A three- dimensional, transient model is developed for studying heat transfer, fluid flow, and mass transfer for the case of a single- pass laser surface alloying process. The coupled momentum, energy, and species conservation equations are solved using a finite volume procedure. Phase change processes are modeled using a fixed-grid enthalpy-porosity technique, which is capable of predicting the continuously evolving solid- liquid interface. The three- dimensional model is able to predict the species concentration distribution inside the molten pool during alloying, as well as in the entire cross section of the solidified alloy. The model is simulated for different values of various significant processing parameters such as laser power, scanning speed, and powder feedrate in order to assess their influences on geometry and dynamics of the pool, cooling rates, as well as species concentration distribution inside the substrate. Effects of incorporating property variations in the numerical model are also discussed.

Turbulent flow computations using upwind RANS technique for a flat plate and an airfoil

Reynolds Averaged Navier Stokes (RANS) equations are solved using third order upwind biased Roe's scheme for the inviscid fluxes and second order central difference scheme for the viscous fluxes. The Baldwin & Lomax turbulence model is employed for Reynolds stresses. The governing equations are solved using finite-volume implicit scheme in body fitted curvilinear coordinate O-grid system. Computations axe reported for a flat plate apart from RAE 2822 and NACA 0012 airfoils. Results for the flat plate at M = 0.3, R-c = 4.0 x 10(6) compare favourably with the analytical solution. Results for the two airfoils are compared with experiment. There is a good agreement in C-p distribution between experiment and computation for both the airfoils. Comparison of C-f distribution with experiment for RAE 2822 airfoil is reasonable.

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.

Filtered density function for large eddy simulation of turbulent reacting flows

A methodology termed the “filtered density function” (FDF) is developed and implemented for large eddy simulation (LES) of chemically reacting turbulent flows. In this methodology, the effects of the unresolved scalar fluctuations are taken into account by considering the probability density function (PDF) of subgrid scale (SGS) scalar quantities. A transport equation is derived for the FDF in which the effect of chemical reactions appears in a closed form. The influences of scalar mixing and convection within the subgrid are modeled. The FDF transport equation is solved numerically via a Lagrangian Monte Carlo scheme in which the solutions of the equivalent stochastic differential equations (SDEs) are obtained. These solutions preserve the Itô-Gikhman nature of the SDEs. The consistency of the FDF approach, the convergence of its Monte Carlo solution and the performance of the closures employed in the FDF transport equation are assessed by comparisons with results obtained by direct numerical simulation (DNS) and by conventional LES procedures in which the first two SGS scalar moments are obtained by a finite difference method (LES-FD). These comparative assessments are conducted by implementations of all three schemes (FDF, DNS and LES-FD) in a temporally developing mixing layer and a spatially developing planar jet under both non-reacting and reacting conditions. In non-reacting flows, the Monte Carlo solution of the FDF yields results similar to those via LES-FD. The advantage of the FDF is demonstrated by its use in reacting flows. In the absence of a closure for the SGS scalar fluctuations, the LES-FD results are significantly different from those based on DNS. The FDF results show a much closer agreement with filtered DNS results. © 1998 American Institute of Physics.

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.