Lessons Learned from Alternative Transportation Fuels: Modeling Transition Dynamics
Welch, C.
2006-02-01
Report focuses on understanding how analytical system modeling and data from AFV experiences could improve our understanding of the dynamic forces governing the transition to a hydrogen future.
Lessons Learned from Alternative Transportation Fuels: Modeling Transition Dynamics
Lessons Learned from Alternative Transportation Fuels: Modeling Transition Dynamics C. Welch Technical Report NREL/TP-540-39446 February 2006 Lessons Learned from Alternative Transportation Fuels: Modeling Transition Dynamics C. Welch Prepared under Task Nos. HS04.2000 and HS06.1002 Technical Report NREL/TP-540-39446 February 2006 National Renewable Energy Laboratory 1617 Cole Boulevard, Golden, Colorado 80401-3393 303-275-3000 * www.nrel.gov Operated for the U.S. Department of Energy Office of
Comparison of bifurcation dynamics of turbulent transport models for the L-H transition
Weymiens, W. Blank, H. J. de; Hogeweij, G. M. D.; Paquay, S.
2014-05-15
In more than three decades, a large amount of models and mechanisms have been proposed to describe a very beneficial feature of magnetically confined fusion plasmas: the L-H transition. Bifurcation theory can be used to compare these different models based on their dynamical transition structure. In this paper, we employ bifurcation theory to distinguish two fundamentally different descriptions of the interaction between turbulence levels and sheared flows. The analytic bifurcation analysis characterises the parameter space structure of the transition dynamics. Herewith, in these models three dynamically different types of transitions are characterised, sharp transitions, oscillatory transitions, and smooth transitions. One of the two models has a very robust transition structure and is therefore likely to be more accurate for such a robust phenomenon as the L-H transition. The other model needs more fine-tuning to get non-oscillatory transitions. These conclusions from the analytic bifurcation analysis are confirmed by dedicated numerical simulations, with the newly developed code Bifurcator.
DYNAMIC MODELING STRATEGY FOR FLOW REGIME TRANSITION IN GAS-LIQUID TWO-PHASE FLOWS
X. Wang; X. Sun; H. Zhao
2011-09-01
In modeling gas-liquid two-phase flows, the concept of flow regime has been used to characterize the global interfacial structure of the flows. Nearly all constitutive relations that provide closures to the interfacial transfers in two-phase flow models, such as the two-fluid model, are often flow regime dependent. Currently, the determination of the flow regimes is primarily based on flow regime maps or transition criteria, which are developed for steady-state, fully-developed flows and widely applied in nuclear reactor system safety analysis codes, such as RELAP5. As two-phase flows are observed to be dynamic in nature (fully-developed two-phase flows generally do not exist in real applications), it is of importance to model the flow regime transition dynamically for more accurate predictions of two-phase flows. The present work aims to develop a dynamic modeling strategy for determining flow regimes in gas-liquid two-phase flows through the introduction of interfacial area transport equations (IATEs) within the framework of a two-fluid model. The IATE is a transport equation that models the interfacial area concentration by considering the creation and destruction of the interfacial area, such as the fluid particle (bubble or liquid droplet) disintegration, boiling and evaporation; and fluid particle coalescence and condensation, respectively. For the flow regimes beyond bubbly flows, a two-group IATE has been proposed, in which bubbles are divided into two groups based on their size and shape (which are correlated), namely small bubbles and large bubbles. A preliminary approach to dynamically identifying the flow regimes is provided, in which discriminators are based on the predicted information, such as the void fraction and interfacial area concentration of small bubble and large bubble groups. This method is expected to be applied to computer codes to improve their predictive capabilities of gas-liquid two-phase flows, in particular for the applications in
Dynamic Modeling Strategy for Flow Regime Transition in Gas-Liquid Two-Phase Flows
Xia Wang; Xiaodong Sun; Benjamin Doup; Haihua Zhao
2012-12-01
In modeling gas-liquid two-phase flows, the concept of flow regimes has been widely used to characterize the global interfacial structure of the flows. Nearly all constitutive relations that provide closures to the interfacial transfers in two-phase flow models, such as the two-fluid model, are flow regime dependent. Current nuclear reactor safety analysis codes, such as RELAP5, classify flow regimes using flow regime maps or transition criteria that were developed for steady-state, fully-developed flows. As twophase flows are dynamic in nature, it is important to model the flow regime transitions dynamically to more accurately predict the two-phase flows. The present work aims to develop a dynamic modeling strategy to determine flow regimes in gas-liquid two-phase flows through introduction of interfacial area transport equations (IATEs) within the framework of a two-fluid model. The IATE is a transport equation that models the interfacial area concentration by considering the creation of the interfacial area, fluid particle (bubble or liquid droplet) disintegration, boiling and evaporation, and the destruction of the interfacial area, fluid particle coalescence and condensation. For flow regimes beyond bubbly flows, a two-group IATE has been proposed, in which bubbles are divided into two groups based on their size and shapes, namely group-1 and group-2 bubbles. A preliminary approach to dynamically identify the flow regimes is discussed, in which discriminator s are based on the predicted information, such as the void fraction and interfacial area concentration. The flow regime predicted with this method shows good agreement with the experimental observations.
Kaper, Tasso J. Kramer, Mark A.; Rotstein, Horacio G.
2013-12-15
Rhythmic neuronal oscillations across a broad range of frequencies, as well as spatiotemporal phenomena, such as waves and bumps, have been observed in various areas of the brain and proposed as critical to brain function. While there is a long and distinguished history of studying rhythms in nerve cells and neuronal networks in healthy organisms, the association and analysis of rhythms to diseases are more recent developments. Indeed, it is now thought that certain aspects of diseases of the nervous system, such as epilepsy, schizophrenia, Parkinson's, and sleep disorders, are associated with transitions or disruptions of neurological rhythms. This focus issue brings together articles presenting modeling, computational, analytical, and experimental perspectives about rhythms and dynamic transitions between them that are associated to various diseases.
Finite temperature spin-dynamics and phase transitions in spin-orbital models
Chen, C.-C.
2010-04-29
We study finite temperature properties of a generic spin-orbital model relevant to transition metal compounds, having coupled quantum Heisenberg-spin and Ising-orbital degrees of freedom. The model system undergoes a phase transition, consistent with that of a 2D Ising model, to an orbitally ordered state at a temperature set by short-range magnetic order. At low temperatures the orbital degrees of freedom freeze-out and the model maps onto a quantum Heisenberg model. The onset of orbital excitations causes a rapid scrambling of the spin spectral weight away from coherent spin-waves, which leads to a sharp increase in uniform magnetic susceptibility just below the phase transition, reminiscent of the observed behavior in the Fe-pnictide materials.
Using System Dynamics to Model the Transition to Biofuels in the United States: Preprint
Bush, B.; Duffy, M.; Sandor, D.; Peterson, S.
2008-06-01
Transitioning to a biofuels industry that is expected to displace about 30% of current U.S. gasoline consumption requires a robust biomass-to-biofuels system-of-systems that operates in concert with the existing markets. This paper discusses employing a system dynamics approach to investigate potential market penetration scenarios for cellulosic ethanol and to help government decision makers focus on areas with greatest potential.
Using System Dynamics to Model the Transition to Biofuels in the United States
Bush, B.; Duffy, M.; Sandor, D.; Peterson, S.
2008-01-01
Today, the U.S. consumes almost 21 million barrels of crude oil per day; approximately 60% of the U.S. demand is supplied by imports. The transportation sector alone accounts for two-thirds of U.S. petroleum use. Biofuels, liquid fuels produced from domestically-grown biomass, have the potential to displace about 30% of current U.S. gasoline consumption. Transitioning to a biofuels industry on this scale will require the creation of a robust biomass-to-biofuels system-of-systems that operates in concert with the existing agriculture, forestry, energy, and transportation markets. The U.S. Department of Energy is employing a system dynamics approach to investigate potential market penetration scenarios for cellulosic ethanol, and to aid decision makers in focusing government actions on the areas with greatest potential to accelerate the deployment of biofuels and ultimately reduce the nationpsilas dependence on imported oil.
Driven Skyrmions and Dynamical Transitions in Chiral Magnets...
Office of Scientific and Technical Information (OSTI)
Driven Skyrmions and Dynamical Transitions in Chiral Magnets Citation Details In-Document Search Title: Driven Skyrmions and Dynamical Transitions in Chiral Magnets Authors: Lin, ...
Energy Transition Model | Open Energy Information
Transition Model Jump to: navigation, search Tool Summary LAUNCH TOOL Name: Energy Transition Model AgencyCompany Organization: Quintel Intelligence Sector: Energy Topics:...
Dynamical transitions in large systems of mean field-coupled...
Office of Scientific and Technical Information (OSTI)
chaos and cluster states Citation Details In-Document Search Title: Dynamical transitions in large systems of mean field-coupled Landau-Stuart oscillators: Extensive chaos and ...
Dynamic Analysis of Fuel Cycle Transitioning
Brent Dixon; Steve Piet; David Shropshire; Gretchen Matthern
2009-09-01
This paper examines the time-dependent dynamics of transitioning from a once-through fuel cycle to a closed fuel cycle. The once-through system involves only Light Water Reactors (LWRs) operating on uranium oxide fuel UOX), while the closed cycle includes both LWRs and fast spectrum reactors (FRs) in either a single-tier system or two-tier fuel system. The single-tier system includes full transuranic recycle in FRs while the two-tier system adds one pass of mixed oxide uranium-plutonium (MOX U-Pu) fuel in the LWR. While the analysis primarily focuses on burner fast reactors, transuranic conversion ratios up to 1.0 are assessed and many of the findings apply to any fuel cycle transitioning from a thermal once-through system to a synergistic thermal-fast recycle system. These findings include uranium requirements for a range of nuclear electricity growth rates, the importance of back end fuel cycle facility timing and magnitude, the impact of employing a range of fast reactor conversion ratios, system sensitivity to used fuel cooling time prior to recycle, impacts on a range of waste management indicators, and projected electricity cost ranges for once-through, single-tier and two-tier systems. The study confirmed that significant waste management benefits can be realized as soon as recycling is initiated, but natural uranium savings are minimal in this century. The use of MOX in LWRs decouples the development of recycle facilities from fast reactor fielding, but also significantly delays and limits fast reactor deployment. In all cases, fast reactor deployment was significantly below than predicted by static equilibrium analyses.
Palo, P.A.; Meggitt, D.J.; Nordell, W.J.
1983-05-01
This paper presents a summary of the development and validation of undersea cable dynamics computer models by the Naval Civil Engineering Laboratory (NCEL) under the sponsorship of the Naval Facilities Engineering Command. These models allow for the analysis of both small displacement (strumming) and large displacement (static and dynamic) deformations of arbitrarily configured cable structures. All of the large displacement models described in this paper are available to the public. This paper does not emphasize the theoretical development of the models (this information is available in other references) but emphasizes the various features of the models, the comparisons between model output and experimental data, and applications for which the models have been used.
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Dynamics Model content top Chemical Supply Chain Analysis Posted by Admin on Mar 1, 2012 in | Comments 0 comments Chemical Supply Chain Analysis NISAC has developed a range of...
Dynamic pathway of the photoinduced magnetic phase transition...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Dynamic pathway of the photoinduced magnetic phase transition of multiferroic TbMnO3 Wednesday, November 25, 2015 - 3:00pm SLAC, Redtail Hawk Conference Room 108A Speaker:...
the-schedule-based-transit-model
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
The Schedule-Based transit model of the Chicago Metropolitan Area Vadim Sokolov Transportation Research and Analysis Computing Center Argonne National Laboratory List of Authors ================ Vadim Sokolov Transportation Research and Analysis Computing Center Argonne National Laboratory 277 International Drive West Chicago, IL 60185 Abstract ========= Usually public transit systems are modeled using so called frequency based approach. In this case transit route times are defined in terms of
UNIVERSALITY OF PHASE TRANSITION DYNAMICS: TOPOLOGICAL DEFECTS FROM SYMMETRY BREAKING
Zurek, Wojciech H.; Del Campo, Adolfo
2014-02-13
In the course of a non-equilibrium continuous phase transition, the dynamics ceases to be adiabatic in the vicinity of the critical point as a result of the critical slowing down (the divergence of the relaxation time in the neighborhood of the critical point). This enforces a local choice of the broken symmetry and can lead to the formation of topological defects. The Kibble-Zurek mechanism (KZM) was developed to describe the associated nonequilibrium dynamics and to estimate the density of defects as a function of the quench rate through the transition. During recent years, several new experiments investigating formation of defects in phase transitions induced by a quench both in classical and quantum mechanical systems were carried out. At the same time, some established results were called into question. We review and analyze the Kibble-Zurek mechanism focusing in particular on this surge of activity, and suggest possible directions for further progress.
Integrated Market Modeling of Hydrogen Transition Scenarios with...
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
Integrated Market Modeling of Hydrogen Transition Scenarios with HyTrans Integrated Market Modeling of Hydrogen Transition Scenarios with HyTrans Presentation by Paul Leiby of Oak ...
Linking the micro and macro: L-H transition dynamics and threshold...
Office of Scientific and Technical Information (OSTI)
Linking the micro and macro: L-H transition dynamics and threshold physics Citation Details In-Document Search Title: Linking the micro and macro: L-H transition dynamics and ...
Dynamical bifurcation as a semiclassical counterpart of a quantum phase transition
Buonsante, P.; Vezzani, A.
2011-12-15
We illustrate how dynamical transitions in nonlinear semiclassical models can be recognized as phase transitions in the corresponding--inherently linear--quantum model, where, in a statistical-mechanics framework, the thermodynamic limit is realized by letting the particle population go to infinity at fixed size. We focus on lattice bosons described by the Bose-Hubbard (BH) model and discrete self-trapping (DST) equations at the quantum and semiclassical levels, respectively. After showing that the Gaussianity of the quantum ground states is broken at the phase transition, we evaluate finite-population effects by introducing a suitable scaling hypothesis; we work out the exact value of the critical exponents and provide numerical evidence confirming our hypothesis. Our analytical results rely on a general scheme obtained from a large-population expansion of the eigenvalue equation of the BH model. In this approach the DST equations resurface as solutions of the zeroth-order problem.
HyPro: Modeling the Hydrogen Transition
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
HyPro: Modeling the H 2 Transition Brian James Directed Technologies, Inc. 9 May 2007 This presentation does not contain any proprietary, confidential, or otherwise restricted information Directed Technologies, Inc. 09 May 2007 2 Outline Model Overview Results Summary Cost & Methodology Assumptions Directed Technologies, Inc. 09 May 2007 3 Project Objectives Overall Create a tool robust enough to test the impact of different assumptions on the development of hydrogen infrastructure and
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Jaiswal, Abhishek; Egami, Takeshi; Zhang, Yang
2015-04-01
The phase behavior of multi-component metallic liquids is exceedingly complex because of the convoluted many-body and many-elemental interactions. Herein, we present systematic studies of the dynamic aspects of such a model ternary metallic liquid Cu40Zr51Al9 using molecular dynamics simulation with embedded atom method. We observed a dynamical crossover from Arrhenius to super-Arrhenius behavior in the transport properties (diffusion coefficient, relaxation times, and shear viscosity) bordered at Tx ~1300K. Unlike in many molecular and macromolecular liquids, this crossover phenomenon occurs in the equilibrium liquid state well above the melting temperature of the system (Tm ~ 900K), and the crossover temperature ismore » roughly twice of the glass-transition temperature (Tg). Below Tx, we found the elemental dynamics decoupled and the Stokes-Einstein relation broke down, indicating the onset of heterogeneous spatially correlated dynamics in the system mediated by dynamic communications among local configurational excitations. To directly characterize and visualize the correlated dynamics, we employed a non-parametric, unsupervised machine learning technique and identified dynamical clusters of atoms with similar atomic mobility. The revealed average dynamical cluster size shows an accelerated increase below Tx and mimics the trend observed in other ensemble averaged quantities that are commonly used to quantify the spatially heterogeneous dynamics such as the non-Gaussian parameter and the four-point correlation function.« less
Space and time renormalization in phase transition dynamics
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Francuz, Anna; Dziarmaga, Jacek; Gardas, Bartłomiej; Zurek, Wojciech H.
2016-02-18
Here, when a system is driven across a quantum critical point at a constant rate, its evolution must become nonadiabatic as the relaxation time τ diverges at the critical point. According to the Kibble-Zurek mechanism (KZM), the emerging post-transition excited state is characterized by a finite correlation length ξˆ set at the time tˆ=τˆ when the critical slowing down makes it impossible for the system to relax to the equilibrium defined by changing parameters. This observation naturally suggests a dynamical scaling similar to renormalization familiar from the equilibrium critical phenomena. We provide evidence for such KZM-inspired spatiotemporal scaling by investigatingmore » an exact solution of the transverse field quantum Ising chain in the thermodynamic limit.« less
First-order electroweak phase transition in the standard model...
Office of Scientific and Technical Information (OSTI)
First-order electroweak phase transition in the standard model with a low cutoff Citation Details In-Document Search Title: First-order electroweak phase transition in the standard ...
Development of a Dynamic DOE Calibration Model
Broader source: Energy.gov [DOE]
A dynamic heavy duty diesel engine model was developed. The model can be applied for calibration and control system optimization.
Pfeffer, A; Das, S; Lawless, D; Ng, B
2006-10-10
Many dynamic systems involve a number of entities that are largely independent of each other but interact with each other via a subset of state variables. We present global/local dynamic models (GLDMs) to capture these kinds of systems. In a GLDM, the state of an entity is decomposed into a globally influenced state that depends on other entities, and a locally influenced state that depends only on the entity itself. We present an inference algorithm for GLDMs called global/local particle filtering, that introduces the principle of reasoning globally about global dynamics and locally about local dynamics. We have applied GLDMs to an asymmetric urban warfare environment, in which enemy units form teams to attack important targets, and the task is to detect such teams as they form. Experimental results for this application show that global/local particle filtering outperforms ordinary particle filtering and factored particle filtering.
Modeling-Computer Simulations At Walker-Lane Transitional Zone...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Walker-Lane Transitional Zone Region (Laney, 2005) Exploration...
Modeling-Computer Simulations At Walker-Lane Transitional Zone...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Walker-Lane Transitional Zone Region (Pritchett, 2004) Exploration...
Modeling-Computer Simulations At Walker-Lane Transitional Zone...
navigation, search GEOTHERMAL ENERGYGeothermal Home Exploration Activity: Modeling-Computer Simulations At Walker-Lane Transitional Zone Region (Biasi, Et Al., 2009) Exploration...
An improved model for the transit entropy of monatomic liquids
Wallace, Duane C; Chisolm, Eric D; Bock, Nicolas
2009-01-01
In the original formulation of V-T theory for monatomic liquid dynamics, the transit contribution to entropy was taken to be a universal constant, calibrated to the constant-volume entropy of melting. This model suffers two deficiencies: (a) it does not account for experimental entropy differences of {+-}2% among elemental liquids, and (b) it implies a value of zero for the transit contribution to internal energy. The purpose of this paper is to correct these deficiencies. To this end, the V-T equation for entropy is fitted to an overall accuracy of {+-}0.1% to the available experimental high temperature entropy data for elemental liquids. The theory contains two nuclear motion contributions: (a) the dominant vibrational contribution S{sub vib}(T/{theta}{sub 0}), where T is temperature and {theta}{sub 0} is the vibrational characteristic temperature, and (b) the transit contribution S{sub tr}(T/{theta}{sub tr}), where {theta}{sub tr} is a scaling temperature for each liquid. The appearance of a common functional form of S{sub tr} for all the liquids studied is a property of the experimental data, when analyzed via the V-T formula. The resulting S{sub tr} implies the correct transit contribution to internal energy. The theoretical entropy of melting is derived, in a single formula applying to normal and anomalous melting alike. An ab initio calculation of {theta}{sub 0}, based on density functional theory, is reported for liquid Na and Cu. Comparison of these calculations with the above analysis of experimental entropy data provides verification of V-T theory. In view of the present results, techniques currently being applied in ab initio simulations of liquid properties can be employed to advantage in the further testing and development of V-T theory.
Jaiswal, Abhishek; Egami, Takeshi; Zhang, Yang
2015-04-01
The phase behavior of multi-component metallic liquids is exceedingly complex because of the convoluted many-body and many-elemental interactions. Herein, we present systematic studies of the dynamic aspects of such a model ternary metallic liquid Cu_{40}Zr_{51}Al_{9} using molecular dynamics simulation with embedded atom method. We observed a dynamical crossover from Arrhenius to super-Arrhenius behavior in the transport properties (diffusion coefficient, relaxation times, and shear viscosity) bordered at T_{x} ~1300K. Unlike in many molecular and macromolecular liquids, this crossover phenomenon occurs in the equilibrium liquid state well above the melting temperature of the system (T_{m} ~ 900K), and the crossover temperature is roughly twice of the glass-transition temperature (T_{g}). Below T_{x}, we found the elemental dynamics decoupled and the Stokes-Einstein relation broke down, indicating the onset of heterogeneous spatially correlated dynamics in the system mediated by dynamic communications among local configurational excitations. To directly characterize and visualize the correlated dynamics, we employed a non-parametric, unsupervised machine learning technique and identified dynamical clusters of atoms with similar atomic mobility. The revealed average dynamical cluster size shows an accelerated increase below T_{x} and mimics the trend observed in other ensemble averaged quantities that are commonly used to quantify the spatially heterogeneous dynamics such as the non-Gaussian parameter and the four-point correlation function.
Linking the micro and macro: L-H transition dynamics and threshold physics
Malkov, M. A. Diamond, P. H.; Miki, K.; Rice, J. E.; Tynan, G. R.
2015-03-15
The links between the microscopic dynamics and macroscopic threshold physics of the L → H transition are elucidated. Emphasis is placed on understanding the physics of power threshold scalings, and especially on understanding the minimum in the power threshold as a function of density P{sub thr} (n). By extending a numerical 1D model to evolve both electron and ion temperatures, including collisional coupling, we find that the decrease in P{sub thr} (n) along the low-density branch is due to the combination of an increase in collisional electron-to-ion energy transfer and an increase in the heating fraction coupled to the ions. Both processes strengthen the edge diamagnetic electric field needed to lock in the mean electric field shear for the L→H transition. The increase in P{sub thr} (n) along the high-density branch is due to the increase with ion collisionality of damping of turbulence-driven shear flows. Turbulence driven shear flows are needed to trigger the transition by extracting energy from the turbulence. Thus, we identify the critical transition physics components of the separatrix ion heat flux and the zonal flow excitation. The model reveals a power threshold minimum in density scans as a crossover between the threshold decrease supported by an increase in heat fraction received by ions (directly or indirectly, from electrons) and a threshold increase, supported by the rise in shear flow damping. The electron/ion heating mix emerges as important to the transition, in that it, together with electron-ion coupling, regulates the edge diamagnetic electric field shear. The importance of possible collisionless electron-ion heat transfer processes is explained.
Hydrogen Transition (HyTRANS) Model
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
optimization of an objective function that reflects private costs and benefits. ... Platform, Requirements & Availability HyTRANS is a dynamic, non- linear optimization, ...
Regional Dynamics Model (REDYN) | Open Energy Information
use the REDYN model to estimate the effects of actions and policies on people and the economy. The REDYN model powers the unique Regional Dynamics Economic Service, an...
Protein-Style Dynamical Transition in a non-Biological Polymer...
Office of Scientific and Technical Information (OSTI)
Protein-Style Dynamical Transition in a non-Biological Polymer and a non-Aqueous Solvent ... This content will become publicly available on January 1, 2017 Title: Protein-Style ...
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Brady, Nathaniel F.; Appavoo, Kannatassen; Seo, Minah; Nag, Joyeeta; Prasankumar, Rohit P.; Haglund, Richard F.; Hilton, David J.
2016-03-02
Here we report on ultrafast optical investigations of the light-induced insulator-to-metal phase transition in vanadium dioxide with controlled disorder generated by substrate mismatch. These results reveal common dynamics of this optically-induced phase transition that are independent of this disorder. Lastly, above the fluence threshold for completing the transition to the rutile crystalline phase, we find a common time scale, independent of sample morphology, of 40.5 ± 2 ps that is consistent with nucleation and growth dynamics of the R phase from the parent M1 ground state.
The electroweak phase transition in the Inert Doublet Model
Blinov, Nikita; Profumo, Stefano; Stefaniak, Tim
2015-07-21
We study the strength of a first-order electroweak phase transition in the Inert Doublet Model (IDM), where particle dark matter (DM) is comprised of the lightest neutral inert Higgs boson. We improve over previous studies in the description and treatment of the finite-temperature effective potential and of the electroweak phase transition. We focus on a set of benchmark models inspired by the key mechanisms in the IDM leading to a viable dark matter particle candidate, and illustrate how to enhance the strength of the electroweak phase transition by adjusting the masses of the yet undiscovered IDM Higgs states. We argue that across a variety of DM masses, obtaining a strong enough first-order phase transition is a generic possibility in the IDM. We find that due to direct dark matter searches and collider constraints, a sufficiently strong transition and a thermal relic density matching the universal DM abundance is possible only in the Higgs funnel regime.
Modeling Temporal Behavior in Large Networks: A Dynamic Mixed-Membership Model
Rossi, R; Gallagher, B; Neville, J; Henderson, K
2011-11-11
Given a large time-evolving network, how can we model and characterize the temporal behaviors of individual nodes (and network states)? How can we model the behavioral transition patterns of nodes? We propose a temporal behavior model that captures the 'roles' of nodes in the graph and how they evolve over time. The proposed dynamic behavioral mixed-membership model (DBMM) is scalable, fully automatic (no user-defined parameters), non-parametric/data-driven (no specific functional form or parameterization), interpretable (identifies explainable patterns), and flexible (applicable to dynamic and streaming networks). Moreover, the interpretable behavioral roles are generalizable, computationally efficient, and natively supports attributes. We applied our model for (a) identifying patterns and trends of nodes and network states based on the temporal behavior, (b) predicting future structural changes, and (c) detecting unusual temporal behavior transitions. We use eight large real-world datasets from different time-evolving settings (dynamic and streaming). In particular, we model the evolving mixed-memberships and the corresponding behavioral transitions of Twitter, Facebook, IP-Traces, Email (University), Internet AS, Enron, Reality, and IMDB. The experiments demonstrate the scalability, flexibility, and effectiveness of our model for identifying interesting patterns, detecting unusual structural transitions, and predicting the future structural changes of the network and individual nodes.
Development of one-equation transition/turbulence models
Edwards, J.R.; Roy, C.J.; Blottner, F.G.; Hassan, H.A.
2000-01-14
This paper reports on the development of a unified one-equation model for the prediction of transitional and turbulent flows. An eddy viscosity--transport equation for nonturbulent fluctuation growth based on that proposed by Warren and Hassan is combined with the Spalart-Allmaras one-equation model for turbulent fluctuation growth. Blending of the two equations is accomplished through a multidimensional intermittency function based on the work of Dhawan and Narasimha. The model predicts both the onset and extent of transition. Low-speed test cases include transitional flow over a flat plate, a single element airfoil, and a multi-element airfoil in landing configuration. High-speed test cases include transitional Mach 3.5 flow over a 5{degree} cone and Mach 6 flow over a flared-cone configuration. Results are compared with experimental data, and the grid-dependence of selected predictions is analyzed.
Quantitative Modeling of High Temperature Magnetization Dynamics
Zhang, Shufeng
2009-03-01
Final Technical Report Project title: Quantitative Modeling of High Temperature Magnetization Dynamics DOE/Office of Science Program Manager Contact: Dr. James Davenport
HyPro: Modeling the Hydrogen Transition
Office of Energy Efficiency and Renewable Energy (EERE)
Presentation by Brian James of Directed Technologies at the Joint Meeting on Hydrogen Delivery Modeling and Analysis, May 8-9, 2007
The Challenges to Coupling Dynamic Geospatial Models
Goldstein, N
2006-06-23
Many applications of modeling spatial dynamic systems focus on a single system and a single process, ignoring the geographic and systemic context of the processes being modeled. A solution to this problem is the coupled modeling of spatial dynamic systems. Coupled modeling is challenging for both technical reasons, as well as conceptual reasons. This paper explores the benefits and challenges to coupling or linking spatial dynamic models, from loose coupling, where information transfer between models is done by hand, to tight coupling, where two (or more) models are merged as one. To illustrate the challenges, a coupled model of Urbanization and Wildfire Risk is presented. This model, called Vesta, was applied to the Santa Barbara, California region (using real geospatial data), where Urbanization and Wildfires occur and recur, respectively. The preliminary results of the model coupling illustrate that coupled modeling can lead to insight into the consequences of processes acting on their own.
Very Large System Dynamics Models - Lessons Learned
Jacob J. Jacobson; Leonard Malczynski
2008-10-01
This paper provides lessons learned from developing several large system dynamics (SD) models. System dynamics modeling practice emphasize the need to keep models small so that they are manageable and understandable. This practice is generally reasonable and prudent; however, there are times that large SD models are necessary. This paper outlines two large SD projects that were done at two Department of Energy National Laboratories, the Idaho National Laboratory and Sandia National Laboratories. This paper summarizes the models and then discusses some of the valuable lessons learned during these two modeling efforts.
Integrated Market Modeling of Hydrogen Transition Scenarios with HyTrans
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
Integrated Market Modeling of Hydrogen Transition Scenarios with HyTrans Paul N. Leiby, David L. Greene and David Bowman Oak Ridge National Laboratory A presentation to the Hydrogen Delivery Analysis Meeting FreedomCAR and Fuels Partnership Delivery, Storage and Hydrogen Pathways Tech Teams May 8-9, 2007 Columbia, MD 2 OAK RIDGE NATIONAL LABORATORY U. S. DEPARTMENT OF ENERGY Drawing from several other DOE models, HyTrans integrates supply and demand in a dynamic non-linear market model to 2050.
Casten, R. F.; Bonatsos, Dennis; McCutchan, E. A.
2009-01-28
Recently, a new signature for quantum phase transitional regions has been discussed. This signature, based on degeneracies of yrast and intrinsic excitations, can distinguish first and second order phase transitions, and is valid not only at or near the analytic critical points described by X(5) and E(5), but along the phase transitional line connecting them as well. In addition, a study of a number of recent analytic solutions to the Bohr Hamiltonian and of the dynamical symmetries of the IBA Hamiltonian has revealed a set of extremely simple and general analytic formulas that describe the energies of 0{sup +} states. For the case of flat-bottomed geometrical potentials, the formula depends solely on the number of relevant dimensions. For the IBA (large boson number limit) a single formula describes all three dynamical symmetries.
Dynamics Modelling of Biolistic Gene Guns
Zhang, M.; Tao, W.; Pianetta, P.A.
2009-06-04
The gene transfer process using biolistic gene guns is a highly dynamic process. To achieve good performance, the process needs to be well understood and controlled. Unfortunately, no dynamic model is available in the open literature for analysing and controlling the process. This paper proposes such a model. Relationships of the penetration depth with the helium pressure, the penetration depth with the acceleration distance, and the penetration depth with the micro-carrier radius are presented. Simulations have also been conducted. The results agree well with experimental results in the open literature. The contribution of this paper includes a dynamic model for improving and manipulating performance of the biolistic gene gun.
Simple Dynamic Gasifier Model That Runs in Aspen Dynamics
Robinson, P.J.; Luyben, W.L.
2008-10-15
Gasification (or partial oxidation) is a vital component of 'clean coal' technology. Sulfur and nitrogen emissions can be reduced, overall energy efficiency is increased, and carbon dioxide recovery and sequestration are facilitated. Gasification units in an electric power generation plant produce a fuel for driving combustion turbines. Gasification units in a chemical plant generate gas, which can be used to produce a wide spectrum of chemical products. Future plants are predicted to be hybrid power/chemical plants with gasification as the key unit operation. The widely used process simulator Aspen Plus provides a library of models that can be used to develop an overall gasifier model that handles solids. So steady-state design and optimization studies of processes with gasifiers can be undertaken. This paper presents a simple approximate method for achieving the objective of having a gasifier model that can be exported into Aspen Dynamics. The basic idea is to use a high molecular weight hydrocarbon that is present in the Aspen library as a pseudofuel. This component should have the same 1:1 hydrogen-to-carbon ratio that is found in coal and biomass. For many plantwide dynamic studies, a rigorous high-fidelity dynamic model of the gasifier is not needed because its dynamics are very fast and the gasifier gas volume is a relatively small fraction of the total volume of the entire plant. The proposed approximate model captures the essential macroscale thermal, flow, composition, and pressure dynamics. This paper does not attempt to optimize the design or control of gasifiers but merely presents an idea of how to dynamically simulate coal gasification in an approximate way.
Development of a One-Equation Transition/Turbulence Model
EDWARDS,JACK R.; ROY,CHRISTOPHER J.; BLOTTNER,FREDERICK G.; HASSAN,HASSAN A.
2000-09-26
This paper reports on the development of a unified one-equation model for the prediction of transitional and turbulent flows. An eddy viscosity - transport equation for non-turbulent fluctuation growth based on that proposed by Warren and Hassan (Journal of Aircraft, Vol. 35, No. 5) is combined with the Spalart-Allmaras one-equation model for turbulent fluctuation growth. Blending of the two equations is accomplished through a multidimensional intermittence function based on the work of Dhawan and Narasimha (Journal of Fluid Mechanics, Vol. 3, No. 4). The model predicts both the onset and extent of transition. Low-speed test cases include transitional flow over a flat plate, a single element airfoil, and a multi-element airfoil in landing configuration. High-speed test cases include transitional Mach 3.5 flow over a 5{degree} cone and Mach 6 flow over a flared-cone configuration. Results are compared with experimental data, and the spatial accuracy of selected predictions is analyzed.
Swarm Intelligence for Urban Dynamics Modelling
Ghnemat, Rawan; Bertelle, Cyrille; Duchamp, Gerard H. E.
2009-04-16
In this paper, we propose swarm intelligence algorithms to deal with dynamical and spatial organization emergence. The goal is to model and simulate the developement of spatial centers using multi-criteria. We combine a decentralized approach based on emergent clustering mixed with spatial constraints or attractions. We propose an extension of the ant nest building algorithm with multi-center and adaptive process. Typically, this model is suitable to analyse and simulate urban dynamics like gentrification or the dynamics of the cultural equipment in urban area.
Dynamics and heterogeneity of a fate determinant during transition towards cell differentiation
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Peláez, Nicolás; Gavalda-Miralles, Arnau; Wang, Bao; Navarro, Heliodoro Tejedor; Gudjonson, Herman; Rebay, Ilaria; Dinner, Aaron R.; Katsaggelos, Aggelos K.; Amaral, Luís AN; Carthew, Richard W.
2015-11-19
Yan is an ETS-domain transcription factor responsible for maintaining Drosophila eye cells in a multipotent state. Yan is at the core of a regulatory network that determines the time and place in which cells transit from multipotency to one of several differentiated lineages. Using a fluorescent reporter for Yan expression, we observed a biphasic distribution of Yan in multipotent cells, with a rapid inductive phase and slow decay phase. Transitions to various differentiated states occurred over the course of this dynamic process, suggesting that Yan expression level does not strongly determine cell potential. Consistent with this conclusion, perturbing Yan expressionmore » by varying gene dosage had no effect on cell fate transitions. However, we observed that as cells transited to differentiation, Yan expression became highly heterogeneous and this heterogeneity was transient. Signals received via the EGF Receptor were necessary for the transience in Yan noise since genetic loss caused sustained noise. Since these signals are essential for eye cells to differentiate, we suggest that dynamic heterogeneity of Yan is a necessary element of the transition process, and cell states are stabilized through noise reduction.« less
Computational Fluid Dynamics Modeling of Diesel Engine Combustion...
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
Computational Fluid Dynamics Modeling of Diesel Engine Combustion and Emissions Computational Fluid Dynamics Modeling of Diesel Engine Combustion and Emissions 2005 Diesel Engine ...
The glass transition temperature of glassy polymers using dynamic mechanical analysis
Rodriguez, E.L.
1994-09-01
Dynamic Mechanical Analysis (DMA) is presented for four glassy polymers. Poly(vinyl acetate), poly(vinyl chloride), poly(styrene), and poly(carbonate) were studied as a function of the heating rate using ramp and step heating programs and a constant frequency of 1 Hz. The effect of frequency on the dynamic mechanical parameters was also examined from 0.01 Hz to 10 Hz. The dynamic elastic storage modulus (E{double_prime}), the dynamic elastic loss modulus (E{double_prime}) and the tan{delta} (E{double_prime}/E{prime}) were affected by both the heating rate and the frequency. Apparent activation energies for the glass transition were also determined for the four polymers which were in the range from 98 of 194 kcal/mol.
Integrated Market Modeling of Hydrogen Transition Scenarios with HyTrans |
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
Department of Energy Integrated Market Modeling of Hydrogen Transition Scenarios with HyTrans Integrated Market Modeling of Hydrogen Transition Scenarios with HyTrans Presentation by Paul Leiby of Oak Ridge National Laboratory at the Joint Meeting on Hydrogen Delivery Modeling and Analysis, May 8-9, 2007 deliv_analysis_leiby.pdf (740.06 KB) More Documents & Publications DOE Hydrogen Transition Analysis Workshop Hydrogen Policy and Analyzing the Transition Hydrogen Transition Study
Protein-style dynamical transition in a non-biological polymer and a non-aqueous solvent
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Mamontov, E.; Sharma, V. K.; Borreguero, J. M.; Tyagi, M.
2016-03-15
Using neutron scattering and molecular dynamics simulation, techniques most often associated with protein dynamical transition studies, we have investigated the microscopic dynamics of one of the most common polymers, polystyrene, which was exposed to toluene vapor, mimicking the process of protein hydration from water vapor. Polystyrene with adsorbed toluene is an example of a solvent-solute system, which, unlike biopolymers, is anhydrous and lacks hydrogen bonding. Nevertheless, it exhibits the essential traits of the dynamical transition in biomolecules, such as a specific dependence of the microscopic dynamics of both solvent and host on the temperature and the amount of solvent adsorbed.more » Ultimately, we conclude that the protein dynamical transition is a manifestation of a universal solvent-solute dynamical relationship, which is not specific to either biomolecules as solute, or aqueous media as solvent, or even a particular type of interactions between solvent and solute.« less
Bettencourt, Luis M. A.
2001-02-15
I consider several Langevin and Fokker-Planck classes of dynamics for scalar field theories in contact with a thermal bath at temperature T. These models have been applied recently in the numerical description of the dynamics of second order phase transitions and associated topological defect formation as well as in other studies of the dynamics of critical phenomena. Closed form solutions of the Fokker-Planck equation are given for a harmonic potential and a dynamical mean-field approximation is developed. These methods allow for an analytical discussion of the behavior of the theories in several circumstances of interest such as critical slowing down at a second order transition and the development of spinodal instabilities.
Model for Aggregated Water Heater Load Using Dynamic Bayesian Networks
Vlachopoulou, Maria; Chin, George; Fuller, Jason C.; Lu, Shuai; Kalsi, Karanjit
2012-07-19
The transition to the new generation power grid, or smart grid, requires novel ways of using and analyzing data collected from the grid infrastructure. Fundamental functionalities like demand response (DR), that the smart grid needs, rely heavily on the ability of the energy providers and distributors to forecast the load behavior of appliances under different DR strategies. This paper presents a new model of aggregated water heater load, based on dynamic Bayesian networks (DBNs). The model has been validated against simulated data from an open source distribution simulation software (GridLAB-D). The results presented in this paper demonstrate that the DBN model accurately tracks the load profile curves of aggregated water heaters under different testing scenarios.
Modeling of Reactor Kinetics and Dynamics
Matthew Johnson; Scott Lucas; Pavel Tsvetkov
2010-09-01
In order to model a full fuel cycle in a nuclear reactor, it is necessary to simulate the short time-scale kinetic behavior of the reactor as well as the long time-scale dynamics that occur with fuel burnup. The former is modeled using the point kinetics equations, while the latter is modeled by coupling fuel burnup equations with the kinetics equations. When the equations are solved simultaneously with a nonlinear equation solver, the end result is a code with the unique capability of modeling transients at any time during a fuel cycle.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Zachara, John M.; Chen, Xingyuan; Murray, Chris; Hammond, Glenn
2016-03-04
A well-field within a uranium (U) plume in the groundwater-surface water transition zone was monitored for a 3 year period for water table elevation and dissolved solutes. The plume discharges to the Columbia River, which displays a dramatic spring stage surge resulting from snowmelt. Groundwater exhibits a low hydrologic gradient and chemical differences with river water. River water intrudes the site in spring. Specific aims were to assess the impacts of river intrusion on dissolved uranium (Uaq), specific conductance (SpC), and other solutes, and to discriminate between transport, geochemical, and source term heterogeneity effects. As a result of water tablemore » elevation fluctuations, river water intrusion, and changes in groundwater flow directions, time series trends for Uaq and SpC were found to be complex and displayed large temporal and well-to-well variability. The wells were clustered into subsets exhibiting common behaviors resulting from the intrusion dynamics of river water and the location of source terms. Hot-spots in Uaq varied in location with increasing water table elevation through the combined effects of advection and source term location. Heuristic reactive transport modeling with PFLOTRAN demonstrated that mobilized Uaq was transported between wells and source terms in complex trajectories, and was diluted as river water entered and exited the groundwater system. Moreover, while Uaq time-series concentration trends varied significantly from year-to-year as a result of climate-caused differences in the spring hydrograph, common and partly predictable response patterns were observed that were driven by water table elevation, and the extent and duration of river water intrusion.« less
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Zachara, John M.; Chen, Xingyuan; Murray, Chris; Hammond, Glenn
2016-03-04
In this study, a well-field within a uranium (U) plume in the groundwater-surface water transition zone was monitored for a 3 year period for water table elevation and dissolved solutes. The plume discharges to the Columbia River, which displays a dramatic spring stage surge resulting from snowmelt. Groundwater exhibits a low hydrologic gradient and chemical differences with river water. River water intrudes the site in spring. Specific aims were to assess the impacts of river intrusion on dissolved uranium (Uaq), specific conductance (SpC), and other solutes, and to discriminate between transport, geochemical, and source term heterogeneity effects. Time series trendsmore » for Uaq and SpC were complex and displayed large temporal and well-to-well variability as a result of water table elevation fluctuations, river water intrusion, and changes in groundwater flow directions. The wells were clustered into subsets exhibiting common behaviors resulting from the intrusion dynamics of river water and the location of source terms. Hot-spots in Uaq varied in location with increasing water table elevation through the combined effects of advection and source term location. Heuristic reactive transport modeling with PFLOTRAN demonstrated that mobilized Uaq was transported between wells and source terms in complex trajectories, and was diluted as river water entered and exited the groundwater system. While Uaq time-series concentration trends varied significantly from year-to-year as a result of climate-caused differences in the spring hydrograph, common and partly predictable response patterns were observed that were driven by water table elevation, and the extent and duration of river water intrusion.« less
Sabin-to-Mahoney Transition Model of Quasispecies Replication
Energy Science and Technology Software Center (OSTI)
2009-05-31
Qspp is an agent-based stochastic simulation model of the Poliovirus Sabin-to-Mahoney transition. This code simulates a cell-to-cell model of Poliovirus replication. The model tracks genotypes (virus genomes) as they are replicated in cells, and as the cells burst and release particles into the medium of a culture dish. An inoculum is then taken from the pool of virions and is used to inoculate cells on a new dish. This process repeats. The Sabin genotype comprisesmore » the initial inoculum. Nucleotide positions that match the Sabin1 (vaccine strain) and Mahoney (wild type) genotypes, as well as the neurovirulent phenotype (from the literature) are enumerated as constants.« less
Phase transitions in a reaction-diffusion model on a line with boundaries
Khorrami, Mohammad Aghamohammadi, Amir
2014-03-15
A one-dimensional model on a line of length L is investigated, which involves particle diffusion as well as single particle annihilation. There are also creation and annihilation at the boundaries. The static and dynamical behaviors of the system are studied. It is seen that the system could exhibit a dynamical phase transition. For small drift velocities, the relaxation time does not depend on the absorption rates at the boundaries. This is the fast phase. For large velocities, the smaller of the absorption rates at boundaries enter the relaxation rate and makes it longer. This is the slow phase. Finally, the effect of a random particle creation in the bulk is also investigated.
Structural system identification: Structural dynamics model validation
Red-Horse, J.R.
1997-04-01
Structural system identification is concerned with the development of systematic procedures and tools for developing predictive analytical models based on a physical structure`s dynamic response characteristics. It is a multidisciplinary process that involves the ability (1) to define high fidelity physics-based analysis models, (2) to acquire accurate test-derived information for physical specimens using diagnostic experiments, (3) to validate the numerical simulation model by reconciling differences that inevitably exist between the analysis model and the experimental data, and (4) to quantify uncertainties in the final system models and subsequent numerical simulations. The goal of this project was to develop structural system identification techniques and software suitable for both research and production applications in code and model validation.
Exact results for models of multichannel quantum nonadiabatic transitions
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Sinitsyn, N. A.
2014-12-11
We consider nonadiabatic transitions in explicitly time-dependent systems with Hamiltonians of the form Hˆ(t)=Aˆ+Bˆt+Cˆ/t, where t is time and Aˆ,Bˆ,Cˆ are Hermitian N × N matrices. We show that in any model of this type, scattering matrix elements satisfy nontrivial exact constraints that follow from the absence of the Stokes phenomenon for solutions with specific conditions at t→–∞. This allows one to continue such solutions analytically to t→+∞, and connect their asymptotic behavior at t→–∞ and t→+∞. This property becomes particularly useful when a model shows additional discrete symmetries. Specifically, we derive a number of simple exact constraints and explicitmore » expressions for scattering probabilities in such systems.« less
Nanoscopic dynamics of phospholipid in unilamellar vesicles: Effect of gel to fluid phase transition
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Sharma, V. K.; Mamontov, E.; Anunciado, D. B.; O’Neill, H.; Urban, V.
2015-03-04
Dynamics of phospholipids in unilamellar vesicles (ULV) is of interest in biology, medical, and food sciences since these molecules are widely used as biocompatible agents and a mimic of cell membrane systems. We have investigated the nanoscopic dynamics of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) phospholipid in ULV as a function of temperature using elastic and quasielastic neutron scattering (QENS). The dependence of the signal on the scattering momentum transfer, which is a critical advantage of neutron scattering techniques, allows the detailed analysis of the lipid motions that cannot be carried out by other means. In agreement with a differential scanning calorimetry measurement, amore » sharp rise in the elastic scattering intensity below ca. 296 K indicates a phase transition from the high-temperature fluid phase to the low-temperature solid gel phase. The microscopic lipid dynamics exhibits qualitative differences between the solid gel phase (in a measurement at 280 K) and the fluid phase (in a measurement at a physiological temperature of 310 K). The data analysis invariably shows the presence of two distinct motions: the whole lipid molecule motion within a monolayer, or lateral diffusion, and the relatively faster internal motion of the DMPC molecule. The lateral diffusion of the whole lipid molecule is found to be Fickian in character, whereas the internal lipid motions are of localized character, consistent with the structure of the vesicles. The lateral motion slows down by an order of magnitude in the solid gel phase, whereas for the internal motion not only the time scale, but also the character of the motion changes upon the phase transition. In the solid gel phase, the lipids are more ordered and undergo uniaxial rotational motion. However, in the fluid phase, the hydrogen atoms of the lipid tails undergo confined translation diffusion rather than uniaxial rotational diffusion. The localized translational diffusion of the hydrogen
Nanoscopic Dynamics of Phospholipid in Unilamellar Vesicles: Effect of Gel to Fluid Phase Transition
Sharma, Veerendra K [ORNL; Mamontov, Eugene [ORNL; Anunciado, Divina B [ORNL; O'Neill, Hugh Michael [ORNL; Urban, Volker S [ORNL
2015-01-01
Dynamics of phospholipids in unilamellar vesicles (ULV) is of interest in biology, medical, and food sciences since these molecules are widely used as biocompatible agents and a mimic of cell membrane systems. We have investigated the nanoscopic dynamics of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) phospholipid in ULV as a function of temperature using elastic and quasielastic neutron scattering (QENS). The dependence of the signal on the scattering momentum transfer, which is a critical advantage of neutron scattering techniques, allows the detailed analysis of the lipid motions that cannot be carried out by other means. In agreement with a differential scanning calorimetry measurement, a sharp rise in the elastic scattering intensity below ca. 296 K indicates a phase transition from the high-temperature fluid phase to the low-temperature solid gel phase. The microscopic lipid dynamics exhibits qualitative differences between the solid gel phase (in a measurement at 280 K) and the fluid phase (in a measurement at a physiological temperature of 310 K). The data analysis invariably shows the presence of two distinct motions: the whole lipid molecule motion within a monolayer, or lateral diffusion, and the relatively faster internal motion of the DMPC molecule. The lateral diffusion of the whole lipid molecule is found to be Fickian in character, whereas the internal lipid motions are of localized character, consistent with the structure of the vesicles. The lateral motion slows down by an order of magnitude in the solid gel phase, whereas for the internal motion not only the time scale, but also the character of the motion changes upon the phase transition. In the solid gel phase, the lipids are more ordered and undergo uniaxial rotational motion. However, in the fluid phase, the hydrogen atoms of the lipid tails undergo confined translation diffusion rather than uniaxial rotational diffusion. The localized translational diffusion of the hydrogen atoms of
DYNAMICAL MODELING OF GALAXY MERGERS USING IDENTIKIT
Privon, G. C.; Evans, A. S.; Barnes, J. E.; Hibbard, J. E.; Yun, M. S.; Mazzarella, J. M.; Armus, L.; Surace, J.
2013-07-10
We present dynamical models of four interacting systems: NGC 5257/8, The Mice, the Antennae, and NGC 2623. The parameter space of the encounters are constrained using the Identikit model-matching and visualization tool. Identikit utilizes hybrid N-body and test particle simulations to enable rapid exploration of the parameter space of galaxy mergers. The Identikit-derived matches of these systems are reproduced with self-consistent collisionless simulations which show very similar results. The models generally reproduce the observed morphology and H I kinematics of the tidal tails in these systems with reasonable properties inferred for the progenitor galaxies. The models presented here are the first to appear in the literature for NGC 5257/8 and NGC 2623, and The Mice and the Antennae are compared with previously published models. Based on the assumed mass model and our derived initial conditions, the models indicate that the four systems are currently being viewed 175-260 Myr after first passage and cover a wide range of merger stages. In some instances there are mismatches between the models and the data (e.g., in the length of a tail); these are likely due to our adoption of a single mass model for all galaxies. Despite the use of a single mass model, these results demonstrate the utility of Identikit in constraining the parameter space for galaxy mergers when applied to real data.
Restoration of the Potosi Dynamic Model 2010
Adushita, Yasmin; Leetaru, Hannes
2014-09-30
In topical Report DOE/FE0002068-1 [2] technical performance evaluations on the Cambrian Potosi Formation were performed through reservoir modeling. The data included formation tops from mud logs, well logs from the VW1 and the CCS1 wells, structural and stratigraphic formation from three dimensional (3D) seismic data, and field data from several waste water injection wells for Potosi Formation. Intention was for two million tons per annum (MTPA) of CO2 to be injected for 20 years. In this Task the 2010 Potosi heterogeneous model (referred to as the "Potosi Dynamic Model 2010" in this report) was re-run using a new injection scenario; 3.2 MTPA for 30 years. The extent of the Potosi Dynamic Model 2010, however, appeared too small for the new injection target. It was not sufficiently large enough to accommodate the evolution of the plume. Also, it might have overestimated the injection capacity by enhancing too much the pressure relief due to the relatively close proximity between the injector and the infinite acting boundaries. The new model, Potosi Dynamic Model 2013a, was built by extending the Potosi Dynamic Model 2010 grid to 30 miles x 30 miles (48 km by 48 km), while preserving all property modeling workflows and layering. This model was retained as the base case. Potosi Dynamic Model 2013.a gives an average CO2 injection rate of 1.4 MTPA and cumulative injection of 43 Mt in 30 years, which corresponds to 45% of the injection target. This implies that according to this preliminary model, a minimum of three (3) wells could be required to achieve the injection target. The injectivity evaluation of the Potosi formation will be revisited in topical Report 15 during which more data will be integrated in the modeling exercise. A vertical flow performance evaluation could be considered for the succeeding task to determine the appropriate tubing size, the required injection tubing head pressure (THP) and to investigate whether the corresponding well injection rate
Dynamical model for light composite fermions
Derman, E.
1981-04-01
A simple dynamical model for the internal structure of the three light lepton and quark generations (..nu../sub e/,e,u,d), (..nu../sub ..mu../,..mu..,c,s), and (..nu../sub tau/,tau,t,b) is proposed. Each generation is constructed of only one fundamental massive generation F=(L-italic/sup 0/,L/sup -/,U,D) with the same (SU/sub 3/)/sub c/ x SU/sub 2/ x U/sub 1/ quantum numbers as the light generations, bound to a core of one or more massive Higgs bosons H, where H is the single physical Higgs boson necessary for spontaneous symmetry breaking in the standard model. For example, e/sup -/=L/sup -/H), ..mu../sup -/=L/sup -/HH), tau/sup -/=L/sup -/HHH). It is shown that the known binding force due to H exchange is attractive and strong enough to produce light bound states. Dynamical calculations for the bound-state composite fermions using the Bethe-Salpeter equation, together with some phenomenological imput, suggest M/sub H/approx.16 TeV and M/sub F/approx.100 GeV. It is likely that such bound states can have properties compatible with the up to now apparently elementary appearance of known fermions, for example, their Dirac magnetic moments and absence of intergeneration radiative decays (such as ..mu -->..e..gamma..). Phenomenological consequences and tests of the model are discussed.
Particle model for skyrmions in metallic chiral magnets: Dynamics...
Office of Scientific and Technical Information (OSTI)
Particle model for skyrmions in metallic chiral magnets: Dynamics, pinning, and creep Citation Details In-Document Search Title: Particle model for skyrmions in metallic chiral ...
Computational social dynamic modeling of group recruitment.
Berry, Nina M.; Lee, Marinna; Pickett, Marc; Turnley, Jessica Glicken; Smrcka, Julianne D.; Ko, Teresa H.; Moy, Timothy David; Wu, Benjamin C.
2004-01-01
The Seldon software toolkit combines concepts from agent-based modeling and social science to create a computationally social dynamic model for group recruitment. The underlying recruitment model is based on a unique three-level hybrid agent-based architecture that contains simple agents (level one), abstract agents (level two), and cognitive agents (level three). This uniqueness of this architecture begins with abstract agents that permit the model to include social concepts (gang) or institutional concepts (school) into a typical software simulation environment. The future addition of cognitive agents to the recruitment model will provide a unique entity that does not exist in any agent-based modeling toolkits to date. We use social networks to provide an integrated mesh within and between the different levels. This Java based toolkit is used to analyze different social concepts based on initialization input from the user. The input alters a set of parameters used to influence the values associated with the simple agents, abstract agents, and the interactions (simple agent-simple agent or simple agent-abstract agent) between these entities. The results of phase-1 Seldon toolkit provide insight into how certain social concepts apply to different scenario development for inner city gang recruitment.
Dynamic density field measurements of an explosively driven ????? phase transition in iron
Hull, L. M.; Gray, G. T.; Warthen, B. J.
2014-07-28
We provide a unique set of observations of the behavior of the ??? phase transition under a complex axially symmetric loading path created by sweeping a detonation wave along the end surface of a cylindrical sample. The primary data sets are the measured mass density distributions acquired at 5 independent times during the sweep of the detonation along the surface. Shocked regions and boundaries are measured, as well as regions and boundaries of elevated density (presumed to be the ??phase iron). The formation and dynamics of these regions were captured and are available for comparisons to material descriptions. We also applied 16 Photon Doppler Velocimetry probes to capture the free surface velocity along a discrete set of radially distributed points in order to compare and correlate the density measurements with previous shock wave studies. The velocimetry data are in nearly exact agreement with previous shock wave studies of the ??? phase transition, the density distributions, while generally in agreement with expectations evolved from the shock wave studies, show that the epsilon phase is generated in regions of high shear stress but at hydrostatic stresses below the typically quoted 13?GPa value. The density field measurements are particularly useful for observing the effects of the forward and reverse transformation kinetics, as well as the reverse transformation hysteresis.
LDRD final report : mesoscale modeling of dynamic loading of heterogeneous
Office of Scientific and Technical Information (OSTI)
materials. (Technical Report) | SciTech Connect LDRD final report : mesoscale modeling of dynamic loading of heterogeneous materials. Citation Details In-Document Search Title: LDRD final report : mesoscale modeling of dynamic loading of heterogeneous materials. Material response to dynamic loading is often dominated by microstructure (grain structure, porosity, inclusions, defects). An example critically important to Sandia's mission is dynamic strength of polycrystalline metals where
Filling dependence of the Mott transition in the degenerate Hubbard model
Koch, Erik; Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801 ; Gunnarsson, Olle; Martin, Richard M.
1999-12-15
Describing the doped Fullerenes using a generalized Hubbard model, we study the Mott transition for different integer fillings of the t{sub 1u} band. We use the opening of the energy-gap E{sub g} as a criterion for the transition. E{sub g} is calculated as a function of the on-site Coulomb interaction U using fixed-node diffusion Monte Carlo. We find that for systems with doping away from half filling the Mott transitions occurs at smaller U than for the half-filled system. We give a simple model for the doping dependence of the Mott transition. (c) 1999 The American Physical Society.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Cao, Gaoqing; He, Lianyi; Zhuang, Pengfei
2014-09-15
It is known that a constant magnetic field is a strong catalyst of dynamical chiral symmetry breaking in 2+1 dimensions, leading to generating dynamical fermion mass even at weakest attraction. In this work we investigate the collective modes associated with the dynamical chiral symmetry breaking in a constant magnetic field in the (2+1)-dimensional Nambu–Jona-Lasinio model with continuous U(1) chiral symmetry. We introduce a self-consistent scheme to evaluate the propagators of the collective modes at the leading order in 1/N. The contributions from the vacuum and from the magnetic field are separated such that we can employ the well-established regularization schememore » for the case of vanishing magnetic field. The same scheme can be applied to the study of the next-to-leading order correction in 1/N. We show that the sigma mode is always a lightly bound state with its mass being twice the dynamical fermion mass for arbitrary strength of the magnetic field. Since the dynamics of the collective modes is always 2+1 dimensional, the finite temperature transition should be of the Kosterlitz-Thouless (KT) type. We determine the KT transition temperature TKT as well as the mass melting temperature T* as a function of the magnetic field. It is found that the pseudogap domain TKT < T < T* is enlarged with increasing strength of the magnetic field. The influence of a chiral imbalance or axial chemical potential μ5 is also studied. We find that even a constant axial chemical potential μ5 can lead to inverse magnetic catalysis of the KT transition temperature in 2+1 dimensions. As a result, the inverse magnetic catalysis behavior is actually the de Haas–van Alphen oscillation induced by the interplay between the magnetic field and the Fermi surface.« less
Davtyan, Aram; Dama, James F.; Voth, Gregory A.; Andersen, Hans C.
2015-04-21
Coarse-grained (CG) models of molecular systems, with fewer mechanical degrees of freedom than an all-atom model, are used extensively in chemical physics. It is generally accepted that a coarse-grained model that accurately describes equilibrium structural properties (as a result of having a well constructed CG potential energy function) does not necessarily exhibit appropriate dynamical behavior when simulated using conservative Hamiltonian dynamics for the CG degrees of freedom on the CG potential energy surface. Attempts to develop accurate CG dynamic models usually focus on replacing Hamiltonian motion by stochastic but Markovian dynamics on that surface, such as Langevin or Brownian dynamics. However, depending on the nature of the system and the extent of the coarse-graining, a Markovian dynamics for the CG degrees of freedom may not be appropriate. In this paper, we consider the problem of constructing dynamic CG models within the context of the Multi-Scale Coarse-graining (MS-CG) method of Voth and coworkers. We propose a method of converting a MS-CG model into a dynamic CG model by adding degrees of freedom to it in the form of a small number of fictitious particles that interact with the CG degrees of freedom in simple ways and that are subject to Langevin forces. The dynamic models are members of a class of nonlinear systems interacting with special heat baths that were studied by Zwanzig [J. Stat. Phys. 9, 215 (1973)]. The properties of the fictitious particles can be inferred from analysis of the dynamics of all-atom simulations of the system of interest. This is analogous to the fact that the MS-CG method generates the CG potential from analysis of equilibrium structures observed in all-atom simulation data. The dynamic models generate a non-Markovian dynamics for the CG degrees of freedom, but they can be easily simulated using standard molecular dynamics programs. We present tests of this method on a series of simple examples that demonstrate that
Dynamic Model Validation with Governor Deadband on the Eastern Interconnection
Kou, Gefei; Hadley, Stanton W; Liu, Yilu
2014-04-01
This report documents the efforts to perform dynamic model validation on the Eastern Interconnection (EI) by modeling governor deadband. An on-peak EI dynamic model is modified to represent governor deadband characteristics. Simulation results are compared with synchrophasor measurements collected by the Frequency Monitoring Network (FNET/GridEye). The comparison shows that by modeling governor deadband the simulated frequency response can closely align with the actual system response.
Vasseur, Romain; Lookman, Turab; Shenoy, Subodh R.
2010-09-01
We show how microstructure can arise in first-order ferroelastic structural transitions, in two and three spatial dimensions, through a local mean-field approximation of their pseudospin Hamiltonians, that include anisotropic elastic interactions. Such transitions have symmetry-selected physical strains as their N{sub OP}-component order parameters, with Landau free energies that have a single zero-strain 'austenite' minimum at high temperatures, and spontaneous-strain 'martensite' minima of N{sub V} structural variants at low temperatures. The total free energy also has gradient terms, and power-law anisotropic effective interactions, induced by 'no-dislocation' St Venant compatibility constraints. In a reduced description, the strains at Landau minima induce temperature dependent, clocklike Z{sub N{sub V+1}} Hamiltonians, with N{sub OP}-component strain-pseudospin vectors S-vector pointing to N{sub V}+1 discrete values (including zero). We study elastic texturing in five such first-order structural transitions through a local mean-field approximation of their pseudospin Hamiltonians, that include the power-law interactions. As a prototype, we consider the two-variant square/rectangle transition, with a one-component pseudospin taking N{sub V}+1=3 values of S=0,{+-}1, as in a generalized Blume-Capel model. We then consider transitions with two-component (N{sub OP}=2) pseudospins: the equilateral to centered rectangle (N{sub V}=3); the square to oblique polygon (N{sub V}=4); the triangle to oblique (N{sub V}=6) transitions; and finally the three-dimensional (3D) cubic to tetragonal transition (N{sub V}=3). The local mean-field solutions in two-dimensional and 3D yield oriented domain-wall patterns as from continuous-variable strain dynamics, showing the discrete-variable models capture the essential ferroelastic texturings. Other related Hamiltonians illustrate that structural transitions in materials science can be the source of interesting spin models in
Dynamic Modeling in Solid-Oxide Fuel Cells Controller Design
Lu, Ning; Li, Qinghe; Sun, Xin; Khaleel, Mohammad A.
2007-06-28
In this paper, a dynamic model of the solid-oxide fuel cell (SOFC) power unit is developed for the purpose of designing a controller to regulate fuel flow rate, fuel temperature, air flow rate, and air temperature to maintain the SOFC stack temperature, fuel utilization rate, and voltage within operation limits. A lumped model is used to consider the thermal dynamics and the electro-chemial dynamics inside an SOFC power unit. The fluid dynamics at the fuel and air inlets are considered by using the in-flow ramp-rates.
Developing a Dynamic Pharmacophore Model for HIV-1 Integrase
Carlson, Heather A.; Masukawa, Keven M.; Rubins, Kathleen; Bushman, Frederic; Jorgensen, William L.; Lins, Roberto; Briggs, James; Mccammon, Andy
2000-05-11
We present the first receptor-based pharmacophore model for HIV-1 integrase. The development of ''dynamic'' pharmacophore models is a new method that accounts for the inherent flexibility of the active site and aims to reduce the entropic penalties associated with binding a ligand. Furthermore, this new drug discovery method overcomes the limitation of an incomplete crystal structure of the target protein. A molecular dynamics (MD) simulation describes the flexibility of the uncomplexed protein. Many conformational models of the protein are saved from the MD simulations and used in a series of multi-unit search for interacting conformers (MUSIC) simulations. MUSIC is a multiple-copy minimization method, available in the BOSS program; it is used to determine binding regions for probe molecules containing functional groups that complement the active site. All protein conformations from the MD are overlaid, and conserved binding regions for the probe molecules are identified. Those conserved binding regions define the dynamic pharmacophore model. Here, the dynamic model is compared to known inhibitors of the integrase as well as a three-point, ligand-based pharmacophore model from the literature. Also, a ''static'' pharmacophore model was determined in the standard fashion, using a single crystal structure. Inhibitors thought to bind in the active site of HIV-1 integrase fit the dynamic model but not the static model. Finally, we have identified a set of compounds from the Available Chemicals Directory that fit the dynamic pharmacophore model, and experimental testing of the compounds has confirmed several new inhibitors.
Dynamic model predicts well bore surge and swab pressures
Bing, Z.; Kaiji, Z.
1996-12-30
A dynamic well control model predicts surge and swab pressures more accurately than a steady-state model, thereby providing better estimates of pressure fluctuations when pipe is tripped. Pressure fluctuations from tripping pipe into a well can contribute to lost circulation, kicks,and well control problems. This dynamic method of predicting surge and swab pressures was verified in a full-scale test well in the Zhong Yuan oil field in China. Both the dynamic model and steady state model were verified through the test data. The test data showed the dynamic model can correctly predict downhole pressures from running or pulling pipe in a well; steady state models may result in relatively large prediction errors, especially in deeper wells.
Computational Fluid Dynamics Modeling of Diesel Engine Combustion and
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
Emissions | Department of Energy Computational Fluid Dynamics Modeling of Diesel Engine Combustion and Emissions Computational Fluid Dynamics Modeling of Diesel Engine Combustion and Emissions 2005 Diesel Engine Emissions Reduction (DEER) Conference Presentations and Posters 2005_deer_reitz.pdf (682.47 KB) More Documents & Publications Experiments and Modeling of Two-Stage Combustion in Low-Emissions Diesel Engines Comparison of Conventional Diesel and Reactivity Controlled Compression
Model based control of dynamic atomic force microscope
Lee, Chibum; Salapaka, Srinivasa M.
2015-04-15
A model-based robust control approach is proposed that significantly improves imaging bandwidth for the dynamic mode atomic force microscopy. A model for cantilever oscillation amplitude and phase dynamics is derived and used for the control design. In particular, the control design is based on a linearized model and robust H{sub ∞} control theory. This design yields a significant improvement when compared to the conventional proportional-integral designs and verified by experiments.
Generic solar photovoltaic system dynamic simulation model specification.
Ellis, Abraham; Behnke, Michael Robert; Elliott, Ryan Thomas
2013-10-01
This document is intended to serve as a specification for generic solar photovoltaic (PV) system positive-sequence dynamic models to be implemented by software developers and approved by the WECC MVWG for use in bulk system dynamic simulations in accordance with NERC MOD standards. Two specific dynamic models are included in the scope of this document. The first, a Central Station PV System model, is intended to capture the most important dynamic characteristics of large scale (> 10 MW) PV systems with a central Point of Interconnection (POI) at the transmission level. The second, a Distributed PV System model, is intended to represent an aggregation of smaller, distribution-connected systems that comprise a portion of a composite load that might be modeled at a transmission load bus.
Modeling ramp compression experiments using large-scale molecular dynamics simulation.
Mattsson, Thomas Kjell Rene; Desjarlais, Michael Paul; Grest, Gary Stephen; Templeton, Jeremy Alan; Thompson, Aidan Patrick; Jones, Reese E.; Zimmerman, Jonathan A.; Baskes, Michael I.; Winey, J. Michael; Gupta, Yogendra Mohan; Lane, J. Matthew D.; Ditmire, Todd; Quevedo, Hernan J.
2011-10-01
Molecular dynamics simulation (MD) is an invaluable tool for studying problems sensitive to atomscale physics such as structural transitions, discontinuous interfaces, non-equilibrium dynamics, and elastic-plastic deformation. In order to apply this method to modeling of ramp-compression experiments, several challenges must be overcome: accuracy of interatomic potentials, length- and time-scales, and extraction of continuum quantities. We have completed a 3 year LDRD project with the goal of developing molecular dynamics simulation capabilities for modeling the response of materials to ramp compression. The techniques we have developed fall in to three categories (i) molecular dynamics methods (ii) interatomic potentials (iii) calculation of continuum variables. Highlights include the development of an accurate interatomic potential describing shock-melting of Beryllium, a scaling technique for modeling slow ramp compression experiments using fast ramp MD simulations, and a technique for extracting plastic strain from MD simulations. All of these methods have been implemented in Sandia's LAMMPS MD code, ensuring their widespread availability to dynamic materials research at Sandia and elsewhere.
Dynamical Models of the Excitations of Nucleon Resonances
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
T. Sato; Lee, T.-S. H.
2009-05-01
The development of a dynamical model for investigating the nucleon resonances using the reactions of meson production frommore » $$\\pi N$$, $$\\gamma N$$, $N(e,e')$, and $$N(\
Developing Generic Dynamic Models for the 2030 Eastern Interconnection Grid
Kou, Gefei; Hadley, Stanton W; Markham, Penn N; Liu, Yilu
2013-12-01
The Eastern Interconnection Planning Collaborative (EIPC) has built three major power flow cases for the 2030 Eastern Interconnection (EI) based on various levels of energy/environmental policy conditions, technology advances, and load growth. Using the power flow cases, this report documents the process of developing the generic 2030 dynamic models using typical dynamic parameters. The constructed model was validated indirectly using the synchronized phasor measurements by removing the wind generation temporarily.
Capacity planning in a transitional economy: What issues? Which models?
Mubayi, V.; Leigh, R.W.; Bright, R.N.
1996-03-01
This paper is devoted to an exploration of the important issues facing the Russian power generation system and its evolution in the foreseeable future and the kinds of modeling approaches that capture those issues. These issues include, for example, (1) trade-offs between investments in upgrading and refurbishment of existing thermal (fossil-fired) capacity and safety enhancements in existing nuclear capacity versus investment in new capacity, (2) trade-offs between investment in completing unfinished (under construction) projects based on their original design versus investment in new capacity with improved design, (3) incorporation of demand-side management options (investments in enhancing end-use efficiency, for example) within the planning framework, (4) consideration of the spatial dimensions of system planning including investments in upgrading electric transmission networks or fuel shipment networks and incorporating hydroelectric generation, (5) incorporation of environmental constraints and (6) assessment of uncertainty and evaluation of downside risk. Models for exploring these issues include low power shutdown (LPS) which are computationally very efficient, though approximate, and can be used to perform extensive sensitivity analyses to more complex models which can provide more detailed answers but are computationally cumbersome and can only deal with limited issues. The paper discusses which models can usefully treat a wide range of issues within the priorities facing decision makers in the Russian power sector and integrate the results with investment decisions in the wider economy.
D’Abramo, Marco; Dipartimento di Chimica, Universitá Sapienza, P.le Aldo Moro, 5, 00185, Rome ; Aschi, Massimiliano; Amadei, Andrea
2014-04-28
Here, we extend a recently introduced theoretical-computational procedure [M. D’Alessandro, M. Aschi, C. Mazzuca, A. Palleschi, and A. Amadei, J. Chem. Phys. 139, 114102 (2013)] to include quantum vibrational transitions in modelling electronic spectra of atomic molecular systems in condensed phase. The method is based on the combination of Molecular Dynamics simulations and quantum chemical calculations within the Perturbed Matrix Method approach. The main aim of the presented approach is to reproduce as much as possible the spectral line shape which results from a subtle combination of environmental and intrinsic (chromophore) mechanical-dynamical features. As a case study, we were able to model the low energy UV-vis transitions of pyrene in liquid acetonitrile in good agreement with the experimental data.
Huse, Nils; Khalil, Munira; Kim, Tae Kyu; Smeigh, Amanda L.; Jamula, Lindsey; McCusker, James K.; Schoenlein, Robert W.
2009-05-24
We report measurements of the photo-induced Fe(II) spin crossover reaction dynamics in solution via time-resolved x-ray absorption spectroscopy. EXAFS measurements reveal that the iron?nitrogen bond lengthens by 0.21+-0.03 Angstrom in the high-spin transient excited state relative to the ground state. XANES measurements at the Fe L-edge show directly the influence of the structural change on the ligand-field splitting of the Fe(II) 3d orbitals associated with the spin transition.
Multichannel Pseudogap Kondo Model: Large-N Solution and Quantum-Critical Dynamics
Vojta, Matthias
2001-08-27
We discuss a multichannel SU(N) Kondo model which displays nontrivial zero-temperature phase transitions due to a conduction electron density of states vanishing with a power law at the Fermi level. In a particular large-N limit, the system is described by coupled integral equations corresponding to a dynamic saddle point. We exactly determine the universal low-energy behavior of spectral densities at the scale-invariant fixed points, obtain anomalous exponents, and compute scaling functions describing the crossover near the quantum-critical points. We argue that our findings are relevant to recent experiments on impurity-doped d -wave superconductors.
Approximate Bisimulation-Based Reduction of Power System Dynamic Models
Stankovic, AM; Dukic, SD; Saric, AT
2015-05-01
In this paper we propose approximate bisimulation relations and functions for reduction of power system dynamic models in differential- algebraic (descriptor) form. The full-size dynamic model is obtained by linearization of the nonlinear transient stability model. We generalize theoretical results on approximate bisimulation relations and bisimulation functions, originally derived for a class of constrained linear systems, to linear systems in descriptor form. An algorithm for transient stability assessment is proposed and used to determine whether the power system is able to maintain the synchronism after a large disturbance. Two benchmark power systems are used to illustrate the proposed algorithm and to evaluate the applicability of approximate bisimulation relations and bisimulation functions for reduction of the power system dynamic models.
Modeling and simulation of consumer response to dynamic pricing.
Valenzuela, J.; Thimmapuram, P.; Kim, J
2012-08-01
Assessing the impacts of dynamic-pricing under the smart grid concept is becoming extremely important for deciding its full deployment. In this paper, we develop a model that represents the response of consumers to dynamic pricing. In the model, consumers use forecasted day-ahead prices to shift daily energy consumption from hours when the price is expected to be high to hours when the price is expected to be low while maintaining the total energy consumption as unchanged. We integrate the consumer response model into the Electricity Market Complex Adaptive System (EMCAS). EMCAS is an agent-based model that simulates restructured electricity markets. We explore the impacts of dynamic-pricing on price spikes, peak demand, consumer energy bills, power supplier profits, and congestion costs. A simulation of an 11-node test network that includes eight generation companies and five aggregated consumers is performed for a period of 1 month. In addition, we simulate the Korean power system.
Langeloth, Michael; Böhm, Michael C.; Müller-Plathe, Florian; Sugii, Taisuke
2015-12-28
We investigate the volumetric glass transition temperature T{sub g} in epoxy thermosets by means of molecular dynamics simulations. The epoxy thermosets consist of the resin bisphenol A diglycidyl ether and the hardener diethylenetriamine. A structure based coarse-grained (CG) force field has been derived using iterative Boltzmann inversion in order to facilitate simulations of larger length scales. We observe that T{sub g} increases clearly with the degree of cross-linking for all-atomistic (AA) and CG simulations. The transition T{sub g} in CG simulations of uncured mixtures is much lower than in AA-simulations due to the soft nature of the CG potentials, but increases all the more with the formation of rigid cross-links. Additional simulations of the CG mixtures in contact with a surface show the existence of an interphase region of about 3 nm thickness in which the network properties deviate significantly from the bulk. In accordance to experimental studies, we observe that T{sub g} is reduced in this interphase region and gradually increases to its bulk value with distance from the surface. The present study shows that the glass transition is a local phenomenon that depends on the network structure in the immediate environment.
Modeling of Endothelial Glyccalyx via Dissipative Particle Dynamics |
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Argonne Leadership Computing Facility Modeling of Endothelial Glyccalyx via Dissipative Particle Dynamics Authors: Mingge, D., Haojun, L., Karniadakis, G We employ Dissipative Particle Dynamics (DPD) to simulate flow in small vessels with the endothelial glycocalyx attached to the wall. Of particular interest is the quantification of the slip velocity at the edge of glycocalyx and of the increased pressure drop at different crafting densities, stiffness and height of the glycocalyx. Results
Merritt, T. R.; Meeker, M. A.; Magill, B. A.; Khodaparast, G. A.; McGill, S.; Tischler, J. G.; Choi, S. G.; Palmstrøm, C. J.
2014-05-21
The excitonic radiative transitions of InAs{sub x}P{sub 1−x} (x = 0.13 and x = 0.40) alloy epitaxial layers were studied through magnetic field and temperature dependent photoluminescence and time-resolved photoluminescence spectroscopy. While the linewidth and lineshape of the exciton transition for x = 0.40 indicate the presence of alloy broadening due to random anion distribution and the existence of localized exciton states, those of x = 0.13 suggest that this type of compositional disorder is absent in x = 0.13. This localization is further supported by the behavior of the exciton transitions at low temperature and high magnetic fields. InAs{sub 0.4}P{sub 0.6} exhibits anomalous “S-shaped” temperature dependence of the excition emission peak below 100 K as well as linewidth broadening at high magnetic fields due to the compression of the excitonic volume amid compositional fluctuations. Finally, photoluminescence decay patterns suggest that the excitons radiatively relax through two channels, a fast and a slow decay. While the lifetime of the fast decay is comparable for both compositions (∼30 ps), that of the slow decay increases from 206 ps to 427 ps as x increases from 0.13 to 0.40, attributable to carrier migration between the localization states of InAs{sub 0.4}P{sub 0.6}.
DETAILED CHEMICAL KINETIC MODELING OF ISO-OCTANE SI-HCCI TRANSITION
Havstad, M A; Aceves, S M; McNenly, M J; Piggott, W T; Edwards, K D; Wagner, R M; Daw, C S; Finney, C A
2009-10-12
The authors describe a CHEMKIN-based multi-zone model that simulates the expected combustion variations in a single-cylinder engine fueled with iso-octane as the engine transitions from spark-ignited (ST) combustion to homogeneous charge compression ignition (HCCI) combustion. The model includes a 63-species reaction mechanism and mass and energy balances for the cylinder and the exhaust flow. For this study they assumed that the SI-to-HCCI transition is implemented by means of increasing the internal exhaust gas recirculation (EGR) at constant engine speed. This transition scneario is consistent with that implemented in previously reported experimental measurements on an experimental engine equipped with variable valve actuation. They find that the model captures many of the important experimental trends, including stable SI combustion at low EGR ({approx} 0.10), a transition to highly unstable combustion at intermediate EGR, and finally stable HCCI combustion at very high EGR ({approx} 0.75). Remaining differences between the predicted and experimental instability patterns indicate that there is further room for model improvement.
Detailed Chemical Kinetic Modeling of Iso-octane SI-HCCI Transition
Havstad, Mark A; Aceves, Salvador M; McNenly, Matthew J; Piggott, William T; Edwards, Kevin Dean; Wagner, Robert M; Daw, C Stuart; FINNEY, Charles E A
2010-01-01
We describe a CHEMKIN-based multi-zone model that simulates the expected combustion variations in a single-cylinder engine fueled with iso-octane as the engine transitions from spark-ignited (SI) combustion to homogenous charge compression ignition (HCCI) combustion. The model includes a 63-species reaction mechanism and mass and energy balances for the cylinder and the exhaust flow. For this study we assumed that the SI-to-HCCI transition is implemented by means of increasing the internal exhaust gas recirculation (EGR) at constant engine speed. This transition scenario is consistent with that implemented in previously reported experimental measurements on an experimental engine equipped with variable valve actuation. We find that the model captures many of the important experimental trends, including stable SI combustion at low EGR (-0.10), a transition to highly unstable combustion at intermediate EGR, and finally stable HCCI combustion at very high EGR (-0.75). Remaining differences between the predicted and experimental instability patterns indicate that there is further room for model improvement.
Magnetic field induced quantum phase transitions in the two-impurity Anderson model
Zhu, Lujun; Zhu, Jian - Xin
2010-11-17
In the two-impurity Anderson model, the inter-impurity spin exchange interaction favors a spin singlet state between two impurities leading to the localization of quasiparticles. We show that a local uniform magnetic field can delocalize the quasiparticies to restore the Kondo resonance. This transition is found to be continuous, accompanied by not only the divergence of the staggered (anti ferromagnetic) susceptibility, but also the divergence of the uniform spin susceptibility. This implies that the magnetic field induced quantum phase transitions in Kondo systems are in favor of the local critical type.
Li, W. C.; Song, X.; Feng, J. J.; Zeng, M.; Gao, X. S.; Qin, M. H.; Jia, X. T.
2015-07-07
In this work, the effects of the random exchange interaction on the phase transitions and phase diagrams of classical frustrated Heisenberg model are investigated by Monte Carlo simulation in order to simulate the chemical doping effect in real materials. It is observed that the antiferromagnetic transitions shift toward low temperature with the increasing magnitude of the random exchange interaction, which can be qualitatively understood from the competitions among local spin states. This study is related to the magnetic properties in the doped iron-based superconductors.
A model of lipid-free Apolipoprotein A-I revealed by iterative molecular dynamics simulation
Zhang, Xing; Lei, Dongsheng; Zhang, Lei; Rames, Matthew; Zhang, Shengli
2015-03-20
Apolipoprotein A-I (apo A-I), the major protein component of high-density lipoprotein, has been proven inversely correlated to cardiovascular risk in past decades. The lipid-free state of apo A-I is the initial stage which binds to lipids forming high-density lipoprotein. Molecular models of lipid-free apo A-I have been reported by methods like X-ray crystallography and chemical cross-linking/mass spectrometry (CCL/MS). Through structural analysis we found that those current models had limited consistency with other experimental results, such as those from hydrogen exchange with mass spectrometry. Through molecular dynamics simulations, we also found those models could not reach a stable equilibrium state. Therefore, by integrating various experimental results, we proposed a new structural model for lipidfree apo A-I, which contains a bundled four-helix N-terminal domain (1–192) that forms a variable hydrophobic groove and a mobile short hairpin C-terminal domain (193–243). This model exhibits an equilibrium state through molecular dynamics simulation and is consistent with most of the experimental results known from CCL/MS on lysine pairs, fluorescence resonance energy transfer and hydrogen exchange. This solution-state lipid-free apo A-I model may elucidate the possible conformational transitions of apo A-I binding with lipids in high-density lipoprotein formation.
A model of lipid-free Apolipoprotein A-I revealed by iterative molecular dynamics simulation
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Zhang, Xing; Lei, Dongsheng; Zhang, Lei; Rames, Matthew; Zhang, Shengli
2015-03-20
Apolipoprotein A-I (apo A-I), the major protein component of high-density lipoprotein, has been proven inversely correlated to cardiovascular risk in past decades. The lipid-free state of apo A-I is the initial stage which binds to lipids forming high-density lipoprotein. Molecular models of lipid-free apo A-I have been reported by methods like X-ray crystallography and chemical cross-linking/mass spectrometry (CCL/MS). Through structural analysis we found that those current models had limited consistency with other experimental results, such as those from hydrogen exchange with mass spectrometry. Through molecular dynamics simulations, we also found those models could not reach a stable equilibrium state. Therefore,more » by integrating various experimental results, we proposed a new structural model for lipidfree apo A-I, which contains a bundled four-helix N-terminal domain (1–192) that forms a variable hydrophobic groove and a mobile short hairpin C-terminal domain (193–243). This model exhibits an equilibrium state through molecular dynamics simulation and is consistent with most of the experimental results known from CCL/MS on lysine pairs, fluorescence resonance energy transfer and hydrogen exchange. This solution-state lipid-free apo A-I model may elucidate the possible conformational transitions of apo A-I binding with lipids in high-density lipoprotein formation.« less
Dynamic crack initiation toughness : experiments and peridynamic modeling.
Foster, John T.
2009-10-01
This is a dissertation on research conducted studying the dynamic crack initiation toughness of a 4340 steel. Researchers have been conducting experimental testing of dynamic crack initiation toughness, K{sub Ic}, for many years, using many experimental techniques with vastly different trends in the results when reporting K{sub Ic} as a function of loading rate. The dissertation describes a novel experimental technique for measuring K{sub Ic} in metals using the Kolsky bar. The method borrows from improvements made in recent years in traditional Kolsky bar testing by using pulse shaping techniques to ensure a constant loading rate applied to the sample before crack initiation. Dynamic crack initiation measurements were reported on a 4340 steel at two different loading rates. The steel was shown to exhibit a rate dependence, with the recorded values of K{sub Ic} being much higher at the higher loading rate. Using the knowledge of this rate dependence as a motivation in attempting to model the fracture events, a viscoplastic constitutive model was implemented into a peridynamic computational mechanics code. Peridynamics is a newly developed theory in solid mechanics that replaces the classical partial differential equations of motion with integral-differential equations which do not require the existence of spatial derivatives in the displacement field. This allows for the straightforward modeling of unguided crack initiation and growth. To date, peridynamic implementations have used severely restricted constitutive models. This research represents the first implementation of a complex material model and its validation. After showing results comparing deformations to experimental Taylor anvil impact for the viscoplastic material model, a novel failure criterion is introduced to model the dynamic crack initiation toughness experiments. The failure model is based on an energy criterion and uses the K{sub Ic} values recorded experimentally as an input. The failure model
Integrated dynamic modeling and management system mission analysis
Lee, A.K.
1994-12-28
This document summarizes the mission analysis performed on the Integrated Dynamic Modeling and Management System (IDMMS). The IDMMS will be developed to provide the modeling and analysis capability required to understand the TWRS system behavior in terms of the identified TWRS performance measures. The IDMMS will be used to demonstrate in a verified and validated manner the satisfactory performance of the TWRS system configuration and assurance that the requirements have been satisfied.
Dynamic Metabolic Modeling of Denitrifying Bacterial Growth: The Cybernetic Approach
Song, Hyun-Seob; Liu, Chongxuan
2015-06-29
Denitrification is a multistage reduction process converting nitrate ultimately to nitrogen gas, carried out mostly by facultative bacteria. Modeling of the denitrification process is challenging due to the complex metabolic regulation that modulates sequential formation and consumption of a series of nitrogen oxide intermediates, which serve as the final electron acceptors for denitrifying bacteria. In this work, we examined the effectiveness and accuracy of the cybernetic modeling framework in simulating the growth dynamics of denitrifying bacteria in comparison with kinetic models. In four different case studies using the literature data, we successfully simulated diauxic and triauxic growth patterns observed in anoxic and aerobic conditions, only by tuning two or three parameters. In order to understand the regulatory structure of the cybernetic model, we systematically analyzed the effect of cybernetic control variables on simulation accuracy. The results showed that the consideration of both enzyme synthesis and activity control through u- and v-variables is necessary and relevant and that uvariables are of greater importance in comparison to v-variables. In contrast, simple kinetic models were unable to accurately capture dynamic metabolic shifts across alternative electron acceptors, unless an inhibition term was additionally incorporated. Therefore, the denitrification process represents a reasonable example highlighting the criticality of considering dynamic regulation for successful metabolic modeling.
Shell model method for Gamow-Teller transitions in heavy, deformed nuclei
Gao Zaochun [Joint Institute for Nuclear Astrophysics and Department of Physics, University of Notre Dame, Notre Dame, Indiana 46556 (United States); Department of Physics and Astronomy and National Superconducting Cyclotron Laboratory, Michigan State University, East Lansing, Michigan 48824 (United States); China Institute of Atomic Energy, P.O. Box 275 (18), Beijing 102413 (China); Sun Yang [Joint Institute for Nuclear Astrophysics and Department of Physics, University of Notre Dame, Notre Dame, Indiana 46556 (United States); Chen, Y.-S. [China Institute of Atomic Energy, P.O. Box 275(18), Beijing 102413 (China); Institute of Theoretical Physics, Academia Sinica, Beijing 100080 (China)
2006-11-15
A method for calculation of Gamow-Teller transition rates is developed by using the concept of the Projected Shell Model (PSM). The shell model basis is constructed by superimposing angular-momentum-projected multiquasiparticle configurations, and nuclear wave functions are obtained by diagonalizing the two-body interactions in these projected states. Calculation of transition matrix elements in the PSM framework is discussed in detail, and the effects caused by the Gamow-Teller residual forces and by configuration-mixing are studied. With this method, it may become possible to perform a state-by-state calculation for {beta}-decay and electron-capture rates in heavy, deformed nuclei at finite temperatures. Our first example indicates that, while experimentally known Gamow-Teller transition rates from the ground state of the parent nucleus are reproduced, stronger transitions from some low-lying excited states are predicted to occur, which may considerably enhance the total decay rates once these nuclei are exposed to hot stellar environments.
Gauge turbulence, topological defect dynamics, and condensation in Higgs models
Gasenzer, Thomas; McLerran, Larry; Pawlowski, Jan M.; Sexty, Dénes
2014-07-28
The real-time dynamics of topological defects and turbulent configurations of gauge fields for electric and magnetic confinement are studied numerically within a 2+1D Abelian Higgs model. It is shown that confinement is appearing in such systems equilibrating after a strong initial quench such as the overpopulation of the infrared modes. While the final equilibrium state does not support confinement, metastable vortex defect configurations appearing in the gauge field are found to be closely related to the appearance of physically observable confined electric and magnetic charges. These phenomena are seen to be intimately related to the approach of a non-thermal fixed point of the far-from-equilibrium dynamical evolution, signaled by universal scaling in the gauge-invariant correlation function of the Higgs field. Even when the parameters of the Higgs action do not support condensate formation in the vacuum, during this approach, transient Higgs condensation is observed. We discuss implications of these results for the far-from-equilibrium dynamics of Yang–Mills fields and potential mechanisms of how confinement and condensation in non-Abelian gauge fields can be understood in terms of the dynamics of Higgs models. These suggest that there is an interesting new class of dynamics of strong coherent turbulent gauge fields with condensates.
Gauge turbulence, topological defect dynamics, and condensation in Higgs models
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Gasenzer, Thomas; McLerran, Larry; Pawlowski, Jan M.; Sexty, Dénes
2014-07-28
The real-time dynamics of topological defects and turbulent configurations of gauge fields for electric and magnetic confinement are studied numerically within a 2+1D Abelian Higgs model. It is shown that confinement is appearing in such systems equilibrating after a strong initial quench such as the overpopulation of the infrared modes. While the final equilibrium state does not support confinement, metastable vortex defect configurations appearing in the gauge field are found to be closely related to the appearance of physically observable confined electric and magnetic charges. These phenomena are seen to be intimately related to the approach of a non-thermal fixedmore » point of the far-from-equilibrium dynamical evolution, signaled by universal scaling in the gauge-invariant correlation function of the Higgs field. Even when the parameters of the Higgs action do not support condensate formation in the vacuum, during this approach, transient Higgs condensation is observed. We discuss implications of these results for the far-from-equilibrium dynamics of Yang–Mills fields and potential mechanisms of how confinement and condensation in non-Abelian gauge fields can be understood in terms of the dynamics of Higgs models. These suggest that there is an interesting new class of dynamics of strong coherent turbulent gauge fields with condensates.« less
Dynamic simulation models and performance of an OTEC power plant
Wormley, D.N.; Carmichael, D.A.; Umans, S.
1983-08-01
In this study, the aspects of plant performance which influence the potential for integration of an OTEC plant into a utility grid are considered. A set of simulation models have been developed for the evaluation of OTEC dynamic plant performance. A detailed nonlinear dynamic model has been forumlated which is useful for the assessment of component performance including heat exchangers, turbines, pumps and control systems. A reduced order linear model has been developed which is useful for studies of plant stability, control system development and transient performance of the plant connected to a utility grid. This model is particularly suitable for transient dynamic studies of an OTEC plant as a unit in a utility grid. A quasi-steady power availability model has also been developed which is useful to determine plant ouput power as a function of ocean thermal gradients so that the influence of daily and seasonal temperature variations may be easily computed. The study has found no fundamental technical barriers which would prohibit the interconnection of an OTEC plant into a utility grid. It has also shown that detailed consideration of turbine nozzle angle control is merited and such a control has the potential to provide superior performance in comparison to turbine bypass valve control.
Mathematical Modeling of Microbial Community Dynamics: A Methodological Review
Song, Hyun-Seob; Cannon, William R.; Beliaev, Alex S.; Konopka, Allan
2014-10-17
Microorganisms in nature form diverse communities that dynamically change in structure and function in response to environmental variations. As a complex adaptive system, microbial communities show higher-order properties that are not present in individual microbes, but arise from their interactions. Predictive mathematical models not only help to understand the underlying principles of the dynamics and emergent properties of natural and synthetic microbial communities, but also provide key knowledge required for engineering them. In this article, we provide an overview of mathematical tools that include not only current mainstream approaches, but also less traditional approaches that, in our opinion, can be potentially useful. We discuss a broad range of methods ranging from low-resolution supra-organismal to high-resolution individual-based modeling. Particularly, we highlight the integrative approaches that synergistically combine disparate methods. In conclusion, we provide our outlook for the key aspects that should be further developed to move microbial community modeling towards greater predictive power.
Melintescu, A.; Galeriu, D.; Diabate, S.; Strack, S.
2015-03-15
The processes involved in tritium transfer in crops are complex and regulated by many feedback mechanisms. A full mechanistic model is difficult to develop due to the complexity of the processes involved in tritium transfer and environmental conditions. First, a review of existing models (ORYZA2000, CROPTRIT and WOFOST) presenting their features and limits, is made. Secondly, the preparatory steps for a robust model are discussed, considering the role of dry matter and photosynthesis contribution to the OBT (Organically Bound Tritium) dynamics in crops.
Greene, David L
2007-02-01
The global energy system faces sweeping changes in the next few decades, with potentially critical implications for the global economy and the global environment. It is important that global institutions have the tools necessary to predict, analyze and plan for such massive change. This report summarizes the proceedings of an international workshop concerning methods of forecasting, analyzing, and planning for global energy transitions and their economic and environmental consequences. A specific case, it focused on the transition from conventional to unconventional oil and other energy sources likely to result from a peak in non-OPEC and/or global production of conventional oil. Leading energy models from around the world in government, academia and the private sector met, reviewed the state-of-the-art of global energy modeling and evaluated its ability to analyze and predict large-scale energy transitions.
Office Of Nuclear Energy Annual Review Meeting Dynamic Simulation Modeling Tool
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
Annual Review Meeting Dynamic Simulation Modeling Tool Lou Qualls ORNL September 16-18, 2014 2 Work Package SR-14OR130108 - Modeling Tools for Dynamic Behavior Simulations of SMRs 2 ü FY14 molten salt cooled model deliverable due. n FY15 web application deliverable due. n FY15 model repository establishment due. n FY15 working collaboration with University partners. n Simplified Dynamic Modeling for Advanced SMRs - Numerous dynamic models are needed to simulate plant behavior
Model of fracture of metal melts and the strength of melts under dynamic conditions
Mayer, P. N. Mayer, A. E.
2015-07-15
The development of a continuum model of deformation and fracture of melts is needed for the description of the behavior of metals in extreme states, in particular, under high-current electron and ultrashort laser irradiation. The model proposed includes the equations of mechanics of a two-phase continuum and the equations of the kinetics of phase transitions. The change (exchange) of the volumes of dispersed and carrier phases and of the number of dispersed particles is described, and the energy and mass exchange between the phases due to phase transitions is taken into account. Molecular dynamic (MD) calculations are carried out with the use of the LAMMPS program. The continuum model is verified by MD, computational, and experimental data. The strength of aluminum, copper, and nickel is determined at various temperatures and strain rates. It is shown that an increase in the strain rate leads to an increase in the strength of a liquid metal, while an increase in temperature leads to a decrease in its strength.
Model of a deterministic detector and dynamical decoherence
Lee, Jae Weon; Shepelyansky, Dima L. [Laboratoire de Physique Theorique, UMR 5152 du CNRS, Univ. P. Sabatier, 31062 Toulouse Cedex 4 (France); Averin, Dmitri V. [Department of Physics, University of Stony Brook, SUNY, Stony Brook, New York 11794 (United States); Benenti, Giuliano [Center for Nonlinear and Complex Systems, Universita degli Studi dell'Insubria and Istituto Nazionale per la Fisica della Materia, Unita di Como, Via Valleggio 11, 22100 Como (Italy)
2005-07-15
We discuss a deterministic model of detector coupled to a two-level system (a qubit). The detector is a quasiclassical object whose dynamics is described by the kicked rotator Hamiltonian. We show that in the regime of quantum chaos the detector acts as a chaotic bath and induces decoherence of the qubit. We discuss the dephasing and relaxation rates and demonstrate that the main features of single-qubit decoherence due to a heat bath can be reproduced by our fully deterministic dynamical model. Moreover, we show that, for strong enough qubit-detector coupling, the dephasing rate is given by the rate of exponential instability of the detector's dynamics, that is, by the Lyapunov exponent of classical motion. Finally, we discuss the measurement in the regimes of strong and weak qubit-detector coupling. For the case of strong coupling the detector performs a measurement of the up/down state of the qubit. In the case of weak coupling, due to chaos, the dynamical evolution of the detector is strongly sensitive to the state of the qubit. However, in this case it is unclear how to extract a signal from any measurement with a coarse-graining in the phase space on a size much larger than the Planck cell.
Modeling Crabbing Dynamics in an Electron-Ion Collider
Castilla, Alejandro; Morozov, Vasiliy S.; Satogata, Todd J.; Delayen, Jean R.
2015-09-01
A local crabbing scheme requires ?/2 (mod ?) horizontal betatron phase advances from an interaction point (IP) to the crab cavities on each side of it. However, realistic phase advances generated by sets of quadrupoles, or Final Focusing Blocks (FFB), between the crab cavities located in the expanded beam regions and the IP differ slightly from ?/2. To understand the effect of crabbing on the beam dynamics in this case, a simple model of the optics of the Medium Energy Electron-Ion Collider (MEIC) including local crabbing was developed using linear matrices and then studied numerically over multiple turns (1000 passes) of both electron and proton bunches. The same model was applied to both local and global crabbing schemes to determine the linear-order dynamical effects of the synchro-betatron coupling induced by crabbing.
Modeling the dynamic crush of impact mitigating materials
Logan, R.W.; McMichael, L.D.
1995-05-12
Crushable materials are commonly utilized in the design of structural components to absorb energy and mitigate shock during the dynamic impact of a complex structure, such as an automobile chassis or drum-type shipping container. The development and application of several finite-element material models which have been developed at various times at LLNL for DYNA3D will be discussed. Between the models, they are able to account for several of the predominant mechanisms which typically influence the dynamic mechanical behavior of crushable materials. One issue we addressed was that no single existing model would account for the entire gambit of constitutive features which are important for crushable materials. Thus, we describe the implementation and use of an additional material model which attempts to provide a more comprehensive model of the mechanics of crushable material behavior. This model combines features of the pre-existing DYNA models and incorporates some new features as well in an invariant large-strain formulation. In addition to examining the behavior of a unit cell in uniaxial compression, two cases were chosen to evaluate the capabilities and accuracy of the various material models in DYNA. In the first case, a model for foam filled box beams was developed and compared to test data from a 4-point bend test. The model was subsequently used to study its effectiveness in energy absorption in an aluminum extrusion, spaceframe, vehicle chassis. The second case examined the response of the AT-400A shipping container and the performance of the overpack material during accident environments selected from 10CFR71 and IAEA regulations.
Erdmann, Thorsten; Albert, Philipp J.; Schwarz, Ulrich S.
2013-11-07
Non-processive molecular motors have to work together in ensembles in order to generate appreciable levels of force or movement. In skeletal muscle, for example, hundreds of myosin II molecules cooperate in thick filaments. In non-muscle cells, by contrast, small groups with few tens of non-muscle myosin II motors contribute to essential cellular processes such as transport, shape changes, or mechanosensing. Here we introduce a detailed and analytically tractable model for this important situation. Using a three-state crossbridge model for the myosin II motor cycle and exploiting the assumptions of fast power stroke kinetics and equal load sharing between motors in equivalent states, we reduce the stochastic reaction network to a one-step master equation for the binding and unbinding dynamics (parallel cluster model) and derive the rules for ensemble movement. We find that for constant external load, ensemble dynamics is strongly shaped by the catch bond character of myosin II, which leads to an increase of the fraction of bound motors under load and thus to firm attachment even for small ensembles. This adaptation to load results in a concave force-velocity relation described by a Hill relation. For external load provided by a linear spring, myosin II ensembles dynamically adjust themselves towards an isometric state with constant average position and load. The dynamics of the ensembles is now determined mainly by the distribution of motors over the different kinds of bound states. For increasing stiffness of the external spring, there is a sharp transition beyond which myosin II can no longer perform the power stroke. Slow unbinding from the pre-power-stroke state protects the ensembles against detachment.
Nuclear Hybrid Energy System Modeling: RELAP5 Dynamic Coupling Capabilities
Piyush Sabharwall; Nolan Anderson; Haihua Zhao; Shannon Bragg-Sitton; George Mesina
2012-09-01
The nuclear hybrid energy systems (NHES) research team is currently developing a dynamic simulation of an integrated hybrid energy system. A detailed simulation of proposed NHES architectures will allow initial computational demonstration of a tightly coupled NHES to identify key reactor subsystem requirements, identify candidate reactor technologies for a hybrid system, and identify key challenges to operation of the coupled system. This work will provide a baseline for later coupling of design-specific reactor models through industry collaboration. The modeling capability addressed in this report focuses on the reactor subsystem simulation.
An opinion-driven behavioral dynamics model for addictive behaviors
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Moore, Thomas W.; Finley, Patrick D.; Apelberg, Benjamin J.; Ambrose, Bridget K.; Brodsky, Nancy S.; Brown, Theresa J.; Husten, Corinne; Glass, Robert J.
2015-04-08
We present a model of behavioral dynamics that combines a social network-based opinion dynamics model with behavioral mapping. The behavioral component is discrete and history-dependent to represent situations in which an individual’s behavior is initially driven by opinion and later constrained by physiological or psychological conditions that serve to maintain the behavior. Additionally, individuals are modeled as nodes in a social network connected by directed edges. Parameter sweeps illustrate model behavior and the effects of individual parameters and parameter interactions on model results. Mapping a continuous opinion variable into a discrete behavioral space induces clustering on directed networks. Clusters providemore » targets of opportunity for influencing the network state; however, the smaller the network the greater the stochasticity and potential variability in outcomes. Furthermore, this has implications both for behaviors that are influenced by close relationships verses those influenced by societal norms and for the effectiveness of strategies for influencing those behaviors.« less
Use dynamic simulation to model HPU reactor depressuring
Ernest, J.B.; Depew, C.A. )
1995-01-01
Dynamic simulation is the best available method for the analysis of hydroprocessing unit (HPU) depressuring. Depressuring is crucial for the safe operation of hydrocracking and other HPUs with catalysts that have hydrocracking activity. Effective design for depressuring is valuable for all types of HPUs, both grass-roots and revamps. Reactor loop depressuring can set design temperatures and pressures for the reactor effluent cooling train and other equipment and piping in an HPU. Unfortunately, usual methods for determining equipment and piping design conditions during depressuring leave much room for improvement because they poorly account for time-dependent temperature and pressure changes. Dynamic simulation makes it practical to more accurately estimate these transient conditions. The paper discusses depressuring design, including the nature of depressuring, the impact of depressuring on design, and depressuring calculation methods. The author then describes modeling of hydroprocessing unit depressuring by discussing the general and particular correspondence of simulation modules to physical equipment using the base case of total electrical power failure. The special data that is required for dynamic simulation is described and typical simulation results are given. Lastly, the advantages of dynamic simulation are summarized.
Colvin, J D; Minich, R W; Kalantar, D H
2007-03-29
The recent diagnostic capability of the Omega laser to study solid-solid phase transitions at pressures greater than 10 GPa and at strain rates exceeding 10{sup 7} s{sup -1} has also provided valuable information on the dynamic elastic-plastic behavior of materials. We have found, for example, that plasticity kinetics modifies the effective loading and thermodynamic paths of the material. In this paper we derive a kinetics equation for the time-dependent plastic response of the material to dynamic loading, and describe the model's implementation in a radiation-hydrodynamics computer code. This model for plasticity kinetics incorporates the Gilman model for dislocation multiplication and saturation. We discuss the application of this model to the simulation of experimental velocity interferometry data for experiments on Omega in which Fe was shock compressed to pressures beyond the {alpha}-to-{var_epsilon} phase transition pressure. The kinetics model is shown to fit the data reasonably well in this high strain rate regime and further allows quantification of the relative contributions of dislocation multiplication and drag. The sensitivity of the observed signatures to the kinetics model parameters is presented.
Eulerian hydrocode modeling of a dynamic tensile extrusion experiment (u)
Burkett, Michael W; Clancy, Sean P
2009-01-01
Eulerian hydrocode simulations utilizing the Mechanical Threshold Stress flow stress model were performed to provide insight into a dynamic extrusion experiment. The dynamic extrusion response of copper (three different grain sizes) and tantalum spheres were simulated with MESA, an explicit, 2-D Eulerian continuum mechanics hydrocode and compared with experimental data. The experimental data consisted of high-speed images of the extrusion process, recovered extruded samples, and post test metallography. The hydrocode was developed to predict large-strain and high-strain-rate loading problems. Some of the features of the features of MESA include a high-order advection algorithm, a material interface tracking scheme and a van Leer monotonic advection-limiting. The Mechanical Threshold Stress (MTS) model was utilized to evolve the flow stress as a function of strain, strain rate and temperature for copper and tantalum. Plastic strains exceeding 300% were predicted in the extrusion of copper at 400 m/s, while plastic strains exceeding 800% were predicted for Ta. Quantitative comparisons between the predicted and measured deformation topologies and extrusion rate were made. Additionally, predictions of the texture evolution (based upon the deformation rate history and the rigid body rotations experienced by the copper during the extrusion process) were compared with the orientation imaging microscopy measurements. Finally, comparisons between the calculated and measured influence of the initial texture on the dynamic extrusion response of tantalum was performed.
Dynamic mesoscale model of dipolar fluids via fluctuating hydrodynamics
Persson, Rasmus A. X.; Chu, Jhih-Wei, E-mail: jwchu@nctu.edu.tw [Institute of Bioinformatics and Systems Biology, National Chiao Tung University, Hsinchu 30068, Taiwan (China); Department of Biological Science and Technology, National Chiao Tung University, Hsinchu 30068, Taiwan (China); Voulgarakis, Nikolaos K. [Department of Mathematics, Washington State University, Richland, Washington 99372 (United States)
2014-11-07
Fluctuating hydrodynamics (FHD) is a general framework of mesoscopic modeling and simulation based on conservational laws and constitutive equations of linear and nonlinear responses. However, explicit representation of electrical forces in FHD has yet to appear. In this work, we devised an Ansatz for the dynamics of dipole moment densities that is linked with the Poisson equation of the electrical potential ? in coupling to the other equations of FHD. The resulting ?-FHD equations then serve as a platform for integrating the essential forces, including electrostatics in addition to hydrodynamics, pressure-volume equation of state, surface tension, and solvent-particle interactions that govern the emergent behaviors of molecular systems at an intermediate scale. This unique merit of ?-FHD is illustrated by showing that the water dielectric function and ion hydration free energies in homogeneous and heterogenous systems can be captured accurately via the mesoscopic simulation. Furthermore, we show that the field variables of ?-FHD can be mapped from the trajectory of an all-atom molecular dynamics simulation such that model development and parametrization can be based on the information obtained at a finer-grained scale. With the aforementioned multiscale capabilities and a spatial resolution as high as 5 , the ?-FHD equations represent a useful semi-explicit solvent model for the modeling and simulation of complex systems, such as biomolecular machines and nanofluidics.
Computational fluid dynamic modeling of fluidized-bed polymerization reactors
Rokkam, Ram
2012-01-01
Polyethylene is one of the most widely used plastics, and over 60 million tons are produced worldwide every year. Polyethylene is obtained by the catalytic polymerization of ethylene in gas and liquid phase reactors. The gas phase processes are more advantageous, and use fluidized-bed reactors for production of polyethylene. Since they operate so close to the melting point of the polymer, agglomeration is an operational concern in all slurry and gas polymerization processes. Electrostatics and hot spot formation are the main factors that contribute to agglomeration in gas-phase processes. Electrostatic charges in gas phase polymerization fluidized bed reactors are known to influence the bed hydrodynamics, particle elutriation, bubble size, bubble shape etc. Accumulation of electrostatic charges in the fluidized-bed can lead to operational issues. In this work a first-principles electrostatic model is developed and coupled with a multi-fluid computational fluid dynamic (CFD) model to understand the effect of electrostatics on the dynamics of a fluidized-bed. The multi-fluid CFD model for gas-particle flow is based on the kinetic theory of granular flows closures. The electrostatic model is developed based on a fixed, size-dependent charge for each type of particle (catalyst, polymer, polymer fines) phase. The combined CFD model is first verified using simple test cases, validated with experiments and applied to a pilot-scale polymerization fluidized-bed reactor. The CFD model reproduced qualitative trends in particle segregation and entrainment due to electrostatic charges observed in experiments. For the scale up of fluidized bed reactor, filtered models are developed and implemented on pilot scale reactor.
Impact of high-order moments on the statistical modeling of transition arrays
Gilleron, Franck; Pain, Jean-Christophe; Bauche, Jacques; Bauche-Arnoult, Claire
2008-02-15
The impact of high-order moments on the statistical modeling of transition arrays in complex spectra is studied. It is shown that a departure from the Gaussian, which is usually employed in such an approach, may be observed even in the shape of unresolved spectra due to the large value of the kurtosis coefficient. The use of a Gaussian shape may also overestimate the width of the spectra in some cases. Therefore, it is proposed to simulate the statistical shape of the transition arrays by the more flexible generalized Gaussian distribution which introduces an additional parameter--the power of the argument in the exponential--that can be constrained by the kurtosis value. The relevance of the statistical line distribution is checked by comparisons with smoothed spectra obtained from detailed line-by-line calculations. The departure from the Gaussian is also confirmed through the analysis of 2p-3d transitions of recent absorption measurements. A numerical fit is proposed for an easy implementation of the statistical profile in atomic-structure codes.
A human breast cell model of pre-invasive to invasive transition
Bissell, Mina J; Rizki, Aylin; Weaver, Valerie M.; Lee, Sun-Young; Rozenberg, Gabriela I.; Chin, Koei; Myers, Connie A.; Bascom, Jamie L.; Mott, Joni D.; Semeiks, Jeremy R.; Grate, Leslie R.; Mian, I. Saira; Borowsky, Alexander D.; Jensen, Roy A.; Idowu, Michael O.; Chen, Fanqing; Chen, David J.; Petersen, Ole W.; Gray, Joe W.; Bissell, Mina J.
2008-03-10
A crucial step in human breast cancer progression is the acquisition of invasiveness. There is a distinct lack of human cell culture models to study the transition from pre-invasive to invasive phenotype as it may occur 'spontaneously' in vivo. To delineate molecular alterations important for this transition, we isolated human breast epithelial cell lines that showed partial loss of tissue polarity in three-dimensional reconstituted-basement membrane cultures. These cells remained non-invasive; however, unlike their non-malignant counterparts, they exhibited a high propensity to acquire invasiveness through basement membrane in culture. The genomic aberrations and gene expression profiles of the cells in this model showed a high degree of similarity to primary breast tumor profiles. The xenograft tumors formed by the cell lines in three different microenvironments in nude mice displayed metaplastic phenotypes, including squamous and basal characteristics, with invasive cells exhibiting features of higher grade tumors. To find functionally significant changes in transition from pre-invasive to invasive phenotype, we performed attribute profile clustering analysis on the list of genes differentially expressed between pre-invasive and invasive cells. We found integral membrane proteins, transcription factors, kinases, transport molecules, and chemokines to be highly represented. In addition, expression of matrix metalloproteinases MMP-9,-13,-15,-17 was up regulated in the invasive cells. Using siRNA based approaches, we found these MMPs to be required for the invasive phenotype. This model provides a new tool for dissection of mechanisms by which pre-invasive breast cells could acquire invasiveness in a metaplastic context.
A self-consistent two-fluid model of a magnetized plasma-wall transition
Gyergyek, T.; Kovačič, J.
2015-09-15
A self-consistent one-dimensional two-fluid model of the magnetized plasma-wall transition is presented. The model includes magnetic field, elastic collisions between ions and electrons, and creation/annihilation of charged particles. Two systems of differential equations are derived. The first system describes the whole magnetized plasma-wall transition region, which consists of the pre-sheath, the magnetized pre-sheath (Chodura layer), and the sheath, which is not neutral, but contains a positive space charge. The second system of equations describes only the neutral part of the plasma-wall transition region—this means only the pre-sheath and the Chodura layer, but not also the sheath. Both systems are solved numerically. The first system of equations has two singularities. The first occurs when ion velocity in the direction perpendicularly to the wall drops below the ion thermal velocity. The second occurs when the electron velocity in the direction perpendicularly to the wall exceeds the electron thermal velocity. The second system of differential equations only has one singularity, which has also been derived analytically. For finite electron to ion mass ratio, the integration of the second system always breaks down before the Bohm criterion is fulfilled. Some properties of the first system of equations are examined. It is shown that the increased collision frequency demagnetizes the plasma. On the other hand, if the magnetic field is so strong that the ion Larmor radius and the Debye length are comparable, the electron velocity in the direction perpendicularly to the wall reaches the electron thermal velocity before the ion velocity in the direction perpendicularly to the wall reaches the ion sound velocity. In this case, the integration of the model equations breaks down before the Bohm criterion is fulfilled and the sheath is formed.
Dark matter and strong electroweak phase transition in a radiative neutrino mass model
Ahriche, Amine; Nasri, Salah E-mail: snasri@uaeu.ac.ae
2013-07-01
We consider an extension of the standard model (SM) with charged singlet scalars and right handed (RH) neutrinos all at the electroweak scale. In this model, the neutrino masses are generated at three loops, which provide an explanation for their smallness, and the lightest RH neutrino, N{sub 1}, is a dark matter candidate. We find that for three generations of RH neutrinos, the model can be consistent with the neutrino oscillation data, lepton flavor violating processes, N{sub 1} can have a relic density in agreement with the recent Planck data, and the electroweak phase transition can be strongly first order. We also show that the charged scalars may enhance the branching ratio h→γγ, where as h→γZ get can get few percent suppression. We also discuss the phenomenological implications of the RH neutrinos at the collider.
Block entropy and quantum phase transition in the anisotropic Kondo necklace model
Mendoza-Arenas, J. J.; Franco, R.; Silva-Valencia, J.
2010-06-15
We study the von Neumann block entropy in the Kondo necklace model for different anisotropies {eta} in the XY interaction between conduction spins using the density matrix renormalization group method. It was found that the block entropy presents a maximum for each {eta} considered, and, comparing it with the results of the quantum criticality of the model based on the behavior of the energy gap, we observe that the maximum block entropy occurs at the quantum critical point between an antiferromagnetic and a Kondo singlet state, so this measure of entanglement is useful for giving information about where a quantum phase transition occurs in this model. We observe that the block entropy also presents a maximum at the quantum critical points that are obtained when an anisotropy {Delta} is included in the Kondo exchange between localized and conduction spins; when {Delta} diminishes for a fixed value of {eta}, the critical point increases, favoring the antiferromagnetic phase.
Best practices for system dynamics model design and construction with powersim studio.
Malczynski, Leonard A.
2011-06-01
This guide addresses software quality in the construction of Powersim{reg_sign} Studio 8 system dynamics simulation models. It is the result of almost ten years of experience with the Powersim suite of system dynamics modeling tools (Constructor and earlier Studio versions). It is a guide that proposes a common look and feel for the construction of Powersim Studio system dynamics models.
Modeling nanoscale hydrodynamics by smoothed dissipative particle dynamics
Lei, Huan; Mundy, Christopher J.; Schenter, Gregory K.; Voulgarakis, Nikolaos
2015-05-21
Thermal fluctuation and hydrophobicity are two hallmarks of fluid hydrodynamics on the nano-scale. It is a challenge to consistently couple the small length and time scale phenomena associated with molecular interaction with larger scale phenomena. The development of this consistency is the essence of mesoscale science. In this study, we develop a nanoscale fluid model based on smoothed dissipative particle dynamics that accounts for the phenomena of associated with density fluctuations and hydrophobicity. We show consistency in the fluctuation spectrum across scales. In doing so, it is necessary to account for finite fluid particle size. Furthermore, we demonstrate that the present model can capture of the void probability and solvation free energy of apolar particles of different sizes. The present fluid model is well suited for a understanding emergent phenomena in nano-scale fluid systems.
Aggregated Residential Load Modeling Using Dynamic Bayesian Networks
Vlachopoulou, Maria; Chin, George; Fuller, Jason C.; Lu, Shuai
2014-09-28
AbstractIt is already obvious that the future power grid will have to address higher demand for power and energy, and to incorporate renewable resources of different energy generation patterns. Demand response (DR) schemes could successfully be used to manage and balance power supply and demand under operating conditions of the future power grid. To achieve that, more advanced tools for DR management of operations and planning are necessary that can estimate the available capacity from DR resources. In this research, a Dynamic Bayesian Network (DBN) is derived, trained, and tested that can model aggregated load of Heating, Ventilation, and Air Conditioning (HVAC) systems. DBNs can provide flexible and powerful tools for both operations and planing, due to their unique analytical capabilities. The DBN model accuracy and flexibility of use is demonstrated by testing the model under different operational scenarios.
Critical dynamics of cluster algorithms in the dilute Ising model
Hennecke, M. Universitaet Karlsruhe ); Heyken, U. )
1993-08-01
Autocorrelation times for thermodynamic quantities at [Tc] are calculated from Monte Carlo simulations of the site-diluted simple cubic Ising model, using the Swendsen-Wand and Wolff cluster algorithms. The results show that for these algorithms the autocorrelation times decrease when reducing the concentration of magnetic sites from 100% down to 40%. This is of crucial importance when estimating static properties of the model, since the variances of these estimators increase with autocorrelation time. The dynamical critical exponents are calculated for both algorithms, observing pronounced finite-size effects in the energy autocorrelation data for the algorithm of Wolff. It is concluded that, when applied to the dilute Ising model, cluster algorithms become even more effective than local algorithms, for which increasing autocorrelation times are expected. 33 refs., 5 figs., 2 tabs.
Vugrin, Eric D.; Rostron, Brian L.; Verzi, Stephen J.; Brodsky, Nancy S.; Brown, Theresa J.; Choiniere, Conrad J.; Coleman, Blair N.; Paredes, Antonio; Apelberg, Benjamin J.
2015-03-27
Background Recent declines in US cigarette smoking prevalence have coincided with increases in use of other tobacco products. Multiple product tobacco models can help assess the population health impacts associated with use of a wide range of tobacco products. Methods and Findings We present a multi-state, dynamical systems population structure model that can be used to assess the effects of tobacco product use behaviors on population health. The model incorporates transition behaviors, such as initiation, cessation, switching, and dual use, related to the use of multiple products. The model tracks product use prevalence and mortality attributable to tobacco use for the overall population and by sex and age group. The model can also be used to estimate differences in these outcomes between scenarios by varying input parameter values. We demonstrate model capabilities by projecting future cigarette smoking prevalence and smoking-attributable mortality and then simulating the effects of introduction of a hypothetical new lower-risk tobacco product under a variety of assumptions about product use. Sensitivity analyses were conducted to examine the range of population impacts that could occur due to differences in input values for product use and risk. We demonstrate that potential benefits from cigarette smokers switching to the lower-risk product can be offset over time through increased initiation of this product. Model results show that population health benefits are particularly sensitive to product risks and initiation, switching, and dual use behaviors. Conclusion Our model incorporates the variety of tobacco use behaviors and risks that occur with multiple products. As such, it can evaluate the population health impacts associated with the introduction of new tobacco products or policies that may result in product switching or dual use. Further model development will include refinement of data inputs for non-cigarette tobacco products and inclusion of health
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Vugrin, Eric D.; Rostron, Brian L.; Verzi, Stephen J.; Brodsky, Nancy S.; Brown, Theresa J.; Choiniere, Conrad J.; Coleman, Blair N.; Paredes, Antonio; Apelberg, Benjamin J.
2015-03-27
Background Recent declines in US cigarette smoking prevalence have coincided with increases in use of other tobacco products. Multiple product tobacco models can help assess the population health impacts associated with use of a wide range of tobacco products. Methods and Findings We present a multi-state, dynamical systems population structure model that can be used to assess the effects of tobacco product use behaviors on population health. The model incorporates transition behaviors, such as initiation, cessation, switching, and dual use, related to the use of multiple products. The model tracks product use prevalence and mortality attributable to tobacco use formore » the overall population and by sex and age group. The model can also be used to estimate differences in these outcomes between scenarios by varying input parameter values. We demonstrate model capabilities by projecting future cigarette smoking prevalence and smoking-attributable mortality and then simulating the effects of introduction of a hypothetical new lower-risk tobacco product under a variety of assumptions about product use. Sensitivity analyses were conducted to examine the range of population impacts that could occur due to differences in input values for product use and risk. We demonstrate that potential benefits from cigarette smokers switching to the lower-risk product can be offset over time through increased initiation of this product. Model results show that population health benefits are particularly sensitive to product risks and initiation, switching, and dual use behaviors. Conclusion Our model incorporates the variety of tobacco use behaviors and risks that occur with multiple products. As such, it can evaluate the population health impacts associated with the introduction of new tobacco products or policies that may result in product switching or dual use. Further model development will include refinement of data inputs for non-cigarette tobacco products and inclusion of
Vugrin, Eric D.; Rostron, Brian L.; Verzi, Stephen J.; Brodsky, Nancy S.; Brown, Theresa J.; Choiniere, Conrad J.; Coleman, Blair N.; Paredes, Antonio; Apelberg, Benjamin J.
2015-03-27
Background Recent declines in US cigarette smoking prevalence have coincided with increases in use of other tobacco products. Multiple product tobacco models can help assess the population health impacts associated with use of a wide range of tobacco products. Methods and Findings We present a multi-state, dynamical systems population structure model that can be used to assess the effects of tobacco product use behaviors on population health. The model incorporates transition behaviors, such as initiation, cessation, switching, and dual use, related to the use of multiple products. The model tracks product use prevalence and mortality attributable to tobacco use for the overall population and by sex and age group. The model can also be used to estimate differences in these outcomes between scenarios by varying input parameter values. We demonstrate model capabilities by projecting future cigarette smoking prevalence and smoking-attributable mortality and then simulating the effects of introduction of a hypothetical new lower-risk tobacco product under a variety of assumptions about product use. Sensitivity analyses were conducted to examine the range of population impacts that could occur due to differences in input values for product use and risk. We demonstrate that potential benefits from cigarette smokers switching to the lower-risk product can be offset over time through increased initiation of this product. Model results show that population health benefits are particularly sensitive to product risks and initiation, switching, and dual use behaviors. Conclusion Our model incorporates the variety of tobacco use behaviors and risks that occur with multiple products. As such, it can evaluate the population health impacts associated with the introduction of new tobacco products or policies that may result in product switching or dual use. Further model development will include refinement of data inputs for non-cigarette tobacco products and inclusion of health
A Dynamic Approach to Modeling Dependence Between Human Failure Events
Boring, Ronald Laurids
2015-09-01
In practice, most HRA methods use direct dependence from THERP—the notion that error be- gets error, and one human failure event (HFE) may increase the likelihood of subsequent HFEs. In this paper, we approach dependence from a simulation perspective in which the effects of human errors are dynamically modeled. There are three key concepts that play into this modeling: (1) Errors are driven by performance shaping factors (PSFs). In this context, the error propagation is not a result of the presence of an HFE yielding overall increases in subsequent HFEs. Rather, it is shared PSFs that cause dependence. (2) PSFs have qualities of lag and latency. These two qualities are not currently considered in HRA methods that use PSFs. Yet, to model the effects of PSFs, it is not simply a matter of identifying the discrete effects of a particular PSF on performance. The effects of PSFs must be considered temporally, as the PSFs will have a range of effects across the event sequence. (3) Finally, there is the concept of error spilling. When PSFs are activated, they not only have temporal effects but also lateral effects on other PSFs, leading to emergent errors. This paper presents the framework for tying together these dynamic dependence concepts.
User Guide for PV Dynamic Model Simulation Written on PSCAD Platform
Muljadi, E.; Singh, M.; Gevorgian, V.
2014-11-01
This document describes the dynamic photovoltaic model developed by the National Renewable Energy Laboratory and is intended as a guide for users of these models.
THE HIGHLY DYNAMIC BEHAVIOR OF THE INNERMOST DUST AND GAS IN THE TRANSITION DISK VARIABLE LRLL 31
Flaherty, K. M.; Rieke, G.; Muzerolle, J.; Gutermuth, R.; Balog, Z.; Herbst, W.; Megeath, S. T.; Kun, M.
2011-05-10
We describe extensive synoptic multi-wavelength observations of the transition disk LRLL 31 in the young cluster IC 348. We combined 4 epochs of IRS spectra, 9 epochs of MIPS photometry, 7 epochs of cold-mission IRAC photometry, and 36 epochs of warm-mission IRAC photometry along with multi-epoch near-infrared spectra, optical spectra, and polarimetry to explore the nature of the rapid variability of this object. We find that the inner disk, as traced by the 2-5 {mu}m excess, stays at the dust sublimation radius while the strength of the excess changes by a factor of eight on weekly timescales, and the 3.6 and 4.5 {mu}m photometry show a drop of 0.35 mag in 1 week followed by a slow 0.5 mag increase over the next 3 weeks. The accretion rate, as measured by Pa{beta} and Br{gamma} emission lines, varies by a factor of five with evidence for a correlation between the accretion rate and the infrared excess. While the gas and dust in the inner disk are fluctuating, the central star stays relatively static. Our observations allow us to put constraints on the physical mechanism responsible for the variability. The variable accretion, and wind, are unlikely to be causes of the variability, but are both effects of the same physical process that disturbs the disk. The lack of periodicity in our infrared monitoring indicates that it is unlikely that there is a companion within {approx}0.4 AU that is perturbing the disk. The most likely explanation is either a companion beyond {approx}0.4 AU or a dynamic interface between the stellar magnetic field and the disk leading to a variable scale height and/or warping of the inner disk.
A criterion for the dynamical to kinematical transition of x-ray diffraction on a bent crystal
Kushnir, V.I.; Macrander, A.T.
1993-09-01
It is well known that the peak reflectivity of a bent crystal, generally speaking, is smaller than that of a plane crystal, and it goes to zero when the crystal curvature goes to infinity. The reason for this is the transition between dynamical and kinematical diffraction that takes place as the crystal curvature increases. The physical explanation is as follows: the deviation from exact Bragg position along the beam changes so fast that the thickness over which the beam is within a Darwin width becomes too small to reflect the beam. Bent crystals are widely used as focusing elements in X-ray optics, and estimation of whether or not a bent crystal is still perfect enough to provide good reflectivity is of great importance. Currently the Advanced Photon Source (APS) is considering a number of bent crystals as focusing elements for future APS beamlines, including a sagittaly focusing monochromator and bent backscattering analyzer for inelastic X-ray scattering experiments. A criterion is given in answer to the question: To what extent is it possible to bend a crystal without loss of X-ray peak reflectivity? An expression based on the work of Chukhovskii, Gabrielyan and Petrashen, is formulated that applies to anisotropic cubic crystal and that can be used not only for conventional asymmetric Bragg diffraction, but also for inclined crystal diffraction. The following special cases are treated as examples: isotropic crystal, standard symmetrical Bragg diffraction, extremely asymmetric diffraction, and backscattering with Bragg angles near 90{degree}. In addition, an asymptotic behavior for high energies is detailed.
Model Studies of the Dynamics of Bacterial Flagellar Motors
Bai, F; Lo, C; Berry, R; Xing, J
2009-03-19
The Bacterial Flagellar Motor is a rotary molecular machine that rotates the helical filaments which propel swimming bacteria. Extensive experimental and theoretical studies exist on the structure, assembly, energy input, power generation and switching mechanism of the motor. In our previous paper, we explained the general physics underneath the observed torque-speed curves with a simple two-state Fokker-Planck model. Here we further analyze this model. In this paper we show (1) the model predicts that the two components of the ion motive force can affect the motor dynamics differently, in agreement with the latest experiment by Lo et al.; (2) with explicit consideration of the stator spring, the model also explains the lack of dependence of the zero-load speed on stator number in the proton motor, recently observed by Yuan and Berg; (3) the model reproduces the stepping behavior of the motor even with the existence of the stator springs and predicts the dwelling time distribution. Predicted stepping behavior of motors with two stators is discussed, and we suggest future experimental verification.
Toward a mechanistic modeling of nitrogen limitation on vegetation dynamics
Xu, Chonggang [Los Alamos National Laboratory (LANL); Fisher, Rosie [National Center for Atmospheric Research (NCAR); Wullschleger, Stan D [ORNL; Wilson, Cathy [Los Alamos National Laboratory (LANL); Cai, Michael [Los Alamos National Laboratory (LANL); McDowell, Nathan [Los Alamos National Laboratory (LANL)
2012-01-01
Nitrogen is a dominant regulator of vegetation dynamics, net primary production, and terrestrial carbon cycles; however, most ecosystem models use a rather simplistic relationship between leaf nitrogen content and photosynthetic capacity. Such an approach does not consider how patterns of nitrogen allocation may change with differences in light intensity, growing-season temperature and CO{sub 2} concentration. To account for this known variability in nitrogen-photosynthesis relationships, we develop a mechanistic nitrogen allocation model based on a trade-off of nitrogen allocated between growth and storage, and an optimization of nitrogen allocated among light capture, electron transport, carboxylation, and respiration. The developed model is able to predict the acclimation of photosynthetic capacity to changes in CO{sub 2} concentration, temperature, and radiation when evaluated against published data of V{sub c,max} (maximum carboxylation rate) and J{sub max} (maximum electron transport rate). A sensitivity analysis of the model for herbaceous plants, deciduous and evergreen trees implies that elevated CO{sub 2} concentrations lead to lower allocation of nitrogen to carboxylation but higher allocation to storage. Higher growing-season temperatures cause lower allocation of nitrogen to carboxylation, due to higher nitrogen requirements for light capture pigments and for storage. Lower levels of radiation have a much stronger effect on allocation of nitrogen to carboxylation for herbaceous plants than for trees, resulting from higher nitrogen requirements for light capture for herbaceous plants. As far as we know, this is the first model of complete nitrogen allocation that simultaneously considers nitrogen allocation to light capture, electron transport, carboxylation, respiration and storage, and the responses of each to altered environmental conditions. We expect this model could potentially improve our confidence in simulations of carbon-nitrogen interactions
Lookman, Turab; Vasseur, Romain
2009-01-01
We obtain the microstructure of ferroelastic transitions in two and three dimensions from the solution of their corresponding discrete pseudo-spin models. In two dimensions we consider two transitions each from the high symmetry square and triangle symmetries: square-to-rectangle (SR), square-to-oblique (SO), triangle-to-centered rectangle (TR) and triangle-to-oblique (TO). In three dimensions we study the corresponding spin model for the cubic to tetragonal transition. The Landau free energies for these transitions result in N+ I states clock models (Z{sub N}) with long range interactions and we derive mean-field self-consistency equations for the clock model Hamiltonians. The microstructures from the mean-field solutions of the models are very similar to those obtained from the original continuum models or Monte Carlo simulations on the spin models (in the SR case), illustrating that these discrete models capture the salient physics. The models, in the presence of disorder, provide the basis for the study of the strain glass phase observed in martensitic alloys.
VISION -- A Dynamic Model of the Nuclear Fuel Cycle
J. J. Jacobson; A. M. Yacout; S. J. Piet; D. E. Shropshire; G. E. Matthern
2006-02-01
The Advanced Fuel Cycle Initiative’s (AFCI) fundamental objective is to provide technology options that – if implemented – would enable long-term growth of nuclear power while improving sustainability and energy security. The AFCI organization structure consists of four areas; Systems Analysis, Fuels, Separations and Transmutations. The Systems Analysis Working Group is tasked with bridging the program technical areas and providing the models, tools, and analyses required to assess the feasibility of design and deploy¬ment options and inform key decision makers. An integral part of the Systems Analysis tool set is the development of a system level model that can be used to examine the implications of the different mixes of reactors, implications of fuel reprocessing, impact of deployment technologies, as well as potential “exit” or “off ramp” approaches to phase out technologies, waste management issues and long-term repository needs. The Verifiable Fuel Cycle Simulation Model (VISION) is a computer-based simulation model that allows performing dynamic simulations of fuel cycles to quantify infrastructure requirements and identify key trade-offs between alternatives. VISION is intended to serve as a broad systems analysis and study tool applicable to work conducted as part of the AFCI (including costs estimates) and Generation IV reactor development studies.
Computational fluid dynamics modeling of proton exchange membrane fuel cells
UM,SUKKEE; WANG,C.Y.; CHEN,KEN S.
2000-02-11
A transient, multi-dimensional model has been developed to simulate proton exchange membrane (PEM) fuel cells. The model accounts simultaneously for electrochemical kinetics, current distribution, hydrodynamics and multi-component transport. A single set of conservation equations valid for flow channels, gas-diffusion electrodes, catalyst layers and the membrane region are developed and numerically solved using a finite-volume-based computational fluid dynamics (CFD) technique. The numerical model is validated against published experimental data with good agreement. Subsequently, the model is applied to explore hydrogen dilution effects in the anode feed. The predicted polarization cubes under hydrogen dilution conditions are found to be in qualitative agreement with recent experiments reported in the literature. The detailed two-dimensional electrochemical and flow/transport simulations further reveal that in the presence of hydrogen dilution in the fuel stream, hydrogen is depleted at the reaction surface resulting in substantial kinetic polarization and hence a lower current density that is limited by hydrogen transport from the fuel stream to the reaction site.
Computational Model of Population Dynamics Based on the Cell Cycle and Local Interactions
Oprisan, Sorinel Adrian; Oprisan, Ana
2005-03-31
Our study bridges cellular (mesoscopic) level interactions and global population (macroscopic) dynamics of carcinoma. The morphological differences and transitions between well and smooth defined benign tumors and tentacular malignat tumors suggest a theoretical analysis of tumor invasion based on the development of mathematical models exhibiting bifurcations of spatial patterns in the density of tumor cells. Our computational model views the most representative and clinically relevant features of oncogenesis as a fight between two distinct sub-systems: the immune system of the host and the neoplastic system. We implemented the neoplastic sub-system using a three-stage cell cycle: active, dormant, and necrosis. The second considered sub-system consists of cytotoxic active (effector) cells -- EC, with a very broad phenotype ranging from NK cells to CTL cells, macrophages, etc. Based on extensive numerical simulations, we correlated the fractal dimensions for carcinoma, which could be obtained from tumor imaging, with the malignat stage. Our computational model was able to also simulate the effects of surgical, chemotherapeutical, and radiotherapeutical treatments.
Relativistic, model-independent, multichannel 2→2 transition amplitudes in a finite volume
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Briceno, Raul A.; Hansen, Maxwell T.
2016-07-13
We derive formalism for determining 2 + J → 2 infinite-volume transition amplitudes from finite-volume matrix elements. Specifically, we present a relativistic, model-independent relation between finite-volume matrix elements of external currents and the physically observable infinite-volume matrix elements involving two-particle asymptotic states. The result presented holds for states composed of two scalar bosons. These can be identical or non-identical and, in the latter case, can be either degenerate or non-degenerate. We further accommodate any number of strongly-coupled two-scalar channels. This formalism will, for example, allow future lattice QCD calculations of themore » $$\\rho$$-meson form factor, in which the unstable nature of the $$\\rho$$ is rigorously accommodated. In conclusion, we also discuss how this work will impact future extractions of nuclear parity and hadronic long-range matrix elements from lattice QCD.« less
Fluid dynamic modeling of nano-thermite reactions
Martirosyan, Karen S.; Zyskin, Maxim; Jenkins, Charles M.; Horie, Yasuyuki
2014-03-14
This paper presents a direct numerical method based on gas dynamic equations to predict pressure evolution during the discharge of nanoenergetic materials. The direct numerical method provides for modeling reflections of the shock waves from the reactor walls that generates pressure-time fluctuations. The results of gas pressure prediction are consistent with the experimental evidence and estimates based on the self-similar solution. Artificial viscosity provides sufficient smoothing of shock wave discontinuity for the numerical procedure. The direct numerical method is more computationally demanding and flexible than self-similar solution, in particular it allows study of a shock wave in its early stage of reaction and allows the investigation of “slower” reactions, which may produce weaker shock waves. Moreover, numerical results indicate that peak pressure is not very sensitive to initial density and reaction time, providing that all the material reacts well before the shock wave arrives at the end of the reactor.
Description of waste pretreatment and interfacing systems dynamic simulation model
Garbrick, D.J.; Zimmerman, B.D.
1995-05-01
The Waste Pretreatment and Interfacing Systems Dynamic Simulation Model was created to investigate the required pretreatment facility processing rates for both high level and low level waste so that the vitrification of tank waste can be completed according to the milestones defined in the Tri-Party Agreement (TPA). In order to achieve this objective, the processes upstream and downstream of the pretreatment facilities must also be included. The simulation model starts with retrieval of tank waste and ends with vitrification for both low level and high level wastes. This report describes the results of three simulation cases: one based on suggested average facility processing rates, one with facility rates determined so that approximately 6 new DSTs are required, and one with facility rates determined so that approximately no new DSTs are required. It appears, based on the simulation results, that reasonable facility processing rates can be selected so that no new DSTs are required by the TWRS program. However, this conclusion must be viewed with respect to the modeling assumptions, described in detail in the report. Also included in the report, in an appendix, are results of two sensitivity cases: one with glass plant water recycle steams recycled versus not recycled, and one employing the TPA SST retrieval schedule versus a more uniform SST retrieval schedule. Both recycling and retrieval schedule appear to have a significant impact on overall tank usage.
Doinikov, Alexander A. Bouakaz, Ayache; Sheeran, Paul S.; Dayton, Paul A.
2014-10-15
Purpose: Perfluorocarbon (PFC) microdroplets, called phase-change contrast agents (PCCAs), are a promising tool in ultrasound imaging and therapy. Interest in PCCAs is motivated by the fact that they can be triggered to transition from the liquid state to the gas state by an externally applied acoustic pulse. This property opens up new approaches to applications in ultrasound medicine. Insight into the physics of vaporization of PFC droplets is vital for effective use of PCCAs and for anticipating bioeffects. PCCAs composed of volatile PFCs (with low boiling point) exhibit complex dynamic behavior: after vaporization by a short acoustic pulse, a PFC droplet turns into a vapor bubble which undergoes overexpansion and damped radial oscillation until settling to a final diameter. This behavior has not been well described theoretically so far. The purpose of our study is to develop an improved theoretical model that describes the vaporization dynamics of volatile PFC droplets and to validate this model by comparison with in vitro experimental data. Methods: The derivation of the model is based on applying the mathematical methods of fluid dynamics and thermodynamics to the process of the acoustic vaporization of PFC droplets. The used approach corrects shortcomings of the existing models. The validation of the model is carried out by comparing simulated results with in vitro experimental data acquired by ultrahigh speed video microscopy for octafluoropropane (OFP) and decafluorobutane (DFB) microdroplets of different sizes. Results: The developed theory allows one to simulate the growth of a vapor bubble inside a PFC droplet until the liquid PFC is completely converted into vapor, and the subsequent overexpansion and damped oscillations of the vapor bubble, including the influence of an externally applied acoustic pulse. To evaluate quantitatively the difference between simulated and experimental results, the L2-norm errors were calculated for all cases where the
Creating dynamic equivalent PV circuit models with impedance spectroscopy for arc-fault modeling.
Johnson, Jay Dean; Kuszmaul, Scott S.; Strauch, Jason E.; Schoenwald, David Alan
2011-06-01
Article 690.11 in the 2011 National Electrical Code{reg_sign} (NEC{reg_sign}) requires new photovoltaic (PV) systems on or penetrating a building to include a listed arc fault protection device. Currently there is little experimental or empirical research into the behavior of the arcing frequencies through PV components despite the potential for modules and other PV components to filter or attenuate arcing signatures that could render the arc detector ineffective. To model AC arcing signal propagation along PV strings, the well-studied DC diode models were found to inadequately capture the behavior of high frequency arcing signals. Instead dynamic equivalent circuit models of PV modules were required to describe the impedance for alternating currents in modules. The nonlinearities present in PV cells resulting from irradiance, temperature, frequency, and bias voltage variations make modeling these systems challenging. Linearized dynamic equivalent circuits were created for multiple PV module manufacturers and module technologies. The equivalent resistances and capacitances for the modules were determined using impedance spectroscopy with no bias voltage and no irradiance. The equivalent circuit model was employed to evaluate modules having irradiance conditions that could not be measured directly with the instrumentation. Although there was a wide range of circuit component values, the complex impedance model does not predict filtering of arc fault frequencies in PV strings for any irradiance level. Experimental results with no irradiance agree with the model and show nearly no attenuation for 1 Hz to 100 kHz input frequencies.
The challenges of modelling antibody repertoire dynamics in HIV infection
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Luo, Shishi; Perelson, Alan S.
2015-07-20
Antibody affinity maturation by somatic hypermutation of B-cell immunoglobulin variable region genes has been studied for decades in various model systems using well-defined antigens. While much is known about the molecular details of the process, our understanding of the selective forces that generate affinity maturation are less well developed, particularly in the case of a co-evolving pathogen such as HIV. Despite this gap in understanding, high-throughput antibody sequence data are increasingly being collected to investigate the evolutionary trajectories of antibody lineages in HIV-infected individuals. Here, we review what is known in controlled experimental systems about the mechanisms underlying antibody selectionmore » and compare this to the observed temporal patterns of antibody evolution in HIV infection. In addition, we describe how our current understanding of antibody selection mechanisms leaves questions about antibody dynamics in HIV infection unanswered. Without a mechanistic understanding of antibody selection in the context of a co-evolving viral population, modelling and analysis of antibody sequences in HIV-infected individuals will be limited in their interpretation and predictive ability.« less
Ultrafast Structural Dynamics in Combustion Relevant Model Systems
Weber, Peter M.
2014-03-31
The research project explored the time resolved structural dynamics of important model reaction system using an array of novel methods that were developed specifically for this purpose. They include time resolved electron diffraction, time resolved relativistic electron diffraction, and time resolved Rydberg fingerprint spectroscopy. Toward the end of the funding period, we also developed time-resolved x-ray diffraction, which uses ultrafast x-ray pulses at LCLS. Those experiments are just now blossoming, as the funding period expired. In the following, the time resolved Rydberg Fingerprint Spectroscopy is discussed in some detail, as it has been a very productive method. The binding energy of an electron in a Rydberg state, that is, the energy difference between the Rydberg level and the ground state of the molecular ion, has been found to be a uniquely powerful tool to characterize the molecular structure. To rationalize the structure sensitivity we invoke a picture from electron diffraction: when it passes the molecular ion core, the Rydberg electron experiences a phase shift compared to an electron in a hydrogen atom. This phase shift requires an adjustment of the binding energy of the electron, which is measurable. As in electron diffraction, the phase shift depends on the molecular, geometrical structure, so that a measurement of the electron binding energy can be interpreted as a measurement of the molecules structure. Building on this insight, we have developed a structurally sensitive spectroscopy: the molecule is first elevated to the Rydberg state, and the binding energy is then measured using photoelectron spectroscopy. The molecules structure is read out as the binding energy spectrum. Since the photoionization can be done with ultrafast laser pulses, the technique is inherently capable of a time resolution in the femtosecond regime. For the purpose of identifying the structures of molecules during chemical reactions, and for the analysis of
Meier, Sebastian B. E-mail: wiebke.sarfert@siemens.com; Hartmann, David; Sarfert, Wiebke E-mail: wiebke.sarfert@siemens.com; Winnacker, Albrecht
2014-09-14
Light-emitting electrochemical cells (LECs) have received increasing attention during recent years due to their simple architecture, based on solely air-stabile materials, and ease of manufacture in ambient atmosphere, using solution-based technologies. The LEC's active layer offers semiconducting, luminescent as well as ionic functionality resulting in device physical processes fundamentally different as compared with organic light-emitting diodes. During operation, electrical double layers (EDLs) form at the electrode interfaces as a consequence of ion accumulation and electrochemical doping sets in leading to the in situ development of a light-emitting p-i-n junction. In this paper, we comment on the use of impedance spectroscopy in combination with complex nonlinear squares fitting to derive key information about the latter events in thin-film ionic transition metal complex-based light-emitting electrochemical cells based on the model compound bis-2-phenylpyridine 6-phenyl-2,2´-bipyridine iridium(III) hexafluoridophosphate ([Ir(ppy)₂(pbpy)][PF₆]). At operating voltages below the bandgap potential of the ionic complex used, we obtain the dielectric constant of the active layer, the conductivity of mobile ions, the transference numbers of electrons and ions, and the thickness of the EDLs, whereas the transient thickness of the p-i-n junction is determined at voltages above the bandgap potential. Most importantly, we find that charge transport is dominated by the ions when carrier injection from the electrodes is prohibited, that ion movement is limited by the presence of transverse internal interfaces and that the width of the intrinsic region constitutes almost 60% of the total active layer thickness in steady state at a low operating voltage.
Nomura, K.; Vretenar, D.; Niksic, T.; Otsuka, T.; Shimizu, N.
2011-07-15
Microscopic energy density functionals have become a standard tool for nuclear structure calculations, providing an accurate global description of nuclear ground states and collective excitations. For spectroscopic applications, this framework has to be extended to account for collective correlations related to restoration of symmetries broken by the static mean field, and for fluctuations of collective variables. In this paper, we compare two approaches to five-dimensional quadrupole dynamics: the collective Hamiltonian for quadrupole vibrations and rotations and the interacting boson model (IBM). The two models are compared in a study of the evolution of nonaxial shapes in Pt isotopes. Starting from the binding energy surfaces of {sup 192,194,196}Pt, calculated with a microscopic energy density functional, we analyze the resulting low-energy collective spectra obtained from the collective Hamiltonian, and the corresponding IBM Hamiltonian. The calculated excitation spectra and transition probabilities for the ground-state bands and the {gamma}-vibration bands are compared to the corresponding sequences of experimental states.
Wu, H. L.; Tian, W. D.; Ma, Y. G.; Cai, X. Z.; Chen, J. G.; Fang, D. Q.; Guo, W.; Wang, H. W.
2010-04-15
Dynamical dipole gamma-ray emission in heavy-ion collisions is explored in the framework of the quantum molecular dynamics model. The studies are focused on systems of {sup 40}Ca bombarding {sup 48}Ca and its isotopes at different incident energies and impact parameters. Yields of gamma rays are calculated and the centroid energy and dynamical dipole emission width of the gamma spectra are extracted to investigate the properties of gamma emission. In addition, sensitivities of dynamical dipole gamma-ray emission to the isospin and the symmetry energy coefficient of the equation of state are studied. The results show that detailed study of dynamical dipole gamma radiation can provide information on the equation of state and the symmetry energy around the normal nuclear density.
Dynamics Modeling and Loads Analysis of an Offshore Floating Wind Turbine
Jonkman, J. M.
2007-12-01
This report describes the development, verification, and application of a comprehensive simulation tool for modeling coupled dynamic responses of offshore floating wind turbines.
THE TRANSITION REGION RESPONSE TO A CORONAL NANOFLARE: FORWARD MODELING AND OBSERVATIONS IN SDO/AIA
Viall, Nicholeen M.; Klimchuk, James A.
2015-01-20
The corona and transition region (TR) are fundamentally coupled through the processes of thermal conduction and mass exchange. It is not possible to understand one without the other. Yet the temperature-dependent emissions from the two locations behave quite differently in the aftermath of an impulsive heating event such as a coronal nanoflare. Whereas the corona cools sequentially, emitting first at higher temperatures and then at lower temperatures, the TR is multithermal and the emission at all temperatures responds in unison. We have previously applied the automated time lag technique of Viall and Klimchuk to disk observations of an active region (AR) made by the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory. Lines of sight passing through coronal plasma show clear evidence for post-nanoflare cooling, while lines of sight intersecting the TR footpoints of coronal strands show zero time lag. In this paper, we use the EBTEL hydrodynamics code to demonstrate that this is precisely the expected behavior when the corona is heated by nanoflares. We also apply the time lag technique for the first time to off-limb observations of an AR. Since TR emission is not present above the limb, the occurrence of zero time lags is greatly diminished, supporting the conclusion that zero time lags measured on the disk are due to TR plasma. Lastly, we show that the ''coronal'' channels in AIA can be dominated by bright TR emission. When defined in a physically meaningful way, the TR reaches a temperature of roughly 60% the peak temperature in a flux tube. The TR resulting from impulsive heating can extend to 3MK and higher, well within the range of the ''coronal'' AIA channels.
Hagos, Samson M.; Leung, Lai-Yung R.
2012-11-01
Cloud resolving model simulations and vector analysis are used to develop a quantitative method of assessing regional variations in the relationships between various large-scale environmental variables and the transition to deep convection. Results of the CRM simulations from three tropical regions are used to cluster environmental conditions under which transition to deep convection does and does not take place. Projections of the large-scale environmental variables on the difference between these two clusters are used to quantify the roles of these variables in the transition to deep convection. While the transition to deep convection is most sensitive to moisture and vertical velocity perturbations, the details of the profiles of the anomalies vary from region to region. In comparison, the transition to deep convection is found to be much less sensitive to temperature anomalies over all three regions. The vector formulation presented in this study represents a simple general framework for quantifying various aspects of how the transition to deep convection is sensitive to environmental conditions.
MATHEMATICAL MODELS OF HYSTERESIS (DYNAMIC PROBLEMS IN HYSTERESIS)
Professor Isaak Mayergoyz
2006-08-21
This research has further advanced the current state of the art in the areas of dynamic aspects of hysteresis and nonlinear large scale magnetization dynamics. The results of this research will find important engineering applications in the areas of magnetic data storage technology and the emerging technology of “spintronics”. Our research efforts have been focused on the following tasks: • Study of fast (pulse) precessional switching of magnetization in magnetic materials. • Analysis of critical fields and critical angles for precessional switching of magnetization. • Development of inverse problem approach to the design of magnetic field pulses for precessional switching of magnetization. • Study of magnetization dynamics induced by spin polarized current injection. • Construction of complete stability diagrams for spin polarized current induced magnetization dynamics. • Development of the averaging technique for the analysis of the slow time scale magnetization dynamics. • Study of thermal effects on magnetization dynamics by using the theory of stochastic processes on graphs.
Integrated dynamic landscape analysis and modeling system (IDLAMS) : installation manual.
Li, Z.; Majerus, K. A.; Sundell, R. C.; Sydelko, P. J.; Vogt, M. C.
1999-02-24
The Integrated Dynamic Landscape Analysis and Modeling System (IDLAMS) is a prototype, integrated land management technology developed through a joint effort between Argonne National Laboratory (ANL) and the US Army Corps of Engineers Construction Engineering Research Laboratories (USACERL). Dr. Ronald C. Sundell, Ms. Pamela J. Sydelko, and Ms. Kimberly A. Majerus were the principal investigators (PIs) for this project. Dr. Zhian Li was the primary software developer. Dr. Jeffrey M. Keisler, Mr. Christopher M. Klaus, and Mr. Michael C. Vogt developed the decision analysis component of this project. It was developed with funding support from the Strategic Environmental Research and Development Program (SERDP), a land/environmental stewardship research program with participation from the US Department of Defense (DoD), the US Department of Energy (DOE), and the US Environmental Protection Agency (EPA). IDLAMS predicts land conditions (e.g., vegetation, wildlife habitats, and erosion status) by simulating changes in military land ecosystems for given training intensities and land management practices. It can be used by military land managers to help predict the future ecological condition for a given land use based on land management scenarios of various levels of training intensity. It also can be used as a tool to help land managers compare different land management practices and further determine a set of land management activities and prescriptions that best suit the needs of a specific military installation.
Scott, Erin L.; Mamajek, Eric E.; Pecaut, Mark J.; Quillen, Alice C.; Moolekamp, Fred; Bell, Cameron P. M.
2014-12-10
We investigate the nature of the unusual eclipsing star OGLE LMC-ECL-11893 (OGLE J05172127-6900558) in the Large Magellanic Cloud recently reported by Dong et al. The eclipse period for this star is 468 days, and the eclipses exhibit a minimum of ?1.4 mag, preceded by a plateau of ?0.8 mag. Spectra and optical/IR photometry are consistent with the eclipsed star being a lightly reddened B9III star of inferred age ?150 Myr and mass ?4 M {sub ?}. The disk appears to have an outer radius of ?0.2 AU with predicted temperatures of ?1100-1400 K. We model the eclipses as being due to either a transiting geometrically thin dust disk or gaseous accretion disk around a secondary object; the debris disk produces a better fit. We speculate on the origin of such a dense circumstellar dust disk structure orbiting a relatively old low-mass companion, and on the similarities of this system to the previously discovered EE Cep.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Schryver, Jack; Nutaro, James; Shankar, Mallikarjun
2015-10-30
An agent-based simulation model hierarchy emulating disease states and behaviors critical to progression of diabetes type 2 was designed and implemented in the DEVS framework. The models are translations of basic elements of an established system dynamics model of diabetes. In this model hierarchy, which mimics diabetes progression over an aggregated U.S. population, was dis-aggregated and reconstructed bottom-up at the individual (agent) level. Four levels of model complexity were defined in order to systematically evaluate which parameters are needed to mimic outputs of the system dynamics model. Moreover, the four estimated models attempted to replicate stock counts representing disease statesmore » in the system dynamics model, while estimating impacts of an elderliness factor, obesity factor and health-related behavioral parameters. Health-related behavior was modeled as a simple realization of the Theory of Planned Behavior, a joint function of individual attitude and diffusion of social norms that spread over each agent s social network. Although the most complex agent-based simulation model contained 31 adjustable parameters, all models were considerably less complex than the system dynamics model which required numerous time series inputs to make its predictions. In all three elaborations of the baseline model provided significantly improved fits to the output of the system dynamics model. The performances of the baseline agent-based model and its extensions illustrate a promising approach to translate complex system dynamics models into agent-based model alternatives that are both conceptually simpler and capable of capturing main effects of complex local agent-agent interactions.« less
Wang, Chi -Jen; Liu, Da -Jiang; Evans, James W.
2015-04-28
Threshold versions of Schloegl’s model on a lattice, which involve autocatalytic creation and spontaneous annihilation of particles, can provide a simple prototype for discontinuous non-equilibrium phase transitions. These models are equivalent to so-called threshold contact processes. A discontinuous transition between populated and vacuum states can occur selecting a threshold of N ≥ 2 for the minimum number, N, of neighboring particles enabling autocatalytic creation at an empty site. Fundamental open questions remain given the lack of a thermodynamic framework for analysis. For a square lattice with N = 2, we show that phase coexistence occurs not at a unique value but for a finite range of particle annihilation rate (the natural control parameter). This generic two-phase coexistence also persists when perturbing the model to allow spontaneous particle creation. Such behavior contrasts both the Gibbs phase rule for thermodynamic systems and also previous analysis for this model. We find metastability near the transition corresponding to a non-zero effective line tension, also contrasting previously suggested critical behavior. As a result, mean-field type analysis, extended to treat spatially heterogeneous states, further elucidates model behavior.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Wang, Chi -Jen; Liu, Da -Jiang; Evans, James W.
2015-04-28
Threshold versions of Schloegl’s model on a lattice, which involve autocatalytic creation and spontaneous annihilation of particles, can provide a simple prototype for discontinuous non-equilibrium phase transitions. These models are equivalent to so-called threshold contact processes. A discontinuous transition between populated and vacuum states can occur selecting a threshold of N ≥ 2 for the minimum number, N, of neighboring particles enabling autocatalytic creation at an empty site. Fundamental open questions remain given the lack of a thermodynamic framework for analysis. For a square lattice with N = 2, we show that phase coexistence occurs not at a unique valuemore » but for a finite range of particle annihilation rate (the natural control parameter). This generic two-phase coexistence also persists when perturbing the model to allow spontaneous particle creation. Such behavior contrasts both the Gibbs phase rule for thermodynamic systems and also previous analysis for this model. We find metastability near the transition corresponding to a non-zero effective line tension, also contrasting previously suggested critical behavior. As a result, mean-field type analysis, extended to treat spatially heterogeneous states, further elucidates model behavior.« less
Zamfir, N.V.; McCutchan, E.A.; Casten, R.F.
2004-09-13
We discuss the empirical evolution of structure in the spherical-deformed phase transition region in the context of the behavior of basic structural signatures in the IBA symmetry triangle. The main signatures of a phase/shape transition are a sharp increase in R{sub 4/2} {identical_to} E(4{sub 1}{sup +})/E(2{sub 1}{sup +}) and minima in E(0{sub 2}{sup +}) and E(2{sub {gamma}}{sup +})
JACKSON VL
2011-08-31
The primary purpose of the tank mixing and sampling demonstration program is to mitigate the technical risks associated with the ability of the Hanford tank farm delivery and celtification systems to measure and deliver a uniformly mixed high-level waste (HLW) feed to the Waste Treatment and Immobilization Plant (WTP) Uniform feed to the WTP is a requirement of 24590-WTP-ICD-MG-01-019, ICD-19 - Interface Control Document for Waste Feed, although the exact definition of uniform is evolving in this context. Computational Fluid Dynamics (CFD) modeling has been used to assist in evaluating scaleup issues, study operational parameters, and predict mixing performance at full-scale.
Greenwald, E. E.; North, S. W.; Georgievskii, Y.; Klippenstein, S. J.; Chemistry; Texas A&M Univ.; SNL
2007-06-28
A two transition state model is applied to the prediction of the isomeric branching in the addition of hydroxyl radical to isoprene. The outer transition state is treated with phase space theory fitted to long-range transition state theory calculations on an electrostatic potential energy surface. High-level quantum chemical estimates are applied to the treatment of the inner transition state. A one-dimensional master equation based on an analytic reduction from two-dimensions for a particular statistical assumption about the rotational part of the energy transfer kernel is employed in the calculation of the pressure dependence of the addition process. We find that an accurate treatment of the two separate transition state regions, at the energy and angular momentum resolved level, is essential to the prediction of the temperature dependence of the addition rate. The transition from a dominant outer transition state to a dominant inner transition state is shown to occur at about 275 K, with significant effects from both transition states over the 30-500 K temperature range. Modest adjustments in the ab initio predicted inner saddle point energies yield predictions that are in quantitative agreement with the available high-pressure limit experimental observations and qualitative agreement with those in the falloff regime. The theoretically predicted capture rate is reproduced to within 10% by the expression [1.71 x 10-10(T/298)-2.58 exp(-608.6/RT) + 5.47 x 10-11(T/298)-1.78 exp(-97.3/RT); with R = 1.987 and T in K] cm3 molecule-1 s-1 over the 30-500 K range. A 300 K branching ratio of 0.67:0.02:0.02:0.29 was determined for formation of the four possible OH-isoprene adduct isomers 1, 2, 3, and 4, respectively, and was found to be relatively insensitive to temperature. An Arrhenius activation energy of -0.77 kcal/mol was determined for the high-pressure addition rate constants around 300 K.
Garcia, L.; Carreras, B. A.; Llerena, I.; Calvo, I. [Departamento de Fisica, Universidad Carlos III de Madrid, 28911 Leganes, Madrid (Spain); Departament d'Algebra i Geometria, Facultat de Matematiques, Universitat de Barcelona, Barcelona (Spain); Laboratorio Nacional de Fusion, Asociacion EURATOM-CIEMAT, 28040 Madrid (Spain)
2009-10-15
For the resistive pressure-gradient-driven turbulence model, the transition from laminar regime to fully developed turbulence is not simple and goes through several phases. For low values of the plasma parameter {beta}, a single quasicoherent structure forms. As {beta} increases, several of these structures may emerge and in turn take the dominant role. Finally, at high {beta}, fully developed turbulence with a broad spectrum is established. A suitable characterization of this transition can be given in terms of topological properties of the flow. Here, we analyze these properties that provide an understanding of the turbulence-induced transport and give a measure of the breaking of the homogeneity of the turbulence. To this end, an approach is developed that allows discriminating between topological properties of plasma turbulence flows that are relevant to the transport dynamics and the ones that are not. This is done using computational homology tools and leads to a faster convergence of numerical results for a fixed level of resolution than previously presented in Phys. Rev. E 78, 066402 (2008)
Dynamic Models for Wind Turbines and Wind Power Plants
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
... Each of these models includes representations of general turbine aerodynamics, the ... 9 1.1.2 Wind power integration and wind turbine modeling ......
Dynamic modelling of a double-pendulum gantry crane system incorporating payload
Ismail, R. M. T. Raja; Ahmad, M. A.; Ramli, M. S.; Ishak, R.; Zawawi, M. A.
2011-06-20
The natural sway of crane payloads is detrimental to safe and efficient operation. Under certain conditions, the problem is complicated when the payloads create a double pendulum effect. This paper presents dynamic modelling of a double-pendulum gantry crane system based on closed-form equations of motion. The Lagrangian method is used to derive the dynamic model of the system. A dynamic model of the system incorporating payload is developed and the effects of payload on the response of the system are discussed. Extensive results that validate the theoretical derivation are presented in the time and frequency domains.
Order, chaos and nuclear dynamics: An introduction
Swiatecki, W.J.
1990-08-01
This is an introductory lecture illustrating by simple examples the anticipated effect on collective nuclear dynamics of a transition from order to chaos in the motions of nucleons inside an idealized nucleus. The destruction of order is paralleled by a transition from a rubber-like to a honey-like behaviour of the independent-particle nuclear model. 10 refs., 6 figs.
A dynamic model for underbalanced drilling with coiled tubing
Rommetveit, R.; Vefring, E.H.; Wang, Z.; Bieseman, T.; Faure, A.M.
1995-11-01
A model for underbalanced drilling with coiled tubing has been developed which takes into account all important factors contributing to the process. This model is a unique tool to plan and execute underbalanced or near balance drilling operations. It is a transient, one-dimensional multi-phase flow model with the following components: Lift gas system model, multiphase hydraulics model, reservoir-wellbore interaction model, drilling model, models for multiphase fluids (lift gas, produced gas, mud, foam, produced gas, oil, water and cuttings). Various alternative geometries for gas injection are modeled as well as all important operations during underbalanced drilling with coiled tubing. The model as well as some simulation results for its use are presented in this paper.
Office of Legacy Management (LM)
Dynamic , and Static , Res.ponse of the Government Oil Shale Mine at ' , . , Rifle, ... AND STATIC RESPONSE 'OF THE GOVERNMENT OIL SHALE MINE A T RIFLE, COLORADO, T O THE, ...
Modelling vegetation dynamics at global scale due to climate changes: Comparison of two approaches
Belotelov, N.V.; Bogatyrev, B.G.; Lobanov, A.I.
1996-12-31
Climate changes will influence vegetation dynamics. One of the ways of forecasting these changes is the creation of mathematical models describing vegetation dynamics. Computer experiments can then be conducted under climate change scenarios. Two main approaches are used to create such models. The first approach is based on a bioclimatic dynamic approach. The second approach is based on modelling the main eco-physiological processes. The bioclimatic dynamic approach consists of hypotheses about vegetation types or biomes, and their interrelationships with climate. In the eco-physiological approach, a detailed description of the processes, such as production, mortality, plants migration and their competition is presented. A number of computer experiments has been conducted for several climatic scenario for Russia and the whole world. A qualitative comparison of the results with the results of an earlier bioclimatic model has been done.
A Mechanical Fluid-Dynamical Model For Ground Movements At Campi...
Mechanical Fluid-Dynamical Model For Ground Movements At Campi Flegrei Caldera Jump to: navigation, search OpenEI Reference LibraryAdd to library Journal Article: A Mechanical...
Dynamical Model for Meson Production off Nucleon and Application to Neutrino-Nucleus Reactions
Nakamura, Satoshi X.
2011-11-23
I explain the Sato-Lee (SL) model and its extension to the neutrino-induced pion production off the nucleon. Then I discuss applications of the SL model to incoherent and coherent pion productions in the neutrino-nucleus scattering. I mention a further extension of this approach with a dynamical coupled-channels model developed in Excited Baryon Analysis Center of JLab.
Wang, Shaobu; Lu, Shuai; Zhou, Ning; Lin, Guang; Elizondo, Marcelo A.; Pai, M. A.
2014-09-04
In interconnected power systems, dynamic model reduction can be applied on generators outside the area of interest to mitigate the computational cost with transient stability studies. This paper presents an approach of deriving the reduced dynamic model of the external area based on dynamic response measurements, which comprises of three steps, dynamic-feature extraction, attribution and reconstruction (DEAR). In the DEAR approach, a feature extraction technique, such as singular value decomposition (SVD), is applied to the measured generator dynamics after a disturbance. Characteristic generators are then identified in the feature attribution step for matching the extracted dynamic features with the highest similarity, forming a suboptimal basis of system dynamics. In the reconstruction step, generator state variables such as rotor angles and voltage magnitudes are approximated with a linear combination of the characteristic generators, resulting in a quasi-nonlinear reduced model of the original external system. Network model is un-changed in the DEAR method. Tests on several IEEE standard systems show that the proposed method gets better reduction ratio and response errors than the traditional coherency aggregation methods.
Chapter 18: Understanding the Developing Cellulosic Biofuels Industry through Dynamic Modeling
Newes, E.; Inman, D.; Bush, B.
2011-01-01
The purpose of this chapter is to discuss a system dynamics model called the Biomass Scenario Model (BSM), which is being developed by the U.S. Department of Energy as a tool to better understand the interaction of complex policies and their potential effects on the burgeoning cellulosic biofuels industry in the United States. The model has also recently been expanded to include advanced conversion technologies and biofuels (i.e., conversion pathways that yield biomass-based gasoline, diesel, jet fuel, and butanol), but we focus on cellulosic ethanol conversion pathways here. The BSM uses a system dynamics modeling approach (Bush et al., 2008) built on the STELLA software platform.
Strnsk, Pavel; Macek, Michal; Cejnar, Pavel
2014-06-15
Quantum systems with a finite number of freedom degrees f develop robust singularities in the energy spectrum of excited states as the systems size increases to infinity. We analyze the general form of these singularities for low f, particularly f=2, clarifying the relation to classical stationary points of the corresponding potential. Signatures in the smoothed energy dependence of the quantum state density and in the flow of energy levels with an arbitrary control parameter are described along with the relevant thermodynamical consequences. The general analysis is illustrated with specific examples of excited-state singularities accompanying the first-order quantum phase transition. -- Highlights: ESQPTs found in infinite-size limit of systems with low numbers of freedom degrees f. ESQPTs related to non-analytical evolutions of classical phasespace properties. ESQPT signatures analyzed for general f, particularly f=2, extending known case f=1. ESQPT signatures identified in smoothened density and flow of energy spectrum. ESQPTs shown to induce a new type of thermodynamic anomalies.
Butlitsky, M. A.; Zelener, B. V.
2014-07-14
A two-component plasma model, which we called a shelf Coulomb model has been developed in this work. A Monte Carlo study has been undertaken to calculate equations of state, pair distribution functions, internal energies, and other thermodynamics properties. A canonical NVT ensemble with periodic boundary conditions was used. The motivation behind the model is also discussed in this work. The shelf Coulomb model can be compared to classical two-component (electron-proton) model where charges with zero size interact via a classical Coulomb law. With important difference for interaction of opposite charges: electrons and protons interact via the Coulomb law for large distances between particles, while interaction potential is cut off on small distances. The cut off distance is defined by an arbitrary ? parameter, which depends on system temperature. All the thermodynamics properties of the model depend on dimensionless parameters ? and ? = ?e{sup 2}n{sup 1/3} (where ? = 1/k{sub B}T, n is the particle's density, k{sub B} is the Boltzmann constant, and T is the temperature) only. In addition, it has been shown that the virial theorem works in this model. All the calculations were carried over a wide range of dimensionless ? and ? parameters in order to find the phase transition region, critical point, spinodal, and binodal lines of a model system. The system is observed to undergo a first order gas-liquid type phase transition with the critical point being in the vicinity of ?{sub crit}?13(T{sub crit}{sup *}?0.076),?{sub crit}?1.8(v{sub crit}{sup *}?0.17),P{sub crit}{sup *}?0.39, where specific volume v* = 1/?{sup 3} and reduced temperature T{sup *} = ?{sup ?1}.
Representation of Dormant and Active Microbial Dynamics for Ecosystem Modeling
Wang, Gangsheng; Mayes, Melanie; Gu, Lianhong; Schadt, Christopher Warren
2014-01-01
Dormancy is an essential strategy for microorganisms to cope with environmental stress. However, global ecosystem models typically ignore microbial dormancy, resulting in notable model uncertainties. To facilitate the consideration of dormancy in these large-scale models, we propose a new microbial physiology component that works for a wide range of substrate availabilities. This new model is based on microbial physiological states and the major parameters are the maximum specific growth and maintenance rates of active microbes and the ratio of dormant to active maintenance rates. A major improvement of our model over extant models is that it can explain the low active microbial fractions commonly observed in undisturbed soils. Our new model shows that the exponentially-increasing respiration from substrate-induced respiration experiments can only be used to determine the maximum specific growth rate and initial active microbial biomass, while the respiration data representing both exponentially-increasing and non-exponentially-increasing phases can robustly determine a range of key parameters including the initial total live biomass, initial active fraction, the maximum specific growth and maintenance rates, and the half-saturation constant. Our new model can be incorporated into existing ecosystem models to account for dormancy in microbially-driven processes and to provide improved estimates of microbial activities.
Applied Dynamic Analysis of the Global Economy (ADAGE) Model...
model capable of examining many types of economic, energy, environmental, climate change mitigation, and trade policies at the international, national, U.S. regional, and U.S....
Modeling icesheets dynamics: forward and inverse problems. (Conference...
Office of Scientific and Technical Information (OSTI)
for Modeling and Scientific Computing MOX Seminar held September 9, 2013 in Milano, Italy. ... Country of Publication: United States Language: English Word Cloud More Like This Full ...
Dynamic Models for Wind Turbines and Wind Power Plants
Singh, M.; Santoso, S.
2011-10-01
The primary objective of this report was to develop universal manufacturer-independent wind turbine and wind power plant models that can be shared, used, and improved without any restrictions by project developers, manufacturers, and engineers. Manufacturer-specific models of wind turbines are favored for use in wind power interconnection studies. While they are detailed and accurate, their usages are limited to the terms of the non-disclosure agreement, thus stifling model sharing. The primary objective of the work proposed is to develop universal manufacturer-independent wind power plant models that can be shared, used, and improved without any restrictions by project developers, manufacturers, and engineers. Each of these models includes representations of general turbine aerodynamics, the mechanical drive-train, and the electrical characteristics of the generator and converter, as well as the control systems typically used. To determine how realistic model performance is, the performance of one of the models (doubly-fed induction generator model) has been validated using real-world wind power plant data. This work also documents selected applications of these models.
Ansong, Charles; Sadler, Natalie C.; Hill, Eric A.; Lewis, Michael P.; Zink, Erika M.; Smith, Richard D.; Beliaev, Alex S.; Konopka, Allan; Wright, Aaron T.
2014-07-03
Protein redox chemistry constitutes a major void in knowledge pertaining to photoautotrophic system regulation and signaling processes. We have employed a chemical biology approach to analyze redox sensitive proteins in live Synechococcus sp. PCC 7002 cells in both light and dark periods, and to understand how cellular redox balance is disrupted during nutrient perturbation. The present work identified several novel putative redox-sensitive proteins that are involved in the generation of reductant, macromolecule synthesis, and carbon flux through central metabolic pathways, and may be involved in cell signaling and response mechanisms. Furthermore, our research suggests that dynamic redox changes in response to specific nutrient limitations contribute to the regulatory changes driven by a shift from light to dark. Taken together, these results contribute to the high-level understanding of post-translational mechanisms regulating flux distributions and therefore present potential metabolic engineering targets for redirecting carbon towards biofuel precursors.
Dynamic Model Validation of PV Inverters Under Short-Circuit Conditions: Preprint
Muljadi, E.; Singh, M.; Bravo, R.; Gevorgian, V.
2013-03-01
Photovoltaic (PV) modules have dramatically decreased in price in the past few years, spurring the expansion of photovoltaic deployment. Residential and commercial rooftop installations are connected to the distribution network; large-scale installation PV power plants (PVPs) have benefited from tax incentives and the low cost of PV modules. As the level penetration of PV generation increases, the impact on power system reliability will also be greater. Utility power system planners must consider the role of PV generation in power systems more realistically by representing PV generation in dynamic stability analyses. Dynamic models of PV inverters have been developed in the positive sequence representation. NREL has developed a PV inverter dynamic model in PSCAD/EMTDC. This paper validates the dynamic model with an actual hardware bench test conducted by Southern California Edison's Distributed Energy Resources laboratory. All the fault combinations -- symmetrical and unsymmetrical -- were performed in the laboratory. We compare the simulation results with the bench test results.
Dynamics of pentaquarks in constituent quark models: recent developments
Stancu, Fl.
2005-06-14
Some recent developments in the study of light and heavy pentaquarks are reviewed, mainly within constituent quark models. Emphasis is made on results obtained in the flavor-spin model where a nearly ideal octet-antidecuplet mixing is obtained. The charmed antisextet is reviewed in the context of an SU(4) classification.
Liese, Eric; Zitney, Stephen E.
2013-01-01
Research in dynamic process simulation for integrated gasification combined cycles (IGCC) with carbon capture has been ongoing at the National Energy Technology Laboratory (NETL), culminating in a full operator training simulator (OTS) and immersive training simulator (ITS) for use in both operator training and research. A derivative work of the IGCC dynamic simulator has been a modification of the combined cycle section to more closely represent a typical natural gas fired combined cycle (NGCC). This paper describes the NGCC dynamic process model and highlights some of the simulator’s current capabilities through a particular startup and shutdown scenario.
Human Performance Modeling for Dynamic Human Reliability Analysis
Boring, Ronald Laurids; Joe, Jeffrey Clark; Mandelli, Diego
2015-08-01
Part of the U.S. Department of Energy’s (DOE’s) Light Water Reac- tor Sustainability (LWRS) Program, the Risk-Informed Safety Margin Charac- terization (RISMC) Pathway develops approaches to estimating and managing safety margins. RISMC simulations pair deterministic plant physics models with probabilistic risk models. As human interactions are an essential element of plant risk, it is necessary to integrate human actions into the RISMC risk framework. In this paper, we review simulation based and non simulation based human reliability analysis (HRA) methods. This paper summarizes the founda- tional information needed to develop a feasible approach to modeling human in- teractions in RISMC simulations.
Advanced Modeling of Renewable Energy Market Dynamics: May 2006
Evans, M.; Little, R.; Lloyd, K.; Malikov, G.; Passolt, G.; Arent, D.; Swezey, B.; Mosey, G.
2007-08-01
This report documents a year-long academic project, presenting selected techniques for analysis of market growth, penetration, and forecasting applicable to renewable energy technologies. Existing mathematical models were modified to incorporate the effects of fiscal policies and were evaluated using available data. The modifications were made based on research and classification of current mathematical models used for predicting market penetration. An analysis of the results was carried out, based on available data. MATLAB versions of existing and new models were developed for research and policy analysis.
Coupled Dynamic Modeling of Floating Wind Turbine Systems: Preprint
Wayman, E. N.; Sclavounos, P. D.; Butterfield, S.; Jonkman, J.; Musial, W.
2006-03-01
This article presents a collaborative research program that the Massachusetts Institute of Technology (MIT) and the National Renewable Energy Laboratory (NREL) have undertaken to develop innovative and cost-effective floating and mooring systems for offshore wind turbines in water depths of 10-200 m. Methods for the coupled structural, hydrodynamic, and aerodynamic analysis of floating wind turbine systems are presented in the frequency domain. This analysis was conducted by coupling the aerodynamics and structural dynamics code FAST [4] developed at NREL with the wave load and response simulation code WAMIT (Wave Analysis at MIT) [15] developed at MIT. Analysis tools were developed to consider coupled interactions between the wind turbine and the floating system. These include the gyroscopic loads of the wind turbine rotor on the tower and floater, the aerodynamic damping introduced by the wind turbine rotor, the hydrodynamic damping introduced by wave-body interactions, and the hydrodynamic forces caused by wave excitation. Analyses were conducted for two floater concepts coupled with the NREL 5-MW Offshore Baseline wind turbine in water depths of 10-200 m: the MIT/NREL Shallow Drafted Barge (SDB) and the MIT/NREL Tension Leg Platform (TLP). These concepts were chosen to represent two different methods of achieving stability to identify differences in performance and cost of the different stability methods. The static and dynamic analyses of these structures evaluate the systems' responses to wave excitation at a range of frequencies, the systems' natural frequencies, and the standard deviations of the systems' motions in each degree of freedom in various wind and wave environments. This article in various wind and wave environments. This article explores the effects of coupling the wind turbine with the floating platform, the effects of water depth, and the effects of wind speed on the systems' performance. An economic feasibility analysis of the two concepts
Modeling Dynamic Ductility: An Equation of State for Porous Metals...
Office of Scientific and Technical Information (OSTI)
with material elastic-plastic response in a 2D hydrocode, and then discuss the modeling of ... Country of Publication: United States Language: English Subject: 36 MATERIALS SCIENCE; 75 ...
Dynamics of holographic vacuum energy in the DGP model
Wu Xing; Zhu Zonghong; Cai Ronggen
2008-02-15
We consider the evolution of the vacuum energy in the Dvali-Gabadadze-Porrati (DGP) model according to the holographic principle under the assumption that the relation linking the IR and UV cutoffs still holds in this scenario. The model is studied when the IR cutoff is chosen to be the Hubble scale H{sup -1}, the particle horizon R{sub ph}, and the future event horizon R{sub eh}, respectively. The two branches of the DGP model are also taken into account. Through numerical analysis, we find that in the cases of H{sup -1} in the (+) branch and R{sub eh} in both branches, the vacuum energy can play the role of dark energy. Moreover, when considering the combination of the vacuum energy and the 5D gravity effect in both branches, the equation of state of the effective dark energy may cross -1, which may lead to the big rip singularity. Besides, we constrain the model with the Type Ia supernovae and baryon oscillation data and find that our model is consistent with current data within 1{sigma}, and that the observations prefer either a pure holographic dark energy or a pure DGP model.
Link, S.O.; Kickert, R.N.; Fayer, M.J.; Gee, G.W.
1993-06-01
This report describes the results of simulation models used to predict soil water storage dynamics at the Field Lysimeter Test Facility (FLTF) weighing lysimeters. The objectives of this research is to develop the capability to predict soil water storage dynamics with plants in support of water infiltration control studies for the Hanford Permanent Isolation Barrier Development Program. It is important to gain confidence in one`s ability to simulate soil water dynamics over long time periods to assess the barrier`s ability to prevent drainage. Two models were compared for their ability to simulate soil water storage dynamics with and without plants in weighing lysimeters, the soil water infiltration and movement (SWIM) and the simulation of production and utilization of rangelands (SPUR-91) models. These models adequately simulated soil water storage dynamics for the weighing lysimeters. The range of root mean square error values for the two models was 7.0 to 19.8. This compares well with the range reported by Fayer et al. (1992) for the bare soil data sets of 8.1 to 22.1. Future research will test the predictive capability of these models for longer term lysimeter data sets and for historical data sets collected in various plant community types.
Jackson, Bret; Nattino, Francesco; Kroes, Geert-Jan
2014-08-07
The dissociative chemisorption of methane on metal surfaces is of great practical and fundamental importance. Not only is it the rate-limiting step in the steam reforming of natural gas, the reaction exhibits interesting mode-selective behavior and a strong dependence on the temperature of the metal. We present a quantum model for this reaction on Ni(100) and Ni(111) surfaces based on the reaction path Hamiltonian. The dissociative sticking probabilities computed using this model agree well with available experimental data with regard to variation with incident energy, substrate temperature, and the vibrational state of the incident molecule. We significantly expand the vibrational basis set relative to earlier studies, which allows reaction probabilities to be calculated for doubly excited initial vibrational states, though it does not lead to appreciable changes in the reaction probabilities for singly excited initial states. Sudden models used to treat the center of mass motion parallel to the surface are compared with results from ab initio molecular dynamics and found to be reasonable. Similar comparisons for molecular rotation suggest that our rotationally adiabatic model is incorrect, and that sudden behavior is closer to reality. Such a model is proposed and tested. A model for predicting mode-selective behavior is tested, with mixed results, though we find it is consistent with experimental studies of normal vs. total (kinetic) energy scaling. Models for energy transfer into lattice vibrations are also examined.
Data Driven Approach for High Resolution Population Distribution and Dynamics Models
Bhaduri, Budhendra L; Bright, Eddie A; Rose, Amy N; Liu, Cheng; Urban, Marie L; Stewart, Robert N
2014-01-01
High resolution population distribution data are vital for successfully addressing critical issues ranging from energy and socio-environmental research to public health to human security. Commonly available population data from Census is constrained both in space and time and does not capture population dynamics as functions of space and time. This imposes a significant limitation on the fidelity of event-based simulation models with sensitive space-time resolution. This paper describes ongoing development of high-resolution population distribution and dynamics models, at Oak Ridge National Laboratory, through spatial data integration and modeling with behavioral or activity-based mobility datasets for representing temporal dynamics of population. The model is resolved at 1 km resolution globally and describes the U.S. population for nighttime and daytime at 90m. Integration of such population data provides the opportunity to develop simulations and applications in critical infrastructure management from local to global scales.
Modeling microbial dynamics in heterogeneous environments: Growth on soil carbon sources
Resat, Haluk; Bailey, Vanessa L.; McCue, Lee Ann; Konopka, Allan
2012-01-01
We have developed a new hybrid model to study how microbial dynamics are affected by the heterogeneity in the physical structure of the environment. The modeling framework can represent porous media such as soil. The individual based biological model can explicitly simulate microbial diversity, and cell metabolism is regulated via optimal allocation of cellular resources to enzyme synthesis, control of growth rate by protein synthesis capacity, and shifts to dormancy. This model was developed to study how microbial community functioning is influenced by local environmental conditions and by the functional attributes of individual microbes. Different strategies for acquisition of carbon from polymeric cellulose were investigated. Bacteria that express membrane-associated hydrolase had different growth and survival dynamics in soil pores than bacteria that release extracellular hydrolases. The kinetic differences may suggest different functional roles for these two classes of microbes in cellulose utilization. Our model predicted an emergent behavior in which co-existence led to higher cellulose utilization efficiency and reduced stochasticity. Microbial community dynamics were simulated at two spatial scales: micro-pores that resemble 6-20 {micro}m size portions of the soil physical structure and in 111 {micro}m size soil aggregates with a random pore structure. Trends in dynamic properties were very similar at these two scales, implying that micro-scale studies can be useful approximations to aggregate scale studies when local effects on microbial dynamics are studied.
Dynamic (G2) Model Design Document, 24590-WTP-MDD-PR-01-002, Rev. 12
Deng, Yueying; Kruger, Albert A.
2013-12-16
The Hanford Tank Waste Treatment and Immobilization Plant (WTP) Statement of Work (Department of Energy Contract DE-AC27-01RV14136, Section C) requires the contractor to develop and use process models for flowsheet analyses and pre-operational planning assessments. The Dynamic (G2) Flowsheet is a discrete-time process model that enables the project to evaluate impacts to throughput from eventdriven activities such as pumping, sampling, storage, recycle, separation, and chemical reactions. The model is developed by the Process Engineering (PE) department, and is based on the Flowsheet Bases, Assumptions, and Requirements Document (24590-WTP-RPT-PT-02-005), commonly called the BARD. The terminologies of Dynamic (G2) Flowsheet and Dynamic (G2) Model are interchangeable in this document. The foundation of this model is a dynamic material balance governed by prescribed initial conditions, boundary conditions, and operating logic. The dynamic material balance is achieved by tracking the storage and material flows within the plant as time increments. The initial conditions include a feed vector that represents the waste compositions and delivery sequence of the Tank Farm batches, and volumes and concentrations of solutions in process equipment before startup. The boundary conditions are the physical limits of the flowsheet design, such as piping, volumes, flowrates, operation efficiencies, and physical and chemical environments that impact separations, phase equilibriums, and reaction extents. The operating logic represents the rules and strategies of running the plant.
Wang, Gangsheng; Post, Wilfred M; Mayes, Melanie
2013-01-01
We developed a Microbial-ENzyme-mediated Decomposition (MEND) model, based on the Michaelis-Menten kinetics, that describes the dynamics of physically defined pools of soil organic matter (SOC). These include particulate, mineral-associated, dissolved organic matter (POC, MOC, and DOC, respectively), microbial biomass, and associated exoenzymes. The ranges and/or distributions of parameters were determined by both analytical steady-state and dynamic analyses with SOC data from the literature. We used an improved multi-objective parameter sensitivity analysis (MOPSA) to identify the most important parameters for the full model: maintenance of microbial biomass, turnover and synthesis of enzymes, and carbon use efficiency (CUE). The model predicted an increase of 2 C (baseline temperature =12 C) caused the pools of POC-Cellulose, MOC, and total SOC to increase with dynamic CUE and decrease with constant CUE, as indicated by the 50% confidence intervals. Regardless of dynamic or constant CUE, the pool sizes of POC, MOC, and total SOC varied from 8% to 8% under +2 C. The scenario analysis using a single parameter set indicates that higher temperature with dynamic CUE might result in greater net increases in both POC-Cellulose and MOC pools. Different dynamics of various SOC pools reflected the catalytic functions of specific enzymes targeting specific substrates and the interactions between microbes, enzymes, and SOC. With the feasible parameter values estimated in this study, models incorporating fundamental principles of microbial-enzyme dynamics can lead to simulation results qualitatively different from traditional models with fast/slow/passive pools.
Gender Transition Guidelines, Workplace
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Gender Transition Guidelines, Workplace
Welsch, Ralph Manthe, Uwe
2014-08-07
The mode-selective chemistry of the title reaction is studied by full-dimensional quantum dynamics simulation on an accurate ab initio potential energy surface for vanishing total angular momentum. Using a rigorous transition state based approach and multi-configurational time-dependent Hartree wave packet propagation, initial state-selected reaction probabilities for many ro-vibrational states of methane are calculated. The theoretical results are compared with experimental trends seen in reactions of methane. An intuitive interpretation of the ro-vibrational control of the chemical reactivity provided by a sudden model based on the quantum transition state concept is discussed.
Comparison of transition densities in the DDHMS model of pre-equilibrium emission
Brito, L.; Carlson, B. V.
2014-11-11
The DDHMS (double differential hybrid Monte Carlo simulation) model treats nucleon-induced pre-equilibrium reactions as a series of particle-particle and particle-hole interactions in the space of energy and angle. This work compares spectra obtained within the model using diferent approximations to the density of accessible states. The calculations are performed with the EMPIRE reaction model code, a modular system containing several nuclear reaction models that permits a fairly complete descritpion of the reaction, from elastic scattering and absorption through the pre-equilbrium stage to the final decay by statistical emission.
Parameter Estimation and Model Validation of Nonlinear Dynamical Networks
Abarbanel, Henry; Gill, Philip
2015-03-31
In the performance period of this work under a DOE contract, the co-PIs, Philip Gill and Henry Abarbanel, developed new methods for statistical data assimilation for problems of DOE interest, including geophysical and biological problems. This included numerical optimization algorithms for variational principles, new parallel processing Monte Carlo routines for performing the path integrals of statistical data assimilation. These results have been summarized in the monograph: “Predicting the Future: Completing Models of Observed Complex Systems” by Henry Abarbanel, published by Spring-Verlag in June 2013. Additional results and details have appeared in the peer reviewed literature.
Multiple dynamical resonances in a discrete neuronal model
Jiang Yu
2005-05-01
The conditions for the occurrence of different multiple resonances in an excitable neuron model are analyzed numerically. It is shown that the excitable system may display both stochastic and coherence resonance, in response to periodic stimuli in the presence of different intensities of additive and parametric noises. It is found that double coherence resonances may take place in the low-amplitude oscillation regimes, and coherence resonance may persists even in the weak oscillatory regimes for control parameters slightly larger than the Hopf bifurcation point, where the system is in the incipient stage of large-amplitude excitation regime.
Modeling aspects of the dynamic response of heterogeneous materials
Ionita, Axinte; Clements, Brad; Mas, Eric
2009-01-01
In numerical simulations of engineering applications involving heterogeneous materials capturing the local response coming from a distribution of heterogeneities can lead to a very large model thus making simulations difficult. The use of homogenization techniques can reduce the size of the problem but will miss the local effects. Homogenization can also be difficult if the constituents obey different types of constitutive laws. Additional complications arise if inelastic deformation. In such cases a two-scale approach is prefened and tills work addresses these issues in the context of a two-scale Finite Element Method (FEM). Examples of using two-scale FEM approaches are presented.
GOLDSCHMIDT, YADIN Y.; LIU, Jin-Tao
2007-08-07
In this paper we use London Langevin molecular dynamics simulations to investigate the vortex matter melting transition in the highly anisotropic high-temperature superconductor material Bi{sub 2}Sr{sub 2}CaCu{sub 2}O{sub 8+{delta}}#14; in the presence of low concentration of columnar defects (CDs). We reproduce with further details our previous results obtained by using Multilevel Monte Carlo simulations that showed that the melting of the nanocrystalline vortex matter occurs in two stages: a first stage melting into nanoliquid vortex matter and a second stage delocalization transition into a homogeneous liquid. Furthermore, we report on new dynamical measurements in the presence of a current that identifies clearly the irreversibility line and the second stage delocalization transition. In addition to CDs aligned along the c-axis we also simulate the case of tilted CDs which are aligned at an angle with respect to the applied magnetic field. Results for CDs tilted by 45{degree} with respect to c-axis show that the locations of the melting and delocalization transitions are not affected by the tilt when the ratio of flux lines to CDs remains constant. On the other hand we argue that some dynamical properties and in particular the position of the irreversibility line should be affected.
Advanced Neutron Source Dynamic Model (ANSDM) code description and user guide
March-Leuba, J.
1995-08-01
A mathematical model is designed that simulates the dynamic behavior of the Advanced Neutron Source (ANS) reactor. Its main objective is to model important characteristics of the ANS systems as they are being designed, updated, and employed; its primary design goal, to aid in the development of safety and control features. During the simulations the model is also found to aid in making design decisions for thermal-hydraulic systems. Model components, empirical correlations, and model parameters are discussed; sample procedures are also given. Modifications are cited, and significant development and application efforts are noted focusing on examination of instrumentation required during and after accidents to ensure adequate monitoring during transient conditions.
Update on Small Modular Reactors Dynamics System Modeling Tool -- Molten Salt Cooled Architecture
Hale, Richard Edward; Cetiner, Sacit M.; Fugate, David L.; Qualls, A L.; Borum, Robert C.; Chaleff, Ethan S.; Rogerson, Doug W.; Batteh, John J.; Tiller, Michael M.
2014-08-01
The Small Modular Reactor (SMR) Dynamic System Modeling Tool project is in the third year of development. The project is designed to support collaborative modeling and study of various advanced SMR (non-light water cooled) concepts, including the use of multiple coupled reactors at a single site. The objective of the project is to provide a common simulation environment and baseline modeling resources to facilitate rapid development of dynamic advanced reactor SMR models, ensure consistency among research products within the Instrumentation, Controls, and Human-Machine Interface (ICHMI) technical area, and leverage cross-cutting capabilities while minimizing duplication of effort. The combined simulation environment and suite of models are identified as the Modular Dynamic SIMulation (MoDSIM) tool. The critical elements of this effort include (1) defining a standardized, common simulation environment that can be applied throughout the program, (2) developing a library of baseline component modules that can be assembled into full plant models using existing geometry and thermal-hydraulic data, (3) defining modeling conventions for interconnecting component models, and (4) establishing user interfaces and support tools to facilitate simulation development (i.e., configuration and parameterization), execution, and results display and capture.
Smith, Kyle K. G.; Poulsen, Jens Aage Nyman, Gunnar; Rossky, Peter J.
2015-06-28
We develop two classes of quasi-classical dynamics that are shown to conserve the initial quantum ensemble when used in combination with the Feynman-Kleinert approximation of the density operator. These dynamics are used to improve the Feynman-Kleinert implementation of the classical Wigner approximation for the evaluation of quantum time correlation functions known as Feynman-Kleinert linearized path-integral. As shown, both classes of dynamics are able to recover the exact classical and high temperature limits of the quantum time correlation function, while a subset is able to recover the exact harmonic limit. A comparison of the approximate quantum time correlation functions obtained from both classes of dynamics is made with the exact results for the challenging model problems of the quartic and double-well potentials. It is found that these dynamics provide a great improvement over the classical Wigner approximation, in which purely classical dynamics are used. In a special case, our first method becomes identical to centroid molecular dynamics.
IMPROVED MODELING OF THE ROSSITER-McLAUGHLIN EFFECT FOR TRANSITING EXOPLANETS
Hirano, Teruyuki; Winn, Joshua N.; Albrecht, Simon; Suto, Yasushi; Taruya, Atsushi; Narita, Norio; Sato, Bun'ei
2011-12-01
We present an improved formula for the anomalous radial velocity of the star during planetary transits due to the Rossiter-McLaughlin (RM) effect. The improvement comes from a more realistic description of the stellar absorption line profiles, taking into account stellar rotation, macroturbulence, thermal broadening, pressure broadening, and instrumental broadening. Although the formula is derived for the case in which radial velocities are measured by cross-correlation, we show through numerical simulations that the formula accurately describes the cases where the radial velocities are measured with the iodine absorption-cell technique. The formula relies on prior knowledge of the parameters describing macroturbulence, instrumental broadening, and other broadening mechanisms, but even 30% errors in those parameters do not significantly change the results in typical circumstances. We show that the new analytic formula agrees with previous ones that had been computed on a case-by-case basis via numerical simulations. Finally, as one application of the new formula, we reassess the impact of the differential rotation on the RM velocity anomaly. We show that differential rotation of a rapidly rotating star may have a significant impact on future RM observations.
A semi-analytic power balance model for low (L) to high (H) mode transition power threshold
Singh, R., E-mail: rsingh129@yahoo.co.in [WCI Center for Fusion Theory, National Fusion Research Institute, Daejeon 305-333 (Korea, Republic of); Institute for Plasma Research, Bhat Gandhinagar 2382 428 (India); Jhang, Hogun [WCI Center for Fusion Theory, National Fusion Research Institute, Daejeon 305-333 (Korea, Republic of); Kaw, P. K. [Institute for Plasma Research, Bhat Gandhinagar 2382 428 (India); Diamond, P. H. [WCI Center for Fusion Theory, National Fusion Research Institute, Daejeon 305-333 (Korea, Republic of); Center for Momentum Transport and Flow Organization, University of California, San Diego, California 92093 (United States); Center for Astrophysics and Space Sciences, University of California, San Diego, 9500 Gilman Dr., La Jolla, California 92093-0424 (United States); Nordman, H. [Department of Earth and Space Sciences, Chalmers University of Technology, SE-412 96 Gteborg (Sweden); Bourdelle, C. [Euratom-CEA Association, CEA/DSM/DRFC, CEA Cadarache F-13108 Saint-Paul-Lez-Durance (France); Loarte, A. [ITER Organization, Route de Vinon Sur Verdon, A. 13115 Saint Paul Lez Durance (France)
2014-06-15
We present a semi-analytic model for low (L) to high (H) mode transition power threshold (P{sub th}). Two main assumptions are made in our study. First, high poloidal mode number drift resistive ballooning modes (high-m DRBM) are assumed to be the dominant turbulence driver in a narrow edge region near to last closed flux surface. Second, the pre-transition edge profile and turbulent diffusivity at the narrow edge region pertain to turbulent equipartition. An edge power balance relation is derived by calculating the dissipated power flux through both turbulent conduction and convection, and radiation in the edge region. P{sub th} is obtained by imposing the turbulence quench rule due to sheared E??B rotation. Evaluation of P{sub th} shows a good agreement with experimental results in existing machines. Increase of P{sub th} at low density (i.e., the existence of roll-over density in P{sub th} vs. density) is shown to originate from the longer scale length of the density profile than that of the temperature profile.
Rudolph, K.W. ); Schlager, W. ); Biddle, K.T. )
1989-05-01
Detailed geologic cross sections, augmented by laboratory calibration of lithology and acoustic impedance, have been used to produce synthetic reflection seismic sections of a carbonate foreslope-to-basin transition. Two areas from the Picco di Vallandro region of the Dolomite Alps were modeled: a progradational section and a retrogradational (backstepping) section. The resulting models show how these complex areas of strata interfingering might be displayed on conventional reflection seismic lines. In the area of progradation, rapid stratigraphic thinning below seismic detectibility, coupled with abrupt impedance changes, produces a reflection discontinuity between steeply dipping reflections of the foreslope and gently dipping paralle reflections of the basin section. This apparent downlap surface marks the toe-of-slope for successive clinoforms but dose not correspond to a discrete stratigraphic surface. In the backstepping example, similar stratigraphic thinning and impedance changes create an apparent onlap surface. Wavelet interference causes complications in both examples. These models indicate how stratigraphic complexity can be simplified by the seismic reflection process and suggest that caution should be exercised when using seismic data to construct general models in areas of complex depositional geometries and rapidly changing facies.
Dynamic Modeling of Adjustable-Speed Pumped Storage Hydropower Plant: Preprint
Muljadi, E.; Singh, M.; Gevorgian, V.; Mohanpurkar, M.; Havsapian, R.; Koritarov, V.
2015-04-06
Hydropower is the largest producer of renewable energy in the U.S. More than 60% of the total renewable generation comes from hydropower. There is also approximately 22 GW of pumped storage hydropower (PSH). Conventional PSH uses a synchronous generator, and thus the rotational speed is constant at synchronous speed. This work details a hydrodynamic model and generator/power converter dynamic model. The optimization of the hydrodynamic model is executed by the hydro-turbine controller, and the electrical output real/reactive power is controlled by the power converter. All essential controllers to perform grid-interface functions and provide ancillary services are included in the model.
Corradini, D.; Rovere, M.; Gallo, P.
2015-09-21
In a previous study [Gallo et al., Nat. Commun. 5, 5806 (2014)], we have shown an important connection between thermodynamic and dynamical properties of water in the supercritical region. In particular, by analyzing the experimental viscosity and the diffusion coefficient obtained in simulations performed using the TIP4P/2005 model, we have found that the line of response function maxima in the one phase region, the Widom line, is connected to a crossover from a liquid-like to a gas-like behavior of the transport coefficients. This is in agreement with recent experiments concerning the dynamics of supercritical simple fluids. We here show how different popular water models (TIP4P/2005, TIP4P, SPC/E, TIP5P, and TIP3P) perform in reproducing thermodynamic and dynamic experimental properties in the supercritical region. In particular, the comparison with experiments shows that all the analyzed models are able to qualitatively predict the dynamical crossover from a liquid-like to a gas-like behavior upon crossing the Widom line. Some of the models perform better in reproducing the pressure-temperature slope of the Widom line of supercritical water once a rigid shift of the phase diagram is applied to bring the critical points to coincide with the experimental ones.
Overview of the synergia 3-D multi-particle dynamics modeling framework
Spentzouris, P.; Amundson, J.F.; Dechow, D.R.; /Tech-X, Boulder
2005-05-01
High precision modeling of space-charge effects is essential for designing future accelerators as well as optimizing the performance of existing machines. Synergia is a high-fidelity parallel beam dynamics simulation package with fully three dimensional space-charge capabilities and a higher-order optics implementation. We describe the Synergia framework and model benchmarks we obtained by comparing to semi-analytic results and other codes. We also present Synergia simulations of the Fermilab Booster accelerator and comparisons with experiment.
Dynamical Coupled-Channel Model of Meson Production Reactions in the Nucleon Resonance Region
T.-S. H. Lee; A. Matsuyama; T. Sato
2006-11-15
A dynamical coupled-channel model is presented for investigating the nucleon resonances (N*) in the meson production reactions induced by pions and photons. Our objective is to extract the N* parameters and to investigate the meson production reaction mechanisms for mapping out the quark-gluon substructure of N* from the data. The model is based on an energy-independent Hamiltonian which is derived from a set of Lagrangians by using a unitary transformation method.
Experimental estimation of dynamic plastic bending moments by plastic hinge models
Sogo, T.; Ujihashi, S.; Matsumoto, H.; Adachi, T.
1995-12-31
In the present paper, the experimental estimation of dynamic plastic bending moments for metallic materials is investigated. The three-point bending, test under impact and static loads is applied to aluminum alloy (JIS A6063S) and mild steel (JIS SS400). It is confirmed that tile dynamic bending deformations in three-point bending test can be modeled as a plastic hinge, tile experimental results show that the consumed energies of the specimens are proportional to the bending angles. The ratio of the consumed energy to the bending angle is approximately equal to the plastic bending moment. In the case of aluminum alloy, the dynamic plastic bending moments for the different average bending angular velocities coincide with the static plastic bending moments. On the other hand, in the case of mild steel, the dynamic plastic bending moments are proportional to the average bending angular velocities. As a result, we confirm that the present method based on the plastic hinge model and the consumed energy is efficient for determining tile dynamic plastic bending moment.
Zhu, H.; Chapman, S. C.; Dendy, R. O.; Itoh, K.
2014-06-15
It is shown that rapid substantial changes in heating rate can induce transitions to improved energy confinement regimes in zero-dimensional models for tokamak plasma phenomenology. We examine for the first time the effect of step changes in heating rate in the models of Kim and Diamond [Phys. Rev. Lett. 90, 185006 (2003)] and Malkov and Diamond [Phys. Plasmas 16, 012504 (2009)], which nonlinearly couple the evolving temperature gradient, micro-turbulence, and a mesoscale flow; and in the extension of Zhu et al. [Phys. Plasmas 20, 042302 (2013)], which couples to a second mesoscale flow component. The temperature gradient rises, as does the confinement time defined by analogy with the fusion context, while micro-turbulence is suppressed. This outcome is robust against variation of heating rise time and against introduction of an additional variable into the model. It is also demonstrated that oscillating changes in heating rate can drive the level of micro-turbulence through a period-doubling path to chaos, where the amplitude of the oscillatory component of the heating rate is the control parameter.
Technical Review of the CENWP Computational Fluid Dynamics Model of the John Day Dam Forebay
Rakowski, Cynthia L.; Serkowski, John A.; Richmond, Marshall C.
2010-12-01
The US Army Corps of Engineers Portland District (CENWP) has developed a computational fluid dynamics (CFD) model of the John Day forebay on the Columbia River to aid in the development and design of alternatives to improve juvenile salmon passage at the John Day Project. At the request of CENWP, Pacific Northwest National Laboratory (PNNL) Hydrology Group has conducted a technical review of CENWP's CFD model run in CFD solver software, STAR-CD. PNNL has extensive experience developing and applying 3D CFD models run in STAR-CD for Columbia River hydroelectric projects. The John Day forebay model developed by CENWP is adequately configured and validated. The model is ready for use simulating forebay hydraulics for structural and operational alternatives. The approach and method are sound, however CENWP has identified some improvements that need to be made for future models and for modifications to this existing model.
Environmental patterns from free trade: Implications from dynamic NAFTA models of Mexico
Gale, L.R.
1994-12-31
Studies of the economic impact on Mexico from joining the North American Free Trade Agreement (NAFTA) point to significant dynamic gains from trade. Few studies, however, have effectively related these changes in economic variables to changes in environmental variables. Using sector-share relationships and projections of income growth from dynamic computable general equilibrium models of Mexico, several possible time paths for carbon dioxide emissions are derived. Various scenarios of trade and investment liberalization and increased fuel efficiency under NAFTA result in simulated pollution paths of carbon dioxide that show a reduction in not only the amount of emissions but also in rate of growth of emissions.
On the characteristics of aerosol indirect effect based on dynamic regimes in global climate models
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Zhang, S.; Wang, M.; Ghan, S. J.; Ding, A.; Wang, H.; Zhang, K.; Neubauer, D.; Lohmann, U.; Ferrachat, S.; Takeamura, T.; et al
2015-09-02
Aerosol-cloud interactions continue to constitute a major source of uncertainty for the estimate of climate radiative forcing. The variation of aerosol indirect effects (AIE) in climate models is investigated across different dynamical regimes, determined by monthly mean 500 hPa vertical pressure velocity (?500), lower-tropospheric stability (LTS) and large-scale surface precipitation rate derived from several global climate models (GCMs), with a focus on liquid water path (LWP) response to cloud condensation nuclei (CCN) concentrations. The LWP sensitivity to aerosol perturbation within dynamic regimes is found to exhibit a large spread among these GCMs. It is in regimes of strong large-scale ascendmore(?500 ?1) and low clouds (stratocumulus and trade wind cumulus) where the models differ most. Shortwave aerosol indirect forcing is also found to differ significantly among different regimes. Shortwave aerosol indirect forcing in ascending regimes is as large as that in stratocumulus regimes, which indicates that regimes with strong large-scale ascend are as important as stratocumulus regimes in studying AIE. It is further shown that shortwave aerosol indirect forcing over regions with high monthly large-scale surface precipitation rate (> 0.1 mm d?1) contributes the most to the total aerosol indirect forcing (from 64 to nearly 100 %). Results show that the uncertainty in AIE is even larger within specific dynamical regimes than that globally, pointing to the need to reduce the uncertainty in AIE in different dynamical regimes.less
Turbulence patterns and neutrino flavor transitions in high-resolution supernova models
Borriello, Enrico; Mirizzi, Alessandro [II. Institut fr Theoretische Physik, Universitt Hamburg, Luruper Chaussee 149, D-22761 Hamburg (Germany); Chakraborty, Sovan [Max-Planck-Institut fr Physik (Werner-Heisenberg-Institut), Fhringer Ring 6, D-80805 Mnchen (Germany); Janka, Hans-Thomas [Max-Planck-Institut fr Astrophysik, Karl-Schwarzschild-Str. 1, 85748 Garching (Germany); Lisi, Eligio, E-mail: enrico.borriello@desy.de, E-mail: sovan@mppmu.mpg.de, E-mail: thj@mpa-garching.mpg.de, E-mail: eligio.lisi@ba.infn.it, E-mail: alessandro.mirizzi@desy.de [INFNSezione di Bari, Via Orabona 4, 70126 Bari (Italy)
2014-11-01
During the shock-wave propagation in a core-collapse supernova (SN), matter turbulence may affect neutrino flavor conversion probabilities. Such effects have been usually studied by adding parametrized small-scale random fluctuations (with arbitrary amplitude) on top of coarse, spherically symmetric matter density profiles. Recently, however, two-dimensional (2D) SN models have reached a space resolution high enough to directly trace anisotropic density profiles, down to scales smaller than the typical neutrino oscillation length. In this context, we analyze the statistical properties of a large set of SN matter density profiles obtained in a high-resolution 2D simulation, focusing on a post-bounce time (2 s) suited to study shock-wave effects on neutrino propagation on scales as small as O(100) km and possibly below. We clearly find the imprint of a broken (Kolmogorov-Kraichnan) power-law structure, as generically expected in 2D turbulence spectra. We then compute the flavor evolution of SN neutrinos along representative realizations of the turbulent matter density profiles, and observe no or modest damping of the neutrino crossing probabilities on their way through the shock wave. In order to check the effect of possibly unresolved fluctuations at scales below O(100) km, we also apply a randomization procedure anchored to the power spectrum calculated from the simulation, and find consistent results within 1? fluctuations. These results show the importance of anchoring turbulence effects on SN neutrinos to realistic, fine-grained SN models.
Using a dynamic point-source percolation model to simulate bubble growth.
Zimmerman, Jonathan A.; Zeigler, David A.; Cowgill, Donald F.
2004-05-01
Accurate modeling of nucleation, growth and clustering of helium bubbles within metal tritide alloys is of high scientific and technological importance. Of interest is the ability to predict both the distribution of these bubbles and the manner in which these bubbles interact at a critical concentration of helium-to-metal atoms to produce an accelerated release of helium gas. One technique that has been used in the past to model these materials, and again revisited in this research, is percolation theory. Previous efforts have used classical percolation theory to qualitatively and quantitatively model the behavior of interstitial helium atoms in a metal tritide lattice; however, higher fidelity models are needed to predict the distribution of helium bubbles and include features that capture the underlying physical mechanisms present in these materials. In this work, we enhance classical percolation theory by developing the dynamic point-source percolation model. This model alters the traditionally binary character of site occupation probabilities by enabling them to vary depending on proximity to existing occupied sites, i.e. nucleated bubbles. This revised model produces characteristics for one and two dimensional systems that are extremely comparable with measurements from three dimensional physical samples. Future directions for continued development of the dynamic model are also outlined.
Li, Ben; He, Feng; Ouyang, Jiting; Duan, Xiaoxi
2015-12-15
Simulation work is very important for understanding the formation of self-organized discharge patterns. Previous works have witnessed different models derived from other systems for simulation of discharge pattern, but most of these models are complicated and time-consuming. In this paper, we introduce a convenient phenomenological dynamic model based on the basic dynamic process of glow discharge and the voltage transfer curve (VTC) to study the dielectric barrier glow discharge (DBGD) pattern. VTC is an important characteristic of DBGD, which plots the change of wall voltage after a discharge as a function of the initial total gap voltage. In the modeling, the combined effect of the discharge conditions is included in VTC, and the activation-inhibition effect is expressed by a spatial interaction term. Besides, the model reduces the dimensionality of the system by just considering the integration effect of current flow. All these greatly facilitate the construction of this model. Numerical simulations turn out to be in good accordance with our previous fluid modeling and experimental result.
Modeling dynamic stall on wind turbine blades under rotationally augmented flow fields
Guntur, S.; Schreck, S.; Sorensen, N. N.; Bergami, L.
2015-04-22
It is well known that airfoils under unsteady flow conditions with a periodically varying angle of attack exhibit aerodynamic characteristics different from those under steady flow conditions, a phenomenon commonly known as dynamic stall. It is also well known that the steady aerodynamic characteristics of airfoils in the inboard region of a rotating blade differ from those under steady two-dimensional (2D) flow conditions, a phenomenon commonly known as rotational augmentation. This paper presents an investigation of these two phenomena together in the inboard parts of wind turbine blades. This analysis is carried out using data from three sources: (1) the National Renewable Energy Laboratory’s Unsteady Aerodynamics Experiment Phase VI experimental data, including constant as well as continuously pitching blade conditions during axial operation, (2) data from unsteady Delayed Detached Eddy Simulations (DDES) carried out using the Technical University of Denmark’s in-house flow solver Ellipsys3D, and (3) data from a simplified model based on the blade element momentum method with a dynamic stall subroutine that uses rotationally augmented steady-state polars obtained from steady Phase VI experimental sequences, instead of the traditional 2D nonrotating data. The aim of this work is twofold. First, the blade loads estimated by the DDES simulations are compared to three select cases of the N sequence experimental data, which serves as a validation of the DDES method. Results show reasonable agreement between the two data in two out of three cases studied. Second, the dynamic time series of the lift and the moment polars obtained from the experiments are compared to those from the dynamic stall subroutine that uses the rotationally augmented steady polars. This allowed the differences between the stall phenomenon on the inboard parts of harmonically pitching blades on a rotating wind turbine and the classic dynamic stall representation in 2D flow to be
Elizondo, Marcelo A.; Tuffner, Francis K.; Schneider, Kevin P.
2016-01-01
Unlike transmission systems, distribution feeders in North America operate under unbalanced conditions at all times, and generally have a single strong voltage source. When a distribution feeder is connected to a strong substation source, the system is dynamically very stable, even for large transients. However if a distribution feeder, or part of the feeder, is separated from the substation and begins to operate as an islanded microgrid, transient dynamics become more of an issue. To assess the impact of transient dynamics at the distribution level, it is not appropriate to use traditional transmission solvers, which generally assume transposed lines and balanced loads. Full electromagnetic solvers capture a high level of detail, but it is difficult to model large systems because of the required detail. This paper proposes an electromechanical transient model of synchronous machine for distribution-level modeling and microgrids. This approach includes not only the machine model, but also its interface with an unbalanced network solver, and a powerflow method to solve unbalanced conditions without a strong reference bus. The presented method is validated against a full electromagnetic transient simulation.
Herron, Jeffrey A.; Scaranto, Jessica; Ferrin, Peter A.; Li, Sha; Mavrikakis, Manos
2014-12-05
We present a first-principles, self-consistent periodic density functional theory (PW91-GGA) study of formic acid (HCOOH) decomposition on model (111) and (100) facets of eight fcc metals (Au, Ag, Cu, Pt, Pd, Ni, Ir, and Rh) and (0001) facets of four hcp (Co, Os, Ru, and Re) metals. The calculated binding energies of key formic acid decomposition intermediates including formate (HCOO), carboxyl (COOH), carbon monoxide (CO), water (H2O), carbon dioxide (CO2), hydroxyl (OH), carbon (C), oxygen (O), and hydrogen (H; H2) are presented. Using these energetics, we develop thermochemical potential energy diagrams for both the carboxyl-mediated and the formate-mediated dehydrogenation mechanisms on each surface. We evaluate the relative stability of COOH, HCOO, and other isomeric intermediates (i.e., CO + OH, CO2 + H, CO + O + H) on these surfaces. These results provide insights into formic acid decomposition selectivity (dehydrogenation versus dehydration), and in conjunction with calculated vibrational frequency modes, the results can assist with the experimental search for the elusive carboxyl (COOH) surface intermediate. Results are compared against experimental reports in the literature.
Sankar, I. V. Chatterjee, Ashok
2014-04-24
The two-dimensional extended Holstein-Hubbard model (EHH) has been considered at strong correlation regime in the non-half-filled band case to understand the self-trapping transition of electrons in strongly correlated electron system. We have used the method of optimized canonical transformations to transform an EHH model into an effective extended Hubbard (EEH) model. In the strong on-site correlation limit an EH model can be transformed into a t-J model which is finally solved using Hartree-Fock approximation (HFA). We found that, for non-half-filled band case, the transition is abrupt in the adiabatic region whereas it is continuous in the anti-adiabatic region.
Joubert-Doriol, Loc; Ryabinkin, Ilya G.; Izmaylov, Artur F.; Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6
2013-12-21
In molecular systems containing conical intersections (CIs), a nontrivial geometric phase (GP) appears in the nuclear and electronic wave functions in the adiabatic representation. We study GP effects in nuclear dynamics of an N-dimensional linear vibronic coupling (LVC) model. The main impact of GP on low-energy nuclear dynamics is reduction of population transfer between the local minima of the LVC lower energy surface. For the LVC model, we proposed an isometric coordinate transformation that confines non-adiabatic effects within a two-dimensional subsystem interacting with an N ? 2 dimensional environment. Since environmental modes do not couple electronic states, all GP effects originate from nuclear dynamics within the subsystem. We explored when the GP affects nuclear dynamics of the isolated subsystem, and how the subsystem-environment interaction can interfere with GP effects. Comparing quantum dynamics with and without GP allowed us to devise simple rules to determine significance of the GP for nuclear dynamics in this model.
A moist aquaplanet variant of the HeldSuarez test for atmospheric model dynamical cores
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Thatcher, D. R.; Jablonowski, C.
2015-09-29
A moist idealized test case (MITC) for atmospheric model dynamical cores is presented. The MITC is based on the HeldSuarez (HS) test that was developed for dry simulations on a flat Earth and replaces the full physical parameterization package with a Newtonian temperature relaxation and Rayleigh damping of the low-level winds. This new variant of the HS test includes moisture and thereby sheds light on the non-linear dynamics-physics moisture feedbacks without the complexity of full physics parameterization packages. In particular, it adds simplified moist processes to the HS forcing to model large-scale condensation, boundary layer mixing, and the exchange ofmorelatent and sensible heat between the atmospheric surface and an ocean-covered planet. Using a variety of dynamical cores of NCAR's Community Atmosphere Model (CAM), this paper demonstrates that the inclusion of the moist idealized physics package leads to climatic states that closely resemble aquaplanet simulations with complex physical parameterizations. This establishes that the MITC approach generates reasonable atmospheric circulations and can be used for a broad range of scientific investigations. This paper provides examples of two application areas. First, the test case reveals the characteristics of the physics-dynamics coupling technique and reproduces coupling issues seen in full-physics simulations. In particular, it is shown that sudden adjustments of the prognostic fields due to moist physics tendencies can trigger undesirable large-scale gravity waves, which can be remedied by a more gradual application of the physical forcing. Second, the moist idealized test case can be used to intercompare dynamical cores. These examples demonstrate the versatility of the MITC approach and suggestions are made for further application areas. The new moist variant of the HS test can be considered a test case of intermediate complexity.less
Broader source: Energy.gov [DOE]
AC Transit (or the Alameda-Contra Costa Transit District) is based in Oakland, California, and provides transportation services to the East Bay of San Francisco. The 360-square-mile service area includes 13 cities and adjacent unincorporated areas in Alameda and Contra Costa counties. AC Transit's approximately 638 vehicles serve more than 65 million annual passengers.
A dynamical model for condition monitoring and fault diagnostics of spur gears
Paya, B.; Esat, I.; Badi, M.N.M.
1996-12-31
The symptoms of condition monitoring and fault diagnostics of machinery based on the dynamic modelling of spur gears are discussed in this paper. The mathematical model presented in the earlier work, assumes two degree of freedom for each gear and the rotor, and also incorporates a varying gear tooth stiffness. This system is assumed to be in good condition (i.e. no fault present). The results obtained from this analytical model are compared with the ones obtained from an experimental model gearbox. This experimental gearbox consists of two meshing spur gears driven by an electric motor. The comparison of the results are encouraging as fundamental (dominant) frequencies of the analytical results correlates very closely to the experimental ones. It is shown that certain vibration frequency of a real gearbox such as the tooth meshing frequencies can be achieved from its mathematical model.
Control method and system for hydraulic machines employing a dynamic joint motion model
Danko, George
2011-11-22
A control method and system for controlling a hydraulically actuated mechanical arm to perform a task, the mechanical arm optionally being a hydraulically actuated excavator arm. The method can include determining a dynamic model of the motion of the hydraulic arm for each hydraulic arm link by relating the input signal vector for each respective link to the output signal vector for the same link. Also the method can include determining an error signal for each link as the weighted sum of the differences between a measured position and a reference position and between the time derivatives of the measured position and the time derivatives of the reference position for each respective link. The weights used in the determination of the error signal can be determined from the constant coefficients of the dynamic model. The error signal can be applied in a closed negative feedback control loop to diminish or eliminate the error signal for each respective link.
LDRD final report : mesoscale modeling of dynamic loading of heterogeneous materials.
Robbins, Joshua; Dingreville, Remi Philippe Michel; Voth, Thomas Eugene; Furnish, Michael David
2013-12-01
Material response to dynamic loading is often dominated by microstructure (grain structure, porosity, inclusions, defects). An example critically important to Sandia's mission is dynamic strength of polycrystalline metals where heterogeneities lead to localization of deformation and loss of shear strength. Microstructural effects are of broad importance to the scientific community and several institutions within DoD and DOE; however, current models rely on inaccurate assumptions about mechanisms at the sub-continuum or mesoscale. Consequently, there is a critical need for accurate and robust methods for modeling heterogeneous material response at this lower length scale. This report summarizes work performed as part of an LDRD effort (FY11 to FY13; project number 151364) to meet these needs.
Depinning transition of a domain wall in ferromagnetic films
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Xi, Bin; Luo, Meng -Bo; Vinokur, Valerii M.; Hu, Xiao
2015-09-14
We report first principle numerical study of domain wall (DW) depinning in two-dimensional magnetic film, which is modeled by 2D random-field Ising system with the dipole-dipole interaction. We observe non-conventional activation-type motion of DW and reveal the fractal structure of DW near the depinning transition. We determine scaling functions describing critical dynamics near the transition and obtain universal exponents establishing connection between thermal softening of pinning potential and critical dynamics. In addition, we observe that tuning the strength of the dipole-dipole interaction switches DW dynamics between two different universality classes, corresponding to two distinct dynamic regimes characterized by non-Arrhenius andmore » conventional Arrhenius-type DW motions.« less
Depinning transition of a domain wall in ferromagnetic films
Xi, Bin; Luo, Meng -Bo; Vinokur, Valerii M.; Hu, Xiao
2015-09-14
We report first principle numerical study of domain wall (DW) depinning in two-dimensional magnetic film, which is modeled by 2D random-field Ising system with the dipole-dipole interaction. We observe non-conventional activation-type motion of DW and reveal the fractal structure of DW near the depinning transition. We determine scaling functions describing critical dynamics near the transition and obtain universal exponents establishing connection between thermal softening of pinning potential and critical dynamics. In addition, we observe that tuning the strength of the dipole-dipole interaction switches DW dynamics between two different universality classes, corresponding to two distinct dynamic regimes characterized by non-Arrhenius and conventional Arrhenius-type DW motions.
Bossard, J.A.; Peck, R.E.; Schmidt, D.K.
1993-03-01
The development of an advanced dynamic model for aeroelastic hypersonic vehicles powered by air breathing engines requires an adequate engine model. This report provides a discussion of some of the more important features of supersonic combustion and their relevance to the analysis and design of supersonic ramjet engines. Of particular interest are those aspects of combustion that impact the control of the process. Furthermore, the report summarizes efforts to enhance the aeropropulsive/aeroelastic dynamic model developed at the Aerospace Research Center of Arizona State University by focusing on combustion and improved modeling of this flow. The expanded supersonic combustor model described here has the capability to model the effects of friction, area change, and mass addition, in addition to the heat addition process. A comparison is made of the results from four cases: (1) heat addition only; (2) heat addition plus friction; (3) heat addition, friction, and area reduction, and (4) heat addition, friction, area reduction, and mass addition. The relative impact of these effects on the Mach number, static temperature, and static pressure distributions within the combustor are then shown. Finally, the effects of frozen versus equilibrium flow conditions within the exhaust plume is discussed.
A Nonlocal Peridynamic Plasticity Model for the Dynamic Flow and Fracture of Concrete.
Vogler, Tracy; Lammi, Christopher James
2014-10-01
A nonlocal, ordinary peridynamic constitutive model is formulated to numerically simulate the pressure-dependent flow and fracture of heterogeneous, quasi-brittle ma- terials, such as concrete. Classical mechanics and traditional computational modeling methods do not accurately model the distributed fracture observed within this family of materials. The peridynamic horizon, or range of influence, provides a characteristic length to the continuum and limits localization of fracture. Scaling laws are derived to relate the parameters of peridynamic constitutive model to the parameters of the classical Drucker-Prager plasticity model. Thermodynamic analysis of associated and non-associated plastic flow is performed. An implicit integration algorithm is formu- lated to calculate the accumulated plastic bond extension and force state. The gov- erning equations are linearized and the simulation of the quasi-static compression of a cylinder is compared to the classical theory. A dissipation-based peridynamic bond failure criteria is implemented to model fracture and the splitting of a concrete cylinder is numerically simulated. Finally, calculation of the impact and spallation of a con- crete structure is performed to assess the suitability of the material and failure models for simulating concrete during dynamic loadings. The peridynamic model is found to accurately simulate the inelastic deformation and fracture behavior of concrete during compression, splitting, and dynamically induced spall. The work expands the types of materials that can be modeled using peridynamics. A multi-scale methodology for simulating concrete to be used in conjunction with the plasticity model is presented. The work was funded by LDRD 158806.
Dynamic and impact contact mechanics of geologic materials: Grain-scale experiments and modeling
Cole, David M.; Hopkins, Mark A.; Ketcham, Stephen A.
2013-06-18
High fidelity treatments of the generation and propagation of seismic waves in naturally occurring granular materials is becoming more practical given recent advancements in our ability to model complex particle shapes and their mechanical interaction. Of particular interest are the grain-scale processes that are activated by impact events and the characteristics of force transmission through grain contacts. To address this issue, we have developed a physics based approach that involves laboratory experiments to quantify the dynamic contact and impact behavior of granular materials and incorporation of the observed behavior indiscrete element models. The dynamic experiments do not involve particle damage and emphasis is placed on measured values of contact stiffness and frictional loss. The normal stiffness observed in dynamic contact experiments at low frequencies (e.g., 10 Hz) are shown to be in good agreement with quasistatic experiments on quartz sand. The results of impact experiments - which involve moderate to extensive levels of particle damage - are presented for several types of naturally occurring granular materials (several quartz sands, magnesite and calcium carbonate ooids). Implementation of the experimental findings in discrete element models is discussed and the results of impact simulations involving up to 5 Multiplication-Sign 105 grains are presented.
Development of a plant-wide dynamic model of an integrated gasification combined cycle (IGCC) plant
Bhattacharyya, D.; Turton, R.; Zitney, S.
2009-01-01
In this presentation, development of a plant-wide dynamic model of an advanced Integrated Gasification Combined Cycle (IGCC) plant with CO2 capture will be discussed. The IGCC reference plant generates 640 MWe of net power using Illinois No.6 coal as the feed. The plant includes an entrained, downflow, General Electric Energy (GEE) gasifier with a radiant syngas cooler (RSC), a two-stage water gas shift (WGS) conversion process, and two advanced 'F' class combustion turbines partially integrated with an elevated-pressure air separation unit (ASU). A subcritical steam cycle is considered for heat recovery steam generation. Syngas is selectively cleaned by a SELEXOL acid gas removal (AGR) process. Sulfur is recovered using a two-train Claus unit with tail gas recycle to the AGR. A multistage intercooled compressor is used for compressing CO2 to the pressure required for sequestration. Using Illinois No.6 coal, the reference plant generates 640 MWe of net power. The plant-wide steady-state and dynamic IGCC simulations have been generated using the Aspen Plus{reg_sign} and Aspen Plus Dynamics{reg_sign} process simulators, respectively. The model is generated based on the Case 2 IGCC configuration detailed in the study available in the NETL website1. The GEE gasifier is represented with a restricted equilibrium reactor model where the temperature approach to equilibrium for individual reactions can be modified based on the experimental data. In this radiant-only configuration, the syngas from the Radiant Syngas Cooler (RSC) is quenched in a scrubber. The blackwater from the scrubber bottom is further cleaned in the blackwater treatment plant. The cleaned water is returned back to the scrubber and also used for slurry preparation. The acid gas from the sour water stripper (SWS) is sent to the Claus plant. The syngas from the scrubber passes through a sour shift process. The WGS reactors are modeled as adiabatic plug flow reactors with rigorous kinetics based on the mid
A spectral transform dynamical core option within the Community Atmosphere Model (CAM4)
Evans, Katherine J; Mahajan, Salil; Branstetter, Marcia L; McClean, Julie L.; Caron, Julie M.; Maltrud, Matthew E.; Hack, James J; Bader, David C; Neale, Rich
2014-01-01
A spectral transform dynamical core with an 85 spectral truncation resolution (T85) within the Community Atmosphere Model (CAM), version 4, is evaluated within the recently released Community Earth System Model, version 1.0 (CESM) global climate model. The spectral dynamical core option provides a well-known base within the climate model community from which to assess climate behavior and statistics, and its relative computational efficiency for smaller computing platforms allows it to be extended to perform climate length simulations using high-resolution configurations in the near term. To establish the characteristics of the CAM4 T85, an ensemble of simulations covering the present day observational period using forced sea surface temperatures and prescribed sea-ice extent are evaluated. Overall, the T85 ensemble attributes and biases are similar to a companion ensemble of simulations using the one degree finite volume (FV1) dynamical core, relative to observed and model derived datasets. Notable improvements with T85 compared to FV1 include the representation of wintertime Arctic sea level pressure and summer precipitation over the Western Indian subcontinent. The mean and spatial patterns of the land surface temperature trends over the AMIP period are generally well simulated with the T85 ensemble relative to observations, however the model is not able to capture the extent nor magnitude of changes in temperature extremes over the boreal summer, where the changes are most dramatic. Biases in the wintertime Arctic surface temperature and annual mean surface stress fields persist with T85 as with the CAM3 version of T85.
Advanced High-Temperature Reactor Dynamic System Model Development: April 2012 Status
Qualls, A L; Cetiner, M S; Wilson, Jr, T L
2012-04-30
The Advanced High-Temperature Reactor (AHTR) is a large-output fluoride-salt-cooled high-temperature reactor (FHR). An early-phase preconceptual design of a 1500 MW(e) power plant was developed in 2011 [Refs. 1 and 2]. An updated version of this plant is shown as Fig. 1. FHRs feature low-pressure liquid fluoride salt cooling, coated-particle fuel, a high-temperature power cycle, and fully passive decay heat rejection. The AHTR is designed to be a “walk away” reactor that requires no action to prevent large off-site releases following even severe reactor accidents. This report describes the development of dynamic system models used to further the AHTR design toward that goal. These models predict system response during warmup, startup, normal operation, and limited off-normal operating conditions. Severe accidents that include a loss-of-fluid inventory are not currently modeled. The scope of the models is limited to the plant power system, including the reactor, the primary and intermediate heat transport systems, the power conversion system, and safety-related or auxiliary heat removal systems. The primary coolant system, the intermediate heat transport system and the reactor building structure surrounding them are shown in Fig. 2. These systems are modeled in the most detail because the passive interaction of the primary system with the surrounding structure and heat removal systems, and ultimately the environment, protects the reactor fuel and the vessel from damage during severe reactor transients. The reactor silo also plays an important role during system warmup. The dynamic system modeling tools predict system performance and response. The goal is to accurately predict temperatures and pressures within the primary, intermediate, and power conversion systems and to study the impacts of design changes on those responses. The models are design tools and are not intended to be used in reactor qualification. The important details to capture in the primary
A graphical interface based model for wind turbine drive train dynamics
Manwell, J.F.; McGowan, J.G.; Abdulwahid, U.; Rogers, A.; McNiff, B.
1996-12-31
This paper presents a summary of a wind turbine drive train dynamics code that has been under development at the University of Massachusetts, under National Renewable Energy Laboratory (NREL) support. The code is intended to be used to assist in the proper design and selection of drive train components. This work summarizes the development of the equations of motion for the model, and discusses the method of solution. In addition, a number of comparisons with analytical solutions and experimental field data are given. The summary includes conclusions and suggestions for future work on the model. 13 refs., 10 figs.
DOE Data Explorer [Office of Scientific and Technical Information (OSTI)]
Users of the VEMAP Portal can access input files of numerical data that include monthly and daily files of geographic data, soil and site files, scenario files, etc. Model results from Phase I, the Equilibrium Response datasets, are available through the NCAR anonymous FTP site at http://www.cgd.ucar.edu/vemap/vresults.html. Phase II, Transient Dynamics, include climate datasets, models results, and analysis tools. Many supplemental files are also available from the main data page at http://www.cgd.ucar.edu/vemap/datasets.html.
DOE Data Explorer [Office of Scientific and Technical Information (OSTI)]
The Vegetation-Ecosystem Modeling and Analysis Project (VEMAP) was a large, collaborative, multi-agency program to simulate and understand ecosystem dynamics for the continental U.S. The project involved the development of common data sets for model input including a high-resolution topographically-adjusted climate history of the U.S. from 1895-1993 on a 0.5? grid, with soils and vegetation cover. The vegetation cover data set includes a detailed agricultural data base based on USDA statistics and remote sensing, as well as natural vegetation (also derived from satellite imagery). Two principal model experiments were run. First, a series of ecosystem models were run from 1895 to 1993 to simulate current ecosystem biogeochemistry. Second, these same models were integrated forward using the output from two climate system models (CCC (Canadian Climate Centre) and Hadley Centre models) using climate results translated into the VEMAP grid and re-adjusted for high-resolution topography for the simulated period 1994-2100.[Quoted from http://www.cgd.ucar.edu/vemap/findings.html] The VEMAP Data Portal is a central collection of files maintained and serviced by the NCAR Data Group. These files (the VEMAP Community Datasets) represent a complete and current collection of VEMAP data files. All data files available through the Data Portal have undergone extensive quality assurance.[Taken from http://www.cgd.ucar.edu/vemap/datasets.html] Users of the VEMAP Portal can access input files of numerical data that include monthly and daily files of geographic data, soil and site files, scenario files, etc. Model results from Phase I, the Equilibrium Response datasets, are available through the NCAR anonymous FTP site at http://www.cgd.ucar.edu/vemap/vresults.html. Phase II, Transient Dynamics, include climate datasets, models results, and analysis tools. Many supplemental files are also available from the main data page at http://www.cgd.ucar.edu/vemap/datasets.html.
Ren, X. D. Liu, R.; Zheng, L. M.; Ren, Y. P.; Hu, Z. Z.; He, H.
2015-10-05
To have a clear insight into nanocrystal growth from graphite to diamond upon high energy pulsed laser irradiation of graphite suspension, synthesis of ultrafine nanocrystalline diamonds with laser energy set up from 0.3 J to 12 J, repetition rate of 10 Hz has been studied. The method allows synthesizing ultrafine nanocrystalline particles continuously at the ambient temperature and normal pressure. The particle size is shown independent of laser energy, which is ultrafine and ranges in 2–6 nm. The theoretical grown size of nano-diamonds is found in well agreement with the experiment results. Four kinds of production were found: nano-diamond, spherical carbon nano-particles, flocculent amorphous carbon, and graphene nano-ribbon rolls. A solid-vapor-plasma-liquid coexistence model describing phase transition from graphite to diamond induced by nanosecond laser processing was proposed. Graphene nano-ribbon rolls might be the intermediate phase in the conversion from graphite to diamond.
Modeling void formation dynamics in fibrous porous media with the lattice Boltzmann method
Spaid, M.A.A.; Phelan, F.R. Jr.
1997-12-31
A novel technique for simulating multiphase fluid flow in the microstructure of a fiber preform is developed, which has the capability of capturing the dynamics of void formation. The model is based on the lattice Boltzmann (LB) method -- a relatively new numerical technique which has rapidly emerged as a powerful tool for simulating multiphase fluid mechanics. The primary benefit of the lattice Boltzmann method is the ability to robustly model the interface between two immiscible fluids without the need for a complex interface tracking algorithm. In a previous paper, it was demonstrated that the lattice Boltzmann method may be modified to solve the Stokes/Brinkman formulation for flow in heterogeneous porous media. Multiphase infiltration of the fiber microstructure is modeled by combining the Stokes/Brinkman LB method, with the multiphase LB algorithm described by Shan and Chen. Numerical results are presented which compare void formation dynamics as a function of the nominal porosity for a model fiber microstructure. In addition, unsaturated permeabilities obtained from the numerical simulations are compared to saturated results for flow in the model porous microstructure.
On the existence of vapor-liquid phase transition in dusty plasmas
Kundu, M.; Sen, A.; Ganesh, R. [Institute for Plasma Research, Bhat, Gandhinagar 382 428, Gujarat (India); Avinash, K. [Department of Physics and Astrophysics, University of Delhi, Delhi 110007 (India)
2014-10-15
The phenomenon of phase transition in a dusty-plasma system (DPS) has attracted some attention in the past. Earlier Farouki and Hamaguchi [J. Chem. Phys. 101, 9876 (1994)] have demonstrated the existence of a liquid to solid transition in DPS where the dust particles interact through a Yukawa potential. However, the question of the existence of a vapor-liquid (VL) transition in such a system remains unanswered and relatively unexplored so far. We have investigated this problem by performing extensive molecular dynamics simulations which show that the VL transition does not have a critical curve in the pressure versus volume diagram for a large range of the Yukawa screening parameter ? and the Coulomb coupling parameter ?. Thus, the VL phase transition is found to be super-critical, meaning that this transition is continuous in the dusty plasma model given by Farouki and Hamaguchi. We provide an approximate analytic explanation of this finding by means of a simple model calculation.