Interpersonal Dynamics in a Simulated Prison: A Methodological Analysis
ERIC Educational Resources Information Center
Banuazizi, Ali; Movahedi, Siamak
1975-01-01
A critical overview is presented of the Stanford Prison Experiment, conducted by Zimbardo and his coinvestigators in which they attempted a structural analysis of the problems of imprisonment. Key assumptions are questioned, primarily on methodological grounds, which casts doubts on the plausibility of the experimenters' final causal inferences.…
Akimov, Alexey V
2016-06-30
The "methodology discovery" library for quantum and classical dynamics simulations is presented. One of the major foci of the code is on nonadiabatic molecular dynamics simulations with model and atomistic Hamiltonians treated on the same footing. The essential aspects of the methodology, design philosophy, and implementation are discussed. The code capabilities are demonstrated on a number of model and atomistic test cases. It is demonstrated how the library can be used to study methodologies for quantum and classical dynamics, as well as a tool for performing detailed atomistic studies of nonadiabatic processes in molecular systems. The source code and additional information are available on the Web at http://www.acsu.buffalo.edu/~alexeyak/libra/index.html. © 2016 Wiley Periodicals, Inc. PMID:27016373
A measurement methodology for dynamic angle of sight errors in hardware-in-the-loop simulation
NASA Astrophysics Data System (ADS)
Zhang, Wen-pan; Wu, Jun-hui; Gan, Lin; Zhao, Hong-peng; Liang, Wei-wei
2015-10-01
In order to precisely measure dynamic angle of sight for hardware-in-the-loop simulation, a dynamic measurement methodology was established and a set of measurement system was built. The errors and drifts, such as synchronization delay, CCD measurement error and drift, laser spot error on diffuse reflection plane and optics axis drift of laser, were measured and analyzed. First, by analyzing and measuring synchronization time between laser and time of controlling data, an error control method was devised and lowered synchronization delay to 21μs. Then, the relationship between CCD device and laser spot position was calibrated precisely and fitted by two-dimension surface fitting. CCD measurement error and drift were controlled below 0.26mrad. Next, angular resolution was calculated, and laser spot error on diffuse reflection plane was estimated to be 0.065mrad. Finally, optics axis drift of laser was analyzed and measured which did not exceed 0.06mrad. The measurement results indicate that the maximum of errors and drifts of the measurement methodology is less than 0.275mrad. The methodology can satisfy the measurement on dynamic angle of sight of higher precision and lager scale.
NASA Astrophysics Data System (ADS)
Mántaras, Daniel A.; Luque, Pablo; Nava, Javier A.; Riva, Paolo; Girón, Pablo; Compadre, Diego; Ferran, Jordi
2013-10-01
A key factor to understand the vehicle dynamic behaviour is to know as accurately as possible the interaction that occurs between the tyre and the road, since it depends on many factors that influence the dynamic response of the vehicle. This paper aims to develop a methodology in order to characterise the tyre-road behaviour, applying it to obtain the tyre-road grip coefficient. This methodology is based on the use of dynamic simulation of a virtual model, integrated into a genetic algorithm that identifies the tyre-road friction coefficient in order to adjust the response obtained by simulation to real data. The numerical model was developed in collaboration with SEAT Technical Centre and it was implemented in multibody dynamic simulation software Adams®, from MSC®.
Dynamic optimization methodology based on subgrid-scale dissipation for large eddy simulation
NASA Astrophysics Data System (ADS)
Yu, Changping; Xiao, Zuoli; Li, Xinliang
2016-01-01
A dynamic procedure based on subgrid-scale dissipation is proposed for large eddy simulation of turbulent flows. In the new method, the model coefficients are determined by minimizing the square error of the resolved dissipation rate based on the Germano identity. A dynamic two-term mixed model is tested and evaluated both a priori and a posteriori in simulations of homogeneous and isotropic turbulence. The new dynamic procedure proves to be more effective to optimize the model coefficients as compared with traditional method. The corresponding dynamic mixed model can predict the physical quantities more accurately than traditional dynamic mixed model.
NASA Astrophysics Data System (ADS)
Tremblay, J. P.; Viau, C. R.
2009-09-01
The paper describes a methodology for characterizing the signatures of targets for Imaging Infrared (IIR) missiles and generating dynamic missile engagement scenarios using MathWorks tools (primarily MATLAB and Simulink). The over-all objective of this work was to develop high fidelity physics-based simulations of the attack of IIR missiles on targets that are using various types of countermeasures for survivability. While the methodology has been implemented in products used for analyses of both ship and main battle tank protection this paper focuses on the ship application. The methodology involves a multi-step process. First the infrared signatures of the objects are characterized using a graphical tool that enables the user to select individual or groups of surfaces on the objects (targets and countermeasures) and specify their surface temperatures and spectral emissivities. Second, a dynamic IR scene generator creates the scene as viewed by the missile's seeker. Then an imaging IR seeker, using the option of several tracking algorithms, discriminates the target. Finally, the inclusion of dynamic models for missile guidance, aerodynamics and propulsion together with signal propagation enable the closing of the loop in the missile's fly-out. The simulation dynamically computes the distance between each surface and the missile seeker and uses the specified atmospheric attenuation profile to produce a simulated IR image at the seeker. This is processed using several optional tracking algorithms to generate steering signals. This process is repeated every time-step of the simulation and determines the trajectory of the missile and the hit or miss of the missile at engagement completion. The paper includes the following topics: characterizing IR signatures, generating dynamic IR scenes, simulating representative close-loop missile fly-out engagements, evaluating performance and running simulation batches.
A novel simulation methodology merging source-sink dynamics and landscape connectivity
Source-sink dynamics are an emergent property of complex species-landscape interactions. This study explores the patterns of source and sink behavior that become established across a large landscape, using a simulation model for the northern spotted owl (Strix occidentalis cauri...
Validation of computational fluid dynamics methodology used for human upper airway flow simulations.
Mylavarapu, Goutham; Murugappan, Shanmugam; Mihaescu, Mihai; Kalra, Maninder; Khosla, Sid; Gutmark, Ephraim
2009-07-22
An anatomically accurate human upper airway model was constructed from multiple magnetic resonance imaging axial scans. This model was used to conduct detailed Computational Fluid Dynamics (CFD) simulations during expiration, to investigate the fluid flow in the airway regions where obstruction could occur. An identical physical model of the same airway was built using stereo lithography. Pressure and velocity measurements were conducted in the physical model. Both simulations and experiments were performed at a peak expiratory flow rate of 200 L/min. Several different numerical approaches within the FLUENT commercial software framework were used in the simulations; unsteady Large Eddy Simulation (LES), steady Reynolds-Averaged Navier-Stokes (RANS) with two-equation turbulence models (i.e. k-epsilon, standard k-omega, and k-omega Shear Stress Transport (SST)) and with one-equation Spalart-Allmaras model. The CFD predictions of the average wall static pressures at different locations along the airway wall were favorably compared with the experimental data. Among all the approaches, standard k-omega turbulence model resulted in the best agreement with the static pressure measurements, with an average error of approximately 20% over all ports. The highest positive pressures were observed in the retroglossal regions below the epiglottis, while the lowest negative pressures were recorded in the retropalatal region. The latter is a result of the airflow acceleration in the narrow retropalatal region. The largest pressure drop was observed at the tip of the soft palate. This location has the smallest cross section of the airway. The good agreement between the computations and the experimental results suggest that CFD simulations can be used to accurately compute aerodynamic flow characteristics of the upper airway. PMID:19501360
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 states 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.
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
Simulation Enabled Safeguards Assessment Methodology
Robert Bean; Trond Bjornard; Thomas Larson
2007-09-01
It is expected that nuclear energy will be a significant component of future supplies. New facilities, operating under a strengthened international nonproliferation regime will be needed. There is good reason to believe virtual engineering applied to the facility design, as well as to the safeguards system design will reduce total project cost and improve efficiency in the design cycle. Simulation Enabled Safeguards Assessment MEthodology (SESAME) has been developed as a software package to provide this capability for nuclear reprocessing facilities. The software architecture is specifically designed for distributed computing, collaborative design efforts, and modular construction to allow step improvements in functionality. Drag and drop wireframe construction allows the user to select the desired components from a component warehouse, render the system for 3D visualization, and, linked to a set of physics libraries and/or computational codes, conduct process evaluations of the system they have designed.
NASA Astrophysics Data System (ADS)
Burresi, E.; Celino, M.
2012-05-01
A single wurtzite phase of cadmium sulfide cluster is investigated by ab-initio molecular dynamics simulations at different temperatures, ranging from 100 K to 600 K. In this study we propose a possible procedure to characterize the CdS quantum dots system by means of molecular dynamics calculations using a standard Car-Parrinello scheme. In order to ensure the accuracy of the numerical approach, preliminary calculations to test pseudopotentials, cutoff and box size on both single atoms systems and Cd-Cd, S-S, Cd-S dimers have been performed. Calculated binding energies and bond lengths are obtained in good agreement with experimental data. Subsequently, an uncapped CdS cluster with size below 2 nm, 48 atoms of cadmium and 48 atoms of sulfur, in a wurtzite geometry was structurally optimized to minimize internal stresses. The CdS cluster has been carefully characterized structurally at several temperatures up to T = 600 K. At the temperature of 340 K atomic diffusion on the surface allows the onset of a new stable atomic configuration.
A methodology for the assessment of manned flight simulator fidelity
NASA Technical Reports Server (NTRS)
Hess, Ronald A.; Malsbury, Terry N.
1989-01-01
A relatively simple analytical methodology for assessing the fidelity of manned flight simulators for specific vehicles and tasks is offered. The methodology is based upon an application of a structural model of the human pilot, including motion cue effects. In particular, predicted pilot/vehicle dynamic characteristics are obtained with and without simulator limitations. A procedure for selecting model parameters can be implemented, given a probable pilot control strategy. In analyzing a pair of piloting tasks for which flight and simulation data are available, the methodology correctly predicted the existence of simulator fidelity problems. The methodology permitted the analytical evaluation of a change in simulator characteristics and indicated that a major source of the fidelity problems was a visual time delay in the simulation.
Real time flight simulation methodology
NASA Technical Reports Server (NTRS)
Parrish, E. A.; Cook, G.; Mcvey, E. S.
1977-01-01
Substitutional methods for digitization, input signal-dependent integrator approximations, and digital autopilot design were developed. The software framework of a simulator design package is described. Included are subroutines for iterative designs of simulation models and a rudimentary graphics package.
Methodology for Validating Building Energy Analysis Simulations
Judkoff, R.; Wortman, D.; O'Doherty, B.; Burch, J.
2008-04-01
The objective of this report was to develop a validation methodology for building energy analysis simulations, collect high-quality, unambiguous empirical data for validation, and apply the validation methodology to the DOE-2.1, BLAST-2MRT, BLAST-3.0, DEROB-3, DEROB-4, and SUNCAT 2.4 computer programs. This report covers background information, literature survey, validation methodology, comparative studies, analytical verification, empirical validation, comparative evaluation of codes, and conclusions.
River system environmental modeling and simulation methodology
Rao, N.B.
1981-01-01
Several computer models have been built to examine pollution in rivers. However, the current state of the art in this field emphasizes problem solving using specific programs. A general methodology for building and simulating models of river systems is lacking. Thus, the purpose of this research was to develop a methodology which can be used to conceptualize, visualize, construct and analyze using simulation, models of pollution in river systems. The conceptualization and visualization of these models was facilitated through a network representation. The implementation of the models was accomplished using the capabilities of an existing simulation language, GASP V. The methodology also provides data management facilities for model outputs through the use of the Simulation Data Language (SDL), and high quality plotting facilities through the use of the graphics package DISSPLA (Display Integrated Software System and Plotting Language). Using this methodology, a river system is modeled as consisting of certain elements, namely reaches, junctions, dams, reservoirs, withdrawals and pollutant sources. All these elements of the river system are described in a standard form which has been implemented on a computer. This model, when executed, produces spatial and temporal distributions of the pollutants in the river system. Furthermore, these outputs can be stored in a database and used to produce high quality plots. The result of this research is a methodology for building, implementing and examining the results of models of pollution in river systems.
NERVA dynamic analysis methodology, SPRVIB
NASA Technical Reports Server (NTRS)
Vronay, D. F.
1972-01-01
The general dynamic computer code called SPRVIB (Spring Vib) developed in support of the NERVA (nuclear engine for rocket vehicle application) program is described. Using normal mode techniques, the program computes kinematical responses of a structure caused by various combinations of harmonic and elliptic forcing functions or base excitations. Provision is made for a graphical type of force or base excitation input to the structure. A description of the required input format and a listing of the program are presented, along with several examples illustrating the use of the program. SPRVIB is written in FORTRAN 4 computer language for use on the CDC 6600 or the IBM 360/75 computers.
Dhingra, Radhika; Jimenez, Violeta; Chang, Howard H; Gambhir, Manoj; Fu, Joshua S; Liu, Yang; Remais, Justin V
2013-09-01
Poikilothermic disease vectors can respond to altered climates through spatial changes in both population size and phenology. Quantitative descriptors to characterize, analyze and visualize these dynamic responses are lacking, particularly across large spatial domains. In order to demonstrate the value of a spatially explicit, dynamic modeling approach, we assessed spatial changes in the population dynamics of Ixodes scapularis, the Lyme disease vector, using a temperature-forced population model simulated across a grid of 4 × 4 km cells covering the eastern United States, using both modeled (Weather Research and Forecasting (WRF) 3.2.1) baseline/current (2001-2004) and projected (Representative Concentration Pathway (RCP) 4.5 and RCP 8.5; 2057-2059) climate data. Ten dynamic population features (DPFs) were derived from simulated populations and analyzed spatially to characterize the regional population response to current and future climate across the domain. Each DPF under the current climate was assessed for its ability to discriminate observed Lyme disease risk and known vector presence/absence, using data from the US Centers for Disease Control and Prevention. Peak vector population and month of peak vector population were the DPFs that performed best as predictors of current Lyme disease risk. When examined under baseline and projected climate scenarios, the spatial and temporal distributions of DPFs shift and the seasonal cycle of key questing life stages is compressed under some scenarios. Our results demonstrate the utility of spatial characterization, analysis and visualization of dynamic population responses-including altered phenology-of disease vectors to altered climate. PMID:24772388
Integrating shadow casting methodology and thermal simulation
Malkawi, A.; Jabi, W.
1996-10-01
This paper describes an experiment that integrates shadow casting methodology and thermal simulation algorithms developed by the authors. The 3D shadow procedures use a polyhedral representation of solids within a Cartesian space that allows for accurate casting of shadows. The algorithm is also capable of calculating surface areas of polygonal shadows of any arbitrary shape and size. The thermal simulation algorithms--using the Transfer Function Method (TFM)--incorporate the shaded area calculations to better predict solar heat gain from glazing based on transmitted, absorbed, and conducted cooling loads. The paper describes the use of a 3D computer model to illustrate the impact of the pattern and area of shading on the visual and thermal properties of building apertures. The paper discusses the objectives of this experiment, the algorithms used, and their integration. Conclusions and findings are drawn.
METHODOLOGICAL NOTES: Strange attractors in rattleback dynamics
NASA Astrophysics Data System (ADS)
Borisov, Aleksei V.; Mamaev, Ivan S.
2003-04-01
This review is dedicated to the dynamics of the rattleback, a phenomenon with curious physical properties that is studied in nonholonomic mechanics. All known analytical results are collected here, and some results of our numerical simulation are presented. In particular, three-dimensional Poincare maps associated with dynamical systems are systematically investigated for the first time. It is shown that the loss of stability of periodic and quasiperiodic solutions, which gives rise to strange attractors, is typical of the three-dimensional maps related to rattleback dynamics. This explains some newly discovered properties of the rattleback related to the transition from regular to chaotic solutions at certain values of the physical parameters.
NASA Technical Reports Server (NTRS)
1987-01-01
The proceedings of the conference are presented. The objective was to provide a forum for the discussion of the structure and status of existing computer programs which are used to simulate the dynamics of a variety of tether applications in space. A major topic was different simulation models and the process of validating them. Guidance on future work in these areas was obtained from a panel discussion; the panel was composed of resource and technical managers and dynamic analysts in the tether field. The conclusions of this panel are also presented.
NASA Technical Reports Server (NTRS)
Joncas, K. P.
1972-01-01
Concepts and techniques for identifying and simulating both the steady state and dynamic characteristics of electrical loads for use during integrated system test and evaluation are discussed. The investigations showed that it is feasible to design and develop interrogation and simulation equipment to perform the desired functions. During the evaluation, actual spacecraft loads were interrogated by stimulating the loads with their normal input voltage and measuring the resultant voltage and current time histories. Elements of the circuits were optimized by an iterative process of selecting element values and comparing the time-domain response of the model with those obtained from the real equipment during interrogation.
Cantera Aerosol Dynamics Simulator
Moffat, Harry
2004-09-01
The Cantera Aerosol Dynamics Simulator (CADS) package is a general library for aerosol modeling to address aerosol general dynamics, including formation from gas phase reactions, surface chemistry (growth and oxidation), bulk particle chemistry, transport by Brownian diffusion, thermophoresis, and diffusiophoresis with linkage to DSMC studies, and thermal radiative transport. The library is based upon Cantera, a C++ Cal Tech code that handles gas phase species transport, reaction, and thermodynamics. The method uses a discontinuous galerkin formulation for the condensation and coagulation operator that conserves particles, elements, and enthalpy up to round-off error. Both O-D and 1-D time dependent applications have been developed with the library. Multiple species in the solid phase are handled as well. The O-D application, called Tdcads (Time Dependent CADS) is distributed with the library. Tdcads can address both constant volume and constant pressure adiabatic homogeneous problems. An extensive set of sample problems for Tdcads is also provided.
Data Systems Dynamic Simulator
NASA Technical Reports Server (NTRS)
Rouff, Christopher; Clark, Melana; Davenport, Bill; Message, Philip
1993-01-01
The Data System Dynamic Simulator (DSDS) is a discrete event simulation tool. It was developed for NASA for the specific purpose of evaluating candidate architectures for data systems of the Space Station era. DSDS provides three methods for meeting this requirement. First, the user has access to a library of standard pre-programmed elements. These elements represent tailorable components of NASA data systems and can be connected in any logical manner. Secondly, DSDS supports the development of additional elements. This allows the more sophisticated DSDS user the option of extending the standard element set. Thirdly, DSDS supports the use of data streams simulation. Data streams is the name given to a technique that ignores packet boundaries, but is sensitive to rate changes. Because rate changes are rare compared to packet arrivals in a typical NASA data system, data stream simulations require a fraction of the CPU run time. Additionally, the data stream technique is considerably more accurate than another commonly-used optimization technique.
Energy Science and Technology Software Center (ESTSC)
2015-10-20
Look-ahead dynamic simulation software system incorporates the high performance parallel computing technologies, significantly reduces the solution time for each transient simulation case, and brings the dynamic simulation analysis into on-line applications to enable more transparency for better reliability and asset utilization. It takes the snapshot of the current power grid status, functions in parallel computing the system dynamic simulation, and outputs the transient response of the power system in real time.
Multibody dynamic simulation of knee contact mechanics
Bei, Yanhong; Fregly, Benjamin J.
2006-01-01
Multibody dynamic musculoskeletal models capable of predicting muscle forces and joint contact pressures simultaneously would be valuable for studying clinical issues related to knee joint degeneration and restoration. Current three-dimensional multi-body knee models are either quasi-static with deformable contact or dynamic with rigid contact. This study proposes a computationally efficient methodology for combining multibody dynamic simulation methods with a deformable contact knee model. The methodology requires preparation of the articular surface geometry, development of efficient methods to calculate distances between contact surfaces, implementation of an efficient contact solver that accounts for the unique characteristics of human joints, and specification of an application programming interface for integration with any multibody dynamic simulation environment. The current implementation accommodates natural or artificial tibiofemoral joint models, small or large strain contact models, and linear or nonlinear material models. Applications are presented for static analysis (via dynamic simulation) of a natural knee model created from MRI and CT data and dynamic simulation of an artificial knee model produced from manufacturer’s CAD data. Small and large strain natural knee static analyses required 1 min of CPU time and predicted similar contact conditions except for peak pressure, which was higher for the large strain model. Linear and nonlinear artificial knee dynamic simulations required 10 min of CPU time and predicted similar contact force and torque but different contact pressures, which were lower for the nonlinear model due to increased contact area. This methodology provides an important step toward the realization of dynamic musculoskeletal models that can predict in vivo knee joint motion and loading simultaneously. PMID:15564115
Simulation Methodology in Nursing Education and Adult Learning Theory
ERIC Educational Resources Information Center
Rutherford-Hemming, Tonya
2012-01-01
Simulation is often used in nursing education as a teaching methodology. Simulation is rooted in adult learning theory. Three learning theories, cognitive, social, and constructivist, explain how learners gain knowledge with simulation experiences. This article takes an in-depth look at each of these three theories as each relates to simulation.…
Molecular dynamics simulations
Alder, B.J.
1985-07-01
The molecular dynamics computer simulation discovery of the slow decay of the velocity autocorrelation function in fluids is briefly reviewed in order to contrast that long time tail with those observed for the stress autocorrelation function in fluids and the velocity autocorrelation function in the Lorentz gas. For a non-localized particle in the Lorentz gas it is made plausible that even if it behaved quantum mechanically its long time tail would be the same as the classical one. The generalization of Fick's law for diffusion for the Lorentz gas, necessary to avoid divergences due to the slow decay of correlations, is presented. For fluids, that generalization has not yet been established, but the region of validity of generalized hydrodynamics is discussed. 20 refs., 5 figs.
Cantera Aerosol Dynamics Simulator
Energy Science and Technology Software Center (ESTSC)
2004-09-01
The Cantera Aerosol Dynamics Simulator (CADS) package is a general library for aerosol modeling to address aerosol general dynamics, including formation from gas phase reactions, surface chemistry (growth and oxidation), bulk particle chemistry, transport by Brownian diffusion, thermophoresis, and diffusiophoresis with linkage to DSMC studies, and thermal radiative transport. The library is based upon Cantera, a C++ Cal Tech code that handles gas phase species transport, reaction, and thermodynamics. The method uses a discontinuous galerkinmore » formulation for the condensation and coagulation operator that conserves particles, elements, and enthalpy up to round-off error. Both O-D and 1-D time dependent applications have been developed with the library. Multiple species in the solid phase are handled as well. The O-D application, called Tdcads (Time Dependent CADS) is distributed with the library. Tdcads can address both constant volume and constant pressure adiabatic homogeneous problems. An extensive set of sample problems for Tdcads is also provided.« less
Application of control theory to dynamic systems simulation
NASA Technical Reports Server (NTRS)
Auslander, D. M.; Spear, R. C.; Young, G. E.
1982-01-01
The application of control theory is applied to dynamic systems simulation. Theory and methodology applicable to controlled ecological life support systems are considered. Spatial effects on system stability, design of control systems with uncertain parameters, and an interactive computing language (PARASOL-II) designed for dynamic system simulation, report quality graphics, data acquisition, and simple real time control are discussed.
Experimental verification of dynamic simulation
NASA Technical Reports Server (NTRS)
Yae, K. Harold; Hwang, Howyoung; Chern, Su-Tai
1989-01-01
The dynamics model here is a backhoe, which is a four degree of freedom manipulator from the dynamics standpoint. Two types of experiment are chosen that can also be simulated by a multibody dynamics simulation program. In the experiment, recorded were the configuration and force histories; that is, velocity and position, and force output and differential pressure change from the hydraulic cylinder, in the time domain. When the experimental force history is used as driving force in the simulation model, the forward dynamics simulation produces a corresponding configuration history. Then, the experimental configuration history is used in the inverse dynamics analysis to generate a corresponding force history. Therefore, two sets of configuration and force histories--one set from experiment, and the other from the simulation that is driven forward and backward with the experimental data--are compared in the time domain. More comparisons are made in regard to the effects of initial conditions, friction, and viscous damping.
Dynamic Shear Modulus of Polymers from Molecular Dynamics Simulations
NASA Astrophysics Data System (ADS)
Byutner, Oleksiy; Smith, Grant
2001-03-01
In this work we describe the methodology for using equilibrium molecular dynamics simulations (MD) simulations to obtain the viscoelastic properties of polymers in the glassy regime. Specifically we show how the time dependent shear stress modulus and frequency dependent complex shear modulus in the high-frequency regime can be determined from the off-diagonal terms of the stress-tensor autocorrelation function obtained from MD trajectories using the Green-Kubo method and appropriate Fourier transforms. In order to test the methodology we have performed MD simulations of a low-molecular-weight polybutadiene system using quantum chemistry based potential functions. Values of the glassy modulus and the maximum loss frequency were found to be in good agreement with experimental data for polybutadiene at 298 K.
Remote manipulator dynamic simulation
NASA Technical Reports Server (NTRS)
Wild, E. C.; Donges, P. K.; Garand, W. A.
1972-01-01
A simulator to generate the real time visual scenes required to perform man in the loop investigations of remote manipulator application and design concepts for the space shuttle is described. The simulated remote manipulator consists of a computed display system that uses a digital computer, the electronic scene generator, an operator's station, and associated interface hardware. A description of the capabilities of the implemented simulation is presented. The mathematical models and programs developed for the simulation are included.
Adaptive LES Methodology for Turbulent Flow Simulations
Oleg V. Vasilyev
2008-06-12
Although turbulent flows are common in the world around us, a solution to the fundamental equations that govern turbulence still eludes the scientific community. Turbulence has often been called one of the last unsolved problem in classical physics, yet it is clear that the need to accurately predict the effect of turbulent flows impacts virtually every field of science and engineering. As an example, a critical step in making modern computational tools useful in designing aircraft is to be able to accurately predict the lift, drag, and other aerodynamic characteristics in numerical simulations in a reasonable amount of time. Simulations that take months to years to complete are much less useful to the design cycle. Much work has been done toward this goal (Lee-Rausch et al. 2003, Jameson 2003) and as cost effective accurate tools for simulating turbulent flows evolve, we will all benefit from new scientific and engineering breakthroughs. The problem of simulating high Reynolds number (Re) turbulent flows of engineering and scientific interest would have been solved with the advent of Direct Numerical Simulation (DNS) techniques if unlimited computing power, memory, and time could be applied to each particular problem. Yet, given the current and near future computational resources that exist and a reasonable limit on the amount of time an engineer or scientist can wait for a result, the DNS technique will not be useful for more than 'unit' problems for the foreseeable future (Moin & Kim 1997, Jimenez & Moin 1991). The high computational cost for the DNS of three dimensional turbulent flows results from the fact that they have eddies of significant energy in a range of scales from the characteristic length scale of the flow all the way down to the Kolmogorov length scale. The actual cost of doing a three dimensional DNS scales as Re{sup 9/4} due to the large disparity in scales that need to be fully resolved. State-of-the-art DNS calculations of isotropic turbulence
NASA Astrophysics Data System (ADS)
Ren, Lei; Zhang, Lin; Tao, Fei; (Luke) Zhang, Xiaolong; Luo, Yongliang; Zhang, Yabin
2012-08-01
Multidisciplinary design of complex products leads to an increasing demand for high performance simulation (HPS) platforms. One great challenge is how to achieve high efficient utilisation of large-scale simulation resources in distributed and heterogeneous environments. This article reports a virtualisation-based methodology to realise a HPS platform. This research is driven by the issues concerning large-scale simulation resources deployment and complex simulation environment construction, efficient and transparent utilisation of fine-grained simulation resources and high reliable simulation with fault tolerance. A framework of virtualisation-based simulation platform (VSIM) is first proposed. Then the article investigates and discusses key approaches in VSIM, including simulation resources modelling, a method to automatically deploying simulation resources for dynamic construction of system environment, and a live migration mechanism in case of faults in run-time simulation. Furthermore, the proposed methodology is applied to a multidisciplinary design system for aircraft virtual prototyping and some experiments are conducted. The experimental results show that the proposed methodology can (1) significantly improve the utilisation of fine-grained simulation resources, (2) result in a great reduction in deployment time and an increased flexibility for simulation environment construction and (3)achieve fault tolerant simulation.
Application of the EPRI simulator qualification methodology to the Bugey simulator
Michel, B.; Le Coq, G. , 75 - Paris )
1990-01-01
The report describes application of the EPRI Simulator Qualification Methodology by EdF to the Bugey Simulator. The major objective of the study was to evaluate the methodology itself rather than to perform a complete qualification of the simulator. The key objectives were achieved. The major advantages and difficulties of using the methodology have been identified and discussed. 7 refs., 21 figs., 11 tabs.
A computer simulator for development of engineering system design methodologies
NASA Technical Reports Server (NTRS)
Padula, S. L.; Sobieszczanski-Sobieski, J.
1987-01-01
A computer program designed to simulate and improve engineering system design methodology is described. The simulator mimics the qualitative behavior and data couplings occurring among the subsystems of a complex engineering system. It eliminates the engineering analyses in the subsystems by replacing them with judiciously chosen analytical functions. With the cost of analysis eliminated, the simulator is used for experimentation with a large variety of candidate algorithms for multilevel design optimization to choose the best ones for the actual application. Thus, the simulator serves as a development tool for multilevel design optimization strategy. The simulator concept, implementation, and status are described and illustrated with examples.
Energy Science and Technology Software Center (ESTSC)
2015-09-14
GridDyn is a part of power grid simulation toolkit. The code is designed using modern object oriented C++ methods utilizing C++11 and recent Boost libraries to ensure compatibility with multiple operating systems and environments.
Dynamic simulation of voltage collapses
Deuse, J.; Stubbe, M. )
1993-08-01
Most of the time the voltage collapse phenomena are studied by means of computer programs designed for the calculation of steady state conditions. But in the real world, the simultaneous occurrences of losses of synchronism, of AVR dynamics or of transformer tap changes call for a full dynamic simulation of voltage phenomena. The present paper shows some examples of dynamic simulations of voltage phenomena using a new general purpose stability program (EUROSTAG), covering in a continuous way the classical fields of transient, mid-term and long-term stability, and also the quasi steady state conditions of a power system.
Methodology for characterizing modeling and discretization uncertainties in computational simulation
ALVIN,KENNETH F.; OBERKAMPF,WILLIAM L.; RUTHERFORD,BRIAN M.; DIEGERT,KATHLEEN V.
2000-03-01
This research effort focuses on methodology for quantifying the effects of model uncertainty and discretization error on computational modeling and simulation. The work is directed towards developing methodologies which treat model form assumptions within an overall framework for uncertainty quantification, for the purpose of developing estimates of total prediction uncertainty. The present effort consists of work in three areas: framework development for sources of uncertainty and error in the modeling and simulation process which impact model structure; model uncertainty assessment and propagation through Bayesian inference methods; and discretization error estimation within the context of non-deterministic analysis.
NASA Astrophysics Data System (ADS)
Grabowski, Wojciech W.
2016-04-01
Formation and growth of cloud and precipitation particles ("cloud microphysics") affect cloud dynamics and such macroscopic cloud field properties as the mean surface rainfall, cloud cover, and liquid/ice water paths. Traditional approaches to investigate the impacts involve parallel simulations with different microphysical schemes or with different scheme parameters (such as the assumed droplet/ice concentration for single-moment bulk schemes or the assumed CCN/IN concentration for double-moment schemes). Such methodologies are not reliable because of the natural variability of a cloud field that is affected by the feedback between cloud microphysics and cloud dynamics. In a nutshell, changing the cloud microphysics leads to a different realization of the cloud-scale flow, and separating dynamical and microphysical impacts is cumbersome. A novel modeling methodology, referred to as the microphysical piggybacking, was recently developed to separate purely microphysical effects from the impact on the dynamics. The main idea is to use two sets of thermodynamic variables driven by two microphysical schemes or by the same scheme with different scheme parameters. One set is coupled to the dynamics and drives the simulation, and the other set piggybacks the simulated flow, that is, it responds to the simulated flow but does not affect it. By switching the sets (i.e., the set driving the simulation becomes the piggybacking one, and vice versa), the impact on the cloud dynamics can be isolated from purely microphysical effects. Application of this methodology to the daytime deep convection development over land based on the observations during the Large-scale Biosphere-Atmosphere (LBA) field project in Amazonia will be discussed applying single-moment and double-moment bulk microphysics schemes. We show that the new methodology documents a small indirect aerosol impact on convective dynamics, and a strong microphysical effect. These results question the postulated strong
Methodology development for evaluation of selective-fidelity rotorcraft simulation
NASA Technical Reports Server (NTRS)
Lewis, William D.; Schrage, D. P.; Prasad, J. V. R.; Wolfe, Daniel
1992-01-01
This paper addressed the initial step toward the goal of establishing performance and handling qualities acceptance criteria for realtime rotorcraft simulators through a planned research effort to quantify the system capabilities of 'selective fidelity' simulators. Within this framework the simulator is then classified based on the required task. The simulator is evaluated by separating the various subsystems (visual, motion, etc.) and applying corresponding fidelity constants based on the specific task. This methodology not only provides an assessment technique, but also provides a technique to determine the required levels of subsystem fidelity for a specific task.
Floating orbital molecular dynamics simulations.
Perlt, Eva; Brüssel, Marc; Kirchner, Barbara
2014-04-21
We introduce an alternative ab initio molecular dynamics simulation as a unification of Hartree-Fock molecular dynamics and the floating orbital approach. The general scheme of the floating orbital molecular dynamics method is presented. Moreover, a simple but sophisticated guess for the orbital centers is provided to reduce the number of electronic structure optimization steps at each molecular dynamics step. The conservation of total energy and angular momentum is investigated in order to validate the floating orbital molecular dynamics approach with and without application of the initial guess. Finally, a water monomer and a water dimer are simulated, and the influence of the orbital floating on certain properties like the dipole moment is investigated. PMID:24600690
Dynamic brain mapping methodology and application.
Itil, T M; Mucci, A; Eralp, E
1991-01-01
Brain mapping has opened important perspectives for the neurophysiological evaluation of patients, for the discrimination of drug effects on the brain and for the study of the relationship between the brain and behavior. Our Dynamic Brain Mapping System is the result of many years of EEG quantification. It was designed as a software-oriented system to favor the largest clinical application and simultaneously stimulate new research objectives. Data collection and analysis procedures are critically important in brain mapping for a good understanding of the results. For clinical use, the maps should answer relevant EEG questions and be interpretable with the consolidated knowledge. Therefore, we have developed a new type of brain mapping technology which is called "Field blending interpolation" mapping offered together with the conventional technology with user-selectable interpolation algorithms. In addition to diagnosis, the use of computer-analyzed EEG and brain mapping can be instrumental in drug monitoring, drug selection and drug discriminations. Prospective studies are, however, required to validate the use of brain mapping in each of these new areas. Spatial analysis is the original goal of brain mapping. The development of a new data collection procedure and analysis will be instrumental in the determination of an adequate time and space resolution. PMID:2010323
A POLLUTION REDUCTION METHODOLOGY FOR CHEMICAL PROCESS SIMULATORS
A pollution minimization methodology was developed for chemical process design using computer simulation. It is based on a pollution balance that at steady state is used to define a pollution index with units of mass of pollution per mass of products. The pollution balance has be...
Tree Modeling and Dynamics Simulation
NASA Astrophysics Data System (ADS)
Tian-shuang, Fu; Yi-bing, Li; Dong-xu, Shen
This paper introduces the theory about tree modeling and dynamic movements simulation in computer graphics. By comparing many methods we choose Geometry-based rendering as our method. The tree is decomposed into branches and leaves, under the rotation and quaternion methods we realize the tree animation and avoid the Gimbals Lock in Euler rotation. We take Orge 3D as render engine, which has good graphics programming ability. By the end we realize the tree modeling and dynamic movements simulation, achieve realistic visual quality with little computation cost.
Helicopter-V/STOL dynamic wind and turbulence design methodology
NASA Technical Reports Server (NTRS)
Bailey, J. Earl
1987-01-01
Aircraft and helicopter accidents due to severe dynamic wind and turbulence continue to present challenging design problems. The development of the current set of design analysis tools for a aircraft wind and turbulence design began in the 1940's and 1950's. The areas of helicopter dynamic wind and turbulence modeling and vehicle response to severe dynamic wind inputs (microburst type phenomena) during takeoff and landing remain as major unsolved design problems from a lack of both environmental data and computational methodology. The development of helicopter and V/STOL dynamic wind and turbulence response computation methology is reviewed, the current state of the design art in industry is outlined, and comments on design methodology are made which may serve to improve future flight vehicle design.
Twisdale, L.A.; Dunn, W.L.
1981-08-01
A probabilistic methodology has been developed to predict the probabilities of tornado-propelled missiles impacting and damaging nuclear power plant structures. Mathematical models of each event in the tornado missile hazard have been developed and sequenced to form an integrated, time-history simulation methodology. The models are data based where feasible. The data include documented records of tornado occurrence, field observations of missile transport, results of wind tunnel experiments, and missile impact tests. Probabilistic Monte Carlo techniques are used to estimate the risk probabilities. The methodology has been encoded in the TORMIS computer code to facilitate numerical analysis and plant-specific tornado missile probability assessments. Sensitivity analyses have been performed on both the individual models and the integrated methodology, and risk has been assessed for a hypothetical nuclear power plant design case study.
Fernández-Arévalo, T; Lizarralde, I; Grau, P; Ayesa, E
2014-09-01
This paper presents a new modelling methodology for dynamically predicting the heat produced or consumed in the transformations of any biological reactor using Hess's law. Starting from a complete description of model components stoichiometry and formation enthalpies, the proposed modelling methodology has integrated successfully the simultaneous calculation of both the conventional mass balances and the enthalpy change of reaction in an expandable multi-phase matrix structure, which facilitates a detailed prediction of the main heat fluxes in the biochemical reactors. The methodology has been implemented in a plant-wide modelling methodology in order to facilitate the dynamic description of mass and heat throughout the plant. After validation with literature data, as illustrative examples of the capability of the methodology, two case studies have been described. In the first one, a predenitrification-nitrification dynamic process has been analysed, with the aim of demonstrating the easy integration of the methodology in any system. In the second case study, the simulation of a thermal model for an ATAD has shown the potential of the proposed methodology for analysing the effect of ventilation and influent characterization. PMID:24852412
Simulation visualization through dynamic instrumentation
Bisset, K.R.
1998-09-01
The goal of the instrument composition system is to allow a simulation user to dynamically create instruments as a simulation executes. Instruments can include graphical displays, data collectors, and debugging aides. Instruments are made up of small building blocks which can be easily combined into larger, more complex instruments. Through the sue of an Attribute Server (a distributed publication/subscription mechanism), the actors and instruments in a simulation can interact without direct knowledge of each other. Instead, each actor publishes the attributes which it has available. An instrument subscribes to the attributes in which it is interested, and is notified whenever the value of one of these attribute changes. An instrument can also publish attributes for use by other instruments. Since the Attribute Server is distributed, the publisher of an attribute need not execute on the same machine as the subscriber. This allows CPU intensive data visualization to execute on separate machines from the simulation, minimizing the impact on the simulation.
VERA Core Simulator Methodology for PWR Cycle Depletion
Kochunas, Brendan; Collins, Benjamin S; Jabaay, Daniel; Kim, Kang Seog; Graham, Aaron; Stimpson, Shane; Wieselquist, William A; Clarno, Kevin T; Palmtag, Scott; Downar, Thomas; Gehin, Jess C
2015-01-01
This paper describes the methodology developed and implemented in MPACT for performing high-fidelity pressurized water reactor (PWR) multi-cycle core physics calculations. MPACT is being developed primarily for application within the Consortium for the Advanced Simulation of Light Water Reactors (CASL) as one of the main components of the VERA Core Simulator, the others being COBRA-TF and ORIGEN. The methods summarized in this paper include a methodology for performing resonance self-shielding and computing macroscopic cross sections, 2-D/1-D transport, nuclide depletion, thermal-hydraulic feedback, and other supporting methods. These methods represent a minimal set needed to simulate high-fidelity models of a realistic nuclear reactor. Results demonstrating this are presented from the simulation of a realistic model of the first cycle of Watts Bar Unit 1. The simulation, which approximates the cycle operation, is observed to be within 50 ppm boron (ppmB) reactivity for all simulated points in the cycle and approximately 15 ppmB for a consistent statepoint. The verification and validation of the PWR cycle depletion capability in MPACT is the focus of two companion papers.
An automated methodology development. [software design for combat simulation
NASA Technical Reports Server (NTRS)
Hawley, L. R.
1985-01-01
The design methodology employed in testing the applicability of Ada in large-scale combat simulations is described. Ada was considered as a substitute for FORTRAN to lower life cycle costs and ease the program development efforts. An object-oriented approach was taken, which featured definitions of military targets, the capability of manipulating their condition in real-time, and one-to-one correlation between the object states and real world states. The simulation design process was automated by the problem statement language (PSL)/problem statement analyzer (PSA). The PSL/PSA system accessed the problem data base directly to enhance the code efficiency by, e.g., eliminating non-used subroutines, and provided for automated report generation, besides allowing for functional and interface descriptions. The ways in which the methodology satisfied the responsiveness, reliability, transportability, modifiability, timeliness and efficiency goals are discussed.
Aggregate Building Simulator (ABS) Methodology Development, Application, and User Manual
Dirks, James A.; Gorrissen, Willy J.
2011-11-30
As the relationship between the national building stock and various global energy issues becomes a greater concern, it has been deemed necessary to develop a system of predicting the energy consumption of large groups of buildings. Ideally this system is to take advantage of the most advanced energy simulation software available, be able to execute runs quickly, and provide concise and useful results at a level of detail that meets the users needs without inundating them with data. The resulting methodology that was developed allows the user to quickly develop and execute energy simulations of many buildings simultaneously, taking advantage of parallel processing to greatly reduce total simulation times. The result of these simulations can then be rapidly condensed and presented in a useful and intuitive manner.
Simulation supported POD: Methodology and HFET validation case
Jenson, F.; Iakovleva, E.; Dominguez, N.
2011-06-23
Structure reliability guaranty requires prior evaluation of non destructive testing methods. The concept of Probability of Detection (POD) is generally used to quantitatively assess performances and reliability of testing operations. Such probabilistic approaches take into account the uncertainties that appear during inspections and that are responsible for the output variability. POD curve determination is based on costly and time consuming experimental campaigns. Increasing demand of NDT configurations requiring POD evaluation makes cost reduction of POD campaigns a major issue. A new trend is to apply simulation in the context of probabilistic approaches in order to replace some of the experimental data required to determine the POD with simulation results. This paper presents results of simulation based POD curves of a high frequency eddy current inspection procedure obtained with the new POD module of CIVA. The methodology used for describing uncertainties on the input simulation parameters is described and comparisons with experimental results are presented and discussed.
Dynamic simulations of membranes with cytoskeletal interactions
NASA Astrophysics Data System (ADS)
Lin, Lawrence C.-L.; Brown, Frank L. H.
2005-07-01
We describe a simulation algorithm for the dynamics of elastic membrane sheets over long length and time scales. Our model includes implicit hydrodynamic coupling between membrane and surrounding solvent and allows for arbitrary external forces acting on the membrane surface. In particular, the methodology is well suited to studying membranes in interaction with cytoskeletal filaments. We present results for the thermal undulations of a lipid bilayer attached to a regular network of spectrin filaments as a model for the red blood cell membrane. The dynamic fluctuations of the bilayer over the spectrin network are quantified and used to predict the macroscopic diffusion constant of band 3 on the surface of the red blood cell. We find that thermal undulations likely play a role in the mobility of band 3 in the plane of the erythrocyte membrane.
Dynamic Simulation Nuclear Power Plants
Energy Science and Technology Software Center (ESTSC)
1992-03-03
DSNP (Dynamic Simulator for Nuclear Power-Plants) is a system of programs and data files by which a nuclear power plant, or part thereof, can be simulated. The acronym DSNP is used interchangeably for the DSNP language, the DSNP libraries, the DSNP precompiler, and the DSNP document generator. The DSNP language is a special-purpose, block-oriented, digital-simulation language developed to facilitate the preparation of dynamic simulations of a large variety of nuclear power plants. It is amore » user-oriented language that permits the user to prepare simulation programs directly from power plant block diagrams and flow charts by recognizing the symbolic DSNP statements for the appropriate physical components and listing these statements in a logical sequence according to the flow of physical properties in the simulated power plant. Physical components of nuclear power plants are represented by functional blocks, or modules. Many of the more complex components are represented by several modules. The nuclear reactor, for example, has a kinetic module, a power distribution module, a feedback module, a thermodynamic module, a hydraulic module, and a radioactive heat decay module. These modules are stored in DSNP libraries in the form of a DSNP subroutine or function, a block of statements, a macro, or a combination of the above. Basic functional blocks such as integrators, pipes, function generators, connectors, and many auxiliary functions representing properties of materials used in nuclear power plants are also available. The DSNP precompiler analyzes the DSNP simulation program, performs the appropriate translations, inserts the requested modules from the library, links these modules together, searches necessary data files, and produces a simulation program in FORTRAN.« less
Accelerated dynamics simulations of nanotubes.
Uberuaga, B. P.; Stuart, S. J.; Voter, A. F.
2002-01-01
We report on the application of accelerated dynamics techniques to the study of carbon nanotubes. We have used the parallel replica method and temperature accelerated dynamics simulations are currently in progress. In the parallel replica study, we have stretched tubes at a rate significantly lower than that used in previous studies. In these preliminary results, we find that there are qualitative differences in the rupture of the nanotubes at different temperatures. We plan on extending this investigation to include nanotubes of various chiralities. We also plan on exploring unique geometries of nanotubes.
Radiation in molecular dynamic simulations
Glosli, J; Graziani, F; More, R; Murillo, M; Streitz, F; Surh, M
2008-10-13
Hot dense radiative (HDR) plasmas common to Inertial Confinement Fusion (ICF) and stellar interiors have high temperature (a few hundred eV to tens of keV), high density (tens to hundreds of g/cc) and high pressure (hundreds of Megabars to thousands of Gigabars). Typically, such plasmas undergo collisional, radiative, atomic and possibly thermonuclear processes. In order to describe HDR plasmas, computational physicists in ICF and astrophysics use atomic-scale microphysical models implemented in various simulation codes. Experimental validation of the models used to describe HDR plasmas are difficult to perform. Direct Numerical Simulation (DNS) of the many-body interactions of plasmas is a promising approach to model validation but, previous work either relies on the collisionless approximation or ignores radiation. We present a new numerical simulation technique to address a currently unsolved problem: the extension of molecular dynamics to collisional plasmas including emission and absorption of radiation. The new technique passes a key test: it relaxes to a blackbody spectrum for a plasma in local thermodynamic equilibrium. This new tool also provides a method for assessing the accuracy of energy and momentum exchange models in hot dense plasmas. As an example, we simulate the evolution of non-equilibrium electron, ion, and radiation temperatures for a hydrogen plasma using the new molecular dynamics simulation capability.
NASA Astrophysics Data System (ADS)
Li, Jun; Fu, Siyao; He, Haibo; Jia, Hongfei; Li, Yanzhong; Guo, Yi
2015-11-01
Large-scale regional evacuation is an important part of national security emergency response plan. Large commercial shopping area, as the typical service system, its emergency evacuation is one of the hot research topics. A systematic methodology based on Cellular Automata with the Dynamic Floor Field and event driven model has been proposed, and the methodology has been examined within context of a case study involving the evacuation within a commercial shopping mall. Pedestrians walking is based on Cellular Automata and event driven model. In this paper, the event driven model is adopted to simulate the pedestrian movement patterns, the simulation process is divided into normal situation and emergency evacuation. The model is composed of four layers: environment layer, customer layer, clerk layer and trajectory layer. For the simulation of movement route of pedestrians, the model takes into account purchase intention of customers and density of pedestrians. Based on evacuation model of Cellular Automata with Dynamic Floor Field and event driven model, we can reflect behavior characteristics of customers and clerks at the situations of normal and emergency evacuation. The distribution of individual evacuation time as a function of initial positions and the dynamics of the evacuation process is studied. Our results indicate that the evacuation model using the combination of Cellular Automata with Dynamic Floor Field and event driven scheduling can be used to simulate the evacuation of pedestrian flows in indoor areas with complicated surroundings and to investigate the layout of shopping mall.
Molecular dynamics simulation of benzene
NASA Astrophysics Data System (ADS)
Trumpakaj, Zygmunt; Linde, Bogumił B. J.
2016-03-01
Intermolecular potentials and a few models of intermolecular interaction in liquid benzene are tested by Molecular Dynamics (MD) simulations. The repulsive part of the Lennard-Jones 12-6 (LJ 12-6) potential is too hard, which yields incorrect results. The exp-6 potential with a too hard repulsive term is also often used. Therefore, we took an expa-6 potential with a small Gaussian correction plus electrostatic interactions. This allows to modify the curvature of the potential. The MD simulations are carried out in the temperature range 280-352 K under normal pressure and at experimental density. The Rayleigh scattering of depolarized light is used for comparison. The results of MD simulations are comparable with the experimental values.
Dynamic simulations of tissue welding
Maitland, D.J.; Eder, D.C.; London, R.A.; Glinsky, M.E.
1996-02-01
The exposure of human skin to near-infrared radiation is numerically simulated using coupled laser, thermal transport and mass transport numerical models. The computer model LATIS is applied in both one-dimensional and two-dimensional geometries. Zones within the skin model are comprised of a topical solder, epidermis, dermis, and fatty tissue. Each skin zone is assigned initial optical, thermal and water density properties consistent with values listed in the literature. The optical properties of each zone (i.e. scattering, absorption and anisotropy coefficients) are modeled as a kinetic function of the temperature. Finally, the water content in each zone is computed from water diffusion where water losses are accounted for by evaporative losses at the air-solder interface. The simulation results show that the inclusion of water transport and evaporative losses in the model are necessary to match experimental observations. Dynamic temperature and damage distributions are presented for the skin simulations.
NASA Astrophysics Data System (ADS)
Rose, D. V.; Welch, D. R.; Genoni, T. C.; Mehlhorn, T. A.; Campbell, R. B.
2008-03-01
Particle-based numerical simulations are required to study the dynamics and evolution of inhomogeneous nonequilibrium multispecies strongly coupled plasmas. Molecular dynamics (MD) and particle-in-cell (PIC) techniques and been compared previously [K. Y. Sanbonmatsu, et al., J. Phys. IV (France) 10, Pr5-259 (2000)], with the PIC methodology demonstrating the capability of improved accuracy over the MD simulations at high resolution. However, the PIC simulations were significantly slower, limiting their utility. Here we explore several schemes to improve the computational speed of such calculations including non-iterative, implicit EM field solvers and subgrid models. The simulations are compared directly with the results of Sanbonmatsu, et al., and a new theoretical analysis of the hypernetted chain model where all inter-species correlations are retained [V. Schwarz, et al., Contrib. Plasma Phys. 47, 324 (2007)].
Gamma ray observatory dynamics simulator in Ada (GRODY)
NASA Technical Reports Server (NTRS)
1990-01-01
This experiment involved the parallel development of dynamics simulators for the Gamma Ray Observatory in both FORTRAN and Ada for the purpose of evaluating the applicability of Ada to the NASA/Goddard Space Flight Center's flight dynamics environment. The experiment successfully demonstrated that Ada is a viable, valuable technology for use in this environment. In addition to building a simulator, the Ada team evaluated training approaches, developed an Ada methodology appropriate to the flight dynamics environment, and established a baseline for evaluating future Ada projects.
The SIMRAND methodology - Simulation of Research and Development Projects
NASA Technical Reports Server (NTRS)
Miles, R. F., Jr.
1984-01-01
In research and development projects, a commonly occurring management decision is concerned with the optimum allocation of resources to achieve the project goals. Because of resource constraints, management has to make a decision regarding the set of proposed systems or tasks which should be undertaken. SIMRAND (Simulation of Research and Development Projects) is a methodology which was developed for aiding management in this decision. Attention is given to a problem description, aspects of model formulation, the reduction phase of the model solution, the simulation phase, and the evaluation phase. The implementation of the considered approach is illustrated with the aid of an example which involves a simplified network of the type used to determine the price of silicon solar cells.
Experimental methodology for computational fluid dynamics code validation
Aeschliman, D.P.; Oberkampf, W.L.
1997-09-01
Validation of Computational Fluid Dynamics (CFD) codes is an essential element of the code development process. Typically, CFD code validation is accomplished through comparison of computed results to previously published experimental data that were obtained for some other purpose, unrelated to code validation. As a result, it is a near certainty that not all of the information required by the code, particularly the boundary conditions, will be available. The common approach is therefore unsatisfactory, and a different method is required. This paper describes a methodology developed specifically for experimental validation of CFD codes. The methodology requires teamwork and cooperation between code developers and experimentalists throughout the validation process, and takes advantage of certain synergisms between CFD and experiment. The methodology employs a novel uncertainty analysis technique which helps to define the experimental plan for code validation wind tunnel experiments, and to distinguish between and quantify various types of experimental error. The methodology is demonstrated with an example of surface pressure measurements over a model of varying geometrical complexity in laminar, hypersonic, near perfect gas, 3-dimensional flow.
Development of test methodology for dynamic mechanical analysis instrumentation
NASA Technical Reports Server (NTRS)
Allen, V. R.
1982-01-01
Dynamic mechanical analysis instrumentation was used for the development of specific test methodology in the determination of engineering parameters of selected materials, esp. plastics and elastomers, over a broad range of temperature with selected environment. The methodology for routine procedures was established with specific attention given to sample geometry, sample size, and mounting techniques. The basic software of the duPont 1090 thermal analyzer was used for data reduction which simplify the theoretical interpretation. Clamps were developed which allowed 'relative' damping during the cure cycle to be measured for the fiber-glass supported resin. The correlation of fracture energy 'toughness' (or impact strength) with the low temperature (glassy) relaxation responses for a 'rubber-modified' epoxy system was negative in result because the low-temperature dispersion mode (-80 C) of the modifier coincided with that of the epoxy matrix, making quantitative comparison unrealistic.
Parametric design methodology for chemical processes using a simulator
Diwekar, U.M.; Rubin, E.S. )
1994-02-01
Parameter design is a method popularized by the Japanese quality expert G. Taguchi, for designing products and manufacturing processes that are robust in the face of uncontrollable variations. At the design stage, the goal of parameter design is to identify design settings that make the product performance less sensitive to the effects of manufacturing and environmental variations and deterioration. Because parameter design reduces performance variation by reducing the influence of the sources of variation rather than by controlling them, it is a cost-effective technique for improving quality. A recent study on the application of parameter design methodology for chemical processes reported that the use of Taguchi's method was not justified and a method based on Monte Carlo simulation combined with optimization was shown to be more effective. However, this method is computationally intensive as a large number of samples are necessary to achieve the given accuracy. Additionally, determination of the number of sample runs required is based on experimentation due to a lack of systematic sampling methods. In an attempt to overcome these problems, the use of a stochastic modeling capability combined with an optimizer is presented in this paper. The objective is that of providing an effective means for application of parameter design methodologies to chemical processes using the ASPEN simulator. This implementation not only presents a generalized tool for use by chemical engineers at large but also provides systematic estimates of the number of sample runs required to attain the specified accuracy. The stochastic model employs the technique of Latin hypercube sampling instead of the traditional Monte Carlo technique and hence has a great potential to reduce the required number of samples. The methodology is illustrated via an example problem of designing a chemical process.
Advanced Methodology for Simulation of Complex Flows Using Structured Grid Systems
NASA Technical Reports Server (NTRS)
Steinthorsson, Erlendur; Modiano, David
1995-01-01
Detailed simulations of viscous flows in complicated geometries pose a significant challenge to current capabilities of Computational Fluid Dynamics (CFD). To enable routine application of CFD to this class of problems, advanced methodologies are required that employ (a) automated grid generation, (b) adaptivity, (c) accurate discretizations and efficient solvers, and (d) advanced software techniques. Each of these ingredients contributes to increased accuracy, efficiency (in terms of human effort and computer time), and/or reliability of CFD software. In the long run, methodologies employing structured grid systems will remain a viable choice for routine simulation of flows in complex geometries only if genuinely automatic grid generation techniques for structured grids can be developed and if adaptivity is employed more routinely. More research in both these areas is urgently needed.
Molecular Dynamics Simulations of Polymers
NASA Astrophysics Data System (ADS)
Han, Jie
1995-01-01
Molecular dynamics (MD) simulations have been undertaken in this work to explore structures and properties of polyethylene (PE), polyisobutylene (PIB), atactic polypropylene (aPP) and atactic polystyrene (aPS). This work has not only demonstrated the reliability of MD simulations by comparing results with available experiments, but more importantly has revealed structure-property relationships on a molecular level for these selected polymers. Structures of these amorphous polymers were characterized by radial distribution functions (RDFs) or scattering profiles, and properties of the polymers studied were pressure-volume -temperature (PVT) equation of state, enthalpy, cohesive energy, the diffusion coefficient of methane in the polymer, and glass transition temperature. Good agreement was found for these structures and properties between simulation and experiment. More importantly, the scientific understanding of structure-property relationships was established on a molecular level. In the order of aPP (PE), PIB and aPS, with the chain surface separation or free volume decreasing, the density increases and the diffusion coefficient decreases. Therefore, the effects of changes or modifications in the chemical structure of monomer molecules (substituting pendent hydrogen with methyl or phenyl) on polymeric materials performance were attributed to the effects of molecular chain structure on packing structure, which, in turn, affects the properties of these polymers. Local chain dynamics and relaxation have been studied for bulk PE and aPS. Cooperative transitions occur at second-neighbor bonds for PE, and first-neighbor bonds for aPS due to the role of side groups. The activation energy is a single torsional barrier for overall conformational transitions, and is single torsional barrier plus locally "trapped" barrier for relaxation. Temperature dependence is Arrhenius for transition time, and is WLF for relaxation time. The mean correlation times derived from
Dynamic risk simulation to assess risk along roads
NASA Astrophysics Data System (ADS)
Voumard, Jérémie; Caspar, Olivier; Derron, Marc-Henri; Jaboyedoff, Michel
2013-04-01
Risk generated through natural hazards on roads is usually calculated with an equation which integrates various parameters of hazard and traffic. These are static variables as hazard frequency and number of vehicles crossing the dangerous section. This traditional methodology cannot take into account the dynamic variations of traffic and interactions between vehicles such as speeds modifications due to the section sinuosity, slowdowns resulting saturated traffic or creation of vehicles columns in front of traffic lights. The influence of traffic dynamics on the risk estimation is not addressed with the standard methodologies. Here we show, with the help of a dynamic traffic simulator specially developed for this project, that the variations of traffic greatly influence the risk results. Several sections of an alpine road in Switzerland were analyzed with the method of dynamic risk and compared with the results of the conventional method of risk calculation. It was possible to demonstrate that risk significantly increases on sinuous sections with the decreasing of vehicles speed. It has been highlighted that well positioned traffic lights, outside the risk area, can prevent a risk increase during a construction site, while a column of vehicles located within the danger zone greatly increases the risk. These results demonstrate the importance to consider the traffic in a dynamic way to assess risk to the closest field reality. Thus, recommendations to reduce risk on the roads can be given using a dynamic traffic simulator, modeling interactions between vehicles. Eventually, dynamic risk assessment can also be integrated into existing methodologies that consider only static parameters.
Methodological aspects of EEG and body dynamics measurements during motion
Reis, Pedro M. R.; Hebenstreit, Felix; Gabsteiger, Florian; von Tscharner, Vinzenz; Lochmann, Matthias
2014-01-01
EEG involves the recording, analysis, and interpretation of voltages recorded on the human scalp which originate from brain gray matter. EEG is one of the most popular methods of studying and understanding the processes that underlie behavior. This is so, because EEG is relatively cheap, easy to wear, light weight and has high temporal resolution. In terms of behavior, this encompasses actions, such as movements that are performed in response to the environment. However, there are methodological difficulties which can occur when recording EEG during movement such as movement artifacts. Thus, most studies about the human brain have examined activations during static conditions. This article attempts to compile and describe relevant methodological solutions that emerged in order to measure body and brain dynamics during motion. These descriptions cover suggestions on how to avoid and reduce motion artifacts, hardware, software and techniques for synchronously recording EEG, EMG, kinematics, kinetics, and eye movements during motion. Additionally, we present various recording systems, EEG electrodes, caps and methods for determinating real/custom electrode positions. In the end we will conclude that it is possible to record and analyze synchronized brain and body dynamics related to movement or exercise tasks. PMID:24715858
Methodological aspects of EEG and body dynamics measurements during motion.
Reis, Pedro M R; Hebenstreit, Felix; Gabsteiger, Florian; von Tscharner, Vinzenz; Lochmann, Matthias
2014-01-01
EEG involves the recording, analysis, and interpretation of voltages recorded on the human scalp which originate from brain gray matter. EEG is one of the most popular methods of studying and understanding the processes that underlie behavior. This is so, because EEG is relatively cheap, easy to wear, light weight and has high temporal resolution. In terms of behavior, this encompasses actions, such as movements that are performed in response to the environment. However, there are methodological difficulties which can occur when recording EEG during movement such as movement artifacts. Thus, most studies about the human brain have examined activations during static conditions. This article attempts to compile and describe relevant methodological solutions that emerged in order to measure body and brain dynamics during motion. These descriptions cover suggestions on how to avoid and reduce motion artifacts, hardware, software and techniques for synchronously recording EEG, EMG, kinematics, kinetics, and eye movements during motion. Additionally, we present various recording systems, EEG electrodes, caps and methods for determinating real/custom electrode positions. In the end we will conclude that it is possible to record and analyze synchronized brain and body dynamics related to movement or exercise tasks. PMID:24715858
A Posteriori Analysis for Hydrodynamic Simulations Using Adjoint Methodologies
Woodward, C S; Estep, D; Sandelin, J; Wang, H
2009-02-26
This report contains results of analysis done during an FY08 feasibility study investigating the use of adjoint methodologies for a posteriori error estimation for hydrodynamics simulations. We developed an approach to adjoint analysis for these systems through use of modified equations and viscosity solutions. Targeting first the 1D Burgers equation, we include a verification of the adjoint operator for the modified equation for the Lax-Friedrichs scheme, then derivations of an a posteriori error analysis for a finite difference scheme and a discontinuous Galerkin scheme applied to this problem. We include some numerical results showing the use of the error estimate. Lastly, we develop a computable a posteriori error estimate for the MAC scheme applied to stationary Navier-Stokes.
Simulation of Fault Zone Dynamics
NASA Astrophysics Data System (ADS)
Mora, P.; Abe, S.; Place, D.
2002-12-01
Particle models such as the discrete element model for granular assemblies and the lattice solid model provide a means to study the dynamics of fault zones. The lattice solid model was developed with the aim of progressively building up the capacity to simulate all relevent physical processes in fault zones. The present implementation of the model is able to simulate the dynamics of a granular lattice consisting of bonded or unbonded circular (2D) or spherical (3D) particles. Thermal effects (frictional hear generation, thermal expansion, heat flow) and pore fluid effects (heat induced pore pressure gradients and the consequent Darcian flow and impact on effective friction) can be modelled. Past work involving both circular particles and non-circular grains constructed as groups of bonded particles have demonstrated that grain shape has a fundamental impact on zero-th order behaviour. When circular particles are used, rolling is the most efficient means to accomodate slip of a simulated fault gouge layer leading to unrealistically low friction, typically around 0.2. This is consistent with laboratory results by Mair and Marone which have demonstrated that gouge consisting entirely of spherical beads shows a lower coefficient of friction than gouge containing irregular shaped particles. Recent work comparing quasi-2D laboratory results using pasta (Marone) with 2D numerical results (Morgan) have confirmed that numerical and laboratory results with circular ``particles'' are in agreement. When irregular grains are modelled at the lowest scale, the friction of simulated gouge layers matches with laboratory observations of rock friction (μ ~ 0.6) and is insentitive to the value used for interparticle friction (Mora et al, 2000). This indicates a self-regulation mechanism is occurring in which the group behaviour of the gouge layer remains constant at around 0.6 by balancing the amount of slip and rolling of grains within the gouge layer. A limitation of these studies
A novel bias correction methodology for climate impact simulations
NASA Astrophysics Data System (ADS)
Sippel, S.; Otto, F. E. L.; Forkel, M.; Allen, M. R.; Guillod, B. P.; Heimann, M.; Reichstein, M.; Seneviratne, S. I.; Thonicke, K.; Mahecha, M. D.
2016-02-01
Understanding, quantifying and attributing the impacts of extreme weather and climate events in the terrestrial biosphere is crucial for societal adaptation in a changing climate. However, climate model simulations generated for this purpose typically exhibit biases in their output that hinder any straightforward assessment of impacts. To overcome this issue, various bias correction strategies are routinely used to alleviate climate model deficiencies, most of which have been criticized for physical inconsistency and the nonpreservation of the multivariate correlation structure. In this study, we introduce a novel, resampling-based bias correction scheme that fully preserves the physical consistency and multivariate correlation structure of the model output. This procedure strongly improves the representation of climatic extremes and variability in a large regional climate model ensemble (HadRM3P, climateprediction.net/weatherathome), which is illustrated for summer extremes in temperature and rainfall over Central Europe. Moreover, we simulate biosphere-atmosphere fluxes of carbon and water using a terrestrial ecosystem model (LPJmL) driven by the bias-corrected climate forcing. The resampling-based bias correction yields strongly improved statistical distributions of carbon and water fluxes, including the extremes. Our results thus highlight the importance of carefully considering statistical moments beyond the mean for climate impact simulations. In conclusion, the present study introduces an approach to alleviate climate model biases in a physically consistent way and demonstrates that this yields strongly improved simulations of climate extremes and associated impacts in the terrestrial biosphere. A wider uptake of our methodology by the climate and impact modelling community therefore seems desirable for accurately quantifying changes in past, current and future extremes.
Dynamic risk simulation to assess risk along roads
NASA Astrophysics Data System (ADS)
Voumard, J.; Caspar, O.; Derron, M.-H.; Jaboyedoff, M.
2013-04-01
Risk generated by natural hazards on roads is usually calculated with equations integrating various parameters related to hazard and traffic. These are static variables, like a rockfall hazard estimation for a road section or the average number of vehicles crossing this section every day. This methodology cannot take into account dynamic variations of traffic and interactions between vehicles such as speed modifications due to the sinuosity, slowdowns resulting of saturated traffic or vehicles columns forming in front of traffic lights. The influence of traffic dynamics on the risk estimation is not assessed with standard methodologies. Here we show, by mean of a dynamic traffic simulator, that the traffic variations may greatly influence the risk estimation over time. The risk is analysed on several alpine road sections in Switzerland using a dynamic vehicles approach and compared with the results of the static methodology. It demonstrates that risk significantly increases on sinuous sections because of the decreasing of vehicles speed. A more realistic risk can be obtained from a dynamic approach especially on mountain roads. A dynamic traffic simulator, modelling interactions between vehicles is a helpful tool to support decision making to reduce risk on roads.
Modeling Multibody Stage Separation Dynamics Using Constraint Force Equation Methodology
NASA Technical Reports Server (NTRS)
Tartabini, Paul V.; Roithmayr, Carlos M.; Toniolo, Matthew D.; Karlgaard, Christopher D.; Pamadi, Bandu N.
2011-01-01
This paper discusses the application of the constraint force equation methodology and its implementation for multibody separation problems using three specially designed test cases. The first test case involves two rigid bodies connected by a fixed joint, the second case involves two rigid bodies connected with a universal joint, and the third test case is that of Mach 7 separation of the X-43A vehicle. For the first two cases, the solutions obtained using the constraint force equation method compare well with those obtained using industry- standard benchmark codes. For the X-43A case, the constraint force equation solutions show reasonable agreement with the flight-test data. Use of the constraint force equation method facilitates the analysis of stage separation in end-to-end simulations of launch vehicle trajectories
A Methodology for Evaluating the Fidelity of Ground-Based Flight Simulators
NASA Technical Reports Server (NTRS)
Zeyada, Y.; Hess, R. A.
1999-01-01
An analytical and experimental investigation was undertaken to model the manner in which pilots perceive and utilize visual, proprioceptive, and vestibular cues in a ground-based flight simulator. The study was part of a larger research effort which has the creation of a methodology for determining flight simulator fidelity requirements as its ultimate goal. The study utilized a closed-loop feedback structure of the pilot/simulator system which included the pilot, the cockpit inceptor, the dynamics of the simulated vehicle and the motion system. With the exception of time delays which accrued in visual scene production in the simulator, visual scene effects were not included in this study. The NASA Ames Vertical Motion Simulator was used in a simple, single-degree of freedom rotorcraft bob-up/down maneuver. Pilot/vehicle analysis and fuzzy-inference identification were employed to study the changes in fidelity which occurred as the characteristics of the motion system were varied over five configurations i The data from three of the five pilots that participated in the experimental study were analyzed in the fuzzy inference identification. Results indicate that both the analytical pilot/vehicle analysis and the fuzzyinference identification can be used to reflect changes in simulator fidelity for the task examined.
A Methodology for Evaluating the Fidelity of Ground-Based Flight Simulators
NASA Technical Reports Server (NTRS)
Zeyada, Y.; Hess, R. A.
1999-01-01
An analytical and experimental investigation was undertaken to model the manner in which pilots perceive and utilize visual, proprioceptive, and vestibular cues in a ground-based flight simulator. The study was part of a larger research effort which has the creation of a methodology for determining flight simulator fidelity requirements as its ultimate goal. The study utilized a closed-loop feedback structure of the pilot/simulator system which included the pilot, the cockpit inceptor, the dynamics of the simulated vehicle and the motion system. With the exception of time delays which accrued in visual scene production in the simulator, visual scene effects were not included in this study. The NASA Ames Vertical Motion Simulator was used in a simple, single-degree of freedom rotorcraft bob-up/down maneuver. Pilot/vehicle analysis and fuzzy-inference identification were employed to study the changes in fidelity which occurred as the characteristics of the motion system were varied over five configurations. The data from three of the five pilots that participated in the experimental study were analyzed in the fuzzy-inference identification. Results indicate that both the analytical pilot/vehicle analysis and the fuzzy-inference identification can be used to reflect changes in simulator fidelity for the task examined.
Dynamics of riboswitches: Molecular simulations.
Sanbonmatsu, Karissa Y
2014-10-01
Riboswitch RNAs play key roles in bacterial metabolism and represent a promising new class of antibiotic targets for treatment of infectious disease. While many studies of riboswitches have been performed, the exact mechanism of riboswitch operation is still not fully understood at the atomistic level of detail. Molecular dynamics simulations are useful for interpreting existing experimental data and producing predictions for new experiments. Here, a wide range of computational studies on riboswitches is reviewed. By elucidating the key principles of riboswitch operation, computation may aid in the effort to design more specific antibiotics with affinities greater than those of the native ligand. Such a detailed understanding may be required to improve efficacy and reduce side effects. These studies are laying the groundwork for understanding the action mechanism of new compounds that inhibit riboswitch activity. Future directions such as magnesium effects, large-scale conformational changes, expression platforms and co-transcriptional folding are also discussed. This article is part of a Special Issue entitled: Riboswitches. PMID:24953187
Methodology for analysis and simulation of large multidisciplinary problems
NASA Technical Reports Server (NTRS)
Russell, William C.; Ikeda, Paul J.; Vos, Robert G.
1989-01-01
The Integrated Structural Modeling (ISM) program is being developed for the Air Force Weapons Laboratory and will be available for Air Force work. Its goal is to provide a design, analysis, and simulation tool intended primarily for directed energy weapons (DEW), kinetic energy weapons (KEW), and surveillance applications. The code is designed to run on DEC (VMS and UNIX), IRIS, Alliant, and Cray hosts. Several technical disciplines are included in ISM, namely structures, controls, optics, thermal, and dynamics. Four topics from the broad ISM goal are discussed. The first is project configuration management and includes two major areas: the software and database arrangement and the system model control. The second is interdisciplinary data transfer and refers to exchange of data between various disciplines such as structures and thermal. Third is a discussion of the integration of component models into one system model, i.e., multiple discipline model synthesis. Last is a presentation of work on a distributed processing computing environment.
A Measurement and Simulation Based Methodology for Cache Performance Modeling and Tuning
NASA Technical Reports Server (NTRS)
Waheed, Abdul; Yan, Jerry; Saini, Subhash (Technical Monitor)
1998-01-01
We present a cache performance modeling methodology that facilitates the tuning of uniprocessor cache performance for applications executing on shared memory multiprocessors by accurately predicting the effects of source code level modifications. Measurements on a single processor are initially used for identifying parts of code where cache utilization improvements may significantly impact the overall performance. Cache simulation based on trace-driven techniques can be carried out without gathering detailed address traces. Minimal runtime information for modeling cache performance of a selected code block includes: base virtual addresses of arrays, virtual addresses of variables, and loop bounds for that code block. Rest of the information is obtained from the source code. We show that the cache performance predictions are as reliable as those obtained through trace-driven simulations. This technique is particularly helpful to the exploration of various "what-if' scenarios regarding the cache performance impact for alternative code structures. We explain and validate this methodology using a simple matrix-matrix multiplication program. We then apply this methodology to predict and tune the cache performance of two realistic scientific applications taken from the Computational Fluid Dynamics (CFD) domain.
NASA Astrophysics Data System (ADS)
Reckinger, Scott J.; Livescu, Daniel; Vasilyev, Oleg V.
2016-05-01
An investigation of compressible Rayleigh-Taylor instability (RTI) using Direct Numerical Simulations (DNS) requires efficient numerical methods, advanced boundary conditions, and consistent initialization in order to capture the wide range of scales and vortex dynamics present in the system, while reducing the computational impact associated with acoustic wave generation and the subsequent interaction with the flow. An advanced computational framework is presented that handles the challenges introduced by considering the compressive nature of RTI systems, which include sharp interfacial density gradients on strongly stratified background states, acoustic wave generation and removal at computational boundaries, and stratification dependent vorticity production. The foundation of the numerical methodology described here is the wavelet-based grid adaptivity of the Parallel Adaptive Wavelet Collocation Method (PAWCM) that maintains symmetry in single-mode RTI systems to extreme late-times. PAWCM is combined with a consistent initialization, which reduces the generation of acoustic disturbances, and effective boundary treatments, which prevent acoustic reflections. A dynamic time integration scheme that can handle highly nonlinear and potentially stiff systems, such as compressible RTI, completes the computational framework. The numerical methodology is used to simulate two-dimensional single-mode RTI to extreme late-times for a wide range of flow compressibility and variable density effects. The results show that flow compressibility acts to reduce the growth of RTI for low Atwood numbers, as predicted from linear stability analysis.
Updating finite element dynamic models using an element-by-element sensitivity methodology
NASA Astrophysics Data System (ADS)
Farhat, Charbel; Hemez, Francois M.
1993-09-01
A sensitivity-based methodology for improving the finite element model of a given structure using test modal data and a few sensors is presented. The proposed method searches for both the location and sources of the mass and stiffness errors and does not interfere with the theory behind the finite element model while correcting these errors. The updating algorithm is derived from the unconstrained minimization of the squared L sub 2 norms of the modal dynamic residuals via an iterative two-step staggered procedure. At each iteration, the measured mode shapes are first expanded assuming that the model is error free, then the model parameters are corrected assuming that the expanded mode shapes are exact. The numerical algorithm is implemented in an element-by-element fashion and is capable of 'zooming' on the detected error locations. Several simulation examples which demonstate the potential of the proposed methodology are discussed.
Molecular dynamics simulations of large macromolecular complexes
Perilla, Juan R.; Goh, Boon Chong; Cassidy, C. Keith; Liu, Bo; Bernardi, Rafael C.; Rudack, Till; Yu, Hang; Wu, Zhe; Schulten, Klaus
2015-01-01
Connecting dynamics to structural data from diverse experimental sources, molecular dynamics simulations permit the exploration of biological phenomena in unparalleled detail. Advances in simulations are moving the atomic resolution descriptions of biological systems into the million-to-billion atom regime, in which numerous cell functions reside. In this opinion, we review the progress, driven by large-scale molecular dynamics simulations, in the study of viruses, ribosomes, bioenergetic systems, and other diverse applications. These examples highlight the utility of molecular dynamics simulations in the critical task of relating atomic detail to the function of supramolecular complexes, a task that cannot be achieved by smaller-scale simulations or existing experimental approaches alone. PMID:25845770
Methodology for computational fluid dynamics code verification/validation
Oberkampf, W.L.; Blottner, F.G.; Aeschliman, D.P.
1995-07-01
The issues of verification, calibration, and validation of computational fluid dynamics (CFD) codes has been receiving increasing levels of attention in the research literature and in engineering technology. Both CFD researchers and users of CFD codes are asking more critical and detailed questions concerning the accuracy, range of applicability, reliability and robustness of CFD codes and their predictions. This is a welcomed trend because it demonstrates that CFD is maturing from a research tool to the world of impacting engineering hardware and system design. In this environment, the broad issue of code quality assurance becomes paramount. However, the philosophy and methodology of building confidence in CFD code predictions has proven to be more difficult than many expected. A wide variety of physical modeling errors and discretization errors are discussed. Here, discretization errors refer to all errors caused by conversion of the original partial differential equations to algebraic equations, and their solution. Boundary conditions for both the partial differential equations and the discretized equations will be discussed. Contrasts are drawn between the assumptions and actual use of numerical method consistency and stability. Comments are also made concerning the existence and uniqueness of solutions for both the partial differential equations and the discrete equations. Various techniques are suggested for the detection and estimation of errors caused by physical modeling and discretization of the partial differential equations.
A dynamic hybrid RANS/LES modeling methodology for turbulent/transitional flow field prediction
NASA Astrophysics Data System (ADS)
Alam, Mohammad Faridul
A dynamic hybrid Reynolds-averaged Navier-Stokes (RANS)-Large Eddy Simulation (LES) modeling framework has been investigated and further developed to improve the Computational Fluid Dynamics (CFD) prediction of turbulent flow features along with laminar-to-turbulent transitional phenomena. In recent years, the use of hybrid RANS/LES (HRL) models has become more common in CFD simulations, since HRL models offer more accuracy than RANS in regions of flow separation at a reduced cost relative to LES in attached boundary layers. The first part of this research includes evaluation and validation of a dynamic HRL (DHRL) model that aims to address issues regarding the RANS-to-LES zonal transition and explicit grid dependence, both of which are inherent to most current HRL models. Simulations of two test cases---flow over a backward facing step and flow over a wing with leading-edge ice accretion---were performed to assess the potential of the DHRL model for predicting turbulent features involved in mainly unsteady separated flow. The DHRL simulation results are compared with experimental data, along with the computational results for other HRL and RANS models. In summary, these comparisons demonstrate that the DHRL framework does address many of the weaknesses inherent in most current HRL models. Although HRL models are widely used in turbulent flow simulations, they have limitations for transitional flow predictions. Most HRL models include a fully turbulent RANS component for attached boundary layer regions. The small number of HRL models that do include transition-sensitive RANS models have issues related to the RANS model itself and to the zonal transition between RANS and LES. In order to address those issues, a new transition-sensitive HRL modeling methodology has been developed that includes the DHRL methodology and a physics-based transition-sensitive RANS model. The feasibility of the transition-sensitive dynamic HRL (TDHRL) model has been investigated by
NASA Technical Reports Server (NTRS)
Prajous, R.; Mazankine, J.; Ippolito, J. C.
1978-01-01
Methods and algorithms used for the simulation of elementary power conditioning units buck, boost, and buck-boost, as well as shunt PWM are described. Definitions are given of similar converters and reduced parameters. The various parts of the simulation to be carried out are dealt with; local stability, corrective network, measurements of input-output impedance and global stability. A simulation example is given.
Enhanced Sampling Techniques in Molecular Dynamics Simulations of Biological Systems
Bernardi, Rafael C.; Melo, Marcelo C. R.; Schulten, Klaus
2014-01-01
Background Molecular Dynamics has emerged as an important research methodology covering systems to the level of millions of atoms. However, insufficient sampling often limits its application. The limitation is due to rough energy landscapes, with many local minima separated by high-energy barriers, which govern the biomolecular motion. Scope of review In the past few decades methods have been developed that address the sampling problem, such as replica-exchange molecular dynamics, metadynamics and simulated annealing. Here we present an overview over theses sampling methods in an attempt to shed light on which should be selected depending on the type of system property studied. Major Conclusions Enhanced sampling methods have been employed for a broad range of biological systems and the choice of a suitable method is connected to biological and physical characteristics of the system, in particular system size. While metadynamics and replica-exchange molecular dynamics are the most adopted sampling methods to study biomolecular dynamics, simulated annealing is well suited to characterize very flexible systems. The use of annealing methods for a long time was restricted to simulation of small proteins; however, a variant of the method, generalized simulated annealing, can be employed at a relatively low computational cost to large macromolecular complexes. General Significance Molecular dynamics trajectories frequently do not reach all relevant conformational substates, for example those connected with biological function, a problem that can be addressed by employing enhanced sampling algorithms. PMID:25450171
NASA Astrophysics Data System (ADS)
Riva, Fabio; Ricci, Paolo; Halpern, Federico D.; Jolliet, Sébastien; Loizu, Joaquim; Mosetto, Annamaria
2014-06-01
Bridging the gap between plasma physics and other scientific domains, in particular, the computational fluid dynamics community, a general, rigorous, and simple-to-apply methodology is presented for both the verification of the correct implementation of the model equations (code verification) and numerical error quantification (solution verification). The proposed code verification procedure consists in using the method of manufactured solutions and executing an order-of-accuracy test, assessing the rate of convergence of the numerical solution to the manufactured one. For the solution verification, the numerical error is quantified by applying the Richardson extrapolation, which provides an approximation of the analytical solution, and by using the grid convergence index to estimate the numerical uncertainty affecting the simulation results. The methodology is applied to verify the correct implementation of the drift-reduced Braginskii equations into the GBS code, and to estimate the numerical error affecting the GBS solutions. The GBS code is successfully verified, and an estimate of the numerical error affecting the simulation results is provided.
Calibration of a new very large eddy simulation (VLES) methodology for turbulent flow simulation
NASA Astrophysics Data System (ADS)
Han, XingSi; Ye, TaoHong; Chen, YiLiang
2012-10-01
Following the idea of Speziale's Very Large Eddy Simulation (VLES) method, a new unified hybrid simulation approach was proposed which can change seamlessly from RANS (Reynolds-Averaged Navier-Stokes) to LES (Large Eddy Simulation) method depending on the numerical resolution. The model constants were calibrated in accordance with other hybrid methods. Besides being able to approach the two limits of RANS and LES, the new model also provides a proper VLES mode between the two limits, and thus can be used for a wide range of mesh resolutions. Also RANS simulation can be recovered near the wall which is similar to the Detached Eddy Simulation (DES) concept. This new methodology was implemented into Wilcox's k-ω model and applications were conducted for fully developed turbulent channel flow at Re τ = 395 and turbulent flow past a square cylinder at Re = 22000. Results were compared with LES predictions and other studies. The new method is found to be quite efficient in resolving large flow structures, and can predict satisfactory results on relative coarse mesh.
Brownian Dynamics Simulations of Dispersed Graphene Sheets
NASA Astrophysics Data System (ADS)
Xu, Yueyi; Green, Micah
2013-03-01
Past simulations of the dynamics of dispersed graphene sheets are limited to static fluids on small timescales, with little attention devoted to flow dynamics. To address this need, we investigated how flow fields affect graphene morphology dynamics using a coarse-grained model; this relatively untouched area is critical given the importance of graphene solution-processing of multifunctional devices and materials. In particular, we developed a Brownian Dynamics (BD) algorithm to study the morphology of sheetlike macromolecules in dilute, flowing solutions. We used a bead-rod lattice to represent the mesoscopic conformation of individual two dimensional sheets. We then analyzed the morphology dynamic modes (stretching, tumbling, crumpling) of these molecules as a function of sheet size, Weissenberg number, and bending stiffness. Our results indicate the model can successfully simulate a range of dynamic modes in a given flow field and yield fundamental insight into the flow processing of graphene sheets.
The "Collisions Cube" Molecular Dynamics Simulator.
ERIC Educational Resources Information Center
Nash, John J.; Smith, Paul E.
1995-01-01
Describes a molecular dynamics simulator that employs ping-pong balls as the atoms or molecules and is suitable for either large lecture halls or small classrooms. Discusses its use in illustrating many of the fundamental concepts related to molecular motion and dynamics and providing a three-dimensional perspective of molecular motion. (JRH)
Simulating Flexible-Spacecraft Dynamics and Control
NASA Technical Reports Server (NTRS)
Fedor, Joseph
1987-01-01
Versatile program applies to many types of spacecraft and dynamical problems. Flexible Spacecraft Dynamics and Control program (FSD) developed to aid in simulation of large class of flexible and rigid spacecraft. Extremely versatile and used in attitude dynamics and control analysis as well as in-orbit support of deployment and control of spacecraft. Applicable to inertially oriented spinning, Earth-oriented, or gravity-gradient-stabilized spacecraft. Written in FORTRAN 77.
Data Mining: A Hybrid Methodology for Complex and Dynamic Research
ERIC Educational Resources Information Center
Lang, Susan; Baehr, Craig
2012-01-01
This article provides an overview of the ways in which data and text mining have potential as research methodologies in composition studies. It introduces data mining in the context of the field of composition studies and discusses ways in which this methodology can complement and extend our existing research practices by blending the best of what…
Educational Validity of Business Gaming Simulation: A Research Methodology Framework
ERIC Educational Resources Information Center
Stainton, Andrew J.; Johnson, Johnnie E.; Borodzicz, Edward P.
2010-01-01
Many past educational validity studies of business gaming simulation, and more specifically total enterprise simulation, have been inconclusive. Studies have focused on the weaknesses of business gaming simulation; which is often regarded as an educational medium that has limitations regarding learning effectiveness. However, no attempts have been…
Molecular Dynamics Simulations of Simple Liquids
ERIC Educational Resources Information Center
Speer, Owner F.; Wengerter, Brian C.; Taylor, Ramona S.
2004-01-01
An experiment, in which students were given the opportunity to perform molecular dynamics simulations on a series of molecular liquids using the Amber suite of programs, is presented. They were introduced to both physical theories underlying classical mechanics simulations and to the atom-atom pair distribution function.
[Methodologic aspects of body water kinetic dynamic studies].
Lobachik, V I; Chupushtanov, S A; Pishchulina, G N; Voronov, S F; Nosovskiĭ, A M
1998-01-01
In studying the level of hydration and liquid phases (LPs) in a space mission there use the stable and/or radioactive isotopes. The investigations are unique, the methods are adequate but not adapted in full measure to the challenges of the problem under study. The methodical approaches to the study of the dynamics of water metabolism are not available. Repeated introduction of the markers for these purposes is not acceptable. Another problem associates with taking the markers orally. In this case, a concentration of the markers will depend on the absorption and excretion processes. Prior to, during and in the readaptation period these functions will be different, making a correlation of the accumulated data difficult. There advances a possible version of solving these problems, namely, to use for the dynamic studies the residual contents of the markers in the LPs after single injection. However, this approach calls for investigating the kinetics of markers in weightlessness or during its ground-based simulation. The kinetics of tritium water has been studied in 6 volunteers under conditions of the 5-day bedrest and in 9 healthy men during free motor activity. There determined the characteristics of marker kinetics in a healthy man during his routine living activities. Under bedrest conditions there have been noted slowing-down of the rate and a decrease in the degree of marker accumulation in the body after its single injection, the shift of a period of relative stabilization of marker content in LP to the more late dates of experiment, slowing-down of the marker excretion rate from the body. PMID:9858979
Observing dynamical SUSY breaking with lattice simulation
Kanamori, Issaku
2008-11-23
On the basis of the recently developed lattice formulation of supersymmetric theories which keeps a part of the supersymmetry, we propose a method of observing dynamical SUSY breaking with lattice simulation. We use Hamiltonian as an order parameter and measure the ground state energy as a zero temperature limit of the finite temperature simulation. Our method provides a way of obtaining a physical result from the lattice simulation for supersymmetric theories.
Simulation of liquid dynamics onboard Sloshsat FLEVO
NASA Astrophysics Data System (ADS)
Vreeburg, J. P. B.
1999-01-01
The Sloshsat FLEVO project has an Investigators Working Group which prepared orbital experiments on the behavior of liquid in spacecraft. These are to be performed with a dedicated small satellite, of about 90 kg empty weight and about 34 kg of water in a 87 litre tank. The spacecraft dynamics are simulated by SMS, the Sloshsat Motion Simulator. SMS predictions and those generated by a CFD simulation are compared for an example.
Molecular dynamics simulations: advances and applications
Hospital, Adam; Goñi, Josep Ramon; Orozco, Modesto; Gelpí, Josep L
2015-01-01
Molecular dynamics simulations have evolved into a mature technique that can be used effectively to understand macromolecular structure-to-function relationships. Present simulation times are close to biologically relevant ones. Information gathered about the dynamic properties of macromolecules is rich enough to shift the usual paradigm of structural bioinformatics from studying single structures to analyze conformational ensembles. Here, we describe the foundations of molecular dynamics and the improvements made in the direction of getting such ensemble. Specific application of the technique to three main issues (allosteric regulation, docking, and structure refinement) is discussed.
Combining molecular dynamics with mesoscopic Green's function reaction dynamics simulations
NASA Astrophysics Data System (ADS)
Vijaykumar, Adithya; Bolhuis, Peter G.; ten Wolde, Pieter Rein
2015-12-01
In many reaction-diffusion processes, ranging from biochemical networks, catalysis, to complex self-assembly, the spatial distribution of the reactants and the stochastic character of their interactions are crucial for the macroscopic behavior. The recently developed mesoscopic Green's Function Reaction Dynamics (GFRD) method enables efficient simulation at the particle level provided the microscopic dynamics can be integrated out. Yet, many processes exhibit non-trivial microscopic dynamics that can qualitatively change the macroscopic behavior, calling for an atomistic, microscopic description. We propose a novel approach that combines GFRD for simulating the system at the mesoscopic scale where particles are far apart, with a microscopic technique such as Langevin dynamics or Molecular Dynamics (MD), for simulating the system at the microscopic scale where reactants are in close proximity. This scheme defines the regions where the particles are close together and simulated with high microscopic resolution and those where they are far apart and simulated with lower mesoscopic resolution, adaptively on the fly. The new multi-scale scheme, called MD-GFRD, is generic and can be used to efficiently simulate reaction-diffusion systems at the particle level.
Structure and dynamics of complex liquid water: Molecular dynamics simulation
NASA Astrophysics Data System (ADS)
S, Indrajith V.; Natesan, Baskaran
2015-06-01
We have carried out detailed structure and dynamical studies of complex liquid water using molecular dynamics simulations. Three different model potentials, namely, TIP3P, TIP4P and SPC-E have been used in the simulations, in order to arrive at the best possible potential function that could reproduce the structure of experimental bulk water. All the simulations were performed in the NVE micro canonical ensemble using LAMMPS. The radial distribution functions, gOO, gOH and gHH and the self diffusion coefficient, Ds, were calculated for all three models. We conclude from our results that the structure and dynamical parameters obtained for SPC-E model matched well with the experimental values, suggesting that among the models studied here, the SPC-E model gives the best structure and dynamics of bulk water.
DUVFEL PHOTOINJECTOR DYNAMICS: MEASUREMENT AND SIMULATION.
GRAVES, W.S.; DIMAURO, L.F.; HEESE, R.; JOHNSON, E.D.; ROSE, J.; RUDATI, J.; SHAFTAN, T.; SHEEHY, B.; YU, L.H.; DOWELL, D.H.
2001-06-18
The DUVFEL photoinjector consists of a 1.6 cell BNL gun IV with copper cathode, variable pulse length Ti:Sapp and solenoid magnet. The beam dynamics and the electromagnetic fields in the photoinjector have been characterized by producing a short electron beam with very low charge that is used as a field probe. Transverse beam size and divergence are measured as a function of initial RF phase and initial spot size and compared with simulations using the code HOMDYN. The electromagnetic fields used in the simulations are produced by SUPERFISH, and have been verified with RF measurements. The simulations and measurements of beam dynamics are presented.
Simulation of wetlands forest vegetation dynamics
Phipps, R.L.
1979-01-01
A computer program, SWAMP, was designed to simulate the effects of flood frequency and depth to water table on southern wetlands forest vegetation dynamics. By incorporating these hydrologic characteristics into the model, forest vegetation and vegetation dynamics can be simulated. The model, based on data from the White River National Wildlife Refuge near De Witt, Arkansas, "grows" individual trees on a 20 x 20-m plot taking into account effects on the tree growth of flooding, depth to water table, shade tolerance, overtopping and crowding, and probability of death and reproduction. A potential application of the model is illustrated with simulations of tree fruit production following flood-control implementation and lumbering. ?? 1979.
Dynamic system simulation of small satellite projects
NASA Astrophysics Data System (ADS)
Raif, Matthias; Walter, Ulrich; Bouwmeester, Jasper
2010-11-01
A prerequisite to accomplish a system simulation is to have a system model holding all necessary project information in a centralized repository that can be accessed and edited by all parties involved. At the Institute of Astronautics of the Technische Universitaet Muenchen a modular approach for modeling and dynamic simulation of satellite systems has been developed called dynamic system simulation (DySyS). DySyS is based on the platform independent description language SysML to model a small satellite project with respect to the system composition and dynamic behavior. A library of specific building blocks and possible relations between these blocks have been developed. From this library a system model of the satellite of interest can be created. A mapping into a C++ simulation allows the creation of an executable system model with which simulations are performed to observe the dynamic behavior of the satellite. In this paper DySyS is used to model and simulate the dynamic behavior of small satellites, because small satellite projects can act as a precursor to demonstrate the feasibility of a system model since they are less complex compared to a large scale satellite project.
Discrete dislocation dynamics simulations in a cylinder
NASA Astrophysics Data System (ADS)
Li, Maosheng; Gao, Chan; Xu, Jianing
2015-02-01
Mechanical properties of material are closely related to the motion of dislocations, and predicting the interactions and resulting collective motion of dislocations is a major task in understanding and modelling plastically deforming materials. A discrete dislocation dynamics model is used to describe the orientation substructure within the microstructure. Discrete dislocation dynamics simulations in three dimensions have been used to examine the role of dislocation multiplication and mobility on the plasticity in small samples under uniaxial compression. In this paper we describe the application of the dislocation dynamics simulations in a cylindrical geometry. The boundary conditions for the simulation were estimated from the distribution of the geometrically necessary dislocation density which was obtained from the orientation map. Numerical studies benchmark could validate the accuracy of the algorithms and the importance of handling the singularity correctly. The results of the simulation explain the formation of the experimentally observed substructure.
Spin dynamics simulations at AGS
Huang, H.; MacKay, W.W.; Meot, F.; Roser, T.
2010-05-23
To preserve proton polarization through acceleration, it is important to have a correct model of the process. It has been known that with the insertion of the two helical partial Siberian snakes in the Alternating Gradient Synchrotron (AGS), the MAD model of AGS can not deal with a field map with offset orbit. The stepwise ray-tracing code Zgoubi provides a tool to represent the real electromagnetic fields in the modeling of the optics and spin dynamics for the AGS. Numerical experiments of resonance crossing, including spin dynamics in presence of the snakes and Q-jump, have been performed in AGS lattice models, using Zgoubi. This contribution reports on various results so obtained.
Causality analysis in business performance measurement system using system dynamics methodology
NASA Astrophysics Data System (ADS)
Yusof, Zainuridah; Yusoff, Wan Fadzilah Wan; Maarof, Faridah
2014-07-01
One of the main components of the Balanced Scorecard (BSC) that differentiates it from any other performance measurement system (PMS) is the Strategy Map with its unidirectional causality feature. Despite its apparent popularity, criticisms on the causality have been rigorously discussed by earlier researchers. In seeking empirical evidence of causality, propositions based on the service profit chain theory were developed and tested using the econometrics analysis, Granger causality test on the 45 data points. However, the insufficiency of well-established causality models was found as only 40% of the causal linkages were supported by the data. Expert knowledge was suggested to be used in the situations of insufficiency of historical data. The Delphi method was selected and conducted in obtaining the consensus of the causality existence among the 15 selected expert persons by utilizing 3 rounds of questionnaires. Study revealed that only 20% of the propositions were not supported. The existences of bidirectional causality which demonstrate significant dynamic environmental complexity through interaction among measures were obtained from both methods. With that, a computer modeling and simulation using System Dynamics (SD) methodology was develop as an experimental platform to identify how policies impacting the business performance in such environments. The reproduction, sensitivity and extreme condition tests were conducted onto developed SD model to ensure their capability in mimic the reality, robustness and validity for causality analysis platform. This study applied a theoretical service management model within the BSC domain to a practical situation using SD methodology where very limited work has been done.
Peptide crystal simulations reveal hidden dynamics
Janowski, Pawel A.; Cerutti, David S.; Holton, James; Case, David A.
2013-01-01
Molecular dynamics simulations of biomolecular crystals at atomic resolution have the potential to recover information on dynamics and heterogeneity hidden in the X-ray diffraction data. We present here 9.6 microseconds of dynamics in a small helical peptide crystal with 36 independent copies of the unit cell. The average simulation structure agrees with experiment to within 0.28 Å backbone and 0.42 Å all-atom rmsd; a model refined against the average simulation density agrees with the experimental structure to within 0.20 Å backbone and 0.33 Å all-atom rmsd. The R-factor between the experimental structure factors and those derived from this unrestrained simulation is 23% to 1.0 Å resolution. The B-factors for most heavy atoms agree well with experiment (Pearson correlation of 0.90), but B-factors obtained by refinement against the average simulation density underestimate the coordinate fluctuations in the underlying simulation where the simulation samples alternate conformations. A dynamic flow of water molecules through channels within the crystal lattice is observed, yet the average water density is in remarkable agreement with experiment. A minor population of unit cells is characterized by reduced water content, 310 helical propensity and a gauche(−) side-chain rotamer for one of the valine residues. Careful examination of the experimental data suggests that transitions of the helices are a simulation artifact, although there is indeed evidence for alternate valine conformers and variable water content. This study highlights the potential for crystal simulations to detect dynamics and heterogeneity in experimental diffraction data, as well as to validate computational chemistry methods. PMID:23631449
Methodology for Uncertainty Analysis of Dynamic Computational Toxicology Models
The task of quantifying the uncertainty in both parameter estimates and model predictions has become more important with the increased use of dynamic computational toxicology models by the EPA. Dynamic toxicological models include physiologically-based pharmacokinetic (PBPK) mode...
Buckybomb: Reactive Molecular Dynamics Simulation.
Chaban, Vitaly V; Fileti, Eudes Eterno; Prezhdo, Oleg V
2015-03-01
Energetic materials, such as explosives, propellants, and pyrotechnics, are widely used in civilian and military applications. Nanoscale explosives represent a special group because of the high density of energetic covalent bonds. The reactive molecular dynamics (ReaxFF) study of nitrofullerene decomposition reported here provides a detailed chemical mechanism of explosion of a nanoscale carbon material. Upon initial heating, C60(NO2)12 disintegrates, increasing temperature and pressure by thousands of Kelvins and bars within tens of picoseconds. The explosion starts with NO2 group isomerization into C-O-N-O, followed by emission of NO molecules and formation of CO groups on the buckyball surface. NO oxidizes into NO2, and C60 falls apart, liberating CO2. At the highest temperatures, CO2 gives rise to diatomic carbon. The study shows that the initiation temperature and released energy depend strongly on the chemical composition and density of the material. PMID:26262672
Dynamic procedure for filtered gyrokinetic simulations
Morel, P.; Banon Navarro, A.; Albrecht-Marc, M.; Carati, D.; Merz, F.; Goerler, T.; Jenko, F.
2012-01-15
Large eddy simulations (LES) of gyrokinetic plasma turbulence are investigated as interesting candidates to decrease the computational cost. A dynamic procedure is implemented in the gene code, allowing for dynamic optimization of the free parameters of the LES models (setting the amplitudes of dissipative terms). Employing such LES methods, one recovers the free energy and heat flux spectra obtained from highly resolved direct numerical simulations. Systematic comparisons are performed for different values of the temperature gradient and magnetic shear, parameters which are of prime importance in ion temperature gradient driven turbulence. Moreover, the degree of anisotropy of the problem, which can vary with parameters, can be adapted dynamically by the method that shows gyrokinetic large eddy simulation to be a serious candidate to reduce numerical cost of gyrokinetic solvers.
Computer simulation of microstructural dynamics
Grest, G.S.; Anderson, M.P.; Srolovitz, D.J.
1985-01-01
Since many of the physical properties of materials are determined by their microstructure, it is important to be able to predict and control microstructural development. A number of approaches have been taken to study this problem, but they assume that the grains can be described as spherical or hexagonal and that growth occurs in an average environment. We have developed a new technique to bridge the gap between the atomistic interactions and the macroscopic scale by discretizing the continuum system such that the microstructure retains its topological connectedness, yet is amenable to computer simulations. Using this technique, we have studied grain growth in polycrystalline aggregates. The temporal evolution and grain morphology of our model are in excellent agreement with experimental results for metals and ceramics.
Dynamic Fracture Simulations of Explosively Loaded Cylinders
Arthur, Carly W.; Goto, D. M.
2015-11-30
This report documents the modeling results of high explosive experiments investigating dynamic fracture of steel (AerMet® 100 alloy) cylinders. The experiments were conducted at Lawrence Livermore National Laboratory (LLNL) during 2007 to 2008 [10]. A principal objective of this study was to gain an understanding of dynamic material failure through the analysis of hydrodynamic computer code simulations. Two-dimensional and three-dimensional computational cylinder models were analyzed using the ALE3D multi-physics computer code.
Dynamic risk simulation to assess natural hazards risk along roads
NASA Astrophysics Data System (ADS)
Voumard, J.; Caspar, O.; Derron, M.-H.; Jaboyedoff, M.
2013-11-01
Risk generated by natural hazards on roads is usually calculated with equations integrating various parameters related to hazard and traffic. These are static variables, like an average number of vehicles crossing this section every day and an average vehicle speed. This methodology cannot take into account dynamic variations of traffic and interactions between vehicles such as speed modifications due to windy roads, slowdowns resulting from saturated traffic or vehicle tailbacks forming in front of traffic lights. Here we show, by means of a dynamic traffic simulator, that traffic variations may greatly influence the risk estimation over time. The risk is analysed on several sections of an Alpine road in Switzerland using a dynamic vehicles approach, and compared with the results of the static methodology. It demonstrates that risk can significantly increase on sinuous sections because of decreasing vehicle speed. For example, along an 800 m-long section of road containing two hairpin bends, the dynamic risk is about 50% higher than the static one. Badly placed signalization, slowing down, or stopping the vehicles in a hazardous area may increase the risk by about 150% (i.e. 2.5 times higher) along a straight road section where vehicles speed is high. A more realistic risk can thus be obtained from a dynamic approach, especially on mountain roads. The dynamic traffic simulator developed for this work appears to be a helpful tool to support decision-making in reducing risk on mountain roads and it shows the importance of keeping the traffic moving as freely as possible.
Molecular dynamics simulations through GPU video games technologies
Loukatou, Styliani; Papageorgiou, Louis; Fakourelis, Paraskevas; Filntisi, Arianna; Polychronidou, Eleftheria; Bassis, Ioannis; Megalooikonomou, Vasileios; Makałowski, Wojciech; Vlachakis, Dimitrios; Kossida, Sophia
2016-01-01
Bioinformatics is the scientific field that focuses on the application of computer technology to the management of biological information. Over the years, bioinformatics applications have been used to store, process and integrate biological and genetic information, using a wide range of methodologies. One of the most de novo techniques used to understand the physical movements of atoms and molecules is molecular dynamics (MD). MD is an in silico method to simulate the physical motions of atoms and molecules under certain conditions. This has become a state strategic technique and now plays a key role in many areas of exact sciences, such as chemistry, biology, physics and medicine. Due to their complexity, MD calculations could require enormous amounts of computer memory and time and therefore their execution has been a big problem. Despite the huge computational cost, molecular dynamics have been implemented using traditional computers with a central memory unit (CPU). A graphics processing unit (GPU) computing technology was first designed with the goal to improve video games, by rapidly creating and displaying images in a frame buffer such as screens. The hybrid GPU-CPU implementation, combined with parallel computing is a novel technology to perform a wide range of calculations. GPUs have been proposed and used to accelerate many scientific computations including MD simulations. Herein, we describe the new methodologies developed initially as video games and how they are now applied in MD simulations. PMID:27525251
Dynamic simulation recalls condensate piping event
Farrell, R.J.; Reneberg, K.O. ); Moy, H.C. )
1994-05-01
This article describes how experience gained from simulating and reconstructing a condensate piping event will be used by Consolidated Edison to analyze control system problems. A cooperative effort by Con Edison and the Chemical Engineering Department at Polytechnic University used modular modeling system to investigate the probable cause of a Con Edison condensate piping event. Con Edison commissioned the work to serve as a case study for the more general problem of control systems analysis using dynamic simulation and MMS.
Molecular dynamics simulation of triclinic lysozyme in a crystal lattice.
Janowski, Pawel A; Liu, Chunmei; Deckman, Jason; Case, David A
2016-01-01
Molecular dynamics simulations of crystals can enlighten interpretation of experimental X-ray crystallography data and elucidate structural dynamics and heterogeneity in biomolecular crystals. Furthermore, because of the direct comparison against experimental data, they can inform assessment of molecular dynamics methods and force fields. We present microsecond scale results for triclinic hen egg-white lysozyme in a supercell consisting of 12 independent unit cells using four contemporary force fields (Amber ff99SB, ff14ipq, ff14SB, and CHARMM 36) in crystalline and solvated states (for ff14SB only). We find the crystal simulations consistent across multiple runs of the same force field and robust to various solvent equilibration schemes. However, convergence is slow compared with solvent simulations. All the tested force fields reproduce experimental structural and dynamic properties well, but Amber ff14SB maintains structure and reproduces fluctuations closest to the experimental model: its average backbone structure differs from the deposited structure by 0.37Å; by contrast, the average backbone structure in solution differs from the deposited by 0.65Å. All the simulations are affected by a small progressive deterioration of the crystal lattice, presumably due to imperfect modeling of hydrogen bonding and other crystal contact interactions; this artifact is smallest in ff14SB, with average lattice positions deviating by 0.20Å from ideal. Side-chain disorder is surprisingly low with fewer than 30% of the nonglycine or alanine residues exhibiting significantly populated alternate rotamers. Our results provide helpful insight into the methodology of biomolecular crystal simulations and indicate directions for future work to obtain more accurate energy models for molecular dynamics. PMID:26013419
Molecular dynamic simulations of ocular tablet dissolution.
Ru, Qian; Fadda, Hala M; Li, Chung; Paul, Daniel; Khaw, Peng T; Brocchini, Steve; Zloh, Mire
2013-11-25
Small tablets for implantation into the subconjunctival space in the eye are being developed to inhibit scarring after glaucoma filtration surgery (GFS). There is a need to evaluate drug dissolution at the molecular level to determine how the chemical structure of the active may correlate with dissolution in the nonsink conditions of the conjunctival space. We conducted molecular dynamics simulations to study the dissolution process of tablets derived from two drugs that can inhibit fibrosis after GFS, 5-fluorouracil (5-FU) and the matrix metalloprotease inhibitor (MMPi), ilomastat. The dissolution was simulated in the presence of simple point charge (SPC) water molecules, and the liquid turnover of the aqueous humor in the subconjunctival space was simulated by removal of the dissolved drug molecules at regular intervals and replacement by new water molecules. At the end of the simulation, the total molecular solvent accessible surface area of 5-FU tablets increased by 60 times more than that of ilomastat as a result of tablet swelling and release of molecules into solution. The tablet dissolution pattern shown in our molecular dynamic simulations tends to correlate with experimental release profiles. This work indicates that a series of molecular dynamic simulations can be used to predict the influence of the molecular properties of a drug on its dissolution profile and could be useful during preformulation where sufficient amounts of the drug are not always available to perform dissolution studies. PMID:24073784
Airborne Simulation of Launch Vehicle Dynamics
NASA Technical Reports Server (NTRS)
Miller, Christopher J.; Orr, Jeb S.; Hanson, Curtis E.; Gilligan, Eric T.
2015-01-01
In this paper we present a technique for approximating the short-period dynamics of an exploration-class launch vehicle during flight test with a high-performance surrogate aircraft in relatively benign endoatmospheric flight conditions. The surrogate vehicle relies upon a nonlinear dynamic inversion scheme with proportional-integral feedback to drive a subset of the aircraft states into coincidence with the states of a time-varying reference model that simulates the unstable rigid body dynamics, servodynamics, and parasitic elastic and sloshing dynamics of the launch vehicle. The surrogate aircraft flies a constant pitch rate trajectory to approximate the boost phase gravity turn ascent, and the aircraft's closed-loop bandwidth is sufficient to simulate the launch vehicle's fundamental lateral bending and sloshing modes by exciting the rigid body dynamics of the aircraft. A novel control allocation scheme is employed to utilize the aircraft's relatively fast control effectors in inducing various failure modes for the purposes of evaluating control system performance. Sufficient dynamic similarity is achieved such that the control system under evaluation is configured for the full-scale vehicle with no changes to its parameters, and pilot-control system interaction studies can be performed to characterize the effects of guidance takeover during boost. High-fidelity simulation and flight-test results are presented that demonstrate the efficacy of the design in simulating the Space Launch System (SLS) launch vehicle dynamics using the National Aeronautics and Space Administration (NASA) Armstrong Flight Research Center Fullscale Advanced Systems Testbed (FAST), a modified F/A-18 airplane (McDonnell Douglas, now The Boeing Company, Chicago, Illinois), over a range of scenarios designed to stress the SLS's Adaptive Augmenting Control (AAC) algorithm.
Airborne Simulation of Launch Vehicle Dynamics
NASA Technical Reports Server (NTRS)
Gilligan, Eric T.; Miller, Christopher J.; Hanson, Curtis E.; Orr, Jeb S.
2014-01-01
In this paper we present a technique for approximating the short-period dynamics of an exploration-class launch vehicle during flight test with a high-performance surrogate aircraft in relatively benign endoatmospheric flight conditions. The surrogate vehicle relies upon a nonlinear dynamic inversion scheme with proportional-integral feedback to drive a subset of the aircraft states into coincidence with the states of a time-varying reference model that simulates the unstable rigid body dynamics, servodynamics, and parasitic elastic and sloshing dynamics of the launch vehicle. The surrogate aircraft flies a constant pitch rate trajectory to approximate the boost phase gravity-turn ascent, and the aircraft's closed-loop bandwidth is sufficient to simulate the launch vehicle's fundamental lateral bending and sloshing modes by exciting the rigid body dynamics of the aircraft. A novel control allocation scheme is employed to utilize the aircraft's relatively fast control effectors in inducing various failure modes for the purposes of evaluating control system performance. Sufficient dynamic similarity is achieved such that the control system under evaluation is optimized for the full-scale vehicle with no changes to its parameters, and pilot-control system interaction studies can be performed to characterize the effects of guidance takeover during boost. High-fidelity simulation and flight test results are presented that demonstrate the efficacy of the design in simulating the Space Launch System (SLS) launch vehicle dynamics using NASA Dryden Flight Research Center's Full-scale Advanced Systems Testbed (FAST), a modified F/A-18 airplane, over a range of scenarios designed to stress the SLS's adaptive augmenting control (AAC) algorithm.
Fully dynamical simulation of central nuclear collisions.
van der Schee, Wilke; Romatschke, Paul; Pratt, Scott
2013-11-27
We present a fully dynamical simulation of central nuclear collisions around midrapidity at LHC energies. Unlike previous treatments, we simulate all phases of the collision, including the equilibration of the system. For the simulation, we use numerical relativity solutions to anti-de Sitter space/conformal field theory for the preequilibrium stage, viscous hydrodynamics for the plasma equilibrium stage, and kinetic theory for the low-density hadronic stage. Our preequilibrium stage provides initial conditions for hydrodynamics, resulting in sizable radial flow. The resulting light particle spectra reproduce the measurements from the ALICE experiment at all transverse momenta. PMID:24329444
Molecular dynamics simulation of ice XII
NASA Astrophysics Data System (ADS)
Borzsák, István; Cummings, Peter T.
1999-02-01
Molecular dynamics simulations have been performed on the newly discovered metastable ice XII. This new crystalline ice phase [C. Lobban, J.L. Finney, W.F. Kuhs, Nature (London) 391 (1998) 268] is proton-disordered. Thus 90 possible configurations of the unit cell can be constructed which differ only in the orientations of the water molecules. The simulation used the TIP4P potential model for water at constant temperature and density. About one-quarter of the initial configurations did not melt in the course of the simulation. This result is supportive of the experimental structure and also demonstrates the ability of this water model to study ice phases.
Hybrid molecular-continuum simulations using smoothed dissipative particle dynamics
Petsev, Nikolai D.; Leal, L. Gary; Shell, M. Scott
2015-01-28
We present a new multiscale simulation methodology for coupling a region with atomistic detail simulated via molecular dynamics (MD) to a numerical solution of the fluctuating Navier-Stokes equations obtained from smoothed dissipative particle dynamics (SDPD). In this approach, chemical potential gradients emerge due to differences in resolution within the total system and are reduced by introducing a pairwise thermodynamic force inside the buffer region between the two domains where particles change from MD to SDPD types. When combined with a multi-resolution SDPD approach, such as the one proposed by Kulkarni et al. [J. Chem. Phys. 138, 234105 (2013)], this method makes it possible to systematically couple atomistic models to arbitrarily coarse continuum domains modeled as SDPD fluids with varying resolution. We test this technique by showing that it correctly reproduces thermodynamic properties across the entire simulation domain for a simple Lennard-Jones fluid. Furthermore, we demonstrate that this approach is also suitable for non-equilibrium problems by applying it to simulations of the start up of shear flow. The robustness of the method is illustrated with two different flow scenarios in which shear forces act in directions parallel and perpendicular to the interface separating the continuum and atomistic domains. In both cases, we obtain the correct transient velocity profile. We also perform a triple-scale shear flow simulation where we include two SDPD regions with different resolutions in addition to a MD domain, illustrating the feasibility of a three-scale coupling.
Hybrid molecular-continuum simulations using smoothed dissipative particle dynamics
NASA Astrophysics Data System (ADS)
Petsev, Nikolai D.; Leal, L. Gary; Shell, M. Scott
2015-01-01
We present a new multiscale simulation methodology for coupling a region with atomistic detail simulated via molecular dynamics (MD) to a numerical solution of the fluctuating Navier-Stokes equations obtained from smoothed dissipative particle dynamics (SDPD). In this approach, chemical potential gradients emerge due to differences in resolution within the total system and are reduced by introducing a pairwise thermodynamic force inside the buffer region between the two domains where particles change from MD to SDPD types. When combined with a multi-resolution SDPD approach, such as the one proposed by Kulkarni et al. [J. Chem. Phys. 138, 234105 (2013)], this method makes it possible to systematically couple atomistic models to arbitrarily coarse continuum domains modeled as SDPD fluids with varying resolution. We test this technique by showing that it correctly reproduces thermodynamic properties across the entire simulation domain for a simple Lennard-Jones fluid. Furthermore, we demonstrate that this approach is also suitable for non-equilibrium problems by applying it to simulations of the start up of shear flow. The robustness of the method is illustrated with two different flow scenarios in which shear forces act in directions parallel and perpendicular to the interface separating the continuum and atomistic domains. In both cases, we obtain the correct transient velocity profile. We also perform a triple-scale shear flow simulation where we include two SDPD regions with different resolutions in addition to a MD domain, illustrating the feasibility of a three-scale coupling.
Development of methodology for horizontal axis wind turbine dynamic analysis
NASA Technical Reports Server (NTRS)
Dugundji, J.
1982-01-01
Horizontal axis wind turbine dynamics were studied. The following findings are summarized: (1) review of the MOSTAS computer programs for dynamic analysis of horizontal axis wind turbines; (2) review of various analysis methods for rotating systems with periodic coefficients; (3) review of structural dynamics analysis tools for large wind turbine; (4) experiments for yaw characteristics of a rotating rotor; (5) development of a finite element model for rotors; (6) development of simple models for aeroelastics; and (7) development of simple models for stability and response of wind turbines on flexible towers.
Digital simulation of stiff linear dynamic systems.
NASA Technical Reports Server (NTRS)
Holland, L. D.; Walsh, J. R., Jr.; Kerr, J. H.
1972-01-01
A method is derived for digital computer simulation of linear time-invariant systems when the insignificant eigenvalues involved in such systems are eliminated by an ALSAP root removal technique. The method is applied to a thirteenth-order dynamic system representing a passive RLC network.
Dynamic simulation of a reverse Brayton refrigerator
Peng, N.; Xiong, L. Y.; Dong, B.; Liu, L. Q.; Lei, L. L.; Tang, J. C.
2014-01-29
A test refrigerator based on the modified Reverse Brayton cycle has been developed in the Chinese Academy of Sciences recently. To study the behaviors of this test refrigerator, a dynamic simulation has been carried out. The numerical model comprises the typical components of the test refrigerator: compressor, valves, heat exchangers, expander and heater. This simulator is based on the oriented-object approach and each component is represented by a set of differential and algebraic equations. The control system of the test refrigerator is also simulated, which can be used to optimize the control strategies. This paper describes all the models and shows the simulation results. Comparisons between simulation results and experimental data are also presented. Experimental validation on the test refrigerator gives satisfactory results.
NASA Technical Reports Server (NTRS)
Pamadi, Bandu N.; Toniolo, Matthew D.; Tartabini, Paul V.; Roithmayr, Carlos M.; Albertson, Cindy W.; Karlgaard, Christopher D.
2016-01-01
The objective of this report is to develop and implement a physics based method for analysis and simulation of multi-body dynamics including launch vehicle stage separation. The constraint force equation (CFE) methodology discussed in this report provides such a framework for modeling constraint forces and moments acting at joints when the vehicles are still connected. Several stand-alone test cases involving various types of joints were developed to validate the CFE methodology. The results were compared with ADAMS(Registered Trademark) and Autolev, two different industry standard benchmark codes for multi-body dynamic analysis and simulations. However, these two codes are not designed for aerospace flight trajectory simulations. After this validation exercise, the CFE algorithm was implemented in Program to Optimize Simulated Trajectories II (POST2) to provide a capability to simulate end-to-end trajectories of launch vehicles including stage separation. The POST2/CFE methodology was applied to the STS-1 Space Shuttle solid rocket booster (SRB) separation and Hyper-X Research Vehicle (HXRV) separation from the Pegasus booster as a further test and validation for its application to launch vehicle stage separation problems. Finally, to demonstrate end-to-end simulation capability, POST2/CFE was applied to the ascent, orbit insertion, and booster return of a reusable two-stage-to-orbit (TSTO) vehicle concept. With these validation exercises, POST2/CFE software can be used for performing conceptual level end-to-end simulations, including launch vehicle stage separation, for problems similar to those discussed in this report.
ADAPTIVE MULTILEVEL SPLITTING IN MOLECULAR DYNAMICS SIMULATIONS*
Aristoff, David; Lelièvre, Tony; Mayne, Christopher G.; Teo, Ivan
2014-01-01
Adaptive Multilevel Splitting (AMS) is a replica-based rare event sampling method that has been used successfully in high-dimensional stochastic simulations to identify trajectories across a high potential barrier separating one metastable state from another, and to estimate the probability of observing such a trajectory. An attractive feature of AMS is that, in the limit of a large number of replicas, it remains valid regardless of the choice of reaction coordinate used to characterize the trajectories. Previous studies have shown AMS to be accurate in Monte Carlo simulations. In this study, we extend the application of AMS to molecular dynamics simulations and demonstrate its effectiveness using a simple test system. Our conclusion paves the way for useful applications, such as molecular dynamics calculations of the characteristic time of drug dissociation from a protein target. PMID:26005670
Test of a flexible spacecraft dynamics simulator
NASA Technical Reports Server (NTRS)
Dichmann, Donald; Sedlak, Joseph
1998-01-01
There are a number of approaches one can take to modeling the dynamics of a flexible body. While one can attempt to capture the full dynamical behavior subject to disturbances from actuators and environmental torques, such a detailed description often is unnecessary. Simplification is possible either by limiting the amplitude of motion to permit linearization of the dynamics equations or by restricting the types of allowed motion. In this work, we study the nonlinear dynamics of bending deformations of wire booms on spinning spacecraft. The theory allows for large amplitude excursions from equilibrium while enforcing constraints on the dynamics to prohibit those modes that are physically less relevant or are expected to damp out fast. These constraints explicitly remove the acoustic modes (i.e., longitudinal sound waves and shear waves) while allowing for arbitrary bending and twisting, motions which typically are of lower frequency. As a test case, a spin axis reorientation maneuver by the Polar Plasma Laboratory (POLAR) spacecraft has been simulated. POLAR was chosen as a representative spacecraft because it has flexible wire antennas that extend to a length of 65 meters. Bending deformations in these antennas could be quite large and have a significant effect on the attitude dynamics of the spacecraft body. Summary results from the simulation are presented along, with a comparison with POLAR flight data.
A Comparative Study of Three Methodologies for Modeling Dynamic Stall
NASA Technical Reports Server (NTRS)
Sankar, L.; Rhee, M.; Tung, C.; ZibiBailly, J.; LeBalleur, J. C.; Blaise, D.; Rouzaud, O.
2002-01-01
During the past two decades, there has been an increased reliance on the use of computational fluid dynamics methods for modeling rotors in high speed forward flight. Computational methods are being developed for modeling the shock induced loads on the advancing side, first-principles based modeling of the trailing wake evolution, and for retreating blade stall. The retreating blade dynamic stall problem has received particular attention, because the large variations in lift and pitching moments encountered in dynamic stall can lead to blade vibrations and pitch link fatigue. Restricting to aerodynamics, the numerical prediction of dynamic stall is still a complex and challenging CFD problem, that, even in two dimensions at low speed, gathers the major difficulties of aerodynamics, such as the grid resolution requirements for the viscous phenomena at leading-edge bubbles or in mixing-layers, the bias of the numerical viscosity, and the major difficulties of the physical modeling, such as the turbulence models, the transition models, whose both determinant influences, already present in static maximal-lift or stall computations, are emphasized by the dynamic aspect of the phenomena.
Two-Speed Gearbox Dynamic Simulation Predictions and Test Validation
NASA Technical Reports Server (NTRS)
Lewicki, David G.; DeSmidt, Hans; Smith, Edward C.; Bauman, Steven W.
2010-01-01
Dynamic simulations and experimental validation tests were performed on a two-stage, two-speed gearbox as part of the drive system research activities of the NASA Fundamental Aeronautics Subsonics Rotary Wing Project. The gearbox was driven by two electromagnetic motors and had two electromagnetic, multi-disk clutches to control output speed. A dynamic model of the system was created which included a direct current electric motor with proportional-integral-derivative (PID) speed control, a two-speed gearbox with dual electromagnetically actuated clutches, and an eddy current dynamometer. A six degree-of-freedom model of the gearbox accounted for the system torsional dynamics and included gear, clutch, shaft, and load inertias as well as shaft flexibilities and a dry clutch stick-slip friction model. Experimental validation tests were performed on the gearbox in the NASA Glenn gear noise test facility. Gearbox output speed and torque as well as drive motor speed and current were compared to those from the analytical predictions. The experiments correlate very well with the predictions, thus validating the dynamic simulation methodologies.
Simulation studies using multibody dynamics code DART
NASA Technical Reports Server (NTRS)
Keat, James E.
1989-01-01
DART is a multibody dynamics code developed by Photon Research Associates for the Air Force Astronautics Laboratory (AFAL). The code is intended primarily to simulate the dynamics of large space structures, particularly during the deployment phase of their missions. DART integrates nonlinear equations of motion numerically. The number of bodies in the system being simulated is arbitrary. The bodies' interconnection joints can have an arbitrary number of degrees of freedom between 0 and 6. Motions across the joints can be large. Provision for simulating on-board control systems is provided. Conservation of energy and momentum, when applicable, are used to evaluate DART's performance. After a brief description of DART, studies made to test the program prior to its delivery to AFAL are described. The first is a large angle reorientating of a flexible spacecraft consisting of a rigid central hub and four flexible booms. Reorientation was accomplished by a single-cycle sine wave shape torque input. In the second study, an appendage, mounted on a spacecraft, was slewed through a large angle. Four closed-loop control systems provided control of this appendage and of the spacecraft's attitude. The third study simulated the deployment of the rim of a bicycle wheel configuration large space structure. This system contained 18 bodies. An interesting and unexpected feature of the dynamics was a pulsing phenomena experienced by the stays whole playout was used to control the deployment. A short description of the current status of DART is given.
Mesoscopic Simulation Methods for Polymer Dynamics
NASA Astrophysics Data System (ADS)
Larson, Ronald
2015-03-01
We assess the accuracy and efficiency of mesoscopic simulation methods, namely Brownian Dynamics (BD), Stochastic Rotation Dynamics (SRD) and Dissipative Particle Dynamics (DPD), for polymers in solution at equilibrium and in flows in microfluidic geometries. Both SRD and DPD use solvent ``particles'' to carry momentum, and so account automatically for hydrodynamic interactions both within isolated polymer coils, and with other polymer molecules and with nearby solid boundaries. We assess quantitatively the effects of artificial particle inertia and fluid compressibility and show that they can be made small with appropriate choice of simulation parameters. We then use these methods to study flow-induced migration of polymer chains produced by: 1) hydrodynamic interactions, 2) streamline curvature or stress-gradients, and 3) convection of wall depletion zones. We show that huge concentration gradients can be produced by these mechanisms in microfluidic geometries that can be exploited for separation of polymers by size in periodic contraction-expansion geometries. We also assess the range of conditions for which BD, SRD or DPD is preferable for mesoscopic simulations. Finally, we show how such methods can be used to simulate quantitatively the swimming of micro-organisms such as E. coli. In collaboration with Lei Jiang and Tongyang Zhao, University of Michigan, Ann Arbor, MI.
NASA Astrophysics Data System (ADS)
Paul, Pijush Kanti
In the fault damage zone modeling study for a field in the Timor Sea, I present a methodology to incorporate geomechanically-based fault damage zones into reservoir simulation. In the studied field, production history suggests that the mismatch between actual production and model prediction is due to preferential fluid flow through the damage zones associated with the reservoir scale faults, which is not included in the baseline petrophysical model. I analyzed well data to estimate stress heterogeneity and fracture distributions in the reservoir. Image logs show that stress orientations are homogenous at the field scale with a strike-slip/normal faulting stress regime and maximum horizontal stress oriented in NE-SW direction. Observed fracture zones in wells are mostly associated with well scale fault and bed boundaries. These zones do not show any anomalies in production logs or well test data, because most of the fractures are not optimally oriented to the present day stress state, and matrix permeability is high enough to mask any small anomalies from the fracture zones. However, I found that fracture density increases towards the reservoir scale faults, indicating high fracture density zones or damage zones close to these faults, which is consistent with the preferred flow direction indicated by interference and tracer test done between the wells. It is well known from geologic studies that there is a concentration of secondary fractures and faults in a damage zone adjacent to larger faults. Because there is usually inadequate data to incorporate damage zone fractures and faults into reservoir simulation models, in this study I utilized the principles of dynamic rupture propagation from earthquake seismology to predict the nature of fractured/damage zones associated with reservoir scale faults. The implemented workflow can be used to more routinely incorporate damage zones into reservoir simulation models. Applying this methodology to a real reservoir utilizing
Dynamic simulation of the mastication muscles
NASA Astrophysics Data System (ADS)
Weingaertner, Tim; Albrecht, Jochen
1998-05-01
The purpose of a simulated operation system in craniofacial surgery is to evaluate and visualize the results of operations on the overall facial shape of the patient and on the functionality of his jaw. This paper presents the analyzation of muscle movements in the mastication system by applying real jaw movements to the simulation. With this method an accurate modeling of the mastication muscles can be performed which is a prerequisite for a realistic simulation and precise intra- operative registration. According to this results a large- scale musculoskeletal model of the mastication system is generated including kinematic and dynamic parameters. By integrating distance sensors in the simulation of a segmented CT (computer tomograph) image of the maxilla and mandible the motions of the masticatory muscles during different kinds of jaw movements have been analyzed. The data for this motions have been recorded by a commercial system (CONDYLOCOMP LR3) on a test person and transformed to the graphical simulation system. This method for the first time allows to observe the dynamics of the mastication muscles and their different parameters like muscle length ratio and velocity. The integration of a kinematic model for the jaw movement makes it possible to analyze non traced movements.
Dynamic Simulation of a Helium Liquefier
Maekawa, R.; Ooba, K.; Mito, T.; Nobutoki, M.
2004-06-23
Dynamic behavior of a helium liquefier has been studied in detail with a Cryogenic Process REal-time SimulaTor (C-PREST) at the National Institute for Fusion Science (NIFS). The C-PREST is being developed to integrate large-scale helium cryogenic plant design, operation and maintenance for optimum process establishment. As a first step of simulations of cooldown to 4.5 K with the helium liquefier model is conducted, which provides a plant-process validation platform. The helium liquefier consists of seven heat exchangers, a liquid-nitrogen (LN2) precooler, two expansion turbines and a liquid-helium (LHe) reservoir. Process simulations are fulfilled with sequence programs, which were implemented with C-PREST based on an existing liquefier operation. The interactions of a JT valve, a JT-bypass valve and a reservoir-return valve have been dynamically simulated. The paper discusses various aspects of refrigeration process simulation, including its difficulties such as a balance between complexity of the adopted models and CPU time.
Methodological advances: using greenhouses to simulate climate change scenarios.
Morales, F; Pascual, I; Sánchez-Díaz, M; Aguirreolea, J; Irigoyen, J J; Goicoechea, N; Antolín, M C; Oyarzun, M; Urdiain, A
2014-09-01
Human activities are increasing atmospheric CO2 concentration and temperature. Related to this global warming, periods of low water availability are also expected to increase. Thus, CO2 concentration, temperature and water availability are three of the main factors related to climate change that potentially may influence crops and ecosystems. In this report, we describe the use of growth chamber - greenhouses (GCG) and temperature gradient greenhouses (TGG) to simulate climate change scenarios and to investigate possible plant responses. In the GCG, CO2 concentration, temperature and water availability are set to act simultaneously, enabling comparison of a current situation with a future one. Other characteristics of the GCG are a relative large space of work, fine control of the relative humidity, plant fertirrigation and the possibility of light supplementation, within the photosynthetic active radiation (PAR) region and/or with ultraviolet-B (UV-B) light. In the TGG, the three above-mentioned factors can act independently or in interaction, enabling more mechanistic studies aimed to elucidate the limiting factor(s) responsible for a given plant response. Examples of experiments, including some aimed to study photosynthetic acclimation, a phenomenon that leads to decreased photosynthetic capacity under long-term exposures to elevated CO2, using GCG and TGG are reported. PMID:25113448
Thermostability of Enzymes from Molecular Dynamics Simulations.
Zeiske, Tim; Stafford, Kate A; Palmer, Arthur G
2016-06-14
Thermodynamic stability is a central requirement for protein function, and one goal of protein engineering is improvement of stability, particularly for applications in biotechnology. Herein, molecular dynamics simulations are used to predict in vitro thermostability of members of the bacterial ribonuclease HI (RNase H) family of endonucleases. The temperature dependence of the generalized order parameter, S, for four RNase H homologues, from psychrotrophic, mesophilic, and thermophilic organisms, is highly correlated with experimentally determined melting temperatures and with calculated free energies of folding at the midpoint temperature of the simulations. This study provides an approach for in silico mutational screens to improve thermostability of biologically and industrially relevant enzymes. PMID:27123810
A New Concurrent Multiscale Methodology for Coupling Molecular Dynamics and Finite Element Analyses
NASA Technical Reports Server (NTRS)
Yamakov, Vesselin; Saether, Erik; Glaessgen, Edward H/.
2008-01-01
The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.
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.
Nanoindentation of Zr by molecular dynamics simulation
NASA Astrophysics Data System (ADS)
Lu (芦子哲), Zizhe; Chernatynskiy, Aleksandr; Noordhoek, Mark J.; Sinnott, Susan B.; Phillpot, Simon R.
2015-12-01
Molecular dynamics simulations of nanoindentation are used to study the deformation behaviors of single crystal Zr for four different surface orientations. The comparison of results for two different potentials, an embedded atom method potential and a charged optimized many body potential, reveals the influence of stable and unstable stacking fault energy on dislocation behaviors under nanoindentation. The load-displacement curve, hardness and deformation behaviors of the various surface orientations Zr are compared and the elastic and plastic deformation behaviors are analyzed.
Numerical Simulations of Ion Cloud Dynamics
NASA Astrophysics Data System (ADS)
Sillitoe, Nicolas; Hilico, Laurent
We explain how to perform accurate numerical simulations of ion cloud dynamics by discussing the relevant orders of magnitude of the characteristic times and frequencies involved in the problem and the computer requirement with respect to the ion cloud size. We then discuss integration algorithms and Coulomb force parallelization. We finally explain how to take into account collisions, cooling laser interaction and chemical reactions in a Monte Carlo approach and discuss how to use random number generators to that end.
Twisdale, L.A.; Dunn, W.L.
1981-08-01
A probabilistic methodology has been developed to predict the probabilities of tornado-propelled missiles impacting and damaging nuclear power plant structures. Mathematical models of each event in the tornado missile hazard have been developed and sequenced to form an integrated, time-history simulation methodology. The models are data based where feasible. The data include documented records of tornado occurrence, field observations of missile transport, results of wind tunnel experiments, and missile impact tests. Probabilistic Monte Carlo techniques are used to estimate the risk probabilities. The methodology has been encoded in the TORMIS computer code to facilitate numerical analysis and plant-specific tornado missile probability assessments.
NASA Technical Reports Server (NTRS)
Chidester, Thomas R.; Kanki, Barbara G.; Helmreich, Robert L.
1989-01-01
The crew-factors research program at NASA Ames has developed a methodology for studying the impact of a variety of variables on the effectiveness of crews flying realistic but high workload simulated trips. The validity of investigations using the methodology is enhanced by careful design of full-mission scenarios, performance assessment using converging sources of data, and recruitment of representative subjects. Recently, portions of this methodology have been adapted for use in assessing the effectiveness of crew coordination among participants in line-oriented flight training.
Simulation of counterflow pedestrian dynamics using spheropolygons
NASA Astrophysics Data System (ADS)
Alonso-Marroquín, Fernando; Busch, Jonathan; Chiew, Coraline; Lozano, Celia; Ramírez-Gómez, Álvaro
2014-12-01
Pedestrian dynamic models are typically designed for comfortable walking or slightly congested conditions and typically use a single disk or combination of three disks for the shape of a pedestrian. Under crowd conditions, a more accurate pedestrian shape has advantages over the traditional single or three-disks model. We developed a method for simulating pedestrian dynamics in a large dense crowd of spheropolygons adapted to the cross section of the chest and arms of a pedestrian. Our numerical model calculates pedestrian motion from Newton's second law, taking into account viscoelastic contact forces, contact friction, and ground-reaction forces. Ground-reaction torque was taken to arise solely from the pedestrians' orientation toward their preferred destination. Simulations of counterflow pedestrians dynamics in corridors were used to gain insight into a tragic incident at the Madrid Arena pavilion in Spain, where five girls were crushed to death. The incident took place at a Halloween Celebration in 2012, in a long, densely crowded hallway used as entrance and exit at the same time. Our simulations reconstruct the mechanism of clogging in the hallway. The hypothetical case of a total evacuation order was also investigated. The results highlights the importance of the pedestrians' density and the effect of counterflow in the onset of avalanches and clogging and provides an estimation of the number of injuries based on a calculation of the contact-force network between the pedestrians.
Simulation of counterflow pedestrian dynamics using spheropolygons.
Alonso-Marroquín, Fernando; Busch, Jonathan; Chiew, Coraline; Lozano, Celia; Ramírez-Gómez, Álvaro
2014-12-01
Pedestrian dynamic models are typically designed for comfortable walking or slightly congested conditions and typically use a single disk or combination of three disks for the shape of a pedestrian. Under crowd conditions, a more accurate pedestrian shape has advantages over the traditional single or three-disks model. We developed a method for simulating pedestrian dynamics in a large dense crowd of spheropolygons adapted to the cross section of the chest and arms of a pedestrian. Our numerical model calculates pedestrian motion from Newton's second law, taking into account viscoelastic contact forces, contact friction, and ground-reaction forces. Ground-reaction torque was taken to arise solely from the pedestrians' orientation toward their preferred destination. Simulations of counterflow pedestrians dynamics in corridors were used to gain insight into a tragic incident at the Madrid Arena pavilion in Spain, where five girls were crushed to death. The incident took place at a Halloween Celebration in 2012, in a long, densely crowded hallway used as entrance and exit at the same time. Our simulations reconstruct the mechanism of clogging in the hallway. The hypothetical case of a total evacuation order was also investigated. The results highlights the importance of the pedestrians' density and the effect of counterflow in the onset of avalanches and clogging and provides an estimation of the number of injuries based on a calculation of the contact-force network between the pedestrians. PMID:25615220
Electronic continuum model for molecular dynamics simulations.
Leontyev, I V; Stuchebrukhov, A A
2009-02-28
A simple model for accounting for electronic polarization in molecular dynamics (MD) simulations is discussed. In this model, called molecular dynamics electronic continuum (MDEC), the electronic polarization is treated explicitly in terms of the electronic continuum (EC) approximation, while the nuclear dynamics is described with a fixed-charge force field. In such a force-field all atomic charges are scaled to reflect the screening effect by the electronic continuum. The MDEC model is rather similar but not equivalent to the standard nonpolarizable force-fields; the differences are discussed. Of our particular interest is the calculation of the electrostatic part of solvation energy using standard nonpolarizable MD simulations. In a low-dielectric environment, such as protein, the standard MD approach produces qualitatively wrong results. The difficulty is in mistreatment of the electronic polarizability. We show how the results can be much improved using the MDEC approach. We also show how the dielectric constant of the medium obtained in a MD simulation with nonpolarizable force-field is related to the static (total) dielectric constant, which includes both the nuclear and electronic relaxation effects. Using the MDEC model, we discuss recent calculations of dielectric constants of alcohols and alkanes, and show that the MDEC results are comparable with those obtained with the polarizable Drude oscillator model. The applicability of the method to calculations of dielectric properties of proteins is discussed. PMID:19256627
Integrated computer simulation on FIR FEL dynamics
Furukawa, H.; Kuruma, S.; Imasaki, K.
1995-12-31
An integrated computer simulation code has been developed to analyze the RF-Linac FEL dynamics. First, the simulation code on the electron beam acceleration and transport processes in RF-Linac: (LUNA) has been developed to analyze the characteristics of the electron beam in RF-Linac and to optimize the parameters of RF-Linac. Second, a space-time dependent 3D FEL simulation code (Shipout) has been developed. The RF-Linac FEL total simulations have been performed by using the electron beam data from LUNA in Shipout. The number of particles using in a RF-Linac FEL total simulation is approximately 1000. The CPU time for the simulation of 1 round trip is about 1.5 minutes. At ILT/ILE, Osaka, a 8.5MeV RF-Linac with a photo-cathode RF-gun is used for FEL oscillation experiments. By using 2 cm wiggler, the FEL oscillation in the wavelength approximately 46 {mu}m are investigated. By the simulations using LUNA with the parameters of an ILT/ILE experiment, the pulse shape and the energy spectra of the electron beam at the end of the linac are estimated. The pulse shape of the electron beam at the end of the linac has sharp rise-up and it slowly decays as a function of time. By the RF-linac FEL total simulations with the parameters of an ILT/ILE experiment, the dependencies of the start up of the FEL oscillations on the pulse shape of the electron beam at the end of the linac are estimated. The coherent spontaneous emission effects and the quick start up of FEL oscillations have been observed by the RF-Linac FEL total simulations.
Molecular dynamics simulation of liquid sulfur dioxide.
Ribeiro, Mauro C C
2006-05-01
A previously proposed model for molecular dynamics (MD) simulation of liquid sulfur dioxide, SO(2), has been reviewed. Thermodynamic, structural, and dynamical properties were calculated for a large range of thermodynamic states. Predicted (P,V,T) of simulated system agrees with an elaborated equation of state recently proposed for liquid SO(2). Calculated heat capacity, expansion coefficient, and isothermal compressibility are also in good agreement with experimental data. Calculated equilibrium structure agrees with X-ray and neutron scattering measurements on liquid SO(2). The model also predicts the same (SO(2))(2) dimer structure as previously determined by ab initio calculations. Detailed analysis of equilibrium structure of liquid SO(2) is provided, indicating that, despite the rather large dipole moment of the SO(2) molecule, the structure is mainly determined by the Lennard-Jones interactions. Both single-particle and collective dynamics are investigated. Temperature dependency of dynamical properties is given. The MD results are compared with previous findings obtained from the analysis of inelastic neutron scattering spectra of liquid SO(2), including wave-vector dependent structural relaxation, tau(k), and viscosity, eta(k). PMID:16640437
Simulating stochastic dynamics using large time steps.
Corradini, O; Faccioli, P; Orland, H
2009-12-01
We present an approach to investigate the long-time stochastic dynamics of multidimensional classical systems, in contact with a heat bath. When the potential energy landscape is rugged, the kinetics displays a decoupling of short- and long-time scales and both molecular dynamics or Monte Carlo (MC) simulations are generally inefficient. Using a field theoretic approach, we perform analytically the average over the short-time stochastic fluctuations. This way, we obtain an effective theory, which generates the same long-time dynamics of the original theory, but has a lower time-resolution power. Such an approach is used to develop an improved version of the MC algorithm, which is particularly suitable to investigate the dynamics of rare conformational transitions. In the specific case of molecular systems at room temperature, we show that elementary integration time steps used to simulate the effective theory can be chosen a factor approximately 100 larger than those used in the original theory. Our results are illustrated and tested on a simple system, characterized by a rugged energy landscape. PMID:20365123
Simulating coronal condensation dynamics in 3D
NASA Astrophysics Data System (ADS)
Moschou, S. P.; Keppens, R.; Xia, C.; Fang, X.
2015-12-01
We present numerical simulations in 3D settings where coronal rain phenomena take place in a magnetic configuration of a quadrupolar arcade system. Our simulation is a magnetohydrodynamic simulation including anisotropic thermal conduction, optically thin radiative losses, and parametrised heating as main thermodynamical features to construct a realistic arcade configuration from chromospheric to coronal heights. The plasma evaporation from chromospheric and transition region heights eventually causes localised runaway condensation events and we witness the formation of plasma blobs due to thermal instability, that evolve dynamically in the heated arcade part and move gradually downwards due to interchange type dynamics. Unlike earlier 2.5D simulations, in this case there is no large scale prominence formation observed, but a continuous coronal rain develops which shows clear indications of Rayleigh-Taylor or interchange instability, that causes the denser plasma located above the transition region to fall down, as the system moves towards a more stable state. Linear stability analysis is used in the non-linear regime for gaining insight and giving a prediction of the system's evolution. After the plasma blobs descend through interchange, they follow the magnetic field topology more closely in the lower coronal regions, where they are guided by the magnetic dips.
Molecular Dynamics Simulation of Shock Induced Detonation
NASA Astrophysics Data System (ADS)
Tomar, Vikas; Zhou, Min
2004-07-01
This research focuses on molecular dynamics (MD) simulation of shock induced detonation in Fe2O3+Al thermite mixtures. A MD model is developed to simulate non-equilibrium stress-induced reactions. The focus is on establishing a criterion for reaction initiation, energy content and rate of energy release as functions of mixture and reinforcement characteristics. A cluster functional potential is proposed for this purpose. The potential uses the electronegativity equalization to account for changes in the charge of different species according to local environment. Parameters in the potential are derived to fit to the properties of Fe, Al, Fe2O3, and Al2O3. NPT MD simulations are carried out to qualitatively check the energetics of the forward (Fe2O3+Al) as well as backward (Al2O3+Fe) thermite reactions. The results show that the potential can account for the energetics of thermite reactions.
Dynamic simulator for PEFC propulsion plant
Hiraide, Masataka; Kaneda, Eiichi; Sato, Takao
1996-12-31
This report covers part of a joint study on a PEFC propulsion system for surface ships, summarized in a presentation to this Seminar, entitled {open_quote}Study on a Polymer Electrolyte Fuel Cell (PEFC) Propulsion System for Surface Ships{close_quotes}, and which envisages application to a 1,500 DWT cargo vessel. The work presented here focuses on a simulation study on PEFC propulsion plant performance, and particularly on the system response to changes in load. Using a dynamic simulator composed of system components including fuel cell, various simulations were executed, to examine the performance of the system as a whole and of the individual system components under quick and large load changes such as occasioned by maneuvering operations and by racing when the propeller emerges above water in heavy sea.
Molecular-dynamics simulations of lead clusters
NASA Astrophysics Data System (ADS)
Hendy, S. C.; Hall, B. D.
2001-08-01
Molecular-dynamics simulations of nanometer-sized lead clusters have been performed using the Lim-Ong-Ercolessi glue potential [Surf. Sci. 269/270, 1109 (1992)]. The binding energies of clusters forming crystalline (fcc), decahedron and icosahedron structures are compared, showing that fcc cuboctahedra are the most energetically favored of these polyhedral model structures. However, simulations of the freezing of liquid droplets produced a characteristic form of surface-reconstructed ``shaved'' icosahedron, in which atoms are absent at the edges and apexes of the polyhedron. This arrangement is energetically favored for 600-4000 atom clusters. Larger clusters favor crystalline structures. Indeed, simulated freezing of a 6525-atom liquid droplet produced an imperfect fcc Wulff particle, containing a number of parallel stacking faults. The effects of temperature on the preferred structure of crystalline clusters below the melting point have been considered. The implications of these results for the interpretation of experimental data is discussed.
Coupling a geodynamic seismic cycling model to rupture dynamic simulations
NASA Astrophysics Data System (ADS)
Gabriel, Alice; van Dinther, Ylona
2014-05-01
The relevance and results of dynamic rupture scenarios are implicitly linked to the geometry and pre-existing stress and strength state on a fault. The absolute stresses stored along faults during interseismic periods, are largely unquantifiable. They are, however, pivotal in defining coseismic rupture styles, near-field ground motion, and macroscopic source properties (Gabriel et al., 2012). Obtaining these in a physically consistent manner requires seismic cycling models, which directly couple long-term deformation processes (over 1000 year periods), the self-consistent development of faults, and the resulting dynamic ruptures. One promising approach to study seismic cycling enables both the generation of spontaneous fault geometries and the development of thermo-mechanically consistent fault stresses. This seismo-thermo-mechanical model has been developed using a methodology similar to that employed to study long-term lithospheric deformation (van Dinther et al., 2013a,b, using I2ELVIS of Gerya and Yuen, 2007). We will innovatively include the absolute stress and strength values along physically consistent evolving non-finite fault zones (regions of strain accumulation) from the geodynamic model into dynamic rupture simulations as an initial condition. The dynamic rupture simulations will be performed using SeisSol, an arbitrary high-order derivative Discontinuous Galerkin (ADER-DG) scheme (Pelties et al., 2012). The dynamic rupture models are able to incorporate the large degree of fault geometry complexity arising in naturally evolving geodynamic models. We focus on subduction zone settings with and without a splay fault. Due to the novelty of the coupling, we first focus on methodological challenges, e.g. the synchronization of both methods regarding the nucleation of events, the localization of fault planes, and the incorporation of similar frictional constitutive relations. We then study the importance of physically consistent fault stress, strength, and
NASA Astrophysics Data System (ADS)
Stringer, David Blake
The overarching objective in this research is the development of a robust, rotor dynamic, physics based model of a helicopter drive train as a foundation for the prognostic modeling for rotary-wing transmissions. Rotorcrafts rely on the integrity of their drive trains for their airworthiness. Drive trains rely on gear technology for their integrity and function. Gears alter the vibration characteristics of a mechanical system and significantly contribute to noise, component fatigue, and personal discomfort prevalent in rotorcraft. This research effort develops methodologies for generating a rotor dynamic model of a rotary-wing transmission based on first principles, through (i) development of a three-dimensional gear-mesh stiffness model for helical and spur gears and integration of this model in a finite element rotor dynamic model, (ii) linear and nonlinear analyses of a geared system for comparison and validation of the gear-mesh model, (iii) development of a modal synthesis technique for potentially providing model reduction and faster analysis capabilities for geared systems, and (iv) extension of the gear-mesh model to bevel and epicyclic configurations. In addition to model construction and validation, faults indigenous to geared systems are presented and discussed. Two faults are selected for analysis and seeded into the transmission model. Diagnostic vibration parameters are presented and used as damage indicators in the analysis. The fault models produce results consistent with damage experienced during experimental testing. The results of this research demonstrate the robustness of the physics-based approach in simulating multiple normal and abnormal conditions. The advantages of this physics-based approach, when combined with contemporary probabilistic and time-series techniques, provide a useful method for improving health monitoring technologies in mechanical systems.
Methodology of the Westinghouse dynamic rod worth measurement technique
Chao, Y.A.; Chapman, D.M.; Easter, M.E.; Hill, D.J.; Hoerner, J.A. ); Kurtz, P.N. )
1992-01-01
During zero-power physics testing, plant operations personnel use one of various techniques to measure the reactivity worth of the control rods to confirm shutdown margin. A simple and fast procedure for measuring rod worths called dynamic rod worth measurement (DRWM) has been developed at Westinghouse. This procedure was tested at the recent startups of Point Beach Nuclear Power Plant Unit 1 cycle 20 and Unit 2 cycle 18. The results of these tests show that DRWM measures rod worths with accuracy comparable to that of both boron dilution and rod bank exchange measurements. The DRWM procedure is a fast process of measuring the reactivity worth of individual banks by inserting and withdrawing the bank continuously at the maximum stepping speed without changing the boron concentration and recording the signals of the ex-core detectors.
INCORPORATING DYNAMIC 3D SIMULATION INTO PRA
Steven R Prescott; Curtis Smith
2011-07-01
provide superior results and insights. We also couple the state model with the dynamic 3D simulation analysis representing events (such as flooding) to determine which (if any) components fail. Not only does the simulation take into account any failed items from the state model, but any failures caused by the simulation are incorporated back into the state model and factored into the overall results. Using this method we incorporate accurate 3D simulation results, eliminate static-based PRA issues, and have time ordered failure information.
Statistical methodologies for the control of dynamic remapping
NASA Technical Reports Server (NTRS)
Saltz, J. H.; Nicol, D. M.
1986-01-01
Following an initial mapping of a problem onto a multiprocessor machine or computer network, system performance often deteriorates with time. In order to maintain high performance, it may be necessary to remap the problem. The decision to remap must take into account measurements of performance deterioration, the cost of remapping, and the estimated benefits achieved by remapping. We examine the tradeoff between the costs and the benefits of remapping two qualitatively different kinds of problems. One problem assumes that performance deteriorates gradually, the other assumes that performance deteriorates suddenly. We consider a variety of policies for governing when to remap. In order to evaluate these policies, statistical models of problem behaviors are developed. Simulation results are presented which compare simple policies with computationally expensive optimal decision policies; these results demonstrate that for each problem type, the proposed simple policies are effective and robust.
A new methodology to simulate subglacial deformation of water saturated granular material
NASA Astrophysics Data System (ADS)
Damsgaard, A.; Egholm, D. L.; Piotrowski, J. A.; Tulaczyk, S.; Larsen, N. K.; Brædstrup, C. F.
2015-07-01
The dynamics of glaciers are to a large degree governed by processes operating at the ice-bed interface, and one of the primary mechanisms of glacier flow over soft unconsolidated sediments is subglacial deformation. However, it has proven difficult to constrain the mechanical response of subglacial sediment to the shear stress of an overriding glacier. In this study, we present a new methodology designed to simulate subglacial deformation using a coupled numerical model for computational experiments on grain-fluid mixtures. The granular phase is simulated on a per-grain basis by the discrete element method. The pore water is modeled as a compressible Newtonian fluid without inertia. The numerical approach allows close monitoring of the internal behavior under a range of conditions. The rheology of a water-saturated granular bed may include both plastic and rate-dependent dilatant hardening or weakening components, depending on the rate of deformation, the material state, clay mineral content, and the hydrological properties of the material. The influence of the fluid phase is negligible when relatively permeable sediment is deformed. However, by reducing the local permeability, fast deformation can cause variations in the pore-fluid pressure. The pressure variations weaken or strengthen the granular phase, and in turn influence the distribution of shear strain with depth. In permeable sediments the strain distribution is governed by the grain-size distribution and effective normal stress and is typically on the order of tens of centimeters. Significant dilatant strengthening in impermeable sediments causes deformation to focus at the hydrologically more stable ice-bed interface, and results in a very shallow cm-to-mm deformational depth. The amount of strengthening felt by the glacier depends on the hydraulic conductivity at the ice-bed interface. Grain-fluid feedbacks can cause complex material properties that vary over time, and which may be of importance for
Monoamine transporters: insights from molecular dynamics simulations
Grouleff, Julie; Ladefoged, Lucy Kate; Koldsø, Heidi; Schiøtt, Birgit
2015-01-01
The human monoamine transporters (MATs) facilitate the reuptake of the neurotransmitters serotonin, dopamine, and norepinephrine from the synaptic cleft. Imbalance in monoaminergic neurotransmission is linked to various diseases including major depression, attention deficit hyperactivity disorder, schizophrenia, and Parkinson’s disease. Inhibition of the MATs is thus an important strategy for treatment of such diseases. The MATs are sodium-coupled transport proteins belonging to the neurotransmitter/Na+ symporter (NSS) family, and the publication of the first high-resolution structure of a NSS family member, the bacterial leucine transporter LeuT, in 2005, proved to be a major stepping stone for understanding this family of transporters. Structural data allows for the use of computational methods to study the MATs, which in turn has led to a number of important discoveries. The process of substrate translocation across the membrane is an intrinsically dynamic process. Molecular dynamics simulations, which can provide atomistic details of molecular motion on ns to ms timescales, are therefore well-suited for studying transport processes. In this review, we outline how molecular dynamics simulations have provided insight into the large scale motions associated with transport of the neurotransmitters, as well as the presence of external and internal gates, the coupling between ion and substrate transport, and differences in the conformational changes induced by substrates and inhibitors. PMID:26528185
Local Refinements in Classical Molecular Dynamics Simulations
NASA Astrophysics Data System (ADS)
Fackeldey, Konstantin; Weber, Marcus
2014-03-01
Quantum mechanics provide a detailed description of the physical and chemical behavior of molecules. However, with increasing size of the system the complexity rises exponentially, which is prohibitive for efficient dynamical simulation. In contrast, classical molecular dynamics procure a coarser description by using less degrees of freedom. Thus, it seems natural to seek for an adequate trade-off between accurateness and computational feasibility in the simulation of molecules. Here, we propose a novel method, which combines classical molecular simulations with quantum mechanics for molecular systems. For this we decompose the state space of the respective molecule into subsets, by employing a meshfree partition of unity. We show, that this partition allows us to localize an empirical force field and to run locally constrained classical trajectories. Within each subset, we compute the energy on the quantum level for a fixed number of spatial states (ab initio points). With these energy values from the ab initio points we have a local scattered data problem, which can be solved by the moving least squares method.
NASA Astrophysics Data System (ADS)
Prasetyo, Niko; Canaval, Lorenz R.; Wijaya, Karna; Armunanto, Ria
2015-01-01
The solvation of Li(I) in liquid ammonia has been investigated by an ab initio quantum mechanical charge-field molecular dynamics (QMCF-MD) simulation. Being the first simulation of a metal cation in liquid ammonia employing this methodology, the work yields a wide range of accurate structural and dynamical data. Li(I) is tetrahedrally coordinated by four ammonia molecules in the first solvation shell at a distance of 2.075 Å. Two ligand exchange attempts have been observed within 12 ps of simulation time. The second solvation shell shows a more labile structure with numerous successful exchanges. The results are in excellent agreement with experiments.
Dynamic transitions in molecular dynamics simulations of supercooled silicon
NASA Astrophysics Data System (ADS)
Mei, Xiaojun; Eapen, Jacob
2013-04-01
Two dynamic transitions or crossovers, one at a low temperature (T* ≈ 1006 K) and the other at a high temperature (T0 ≈ 1384 K), are shown to emerge in supercooled liquid silicon using molecular dynamics simulations. The high-temperature transition (T0) marks the decoupling of stress, density, and energy relaxation mechanisms. At the low-temperature transition (T*), depending on the cooling rate, supercooled silicon can either undergo a high-density-liquid to low-density-liquid (HDL-LDL) phase transition or experience an HDL-HDL crossover. Dynamically heterogeneous domains that emerge with supercooling become prominent across the HDL-HDL transition at 1006 K, with well-separated mobile and immobile regions. Interestingly, across the HDL-LDL transition, the most mobile atoms form large prominent aggregates while the least mobile atoms get spatially dispersed akin to that in a crystalline state. The attendant partial return to spatial uniformity with the HDL-LDL phase transition indicates a dynamic mechanism for relieving the frustration in supercooled states.
Dynamical simulation of dipolar Janus colloids: Dynamical properties
NASA Astrophysics Data System (ADS)
Hagy, Matthew C.; Hernandez, Rigoberto
2013-05-01
The dynamical properties of dipolar Janus particles are studied through simulation using our previously-developed detailed pointwise (PW) model and an isotropically coarse-grained (CG) model [M. C. Hagy and R. Hernandez, J. Chem. Phys. 137, 044505 (2012), 10.1063/1.4737432]. The CG model is found to have accelerated dynamics relative to the PW model over a range of conditions for which both models have near identical static equilibrium properties. Physically, this suggests dipolar Janus particles have slower transport properties (such as diffusion) in comparison to isotropically attractive particles. Time rescaling and damping with Langevin friction are explored to map the dynamics of the CG model to that of the PW model. Both methods map the diffusion constant successfully and improve the velocity autocorrelation function and the mean squared displacement of the CG model. Neither method improves the distribution of reversible bond durations f(tb) observed in the CG model, which is found to lack the longer duration reversible bonds observed in the PW model. We attribute these differences in f(tb) to changes in the energetics of multiple rearrangement mechanisms. This suggests a need for new methods that map the coarse-grained dynamics of such systems to the true time scale.
Molecular Dynamics Simulations of Graphene Oxide Frameworks
Zhu, Pan; Sumpter, Bobby G; Meunier, V.; Nicolai, Adrien
2013-01-01
We use quantum mechanical calculations to develop a full set of force field parameters in order to perform molecular dynamics simulations to understand and optimize the molecular storage properties inside Graphene Oxide Frameworks (GOFs). A set of boron-related parameters for commonly used empirical force fields is determined to describe the non-bonded and bonded interactions between linear boronic acid linkers and graphene sheets of GOF materials. The transferability of the parameters is discussed and their validity is quantified by comparing quantum mechanical and molecular mechanical structural and vibrational properties. The application of the model to the dynamics of water inside the GOFs reveals significant variations in structural flexibility of GOF depending on the linker density, which is shown to be usable as a tuning parameter for desired diffusion properties.
Finite element simulation of pipe dynamic response
Slagis, G.C.; Litton, R.W.
1996-12-01
Nonlinear finite element dynamic analyses of the response of a pipe span to controlled-displacement, sinusoidal vibration have been performed. The objective of this preliminary study is to compare strain and acceleration response data to those generated by Beaney in the Berkeley Nuclear Laboratories experiments. Results for an unpressurized, 5 Hz, carbon steel pipe are in good agreement with the experiments. Hence, it appears that analytical simulation will be useful to assess seismic margins. Recommendations for additional studies are provided. The analyses confirm the test results--dynamic response is greatly attenuated by material plasticity. Analytical strains and accelerations are about 30% higher than test data. There are several possible explanations for the differences. To assess the effect of frequency on response, the length of the pipe span was increased. Analysis of the longer, 2 Hz, pipe span shows significantly greater cyclic strains than the 5 Hz span at the same input excitation levels.
Molecular dynamics simulations of weak detonations.
Am-Shallem, Morag; Zeiri, Yehuda; Zybin, Sergey V; Kosloff, Ronnie
2011-12-01
Detonation of a three-dimensional reactive nonisotropic molecular crystal is modeled using molecular dynamics simulations. The detonation process is initiated by an impulse, followed by the creation of a stable fast reactive shock wave. The terminal shock velocity is independent of the initiation conditions. Further analysis shows supersonic propagation decoupled from the dynamics of the decomposed material left behind the shock front. The dependence of the shock velocity on crystal nonlinear compressibility resembles solitary behavior. These properties categorize the phenomena as a weak detonation. The dependence of the detonation wave on microscopic potential parameters was investigated. An increase in detonation velocity with the reaction exothermicity reaching a saturation value is observed. In all other respects the model crystal exhibits typical properties of a molecular crystal. PMID:22304055
Osmosis : a molecular dynamics computer simulation study
NASA Astrophysics Data System (ADS)
Lion, Thomas
Osmosis is a phenomenon of critical importance in a variety of processes ranging from the transport of ions across cell membranes and the regulation of blood salt levels by the kidneys to the desalination of water and the production of clean energy using potential osmotic power plants. However, despite its importance and over one hundred years of study, there is an ongoing confusion concerning the nature of the microscopic dynamics of the solvent particles in their transfer across the membrane. In this thesis the microscopic dynamical processes underlying osmotic pressure and concentration gradients are investigated using molecular dynamics (MD) simulations. I first present a new derivation for the local pressure that can be used for determining osmotic pressure gradients. Using this result, the steady-state osmotic pressure is studied in a minimal model for an osmotic system and the steady-state density gradients are explained using a simple mechanistic hopping model for the solvent particles. The simulation setup is then modified, allowing us to explore the timescales involved in the relaxation dynamics of the system in the period preceding the steady state. Further consideration is also given to the relative roles of diffusive and non-diffusive solvent transport in this period. Finally, in a novel modification to the classic osmosis experiment, the solute particles are driven out-of-equilibrium by the input of energy. The effect of this modification on the osmotic pressure and the osmotic ow is studied and we find that active solute particles can cause reverse osmosis to occur. The possibility of defining a new "osmotic effective temperature" is also considered and compared to the results of diffusive and kinetic temperatures..
Isotropic MD simulations of dynamic brittle fracture
Espanol, P.; Rubio, M.A.; Zuniga, I.
1996-12-01
The authors present results obtained by molecular dynamics simulations on the propagation of fast cracks in triangular 2D lattices. Their aim is to simulate Mode 1 fracture of brittle isotropic materials. They propose a force law that respects the isotropy of the material. The code yields the correct imposed sound c{sub {parallel}}, shear c{sub {perpendicular}} and surface V{sub R} wave speeds. Different notch lengths are systematically studied. They observed that initially the cracks are linear and always branch at a particular critical velocity c* {approx} 0.8V{sub R} and that this occurs when the crack tip reaches the position of a front emitted from the initial crack tip and propagating at a speed c = 0.68V{sub R}.
Dynamic Curving Simulation of Tilting Train
NASA Astrophysics Data System (ADS)
Zeng, Jing; Luo, Ren
The application of carbody tilting technology is the most efficient way to raise train speed during curve negotiations. This paper mainly deals with the dynamic performance simulation of the tilting train. Through the establishment of the nonlinear mathematical model for the titling train electromechanical coupled system, the carbody tilting control law, bogie radial steering mechanism, and titling train curving performance are investigated. The effect of time delay caused by the sensing and control system on the tilting performance of the train is analyzed, and the compensation methods for the time delay effect are studied.
Molecular dynamics simulations of dense plasmas
Collins, L.A.; Kress, J.D.; Kwon, I.; Lynch, D.L.; Troullier, N.
1993-12-31
We have performed quantum molecular dynamics simulations of hot, dense plasmas of hydrogen over a range of temperatures(0.1-5eV) and densities(0.0625-5g/cc). We determine the forces quantum mechanically from density functional, extended Huckel, and tight binding techniques and move the nuclei according to the classical equations of motion. We determine pair-correlation functions, diffusion coefficients, and electrical conductivities. We find that many-body effects predominate in this regime. We begin to obtain agreement with the OCP and Thomas-Fermi models only at the higher temperatures and densities.
Dynamic Deployment Simulations of Inflatable Space Structures
NASA Technical Reports Server (NTRS)
Wang, John T.
2005-01-01
The feasibility of using Control Volume (CV) method and the Arbitrary Lagrangian Eulerian (ALE) method in LSDYNA to simulate the dynamic deployment of inflatable space structures is investigated. The CV and ALE methods were used to predict the inflation deployments of three folded tube configurations. The CV method was found to be a simple and computationally efficient method that may be adequate for modeling slow inflation deployment sine the inertia of the inflation gas can be neglected. The ALE method was found to be very computationally intensive since it involves the solving of three conservative equations of fluid as well as dealing with complex fluid structure interactions.
Real-Time TD-DFT with Classical Ion Dynamics: Methodology and Applications.
Kolesov, Grigory; Grånäs, Oscar; Hoyt, Robert; Vinichenko, Dmitry; Kaxiras, Efthimios
2016-02-01
We present a method for real-time propagation of electronic wave functions, within time-dependent density functional theory (RT-TDDFT), coupled to ionic motion through mean-field classical dynamics. The goal of our method is to treat large systems and complex processes, in particular photocatalytic reactions and electron transfer events on surfaces and thin films. Due to the complexity of these processes, computational approaches are needed to provide insight into the underlying physical mechanisms and are therefore crucial for the rational design of new materials. Because of the short time step required for electron propagation (of order ∼10 attoseconds), these simulations are computationally very demanding. Our methodology is based on numerical atomic-orbital-basis sets for computational efficiency. In the computational package, to which we refer as TDAP-2.0 (Time-evolving Deterministic Atom Propagator), we have implemented a number of important features and analysis tools for more accurate and efficient treatment of large, complex systems and time scales that reach into a fraction of a picosecond. We showcase the capabilities of our method using four different examples: (i) photodissociation into radicals of opposite spin, (ii) hydrogen adsorption on aluminum surfaces, (iii) optical absorption of spin-polarized organic molecule containing a metal ion, and (iv) electron transfer in a prototypical dye-sensitized solar cell. PMID:26680129
Computational strategies in the dynamic simulation of constrained flexible MBS
NASA Technical Reports Server (NTRS)
Amirouche, F. M. L.; Xie, M.
1993-01-01
This research focuses on the computational dynamics of flexible constrained multibody systems. At first a recursive mapping formulation of the kinematical expressions in a minimum dimension as well as the matrix representation of the equations of motion are presented. The method employs Kane's equation, FEM, and concepts of continuum mechanics. The generalized active forces are extended to include the effects of high temperature conditions, such as creep, thermal stress, and elastic-plastic deformation. The time variant constraint relations for rolling/contact conditions between two flexible bodies are also studied. The constraints for validation of MBS simulation of gear meshing contact using a modified Timoshenko beam theory are also presented. The last part deals with minimization of vibration/deformation of the elastic beam in multibody systems making use of time variant boundary conditions. The above methodologies and computational procedures developed are being implemented in a program called DYAMUS.
A novel methodology for identifying emerging themes in small group dynamics.
Neuman, Yair; Assaf, Dan; Cohen, Yohai
2012-01-01
An understanding of group dynamics, such as the dynamics in family psychotherapy, is of great importance to mental health practitioners and other experts who seek to understand the group as a working whole. In this context, identifying themes emerging from the dynamics is of prime importance. However, the lack of formal tools makes it difficult to validly identify emerging themes. The current article presents a novel methodology for identifying themes emerging from group dynamics. With this methodology, the verbal utterances of group members are automatically analyzed to produce a "group matrix." Motifs that emerge in the complex network of signs that is generated by the group are analyzed by means of a computer program, and the explanatory value of the identified motifs is elaborated. The methodology and its benefits are presented and illustrated through the analysis of (1) the family dynamics in a literary piece, Tennessee Williams's The Glass Menagerie and (2) the group dynamics of Israeli and Palestinian participants discussing the notion of forgiveness. PMID:22409206
Methodology for testing infrared focal plane arrays in simulated nuclear radiation environments
NASA Astrophysics Data System (ADS)
Divita, E. L.; Mills, R. E.; Koch, T. L.; Gordon, M. J.; Wilcox, R. A.; Williams, R. E.
1992-07-01
This paper summarizes test methodology for focal plane array (FPA) testing that can be used for benign (clear) and radiation environments, and describes the use of custom dewars and integrated test equipment in an example environment. The test methodology, consistent with American Society for Testing Materials (ASTM) standards, is presented for the total accumulated gamma dose, transient dose rate, gamma flux, and neutron fluence environments. The merits and limitations of using Cobalt 60 for gamma environment simulations and of using various fast-neutron reactors and neutron sources for neutron simulations are presented. Test result examples are presented to demonstrate test data acquisition and FPA parameter performance under different measurement conditions and environmental simulations.
Visualizing Functional Motions of Membrane Transporters with Molecular Dynamics Simulations
2013-01-01
Computational modeling and molecular simulation techniques have become an integral part of modern molecular research. Various areas of molecular sciences continue to benefit from, indeed rely on, the unparalleled spatial and temporal resolutions offered by these technologies, to provide a more complete picture of the molecular problems at hand. Because of the continuous development of more efficient algorithms harvesting ever-expanding computational resources, and the emergence of more advanced and novel theories and methodologies, the scope of computational studies has expanded significantly over the past decade, now including much larger molecular systems and far more complex molecular phenomena. Among the various computer modeling techniques, the application of molecular dynamics (MD) simulation and related techniques has particularly drawn attention in biomolecular research, because of the ability of the method to describe the dynamical nature of the molecular systems and thereby to provide a more realistic representation, which is often needed for understanding fundamental molecular properties. The method has proven to be remarkably successful in capturing molecular events and structural transitions highly relevant to the function and/or physicochemical properties of biomolecular systems. Herein, after a brief introduction to the method of MD, we use a number of membrane transport proteins studied in our laboratory as examples to showcase the scope and applicability of the method and its power in characterizing molecular motions of various magnitudes and time scales that are involved in the function of this important class of membrane proteins. PMID:23298176
Dynamic stiffness removal for direct numerical simulations
Lu, Tianfeng; Law, Chung K.; Yoo, Chun Sang; Chen, Jacqueline H.
2009-08-15
A systematic approach was developed to derive non-stiff reduced mechanisms for direct numerical simulations (DNS) with explicit integration solvers. The stiffness reduction was achieved through on-the-fly elimination of short time-scales induced by two features of fast chemical reactivity, namely quasi-steady-state (QSS) species and partial-equilibrium (PE) reactions. The sparse algebraic equations resulting from QSS and PE approximations were utilized such that the efficiency of the dynamic stiffness reduction is high compared with general methods of time-scale reduction based on Jacobian decomposition. Using the dimension reduction strategies developed in our previous work, a reduced mechanism with 52 species was first derived from a detailed mechanism with 561 species. The reduced mechanism was validated for ignition and extinction applications over the parameter range of equivalence ratio between 0.5 and 1.5, pressure between 10 and 50 atm, and initial temperature between 700 and 1600 K for ignition, and worst-case errors of approximately 30% were observed. The reduced mechanism with dynamic stiffness removal was then applied in homogeneous and 1-D ignition applications, as well as a 2-D direct numerical simulation of ignition with temperature inhomogeneities at constant volume with integration time-steps of 5-10 ns. The integration was numerically stable and good accuracy was achieved. (author)
NASA Astrophysics Data System (ADS)
Jørgensen, Claus; Mark, Ole; Djordjevic, Slobodan; Hammond, Michael; Khan, David M.; Erichsen, Anders; Dorrit Enevoldsen, Ann; Heinicke, Gerald; Helwigh, Birgitte
2015-04-01
flood water, based on either measured waste water pathogen concentrations or on assumptions regarding the prevalence of infections in the population. The exposure (dosage) to pathogens was estimated by multiplying the concentration with literature values for the ingestion of water for different exposure groups (e.g. children, adults). The probability of infection was determined by applying dose response relations and MonteCarlo simulation. The methodology is demonstrated on two cases, i.e one case from a developing country with poor sanitation and one case from a developed country, where climate adaptation is the main issue: The risk of cholera in the City of Dhaka, Bangladesh during a flood event 2004, and the risk of bacterial and viral infections of during a flood event in Copenhagen, Denmark in 2011. Results PIC The historical flood events in Dhaka (2004) and Copenhagen (2011) were successfully modelled. The urban flood model was successfully coupled to QMRA. An example of the results of the quantitative microbial risk assessment given as the average estimated risk of cholera infection for children below 5 years living in slum areas in Dhaka is shown in the figure. Similarly, the risk of infection during the flood event in Copenhagen will be presented in the article. Conclusions We have developed a methodology for the dynamic modeling of the risk of infection during waste water influenced urban flooding. The outcome of the modelling exercise indicates that direct contact with polluted flood water is a likely route of transmission of cholera in Dhaka, and bacterial and viral infectious diseases in Copenhagen. It demonstrates the applicability and the potential for linking urban flood models with QMRA in order to identify interventions to reduce the burden of disease on the population in Dhaka City and Copenhagen.
Marshall, Deborah A; Burgos-Liz, Lina; IJzerman, Maarten J; Crown, William; Padula, William V; Wong, Peter K; Pasupathy, Kalyan S; Higashi, Mitchell K; Osgood, Nathaniel D
2015-03-01
In a previous report, the ISPOR Task Force on Dynamic Simulation Modeling Applications in Health Care Delivery Research Emerging Good Practices introduced the fundamentals of dynamic simulation modeling and identified the types of health care delivery problems for which dynamic simulation modeling can be used more effectively than other modeling methods. The hierarchical relationship between the health care delivery system, providers, patients, and other stakeholders exhibits a level of complexity that ought to be captured using dynamic simulation modeling methods. As a tool to help researchers decide whether dynamic simulation modeling is an appropriate method for modeling the effects of an intervention on a health care system, we presented the System, Interactions, Multilevel, Understanding, Loops, Agents, Time, Emergence (SIMULATE) checklist consisting of eight elements. This report builds on the previous work, systematically comparing each of the three most commonly used dynamic simulation modeling methods-system dynamics, discrete-event simulation, and agent-based modeling. We review criteria for selecting the most suitable method depending on 1) the purpose-type of problem and research questions being investigated, 2) the object-scope of the model, and 3) the method to model the object to achieve the purpose. Finally, we provide guidance for emerging good practices for dynamic simulation modeling in the health sector, covering all aspects, from the engagement of decision makers in the model design through model maintenance and upkeep. We conclude by providing some recommendations about the application of these methods to add value to informed decision making, with an emphasis on stakeholder engagement, starting with the problem definition. Finally, we identify areas in which further methodological development will likely occur given the growing "volume, velocity and variety" and availability of "big data" to provide empirical evidence and techniques
NASA Technical Reports Server (NTRS)
Murphy, Patrick C.; Klein, Vladislav
2003-01-01
A basic problem in flight dynamics is the mathematical formulation of the aerodynamic model for aircraft. This study is part of an ongoing effort at NASA Langley to develop a more general formulation of the aerodynamic model for aircraft that includes nonlinear unsteady aerodynamics and to develop appropriate test techniques that facilitate identification of these models. A methodology for modeling and testing using wide-band inputs to estimate the unsteady form of the aircraft aerodynamic model was developed previously but advanced test facilities were not available at that time to allow complete validation of the methodology. The new model formulation retained the conventional static and rotary dynamic terms but replaced conventional acceleration terms with more general indicial functions. In this study advanced testing techniques were utilized to validate the new methodology for modeling. Results of static, conventional forced oscillation, wide-band forced oscillation, oscillatory coning, and ramp tests are presented.
ERIC Educational Resources Information Center
Lenk, Hans
This document contains nine essays on the sociology and social psychology of team dynamics, including methodological and epistemological issues involved in such study. Essay titles are: (1) Conflict and Achievement in Top Athletic Teams--Sociometric Structures of Racing Eight Oar Crews; (2) Top Performance Despite Internal Conflict--An Antithesis…
ERIC Educational Resources Information Center
Hunleth, Jean
2011-01-01
By taking a reflexive approach to research methodology, this article contributes to discussions on power dynamics and knowledge production in the social studies of children. The author describes and analyzes three research methods that she used with children--drawing, child-led tape-recording and focus group discussions. These methods were carried…
Multiscale Simulation of Microbe Structure and Dynamics
Joshi, Harshad; Singharoy, Abhishek; Sereda, Yuriy V.; Cheluvaraja, Srinath C.; Ortoleva, Peter J.
2012-01-01
A multiscale mathematical and computational approach is developed that captures the hierarchical organization of a microbe. It is found that a natural perspective for understanding a microbe is in terms of a hierarchy of variables at various levels of resolution. This hierarchy starts with the N -atom description and terminates with order parameters characterizing a whole microbe. This conceptual framework is used to guide the analysis of the Liouville equation for the probability density of the positions and momenta of the N atoms constituting the microbe and its environment. Using multiscale mathematical techniques, we derive equations for the co-evolution of the order parameters and the probability density of the N-atom state. This approach yields a rigorous way to transfer information between variables on different space-time scales. It elucidates the interplay between equilibrium and far-from-equilibrium processes underlying microbial behavior. It also provides framework for using coarse-grained nanocharacterization data to guide microbial simulation. It enables a methodical search for free-energy minimizing structures, many of which are typically supported by the set of macromolecules and membranes constituting a given microbe. This suite of capabilities provides a natural framework for arriving at a fundamental understanding of microbial behavior, the analysis of nanocharacterization data, and the computer-aided design of nanostructures for biotechnical and medical purposes. Selected features of the methodology are demonstrated using our multiscale bionanosystem simulator DeductiveMultiscaleSimulator. Systems used to demonstrate the approach are structural transitions in the cowpea chlorotic mosaic virus, RNA of satellite tobacco mosaic virus, virus-like particles related to human papillomavirus, and iron-binding protein lactoferrin. PMID:21802438
Molecular dynamics simulation of amorphous indomethacin.
Xiang, Tian-Xiang; Anderson, Bradley D
2013-01-01
Molecular dynamics (MD) simulations have been conducted using an assembly consisting of 105 indomethacin (IMC) molecules and 12 water molecules to investigate the underlying dynamic (e.g., rotational and translational diffusivities and conformation relaxation rates) and structural properties (e.g., conformation, hydrogen-bonding distributions, and interactions of water with IMC) of amorphous IMC. These properties may be important in predicting physical stability of this metastable material. The IMC model was constructed using X-ray diffraction data with the force-field parameters mostly assigned by analogy with similar groups in Amber-ff03 and atomic charges calculated with the B3LYP/ccpVTZ30, IEFPCM, and RESP models. The assemblies were initially equilibrated in their molten state and cooled through the glass transition temperature to form amorphous solids. Constant temperature dynamic runs were then carried out above and below the T(g) (i.e., at 600 K (10 ns), 400 K (350 ns), and 298 K (240 ns)). The density (1.312 ± 0.003 g/cm(3)) of the simulated amorphous solid at 298 K was close to the experimental value (1.32 g/cm(3)) while the estimated T(g) (384 K) was ~64 degrees higher than the experimental value (320 K) due to the faster cooling rate. Due to the hindered rotation of its amide bond, IMC can exist in different diastereomeric states. Different IMC conformations were sufficiently sampled in the IMC melt or vapor, but transitions occurred rarely in the glass. The hydrogen-bonding patterns in amorphous IMC are more complex in the amorphous state than in the crystalline polymorphs. Carboxylic dimers that are dominant in α- and γ-crystals were found to occur at a much lower probability in the simulated IMC glasses while hydrogen-bonded IMC chains were more easily identified patterns in the simulated amorphous solids. To determine molecular diffusivity, a novel analytical method is proposed to deal with the non-Einsteinian behavior, in which the temporal
Photodynamics of oxybenzone sunscreen: Nonadiabatic dynamics simulations.
Li, Chun-Xiang; Guo, Wei-Wei; Xie, Bin-Bin; Cui, Ganglong
2016-08-21
Herein we have used combined static electronic structure calculations and "on-the-fly" global-switching trajectory surface-hopping dynamics simulations to explore the photochemical mechanism of oxybenzone sunscreen. We have first employed the multi-configurational CASSCF method to optimize minima, conical intersections, and minimum-energy reaction paths related to excited-state intramolecular proton transfer (ESIPT) and excited-state decays in the (1)ππ(∗), (1)nπ(∗), and S0 states (energies are refined at the higher MS-CASPT2 level). According to the mapped potential energy profiles, we have identified two ultrafast excited-state deactivation pathways for the initially populated (1)ππ(∗) system. The first is the diabatic ESIPT process along the (1)ππ(∗) potential energy profile. The generated (1)ππ(∗) keto species then decays to the S0 state via the keto (1)ππ(∗)/gs conical intersection. The second is internal conversion to the dark (1)nπ(∗) state near the (1)ππ(∗) /(1)nπ(∗) crossing point in the course of the diabatic (1)ππ(∗) ESIPT process. Our following dynamics simulations have shown that the ESIPT and (1)ππ(∗) → S0 internal conversion times are 104 and 286 fs, respectively. Finally, our present work demonstrates that in addition to the ESIPT process and the (1)ππ(∗) → S0 internal conversion in the keto region, the (1)ππ(∗) → (1)nπ(∗) internal conversion in the enol region plays as well an important role for the excited-state relaxation dynamics of oxybenzone. PMID:27544106
Atomistic molecular dynamic simulations of multiferroics.
Wang, Dawei; Weerasinghe, Jeevaka; Bellaiche, L
2012-08-10
A first-principles-based approach is developed to simulate dynamical properties, including complex permittivity and permeability in the GHz-THz range, of multiferroics at finite temperatures. It includes both structural degrees of freedom and magnetic moments as dynamic variables in Newtonian and Landau-Lifshitz-Gilbert (LLG) equations within molecular dynamics, respectively, with the couplings between these variables being incorporated. The use of a damping coefficient and of the fluctuation field in the LLG equations is required to obtain equilibrated magnetic properties at any temperature. No electromagnon is found in the spin-canted structure of BiFeO3. On the other hand, two magnons with very different frequencies are predicted via the use of this method. The smallest-in-frequency magnon corresponds to oscillations of the weak ferromagnetic vector in the basal plane being perpendicular to the polarization while the second magnon corresponds to magnetic dipoles going in and out of this basal plane. The large value of the frequency of this second magnon is caused by static couplings between magnetic dipoles with electric dipoles and oxygen octahedra tiltings. PMID:23006300
Atomistic Molecular Dynamic Simulations of Multiferroics
NASA Astrophysics Data System (ADS)
Wang, Dawei; Weerasinghe, Jeevaka; Bellaiche, L.
2012-08-01
A first-principles-based approach is developed to simulate dynamical properties, including complex permittivity and permeability in the GHz-THz range, of multiferroics at finite temperatures. It includes both structural degrees of freedom and magnetic moments as dynamic variables in Newtonian and Landau-Lifshitz-Gilbert (LLG) equations within molecular dynamics, respectively, with the couplings between these variables being incorporated. The use of a damping coefficient and of the fluctuation field in the LLG equations is required to obtain equilibrated magnetic properties at any temperature. No electromagnon is found in the spin-canted structure of BiFeO3. On the other hand, two magnons with very different frequencies are predicted via the use of this method. The smallest-in-frequency magnon corresponds to oscillations of the weak ferromagnetic vector in the basal plane being perpendicular to the polarization while the second magnon corresponds to magnetic dipoles going in and out of this basal plane. The large value of the frequency of this second magnon is caused by static couplings between magnetic dipoles with electric dipoles and oxygen octahedra tiltings.
Dynamical simulations of vesicle growth and division
NASA Astrophysics Data System (ADS)
Ruiz-Herrero, Teresa; Mahadevan, L.
2015-03-01
Prebiotic cells constitute a beautiful and intriguing example of self-replicating vesicles. How these cells managed to grow and divide without sophisticated machinery is still an open question. The properties of these primitive vesicles can shed light on the ways modern cells have evolved by exploiting those characteristics to develop their replication mechanisms. The equilibrium configurations of elastic shells are well understood, however the dynamical behavior during growth still lacks of a deep theoretical understanding. To study vesicle growth from a general perspective, we have developed a minimal generic model where vesicles are represented by a 2D spring network and characterized by a minimum set of magnitudes: growth rate, permeability, bending stiffness, viscosity and temperature. We have performed hybrid molecuar dynamic simulations as a function of a reduced set of dimensionless parameters. Three main outcomes were observed: vesicles that grow without division, vesicles that divide symmetrically, and vesicles that act as generators of daughter vesicles. The type of outcome depends on the system parameters and specifically on its dynamics via two timescales. Furthermore, we found sets of parameters where the system shows size homeostasis. TRH was supported by Ramon Areces Foundation.
Numerical Simulations of Droplet Dynamics in PEM Fuel Cell Microchannels
NASA Astrophysics Data System (ADS)
Cauble, Eric; Owkes, Mark
2015-11-01
Proton exchange membrane (PEM) fuel cells are of beneficial interest due to their capability of producing clean energy with zero emissions. An important design challenge hindering the performance of fuel cells is controlling water removal to maintain a hydrated membrane while avoiding excess water that may lead to channel blockage. Fuel cell water management requires a detailed knowledge of multiphase flow dynamics within microchannels. Direct observation of gas-liquid flows is difficult due to the small scale and viewing obstructions of the channels within the fuel cell. Instead, this work uses a CFD approach to compute the formation and dynamics of droplets in fuel cell channels. The method leverages a conservative volume-of-fluid (VOF) formulation coupled with a novel methodology to track dynamic contact angles. We present details of the numerical approach and simulation results relevant to water management in PEM fuel cells. In particular, it is shown that variation of the contact hysteresis angle influences the wetting properties of the droplet and significantly impacts water transport throughout the a fuel cell channel.
Quantum molecular dynamics simulations of dense matter
Collins, L.; Kress, J.; Troullier, N.; Lenosky, T.; Kwon, I.
1997-12-31
The authors have developed a quantum molecular dynamics (QMD) simulation method for investigating the properties of dense matter in a variety of environments. The technique treats a periodically-replicated reference cell containing N atoms in which the nuclei move according to the classical equations-of-motion. The interatomic forces are generated from the quantum mechanical interactions of the (between?) electrons and nuclei. To generate these forces, the authors employ several methods of varying sophistication from the tight-binding (TB) to elaborate density functional (DF) schemes. In the latter case, lengthy simulations on the order of 200 atoms are routinely performed, while for the TB, which requires no self-consistency, upwards to 1000 atoms are systematically treated. The QMD method has been applied to a variety cases: (1) fluid/plasma Hydrogen from liquid density to 20 times volume-compressed for temperatures of a thousand to a million degrees Kelvin; (2) isotopic hydrogenic mixtures, (3) liquid metals (Li, Na, K); (4) impurities such as Argon in dense hydrogen plasmas; and (5) metal/insulator transitions in rare gas systems (Ar,Kr) under high compressions. The advent of parallel versions of the methods, especially for fast eigensolvers, presage LDA simulations in the range of 500--1000 atoms and TB runs for tens of thousands of particles. This leap should allow treatment of shock chemistry as well as large-scale mixtures of species in highly transient environments.
Digital system for structural dynamics simulation
NASA Technical Reports Server (NTRS)
Krauter, A. I.; Lagace, L. J.; Wojnar, M. K.; Glor, C.
1982-01-01
State-of-the-art digital hardware and software for the simulation of complex structural dynamic interactions, such as those which occur in rotating structures (engine systems). System were incorporated in a designed to use an array of processors in which the computation for each physical subelement or functional subsystem would be assigned to a single specific processor in the simulator. These node processors are microprogrammed bit-slice microcomputers which function autonomously and can communicate with each other and a central control minicomputer over parallel digital lines. Inter-processor nearest neighbor communications busses pass the constants which represent physical constraints and boundary conditions. The node processors are connected to the six nearest neighbor node processors to simulate the actual physical interface of real substructures. Computer generated finite element mesh and force models can be developed with the aid of the central control minicomputer. The control computer also oversees the animation of a graphics display system, disk-based mass storage along with the individual processing elements.
Dynamics simulations for engineering macromolecular interactions
NASA Astrophysics Data System (ADS)
Robinson-Mosher, Avi; Shinar, Tamar; Silver, Pamela A.; Way, Jeffrey
2013-06-01
The predictable engineering of well-behaved transcriptional circuits is a central goal of synthetic biology. The artificial attachment of promoters to transcription factor genes usually results in noisy or chaotic behaviors, and such systems are unlikely to be useful in practical applications. Natural transcriptional regulation relies extensively on protein-protein interactions to insure tightly controlled behavior, but such tight control has been elusive in engineered systems. To help engineer protein-protein interactions, we have developed a molecular dynamics simulation framework that simplifies features of proteins moving by constrained Brownian motion, with the goal of performing long simulations. The behavior of a simulated protein system is determined by summation of forces that include a Brownian force, a drag force, excluded volume constraints, relative position constraints, and binding constraints that relate to experimentally determined on-rates and off-rates for chosen protein elements in a system. Proteins are abstracted as spheres. Binding surfaces are defined radially within a protein. Peptide linkers are abstracted as small protein-like spheres with rigid connections. To address whether our framework could generate useful predictions, we simulated the behavior of an engineered fusion protein consisting of two 20 000 Da proteins attached by flexible glycine/serine-type linkers. The two protein elements remained closely associated, as if constrained by a random walk in three dimensions of the peptide linker, as opposed to showing a distribution of distances expected if movement were dominated by Brownian motion of the protein domains only. We also simulated the behavior of fluorescent proteins tethered by a linker of varying length, compared the predicted Förster resonance energy transfer with previous experimental observations, and obtained a good correspondence. Finally, we simulated the binding behavior of a fusion of two ligands that could
Dynamics simulations for engineering macromolecular interactions
Robinson-Mosher, Avi; Shinar, Tamar; Silver, Pamela A.; Way, Jeffrey
2013-01-01
The predictable engineering of well-behaved transcriptional circuits is a central goal of synthetic biology. The artificial attachment of promoters to transcription factor genes usually results in noisy or chaotic behaviors, and such systems are unlikely to be useful in practical applications. Natural transcriptional regulation relies extensively on protein-protein interactions to insure tightly controlled behavior, but such tight control has been elusive in engineered systems. To help engineer protein-protein interactions, we have developed a molecular dynamics simulation framework that simplifies features of proteins moving by constrained Brownian motion, with the goal of performing long simulations. The behavior of a simulated protein system is determined by summation of forces that include a Brownian force, a drag force, excluded volume constraints, relative position constraints, and binding constraints that relate to experimentally determined on-rates and off-rates for chosen protein elements in a system. Proteins are abstracted as spheres. Binding surfaces are defined radially within a protein. Peptide linkers are abstracted as small protein-like spheres with rigid connections. To address whether our framework could generate useful predictions, we simulated the behavior of an engineered fusion protein consisting of two 20 000 Da proteins attached by flexible glycine/serine-type linkers. The two protein elements remained closely associated, as if constrained by a random walk in three dimensions of the peptide linker, as opposed to showing a distribution of distances expected if movement were dominated by Brownian motion of the protein domains only. We also simulated the behavior of fluorescent proteins tethered by a linker of varying length, compared the predicted Förster resonance energy transfer with previous experimental observations, and obtained a good correspondence. Finally, we simulated the binding behavior of a fusion of two ligands that could
A nonlinear dynamic finite element approach for simulating muscular hydrostats.
Vavourakis, V; Kazakidi, A; Tsakiris, D P; Ekaterinaris, J A
2014-01-01
An implicit nonlinear finite element model for simulating biological muscle mechanics is developed. The numerical method is suitable for dynamic simulations of three-dimensional, nonlinear, nearly incompressible, hyperelastic materials that undergo large deformations. These features characterise biological muscles, which consist of fibres and connective tissues. It can be assumed that the stress distribution inside the muscles is the superposition of stresses along the fibres and the connective tissues. The mechanical behaviour of the surrounding tissues is determined by adopting a Mooney-Rivlin constitutive model, while the mechanical description of fibres is considered to be the sum of active and passive stresses. Due to the nonlinear nature of the problem, evaluation of the Jacobian matrix is carried out in order to subsequently utilise the standard Newton-Raphson iterative procedure and to carry out time integration with an implicit scheme. The proposed methodology is implemented into our in-house, open source, finite element software, which is validated by comparing numerical results with experimental measurements and other numerical results. Finally, the numerical procedure is utilised to simulate primitive octopus arm manoeuvres, such as bending and reaching. PMID:23025686
NASA Astrophysics Data System (ADS)
Louri, Ahmed; Major, Michael C.
1995-07-01
Research in the field of free-space optical interconnection networks has reached a point where simula-tors and other design tools are desirable for reducing development costs and for improving design time. Previously proposed methodologies have only been applicable to simple systems. Our goal was to develop a simulation methodology capable of evaluating the performance characteristics for a variety of different free-space networks under a range of different configurations and operating states. The proposed methodology operates by first establishing the optical signal powers at various locations in the network. These powers are developed through the simulation by diffraction analysis of the light propagation through the network. After this evaluation, characteristics such as bit-error rate, signal-to-noise ratio, and system bandwidth are calculated. Further, the simultaneous evaluation of this process for a set of component misalignments provides a measure of the alignment tolerance of a design. We discuss this simulation process in detail as well as provide models for different optical interconnection network components.
A methodology for dynamic task allocation in a man-machine system
Parker, L.E.; Pin, F.G.
1987-01-01
This paper presents a methodological approach to the dynamic allocation of tasks in a man-machine symbiotic system in the context of dexterous manipulation and teleoperation. This paper addresses symbiosis containing two symbiotic partners which work toward controlling a single manipulator arm for the execution of a series of sequential manipulation tasks. The proposed automated task allocator uses knowledge about the constraints/criteria of the problem, the available resources, the tasks to be performed, and the environment to dynamically allocate tasks to the man and the machine. The presentation of the methodology includes discussions concerning the characteristics of the man-machine symbiotic system, the interaction of the knowledge areas, the flow of execution, and the dynamic nature of the task allocation.
Aldemir, Tunc; Denning, Richard; Catalyurek, Umit; Unwin, Stephen
2015-01-23
Reduction in safety margin can be expected as passive structures and components undergo degradation with time. Limitations in the traditional probabilistic risk assessment (PRA) methodology constrain its value as an effective tool to address the impact of aging effects on risk and for quantifying the impact of aging management strategies in maintaining safety margins. A methodology has been developed to address multiple aging mechanisms involving large numbers of components (with possibly statistically dependent failures) within the PRA framework in a computationally feasible manner when the sequencing of events is conditioned on the physical conditions predicted in a simulation environment, such as the New Generation System Code (NGSC) concept. Both epistemic and aleatory uncertainties can be accounted for within the same phenomenological framework and maintenance can be accounted for in a coherent fashion. The framework accommodates the prospective impacts of various intervention strategies such as testing, maintenance, and refurbishment. The methodology is illustrated with several examples.
Magnetotail dynamics: MHD simulations of driven and spontaneous dynamic changes
Birn, J.; Schindler, K.; Hesse, M.
1994-05-01
The dynamic evolution of the magnetotail during growth phase and expansion phase of a substorm is studied through threedimensional time-dependent MHD simulations. To model growth phase effects, an external electric field with an equatorward inflow is applied at the boundaries over a finite time period. This leads to the formation of a thin current sheet with greatly enhanced current density in the near tail, embedded in the wider plasma/current sheet, which becomes diminished in strength. A faster, spontaneous current sheet formation occurs when entropy conservation is released in an isobaric model, while the ideal MHD constraint persists. This may be a suitable model for the late, explosive part of the growth phase. The transition to the substorm expansive phase is modeled by an increase in anomalous resistivity, using either uniform resistivity or a current density dependent resistivity which is turned on when the current density exceeds a certain threshold. In both cases the violation of ideal MHD leads to resistive instability and the formation of a near-Earth neutral line, fast flow, and plasmoid ejection, together with the dipolarization and current reduction in the region further earthward. The spontaneous increase in total region 1 type field-aligned currents associated with the disruptions of the thin current sheets is less significant than that found in earlier simulations of the disruption of a wider current sheet, whereas the driven increase in the region 1 type current is substantial. The results demonstrate that the same dynamic process which appears spontaneous in the behavior of some quantities might be interpreted as entirely driven from the observation of others.
Assessing Molecular Dynamics Simulations with Solvatochromism Modeling.
Schwabe, Tobias
2015-08-20
For the modeling of solvatochromism with an explicit representation of the solvent molecules, the quality of preceding molecular dynamics simulations is crucial. Therefore, the possibility to apply force fields which are derived with as little empiricism as possible seems desirable. Such an approach is tested here by exploiting the sensitive solvatochromism of p-nitroaniline, and the use of reliable excitation energies based on approximate second-order coupled cluster results within a polarizable embedding scheme. The quality of the various MD settings for four different solvents, water, methanol, ethanol, and dichloromethane, is assessed. In general, good agreement with the experiment is observed when polarizable force fields and special treatment of hydrogen bonding are applied. PMID:26220273
CADS:Cantera Aerosol Dynamics Simulator.
Moffat, Harry K.
2007-07-01
This manual describes a library for aerosol kinetics and transport, called CADS (Cantera Aerosol Dynamics Simulator), which employs a section-based approach for describing the particle size distributions. CADS is based upon Cantera, a set of C++ libraries and applications that handles gas phase species transport and reactions. The method uses a discontinuous Galerkin formulation to represent the particle distributions within each section and to solve for changes to the aerosol particle distributions due to condensation, coagulation, and nucleation processes. CADS conserves particles, elements, and total enthalpy up to numerical round-off error, in all of its formulations. Both 0-D time dependent and 1-D steady state applications (an opposing-flow flame application) have been developed with CADS, with the initial emphasis on developing fundamental mechanisms for soot formation within fires. This report also describes the 0-D application, TDcads, which models a time-dependent perfectly stirred reactor.
Molecular Dynamics Simulations of Water Evaporation
NASA Astrophysics Data System (ADS)
Wen, Chengyuan; Grest, Gary; Cheng, Shengfeng
2015-03-01
The evaporation of water from the liquid/vapor interface is studied via large-scale molecular dynamics simulations for systems of more than a million atoms at 550K and 600K. The TIP4P-2005 water model whose liquid/vapor surface tension is in excellent agreement with experiments is used. Evaporative cooling at the interface is observed from temperature profiles determined from both translational and rotational kinetic energy. During evaporation, the density of water is slightly enhanced near the liquid-vapor interface. The velocity distribution of water molecules in the vapor phase during evaporation at various distances relative to the interface fit a Maxwell-Boltzmann distribution. While our results indicate an imbalance between evaporating and condensing water molecules, local thermal equilibrium is found to hold in addition to mechanical equilibrium. Department of Physics, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA.
Mathematical simulation of Earth system dynamics
NASA Astrophysics Data System (ADS)
Dymnikov, V. P.; Lykosov, V. N.; Volodin, E. M.
2015-05-01
The mathematical simulation of the Earth system, the dynamics of which depends on physical, chemical, biological, and other processes and which requires interdisciplinary approaches to studying this problem, is considered. The term "the Earth system" extends the concept "the climatic system," since additional geospheres (lithosphere, heliosphere, etc.) are taken into account and a wider range of physical, chemical, biological, and social interactions is described. The present-day level of climate modeling is discussed, and some data obtained at the Institute of Numerical Mathematics, Russian Academy of Sciences (INM RAS), are presented for this purpose. The prospects for further development of climate models toward the creation of the Earth system models based on a seamless approach, according to which a unified model is used to make short-term (several days) and long-term (climatic) prediction, are considered.
A new methodology to simulate subglacial deformation of water-saturated granular material
NASA Astrophysics Data System (ADS)
Damsgaard, A.; Egholm, D. L.; Piotrowski, J. A.; Tulaczyk, S.; Larsen, N. K.; Brædstrup, C. F.
2015-11-01
The dynamics of glaciers are to a large degree governed by processes operating at the ice-bed interface, and one of the primary mechanisms of glacier flow over soft unconsolidated sediments is subglacial deformation. However, it has proven difficult to constrain the mechanical response of subglacial sediment to the shear stress of an overriding glacier. In this study, we present a new methodology designed to simulate subglacial deformation using a coupled numerical model for computational experiments on grain-fluid mixtures. The granular phase is simulated on a per-grain basis by the discrete element method. The pore water is modeled as a compressible Newtonian fluid without inertia. The numerical approach allows close monitoring of the internal behavior under a range of conditions. Our computational experiments support the findings of previous studies where the rheology of a slowly deforming water-saturated granular bed in the steady state generally conforms to the rate-independent plastic rheology. Before this so-called critical state, deformation is in many cases accompanied by volumetric changes as grain rearrangement in active shear zones changes the local porosity. For previously consolidated beds porosity increases can cause local pore-pressure decline, dependent on till permeability and shear rate. We observe that the pore-water pressure reduction strengthens inter-granular contacts, which results in increased shear strength of the granular material. In contrast, weakening takes place when shear deformation causes consolidation of dilated sediments or during rapid fabric development. Both processes of strengthening and weakening depend inversely on the sediment permeability and are transient phenomena tied to the porosity changes during the early stages of shear. We find that the transient strengthening and weakening in turn influences the distribution of shear strain in the granular bed. Dilatant strengthening has the ability to distribute strain during
Reliability analysis of repairable systems using system dynamics modeling and simulation
NASA Astrophysics Data System (ADS)
Srinivasa Rao, M.; Naikan, V. N. A.
2014-07-01
Repairable standby system's study and analysis is an important topic in reliability. Analytical techniques become very complicated and unrealistic especially for modern complex systems. There have been attempts in the literature to evolve more realistic techniques using simulation approach for reliability analysis of systems. This paper proposes a hybrid approach called as Markov system dynamics (MSD) approach which combines the Markov approach with system dynamics simulation approach for reliability analysis and to study the dynamic behavior of systems. This approach will have the advantages of both Markov as well as system dynamics methodologies. The proposed framework is illustrated for a standby system with repair. The results of the simulation when compared with that obtained by traditional Markov analysis clearly validate the MSD approach as an alternative approach for reliability analysis.
Fiber lubrication: A molecular dynamics simulation study
NASA Astrophysics Data System (ADS)
Liu, Hongyi
Molecular and mesoscopic level description of friction and lubrication remains a challenge because of difficulties in the phenomenological understanding of to the behaviors of solid-liquid interfaces during sliding. Fortunately, there is the computational simulation approach opens an opportunity to predict and analyze interfacial phenomena, which were studied with molecular dynamics (MD) and mesoscopic dynamics (MesoDyn) simulations. Polypropylene (PP) and cellulose are two of most common polymers in textile fibers. Confined amorphous surface layers of PP and cellulose were built successfully with xenon crystals which were used to compact the polymers. The physical and surface properties of the PP and cellulose surface layers were investigated by MD simulations, including the density, cohesive energy, volumetric thermal expansion, and contact angle with water. The topology method was employed to predict the properties of poly(alkylene glycol) (PAG) diblock copolymers and Pluronic triblock copolymers used as lubricants on surfaces. Density, zero shear viscosity, shear module, cohesive energy and solubility parameter were predicted with each block copolymer. Molecular dynamics simulations were used to study the interaction energy per unit contact area of block copolymer melts with PP and cellulose surfaces. The interaction energy is defined as the ratio of interfacial interaction energy to the contact area. Both poly(proplene oxide) (PPO) and poly(ethylene oxide) (PEO) segments provided a lipophilic character to both PP and cellulose surfaces. The PPO/PEO ratio and the molecular weight were found to impact the interaction energy on both PP and cellulose surfaces. In aqueous solutions, the interaction energy is complicated due to the presence of water and the cross interactions between the multiple molecular components. The polymer-water-surface (PWS) calculation method was proposed to calculate such complex systems. In a contrast with a vacuum condition, the presence
In silico FRET from simulated dye dynamics
NASA Astrophysics Data System (ADS)
Hoefling, Martin; Grubmüller, Helmut
2013-03-01
Single molecule fluorescence resonance energy transfer (smFRET) experiments probe molecular distances on the nanometer scale. In such experiments, distances are recorded from FRET transfer efficiencies via the Förster formula, E=1/(1+(). The energy transfer however also depends on the mutual orientation of the two dyes used as distance reporter. Since this information is typically inaccessible in FRET experiments, one has to rely on approximations, which reduce the accuracy of these distance measurements. A common approximation is an isotropic and uncorrelated dye orientation distribution. To assess the impact of such approximations, we present the algorithms and implementation of a computational toolkit for the simulation of smFRET on the basis of molecular dynamics (MD) trajectory ensembles. In this study, the dye orientation dynamics, which are used to determine dynamic FRET efficiencies, are extracted from MD simulations. In a subsequent step, photons and bursts are generated using a Monte Carlo algorithm. The application of the developed toolkit on a poly-proline system demonstrated good agreement between smFRET simulations and experimental results and therefore confirms our computational method. Furthermore, it enabled the identification of the structural basis of measured heterogeneity. The presented computational toolkit is written in Python, available as open-source, applicable to arbitrary systems and can easily be extended and adapted to further problems. Catalogue identifier: AENV_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AENV_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GPLv3, the bundled SIMD friendly Mersenne twister implementation [1] is provided under the SFMT-License. No. of lines in distributed program, including test data, etc.: 317880 No. of bytes in distributed program, including test data, etc.: 54774217 Distribution format: tar.gz Programming language
Building Energy Simulation Test for Existing Homes (BESTEST-EX) Methodology: Preprint
Judkoff, R.; Polly, B.; Bianchi, M.; Neymark, J.
2011-11-01
The test suite represents a set of cases applying the new Building Energy Simulation Test for Existing Homes (BESTEST-EX) Methodology developed by NREL. (Judkoff et al. 2010a). The NREL team developed the test cases in consultation with the home retrofit industry (BESTEST-EX Working Group 2009), and adjusted the test specifications in accordance with information supplied by a participant with access to large utility bill datasets (Blasnik 2009).
Numerical simulations of Modified Newtonian Dynamics
NASA Astrophysics Data System (ADS)
Candlish, G. N.; Smith, R.; Fellhauer, M.
2016-05-01
The ΛCDM standard cosmological model is strongly supported by multiple lines of evidence, particularly from observations at large scales such as the CMB and large scale structure. There are some indications, however, of problems at smaller scales. An alternative to the CDM approach is to modify the gravitational force, as exemplified by the MOdified Newtonian Dynamics (MOND) idea. While evidence suggests MOND cannot account for dynamics at all scales without dark matter, it has been successful at galactic scales. Due to the complexity of the theory, however, most tests of MOND have extended no further than using a simple scaling relation to determine rotation curves or velocity dispersions. Therefore, to test the concept more thoroughly we require numerical simulations. We discuss the development and testing of a new N-body solver, using two distinct formulations of MOND, that is incorporated into the RAMSES code. The theory of MOND as a modification of Newtonian gravity is briefly summarised. We then show how it is implemented in the code, providing an example of an idealised test case and future applications.
Numerical simulation of tulip flame dynamics
Cloutman, L.D.
1991-11-30
A finite difference reactive flow hydrodynamics program based on the full Navier-Stokes equations was used to simulate the combustion process in a homogeneous-charge, constant-volume combustion bomb in which an oddly shaped flame, known as a ``tulip flame`` in the literature, occurred. The ``tulip flame`` was readily reproduced in the numerical simulations, producing good agreement with the experimental flame shapes and positions at various times. The calculations provide sufficient detail about the dynamics of the experiment to provide some insight into the physical mechanisms responsible for the peculiar flame shape. Several factors seem to contribute to the tulip formation. The most important process is the baroclinic production of vorticity by the flame front, and this rate of production appears to be dramatically increased by the nonaxial flow generated when the initial semicircular flame front burns out along the sides of the chamber. The vorticity produces a pair of vortices behind the flame that advects the flame into the tulip shape. Boundary layer effects contribute to the details of the flame shape next to the walls of the chamber, but are otherwise not important. 24 refs.
Numerical simulation of tulip flame dynamics
Cloutman, L.D.
1991-11-30
A finite difference reactive flow hydrodynamics program based on the full Navier-Stokes equations was used to simulate the combustion process in a homogeneous-charge, constant-volume combustion bomb in which an oddly shaped flame, known as a tulip flame'' in the literature, occurred. The tulip flame'' was readily reproduced in the numerical simulations, producing good agreement with the experimental flame shapes and positions at various times. The calculations provide sufficient detail about the dynamics of the experiment to provide some insight into the physical mechanisms responsible for the peculiar flame shape. Several factors seem to contribute to the tulip formation. The most important process is the baroclinic production of vorticity by the flame front, and this rate of production appears to be dramatically increased by the nonaxial flow generated when the initial semicircular flame front burns out along the sides of the chamber. The vorticity produces a pair of vortices behind the flame that advects the flame into the tulip shape. Boundary layer effects contribute to the details of the flame shape next to the walls of the chamber, but are otherwise not important. 24 refs.
Molecular dynamics simulations of microscale fluid transport
Wong, C.C.; Lopez, A.R.; Stevens, M.J.; Plimpton, S.J.
1998-02-01
Recent advances in micro-science and technology, like Micro-Electro-Mechanical Systems (MEMS), have generated a group of unique liquid flow problems that involve characteristic length scales of a Micron. Also, in manufacturing processes such as coatings, current continuum models are unable to predict microscale physical phenomena that appear in these non-equilibrium systems. It is suspected that in these systems, molecular-level processes can control the interfacial energy and viscoelastic properties at the liquid/solid boundary. A massively parallel molecular dynamics (MD) code has been developed to better understand microscale transport mechanisms, fluid-structure interactions, and scale effects in micro-domains. Specifically, this MD code has been used to analyze liquid channel flow problems for a variety of channel widths, e.g. 0.005-0.05 microns. This report presents results from MD simulations of Poiseuille flow and Couette flow problems and addresses both scaling and modeling issues. For Poiseuille flow, the numerical predictions are compared with existing data to investigate the variation of the friction factor with channel width. For Couette flow, the numerical predictions are used to determine the degree of slip at the liquid/solid boundary. Finally, the results also indicate that shear direction with respect to the wall lattice orientation can be very important. Simulation results of microscale Couette flow and microscale Poiseuille flow for two different surface structures and two different shear directions will be presented.
Lewis, F.M.; Voss, C.I.; Rubin, J.
1987-01-01
Methodologies that account for specific types of chemical reactions in the simulation of solute transport can be developed so they are compatible with solution algorithms employed in existing transport codes. This enables the simulation of reactive transport in complex multidimensional flow regimes, and provides a means for existing codes to account for some of the fundamental chemical processes that occur among transported solutes. Two equilibrium-controlled reaction systems demonstrate a methodology for accommodating chemical interaction into models of solute transport. One system involves the sorption of a given chemical species, as well as two aqueous complexations in which the sorbing species is a participant. The other reaction set involves binary ion exchange coupled with aqueous complexation involving one of the exchanging species. The methodology accommodates these reaction systems through the addition of nonlinear terms to the transport equations for the sorbing species. Example simulation results show (1) the effect equilibrium chemical parameters have on the spatial distributions of concentration for complexing solutes; (2) that an interrelationship exists between mechanical dispersion and the various reaction processes; (3) that dispersive parameters of the porous media cannot be determined from reactive concentration distributions unless the reaction is accounted for or the influence of the reaction is negligible; (4) how the concentration of a chemical species may be significantly affected by its participation in an aqueous complex with a second species which also sorbs; and (5) that these coupled chemical processes influencing reactive transport can be demonstrated in two-dimensional flow regimes. ?? 1987.
Consequence modeling using the fire dynamics simulator.
Ryder, Noah L; Sutula, Jason A; Schemel, Christopher F; Hamer, Andrew J; Van Brunt, Vincent
2004-11-11
The use of Computational Fluid Dynamics (CFD) and in particular Large Eddy Simulation (LES) codes to model fires provides an efficient tool for the prediction of large-scale effects that include plume characteristics, combustion product dispersion, and heat effects to adjacent objects. This paper illustrates the strengths of the Fire Dynamics Simulator (FDS), an LES code developed by the National Institute of Standards and Technology (NIST), through several small and large-scale validation runs and process safety applications. The paper presents two fire experiments--a small room fire and a large (15 m diameter) pool fire. The model results are compared to experimental data and demonstrate good agreement between the models and data. The validation work is then extended to demonstrate applicability to process safety concerns by detailing a model of a tank farm fire and a model of the ignition of a gaseous fuel in a confined space. In this simulation, a room was filled with propane, given time to disperse, and was then ignited. The model yields accurate results of the dispersion of the gas throughout the space. This information can be used to determine flammability and explosive limits in a space and can be used in subsequent models to determine the pressure and temperature waves that would result from an explosion. The model dispersion results were compared to an experiment performed by Factory Mutual. Using the above examples, this paper will demonstrate that FDS is ideally suited to build realistic models of process geometries in which large scale explosion and fire failure risks can be evaluated with several distinct advantages over more traditional CFD codes. Namely transient solutions to fire and explosion growth can be produced with less sophisticated hardware (lower cost) than needed for traditional CFD codes (PC type computer verses UNIX workstation) and can be solved for longer time histories (on the order of hundreds of seconds of computed time) with
NASA Technical Reports Server (NTRS)
Lowrie, J. W.; Fermelia, A. J.; Haley, D. C.; Gremban, K. D.; Vanbaalen, J.; Walsh, R. W.
1982-01-01
The derivation of the equations is presented, the rate control algorithm described, and simulation methodologies summarized. A set of dynamics equations that can be used recursively to calculate forces and torques acting at the joints of an n link manipulator given the manipulator joint rates are derived. The equations are valid for any n link manipulator system with any kind of joints connected in any sequence. The equations of motion for the class of manipulators consisting of n rigid links interconnected by rotary joints are derived. A technique is outlined for reducing the system of equations to eliminate contraint torques. The linearized dynamics equations for an n link manipulator system are derived. The general n link linearized equations are then applied to a two link configuration. The coordinated rate control algorithm used to compute individual joint rates when given end effector rates is described. A short discussion of simulation methodologies is presented.
Dynamic Shade and Irradiance Simulation of Aquatic Landscapes and Watersheds
Penumbra is a landscape shade and irradiance simulation model that simulates how solar energy spatially and temporally interacts within dynamic ecosystems such as riparian zones, forests, and other terrain that cast topological shadows. Direct and indirect solar energy accumulate...
Validation of vehicle dynamics simulation models - a review
NASA Astrophysics Data System (ADS)
Kutluay, Emir; Winner, Hermann
2014-02-01
In this work, a literature survey on the validation of vehicle dynamics simulation models is presented. Estimating the dynamic responses of existing or proposed vehicles has a wide array of applications in the development of vehicle technologies, e.g. active suspensions, controller design, driver assistance systems, etc. Although simulation environments, measurement tools and mathematical theories on vehicle dynamics are well established, the methodical link between the experimental test data and validity analysis of the simulation model is still lacking. This report presents different views on the definition of validation, and its usage in vehicle dynamics simulation models.
Studying Interactions by Molecular Dynamics Simulations at High Concentration
Fogolari, Federico; Corazza, Alessandra; Toppo, Stefano; Tosatto, Silvio C. E.; Viglino, Paolo; Ursini, Fulvio; Esposito, Gennaro
2012-01-01
Molecular dynamics simulations have been used to study molecular encounters and recognition. In recent works, simulations using high concentration of interacting molecules have been performed. In this paper, we consider the practical problems for setting up the simulation and to analyse the results of the simulation. The simulation of beta 2-microglobulin association and the simulation of the binding of hydrogen peroxide by glutathione peroxidase are provided as examples. PMID:22500085
Analysis & Simulation of Dynamics in Supercooled Liquids
NASA Astrophysics Data System (ADS)
Elmatad, Yael Sarah
2011-12-01
The nature of supercooled liquids and the glass transition has been debated by many scientists. Several theories have been put forth to describe the remarkable properties of this out-of-equilibrium material. Each of these theories makes specific predictions as to how the scaling of various transport properties in supercooled materials should behave. Given access to a large pool of high-quality supercooled liquid data we seek to compare these theories to one another. Moreover, we explore properties of a pair of models which are the basis for one particularly attractive theory---Chandler-Garrahan theory---and discuss the models' behavior in space-time and possible implications to the behavior of experimental supercooled liquids. Here we investigate the nature of dynamics in supercooled liquids using a two pronged approach. First we analyze the transport properties found in experiments and simulations of supercooled liquids. Then, we analyze simulation trajectories for lattice models which reproduce many of the interesting properties of supercooled liquids. In doing so, we illuminate several glass universalities, common properties of a wide variety of glass formers. By analyzing relaxation time and viscosity data for over 50 data sets and 1200 points, we find that relaxation time can be collapsed onto a single, parabolic curve. This collapse supports a theory of universal glass behavior based on facilitated models proposed by David Chandler and Juan Garrahan in 2003. We then show that the parabolic fit parameters for any particular liquid are a material property: they converge fast and are capable of predicting behavior in regions beyond the included data sets. We compare this property to other popular fitting schemes such as the Vogel-Fulcher, double exponential, and fractional exponential forms and conclude that these three forms result in parameters which are non predictive and therefore not material properties. Additionally, we examine the role of attractive
Interest of the Theory of Uncertain in the Dynamic LCA- Fire Methodology to Assess Fire Effects
NASA Astrophysics Data System (ADS)
Chettouh, Samia; Hamzi, Rachida; Innal, Fares; Haddad, Djamel
Life Cycle Impact Assessment (LCIA) is the third phase of Life Cycle Assessment (LCA) described in ISO 14042. The purpose of LCIA is to assess a product system's life cycle inventory analysis (LCI) in order to better understand its environmental significance. However, LCIA typically excludes spatial, temporal, threshold and dose-response information, and combines emissions or activities over space and/or time. This may diminish the environmental relevance of the indicator result. The methodology, Dynamic LCA -Fire proposed in this paper to complete the International Standard ISO 14042 in the fire field, combines the LCA - Fire method with the Dispersion Numerical Model. It is based on the use of the plume model used to assess pollutant concentrations and thermal effects from fire accident scenarios. In this study, The Dynamic LCA - Fire methodology is applied to a case study for petroleum production process management.
Probabilistic dynamics of mistuned bladed disc systems using subset simulation
NASA Astrophysics Data System (ADS)
Yuan, Jie; Allegri, Giuliano; Scarpa, Fabrizio; Rajasekaran, Ramesh; Patsias, Sophoclis
2015-08-01
The work describes an assessment of subset simulation (SubSim) techniques to increase the computational efficiency for the predictions of probabilistic dynamic behaviour in mistuned bladed disc systems. SubSim is an adaptive stochastic procedure to efficiently compute small failure probabilities, which are expressed as a product of large conditional failures probabilities by introducing intermediate failure events. The original version of SubSim with a classical modified Markov chain Monte Carlo (MCMC) method is used in this work to generate samples related to intermediate failure events. A 2-DOFs model with lumped parameters identified from a high-fidelity finite element model is used to represent a bladed disc. The statistics associated to the maximum forced frequency response amplitudes are evaluated from different levels of the blade mistuning using stiffness perturbations of the blades. Direct Monte Carlo simulations (MCS) are used to benchmark the results from the SubSim. The proposed methodology is shown to capture efficiently the statistical properties of the mistuned blades with less than 5% samples compared to the direct MCS method. Trade-off parametric studies of the SubSim method indicate that 2000 samples at each level yield an overall good computational efficiency and accuracy for the bladed disk system considered in this work. The study confirms that SubSim techniques can be effectively used in stochastic analysis of bladed disc systems with uncertainty related to the blade configurations.
Rotational Brownian Dynamics simulations of clathrin cage formation
Ilie, Ioana M.; Briels, Wim J.; Otter, Wouter K. den
2014-08-14
The self-assembly of nearly rigid proteins into ordered aggregates is well suited for modeling by the patchy particle approach. Patchy particles are traditionally simulated using Monte Carlo methods, to study the phase diagram, while Brownian Dynamics simulations would reveal insights into the assembly dynamics. However, Brownian Dynamics of rotating anisotropic particles gives rise to a number of complications not encountered in translational Brownian Dynamics. We thoroughly test the Rotational Brownian Dynamics scheme proposed by Naess and Elsgaeter [Macromol. Theory Simul. 13, 419 (2004); Naess and Elsgaeter Macromol. Theory Simul. 14, 300 (2005)], confirming its validity. We then apply the algorithm to simulate a patchy particle model of clathrin, a three-legged protein involved in vesicle production from lipid membranes during endocytosis. Using this algorithm we recover time scales for cage assembly comparable to those from experiments. We also briefly discuss the undulatory dynamics of the polyhedral cage.
Mathematical model of marine diesel engine simulator for a new methodology of self propulsion tests
NASA Astrophysics Data System (ADS)
Izzuddin, Nur; Sunarsih, Priyanto, Agoes
2015-05-01
As a vessel operates in the open seas, a marine diesel engine simulator whose engine rotation is controlled to transmit through propeller shaft is a new methodology for the self propulsion tests to track the fuel saving in a real time. Considering the circumstance, this paper presents the real time of marine diesel engine simulator system to track the real performance of a ship through a computer-simulated model. A mathematical model of marine diesel engine and the propeller are used in the simulation to estimate fuel rate, engine rotating speed, thrust and torque of the propeller thus achieve the target vessel's speed. The input and output are a real time control system of fuel saving rate and propeller rotating speed representing the marine diesel engine characteristics. The self-propulsion tests in calm waters were conducted using a vessel model to validate the marine diesel engine simulator. The simulator then was used to evaluate the fuel saving by employing a new mathematical model of turbochargers for the marine diesel engine simulator. The control system developed will be beneficial for users as to analyze different condition of vessel's speed to obtain better characteristics and hence optimize the fuel saving rate.
Nanosculpt: A methodology for generating complex realistic configurations for atomistic simulations
Prakash, A.; Hummel, M.; Schmauder, S.; Bitzek, E.
2016-01-01
Atomistic simulations have now become commonplace in the study of the deformation and failure of materials. Increase in computing power in recent years has made large-scale simulations with billions, or even trillions, of atoms a possibility. Most simulations to-date, however, are still performed with quasi-2D geometries or rather simplistic 3D setups. Although controlled studies on such well-defined structures are often required to obtain quantitative information from atomistic simulations, for qualitative studies focusing on e.g. the identification of mechanisms, researchers would greatly benefit from a methodology that helps realize more realistic configurations. The ideal scenario would be a one-on-one reconstruction of experimentally observed structures. To this end, we propose a new method and software tool called nanosculpt with the following features:•The method allows for easy sample generation for atomistic simulations from any arbitrarily shaped 3D enclosed volume.•The tool can be used to build atomistic samples from artificial geometries, including CAD geometries and structures obtained from simulation methods other than atomistic simulations.•The tool enables the generation of experimentally informed atomistic samples, by e.g. digitization of micrographs or usage of tomography data. PMID:27054098
Mathematical model of marine diesel engine simulator for a new methodology of self propulsion tests
Izzuddin, Nur; Sunarsih,; Priyanto, Agoes
2015-05-15
As a vessel operates in the open seas, a marine diesel engine simulator whose engine rotation is controlled to transmit through propeller shaft is a new methodology for the self propulsion tests to track the fuel saving in a real time. Considering the circumstance, this paper presents the real time of marine diesel engine simulator system to track the real performance of a ship through a computer-simulated model. A mathematical model of marine diesel engine and the propeller are used in the simulation to estimate fuel rate, engine rotating speed, thrust and torque of the propeller thus achieve the target vessel’s speed. The input and output are a real time control system of fuel saving rate and propeller rotating speed representing the marine diesel engine characteristics. The self-propulsion tests in calm waters were conducted using a vessel model to validate the marine diesel engine simulator. The simulator then was used to evaluate the fuel saving by employing a new mathematical model of turbochargers for the marine diesel engine simulator. The control system developed will be beneficial for users as to analyze different condition of vessel’s speed to obtain better characteristics and hence optimize the fuel saving rate.
Prakash, A; Hummel, M; Schmauder, S; Bitzek, E
2016-01-01
Atomistic simulations have now become commonplace in the study of the deformation and failure of materials. Increase in computing power in recent years has made large-scale simulations with billions, or even trillions, of atoms a possibility. Most simulations to-date, however, are still performed with quasi-2D geometries or rather simplistic 3D setups. Although controlled studies on such well-defined structures are often required to obtain quantitative information from atomistic simulations, for qualitative studies focusing on e.g. the identification of mechanisms, researchers would greatly benefit from a methodology that helps realize more realistic configurations. The ideal scenario would be a one-on-one reconstruction of experimentally observed structures. To this end, we propose a new method and software tool called nano sculpt with the following features:•The method allows for easy sample generation for atomistic simulations from any arbitrarily shaped 3D enclosed volume.•The tool can be used to build atomistic samples from artificial geometries, including CAD geometries and structures obtained from simulation methods other than atomistic simulations.•The tool enables the generation of experimentally informed atomistic samples, by e.g. digitization of micrographs or usage of tomography data. PMID:27054098
High frequency dynamic engine simulation. [TF-30 engine
NASA Technical Reports Server (NTRS)
Schuerman, J. A.; Fischer, K. E.; Mclaughlin, P. W.
1977-01-01
A digital computer simulation of a mixed flow, twin spool turbofan engine was assembled to evaluate and improve the dynamic characteristics of the engine simulation to disturbance frequencies of at least 100 Hz. One dimensional forms of the dynamic mass, momentum and energy equations were used to model the engine. A TF30 engine was simulated so that dynamic characteristics could be evaluated against results obtained from testing of the TF30 engine at the NASA Lewis Research Center. Dynamic characteristics of the engine simulation were improved by modifying the compression system model. Modifications to the compression system model were established by investigating the influence of size and number of finite dynamic elements. Based on the results of this program, high frequency engine simulations using finite dynamic elements can be assembled so that the engine dynamic configuration is optimum with respect to dynamic characteristics and computer execution time. Resizing of the compression systems finite elements improved the dynamic characteristics of the engine simulation but showed that additional refinements are required to obtain close agreement simulation and actual engine dynamic characteristics.
Nanoscale deicing by molecular dynamics simulation
NASA Astrophysics Data System (ADS)
Xiao, Senbo; He, Jianying; Zhang, Zhiliang
2016-07-01
Deicing is important to human activities in low-temperature circumstances, and is critical for combating the damage caused by excessive accumulation of ice. The aim of creating anti-icing materials, surfaces and applications relies on the understanding of fundamental nanoscale ice adhesion mechanics. Here in this study, we employ all-atom modeling and molecular dynamics simulation to investigate ice adhesion. We apply force to detach and shear nano-sized ice cubes for probing the determinants of atomistic adhesion mechanics, and at the same time investigate the mechanical effect of a sandwiched aqueous water layer between ice and substrates. We observe that high interfacial energy restricts ice mobility and increases both ice detaching and shearing stresses. We quantify up to a 60% decrease in ice adhesion strength by an aqueous water layer, and provide atomistic details that support previous experimental studies. Our results contribute quantitative comparison of nanoscale adhesion strength of ice on hydrophobic and hydrophilic surfaces, and supply for the first time theoretical references for understanding the mechanics at the atomistic origins of macroscale ice adhesion.Deicing is important to human activities in low-temperature circumstances, and is critical for combating the damage caused by excessive accumulation of ice. The aim of creating anti-icing materials, surfaces and applications relies on the understanding of fundamental nanoscale ice adhesion mechanics. Here in this study, we employ all-atom modeling and molecular dynamics simulation to investigate ice adhesion. We apply force to detach and shear nano-sized ice cubes for probing the determinants of atomistic adhesion mechanics, and at the same time investigate the mechanical effect of a sandwiched aqueous water layer between ice and substrates. We observe that high interfacial energy restricts ice mobility and increases both ice detaching and shearing stresses. We quantify up to a 60% decrease in ice
Molecular Dynamics Simulation of Disordered Zircon
Devanathan, Ram; Corrales, Louis R.; Weber, William J.; Chartier, Alain; Meis, Constantin
2004-02-27
The melting of zircon and the amorphous state produced by quenching from the melt were simulated by molecular dynamics using a new partial charge model combined with the Ziegler-Biersack-Littmark potential. The model has been established for the description of the crystalline and aperiodic structures of zircon in order to be used for the simulation of displacement cascades. It provides an excellent fit to the structure, and accounts with convenient precision the mechanical and thermodynamic properties of zircon. The calculated melting temperature is about 2100 K. The activation energy for self-diffusion of ions in the liquid state was determined to be 190-200 kJ/mole. Melt quenching was employed to produce two different disordered states with distinct densities and structures. In the high density disordered state, the zircon structure is intact but the bond angle distributions are broader, 4% of the Si units are polymerized, and the volume swelling is about 8%. In the low density amorphous state, the Zr and Si coordination numbers are lower, and the Zr-O and Si-O bond lengths are shorter than corresponding values for the crystal. In addition, a highly polymerized Si network, with average connectivity of two, is observed in the low density amorphous state. These features have all been experimentally observed in natural metamict zircon. The present findings, when considered in light of experimental radiation effects studies, suggest that the swelling in zircon arises initially from disorder in the zircon crystal, and at high doses the disordered crystal is unable to accommodate the volume expansion and transforms to the amorphous state.
Molecular dynamics simulations of supramolecular polymer rheology
NASA Astrophysics Data System (ADS)
Li, Zhenlong; Djohari, Hadrian; Dormidontova, Elena E.
2010-11-01
Using equilibrium and nonequilibrium molecular dynamics simulations, we studied the equilibrium and rheological properties of dilute and semidilute solutions of head-to-tail associating polymers. In our simulation model, a spontaneous complementary reversible association between the donor and the acceptor groups at the ends of oligomers was achieved by introducing a combination of truncated pseudo-Coulombic attractive potential and Lennard Jones repulsive potential between donor, acceptor, and neighboring groups. We have calculated the equilibrium properties of supramolecular polymers, such as the ring/chain equilibrium, average molecular weight, and molecular weight distribution of self-assembled chains and rings, which all agree well with previous analytical and computer modeling results. We have investigated shear thinning of solutions of 8- and 20-bead associating oligomers with different association energies at different temperatures and oligomer volume fractions. All reduced viscosity data for a given oligomer length can be collapsed into one master curve, exhibiting two power-law regions of shear-thinning behavior with an exponent of -0.55 at intermediate ranges of the reduced shear rate β and -0.8 (or -0.9) at larger shear rates. The equilibrium viscosity of supramolecular solutions with different oligomer lengths and associating energies is found to obey a power-law scaling dependence on oligomer volume fraction with an exponent of 1.5, in agreement with the experimental observations for several dilute or semidilute solutions of supramolecular polymers. This implies that dilute and semidilute supramolecular polymer solutions exhibit high polydispersity but may not be sufficiently entangled to follow the reptation mechanism of relaxation.
Numerical simulation of magma chamber dynamics.
NASA Astrophysics Data System (ADS)
Longo, Antonella; Papale, Paolo; Montagna, Chiara Paola; Vassalli, Melissa; Giudice, Salvatore; Cassioli, Andrea
2010-05-01
Magma chambers are characterized by periodic arrivals of deep magma batches that give origin to complex patterns of magma convection and mixing, and modify the distribution of physical quantities inside the chamber. We simulate the transient, 2D, multi-component homogeneous dynamics in geometrically complex dyke+chamber systems, by means of GALES, a finite element parallel C++ code solving mass, momentum and energy equations for multi-component homogeneous gas-liquid (± crystals) mixtures in compressible-to-incompressible flow conditions. Code validation analysis includes several cases from the classical engineering literature, corresponding to a variety of subsonic to supersonic gas-liquid flow regimes (see http://www.pi.ingv.it/~longo/gales/gales.html). The model allows specification of the composition of the different magmas in the domain, in terms of ten major oxides plus the two volatile species H2O and CO2. Gas-liquid thermodynamics are modeled by using the compositional dependent, non-ideal model in Papale et al. (Chem.. Geol., 2006). Magma properties are defined in terms of local pressure, temperature, and composition including volatiles. Several applications are performed within domains characterized by the presence of one or more magma chambers and one or more dykes, with different geometries and characteristic size from hundreds of m to several km. In most simulations an initial compositional interface is placed at the top of a feeding dyke, or at larger depth, with the deeper magma having a lower density as a consequence of larger volatile content. The numerical results show complex patterns of magma refilling in the chamber, with alternating phases of magma ingression and magma sinking from the chamber into the feeding dyke. Intense mixing takes place in feeding dykes, so that the new magma entering the chamber is always a mixture of the deep and the initially resident magma. Buoyant plume rise occurs through the formation of complex convective
Constraint Force Equation Methodology for Modeling Multi-Body Stage Separation Dynamics
NASA Technical Reports Server (NTRS)
Toniolo, Matthew D.; Tartabini, Paul V.; Pamadi, Bandu N.; Hotchko, Nathaniel
2008-01-01
This paper discusses a generalized approach to the multi-body separation problems in a launch vehicle staging environment based on constraint force methodology and its implementation into the Program to Optimize Simulated Trajectories II (POST2), a widely used trajectory design and optimization tool. This development facilitates the inclusion of stage separation analysis into POST2 for seamless end-to-end simulations of launch vehicle trajectories, thus simplifying the overall implementation and providing a range of modeling and optimization capabilities that are standard features in POST2. Analysis and results are presented for two test cases that validate the constraint force equation methodology in a stand-alone mode and its implementation in POST2.
Simulation of chemical isomerization reaction dynamics on a NMR quantum simulator.
Lu, Dawei; Xu, Nanyang; Xu, Ruixue; Chen, Hongwei; Gong, Jiangbin; Peng, Xinhua; Du, Jiangfeng
2011-07-01
Quantum simulation can beat current classical computers with minimally a few tens of qubits. Here we report an experimental demonstration that a small nuclear-magnetic-resonance quantum simulator is already able to simulate the dynamics of a prototype laser-driven isomerization reaction using engineered quantum control pulses. The experimental results agree well with classical simulations. We conclude that the quantum simulation of chemical reaction dynamics not computable on current classical computers is feasible in the near future. PMID:21797586
A Computational Methodology for Simulating Thermal Loss Testing of the Advanced Stirling Convertor
NASA Technical Reports Server (NTRS)
Reid, Terry V.; Wilson, Scott D.; Schifer, Nicholas A.; Briggs, Maxwell H.
2011-01-01
Convertor and generator testing is carried out in tests designed to characterize convertor performance when subjected to environments intended to simulate launch and space conditions. The value of net heat input must be known in order to calculate convertor efficiency and to validate convertor performance. Specially designed test hardware was used to verify and validate a two step methodology for the prediction of net heat input. This lessons learned from these simulations have been applied to previous convertor simulations. As heat is supplied to the convertors, electric power is produced and measured. Net heat input to the convertor is one parameter that will contribute to the calculation of efficiency. This parameter is not measured directly. Insulation Loss. Determine the current status of the thermal conductivity of the micro-porous insulation. - Match heat source and hot end temperatures. - Match temperature difference across Kaowool insulation
Staggered solution procedures for multibody dynamics simulation
NASA Technical Reports Server (NTRS)
Park, K. C.; Chiou, J. C.; Downer, J. D.
1990-01-01
The numerical solution procedure for multibody dynamics (MBD) systems is termed a staggered MBD solution procedure that solves the generalized coordinates in a separate module from that for the constraint force. This requires a reformulation of the constraint conditions so that the constraint forces can also be integrated in time. A major advantage of such a partitioned solution procedure is that additional analysis capabilities such as active controller and design optimization modules can be easily interfaced without embedding them into a monolithic program. After introducing the basic equations of motion for MBD system in the second section, Section 3 briefly reviews some constraint handling techniques and introduces the staggered stabilized technique for the solution of the constraint forces as independent variables. The numerical direct time integration of the equations of motion is described in Section 4. As accurate damping treatment is important for the dynamics of space structures, we have employed the central difference method and the mid-point form of the trapezoidal rule since they engender no numerical damping. This is in contrast to the current practice in dynamic simulations of ground vehicles by employing a set of backward difference formulas. First, the equations of motion are partitioned according to the translational and the rotational coordinates. This sets the stage for an efficient treatment of the rotational motions via the singularity-free Euler parameters. The resulting partitioned equations of motion are then integrated via a two-stage explicit stabilized algorithm for updating both the translational coordinates and angular velocities. Once the angular velocities are obtained, the angular orientations are updated via the mid-point implicit formula employing the Euler parameters. When the two algorithms, namely, the two-stage explicit algorithm for the generalized coordinates and the implicit staggered procedure for the constraint Lagrange
Continuum and molecular-dynamics simulation of nanodroplet collisions
NASA Astrophysics Data System (ADS)
Bardia, Raunak; Liang, Zhi; Keblinski, Pawel; Trujillo, Mario F.
2016-05-01
The extent to which the continuum treatment holds in binary droplet collisions is examined in the present work by using a continuum-based implicit surface capturing strategy (volume-of-fluid coupled to Navier-Stokes) and a molecular dynamics methodology. The droplet pairs are arranged in a head-on-collision configuration with an initial separation distance of 5.3 nm and a velocity of 3 ms-1. The size of droplets ranges from 10-50 nm. Inspecting the results, the collision process can be described as consisting of two periods: a preimpact phase that ends with the initial contact of both droplets, and a postimpact phase characterized by the merging, deformation, and coalescence of the droplets. The largest difference between the continuum and molecular dynamics (MD) predictions is observed in the preimpact period, where the continuum-based viscous and pressure drag forces significantly overestimate the MD predictions. Due to large value of Knudsen number in the gas (Kngas=1.972 ), this behavior is expected. Besides the differences between continuum and MD, it is also observed that the continuum simulations do not converge for the set of grid sizes considered. This is shown to be directly related to the initial velocity profile and the minute size of the nanodroplets. For instance, for micrometer-size droplets, this numerical sensitivity is not an issue. During the postimpact period, both MD and continuum-based simulations are strikingly similar, with only a moderate difference in the peak kinetic energy recorded during the collision process. With values for the Knudsen number in the liquid (Knliquid=0.01 for D =36 nm ) much closer to the continuum regime, this behavior is expected. The 50 nm droplet case is sufficiently large to be predicted reasonably well with the continuum treatment. However, for droplets smaller than approximately 36 nm, the departure from continuum behavior becomes noticeably pronounced, and becomes drastically different for the 10 nm droplets.
Fogarty, Aoife C. Potestio, Raffaello Kremer, Kurt
2015-05-21
A fully atomistic modelling of many biophysical and biochemical processes at biologically relevant length- and time scales is beyond our reach with current computational resources, and one approach to overcome this difficulty is the use of multiscale simulation techniques. In such simulations, when system properties necessitate a boundary between resolutions that falls within the solvent region, one can use an approach such as the Adaptive Resolution Scheme (AdResS), in which solvent particles change their resolution on the fly during the simulation. Here, we apply the existing AdResS methodology to biomolecular systems, simulating a fully atomistic protein with an atomistic hydration shell, solvated in a coarse-grained particle reservoir and heat bath. Using as a test case an aqueous solution of the regulatory protein ubiquitin, we first confirm the validity of the AdResS approach for such systems, via an examination of protein and solvent structural and dynamical properties. We then demonstrate how, in addition to providing a computational speedup, such a multiscale AdResS approach can yield otherwise inaccessible physical insights into biomolecular function. We use our methodology to show that protein structure and dynamics can still be correctly modelled using only a few shells of atomistic water molecules. We also discuss aspects of the AdResS methodology peculiar to biomolecular simulations.
System dynamic simulation of precision segmented reflector
NASA Technical Reports Server (NTRS)
Shih, Choon-Foo; Lou, Michael C.
1991-01-01
A joint effort was undertaken on a Precision Segmented Reflector (PSR) Project. The missions in which the PSR is to be used will use large (up to 20 m in diameter) telescopes. The essential requirement for the telescopes is that the reflector surface of the primary mirror must be made extremely precise to allow no more than a few microns of errors and, additionally, this high surface precision must be maintained when the telescope is subjected to on-orbital mechanical and thermal disturbances. Based on the mass, size, and stability considerations, reflector surface formed by segmented, probably actively or passively controlled, composite panels are regarded as most suitable for future space based astronomical telescope applications. In addition to the design and fabrication of composite panels with a surface error of less than 3 microns RMS, PSR also develops related reflector structures, materials, control, and sensing technologies. As part of the planning effort for PSR Technology Demonstration, a system model which couples the reflector, consisting of panels, support truss and actuators, and the optical bench was assembled for dynamic simulations. Random vibration analyses using seismic data obtained from actual measurements at the test site designated for PSR Technology Demonstration are described.
Kinetic simulations of plasmoid chain dynamics
Markidis, S.; Henri, P.; Lapenta, G.; Divin, A.; Goldman, M.; Newman, D.; Laure, E.
2013-08-15
The dynamics of a plasmoid chain is studied with three dimensional Particle-in-Cell simulations. The evolution of the system with and without a uniform guide field, whose strength is 1/3 the asymptotic magnetic field, is investigated. The plasmoid chain forms by spontaneous magnetic reconnection: the tearing instability rapidly disrupts the initial current sheet generating several small-scale plasmoids that rapidly grow in size coalescing and kinking. The plasmoid kink is mainly driven by the coalescence process. It is found that the presence of guide field strongly influences the evolution of the plasmoid chain. Without a guide field, a main reconnection site dominates and smaller reconnection regions are included in larger ones, leading to an hierarchical structure of the plasmoid-dominated current sheet. On the contrary in presence of a guide field, plasmoids have approximately the same size and the hierarchical structure does not emerge, a strong core magnetic field develops in the center of the plasmoid in the direction of the existing guide field, and bump-on-tail instability, leading to the formation of electron holes, is detected in proximity of the plasmoids.
Molecular Dynamics Simulations of Coulomb Explosion
Bringa, E M
2002-05-17
A swift ion creates a track of electronic excitations in the target material. A net repulsion inside the track can cause a ''Coulomb Explosion'', which can lead to damage and sputtering of the material. Here we report results from molecular-dynamics (MD) simulations of Coulomb explosion for a cylindrical track as a function of charge density and neutralization/quenching time, {tau}. Screening by the free electrons is accounted for using a screened Coulomb potential for the interaction among charges. The yield exhibits a prompt component from the track core and a component, which dominates at higher excitation density, from the heated region produced. For the cases studied, the number of atoms ejected per incident ion, i.e. the sputtering yield Y, is quadratic with charge density along the track as suggested by simple models. Y({tau} = 0.2 Debye periods) is nearly 20% of the yield when there is no neutralization ({tau} {yields} {infinity}). The connections between ''Coulomb explosions'', thermal spikes and measurements of electronic sputtering are discussed.
Annual Report 1999 Environmental Dynamics and Simulation
NS Foster-Mills
2000-06-28
This annual report describes selected 1999 research accomplishments for the Environmental Dynamics and Simulation (ED and S) directorate, one of six research organizations in the William R. Wiley Environmental Molecular Sciences Laboratory (EMSL). These accomplishments are representative of the different lines of research underway in the ED and S directorate. EMSL is one of US Department of Energy's (DOE) national scientific user facilities and is the centerpiece of DOE's commitment to providing world-class experimental, theoretical, and computational capabilities for solving the nation's environmental problems. Capabilities in the EMSL include over 100 major instrument systems for use by the resident research staff, their collaborators, and users of the EMSL. These capabilities are used to address the fundamental science that will be the basis for finding solutions to national environmental issues such as cleaning up contamianted areas at DOE sites across the country and developing green technologies that will reduce or eliminate future pollution production. The capabilities are also used to further the understanding of global climate change and environmental issues relevant to energy production and use and health effects resulting from exposure to contaminated environments.
Expansion techniques for collisionless stellar dynamical simulations
Meiron, Yohai; Li, Baile; Holley-Bockelmann, Kelly; Spurzem, Rainer
2014-09-10
We present graphics processing unit (GPU) implementations of two fast force calculation methods based on series expansions of the Poisson equation. One method is the self-consistent field (SCF) method, which is a Fourier-like expansion of the density field in some basis set; the other method is the multipole expansion (MEX) method, which is a Taylor-like expansion of the Green's function. MEX, which has been advocated in the past, has not gained as much popularity as SCF. Both are particle-field methods and optimized for collisionless galactic dynamics, but while SCF is a 'pure' expansion, MEX is an expansion in just the angular part; thus, MEX is capable of capturing radial structure easily, while SCF needs a large number of radial terms. We show that despite the expansion bias, these methods are more accurate than direct techniques for the same number of particles. The performance of our GPU code, which we call ETICS, is profiled and compared to a CPU implementation. On the tested GPU hardware, a full force calculation for one million particles took ∼0.1 s (depending on expansion cutoff), making simulations with as many as 10{sup 8} particles fast for a comparatively small number of nodes.
Nanoscale deicing by molecular dynamics simulation.
Xiao, Senbo; He, Jianying; Zhang, Zhiliang
2016-08-14
Deicing is important to human activities in low-temperature circumstances, and is critical for combating the damage caused by excessive accumulation of ice. The aim of creating anti-icing materials, surfaces and applications relies on the understanding of fundamental nanoscale ice adhesion mechanics. Here in this study, we employ all-atom modeling and molecular dynamics simulation to investigate ice adhesion. We apply force to detach and shear nano-sized ice cubes for probing the determinants of atomistic adhesion mechanics, and at the same time investigate the mechanical effect of a sandwiched aqueous water layer between ice and substrates. We observe that high interfacial energy restricts ice mobility and increases both ice detaching and shearing stresses. We quantify up to a 60% decrease in ice adhesion strength by an aqueous water layer, and provide atomistic details that support previous experimental studies. Our results contribute quantitative comparison of nanoscale adhesion strength of ice on hydrophobic and hydrophilic surfaces, and supply for the first time theoretical references for understanding the mechanics at the atomistic origins of macroscale ice adhesion. PMID:27431975
NASA Technical Reports Server (NTRS)
Przekwas, A. J.; Athavale, M. M.; Hendricks, R. C.; Steinetz, B. M.
2006-01-01
Detailed information of the flow-fields in the secondary flowpaths and their interaction with the primary flows in gas turbine engines is necessary for successful designs with optimized secondary flow streams. Present work is focused on the development of a simulation methodology for coupled time-accurate solutions of the two flowpaths. The secondary flowstream is treated using SCISEAL, an unstructured adaptive Cartesian grid code developed for secondary flows and seals, while the mainpath flow is solved using TURBO, a density based code with capability of resolving rotor-stator interaction in multi-stage machines. An interface is being tested that links the two codes at the rim seal to allow data exchange between the two codes for parallel, coupled execution. A description of the coupling methodology and the current status of the interface development is presented. Representative steady-state solutions of the secondary flow in the UTRC HP Rig disc cavity are also presented.
Dynamic Simulation over Long Time Periods with 100% Solar Generation.
Concepcion, Ricky James; Elliott, Ryan Thomas
2015-12-01
This project aimed to identify the path forward for dynamic simulation tools to accommodate these needs by characterizing the properties of power systems (with high PV penetration), analyzing how these properties affect dynamic simulation software, and offering solutions for potential problems.
Sezer, Deniz; Freed, Jack H; Roux, Benoit
2008-09-01
Simulating electron spin resonance (ESR) spectra directly from molecular dynamics simulations of a spin-labeled protein necessitates a large number (hundreds or thousands) of relatively long (hundreds of nanoseconds) trajectories. To meet this challenge, we explore the possibility of constructing accurate stochastic models of the spin label dynamics from atomistic trajectories. A systematic, two-step procedure, based on the probabilistic framework of hidden Markov models, is developed to build a discrete-time Markov chain process that faithfully captures the internal spin label dynamics on time scales longer than about 150 ps. The constructed Markov model is used both to gain insight into the long-lived conformations of the spin label and to generate the stochastic trajectories required for the simulation of ESR spectra. The methodology is illustrated with an application to the case of a spin-labeled poly alanine alpha helix in explicit solvent. PMID:18698714