Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical 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 OSTI2015

Physics-based multiscale coupling for full core nuclear reactor simulation

Numerical simulation of nuclear reactors is a key technology in the quest for improvements in efficiency, safety, and reliability of both existing and future reactor designs. Historically, simulation of an entire reactor was accomplished by linking together multiple existing codes that each simulated a subset of the relevant multiphysics phenomena. Recent advances in the MOOSE (Multiphysics Object Oriented Simulation Environment) framework have enabled a new approach: multiple domain-specific applications, all built on the same software framework, are efficiently linked to create a cohesive application. This is accomplished with a flexible coupling capability that allows for a variety of different data exchanges to occur simultaneously on high performance parallel computational hardware. Examples based on the KAIST-3A benchmark core, as well as a simplified Westinghouse AP-1000 configuration, demonstrate the power of this new framework for tackling—in a coupled, multiscale manner—crucial reactor phenomena such as CRUD-induced power shift and fuel shuffle. 2014 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-SA license

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Analysis of rapid increase in the plasma density during the ramp-up phase in a radio frequency negative ion source by large-scale particle simulation

Numerical simulations become useful for the developing RF-ICP (Radio Frequency Inductively Coupled Plasma) negative ion sources. We are developing and parallelizing a two-dimensional three velocity electromagnetic Particle-In-Cell code. The result shows rapid increase in the electron density during the density ramp-up phase. A radial electric field due to the space charge is produced with increase in the electron density and the electron transport in the radial direction is suppressed. As a result, electrons stay for a long period in the region where the inductive electric field is strong, and this leads efficient electron acceleration and a rapid increasing of the electron density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Numerical experiments on the detailed energy conversion and spectrum studies in a corona current sheet

In this paper, we study the energy conversion and spectra in a corona current sheet (CS) by 2.5 dimensional MHD numerical simulations. Numerical results show that many Petschek-like fine structures with slow-mode shocks mediated by plasmoid instabilities develop during the magnetic reconnection process. The termination shocks can also be formed above the primary magnetic island and at the head of secondary islands. These shocks play important roles in generating thermal energy in a corona CS. For a numerical simulation with initial conditions close to the solar corona environment, the ratio of the generated thermal energy to the total dissipated magnetic energy is around 1/5 before secondary islands appear. After secondary islands appear, the generated thermal energy starts to increase sharply and this ratio can reach a value of about 3/5. In an environment with a relatively lower plasma density and plasma β, the plasma can be heated to a much higher temperature. After secondary islands appear, the one-dimensional energy spectra along the CS do not behave as a simple power law and the spectrum index increases with the wave number. The average spectrum index for the magnetic energy spectrum along the CS is about 1.8. The two-dimensional spectra intuitively show that part of the high energy is cascaded to large kx and ky space after secondary islands appear. The plasmoid distribution function calculated from numerical simulations behaves as a power law closer to f(ψ)∼ψ{sup −1} in the intermediate ψ regime. By using η{sub eff}=v{sub inflow}⋅L, the effective magnetic diffusivity is estimated to be about 10{sup 11} ∼ 10{sup 12} m{sup 2} s{sup −1}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

OBSERVATIONAL IMPLICATIONS OF GAMMA-RAY BURST AFTERGLOW JET SIMULATIONS AND NUMERICAL LIGHT CURVE CALCULATIONS

We discuss jet dynamics for narrow and wide gamma-ray burst (GRB) afterglow jets and the observational implications of numerical simulations of relativistic jets in two dimensions. We confirm earlier numerical results that sideways expansion of relativistic jets during the bulk of the afterglow emission phase is logarithmic in time and find that this also applies to narrow jets with half opening angle of 0.05 rad. As a result, afterglow jets remain highly nonspherical until after they have become nonrelativistic. Although sideways expansion steepens the afterglow light curve after the jet break, the jet edges becoming visible dominates the jet break, which means that the jet break is sensitive to the observer angle even for narrow jets. Failure to take the observer angle into account can lead to an overestimation of the jet energy by up to a factor of four. This weakens the challenge posed to the magneter energy limit by extreme events such as GRB090926A. Late-time radio calorimetry based on a spherical nonrelativistic outflow model remains relevant when the observer is approximately on-axis and where differences of a few in flux level between the model and the simulation are acceptable. However, this does not imply sphericity of the outflow and therefore does not translate to high observer angles relevant to orphan afterglows. For more accurate calorimetry and in order to model significant late-time features such as the rise of the counterjet, detailed jet simulations remain indispensable.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Simulating a numerical UV completion of quartic Galileons

The Galileon theory is a prototypical effective field theory that incorporates the Vainshtein screening mechanism—a feature that arises in some extensions of general relativity, such as massive gravity. The Vainshtein effect requires that the theory contain higher order derivative interactions, which results in Galileons, and theories like them, failing to be technically well posed. While this is not a fundamental issue when the theory is correctly treated as an effective field theory, it nevertheless poses significant practical problems when numerically simulating this model. These problems can be tamed using a number of different approaches: introducing an active low-pass filter and/or constructing a UV completion at the level of the equations of motion, which controls the high momentum modes. These methods have been tested on cubic Galileon interactions, and have been shown to reproduce the correct low-energy behavior. Here we show how the numerical UV-completion method can be applied to quartic Galileon interactions, and present the first simulations of the quartic Galileon model using this technique. We demonstrate that our approach can probe physics in the regime of the effective field theory in which the quartic term dominates, while successfully reproducing the known results for cubic interactions. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2011

Process Simulation of Cold Pressing and Sintering of Armstrong CP-Ti Powders

A computational methodology is presented for the process simulation of cold pressing and sintering of Armstrong CP-Ti powders. Since the powder consolidation is governed by specific pressure-dependent constitutive equations, solution algorithms were developed for the ABAQUS user material subroutine, UMAT, for computing the plastic strain increments based on an implicit integration of the nonlinear yield function, flow rule, and hardening equations. Sintering was simulated using a model based on diffusional creep using the user subroutine CREEP. The initial mesh, stress, and density for the simulation of sintering were obtained from the results of the cold pressing simulation, minimizing the errors from decoupling the cold pressing and sintering simulations. Numerical simulation results are presented for the cold compaction followed by a sintering step of the Ti powders. The numerical simulation results for the relative density were compared to those measured from experiments before and after sintering, showing that the relative density can be accurately predicted. Notice: This manuscript has been authored by UT-Battelle, LLC, under Contract No. DE-AC05-00OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. ACKNOWLEDGEMENTS This research was sponsored by the U.S. DOE, and carried out at ORNL, under Contract DE-AC05-00OR22725 with UT-Battelle, LLC. This research was sponsored by the U.S. DOE, EERE Industrial Technology Program Office under CPS Agreement # 17881.

36 MATERIALS SCIENCE↗
DOE OSTI2010

Simulation of substrate degradation in composting of sewage sludge

To simulate the substrate degradation kinetics of the composting process, this paper develops a mathematical model with a first-order reaction assumption and heat/mass balance equations. A pilot-scale composting test with a mixture of sewage sludge and wheat straw was conducted in an insulated reactor. The BVS (biodegradable volatile solids) degradation process, matrix mass, MC (moisture content), DM (dry matter) and VS (volatile solid) were simulated numerically by the model and experimental data. The numerical simulation offered a method for simulating k (the first-order rate constant) and estimating k{sub 20} (the first-order rate constant at 20 {sup o}C). After comparison with experimental values, the relative error of the simulation value of the mass of the compost at maturity was 0.22%, MC 2.9%, DM 4.9% and VS 5.2%, which mean that the simulation is a good fit. The k of sewage sludge was simulated, and k{sub 20}, k{sub 20s} (first-order rate coefficient of slow fraction of BVS at 20 {sup o}C) of the sewage sludge were estimated as 0.082 and 0.015 d{sup -1}, respectively.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Monitoring of corrosion effects in pipes with multi-mode acoustic signals

Current commercial long-range guided wave inspection and monitoring technologies for pipes usually require establishment of certain specific acoustic modes for wave propagation and use a large number of sensors in a certain geometric arrangement (for example, uniformly-spaced arrays of sensors on a collar, placed around a pipe). This paper, in contrast, presents a long-range corrosion detection method in pipes based on multi-mode guided acoustic waves, generated and received using simple, single-element sensors. In this technique, available guided wave modes in any interrogated zone of the structure are automatically down-selected by the medium, and changes in energy distribution over the selected modes are analyzed through signal processing to monitor material loss or addition. Here, to investigate the feasibility of the proposed method, we have performed both experimental and numerical simulation studies on 50.8-mm (2-in)-diameter schedule-40 carbon steel pipes. In these studies, a linear chirp signal is injected in compression mode, perpendicular to the surface of the pipe, generating multi-mode, multi-directional guided waves that travel within the pipe wall. Corrosion effects in the pipes are simulated numerically and experimentally. To reduce errors in the measurement technique due to changes in the local or global temperature of the measured media, a universal temperature compensation algorithm is introduced. The experimental results show good agreement with simulation results, indicating that the proposed multi-mode wave generation and detection in a pipe could be used for long-range detection and monitoring of corrosion effects, such as pipe wall losses or solid material deposits.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Verification of numerical models for seismic fluid-structure interaction analysis of internal components in liquid-filled advanced reactors

Earthquake shaking of a liquid-filled advanced reactor induces fluid-structure interaction (FSI) between the reactor vessel (tank), its internal components, and the contained liquid. Seismic design, qualification, and risk assessment of such reactors must consider fluid-structure responses, their geometries and support conditions, and three-directional seismic inputs, all of which require numerical simulations. Numerical models used for analysis of safety-related nuclear equipment must be verified and validated. Here, this paper verifies numerical models of submerged components using analytical solutions. Seismic FSI analysis of the numerical models is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Prior analytical solutions are reworked in this paper, and calculation errors are identified and corrected. The solutions address frequencies of two concentric cylindrical pipes filled with liquid: the inner pipe is a submerged component, and the outer pipe contains the liquid. A unitless frequency coefficient is defined and presented for pipes of different materials, filled with different liquids, and with a range of dimensions for application to advanced reactors. The numerical models are verified here by comparing the lateral frequencies of submerged components with those calculated using the corrected analytical solutions. Recommendations for verification of numerical models of internal components in advanced reactors for seismic FSI analysis are provided. Although the unitless frequency coefficients and verification procedures are developed for application to advanced nuclear reactors, they are broadly applicable to FSI analysis of submerged components in liquid-filled vessels such as storage tanks, boilers, and steam generators.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Measurements of the populations of metastable and resonance levels in the plasma of an RF capacitive discharge in argon

The behavior of the populations of two metastable and two lower resonance levels of argon atoms in the plasma of an RF capacitive discharge was studied. The populations were measured by two methods: the method of emission self-absorption and the method based on measurements of the intensity ratios of spectral lines. It is shown that the populations of resonance levels increase with increasing power deposited in the discharge, whereas the populations of metastable levels is independent of the RF power. The distribution of the populations over energy levels is not equilibrium under these conditions. The population kinetics of argon atomic levels in the discharge plasma is simulated numerically. The distribution function of plasma electrons recovered from the measured populations of atomic levels and numerical simulations is found to be non-Maxwellian.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

History Matching and Performance Prediction of a Polymer Flood Pilot in Heavy Oil Reservoir on Alaska North Slope

The first-ever polymer flood pilot to enhance heavy oil recovery on Alaska North Slope (ANS) is ongoing. After more than 2.5 years of polymer injection, significant benefit has been observed from the decrease in water cut from 65% to less than 15% in the project producers. The primary objective of this study is to develop a robust history-matched reservoir simulation model capable of predicting future polymer flood performance. In this work, the reservoir simulation model has been developed based on the geological model and available reservoir and fluid data. In particular, four high transmissibility strips were introduced to connect the injector-producer well pairs, simulating short-circuiting flow behavior that can be explained by viscous fingering and reproducing the water cut history. The strip transmissibilities were manually tuned to improve the history matching results during the waterflooding and polymer flooding periods, respectively. It has been found that higher strip transmissibilities match the sharp water cut increase very well in the waterflooding period. Then the strip transmissibilities need to be reduced with time to match the significant water cut reduction. The viscous fingering effect in the reservoir during waterflooding and the restoration of injection conformance during polymer flooding have been effectively represented. Based on the validated simulation model, numerical simulation tests have been conducted to investigate the oil recovery performance under different development strategies, with consideration for sensitivity to polymer parameter uncertainties. The oil recovery factor with polymer flooding can reach about 39% in 30 years, twice as much as forecasted with continued waterflooding. Besides, the updated reservoir model has been successfully employed to forecast polymer utilization, a valuable parameter to evaluate the pilot test’s economic efficiency. All the investigated development strategies indicate polymer utilization lower than 3.5 lbs/bbl in 30 years, which is economically attractive.

Wang, Xindan↗
DOE OSTI2015

Computational hydrodynamics and optical performance of inductively-coupled plasma adaptive lenses

This study addresses the optical performance of a plasma adaptive lens for aero-optical applications by using both axisymmetric and three-dimensional numerical simulations. Plasma adaptive lenses are based on the effects of free electrons on the phase velocity of incident light, which, in theory, can be used as a phase-conjugation mechanism. A closed cylindrical chamber filled with Argon plasma is used as a model lens into which a beam of light is launched. The plasma is sustained by applying a radio-frequency electric current through a coil that envelops the chamber. Four different operating conditions, ranging from low to high powers and induction frequencies, are employed in the simulations. The numerical simulations reveal complex hydrodynamic phenomena related to buoyant and electromagnetic laminar transport, which generate, respectively, large recirculating cells and wall-normal compression stresses in the form of local stagnation-point flows. In the axisymmetric simulations, the plasma motion is coupled with near-wall axial striations in the electron-density field, some of which propagate in the form of low-frequency traveling disturbances adjacent to vortical quadrupoles that are reminiscent of Taylor-Görtler flow structures in centrifugally unstable flows. Although the refractive-index fields obtained from axisymmetric simulations lead to smooth beam wavefronts, they are found to be unstable to azimuthal disturbances in three of the four three-dimensional cases considered. The azimuthal striations are optically detrimental, since they produce high-order angular aberrations that account for most of the beam wavefront error. A fourth case is computed at high input power and high induction frequency, which displays the best optical properties among all the three-dimensional simulations considered. In particular, the increase in induction frequency prevents local thermalization and leads to an axisymmetric distribution of electrons even after introduction of spatial disturbances. The results highlight the importance of accounting for spatial effects in the numerical computations when optical analyses of plasma lenses are pursued in this range of operating conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Behavior of lithium ions in the turbulent near-wall tokamak plasma under heating of ions and electrons of the main plasma

Turbulent dynamics of the near-wall tokamak plasma is simulated by numerically solving the nonlinear reduced Braginskii magnetohydrodynamic equations with allowance for a lithium ion admixture. The effects of turbulence and radiation of the admixture are analyzed in the framework of a self-consistent approach. The radial distributions of the radiative loss power and the density of Li{sup 0} atoms and Li{sup +1} ions are obtained as functions of the electron and ion temperatures of the main plasma in the near-wall layer. The results of numerical simulations show that supply of lithium ions into the low-temperature near-wall plasma substantially depends on whether the additional power is deposited into the electron or ion component of the main plasma. If the electron temperature in the layer increases (ECR heating), then the ion density drops. At the same time, an increase in the temperature of the main ions (ICR heating) leads to an increase in the density of Li{sup +1} ions. The results of numerical simulations are explained by the different influence of the electron and ion temperatures on the atomic processes governing the accumulation and loss of particles in the balance equations for neutral Li{sup 0} atoms and Li{sup +1} ions in the admixture. The radial profile of the electron temperature and the corresponding distribution of the radiative loss power for different densities of neutral Li{sup 0} atoms on the wall are obtained. The calculations show that the presence of Li{sup +1} ions affects turbulent transport of the main ions. In this case, the electron heat flux increases by 20–30% with increasing Li{sup +1} density, whereas the flux of the main ions drops by nearly the same amount. The radial profile of the turbulent flux of lithium ions is obtained. It is demonstrated that the appearance of the pinch effect is related to the positive density gradient of lithium ions across the calculation layer. For the parameters of the T-10 tokamak, the effect of radiative cooling of the near-wall plasma layer becomes appreciable when the near-wall density of neutral lithium atoms exceeds 7 × 10{sup 11} cm{sup −3}. In this case, the density of radiative loss power in the center of the layer is estimated to be about 500–600 kW/m{sup 3}.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Extracting Non-propagating Oscillatory Fields in Concrete to Detect Distributed Cracking

We introduce work to find and locate distributed subsurface cracks in concrete by extracting non-propagating oscillatory fields. The medium of interest is concrete, but the approach also applies to other types of inhomogeneous media. The theoretical basis of the work is first presented through a one-dimensional point-scatterer model that considers the wavefield set up by multiple distinct scatterers. More complex scattering scenarios are then investigated using numerical simulation. The numerical models consider two types of scatterers: elliptic large-scale particles distributed throughout a medium and small-sized cracks localized within a damage zone. The theoretical and numerical analyses show that forward propagating waves undergo distinct scattering behavior within the crack damaged zone: non-propagating resonance-like oscillatory fields are set up within the cracked zone that are distinct from the scatter caused by the large-scale particles. We employ frequency-wavenumber (f-k) domain analysis to extract the energy of non-propagating oscillatory fields and thus to detect and locate zones of distributed cracking. The proposed approach is evaluated using numerical simulation and experimental data collected from concrete specimens that contain simulated distributed cracks. The results reflect the location of distributed crack zones in discrete random media such as concrete can be successfully detected.

36 MATERIALS SCIENCE↗
DOE OSTI2020

DECOVALEX-2019 (Task G Final Report)

The DECOVALEX Project is an on-going international research collaboration, established in 1992, to advance the understanding and modeling of coupled Thermal (T), Hydrological (H), Mechanical (M) and Chemical (C) processes in geological in geological systems. This document is the final report of Task G which was proposed by the Swedish Radiation Safety Authority (SSM) and coordinated by geomecon GmbH, presenting the technical definitions of the problems studied, approaches applied, achievements made and outstanding issues for future research. Task G in DECOVALEX-2019 is primarily concerned with the evolution of transmissivity throughout the lifetime of a repository for radioactive waste, and in particular of spent nuclear fuel, in sparsely fractured and competent crystalline rock masses, and the representation of these processes via numerical simulation. The numerical simulations were to be evaluated for their suitability, and to be validated. Over the course of the work presented here, not only were transmissivity changes and the best approaches to simulate those change examined, but also strategies were to be developed on how to monitor transmissivity change. Guidelines on what a control program and a monitoring system for a repository should include and be able to measure, and how this can be done in practice, was a key objective of the research project.

58 GEOSCIENCES↗
DOE OSTI2016

RELATIVISTIC ELECTRONS PRODUCED BY RECONNECTING ELECTRIC FIELDS IN A LASER-DRIVEN BENCH-TOP SOLAR FLARE

Laboratory experiments have been carried out to model the magnetic reconnection process in a solar flare with powerful lasers. Relativistic electrons with energy up to megaelectronvolts are detected along the magnetic separatrices bounding the reconnection outflow, which exhibit a kappa-like distribution with an effective temperature of ∼10{sup 9} K. The acceleration of non-thermal electrons is found to be more efficient in the case with a guide magnetic field (a component of a magnetic field along the reconnection-induced electric field) than in the case without a guide field. Hardening of the spectrum at energies ≥500 keV is observed in both cases, which remarkably resembles the hardening of hard X-ray and γ -ray spectra observed in many solar flares. This supports a recent proposal that the hardening in the hard X-ray and γ -ray emissions of solar flares is due to a hardening of the source-electron spectrum. We also performed numerical simulations that help examine behaviors of electrons in the reconnection process with the electromagnetic field configurations occurring in the experiments. The trajectories of non-thermal electrons observed in the experiments were well duplicated in the simulations. Our numerical simulations generally reproduce the electron energy spectrum as well, except for the hardening of the electron spectrum. This suggests that other mechanisms such as shock or turbulence may play an important role in the production of the observed energetic electrons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Neutrino transport in general relativistic neutron star merger simulations

Abstract Numerical simulations of neutron star–neutron star and neutron star–black hole binaries play an important role in our ability to model gravitational-wave and electromagnetic signals powered by these systems. These simulations have to take into account a wide range of physical processes including general relativity, magnetohydrodynamics, and neutrino radiation transport. The latter is particularly important in order to understand the properties of the matter ejected by many mergers, the optical/infrared signals powered by nuclear reactions in the ejecta, and the contribution of that ejecta to astrophysical nucleosynthesis. However, accurate evolutions of the neutrino transport equations that include all relevant physical processes remain beyond our current reach. In this review, I will discuss the current state of neutrino modeling in general relativistic simulations of neutron star mergers and of their post-merger remnants. I will focus on the three main types of algorithms used in simulations so far: leakage, moments, and Monte-Carlo scheme. I will review the advantages and limitations of each scheme, as well as the various neutrino–matter interactions that should be included in simulations. We will see that the quality of the treatment of neutrinos in merger simulations has greatly increased over the last decade, but also that many potentially important interactions remain difficult to take into account in simulations (pair annihilation, oscillations, inelastic scattering).

Foucart, Francois (ORCID:0000000346174738)↗
DOE OSTI2012

Evidence of critical balance in kinetic Alfven wave turbulence simulations

A numerical simulation of kinetic plasma turbulence is performed to assess the applicability of critical balance to kinetic, dissipation scale turbulence. The analysis is performed in the frequency domain to obviate complications inherent in performing a local analysis of turbulence. A theoretical model of dissipation scale critical balance is constructed and compared to simulation results, and excellent agreement is found. This result constitutes the first evidence of critical balance in a kinetic turbulence simulation and provides evidence of an anisotropic turbulence cascade extending into the dissipation range. We also perform an Eulerian frequency analysis of the simulation data and compare it to the results of a previous study of magnetohydrodynamic turbulence simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗