Two-phase flow occurs during the production of oil and gas in the wellbores. Modeling this phenomenon is important for monitoring well productivity and designing surface facilities. Since the transient time period in the wellbore is usually shorter than reservoir time steps, stabilized flow is assumed in the wellbore. As such, semi-steady state models are used for modeling wellbore flow dynamics. However, in the case that flow variations happen in a short period of time (i.e., a gas kick during drilling) the use of a transient two-phase model is crucial. Over the last few years, a number of numerical and analytical wellbore simulators have been developed to mimic wellbore-reservoir interaction. However, some issues still remain a concern in these studies. The main issues surrounding a comprehensive wellbore model consist of fluid property calculations, such as black-oil or compositional models, governing equations, such as mechanistic or correlation-based models, effect of temperature variation and non-isothermal assumption, and methods for coupling the wellbore to the reservoir. In most cases, only standalone wellbore models for blackoil have been used to simulate reservoir and wellbore dynamic interactions. Those models are based on simplified assumptions that lead to an unrealistic estimation of pressure and temperature distributions inside the well. In addition, most reservoir simulators use rough estimates for the perforation pressure as a coupling condition between the wellbore and the reservoir, neglecting pressure drops in the horizontal section. In this study, we present an implementation of a compositional, pseudo steady-state, non-isothermal, coupled wellbore-reservoir simulator for fluid flow in wellbores with a vertical section and a horizontal section embedded on the producing reservoir. In addition, we present the implementation of a pseudo-compositional, fully implicit, transient two-fluid model for two-phase flow in wellbores. In this model, we solve gas/liquid mass balance, gas/liquid momentum balance, and two-phase energy equations in order to obtain the five primary variables: liquid velocity, gas velocity, pressure, holdup and temperature. In our simulation, we compared stratified, bubbly, intermittent flow effects on pressure and temperature distributions in either a transient or steady-state condition. We found that flow geometry variation in different regimes can significantly affect the flow parameters. We also observed that there are significant differences in flow rate prediction between a coupled wellbore-reservoir simulator and a stand-alone reservoir simulator, at the early stages of production. The outcome of this research leads to a more accurate and reliable simulation of multiphase flow in the wellbore, which can be applied to surface facility design, well performance optimization, and wellbore damage estimation. / text

A steady-state model for the calculation of temperature and pressure distributions, and heat and work balance for the STAR-LM and the STAR-H2 systems was developed. The STAR-LM system is designed for electricity production and consists of the lead cooled reactor on natural circulation and the supercritical carbon dioxide Brayton cycle. The STAR-H2 system uses the same reactor which is coupled to the hydrogen production plant, the Brayton cycle, and the water desalination plant. The Brayton cycle produces electricity for the on-site needs. Realistic modules for each system component were developed. The model also performs design calculations for the turbine and compressors for the CO2 Brayton cycle. The model was used to optimize the performance of the entire system as well as every system component. The size of each component was calculated. For the 400 MWt reactor power the STAR-LM produces 174.4 MWe (44% efficiency) and the STAR-H2 system produces 7450 kg H2/hr. The steady state model was used to conduct quasi-static passive load follow analysis. The control strategy was developed for each system; no control action on the reactor is required. As a main safety criterion, the peak cladding temperature is used. It was demonstrated that this temperature remains below the safety limit during both normal operation and load follow.

Electrophysiologic responses are used to estimate hearing thresholds and fit hearing aids in young infants, but these estimates are not exact. An objective test of speech encoding could be used to validate infant fittings by showing that speech has been registered in the central auditory system. Such a test could also show the effects of auditory processing problems on the neural representation of speech. This thesis describes techniques for recording electrophysiologic responses to natural speech stimuli from the brainstem and auditory cortex. The first technique uses a Fourier analyzer to measure steady-state brainstem responses to periodicities and envelope changes in vowels, and the second uses a windowed cross-correlation procedure to measure cortical responses to the envelopes of sentences. Two studies were conducted with the Fourier analyzer. The first measured responses to natural vowels with steady and changing fundamentals, and changing formants. Significant responses to the fundamental were detected for all of the vowels, in all of the subjects, in 19 – 73 s (on average). The second study recorded responses to a vowel fundamental and harmonics. Vowels were presented in opposite polarities to distinguish envelope responses from responses to the spectrum. Significant envelope responses were detected in all subjects at the fundamental. Significant spectral responses were detected in most subjects at harmonics near formant peaks. The third study used cross-correlation to measure cortical responses to sentences. Significant envelope responses were detected to all sentences, at delays of roughly 180 ms. Responses were localized to the posterior auditory cortices. A model based on a series of overlapping transient responses to envelope changes could also account for the results, suggesting that the cortex either directly follows the speech envelope or consistently reacts to changes in this envelope. The strengths and weaknesses of both techniques are discussed in relation to their potential clinical applications.

No description available.

Discrete event simulation is well known to be a powerful approach to investigate behaviour of complex dynamic stochastic systems, especially when the system is analytically not tractable. The estimation of mean values has traditionally been the main goal of simulation output analysis, even though it provides limited information about the analysed system's performance. Because of its complexity, quantile analysis is not as frequently applied, despite its ability to provide much deeper insights into the system of interest. A set of quantiles can be used to approximate a cumulative distribution function, providing fuller information about a given performance characteristic of the simulated system. This thesis employs the distributed computing power of multiple computers by proposing new methods for sequential and automated analysis of quantile-based performance measures of such dynamic systems. These new methods estimate steady state quantiles based on replicating simulations on clusters of workstations as simulation engines. A general contribution to the problem of the length of the initial transient is made by considering steady state in terms of the underlying probability distribution. Our research focuses on sequential and automated methods to guarantee a satisfactory level of confidence of the final results. The correctness of the proposed methods has been exhaustively studied by means of sequential coverage analysis. Quantile estimates are used to investigate underlying probability distributions. We demonstrate that synchronous replications greatly assist this kind of analysis.

This thesis describes a new approach to accelerate the simulation of the steady-state response of nonlinear circuits using the Harmonic Balance (HB) technique. The approach presented in this work focuses on direct factorization of the sparse Jacobian matrix of the HB nonlinear equations using a Graphics Processing Unit (GPU) platform. This approach exploits the heterogeneous structure of the Jacobian matrix. The computational core of the proposed approach is based on developing a block-wise version of the KLU factorization algorithm, where scalar arithmetic operations are replaced by block-aware matrix operations. For a large number of harmonics, or excitation tones, or both the Block-KLU (BKLU) approach effectively raises the ratio of floating-point operations to other operations and, therefore, becomes an ideal vehicle for implementation on a GPU-based platform. Motivated by this fact, a GPU-based Hybrid Block KLU framework is developed to implement the BKLU. The proposed approach in this thesis is named Hybrid-BKLU. The Hybrid-BKLU is implemented in two parts, on the host CPU and on the graphic card’s GPU, using the OpenCL heterogeneous parallel programming language. To show the efficiency of the Hybrid-BKLU approach, its performance is compared with BKLU approach performing HB analysis on several test circuits. The Hybrid-BKLU approach yields speedup by up to 89 times over conventional BKLU on CPU.

Electrophysiologic responses are used to estimate hearing thresholds and fit hearing aids in young infants, but these estimates are not exact. An objective test of speech encoding could be used to validate infant fittings by showing that speech has been registered in the central auditory system. Such a test could also show the effects of auditory processing problems on the neural representation of speech. This thesis describes techniques for recording electrophysiologic responses to natural speech stimuli from the brainstem and auditory cortex. The first technique uses a Fourier analyzer to measure steady-state brainstem responses to periodicities and envelope changes in vowels, and the second uses a windowed cross-correlation procedure to measure cortical responses to the envelopes of sentences. Two studies were conducted with the Fourier analyzer. The first measured responses to natural vowels with steady and changing fundamentals, and changing formants. Significant responses to the fundamental were detected for all of the vowels, in all of the subjects, in 19 – 73 s (on average). The second study recorded responses to a vowel fundamental and harmonics. Vowels were presented in opposite polarities to distinguish envelope responses from responses to the spectrum. Significant envelope responses were detected in all subjects at the fundamental. Significant spectral responses were detected in most subjects at harmonics near formant peaks. The third study used cross-correlation to measure cortical responses to sentences. Significant envelope responses were detected to all sentences, at delays of roughly 180 ms. Responses were localized to the posterior auditory cortices. A model based on a series of overlapping transient responses to envelope changes could also account for the results, suggesting that the cortex either directly follows the speech envelope or consistently reacts to changes in this envelope. The strengths and weaknesses of both techniques are discussed in relation to their potential clinical applications.

An attempt was made to calibrate a steady state activated sludge model (ASM2d) for the biological nutrient removal process at the Gold bar wastewater treatment plant. This calibrated model could be used on a regular basis to test various operational strategies and predict effluent quality under different scenario. To achieve this historic data from the plant database was collected based on 24 composite samples. A trial and error method of wastewater characterization of the primary effluent was attempted using the influent advisor module of the GPS-[X] software. Sensitivity analysis of kinetic parameters was carried out and the most important ones identified were calibrated (default values were modified) based on literature. After calibration it was observed that the model was overestimating the concentrations of carbonaceous biological oxygen demand, total suspended solids and orthophosphate in the effluent, compared to the actual value measured at the plant. Similarly the effluent ammonia concentration was underestimated for most days along with the nitrate and nitrite concentration. This clearly indicated the need for a more accurate calibration based on experimental data to improve prediction capabilities and the reliability of the model. / Environmental Science

Incremental sheet forming process is a relatively new method of forming which is increasingly being used in the industry. Complex shapes can be manufactured using this method and the forming operation doesn't require any dies. High strains of over 300 % can also be achieved. Incremental sheet forming method is used to manufacture many different components presently. Prototype examples include car headlights, tubs, train body panels and medical products. The work done in the thesis deals with the prediction of the steady state forces acting on the tool during forming. Prediction of forces generated would help to design the machine against excessive vibrations. It would help the user to protect the tool and the material blank from failure. An efficient design ensures that the tool would not get deflected out of its path while forming, improving the accuracy of the finished part. To study the forces, experiments were conducted by forming pyramid and cone shapes. An experimental arrangement was set up and experimental data was collected using a data acquisition system. The effect that the various process parameters, like the thickness of the sheet, wall angle of the part and tool diameter had on the steady state force were studied. A three dimensional model was developed using commercial finite element software ABAQUS using a new modeling technique to simulate the deformation of the sheet metal blank during incremental sheet forming. The steady state forces generated for any shape, with any set of parameters used, could be predicted using the numerical model. The advantage of having a numerical model is that the forces can be predicted without doing experiments. The model was used to predict the steady state forces developed during forming of pyramid and cone shapes. The results were compared and were seen to be reasonably close to the experimental results. Later, the numerical model was validated by forming arbitrary shapes and comparing the value obtained from simulations to the value of the measured steady state forces. The results obtained from the numerical model were seen to match very well with the experimental forces for the new shapes. The numerical model developed using the new technique was seen to predict forces to a reasonable extent with less computational time as compared to the models currently available.

This thesis describes a new approach to accelerate the simulation of the steady-state response of nonlinear circuits using the Harmonic Balance (HB) technique. The approach presented in this work focuses on direct factorization of the sparse Jacobian matrix of the HB nonlinear equations using a Graphics Processing Unit (GPU) platform. This approach exploits the heterogeneous structure of the Jacobian matrix. The computational core of the proposed approach is based on developing a block-wise version of the KLU factorization algorithm, where scalar arithmetic operations are replaced by block-aware matrix operations. For a large number of harmonics, or excitation tones, or both the Block-KLU (BKLU) approach effectively raises the ratio of floating-point operations to other operations and, therefore, becomes an ideal vehicle for implementation on a GPU-based platform. Motivated by this fact, a GPU-based Hybrid Block KLU framework is developed to implement the BKLU. The proposed approach in this thesis is named Hybrid-BKLU. The Hybrid-BKLU is implemented in two parts, on the host CPU and on the graphic card’s GPU, using the OpenCL heterogeneous parallel programming language. To show the efficiency of the Hybrid-BKLU approach, its performance is compared with BKLU approach performing HB analysis on several test circuits. The Hybrid-BKLU approach yields speedup by up to 89 times over conventional BKLU on CPU.