Gas-solid two-phase flow in the riser of circulating fluidized beds: mathematical modelling and numerical simulation

A mathematical model is developed for gas-solids flows in circulating fluidized beds. An Eulerian formulation is followed based on the two-fluids model approach where both the fluid and the particulate phases are treated as a continuum. The physical modelling is discussed, including the formulation of boundary conditions and the description of the numerical methodology. Results of numerical simulation are presented and discussed. The model is validated through comparison to experiment, and simulation is performed to investigate the effects on the flow hydrodynamics of the solids viscosity.

Chaos and Unpredictability in the Vibration of an Elasto-Plastic Beam

This contribution discusses the nonlinear dynamics of a pin-ended elasto-plastic beam with both kinematic and isotropic hardening. An iterative numerical procedure based on the operator split technique is developed in order to deal with the nonlinearities in the equations of motion. Free and forced responses for harmonic sinusoidal and square wave excitations are investigated. Numerical simulations present many interesting behaviors such as jump phenomena, sensitivity to initial conditions, chaos and transient chaos. These results indicate that there are practical problems in predicting the response of the beam even when periodic steady state response is expected.

Covalidation of Integral Transforms and Method of Lines in Nonlinear Convection-Diffusion with Mathematica

The Mathematica system (version 4.0) is employed in the solution of nonlinear difusion and convection-difusion problems, formulated as transient one-dimensional partial diferential equations with potential dependent equation coefficients. The Generalized Integral Transform Technique (GITT) is first implemented for the hybrid numerical-analytical solution of such classes of problems, through the symbolic integral transformation and elimination of the space variable, followed by the utilization of the built-in Mathematica function NDSolve for handling the resulting transformed ODE system. This approach ofers an error-controlled final numerical solution, through the simultaneous control of local errors in this reliable ODE's solver and of the proposed eigenfunction expansion truncation order. For covalidation purposes, the same built-in function NDSolve is employed in the direct solution of these partial diferential equations, as made possible by the algorithms implemented in Mathematica (versions 3.0 and up), based on application of the method of lines. Various numerical experiments are performed and relative merits of each approach are critically pointed out.

Numerical Simulation of the Cavitation in the Hydrodynamic Lubrication of Journal Bearings: a Parallel Algorithm

In this paper we present an algorithm for the numerical simulation of the cavitation in the hydrodynamic lubrication of journal bearings. Despite the fact that this physical process is usually modelled as a free boundary problem, we adopted the equivalent variational inequality formulation. We propose a two-level iterative algorithm, where the outer iteration is associated to the penalty method, used to transform the variational inequality into a variational equation, and the inner iteration is associated to the conjugate gradient method, used to solve the linear system generated by applying the finite element method to the variational equation. This inner part was implemented using the element by element strategy, which is easily parallelized. We analyse the behavior of two physical parameters and discuss some numerical results. Also, we analyse some results related to the performance of a parallel implementation of the algorithm.

A Well Stated Time Domain Integral Representation for Elastodynamic Analysis and Applications

This article discusses three possible ways to derive time domain boundary integral representations for elastodynamics. This discussion points out possible difficulties found when using those formulations to deal with practical applications. The discussion points out recommendations to select the convenient integral representation to deal with elastodynamic problems and opens the possibility of deriving simplified schemes. The proper way to take into account initial conditions applied to the body is an interesting topict shown. It illustrates the main differences between the discussed boundary integral representation expressions, their singularities and possible numerical problems. The correct way to use collocation points outside the analyzed domain is carefully described. Some applications are shown at the end of the paper, in order to demonstrate the capabilities of the technique when properly used.

An Inspection System to Identify Fatigue Damage on Steel-Bridge Structures

Crack formation and growth in steel bridge structural elements may be due to loading oscillations. The welded elements are liable to internal discontinuities along welded joints and sensible to stress variations. The evaluation of the remaining life of a bridge is needed to make cost-effective decisions regarding inspection, repair, rehabilitation, and replacement. A steel beam model has been proposed to simulate crack openings due to cyclic loads. Two possible alternatives have been considered to model crack propagation, which the initial phase is based on the linear fracture mechanics. Then, the model is extended to take into account the elastoplastic fracture mechanic concepts. The natural frequency changes are directly related to moment of inertia variation and consequently to a reduction in the flexural stiffness of a steel beam. Thus, it is possible to adopt a nondestructive technique during steel bridge inspection to quantify the structure eigenvalue variation that will be used to localize the grown fracture. A damage detection algorithm is developed for the proposed model and the numerical results are compared with the solutions achieved by using another well know computer code.

Numerical study of wedge supported oblique shock wave-oblique detonation wave transitions

The results of a numerical study of premixed Hydrogen-air flows ignition by an oblique shock wave (OSW) stabilized by a wedge are presented, in situations when initial and boundary conditions are such that transition between the initial OSW and an oblique detonation wave (ODW) is observed. More precisely, the objectives of the paper are: (i) to identify the different possible structures of the transition region that exist between the initial OSW and the resulting ODW and (ii) to evidence the effect on the ODW of an abrupt decrease of the wedge angle in such a way that the final part of the wedge surface becomes parallel to the initial flow. For such a geometrical configuration and for the initial and boundary conditions considered, the overdriven detonation supported by the initial wedge angle is found to relax towards a Chapman-Jouguet detonation in the region where the wedge surface is parallel to the initial flow. Computations are performed using an adaptive, unstructured grid, finite volume computer code previously developed for the sake of the computations of high speed, compressible flows of reactive gas mixtures. Physico-chemical properties are functions of the local mixture composition, temperature and pressure, and they are computed using the CHEMKIN-II subroutines.

On state and parameter estimation in chaotic systems

State-of-the-art predictions of atmospheric states rely on large-scale numerical models of chaotic systems. This dissertation studies numerical methods for state and parameter estimation in such systems. The motivation comes from weather and climate models and a methodological perspective is adopted. The dissertation comprises three sections: state estimation, parameter estimation and chemical data assimilation with real atmospheric satellite data. In the state estimation part of this dissertation, a new filtering technique based on a combination of ensemble and variational Kalman filtering approaches, is presented, experimented and discussed. This new filter is developed for large-scale Kalman filtering applications. In the parameter estimation part, three different techniques for parameter estimation in chaotic systems are considered. The methods are studied using the parameterized Lorenz 95 system, which is a benchmark model for data assimilation. In addition, a dilemma related to the uniqueness of weather and climate model closure parameters is discussed. In the data-oriented part of this dissertation, data from the Global Ozone Monitoring by Occultation of Stars (GOMOS) satellite instrument are considered and an alternative algorithm to retrieve atmospheric parameters from the measurements is presented. The validation study presents first global comparisons between two unique satellite-borne datasets of vertical profiles of nitrogen trioxide (NO3), retrieved using GOMOS and Stratospheric Aerosol and Gas Experiment III (SAGE III) satellite instruments. The GOMOS NO3 observations are also considered in a chemical state estimation study in order to retrieve stratospheric temperature profiles. The main result of this dissertation is the consideration of likelihood calculations via Kalman filtering outputs. The concept has previously been used together with stochastic differential equations and in time series analysis. In this work, the concept is applied to chaotic dynamical systems and used together with Markov chain Monte Carlo (MCMC) methods for statistical analysis. In particular, this methodology is advocated for use in numerical weather prediction (NWP) and climate model applications. In addition, the concept is shown to be useful in estimating the filter-specific parameters related, e.g., to model error covariance matrix parameters.

Homogeneous assay technologies in drug screening: Quenching Resonance Energy Transfer (QRET) technique

Information gained from the human genome project and improvements in compound synthesizing have increased the number of both therapeutic targets and potential lead compounds. This has evolved a need for better screening techniques to have a capacity to screen number of compound libraries against increasing amount of targets. Radioactivity based assays have been traditionally used in drug screening but the fluorescence based assays have become more popular in high throughput screening (HTS) as they avoid safety and waste problems confronted with radioactivity. In comparison to conventional fluorescence more sensitive detection is obtained with time-resolved luminescence which has increased the popularity of time-resolved fluorescence resonance energy transfer (TR-FRET) based assays. To simplify the current TR-FRET based assay concept the luminometric homogeneous single-label utilizing assay technique, Quenching Resonance Energy Transfer (QRET), was developed. The technique utilizes soluble quencher to quench non-specifically the signal of unbound fraction of lanthanide labeled ligand. One labeling procedure and fewer manipulation steps in the assay concept are saving resources. The QRET technique is suitable for both biochemical and cell-based assays as indicated in four studies:1) ligand screening study of β2 -adrenergic receptor (cell-based), 2) activation study of Gs-/Gi-protein coupled receptors by measuring intracellular concentration of cyclic adenosine monophosphate (cell-based), 3) activation study of G-protein coupled receptors by observing the binding of guanosine-5’-triphosphate (cell membranes), and 4) activation study of small GTP binding protein Ras (biochemical). Signal-to-background ratios were between 2.4 to 10 and coefficient of variation varied from 0.5 to 17% indicating their suitability to HTS use.

Numerical Simulation of Stochastic Di erential Equations

Stochastic differential equation (SDE) is a differential equation in which some of the terms and its solution are stochastic processes. SDEs play a central role in modeling physical systems like finance, Biology, Engineering, to mention some. In modeling process, the computation of the trajectories (sample paths) of solutions to SDEs is very important. However, the exact solution to a SDE is generally difficult to obtain due to non-differentiability character of realizations of the Brownian motion. There exist approximation methods of solutions of SDE. The solutions will be continuous stochastic processes that represent diffusive dynamics, a common modeling assumption for financial, Biology, physical, environmental systems. This Masters' thesis is an introduction and survey of numerical solution methods for stochastic differential equations. Standard numerical methods, local linearization methods and filtering methods are well described. We compute the root mean square errors for each method from which we propose a better numerical scheme. Stochastic differential equations can be formulated from a given ordinary differential equations. In this thesis, we describe two kind of formulations: parametric and non-parametric techniques. The formulation is based on epidemiological SEIR model. This methods have a tendency of increasing parameters in the constructed SDEs, hence, it requires more data. We compare the two techniques numerically.

Genetic diversity among proso millet (Panicum miliaceum) biotypes assessed by AFLP technique

The Amplified Fragment Length Polymorphism (AFLP) technique was used to access genetic diversity between three domestic and nine wild proso millet biotypes from the United States and Canada. Eight primer combinations detected 39 polymorphic DNA fragments, with the genetic distance estimates among biotypes ranging from 0.02 to 0.04. Colorado-Weld County black seeded and Wyoming-Platte County were the most distinct biotypes according to the dissimilarity level. A UPGMA cluster analysis revealed two distinct groups of proso millet without any geographic association. Six weed biotypes exhibiting some characters of cultivated plants were grouped together with domesticated biotypes of proso millet while the three typical wild phenotypes were clearly clustered into another group according to AFLP markers.

Eletrolite leakage as a technique to diagnose euphorbia heterophylla biotypes resistant to ppo-inhibitors herbicides

The action of herbicides that affect the integrity of cell membranes and cause leakage, like PPO-inhibitors, can be detected by measuring the electric conductivity (EC) of a solution in which the plant tissue target is incubated in the presence of herbicide. The objectives of this work were to confirm PPO resistance in a new Euphorbia heterophylla (EPHHL) biotype, and to compare the electrolyte leakage from R and S to PPO-inhibitors biotypes, using two different methods of incubation in a solution containing herbicides. One experiment was carried in greenhouse and three in laboratory, with a completely randomized design. In the greenhouse experiment, four biotypes of EPHHL were sprayed with seven rates of fomesafen to confirm resistance in suspected biotypes. Leaf disks from R and S EPHHL biotypes in the second and the third experiments and entire leaves in the fourth experiment were incubated in a solution containing PPO-inhibitors to subsequently determine EC of solution. The study confirmed the resistance to PPO-inhibitors in two EPHHL biotypes. There were no significant differences between S and R biotypes in the experiments with the incubation of leaf disks, but incubation of entire leaves of EPHHL S biotype showed higher EC when in a solution with fomesafen, in comparison to the R biotype. The results of this work are an indirect evidence that resistance to PPO-inhibitors is related to lower absorption of herbicide by the shoots and also to some kind of mechanism to cope with oxidative stress.

Optimization and validation of the solid-liquid extraction technique for determination of picloram in soils by high performance liquid chromatography

The objective of this study was to optimize and validate the solid-liquid extraction (ESL) technique for determination of picloram residues in soil samples. At the optimization stage, the optimal conditions for extraction of soil samples were determined using univariate analysis. Ratio soil/solution extraction, type and time of agitation, ionic strength and pH of extraction solution were evaluated. Based on the optimized parameters, the following method of extraction and analysis of picloram was developed: weigh 2.00 g of soil dried and sieved through a sieve mesh of 2.0 mm pore, add 20.0 mL of KCl concentration of 0.5 mol L-1, shake the bottle in the vortex for 10 seconds to form suspension and adjust to pH 7.00, with alkaline KOH 0.1 mol L-1. Homogenate the system in a shaker system for 60 minutes and then let it stand for 10 minutes. The bottles are centrifuged for 10 minutes at 3,500 rpm. After the settlement of the soil particles and cleaning of the supernatant extract, an aliquot is withdrawn and analyzed by high performance liquid chromatography. The optimized method was validated by determining the selectivity, linearity, detection and quantification limits, precision and accuracy. The ESL methodology was efficient for analysis of residues of the pesticides studied, with percentages of recovery above 90%. The limits of detection and quantification were 20.0 and 66.0 mg kg-1 soil for the PVA, and 40.0 and 132.0 mg kg-1 soil for the VLA. The coefficients of variation (CV) were equal to 2.32 and 2.69 for PVA and TH soils, respectively. The methodology resulted in low organic solvent consumption and cleaner extracts, as well as no purification steps for chromatographic analysis were required. The parameters evaluated in the validation process indicated that the ESL methodology is efficient for the extraction of picloram residues in soils, with low limits of detection and quantification.

The Combined treatment of UV Light and Titanium Dioxide (TiO2) in Grease Filtration Technique

The aim of this thesis was to identify the best grease removal technique with the application of low power of UV light to TiO2 coated grease filters. The treatment with various power series of ozone generating and ozone free lamps to normal grease filters and TiO2 coated grease filters were examined and the obtained results are compared to each other in this paper. The effect of ozone reaction was observed and compared with the effect of TiO2. The experiments were solely based on the photo oxidation and photo catalytic oxidation reactions. TiO2 is a green catalyst used in the photocatalytic reaction. Sunflower oil was used for grease production and tetracholoroethylene as a solvent. Grease samples were collected from the ventilation duct connected to the cooking hood system. Sample extraction was done in ultrasonic bath with the principle of sonication. The sample analysis was done by FTIR machine. The result determining the concentration of grease was the quantification of saturated C-H bonds in the chosen peak group of the spectrum. A very low power of UVC light functions perfectly with the Titanium dioxide. The experimental results have shown the combined treatment of titanium dioxide and UV light is an effective method in grease removal process. The photocatalytic reaction with titanium dioxide is better than photo oxidation reaction with ozone treatment. Photocatalytic reaction is environmentally friendly, energy efficient and economical.

Transient current technique simulation of semiconductor detectors

Nowadays advanced simulation technologies of semiconductor devices occupies an important place in microelectronics production process. Simulation helps to understand devices internal processes physics, detect new effects and find directions for optimization. Computer calculation reduces manufacturing costs and time. Modern simulation suits such as Silcaco TCAD allow simulating not only individual semiconductor structures, but also these structures in the circuit. For that purpose TCAD include MixedMode tool. That tool can simulate circuits using compact circuit models including semiconductor structures with their physical models. In this work, MixedMode is used for simulating transient current technique setup, which include detector and supporting electrical circuit. This technique was developed by RD39 collaboration project for investigation radiation detectors radiation hard properties.