Non-linear boundary element analysis of floor slabs reinforced with rectangular beams

In this work, a numerical model to perform non-linear analysis of building floor structures is proposed. The presented model is derived from the Kirchhoff-s plate bending formulation of the boundary element method (BENI) for zoned domains, in which the plate stiffness is modified by the presence of membrane effects. In this model, no approximation of the generalized forces along the interface is required and the compatibility and equilibrium conditions along interfaces are imposed at the integral equation level. In order to reduce the number of degrees of freedom, the Navier Bernoulli hypothesis is assumed to simplify the strain field for the thin sub-regions (rectangular beams). The non-linear formulation is obtained from the linear formulation by incorporating initial internal force fields, which are approximated by using the well-known cell sub-division. Then, the non-linear solution of algebraic equations is obtained by using the concept of the consistent tangent operator. The Von Mises criterion is adopted to govern the elasto-plastic material behaviour checked at points along the plate thickness and along the rectangular beam element axes. The numerical representations are accurately obtained by either computing analytically the element integrals or performing the numerical integration accurately using an appropriate sub-elementation scheme. (C) 2007 Elsevier Ltd. All rights reserved.

A BEM formulation based on Reissner's theory to perform simple bending analysis of plates reinforced by rectangular beams

In this work, the plate bending formulation of the boundary element method (BEM), based on the Reissner's hypothesis, is extended to the analysis of plates reinforced by rectangular beams. This composed structure is modelled by a zoned plate, being the beams represented by narrow sub-regions with larger thickness. The integral equations are derived by applying the weighted residual method to each sub-region, and summing them to get the equation for the whole plate. Equilibrium and compatibility conditions are automatically imposed by the integral equations, which treat this composed structure as a single body. In order to decrease the number of degrees of freedom, some approximations are considered for both displacements and tractions along the beam width. The accuracy of the proposed model is illustrated by simple examples whose exact solution are known as well as by more complex examples whose numerical results are compared with a well-known finite element code.

A steady-state model for heat pipes of circular or rectangular cross-sections

The objective of this paper is to present a generalized analytical-numerical model of the internal flow in heat pipes. The model formulation is based on two-dimensional formulation of the energy and momentum equations in the vapour and liquid regions and also in the metallic tube. The numerical solution of the model is obtained by using the descretization scheme LOAD and the SIMPLE numerical code. The flow fields, as well as the pressure fields, for different geometries were obtained and discussed. Copyright © 1996 Elsevier Science Ltd.

Study on transient ice formation of laminar flow inside externally supercooled rectangular duct

The transient process of solidification of laminar liquid flow (water) submitted to super-cooling was investigated both theoretically and experimentally. In this study an alternative analytical formulation and numerical approach were adopted resulting in the unsteady model with temperature dependent thermophysical properties in the solid region. The proposed model is based upon the fundamental equations of energy balance in the solid and liquid regions as well as across the solidification front. The basic equations and the associated boundary and initial conditions were made dimensionless by using the Landau transformation to immobilize the moving front and render the problem to a fixed plane type problem. A laminar velocity profile is admitted in the liquid domain and the resulting equations were discretized using the finite difference approach. The numerical predictions obtained were compared with the available results based on other models and concepts such as Neumann analytical model, the apparent thermal capacity model due to Bonacina and the conventional fixed grid energy model due to Goodrich. To obtain further comparisons and more validation of the model and the numerical solution, an experimental rig was constructed and instrumented permitting very well controlled experimental measurements. The numerical predictions were compared with the experimental results and the agreement was found satisfactory.

Determination of resonant frequencies of triangular and rectangular microstrip antennas, using artificial neural networks

This paper presents models that can be used in the design of microstrip antennas for mobile communications. The antennas can be triangular or rectangular. The presented models are compared with deterministic and empirical models based on artificial neural networks (ANN) presented in the literature. The models are based on Perceptron Multilayer (PML) and Radial Basis Function (RBF) ANN. RBF based models presented the best results. Also, the models can be embedded in CAD systems, in order to design microstrip antennas for mobile communications.

New contributions to the determination of resonance frequencies of rectangular microstrip antennas by neural networks

This paper uses artificial neural networks (ANN) to compute the resonance frequencies of rectangular microstrip antennas (MSA), used in mobile communications. Perceptron Multi-layers (PML) networks were used, with the Quasi-Newton method proposed by Broyden, Fletcher, Goldfarb and Shanno (BFGS). Due to the nature of the problem, two hundred and fifty networks were trained, and the resonance frequency for each test antenna was calculated by statistical methods. The estimate resonance frequencies for six test antennas were compared with others results obtained by deterministic and ANN based empirical models from the literature, and presented a better agreement with the experimental values.

A quantum theory for the excitation spectrum of a rectangular Andreev billiard

We consider a N - S box system consisting of a rectangular conductor coupled to a superconductor. The Green functions are constructed by solving the Bogoliubov-de Gennes equations at each side of the interface, with the pairing potential described by a step-like function. Taking into account the mismatch in the Fermi wave number and the effective masses of the normal metal - superconductor and the tunnel barrier at the interface, we use the quantum section method in order to find the exact energy Green function yielding accurate computed eigenvalues and the density of states. Furthermore, this procedure allow us to analyze in detail the nontrivial semiclassical limit and examine the range of applicability of the Bohr-Sommerfeld quantization method.

Non-intrusive pressure measurement in microchannels

A non-intrusive interferometric measurement technique has been successfully developed to measure fluid compressibility in both gas and liquid phases via refractive index (RI) changes. The technique, consisting of an unfocused laser beam impinging a glass channel, can be used to separate and quantify cell deflection, fluid flow rates, and pressure variations in microchannels. Currently in fields such as microfluidics, pressure and flow rate measurement devices are orders of magnitude larger than the channel cross-sections making direct pressure and fluid flow rate measurements impossible. Due to the non-intrusive nature of this technique, such measurements are now possible, opening the door for a myriad of new scientific research and experimentation. This technique, adapted from the concept of Micro Interferometric Backscatter Detection (MIBD), boasts the ability to provide comparable sensitivities in a variety of channel types and provides quantification capability not previously demonstrated in backscatter detection techniques. Measurement sensitivity depends heavily on experimental parameters such as beam impingement angle, fluid volume, photodetector sensitivity, and a channel’s dimensional tolerances. The current apparatus readily quantifies fluid RI changes of 10-5 refractive index units (RIU) corresponding to pressures of approximately 14 psi and 1 psi in water and air, respectively. MIBD reports detection capability as low as 10-9 RIU and the newly adapted technique has the potential to meet and exceed this limit providing quantification in the place of detection. Specific device sensitivities are discussed and suggestions are provided on how the technique may be refined to provide optimal quantification capabilities based on experimental conditions.

Experimental investigation of regular fluids and nanofluids during flow boiling in a single microchannel at different heat fluxes and mass fluxes

The dissipation of high heat flux from integrated circuit chips and the maintenance of acceptable junction temperatures in high powered electronics require advanced cooling technologies. One such technology is two-phase cooling in microchannels under confined flow boiling conditions. In macroscale flow boiling bubbles will nucleate on the channel walls, grow, and depart from the surface. In microscale flow boiling bubbles can fill the channel diameter before the liquid drag force has a chance to sweep them off the channel wall. As a confined bubble elongates in a microchannel, it traps thin liquid films between the heated wall and the vapor core that are subject to large temperature gradients. The thin films evaporate rapidly, sometimes faster than the incoming mass flux can replenish bulk fluid in the microchannel. When the local vapor pressure spike exceeds the inlet pressure, it forces the upstream interface to travel back into the inlet plenum and create flow boiling instabilities. Flow boiling instabilities reduce the temperature at which critical heat flux occurs and create channel dryout. Dryout causes high surface temperatures that can destroy the electronic circuits that use two-phase micro heat exchangers for cooling. Flow boiling instability is characterized by periodic oscillation of flow regimes which induce oscillations in fluid temperature, wall temperatures, pressure drop, and mass flux. When nanofluids are used in flow boiling, the nanoparticles become deposited on the heated surface and change its thermal conductivity, roughness, capillarity, wettability, and nucleation site density. It also affects heat transfer by changing bubble departure diameter, bubble departure frequency, and the evaporation of the micro and macrolayer beneath the growing bubbles. Flow boiling was investigated in this study using degassed, deionized water, and 0.001 vol% aluminum oxide nanofluids in a single rectangular brass microchannel with a hydraulic diameter of 229 µm for one inlet fluid temperature of 63°C and two constant flow rates of 0.41 ml/min and 0.82 ml/min. The power input was adjusted for two average surface temperatures of 103°C and 119°C at each flow rate. High speed images were taken periodically for water and nanofluid flow boiling after durations of 25, 75, and 125 minutes from the start of flow. The change in regime timing revealed the effect of nanoparticle suspension and deposition on the Onset of Nucelate Boiling (ONB) and the Onset of Bubble Elongation (OBE). Cycle duration and bubble frequencies are reported for different nanofluid flow boiling durations. The addition of nanoparticles was found to stabilize bubble nucleation and growth and limit the recession rate of the upstream and downstream interfaces, mitigating the spreading of dry spots and elongating the thin film regions to increase thin film evaporation.

Thickness measurement of dynamic thin liquid films generated by plug-annular flow in non-wetting microchannels

Surface tension forces are significant at millimeter length-scales, causing profoundly different flow morphologies in microchannels than in macroscale flows. The existence and morphology of thin liquid films is particularly relevant for predicting performance and operational stability of devices containing microscale two phase flows. Analytical, computational, and experimental methods previously employed in the study of thin liquid films are discussed. Thicknesses before and after a novel film morphology, referred to as a `shock,' are measured with a novel film thickness measurement technique that uses confocal microscopy. Film thicknesses predicted by previous work are compared to experimental results. Methods for increasing the accuracy of the confocal film thickness measurement technique are discussed.

Two-Phase Flow In Microchannels: Morphology And Interface Phenomena

The existence and morphology, as well as the dynamics of micro-scale gas-liquid interfaces is investigated numerically and experimentally. These studies can be used to assess liquid management issues in microsystems such as PEMFC gas flow channels, and are meant to open new research perspectives in two-phase flow, particularly in film deposition on non-wetting surfaces. For example the critical plug volume data can be used to deliver desired length plugs, or to determine the plug formation frequency. The dynamics of gas-liquid interfaces, of interest for applications involving small passages (e.g. heat exchangers, phase separators and filtration systems), was investigated using high-speed microscopy - a method that also proved useful for the study of film deposition processes. The existence limit for a liquid plug forming in a mixed wetting channel is determined by numerical simulations using Surface Evolver. The plug model simulate actual conditions in the gas flow channels of PEM fuel cells, the wetting of the gas diffusion layer (GDL) side of the channel being different from the wetting of the bipolar plate walls. The minimum plug volume, denoted as critical volume is computed for a series of GDL and bipolar plate wetting properties. Critical volume data is meant to assist in the water management of PEMFC, when corroborated with experimental data. The effect of cross section geometry is assessed by computing the critical volume in square and trapezoidal channels. Droplet simulations show that water can be passively removed from the GDL surface towards the bipolar plate if we take advantage on differing wetting properties between the two surfaces, to possibly avoid the gas transport blockage through the GDL. High speed microscopy was employed in two-phase and film deposition experiments with water in round and square capillary tubes. Periodic interface destabilization was observed and the existence of compression waves in the gas phase is discussed by taking into consideration a naturally occurring convergent-divergent nozzle formed by the flowing liquid phase. The effect of channel geometry and wetting properties was investigated through two-phase water-air flow in square and round microchannels, having three static contact angles of 20, 80 and 105 degrees. Four different flow regimes are observed for a fixed flow rate, this being thought to be caused by the wetting behavior of liquid flowing in the corners as well as the liquid film stability. Film deposition experiments in wetting and non-wetting round microchannels show that a thicker film is deposited for wetting conditions departing from the ideal 0 degrees contact angle. A film thickness dependence with the contact angle theta as well as the Capillary number, in the form h_R ~ Ca^(2/3)/ cos(theta) is inferred from scaling arguments, for contact angles smaller than 36 degrees. Non-wetting film deposition experiments reveal that a film significantly thicker than the wetting Bretherton film is deposited. A hydraulic jump occurs if critical conditions are met, as given by a proposed nondimensional parameter similar to the Froude number. Film thickness correlations are also found by matching the measured and the proposed velocity derived in the shock theory. The surface wetting as well as the presence of the shock cause morphological changes in the Taylor bubble flow.

Analytical form of the Probability Density Function of travel times in rectangular zone with Tchebyshev’s metric

Geometrical dependencies are being researched for analytical representation of the probability density function (pdf) for the travel time between a random, and a known or another random point in Tchebyshev’s metric. In the most popular case - a rectangular area of service - the pdf of this random variable depends directly on the position of the server. Two approaches have been introduced for the exact analytical calculation of the pdf: Ad-hoc approach – useful for a ‘manual’ solving of a specific case; by superposition – an algorithmic approach for the general case. The main concept of each approach is explained, and a short comparison is done to prove the faithfulness.