Convective heat dissipation with lattice-frame materials

This paper presents experimental results on heat transfer and pressure drop for a compact heat sink made of fully triangulated, lightweight (porosity∼0.938), aluminum lattice-frame materials (LFMs). Due to the inherent structural anisotropy of the LFMs, two mutually perpendicular orientations were selected for the measurements. Constant heat flux was applied to the heat sink under steady state conditions, and dissipated by forced air convection. The experimental data were compared with those predicted from an analytical model based on fin analogy. The experimental results revealed that pressure drop is strongly dependent upon the orientation of the structure, due mainly to the flow blockage effect. For heat transfer measurements, typical local temperature distributions on the substrate under constant heat flux conditions were captured with infrared camera. The thermal behavior of LFMs was found to follow closely that of cylinder banks, with early transition Reynolds number (based on strut diameter) equal to about 300. The Nusselt number prediction from the fin-analogy correlates well with experimental measurements, except at low Reynolds numbers where a slightly underestimation is observed. Comparisons with empty channels and commonly used heat exchanger media show that the present LFM heat sink can remove heat approximately seven times more efficient than an empty channel and as efficient as a bank of cylinders at the same porosity level. The aluminum LFMs are extremely stiff and strong, making them ideal candidates for multifunctional structures requiring both heat dissipation and mechanical load carrying capabilities. © 2003 Elsevier Ltd. All rights reserved.

Properties of a-C:H films deposited from a methane electron cyclotron wave resonant plasma

An electron cyclotron wave resonant methane plasma discharge was used for the high rate deposition of hydrogenated amorphous carbon (a-C:H). Deposition rates of up to ∼400 Å/min were obtained over substrates up to 2.5 in. in diameter with a film thickness uniformity of ∼±10%. The deposited films were characterised in terms of their mass density, sp3 and hydrogen contents, C-H bonding, intrinsic stress, scratch resistance and friction properties. The deposited films possessed an average sp3 content, mass density and refractive index of ∼58%, 1.76 g/cm3 and 2.035 respectively.Mechanical characterisation indicated that the films possessed very low steady-state coefficients of friction (ca. 0.06) and a moderate shear strength of ∼141 MPa. Nano-indentation measurements also indicated a hardness and elastic modulus of ∼16.1 and 160 GPa respectively. The critical loads required to induce coating failure were also observed to increase with ion energy as a consequence of the increase in degree of ion mixing at the interface. Furthermore, coating failure under scratch test conditions was observed to take place via fracture within the silicon substrate itself, rather than either in the coating or at the film/substrate interface. © 2003 Elsevier B.V. All rights reserved.

Characterization of the Fracture Work for Ductile Film Undergoing the Micro-Scratch

The interface adhesion strength (or interface toughness) of a thin film/substrate system is often assessed by the micro-scratch test. For a brittle film material, the interface adhesion strength is easily obtained through measuring the scratch driving forces. However, to measure the interface adhesion strength (or interface toughness) for a metal thin film material (the ductile material) by the microscratch test is very difficult, because intense plastic deformation is involved and the problem is a three-dimensional elastic-plastic one. In the present research, using a double-cohesive zone model, the failure characteristics of the thin film/substrate system can be described and further simulated. For a steady-state scratching process, a three-dimensional elastic-plastic finite element method based on the double cohesive zone model is developed and adopted, and the steady-state fracture work of the total system is calculated. The parameter relations between the horizontal driving forces (or energy release rate of the scratching process) and the separation strength of thin film/substrate interface, and the material shear strength, as well as the material parameters are developed. Furthermore, a scratch experiment for the Al/Si film/substrate system is carried out and the failure mechanisms are explored. Finally, the prediction results are applied to a scratch experiment for the Pt/NiO material system given in the literature.

Robust approach to dynamic statistical process control

Statistical Process Control (SPC) technique are well established across a wide range of industries. In particular, the plotting of key steady state variables with their statistical limit against time (Shewart charting) is a common approach for monitoring the normality of production. This paper aims with extending Shewart charting techniques to the quality monitoring of variables driven by uncertain dynamic processes, which has particular application in the process industries where it is desirable to monitor process variables on-line as well as final product. The robust approach to dynamic SPC is based on previous work on guaranteed cost filtering for linear systems and is intended to provide a basis for both a wide application of SPC monitoring and also motivate unstructured fault detection.

Determination Ofinterfacial Mechanical Parameters For An Al/Epoxy/Al2O3 System Byusing Peel Test Simulations

Peel test measurements and simulations of the interfacial mechanical parameters for the Al/Epoxy/Al2O3 system are performed in the present investigation. A series of Al film thicknesses between 20 and 250 microns and three peel angles of 90, 135 and 180 degrees are considered. Two types of epoxy adhesives are adopted to obtain both strong and weak interface adhesions. A finite element model with cohesive zone elements is used to identify the interfacial parameters and simulate the peel test process. By simulating and recording normal stress near the crack tip, the separation strength is obtained. Furthermore, the cohesive energy is identified by comparing the simulated steady-state peel force and the experimental result. It is found from the research that both the cohesive energy and the separation strength can be taken as the intrinsic interfacial parameters which are dependent on the thickness of the adhesive layer and independent of the film thickness and peel angle.

Peeling Experiments Of Ductile Thin Films Along Ceramic Substrates - Critical Assessment Of Analytical Models

Two types of peeling experiments are performed in the present research. One is for the Al film/Al2O3 substrate system with an adhesive layer between the film and the substrate. The other one is for the Cu film/Al2O3 substrate system without adhesive layer between the film and the substrate, and the Cu films are electroplated onto the Al2O3 substrates. For the case with adhesive layer, two kinds of adhesives are selected, which are all the mixtures of epoxy and polyimide with mass ratios 1:1.5 and 1:1, respectively. The relationships between energy release rate, the film thickness and the adhesive layer thickness are measured during the steady-state peeling process. The effects of the adhesive layer on the energy release rate are analyzed. Using the experimental results, several analytical criteria for the steady-state peeling based on the bending model and on the two-dimensional finite element analysis model are critically assessed. Through assessment of analytical models, we find that the cohesive zone criterion based on the beam bend model is suitable for a weak interface strength case and it describes a macroscale fracture process zone case, while the two-dimensional finite element model is effective to both the strong interface and weak interface, and it describes a small-scale fracture process zone case. (C) 2007 Elsevier Ltd. All rights reserved.

Shock wave reflection from a wedge in a dusty gas

The present paper contains a detailed study of shock wave reflection from a wedge placed in various suspensions. In past works, the incident shock propagated initially in pure gas and the suspension started only at the leading edge of the deflecting wedge. However, in the present case the entire flow field is filled with a gas-dust suspension and the initial shock wave has steady-state structure relative to the shock front. In former studies the transmitted shock wave starts its propagation into the suspension and is reflected from the wedge at the same time. It is therefore obvious that the two unrelated processes of (2D) reflection and (1D) "transitional" relaxation occur simultaneously. In the present case the suspension behind the incident shock wave has reached steady state (i.e., it is a traveling wave) before the shock reaches the wedge leading edge. The reflection process from the deflecting wedge is studied for different dust mass loadings and different dust-particle diameter. It is shown that when the dust loading is low and the dust particle diameter is small the wave reflection pattern is similar to that observed in a similar pure gas case. In addition, an equilibrium state is reached, behind the evolved waves, very quickly. On the other hand, when the dust loading is relatively high and/or the dust particle diameter is relatively large, the observed reflection wave pattern is very different from that seen in a similar pure gas case. In such cases it takes much longer time to reach an equilibrium state behind the reflecting waves. It is also shown that the dust presence significantly affects the (gas) pressure on the wedge surface. The higher the dust loading is, the higher the pressure on the wedge surface. Suspensions composed of solid particle of different size, but having the same dust mass loading, will approach the same equilibrium pressure. However, it will take longer time to reach an equilibrium state for suspensions having large diameter particles. (C) 2004 Elsevier Ltd. All rights reserved.

Numerical simulation of Richtmyer-Meshkov instability

The compressible Navier-Stokes equations discretized with a fourth order accurate compact finite difference scheme with group velocity control are used to simulate the Richtmyer-Meshkov (R-M) instability problem produced by cylindrical shock-cylindrical material interface with shock Mach number Ms = 1.2 and density ratio 1:20 (interior density/outer density). Effect of shock refraction, reflection, interaction of the reflected shock with the material interface, and effect of initial perturbation modes on R-M instability are investigated numerically. It is noted that the shock refraction is a main physical mechanism of the initial phase changing of the material surface. The multiple interactions of the reflected shock from the origin with the interface and the R-M instability near the material interface are the reason for formation of the spike-bubble structures. Different viscosities lead to different spike-bubble structure characteristics. The vortex pairing phenomenon is found in the initial double mode simulation. The mode interaction is the main factor of small structures production near the interface.

Internal resonance of an L-shaped beam with a limit stop .1. Free vibration

Real-life structures often possess piecewise stiffness because of clearances or interference between subassemblies. Such an aspect can alter a system's fundamental free vibration response and leads to complex mode interaction. The free vibration behaviour of an L-shaped beam with a limit stop is analyzed by using the frequency response function and the incremental harmonic balance method. The presence of multiple internal resonances, which involve interactions among the first five modes and are extremely complex, have been discovered by including higher harmonics in the analysis. The results show that mode interaction may occur if the higher harmonics of a vibration mode are close to the natural frequency of a higher mode. The conditions for the existence of internal resonance are explored, and it is shown that a prerequisite is the presence of bifurcation points in the form of intersecting backbone curves. A method to compute such intersections by using only one harmonic in the free vibration solution is proposed. (C) 1996 Academic Press Limited

Concentration distribution in crystallization from solution under microgravity

Concentration distribution in crystallization from solution under microgravity is numerically studied. A quasi-steady state growth and dissolution in a 2D rectangular enclosure filled with sodium chlorate (NaClO3) aqueous solution, in which one wall is the growth surface of the crystal and the opposite one is the dissolution surface, is considered. The solute transport process at the growth surface is described by the diffusion-reaction theory with finite interface kinetics coefficient. The results show that the concentration at the growth surface is supersaturated and the supersaturation distribution is of non-uniformity, i.e. the supersaturation in a region facing an incoming flow is high. On the other hand, the non-uniformity of supersaturation at the growth surface is closely related to the gravity level even under microgravity, it exponentially increases as the thermal Rayleigh number on behalf of the gravity level rises.

Kinetic theory for particle concentration distribution in two phase flow

The results of experiments in open channels and closed pipelines show two kinds of patterns for the vertical distribution of particle concentration (i.e., pattern I and pattern II). The former shows a pattern of maximum concentration at some location above the bottom and the downward decay of the concentration below the location. The latter always shows an increase of the particle concentration downward over the whole vertical, with the maximum value at the bottom. Many investigations were made on the pattern II, but few were made on pattern I. In this paper, a particle velocity distribution function is first obtained in the equilibrium state or in dilute steady state for the particle in two-phase flows, then a theoretical model for the particle concentration distribution is derived from the kinetic theory. More attention is paid to the predictions of the concentration distribution of pattern I and comparisons of the present model are made with the data measured by means of laser doppler anemometry (LDA). Very good agreements are obtained between the measured and calculated results.

Numerical experiments on one-dimensional model of turbulence

The initial-value problem of a forced Burgers equation is numerically solved by the Fourier expansion method. It is found that its solutions finally reach a steady state of 'laminar flow' which has no randomness and is stable to disturbances. Hence, strictly speaking, the so-called Burgers turbulence is not a turbulence. A new one-dimensional model is proposed to simulate the Navier-Stokes turbulence. A series of numerical experiments on this one-dimensional turbulence is made and is successful in obtaining Kolmogorov's (1941) k exp(-5/3) inertial-range spectrum. The (one-dimensional) Kolmogorov constant ranges from 0.5 to 0.65.

Fiscal Policy and Economic Stability: Does PIGS stand for Procyclicality In Government Spending?

The Financial Crisis has hit particularly hard countries like Ireland or Spain. Procyclical fiscal policy has contributed to a boom-bust cycle that undermined fiscal positions and deepened current account deficits during the boom. We set up an RBC model of a small open economy, following Mendoza (1991), and introduce the effect of fiscal policy decisions that change over the cycle. We calibrate the model on data for Ireland, and simulate the effect of different spending policies in response to supply shocks. Procyclical fiscal policy distorts intertemporal allocation decisions. Temporary spending boosts in booms spur investment, and hence the need for external finance, and so generates very volatile cycles in investment and the current account. This economic instability is also harmful for the steady state level of output. Our model is able to replicate the relation between the degree of cyclicality of fiscal policy, and the volatility of consumption, investment and the current account observed in OECD countries.

Reference Points and Optimal Management in Stochastic Age-Structured Fisheries Models

The purpose of this article is to characterize dynamic optimal harvesting trajectories that maximize discounted utility assuming an age-structured population model, in the same line as Tahvonen (2009). The main novelty of our study is that uses as an age-structured population model the standard stochastic cohort framework applied in Virtual Population Analysis for fish stock assessment. This allows us to compare optimal harvesting in a discounted economic context with standard reference points used by fisheries agencies for long term management plans (e.g. Fmsy). Our main findings are the following. First, optimal steady state is characterized and sufficient conditions that guarantees its existence and uniqueness for the general case of n cohorts are shown. It is also proved that the optimal steady state coincides with the traditional target Fmsy when the utility function to be maximized is the yield and the discount rate is zero. Second, an algorithm to calculate the optimal path that easily drives the resource to the steady state is developed. And third, the algorithm is applied to the Northern Stock of hake. Results show that management plans based exclusively on traditional reference targets as Fmsy may drive fishery economic results far from the optimal.

Growth in Overlapping Generation Economies with Non-Renewable Resources

Published as an article in: Journal of Environmental Economics and Management, 2005, vol. 50, issue 2, pages 387-407.