3-dimensional Oil Drift Simulations
NASA Astrophysics Data System (ADS)
Wettre, C.; Reistad, M.; Hjøllo, B.Å.
Simulation of oil drift has been an ongoing activity at the Norwegian Meteorological Institute since the 1970's. The Marine Forecasting Centre provides a 24-hour service for the Norwegian Pollution Control Authority and the oil companies operating in the Norwegian sector. The response time is 30 minutes. From 2002 the service is extended to simulation of oil drift from oil spills in deep water, using the DeepBlow model developed by SINTEF Applied Chemistry. The oil drift model can be applied both for instantaneous and continuous releases. The changes in the mass of oil and emulsion as a result of evaporation and emulsion are computed. For oil spill at deep water, hydrate formation and gas dissolution are taken into account. The properties of the oil depend on the oil type, and in the present version 64 different types of oil can be simulated. For accurate oil drift simulations it is important to have the best possible data on the atmospheric and oceanic conditions. The oil drift simulations at the Norwegian Meteorological Institute are always based on the most updated data from numerical models of the atmosphere and the ocean. The drift of the surface oil is computed from the vectorial sum of the surface current from the ocean model and the wave induced Stokes drift computed from wave energy spectra from the wave prediction model. In the new model the current distribution with depth is taken into account when calculating the drift of the dispersed oil droplets. Salinity and temperature profiles from the ocean model are needed in the DeepBlow model. The result of the oil drift simulations can be plotted on sea charts used for navigation, either as trajectory plots or particle plots showing the situation at a given time. The results can also be sent as data files to be included in the user's own GIS system.
MacOil--An Oil Exploration Simulation.
ERIC Educational Resources Information Center
Burger, H. Robert
1989-01-01
Described is a simulation program designed to teach the geology of oil exploration in a game format. Discussed are information and tools provided, information management, data sets, and availability of the program. (CW)
Time accurate simulations of compressible shear flows
NASA Technical Reports Server (NTRS)
Givi, Peyman; Steinberger, Craig J.; Vidoni, Thomas J.; Madnia, Cyrus K.
1993-01-01
The objectives of this research are to employ direct numerical simulation (DNS) to study the phenomenon of mixing (or lack thereof) in compressible free shear flows and to suggest new means of enhancing mixing in such flows. The shear flow configurations under investigation are those of parallel mixing layers and planar jets under both non-reacting and reacting nonpremixed conditions. During the three-years of this research program, several important issues regarding mixing and chemical reactions in compressible shear flows were investigated.
Accurate simulation of optical properties in dyes.
Jacquemin, Denis; Perpète, Eric A; Ciofini, Ilaria; Adamo, Carlo
2009-02-17
Since Antiquity, humans have produced and commercialized dyes. To this day, extraction of natural dyes often requires lengthy and costly procedures. In the 19th century, global markets and new industrial products drove a significant effort to synthesize artificial dyes, characterized by low production costs, huge quantities, and new optical properties (colors). Dyes that encompass classes of molecules absorbing in the UV-visible part of the electromagnetic spectrum now have a wider range of applications, including coloring (textiles, food, paintings), energy production (photovoltaic cells, OLEDs), or pharmaceuticals (diagnostics, drugs). Parallel to the growth in dye applications, researchers have increased their efforts to design and synthesize new dyes to customize absorption and emission properties. In particular, dyes containing one or more metallic centers allow for the construction of fairly sophisticated systems capable of selectively reacting to light of a given wavelength and behaving as molecular devices (photochemical molecular devices, PMDs).Theoretical tools able to predict and interpret the excited-state properties of organic and inorganic dyes allow for an efficient screening of photochemical centers. In this Account, we report recent developments defining a quantitative ab initio protocol (based on time-dependent density functional theory) for modeling dye spectral properties. In particular, we discuss the importance of several parameters, such as the methods used for electronic structure calculations, solvent effects, and statistical treatments. In addition, we illustrate the performance of such simulation tools through case studies. We also comment on current weak points of these methods and ways to improve them. PMID:19113946
Next generation oil reservoir simulations
Joubert, W.
1996-04-01
This paper describes a collaborative effort between Amoco Production Company, Los Alamos National Laboratory and Cray Research Inc. to develop a next-generation massively parallel oil reservoir simulation code. The simulator, code-named Falcon, enables highly detailed simulations to be performed on a range of platforms such as the Cray T3D and T3E. The code is currently being used by Amoco to perform a sophisticated field study using multiple geostatistical realizations on a scale of 2-5 million grid blocks and 1000-2000 wells. In this paper we discuss the nature of this collaborative effort, the software design and engineering aspects of the code, parallelization experiences, and performance studies. The code will be marketed to the oil industry by a third-party independent software vendor in mid-1996.
RETORT. Oil Shale Retorting Simulation
Eyberger, L.R.
1992-02-26
RETORT is a one-dimensional mathematical model for simulating the chemical and physical processes involved in the vertical retorting of a fixed or moving rubbled bed of oil shale. The model includes those processes believed to have the most important effects in either the hot-gas retorting mode or the forward combustion retorting mode. The physical processes are: axial convective transport of heat and mass, axial thermal dispersion, axial pressure drop, gas-solid heat transfer, intraparticle thermal conductivity, water evaporation and condensation, wall heat loss, and movement of shale countercurrent to flow of gas. The chemical reactions within the shale particles are: release of bound water, pyrolysis of kerogen, coking of oil, pyrolysis of char, decomposition of carbonate minerals, and gasification of residual organic carbon with CO2, H2O, and O2. The chemical reactions in the bulk-gas stream are: combustion and cracking of oil vapor, combustion of H2, CH4, CHx, and CO, and the water-gas shift. The RETORT model is meant to simulate adiabatic laboratory retorts and in situ retorts that have been prepared with fairly uniform lateral distribution of shale particle sizes, void volume, and permeability. The model`s main role is to calculate, as a function of time and axial location in the retort, the flow rate of the bulk-gas stream and the composition and temperature of both the fluid stream and the shale particles.
RETORT. Oil Shale Retorting Simulation
Braun, R.L.
1992-02-26
RETORT is a one-dimensional mathematical model for simulating the chemical and physical processes involved in the vertical retorting of a fixed or moving rubbled bed of oil shale. The model includes those processes believed to have the most important effects in either the hot-gas retorting mode or the forward combustion retorting mode. The physical processes are: axial convective transport of heat and mass, axial thermal dispersion, axial pressure drop, gas-solid heat transfer, intraparticle thermal conductivity, water evaporation and condensation, wall heat loss, and movement of shale countercurrent to flow of gas. The chemical reactions within the shale particles are: release of bound water, pyrolysis of kerogen, coking of oil, pyrolysis of char, decomposition of carbonate minerals, and gasification of residual organic carbon with CO2, H2O, and O2. The chemical reactions in the bulk-gas stream are: combustion and cracking of oil vapor, combustion of H2, CH4, CHx, and CO, and the water- gas shift. The RETORT model is meant to simulate adiabatic laboratory retorts and in situ retorts that have been prepared with fairly uniform lateral distribution of shale particle sizes, void volume, and permeability. The model`s main role is to calculate, as a function of time and axial location in the retort, the flow rate of the bulk-gas stream and the composition and temperature of both the fluid stream and the shale particles.
D-BRAIN: Anatomically Accurate Simulated Diffusion MRI Brain Data.
Perrone, Daniele; Jeurissen, Ben; Aelterman, Jan; Roine, Timo; Sijbers, Jan; Pizurica, Aleksandra; Leemans, Alexander; Philips, Wilfried
2016-01-01
Diffusion Weighted (DW) MRI allows for the non-invasive study of water diffusion inside living tissues. As such, it is useful for the investigation of human brain white matter (WM) connectivity in vivo through fiber tractography (FT) algorithms. Many DW-MRI tailored restoration techniques and FT algorithms have been developed. However, it is not clear how accurately these methods reproduce the WM bundle characteristics in real-world conditions, such as in the presence of noise, partial volume effect, and a limited spatial and angular resolution. The difficulty lies in the lack of a realistic brain phantom on the one hand, and a sufficiently accurate way of modeling the acquisition-related degradation on the other. This paper proposes a software phantom that approximates a human brain to a high degree of realism and that can incorporate complex brain-like structural features. We refer to it as a Diffusion BRAIN (D-BRAIN) phantom. Also, we propose an accurate model of a (DW) MRI acquisition protocol to allow for validation of methods in realistic conditions with data imperfections. The phantom model simulates anatomical and diffusion properties for multiple brain tissue components, and can serve as a ground-truth to evaluate FT algorithms, among others. The simulation of the acquisition process allows one to include noise, partial volume effects, and limited spatial and angular resolution in the images. In this way, the effect of image artifacts on, for instance, fiber tractography can be investigated with great detail. The proposed framework enables reliable and quantitative evaluation of DW-MR image processing and FT algorithms at the level of large-scale WM structures. The effect of noise levels and other data characteristics on cortico-cortical connectivity and tractography-based grey matter parcellation can be investigated as well. PMID:26930054
D-BRAIN: Anatomically Accurate Simulated Diffusion MRI Brain Data
Perrone, Daniele; Jeurissen, Ben; Aelterman, Jan; Roine, Timo; Sijbers, Jan; Pizurica, Aleksandra; Leemans, Alexander; Philips, Wilfried
2016-01-01
Diffusion Weighted (DW) MRI allows for the non-invasive study of water diffusion inside living tissues. As such, it is useful for the investigation of human brain white matter (WM) connectivity in vivo through fiber tractography (FT) algorithms. Many DW-MRI tailored restoration techniques and FT algorithms have been developed. However, it is not clear how accurately these methods reproduce the WM bundle characteristics in real-world conditions, such as in the presence of noise, partial volume effect, and a limited spatial and angular resolution. The difficulty lies in the lack of a realistic brain phantom on the one hand, and a sufficiently accurate way of modeling the acquisition-related degradation on the other. This paper proposes a software phantom that approximates a human brain to a high degree of realism and that can incorporate complex brain-like structural features. We refer to it as a Diffusion BRAIN (D-BRAIN) phantom. Also, we propose an accurate model of a (DW) MRI acquisition protocol to allow for validation of methods in realistic conditions with data imperfections. The phantom model simulates anatomical and diffusion properties for multiple brain tissue components, and can serve as a ground-truth to evaluate FT algorithms, among others. The simulation of the acquisition process allows one to include noise, partial volume effects, and limited spatial and angular resolution in the images. In this way, the effect of image artifacts on, for instance, fiber tractography can be investigated with great detail. The proposed framework enables reliable and quantitative evaluation of DW-MR image processing and FT algorithms at the level of large-scale WM structures. The effect of noise levels and other data characteristics on cortico-cortical connectivity and tractography-based grey matter parcellation can be investigated as well. PMID:26930054
Accurate Langevin approaches to simulate Markovian channel dynamics
NASA Astrophysics Data System (ADS)
Huang, Yandong; Rüdiger, Sten; Shuai, Jianwei
2015-12-01
The stochasticity of ion-channels dynamic is significant for physiological processes on neuronal cell membranes. Microscopic simulations of the ion-channel gating with Markov chains can be considered to be an accurate standard. However, such Markovian simulations are computationally demanding for membrane areas of physiologically relevant sizes, which makes the noise-approximating or Langevin equation methods advantageous in many cases. In this review, we discuss the Langevin-like approaches, including the channel-based and simplified subunit-based stochastic differential equations proposed by Fox and Lu, and the effective Langevin approaches in which colored noise is added to deterministic differential equations. In the framework of Fox and Lu’s classical models, several variants of numerical algorithms, which have been recently developed to improve accuracy as well as efficiency, are also discussed. Through the comparison of different simulation algorithms of ion-channel noise with the standard Markovian simulation, we aim to reveal the extent to which the existing Langevin-like methods approximate results using Markovian methods. Open questions for future studies are also discussed.
Progress in fast, accurate multi-scale climate simulations
Collins, W. D.; Johansen, H.; Evans, K. J.; Woodward, C. S.; Caldwell, P. M.
2015-06-01
We present a survey of physical and computational techniques that have the potential to contribute to the next generation of high-fidelity, multi-scale climate simulations. Examples of the climate science problems that can be investigated with more depth with these computational improvements include the capture of remote forcings of localized hydrological extreme events, an accurate representation of cloud features over a range of spatial and temporal scales, and parallel, large ensembles of simulations to more effectively explore model sensitivities and uncertainties. Numerical techniques, such as adaptive mesh refinement, implicit time integration, and separate treatment of fast physical time scales are enablingmore » improved accuracy and fidelity in simulation of dynamics and allowing more complete representations of climate features at the global scale. At the same time, partnerships with computer science teams have focused on taking advantage of evolving computer architectures such as many-core processors and GPUs. As a result, approaches which were previously considered prohibitively costly have become both more efficient and scalable. In combination, progress in these three critical areas is poised to transform climate modeling in the coming decades.« less
Progress in fast, accurate multi-scale climate simulations
Collins, W. D.; Johansen, H.; Evans, K. J.; Woodward, C. S.; Caldwell, P. M.
2015-06-01
We present a survey of physical and computational techniques that have the potential to contribute to the next generation of high-fidelity, multi-scale climate simulations. Examples of the climate science problems that can be investigated with more depth with these computational improvements include the capture of remote forcings of localized hydrological extreme events, an accurate representation of cloud features over a range of spatial and temporal scales, and parallel, large ensembles of simulations to more effectively explore model sensitivities and uncertainties. Numerical techniques, such as adaptive mesh refinement, implicit time integration, and separate treatment of fast physical time scales are enabling improved accuracy and fidelity in simulation of dynamics and allowing more complete representations of climate features at the global scale. At the same time, partnerships with computer science teams have focused on taking advantage of evolving computer architectures such as many-core processors and GPUs. As a result, approaches which were previously considered prohibitively costly have become both more efficient and scalable. In combination, progress in these three critical areas is poised to transform climate modeling in the coming decades.
How Accurate Are Transition States from Simulations of Enzymatic Reactions?
2015-01-01
The rate expression of traditional transition state theory (TST) assumes no recrossing of the transition state (TS) and thermal quasi-equilibrium between the ground state and the TS. Currently, it is not well understood to what extent these assumptions influence the nature of the activated complex obtained in traditional TST-based simulations of processes in the condensed phase in general and in enzymes in particular. Here we scrutinize these assumptions by characterizing the TSs for hydride transfer catalyzed by the enzyme Escherichia coli dihydrofolate reductase obtained using various simulation approaches. Specifically, we compare the TSs obtained with common TST-based methods and a dynamics-based method. Using a recently developed accurate hybrid quantum mechanics/molecular mechanics potential, we find that the TST-based and dynamics-based methods give considerably different TS ensembles. This discrepancy, which could be due equilibrium solvation effects and the nature of the reaction coordinate employed and its motion, raises major questions about how to interpret the TSs determined by common simulation methods. We conclude that further investigation is needed to characterize the impact of various TST assumptions on the TS phase-space ensemble and on the reaction kinetics. PMID:24860275
Progress in Fast, Accurate Multi-scale Climate Simulations
Collins, William D; Johansen, Hans; Evans, Katherine J; Woodward, Carol S.; Caldwell, Peter
2015-01-01
We present a survey of physical and computational techniques that have the potential to con- tribute to the next generation of high-fidelity, multi-scale climate simulations. Examples of the climate science problems that can be investigated with more depth include the capture of remote forcings of localized hydrological extreme events, an accurate representation of cloud features over a range of spatial and temporal scales, and parallel, large ensembles of simulations to more effectively explore model sensitivities and uncertainties. Numerical techniques, such as adaptive mesh refinement, implicit time integration, and separate treatment of fast physical time scales are enabling improved accuracy and fidelity in simulation of dynamics and allow more complete representations of climate features at the global scale. At the same time, part- nerships with computer science teams have focused on taking advantage of evolving computer architectures, such as many-core processors and GPUs, so that these approaches which were previously considered prohibitively costly have become both more efficient and scalable. In combination, progress in these three critical areas is poised to transform climate modeling in the coming decades.
Anisotropic Turbulence Modeling for Accurate Rod Bundle Simulations
Baglietto, Emilio
2006-07-01
An improved anisotropic eddy viscosity model has been developed for accurate predictions of the thermal hydraulic performances of nuclear reactor fuel assemblies. The proposed model adopts a non-linear formulation of the stress-strain relationship in order to include the reproduction of the anisotropic phenomena, and in combination with an optimized low-Reynolds-number formulation based on Direct Numerical Simulation (DNS) to produce correct damping of the turbulent viscosity in the near wall region. This work underlines the importance of accurate anisotropic modeling to faithfully reproduce the scale of the turbulence driven secondary flows inside the bundle subchannels, by comparison with various isothermal and heated experimental cases. The very low scale secondary motion is responsible for the increased turbulence transport which produces a noticeable homogenization of the velocity distribution and consequently of the circumferential cladding temperature distribution, which is of main interest in bundle design. Various fully developed bare bundles test cases are shown for different geometrical and flow conditions, where the proposed model shows clearly improved predictions, in close agreement with experimental findings, for regular as well as distorted geometries. Finally the applicability of the model for practical bundle calculations is evaluated through its application in the high-Reynolds form on coarse grids, with excellent results. (author)
Symphony: a framework for accurate and holistic WSN simulation.
Riliskis, Laurynas; Osipov, Evgeny
2015-01-01
Research on wireless sensor networks has progressed rapidly over the last decade, and these technologies have been widely adopted for both industrial and domestic uses. Several operating systems have been developed, along with a multitude of network protocols for all layers of the communication stack. Industrial Wireless Sensor Network (WSN) systems must satisfy strict criteria and are typically more complex and larger in scale than domestic systems. Together with the non-deterministic behavior of network hardware in real settings, this greatly complicates the debugging and testing of WSN functionality. To facilitate the testing, validation, and debugging of large-scale WSN systems, we have developed a simulation framework that accurately reproduces the processes that occur inside real equipment, including both hardware- and software-induced delays. The core of the framework consists of a virtualized operating system and an emulated hardware platform that is integrated with the general purpose network simulator ns-3. Our framework enables the user to adjust the real code base as would be done in real deployments and also to test the boundary effects of different hardware components on the performance of distributed applications and protocols. Additionally we have developed a clock emulator with several different skew models and a component that handles sensory data feeds. The new framework should substantially shorten WSN application development cycles. PMID:25723144
Symphony: A Framework for Accurate and Holistic WSN Simulation
Riliskis, Laurynas; Osipov, Evgeny
2015-01-01
Research on wireless sensor networks has progressed rapidly over the last decade, and these technologies have been widely adopted for both industrial and domestic uses. Several operating systems have been developed, along with a multitude of network protocols for all layers of the communication stack. Industrial Wireless Sensor Network (WSN) systems must satisfy strict criteria and are typically more complex and larger in scale than domestic systems. Together with the non-deterministic behavior of network hardware in real settings, this greatly complicates the debugging and testing of WSN functionality. To facilitate the testing, validation, and debugging of large-scale WSN systems, we have developed a simulation framework that accurately reproduces the processes that occur inside real equipment, including both hardware- and software-induced delays. The core of the framework consists of a virtualized operating system and an emulated hardware platform that is integrated with the general purpose network simulator ns-3. Our framework enables the user to adjust the real code base as would be done in real deployments and also to test the boundary effects of different hardware components on the performance of distributed applications and protocols. Additionally we have developed a clock emulator with several different skew models and a component that handles sensory data feeds. The new framework should substantially shorten WSN application development cycles. PMID:25723144
Accurate bs and w testing important for crude-oil custody transfer
Williams, J. )
1990-11-12
This paper discusses how monitoring crude-oil sediment and water content at the field production site is essential in accurate crude-oil custody transfer operations. This is accomplished by manual methods, or on-line devices like capacitance, density, or energy-absorption analyzers. For custody-transfer purposes, sediment and water is determined by a test which follows one of the API manuals of petroleum measurement standards (MPMS). Typically, this test is conducted in the field by the field centrifuge method which, if performed properly, yields very accurate results. Laboratory tests can be performed, but sample handling becomes even more critical.
Massively Parallel Processing for Fast and Accurate Stamping Simulations
NASA Astrophysics Data System (ADS)
Gress, Jeffrey J.; Xu, Siguang; Joshi, Ramesh; Wang, Chuan-tao; Paul, Sabu
2005-08-01
The competitive automotive market drives automotive manufacturers to speed up the vehicle development cycles and reduce the lead-time. Fast tooling development is one of the key areas to support fast and short vehicle development programs (VDP). In the past ten years, the stamping simulation has become the most effective validation tool in predicting and resolving all potential formability and quality problems before the dies are physically made. The stamping simulation and formability analysis has become an critical business segment in GM math-based die engineering process. As the simulation becomes as one of the major production tools in engineering factory, the simulation speed and accuracy are the two of the most important measures for stamping simulation technology. The speed and time-in-system of forming analysis becomes an even more critical to support the fast VDP and tooling readiness. Since 1997, General Motors Die Center has been working jointly with our software vendor to develop and implement a parallel version of simulation software for mass production analysis applications. By 2001, this technology was matured in the form of distributed memory processing (DMP) of draw die simulations in a networked distributed memory computing environment. In 2004, this technology was refined to massively parallel processing (MPP) and extended to line die forming analysis (draw, trim, flange, and associated spring-back) running on a dedicated computing environment. The evolution of this technology and the insight gained through the implementation of DM0P/MPP technology as well as performance benchmarks are discussed in this publication.
The FLUKA Code: An Accurate Simulation Tool for Particle Therapy
Battistoni, Giuseppe; Bauer, Julia; Boehlen, Till T.; Cerutti, Francesco; Chin, Mary P. W.; Dos Santos Augusto, Ricardo; Ferrari, Alfredo; Ortega, Pablo G.; Kozłowska, Wioletta; Magro, Giuseppe; Mairani, Andrea; Parodi, Katia; Sala, Paola R.; Schoofs, Philippe; Tessonnier, Thomas; Vlachoudis, Vasilis
2016-01-01
Monte Carlo (MC) codes are increasingly spreading in the hadrontherapy community due to their detailed description of radiation transport and interaction with matter. The suitability of a MC code for application to hadrontherapy demands accurate and reliable physical models capable of handling all components of the expected radiation field. This becomes extremely important for correctly performing not only physical but also biologically based dose calculations, especially in cases where ions heavier than protons are involved. In addition, accurate prediction of emerging secondary radiation is of utmost importance in innovative areas of research aiming at in vivo treatment verification. This contribution will address the recent developments of the FLUKA MC code and its practical applications in this field. Refinements of the FLUKA nuclear models in the therapeutic energy interval lead to an improved description of the mixed radiation field as shown in the presented benchmarks against experimental data with both 4He and 12C ion beams. Accurate description of ionization energy losses and of particle scattering and interactions lead to the excellent agreement of calculated depth–dose profiles with those measured at leading European hadron therapy centers, both with proton and ion beams. In order to support the application of FLUKA in hospital-based environments, Flair, the FLUKA graphical interface, has been enhanced with the capability of translating CT DICOM images into voxel-based computational phantoms in a fast and well-structured way. The interface is capable of importing also radiotherapy treatment data described in DICOM RT standard. In addition, the interface is equipped with an intuitive PET scanner geometry generator and automatic recording of coincidence events. Clinically, similar cases will be presented both in terms of absorbed dose and biological dose calculations describing the various available features. PMID:27242956
The FLUKA Code: An Accurate Simulation Tool for Particle Therapy.
Battistoni, Giuseppe; Bauer, Julia; Boehlen, Till T; Cerutti, Francesco; Chin, Mary P W; Dos Santos Augusto, Ricardo; Ferrari, Alfredo; Ortega, Pablo G; Kozłowska, Wioletta; Magro, Giuseppe; Mairani, Andrea; Parodi, Katia; Sala, Paola R; Schoofs, Philippe; Tessonnier, Thomas; Vlachoudis, Vasilis
2016-01-01
Monte Carlo (MC) codes are increasingly spreading in the hadrontherapy community due to their detailed description of radiation transport and interaction with matter. The suitability of a MC code for application to hadrontherapy demands accurate and reliable physical models capable of handling all components of the expected radiation field. This becomes extremely important for correctly performing not only physical but also biologically based dose calculations, especially in cases where ions heavier than protons are involved. In addition, accurate prediction of emerging secondary radiation is of utmost importance in innovative areas of research aiming at in vivo treatment verification. This contribution will address the recent developments of the FLUKA MC code and its practical applications in this field. Refinements of the FLUKA nuclear models in the therapeutic energy interval lead to an improved description of the mixed radiation field as shown in the presented benchmarks against experimental data with both (4)He and (12)C ion beams. Accurate description of ionization energy losses and of particle scattering and interactions lead to the excellent agreement of calculated depth-dose profiles with those measured at leading European hadron therapy centers, both with proton and ion beams. In order to support the application of FLUKA in hospital-based environments, Flair, the FLUKA graphical interface, has been enhanced with the capability of translating CT DICOM images into voxel-based computational phantoms in a fast and well-structured way. The interface is capable of importing also radiotherapy treatment data described in DICOM RT standard. In addition, the interface is equipped with an intuitive PET scanner geometry generator and automatic recording of coincidence events. Clinically, similar cases will be presented both in terms of absorbed dose and biological dose calculations describing the various available features. PMID:27242956
Exploring accurate Poisson–Boltzmann methods for biomolecular simulations
Wang, Changhao; Wang, Jun; Cai, Qin; Li, Zhilin; Zhao, Hong-Kai; Luo, Ray
2013-01-01
Accurate and efficient treatment of electrostatics is a crucial step in computational analyses of biomolecular structures and dynamics. In this study, we have explored a second-order finite-difference numerical method to solve the widely used Poisson–Boltzmann equation for electrostatic analyses of realistic bio-molecules. The so-called immersed interface method was first validated and found to be consistent with the classical weighted harmonic averaging method for a diversified set of test biomolecules. The numerical accuracy and convergence behaviors of the new method were next analyzed in its computation of numerical reaction field grid potentials, energies, and atomic solvation forces. Overall similar convergence behaviors were observed as those by the classical method. Interestingly, the new method was found to deliver more accurate and better-converged grid potentials than the classical method on or nearby the molecular surface, though the numerical advantage of the new method is reduced when grid potentials are extrapolated to the molecular surface. Our exploratory study indicates the need for further improving interpolation/extrapolation schemes in addition to the developments of higher-order numerical methods that have attracted most attention in the field. PMID:24443709
Developing accurate simulations for high-speed fiber links
NASA Astrophysics Data System (ADS)
Searcy, Steven; Stark, Andrew; Hsueh, Yu-Ting; Detwiler, Thomas; Tibuleac, Sorin; Chang, GK; Ralph, Stephen E.
2011-01-01
Reliable simulations of high-speed fiber optic links are necessary to understand, design, and deploy fiber networks. Laboratory experiments cannot explore all possible component variations and fiber environments that are found in today's deployed systems. Simulations typically depict relative penalties compared to a reference link. However, absolute performance metrics are required to assess actual deployment configurations. Here we detail the efforts within the Georgia Tech 100G Consortium towards achieving high absolute accuracy between simulation and experimental performance with a goal of +/-0.25 dB for back-to-back configuration, and +/-0.5 dB for transmission over multiple spans with different dispersion maps. We measure all possible component parameters including fiber length, loss, and dispersion for use in simulation. We also validate experimental methods of performance evaluation including OSNR assessment and DSP-based demodulation. We investigate a wide range of parameters including modulator chirp, polarization state, polarization dependent loss, transmit spectrum, laser linewidth, and fiber nonlinearity. We evaluate 56 Gb/s (single-polarization) and 112 Gb/s (dual-polarization) DQPSK and coherent QPSK within a 50 GHz DWDM environment with 10 Gb/s OOK adjacent channels for worst-case XPM effects. We demonstrate good simulation accuracy within linear and some nonlinear regimes for a wide range of OSNR in both back-to-back configuration and up to eight spans, over a range of launch powers. This allows us to explore a wide range of environments not available in the lab, including different fiber types, ROADM passbands, and levels of crosstalk. Continued exploration is required to validate robustness over various demodulation algorithms.
Accurate, practical simulation of satellite infrared radiometer spectral data
Sullivan, T.J.
1982-09-01
This study's purpose is to determine whether a relatively simple random band model formulation of atmospheric radiation transfer in the infrared region can provide valid simulations of narrow interval satellite-borne infrared sounder system data. Detailed ozonesondes provide the pertinent atmospheric information and sets of calibrated satellite measurements provide the validation. High resolution line-by-line model calculations are included to complete the evaluation.
Accurate simulation of terahertz transmission through doped silicon junctions
NASA Astrophysics Data System (ADS)
Jen, Chih-Yu; Richter, Christiaan
2015-03-01
In the previous work we presented results demonstrating the ability of transmission mode terahertz time domain spectroscopy (THz-TDS) to detect doping profile differences and deviations in silicon. This capability is potentially useful for quality control in the semiconductor and photovoltaic industry. We shared subsequent experimental results revealing that terahertz interactions with both electrons and holes are strong enough to recognize both n- and p-type doping profile changes. We also displayed that the relatively long wavelength (~ 1 mm) of THz radiation allows this approach to be compatible with surface treatments like for instance the texturing (scattering layer) typically used in the solar industry. In this work we continuously demonstrate the accuracy with which current terahertz optical models can simulate the power spectrum of terahertz radiation transmitted through junctions with known doping profiles (as determined with SIMS). We conclude that current optical models predict the terahertz transmission and absorption in silicon junctions well.
Accurate Position Sensing of Defocused Beams Using Simulated Beam Templates
Awwal, A; Candy, J; Haynam, C; Widmayer, C; Bliss, E; Burkhart, S
2004-09-29
In position detection using matched filtering one is faced with the challenge of determining the best position in the presence of distortions such as defocus and diffraction noise. This work evaluates the performance of simulated defocused images as the template against the real defocused beam. It was found that an amplitude modulated phase-only filter is better equipped to deal with real defocused images that suffer from diffraction noise effects resulting in a textured spot intensity pattern. It is shown that the there is a tradeoff of performance dependent upon the type and size of the defocused image. A novel automated system was developed that can automatically select the right template type and size. Results of this automation for real defocused images are presented.
Accurate direct Eulerian simulation of dynamic elastic-plastic flow
Kamm, James R; Walter, John W
2009-01-01
The simulation of dynamic, large strain deformation is an important, difficult, and unsolved computational challenge. Existing Eulerian schemes for dynamic material response are plagued by unresolved issues. We present a new scheme for the first-order system of elasto-plasticity equations in the Eulerian frame. This system has an intrinsic constraint on the inverse deformation gradient. Standard Godunov schemes do not satisfy this constraint. The method of Flux Distributions (FD) was devised to discretely enforce such constraints for numerical schemes with cell-centered variables. We describe a Flux Distribution approach that enforces the inverse deformation gradient constraint. As this approach is new and novel, we do not yet have numerical results to validate our claims. This paper is the first installment of our program to develop this new method.
Toward the Accurate Simulation of Two-Dimensional Electronic Spectra
NASA Astrophysics Data System (ADS)
Giussani, Angelo; Nenov, Artur; Segarra-Martí, Javier; Jaiswal, Vishal K.; Rivalta, Ivan; Dumont, Elise; Mukamel, Shaul; Garavelli, Marco
2015-06-01
Two-dimensional pump-probe electronic spectroscopy is a powerful technique able to provide both high spectral and temporal resolution, allowing the analysis of ultrafast complex reactions occurring via complementary pathways by the identification of decay-specific fingerprints. [1-2] The understanding of the origin of the experimentally recorded signals in a two-dimensional electronic spectrum requires the characterization of the electronic states involved in the electronic transitions photoinduced by the pump/probe pulses in the experiment. Such a goal constitutes a considerable computational challenge, since up to 100 states need to be described, for which state-of-the-art methods as RASSCF and RASPT2 have to be wisely employed. [3] With the present contribution, the main features and potentialities of two-dimensional electronic spectroscopy are presented, together with the machinery in continuous development in our groups in order to compute two-dimensional electronic spectra. The results obtained using different level of theory and simulations are shown, bringing as examples the computed two-dimensional electronic spectra for some specific cases studied. [2-4] [1] Rivalta I, Nenov A, Cerullo G, Mukamel S, Garavelli M, Int. J. Quantum Chem., 2014, 114, 85 [2] Nenov A, Segarra-Martí J, Giussani A, Conti I, Rivalta I, Dumont E, Jaiswal V K, Altavilla S, Mukamel S, Garavelli M, Faraday Discuss. 2015, DOI: 10.1039/C4FD00175C [3] Nenov A, Giussani A, Segarra-Martí J, Jaiswal V K, Rivalta I, Cerullo G, Mukamel S, Garavelli M, J. Chem. Phys. submitted [4] Nenov A, Giussani A, Fingerhut B P, Rivalta I, Dumont E, Mukamel S, Garavelli M, Phys. Chem. Chem. Phys. Submitted [5] Krebs N, Pugliesi I, Hauer J, Riedle E, New J. Phys., 2013,15, 08501
Chen, Chun-Chi; Chen, Chao-Lieh; Lin, You-Ting
2016-01-01
This study proposes a new behavioral simulator that uses SIMULINK for all-digital CMOS time-domain smart temperature sensors (TDSTSs) for performing rapid and accurate simulations. Inverter-based TDSTSs offer the benefits of low cost and simple structure for temperature-to-digital conversion and have been developed. Typically, electronic design automation tools, such as HSPICE, are used to simulate TDSTSs for performance evaluations. However, such tools require extremely long simulation time and complex procedures to analyze the results and generate figures. In this paper, we organize simple but accurate equations into a temperature-dependent model (TDM) by which the TDSTSs evaluate temperature behavior. Furthermore, temperature-sensing models of a single CMOS NOT gate were devised using HSPICE simulations. Using the TDM and these temperature-sensing models, a novel simulator in SIMULINK environment was developed to substantially accelerate the simulation and simplify the evaluation procedures. Experiments demonstrated that the simulation results of the proposed simulator have favorable agreement with those obtained from HSPICE simulations, showing that the proposed simulator functions successfully. This is the first behavioral simulator addressing the rapid simulation of TDSTSs. PMID:27509507
The report describes procedures developed at the Oil and Hazardous Material Simulated Environmental Test Tank (OHMSETT) for sampling oil and water mixtures. Two procedures for sampling in containers are discussed: grab and stratified sampling. Both of these techniques require str...
MERCURY MASS DISTRIBUTION DURING LABORATORY AND SIMULATED IN-SITU OIL SHALE RETORTING
Total mercury mass in oil shale retort offgas was quantified in a series of laboratory retorting experiments and in a simulated modified in-situ (MIS) retorting experiment. Accurate quantitative determinations of offgas Hg mass were made possible by the use of a continuous on-lin...
Material Models for Accurate Simulation of Sheet Metal Forming and Springback
NASA Astrophysics Data System (ADS)
Yoshida, Fusahito
2010-06-01
For anisotropic sheet metals, modeling of anisotropy and the Bauschinger effect is discussed in the framework of Yoshida-Uemori kinematic hardening model incorporating with anisotropic yield functions. The performances of the models in predicting yield loci, cyclic stress-strain responses on several types of steel and aluminum sheets are demonstrated by comparing the numerical simulation results with the corresponding experimental observations. From some examples of FE simulation of sheet metal forming and springback, it is concluded that modeling of both the anisotropy and the Bauschinger effect is essential for the accurate numerical simulation.
NASA Technical Reports Server (NTRS)
Alter, Stephen J.; Brauckmann, Gregory J.; Kleb, Bil; Streett, Craig L; Glass, Christopher E.; Schuster, David M.
2015-01-01
Using the Fully Unstructured Three-Dimensional (FUN3D) computational fluid dynamics code, an unsteady, time-accurate flow field about a Space Launch System configuration was simulated at a transonic wind tunnel condition (Mach = 0.9). Delayed detached eddy simulation combined with Reynolds Averaged Naiver-Stokes and a Spallart-Almaras turbulence model were employed for the simulation. Second order accurate time evolution scheme was used to simulate the flow field, with a minimum of 0.2 seconds of simulated time to as much as 1.4 seconds. Data was collected at 480 pressure taps at locations, 139 of which matched a 3% wind tunnel model, tested in the Transonic Dynamic Tunnel (TDT) facility at NASA Langley Research Center. Comparisons between computation and experiment showed agreement within 5% in terms of location for peak RMS levels, and 20% for frequency and magnitude of power spectral densities. Grid resolution and time step sensitivity studies were performed to identify methods for improved accuracy comparisons to wind tunnel data. With limited computational resources, accurate trends for reduced vibratory loads on the vehicle were observed. Exploratory methods such as determining minimized computed errors based on CFL number and sub-iterations, as well as evaluating frequency content of the unsteady pressures and evaluation of oscillatory shock structures were used in this study to enhance computational efficiency and solution accuracy. These techniques enabled development of a set of best practices, for the evaluation of future flight vehicle designs in terms of vibratory loads.
Liquid propellant rocket engine combustion simulation with a time-accurate CFD method
NASA Technical Reports Server (NTRS)
Chen, Y. S.; Shang, H. M.; Liaw, Paul; Hutt, J.
1993-01-01
Time-accurate computational fluid dynamics (CFD) algorithms are among the basic requirements as an engineering or research tool for realistic simulations of transient combustion phenomena, such as combustion instability, transient start-up, etc., inside the rocket engine combustion chamber. A time-accurate pressure based method is employed in the FDNS code for combustion model development. This is in connection with other program development activities such as spray combustion model development and efficient finite-rate chemistry solution method implementation. In the present study, a second-order time-accurate time-marching scheme is employed. For better spatial resolutions near discontinuities (e.g., shocks, contact discontinuities), a 3rd-order accurate TVD scheme for modeling the convection terms is implemented in the FDNS code. Necessary modification to the predictor/multi-corrector solution algorithm in order to maintain time-accurate wave propagation is also investigated. Benchmark 1-D and multidimensional test cases, which include the classical shock tube wave propagation problems, resonant pipe test case, unsteady flow development of a blast tube test case, and H2/O2 rocket engine chamber combustion start-up transient simulation, etc., are investigated to validate and demonstrate the accuracy and robustness of the present numerical scheme and solution algorithm.
NASA Astrophysics Data System (ADS)
Dinpajooh, Mohammadhasan; Bai, Peng; Allan, Douglas A.; Siepmann, J. Ilja
2015-09-01
Since the seminal paper by Panagiotopoulos [Mol. Phys. 61, 813 (1997)], the Gibbs ensemble Monte Carlo (GEMC) method has been the most popular particle-based simulation approach for the computation of vapor-liquid phase equilibria. However, the validity of GEMC simulations in the near-critical region has been questioned because rigorous finite-size scaling approaches cannot be applied to simulations with fluctuating volume. Valleau [Mol. Simul. 29, 627 (2003)] has argued that GEMC simulations would lead to a spurious overestimation of the critical temperature. More recently, Patel et al. [J. Chem. Phys. 134, 024101 (2011)] opined that the use of analytical tail corrections would be problematic in the near-critical region. To address these issues, we perform extensive GEMC simulations for Lennard-Jones particles in the near-critical region varying the system size, the overall system density, and the cutoff distance. For a system with N = 5500 particles, potential truncation at 8σ and analytical tail corrections, an extrapolation of GEMC simulation data at temperatures in the range from 1.27 to 1.305 yields Tc = 1.3128 ± 0.0016, ρc = 0.316 ± 0.004, and pc = 0.1274 ± 0.0013 in excellent agreement with the thermodynamic limit determined by Potoff and Panagiotopoulos [J. Chem. Phys. 109, 10914 (1998)] using grand canonical Monte Carlo simulations and finite-size scaling. Critical properties estimated using GEMC simulations with different overall system densities (0.296 ≤ ρt ≤ 0.336) agree to within the statistical uncertainties. For simulations with tail corrections, data obtained using rcut = 3.5σ yield Tc and pc that are higher by 0.2% and 1.4% than simulations with rcut = 5 and 8σ but still with overlapping 95% confidence intervals. In contrast, GEMC simulations with a truncated and shifted potential show that rcut = 8σ is insufficient to obtain accurate results. Additional GEMC simulations for hard-core square-well particles with various ranges of the
Dinpajooh, Mohammadhasan; Bai, Peng; Allan, Douglas A.; Siepmann, J. Ilja
2015-09-21
Since the seminal paper by Panagiotopoulos [Mol. Phys. 61, 813 (1997)], the Gibbs ensemble Monte Carlo (GEMC) method has been the most popular particle-based simulation approach for the computation of vapor–liquid phase equilibria. However, the validity of GEMC simulations in the near-critical region has been questioned because rigorous finite-size scaling approaches cannot be applied to simulations with fluctuating volume. Valleau [Mol. Simul. 29, 627 (2003)] has argued that GEMC simulations would lead to a spurious overestimation of the critical temperature. More recently, Patel et al. [J. Chem. Phys. 134, 024101 (2011)] opined that the use of analytical tail corrections would be problematic in the near-critical region. To address these issues, we perform extensive GEMC simulations for Lennard-Jones particles in the near-critical region varying the system size, the overall system density, and the cutoff distance. For a system with N = 5500 particles, potential truncation at 8σ and analytical tail corrections, an extrapolation of GEMC simulation data at temperatures in the range from 1.27 to 1.305 yields T{sub c} = 1.3128 ± 0.0016, ρ{sub c} = 0.316 ± 0.004, and p{sub c} = 0.1274 ± 0.0013 in excellent agreement with the thermodynamic limit determined by Potoff and Panagiotopoulos [J. Chem. Phys. 109, 10914 (1998)] using grand canonical Monte Carlo simulations and finite-size scaling. Critical properties estimated using GEMC simulations with different overall system densities (0.296 ≤ ρ{sub t} ≤ 0.336) agree to within the statistical uncertainties. For simulations with tail corrections, data obtained using r{sub cut} = 3.5σ yield T{sub c} and p{sub c} that are higher by 0.2% and 1.4% than simulations with r{sub cut} = 5 and 8σ but still with overlapping 95% confidence intervals. In contrast, GEMC simulations with a truncated and shifted potential show that r{sub cut} = 8σ is insufficient to obtain accurate results. Additional GEMC simulations for hard
CgWind: A high-order accurate simulation tool for wind turbines and wind farms
Chand, K K; Henshaw, W D; Lundquist, K A; Singer, M A
2010-02-22
CgWind is a high-fidelity large eddy simulation (LES) tool designed to meet the modeling needs of wind turbine and wind park engineers. This tool combines several advanced computational technologies in order to model accurately the complex and dynamic nature of wind energy applications. The composite grid approach provides high-quality structured grids for the efficient implementation of high-order accurate discretizations of the incompressible Navier-Stokes equations. Composite grids also provide a natural mechanism for modeling bodies in relative motion and complex geometry. Advanced algorithms such as matrix-free multigrid, compact discretizations and approximate factorization will allow CgWind to perform highly resolved calculations efficiently on a wide class of computing resources. Also in development are nonlinear LES subgrid-scale models required to simulate the many interacting scales present in large wind turbine applications. This paper outlines our approach, the current status of CgWind and future development plans.
Margot Gerritsen
2008-10-31
Gas-injection processes are widely and increasingly used for enhanced oil recovery (EOR). In the United States, for example, EOR production by gas injection accounts for approximately 45% of total EOR production and has tripled since 1986. The understanding of the multiphase, multicomponent flow taking place in any displacement process is essential for successful design of gas-injection projects. Due to complex reservoir geometry, reservoir fluid properties and phase behavior, the design of accurate and efficient numerical simulations for the multiphase, multicomponent flow governing these processes is nontrivial. In this work, we developed, implemented and tested a streamline based solver for gas injection processes that is computationally very attractive: as compared to traditional Eulerian solvers in use by industry it computes solutions with a computational speed orders of magnitude higher and a comparable accuracy provided that cross-flow effects do not dominate. We contributed to the development of compositional streamline solvers in three significant ways: improvement of the overall framework allowing improved streamline coverage and partial streamline tracing, amongst others; parallelization of the streamline code, which significantly improves wall clock time; and development of new compositional solvers that can be implemented along streamlines as well as in existing Eulerian codes used by industry. We designed several novel ideas in the streamline framework. First, we developed an adaptive streamline coverage algorithm. Adding streamlines locally can reduce computational costs by concentrating computational efforts where needed, and reduce mapping errors. Adapting streamline coverage effectively controls mass balance errors that mostly result from the mapping from streamlines to pressure grid. We also introduced the concept of partial streamlines: streamlines that do not necessarily start and/or end at wells. This allows more efficient coverage and avoids
Improving light propagation Monte Carlo simulations with accurate 3D modeling of skin tissue
Paquit, Vincent C; Price, Jeffery R; Meriaudeau, Fabrice; Tobin Jr, Kenneth William
2008-01-01
In this paper, we present a 3D light propagation model to simulate multispectral reflectance images of large skin surface areas. In particular, we aim to simulate more accurately the effects of various physiological properties of the skin in the case of subcutaneous vein imaging compared to existing models. Our method combines a Monte Carlo light propagation model, a realistic three-dimensional model of the skin using parametric surfaces and a vision system for data acquisition. We describe our model in detail, present results from the Monte Carlo modeling and compare our results with those obtained with a well established Monte Carlo model and with real skin reflectance images.
Time-Accurate Simulations and Acoustic Analysis of Slat Free-Shear-Layer. Part 2
NASA Technical Reports Server (NTRS)
Khorrami, Mehdi R.; Singer, Bart A.; Lockard, David P.
2002-01-01
Unsteady computational simulations of a multi-element, high-lift configuration are performed. Emphasis is placed on accurate spatiotemporal resolution of the free shear layer in the slat-cove region. The excessive dissipative effects of the turbulence model, so prevalent in previous simulations, are circumvented by switching off the turbulence-production term in the slat cove region. The justifications and physical arguments for taking such a step are explained in detail. The removal of this excess damping allows the shear layer to amplify large-scale structures, to achieve a proper non-linear saturation state, and to permit vortex merging. The large-scale disturbances are self-excited, and unlike our prior fully turbulent simulations, no external forcing of the shear layer is required. To obtain the farfield acoustics, the Ffowcs Williams and Hawkings equation is evaluated numerically using the simulated time-accurate flow data. The present comparison between the computed and measured farfield acoustic spectra shows much better agreement for the amplitude and frequency content than past calculations. The effect of the angle-of-attack on the slat's flow features radiated acoustic field are also simulated presented.
High-Resolution Tsunami Inundation Simulations Based on Accurate Estimations of Coastal Waveforms
NASA Astrophysics Data System (ADS)
Oishi, Y.; Imamura, F.; Sugawara, D.; Furumura, T.
2015-12-01
We evaluate the accuracy of high-resolution tsunami inundation simulations in detail using the actual observational data of the 2011 Tohoku-Oki earthquake (Mw9.0) and investigate the methodologies to improve the simulation accuracy.Due to the recent development of parallel computing technologies, high-resolution tsunami inundation simulations are conducted more commonly than before. To evaluate how accurately these simulations can reproduce inundation processes, we test several types of simulation configurations on a parallel computer, where we can utilize the observational data (e.g., offshore and coastal waveforms and inundation properties) that are recorded during the Tohoku-Oki earthquake.Before discussing the accuracy of inundation processes on land, the incident waves at coastal sites must be accurately estimated. However, for megathrust earthquakes, it is difficult to find the tsunami source that can provide accurate estimations of tsunami waveforms at every coastal site because of the complex spatiotemporal distribution of the source and the limitation of observation. To overcome this issue, we employ a site-specific source inversion approach that increases the estimation accuracy within a specific coastal site by applying appropriate weighting to the observational data in the inversion process.We applied our source inversion technique to the Tohoku tsunami and conducted inundation simulations using 5-m resolution digital elevation model data (DEM) for the coastal area around Miyako Bay and Sendai Bay. The estimated waveforms at the coastal wave gauges of these bays successfully agree with the observed waveforms. However, the simulations overestimate the inundation extent indicating the necessity to improve the inundation model. We find that the value of Manning's roughness coefficient should be modified from the often-used value of n = 0.025 to n = 0.033 to obtain proper results at both cities.In this presentation, the simulation results with several
NASA Astrophysics Data System (ADS)
Heidari, M.; Cortes-Huerto, R.; Donadio, D.; Potestio, R.
2016-07-01
In adaptive resolution simulations the same system is concurrently modeled with different resolution in different subdomains of the simulation box, thereby enabling an accurate description in a small but relevant region, while the rest is treated with a computationally parsimonious model. In this framework, electrostatic interaction, whose accurate treatment is a crucial aspect in the realistic modeling of soft matter and biological systems, represents a particularly acute problem due to the intrinsic long-range nature of Coulomb potential. In the present work we propose and validate the usage of a short-range modification of Coulomb potential, the Damped shifted force (DSF) model, in the context of the Hamiltonian adaptive resolution simulation (H-AdResS) scheme. This approach, which is here validated on bulk water, ensures a reliable reproduction of the structural and dynamical properties of the liquid, and enables a seamless embedding in the H-AdResS framework. The resulting dual-resolution setup is implemented in the LAMMPS simulation package, and its customized version employed in the present work is made publicly available.
Turmoil: A Simulation Game Dealing With International Oil Trade
ERIC Educational Resources Information Center
Kelly, Robert
1976-01-01
This simulation game is intended to help secondary students understand the complexities of the international oil trade. Students represent nations involved in trading oil and other commodities. The game takes about five classroom periods to teach. The article includes all essential materials. (Author/RM)
Time Accurate CFD Simulations of the Orion Launch Abort Vehicle in the Transonic Regime
NASA Technical Reports Server (NTRS)
Ruf, Joseph; Rojahn, Josh
2011-01-01
Significant asymmetries in the fluid dynamics were calculated for some cases in the CFD simulations of the Orion Launch Abort Vehicle through its abort trajectories. The CFD simulations were performed steady state with symmetric boundary conditions and geometries. The trajectory points at issue were in the transonic regime, at 0 and 5 angles of attack with the Abort Motors with and without the Attitude Control Motors (ACM) firing. In some of the cases the asymmetric fluid dynamics resulted in aerodynamic side forces that were large enough that would overcome the control authority of the ACMs. MSFC s Fluid Dynamics Group supported the investigation into the cause of the flow asymmetries with time accurate CFD simulations, utilizing a hybrid RANS-LES turbulence model. The results show that the flow over the vehicle and the subsequent interaction with the AB and ACM motor plumes were unsteady. The resulting instantaneous aerodynamic forces were oscillatory with fairly large magnitudes. Time averaged aerodynamic forces were essentially symmetric.
Accurate and efficient halo-based galaxy clustering modelling with simulations
NASA Astrophysics Data System (ADS)
Zheng, Zheng; Guo, Hong
2016-06-01
Small- and intermediate-scale galaxy clustering can be used to establish the galaxy-halo connection to study galaxy formation and evolution and to tighten constraints on cosmological parameters. With the increasing precision of galaxy clustering measurements from ongoing and forthcoming large galaxy surveys, accurate models are required to interpret the data and extract relevant information. We introduce a method based on high-resolution N-body simulations to accurately and efficiently model the galaxy two-point correlation functions (2PCFs) in projected and redshift spaces. The basic idea is to tabulate all information of haloes in the simulations necessary for computing the galaxy 2PCFs within the framework of halo occupation distribution or conditional luminosity function. It is equivalent to populating galaxies to dark matter haloes and using the mock 2PCF measurements as the model predictions. Besides the accurate 2PCF calculations, the method is also fast and therefore enables an efficient exploration of the parameter space. As an example of the method, we decompose the redshift-space galaxy 2PCF into different components based on the type of galaxy pairs and show the redshift-space distortion effect in each component. The generalizations and limitations of the method are discussed.
Reducing US Oil Dependence Using Simulation
NASA Technical Reports Server (NTRS)
Ayoub, Fadi; Arnaout, Georges M.
2011-01-01
People across the world are addicted to oil; as a result, the instability of oil prices and the shortage of oil reserves have influenced human behaviors and global businesses. Today, the United States makes up only 5% of the global population but consumes 25% of the. world total energy. Most of this energy is generated from fossil fuels in the form of electricity. The contribution of this paper is to examine the possibilities of replacing fossil fuel with renewable energies to generate electricity as well as to examine other methods to reduce oil and gas consumption. We propose a system dynamics model in an attempt to predict the future US dependence on fossil fuels by using renewable energy resources such as, nuclear, wind, solar, and hydro powers. Based on the findings of our model, the study expects to provide insights towards promising solutions of the oil dependency problem.
Wijma, Hein J; Marrink, Siewert J; Janssen, Dick B
2014-07-28
Computational approaches could decrease the need for the laborious high-throughput experimental screening that is often required to improve enzymes by mutagenesis. Here, we report that using multiple short molecular dynamics (MD) simulations makes it possible to accurately model enantioselectivity for large numbers of enzyme-substrate combinations at low computational costs. We chose four different haloalkane dehalogenases as model systems because of the availability of a large set of experimental data on the enantioselective conversion of 45 different substrates. To model the enantioselectivity, we quantified the frequency of occurrence of catalytically productive conformations (near attack conformations) for pairs of enantiomers during MD simulations. We found that the angle of nucleophilic attack that leads to carbon-halogen bond cleavage was a critical variable that limited the occurrence of productive conformations; enantiomers for which this angle reached values close to 180° were preferentially converted. A cluster of 20-40 very short (10 ps) MD simulations allowed adequate conformational sampling and resulted in much better agreement to experimental enantioselectivities than single long MD simulations (22 ns), while the computational costs were 50-100 fold lower. With single long MD simulations, the dynamics of enzyme-substrate complexes remained confined to a conformational subspace that rarely changed significantly, whereas with multiple short MD simulations a larger diversity of conformations of enzyme-substrate complexes was observed. PMID:24916632
NASA Technical Reports Server (NTRS)
VanZante, Dale E.; Strazisar, Anthony J.; Wood, Jerry R,; Hathaway, Michael D.; Okiishi, Theodore H.
2000-01-01
The tip clearance flows of transonic compressor rotors are important because they have a significant impact on rotor and stage performance. While numerical simulations of these flows are quite sophisticated. they are seldom verified through rigorous comparisons of numerical and measured data because these kinds of measurements are rare in the detail necessary to be useful in high-speed machines. In this paper we compare measured tip clearance flow details (e.g. trajectory and radial extent) with corresponding data obtained from a numerical simulation. Recommendations for achieving accurate numerical simulation of tip clearance flows are presented based on this comparison. Laser Doppler Velocimeter (LDV) measurements acquired in a transonic compressor rotor, NASA Rotor 35, are used. The tip clearance flow field of this transonic rotor was simulated using a Navier-Stokes turbomachinery solver that incorporates an advanced k-epsilon turbulence model derived for flows that are not in local equilibrium. Comparison between measured and simulated results indicates that simulation accuracy is primarily dependent upon the ability of the numerical code to resolve important details of a wall-bounded shear layer formed by the relative motion between the over-tip leakage flow and the shroud wall. A simple method is presented for determining the strength of this shear layer.
NASA Astrophysics Data System (ADS)
Garrison, Stephen L.
2005-07-01
The combination of molecular simulations and potentials obtained from quantum chemistry is shown to be able to provide reasonably accurate thermodynamic property predictions. Gibbs ensemble Monte Carlo simulations are used to understand the effects of small perturbations to various regions of the model Lennard-Jones 12-6 potential. However, when the phase behavior and second virial coefficient are scaled by the critical properties calculated for each potential, the results obey a corresponding states relation suggesting a non-uniqueness problem for interaction potentials fit to experimental phase behavior. Several variations of a procedure collectively referred to as quantum mechanical Hybrid Methods for Interaction Energies (HM-IE) are developed and used to accurately estimate interaction energies from CCSD(T) calculations with a large basis set in a computationally efficient manner for the neon-neon, acetylene-acetylene, and nitrogen-benzene systems. Using these results and methods, an ab initio, pairwise-additive, site-site potential for acetylene is determined and then improved using results from molecular simulations using this initial potential. The initial simulation results also indicate that a limited range of energies important for accurate phase behavior predictions. Second virial coefficients calculated from the improved potential indicate that one set of experimental data in the literature is likely erroneous. This prescription is then applied to methanethiol. Difficulties in modeling the effects of the lone pair electrons suggest that charges on the lone pair sites negatively impact the ability of the intermolecular potential to describe certain orientations, but that the lone pair sites may be necessary to reasonably duplicate the interaction energies for several orientations. Two possible methods for incorporating the effects of three-body interactions into simulations within the pairwise-additivity formulation are also developed. A low density
Simulation of residual oil displacement in a sinusoidal channel with the lattice Boltzmann method
NASA Astrophysics Data System (ADS)
Otomo, Hiroshi; Fan, Hongli; Hazlett, Randy; Li, Yong; Staroselsky, Ilya; Zhang, Raoyang; Chen, Hudong
2015-10-01
We simulate oil slug displacement in a sinusoidal channel in order to validate computational models and algorithms for multi-component flow. This case fits in the gap between fully realistic cases characterized by complicated geometry and academic cases with simplistic geometry. Our computational model is based on the lattice Boltzmann method and allows for variation of physical parameters such as wettability and viscosity. The effect of variation of model parameters is analyzed, in particular via comparison with analytical solutions. We discuss the requirements for accurate solution of the oil slug displacement problem.
Application of the multicomponent lattice Boltzmann simulation method to oil/water dispersions
NASA Astrophysics Data System (ADS)
Seaton, M. A.; Halliday, I.; Masters, A. J.
2011-03-01
We study the propagation of acoustic fields in bounded, two-dimensional, mono-disperse oil/water emulsions using a carefully modified and appropriately calibrated single relaxation time multicomponent lattice Boltzmann equation simulation. Our model is a variant of an algorithm applying both interface forces based on macroscopic surface tensions and a kinematic condition for phase separation, adapted to allow sonic speed variations between its oil and water components. Appropriate second-order accurate acoustic boundary conditions are obtained from a node-based lattice closure with local mass conservation and applicability for varying fluid viscosities. Data from an example simulation of a single oil drop in water interacting with a generated standing acoustic wave are presented and, where appropriate, compared with empirical theories and analogous calculations for a solid object.
Development of modified cable models to simulate accurate neuronal active behaviors
2014-01-01
In large network and single three-dimensional (3-D) neuron simulations, high computing speed dictates using reduced cable models to simulate neuronal firing behaviors. However, these models are unwarranted under active conditions and lack accurate representation of dendritic active conductances that greatly shape neuronal firing. Here, realistic 3-D (R3D) models (which contain full anatomical details of dendrites) of spinal motoneurons were systematically compared with their reduced single unbranched cable (SUC, which reduces the dendrites to a single electrically equivalent cable) counterpart under passive and active conditions. The SUC models matched the R3D model's passive properties but failed to match key active properties, especially active behaviors originating from dendrites. For instance, persistent inward currents (PIC) hysteresis, frequency-current (FI) relationship secondary range slope, firing hysteresis, plateau potential partial deactivation, staircase currents, synaptic current transfer ratio, and regional FI relationships were not accurately reproduced by the SUC models. The dendritic morphology oversimplification and lack of dendritic active conductances spatial segregation in the SUC models caused significant underestimation of those behaviors. Next, SUC models were modified by adding key branching features in an attempt to restore their active behaviors. The addition of primary dendritic branching only partially restored some active behaviors, whereas the addition of secondary dendritic branching restored most behaviors. Importantly, the proposed modified models successfully replicated the active properties without sacrificing model simplicity, making them attractive candidates for running R3D single neuron and network simulations with accurate firing behaviors. The present results indicate that using reduced models to examine PIC behaviors in spinal motoneurons is unwarranted. PMID:25277743
Direct Simulations of Transition and Turbulence Using High-Order Accurate Finite-Difference Schemes
NASA Technical Reports Server (NTRS)
Rai, Man Mohan
1997-01-01
In recent years the techniques of computational fluid dynamics (CFD) have been used to compute flows associated with geometrically complex configurations. However, success in terms of accuracy and reliability has been limited to cases where the effects of turbulence and transition could be modeled in a straightforward manner. Even in simple flows, the accurate computation of skin friction and heat transfer using existing turbulence models has proved to be a difficult task, one that has required extensive fine-tuning of the turbulence models used. In more complex flows (for example, in turbomachinery flows in which vortices and wakes impinge on airfoil surfaces causing periodic transitions from laminar to turbulent flow) the development of a model that accounts for all scales of turbulence and predicts the onset of transition may prove to be impractical. Fortunately, current trends in computing suggest that it may be possible to perform direct simulations of turbulence and transition at moderate Reynolds numbers in some complex cases in the near future. This seminar will focus on direct simulations of transition and turbulence using high-order accurate finite-difference methods. The advantage of the finite-difference approach over spectral methods is that complex geometries can be treated in a straightforward manner. Additionally, finite-difference techniques are the prevailing methods in existing application codes. In this seminar high-order-accurate finite-difference methods for the compressible and incompressible formulations of the unsteady Navier-Stokes equations and their applications to direct simulations of turbulence and transition will be presented.
NASA Astrophysics Data System (ADS)
Rensonnet, Gaëtan; Jacobs, Damien; Macq, Benoît.; Taquet, Maxime
2016-03-01
Diffusion-weighted magnetic resonance imaging (DW-MRI) is a powerful tool to probe the diffusion of water through tissues. Through the application of magnetic gradients of appropriate direction, intensity and duration constituting the acquisition parameters, information can be retrieved about the underlying microstructural organization of the brain. In this context, an important and open question is to determine an optimal sequence of such acquisition parameters for a specific purpose. The use of simulated DW-MRI data for a given microstructural configuration provides a convenient and efficient way to address this problem. We first present a novel hybrid method for the synthetic simulation of DW-MRI signals that combines analytic expressions in simple geometries such as spheres and cylinders and Monte Carlo (MC) simulations elsewhere. Our hybrid method remains valid for any acquisition parameters and provides identical levels of accuracy with a computational time that is 90% shorter than that required by MC simulations for commonly-encountered microstructural configurations. We apply our novel simulation technique to estimate the radius of axons under various noise levels with different acquisition protocols commonly used in the literature. The results of our comparison suggest that protocols favoring a large number of gradient intensities such as a Cube and Sphere (CUSP) imaging provide more accurate radius estimation than conventional single-shell HARDI acquisitions for an identical acquisition time.
A fast and accurate simulator for the design of birdcage coils in MRI.
Giovannetti, Giulio; Landini, Luigi; Santarelli, Maria Filomena; Positano, Vincenzo
2002-11-01
The birdcage coils are extensively used in MRI systems since they introduce a high signal to noise ratio and a high radiofrequency magnetic field homogeneity that guarantee a large field of view. The present article describes the implementation of a birdcage coil simulator, operating in high-pass and low-pass modes, using magnetostatic analysis of the coil. Respect to other simulators described in literature, our simulator allows to obtain in short time not only the dominant frequency mode, but also the complete resonant frequency spectrum and the relevant magnetic field pattern with high accuracy. Our simulator accounts for all the inductances including the mutual inductances between conductors. Moreover, the inductance calculation includes an accurately birdcage geometry description and the effect of a radiofrequency shield. The knowledge of all the resonance modes introduced by a birdcage coil is twofold useful during birdcage coil design: --higher order modes should be pushed far from the fundamental one, --for particular applications, it is necessary to localize other resonant modes (as the Helmholtz mode) jointly to the dominant mode. The knowledge of the magnetic field pattern allows to a priori verify the field homogeneity created inside the coil, when varying the coil dimension and mainly the number of the coil legs. The coil is analyzed using equivalent circuit method. Finally, the simulator is validated by implementing a low-pass birdcage coil and comparing our data with the literature. PMID:12413563
Application of the G-JF discrete-time thermostat for fast and accurate molecular simulations
NASA Astrophysics Data System (ADS)
Grønbech-Jensen, Niels; Hayre, Natha Robert; Farago, Oded
2014-02-01
A new Langevin-Verlet thermostat that preserves the fluctuation-dissipation relationship for discrete time steps is applied to molecular modeling and tested against several popular suites (AMBER, GROMACS, LAMMPS) using a small molecule as an example that can be easily simulated by all three packages. Contrary to existing methods, the new thermostat exhibits no detectable changes in the sampling statistics as the time step is varied in the entire numerical stability range. The simple form of the method, which we express in the three common forms (Velocity-Explicit, Störmer-Verlet, and Leap-Frog), allows for easy implementation within existing molecular simulation packages to achieve faster and more accurate results with no cost in either computing time or programming complexity.
Maudlin, P.J.; Stout, M.G.
1996-09-01
Strength and fracture constitutive relationships containing strain rate dependence and thermal softening are important for accurate simulation of metal cutting. The mechanical behavior of a hardened 4340 steel was characterized using the von Mises yield function, the Mechanical Threshold Stress model and the Johnson- Cook fracture model. This constitutive description was implemented into the explicit Lagrangian FEM continuum-mechanics code EPIC, and orthogonal plane-strain metal cutting calculations were performed. Heat conduction and friction at the toolwork-piece interface were included in the simulations. These transient calculations were advanced in time until steady state machining behavior (force) was realized. Experimental cutting force data (cutting and thrust forces) were measured for a planning operation and compared to the calculations. 13 refs., 6 figs.
Nielsen, Jens; D’Avezac, Mayeul; Hetherington, James; Stamatakis, Michail
2013-12-14
Ab initio kinetic Monte Carlo (KMC) simulations have been successfully applied for over two decades to elucidate the underlying physico-chemical phenomena on the surfaces of heterogeneous catalysts. These simulations necessitate detailed knowledge of the kinetics of elementary reactions constituting the reaction mechanism, and the energetics of the species participating in the chemistry. The information about the energetics is encoded in the formation energies of gas and surface-bound species, and the lateral interactions between adsorbates on the catalytic surface, which can be modeled at different levels of detail. The majority of previous works accounted for only pairwise-additive first nearest-neighbor interactions. More recently, cluster-expansion Hamiltonians incorporating long-range interactions and many-body terms have been used for detailed estimations of catalytic rate [C. Wu, D. J. Schmidt, C. Wolverton, and W. F. Schneider, J. Catal. 286, 88 (2012)]. In view of the increasing interest in accurate predictions of catalytic performance, there is a need for general-purpose KMC approaches incorporating detailed cluster expansion models for the adlayer energetics. We have addressed this need by building on the previously introduced graph-theoretical KMC framework, and we have developed Zacros, a FORTRAN2003 KMC package for simulating catalytic chemistries. To tackle the high computational cost in the presence of long-range interactions we introduce parallelization with OpenMP. We further benchmark our framework by simulating a KMC analogue of the NO oxidation system established by Schneider and co-workers [J. Catal. 286, 88 (2012)]. We show that taking into account only first nearest-neighbor interactions may lead to large errors in the prediction of the catalytic rate, whereas for accurate estimates thereof, one needs to include long-range terms in the cluster expansion.
Numerical Simulation of the 2004 Indian Ocean Tsunami: Accurate Flooding and drying in Banda Aceh
NASA Astrophysics Data System (ADS)
Cui, Haiyang; Pietrzak, Julie; Stelling, Guus; Androsov, Alexey; Harig, Sven
2010-05-01
The Indian Ocean Tsunami on December 26, 2004 caused one of the largest tsunamis in recent times and led to widespread devastation and loss of life. One of the worst hit regions was Banda Aceh, which is the capital of the Aceh province, located in the northern part of Sumatra, 150km from the source of the earthquake. A German-Indonesian Tsunami Early Warning System (GITEWS) (www.gitews.de) is currently under active development. The work presented here is carried out within the GITEWS framework. One of the aims of this project is the development of accurate models with which to simulate the propagation, flooding and drying, and run-up of a tsunami. In this context, TsunAWI has been developed by the Alfred Wegener Institute; it is an explicit, () finite element model. However, the accurate numerical simulation of flooding and drying requires the conservation of mass and momentum. This is not possible in the current version of TsunAWi. The P1NC - P1element guarantees mass conservation in a global sense, yet as we show here it is important to guarantee mass conservation at the local level, that is within each individual cell. Here an unstructured grid, finite volume ocean model is presented. It is derived from the P1NC - P1 element, and is shown to be mass and momentum conserving. Then a number of simulations are presented, including dam break problems flooding over both a wet and a dry bed. Excellent agreement is found. Then we present simulations for Banda Aceh, and compare the results to on-site survey data, as well as to results from the original TsunAWI code.
De Vos, Maarten; De Lathauwer, Lieven; Vanrumste, Bart; Van Huffel, Sabine; Van Paesschen, W.
2007-01-01
Long-term electroencephalographic (EEG) recordings are important in the presurgical evaluation of refractory partial epilepsy for the delineation of the ictal onset zones. In this paper, we introduce a new concept for an automatic, fast, and objective localisation of the ictal onset zone in ictal EEG recordings. Canonical decomposition of ictal EEG decomposes the EEG in atoms. One or more atoms are related to the seizure activity. A single dipole was then fitted to model the potential distribution of each epileptic atom. In this study, we performed a simulation study in order to estimate the dipole localisation error. Ictal dipole localisation was very accurate, even at low signal-to-noise ratios, was not affected by seizure activity frequency or frequency changes, and was minimally affected by the waveform and depth of the ictal onset zone location. Ictal dipole localisation error using 21 electrodes was around 10.0 mm and improved more than tenfold in the range of 0.5–1.0 mm using 148 channels. In conclusion, our simulation study of canonical decomposition of ictal scalp EEG allowed a robust and accurate localisation of the ictal onset zone. PMID:18301715
Recommendations for accurate numerical blood flow simulations of stented intracranial aneurysms.
Janiga, Gábor; Berg, Philipp; Beuing, Oliver; Neugebauer, Mathias; Gasteiger, Rocco; Preim, Bernhard; Rose, Georg; Skalej, Martin; Thévenin, Dominique
2013-06-01
The number of scientific publications dealing with stented intracranial aneurysms is rapidly increasing. Powerful computational facilities are now available; an accurate computational modeling of hemodynamics in patient-specific configurations is, however, still being sought. Furthermore, there is still no general agreement on the quantities that should be computed and on the most adequate analysis for intervention support. In this article, the accurate representation of patient geometry is first discussed, involving successive improvements. Concerning the second step, the mesh required for the numerical simulation is especially challenging when deploying a stent with very fine wire structures. Third, the description of the fluid properties is a major challenge. Finally, a founded quantitative analysis of the simulation results is obviously needed to support interventional decisions. In the present work, an attempt has been made to review the most important steps for a high-quality computational fluid dynamics computation of virtually stented intracranial aneurysms. In consequence, this leads to concrete recommendations, whereby the obtained results are not discussed for their medical relevance but for the evaluation of their quality. This investigation might hopefully be helpful for further studies considering stent deployment in patient-specific geometries, in particular regarding the generation of the most appropriate computational model. PMID:23729530
NASA Technical Reports Server (NTRS)
Alter, Stephen J.; Brauckmann, Gregory J.; Kleb, William L.; Glass, Christopher E.; Streett, Craig L.; Schuster, David M.
2015-01-01
A transonic flow field about a Space Launch System (SLS) configuration was simulated with the Fully Unstructured Three-Dimensional (FUN3D) computational fluid dynamics (CFD) code at wind tunnel conditions. Unsteady, time-accurate computations were performed using second-order Delayed Detached Eddy Simulation (DDES) for up to 1.5 physical seconds. The surface pressure time history was collected at 619 locations, 169 of which matched locations on a 2.5 percent wind tunnel model that was tested in the 11 ft. x 11 ft. test section of the NASA Ames Research Center's Unitary Plan Wind Tunnel. Comparisons between computation and experiment showed that the peak surface pressure RMS level occurs behind the forward attach hardware, and good agreement for frequency and power was obtained in this region. Computational domain, grid resolution, and time step sensitivity studies were performed. These included an investigation of pseudo-time sub-iteration convergence. Using these sensitivity studies and experimental data comparisons, a set of best practices to date have been established for FUN3D simulations for SLS launch vehicle analysis. To the author's knowledge, this is the first time DDES has been used in a systematic approach and establish simulation time needed, to analyze unsteady pressure loads on a space launch vehicle such as the NASA SLS.
Oil Droplet Size Distribution and Optical Properties During Wave Tank Simulated Oil Spills
NASA Astrophysics Data System (ADS)
Conmy, R. N.; Venosa, A.; Courtenay, S.; King, T.; Robinson, B.; Ryan, S.
2013-12-01
Fate and transport of spilled petroleum oils in aquatic environments is highly dependent upon oil droplet behavior which is a function of chemical composition, dispersibility (natural and chemically-enhanced) and droplet size distribution (DSD) of the oil. DSD is influenced by mixing energy, temperature, salinity, pressure, presence of dissolved and particulate materials, flow rate of release, and application of dispersants. To better understand DSD and droplet behavior under varying physical conditions, flask-scale experiments are often insufficient. Rather, wave tank simulations allow for scaling to field conditions. Presented here are experiment results from the Bedford Institute of Oceanography wave tank facility, where chemically-dispersed (Corexit 9500; DOR = 1:20) Louisiana Sweet crude, IFO-120 and ANS crude oil were exposed to mixing energies to achieve dispersant effectiveness observed in the field. Oil plumes were simulated, both surface and subsea releases with varying water temperature and flow rate. Fluorometers (Chelsea Technologies Group AQUATracka, Turner Designs Cyclops, WET Labs Inc ECO) and particle size analyzers (Sequoia LISST) were used to track the dispersed plumes in the tank and characterize oil droplets. Sensors were validated with known oil volumes (down to 300 ppb) and measured Total Petroleum Hydrocarbons (TPH) and Benzene-Toluene-Ethylbenzene-Xylene (BTEX) values. This work has large implications for tracking surface and deep sea oil plumes with fluorescence and particle size analyzers, improved weathering and biodegradation estimates, and understanding the fate and transport of spill oil.
Unfitted Two-Phase Flow Simulations in Pore-Geometries with Accurate
NASA Astrophysics Data System (ADS)
Heimann, Felix; Engwer, Christian; Ippisch, Olaf; Bastian, Peter
2013-04-01
The development of better macro scale models for multi-phase flow in porous media is still impeded by the lack of suitable methods for the simulation of such flow regimes on the pore scale. The highly complicated geometry of natural porous media imposes requirements with regard to stability and computational efficiency which current numerical methods fail to meet. Therefore, current simulation environments are still unable to provide a thorough understanding of porous media in multi-phase regimes and still fail to reproduce well known effects like hysteresis or the more peculiar dynamics of the capillary fringe with satisfying accuracy. Although flow simulations in pore geometries were initially the domain of Lattice-Boltzmann and other particle methods, the development of Galerkin methods for such applications is important as they complement the range of feasible flow and parameter regimes. In the recent past, it has been shown that unfitted Galerkin methods can be applied efficiently to topologically demanding geometries. However, in the context of two-phase flows, the interface of the two immiscible fluids effectively separates the domain in two sub-domains. The exact representation of such setups with multiple independent and time depending geometries exceeds the functionality of common unfitted methods. We present a new approach to pore scale simulations with an unfitted discontinuous Galerkin (UDG) method. Utilizing a recursive sub-triangulation algorithm, we extent the UDG method to setups with multiple independent geometries. This approach allows an accurate representation of the moving contact line and the interface conditions, i.e. the pressure jump across the interface. Example simulations in two and three dimensions illustrate and verify the stability and accuracy of this approach.
Lee, M.W.; Meuwly, M.
2013-01-01
The evaluation of hydration free energies is a sensitive test to assess force fields used in atomistic simulations. We showed recently that the vibrational relaxation times, 1D- and 2D-infrared spectroscopies for CN(-) in water can be quantitatively described from molecular dynamics (MD) simulations with multipolar force fields and slightly enlarged van der Waals radii for the C- and N-atoms. To validate such an approach, the present work investigates the solvation free energy of cyanide in water using MD simulations with accurate multipolar electrostatics. It is found that larger van der Waals radii are indeed necessary to obtain results close to the experimental values when a multipolar force field is used. For CN(-), the van der Waals ranges refined in our previous work yield hydration free energy between -72.0 and -77.2 kcal mol(-1), which is in excellent agreement with the experimental data. In addition to the cyanide ion, we also study the hydroxide ion to show that the method used here is readily applicable to similar systems. Hydration free energies are found to sensitively depend on the intermolecular interactions, while bonded interactions are less important, as expected. We also investigate in the present work the possibility of applying the multipolar force field in scoring trajectories generated using computationally inexpensive methods, which should be useful in broader parametrization studies with reduced computational resources, as scoring is much faster than the generation of the trajectories.
Linaro, Daniele; Storace, Marco; Giugliano, Michele
2011-01-01
Stochastic channel gating is the major source of intrinsic neuronal noise whose functional consequences at the microcircuit- and network-levels have been only partly explored. A systematic study of this channel noise in large ensembles of biophysically detailed model neurons calls for the availability of fast numerical methods. In fact, exact techniques employ the microscopic simulation of the random opening and closing of individual ion channels, usually based on Markov models, whose computational loads are prohibitive for next generation massive computer models of the brain. In this work, we operatively define a procedure for translating any Markov model describing voltage- or ligand-gated membrane ion-conductances into an effective stochastic version, whose computer simulation is efficient, without compromising accuracy. Our approximation is based on an improved Langevin-like approach, which employs stochastic differential equations and no Montecarlo methods. As opposed to an earlier proposal recently debated in the literature, our approximation reproduces accurately the statistical properties of the exact microscopic simulations, under a variety of conditions, from spontaneous to evoked response features. In addition, our method is not restricted to the Hodgkin-Huxley sodium and potassium currents and is general for a variety of voltage- and ligand-gated ion currents. As a by-product, the analysis of the properties emerging in exact Markov schemes by standard probability calculus enables us for the first time to analytically identify the sources of inaccuracy of the previous proposal, while providing solid ground for its modification and improvement we present here. PMID:21423712
Lee, Myung Won; Meuwly, Markus
2013-12-14
The evaluation of hydration free energies is a sensitive test to assess force fields used in atomistic simulations. We showed recently that the vibrational relaxation times, 1D- and 2D-infrared spectroscopies for CN(-) in water can be quantitatively described from molecular dynamics (MD) simulations with multipolar force fields and slightly enlarged van der Waals radii for the C- and N-atoms. To validate such an approach, the present work investigates the solvation free energy of cyanide in water using MD simulations with accurate multipolar electrostatics. It is found that larger van der Waals radii are indeed necessary to obtain results close to the experimental values when a multipolar force field is used. For CN(-), the van der Waals ranges refined in our previous work yield hydration free energy between -72.0 and -77.2 kcal mol(-1), which is in excellent agreement with the experimental data. In addition to the cyanide ion, we also study the hydroxide ion to show that the method used here is readily applicable to similar systems. Hydration free energies are found to sensitively depend on the intermolecular interactions, while bonded interactions are less important, as expected. We also investigate in the present work the possibility of applying the multipolar force field in scoring trajectories generated using computationally inexpensive methods, which should be useful in broader parametrization studies with reduced computational resources, as scoring is much faster than the generation of the trajectories. PMID:24170171
Meek, Garrett A; Levine, Benjamin G
2014-07-01
Spikes in the time-derivative coupling (TDC) near surface crossings make the accurate integration of the time-dependent Schrödinger equation in nonadiabatic molecular dynamics simulations a challenge. To address this issue, we present an approximation to the TDC based on a norm-preserving interpolation (NPI) of the adiabatic electronic wave functions within each time step. We apply NPI and two other schemes for computing the TDC in numerical simulations of the Landau-Zener model, comparing the simulated transfer probabilities to the exact solution. Though NPI does not require the analytical calculation of nonadiabatic coupling matrix elements, it consistently yields unsigned population transfer probability errors of ∼0.001, whereas analytical calculation of the TDC yields errors of 0.0-1.0 depending on the time step, the offset of the maximum in the TDC from the beginning of the time step, and the coupling strength. The approximation of Hammes-Schiffer and Tully yields errors intermediate between NPI and the analytical scheme. PMID:26279558
Time Accurate CFD Simulations of the Orion Launch Abort Vehicle in the Transonic Regime
NASA Technical Reports Server (NTRS)
Rojahn, Josh; Ruf, Joe
2011-01-01
Significant asymmetries in the fluid dynamics were calculated for some cases in the CFD simulations of the Orion Launch Abort Vehicle through its abort trajectories. The CFD simulations were performed steady state and in three dimensions with symmetric geometries, no freestream sideslip angle, and motors firing. The trajectory points at issue were in the transonic regime, at 0 and +/- 5 angles of attack with the Abort Motors with and without the Attitude Control Motors (ACM) firing. In some of the cases the asymmetric fluid dynamics resulted in aerodynamic side forces that were large enough that would overcome the control authority of the ACMs. MSFC's Fluid Dynamics Group supported the investigation into the cause of the flow asymmetries with time accurate CFD simulations, utilizing a hybrid RANS-LES turbulence model. The results show that the flow over the vehicle and the subsequent interaction with the AB and ACM motor plumes were unsteady. The resulting instantaneous aerodynamic forces were oscillatory with fairly large magnitudes. Time averaged aerodynamic forces were essentially symmetric.
Cartesian Off-Body Grid Adaption for Viscous Time- Accurate Flow Simulation
NASA Technical Reports Server (NTRS)
Buning, Pieter G.; Pulliam, Thomas H.
2011-01-01
An improved solution adaption capability has been implemented in the OVERFLOW overset grid CFD code. Building on the Cartesian off-body approach inherent in OVERFLOW and the original adaptive refinement method developed by Meakin, the new scheme provides for automated creation of multiple levels of finer Cartesian grids. Refinement can be based on the undivided second-difference of the flow solution variables, or on a specific flow quantity such as vorticity. Coupled with load-balancing and an inmemory solution interpolation procedure, the adaption process provides very good performance for time-accurate simulations on parallel compute platforms. A method of using refined, thin body-fitted grids combined with adaption in the off-body grids is presented, which maximizes the part of the domain subject to adaption. Two- and three-dimensional examples are used to illustrate the effectiveness and performance of the adaption scheme.
NASA Technical Reports Server (NTRS)
Przekwas, A. J.; Athavale, M. M.; Hendricks, R. C.; Steinetz, B. M.
2006-01-01
Detailed information of the flow-fields in the secondary flowpaths and their interaction with the primary flows in gas turbine engines is necessary for successful designs with optimized secondary flow streams. Present work is focused on the development of a simulation methodology for coupled time-accurate solutions of the two flowpaths. The secondary flowstream is treated using SCISEAL, an unstructured adaptive Cartesian grid code developed for secondary flows and seals, while the mainpath flow is solved using TURBO, a density based code with capability of resolving rotor-stator interaction in multi-stage machines. An interface is being tested that links the two codes at the rim seal to allow data exchange between the two codes for parallel, coupled execution. A description of the coupling methodology and the current status of the interface development is presented. Representative steady-state solutions of the secondary flow in the UTRC HP Rig disc cavity are also presented.
Pacheco, P; Miller, P; Kim, J; Leese, T; Zabiyaka, Y
2003-05-07
Object-oriented NeuroSys (ooNeuroSys) is a collection of programs for simulating very large networks of biologically accurate neurons on distributed memory parallel computers. It includes two principle programs: ooNeuroSys, a parallel program for solving the large systems of ordinary differential equations arising from the interconnected neurons, and Neurondiz, a parallel program for visualizing the results of ooNeuroSys. Both programs are designed to be run on clusters and use the MPI library to obtain parallelism. ooNeuroSys also includes an easy-to-use Python interface. This interface allows neuroscientists to quickly develop and test complex neuron models. Both ooNeuroSys and Neurondiz have a design that allows for both high performance and relative ease of maintenance.
Time-Accurate Computational Fluid Dynamics Simulation of a Pair of Moving Solid Rocket Boosters
NASA Technical Reports Server (NTRS)
Strutzenberg, Louise L.; Williams, Brandon R.
2011-01-01
Since the Columbia accident, the threat to the Shuttle launch vehicle from debris during the liftoff timeframe has been assessed by the Liftoff Debris Team at NASA/MSFC. In addition to engineering methods of analysis, CFD-generated flow fields during the liftoff timeframe have been used in conjunction with 3-DOF debris transport methods to predict the motion of liftoff debris. Early models made use of a quasi-steady flow field approximation with the vehicle positioned at a fixed location relative to the ground; however, a moving overset mesh capability has recently been developed for the Loci/CHEM CFD software which enables higher-fidelity simulation of the Shuttle transient plume startup and liftoff environment. The present work details the simulation of the launch pad and mobile launch platform (MLP) with truncated solid rocket boosters (SRBs) moving in a prescribed liftoff trajectory derived from Shuttle flight measurements. Using Loci/CHEM, time-accurate RANS and hybrid RANS/LES simulations were performed for the timeframe T0+0 to T0+3.5 seconds, which consists of SRB startup to a vehicle altitude of approximately 90 feet above the MLP. Analysis of the transient flowfield focuses on the evolution of the SRB plumes in the MLP plume holes and the flame trench, impingement on the flame deflector, and especially impingment on the MLP deck resulting in upward flow which is a transport mechanism for debris. The results show excellent qualitative agreement with the visual record from past Shuttle flights, and comparisons to pressure measurements in the flame trench and on the MLP provide confidence in these simulation capabilities.
Achieving accurate simulations of urban impacts on ozone at high resolution
NASA Astrophysics Data System (ADS)
Li, J.; Georgescu, M.; Hyde, P.; Mahalov, A.; Moustaoui, M.
2014-11-01
The effects of urbanization on ozone levels have been widely investigated over cities primarily located in temperate and/or humid regions. In this study, nested WRF-Chem simulations with a finest grid resolution of 1 km are conducted to investigate ozone concentrations [O3] due to urbanization within cities in arid/semi-arid environments. First, a method based on a shape preserving Monotonic Cubic Interpolation (MCI) is developed and used to downscale anthropogenic emissions from the 4 km resolution 2005 National Emissions Inventory (NEI05) to the finest model resolution of 1 km. Using the rapidly expanding Phoenix metropolitan region as the area of focus, we demonstrate the proposed MCI method achieves ozone simulation results with appreciably improved correspondence to observations relative to the default interpolation method of the WRF-Chem system. Next, two additional sets of experiments are conducted, with the recommended MCI approach, to examine impacts of urbanization on ozone production: (1) the urban land cover is included (i.e., urbanization experiments) and, (2) the urban land cover is replaced with the region’s native shrubland. Impacts due to the presence of the built environment on [O3] are highly heterogeneous across the metropolitan area. Increased near surface [O3] due to urbanization of 10-20 ppb is predominantly a nighttime phenomenon while simulated impacts during daytime are negligible. Urbanization narrows the daily [O3] range (by virtue of increasing nighttime minima), an impact largely due to the region’s urban heat island. Our results demonstrate the importance of the MCI method for accurate representation of the diurnal profile of ozone, and highlight its utility for high-resolution air quality simulations for urban areas.
At the request of the US EPA Oil Program Center, ERD is developing an oil spill model that focuses on fate and transport of oil components under various response scenarios. This model includes various simulation options, including the use of chemical dispersing agents on oil sli...
Finite domain simulations with adaptive boundaries: accurate potentials and nonequilibrium movesets.
Wagoner, Jason A; Pande, Vijay S
2013-12-21
We extend the theory of hybrid explicit/implicit solvent models to include an explicit domain that grows and shrinks in response to a solute's evolving configuration. The goal of this model is to provide an appropriate but not excessive amount of solvent detail, and the inclusion of an adjustable boundary provides a significant computational advantage for solutes that explore a range of configurations. In addition to the theoretical development, a successful implementation of this method requires (1) an efficient moveset that propagates the boundary as a new coordinate of the system, and (2) an accurate continuum solvent model with parameters that are transferable to an explicit domain of any size. We address these challenges and develop boundary updates using Monte Carlo moves biased by nonequilibrium paths. We obtain the desired level of accuracy using a "decoupling interface" that we have previously shown to remove boundary artifacts common to hybrid solvent models. Using an uncharged, coarse-grained solvent model, we then study the efficiency of nonequilibrium paths that a simulation takes by quantifying the dissipation. In the spirit of optimization, we study this quantity over a range of simulation parameters. PMID:24359359
An accurate and efficient 3-D micromagnetic simulation of metal evaporated tape
NASA Astrophysics Data System (ADS)
Jones, M.; Miles, J. J.
1997-07-01
Metal evaporated tape (MET) has a complex column-like structure in which magnetic domains are arranged randomly. In order to accurately simulate the behaviour of MET it is important to capture these aspects of the material in a high-resolution 3-D micromagnetic model. The scale of this problem prohibits the use of traditional scalar computers and leads us to develop algorithms for a vector processor architecture. We demonstrate that despite the materials highly non-uniform structure, it is possible to develop fast vector algorithms for the computation of the magnetostatic interaction field. We do this by splitting the field calculation into near and far components. The near field component is calculated exactly using an efficient vector algorithm, whereas the far field is calculated approximately using a novel fast Fourier transform (FFT) technique. Results are presented which demonstrate that, in practice, the algorithms require sub-O( N log( N)) computation time. In addition results of highly realistic simulation of hysteresis in MET are presented.
How to obtain accurate resist simulations in very low-k1 era?
NASA Astrophysics Data System (ADS)
Chiou, Tsann-Bim; Park, Chan-Ha; Choi, Jae-Seung; Min, Young-Hong; Hansen, Steve; Tseng, Shih-En; Chen, Alek C.; Yim, Donggyu
2006-03-01
A procedure for calibrating a resist model iteratively adjusts appropriate parameters until the simulations of the model match the experimental data. The tunable parameters may include the shape of the illuminator, the geometry and transmittance/phase of the mask, light source and scanner-related parameters that affect imaging quality, resist process control and most importantly the physical/chemical factors in the resist model. The resist model can be accurately calibrated by measuring critical dimensions (CD) of a focus-exposure matrix (FEM) and the technique has been demonstrated to be very successful in predicting lithographic performance. However, resist model calibration is more challenging in the low k1 (<0.3) regime because numerous uncertainties, such as mask and resist CD metrology errors, are becoming too large to be ignored. This study demonstrates a resist model calibration procedure for a 0.29 k1 process using a 6% halftone mask containing 2D brickwall patterns. The influence of different scanning electron microscopes (SEM) and their wafer metrology signal analysis algorithms on the accuracy of the resist model is evaluated. As an example of the metrology issue of the resist pattern, the treatment of a sidewall angle is demonstrated for the resist line ends where the contrast is relatively low. Additionally, the mask optical proximity correction (OPC) and corner rounding are considered in the calibration procedure that is based on captured SEM images. Accordingly, the average root-mean-square (RMS) error, which is the difference between simulated and experimental CDs, can be improved by considering the metrological issues. Moreover, a weighting method and a measured CD tolerance are proposed to handle the different CD variations of the various edge points of the wafer resist pattern. After the weighting method is implemented and the CD selection criteria applied, the RMS error can be further suppressed. Therefore, the resist CD and process window can
Numerical simulation of oil pool boundary evolution
NASA Astrophysics Data System (ADS)
Khudobina, Yulia; Bubenchikov, Aleksey; Bubenchikov, Mikhail; Matvienko, Oleg; Libin, Eduard
2016-01-01
The study of spatial distribution of hydrocarbon resources and forecasting their geographical location is of great importance for the most complete recovery of hydrocarbons from deposits. The present study gives new mathematical results in the theory of stratified fluid flow in a porous medium. This paper analyzes the evolution of oil pool boundary basing on vortex numerical model for movement of the boundary separating fluids of different densities. It presents the investigation of how the location of light fluid regarding the heavier fluid influences on the changes in the boundary between two media in case of various shifting of the well.
Stochastic simulation model of oil spill fate and transfer
Al-Rabeh, A.H.; Cekirge, H.M.; Gunay, N. )
1989-06-01
Over the past few years, considerable research has been directed toward the development of mathematical models to describe the behavior of oil spills. A successful model would be of great value in selecting locations for the deployment of containment and collection systems to mitigate the effects of the pollutant on the environment. In this study, a comprehensive stochastic model is formulated to simulate the fate and transport of oil spills. The model consists of a set of algorithms describing the processes of advection, turbulent diffusion, surface spreading, vertical mechanical dispersion, emulsification, and evaporation. Each algorithm is developed separately and is linked to related processes and to environmental and other parameters. The model requires as input the velocity field of the transporting medium. This can be obtained from a three-dimensional hydrodynamic model for tidal and wind-driven currents for the region of interest. The oil spill fate and transport model is used to simulate a surface oil spill in the Abu Ali region on the western side of the Arabian Gulf. The simulation results indicate that the model can predict the fate and transport of oil slicks with reasonable accuracy. 45 refs., 10 figs.
Experiences with linear solvers for oil reservoir simulation problems
Joubert, W.; Janardhan, R.; Biswas, D.; Carey, G.
1996-12-31
This talk will focus on practical experiences with iterative linear solver algorithms used in conjunction with Amoco Production Company`s Falcon oil reservoir simulation code. The goal of this study is to determine the best linear solver algorithms for these types of problems. The results of numerical experiments will be presented.
A Three Dimensional Parallel Time Accurate Turbopump Simulation Procedure Using Overset Grid Systems
NASA Technical Reports Server (NTRS)
Kiris, Cetin; Chan, William; Kwak, Dochan
2001-01-01
The objective of the current effort is to provide a computational framework for design and analysis of the entire fuel supply system of a liquid rocket engine, including high-fidelity unsteady turbopump flow analysis. This capability is needed to support the design of pump sub-systems for advanced space transportation vehicles that are likely to involve liquid propulsion systems. To date, computational tools for design/analysis of turbopump flows are based on relatively lower fidelity methods. An unsteady, three-dimensional viscous flow analysis tool involving stationary and rotational components for the entire turbopump assembly has not been available for real-world engineering applications. The present effort provides developers with information such as transient flow phenomena at start up, and non-uniform inflows, and will eventually impact on system vibration and structures. In the proposed paper, the progress toward the capability of complete simulation of the turbo-pump for a liquid rocket engine is reported. The Space Shuttle Main Engine (SSME) turbo-pump is used as a test case for evaluation of the hybrid MPI/Open-MP and MLP versions of the INS3D code. CAD to solution auto-scripting capability is being developed for turbopump applications. The relative motion of the grid systems for the rotor-stator interaction was obtained using overset grid techniques. Unsteady computations for the SSME turbo-pump, which contains 114 zones with 34.5 million grid points, are carried out on Origin 3000 systems at NASA Ames Research Center. Results from these time-accurate simulations with moving boundary capability will be presented along with the performance of parallel versions of the code.
A Three-Dimensional Parallel Time-Accurate Turbopump Simulation Procedure Using Overset Grid System
NASA Technical Reports Server (NTRS)
Kiris, Cetin; Chan, William; Kwak, Dochan
2002-01-01
The objective of the current effort is to provide a computational framework for design and analysis of the entire fuel supply system of a liquid rocket engine, including high-fidelity unsteady turbopump flow analysis. This capability is needed to support the design of pump sub-systems for advanced space transportation vehicles that are likely to involve liquid propulsion systems. To date, computational tools for design/analysis of turbopump flows are based on relatively lower fidelity methods. An unsteady, three-dimensional viscous flow analysis tool involving stationary and rotational components for the entire turbopump assembly has not been available for real-world engineering applications. The present effort provides developers with information such as transient flow phenomena at start up, and nonuniform inflows, and will eventually impact on system vibration and structures. In the proposed paper, the progress toward the capability of complete simulation of the turbo-pump for a liquid rocket engine is reported. The Space Shuttle Main Engine (SSME) turbo-pump is used as a test case for evaluation of the hybrid MPI/Open-MP and MLP versions of the INS3D code. CAD to solution auto-scripting capability is being developed for turbopump applications. The relative motion of the grid systems for the rotor-stator interaction was obtained using overset grid techniques. Unsteady computations for the SSME turbo-pump, which contains 114 zones with 34.5 million grid points, are carried out on Origin 3000 systems at NASA Ames Research Center. Results from these time-accurate simulations with moving boundary capability are presented along with the performance of parallel versions of the code.
Romanov, V N; Cygan, R T; Myshakin, E M
2012-06-21
Naturally occurring clay minerals provide a distinctive material for carbon capture and carbon dioxide sequestration. Swelling clay minerals, such as the smectite variety, possess an aluminosilicate structure that is controlled by low-charge layers that readily expand to accommodate water molecules and, potentially, CO2. Recent experimental studies have demonstrated the efficacy of intercalating CO2 in the interlayer of layered clays, but little is known about the molecular mechanisms of the process and the extent of carbon capture as a function of clay charge and structure. A series of molecular dynamics simulations and vibrational analyses have been completed to assess the molecular interactions associated with incorporation of CO2 and H2O in the interlayer of montmorillonite clay and to help validate the models with experimental observation. An accurate and fully flexible set of interatomic potentials for CO2 is developed and combined with Clayff potentials to help evaluate the intercalation mechanism and examine the effect of molecular flexibility onthe diffusion rate of CO2 in water.
Candel, A.; Kabel, A.; Lee, L.; Li, Z.; Limborg, C.; Ng, C.; Prudencio, E.; Schussman, G.; Uplenchwar, R.; Ko, K.; /SLAC
2009-06-19
Over the past years, SLAC's Advanced Computations Department (ACD), under SciDAC sponsorship, has developed a suite of 3D (2D) parallel higher-order finite element (FE) codes, T3P (T2P) and Pic3P (Pic2P), aimed at accurate, large-scale simulation of wakefields and particle-field interactions in radio-frequency (RF) cavities of complex shape. The codes are built on the FE infrastructure that supports SLAC's frequency domain codes, Omega3P and S3P, to utilize conformal tetrahedral (triangular)meshes, higher-order basis functions and quadratic geometry approximation. For time integration, they adopt an unconditionally stable implicit scheme. Pic3P (Pic2P) extends T3P (T2P) to treat charged-particle dynamics self-consistently using the PIC (particle-in-cell) approach, the first such implementation on a conformal, unstructured grid using Whitney basis functions. Examples from applications to the International Linear Collider (ILC), Positron Electron Project-II (PEP-II), Linac Coherent Light Source (LCLS) and other accelerators will be presented to compare the accuracy and computational efficiency of these codes versus their counterparts using structured grids.
Joldes, Grand Roman; Wittek, Adam; Miller, Karol
2008-01-01
Real time computation of soft tissue deformation is important for the use of augmented reality devices and for providing haptic feedback during operation or surgeon training. This requires algorithms that are fast, accurate and can handle material nonlinearities and large deformations. A set of such algorithms is presented in this paper, starting with the finite element formulation and the integration scheme used and addressing common problems such as hourglass control and locking. The computation examples presented prove that by using these algorithms, real time computations become possible without sacrificing the accuracy of the results. For a brain model having more than 7000 degrees of freedom, we computed the reaction forces due to indentation with frequency of around 1000 Hz using a standard dual core PC. Similarly, we conducted simulation of brain shift using a model with more than 50 000 degrees of freedom in less than a minute. The speed benefits of our models results from combining the Total Lagrangian formulation with explicit time integration and low order finite elements. PMID:19152791
A mechanistic approach for accurate simulation of village scale malaria transmission
Bomblies, Arne; Duchemin, Jean-Bernard; Eltahir, Elfatih AB
2009-01-01
Background Malaria transmission models commonly incorporate spatial environmental and climate variability for making regional predictions of disease risk. However, a mismatch of these models' typical spatial resolutions and the characteristic scale of malaria vector population dynamics may confound disease risk predictions in areas of high spatial hydrological variability such as the Sahel region of Africa. Methods Field observations spanning two years from two Niger villages are compared. The two villages are separated by only 30 km but exhibit a ten-fold difference in anopheles mosquito density. These two villages would be covered by a single grid cell in many malaria models, yet their entomological activity differs greatly. Environmental conditions and associated entomological activity are simulated at high spatial- and temporal resolution using a mechanistic approach that couples a distributed hydrology scheme and an entomological model. Model results are compared to regular field observations of Anopheles gambiae sensu lato mosquito populations and local hydrology. The model resolves the formation and persistence of individual pools that facilitate mosquito breeding and predicts spatio-temporal mosquito population variability at high resolution using an agent-based modeling approach. Results Observations of soil moisture, pool size, and pool persistence are reproduced by the model. The resulting breeding of mosquitoes in the simulated pools yields time-integrated seasonal mosquito population dynamics that closely follow observations from captured mosquito abundance. Interannual difference in mosquito abundance is simulated, and the inter-village difference in mosquito population is reproduced for two years of observations. These modeling results emulate the known focal nature of malaria in Niger Sahel villages. Conclusion Hydrological variability must be represented at high spatial and temporal resolution to achieve accurate predictive ability of malaria risk
Simulation of oil pollution in the Persian Gulf near Assaluyeh oil terminal.
Faghihifard, M; Badri, M A
2016-04-15
Numerical simulation of oil slick movement with respect to tidal factors and wind effects was performed in order to counteract oil pollution in the Persian Gulf. First, a flow model was invoked with respect to water level fluctuations. The main tidal constituents were applied to the model using the initial conditions of water level variations in the Hormuz Strait near the Hangam Island. The movement of oil pollution was determined due to wind, tide and temperature effects and confirmed by applying a verified field results. Simulations were focused near an important terminal in the Persian Gulf, Assaluyeh Port. The results were led to preparing a risk-taking map in a parallel research for the Persian Gulf. PMID:26906497
NASA Astrophysics Data System (ADS)
Shukla, Ratnesh K.
2014-11-01
Single fluid schemes that rely on an interface function for phase identification in multicomponent compressible flows are widely used to study hydrodynamic flow phenomena in several diverse applications. Simulations based on standard numerical implementation of these schemes suffer from an artificial increase in the width of the interface function owing to the numerical dissipation introduced by an upwind discretization of the governing equations. In addition, monotonicity requirements which ensure that the sharp interface function remains bounded at all times necessitate use of low-order accurate discretization strategies. This results in a significant reduction in accuracy along with a loss of intricate flow features. In this paper we develop a nonlinear transformation based interface capturing method which achieves superior accuracy without compromising the simplicity, computational efficiency and robustness of the original flow solver. A nonlinear map from the signed distance function to the sigmoid type interface function is used to effectively couple a standard single fluid shock and interface capturing scheme with a high-order accurate constrained level set reinitialization method in a way that allows for oscillation-free transport of the sharp material interface. Imposition of a maximum principle, which ensures that the multidimensional preconditioned interface capturing method does not produce new maxima or minima even in the extreme events of interface merger or breakup, allows for an explicit determination of the interface thickness in terms of the grid spacing. A narrow band method is formulated in order to localize computations pertinent to the preconditioned interface capturing method. Numerical tests in one dimension reveal a significant improvement in accuracy and convergence; in stark contrast to the conventional scheme, the proposed method retains its accuracy and convergence characteristics in a shifted reference frame. Results from the test
A time-accurate adaptive grid method and the numerical simulation of a shock-vortex interaction
NASA Technical Reports Server (NTRS)
Bockelie, Michael J.; Eiseman, Peter R.
1990-01-01
A time accurate, general purpose, adaptive grid method is developed that is suitable for multidimensional steady and unsteady numerical simulations. The grid point movement is performed in a manner that generates smooth grids which resolve the severe solution gradients and the sharp transitions in the solution gradients. The temporal coupling of the adaptive grid and the PDE solver is performed with a grid prediction correction method that is simple to implement and ensures the time accuracy of the grid. Time accurate solutions of the 2-D Euler equations for an unsteady shock vortex interaction demonstrate the ability of the adaptive method to accurately adapt the grid to multiple solution features.
NASA Astrophysics Data System (ADS)
Wosnik, M.; Bachant, P.
2014-12-01
Cross-flow turbines, often referred to as vertical-axis turbines, show potential for success in marine hydrokinetic (MHK) and wind energy applications, ranging from small- to utility-scale installations in tidal/ocean currents and offshore wind. As turbine designs mature, the research focus is shifting from individual devices to the optimization of turbine arrays. It would be expensive and time-consuming to conduct physical model studies of large arrays at large model scales (to achieve sufficiently high Reynolds numbers), and hence numerical techniques are generally better suited to explore the array design parameter space. However, since the computing power available today is not sufficient to conduct simulations of the flow in and around large arrays of turbines with fully resolved turbine geometries (e.g., grid resolution into the viscous sublayer on turbine blades), the turbines' interaction with the energy resource (water current or wind) needs to be parameterized, or modeled. Models used today--a common model is the actuator disk concept--are not able to predict the unique wake structure generated by cross-flow turbines. This wake structure has been shown to create "constructive" interference in some cases, improving turbine performance in array configurations, in contrast with axial-flow, or horizontal axis devices. Towards a more accurate parameterization of cross-flow turbines, an extensive experimental study was carried out using a high-resolution turbine test bed with wake measurement capability in a large cross-section tow tank. The experimental results were then "interpolated" using high-fidelity Navier--Stokes simulations, to gain insight into the turbine's near-wake. The study was designed to achieve sufficiently high Reynolds numbers for the results to be Reynolds number independent with respect to turbine performance and wake statistics, such that they can be reliably extrapolated to full scale and used for model validation. The end product of
Thermal numerical simulator for laboratory evaluation of steamflood oil recovery
Sarathi, P.
1991-04-01
A thermal numerical simulator running on an IBM AT compatible personal computer is described. The simulator was designed to assist laboratory design and evaluation of steamflood oil recovery. An overview of the historical evolution of numerical thermal simulation, NIPER's approach to solving these problems with a desk top computer, the derivation of equations and a description of approaches used to solve these equations, and verification of the simulator using published data sets and sensitivity analysis are presented. The developed model is a three-phase, two-dimensional multicomponent simulator capable of being run in one or two dimensions. Mass transfer among the phases and components is dictated by pressure- and temperature-dependent vapor-liquid equilibria. Gravity and capillary pressure phenomena were included. Energy is transferred by conduction, convection, vaporization and condensation. The model employs a block centered grid system with a five-point discretization scheme. Both areal and vertical cross-sectional simulations are possible. A sequential solution technique is employed to solve the finite difference equations. The study clearly indicated the importance of heat loss, injected steam quality, and injection rate to the process. Dependence of overall recovery on oil volatility and viscosity is emphasized. The process is very sensitive to relative permeability values. Time-step sensitivity runs indicted that the current version is time-step sensitive and exhibits conditional stability. 75 refs., 19 figs., 19 tabs.
Conmy, Robyn N; Coble, Paula G; Farr, James; Wood, A Michelle; Lee, Kenneth; Pegau, W Scott; Walsh, Ian D; Koch, Corey R; Abercrombie, Mary I; Miles, M Scott; Lewis, Marlon R; Ryan, Scott A; Robinson, Brian J; King, Thomas L; Kelble, Christopher R; Lacoste, Jordanna
2014-01-01
In situ fluorometers were deployed during the Deepwater Horizon (DWH) Gulf of Mexico oil spill to track the subsea oil plume. Uncertainties regarding instrument specifications and capabilities necessitated performance testing of sensors exposed to simulated, dispersed oil plumes. Dynamic ranges of the Chelsea Technologies Group AQUAtracka, Turner Designs Cyclops, Satlantic SUNA and WET Labs, Inc. ECO, exposed to fresh and artificially weathered crude oil, were determined. Sensors were standardized against known oil volumes and total petroleum hydrocarbons and benzene-toluene-ethylbenzene-xylene measurements-both collected during spills, providing oil estimates during wave tank dilution experiments. All sensors estimated oil concentrations down to 300 ppb oil, refuting previous reports. Sensor performance results assist interpretation of DWH oil spill data and formulating future protocols. PMID:24377909
From oil shortage to oil glut: simulation of growth prospects in the Nigerian economy
Olofin, S.; Iyaniwura, J.O.
1983-11-01
During the 1970s, the economy of Nigeria provided one of the most interesting cases of development financed through oil revenue. Between 1970 and 1980, the country's GNP grew at an outstanding rate, but after the transition from oil shortage to oil glut, the economy of Nigeria ran into dramatic financial difficulties, which are now placing major constraints to its development. To investigate the transition from an oil-based economy to a stage characterized by greater diversification of exports and more balanced sectoral growth, a model has been built by the University of Ibadan in Nigeria. The model was developed in association with Project LINK staffing for the future inclusion in the Project. According to the finding presented in the study, the annual growth rate of GDP of Nigeria between 1980 and 1988 will be around 2.5%. To compensate the drop of the foreign-exchange earnings caused by the contraction of oil prices and demand, a vigorous export drive of agricultural products is simulated. 8 references, 7 figures, 4 tables.
NASA Astrophysics Data System (ADS)
Sellers, Michael; Lisal, Martin; Brennan, John
2015-06-01
Investigating the ability of a molecular model to accurately represent a real material is crucial to model development and use. When the model simulates materials in extreme conditions, one such property worth evaluating is the phase transition point. However, phase transitions are often overlooked or approximated because of difficulty or inaccuracy when simulating them. Techniques such as super-heating or super-squeezing a material to induce a phase change suffer from inherent timescale limitations leading to ``over-driving,'' and dual-phase simulations require many long-time runs to seek out what frequently results in an inexact location of phase-coexistence. We present a compilation of methods for the determination of solid-solid and solid-liquid phase transition points through the accurate calculation of the chemical potential. The methods are applied to the Smith-Bharadwaj atomistic potential's representation of cyclotrimethylene trinitramine (RDX) to accurately determine its melting point (Tm) and the alpha to gamma solid phase transition pressure. We also determine Tm for a coarse-grain model of RDX, and compare its value to experiment and atomistic counterpart. All methods are employed via the LAMMPS simulator, resulting in 60-70 simulations that total 30-50 ns. Approved for public release. Distribution is unlimited.
NASA Astrophysics Data System (ADS)
Moore, Christopher; Hopkins, Matthew; Moore, Stan; Boerner, Jeremiah; Cartwright, Keith
2015-09-01
Simulation of breakdown is important for understanding and designing a variety of applications such as mitigating undesirable discharge events. Such simulations need to be accurate through early time arc initiation to late time stable arc behavior. Here we examine constraints on the timestep and mesh size required for arc simulations using the particle-in-cell (PIC) method with direct simulation Monte Carlo (DMSC) collisions. Accurate simulation of electron avalanche across a fixed voltage drop and constant neutral density (reduced field of 1000 Td) was found to require a timestep ~ 1/100 of the mean time between collisions and a mesh size ~ 1/25 the mean free path. These constraints are much smaller than the typical PIC-DSMC requirements for timestep and mesh size. Both constraints are related to the fact that charged particles are accelerated by the external field. Thus gradients in the electron energy distribution function can exist at scales smaller than the mean free path and these must be resolved by the mesh size for accurate collision rates. Additionally, the timestep must be small enough that the particle energy change due to the fields be small in order to capture gradients in the cross sections versus energy. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. DOE's National Nuclear Security Administration under Contract DE-AC04-94AL85000.
NASA Astrophysics Data System (ADS)
Grasso, Robert J.; Russo, Leonard P.; Barrett, John L.; Odhner, Jefferson E.; Egbert, Paul I.
2007-09-01
BAE Systems presents the results of a program to model the performance of Raman LIDAR systems for the remote detection of atmospheric gases, air polluting hydrocarbons, chemical and biological weapons, and other molecular species of interest. Our model, which integrates remote Raman spectroscopy, 2D and 3D LADAR, and USAF atmospheric propagation codes permits accurate determination of the performance of a Raman LIDAR system. The very high predictive performance accuracy of our model is due to the very accurate calculation of the differential scattering cross section for the specie of interest at user selected wavelengths. We show excellent correlation of our calculated cross section data, used in our model, with experimental data obtained from both laboratory measurements and the published literature. In addition, the use of standard USAF atmospheric models provides very accurate determination of the atmospheric extinction at both the excitation and Raman shifted wavelengths.
Simulating Oil Production from Fractured/Faulted Basement Reservoirs
NASA Astrophysics Data System (ADS)
Wang, H.; Forster, C.; Fu, Y.; Huang, C.; Yang, Y.; Deo, M.
2006-12-01
A fully-implicit, three-dimensional (3D), three-phase, discrete fault/fracture, black oil simulator provides new insight and understanding of oil production from reservoirs in fractured, low-permeability basement rocks. Results obtained with a controlled volume finite element (CVFE) method compare favorably to those obtained using both single- and dual-porosity finite difference methods (e.g., ECLIPSE). A regularized network of 30 orthogonal faults within a 1000 by 1000 by 200 feet model domain is used to compare the simulator results and to explore the implications of grid sensitivity. In this simple reservoir, cumulative oil recoveries over 900 days of production are similar for CVFE, single-porosity and dual-porosity approaches. CVFE is used to simulate a complex network of intersecting faults that mimic a more realistic basement reservoir with the same fault surface area and fault volume as the regularized network. Cumulative oil production at 900 days is about 3% lower than for the regularized network. The CVFE method provides a much improved ability to represent complex fracture/fault geometries and spatial variations in the internal properties of faults. CVFE simulations of the realistic network illustrate the possible consequences of uncertainty in knowing fracture/fault properties (e.g., porosity, permeability, thickness, dip orientation, connectivity and flow transmissibility). For example, introducing spatial variability in permeability within the fault planes (using spatially randomized patterns of 10, 100 and 1000 md), while retaining a constant geometric mean permeability of 100 md, yields enhanced oil production due to the high-permeability pathways. A 50:50 mix of 10 and 1000 md elements yields 36%OOIP while a 33:33:33 mix of 10, 100 and 1000 md yields 24%OOIP. These results are 26% and 14% greater, respectively, than that obtained for the uniform 100 md case (11%OOIP). This inherent variability, combined with uncertainty in knowing the detailed
Accurate simulation of the electron cloud in the Fermilab Main Injector with VORPAL
Lebrun, Paul L.G.; Spentzouris, Panagiotis; Cary, John R.; Stoltz, Peter; Veitzer, Seth A.; /Tech-X, Boulder
2010-05-01
Precision simulations of the electron cloud at the Fermilab Main Injector have been studied using the plasma simulation code VORPAL. Fully 3D and self consistent solutions that includes E.M. field maps generated by the cloud and the proton bunches have been obtained, as well detailed distributions of the electron's 6D phase space. We plan to include such maps in the ongoing simulation of the space charge effects in the Main Injector. Simulations of the response of beam position monitors, retarding field analyzers and microwave transmission experiments are ongoing.
A CUDA based parallel multi-phase oil reservoir simulator
NASA Astrophysics Data System (ADS)
Zaza, Ayham; Awotunde, Abeeb A.; Fairag, Faisal A.; Al-Mouhamed, Mayez A.
2016-09-01
Forward Reservoir Simulation (FRS) is a challenging process that models fluid flow and mass transfer in porous media to draw conclusions about the behavior of certain flow variables and well responses. Besides the operational cost associated with matrix assembly, FRS repeatedly solves huge and computationally expensive sparse, ill-conditioned and unsymmetrical linear system. Moreover, as the computation for practical reservoir dimensions lasts for long times, speeding up the process by taking advantage of parallel platforms is indispensable. By considering the state of art advances in massively parallel computing and the accompanying parallel architecture, this work aims primarily at developing a CUDA-based parallel simulator for oil reservoir. In addition to the initial reported 33 times speed gain compared to the serial version, running experiments showed that BiCGSTAB is a stable and fast solver which could be incorporated in such simulations instead of the more expensive, storage demanding and usually utilized GMRES.
Advanced material testing in support of accurate sheet metal forming simulations
NASA Astrophysics Data System (ADS)
Kuwabara, Toshihiko
2013-05-01
This presentation is a review of experimental methods for accurately measuring and modeling the anisotropic plastic deformation behavior of metal sheets under a variety of loading paths: biaxial compression test, hydraulic bulge test, biaxial tension test using a cruciform specimen, multiaxial tube expansion test using a closed-loop electrohydraulic testing machine for the measurement of forming limit strains and stresses, combined tension-shear test, and in-plane stress reversal test. Observed material responses are compared with predictions using phenomenological plasticity models to highlight the importance of accurate material testing. Special attention is paid to the plastic deformation behavior of sheet metals commonly used in industry, and to verifying the validity of constitutive models based on anisotropic yield functions at a large plastic strain range. The effects of using appropriate material models on the improvement of predictive accuracy for forming defects, such as springback and fracture, are also presented.
A Variable Coefficient Method for Accurate Monte Carlo Simulation of Dynamic Asset Price
NASA Astrophysics Data System (ADS)
Li, Yiming; Hung, Chih-Young; Yu, Shao-Ming; Chiang, Su-Yun; Chiang, Yi-Hui; Cheng, Hui-Wen
2007-07-01
In this work, we propose an adaptive Monte Carlo (MC) simulation technique to compute the sample paths for the dynamical asset price. In contrast to conventional MC simulation with constant drift and volatility (μ,σ), our MC simulation is performed with variable coefficient methods for (μ,σ) in the solution scheme, where the explored dynamic asset pricing model starts from the formulation of geometric Brownian motion. With the method of simultaneously updated (μ,σ), more than 5,000 runs of MC simulation are performed to fulfills basic accuracy of the large-scale computation and suppresses statistical variance. Daily changes of stock market index in Taiwan and Japan are investigated and analyzed.
NASA Technical Reports Server (NTRS)
Venkatachari, Balaji Shankar; Streett, Craig L.; Chang, Chau-Lyan; Friedlander, David J.; Wang, Xiao-Yen; Chang, Sin-Chung
2016-01-01
Despite decades of development of unstructured mesh methods, high-fidelity time-accurate simulations are still predominantly carried out on structured, or unstructured hexahedral meshes by using high-order finite-difference, weighted essentially non-oscillatory (WENO), or hybrid schemes formed by their combinations. In this work, the space-time conservation element solution element (CESE) method is used to simulate several flow problems including supersonic jet/shock interaction and its impact on launch vehicle acoustics, and direct numerical simulations of turbulent flows using tetrahedral meshes. This paper provides a status report for the continuing development of the space-time conservation element solution element (CESE) numerical and software framework under the Revolutionary Computational Aerosciences (RCA) project. Solution accuracy and large-scale parallel performance of the numerical framework is assessed with the goal of providing a viable paradigm for future high-fidelity flow physics simulations.
NASA Technical Reports Server (NTRS)
Vatsa, Veer N.; Singer, Bart A.
2003-01-01
We evaluate the applicability of a production computational fluid dynamics code for conducting detached eddy simulation for unsteady flows. A second-order accurate Navier-Stokes code developed at NASA Langley Research Center, known as TLNS3D, is used for these simulations. We focus our attention on high Reynolds number flow (Re = 5 x 10(sup 4) - 1.4 x 10(sup 5)) past a circular cylinder to simulate flows with large-scale separations. We consider two types of flow situations: one in which the flow at the separation point is laminar, and the other in which the flow is already turbulent when it detaches from the surface of the cylinder. Solutions are presented for two- and three-dimensional calculations using both the unsteady Reynolds-averaged Navier-Stokes paradigm and the detached eddy simulation treatment. All calculations use the standard Spalart-Allmaras turbulence model as the base model.
Accurate time delay technology in simulated test for high precision laser range finder
NASA Astrophysics Data System (ADS)
Chen, Zhibin; Xiao, Wenjian; Wang, Weiming; Xue, Mingxi
2015-10-01
With the continuous development of technology, the ranging accuracy of pulsed laser range finder (LRF) is higher and higher, so the maintenance demand of LRF is also rising. According to the dominant ideology of "time analog spatial distance" in simulated test for pulsed range finder, the key of distance simulation precision lies in the adjustable time delay. By analyzing and comparing the advantages and disadvantages of fiber and circuit delay, a method was proposed to improve the accuracy of the circuit delay without increasing the count frequency of the circuit. A high precision controllable delay circuit was designed by combining the internal delay circuit and external delay circuit which could compensate the delay error in real time. And then the circuit delay accuracy could be increased. The accuracy of the novel circuit delay methods proposed in this paper was actually measured by a high sampling rate oscilloscope actual measurement. The measurement result shows that the accuracy of the distance simulated by the circuit delay is increased from +/- 0.75m up to +/- 0.15m. The accuracy of the simulated distance is greatly improved in simulated test for high precision pulsed range finder.
Enabling R&D for accurate simulation of non-ideal explosives.
Aidun, John Bahram; Thompson, Aidan Patrick; Schmitt, Robert Gerard
2010-09-01
We implemented two numerical simulation capabilities essential to reliably predicting the effect of non-ideal explosives (NXs). To begin to be able to treat the multiple, competing, multi-step reaction paths and slower kinetics of NXs, Sandia's CTH shock physics code was extended to include the TIGER thermochemical equilibrium solver as an in-line routine. To facilitate efficient exploration of reaction pathways that need to be identified for the CTH simulations, we implemented in Sandia's LAMMPS molecular dynamics code the MSST method, which is a reactive molecular dynamics technique for simulating steady shock wave response. Our preliminary demonstrations of these two capabilities serve several purposes: (i) they demonstrate proof-of-principle for our approach; (ii) they provide illustration of the applicability of the new functionality; and (iii) they begin to characterize the use of the new functionality and identify where improvements will be needed for the ultimate capability to meet national security needs. Next steps are discussed.
Offner, Stella S. R.; Robitaille, Thomas P.; Hansen, Charles E.; Klein, Richard I.; McKee, Christopher F.
2012-07-10
The properties of unresolved protostars and their local environment are frequently inferred from spectral energy distributions (SEDs) using radiative transfer modeling. In this paper, we use synthetic observations of realistic star formation simulations to evaluate the accuracy of properties inferred from fitting model SEDs to observations. We use ORION, an adaptive mesh refinement (AMR) three-dimensional gravito-radiation-hydrodynamics code, to simulate low-mass star formation in a turbulent molecular cloud including the effects of protostellar outflows. To obtain the dust temperature distribution and SEDs of the forming protostars, we post-process the simulations using HYPERION, a state-of-the-art Monte Carlo radiative transfer code. We find that the ORION and HYPERION dust temperatures typically agree within a factor of two. We compare synthetic SEDs of embedded protostars for a range of evolutionary times, simulation resolutions, aperture sizes, and viewing angles. We demonstrate that complex, asymmetric gas morphology leads to a variety of classifications for individual objects as a function of viewing angle. We derive best-fit source parameters for each SED through comparison with a pre-computed grid of radiative transfer models. While the SED models correctly identify the evolutionary stage of the synthetic sources as embedded protostars, we show that the disk and stellar parameters can be very discrepant from the simulated values, which is expected since the disk and central source are obscured by the protostellar envelope. Parameters such as the stellar accretion rate, stellar mass, and disk mass show better agreement, but can still deviate significantly, and the agreement may in some cases be artificially good due to the limited range of parameters in the set of model SEDs. Lack of correlation between the model and simulation properties in many individual instances cautions against overinterpreting properties inferred from SEDs for unresolved protostellar
A hybrid method for efficient and accurate simulations of diffusion compartment imaging signals
NASA Astrophysics Data System (ADS)
Rensonnet, Gaëtan; Jacobs, Damien; Macq, Benoît; Taquet, Maxime
2015-12-01
Diffusion-weighted imaging is sensitive to the movement of water molecules through the tissue microstructure and can therefore be used to gain insight into the tissue cellular architecture. While the diffusion signal arising from simple geometrical microstructure is known analytically, it remains unclear what diffusion signal arises from complex microstructural configurations. Such knowledge is important to design optimal acquisition sequences, to understand the limitations of diffusion-weighted imaging and to validate novel models of the brain microstructure. We present a novel framework for the efficient simulation of high-quality DW-MRI signals based on the hybrid combination of exact analytic expressions in simple geometric compartments such as cylinders and spheres and Monte Carlo simulations in more complex geometries. We validate our approach on synthetic arrangements of parallel cylinders representing the geometry of white matter fascicles, by comparing it to complete, all-out Monte Carlo simulations commonly used in the literature. For typical configurations, equal levels of accuracy are obtained with our hybrid method in less than one fifth of the computational time required for Monte Carlo simulations.
NASA Astrophysics Data System (ADS)
Yi, Sha-Sha; Pan, Cong; Hu, Zhong-Han
2015-12-01
Modern computer simulations of biological systems often involve an explicit treatment of the complex interactions among a large number of molecules. While it is straightforward to compute the short-ranged Van der Waals interaction in classical molecular dynamics simulations, it has been a long-lasting issue to develop accurate methods for the longranged Coulomb interaction. In this short review, we discuss three types of methodologies for the accurate treatment of electrostatics in simulations of explicit molecules: truncation-type methods, Ewald-type methods, and mean-field-type methods. Throughout the discussion, we brief the formulations and developments of these methods, emphasize the intrinsic connections among the three types of methods, and focus on the existing problems which are often associated with the boundary conditions of electrostatics. This brief survey is summarized with a short perspective on future trends along the method developments and applications in the field of biological simulations. Project supported by the National Natural Science Foundation of China (Grant Nos. 91127015 and 21522304) and the Open Project from the State Key Laboratory of Theoretical Physics, and the Innovation Project from the State Key Laboratory of Supramolecular Structure and Materials.
Zhao, A P; Cvetkovic, S R
1994-08-20
An efficient, accurate, and automated vectorial finite-element software package (named WAVEGIDE), which is implemented within a PDE/Protran problem-solving environment, has been extended to general multilayer anisotropic waveguides. With our system, through an interactive question-and-answer session, the problem can be simply defined with high-level PDE/Protran commands. The problem can then be solved easily and quickly by the main processor within this intelligent environment. In particular, in our system the eigenvalue of waveguide problems may be either a propagation constant (β) or an operated light frequency (F). Furthermore, the cutoff frequencies of propagation modes in waveguides can be calculated. As an application of this approach, numerical results for both scalar and hybrid modes in multilayer anisotropic waveguides are presented and are also compared with results obtained with the domain-integral method. These results clearly illustrate the unique flexibility, accuracy, and the ease of use f the WAVEGIDE program. PMID:20935964
Chemically Accurate Simulation of a Polyatomic Molecule-Metal Surface Reaction.
Nattino, Francesco; Migliorini, Davide; Kroes, Geert-Jan; Dombrowski, Eric; High, Eric A; Killelea, Daniel R; Utz, Arthur L
2016-07-01
Although important to heterogeneous catalysis, the ability to accurately model reactions of polyatomic molecules with metal surfaces has not kept pace with developments in gas phase dynamics. Partnering the specific reaction parameter (SRP) approach to density functional theory with ab initio molecular dynamics (AIMD) extends our ability to model reactions with metals with quantitative accuracy from only the lightest reactant, H2, to essentially all molecules. This is demonstrated with AIMD calculations on CHD3 + Ni(111) in which the SRP functional is fitted to supersonic beam experiments, and validated by showing that AIMD with the resulting functional reproduces initial-state selected sticking measurements with chemical accuracy (4.2 kJ/mol ≈ 1 kcal/mol). The need for only semilocal exchange makes our scheme computationally tractable for dissociation on transition metals. PMID:27284787
Chemically Accurate Simulation of a Polyatomic Molecule-Metal Surface Reaction
2016-01-01
Although important to heterogeneous catalysis, the ability to accurately model reactions of polyatomic molecules with metal surfaces has not kept pace with developments in gas phase dynamics. Partnering the specific reaction parameter (SRP) approach to density functional theory with ab initio molecular dynamics (AIMD) extends our ability to model reactions with metals with quantitative accuracy from only the lightest reactant, H2, to essentially all molecules. This is demonstrated with AIMD calculations on CHD3 + Ni(111) in which the SRP functional is fitted to supersonic beam experiments, and validated by showing that AIMD with the resulting functional reproduces initial-state selected sticking measurements with chemical accuracy (4.2 kJ/mol ≈ 1 kcal/mol). The need for only semilocal exchange makes our scheme computationally tractable for dissociation on transition metals. PMID:27284787
Differential-equation-based representation of truncation errors for accurate numerical simulation
NASA Astrophysics Data System (ADS)
MacKinnon, Robert J.; Johnson, Richard W.
1991-09-01
High-order compact finite difference schemes for 2D convection-diffusion-type differential equations with constant and variable convection coefficients are derived. The governing equations are employed to represent leading truncation terms, including cross-derivatives, making the overall O(h super 4) schemes conform to a 3 x 3 stencil. It is shown that the two-dimensional constant coefficient scheme collapses to the optimal scheme for the one-dimensional case wherein the finite difference equation yields nodally exact results. The two-dimensional schemes are tested against standard model problems, including a Navier-Stokes application. Results show that the two schemes are generally more accurate, on comparable grids, than O(h super 2) centered differencing and commonly used O(h) and O(h super 3) upwinding schemes.
A new class of accurate, mesh-free hydrodynamic simulation methods
NASA Astrophysics Data System (ADS)
Hopkins, Philip F.
2015-06-01
We present two new Lagrangian methods for hydrodynamics, in a systematic comparison with moving-mesh, smoothed particle hydrodynamics (SPH), and stationary (non-moving) grid methods. The new methods are designed to simultaneously capture advantages of both SPH and grid-based/adaptive mesh refinement (AMR) schemes. They are based on a kernel discretization of the volume coupled to a high-order matrix gradient estimator and a Riemann solver acting over the volume `overlap'. We implement and test a parallel, second-order version of the method with self-gravity and cosmological integration, in the code GIZMO:1 this maintains exact mass, energy and momentum conservation; exhibits superior angular momentum conservation compared to all other methods we study; does not require `artificial diffusion' terms; and allows the fluid elements to move with the flow, so resolution is automatically adaptive. We consider a large suite of test problems, and find that on all problems the new methods appear competitive with moving-mesh schemes, with some advantages (particularly in angular momentum conservation), at the cost of enhanced noise. The new methods have many advantages versus SPH: proper convergence, good capturing of fluid-mixing instabilities, dramatically reduced `particle noise' and numerical viscosity, more accurate sub-sonic flow evolution, and sharp shock-capturing. Advantages versus non-moving meshes include: automatic adaptivity, dramatically reduced advection errors and numerical overmixing, velocity-independent errors, accurate coupling to gravity, good angular momentum conservation and elimination of `grid alignment' effects. We can, for example, follow hundreds of orbits of gaseous discs, while AMR and SPH methods break down in a few orbits. However, fixed meshes minimize `grid noise'. These differences are important for a range of astrophysical problems.
Hepburn, I; Chen, W; De Schutter, E
2016-08-01
Spatial stochastic molecular simulations in biology are limited by the intense computation required to track molecules in space either in a discrete time or discrete space framework, which has led to the development of parallel methods that can take advantage of the power of modern supercomputers in recent years. We systematically test suggested components of stochastic reaction-diffusion operator splitting in the literature and discuss their effects on accuracy. We introduce an operator splitting implementation for irregular meshes that enhances accuracy with minimal performance cost. We test a range of models in small-scale MPI simulations from simple diffusion models to realistic biological models and find that multi-dimensional geometry partitioning is an important consideration for optimum performance. We demonstrate performance gains of 1-3 orders of magnitude in the parallel implementation, with peak performance strongly dependent on model specification. PMID:27497550
NASA Astrophysics Data System (ADS)
Hepburn, I.; Chen, W.; De Schutter, E.
2016-08-01
Spatial stochastic molecular simulations in biology are limited by the intense computation required to track molecules in space either in a discrete time or discrete space framework, which has led to the development of parallel methods that can take advantage of the power of modern supercomputers in recent years. We systematically test suggested components of stochastic reaction-diffusion operator splitting in the literature and discuss their effects on accuracy. We introduce an operator splitting implementation for irregular meshes that enhances accuracy with minimal performance cost. We test a range of models in small-scale MPI simulations from simple diffusion models to realistic biological models and find that multi-dimensional geometry partitioning is an important consideration for optimum performance. We demonstrate performance gains of 1-3 orders of magnitude in the parallel implementation, with peak performance strongly dependent on model specification.
H2 Adsorption in a Porous Crystal: Accurate First-Principles Quantum Simulation.
D'Arcy, Jordan H; Jordan, Meredith J T; Frankcombe, Terry J; Collins, Michael A
2015-12-17
A general method is presented for constructing, from ab initio quantum chemistry calculations, the potential energy surface (PES) for H2 absorbed in a porous crystalline material. The method is illustrated for the metal-organic framework material MOF-5. Rigid body quantum diffusion Monte Carlo simulations are used in the construction of the PES and to evaluate the quantum ground state of H2 in MOF-5, the zero-point energy, and the enthalpy of adsorption at 0 K. PMID:26322374
Schwörer, Magnus; Lorenzen, Konstantin; Mathias, Gerald; Tavan, Paul
2015-03-14
Recently, a novel approach to hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations has been suggested [Schwörer et al., J. Chem. Phys. 138, 244103 (2013)]. Here, the forces acting on the atoms are calculated by grid-based density functional theory (DFT) for a solute molecule and by a polarizable molecular mechanics (PMM) force field for a large solvent environment composed of several 10(3)-10(5) molecules as negative gradients of a DFT/PMM hybrid Hamiltonian. The electrostatic interactions are efficiently described by a hierarchical fast multipole method (FMM). Adopting recent progress of this FMM technique [Lorenzen et al., J. Chem. Theory Comput. 10, 3244 (2014)], which particularly entails a strictly linear scaling of the computational effort with the system size, and adapting this revised FMM approach to the computation of the interactions between the DFT and PMM fragments of a simulation system, here, we show how one can further enhance the efficiency and accuracy of such DFT/PMM-MD simulations. The resulting gain of total performance, as measured for alanine dipeptide (DFT) embedded in water (PMM) by the product of the gains in efficiency and accuracy, amounts to about one order of magnitude. We also demonstrate that the jointly parallelized implementation of the DFT and PMM-MD parts of the computation enables the efficient use of high-performance computing systems. The associated software is available online. PMID:25770527
Accurate simulation of the electron cloud in the Fermilab Main Injector with VORPAL
Lebrun, Paul L.G.; Spentzouris, Panagiotis; Cary, John R.; Stoltz, Peter; Veitzer, Seth A.; /Tech-X, Boulder
2011-01-01
We present results from a precision simulation of the electron cloud (EC) in the Fermilab Main Injector using the code VORPAL. This is a fully 3d and self consistent treatment of the EC. Both distributions of electrons in 6D phase-space and E.M. field maps have been generated. This has been done for various configurations of the magnetic fields found around the machine have been studied. Plasma waves associated to the fluctuation density of the cloud have been analyzed. Our results are compared with those obtained with the POSINST code. The response of a Retarding Field Analyzer (RFA) to the EC has been simulated, as well as the more challenging microwave absorption experiment. Definite predictions of their exact response are difficult to obtain,mostly because of the uncertainties in the secondary emission yield and, in the case of the RFA, because of the sensitivity of the electron collection efficiency to unknown stray magnetic fields. Nonetheless, our simulations do provide guidance to the experimental program.
Pre-Stall Behavior of a Transonic Axial Compressor Stage via Time-Accurate Numerical Simulation
NASA Technical Reports Server (NTRS)
Chen, Jen-Ping; Hathaway, Michael D.; Herrick, Gregory P.
2008-01-01
CFD calculations using high-performance parallel computing were conducted to simulate the pre-stall flow of a transonic compressor stage, NASA compressor Stage 35. The simulations were run with a full-annulus grid that models the 3D, viscous, unsteady blade row interaction without the need for an artificial inlet distortion to induce stall. The simulation demonstrates the development of the rotating stall from the growth of instabilities. Pressure-rise performance and pressure traces are compared with published experimental data before the study of flow evolution prior to the rotating stall. Spatial FFT analysis of the flow indicates a rotating long-length disturbance of one rotor circumference, which is followed by a spike-type breakdown. The analysis also links the long-length wave disturbance with the initiation of the spike inception. The spike instabilities occur when the trajectory of the tip clearance flow becomes perpendicular to the axial direction. When approaching stall, the passage shock changes from a single oblique shock to a dual-shock, which distorts the perpendicular trajectory of the tip clearance vortex but shows no evidence of flow separation that may contribute to stall.
Schwörer, Magnus; Lorenzen, Konstantin; Mathias, Gerald; Tavan, Paul
2015-03-14
Recently, a novel approach to hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations has been suggested [Schwörer et al., J. Chem. Phys. 138, 244103 (2013)]. Here, the forces acting on the atoms are calculated by grid-based density functional theory (DFT) for a solute molecule and by a polarizable molecular mechanics (PMM) force field for a large solvent environment composed of several 10{sup 3}-10{sup 5} molecules as negative gradients of a DFT/PMM hybrid Hamiltonian. The electrostatic interactions are efficiently described by a hierarchical fast multipole method (FMM). Adopting recent progress of this FMM technique [Lorenzen et al., J. Chem. Theory Comput. 10, 3244 (2014)], which particularly entails a strictly linear scaling of the computational effort with the system size, and adapting this revised FMM approach to the computation of the interactions between the DFT and PMM fragments of a simulation system, here, we show how one can further enhance the efficiency and accuracy of such DFT/PMM-MD simulations. The resulting gain of total performance, as measured for alanine dipeptide (DFT) embedded in water (PMM) by the product of the gains in efficiency and accuracy, amounts to about one order of magnitude. We also demonstrate that the jointly parallelized implementation of the DFT and PMM-MD parts of the computation enables the efficient use of high-performance computing systems. The associated software is available online.
NASA Astrophysics Data System (ADS)
Schwörer, Magnus; Lorenzen, Konstantin; Mathias, Gerald; Tavan, Paul
2015-03-01
Recently, a novel approach to hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations has been suggested [Schwörer et al., J. Chem. Phys. 138, 244103 (2013)]. Here, the forces acting on the atoms are calculated by grid-based density functional theory (DFT) for a solute molecule and by a polarizable molecular mechanics (PMM) force field for a large solvent environment composed of several 103-105 molecules as negative gradients of a DFT/PMM hybrid Hamiltonian. The electrostatic interactions are efficiently described by a hierarchical fast multipole method (FMM). Adopting recent progress of this FMM technique [Lorenzen et al., J. Chem. Theory Comput. 10, 3244 (2014)], which particularly entails a strictly linear scaling of the computational effort with the system size, and adapting this revised FMM approach to the computation of the interactions between the DFT and PMM fragments of a simulation system, here, we show how one can further enhance the efficiency and accuracy of such DFT/PMM-MD simulations. The resulting gain of total performance, as measured for alanine dipeptide (DFT) embedded in water (PMM) by the product of the gains in efficiency and accuracy, amounts to about one order of magnitude. We also demonstrate that the jointly parallelized implementation of the DFT and PMM-MD parts of the computation enables the efficient use of high-performance computing systems. The associated software is available online.
Subramanian, Swetha; Mast, T Douglas
2015-10-01
Computational finite element models are commonly used for the simulation of radiofrequency ablation (RFA) treatments. However, the accuracy of these simulations is limited by the lack of precise knowledge of tissue parameters. In this technical note, an inverse solver based on the unscented Kalman filter (UKF) is proposed to optimize values for specific heat, thermal conductivity, and electrical conductivity resulting in accurately simulated temperature elevations. A total of 15 RFA treatments were performed on ex vivo bovine liver tissue. For each RFA treatment, 15 finite-element simulations were performed using a set of deterministically chosen tissue parameters to estimate the mean and variance of the resulting tissue ablation. The UKF was implemented as an inverse solver to recover the specific heat, thermal conductivity, and electrical conductivity corresponding to the measured area of the ablated tissue region, as determined from gross tissue histology. These tissue parameters were then employed in the finite element model to simulate the position- and time-dependent tissue temperature. Results show good agreement between simulated and measured temperature. PMID:26352462
NASA Astrophysics Data System (ADS)
Subramanian, Swetha; Mast, T. Douglas
2015-09-01
Computational finite element models are commonly used for the simulation of radiofrequency ablation (RFA) treatments. However, the accuracy of these simulations is limited by the lack of precise knowledge of tissue parameters. In this technical note, an inverse solver based on the unscented Kalman filter (UKF) is proposed to optimize values for specific heat, thermal conductivity, and electrical conductivity resulting in accurately simulated temperature elevations. A total of 15 RFA treatments were performed on ex vivo bovine liver tissue. For each RFA treatment, 15 finite-element simulations were performed using a set of deterministically chosen tissue parameters to estimate the mean and variance of the resulting tissue ablation. The UKF was implemented as an inverse solver to recover the specific heat, thermal conductivity, and electrical conductivity corresponding to the measured area of the ablated tissue region, as determined from gross tissue histology. These tissue parameters were then employed in the finite element model to simulate the position- and time-dependent tissue temperature. Results show good agreement between simulated and measured temperature.
Accurate simulation of MPPT methods performance when applied to commercial photovoltaic panels.
Cubas, Javier; Pindado, Santiago; Sanz-Andrés, Ángel
2015-01-01
A new, simple, and quick-calculation methodology to obtain a solar panel model, based on the manufacturers' datasheet, to perform MPPT simulations, is described. The method takes into account variations on the ambient conditions (sun irradiation and solar cells temperature) and allows fast MPPT methods comparison or their performance prediction when applied to a particular solar panel. The feasibility of the described methodology is checked with four different MPPT methods applied to a commercial solar panel, within a day, and under realistic ambient conditions. PMID:25874262
Accurate Simulation of MPPT Methods Performance When Applied to Commercial Photovoltaic Panels
2015-01-01
A new, simple, and quick-calculation methodology to obtain a solar panel model, based on the manufacturers' datasheet, to perform MPPT simulations, is described. The method takes into account variations on the ambient conditions (sun irradiation and solar cells temperature) and allows fast MPPT methods comparison or their performance prediction when applied to a particular solar panel. The feasibility of the described methodology is checked with four different MPPT methods applied to a commercial solar panel, within a day, and under realistic ambient conditions. PMID:25874262
High-order accurate multi-phase simulations: building blocks and whats tricky about them
NASA Astrophysics Data System (ADS)
Kummer, Florian
2015-11-01
We are going to present a high-order numerical method for multi-phase flow problems, which employs a sharp interface representation by a level-set and an extended discontinuous Galerkin (XDG) discretization for the flow properties. The shape of the XDG basis functions is dynamically adapted to the position of the fluid interface, so that the spatial approximation space can represent jumps in pressure and kinks in velocity accurately. By this approach, the `hp-convergence' property of the classical discontinuous Galerkin (DG) method can be preserved for the low-regularity, discontinuous solutions, such as those appearing in multi-phase flows. Within the past years, several building blocks of such a method were presented: this includes numerical integration on cut-cells, the spatial discretization by the XDG method, precise evaluation of curvature and level-set algorithms tailored to the special requirements of XDG-methods. The presentation covers a short review on these building-block and their integration into a full multi-phase solver. A special emphasis is put on the discussion of the several pitfalls one may expire in the formulation of such a solver. German Research Foundation.
Zhou, Nengji; Chen, Lipeng; Huang, Zhongkai; Sun, Kewei; Tanimura, Yoshitaka; Zhao, Yang
2016-03-10
By employing the Dirac-Frenkel time-dependent variational principle, we study the dynamical properties of the Holstein molecular crystal model with diagonal and off-diagonal exciton-phonon coupling. A linear combination of the Davydov D1 (D2) ansatz, referred to as the "multi-D1 ansatz" ("multi-D2 ansatz"), is used as the trial state with enhanced accuracy but without sacrificing efficiency. The time evolution of the exciton probability is found to be in perfect agreement with that of the hierarchy equations of motion, demonstrating the promise the multiple Davydov trial states hold as an efficient, robust description of dynamics of complex quantum systems. In addition to the linear absorption spectra computed for both diagonal and off-diagonal cases, for the first time, 2D spectra have been calculated for systems with off-diagonal exciton-phonon coupling by employing the multiple D2 ansatz to compute the nonlinear response function, testifying to the great potential of the multiple D2 ansatz for fast, accurate implementation of multidimensional spectroscopy. It is found that the signal exhibits a single peak for weak off-diagonal coupling, while a vibronic multipeak structure appears for strong off-diagonal coupling. PMID:26871592
Consistent Multigroup Theory Enabling Accurate Course-Group Simulation of Gen IV Reactors
Rahnema, Farzad; Haghighat, Alireza; Ougouag, Abderrafi
2013-11-29
The objective of this proposal is the development of a consistent multi-group theory that accurately accounts for the energy-angle coupling associated with collapsed-group cross sections. This will allow for coarse-group transport and diffusion theory calculations that exhibit continuous energy accuracy and implicitly treat cross- section resonances. This is of particular importance when considering the highly heterogeneous and optically thin reactor designs within the Next Generation Nuclear Plant (NGNP) framework. In such reactors, ignoring the influence of anisotropy in the angular flux on the collapsed cross section, especially at the interface between core and reflector near which control rods are located, results in inaccurate estimates of the rod worth, a serious safety concern. The scope of this project will include the development and verification of a new multi-group theory enabling high-fidelity transport and diffusion calculations in coarse groups, as well as a methodology for the implementation of this method in existing codes. This will allow for a higher accuracy solution of reactor problems while using fewer groups and will reduce the computational expense. The proposed research represents a fundamental advancement in the understanding and improvement of multi- group theory for reactor analysis.
How accurate are volcanic ash simulations of the 2010 Eyjafjallajökull eruption?
NASA Astrophysics Data System (ADS)
Dacre, Helen; Harvey, Natalie; Webley, Peter; Morton, Don
2016-04-01
In the event of a volcanic eruption the decision to close airspace is based on forecast ash maps, produced using volcanic ash transport and dispersion models. In this paper we quantitatively evaluate the spatial skill of volcanic ash simulations using satellite retrievals of ash from the Eyjafjallajökull eruption during the period from 7-16 May 2010. We find that at the start of this period, 7-10 May, the model (FLEXPART) has excellent skill and can predict the spatial distribution of the satellite retrieved ash to within 0.5°× 0.5° lat/lon. However, on the 10 May there is a decrease in the spatial accuracy of the model, to 2.5°× 2.5° lat/lon, and between 11-12 May the simulated ash location errors grow rapidly. On the 11 May ash is located close to a bifurcation point in the atmosphere, resulting in a rapid divergence in the modeled and satellite ash locations. In general, the model skill reduces as the residence time of ash increases. However, the error growth is not always steady. Rapid increases in error growth are linked to critical points in the ash trajectories. Ensemble modeling using perturbed meteorological data would help to represent this uncertainty and assimilation of satellite ash data would help to reduce uncertainty in volcanic ash forecasts.
How accurate are volcanic ash simulations of the 2010 Eyjafjallajökull eruption?
NASA Astrophysics Data System (ADS)
Dacre, H. F.; Harvey, N. J.; Webley, P. W.; Morton, D.
2016-04-01
In the event of a volcanic eruption the decision to close airspace is based on forecast ash maps, produced using volcanic ash transport and dispersion models. In this paper we quantitatively evaluate the spatial skill of volcanic ash simulations using satellite retrievals of ash from the Eyjafjallajökull eruption during the period from 7 to 16 May 2010. We find that at the start of this period, 7-10 May, the model (FLEXible PARTicle) has excellent skill and can predict the spatial distribution of the satellite-retrieved ash to within 0.5° × 0.5° latitude/longitude. However, on 10 May there is a decrease in the spatial accuracy of the model to 2.5°× 2.5° latitude/longitude, and between 11 and 12 May the simulated ash location errors grow rapidly. On 11 May ash is located close to a bifurcation point in the atmosphere, resulting in a rapid divergence in the modeled and satellite ash locations. In general, the model skill reduces as the residence time of ash increases. However, the error growth is not always steady. Rapid increases in error growth are linked to key points in the ash trajectories. Ensemble modeling using perturbed meteorological data would help to represent this uncertainty, and assimilation of satellite ash data would help to reduce uncertainty in volcanic ash forecasts.
Lippert, Ross A; Predescu, Cristian; Ierardi, Douglas J; Mackenzie, Kenneth M; Eastwood, Michael P; Dror, Ron O; Shaw, David E
2013-10-28
In molecular dynamics simulations, control over temperature and pressure is typically achieved by augmenting the original system with additional dynamical variables to create a thermostat and a barostat, respectively. These variables generally evolve on timescales much longer than those of particle motion, but typical integrator implementations update the additional variables along with the particle positions and momenta at each time step. We present a framework that replaces the traditional integration procedure with separate barostat, thermostat, and Newtonian particle motion updates, allowing thermostat and barostat updates to be applied infrequently. Such infrequent updates provide a particularly substantial performance advantage for simulations parallelized across many computer processors, because thermostat and barostat updates typically require communication among all processors. Infrequent updates can also improve accuracy by alleviating certain sources of error associated with limited-precision arithmetic. In addition, separating the barostat, thermostat, and particle motion update steps reduces certain truncation errors, bringing the time-average pressure closer to its target value. Finally, this framework, which we have implemented on both general-purpose and special-purpose hardware, reduces software complexity and improves software modularity. PMID:24182003
A novel fast and accurate pseudo-analytical simulation approach for MOAO
NASA Astrophysics Data System (ADS)
Gendron, É.; Charara, A.; Abdelfattah, A.; Gratadour, D.; Keyes, D.; Ltaief, H.; Morel, C.; Vidal, F.; Sevin, A.; Rousset, G.
2014-08-01
Multi-object adaptive optics (MOAO) is a novel adaptive optics (AO) technique for wide-field multi-object spectrographs (MOS). MOAO aims at applying dedicated wavefront corrections to numerous separated tiny patches spread over a large field of view (FOV), limited only by that of the telescope. The control of each deformable mirror (DM) is done individually using a tomographic reconstruction of the phase based on measurements from a number of wavefront sensors (WFS) pointing at natural and artificial guide stars in the field. We have developed a novel hybrid, pseudo-analytical simulation scheme, somewhere in between the end-to- end and purely analytical approaches, that allows us to simulate in detail the tomographic problem as well as noise and aliasing with a high fidelity, and including fitting and bandwidth errors thanks to a Fourier-based code. Our tomographic approach is based on the computation of the minimum mean square error (MMSE) reconstructor, from which we derive numerically the covariance matrix of the tomographic error, including aliasing and propagated noise. We are then able to simulate the point-spread function (PSF) associated to this covariance matrix of the residuals, like in PSF reconstruction algorithms. The advantage of our approach is that we compute the same tomographic reconstructor that would be computed when operating the real instrument, so that our developments open the way for a future on-sky implementation of the tomographic control, plus the joint PSF and performance estimation. The main challenge resides in the computation of the tomographic reconstructor which involves the inversion of a large matrix (typically 40 000 × 40 000 elements). To perform this computation efficiently, we chose an optimized approach based on the use of GPUs as accelerators and using an optimized linear algebra library: MORSE providing a significant speedup against standard CPU oriented libraries such as Intel MKL. Because the covariance matrix is
Cohen, Trevor; Blatter, Brett; Patel, Vimla
2008-01-01
Cognitive studies reveal that less-than-expert clinicians are less able to recognize meaningful patterns of data in clinical narratives. Accordingly, psychiatric residents early in training fail to attend to information that is relevant to diagnosis and the assessment of dangerousness. This manuscript presents cognitively motivated methodology for the simulation of expert ability to organize relevant findings supporting intermediate diagnostic hypotheses. Latent Semantic Analysis is used to generate a semantic space from which meaningful associations between psychiatric terms are derived. Diagnostically meaningful clusters are modeled as geometric structures within this space and compared to elements of psychiatric narrative text using semantic distance measures. A learning algorithm is defined that alters components of these geometric structures in response to labeled training data. Extraction and classification of relevant text segments is evaluated against expert annotation, with system-rater agreement approximating rater-rater agreement. A range of biomedical informatics applications for these methods are suggested. PMID:18455483
TRIM—3D: a three-dimensional model for accurate simulation of shallow water flow
Casulli, Vincenzo; Bertolazzi, Enrico; Cheng, Ralph T.
1993-01-01
A semi-implicit finite difference formulation for the numerical solution of three-dimensional tidal circulation is discussed. The governing equations are the three-dimensional Reynolds equations in which the pressure is assumed to be hydrostatic. A minimal degree of implicitness has been introduced in the finite difference formula so that the resulting algorithm permits the use of large time steps at a minimal computational cost. This formulation includes the simulation of flooding and drying of tidal flats, and is fully vectorizable for an efficient implementation on modern vector computers. The high computational efficiency of this method has made it possible to provide the fine details of circulation structure in complex regions that previous studies were unable to obtain. For proper interpretation of the model results suitable interactive graphics is also an essential tool.
Evaluation of the EURO-CORDEX RCMs to accurately simulate the Etesian wind system
NASA Astrophysics Data System (ADS)
Dafka, Stella; Xoplaki, Elena; Toreti, Andrea; Zanis, Prodromos; Tyrlis, Evangelos; Luterbacher, Jürg
2016-04-01
The Etesians are among the most persistent regional scale wind systems in the lower troposphere that blow over the Aegean Sea during the extended summer season. ΑAn evaluation of the high spatial resolution, EURO-CORDEX Regional Climate Models (RCMs) is here presented. The study documents the performance of the individual models in representing the basic spatiotemporal pattern of the Etesian wind system for the period 1989-2004. The analysis is mainly focused on evaluating the abilities of the RCMs in simulating the surface wind over the Aegean Sea and the associated large scale atmospheric circulation. Mean Sea Level Pressure (SLP), wind speed and geopotential height at 500 hPa are used. The simulated results are validated against reanalysis datasets (20CR-v2c and ERA20-C) and daily observational measurements (12:00 UTC) from the mainland Greece and Aegean Sea. The analysis highlights the general ability of the RCMs to capture the basic features of the Etesians, but also indicates considerable deficiencies for selected metrics, regions and subperiods. Some of these deficiencies include the significant underestimation (overestimation) of the mean SLP in the northeastern part of the analysis domain in all subperiods (for May and June) when compared to 20CR-v2c (ERA20-C), the significant overestimation of the anomalous ridge over the Balkans and central Europe and the underestimation of the wind speed over the Aegean Sea. Future work will include an assessment of the Etesians for the next decades using EURO-CORDEX projections under different RCP scenarios and estimate the future potential for wind energy production.
NASA Astrophysics Data System (ADS)
Ren, Zijian; Ma, Chunyong; Chen, Lu; Chen, Ge
2016-05-01
A high-efficiency anisotropic model for bidirectional reflectance distribution function (BRDF) of seawater covered by oil slicks (SWCOS) was proposed. This model was set by combining a BRDF model for anisotropic rough sea surface whose slopes follow Gaussian distribution and the two-beam inference theory of a thin film. We have simulated the BRDFs of oil slicks by using the above model and the measured complex refractive index data of Romashkino crude oil. In addition, the relationships between the BRDF of oil slicks and the wind speed of sea surface, thickness of oil slick, complex refractive index of crude oil and the incident zenith angle were analyzed. Also, the differences between optical characteristics of clean water and of polluted water were discussed in the context of the optical contrast of SWCOS. With high simulation speed and reliable simulation precision, this model provides a theoretical basis for rapid detection of oil spill.
A novel approach for accurate radiative transfer in cosmological hydrodynamic simulations
NASA Astrophysics Data System (ADS)
Petkova, Margarita; Springel, Volker
2011-08-01
accurately deal with non-equilibrium effects. We discuss several tests of the new method, including shadowing configurations in two and three dimensions, ionized sphere expansion in static and dynamic density fields and the ionization of a cosmological density field. The tests agree favourably with analytical expectations and results based on other numerical radiative transfer approximations.
Accurate simulation of near-wall turbulence over a compliant tensegrity fabric
NASA Astrophysics Data System (ADS)
Luo, Haoxiang; Bewley, Thomas R.
2005-05-01
This paper presents a new class of compliant surfaces, dubbed tensegrity fabrics, for the problem of reducing the drag induced by near-wall turbulent flows. The substructure upon which this compliant surface is built is based on the "tensegrity" structural paradigm, and is formed as a stable pretensioned network of compressive members ("bars") interconnected by tensile members ("tendons"). Compared with existing compliant surface studies, most of which are based on spring-supported plates or membranes, tensegrity fabrics appear to be better configured to respond to the shear stress fluctuations (in addition to the pressure fluctuations) generated by near-wall turbulence. As a result, once the several parameters affecting the compliance characteristics of the structure are tuned appropriately, the tensegrity fabric might exhibit an improved capacity for dampening the fluctuations of near-wall turbulence, thereby reducing drag. This paper improves our previous work (SPIE Paper 5049-57) and uses a 3D time-dependent coordinate transformation in the flow simulations to account for the motion of the channel walls, and the Cartesian components of the velocity are used as the flow variables. For the spatial discretization, a dealiased pseudospectral scheme is used in the homogeneous directions and a second-order finite difference scheme is used in the wall-normal direction. The code is first validated with several benchmark results that are available in the published literature for flows past both stationary and nonstationary walls. Direct numerical simulations of turbulent flows at Re_tau=150 over the compliant tensegrity fabric are then presented. It is found that, when the stiffness, mass, damping, and orientation of the members of the the unit cell defining the tensegrity fabric are selected appropriately, the near-wall statistics of the turbulence are altered significantly. The flow/structure interface is found to form streamwise-travelling waves reminiscent of those
Midgley, David J.; Greenfield, Paul; Shaw, Janet M.; Oytam, Yalchin; Li, Dongmei; Kerr, Caroline A.; Hendry, Philip
2012-01-01
The second generation (G2) PhyloChip is designed to detect over 8700 bacteria and archaeal and has been used over 50 publications and conference presentations. Many of those publications reveal that the PhyloChip measures of species richness greatly exceed statistical estimates of richness based on other methods. An examination of probes downloaded from Greengenes suggested that the system may have the potential to distort the observed community structure. This may be due to the sharing of probes by taxa; more than 21% of the taxa in that downloaded data have no unique probes. In-silico simulations using these data showed that a population of 64 taxa representing a typical anaerobic subterranean community returned 96 different taxa, including 15 families incorrectly called present and 19 families incorrectly called absent. A study of nasal and oropharyngeal microbial communities by Lemon et al (2010) found some 1325 taxa using the G2 PhyloChip, however, about 950 of these taxa have, in the downloaded data, no unique probes and cannot be definitively called present. Finally, data from Brodie et al (2007), when re-examined, indicate that the abundance of the majority of detected taxa, are highly correlated with one another, suggesting that many probe sets do not act independently. Based on our analyses of downloaded data, we conclude that outputs from the G2 PhyloChip should be treated with some caution, and that the presence of taxa represented solely by non-unique probes be independently verified. PMID:22457798
Gray, Alan; Harlen, Oliver G.; Harris, Sarah A.; Khalid, Syma; Leung, Yuk Ming; Lonsdale, Richard; Mulholland, Adrian J.; Pearson, Arwen R.; Read, Daniel J.; Richardson, Robin A.
2015-01-01
Despite huge advances in the computational techniques available for simulating biomolecules at the quantum-mechanical, atomistic and coarse-grained levels, there is still a widespread perception amongst the experimental community that these calculations are highly specialist and are not generally applicable by researchers outside the theoretical community. In this article, the successes and limitations of biomolecular simulation and the further developments that are likely in the near future are discussed. A brief overview is also provided of the experimental biophysical methods that are commonly used to probe biomolecular structure and dynamics, and the accuracy of the information that can be obtained from each is compared with that from modelling. It is concluded that progress towards an accurate spatial and temporal model of biomacromolecules requires a combination of all of these biophysical techniques, both experimental and computational. PMID:25615870
Dissipative particle dynamics simulation on the rheological properties of heavy crude oil
NASA Astrophysics Data System (ADS)
Wang, Sibo; Xu, Junbo; Wen, Hao
2015-11-01
The rheological properties of heavy crude oil have a significant impact on the production, refining and transportation. In this paper, dissipative particle dynamics (DPD) simulations were performed to study the effects of the addition of light crude oil and emulsification on the rheological properties of heavy crude oil. The simulation results reflected that the addition of light crude oil reduced the viscosity effectively. The shear thinning behaviour of crude oil mixtures were becoming less distinct as the increase of the mass fraction of light crude oil. According to the statistics, the shear had an influence on the aggregation and spatial orientation of asphaltene molecules. In addition, the relationship between the viscosity and the oil mass fraction was investigated in the simulations of emulsion systems. The viscosity increased with the oil mass fraction slowly in oil-in-water emulsions. When the oil mass fraction was higher than 50%, the increase became much faster since systems had been converted into water-in-oil emulsions. The equilibrated morphologies of emulsion systems were shown to illustrate the phase inversion. The surfactant-like feature of asphaltenes was also studied in the simulations.
NASA Astrophysics Data System (ADS)
Klostermann, U. K.; Mülders, T.; Schmöller, T.; Lorusso, G. F.; Hendrickx, E.
2010-04-01
In this paper, we discuss the performance of EUV resist models in terms of predictive accuracy, and we assess the readiness of the corresponding model calibration methodology. The study is done on an extensive OPC data set collected at IMEC for the ShinEtsu resist SEVR-59 on the ASML EUV Alpha Demo Tool (ADT), with the data set including more than thousand CD values. We address practical aspects such as the speed of calibration and selection of calibration patterns. The model is calibrated on 12 process window data series varying in pattern width (32, 36, 40 nm), orientation (H, V) and pitch (dense, isolated). The minimum measured feature size at nominal process condition is a 32 nm CD at a dense pitch of 64 nm. Mask metrology is applied to verify and eventually correct nominal width of the drawn CD. Cross-sectional SEM information is included in the calibration to tune the simulated resist loss and sidewall angle. The achieved calibration RMS is ~ 1.0 nm. We show what elements are important to obtain a well calibrated model. We discuss the impact of 3D mask effects on the Bossung tilt. We demonstrate that a correct representation of the flare level during the calibration is important to achieve a high predictability at various flare conditions. Although the model calibration is performed on a limited subset of the measurement data (one dimensional structures only), its accuracy is validated based on a large number of OPC patterns (at nominal dose and focus conditions) not included in the calibration; validation RMS results as small as 1 nm can be reached. Furthermore, we study the model's extendibility to two-dimensional end of line (EOL) structures. Finally, we correlate the experimentally observed fingerprint of the CD uniformity to a model, where EUV tool specific signatures are taken into account.
NASA Astrophysics Data System (ADS)
Zhang, Xiang; Vu-Quoc, Loc
2007-07-01
We present in this paper the displacement-driven version of a tangential force-displacement (TFD) model that accounts for both elastic and plastic deformations together with interfacial friction occurring in collisions of spherical particles. This elasto-plastic frictional TFD model, with its force-driven version presented in [L. Vu-Quoc, L. Lesburg, X. Zhang. An accurate tangential force-displacement model for granular-flow simulations: contacting spheres with plastic deformation, force-driven formulation, Journal of Computational Physics 196(1) (2004) 298-326], is consistent with the elasto-plastic frictional normal force-displacement (NFD) model presented in [L. Vu-Quoc, X. Zhang. An elasto-plastic contact force-displacement model in the normal direction: displacement-driven version, Proceedings of the Royal Society of London, Series A 455 (1991) 4013-4044]. Both the NFD model and the present TFD model are based on the concept of additive decomposition of the radius of contact area into an elastic part and a plastic part. The effect of permanent indentation after impact is represented by a correction to the radius of curvature. The effect of material softening due to plastic flow is represented by a correction to the elastic moduli. The proposed TFD model is accurate, and is validated against nonlinear finite element analyses involving plastic flows in both the loading and unloading conditions. The proposed consistent displacement-driven, elasto-plastic NFD and TFD models are designed for implementation in computer codes using the discrete-element method (DEM) for granular-flow simulations. The model is shown to be accurate and is validated against nonlinear elasto-plastic finite-element analysis.
NASA Astrophysics Data System (ADS)
Moreira, António H. J.; Queirós, Sandro; Morais, Pedro; Rodrigues, Nuno F.; Correia, André Ricardo; Fernandes, Valter; Pinho, A. C. M.; Fonseca, Jaime C.; Vilaça, João. L.
2015-03-01
The success of dental implant-supported prosthesis is directly linked to the accuracy obtained during implant's pose estimation (position and orientation). Although traditional impression techniques and recent digital acquisition methods are acceptably accurate, a simultaneously fast, accurate and operator-independent methodology is still lacking. Hereto, an image-based framework is proposed to estimate the patient-specific implant's pose using cone-beam computed tomography (CBCT) and prior knowledge of implanted model. The pose estimation is accomplished in a threestep approach: (1) a region-of-interest is extracted from the CBCT data using 2 operator-defined points at the implant's main axis; (2) a simulated CBCT volume of the known implanted model is generated through Feldkamp-Davis-Kress reconstruction and coarsely aligned to the defined axis; and (3) a voxel-based rigid registration is performed to optimally align both patient and simulated CBCT data, extracting the implant's pose from the optimal transformation. Three experiments were performed to evaluate the framework: (1) an in silico study using 48 implants distributed through 12 tridimensional synthetic mandibular models; (2) an in vitro study using an artificial mandible with 2 dental implants acquired with an i-CAT system; and (3) two clinical case studies. The results shown positional errors of 67+/-34μm and 108μm, and angular misfits of 0.15+/-0.08° and 1.4°, for experiment 1 and 2, respectively. Moreover, in experiment 3, visual assessment of clinical data results shown a coherent alignment of the reference implant. Overall, a novel image-based framework for implants' pose estimation from CBCT data was proposed, showing accurate results in agreement with dental prosthesis modelling requirements.
A specific PFT and sub-canopy structure for simulating oil palm in the Community Land Model
NASA Astrophysics Data System (ADS)
Fan, Y.; Knohl, A.; Roupsard, O.; Bernoux, M.; LE Maire, G.; Panferov, O.; Kotowska, M.; Meijide, A.
2015-12-01
Towards an effort to quantify the effects of rainforests to oil palm conversion on land-atmosphere carbon, water and energy fluxes, a specific plant functional type (PFT) and sub-canopy structure are developed for simulating oil palm within the Community Land Model (CLM4.5). Current global land surface models only simulate annual crops beside natural vegetation. In this study, a multilayer oil palm subroutine is developed in CLM4.5 for simulating oil palm's phenology and carbon and nitrogen allocation. The oil palm has monopodial morphology and sequential phenology of around 40 stacked phytomers, each carrying a large leaf and a fruit bunch, forming a natural multilayer canopy. A sub-canopy phenological and physiological parameterization is thus introduced, so that multiple phytomer components develop simultaneously but according to their different phenological steps (growth, yield and senescence) at different canopy layers. This specific multilayer structure was proved useful for simulating canopy development in terms of leaf area index (LAI) and fruit yield in terms of carbon and nitrogen outputs in Jambi, Sumatra (Fan et al. 2015). The study supports that species-specific traits, such as palm's monopodial morphology and sequential phenology, are necessary representations in terrestrial biosphere models in order to accurately simulate vegetation dynamics and feedbacks to climate. Further, oil palm's multilayer structure allows adding all canopy-level calculations of radiation, photosynthesis, stomatal conductance and respiration, beside phenology, also to the sub-canopy level, so as to eliminate scale mismatch problem among different processes. A series of adaptations are made to the CLM model. Initial results show that the adapted multilayer radiative transfer scheme and the explicit represention of oil palm's canopy structure improve on simulating photosynthesis-light response curve. The explicit photosynthesis and dynamic leaf nitrogen calculations per canopy
NASA Astrophysics Data System (ADS)
Farah, A.
The Ionospheric delay is still one of the largest sources of error that affects the positioning accuracy of any satellite positioning system. This problem could be solved due to the dispersive nature of the Ionosphere by combining simultaneous measurements of signals at two different frequencies but it is still there for single- frequency users. Much effort has been made in establishing models for single- frequency users to make this effect as small as possible. These models vary in accuracy, input data and computational complexity, so the choice between the different models depends on the individual circumstances of the user. From the simulation point of view, the model needed should be accurate with a global coverage and good description to the Ionosphere's variable nature with both time and location. The author reviews some of these established models, starting with the BENT model, the Klobuchar model and the IRI (International Reference Ionosphere) model. Since quiet a long time, Klobuchar model considers the most widely used model ever in this field, due to its simplicity and time saving. Any GPS user could find Klobuchar model's coefficients in the broadcast navigation message. CODE, Centre for Orbit Determination in Europe provides a new set of coefficients for Klobuchar model, which gives more accurate results for the Ionospheric delay computation. IGS (International GPS Service) services include providing GPS community with a global Ionospheric maps in IONEX-format (IONosphere Map Exchange format) which enables the computation of the Ionospheric delay at the desired location and time. The study was undertaken from GPS-data simulation point of view. The aim was to select a model for the simulation of GPS data that gives a good description of the Ionosphere's nature with a high degree of accuracy in computing the Ionospheric delay that yields to better-simulated data. A new model developed by the author based on IGS global Ionospheric maps. A comparison
DISPERSANT EFFECTIVENESS ON THREE OILS UNDER VARIOUS SIMULATED ENVIRONMENTAL CONDITIONS
The complexity of chemical and physical interactions between spilled oils, dispersants and the sea, necessitates an empirical approach for describing the interaction between the dispersant and oil slick which may provide a guide to dispersant effects on oil slicks. Recently, US ...
NASA Astrophysics Data System (ADS)
Shauly, Eitan; Rotstein, Israel; Peltinov, Ram; Latinski, Sergei; Adan, Ofer; Levi, Shimon; Menadeva, Ovadya
2009-03-01
The continues transistors scaling efforts, for smaller devices, similar (or larger) drive current/um and faster devices, increase the challenge to predict and to control the transistor off-state current. Typically, electrical simulators like SPICE, are using the design intent (as-drawn GDS data). At more sophisticated cases, the simulators are fed with the pattern after lithography and etch process simulations. As the importance of electrical simulation accuracy is increasing and leakage is becoming more dominant, there is a need to feed these simulators, with more accurate information extracted from physical on-silicon transistors. Our methodology to predict changes in device performances due to systematic lithography and etch effects was used in this paper. In general, the methodology consists on using the OPCCmaxTM for systematic Edge-Contour-Extraction (ECE) from transistors, taking along the manufacturing and includes any image distortions like line-end shortening, corner rounding and line-edge roughness. These measurements are used for SPICE modeling. Possible application of this new metrology is to provide a-head of time, physical and electrical statistical data improving time to market. In this work, we applied our methodology to analyze a small and large array's of 2.14um2 6T-SRAM, manufactured using Tower Standard Logic for General Purposes Platform. 4 out of the 6 transistors used "U-Shape AA", known to have higher variability. The predicted electrical performances of the transistors drive current and leakage current, in terms of nominal values and variability are presented. We also used the methodology to analyze an entire SRAM Block array. Study of an isolation leakage and variability are presented.
Gray, Alan; Harlen, Oliver G.; Harris, Sarah A.; Khalid, Syma; Leung, Yuk Ming; Lonsdale, Richard; Mulholland, Adrian J.; Pearson, Arwen R.; Read, Daniel J.; Richardson, Robin A.
2015-01-01
The current computational techniques available for biomolecular simulation are described, and the successes and limitations of each with reference to the experimental biophysical methods that they complement are presented. Despite huge advances in the computational techniques available for simulating biomolecules at the quantum-mechanical, atomistic and coarse-grained levels, there is still a widespread perception amongst the experimental community that these calculations are highly specialist and are not generally applicable by researchers outside the theoretical community. In this article, the successes and limitations of biomolecular simulation and the further developments that are likely in the near future are discussed. A brief overview is also provided of the experimental biophysical methods that are commonly used to probe biomolecular structure and dynamics, and the accuracy of the information that can be obtained from each is compared with that from modelling. It is concluded that progress towards an accurate spatial and temporal model of biomacromolecules requires a combination of all of these biophysical techniques, both experimental and computational.
NASA Astrophysics Data System (ADS)
Jolivet, L.; Cohen, M.; Ruas, A.
2015-08-01
Landscape influences fauna movement at different levels, from habitat selection to choices of movements' direction. Our goal is to provide a development frame in order to test simulation functions for animal's movement. We describe our approach for such simulations and we compare two types of functions to calculate trajectories. To do so, we first modelled the role of landscape elements to differentiate between elements that facilitate movements and the ones being hindrances. Different influences are identified depending on landscape elements and on animal species. Knowledge were gathered from ecologists, literature and observation datasets. Second, we analysed the description of animal movement recorded with GPS at fine scale, corresponding to high temporal frequency and good location accuracy. Analysing this type of data provides information on the relation between landscape features and movements. We implemented an agent-based simulation approach to calculate potential trajectories constrained by the spatial environment and individual's behaviour. We tested two functions that consider space differently: one function takes into account the geometry and the types of landscape elements and one cost function sums up the spatial surroundings of an individual. Results highlight the fact that the cost function exaggerates the distances travelled by an individual and simplifies movement patterns. The geometry accurate function represents a good bottom-up approach for discovering interesting areas or obstacles for movements.
Multigrid Methods for Fully Implicit Oil Reservoir Simulation
NASA Technical Reports Server (NTRS)
Molenaar, J.
1996-01-01
In this paper we consider the simultaneous flow of oil and water in reservoir rock. This displacement process is modeled by two basic equations: the material balance or continuity equations and the equation of motion (Darcy's law). For the numerical solution of this system of nonlinear partial differential equations there are two approaches: the fully implicit or simultaneous solution method and the sequential solution method. In the sequential solution method the system of partial differential equations is manipulated to give an elliptic pressure equation and a hyperbolic (or parabolic) saturation equation. In the IMPES approach the pressure equation is first solved, using values for the saturation from the previous time level. Next the saturations are updated by some explicit time stepping method; this implies that the method is only conditionally stable. For the numerical solution of the linear, elliptic pressure equation multigrid methods have become an accepted technique. On the other hand, the fully implicit method is unconditionally stable, but it has the disadvantage that in every time step a large system of nonlinear algebraic equations has to be solved. The most time-consuming part of any fully implicit reservoir simulator is the solution of this large system of equations. Usually this is done by Newton's method. The resulting systems of linear equations are then either solved by a direct method or by some conjugate gradient type method. In this paper we consider the possibility of applying multigrid methods for the iterative solution of the systems of nonlinear equations. There are two ways of using multigrid for this job: either we use a nonlinear multigrid method or we use a linear multigrid method to deal with the linear systems that arise in Newton's method. So far only a few authors have reported on the use of multigrid methods for fully implicit simulations. Two-level FAS algorithm is presented for the black-oil equations, and linear multigrid for
Reactive plasma upgrade of squalane - a heavy oil simulant
Kong, P.C.; Watkins, A.D.; Detering, B.A.; Thomas, C.P.
1995-10-01
U.S. light crude oil production has steadily declined over the last two decades. However, huge known heavy oil deposits in the North American continent remain largely untapped. In the past 10 years, the API gravity of crude oils has been decreasing by about 0.17% per year, and the sulfur content has been increasing by about 0.027% per year. As the API gravity of crude oil decreases, there will be an urgent need for economically viable new technologies to upgrade the heavy oil to a high API gravity feed stock for the refineries. The Idaho National Engineering Laboratory is investigating an innovative plasma process to upgrade heavy oil and refinery residuum. This paper will present some of the results and the implications of this technology for heavy oil upgrade and conversion.
Boriskina, Svetlana V; Sewell, Phillip; Benson, Trevor M; Nosich, Alexander I
2004-03-01
A fast and accurate method is developed to compute the natural frequencies and scattering characteristics of arbitrary-shape two-dimensional dielectric resonators. The problem is formulated in terms of a uniquely solvable set of second-kind boundary integral equations and discretized by the Galerkin method with angular exponents as global test and trial functions. The log-singular term is extracted from one of the kernels, and closed-form expressions are derived for the main parts of all the integral operators. The resulting discrete scheme has a very high convergence rate. The method is used in the simulation of several optical microcavities for modern dense wavelength-division-multiplexed systems. PMID:15005404
Loco, Daniele; Jurinovich, Sandro; Di Bari, Lorenzo; Mennucci, Benedetta
2016-01-14
We present and discuss a simple and fast computational approach to the calculation of electronic circular dichroism spectra of nucleic acids. It is based on a exciton model in which the couplings are obtained in terms of the full transition-charge distributions, as resulting from TDDFT methods applied on the individual nucleobases. We validated the method on two systems, a DNA G-quadruplex and a RNA β-hairpin whose solution structures have been accurately determined by means of NMR. We have shown that the different characteristics of composition and structure of the two systems can lead to quite important differences in the dependence of the accuracy of the simulation on the excitonic parameters. The accurate reproduction of the CD spectra together with their interpretation in terms of the excitonic composition suggest that this method may lend itself as a general computational tool to both predict the spectra of hypothetic structures and define clear relationships between structural and ECD properties. PMID:26646952
NASA Astrophysics Data System (ADS)
Grimminck, Dennis L. A. G.; Polman, Ben J. W.; Kentgens, Arno P. M.; Leo Meerts, W.
2011-08-01
A fast and accurate fit program is presented for deconvolution of one-dimensional solid-state quadrupolar NMR spectra of powdered materials. Computational costs of the synthesis of theoretical spectra are reduced by the use of libraries containing simulated time/frequency domain data. These libraries are calculated once and with the use of second-party simulation software readily available in the NMR community, to ensure a maximum flexibility and accuracy with respect to experimental conditions. EASY-GOING deconvolution ( EGdeconv) is equipped with evolutionary algorithms that provide robust many-parameter fitting and offers efficient parallellised computing. The program supports quantification of relative chemical site abundances and (dis)order in the solid-state by incorporation of (extended) Czjzek and order parameter models. To illustrate EGdeconv's current capabilities, we provide three case studies. Given the program's simple concept it allows a straightforward extension to include other NMR interactions. The program is available as is for 64-bit Linux operating systems.
NASA Astrophysics Data System (ADS)
Sepp Neves, Antonio Augusto; Pinardi, Nadia; martins, Flavio
2016-04-01
Every year, 270,000 tonnes of oil are estimated to be spilled in the ocean by vessel operations (e.g. tank washing, leakage of lubricants) and the so called operational spills are typically associated with small volumes and high occurrence rate. Vessel-related accidental spills (e.g. collisions, explosions) seldom occur and usually involve high volumes of oil, accounting for about 100,000 tonnes/year. The occurrence of accidental spills and their impacts have been well documented in the available literature. On the other hand, occurrence rates of operational spills and the effects they have on the marine and coastal environments remain very uncertain due to insufficient sampling effort and methodological limitations. Trying to foresee when and where an oil spill will occur in a certain area, its characteristics and impacts is, at present, impossible. Oil spill risk assessments (OSRAs) have been employed in several parts of the globe in order to deal with such uncertainties and protect the marine environment. In the present work, we computed the oil spill risk applying ensemble oil spill simulations following an ISO-31000 compliant OSRA methodology (Sepp Neves et al. , 2015). The ensemble experiment was carried out for the Algarve coast (southern Portugal) generating a unique data set of 51,200 numerical oil spill simulations covering the main sources of uncertainties (i.e. where and when the spill will happen and oil spill model configuration). From the generated data set, the risk due to accidental and operational spills was mapped for the Algarve municipalities based on the frequency and magnitude (i.e. concentrations) of beaching events and the main sources of risk were identified. The socioeconomic and environmental dimensions of the risk were treated separately. Seasonal changes in the risk index proposed due to the variability of meteo-oceanographic variables (i.e. currents and waves) were also quantified.
A scalable parallel black oil simulator on distributed memory parallel computers
NASA Astrophysics Data System (ADS)
Wang, Kun; Liu, Hui; Chen, Zhangxin
2015-11-01
This paper presents our work on developing a parallel black oil simulator for distributed memory computers based on our in-house parallel platform. The parallel simulator is designed to overcome the performance issues of common simulators that are implemented for personal computers and workstations. The finite difference method is applied to discretize the black oil model. In addition, some advanced techniques are employed to strengthen the robustness and parallel scalability of the simulator, including an inexact Newton method, matrix decoupling methods, and algebraic multigrid methods. A new multi-stage preconditioner is proposed to accelerate the solution of linear systems from the Newton methods. Numerical experiments show that our simulator is scalable and efficient, and is capable of simulating extremely large-scale black oil problems with tens of millions of grid blocks using thousands of MPI processes on parallel computers.
Models, Simulators, and Data-driven Resources for Oil and Natural Gas Research
NETL provides a number of analytical tools to assist in conducting oil and natural gas research. Software, developed under various DOE/NETL projects, includes numerical simulators, analytical models, databases, and documentation.[copied from http://www.netl.doe.gov/technologies/oil-gas/Software/Software_main.html] Links lead users to methane hydrates models, preedictive models, simulators, databases, and other software tools or resources.
NASA Astrophysics Data System (ADS)
Xiong, Yi
solutions and results of a commercial simulator before conducting numerical studies. The numerical studies demonstrate the effect of capillary pressure on VLE, and further on production performance. The significant effect of capillary pressure on VLE leads to the suppression of bubble-point pressure and more light components dissolved in the oil phase. Consequently it is observed that there is smaller gas saturation, larger mole fractions of light components, and faster pressure decreasing at reservoir conditions; meanwhile less gas and more oil are produced at surface. The substantial decrease in reservoir pore pressure results in a large increase of effective stress, which induces the changes of rock properties and influences the production performance. The stress-induced degradation of permeability undermines the production performance, and the geomechanical effect on the permeability of natural fractures is mainly responsible for the undermined production performance. The reduction of pore size due to the geomechanical effect could increase the capillary pressure, which enlarges the influence of capillarity on VLE and further suppresses bubble-point pressure. On the other hand, the effect of capillary pressure on VLE influences the fluid flow and therefore influences the effective stress through the flow-stress coupling process. Thus the interaction between pore confinement and rock compaction can be modeled with MSFLOW_COM, and illustrated through numerical studies. This research provides a three-dimensional numerical tool for accurately modeling porous and fractured tight oil reservoirs. The developed simulator is able to assist scientists and engineers to study and understand the complex multiphase, multi-component fluid flow behaviors in tight oil reservoirs.
NASA Astrophysics Data System (ADS)
Jiang, Xikai; Karpeev, Dmitry; Li, Jiyuan; de Pablo, Juan; Hernandez-Ortiz, Juan; Heinonen, Olle
Boundary integrals arise in many electrostatic and magnetostatic problems. In computational modeling of these problems, although the integral is performed only on the boundary of a domain, its direct evaluation needs O(N2) operations, where N is number of unknowns on the boundary. The O(N2) scaling impedes a wider usage of the boundary integral method in scientific and engineering communities. We have developed a parallel computational approach that utilize the Fast Multipole Method to evaluate the boundary integral in O(N) operations. To demonstrate the accuracy, efficiency, and scalability of our approach, we consider two test cases. In the first case, we solve a boundary value problem for a ferroelectric/ferromagnetic volume in free space using a hybrid finite element-boundary integral method. In the second case, we solve an electrostatic problem involving the polarization of dielectric objects in free space using the boundary element method. The results from test cases show that our parallel approach can enable highly efficient and accurate simulations of mesoscale electrostatic/magnetostatic problems. Computing resources was provided by Blues, a high-performance cluster operated by the Laboratory Computing Resource Center at Argonne National Laboratory. Work at Argonne was supported by U. S. DOE, Office of Science under Contract No. DE-AC02-06CH11357.
Accurate path integral molecular dynamics simulation of ab-initio water at near-zero added cost
NASA Astrophysics Data System (ADS)
Elton, Daniel; Fritz, Michelle; Soler, José; Fernandez-Serra, Marivi
It is now established that nuclear quantum motion plays an important role in determining water's structure and dynamics. These effects are important to consider when evaluating DFT functionals and attempting to develop better ones for water. The standard way of treating nuclear quantum effects, path integral molecular dynamics (PIMD), multiplies the number of energy/force calculations by the number of beads, which is typically 32. Here we introduce a method whereby PIMD can be incorporated into a DFT molecular dynamics simulation at virtually zero cost. The method is based on the cluster (many body) expansion of the energy. We first subtract the DFT monomer energies, using a custom DFT-based monomer potential energy surface. The evolution of the PIMD beads is then performed using only the more-accurate Partridge-Schwenke monomer energy surface. The DFT calculations are done using the centroid positions. Various bead thermostats can be employed to speed up the sampling of the quantum ensemble. The method bears some resemblance to multiple timestep algorithms and other schemes used to speed up PIMD with classical force fields. We show that our method correctly captures some of key effects of nuclear quantum motion on both the structure and dynamics of water. We acknowledge support from DOE Award No. DE-FG02-09ER16052 (D.E.) and DOE Early Career Award No. DE-SC0003871 (M.V.F.S.).
Gao, Yue; Jia, Yu-Yan; Gao, Bao-Yu; Cao, Bai-Chuan; Zhang, Yong-Qiang; Lu, Lei
2010-10-01
A dithiocarbamate flocculant, DTC (T403), was prepared by the reaction of amine-terminated polyoxypropane-ether compound known as Jeffamine-T403 and carbon disulfide in alkaline solution. The oil removal efficiency of DTC (T403) for simulated produced wastewater from polymer flooding in oil production was studied by Jar-test. The effect of the dosage of DTC (T403), hydrolyzed polyacrylamide (HPAM), Fe2+ and Fe3+ ions, and pH on the oil removal efficiency of DTC (T403) was investigated. The results showed that the chelate polymer formed by DTC (T403) and Fe2+ ion has good oil removal performance by net capturing mechanism. HPAM had a negative effect on oil removal efficiency of DTC (T403). For the treatment of the simulated wastewater containing 0-900 mg/L of HPAM and 300 mg/L of oil, the residual oil concentrations in water samples decreased below 10 mg/L when the dosage of Fe2+ and DTC (T403) was 10 mg/L and 25 mg/L, respectively. The oil removal efficiency of DTC (T403) was affected by pH and good oil removal efficiency was obtained when the pH was below 7.5. DTC (T403) is appropriate for the treatment of oily wastewater containing Fe2+ ion. PMID:21229744
NASA Astrophysics Data System (ADS)
Jamil, I. E.; Al-Kayiem, H. H.
2015-12-01
The maturing oil fields with increasing water production can pose a challenge in terms of produced water handling and disposal issues. This paper presents the modelling and simulation procedure of the two-phase flow of water/oil in a downhole using ANSYS- FLUENT 14 software. The developed procedure successfully simulated the production zone and the interaction of the two fluids in a natural environment where the reservoir pressure is the sole driving force. The results show significant difference of volume distribution in the flows with different oil content. The model can become an essential tool to assist in prediction of the behavior of oil/water mixture flow in the wellbores, and to serve in designing downhole oil/water separators.
Liu, Mengge; Chen, Guang; Guo, Hailong; Fan, Baolei; Liu, Jianjun; Fu, Qiang; Li, Xiu; Lu, Xiaomin; Zhao, Xianen; Li, Guoliang; Sun, Zhiwei; Xia, Lian; Zhu, Shuyun; Yang, Daoshan; Cao, Ziping; Wang, Hua; Suo, Yourui; You, Jinmao
2015-09-16
Determination of plant growth regulators (PGRs) in a signal transduction system (STS) is significant for transgenic food safety, but may be challenged by poor accuracy and analyte instability. In this work, a microwave-assisted extraction-derivatization (MAED) method is developed for six acidic PGRs in oil samples, allowing an efficient (<1.5 h) and facile (one step) pretreatment. Accuracies are greatly improved, particularly for gibberellin A3 (-2.72 to -0.65%) as compared with those reported (-22 to -2%). Excellent selectivity and quite low detection limits (0.37-1.36 ng mL(-1)) are enabled by fluorescence detection-mass spectrum monitoring. Results show the significant differences in acidic PGRs between transgenic and nontransgenic oils, particularly 1-naphthaleneacetic acid (1-NAA), implying the PGRs induced variations of components and genes. This study provides, for the first time, an accurate and efficient determination for labile PGRs involved in STS and a promising concept for objectively evaluating the safety of transgenic foods. PMID:26309068
NASA Astrophysics Data System (ADS)
Dougherty, N. S.; Burnette, D. W.; Holt, J. B.; Matienzo, Jose
1993-07-01
Time-accurate unsteady flow simulations are being performed supporting the SRM T+68sec pressure 'spike' anomaly investigation. The anomaly occurred in the RH SRM during the STS-54 flight (STS-54B) but not in the LH SRM (STS-54A) causing a momentary thrust mismatch approaching the allowable limit at that time into the flight. Full-motor internal flow simulations using the USA-2D axisymmetric code are in progress for the nominal propellant burn-back geometry and flow conditions at T+68-sec--Pc = 630 psi, gamma = 1.1381, T(sub c) = 6200 R, perfect gas without aluminum particulate. In a cooperative effort with other investigation team members, CFD-derived pressure loading on the NBR and castable inhibitors was used iteratively to obtain nominal deformed geometry of each inhibitor, and the deformed (bent back) inhibitor geometry was entered into this model. Deformed geometry was computed using structural finite-element models. A solution for the unsteady flow has been obtained for the nominal flow conditions (existing prior to the occurrence of the anomaly) showing sustained standing pressure oscillations at nominally 14.5 Hz in the motor IL acoustic mode that flight and static test data confirm to be normally present at this time. Average mass flow discharged from the nozzle was confirmed to be the nominal expected (9550 lbm/sec). The local inlet boundary condition is being perturbed at the location of the presumed reconstructed anomaly as identified by interior ballistics performance specialist team members. A time variation in local mass flow is used to simulate sudden increase in burning area due to localized propellant grain cracks. The solution will proceed to develop a pressure rise (proportional to total mass flow rate change squared). The volume-filling time constant (equivalent to 0.5 Hz) comes into play in shaping the rise rate of the developing pressure 'spike' as it propagates at the speed of sound in both directions to the motor head end and nozzle. The
NASA Technical Reports Server (NTRS)
Dougherty, N. S.; Burnette, D. W.; Holt, J. B.; Matienzo, Jose
1993-01-01
Time-accurate unsteady flow simulations are being performed supporting the SRM T+68sec pressure 'spike' anomaly investigation. The anomaly occurred in the RH SRM during the STS-54 flight (STS-54B) but not in the LH SRM (STS-54A) causing a momentary thrust mismatch approaching the allowable limit at that time into the flight. Full-motor internal flow simulations using the USA-2D axisymmetric code are in progress for the nominal propellant burn-back geometry and flow conditions at T+68-sec--Pc = 630 psi, gamma = 1.1381, T(sub c) = 6200 R, perfect gas without aluminum particulate. In a cooperative effort with other investigation team members, CFD-derived pressure loading on the NBR and castable inhibitors was used iteratively to obtain nominal deformed geometry of each inhibitor, and the deformed (bent back) inhibitor geometry was entered into this model. Deformed geometry was computed using structural finite-element models. A solution for the unsteady flow has been obtained for the nominal flow conditions (existing prior to the occurrence of the anomaly) showing sustained standing pressure oscillations at nominally 14.5 Hz in the motor IL acoustic mode that flight and static test data confirm to be normally present at this time. Average mass flow discharged from the nozzle was confirmed to be the nominal expected (9550 lbm/sec). The local inlet boundary condition is being perturbed at the location of the presumed reconstructed anomaly as identified by interior ballistics performance specialist team members. A time variation in local mass flow is used to simulate sudden increase in burning area due to localized propellant grain cracks. The solution will proceed to develop a pressure rise (proportional to total mass flow rate change squared). The volume-filling time constant (equivalent to 0.5 Hz) comes into play in shaping the rise rate of the developing pressure 'spike' as it propagates at the speed of sound in both directions to the motor head end and nozzle. The
Oil spills and AI: How to manage resources through simulation
Giribone, P.; Bruzzone, A.G.; Caddeo, S.
1995-12-31
Today, in the Mediterranean theater of the Upper Tyrrhenian, the ecological risk involving oil installations is still quite high. This is due to the fact that valuable environmental and tourist areas exist together with large industrial and port structures; in particular, recent events have demonstrated the danger involving oil spills along the Ligurian coastline. This study proposes an approach to plan the operations that should be performed when accidents occur, based on the use of AI techniques.
Calculation of high-temperature crude oil/water/vapor separations using simulated distillation data
Langhoff, J.A.; Wu, C.H.
1984-09-01
High-temperature crude oil-water-vapor separation takes place in steam flooding and in-situ combustion processes. It also takes place in hydrocarbon recovery from deep volatile oil reservoirs. A practical procedure using the Holland and Welch method and simulated distillation data was developed to calculate crude oil-water-vapor separations at 387/sup 0/F (197/sup 0/C) and 456/sup 0/F (235/sup 0/C). The overhead yields obtained from the calculations were expressed as a function of the steam distillation factor, V /SUB w/ /V /SUB oi/ . The results were compared with laboratory crude oil steam distillation data. The approach satisfactorily predicted the overhead yields of thirteen out of sixteen crude oils with an average error of 11%. This is within experimental error of crude oil steam distillation. Twelve pseudocomponents of crude oils were selected and characterized using simulated distillation data for the calculations. The physical properties of the pseudocomponents were determined from existing correlations and from matching laboratory steam distillation data. The use of simulated distillation data eliminates the uncertainty and assumptions normally involved in the selection of crude oil pseudocomponents using Bureau of Mines distillation data, and thus, improves the reliability of the proposed computational approach. The proposed approach has many advantages. It eliminates the need of conducting experimental steam distillation tests if crude oil simulated distillation data are available. It is easy and fast to calculate the overhead yields and densities without using the equation of state and uncertain pseudocomponent critical properties and interaction parameters. The proposed approach will provide useful information for designing and operating thermal recovery processes, and for predicting hydrocarbon recovery from high-temperature volatile oil reservoirs.
Underwater Oil Plume Intrusion from Deepwater Blowouts - A Large-Eddy Simulation Study
NASA Astrophysics Data System (ADS)
Yang, D.; Chen, B.; Chamecki, M.; Meneveau, C. V.
2015-12-01
The interaction of buoyancy-driven hydrocarbon plumes with the stably stratified deep-ocean environment plays a crucial role in the formation of underwater oil intrusions. As gas bubbles and oil droplets are released from an underwater oil well blowout, they induce a strong buoyancy flux that lifts entrained sea water to form an upward plume. Towards higher elevations, the stratification-induced negative buoyancy increases and eventually exceeds the gas/oil-induced buoyancy, causing the plume to decelerate and a large fraction of entrained sea water to peel off from the rising plume to form a fountain-like downward outer plume. During this peeling process, weakly buoyant particles (e.g. small oil droplets) are trapped and fall together with the detrained fluid, and then migrate horizontally at the equilibrium buoyancy depth, forming underwater oil intrusion layers. In this study, the complex plume dynamics and oil intrusion are studied using a large-eddy simulation (LES) model. The LES model captures the essential characteristics of the plume structure and the peeling/intrusion processes, and yields good agreement with prior laboratory experiments. Applying to the Deepwater Horizon oil well blowout condition, the LES model shows considerable underwater trapping and intrusion of oil droplets under various conditions, with the trapping rate significantly affected by the diameter of the oil droplet. This study is supported by Gulf of Mexico Research Initiative RFP-II research grant.
Experimental and numerical simulation study of microbial enhanced oil recovery using bio-surfactants
NASA Astrophysics Data System (ADS)
Maudgalya, Saikrishna
An experimental and numerical study were conducted to investigate the ability of bio-surfactant produced by the microbe Bacillus mojavensis strain JF-2 to recover residual oil from consolidated porous media. Experiments showed that the bio-surfactant at concentrations as low as 40.0 ppm. (0.04 mg/scc) and viscosified with 1000.0 ppm of polymer could recover 10.0 % to 40.0 % of residual oil when injected through sandstone cores at typical field rates. A 2-phase, 10-component microbial enhanced oil recovery numerical simulator was modified to include reservoir salinity and facilitate surfactant and polymer injection. The effects of reservoir brine salinity and divalent ion effects on bio-surfactant and polymer adsorption, polymer retention, polymer viscosity, bio-surfactant interfacial tension and the shear rate effect on polymer viscosity were added to the simulator. Core flood experiments where JF-2 bio-surfactant viscosified with partially hydrolyzed polyacrylamide was injected into Berea cores at waterflood residual oil saturation were simulated. The effects of brine salinity and hardness on surfactant and polymer behavior were tested and the core flood simulation results compared with the experimental results. After the laboratory and simulation studies, a residual oil recovery method based on non-aqueous phase liquid (NAPL) contaminant removal from aquifers is discussed and functional form of the transport equation presented. In this method, residual oil is treated as another chemical species dispersed in porous media instead of a phase that is uniformly distributed across the media.
NASA Astrophysics Data System (ADS)
Bozinoski, Radoslav
Significant research has been performed over the last several years on understanding the unsteady aerodynamics of various fluid flows. Much of this work has focused on quantifying the unsteady, three-dimensional flow field effects which have proven vital to the accurate prediction of many fluid and aerodynamic problems. Up until recently, engineers have predominantly relied on steady-state simulations to analyze the inherently three-dimensional ow structures that are prevalent in many of today's "real-world" problems. Increases in computational capacity and the development of efficient numerical methods can change this and allow for the solution of the unsteady Reynolds-Averaged Navier-Stokes (RANS) equations for practical three-dimensional aerodynamic applications. An integral part of this capability has been the performance and accuracy of the turbulence models coupled with advanced parallel computing techniques. This report begins with a brief literature survey of the role fully three-dimensional, unsteady, Navier-Stokes solvers have on the current state of numerical analysis. Next, the process of creating a baseline three-dimensional Multi-Block FLOw procedure called MBFLO3 is presented. Solutions for an inviscid circular arc bump, laminar at plate, laminar cylinder, and turbulent at plate are then presented. Results show good agreement with available experimental, numerical, and theoretical data. Scalability data for the parallel version of MBFLO3 is presented and shows efficiencies of 90% and higher for processes of no less than 100,000 computational grid points. Next, the description and implementation techniques used for several turbulence models are presented. Following the successful implementation of the URANS and DES procedures, the validation data for separated, non-reattaching flows over a NACA 0012 airfoil, wall-mounted hump, and a wing-body junction geometry are presented. Results for the NACA 0012 showed significant improvement in flow predictions
NASA Astrophysics Data System (ADS)
Mehmani, Yashar; Oostrom, Mart; Balhoff, Matthew T.
2014-03-01
Several approaches have been developed in the literature for solving flow and transport at the pore scale. Some authors use a direct modeling approach where the fundamental flow and transport equations are solved on the actual pore-space geometry. Such direct modeling, while very accurate, comes at a great computational cost. Network models are computationally more efficient because the pore-space morphology is approximated. Typically, a mixed cell method (MCM) is employed for solving the flow and transport system which assumes pore-level perfect mixing. This assumption is invalid at moderate to high Peclet regimes. In this work, a novel Eulerian perspective on modeling flow and transport at the pore scale is developed. The new streamline splitting method (SSM) allows for circumventing the pore-level perfect-mixing assumption, while maintaining the computational efficiency of pore-network models. SSM was verified with direct simulations and validated against micromodel experiments; excellent matches were obtained across a wide range of pore-structure and fluid-flow parameters. The increase in the computational cost from MCM to SSM is shown to be minimal, while the accuracy of SSM is much higher than that of MCM and comparable to direct modeling approaches. Therefore, SSM can be regarded as an appropriate balance between incorporating detailed physics and controlling computational cost. The truly predictive capability of the model allows for the study of pore-level interactions of fluid flow and transport in different porous materials. In this paper, we apply SSM and MCM to study the effects of pore-level mixing on transverse dispersion in 3-D disordered granular media.
NASA Astrophysics Data System (ADS)
Sagui, Celeste; Pedersen, Lee G.; Darden, Thomas A.
2004-01-01
The accurate simulation of biologically active macromolecules faces serious limitations that originate in the treatment of electrostatics in the empirical force fields. The current use of "partial charges" is a significant source of errors, since these vary widely with different conformations. By contrast, the molecular electrostatic potential (MEP) obtained through the use of a distributed multipole moment description, has been shown to converge to the quantum MEP outside the van der Waals surface, when higher order multipoles are used. However, in spite of the considerable improvement to the representation of the electronic cloud, higher order multipoles are not part of current classical biomolecular force fields due to the excessive computational cost. In this paper we present an efficient formalism for the treatment of higher order multipoles in Cartesian tensor formalism. The Ewald "direct sum" is evaluated through a McMurchie-Davidson formalism [L. McMurchie and E. Davidson, J. Comput. Phys. 26, 218 (1978)]. The "reciprocal sum" has been implemented in three different ways: using an Ewald scheme, a particle mesh Ewald (PME) method, and a multigrid-based approach. We find that even though the use of the McMurchie-Davidson formalism considerably reduces the cost of the calculation with respect to the standard matrix implementation of multipole interactions, the calculation in direct space remains expensive. When most of the calculation is moved to reciprocal space via the PME method, the cost of a calculation where all multipolar interactions (up to hexadecapole-hexadecapole) are included is only about 8.5 times more expensive than a regular AMBER 7 [D. A. Pearlman et al., Comput. Phys. Commun. 91, 1 (1995)] implementation with only charge-charge interactions. The multigrid implementation is slower but shows very promising results for parallelization. It provides a natural way to interface with continuous, Gaussian-based electrostatics in the future. It is
Mehmani, Yashar; Oostrom, Martinus; Balhoff, Matthew
2014-03-20
Several approaches have been developed in the literature for solving flow and transport at the pore-scale. Some authors use a direct modeling approach where the fundamental flow and transport equations are solved on the actual pore-space geometry. Such direct modeling, while very accurate, comes at a great computational cost. Network models are computationally more efficient because the pore-space morphology is approximated. Typically, a mixed cell method (MCM) is employed for solving the flow and transport system which assumes pore-level perfect mixing. This assumption is invalid at moderate to high Peclet regimes. In this work, a novel Eulerian perspective on modeling flow and transport at the pore-scale is developed. The new streamline splitting method (SSM) allows for circumventing the pore-level perfect mixing assumption, while maintaining the computational efficiency of pore-network models. SSM was verified with direct simulations and excellent matches were obtained against micromodel experiments across a wide range of pore-structure and fluid-flow parameters. The increase in the computational cost from MCM to SSM is shown to be minimal, while the accuracy of SSM is much higher than that of MCM and comparable to direct modeling approaches. Therefore, SSM can be regarded as an appropriate balance between incorporating detailed physics and controlling computational cost. The truly predictive capability of the model allows for the study of pore-level interactions of fluid flow and transport in different porous materials. In this paper, we apply SSM and MCM to study the effects of pore-level mixing on transverse dispersion in 3D disordered granular media.
Simulation of impact of oil spill in the ocean--a case study of Arabian Gulf.
Verma, Parikshit; Wate, Satish R; Devotta, Sukumar
2008-11-01
To meet the growing energy demand worldwide, oil and gas exploration and production activities have increased rapidly both in onshore and offshore areas. The produced oil from the ocean bed is transported onshore either by ship or pipeline. This has increased the risk of oil spill in the coastal area. In order to prepare an emergency preparedness plan and to assess the magnitude of risk involved in transporting and offloading oil, oil spill simulation studies play an important role. This paper describes a simulation of oil spill in coastal bay of Arabian Gulf where new developments are taking place using MIKE 21 model. The developments include a diesel based thermal power plant near Sir Baniyas Island, which is an ecological fragile area. Based on the project activity, two probable scenarios, one for diesel leak (250 m3/h) for 1 h and the other for instantaneous spill (500 m3) are considered. The MIKE 21 model was calibrated for hydrodynamics using measured field data followed by diesel-spill simulation to track its movement in the Arabian Gulf. The results for both leak and instantaneous spill indicate that spilled diesel will not move towards the Sir Banyas Island and more than 45% of the diesel will be evaporated within 48 h of oil spill. Based on the results, a clean up and contingency plan is proposed to mitigate the adverse impacts arising due to diesel spill in the study area. PMID:18095178
Simulation of Anomalous Oil Filtration in a Porous Bed
NASA Astrophysics Data System (ADS)
Kelbaliev, G. I.; Rzaev, Ab. G.; Rasulov, S. R.; Guseinova, L. V.
2015-03-01
The problems of modeling the processes of filtration of anomalous structurized oils with coagulation structures present in a porous bed are considered. An equation for the filtration of Bingham fluids in a carrier bed that accounts for the dependence of the filtration rate on the pressure gradient and shear stress has been derived. Models for calculating the effective viscosity and mobility of oil depending on the change of the pressure gradient in the bed and of the concentration of particles in it have been developed. A comparison of these models with the experimental data available in the literature for various wells yielded satisfactory results.
NASA Astrophysics Data System (ADS)
Sepp Neves, Antonio Augusto; Pinardi, Nadia; Martins, Flavio
2016-08-01
Oil Spill Risk Assessments (OSRAs) are widely employed to support decision making regarding oil spill risks. This article adapts the ISO-compliant OSRA framework developed by Sepp Neves et al. (J Environ Manag 159:158-168, 2015) to estimate risks in a complex scenario where uncertainties related to the meteo-oceanographic conditions, where and how a spill could happen exist and the risk computation methodology is not yet well established (ensemble oil spill modeling). The improved method was applied to the Algarve coast, Portugal. Over 50,000 simulations were performed in 2 ensemble experiments to estimate the risks due to operational and accidental spill scenarios associated with maritime traffic. The level of risk was found to be important for both types of scenarios, with significant seasonal variations due to the the currents and waves variability. Higher frequency variability in the meteo-oceanographic variables were also found to contribute to the level of risk. The ensemble results show that the distribution of oil concentrations found on the coast is not Gaussian, opening up new fields of research on how to deal with oil spill risks and related uncertainties.
NASA Astrophysics Data System (ADS)
Sepp Neves, Antonio Augusto; Pinardi, Nadia; Martins, Flavio
2016-06-01
Oil Spill Risk Assessments (OSRAs) are widely employed to support decision making regarding oil spill risks. This article adapts the ISO-compliant OSRA framework developed by Sepp Neves et al. (J Environ Manag 159:158-168, 2015) to estimate risks in a complex scenario where uncertainties related to the meteo-oceanographic conditions, where and how a spill could happen exist and the risk computation methodology is not yet well established (ensemble oil spill modeling). The improved method was applied to the Algarve coast, Portugal. Over 50,000 simulations were performed in 2 ensemble experiments to estimate the risks due to operational and accidental spill scenarios associated with maritime traffic. The level of risk was found to be important for both types of scenarios, with significant seasonal variations due to the the currents and waves variability. Higher frequency variability in the meteo-oceanographic variables were also found to contribute to the level of risk. The ensemble results show that the distribution of oil concentrations found on the coast is not Gaussian, opening up new fields of research on how to deal with oil spill risks and related uncertainties.
Modeling and simulation of an enzymatic reactor for hydrolysis of palm oil.
Bhatia, S; Naidu, A D; Kamaruddin, A H
1999-01-01
Hydrolysis of palm oil has become an important process in Oleochemical industries. Therefore, an investigation was carried out for hydrolysis of palm oil to fatty acid and glycerol using immobilized lipase in packed bed reactor. The conversion vs. residence time data were used in Michaelis-Menten rate equation to evaluate the kinetic parameters. A mathematical model for the rate of palm oil hydrolysis was proposed incorporating role of external mass transfer and pore diffusion. The model was simulated for steady-state isothermal operation of immobilized lipase packed bed reactor. The experimental data were compared with the simulated results. External mass transfer was found to affect the rate of palm oil hydrolysis at higher residence time. PMID:10595445
Computer simulation of the probability that endangered whales will interact with oil spills
Reed, M.; Jayko, K.; Bowles, A.; Anderson, E.; Leatherwood, S.
1987-03-01
A numerical model system was developed to assess quantitatively the probability that endangered bowhead and gray whales will encounter spilled oil in Alaskan waters. Bowhead and gray whale migration and diving-surfacing models, and an oil-spill trajectory model comprise the system. The migration models were developed from conceptual considerations, then calibrated with and tested against observations. The movement of a whale point is governed by a random walk algorithm which stochastically follows a migratory pathway. The oil-spill model, developed under a series of other contracts, accounts for transport and spreading behavior in open water and in the presence of sea ice. Historical wind records and heavy, normal, or light ice cover data sets are selected at random to provide stochastic oil-spill scenarios for whale-oil interaction simulations.
BOAST-VHS. 3D 3-Phase Black Oil Applied Simulator
Chang, M.; Sarathi, P.; Heemstra, R.J.; Cheng, A.M.; Pautz, J.F.
1992-01-01
BOAST-VHS is a three-dimensional, three-phase, finite-difference black oil simulator developed for use on a personal computer. It simulates isothermal, Darcy flow in three dimensions. The simulator assumes that the reservoir fluids can be described by three fluid phases (oil, gas, and water) of constant composition whose properties are functions of pressure only. BOAST-VHS can simulate oil and/or gas recovery by fluid expansion, displacement, gravity drainage, and imbibition mechanisms. BOAST-VHS is recommended as a cost-effective reservoir simulation tool for the study of such problems as primary depletion, pressure maintenance (by water and/or gas injection) and basic secondary recovery operations (such as waterflooding) in a black oil reservoir using slanted or horizontal wells, in addition to conventional vertical wells. The well model in BOAST-VHS permits specification of rate or pressure constraints on well performance. The model also allows the user to add or recomplete wells during the period represented by the simulation. BOAST-VHS has flexible initialization capabilities, a bubble point tracking scheme, an automatic time-step control method, a zero transmissibility option (inactive grid blocks), and a material balance check on solution stability.
A dynamic model to simulate foamy oil flow in porous media
Sheng, J.J.; Hayes, R.E.; Maini, B.B.; Tortike, W.S.
1996-12-31
A number of heavy oil reservoirs under solution gas drive show anomalously good primary performance. Foamy oil behavior is believed to be one of the reasons. Currently, numerical simulation of primary depletion in foamy oil reservoirs is based primarily on empirical adjustments to the conventional solution gas drive models. This paper presents a numerical model including the rate processes related to the nucleation, bubble growth and coalescence. The rate of nucleation is assumed to be instantaneous; the rate of bubble growth is a function of supersaturation and time; and the rate of coalescence is proportional to the amount of gas bubbles dispersed in the oil phase. The two-phase flow of oil and gas is modelled with the normal two-phase relative permeability-saturation relationship. The dispersed gas is assumed to flow with the oil as if it was a part of the liquid phase. The model satisfactorily matched the primary depletion tests in a sand-pack. It was observed that the solution gas drive recovery factor increased dramatically as the rate of pressure decline was increased. This model adequately accounts for the rate processes under solution gas drive in heavy oil reservoirs. It can be extended to investigate the effects of various process parameters on oil recovery and it may provide more reliable prediction of field performance. It also provides a tool to evaluate the significance of dynamic processes under various operation conditions.
Simulating oil droplet dispersal from the Deepwater Horizon spill with a Lagrangian approach
North, Elizabeth W.; Adams, E. Eric; Sherwood, Christopher R.; He, Ruoying; Hyun, Kyung Hoon; Socolofsky, Scott A.; Schlag, Zachary
2011-01-01
An analytical multiphase plume model, combined with time-varying flow and hydrographic fields generated by the 3-D South Atlantic Bight and Gulf of Mexico model (SABGOM) hydrodynamic model, were used as input to a Lagrangian transport model (LTRANS), to simulate transport of oil droplets dispersed at depth from the recent Deepwater Horizon MC 252 oil spill. The plume model predicts a stratification-dominated near field, in which small oil droplets detrain from the central plume containing faster rising large oil droplets and gas bubbles and become trapped by density stratification. Simulated intrusion (trap) heights of ∼ 310–370 m agree well with the midrange of conductivity-temperature-depth observations, though the simulated variation in trap height was lower than observed, presumably in part due to unresolved variability in source composition (percentage oil versus gas) and location (multiple leaks during first half of spill). Simulated droplet trajectories by the SABGOM-LTRANS modeling system showed that droplets with diameters between 10 and 50 μm formed a distinct subsurface plume, which was transported horizontally and remained in the subsurface for >1 month. In contrast, droplets with diameters ≥90 μm rose rapidly to the surface. Simulated trajectories of droplets ≤50 μm in diameter were found to be consistent with field observations of a southwest-tending subsurface plume in late June 2010 reported by Camilli et al. [2010]. Model results suggest that the subsurface plume looped around to the east, with potential subsurface oil transport to the northeast and southeast. Ongoing work is focusing on adding degradation processes to the model to constrain droplet dispersal.
Numerical Simulation and Structural Optimization of the Inclined Oil/Water Separator
Wu, Shijuan; Lu, Hongfang; Zhao, Lijie
2015-01-01
Improving the separation efficiency of the inclined oil/water separator, a new type of gravity separation equipment, is of great importance. In order to obtain a comprehensive understanding of the internal flow field of the separation process of oil and water within this separator, a numerical simulation based on Euler multiphase flow analysis and the realizable k-ε two equation turbulence model was executed using Fluent software. The optimal value ranges of the separator’s various structural parameters used in the numerical simulation were selected through orthogonal array experiments. A field experiment on the separator was conducted with optimized structural parameters in order to validate the reliability of the numerical simulation results. The research results indicated that the horizontal position of the dispenser, the hole number, and the diameter had significant effects on the oil/water separation efficiency, and that the longitudinal position of the dispenser and the position of the weir plate had insignificant effects on the oil/water separation efficiency. The optimal structural parameters obtained through the orthogonal array experiments resulted in an oil/water separation efficiency of up to 95%, which was 4.996% greater than that realized by the original structural parameters. PMID:25874861
Oil and air dispersion in a simulated fermentation broth as a function of mycelial morphology.
Lucatero, Savidra; Larralde-Corona, Claudia Patricia; Corkidi, Gabriel; Galindo, Enrique
2003-01-01
The culture conditions of a multiphase fermentation involving morphologically complex mycelia were simulated in order to investigate the influence of mycelial morphology (Trichoderma harzianum) on castor oil and air dispersion. Measurements of oil drops and air bubbles were obtained using an image analysis system coupled to a mixing tank. Complex interactions of the phases involved could be clearly observed. The Sauter diameter and the size distributions of drops and bubbles were affected by the morphological type of biomass (pellets or dispersed mycelia) added to the system. Larger oil drop sizes were obtained with dispersed mycelia than with pellets, as a result of the high apparent viscosity of the broth, which caused a drop in the power drawn, reducing oil drop break-up. Unexpectedly, bubble sizes observed with dispersed mycelia were smaller than with pellets, a phenomenon which can be explained by the segregation occurring at high biomass concentrations with the dispersed mycelia. Very complex oil drops were produced, containing air bubbles and a high number of structures likely consisting of small water droplets. Bubble location was influenced by biomass morphology. The percentage (in volume) of oil-trapped bubbles increased (from 32 to 80%) as dispersed mycelia concentration increased. A practically constant (32%) percentage of oil-trapped bubbles was observed with pelleted morphology at all biomass concentrations. The results evidenced the high complexity of phases interactions and the importance of mycelial morphology in such processes. PMID:12675561
A framework to quantify uncertainty in simulations of oil transport in the ocean
NASA Astrophysics Data System (ADS)
Gonçalves, Rafael C.; Iskandarani, Mohamed; Srinivasan, Ashwanth; Thacker, W. Carlisle; Chassignet, Eric; Knio, Omar M.
2016-04-01
An uncertainty quantification framework is developed for the DeepC Oil Model based on a nonintrusive polynomial chaos method. This allows the model's output to be presented in a probabilistic framework so that the model's predictions reflect the uncertainty in the model's input data. The new capability is illustrated by simulating the far-field dispersal of oil in a Deepwater Horizon blowout scenario. The uncertain input consisted of ocean current and oil droplet size data and the main model output analyzed is the ensuing oil concentration in the Gulf of Mexico. A 1331 member ensemble was used to construct a surrogate for the model which was then mined for statistical information. The mean and standard deviations in the oil concentration were calculated for up to 30 days, and the total contribution of each input parameter to the model's uncertainty was quantified at different depths. Also, probability density functions of oil concentration were constructed by sampling the surrogate and used to elaborate probabilistic hazard maps of oil impact. The performance of the surrogate was constantly monitored in order to demarcate the space-time zones where its estimates are reliable.
Pulido-Mayoral, Nancy; Galindo, Enrique
2004-01-01
The sizes of air bubbles and castor oil drops were studied by image analysis as a function of the concentration of soluble protein (bovine serum albumin [BSA] and lipase, as model proteins) in a three-phase system using a simulated fermentation medium (aqueous salt solution, castor oil, and air). Small amounts of proteins (<0.02 g/L) caused an important decrease in oil drops and bubbles sizes, together with a pronounced decrease in surface tension. The extent and profiles of this decrease seem to be determined by the conformation of the protein at the interface. The kLa value increased considerably for increasing concentration (up to 0.02 g/L) of the two proteins but was very different (2-fold higher for the lipase) at the highest concentrations tested (0.5 g/L), a phenomenon that can be caused by the extent to which bubbles are trapped within oil drops. PMID:15458353
Chocholousová, Jana; Feig, Michael
2006-04-30
Different integrator time steps in NVT and NVE simulations of protein and nucleic acid systems are tested with the GBMV (Generalized Born using Molecular Volume) and GBSW (Generalized Born with simple SWitching) methods. The simulation stability and energy conservation is investigated in relation to the agreement with the Poisson theory. It is found that very close agreement between generalized Born methods and the Poisson theory based on the commonly used sharp molecular surface definition results in energy drift and simulation artifacts in molecular dynamics simulation protocols with standard 2-fs time steps. New parameters are proposed for the GBMV method, which maintains very good agreement with the Poisson theory while providing energy conservation and stable simulations at time steps of 1 to 1.5 fs. PMID:16518883
An evaluation of the effectiveness of the Texas Oil Spill Control, Inc. (TOSCON) weir skimmer and gravity differential separator was conducted at the U.S. Environmental Protection Agency's Oil and Hazardous Materials Simulated Environmental Test Tank (OHMSETT) facility in October...
Changes in the chemical components of light crude oil during simulated short term weathering
NASA Astrophysics Data System (ADS)
Ma, Qimin; Luo, Xia; Yu, Zhigang
2007-07-01
To unambiguously identify spilled oils and to link them to the known sources are extremely important in settling questions of environmental impact and legal liability. The fate and behavior of spilled oils in the environment depend on a number of physicochemical and biological factors. This paper presents the results regarding changes in chemical composition of light crude oil during simulated short-term weathering based on natural environmental conditions. The results show that the saturated hydrocarbons of the light crude oil mainly distribute between n-C8 and n-C23 and the most abundant n-alkanes are found in the n-C10 to n-C16. The main chemical components of the light crude oil are n-alkanes and isoprenoids. The aromatic compounds are subordinate chemical components. Under the conditions of the weathering simulation experiment, n-alkanes less than n-C12, toluene and 1,3-dimethyl benzene are lost after 1 d weathering, the n-C13, n-C14, naphthalene and 2-methyl-naphthalene are lost on the fifth day of weathering, and n-C15 alkane components show certain weatherproof capability. The ratios n-C17/pristane and n-C18/phytane are unaltered and can be used to identify the source of the light crude oil during the first 8 d of weathering. After 21 d, the ratio pristine/phytane can not provide much information on the source of the spilled light crude oil. Triterpanes (m/z 191) as biomarker compounds of light crude oil are more valuable.
NASA Astrophysics Data System (ADS)
González, J.; Figueiras, F. G.; Aranguren-Gassis, M.; Crespo, B. G.; Fernández, E.; Morán, X. A. G.; Nieto-Cid, M.
2009-07-01
Two microcosm experiments were carried out to simulate the effect of sporadic oil spills derived from tanker accidents on oceanic and coastal marine phytoplankton assemblages. Treatments were designed to reproduce the spill from the Prestige, which took place in Galician coastal waters (NW Iberia) in November 2002. Two different concentrations of the water soluble fraction of oil were used: low (8.6 ± 0.7 μg l -1 of chrysene equivalents) and high (23 ± 5 μg l -1 of chrysene equivalents l -1). Photosynthetic activity and chlorophyll a concentration decreased in both assemblages after 24-72 h of exposure to the two oil concentrations, even though the effect was more severe on the oceanic assemblage. These variables progressively recovered up to values close or higher than those in the controls, but the short-term negative effect of oil, which was generally stronger at the high concentration, also induced changes in the structure of the plankton community. While the biomass of nanoflagellates increased in both assemblages, oceanic picophytoplankton was drastically reduced by the addition of oil. Effects on diatoms were also observed, particularly in the coastal assemblage. The response of coastal diatoms to oil addition showed a clear dependence on size. Small diatoms (<20 μm) were apparently stimulated by oil, whereas diatoms >20 μm were only negatively affected by the high oil concentration. These differences, which could be partially due to indirect trophic interactions, might also be related to different sensitivity of species to PAHs. These results, in agreement with previous observations, additionally show that the negative effect of the water soluble fraction of oil on oceanic phytoplankton was stronger than on coastal phytoplankton.
NASA Astrophysics Data System (ADS)
Pokorna, Lucie; Kliegrova, Stanislava; Huth, Radan; Farda, Ales; Stepanek, Petr
2014-05-01
Regional climate models (RCM) are a useful tool for a simulation of surface climate with respect to conditions of individual regions. The need of the realistic representation of surface elements at the local scale is important particularly in terrain with complex orography. The Czech Republic with the mountain chains along its border and highlands as well as lowlands in the inland seems to be a good representation of such region. A good performance of the models in reproducing recent temporal and spatial distribution of temperature and precipitation can enhance our confidence in the changes projected for future climate conditions. In this study, we compare two versions of the RCM ALARO covering a 30-year climate period (1961-1990); a simulation with a common resolution 25-km and a simulation with a very high resolution 6-km. The ALARO-Climate RCM has been developed in recent years in the Czech Hydrometeorological Institute on the basis of the numerical weather prediction model ALADIN and is already operated at other five national meteorological services. Both presented simulations are driven by the ERA-40 reanalysis and run on the large pan-European integration domain ("ENSEMBLES / Euro-Cordex domain"). As the reference dataset we use technical homogenized series based on time series from stations in the Czech Republic interpolated to the same network as both model simulations but with real altitude of the grid points (GriSt). The seasonal and monthly values of mean, maximum and minimum temperature as well as precipitation amounts are examined. We display a spatial distribution of biases of seasonal means and the temporal distribution of biases based on monthly values with respect to the altitude for both simulations. The results indicate that a higher resolution of model tends to improve the simulation of present day climate, with larger improvements in areas affected by mountains.
Molecular dynamic simulation of asphaltene co-aggregation with humic acid during oil spill.
Zhu, Xinzhe; Chen, Daoyi; Wu, Guozhong
2015-11-01
Humic acid in water and sediment plays a key role in the fate and transport of the spilled oil, but little is known about its influence on the aggregation of heavy oil asphaltenes which is adverse for remediation. Molecular dynamic simulation was performed to characterize the co-aggregation of asphaltenes (continental model and Violanthrone-79 model) with Leonardite humic acid (LHA) at the toluene-water interface and in bulk water, respectively, to simulate the transport of asphaltenes from oil to water. At the toluene-water interface, a LHA layer tended to form and bind to the water by hydrogen bonding which provided a surface for the accumulation of asphaltenes by parallel or T-shape stacking. After entering the bulk water, asphaltene aggregates stacked in parallel were tightly sequestrated inside the inner cavity of LHA aggregates following surface adsorption and structure deformation. Asphaltene aggregation in water was 2-fold higher than at the toluene-water interface. The presence of LHA increased the intensity of asphaltene aggregation by up to 83% in bulk water while relatively less influence was observed at the toluene-water interface. Overall results suggested that the co-aggregation of asphaltene with humic acid should be incorporated to the current oil spill models for better interpreting the overall environmental risks of oil spill. PMID:26149857
NASA Astrophysics Data System (ADS)
Chamecki, M.; Yang, D.; Meneveau, C. V.
2013-12-01
Deep water blowouts generate plumes of oil droplets that rise through, and interact with various layers of the ocean. When plumes reach the ocean mixed layer (OML), the interactions among oil droplet plume, Ekman Spiral and Langmuir turbulence strongly affect the final rates of dilution and bio-degradation. The present study aims at developing a large-eddy simulation (LES) capability for the study of the physical distribution and dispersion of oil droplets under the action of physical oceanographic processes in the OML. In the current LES approach, the velocity and temperature fields are simulated using a hybrid pseudo-spectral and finite-difference scheme; the oil field is described by an Eulerian concentration field and it is simulated using a bounded finite-volume scheme. Fluid accelerations induced by buoyancy of the oil plume are included, and a number of subgrid-scale models for the flow solver are implemented and tested. The LES capability is then applied to the simulation of oil plume dispersion in the OML. Graphical visualization of the LES results shows surface oil slick distribution consistent with the satellite and aerial images of surface oil slicks reported in the literature. Different combinations of Lamgmuir turbulence and droplet size lead to different oil slick patterns at the surface and significantly impact oil concentration. Possible effects for bio-degradation are also discussed. Funding from the GoMRI RFP-II is gratefully acknowledged.
Atomistic Molecular Dynamics Simulations of Crude Oil/Brine Displacement in Calcite Mesopores.
Sedghi, Mohammad; Piri, Mohammad; Goual, Lamia
2016-04-12
Unconventional reservoirs such as hydrocarbon-bearing shale formations and ultratight carbonates generate a large fraction of oil and gas production in North America. The characteristic feature of these reservoirs is their nanoscale porosity that provides significant surface areas between the pore walls and the occupying fluids. To better assess hydrocarbon recovery from these formations, it is crucial to develop an improved insight into the effects of wall-fluid interactions on the interfacial phenomena in these nanoscale confinements. One of the important properties that controls the displacement of fluids inside the pores is the threshold capillary pressure. In this study, we present the results of an integrated series of large-scale molecular dynamics (MD) simulations performed to investigate the effects of wall-fluid interactions on the threshold capillary pressures of oil-water/brine displacements in a calcite nanopore with a square cross section. Fully atomistic models are utilized to represent crude oil, brine, and calcite in order to accommodate electrostatic interactions and H-bonding between the polar molecules and the calcite surface. To this end, we create mixtures of various polar and nonpolar organic molecules to better represent the crude oil. The interfacial tension between oil and water/brine and their contact angle on calcite surface are simulated. We study the effects of oil composition, water salinity, and temperature and pressure conditions on these properties. The threshold capillary pressure values are also obtained from the MD simulations for the calcite nanopore. We then compare the MD results against those generated using the Mayer-Stowe-Princen (MSP) method and explain the differences. PMID:27010399
Sharapov, Vladimir A; Mandelshtam, Vladimir A
2007-10-18
We consider systems undergoing very-low-temperature solid-solid transitions associated with minima of similar energy but different symmetry, and separated by a high potential barrier. In such cases the well-known "broken-ergodicity" problem is often difficult to overcome, even using the most advanced Monte Carlo (MC) techniques, including the replica exchange method (REM). The methodology that we develop in this paper is suitable for the above specified cases and is numerically accurate and efficient. It is based on a new MC move implemented within the REM framework, in which trial points are generated analytically using an auxiliary harmonic superposition system that mimics well the true system at low temperatures. Due to the new move, the low-temperature random walks are able to frequently switch the relevant potential energy funnels leading to an efficient sampling. Numerically accurate results are obtained for a number of Lennard-Jones clusters, including those that have so far been treated only by the harmonic superposition approximation (HSA). The latter is believed to provide good estimates for low-temperature equilibrium properties but is manifestly uncontrollable and is difficult to validate. The present results provide a good test for the HSA and demonstrate its reliability, particularly for estimation of the solid-solid transition temperatures in most cases considered. PMID:17685597
3D simulation of floral oil storage in the scopa of South American insects
NASA Astrophysics Data System (ADS)
Ruettgers, Alexander; Griebel, Michael; Pastrik, Lars; Schmied, Heiko; Wittmann, Dieter; Scherrieble, Andreas; Dinkelmann, Albrecht; Stegmaier, Thomas; InstituteNumerical Simulation Team; Institute of Crop Science; Resource Conservation Team; Institute of Textile Technology; Process Engineering Team
2014-11-01
Several species of bees in South America possess structures to store and transport floral oils. By using closely spaced hairs at their back legs, the so called scopa, these bees can absorb and release oil droplets without loss. The high efficiency of this process is a matter of ongoing research. Basing on recent x-ray microtomography scans from the scopa of these bees at the Institute of Textile Technology and Process Engineering Denkendorf, we build a three-dimensional computer model. Using NaSt3DGPF, a two-phase flow solver developed at the Institute for Numerical Simulation of the University of Bonn, we perform massively parallel flow simulations with the complex micro-CT data. In this talk, we discuss the results of our simulations and the transfer of the x-ray measurement into a computer model. This research was funded under GR 1144/18-1 by the Deutsche Forschungsgemeinschaft (DFG).
A New High-Speed Oil-Free Turbine Engine Rotordynamic Simulator Test Rig
NASA Technical Reports Server (NTRS)
Howard, Samuel A.
2007-01-01
A new test rig has been developed for simulating high-speed turbomachinery rotor systems using Oil-Free foil air bearing technology. Foil air bearings have been used in turbomachinery, primarily air cycle machines, for the past four decades to eliminate the need for oil lubrication. The goal of applying this bearing technology to other classes of turbomachinery has prompted the fabrication of this test rig. The facility gives bearing designers the capability to test potential bearing designs with shafts that simulate the rotating components of a target machine without the high cost of building "make-and-break" hardware. The data collected from this rig can be used to make design changes to the shaft and bearings in subsequent design iterations. This paper describes the new test rig and demonstrates its capabilities through the initial run with a simulated shaft system.
NASA Astrophysics Data System (ADS)
Ou, Yihong; Du, Yang; Jiang, Xingsheng; Wang, Dong; Liang, Jianjun
2010-04-01
The study on the special phenomenon, occurrence process and control mechanism of gasoline-air mixture thermal ignition in underground oil depots is of important academic and applied value for enriching scientific theories of explosion safety, developing protective technology against fire and decreasing the number of fire accidents. In this paper, the research on thermal ignition process of gasoline-air mixture in model underground oil depots tunnel has been carried out by using experiment and numerical simulation methods. The calculation result has been demonstrated by the experiment data. The five stages of thermal ignition course, which are slow oxidation stage, rapid oxidation stage, fire stage, flameout stage and quench stage, have been firstly defined and accurately descried. According to the magnitude order of concentration, the species have been divided into six categories, which lay the foundation for explosion-proof design based on the role of different species. The influence of space scale on thermal ignition in small-scale space has been found, and the mechanism for not easy to fire is that the wall reflection causes the reflux of fluids and changes the distribution of heat and mass, so that the progress of chemical reactions in the whole space are also changed. The novel mathematical model on the basis of unification chemical kinetics and thermodynamics established in this paper provides supplementary means for the analysis of process and mechanism of thermal ignition.
NASA Technical Reports Server (NTRS)
Buggele, A. E.
1973-01-01
The complex problem why large space simulation chambers do not realize the true ultimate vacuum was investigated. Some contaminating factors affecting diffusion pump performance have been identified and some advances in vacuum/distillation/fractionation technology have been achieved which resulted in a two decade or more lower ultimate pressure. Data are presented to show the overall or individual contaminating effect of commonly used phthalate ester plasticizers of 390 to 530 molecular weight on diffusion pump performance. Methods for removing contaminants from diffusion pump silicone oil during operation and reclaiming contaminated oil by high vacuum molecular distillation are described.
Predictive data-derived Bayesian statistic-transport model and simulator of sunken oil mass
NASA Astrophysics Data System (ADS)
Echavarria Gregory, Maria Angelica
Sunken oil is difficult to locate because remote sensing techniques cannot as yet provide views of sunken oil over large areas. Moreover, the oil may re-suspend and sink with changes in salinity, sediment load, and temperature, making deterministic fate models difficult to deploy and calibrate when even the presence of sunken oil is difficult to assess. For these reasons, together with the expense of field data collection, there is a need for a statistical technique integrating limited data collection with stochastic transport modeling. Predictive Bayesian modeling techniques have been developed and demonstrated for exploiting limited information for decision support in many other applications. These techniques brought to a multi-modal Lagrangian modeling framework, representing a near-real time approach to locating and tracking sunken oil driven by intrinsic physical properties of field data collected following a spill after oil has begun collecting on a relatively flat bay bottom. Methods include (1) development of the conceptual predictive Bayesian model and multi-modal Gaussian computational approach based on theory and literature review; (2) development of an object-oriented programming and combinatorial structure capable of managing data, integration and computation over an uncertain and highly dimensional parameter space; (3) creating a new bi-dimensional approach of the method of images to account for curved shoreline boundaries; (4) confirmation of model capability for locating sunken oil patches using available (partial) real field data and capability for temporal projections near curved boundaries using simulated field data; and (5) development of a stand-alone open-source computer application with graphical user interface capable of calibrating instantaneous oil spill scenarios, obtaining sets maps of relative probability profiles at different prediction times and user-selected geographic areas and resolution, and capable of performing post
Crowley, Jason M; Tahir-Kheli, Jamil; Goddard, William A
2015-10-01
It has been established experimentally that Bi2Te3 and Bi2Se3 are topological insulators, with zero band gap surface states exhibiting linear dispersion at the Fermi energy. Standard density functional theory (DFT) methods such as PBE lead to large errors in the band gaps for such strongly correlated systems, while more accurate GW methods are too expensive computationally to apply to the thin films studied experimentally. We show here that the hybrid B3PW91 density functional yields GW-quality results for these systems at a computational cost comparable to PBE. The efficiency of our approach stems from the use of Gaussian basis functions instead of plane waves or augmented plane waves. This remarkable success without empirical corrections of any kind opens the door to computational studies of real chemistry involving the topological surface state, and our approach is expected to be applicable to other semiconductors with strong spin-orbit coupling. PMID:26722872
NASA Astrophysics Data System (ADS)
Pau, George Shu Heng; Shen, Chaopeng; Riley, William J.; Liu, Yaning
2016-02-01
The topography, and the biotic and abiotic parameters are typically upscaled to make watershed-scale hydrologic-biogeochemical models computationally tractable. However, upscaling procedure can produce biases when nonlinear interactions between different processes are not fully captured at coarse resolutions. Here we applied the Proper Orthogonal Decomposition Mapping Method (PODMM) to downscale the field solutions from a coarse (7 km) resolution grid to a fine (220 m) resolution grid. PODMM trains a reduced-order model (ROM) with coarse-resolution and fine-resolution solutions, here obtained using PAWS+CLM, a quasi-3-D watershed processes model that has been validated for many temperate watersheds. Subsequent fine-resolution solutions were approximated based only on coarse-resolution solutions and the ROM. The approximation errors were efficiently quantified using an error estimator. By jointly estimating correlated variables and temporally varying the ROM parameters, we further reduced the approximation errors by up to 20%. We also improved the method's robustness by constructing multiple ROMs using different set of variables, and selecting the best approximation based on the error estimator. The ROMs produced accurate downscaling of soil moisture, latent heat flux, and net primary production with O(1000) reduction in computational cost. The subgrid distributions were also nearly indistinguishable from the ones obtained using the fine-resolution model. Compared to coarse-resolution solutions, biases in upscaled ROM solutions were reduced by up to 80%. This method has the potential to help address the long-standing spatial scaling problem in hydrology and enable long-time integration, parameter estimation, and stochastic uncertainty analysis while accurately representing the heterogeneities.
NASA Astrophysics Data System (ADS)
Adidharma, Hertanto; Tan, Sugata P.
2016-07-01
Canonical Monte Carlo simulations on face-centered cubic (FCC) and hexagonal closed packed (HCP) Lennard-Jones (LJ) solids are conducted at very low temperatures (0.10 ≤ T∗ ≤ 1.20) and high densities (0.96 ≤ ρ∗ ≤ 1.30). A simple and robust method is introduced to determine whether or not the cutoff distance used in the simulation is large enough to provide accurate thermodynamic properties, which enables us to distinguish the properties of FCC from that of HCP LJ solids with confidence, despite their close similarities. Free-energy expressions derived from the simulation results are also proposed, not only to describe the properties of those individual structures but also the FCC-liquid, FCC-vapor, and FCC-HCP solid phase equilibria.
Lobato, I; Van Dyck, D
2015-09-01
The main features and the GPU implementation of the MULTEM program are presented and described. This new program performs accurate and fast multislice simulations by including higher order expansion of the multislice solution of the high energy Schrödinger equation, the correct subslicing of the three-dimensional potential and top-bottom surfaces. The program implements different kinds of simulation for CTEM, STEM, ED, PED, CBED, ADF-TEM and ABF-HC with proper treatment of the spatial and temporal incoherences. The multislice approach described here treats the specimen as amorphous material which allows a straightforward implementation of the frozen phonon approximation. The generalized transmission function for each slice is calculated when is needed and then discarded. This allows us to perform large simulations that can include millions of atoms and keep the computer memory requirements to a reasonable level. PMID:25965576
Ovchinnikov, Victor; Nam, Kwangho; Karplus, Martin
2016-08-25
A method is developed to obtain simultaneously free energy profiles and diffusion constants from restrained molecular simulations in diffusive systems. The method is based on low-order expansions of the free energy and diffusivity as functions of the reaction coordinate. These expansions lead to simple analytical relationships between simulation statistics and model parameters. The method is tested on 1D and 2D model systems; its accuracy is found to be comparable to or better than that of the existing alternatives, which are briefly discussed. An important aspect of the method is that the free energy is constructed by integrating its derivatives, which can be computed without need for overlapping sampling windows. The implementation of the method in any molecular simulation program that supports external umbrella potentials (e.g., CHARMM) requires modification of only a few lines of code. As a demonstration of its applicability to realistic biomolecular systems, the method is applied to model the α-helix ↔ β-sheet transition in a 16-residue peptide in implicit solvent, with the reaction coordinate provided by the string method. Possible modifications of the method are briefly discussed; they include generalization to multidimensional reaction coordinates [in the spirit of the model of Ermak and McCammon (Ermak, D. L.; McCammon, J. A. J. Chem. Phys. 1978, 69, 1352-1360)], a higher-order expansion of the free energy surface, applicability in nonequilibrium systems, and a simple test for Markovianity. In view of the small overhead of the method relative to standard umbrella sampling, we suggest its routine application in the cases where umbrella potential simulations are appropriate. PMID:27135391
Castaneda, Jaime N.; Shollenberger, Kim Ann; Torczynski, John Robert; Cote, Raymond O.; Barney, Jeremy; O'Hern, Timothy John
2003-10-01
An experimental program is being conducted to study a proposed approach for oil reintroduction in the Strategic Petroleum Reserve (SPR). The goal is to assess whether useful oil is rendered unusable through formation of a stable oil-brine emulsion during reintroduction of degassed oil into the brine layer in storage caverns. This report documents the first stage of the program, in which simulant liquids are used to characterize the buoyant plume that is produced when a jet of crude oil is injected downward from a tube into brine. The experiment consists of a large transparent vessel that is a scale model of the proposed oil injection process at the SPR. An oil layer is floated on top of a brine layer. Silicon oil (Dow Corning 200{reg_sign} Fluid, 5 cSt) is used as the simulant for crude oil to allow visualization of the flow and to avoid flammability and related concerns. Sodium nitrate solution is used as the simulant for brine because it is not corrosive and it can match the density ratio between brine and crude oil. The oil is injected downward through a tube into the brine at a prescribed depth below the oil-brine interface. Flow rates are determined by scaling to match the ratio of buoyancy to momentum between the experiment and the SPR. Initially, the momentum of the flow produces a downward jet of oil below the tube end. Subsequently, the oil breaks up into droplets due to shear forces, buoyancy dominates the flow, and a plume of oil droplets rises to the interface. The interface is deflected upward by the impinging oil-brine plume. Two different diameter injection tubes were used (1/2-inch and 1-inch OD) to vary the scaling. Use of the 1-inch injection tube also assured that turbulent pipe flow was achieved, which was questionable for lower flow rates in the 1/2-inch tube. In addition, a 1/2-inch J-tube was used to direct the buoyant jet upwards rather than downwards to determine whether flow redirection could substantially reduce the oil-plume size and the
Assaraf, Roland; Caffarel, Michel; Kollias, A C
2011-04-15
We present a method to efficiently evaluate small energy differences of two close N-body systems by employing stochastic processes having a stability versus chaos property. By using the same random noise, energy differences are computed from close trajectories without reweighting procedures. The approach is presented for quantum systems but can be applied to classical N-body systems as well. It is exemplified with diffusion Monte Carlo simulations for long chains of hydrogen atoms and molecules for which it is shown that the long-standing problem of computing energy derivatives is solved. PMID:21568537
NASA Astrophysics Data System (ADS)
Kairn, T.; Crowe, S. B.; Charles, P. H.; Trapp, J. V.
2014-03-01
This study investigates the variation of photon field penumbra shape with initial electron beam diameter, for very narrow beams. A Varian Millenium MLC (Varian Medical Systems, Palo Alto, USA) and a Brainlab m3 microMLC (Brainlab AB. Feldkirchen, Germany) were used, with one Varian iX linear accelerator, to produce fields that were (nominally) 0.20 cm across. Dose profiles for these fields were measured using radiochromic film and compared with the results of simulations completed using BEAMnrc and DOSXYZnrc, where the initial electron beam was set to FWHM = 0.02, 0.10, 0.12, 0.15, 0.20 and 0.50 cm. Increasing the electron-beam FWHM produced increasing occlusion of the photon source by the closely spaced collimator leaves and resulted in blurring of the simulated profile widths from 0.24 to 0.58 cm, for the MLC, from 0.11 to 0.40 cm, for the microMLC. Comparison with measurement data suggested that the electron spot size in the clinical linear accelerator was between FWHM = 0.10 and 0.15 cm, encompassing the result of our previous output-factor based work, which identified a FWHM of 0.12 cm. Investigation of narrow-beam penumbra variation has been found to be a useful procedure, with results varying noticeably with linear accelerator spot size and allowing FWHM estimates obtained using other methods to be verified.
Molecular Dynamics Simulation of β-Lactoglobulin at Different Oil/Water Interfaces.
Zare, Davoud; Allison, Jane R; McGrath, Kathryn M
2016-05-01
Controlling and manipulating protein behavior at an interface is of immense relevance to a broad range of physicochemical and biological phenomena and technological processes. Although many experimental studies have contributed to rapid progress in the fundamental knowledge of protein behavior at interfaces, detailed molecular-level understanding of the mechanism of protein adsorption at an interface is still remarkably lacking. In this study, atomistic molecular dynamics simulations were used to characterize the adsorption of β-lactoglobulin at two different oil/water (O/W) interfaces, where the oil was either the marginally hydrophilic octanol or the more hydrophilic triolein, and the results were compared to those of a previous study utilizing the hydrophobic oil decane. Both the approach to the surface and the mechanism of adsorption depend upon the hydrophilicity of the oil and the interfacial tension of the O/W interface, with the nature of the adsorption, the accompanying structural changes, and the energetic driving force differing markedly between the different oils. Intriguingly, the behavior of the protein resembles that predicted for a soft spherical particle at an O/W interface. The results are also in agreement with key experimental findings, particularly the observation that proteins undergo more conformational change upon adsorption to hydrophobic surfaces, flattening out to expose hydrophobic interior residues to the surface, and that a thicker layer of native-like adsorbed protein forms at hydrophilic surfaces, and reveal structural and mechanistic detail behind each mechanism of adsorption. PMID:27075297
Simulation of EOR (enhanced oil recovery) processes in stochastically generated permeable media
Waggoner, J.R.; Castillo, J.L.; Lake, L.W. . Dept. of Petroleum Engineering)
1990-01-01
Many enhanced oil recovery (EOR) processes involve injecting an agent, such as steam or CO{sub 2}, that is much more mobile than the resident oil. Other EOR processes attempt to improve sweep efficiency by adding polymer or surfactant to the injected water to create a favorable mobility ratio. This study examines the effect of statistically generated heterogeneity on miscible displacements at unfavorable and favorable mobility ratios. The principal goal is to delineate the effects of fingering, dispersion and channeling on volumetric sweep efficiency. Two-dimensional heterogeneous permeability fields are generated with variability (heterogeneity) and spatial correlation as characterizing parameters. Four levels of correlation and three of variability make up a 12 element matrix. At each element of the matrix, a miscible displacement simulation at unit mobility ratio shows the effect of the heterogeneity, and simulations at mobility ratios of 10 and 0.5 show the effect of viscous force differences combined with heterogeneity. 20 refs., 7 figs., 3 tabs.
Mesoscopic simulation study on phase diagram of the system oil/water/aerosol OT
NASA Astrophysics Data System (ADS)
Yuan, Shi-Ling; Cai, Zheng-Ting; Xu, Gui-Ying; Jiang, Yuan-Sheng
2002-10-01
A simple model, i.e., sodium bis(2-ethylhexyl) sulfosuccinate (AOT) represented by one-head and two-tail beads tied together by a harmonic spring and water or iso-octane by one bead, was put forward via dissipative particles dynamics (DPD) simulation method. Using the changes of interfacial tension between water and oil phase, a ternary phase diagram of AOT/water/iso-octane system was drawn.
A dual-porosity model for simulating solute transport in oil shale
Glover, K.C.
1987-01-01
A model is described for simulating three-dimensional groundwater flow and solute transport in oil shale and associated geohydrologic units. The model treats oil shale as a dual-porosity medium by simulating flow and transport within fractures using the finite-element method. Diffusion of solute between fractures and the essentially static water of the shale matrix is simulated by including an analytical solution that acts as a source-sink term to the differential equation of solute transport. While knowledge of fracture orientation and spacing is needed to effectively use the model, it is not necessary to map the locations of individual fractures. The computer program listed in the report incorporates many of the features of previous dual-porosity models while retaining a practical approach to solving field problems. As a result the theory of solute transport is not extended in any appreciable way. The emphasis is on bringing together various aspects of solute transport theory in a manner that is particularly suited to the unusual groundwater flow and solute transport characteristics of oil shale systems. (Author 's abstract)
NASA Astrophysics Data System (ADS)
Baiardi, Alberto; Barone, Vincenzo; Biczysko, Malgorzata; Bloino, Julien
2014-06-01
Two parallel theories including Franck-Condon, Herzberg-Teller and Duschinsky (i.e., mode mixing) effects, allowing different approximations for the description of excited state PES have been developed in order to simulate realistic, asymmetric, electronic spectra line-shapes taking into account the vibrational structure: the so-called sum-over-states or time-independent (TI) method and the alternative time-dependent (TD) approach, which exploits the properties of the Fourier transform. The integrated TI-TD procedure included within a general purpose QM code [1,2], allows to compute one photon absorption, fluorescence, phosphorescence, electronic circular dichroism, circularly polarized luminescence and resonance Raman spectra. Combining both approaches, which use a single set of starting data, permits to profit from their respective advantages and minimize their respective limits: the time-dependent route automatically includes all vibrational states and, possibly, temperature effects, while the time-independent route allows to identify and assign single vibronic transitions. Interpretation, analysis and assignment of experimental spectra based on integrated TI-TD vibronic computations will be illustrated for challenging cases of medium-sized open-shell systems in the gas and condensed phases with inclusion of leading anharmonic effects. 1. V. Barone, A. Baiardi, M. Biczysko, J. Bloino, C. Cappelli, F. Lipparini Phys. Chem. Chem. Phys, 14, 12404, (2012) 2. A. Baiardi, V. Barone, J. Bloino J. Chem. Theory Comput., 9, 4097-4115 (2013)
NASA Astrophysics Data System (ADS)
Lemoine, X.; Sriram, S.; Kergen, R.
2011-05-01
ArcelorMittal continuously develops new steel grades (AHSS) with high performance for the automotive industry to improve the weight reduction and the passive safety. The wide market introduction of AHSS raises a new challenge for manufacturers in terms of material models in the prediction of forming—especially formability and springback. The relatively low uniform elongation, the high UTS and the low forming limit curve of these AHSS may cause difficulties in forming simulations. One of these difficulties is the consequence of the relatively low uniform elongation on the parameters identification of isotropic hardening model. Different experimental tests allow to reach large plastic strain levels (hydraulic bulge test, stack compression test, shear test…). After a description on how to determine the flow curve in these experimental tests, a comparison of the different flow curves is made for different steel grades. The ArcelorMittal identification protocol for hardening models is only based on stress-strain curves determined in uniaxial tension. Experimental tests where large plastic strain levels are reached are used to validate our identification protocol and to recommend some hardening models. Finally, the influence of isotropic hardening models and yield loci in forming prediction for AHSS steels will be presented.
NASA Astrophysics Data System (ADS)
Yu, D. O.; Kwon, O. J.
2014-06-01
In the present study, aeroelastic simulations of horizontal-axis wind turbine rotor blades were conducted using a coupled CFD-CSD method. The unsteady blade aerodynamic loads and the dynamic blade response due to yaw misalignment and non-uniform sheared wind were investigated. For this purpose, a CFD code solving the RANS equations on unstructured meshes and a FEM-based CSD beam solver were used. The coupling of the CFD and CSD solvers was made by exchanging the data between the two solvers in a loosely coupled manner. The present coupled CFD-CSD method was applied to the NREL 5MW reference wind turbine rotor, and the results were compared with those of CFD-alone rigid blade calculations. It was found that aeroelastic blade deformation leads to a significant reduction of blade aerodynamic loads, and alters the unsteady load behaviours, mainly due to the torsional deformation. The reduction of blade aerodynamic loads is particularly significant at the advancing rotor blade side for yawed flow conditions, and at the upper half of rotor disk where wind velocity is higher due to wind shear.
NASA Astrophysics Data System (ADS)
Geng, Xiaolong; Pan, Zhong; Boufadel, Michel C.; Ozgokmen, Tamay; Lee, Kenneth; Zhao, Lin
2016-04-01
Numerical experiments of oil bioremediation of tidally influenced beach were simulated using the model BIOMARUN. Nutrient and dissolved oxygen were assumed present in a solution applied on the exposed beach face, and the concentration of these amendments was tracked throughout the beach for up to 6 months. It was found that, in comparison to natural attenuation, bioremediation increased the removal efficiency by 76% and 65% for alkanes and aromatics, respectively. Increasing the nutrient concentration in the applied solution did not always enhance biodegradation as oxygen became limiting even when the beach was originally oxygen-rich. Therefore, replenishment of oxygen to oil-contaminated zone was also essential. Stimulation of oil biodegradation was more evident in the upper and midintertidal zone of the beach, and less in the lower intertidal zone. This was due to reduced nutrient and oxygen replenishment, as very little of the amendment solution reached that zone. It was found that under continual application, most of the oil biodegraded within 2 months, while it persisted for 6 months under natural conditions. While the difference in duration suggests minimal long-term effects, there are situations where the beach would need to be cleaned for major ecological functions, such as temporary nesting or feeding for migratory birds. Biochemical retention time map (BRTM) showed that the duration of solution application was dependent upon the stimulated oil biodegradation rate. By contrast, the application rate of the amendment solution was dependent upon the subsurface extent of the oil-contaminated zone. Delivery of nutrient and oxygen into coastal beach involved complex interaction among amendment solution, groundwater, and seawater. Therefore, approaches that ignore the hydrodynamics due to tide are unlikely to provide the optimal solutions for shoreline bioremediation.
NASA Astrophysics Data System (ADS)
Reinhardt, Colin N.; Ritcey, James A.
2015-09-01
We present a novel method for efficient and physically-accurate modeling & simulation of anisoplanatic imaging through the atmosphere; in particular we present a new space-variant volumetric image blur algorithm. The method is based on the use of physical atmospheric meteorology models, such as vertical turbulence profiles and aerosol/molecular profiles which can be in general fully spatially-varying in 3 dimensions and also evolving in time. The space-variant modeling method relies on the metadata provided by 3D computer graphics modeling and rendering systems to decompose the image into a set of slices which can be treated in an independent but physically consistent manner to achieve simulated image blur effects which are more accurate and realistic than the homogeneous and stationary blurring methods which are commonly used today. We also present a simple illustrative example of the application of our algorithm, and show its results and performance are in agreement with the expected relative trends and behavior of the prescribed turbulence profile physical model used to define the initial spatially-varying environmental scenario conditions. We present the details of an efficient Fourier-transform-domain formulation of the SV volumetric blur algorithm and detailed algorithm pseudocode description of the method implementation and clarification of some nonobvious technical details.
NASA Astrophysics Data System (ADS)
Chang, Y.; Huang, W.
2007-12-01
Evolution of fluorescence color of inclusion oils has been simulated by measuring in-situ the fluorescence of `live' oils generated from thirteen oil-prone kerogens from different depositional environments during a closed system pyrolysis in a diamond anvil cell at three heating rates (3, 8, and 25°C /min) up to 600°C. The measured fluorescence intensity of samples increases considerably within maturation intervals close to oil windows, while the fluorescence spectra of oils generated from all studied kerogens exhibit exclusively a progressive blue-shift of peak wavelengths (λmax) and red/green quotient (Q) upon increasing maturity. The observation is consistent with the maturity dependence of spectral shift trend widely recognized in natural hydrocarbon inclusions or crude oils. This study furthermore reveals that the acclaimed direction of spectral shift for inclusion oils is mostly independent of sources of their parental kerogens, implying that some reverse or anomalous trends reported in inclusion oils may be attributed to other processes subsequent to their generation, which significantly altered the fluorescence properties of oils. However, the experimental maturity corresponding to each color (λmax or Q) of oils can vary significantly (± 0.2 %Ro) among their sourced kerogens, suggesting that single fluorescence color of crude or inclusion oil is both maturity- and source-dependent and therefore may not be a good indication of its maturity. In addition, the blue-shift of cumulative oils generated from all kerogens approaches similar minima λmax around 564 nm or Q- value around 0.6 at maturity close to the middle or late stage of oil generation, implying that most late cumulative oils may exhibit similar colors. The oils generated in a maturity interval in late stage, however, can exhibit color of shorter wavelength less than the minimum.
NASA Astrophysics Data System (ADS)
Byun, Jaeseung; Bodony, Daniel; Pantano, Carlos
2014-11-01
Improved order-of-accuracy discretizations often require careful consideration of their numerical stability. We report on new high-order finite difference schemes using Summation-By-Parts (SBP) operators along with the Simultaneous-Approximation-Terms (SAT) boundary condition treatment for first and second-order spatial derivatives with variable coefficients. In particular, we present a highly accurate operator for SBP-SAT-based approximations of second-order derivatives with variable coefficients for Dirichlet and Neumann boundary conditions. These terms are responsible for approximating the physical dissipation of kinetic and thermal energy in a simulation, and contain grid metrics when the grid is curvilinear. Analysis using the Laplace transform method shows that strong stability is ensured with Dirichlet boundary conditions while weaker stability is obtained for Neumann boundary conditions. Furthermore, the benefits of the scheme is shown in the direct numerical simulation (DNS) of a Mach 1.5 compressible turbulent supersonic jet using curvilinear grids and skew-symmetric discretization. Particularly, we show that the improved methods allow minimization of the numerical filter often employed in these simulations and we discuss the qualities of the simulation.
Reeves, T.L.; Turner, J.P.; Hasfurther, V.R.; Skinner, Q.D.
1992-06-01
The scope of this program is to study interacting hydrologic, geotechnical, and chemical factors affecting the behavior and disposal of combusted processed oil shale. The research combines bench-scale testing with large scale research sufficient to describe commercial scale embankment behavior. The large scale approach was accomplished by establishing five lysimeters, each 7.3 {times} 3.0 {times} 3.0 m deep, filled with processed oil shale that has been retorted and combusted by the Lurgi-Ruhrgas (Lurgi) process. Approximately 400 tons of Lurgi processed oil shale waste was provided by RBOSC to carry out this study. Research objectives were designed to evaluate hydrologic, geotechnical, and chemical properties and conditions which would affect the design and performance of large-scale embankments. The objectives of this research are: assess the unsaturated movement and redistribution of water and the development of potential saturated zones and drainage in disposed processed oil shale under natural and simulated climatic conditions; assess the unsaturated movement of solubles and major chemical constituents in disposed processed oil shale under natural and simulated climatic conditions; assess the physical and constitutive properties of the processed oil shale and determine potential changes in these properties caused by disposal and weathering by natural and simulated climatic conditions; assess the use of previously developed computer model(s) to describe the infiltration, unsaturated movement, redistribution, and drainage of water in disposed processed oil shale; evaluate the stability of field scale processed oil shale solid waste embankments using computer models.
Reeves, T.L.; Turner, J.P.; Hasfurther, V.R.; Skinner, Q.D.
1992-06-01
The scope of this program is to study interacting hydrologic, geotechnical, and chemical factors affecting the behavior and disposal of combusted processed oil shale. The research combines bench-scale testing with large scale research sufficient to describe commercial scale embankment behavior. The large scale approach was accomplished by establishing five lysimeters, each 7.3 [times] 3.0 [times] 3.0 m deep, filled with processed oil shale that has been retorted and combusted by the Lurgi-Ruhrgas (Lurgi) process. Approximately 400 tons of Lurgi processed oil shale waste was provided by RBOSC to carry out this study. Research objectives were designed to evaluate hydrologic, geotechnical, and chemical properties and conditions which would affect the design and performance of large-scale embankments. The objectives of this research are: assess the unsaturated movement and redistribution of water and the development of potential saturated zones and drainage in disposed processed oil shale under natural and simulated climatic conditions; assess the unsaturated movement of solubles and major chemical constituents in disposed processed oil shale under natural and simulated climatic conditions; assess the physical and constitutive properties of the processed oil shale and determine potential changes in these properties caused by disposal and weathering by natural and simulated climatic conditions; assess the use of previously developed computer model(s) to describe the infiltration, unsaturated movement, redistribution, and drainage of water in disposed processed oil shale; evaluate the stability of field scale processed oil shale solid waste embankments using computer models.
Simulation of the landfall of the Deepwater Horizon oil on the shorelines of the Gulf of Mexico.
Boufadel, Michel C; Abdollahi-Nasab, Ali; Geng, Xiaolong; Galt, Jerry; Torlapati, Jagadish
2014-08-19
We conducted simulations of oil transport from the footprint of the Macondo Well on the water surface throughout the Gulf of Mexico, including deposition on the shorelines. We used the U.S. National Oceanic Atmospheric Administration (NOAA) model General NOAA Operational Modeling Environment (GNOME) and the same parameter values and input adopted by NOAA following the Deepwater Horizon (DWH) blowout. We found that the disappearance rate of oil off the water surface was most likely around 20% per day based on satellite-based observations of the disappearance rate of oil detected on the sea surface after the DWH wellhead was capped. The simulations and oil mass estimates suggest that the mass of oil that reached the shorelines was between 10,000 and 30,000 tons, with an expected value of 22,000 tons. More than 90% of the oil deposition occurred on the Louisiana shorelines, and it occurred in two batches. Simulations revealed that capping the well after 2 weeks would have resulted in only 30% of the total oil depositing on the shorelines, while capping after 3 weeks would have resulted in 60% deposition. Additional delay in capping after 3 weeks would have averted little additional shoreline oiling over the ensuing 4 weeks. PMID:25068902
Kirtland, D.A.
1986-01-01
The stimulus for this research is the belief that the International Energy Agency (IEA) would be better able to plant for and implement its Emergency Oil Sharing Plan (EOSP) if it modeled important aspects of the world oil economy from an operational point of view. A tanker simulation (TANKSIM) model is offered as a prototype to produce performance measures that may be used in investigating the adequacy and responsiveness of the world's ports and tankers to various patterns of trade. TANKSIM is a SLAM-based, discrete-event simulation model that tracks the movements of tankers worldwide and compiles information on the fleet disposition and status and port and vessel operating characteristics associated with an input trade pattern. Any pattern of crude oil trade whether historical, prescribed, or predicted, may be used to drive the model. A flexible model structure enables a user to investigate the behavior of any configuration of nations, ports, and tanker types. Model use is demonstrated with a pair of hypothetical runs, a Base Case and Scenario. The Base Case considers a portion of the world's crude-hauling tankers operating among IEA nations and members of the Organization of Petroleum Exporting Countries (OPEC) for one year. The Scenario alters the Base Case trade pattern by representing a hypothetical cessation of Iranian exports due to the Iran/Iraq war and an IEA-EOSP sponsored increase in Saudi Arabian exports in response.
Chakravarty, A.; Emanuel, A.S.; Bernath, J.A.
1997-08-01
The application of streamtube models for reservoir simulation has an extensive history in the oil industry. Although these models are strictly applicable only to fields under voidage balance, they have proved to be useful in a large number of fields provided that there is no solution gas evolution and production. These models combine the benefit of very fast computational time with the practical ability to model a large reservoir over the course of its history. These models do not, however, directly incorporate the detailed geological information that recent experience has taught is important. This paper presents a technique for mapping the saturation information contained in a history matched streamtube model onto a detailed geostatistically derived finite difference grid. With this technique, the saturation information in a streamtube model, data that is actually statistical in nature, can be identified with actual physical locations in a field and a picture of the remaining oil saturation can be determined. Alternatively, the streamtube model can be used to simulate the early development history of a field and the saturation data then used to initialize detailed late time finite difference models. The proposed method is presented through an example application to the Ninian reservoir. This reservoir, located in the North Sea (UK), is a heterogeneous sandstone characterized by a line drive waterflood, with about 160 wells, and a 16 year history. The reservoir was satisfactorily history matched and mapped for remaining oil saturation. A comparison to 3-D seismic survey and recently drilled wells have provided preliminary verification.
Paris, Claire B; Hénaff, Matthieu Le; Aman, Zachary M; Subramaniam, Ajit; Helgers, Judith; Wang, Dong-Ping; Kourafalou, Vassiliki H; Srinivasan, Ashwanth
2012-12-18
During the Deepwater Horizon incident, crude oil flowed into the Gulf of Mexico from 1522 m underwater. In an effort to prevent the oil from rising to the surface, synthetic dispersants were applied at the wellhead. However, uncertainties in the formation of oil droplets and difficulties in measuring their size in the water column, complicated further assessment of the potential effect of the dispersant on the subsea-to-surface oil partition. We adapted a coupled hydrodynamic and stochastic buoyant particle-tracking model to the transport and fate of hydrocarbon fractions and simulated the far-field transport of the oil from the intrusion depth. The evaluated model represented a baseline for numerical experiments where we varied the distributions of particle sizes and thus oil mass. The experiments allowed to quantify the relative effects of chemical dispersion, vertical currents, and inertial buoyancy motion on oil rise velocities. We present a plausible model scenario, where some oil is trapped at depth through shear emulsification due to the particular conditions of the Macondo blowout. Assuming effective mixing of the synthetic dispersants at the wellhead, the model indicates that the submerged oil mass is shifted deeper, decreasing only marginally the amount of oil surfacing. In this scenario, the oil rises slowly to the surface or stays immersed. This suggests that other mechanisms may have contributed to the rapid surfacing of oil-gas mixture observed initially. The study also reveals local topographic and hydrodynamic processes that influence the oil transport in eddies and multiple layers. This numerical approach provides novel insights on oil transport mechanisms from deep blowouts and on gauging the subsea use of synthetic dispersant in mitigating coastal damage. PMID:23167517
Global simulations of smoke from Kuwaiti oil fires and possible effects on climate
Glatzmaier, G.A.; Malone, R.C.; Kao, C.Y.J.
1991-12-31
The Los Alamos Global Climate Model has bee used to simulate the global evolution of the Kuwaiti oil fire smoke and its potential effects on the climate. The initial simulations were done shortly before the fires were lit in January 1991. They indicated that such an event would not result in a ``Mini Nuclear Winter`` as some people were suggesting. Further simulations during the year suggested that the smoke could be responsible for subtle regional climate changes in the spring such as a 5 degree centigrade decrease in the surface temperature in Kuwait, a 10% decrease in precipitation in Saudi Arabia and a 10% increase in precipitation in the Tibetan Plateau region. These results are in qualitative agreement with the observations this year.
A fortran program for Monte Carlo simulation of oil-field discovery sequences
Bohling, G.C.; Davis, J.C.
1993-01-01
We have developed a program for performing Monte Carlo simulation of oil-field discovery histories. A synthetic parent population of fields is generated as a finite sample from a distribution of specified form. The discovery sequence then is simulated by sampling without replacement from this parent population in accordance with a probabilistic discovery process model. The program computes a chi-squared deviation between synthetic and actual discovery sequences as a function of the parameters of the discovery process model, the number of fields in the parent population, and the distributional parameters of the parent population. The program employs the three-parameter log gamma model for the distribution of field sizes and employs a two-parameter discovery process model, allowing the simulation of a wide range of scenarios. ?? 1993.
Hsieh, Paul
2010-01-01
This report describes the application of a computer model to simulate reservoir depletion and oil flow from the Macondo well following the Deepwater Horizon blowout. Reservoir and fluid data used for model development are based on (1) information released in BP's investigation report of the incident, (2) information provided by BP personnel during meetings in Houston, Texas, and (3) calibration by history matching to shut-in pressures measured in the capping stack during the Well Integrity Test. The model is able to closely match the measured shut-in pressures. In the simulation of the 86-day period from the blowout to shut in, the simulated reservoir pressure at the well face declines from the initial reservoir pressure of 11,850 pounds per square inch (psi) to 9,400 psi. After shut in, the simulated reservoir pressure recovers to a final value of 10,300 psi. The pressure does not recover back to the initial pressure owing to reservoir depletion caused by 86 days of oil discharge. The simulated oil flow rate declines from 63,600 stock tank barrels per day just after the Deepwater Horizon blowout to 52,600 stock tank barrels per day just prior to shut in. The simulated total volume of oil discharged is 4.92 million stock tank barrels. The overall uncertainty in the simulated flow rates and total volume of oil discharged is estimated to be + or - 10 percent.
Integration of seismic methods with reservoir simulation, Pikes Peak heavy oil field, Saskatchewan
NASA Astrophysics Data System (ADS)
Zou, Ying
The Pikes Peak heavy oil field has been operated by Husky Energy Ltd since 1981. Steam injection has been successfully employed to increase production. Efforts in geophysics and reservoir engineering have been made to improve interpretations in the mapping of reservoir conditions. This dissertation developed tools and a working flow for integrating the analysis of time-lapse seismic surveys with reservoir simulation, and applied them to the Pikes Peak field. Two time-lapse 2D seismic lines acquired in February 1991 and March 2000 in the eastern part of the field were carefully processed to produce wavelet and structure matched final sections. Reservoir simulation based on the field reservoir production history was carried out. It provided independent complementary information for the time-lapse seismic analysis. A rock physics procedure based on Gassmann's equation and Batzle and Wang's empirical relationship successfully linked the reservoir engineering to the seismic method. Based on the resultant seismic models, synthetic seismic sections were generated as the analogy of field seismic sections. The integrated interpretation for the Pikes Peak reservoir drew the following conclusions: The areas with a gas saturation difference, between two compared time steps, have seismic differences. Thicker gas zones correspond with large reflectivity changes on the top of the reservoir and larger traveltime delays in the seismic section. The thin gas zones only induce large reflectivity changes on the top of the reservoir, and do not have large time delays below the reservoir zone. High temperature regions also correlate with areas having large seismic energy differences. High temperature with thick gas (steam and methane) zones may be evidence for steam existence. The seismic differences at locations far from the production zone are due to the lower pressure that causes solution gas to evolve from the oil. Pressure changes propagate much faster (˜20 m in one month) than
The sequential method for the black-oil reservoir simulation on unstructured grids
NASA Astrophysics Data System (ADS)
Li, Baoyan; Chen, Zhangxin; Huan, Guanren
2003-11-01
This paper presents new results for applying the sequential solution method to the black-oil reservoir simulation with unstructured grids. The fully implicit solution method has been successfully applied to reservoir simulation with unstructured grids. However, the complexity of the fully implicit method and the irregularity of the grids result in a very complicated structure of linear equation systems (LESs) and in high computational cost to solve them. To tackle this problem, the sequential method is applied to reduce the size of the LESs. To deal with instable problems caused by the low implicit degree of this method, some practical techniques are introduced to control convergence of Newton-Raphson's iterations which are exploited in the linearization of the governing equations of the black-oil model. These techniques are tested with the benchmark problem of the ninth comparative solution project (CSP) organized by the society of petroleum engineers (SPE) and applied to field-scale models of both saturated and undersaturated reservoirs. The simulation results show that the sequential method uses as little as 20.01% of the memory for solving the LESs and 23.89% of the total computational time of the fully implicit method to reach the same precision for the undersaturated reservoirs, when the same iteration control parameters are used for both solution methods. However, for the saturated reservoirs, the sequential method must use stricter iteration control parameters to reach the same precision as the fully implicit method.
NASA Astrophysics Data System (ADS)
Iskandarani, Mohamed; Wang, Shitao; Srinivasan, Ashwanth; Carlisle Thacker, W.; Winokur, Justin; Knio, Omar M.
2016-04-01
We give an overview of four different ensemble-based techniques for uncertainty quantification and illustrate their application in the context of oil plume simulations. These techniques share the common paradigm of constructing a model proxy that efficiently captures the functional dependence of the model output on uncertain model inputs. This proxy is then used to explore the space of uncertain inputs using a large number of samples, so that reliable estimates of the model's output statistics can be calculated. Three of these techniques use polynomial chaos (PC) expansions to construct the model proxy, but they differ in their approach to determining the expansions' coefficients; the fourth technique uses Gaussian Process Regression (GPR). An integral plume model for simulating the Deepwater Horizon oil-gas blowout provides examples for illustrating the different techniques. A Monte Carlo ensemble of 50,000 model simulations is used for gauging the performance of the different proxies. The examples illustrate how regression-based techniques can outperform projection-based techniques when the model output is noisy. They also demonstrate that robust uncertainty analysis can be performed at a fraction of the cost of the Monte Carlo calculation.
NASA Astrophysics Data System (ADS)
Vaganova, N. A.; Filimonov, M. Yu.
2015-11-01
A mathematical model, numerical algorithm and program code for simulation and long-term forecasting of changes in permafrost as a result of operation of a multiple well pad of northern oil and gas field are presented. In the model the most significant climatic and physical factors are taken into account such as solar radiation, determined by specific geographical location, heterogeneous structure of frozen soil, thermal stabilization of soil, possible insulation of the objects, seasonal fluctuations in air temperature, and freezing and thawing of the upper soil layer. Results of computing are presented.
Vaganova, N. A.; Filimonov, M. Yu.
2015-11-30
A mathematical model, numerical algorithm and program code for simulation and long-term forecasting of changes in permafrost as a result of operation of a multiple well pad of northern oil and gas field are presented. In the model the most significant climatic and physical factors are taken into account such as solar radiation, determined by specific geographical location, heterogeneous structure of frozen soil, thermal stabilization of soil, possible insulation of the objects, seasonal fluctuations in air temperature, and freezing and thawing of the upper soil layer. Results of computing are presented.
NASA Astrophysics Data System (ADS)
Zhang, Na; Yao, Jun; Huang, Zhaoqin; Wang, Yueying
2013-06-01
Numerical simulation in naturally fractured media is challenging because of the coexistence of porous media and fractures on multiple scales that need to be coupled. We present a new approach to reservoir simulation that gives accurate resolution of both large-scale and fine-scale flow patterns. Multiscale methods are suitable for this type of modeling, because it enables capturing the large scale behavior of the solution without solving all the small features. Dual-porosity models in view of their strength and simplicity can be mainly used for sugar-cube representation of fractured media. In such a representation, the transfer function between the fracture and the matrix block can be readily calculated for water-wet media. For a mixed-wet system, the evaluation of the transfer function becomes complicated due to the effect of gravity. In this work, we use a multiscale finite element method (MsFEM) for two-phase flow in fractured media using the discrete-fracture model. By combining MsFEM with the discrete-fracture model, we aim towards a numerical scheme that facilitates fractured reservoir simulation without upscaling. MsFEM uses a standard Darcy model to approximate the pressure and saturation on a coarse grid, whereas fine scale effects are captured through basis functions constructed by solving local flow problems using the discrete-fracture model. The accuracy and the robustness of MsFEM are shown through several examples. In the first example, we consider several small fractures in a matrix and then compare the results solved by the finite element method. Then, we use the MsFEM in more complex models. The results indicate that the MsFEM is a promising path toward direct simulation of highly resolution geomodels.
NASA Astrophysics Data System (ADS)
Ahmed, Mahmoud; Eslamian, Morteza
2015-07-01
Laminar natural convection in differentially heated ( β = 0°, where β is the inclination angle), inclined ( β = 30° and 60°), and bottom-heated ( β = 90°) square enclosures filled with a nanofluid is investigated, using a two-phase lattice Boltzmann simulation approach. The effects of the inclination angle on Nu number and convection heat transfer coefficient are studied. The effects of thermophoresis and Brownian forces which create a relative drift or slip velocity between the particles and the base fluid are included in the simulation. The effect of thermophoresis is considered using an accurate and quantitative formula proposed by the authors. Some of the existing results on natural convection are erroneous due to using wrong thermophoresis models or simply ignoring the effect. Here we show that thermophoresis has a considerable effect on heat transfer augmentation in laminar natural convection. Our non-homogenous modeling approach shows that heat transfer in nanofluids is a function of the inclination angle and Ra number. It also reveals some details of flow behavior which cannot be captured by single-phase models. The minimum heat transfer rate is associated with β = 90° (bottom-heated) and the maximum heat transfer rate occurs in an inclination angle which varies with the Ra number.
Ahmed, Mahmoud; Eslamian, Morteza
2015-12-01
Laminar natural convection in differentially heated (β = 0°, where β is the inclination angle), inclined (β = 30° and 60°), and bottom-heated (β = 90°) square enclosures filled with a nanofluid is investigated, using a two-phase lattice Boltzmann simulation approach. The effects of the inclination angle on Nu number and convection heat transfer coefficient are studied. The effects of thermophoresis and Brownian forces which create a relative drift or slip velocity between the particles and the base fluid are included in the simulation. The effect of thermophoresis is considered using an accurate and quantitative formula proposed by the authors. Some of the existing results on natural convection are erroneous due to using wrong thermophoresis models or simply ignoring the effect. Here we show that thermophoresis has a considerable effect on heat transfer augmentation in laminar natural convection. Our non-homogenous modeling approach shows that heat transfer in nanofluids is a function of the inclination angle and Ra number. It also reveals some details of flow behavior which cannot be captured by single-phase models. The minimum heat transfer rate is associated with β = 90° (bottom-heated) and the maximum heat transfer rate occurs in an inclination angle which varies with the Ra number. PMID:26183389
NASA Astrophysics Data System (ADS)
Eid, Mohamed El Gohary
This study is combining two important and complicated processes; Enhanced Oil Recovery, EOR, from the oil rim and Enhanced Gas Recovery, EGR from the gas cap using nonhydrocarbon injection gases. EOR is proven technology that is continuously evolving to meet increased demand and oil production and desire to augment oil reserves. On the other hand, the rapid growth of the industrial and urban development has generated an unprecedented power demand, particularly during summer months. The required gas supplies to meet this demand are being stretched. To free up gas supply, alternative injectants to hydrocarbon gas are being reviewed to support reservoir pressure and maximize oil and gas recovery in oil rim reservoirs. In this study, a multi layered heterogeneous gas reservoir with an oil rim was selected to identify the most optimized development plan for maximum oil and gas recovery. The integrated reservoir characterization model and the pertinent transformed reservoir simulation history matched model were quality assured and quality checked. The development scheme is identified, in which the pattern and completion of the wells are optimized to best adapt to the heterogeneity of the reservoir. Lateral and maximum block contact holes will be investigated. The non-hydrocarbon gases considered for this study are hydrogen sulphide, carbon dioxide and nitrogen, utilized to investigate miscible and immiscible EOR processes. In November 2010, re-vaporization study, was completed successfully, the first in the UAE, with an ultimate objective is to examine the gas and condensate production in gas reservoir using non hydrocarbon gases. Field development options and proces schemes as well as reservoir management and long term business plans including phases of implementation will be identified and assured. The development option that maximizes the ultimate recovery factor will be evaluated and selected. The study achieved satisfactory results in integrating gas and oil
NASA Astrophysics Data System (ADS)
Yang, Pengju; Guo, Lixin; Jia, Chungang
2016-01-01
Based on the model of Lombardini et al. [J. Atmos. Ocean. Technol.
Oldenburg, Curtis M; Freifeld, Barry M; Pruess, Karsten; Pan, Lehua; Finsterle, Stefan; Moridis, George J
2012-12-11
In response to the urgent need for estimates of the oil and gas flow rate from the Macondo well MC252-1 blowout, we assembled a small team and carried out oil and gas flow simulations using the TOUGH2 codes over two weeks in mid-2010. The conceptual model included the oil reservoir and the well with a top boundary condition located at the bottom of the blowout preventer. We developed a fluid properties module (Eoil) applicable to a simple two-phase and two-component oil-gas system. The flow of oil and gas was simulated using T2Well, a coupled reservoir-wellbore flow model, along with iTOUGH2 for sensitivity analysis and uncertainty quantification. The most likely oil flow rate estimated from simulations based on the data available in early June 2010 was about 100,000 bbl/d (barrels per day) with a corresponding gas flow rate of 300 MMscf/d (million standard cubic feet per day) assuming the well was open to the reservoir over 30 m of thickness. A Monte Carlo analysis of reservoir and fluid properties provided an uncertainty distribution with a long tail extending down to 60,000 bbl/d of oil (170 MMscf/d of gas). The flow rate was most strongly sensitive to reservoir permeability. Conceptual model uncertainty was also significant, particularly with regard to the length of the well that was open to the reservoir. For fluid-entry interval length of 1.5 m, the oil flow rate was about 56,000 bbl/d. Sensitivity analyses showed that flow rate was not very sensitive to pressure-drop across the blowout preventer due to the interplay between gas exsolution and oil flow rate. PMID:21730177
Oldenburg, Curtis M.; Freifeld, Barry M.; Pruess, Karsten; Pan, Lehua; Finsterle, Stefan; Moridis, George J.
2012-01-01
In response to the urgent need for estimates of the oil and gas flow rate from the Macondo well MC252-1 blowout, we assembled a small team and carried out oil and gas flow simulations using the TOUGH2 codes over two weeks in mid-2010. The conceptual model included the oil reservoir and the well with a top boundary condition located at the bottom of the blowout preventer. We developed a fluid properties module (Eoil) applicable to a simple two-phase and two-component oil-gas system. The flow of oil and gas was simulated using T2Well, a coupled reservoir-wellbore flow model, along with iTOUGH2 for sensitivity analysis and uncertainty quantification. The most likely oil flow rate estimated from simulations based on the data available in early June 2010 was about 100,000 bbl/d (barrels per day) with a corresponding gas flow rate of 300 MMscf/d (million standard cubic feet per day) assuming the well was open to the reservoir over 30 m of thickness. A Monte Carlo analysis of reservoir and fluid properties provided an uncertainty distribution with a long tail extending down to 60,000 bbl/d of oil (170 MMscf/d of gas). The flow rate was most strongly sensitive to reservoir permeability. Conceptual model uncertainty was also significant, particularly with regard to the length of the well that was open to the reservoir. For fluid-entry interval length of 1.5 m, the oil flow rate was about 56,000 bbl/d. Sensitivity analyses showed that flow rate was not very sensitive to pressure-drop across the blowout preventer due to the interplay between gas exsolution and oil flow rate. PMID:21730177
Johnson, B.T.; Petty, J.D.; Huckins, J.N.; Lee, Kenneth; Gauthier, J.
2004-01-01
Phytoremediation in a simulated crude oil spill was studied with a "minimalistic" approach. The SPMD-TOX paradigm - a miniature passive sorptive device to collect and concentrate chemicals and microscale tests to detect toxicity - was used to monitor over time the bioavailability and potential toxicity of an oil spill. A simulated crude oil spill was initiated on an intertidal freshwater grass-wetland along the St. Lawrence River southwest of Quebec City, Quebec, Canada. Several phytoremediation treatments were investigated; to dissipate and ameliorate the spill, treatments included nutrient amendments with inorganic nitrogen sources (ammonium nitrate and sodium nitrate) and phosphate (super triple phosphate) with and without cut plants, with natural attenuation (no phytoremedial treatment) as a control. Sequestered oil residues were bioavailable in all oil-treated plots in Weeks 1 and 2. Interestingly, the samples were colored and fluoresced under ultraviolet light. In addition, microscale tests showed that sequestered residues were acutely toxic and genotoxic, as well as that they induced hepatic P450 enzymes. Analysis of these data suggested that polycyclic aromatic hydrocarbons were among the bioavailable residues sequestered. In addition, these findings suggested that the toxic bioavailable fractions of the oil spill and degradation products dissipated rapidly over time because after the second week the water column contained no oil or detectable degradation products in this riverine intertidal wetland. SPMD-TOX revealed no evidence of bioavailable oil products in Weeks 4, 6, 8, and 12. All phytoremediation efforts appeared to be ineffective in changing either the dissipation rate or the ability to ameliorate the oil toxicity. SPMD-TOX analysis of the water columns from these riverine experimental plots profiled the occurrence, dissipation, and influence of phytoremediation on the bioavailability and toxicity of oil products (parent or degradation products
M.J. McInerney; S.K. Maudgalya; R. Knapp; M. Folmsbee
2004-05-31
Current technology recovers only one-third to one-half of the oil that is originally present in an oil reservoir. Entrapment of petroleum hydrocarbons by capillary forces is a major factor that limits oil recovery (1, 3, 4). Hydrocarbon displacement can occur if interfacial tension (IFT) between the hydrocarbon and aqueous phases is reduced by several orders of magnitude. Microbially-produced biosurfactants may be an economical method to recover residual hydrocarbons since they are effective at low concentrations. Previously, we showed that substantial mobilization of residual hydrocarbon from a model porous system occurs at biosurfactant concentrations made naturally by B. mojavensis strain JF-1 if a polymer and 2,3-butanediol were present (2). In this report, we include data on oil recovery from Berea sandstone experiments along with our previous data from sand pack columns in order to relate biosurfactant concentration to the fraction of oil recovered. We also investigate the effect that the JF-2 biosurfactant has on interfacial tension (IFT). The presence of a co-surfactant, 2,3-butanediol, was shown to improve oil recoveries possibly by changing the optimal salinity concentration of the formulation. The JF-2 biosurfactant lowered IFT by nearly 2 orders of magnitude compared to typical values of 28-29 mN/m. Increasing the salinity increased the IFT with or without 2,3-butanediol present. The lowest interfacial tension observed was 0.1 mN/m. Tertiary oil recovery experiments showed that biosurfactant solutions with concentrations ranging from 10 to 60 mg/l in the presence of 0.1 mM 2,3-butanediol and 1 g/l of partially hydrolyzed polyacrylamide (PHPA) recovered 10-40% of the residual oil present in Berea sandstone cores. When PHPA was used alone, about 10% of the residual oil was recovered. Thus, about 10% of the residual oil recovered in these experiments was due to the increase in viscosity of the displacing fluid. Little or no oil was recovered at
Wu, Guozhong; He, Lin; Chen, Daoyi
2013-09-01
The molecular scale sorption, diffusion and distribution of asphaltene, resin, aromatic and saturate fractions of heavy crude oil on quartz surface were studied using molecular dynamic simulation. Sorption of saturates on quartz decreased by 31% when temperature increased from 298 to 398K while opposite trend was observed for resins, but insignificant changes were found in asphaltenes and aromatics. Despite of this variety, the main contribution of interactions was van der Waals energy (>90%) irrespective of molecular components and temperatures. The diffusion coefficient of saturates was predicted as 10.8×10(-10)m(2)s(-1) while that of the remaining fractions was about 4×10(-10)m(2)s(-1) at 298K. The most likely oil distribution on quartz surface was that aromatics and saturates transported randomly into and out of the complex consisting of asphaltenes surrounded by resins, which was influenced by temperature. Overall, the knowledge on quartz-oil and oil-oil interactions gained in this study is essential for future risk assessment and remediation activities as previous studies on soil remediation either limited to light oil fractions with <40 carbons or treated the heavy crude oil as a single pseudo entity ignoring the interactions between oil fractions. PMID:23632245
Sramala, Issara; Pinket, Wichchunee; Pongwan, Pawinee; Jarussophon, Suwatchai; Kasemwong, Kittiwut
2016-01-01
In this study, tea (Camellia oleifera) seed oil was formulated into self-emulsifying oil formulations (SEOF) to enhance the aqueous dispersibility and intestinal retention to achieve higher bioavailability. Self-emulsifying tea seed oils were developed by using different concentrations of lecithin in combination with surfactant blends (Span(®)80 and Tween(®)80). The lecithin/surfactant systems were able to provide clear and stable liquid formulations. The SEOF were investigated for physicochemical properties including appearance, emulsion droplets size, PDI and zeta potential. The chemical compositions of tea seed oil and SEOF were compared using GC-MS techniques. In addition, the oil adsorption measurement on artificial membranes was performed using a Franz cell apparatus and colorimetric analysis. The microscopic structure of membranes was observed with scanning electron microscopy (SEM). After aqueous dilution with fed-state simulated gastric fluid (FeSSGF), the droplet size of all SEOF was close to 200 nm with low PDI values and the zeta potential was negative. GC-MS chromatograms revealed that the chemical compositions of SEOF were not significantly different from that of the original tea seed oil. The morphological study showed that only the SEOF could form film layers. The oil droplets were extracted both from membrane treated with tea seed oil and the SEOF in order to evaluate the chemical compositions by GC-MS. PMID:27136528
Diwan, Utpal; Kovscek, Anthony R.
1999-08-09
In this investigation, existing analytical models for cyclic steam injection and oil recovery are reviewed and a new model is proposed that is applicable to horizontal wells. A new flow equation is developed for oil production during cyclic steaming of horizontal wells. The model accounts for the gravity-drainage of oil along the steam-oil interface and through the steam zone. Oil viscosity, effective permeability, geometry of the heated zone, porosity, mobile oil saturation, and thermal diffusivity of the reservoir influence the flow rate of oil in the model. The change in reservoir temperature with time is also modeled, and it results in the expected decline in oil production rate during the production cycle as the reservoir cools. Wherever appropriate, correlations and incorporated to minimize data requirements. A limited comparison to numerical simulation results agrees well, indicating that essential physics are successfully captured. Cyclic steaming appears to be a systematic met hod for heating a cold reservoir provided that a relatively uniform distribution of steam is obtained along the horizontal well during injection. A sensitivity analysis shows that the process is robust over the range of expected physical parameters.
NASA Astrophysics Data System (ADS)
Tolpeshta, I. I.; Trofimov, S. Ya.; Erkenova, M. I.; Sokolova, T. A.; Stepanov, A. L.; Lysak, L. V.; Lobanenkov, A. M.
2015-03-01
A model experiment has been performed on the successive aerobic and anaerobic degradation of oil products in samples of oil-contaminated peat sampled from a pine-subshrub-sphagnum bog near the Sutormin oilfield pipeline in the Yamal-Nenets autonomous district. During the incubation of oil-contaminated peat with lime and mineral fertilizers under complete flooding, favorable conditions are created for the aerobic oxidation of oil products at the beginning of the experiment and, as the redox potential decreases, for the anaerobic degradation of oil products conjugated with the reduction of N5+ and S+6 and methanogenesis. From the experimental data on the dynamics of the pH; Eh; and the NO{3/-}, NO{2/-}, and SO{4/2-} concentrations in the liquid phase of the samples, it has been found that denitrifiers significantly contributed to the biodegradation of oil products under the experimental conditions. After the end of the experiment, the content of oil products in the contaminated samples decreased by 21-26%.
NASA Astrophysics Data System (ADS)
Cao, Liyuan
In order to develop and apply nanotechnology in oil industry, nanoparticles transport in porous media has been studied in the past few years. Theoretical modeling were carried out to evaluate nanoparticle mobility and investigate nanoparticle retention mechanism. In this study, a simulator based on Ju and Fan's mathematical model was used to study nanoparticles transport in porous media on a reservoir scale. The simulator was verified with two simulation software, Eclipse from Schlumberger and MNM1D (Micro- and Nanoparticle transport Model in porous media in 1D geometry) developed by Tosco et al. Different injection scenarios were simulated: continuous injection, slug injection, and postflush. The effect of injection time, injection rate, and slug size on oil recovery were studied. The result discovered that when nanofluids flooding is used after water flooding as tertiary recovery method, early nanofluids injection will lead to higher oil recovery, but with more nanoparticle loss. Higher injection rate of nanofluids could help improve the flooding efficiency, but not the ultimate oil recovery for field development. Also, it can cause more nanoparticle loss. Brine water postflush is recommended when doing nanoflooding. It can significantly improve the recovery of nanoparticles, and for a homogeneous or heterogeneous reservoir, oil recovery is better compared to water flooding.
Hamidi, Hossein; Mohammadian, Erfan; Junin, Radzuan; Rafati, Roozbeh; Manan, Mohammad; Azdarpour, Amin; Junid, Mundzir
2014-02-01
Theoretically, Ultrasound method is an economical and environmentally friendly or "green" technology, which has been of interest for more than six decades for the purpose of enhancement of oil/heavy-oil production. However, in spite of many studies, questions about the effective mechanisms causing increase in oil recovery still existed. In addition, the majority of the mechanisms mentioned in the previous studies are theoretical or speculative. One of the changes that could be recognized in the fluid properties is viscosity reduction due to radiation of ultrasound waves. In this study, a technique was developed to investigate directly the effect of ultrasonic waves (different frequencies of 25, 40, 68 kHz and powers of 100, 250, 500 W) on viscosity changes of three types of oil (Paraffin oil, Synthetic oil, and Kerosene) and a Brine sample. The viscosity calculations in the smooth capillary tube were based on the mathematical models developed from the Poiseuille's equation. The experiments were carried out for uncontrolled and controlled temperature conditions. It was observed that the viscosity of all the liquids was decreased under ultrasound in all the experiments. This reduction was more significant for uncontrolled temperature condition cases. However, the reduction in viscosity under ultrasound was higher for lighter liquids compare to heavier ones. Pressure difference was diminished by decreasing in the fluid viscosity in all the cases which increases fluid flow ability, which in turn aids to higher oil recovery in enhanced oil recovery (EOR) operations. Higher ultrasound power showed higher liquid viscosity reduction in all the cases. Higher ultrasound frequency revealed higher and lower viscosity reduction for uncontrolled and controlled temperature condition experiments, respectively. In other words, the reduction in viscosity was inversely proportional to increasing the frequency in temperature controlled experiments. It was concluded that cavitation