SEARCH · Engineering Papers

Results for “core simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,207 records · Page 67

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Investigation of the OECD/NEA PWR MOX/UO{sub 2} core transient benchmark using a coupled whole core pin-by-pin route in WIMS

This paper investigates the new CAMELOT coupling route developed in WIMS, which uses prepared tabulated cross-sections combined with the flux solver MERLIN module and the integrated sub-channel thermal hydraulics solver ARTHUR module to solve for the coupled core state. This coupling route is applied to the OECD/NEA pressurized water reactor (PWR) core transient benchmark exercise, involving 3-D neutronics and thermal-hydraulics analysis of a full-sized PWR core featuring MOX fuel at various core states, as well as simulating a control rod ejection transient. Analysis of the results for each part of the benchmark demonstrate that the high-resolution pin-by-pin analysis used in WIMS closely matches the other participants in the benchmark, with differences of between 20-200 pcm for k-effective, 1-2.5% RMS difference for assembly-averaged powers, 2-8% difference in rod worth calculations, < 200 ppm for calculations of critical boron concentration, zero difference in delayed neutron fraction calculation and a very close match for power response to the rod ejection transient. This is the first application of the WIMS-CAMELOT approach for the transient analysis of a full-sized reactor and as such further improvements such as reducing computational cost through parallelization, optimisation and memory reduction, are currently in active development. This paper demonstrates the accuracy and flexibility of the CAMELOT coupling route, providing a straightforward and consistent user image that can be easily used for a variety of modelling problems involving coupled neutronics and thermal-hydraulics. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2025

Wellbore Stability and Mud Loss Management in Geothermal Drilling: Optimizing Mud Weight to Mitigate Tensile Wellbore Fracturing at The Geysers, California

As part of a U.S. Department of Energy (DOE) Geothermal Technologies Office-funded initiative, Geysers Power Company, LLC, a subsidiary of Calpine Corporation, has been working to enhance drilling performance at the world’s largest geothermal field, The Geysers, in northern California. In a recent drilling operation of the GDC-36 well, excessive mud losses were encountered, initially addressed through repeated but largely ineffective cement plugging. Ultimately, the most effective strategy was to drill blind through the loss zones, made feasible by the high rate of penetration (ROP) achieved with PDC bits, allowing significant progress before the mud tanks were depleted and water-sensitive argillic formation layers could collapse. In response to these challenges, the project team explored alternative methods to minimize downtime and risks associated with cement plugging and continuous mud loss and to contemplate the driving mechanisms for the losses. Wellbore imaging using Formation MicroImager (FMI) and Ultrasonic Borehole Imager (UBI) tools revealed longitudinal tensile fractures, which were attributed to mud weights exceeding the minimum circumferential stress resulting from the native stress field and formation pressure. This study examines the mud losses encountered and leverages wellbore imaging data to understand the mechanisms behind mud induced tensile fracturing in specific rock facies. Understanding fracture behavior across different lithologies is crucial, as fractures within the reservoir can enhance steam migration throughout the system. The reservoir at The Geysers lies within the Mesozoic Franciscan Assemblage, a tectonic mélange formed by subduction. It consists of metamorphosed turbidite sandstone (greywacke) and mudstone (argillite), oceanic upper crust (including greenstone and chert), and serpentinized ultramafic rocks - each exhibiting distinct geomechanical fracturing properties. The structural fabric of the Franciscan Assemblage was shaped by low-angle Mesozoic thrust faulting and later overprinted by sub-vertical strike-slip structures related to the Pacific-North American plate boundary. A wellbore stability model was developed using core measurements and logs to simulate fracturing scenarios during drilling under varying stress conditions. These simulations guided the development of an optimized mud weight management strategy that should enable adaptive adjustments during drilling, reducing the likelihood of tensile fracturing and mud losses, ultimately improving operational efficiency.

15 GEOTHERMAL ENERGY↗

DOE OSTI2025

$S^5$: Tidal Disruption in Crater 2 and Formation of Diffuse Dwarf Galaxies in the Local Group

We present results of a spectroscopic campaign around the diffuse dwarf galaxy Crater 2 (Cra2) and its tidal tails as part of the Southern Stellar Stream Spectroscopic Survey ($S^5$). Cra2 is a Milky Way dwarf spheroidal satellite with extremely cold kinematics, but a huge size similar to the Small Magellanic Cloud, which may be difficult to explain within collisionless cold dark matter. We identify 143 Cra2 members, of which 114 belong to the galaxy's main body and 29 are deemed part of its stellar stream. We confirm that Cra2 is dynamically cold (central velocity dispersion $2.51^{+0.33}_{-0.30}\,{\rm km\,s^{-1}}$) and also discover a $\approx$7$Ï$ velocity gradient consistent with its tidal debris track. We separately estimate the stream velocity dispersion to be $5.74^{+0.98}_{-0.83}\,{\rm km\,s^{-1}}$. We develop a suite of $N$-body simulations with both cuspy and cored density profiles on a realistic Cra2 orbit to compare with $S^5$ observations. We find that the velocity dispersion ratio between Cra2 stream and galaxy ($2.30^{+0.41}_{-0.35}$) is difficult to reconcile with a cuspy halo with fiducial concentration and an initial mass predicted by standard stellar mass$-$halo mass relationships. Instead, either a cored halo with relatively small core radius or a low-concentration cuspy model can reproduce this ratio. Despite tidal mass loss, Cra2 is metal-poor ($\langle \rm[Fe/H]\rangle=-2.16\pm0.04$) compared to the stellar mass$-$metallicity relation for its luminosity. Other diffuse dwarf galaxies similar to Cra2 in the Local Group (Antlia 2 and Andromeda 19) also challenge galaxy formation models. Finally, we discuss possible formation scenarios for Cra2, including ram-pressure stripping of a gas-rich progenitor combined with tides.

Limberg, Guilherme [Chicago U., KICP; Chicago U.] ↗

DOE OSTI2014

State space modeling of reactor core in a pressurized water reactor

The power control system of a nuclear reactor is the key system that ensures a safe operation for a nuclear power plant. However, a mathematical model of a nuclear power plant is in the form of nonlinear process and time dependent that give very hard to be described. One of the important components of a Pressurized Water Reactor is the Reactor core. The aim of this study is to analyze the performance of power produced from a reactor core using temperature of the moderator as an input. Mathematical representation of the state space model of the reactor core control system is presented and analyzed in this paper. The data and parameters are taken from a real time VVER-type Pressurized Water Reactor and will be verified using Matlab and Simulink. Based on the simulation conducted, the results show that the temperature of the moderator plays an important role in determining the power of reactor core.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Multi-scale gyrokinetic simulations of an Alcator C-Mod, ELM-y H-mode plasma

High fidelity, multi-scale gyrokinetic simulations capable of capturing both ion (${k}_{\theta }{\rho }_{s}\sim { \mathcal O }(1.0)$) and electron-scale (${k}_{\theta }{\rho }_{e}\sim { \mathcal O }(1.0)$) turbulence were performed in the core of an Alcator C-Mod ELM-y H-mode discharge which exhibits reactor-relevant characteristics. These simulations, performed with all experimental inputs and realistic ion to electron mass ratio (${({m}_{i}/{m}_{e})}^{1/2}=60.0$) provide insight into the physics fidelity that may be needed for accurate simulation of the core of fusion reactor discharges. Three multi-scale simulations and series of separate ion and electron-scale simulations performed using the GYRO code are presented. As with earlier multi-scale results in L-mode conditions, both ion and multi-scale simulations results are compared with experimentally inferred ion and electron heat fluxes, as well as the measured values of electron incremental thermal diffusivities—indicative of the experimental electron temperature profile stiffness. Consistent with the L-mode results, cross-scale coupling is found to play an important role in the simulation of these H-mode conditions. Extremely stiff ion-scale transport is observed in these high-performance conditions which is shown to likely play and important role in the reproduction of measurements of perturbative transport. In conclusion, these results provide important insight into the role of multi-scale plasma turbulence in the core of reactor-relevant plasmas and establish important constraints on the fidelity of models needed for predictive simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

SDSS-IV MaNGA: Cannibalism Caught in the Act—On the Frequency of Occurrence of Multiple Cores in Brightest Cluster Galaxies

Although it is generally accepted that massive galaxies form in a two-phased fashion, beginning with a rapid mass buildup through intense starburst activities followed by primarily dry mergers that mainly deposit stellar mass at outskirts, the late time stellar mass growth of brightest cluster galaxies (BCGs), the most massive galaxies in the universe, is still not well understood. Several independent measurements have indicated a slower mass growth rate than predictions from theoretical models. We attempt to resolve the discrepancy by measuring the frequency of BCGs with multiple cores, which serve as a proxy of the merger rates in the central region and facilitate a more direct comparison with theoretical predictions. Using 79 BCGs at z = 0.06–0.15 with integral field spectroscopic data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) project, we obtain a multiple-core fraction of 0.11 ± 0.04 at z ≈ 0.1 within an 18 kpc radius from the center, which is comparable to the value of 0.08 ± 0.04 derived from mock observations of 218 simulated BCGs from the cosmological hydrodynamical simulation IllustrisTNG. We find that most cores that appear close to the BCGs from imaging data turn out to be physically associated systems. Anchoring on the similarity in the multiple-core frequency between the MaNGA and IllustrisTNG, we discuss the mass growth rate of BCGs over the past 4.5 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Integrating static PRA information with risk informed safety margin characterization (RISMC) simulation methods

The overall objective of the project was to develop a computationally feasible and user-friendly mechanized process to integrate traditional probabilistic risk assessment (PRA) and dynamic PRA (DPRA) results. Starting with the systematic identification of items in an existing PRA that need dynamic augmentation, the project used a generic 4-loop pressurized reactor (PWR) and 3-loop PWR as example plants. Station blackout (SBO) and large break loss of coolant accident SBLOCA) were selected as the example initiating events. Using the traditional event-tree (ET)/fault-tree (FT) methodology augmented by dynamic evet tree approach, the potential consequences of the initiating events were simulated with RELAP-3D and MELCOR/RASCAL codes to cover Level 1 through Level 3 of PRA. RAVEN and ADAPT software were used to generate Level 1 simulations with RELAP-3D and Level 2/3 simulations with MELCOR (Level 2)/RASCAL (Level 3), respectively. Example branching conditions (BCs) for SBO included AC power recovery time, valve repair failure time, reactor coolant pump leak time/break size and emergency power supply duration to a total of 9. Example BCs for LOCA included off-site power recovery time, diesel generator power recovery time, auxiliary feed water system operation time, safety relief valve failure to open upon demand, reactor coolant pump seal break time and size to a total of 21. Each RELAP-3D simulation (9,587 scenarios) was labelled OK or Core Damage based on the maximum allowed peak clad temperature (2,100oF). Each MELCOR simulation (4610 scenarios) was labeled as Bin over 10rem or Bin 0-10rem based on the dose at the site boundary. The scenarios were clustered based on the criteria above using the mean shift methodology. Classical PRA (CPRA) and DPRA results were compared to identify the ET sequences that need DPRA augmentation. Several approaches were proposed for the incorporation of these sequences into CPRA using clustering with the mean shift methodology, restructuring the CPRA ETs by adding new BCs/sequences, and using the concept of a limit surface. Procedures for decision making regarding the possible consequences of an initiating event (e.g. core damage or not, site evacuation or not) were developed using a convolutional neural network (CNN), a recurrent neural network (RNN) and a transformer neural network (TNN). The project has led to two PhD degrees, three archival journal papers and five refereed conference proceedings.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS2012

Next Generation Simulation Framework for Robotic and Human Space Missions

The Dartslab team at NASA's Jet Propulsion Laboratory (JPL) has a long history of developing physics-based simulations based on the Darts/Dshell simulation framework that have been used to simulate many planetary robotic missions, such as the Cassini spacecraft and the rovers that are currently driving on Mars. Recent collaboration efforts between the Dartslab team at JPL and the Mission Operations Directorate (MOD) at NASA Johnson Space Center (JSC) have led to significant enhancements to the Dartslab DSENDS (Dynamics Simulator for Entry, Descent and Surface landing) software framework. The new version of DSENDS is now being used for new planetary mission simulations at JPL. JSC is using DSENDS as the foundation for a suite of software known as COMPASS (Core Operations, Mission Planning, and Analysis Spacecraft Simulation) that is the basis for their new human space mission simulations and analysis. In this paper, we will describe the collaborative process with the JPL Dartslab and the JSC MOD team that resulted in the redesign and enhancement of the DSENDS software. We will outline the improvements in DSENDS that simplify creation of new high-fidelity robotic/spacecraft simulations. We will illustrate how DSENDS simulations are assembled and show results from several mission simulations.

Dartslab↗

DOE OSTI2018

A numerical simulation of a dry-out process for CO 2 sequestration in heterogeneous deep saline aquifers: Original Research Article: Dry-out problem: Multifactorial numerical simulation

Salt precipitation can be generated near the injection well when dry supercritical CO 2 (scCO 2 ) is injected into deep saline aquifers. Traditional laboratory experiments and numerical simulations focus on the generation of salt precipitation and its impact on core permeability. Current laboratory experimental studies have shown that the heterogeneity of the core and water content on the surface of the particles plays an important role in the formation of salt precipitation. The aim of this study is to explore the effects of brine back-flow, heterogeneity, and low water content on salt precipitation. Numerical simulations were performed using core-scale and site-scale models. Simulation results showed the following: 1. Back-flow plays an important role in the salt accumulation near the well – an homogeneous laboratory experiment used with a small core may not provide a full picture of the back-flow phenomenon. 2. The continuous dry-out process reduces the water content near the injection well to zero. As a result, the full range of saturation flow equations should be used to simulate the dry-out process accurately. 3. The heterogeneity of the matrix not only affects the distribution of salt precipitation but also increases the amount of salt precipitation near the well. Homogeneous simulation underestimates the amount of salt precipitation. 4. The reduction in porosity due to salt precipitation could extend the dry-out process in the low permeability zone, further expanding the salt precipitation area.

58 GEOSCIENCES↗

DOE OSTI2013

Integrated Simulation Development and Decision Support Tool-Set for Utility Market and Distributed Solar Power Generation

This project will enable utilities to develop long-term strategic plans that integrate high levels of renewable energy generation, and to better plan power system operations under high renewable penetration. The program developed forecast data streams for decision support and effective integration of centralized and distributed solar power generation in utility operations. This toolset focused on real time simulation of distributed power generation within utility grids with the emphasis on potential applications in day ahead (market) and real time (reliability) utility operations. The project team developed and demonstrated methodologies for quantifying the impact of distributed solar generation on core utility operations, identified protocols for internal data communication requirements, and worked with utility personnel to adapt the new distributed generation (DG) forecasts seamlessly within existing Load and Generation procedures through a sophisticated DMS. This project supported the objectives of the SunShot Initiative and SUNRISE by enabling core utility operations to enhance their simulation capability to analyze and prepare for the impacts of high penetrations of solar on the power grid. The impact of high penetration solar PV on utility operations is not only limited to control centers, but across many core operations. Benefits of an enhanced DMS using state-of-the-art solar forecast data were demonstrated within this project and have had an immediate direct operational cost savings for Energy Marketing for Day Ahead generation commitments, Real Time Operations, Load Forecasting (at an aggregate system level for Day Ahead), Demand Response, Long term Planning (asset management), Distribution Operations, and core ancillary services as required for balancing and reliability. This provided power system operators with the necessary tools and processes to operate the grid in a reliable manner under high renewable penetration.

14 SOLAR ENERGY↗

DOE OSTI2014

Final Technical Report for Contract No. DE-EE0006332, "Integrated Simulation Development and Decision Support Tool-Set for Utility Market and Distributed Solar Power Generation"

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2015

Non-universality of dark-matter halos: Cusps, cores, and the central potential

Dark-matter halos grown in cosmological simulations appear to have central Navarro–Frenk–White-like density cusps with mean values of dlogρ/dlogr≈−1, and some dispersion, which is generally parametrized by the varying index α in the Einasto density profile fitting function. Non-universality in profile shapes is also seen in observed galaxy clusters and possibly dwarf galaxies. Here we show that non-universality, at any given mass scale, is an intrinsic property of DARKexp, a theoretically derived model for collisionless self-gravitating systems. We demonstrate that DARKexp—which has only one shape parameter, ϕ{sub 0}—fits the dispersion in profile shapes of massive simulated halos as well as observed clusters very well. DARKexp also allows for cored dark-matter profiles, such as those found for dwarf spheroidal galaxies. We provide approximate analytical relations between DARKexp ϕ{sub 0}, Einasto α, or the central logarithmic slope in the Dehnen–Tremaine analytical γ-models. The range in halo parameters reflects a substantial variation in the binding energies per unit mass of dark-matter halos.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Containers for Massive Ensemble of I/O Bound Hierarchical Coupled Simulations

We present our experience using containers to scale up a massive ensemble of coupled I/O bound workloads on the NERSC Cori supercomputer. We describe the design of a hierarchical simulation structure using the Integrated Plasma Simulator (IPS) that enables the flexible execution of coupled simulations at the system, node, and core level using the same coupling abstraction and API. The hierarchical design allows for the node-level execution to be efficiently executed using containers while not impacting the structure of the simulation at the system level. We demonstrate the viability of the approach by presenting experimental results from applications in coupled fusion plasma simulations that illustrate the performance impact of using containers to deploy the node-level workloads, in conjunction with the user mountable XFS file systems to ameliorate the load on the Lustre parallel file system. We also present results from production runs showing the ability of the ensemble simulations to scale to hundreds of Cori Haswell nodes, with little or no overhead.

Elwasif, Wael↗

DOE OSTI2016

Severe Accident Modeling for Cyber Scenarios

A generic severe accident model of a PWR is developed to support analyses of cyber-fault scenarios that may lead to extensive core damage. For this purpose, a MELCOR 2.1 model of a generic PWR is created that is capable of simulating primary and secondary thermal-hydraulics, core heatup and fuel degradation, RPV failure, containment response, and radionuclide release and transport to the environment. The size and detail of the severe accident model are intended to facilitate scoping analyses of various cyber-based scenarios in an automated, dynamic framework. A dynamic job scheduler such as ADAPT will be used to evaluate many accident scenarios in an automated fashion utilizing parallel computing. Individual test calculations are first conducted of a basic unmitigated scenario (i.e. no significant operator intervention) with a cyber-based accident initiator. This preliminary assessment is used to gauge modelling capabilities for further, expanded cyber accident scenarios in conjunction with ADAPT. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2016

THE LAST MINUTES OF OXYGEN SHELL BURNING IN A MASSIVE STAR

We present the first 4 π– three-dimensional (3D) simulation of the last minutes of oxygen shell burning in an 18 M {sub ⊙} supernova progenitor up to the onset of core collapse. A moving inner boundary is used to accurately model the contraction of the silicon and iron core according to a one-dimensional stellar evolution model with a self-consistent treatment of core deleptonization and nuclear quasi-equilibrium. The simulation covers the full solid angle to allow the emergence of large-scale convective modes. Due to core contraction and the concomitant acceleration of nuclear burning, the convective Mach number increases to ∼0.1 at collapse, and an ℓ = 2 mode emerges shortly before the end of the simulation. Aside from a growth of the oxygen shell from 0.51 M {sub ⊙} to 0.56 M {sub ⊙} due to entrainment from the carbon shell, the convective flow is reasonably well described by mixing-length theory, and the dominant scales are compatible with estimates from linear stability analysis. We deduce that artificial changes in the physics, such as accelerated core contraction, can have precarious consequences for the state of convection at collapse. We argue that scaling laws for the convective velocities and eddy sizes furnish good estimates for the state of shell convection at collapse and develop a simple analytic theory for the impact of convective seed perturbations on shock revival in the ensuing supernova. We predict a reduction of the critical luminosity for explosion by 12% – 24% due to seed asphericities for our 3D progenitor model relative to the case without large seed perturbations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Enhanced Severe Transient Analysis for Prevention Technical Program Plan

This document outlines the development of a high fidelity, best estimate nuclear power plant severe transient simulation capability that will complement or enhance the integral system codes historically used for licensing and analysis of severe accidents. As with other tools in the Risk Informed Safety Margin Characterization (RISMC) Toolkit, the ultimate user of Enhanced Severe Transient Analysis and Prevention (ESTAP) capability is the plant decision-maker; the deliverable to that customer is a modern, simulation-based safety analysis capability, applicable to a much broader class of safety issues than is traditional Light Water Reactor (LWR) licensing analysis. Currently, the RISMC pathway’s major emphasis is placed on developing RELAP-7, a next-generation safety analysis code, and on showing how to use RELAP-7 to analyze margin from a modern point of view: that is, by characterizing margin in terms of the probabilistic spectra of the “loads” applied to systems, structures, and components (SSCs), and the “capacity” of those SSCs to resist those loads without failing. The first objective of the ESTAP task, and the focus of one task of this effort, is to augment RELAP-7 analyses with user-selected multi-dimensional, multi-phase models of specific plant components to simulate complex phenomena that may lead to, or exacerbate, severe transients and core damage. Such phenomena include: coolant crossflow between PWR assemblies during a severe reactivity transient, stratified single or two-phase coolant flow in primary coolant piping, inhomogeneous mixing of emergency coolant water or boric acid with hot primary coolant, and water hammer. These are well-documented phenomena associated with plant transients but that are generally not captured in system codes. They are, however, generally limited to specific components, structures, and operating conditions. The second ESTAP task is to similarly augment a severe (post-core damage) accident integral analyses code with high fidelity simulations that would allow investigation of multi-dimensional, multi-phase containment phenomena that are only treated approximately in established codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2012

Benchmarking of software and methods for use in transient multidimensional fuel performance with spatial reactor kinetics

The key physics involved in accurate prediction of reactor-fuel-element behavior includes neutron transport and thermal hydraulics. The thermal hydraulic feedback mechanism is primarily provided through cross sections to the neutron transport that are temperature and density dependent. Historically, this coupling was primarily seen only in reactor simulators, which are well suited to model the reactor core, giving only a coarse treatment to individual fuel pins as well as simple models for thermal distribution calculations. This poor resolution on the primary coupling mechanisms can lead to conservatisms that should be removed to improve fuel design and performance. This work seeks to address the resolution of space-time-dependent neutron kinetics with thermal feedback within the fuel pin scale in the multi-physics framework. The specific application of this new capability is transient performance analysis of space-time-dependent temperature distribution of fuel elements. The coupling between the neutron transport and the thermal feedback is extremely important in this highly coupled problem, primarily applicable to reactivity-initiated- accidents (RIAs) and loss-of-coolant-accidents (LOCAs). The capability developed will include the coupling of the time-dependent neutron transport with the time-dependent thermal diffusion capability. An improvement in resolution and coupling is proposed by developing neutron transport models that are internally coupled with high fidelity within fuel pin thermal calculations in a multi-physics framework. Good agreement is shown with benchmarks and problems from the literature of RIAs and LOCAs for the tools used. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2015

A molecular dynamics study of phase transition in strongly coupled pair-ion plasmas

Existence of phase transition in strongly coupled pair-ion plasmas with soft core is investigated. Extensive Molecular Dynamics (MD) simulations are performed in the canonical ensemble, for such plasmas, at different temperatures, to analyze phase stability. Our studies show interesting phase co-existence between liquid-like and vapor-like phases. The different phases are identified by calculating the ensemble averaged density. This and the corresponding critical properties are calculated directly from MD simulation. The critical temperature of vapor-liquid coexistence is obtained, and the corresponding critical value of density is also estimated for different sizes of the soft core. We have used a novel method that allows the location of phase coexistence through a constant density simulation in which the temperature is changed in a single time-step (quenching) in order to place the system in a thermodynamically and mechanically unstable state, resulting in spontaneous separation of two coexisting phases. The results obtained from this temperature quench MD method also show the coexistence of vapor-liquid phase in pair-ion plasmas. The critical exponents obtained directly from MD simulation are found to be in close agreement with the values predicted by a mean-field theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗