Computational modeling of optical trapping of vaterite microspheres

The ability to grow microscopic spherical birefringent crystals of vaterite, a calcium carbonate mineral, has allowed the development of an optical microrheometer based on optical tweezers. However, since these crystals are birefringent, and worse, are expected to have non-uniform birefringence, computational modeling of the microrheometer is a highly challenging task. Modeling the microrheometer - and optical tweezers in general - typically requires large numbers of repeated calculations for the same trapped particle. This places strong demands on the efficiency of computational methods used. While our usual method of choice for computational modelling of optical tweezers - the T-matrix method - meets this requirement of efficiency, it is restricted to homogeneous isotropic particles. General methods that can model complex structures such as the vaterite particles, such as finite-difference time-domain (FDTD) or finite-difference frequency-domain (FDFD) methods, are inefficient. Therefore, we have developed a hybrid FDFD/T-matrix method that combines the generality of volume-discretisation methods such as FDFD with the efficiency of the T-matrix method. We have used this hybrid method to calculate optical forces and torques on model vaterite spheres in optical traps. We present and compare the results of computational modelling and experimental measurements.

The programming language python in earth system simulations

-scale vary from a planetary scale and million years for convection problems to 100km and 10 years for fault systems simulations. Various techniques are in use to deal with the time dependency (e.g. Crank-Nicholson), with the non-linearity (e.g. Newton-Raphson) and weakly coupled equations (e.g. non-linear Gauss-Seidel). Besides these high-level solution algorithms discretization methods (e.g. finite element method (FEM), boundary element method (BEM)) are used to deal with spatial derivatives. Typically, large-scale, three dimensional meshes are required to resolve geometrical complexity (e.g. in the case of fault systems) or features in the solution (e.g. in mantel convection simulations). The modelling environment escript allows the rapid implementation of new physics as required for the development of simulation codes in earth sciences. Its main object is to provide a programming language, where the user can define new models and rapidly develop high-level solution algorithms. The current implementation is linked with the finite element package finley as a PDE solver. However, the design is open and other discretization technologies such as finite differences and boundary element methods could be included. escript is implemented as an extension of the interactive programming environment python (see www.python.org). Key concepts introduced are Data objects, which are holding values on nodes or elements of the finite element mesh, and linearPDE objects, which are defining linear partial differential equations to be solved by the underlying discretization technology. In this paper we will show the basic concepts of escript and will show how escript is used to implement a simulation code for interacting fault systems. We will show some results of large-scale, parallel simulations on an SGI Altix system. Acknowledgements: Project work is supported by Australian Commonwealth Government through the Australian Computational Earth Systems Simulator Major National Research Facility, Queensland State Government Smart State Research Facility Fund, The University of Queensland and SGI.

The fallacy of integration:work and non-work in professional services

Many organisations are encouraging their staff to integrate work and non-work, but a qualitative study of young professionals found that many crave greater segregation rather than more integration. Most wished to build boundaries to separate the two and simplify a complex world. Where working practices render traditional boundaries of time and space ineffective, this population seems to create new idiosyncratic boundaries to segregate work from non-work. These idiosyncratic boundaries depended on age, culture and life-stage though for most of this population there was no appreciable gender difference in attitudes to segregating work and non-work. Gender differences only became noticeable for parents. A matrix defining the dimensions to these boundaries is proposed that may advance understanding of how individuals separate their work and personal lives. In turn, this may facilitate the development of policies and practices to integrate work and non-work that meet individual as well as organisational needs.

Molecular phase space transport in water:non-stationary random walk model

Molecular transport in phase space is crucial for chemical reactions because it defines how pre-reactive molecular configurations are found during the time evolution of the system. Using Molecular Dynamics (MD) simulated atomistic trajectories we test the assumption of the normal diffusion in the phase space for bulk water at ambient conditions by checking the equivalence of the transport to the random walk model. Contrary to common expectations we have found that some statistical features of the transport in the phase space differ from those of the normal diffusion models. This implies a non-random character of the path search process by the reacting complexes in water solutions. Our further numerical experiments show that a significant long period of non-stationarity in the transition probabilities of the segments of molecular trajectories can account for the observed non-uniform filling of the phase space. Surprisingly, the characteristic periods in the model non-stationarity constitute hundreds of nanoseconds, that is much longer time scales compared to typical lifetime of known liquid water molecular structures (several picoseconds).

Low Volatility Options and Numerical Diffusion of Finite Difference Schemes

2000 Mathematics Subject Classification: 65M06, 65M12.

Efficiency improvement of non-uniformly-aged PV arrays

The utilization of solar energy by photovoltaic (PV) systems have received much research and development (R&D) attention across the globe. In the past decades, a large number of PV array have been installed. Since the installed PV arrays often operate in harsh environments, non-uniform aging can occur and impact adversely on the performance of PV systems, especially in the middle and late periods of their service life. Due to the high cost of replacing aged PV modules by new modules, it is appealing to improve energy efficiency of aged PV systems. For this purpose, this paper presents a PV module reconfiguration strategy to achieve the maximum power generation from non-uniformly aged PV arrays without significant investment. The proposed reconfiguration strategy is based on the cell-unit structure of PV modules, the operating voltage limit of gird-connected converter, and the resulted bucket-effect of the maximum short circuit current. The objectives are to analyze all the potential reorganization options of the PV modules, find the maximum power point and express it in a proposition. This proposition is further developed into a novel implementable algorithm to calculate the maximum power generation and the corresponding reconfiguration of the PV modules. The immediate benefits from this reconfiguration are the increased total power output and maximum power point voltage information for global maximum power point tracking (MPPT). A PV array simulation model is used to illustrate the proposed method under three different cases. Furthermore, an experimental rig is built to verify the effectiveness of the proposed method. The proposed method will open an effective approach for condition-based maintenance of emerging aging PV arrays.

Advanced thermal analysis of microelectronics using spreading resistance models

Thermal analysis of electronic devices is one of the most important steps for designing of modern devices. Precise thermal analysis is essential for designing an effective thermal management system of modern electronic devices such as batteries, LEDs, microelectronics, ICs, circuit boards, semiconductors and heat spreaders. For having a precise thermal analysis, the temperature profile and thermal spreading resistance of the device should be calculated by considering the geometry, property and boundary conditions. Thermal spreading resistance occurs when heat enters through a portion of a surface and flows by conduction. It is the primary source of thermal resistance when heat flows from a tiny heat source to a thin and wide heat spreader. In this thesis, analytical models for modeling the temperature behavior and thermal resistance in some common geometries of microelectronic devices such as heat channels and heat tubes are investigated. Different boundary conditions for the system are considered. Along the source plane, a combination of discretely specified heat flux, specified temperatures and adiabatic condition are studied. Along the walls of the system, adiabatic or convective cooling boundary conditions are assumed. Along the sink plane, convective cooling with constant or variable heat transfer coefficient are considered. Also, the effect of orthotropic properties is discussed. This thesis contains nine chapters. Chapter one is the introduction and shows the concepts of thermal spreading resistance besides the originality and importance of the work. Chapter two reviews the literatures on the thermal spreading resistance in the past fifty years with a focus on the recent advances. In chapters three and four, thermal resistance of a twodimensional flux channel with non-uniform convection coefficient in the heat sink plane is studied. The non-uniform convection is modeled by using two functions than can simulate a wide variety of different heat sink configurations. In chapter five, a non-symmetrical flux channel with different heat transfer coefficient along the right and left edges and sink plane is analytically modeled. Due to the edge cooling and non-symmetry, the eigenvalues of the system are defined using the heat transfer coefficient on both edges and for satisfying the orthogonality condition, a normalized function is calculated. In chapter six, thermal behavior of two-dimensional rectangular flux channel with arbitrary boundary conditions on the source plane is presented. The boundary condition along the source plane can be a combination of the first kind boundary condition (Dirichlet or prescribed temperature) and the second kind boundary condition (Neumann or prescribed heat flux). The proposed solution can be used for modeling the flux channels with numerous different source plane boundary conditions without any limitations in the number and position of heat sources. In chapter seven, temperature profile of a circular flux tube with discretely specified boundary conditions along the source plane is presented. Also, the effect of orthotropic properties are discussed. In chapter 8, a three-dimensional rectangular flux channel with a non-uniform heat convection along the heat sink plane is analytically modeled. In chapter nine, a summary of the achievements is presented and some systems are proposed for the future studies. It is worth mentioning that all the models and case studies in the thesis are compared with the Finite Element Method (FEM).

Principal Component Analysis of Results Obtained from Finite-Difference Time-Domain Algorithms

Finite-Differences Time-Domain (FDTD) algorithms are well established tools of computational electromagnetism. Because of their practical implementation as computer codes, they are affected by many numerical artefact and noise. In order to obtain better results we propose using Principal Component Analysis (PCA) based on multivariate statistical techniques. The PCA has been successfully used for the analysis of noise and spatial temporal structure in a sequence of images. It allows a straightforward discrimination between the numerical noise and the actual electromagnetic variables, and the quantitative estimation of their respective contributions. Besides, The GDTD results can be filtered to clean the effect of the noise. In this contribution we will show how the method can be applied to several FDTD simulations: the propagation of a pulse in vacuum, the analysis of two-dimensional photonic crystals. In this last case, PCA has revealed hidden electromagnetic structures related to actual modes of the photonic crystal.

Effect of web holes on web crippling strength of cold-formed steel channel sections under end-one-flange loading condition – Part I: Tests and finite element analysis

Web openings could be used in cold-formed steel beam members, such as wall studs or floor joints, to facilitate ease of services in buildings. In this paper, a combination of tests and non-linear finite element analyses is used to investigate the effect of such holes on web crippling under end-one-flange (EOF) loading condition; the cases of both flanges fastened and unfastened to the bearing plates are considered. The results of 74 web crippling tests are presented, with 22 tests conducted on channel sections without web openings and 52 tests conducted on channel sections with web openings. In the case of the tests with web openings, the hole was either located centred above the bearing plates or having a horizontal clear distance to the near edge of the bearing plates. A good agreement between the tests and finite element analyses was obtained in term of both strength and failure modes.

Influência de feixes de condutores em linhas de transmissão UAT CA

Dissertação (mestrado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Elétrica, 2015.

Finite difference calculations of permeability in large domains in a wide porosity range.

Determining effective hydraulic, thermal, mechanical and electrical properties of porous materials by means of classical physical experiments is often time-consuming and expensive. Thus, accurate numerical calculations of material properties are of increasing interest in geophysical, manufacturing, bio-mechanical and environmental applications, among other fields. Characteristic material properties (e.g. intrinsic permeability, thermal conductivity and elastic moduli) depend on morphological details on the porescale such as shape and size of pores and pore throats or cracks. To obtain reliable predictions of these properties it is necessary to perform numerical analyses of sufficiently large unit cells. Such representative volume elements require optimized numerical simulation techniques. Current state-of-the-art simulation tools to calculate effective permeabilities of porous materials are based on various methods, e.g. lattice Boltzmann, finite volumes or explicit jump Stokes methods. All approaches still have limitations in the maximum size of the simulation domain. In response to these deficits of the well-established methods we propose an efficient and reliable numerical method which allows to calculate intrinsic permeabilities directly from voxel-based data obtained from 3D imaging techniques like X-ray microtomography. We present a modelling framework based on a parallel finite differences solver, allowing the calculation of large domains with relative low computing requirements (i.e. desktop computers). The presented method is validated in a diverse selection of materials, obtaining accurate results for a large range of porosities, wider than the ranges previously reported. Ongoing work includes the estimation of other effective properties of porous media.

numerical problems in the calculation of electric field and charge density profiles for mass impregnated non-draining hvdc cables subjected to fast polarity reversals

The goal of this simulation thesis is to present a tool for studying and eliminating various numerical problems observed while analyzing the behavior of the MIND cable during fast voltage polarity reversal. The tool is built on the MATLAB environment, where several simulations were run to achieve oscillation-free results. This thesis will add to earlier research on HVDC cables subjected to polarity reversals. Initially, the code does numerical simulations to analyze the electric field and charge density behavior of a MIND cable for certain scenarios such as before, during, and after polarity reversal. However, the primary goal is to reduce numerical oscillations from the charge density profile. The generated code is notable for its usage of the Arithmetic Mean Approach and the Non-Uniform Field Approach for filtering and minimizing oscillations even under time and temperature variations.

Bile And Liver Metallothionein Behavior In Copper-exposed Fish.

The present study analyzed metallothionein (MT) excretion from liver to bile in Nile Tilapia (Oreochromis niloticus) exposed to sub-lethal copper concentrations (2mgL(-1)) in a laboratory setting. MTs in liver and bile were quantified by spectrophotometry after thermal incubation and MT metal-binding profiles were characterized by size exclusion high performance liquid chromatography coupled to ICP-MS (SEC-HPLC-ICP-MS). Results show that liver MT is present in approximately 250-fold higher concentrations than bile MT in non-exposed fish. Differences between the MT profiles from the control and exposed group were observed for both matrices, indicating differential metal-binding behavior when comparing liver and bile MT. This is novel data regarding intra-organ MT comparisons, since differences between organs are usually present only with regard to quantification, not metal-binding behavior. Bile MT showed statistically significant differences between the control and exposed group, while the same did not occur with liver MT. This indicates that MTs synthesized in the liver accumulate more slowly than MTs excreted from liver to bile, since the same fish presented significantly higher MT levels in liver when compared to bile. We postulate that bile, although excreted in the intestine and partially reabsorbed by the same returning to the liver, may also release MT-bound metals more rapidly and efficiently, which may indicate an efficient detoxification route. Thus, we propose that the analysis of bile MTs to observe recent metal exposure may be more adequate than the analysis of liver MTs, since organism responses to metals are more quickly observed in bile, although further studies are necessary.

Modelagem matemática e simulação numérica do resfriamento rápido de morango com ar forçado

This work approaches the forced air cooling of strawberry by numerical simulation. The mathematical model that was used describes the process of heat transfer, based on the Fourier's law, in spherical coordinates and simplified to describe the one-dimensional process. For the resolution of the equation expressed for the mathematical model, an algorithm was developed based on the explicit scheme of the numerical method of the finite differences and implemented in the scientific computation program MATLAB 6.1. The validation of the mathematical model was made by the comparison between theoretical and experimental data, where strawberries had been cooled with forced air. The results showed to be possible the determination of the convective heat transfer coefficient by fitting the numerical and experimental data. The methodology of the numerical simulations was showed like a promising tool in the support of the decision to use or to develop equipment in the area of cooling process with forced air of spherical fruits.

A hybrid-Trefftz formulation for plane elasticity with selective enrichment of the approximations

This work is related to the so-called non-conventional finite element formulations. Essentially, a methodology for the enrichment of the initial approximation which is typical of the meshless methods and based on the clouds concept is introduced in the hybrid-Trefftz formulation for plane elasticity. The formulation presented allows for the approximation and direct enrichment of two independent fields: stresses in the domains and displacements on the boundaries of the elements. Defined by a set of elements and interior boundaries sharing a common node, the cloud notion is employed to select the enrichment support for the approximation fields. The numerical analysis performed reveals an excellent performance of the resulting formulation, characterized by the good approximation ability and a reduced computational effort. Copyright (C) 2009 John Wiley & Sons, Ltd.