Engineering Papers⌕ Search

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.

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

At least 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

DEUTERIUM BURNING IN MASSIVE GIANT PLANETS AND LOW-MASS BROWN DWARFS FORMED BY CORE-NUCLEATED ACCRETION

Using detailed numerical simulations, we study the formation of bodies near the deuterium-burning limit according to the core-nucleated giant planet accretion scenario. The objects, with heavy-element cores in the range 5-30 M{sub Circled-Plus }, are assumed to accrete gas up to final masses of 10-15 Jupiter masses (M{sub Jup}). After the formation process, which lasts 1-5 Myr and which ends with a ''cold-start'', low-entropy configuration, the bodies evolve at constant mass up to an age of several Gyr. Deuterium burning via proton capture is included in the calculation, and we determined the mass, M{sub 50}, above which more than 50% of the initial deuterium is burned. This often-quoted borderline between giant planets and brown dwarfs is found to depend only slightly on parameters, such as core mass, stellar mass, formation location, solid surface density in the protoplanetary disk, disk viscosity, and dust opacity. The values for M{sub 50} fall in the range 11.6-13.6 M{sub Jup}, in agreement with previous determinations that do not take the formation process into account. For a given opacity law during the formation process, objects with higher core masses form more quickly. The result is higher entropy in the envelope at the completion of accretion, yielding lower values of M{sub 50}. For masses above M{sub 50}, during the deuterium-burning phase, objects expand and increase in luminosity by one to three orders of magnitude. Evolutionary tracks in the luminosity versus time diagram are compared with the observed position of the companion to Beta Pictoris.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A multiphysics model of the versatile test reactor based on the MOOSE framework

The traditional modeling approach for sodium fast reactor cores relies on separate physics models, where the fuel performance, thermal–hydraulics, and neutronics calculations required to predict the core physics characteristics for nominal conditions are decoupled by relying on user-imposed boundary conditions. Here, this paper aims at evaluating the impact of multiphysics simulations for predicting the core characteristics of the Versatile Test Reactor, which is being designed as a 300-MWt sodium-cooled fast reactor. The purpose of the Versatile Test Reactor is to accelerate the testing of advanced nuclear materials in the United States. The proposed multiphysics model relies on the Griffin reactor physics code, the SAM thermal–hydraulic system code, the BISON fuel performance code, as well as generic Multiphysics Object-Oriented Simulation Environment capabilities implemented in the open-source tensor mechanics module. For k eff calculations, the introduction of a tight coupling between the neutronics, thermo-mechanical and thermal–hydraulics models induces a change of around 543 pcm in the eigenvalue, compared to the traditional standalone neutronics calculation where approximate temperature profiles are used. The multiphysics model is then employed for quantifying the impact of the thermal conductivity uncertainties on some of the key figures of merit, such as the fuel centerline temperature, assembly powers, and keff for nominal core conditions. As anticipated, uncertainties on fuel thermal conductivity mostly impact the fuel centerline temperature, and to a lesser extend the k eff .

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Constraints for stellar electron-capture rates on 86 Kr via the 86 Kr(t, 3 He +γ) 86 Br reaction and the implications for core-collapse supernovae

Background: In the late stages of stellar core collapse just prior to core bounce, electron captures on medium-heavy nuclei drive deleptonization. Therefore, simulations require the use of accurate reaction rates. Nuclei with neutron number near N=50 above atomic number Z=28 play an important role. Rates presently used in astrophysical simulations rely primarily on a relatively simple single-state approximation. In order to improve the accuracy of the astrophysical simulations, experimental data are needed to test the electron-capture rates and to guide the development of better theoretical models and astrophysical simulations. Purpose: The purpose of the present work was to measure the Gamow-Teller transition strength from 86 Kr to 86 Br, to derive the stellar electron-capture rates based on the extracted strengths, and to compare the derived rates with rates based on shell-model and quasiparticle random-phase approximation (QRPA) Gamow-Teller strengths calculations, as well as the single-state approximation. An additional purpose was to test the impact of using improved electron-capture rates on the late evolution of core-collapse supernovae. Method: The Gamow-Teller strengths from 86 Kr were extracted from the 86 Kr(t, 3 He +γ) charge-exchange reaction at 115MeV/u. The electron-capture rates were calculated as a function of stellar density and temperature. Besides the case of 86 Kr, the electron-capture rates based on the QRPA calculations were calculated for 78 additional isotopes near N=50 above Z=28. The impact of using these rates instead of those based on the single-state approximation is studied in a spherically symmetrical simulation of core collapse just prior to bounce. Results: The derived electron-capture rates on 86 Kr from the experimental Gamow-Teller strength distribution are much smaller than the rates estimated based on the single-state approximation. Rates based on Gamow-Teller strengths estimated in shell-model and QRPA calculations are more accurate. The core-collapse supernova simulation with electron-capture rates based on the QRPA calculations indicate a significant reduction in the deleptonization during the collapse phase. Conclusions: It is important to utilize microscopic theoretical models that are tested by experimental data to constrain and estimate Gamow-Teller strengths and derived electron-capture rates for nuclei near N=50 that are inputs for astrophysical simulations of core-collapse supernovae and their multimessenger signals, such as the emission of neutrinos and gravitational waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Effect of Hydrophobic Core Topology and Composition on the Structure and Kinetics of Star Polymers: A Molecular Dynamics Study

We present a molecular dynamics study of the effect of core chemistry on star polymer structural and kinetic properties. This work serves to validate the choice of a model adamantane core used in previous simulations to represent larger star polymeric systems in an aqueous environment, as well as to explore how the choice of size and core chemistry using a dendrimer or nanogel core may affect these polymeric nanoparticle systems, particularly with respect to thermosensitivity and solvation properties that are relevant for applications in drug loading and delivery.

Carr, Amber C.↗
DOE OSTI2011

Simulation-based Performance Analysis and Tuning for a Two-level Directly Connected System

Hardware and software co-design is becoming increasingly important due to complexities in supercomputing architectures. Simulating applications before there is access to the real hardware can assist machine architects in making better design decisions that can optimize application performance. At the same time, the application and runtime can be optimized and tuned beforehand. BigSim is a simulation-based performance prediction framework designed for these purposes. It can be used to perform packet-level network simulations of parallel applications using existing parallel machines. In this paper, we demonstrate the utility of BigSim in analyzing and optimizing parallel application performance for future systems based on the PERCS network. We present simulation studies using benchmarks and real applications expected to run on future supercomputers. Future petascale systems will have more than 100,000 cores, and we present simulations at that scale.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Utilizing many-core accelerators for halo and center finding within a cosmology simulation

Efficiently finding and computing statistics about “halos” (regions of high density) are essential analysis steps for N-body cosmology simulations. However, in state-of-the-art simulation codes, these analysis operators do not currently take advantage of the shared-memory data-parallelism available on multi-core and many-core architectures. The Hybrid / Hardware Accelerated Cosmology Code (HACC) is designed as an MPI+X code, but the analysis operators are parallelized only among MPI ranks, because of the difficulty in porting different X implementations (e.g., OpenMP, CUDA) across all architectures on which it is run. In this paper, we present portable data-parallel algorithms for several variations of halo finding and halo center finding algorithms. These are implemented with the PISTON component of the VTK-m framework, which uses Nvidia's Thrust library to construct data-parallel algorithms that allow a single implementation to be compiled to multiple backends to target a variety of multi-core and many-core architectures. Finally, we compare the performance of our halo and center finding algorithms against the original HACC implementations on the Moonlight, Stampede, and Titan supercomputers. The portability of Thrust allowed the same code to run efficiently on each of these architectures. On Titan, the performance improvements using our code have enabled halo analysis to be performed on a very large data set (81923 particles across 16,384 nodes of Titan) for which analysis using only the existing CPU algorithms was not feasible.

Sewell, Christopher↗
DOE OSTI2023

Strong Dark Matter Self-interactions Diversify Halo Populations within and surrounding the Milky Way

Abstract We perform a high-resolution cosmological zoom-in simulation of a Milky Way (MW)–like system, which includes a realistic Large Magellanic Cloud analog, using a large differential elastic dark matter self-interaction cross section that reaches ≈100 cm 2 g −1 at relative velocities of ≈10 km s −1 , motivated by the diverse and orbitally dependent central densities of dwarf galaxies within and surrounding the MW. We explore the effects of dark matter self-interactions on satellite, splashback, and isolated halos through their abundance, central densities, maximum circular velocities, orbital parameters, and correlations between these variables. We use an effective constant cross section model to analytically predict the stages of our simulated halos’ gravothermal evolution, demonstrating that deviations from the collisionless R max – V max relation can be used to select deeply core-collapsed halos, where V max is a halo’s maximum circular velocity, and R max is the radius at which it occurs. We predict that a sizable fraction (≈20%) of subhalos with masses down to ≈10 8 M ⊙ is deeply core collapsed in our SIDM model. Core-collapsed systems form ≈10% of the isolated halo population down to the same mass; these isolated, core-collapsed halos would host faint dwarf field galaxies with extremely steep central density profiles. Finally, most halos with masses above ≈10 9 M ⊙ are core-forming in our simulation. Our study thus demonstrates how self-interactions diversify halo populations in an environmentally dependent fashion within and surrounding MW-mass hosts, providing a compelling avenue to address the diverse dark matter distributions of observed dwarf galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Interim MELCOR Simulation of the Fukushima Daiichi Unit 2 Accident Reactor Core Isolation Cooling Operation

Data, a brief description of key boundary conditions, and results of Sandia National Laboratories’ ongoing MELCOR analysis of the Fukushima Unit 2 accident are given for the reactor core isolation cooling (RCIC) system. Important assumptions and related boundary conditions in the current analysis additional to or different than what was assumed/imposed in the work of SAND2012-6173 are identified. This work is for the U.S. Department of Energy’s Nuclear Energy University Programs fiscal year 2014 Reactor Safety Technologies Research and Development Program RC-7: RCIC Performance under Severe Accident Conditions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Core-corona hadronization model and muon puzzle in extensive air showers

The discrepancy between models and measurements concerning the muon content in air showers produced by ultra-high energy cosmic rays, the so-called muon puzzle, remains unsolved. String fragmentation models currently used in hadronic interactions fail to reproduce the observations, while recent measurements at the LHC hint towards the existence of production mechanisms, such as collective statistical hadronization, which lead to an increase in the muon production after a hadronic interaction. The core-corona model of heavy ions implements both production mechanisms, where the large-density region of interaction (core) hadronizes statistically, while the low-density (corona) hadronizes following string fragmentation. In this contribution, we present a modification of the CONEX framework to implement the core-corona model in air shower simulations. We demonstrate a significant impact of the core effect, as observed at the LHC, on the muon content in air showers produced by ultra-high energy cosmic rays.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1997

The Universe at Moderate Redshift

The report covers the work done in the past year and a wide range of fields including properties of clusters of galaxies; topological properties of galaxy distributions in terms of galaxy types; patterns of gravitational nonlinear clustering process; development of a ray tracing algorithm to study the gravitational lensing phenomenon by galaxies, clusters and large-scale structure, one of whose applications being the effects of weak gravitational lensing by large-scale structure on the determination of q(0); the origin of magnetic fields on the galactic and cluster scales; the topological properties of Ly(alpha) clouds the Ly(alpha) optical depth distribution; clustering properties of Ly(alpha) clouds; and a determination (lower bound) of Omega(b) based on the observed Ly(alpha) forest flux distribution. In the coming year, we plan to continue the investigation of Ly(alpha) clouds using larger dynamic range (about a factor of two) and better simulations (with more input physics included) than what we have now. We will study the properties of galaxies on 1 - 100h(sup -1) Mpc scales using our state-of-the-art large scale galaxy formation simulations of various cosmological models, which will have a resolution about a factor of 5 (in each dimension) better than our current, best simulations. We will plan to study the properties of X-ray clusters using unprecedented, very high dynamic range (20,000) simulations which will enable us to resolve the cores of clusters while keeping the simulation volume sufficiently large to ensure a statistically fair sample of the objects of interest. The details of the last year's works are now described.

Cen, Renyue↗
NASA NTRS2011

Computational Experiments for Science and Engineering Education

How to integrate simulation-based engineering and science (SBES) into the science curriculum smoothly is a challenging question. For the importance of SBES to be appreciated, the core value of simulations-that they help people understand natural phenomena and solve engineering problems-must be taught. A strategy to achieve this goal is to introduce computational experiments to the science curriculum to replace or supplement textbook illustrations and exercises and to complement or frame hands-on or wet lab experiments. In this way, students will have an opportunity to learn about SBES without compromising other learning goals required by the standards and teachers will welcome these tools as they strengthen what they are already teaching. This paper demonstrates this idea using a number of examples in physics, chemistry, and engineering. These exemplary computational experiments show that it is possible to create a curriculum that is both deeper and wider.

Xie, Charles↗
NASA NTRS2016

Post-Buckling Analysis of Curved Honeycomb Sandwich Panels Containing Interfacial Disbonds

A numerical study on the effect of facesheet-core disbonds on the post-buckling response of curved honeycomb sandwich panels is presented herein. This work was conducted as part of the development of a damage tolerance plan for the next-generation Space Launch System heavy lift launch vehicle payload fairing. As such, the study utilized full-scale fairing barrel segments as the structure of interest. The panels were composed of carbon fiber reinforced polymer facesheets and aluminum honeycomb core. The panels were analyzed numerically using the finite element method incorporating geometric nonlinearity. In a predetermined circular region, facesheet and core nodes were detached to simulate a disbond, between the outer mold line facesheet and honeycomb core, induced via low-speed impact. Surface-to-surface contact in the disbonded region was invoked to prevent interpenetration of the facesheet and core elements and obtain realistic stresses in the core. The diameter of this disbonded region was varied and the effect of the size of the disbond on the post-buckling response was observed. Significant changes in the slope of the edge load-deflection response were used to determine the onset of global buckling and corresponding buckling load. Finally, several studies were conducted to determine the sensitivity of the numerical predictions to refinement in the finite element mesh.

elastic buckling↗
DOE OSTI2022

Data assimilation for burnup distribution of PWR with three-dimensional variational algorithm and artificial neutral network

In this paper, a data-assimilation method has been proposed and applied for the burnup distribution of PWR. The burnup distribution is significant to the safety and economy of the reactor, as it is essential for the fuel-reloading design and optimization. Due to the burnup distribution cannot be measured directly during the reactor operation, the numerical simulation is widely applied to determine the burnup distribution. However, there is a deviation between the numerical simulation and the actual core due to some unavoidable factors, such as component manufacturing deviation, uneven flow distribution and so on. These differences would induce the errors to the simulation values of power distributions and hence to the burnup distribution. To address this problem, a data-assimilation method for the burnup distribution has been proposed with the application of power-distribution measurements. In our research, the three-dimensional variational (3DVAR) algorithm was applied for burnup-distribution calibration and the artificial neutral network (ANN) was applied to establish the relation between power distribution and corresponding burnup distribution. As engineering verification, the proposed data-assimilation method has been applied to the CNP1000 PWR operated in China. The numerical results indicated that the burnup-distribution errors can be reduced notably, as the maximum value of relative errors for power distribution can be reduced from 5.35% to 3.96%. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2011

Dynamic Simulation of Human Thermoregulation and Heat Transfer for Spaceflight Applications

Models of human thermoregulation and heat transfer date from the early 1970s and have been developed for applications ranging from evaluating thermal comfort in spacecraft and aircraft cabin environments to predicting heat stress during EVAs. Most lumped or compartment models represent the body as an assemblage cylindrical and spherical elements which may be subdivided into layers to describe tissue heterogeneity. Many existing models are of limited usefulness in asymmetric thermal environments, such as may be encountered during an EVA. Conventional whole-body clothing models also limit the ability to describe local surface thermal and evaporation effects in sufficient detail. A further limitation is that models based on a standard man model are not readily scalable to represent large or small subjects. This work describes development of a new human thermal model derived from the 41-node man model. Each segment is divided into four concentric, constant thickness cylinders made up of a central core surrounded by muscle, fat, and skin, respectively. These cylinders are connected by the flow of blood from a central blood pool to each part. The central blood pool is updated at each time step, based on a whole-body energy balance. Results show the model simulates core and surface temperature histories, sweat evaporation and metabolic rates which generally are consistent with controlled exposures of human subjects. Scaling rules are developed to enable simulation of small and large subjects (5th percentile and 95th percentile). Future refinements will include a clothing model that addresses local surface insulation and permeation effects and developing control equations to describe thermoregulatory effects such as may occur with prolonged weightlessness or with aging.

Miller, Thomas R.↗
DOE OSTI2016

Structure, dynamics, and deuterium fractionation of massive pre-stellar cores

High levels of deuterium fraction in N{sub 2}H{sup +} are observed in some pre-stellar cores. Single-zone chemical models find that the timescale required to reach observed values (D{sub frac}{sup N{sub 2}H{sup +}}≡N{sub 2}D{sup +}/N{sub 2}H{sup +}≳0.1) is longer than the free-fall time, possibly 10 times longer. Here, we explore the deuteration of turbulent, magnetized cores with 3D magnetohydrodynamics simulations. We use an approximate chemical model to follow the growth in abundances of N{sub 2}H{sup +} and N{sub 2}D{sup +}. We then examine the dynamics of the core using each tracer for comparison to observations. We find that the velocity dispersion of the core as traced by N{sub 2}D{sup +} appears slightly sub-virial compared to predictions of the Turbulent Core Model of McKee and Tan, except at late times just before the onset of protostar formation. By varying the initial mass surface density, the magnetic energy, the chemical age, and the ortho-to-para ratio of H{sub 2}, we also determine the physical and temporal properties required for high deuteration. We find that low initial ortho-to-para ratios (≲0.01) and/or multiple free-fall times (≳3) of prior chemical evolution are necessary to reach the observed values of deuterium fraction in pre-stellar cores.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Induced core formation time in subcritical magnetic clouds by large-scale trans-Alfvénic flows

We clarify the mechanism of accelerated core formation by large-scale nonlinear flows in subcritical magnetic clouds by finding a semi-analytical formula for the core formation time and describing the physical processes that lead to them. Recent numerical simulations show that nonlinear flows induce rapid ambipolar diffusion that leads to localized supercritical regions that can collapse. Here, we employ non-ideal magnetohydrodynamic simulations including ambipolar diffusion for gravitationally stratified sheets threaded by vertical magnetic fields. One of the horizontal dimensions is eliminated, resulting in a simpler two-dimensional simulation that can clarify the basic process of accelerated core formation. A parameter study of simulations shows that the core formation time is inversely proportional to the square of the flow speed when the flow speed is greater than the Alfvén speed. We find a semi-analytical formula that explains this numerical result. The formula also predicts that the core formation time is about three times shorter than that with no turbulence, when the turbulent speed is comparable to the Alfvén speed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Demonstration of Load-Follow Operations with Standalone BISON through VERA

As part of the Consortium for Advanced Simulation of Light Water Reactors (CASL), the Virtual Environment for Reactor Applications (VERA) is being developed to provide high-fidelity multiphysics simulations of nuclear reactor cores. VERA contains a number of single physics packages that are coupled together. MPACT is the primary deterministic neutron transport solver, CTF provides thermal hydraulics solutions, and BISON - which is built off of the MOOSE framework - is used for fuel performance simulations. Over the past year, a standalone BISON capability has been developed using a one-way coupling between MPACT/CTF coupling calculations and BISON. This has allowed CASL contributors to gain insight into the fuel performance characteristics while a more tightly coupled methodology between MPACT, CTF, and BISON has been under development. The standalone capability has been used to simulate Watts Bar Unit 1 (WBN1) Cycles 1-3 and to screen for quantities of interest to pellet-clad interaction in Cycles 6-7. The simulation and screening results are presented here. CASL recently partnered with Exelon and the University of Illinois at Urbana Champaign (UIUC) to begin simulating some of Exelon's plants, with particular focus on load-follow operations. In such operations, the power output from the plant will vary based on the demand or anticipated demand at the time. For example, demand is much lower late at night than during the day when people are awake, so it operating the plant at lower power in the evenings could be advantageous. Historically, nuclear power generation has been among the cheapest sources of energy, so most plants in the United States have been operated at full power except for outages and various small scale events. However, in the current economic climate, priorities have changed. A number of questions arise when load-follow operations are considered, particularly with respect to VERA's tools and capabilities. Can the coupled neutronics and thermal hydraulics simulations with MPACT and CTF perform well without substantial convergence issues? Can BISON handle the somewhat rapid power changes present with load-follow? Will clad hoop stresses be alarmingly high as a result of these power changes? This work addresses these questions and demonstrates operations to assess feasibility and robustness before Exelon and UIUC proceed with their simulations using VERA. To demonstrate load-follow capability, two test problems were developed. The first is a single rod that experiences nominal operations during Cycle 1, and then it simulates 18 months of representative load-follow operations in Cycle 2. This single-rod case was used to test the robustness of MPACT, CTF, and BISON on a typical rod. The second problem is a quarter-core model based on WBN1, Cycle 3, with a single month of load-follow operations inserted into the middle of the cycle. This provides better understanding of the stresses to the clad during the numerous power changes. Sections of this article are adapted from the original, full technical report on this capability.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

2D properties of core turbulence on DIII-D and comparison to gyrokinetic simulations

Quantitative 2D characteristics of localized density fluctuations are presented over the range of 0.3 < r/a < 0.9 in L-mode plasmas on DIII-D [J. L. Luxon, Nucl. Fusion 42, 614 (2002)]. Broadband density fluctuations increase in amplitude from (n) over tilde/n < 0.5% in the deep core to (n) over tilde/n similar to 2.5% near the outer region. The observed Doppler-shift due to the E x B velocity matches well with the measured turbulence group and phase velocities (in toroidally rotating neutral beam heated plasmas). Turbulence decorrelation rates are found to be similar to 200 kHz at the edge and to decrease toward the core (0.45 < r/a < 0.9) where they approach the E x B shearing rate (similar to 50 kHz). Radial and poloidal correlation lengths are found to scale with the ion gyroradius and exhibit an asymmetric poloidally elongated eddy structure. The ensemble-averaged turbulent eddy structure changes its tilt with respect to the radial-poloidal coordinates in the core, consistent with an E x B shear mechanism. The 2D spatial correlation and wavenumber spectra [S(k(r); k(theta))] are presented and compared to nonlinear flux-tube GYRO simulations at two radii, r/a = 0.5 and r/a = 0.75, showing reasonable overall agreement, but the GYRO spectrum exhibits a peak at finite kr for r/a = 0.75 that is not observed experimentally; E x B shear may cause this discrepancy. (C) 2012 American Institute of Physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗