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 1,405 records · Page 78

Evolution of Photon and Particle Spectra in Compact, Luminous Objects

Physical conditions in the radiating plasma in the cores of radio-strong quasars and active galactic nuclei cannot be derived from observations until the effects of relativistic aberration are understood. This requires determining both the bulk flow speeds and any wave or signal speed in the parsec-scale nuclear jets. In this project we studied several aspects of such waves. We considered constraints on jet deceleration by mass pickup, and found that bolometric luminosities of the active nuclei cannot constrain core jet speeds usefully. We also simulated observations of ballistic, helical trajectories and helical waves moving directly outwards along the jet. We found that ballistic trajectories are not allowed by the data; the helical features seen are very likely to be helical waves. We believe these are waves propagating in the jet plasma. To this end, we studied waves propagating in relativistic pair plasma jets. In particular, we undertook a program whose goal was to determine the nature of waves which can propagate in relativistic pair plasmas, and how such waves propagating in streaming jet plasma would be observed by an external observer. We developed the possibility of using pulsars as test cases for our models; this takes advantage of new technology in pulsar observations, and the similarity of the physical conditions in the pulsar magnetosphere to the dense, relativistic pair plasmas which exist in radio-strong quasars.

Eilek, Jean A.↗

Quasar Host Galaxies/Neptune Rotation/Galaxy Building Blocks/Hubble Deep Field/Saturn Storm

Computerized animations simulate a quasar erupting in the core of a normal spiral galaxy, the collision of two interacting galaxies, and the evolution of the universe. Hubble Space Telescope (HST) images show six quasars' host galaxies (including spirals, ellipticals, and colliding galaxies) and six clumps of galaxies approximately 11 billion light years away. A false color time lapse movie of Neptune displays the planet's 16-hour rotation, and the evolution of a storm on Saturn is seen though a video of the planet's rotation. A zoom sequence starts with a ground-based image of the constellation Ursa major and ends with the Hubble Deep Field through progressively narrower and deeper views.

Source record↗

COR1 Engineering Test Unit Measurements at the NCAR/HAO Vacuum Tunnel Facility, October-November 2002

The Engineering Test Unit (ETU) of COR1 was made in two configurations. The first configuration, ETU-1, was for vibration testing, while the second, ETU-2, was for optical testing. This is a report on the optical testing performed on ETU-2 at the NCAR/HAO Vacuum Tunnel Facility during the months of October and November, 2002. This was the same facility used to test the two previous breadboard models. In both configurations, the first two tube sections were complete, with all optical elements aligned. The vibration model ETU-1 had the remaining tube sections attached, with mass models for the remaining optics, for the various mechanisms, and for the focal plane assembly. It was then converted into the optical model ETU-2 by removing tube sections 3 to 5, and mounting the remaining optics on commercial mounts. (The bandpass filter was also installed into tube 2, which had been replaced in ETU-1 by a mass model, so that pre- and post-vibration optical measurements could be made.) Doublet 2 was installed in a Newport LP-2 carrier, and aligned to the other optics in the first two tube sections. The LP-2 adjustment screws were then uralened so that the alignment could be maintained during shipping. Because neither the flight polarizer nor Hollow Core Motor were available, they were simulated by a commercial polarizer and rotational mount, both from Oriel corporation. The Oriel rotational stage was not designed for vacuum use, but it was determined after consultation with the company, and lab testing, that the stage could be used in the moderate vacuum conditions at the NCAR/HAO facility. The shutter and focal plane assembly were simulated with the same camera used for the previous two breadboard tests. The focal plane mask was simulated with a plane of BK7 glass with a mask glued on, using the same procedure as for the Lyot spot on Doublet 1, and mounted in an adjustable LP-2 carrier. Two masks were made, one made to the precise specifications of the optical design, the other slightly bigger to make alignment easier.

Thompson, William↗

Ndarts

Ndarts software provides algorithms for computing quantities associated with the dynamics of articulated, rigid-link, multibody systems. It is designed as a general-purpose dynamics library that can be used for the modeling of robotic platforms, space vehicles, molecular dynamics, and other such applications. The architecture and algorithms in Ndarts are based on the Spatial Operator Algebra (SOA) theory for computational multibody and robot dynamics developed at JPL. It uses minimal, internal coordinate models. The algorithms are low-order, recursive scatter/ gather algorithms. In comparison with the earlier Darts++ software, this version has a more general and cleaner design needed to support a larger class of computational dynamics needs. It includes a frames infrastructure, allows algorithms to operate on subgraphs of the system, and implements lazy and deferred computation for better efficiency. Dynamics modeling modules such as Ndarts are core building blocks of control and simulation software for space, robotic, mechanism, bio-molecular, and material systems modeling.

Jain, Abhinandan↗

Time-Accurate Local Time Stepping and High-Order Time CESE Methods for Multi-Dimensional Flows Using Unstructured Meshes

With the wide availability of affordable multiple-core parallel supercomputers, next generation numerical simulations of flow physics are being focused on unsteady computations for problems involving multiple time scales and multiple physics. These simulations require higher solution accuracy than most algorithms and computational fluid dynamics codes currently available. This paper focuses on the developmental effort for high-fidelity multi-dimensional, unstructured-mesh flow solvers using the space-time conservation element, solution element (CESE) framework. Two approaches have been investigated in this research in order to provide high-accuracy, cross-cutting numerical simulations for a variety of flow regimes: 1) time-accurate local time stepping and 2) highorder CESE method. The first approach utilizes consistent numerical formulations in the space-time flux integration to preserve temporal conservation across the cells with different marching time steps. Such approach relieves the stringent time step constraint associated with the smallest time step in the computational domain while preserving temporal accuracy for all the cells. For flows involving multiple scales, both numerical accuracy and efficiency can be significantly enhanced. The second approach extends the current CESE solver to higher-order accuracy. Unlike other existing explicit high-order methods for unstructured meshes, the CESE framework maintains a CFL condition of one for arbitrarily high-order formulations while retaining the same compact stencil as its second-order counterpart. For large-scale unsteady computations, this feature substantially enhances numerical efficiency. Numerical formulations and validations using benchmark problems are discussed in this paper along with realistic examples.

Chang, Chau-Lyan↗

Reverse Computation for Rollback-based Fault Tolerance in Large Parallel Systems

Reverse computation is presented here as an important future direction in addressing the challenge of fault tolerant execution on very large cluster platforms for parallel computing. As the scale of parallel jobs increases, traditional checkpointing approaches suffer scalability problems ranging from computational slowdowns to high congestion at the persistent stores for checkpoints. Reverse computation can overcome such problems and is also better suited for parallel computing on newer architectures with smaller, cheaper or energy-efficient memories and file systems. Initial evidence for the feasibility of reverse computation in large systems is presented with detailed performance data from a particle simulation scaling to 65,536 processor cores and 950 accelerators (GPUs). Reverse computation is observed to deliver very large gains relative to checkpointing schemes when nodes rely on their host processors/memory to tolerate faults at their accelerators. A comparison between reverse computation and checkpointing with measurements such as cache miss ratios, TLB misses and memory usage indicates that reverse computation is hard to ignore as a future alternative to be pursued in emerging architectures.

Perumalla, Kalyan S↗

Prediction of number of breached rods following a LBLOCA of Candu plants using a BEPU approach

Radioactive doses following design basis accidents (DBA) have been important safety criteria of Candu nuclear power plant and they have been predicted in terms of the number of breached fuel rods. To support the licensing review on this concern, an analysis of LBLOCA has been conducted by using the BEPU method of KINS, KINS-REM. Number of Breached Rods (NBR) following a LBLOCA was predicted at 95 percentile probabilistic upper level in 95 percentile confidence level. Peak Cladding Temperatures (PCT) of the 84 bundles in the core pass 4 were calculated from the 124 MARS code runs in which the uncertainties of 10 major parameters including fuel thermal conductivity and break flow model were implemented. The fuel rod breaching criteria, PCT>1477 K, was used to determine the NBR 95/95. From the calculation, the predicted NBR 95/95 was 1591 rods and the calculated maximum NBR was lower than 2000 rods. Through the further improvements in feedback of the channel power behavior to thermalhydraulic calculation and in channel group modeling, NBR in more reliable level can be expected. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Progress on ATF development of SiC for LWR

Since 2011, Toshiba has been participating in a study whose objective is the commercial utilization of SiC instead of Zircaloy as LWR fuel material. During the period (Phase-1) from 2012 to 2015, feasibility studies on material design and fabrication were performed. This paper reports the principal results of a comprehensive feasibility study on SiC involving experiments and trial fabrication in terms of material properties, core and fuel design, and analytical simulation for safety analysis. The principal results are summarized as follows. First, the fabrication trials for the reduced scales of channel box and cladding have been successfully performed. Secondly, with regard to fuel rod performance analysis, hoop stress during normal operation may be a potential issue in fuel mechanical design as it is relatively high. Thirdly, with regard to safety analysis, the hydrogen generation of SiC fuel tends to decrease significantly compared with that of Zircaloy fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Time-resolved measurement of subchannel void distribution for boiling two-phase flow in 5 x 5 rod bundle at high-pressure and high-temperature

A subchannel void sensor (SCVS) was improved to acquire time-resolve boiling two-phase flow in a 5-by-5 heated rod bundle at high-pressure and high-temperature. The SCVS consists of 6-by-6 wire electrodes and 5-by-5 rod electrodes. The SCVS can acquire a time series data of local void fraction at 32 points of central subchannel region in addition to 100 points of rod surface regions between rods. The devised sensors are installed in the heated rod bundle at eight height levels to acquire two-phase flow dynamics and its development. The axial and radial power profile of the heated rod bundle are uniform, and eight pairs of sheath thermocouples are embedded on the heated rod to monitor its surface temperature distribution. The paper addresses the axial and cross-sectional distributions of void fraction at various pressure and rod bundle power conditions, which simulates an accidental condition in reactor core and in a spent fuel pool of a boiling water reactor (BWR). The experimental data are suitable to validate the computational multi-fluid dynamics (CMFD) and subchannel codes. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fuel cycle depletion validation and code-to-code verification studies for High Flux Isotope Reactor highly and low-enriched uranium fuel designs

Here, this paper documents fuel cycle depletion validation and code-to-code verification studies for the High Flux Isotope Reactor (HFIR) highly enriched uranium (HEU) and proposed low-enriched uranium (LEU) fuel designs. In support of HFIR’s world-leading performance, transport and depletion simulations are performed to ensure safe operations, design and qualify irradiation experiments, enhance core components and irradiation facilities, and design and characterize LEU fuel designs. Identifying well-validated, computationally efficient codes is required for the success of these efforts. The HFIR Controller, Shift, and VESTA codes were deployed to simulate HEU uranium–oxide dispersion fuel cycles at 85, 95, and 100 MW operations, as well as LEU fuel cycles operating at 95 MW with uranium–silicide dispersion and uranium–molybdenum monolithic alloy fuel forms. Excellent agreement between the codes and with experimentally obtained 235 U enrichment distributions provides increased confidence in the ability of these codes to model and simulate HFIR’s unique core design.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Numerical challenges in modeling gravothermal collapse in Self-Interacting Dark Matter halos

When dark matter has a large cross section for self scattering, halos can undergo a process known as gravothermal core collapse, where the inner core rapidly increases in density and temperature. To date, several methods have been used to implement Self-Interacting Dark Matter (SIDM) in N-body codes, but there has been no systematic study of these different methods or their accuracy in the core-collapse phase. In this paper, we compare three different numerical implementations of SIDM, including the standard methods from the GIZMO and Arepo codes, by simulating idealized dwarf halos undergoing significant dark matter self interactions (σ/m = 50 cm 2 /g). When simulating these halos, we also vary the mass resolution, time-stepping criteria, and gravitational force-softening scheme. The various SIDM methods lead to distinct differences in a halo's evolution during the core-collapse phase, as each results in spurious scattering rate differences and energy gains/losses. The use of adaptive force softening for gravity can lead to numerical heating that artificially accelerates core collapse, while an insufficiently small simulation time step can cause core evolution to stall or completely reverse. Additionally, particle numbers must be large enough to ensure that the simulated halos are not sensitive to noise in the initial conditions. Even for the highest-resolution simulations tested in this study (10 6 particles per halo), we find that variations of order 10% in collapse time are still present. Here, the results of this work underscore the sensitivity of SIDM modeling on the choice of numerical implementation and motivate a careful study of how these results generalize to halos in a cosmological context.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Prometheus: Scalable and Accurate Emulation of Task-Based Applications on Many-Core Systems.

Modeling the performance of non-deterministic parallel applications on future many-core systems requires the development of novel simulation and emulation techniques and tools. We present “Prometheus”, a fast, accurate and modular emulation framework for task-based applications. By raising the level of abstraction and focusing on runtime synchronization, Prometheus can accurately predict applications’ performance on very large many-core systems. We validate our emulation framework against two real platforms (AMD Interlagos and Intel MIC) and report error rates generally below 4%. We, then, evaluate Prometheus’ performance and scalability: our results show that Prometheus can emulate a task-based application on a system with 512K cores in 11.5 hours. We present two test cases that show how Prometheus can be used to study the performance and behavior of systems that present some of the characteristics expected from exascale supercomputer nodes, such as active power management and processors with a high number of cores but reduced cache per core.

exascale systems↗

Main results of Phenix final core physics tests

The French sodium cooled fast reactor Phenix was shut down in 2009 after 35 years of operation. Before decommissioning a final set of tests was performed. This paper focuses on the following core physics tests: - measurement of control rod reactivity worth by several methods (sub-critical, critical, rod-drop methods), - control rod shifting during a full power operation state, - measurement of individual subassembly reactivity worth (fresh and burned fuel and fertile, sodium hole), - simulation of a gas bubble crossing the core. The control rod measurement test has shown some discrepancies between the different measurement methods, mainly for the rod bank worth. Considering a macroscopic parameter (the reactivity loss estimation), NSMM method seems to produce the better results. The control rod shifting test has highlighted the impact of the spatial effects generated by control rods movements over the power map. These spatial effects, known as shadowing effects, can modify up to 15% the individual control rod worth. The tests concerning individual subassembly worth and gas bubble have permit to evaluate the impact of local perturbations on the reactivity. All these tests can be considered as successful and the ability of the European neutronic code for fast reactors, ERANOS 2.2, to reproduce complex and atypical configuration has been demonstrated by good agreement between measured and computed values. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The impact of feedback during massive star formation by core accretion

We study feedback during massive star formation using semi-analytic methods, considering the effects of disk winds, radiation pressure, photoevaporation, and stellar winds, while following protostellar evolution in collapsing massive gas cores. We find that disk winds are the dominant feedback mechanism setting star formation efficiencies (SFEs) from initial cores of ∼0.3–0.5. However, radiation pressure is also significant to widen the outflow cavity causing reductions of SFE compared to the disk-wind only case, especially for >100 M{sub ⊙} star formation at clump mass surface densities Σ{sub cl}≲0.3 g cm{sup −2}. Photoevaporation is of relatively minor importance due to dust attenuation of ionizing photons. Stellar winds have even smaller effects during the accretion stage. For core masses M{sub c}≃10–1000 M{sub ⊙} and Σ{sub cl}≃0.1–3 g cm{sup −2}, we find the overall SFE to be ε-bar {sub ∗f}=0.31(R{sub c}/0.1pc){sup −0.39}, potentially a useful sub-grid star formation model in simulations that can resolve pre-stellar core radii, R{sub c}=0.057(M{sub c}/60M{sub ⊙}){sup 1/2}(Σ{sub cl}/gcm{sup −2}){sup −1/2} pc. The decline of SFE with M {sub c} is gradual with no evidence for a maximum stellar-mass set by feedback processes up to stellar masses of m{sub ∗}∼300 M{sub ⊙}. We thus conclude that the observed truncation of the high-mass end of the IMF is shaped mostly by the pre-stellar core mass function or internal stellar processes. To form massive stars with the observed maximum masses of ∼150–300M{sub ⊙}, initial core masses need to be ≳500–1000 M{sub ⊙}. We also apply our feedback model to zero-metallicity primordial star formation, showing that, in the absence of dust, photoevaporation staunches accretion at ∼50 M{sub ⊙}. Our model implies radiative feedback is most significant at metallicities ∼10{sup −2}Z{sub ⊙}, since both radiation pressure and photoevaporation are effective in this regime.

79 ASTRONOMY AND ASTROPHYSICS↗

Synergistic cross-scale coupling of turbulence in a tokamak plasma

For the first time, nonlinear gyrokinetic simulations spanning both the ion and electron spatio-temporal scales have been performed with realistic electron mass ratio ((m D /m e ) 1/2 = 60.0), realistic geometry, and all experimental inputs, demonstrating the coexistence and synergy of ion ( ( k θ ρ s ∼ O ( 1.0 ) ) ) and electron-scale ( ( k θ ρ e ∼ O ( 1.0 ) ) ) turbulence in the core of a tokamak plasma. All multi-scale simulations utilized the GYRO code [J. Candy and R. E. Waltz, J. Comput. Phys. 186, 545 (2003)] to study the coupling of ion and electron-scale turbulence in the core (r/a = 0.6) of an Alcator C-Mod L-mode discharge shown previously to exhibit an under-prediction of the electron heat flux when using simulations only including ion-scale turbulence. Electron-scale turbulence is found to play a dominant role in setting the electron heat flux level and radially elongated (k r ≪ k θ ) “streamers” are found to coexist with ion-scale eddies in experimental plasma conditions. Inclusion of electron-scale turbulence in these simulations is found to increase both ion and electron heat flux levels by enhancing the transport at the ion-scale while also driving electron heat flux at sub-ρ i scales. The combined increases in the low and high-k driven electron heat flux may explain previously observed discrepancies between simulated and experimental electron heat fluxes and indicates a complex interaction of short and long wavelength turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Black Hole Supernovae, Their Equation of State Dependence, and Ejecta Composition

Abstract Recent literature on core-collapse supernovae suggests that a black hole (BH) can form within ∼1 s of shock revival, while still culminating in a successful supernova. We refer to these as BH supernovae, as they are distinct from other BH formation channels in both timescale and impact on the explosion. We simulate these events self-consistently from core collapse until 20–50 days after collapse using three axisymmetric models of a 60 M ⊙ zero-age main-sequence progenitor star and investigate how the composition of the ejecta is impacted by the BH formation. We employ Skyrme-type equations of state (EOSs) and vary the uncertain nucleonic effective mass, which affects the pressure inside the proto–neutron star through the thermal part of the EOS. This results in different BH formation times and explosion energies at BH formation, yielding final explosion energies between 0.06 and 0.72 × 10 51 erg with 21.8–23.3 M ⊙ of ejecta, of which 0–0.018 M ⊙ is 56 Ni. Compared to expectations from 1D simulations, we find more nuanced EOS dependences of the explosion dynamics, the mass of the BH remnant, and the elemental composition of the ejecta. We investigate why the explosions survive despite the massive overburden and link the shape of the diagnostic energy curve and character of the ejecta evolution to the progenitor structure.

Eggenberger Andersen, Oliver (ORCID:00000002966079↗

Warm-Core Intensification Through Horizontal Eddy Heat Transports into the Eye

A simulation of Hurricane Bob (1991) using the PSU/NCAR MM5 mesoscale model with a finest mesh spacing of 1.3 km is used to diagnose the heat budget of the hurricane. Heat budget terms, including latent and radiative heating, boundary layer forcing, and advection terms were output directly from the model for a 6-h period with 2-min frequency. Previous studies of warm core formation have emphasized the warming associated with gentle subsidence within the eye. The simulation of Hurricane Bob confirms subsidence warming as a major factor for eye warming, but also shows a significant contribution from horizontal advective terms. When averaged over the area of the eye, subsidence is found to strongly warm the mid-troposphere (2-9 km) while horizontal advection warms the mid to upper troposphere (5-13 km) with about equal magnitude. Partitioning of the horizontal advective terms into azimuthal mean and eddy components shows that the mean radial circulation does not, as expected, generally contribute to this warming, but that it is produced almost entirely by the horizontal eddy transport of heat into the eye. A further breakdown of the eddy components into azimuthal wave numbers 1, 2, and higher indicates that the warming is dominated by wave number 1 asymmetries, with smaller coming from higher wave numbers. Warming by horizontal eddy transport is consistent with idealized modeling of vortex Rossby waves and work is in progress to identify and clarify the role of vortex Rossby waves in warm-core intensification in both the full-physics model and idealized models.

Braun, Scott A.↗

The TAVERNS emulator: An Ada simulation of the space station data communications network and software development environment

The Space Station DMS (Data Management System) is the onboard component of the Space Station Information System (SSIS) that includes the computers, networks and software that support the various core and payload subsystems of the Space Station. TAVERNS (Test And Validation Environment for Remote Networked Systems) is a distributed approach for development and validation of application software for Space Station. The TAVERNS concept assumes that the different subsystems will be developed by different contractors who may be geographically separated. The TAVERNS Emulator is an Ada simulation of a TAVERNS on the ASD VAX. The software services described in the DMS Test Bed User's Manual are being emulated on the VAX together with simulations of some of the core subsystems and a simulation of the DCN. The TAVERNS Emulator will be accessible remotely from any VAX that can communicate with the ASD VAX.

Howes, Norman R.↗