Hybrid multi-objective optimization and hybridized self-organizing response surface method

Numerical optimization is a technique where a computer is used to explore design parameter combinations to find extremes in performance factors. In multi-objective optimization several performance factors can be optimized simultaneously. The solution to multi-objective optimization problems is not a single design, but a family of optimized designs referred to as the Pareto frontier. The Pareto frontier is a trade-off curve in the objective function space composed of solutions where performance in one objective function is traded for performance in others. A Multi-Objective Hybridized Optimizer (MOHO) was created for the purpose of solving multi-objective optimization problems by utilizing a set of constituent optimization algorithms. MOHO tracks the progress of the Pareto frontier approximation development and automatically switches amongst those constituent evolutionary optimization algorithms to speed the formation of an accurate Pareto frontier approximation. Aerodynamic shape optimization is one of the oldest applications of numerical optimization. MOHO was used to perform shape optimization on a 0.5-inch ballistic penetrator traveling at Mach number 2.5. Two objectives were simultaneously optimized: minimize aerodynamic drag and maximize penetrator volume. This problem was solved twice. The first time the problem was solved by using Modified Newton Impact Theory (MNIT) to determine the pressure drag on the penetrator. In the second solution, a Parabolized Navier-Stokes (PNS) solver that includes viscosity was used to evaluate the drag on the penetrator. The studies show the difference in the optimized penetrator shapes when viscosity is absent and present in the optimization. In modern optimization problems, objective function evaluations may require many hours on a computer cluster to perform these types of analysis. One solution is to create a response surface that models the behavior of the objective function. Once enough data about the behavior of the objective function has been collected, a response surface can be used to represent the actual objective function in the optimization process. The Hybrid Self-Organizing Response Surface Method (HYBSORSM) algorithm was developed and used to make response surfaces of objective functions. HYBSORSM was evaluated using a suite of 295 non-linear functions. These functions involve from 2 to 100 variables demonstrating robustness and accuracy of HYBSORSM.

Lattice Boltzmann method for flow and heat transfer in microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Method development for the analysis of smokeless powders and organic gunshot residue by ultra performance liquid chromatography with tandem mass spectrometry

The goal of this project was to develop a rapid separation and detection method for analyzing organic compounds in smokeless powders and then test its applicability on gunshot residue (GSR) samples. In this project, a total of 20 common smokeless powder additives and their decomposition products were separated by ultra performance liquid chromatography (UPLC) and confirmed by tandem mass spectrometry (MS/MS) using multiple reaction monitoring mode (MRM). Some of the targeted compounds included diphenylamines, centralites, nitrotoluenes, nitroglycerin, and various phthalates. The compounds were ionized in the MS source using simultaneous positive and negative electrospray ionization (ESI) with negative atmospheric pressure chemical ionization (APCI) in order to detect all compounds in a single analysis. The developed UPLC/MS/MS method was applied to commercially available smokeless powders and gunshot residue samples recovered from the hands of shooters, spent cartridges, and smokeless powder retrieved from unfired cartridges. Distinct compositions were identified for smokeless powders from different manufacturers and from separate manufacturing lots. The procedure also produced specific chemical profiles when tested on gunshot residues from different manufacturers. Overall, this thesis represents the development of a rapid and reproducible procedure capable of simultaneously detecting the widest possible range of components present in organic gunshot residue.^

Lattice Boltzmann Method for Flow and Heat Transfer in Microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Environmental analysis of polar herbicides in complex organic-rich matrices by high performance liquid chromatography atmospheric pressure ionization mass spectrometry (HPLC-API-MS)

A comprehensive forensic investigation of sensitive ecosystems in the Everglades Area is presented. Assessing the background levels of contamination in these ecosystems represents a vital resource to build up forensic evidence required to enforce future environmental crimes within the studied areas. This investigation presents the development and validation of a fractionation and isolation method for two families of herbicides commonly applied in the vicinity of the study area, including phenoxy acids like 2,4-D, MCPA, and silvex; as well as the most common triazine-based herbicides like atrazine, prometyne, simazine and related metabolites like DIA and DEA. Accelerated solvent extraction (ASE) and solid phase extraction (SPE) were used to isolate the analytes from abiotic matrices containing large amounts of organic material. Atmospheric-pressure ionization (API) with electrospray ionization in negative mode (ESP-), and Chemical Ionization in the positive mode (APCI+) were used to perform the characterization of the herbicides of interest.

Advanced control of managed pressure drilling

Automation of managed pressure drilling (MPD) enhances the safety and increases efficiency of drilling and that drives the development of controllers and observers for MPD. The objective is to maintain the bottom hole pressure (BHP) within the pressure window formed by the reservoir pressure and fracture pressure and also to reject kicks. Practical MPD automation solutions must address the nonlinearities and uncertainties caused by the variations in mud flow rate, choke opening, friction factor, mud density, etc. It is also desired that if pressure constraints are violated the controller must take appropriate actions to reject the ensuing kick. The objectives are addressed by developing two controllers: a gain switching robust controller and a nonlinear model predictive controller (NMPC). The robust gain switching controller is designed using H1 loop shaping technique, which was implemented using high gain bumpless transfer and 2D look up table. Six candidate controllers were designed in such a way they preserve robustness and performance for different choke openings and flow rates. It is demonstrated that uniform performance is maintained under different operating conditions and the controllers are able to reject kicks using pressure control and maintain BHP during drill pipe extension. The NMPC was designed to regulate the BHP and contain the outlet flow rate within certain tunable threshold. The important feature of that controller is that it can reject kicks without requiring any switching and thus there is no scope for shattering due to switching between pressure and flow control. That is achieved by exploiting the constraint handling capability of NMPC. Active set method was used for computing control inputs. It is demonstrated that NMPC is able to contain kicks and maintain BHP during drill pipe extension.

Polarimetric sensitivity to hydrostatic pressure and temperature in birefringent dual-core microstructured polymer fiber

We experimentally characterized a birefringent microstructured polymer fiber of specific construction, which allows for single mode propagation in two cores separated by a pair of large holes. The fiber exhibits high birefringence in each of the cores as well as relatively weak coupling between the cores. Spectral dependence of the group and the phase modal birefringence was measured using an interferometric method. We have also measured the sensing characteristics of the fiber such as the polarimetric sensitivity to hydrostatic pressure and temperature. © 2010 SPIE.

Ultra High Pressure Hydrogen Studies

Hydrogen has been called the fuel of the future, and as it’s non- renewable counterparts become scarce the economic viability of hydrogen gains traction. The potential of hydrogen is marked by its high mass specific energy density and wide applicability as a fuel in fuel cell vehicles and homes. However hydrogen’s volume must be reduced via pressurization or liquefaction in order to make it more transportable and volume efficient. Currently the vast majority of industrially produced hydrogen comes from steam reforming of natural gas. This practice yields low-pressure gas which must then be compressed at considerable cost and uses fossil fuels as a feedstock leaving behind harmful CO and CO2 gases as a by-product. The second method used by industry to produce hydrogen gas is low pressure electrolysis. In comparison the electrolysis of water at low pressure can produce pure hydrogen and oxygen gas with no harmful by-products using only water as a feedstock, but it will still need to be compressed before use. Multiple theoretical works agree that high pressure electrolysis could reduce the energy losses due to product gas compression. However these works openly admit that their projected gains are purely theoretical and ignore the practical limitations and resistances of a real life high pressure system. The goal of this work is to experimentally confirm the proposed thermodynamic gains of ultra-high pressure electrolysis in alkaline solution and characterize the behavior of a real life high pressure system.

A New Method for Modeling Free Surface Flows and Fluid-structure Interaction with Ocean Applications

The computational modeling of ocean waves and ocean-faring devices poses numerous challenges. Among these are the need to stably and accurately represent both the fluid-fluid interface between water and air as well as the fluid-structure interfaces arising between solid devices and one or more fluids. As techniques are developed to stably and accurately balance the interactions between fluid and structural solvers at these boundaries, a similarly pressing challenge is the development of algorithms that are massively scalable and capable of performing large-scale three-dimensional simulations on reasonable time scales. This dissertation introduces two separate methods for approaching this problem, with the first focusing on the development of sophisticated fluid-fluid interface representations and the second focusing primarily on scalability and extensibility to higher-order methods.

We begin by introducing the narrow-band gradient-augmented level set method (GALSM) for incompressible multiphase Navier-Stokes flow. This is the first use of the high-order GALSM for a fluid flow application, and its reliability and accuracy in modeling ocean environments is tested extensively. The method demonstrates numerous advantages over the traditional level set method, among these a heightened conservation of fluid volume and the representation of subgrid structures.

Next, we present a finite-volume algorithm for solving the incompressible Euler equations in two and three dimensions in the presence of a flow-driven free surface and a dynamic rigid body. In this development, the chief concerns are efficiency, scalability, and extensibility (to higher-order and truly conservative methods). These priorities informed a number of important choices: The air phase is substituted by a pressure boundary condition in order to greatly reduce the size of the computational domain, a cut-cell finite-volume approach is chosen in order to minimize fluid volume loss and open the door to higher-order methods, and adaptive mesh refinement (AMR) is employed to focus computational effort and make large-scale 3D simulations possible. This algorithm is shown to produce robust and accurate results that are well-suited for the study of ocean waves and the development of wave energy conversion (WEC) devices.

Nonlinear soil-structure interaction for buried pressure pipes under differential ground motion

Pipelines extend thousands of kilometers across wide geographic areas as a network to provide essential services for modern life. It is inevitable that pipelines must pass through unfavorable ground conditions, which are susceptible to natural disasters. This thesis investigates the behaviour of buried pressure pipelines experiencing ground distortions induced by normal faulting. A recent large database of physical modelling observations on buried pipes of different stiffness relative to the surrounding soil subjected to normal faults provided a unique opportunity to calibrate numerical tools. Three-dimensional finite element models were developed to enable the complex soil-structure interaction phenomena to be further understood, especially on the subjects of gap formation beneath the pipe and the trench effect associated with the interaction between backfill and native soils. Benchmarked numerical tools were then used to perform parametric analysis regarding project geometry, backfill material, relative pipe-soil stiffness and pipe diameter. Seismic loading produces a soil displacement profile that can be expressed by isoil, the distance between the peak curvature and the point of contraflexure. A simplified design framework based on this length scale (i.e., the Kappa method) was developed, which features estimates of longitudinal bending moments of buried pipes using a characteristic length, ipipe, the distance from peak to zero curvature. Recent studies indicated that empirical soil springs that were calibrated against rigid pipes are not suitable for analyzing flexible pipes, since they lead to excessive conservatism (for design). A large-scale split-box normal fault simulator was therefore assembled to produce experimental data for flexible PVC pipe responses to a normal fault. Digital image correlation (DIC) was employed to analyze the soil displacement field, and both optical fibres and conventional strain gauges were used to measure pipe strains. A refinement to the Kappa method was introduced to enable the calculation of axial strains as a function of pipe elongation induced by flexure and an approximation of the longitudinal ground deformations. A closed-form Winkler solution of flexural response was also derived to account for the distributed normal fault pattern. Finally, these two analytical solutions were evaluated against the pipe responses observed in the large-scale laboratory tests.

Simulation of Oil Shale Semi-Coke Particle Cold Transportation in a Spouted Bed Using CPFD Method

The spouted bed was widely used due to its good mixing of particles and effective phase transferability between the gas and solid phase. In this paper, the transportation process of particles in a 3D spouted bed was studied using the Computational Particle Fluid Dynamics (CPFD) numerical method. Experiments were conducted to verify the validity of the simulation results. Distributions of the pressure, velocities and particle concentration of transportation devices were investigated. The motion state and characteristics of multiphase flows in the transportation device were demonstrated under various operating conditions. The results showed that a good consistency was obtained between the simulated results and the experimental results. The motion characteristics of the gas-solid two-phase flow in the device was effectively predicted, which could assist the optimal operating condition estimation for the spouted transportation process.

The Use of a Pressure-Indicating Sensor Film to Provide Feedback upon Hydrogel-Forming Microneedle Array Self-Application In Vivo

PURPOSE:
To evaluate the combination of a pressure-indicating sensor film with hydrogel-forming microneedle arrays, as a method of feedback to confirm MN insertion in vivo.
METHODS:
Pilot in vitro insertion studies were conducted using a Texture Analyser to insert MN arrays, coupled with a pressure-indicating sensor film, at varying forces into excised neonatal porcine skin. In vivo studies involved twenty human volunteers, who self-applied two hydrogel-forming MN arrays, one with a pressure-indicating sensor film incorporated and one without. Optical coherence tomography was employed to measure the resulting penetration depth and colorimetric analysis to investigate the associated colour change of the pressure-indicating sensor film.
RESULTS:
Microneedle insertion was achieved in vitro at three different forces, demonstrating the colour change of the pressure-indicating sensor film upon application of increasing pressure. When self-applied in vivo, there was no significant difference in the microneedle penetration depth resulting from each type of array, with a mean depth of 237 μm recorded. When the pressure-indicating sensor film was present, a colour change occurred upon each application, providing evidence of insertion.
CONCLUSIONS:
For the first time, this study shows how the incorporation of a simple, low-cost pressure-indicating sensor film can indicate microneedle insertion in vitro and in vivo, providing visual feedback to assure the user of correct application. Such a strategy may enhance usability of a microneedle device and, hence, assist in the future translation of the technology to widespread clinical use.

High Pressure Speed of Sound and Related Thermodynamic Properties of 1-Alkyl-3-methylimidazolium Bis[(trifluoromethyl)sulfonyl]imides (from 1-Propyl- to 1-Hexyl-)

The knowledge of thermodynamic high-pressure speed of sound in ionic liquids (ILs) is a crucial way either to study the nature of the molecular interactions, structure and packing effects or to determine other key thermodynamic properties of ILs essential for their applications in any chemical and industrial processes. Herein, we report the speed of sound as a function temperature at pressures up to 101 MPa in four ultrapure ILs: 1-propyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide, 1-butyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide, 1-pentyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide, and 1-hexyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide, taking into consideration their relaxation behavior. Additionally, to further improve the reliability of the speed of sound results, the density, isentropic compressibility, and isobaric heat capacity as a function of temperature and pressure are calculated using an acoustic method.

New Method Based on the UNIFAC–VISCO Model for the Estimation of Ionic Liquids Viscosity Using the Experimental Data Recommended by Mathematical Gnostics

The viscosity of ionic liquids (ILs) has been modeled as a function of temperature and at atmospheric pressure using a new method based on the UNIFAC–VISCO method. This model extends the calculations previously reported by our group (see Zhao et al. J. Chem. Eng. Data 2016, 61, 2160–2169) which used 154 experimental viscosity data points of 25 ionic liquids for regression of a set of binary interaction parameters and ion Vogel–Fulcher–Tammann (VFT) parameters. Discrepancies in the experimental data of the same IL affect the quality of the correlation and thus the development of the predictive method. In this work, mathematical gnostics was used to analyze the experimental data from different sources and recommend one set of reliable data for each IL. These recommended data (totally 819 data points) for 70 ILs were correlated using this model to obtain an extended set of binary interaction parameters and ion VFT parameters, with a regression accuracy of 1.4%. In addition, 966 experimental viscosity data points for 11 binary mixtures of ILs were collected from literature to establish this model. All the binary data consist of 128 training data points used for the optimization of binary interaction parameters and 838 test data points used for the comparison of the pure evaluated values. The relative average absolute deviation (RAAD) for training and test is 2.9% and 3.9%, respectively.

Thomson scattering on non-thermal atmospheric pressure plasma jets

To characterize non-thermal atmospheric pressure plasmas experimentally, a large variety of methods and techniques is available, each having its own specific possibilities and limitations. A rewarding method to investigate these plasma sources is laser Thomson scattering. However, that is challenging. Non-thermal atmospheric pressure plasmas (gas temperatures close to room temperature and electron temperatures of a few eV) have usually small dimensions (below 1 mm) and a low degree of ionization (below 10^-4). Here an overview is presented of how Thomson scattering can be applied to such plasmas and used to measure directly spatially and temporally resolved the electron density and energy distribution. A general description of the scattering of photons and the guidelines for an experimental setup of this active diagnostic are provided. Special attention is given to the design concepts required to achieve the maximum signal photon flux with a minimum of unwanted signals. Recent results from the literature are also presented and discussed.