Search NASA⌕ Search

SEARCH · Search NASA

Results for “boundary conditions”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

An Eulerian Vlasov-Fokker–Planck algorithm for spherical implosion simulations of inertial confinement fusion capsules

Here, we present a numerical algorithm that enables a phase-space adaptive Eulerian Vlasov–Fokker–Planck (VFP) simulation of inertial confinement fusion (ICF) capsule implosions. The approach relies on extending a recent mass, momentum, and energy conserving phase-space moving-mesh adaptivity strategy to spherical geometry. In configuration space, we employ a mesh motion partial differential equation (MMPDE) strategy while, in velocity space, the mesh is expanded/contracted and shifted with the plasma’s evolving temperature and drift velocity. The mesh motion is dealt with by transforming the underlying VFP equations into a computational (logical) coordinate, with the resulting inertial terms carefully discretized to ensure conservation. To deal with the spatial and temporally varying dynamics in a spherically imploding system, we have developed a novel nonlinear stabilization strategy for MMPDE in the configuration space. The strategy relies on a nonlinear optimization procedure that optimizes between mesh quality and the volumetric rate change of the mesh to ensure both accuracy and stability of the solution. Implosions of ICF capsules are driven by several boundary conditions: (1) an elastic moving wall boundary; (2) a time-dependent Maxwellian Dirichlet boundary; and (3) a pressure-driven Lagrangian boundary. Of these, the pressure-driven Lagrangian boundary driver is new to our knowledge. The implementation of our strategy is verified through a set of test problems, including the Guderley and Van-Dyke implosion problems — the first-ever reported using a Vlasov–Fokker–Planck model.

1D2V↗
DOE OSTI2022

Laminar Flow Heat Transfer In Helical Oval-Twisted Tube for Heat Exchanger Applications

The heat transfer performance of a novel tube configuration that combines the swirling velocity induced by oval-twisting and secondary flow generated by a helical geometrical flow path is presented. The Nusselt number (Nu) and friction factor (f) are compared for the laminar flow regime (Re = 250- 2000) under isothermal wall conditions. Under the same flow and boundary conditions, the oval-twisted helical tube increased the Nu and slightly increased the f over the circular helical tube. The best performance with the highest Nuand lowest f occurs at the coil curvature ratio (d h /D) of 0.17. The quantified enhancement performance factor (h) shows a 46 - 56% increase for the oval-twisted helical tube over the circular helical tube with d h /D = 0.17. Correlations for the Nu and f are developed and agreed with the numerical results: ± 6% and ± 5% for the Nu and f, respectively.

42 ENGINEERING↗
DOE OSTI2024

Nonlinear elasticity with the Shifted Boundary Method

Here, we propose a new unfitted/immersed computational framework for nonlinear solid mechanics, which bypasses the complexities associated with the generation of CAD representations and subsequent body-fitted meshing. This approach allows to speed up the cycle of design and analysis in complex geometry and requires relatively simple computer graphics representations of the surface geometries to be simulated, such as the Standard Tessellation Language (STL format). Complex data structures and integration on cut elements are avoided by means of an approximate boundary representation and a modification (shifting) of the boundary conditions to maintain optimal accuracy. An extensive set of computational experiments in two and three dimensions is included.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Poisson equation method for prescribing fully developed non-Newtonian inlet conditions for computational fluid dynamics simulations in models of arbitrary cross-section

Prescribing inlet boundary conditions for computational fluid dynamics (CFD) simulations of internal flow in complex geometries such as anatomical vascular models is challenging. In the absence of patient-specific inlet velocity data, a common approach for long blood vessels is to assume that the inlet flow is fully developed. In vessels of irregular cross section, however, prescribing fully developed conditions is complicated due to the lack of a general closed-form analytical solution. In this study, we develop a simple Poisson equation method for prescribing fully developed inlet conditions for the flow of either Newtonian or non-Newtonian fluids in CFD models of arbitrary cross-section. We first derive the generalized Poisson equation for fully developed flow of a non-Newtonian fluid and we then develop and verify a methodology for numerically computing the solution on any planar boundary domain. In addition, we develop a simple extension of the method for prescribing a non-orthogonal inlet velocity that represents fully developed flow from an upstream tube that is connected to the CFD inlet at a non-orthogonal angle. This may be used to investigate a common source of uncertainty in CFD simulations of internal flow that is due to a lack of information concerning the exact streamwise flow direction at the inlets. Comparison to several Newtonian and non-Newtonian benchmark verification solutions shows the method to be extremely accurate. As a practical demonstration case, we use the method to prescribe fully developed conditions on multiple non-circular inlets for the non-Newtonian flow of blood in a patient-specific model of the inferior vena cava (IVC). Finally, we further demonstrate the utility of the method by performing a sensitivity study using the patient-specific IVC model, wherein we investigate the influence of inlet velocity flow direction on the non-Newtonian IVC hemodynamics. Given its simplicity and computational efficiency, the method is shown to be far superior to alternative approaches for prescribing fully developed inlet conditions in such complicated geometries. In conclusion, to facilitate the adoption of our Poisson equation method, we have distributed our OpenFOAM source code and the associated test cases from this study as open-source software.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Propagation Pattern for Moment Representation of the Lattice Boltzmann Method

A propagation pattern for the moment representation of the regularized lattice Boltzmann method (LBM) in three dimensions is presented. Using effectively lossless compression, the simulation state is stored as a set of moments of the lattice Boltzmann distribution function, instead of the distribution function itself. An efficient cache-aware propagation pattern for this moment representation has the effect of substantially reducing both the storage and memory bandwidth required for LBM simulations. This article extends recent work with the moment representation by expanding the performance analysis on central processing unit (CPU) architectures, considering how boundary conditions are implemented, and demonstrating the effectiveness of the moment representation on a graphics processing unit (GPU) architecture.

42 ENGINEERING↗
DOE OSTI2020

Parallel three-dimensional simulations of quasi-static elastoplastic solids

Hypo-elastoplasticity is a flexible framework for modeling the mechanics of many hard materials under small elastic deformation and large plastic deformation. Under typical loading rates, most laboratory tests of these materials happen in the quasi-static limit, but there are few existing numerical methods tailor-made for this physical regime. Here, we extend to three dimensions a recent projection method for simulating quasi-static hypo-elastoplastic materials. The method is based on a mathematical correspondence to the incompressible Navier–Stokes equations, where the projection method of Chorin (1968) is an established numerical technique. We develop and utilize a three-dimensional parallel geometric multigrid solver employed to solve a linear system for the quasi-static projection. Our method is tested through simulation of three-dimensional shear band nucleation and growth, a precursor to failure in many materials. As an example system, we employ a physical model of a bulk metallic glass based on the shear transformation zone theory, but the method can be applied to any elastoplasticity model. We consider several examples of three-dimensional shear banding, and examine shear band formation in physically realistic materials with heterogeneous initial conditions under both simple shear deformation and boundary conditions inspired by friction welding.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Systematic Way to Extend the Debye–Hückel Theory beyond Dilute Electrolyte Solutions

An extended Debye–Hückel theory with fourth order gradient term is developed for electrolyte solutions; namely, the electric potential φ(r) of the bulk electrolyte solution can be described by ∇ 2 φ(r) = κ 2 φ(r) + L Q 2 ∇ 4 φ(r), where the parameters κ and L Q are chosen to reproduce the first two roots of the dielectric response function of the bulk solution. Three boundary conditions for solving the electric potential problem are proposed based upon the continuity conditions of involving functions at the dielectric boundary, with which a boundary element method for the electric potential of a solute with a general geometrical shape and charge distribution is derived. Solutions for the electric potential of a spherical ion and a diatomic molecule are found and used to calculate their electrostatic solvation energies. In this work, the validity of the theory is successfully demonstrated when applied to binary as well as multicomponent primitive models of electrolyte solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

A Machine Learning Correction Model of the Winter Clear-Sky Temperature Bias over the Arctic Sea Ice in Atmospheric Reanalyses

Atmospheric reanalyses are widely used to estimate the past atmospheric near-surface state over sea ice. They provide boundary conditions for sea ice and ocean numerical simulations and relevant information for studying polar variability and anthropogenic climate change. Previous research revealed the existence of large near-surface temperature biases (mostly warm) over the Arctic sea ice in the current generation of atmospheric reanalyses, which is linked to a poor representation of the snow over the sea ice and the stably stratified boundary layer in the forecast models used to produce the reanalyses. These errors can compromise the employment of reanalysis products in support of polar research. Here, we train a fully connected neural network that learns from remote sensing infrared temperature observations to correct the existing generation of uncoupled atmospheric reanalyses (ERA5, JRA-55) based on a set of sea ice and atmospheric predictors, which are themselves reanalysis products. The advantages of the proposed correction scheme over previous calibration attempts are the consideration of the synoptic weather and cloud state, compatibility of the predictors with the mechanism responsible for the bias, and a self-emerging seasonality and multidecadal trend consistent with the declining sea ice state in the Arctic. The correction leads on average to a 27% temperature bias reduction for ERA5 and 7% for JRA-55 if compared to independent in situ observations from the MOSAiC campaign (respectively, 32% and 10% under clear-sky conditions). These improvements can be beneficial for forced sea ice and ocean simulations, which rely on reanalyses surface fields as boundary conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Falling-film thermosyphons: Application to water harvesting from humid gas streams

Water harvesting from humid gas streams is central to fresh-water production, desalination, and particulate removal from combustion gas streams. Two-phase thermosyphons are well suited for these applications due to their very low thermal resistance, but have several operating limits at the high heat flows required for such applications. This work introduces the falling-film thermosyphon for large-scale condensation applications. Working fluid is pumped to the top of the evaporator to provide the evaporating liquid film on the inner tube wall, rather than by vapor condensation. The flooding, dry-out and pool-boiling limits are eliminated, resulting in significantly higher maximum heat fluxes for the same physical evaporator size. Furthermore, the condenser no longer needs to be located vertically above the evaporator, which allows for a standard steam condenser to be used. A model for the condensation process in humid air was developed that estimates the fluid and heat transfer and is confirmed experimentally. Further benefits include the use of a high-thermal-conductivity polymer material for the evaporator section to minimize corrosion, and the ability to impose a temperature boundary condition on the evaporator, which is made possible due to the elimination of a liquid pool.

42 ENGINEERING↗
DOE OSTI2020

Experiences with Mobile Units for Liquid RAW Management in Slovak and Czech Nuclear Power Plants - 20089

In Nuclear Power Plants VVER 440, Russian type of reactor there was designed auxiliary building where according to basic design all operational liquid radioactive waste (RAW) should be stored in large capacity tanks during the whole operational period. That's why there are in the auxiliary building several types of large storage tanks with different internal volume from 110 m{sup 3}, 460 m{sup 3} to 550 m{sup 3}. Storage tanks are usually situated in rooms with limited access and minimal handling space. Operational waste was here transported by the system of pipes from the reactor building. There is possibility to handle with the waste and set the pipeline in to the chosen storage tanks. Nobody takes care about the quality of waste in the storage tanks and waste waits for the period of decommissioning when it should be treated. The storage capacity was designed for the whole lifetime of nuclear power plant. Currently the life time of Nuclear Power Plant's (NPP) is extended and there is necessity to retrieve and treat the RAW in advance. Another reason why the stakeholder, who operates the NPPs, decides to treat the operational RAW during the operational period of NPP, is that there are still new requirements from safety point of view which need to have free available space for example because of potential accident. So treating the operational RAW during the life time of NPP is present issue. Because in basic design nobody cares how the RAW will be retrieved and treated in the future, there is usually no technology for that in existing NPP. Some of the NPPs build an additional building with suitable technologies for treatment or hire the company which offer mobile technology for retrieving the waste from the storage tanks. Wood Nuclear Slovakia PLC provides all activities within liquid RAW management projects by its own staff including management, design, construction and deployment of equipment, engineering support, as well as all RAW retrieval and solidification activities. We offer RAW management as a service without any significant changes in existing NPP's buildings and technology. It is always a technical challenge to conduct successful retrieval, pumping and treatment of RAW. The complex RAW retrieval system can be designed as tailored solution for any RAW storage according to known boundary conditions and specific needs. These systems are based on years of experience with RAW manipulation, processing and treatment and the whole RAW management knowledge by Wood Nuclear Slovakia LTD. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

The repulsive Casimir-type forces in binary mixture of weakly interacting Bose–Einstein condensates gases

Highlights: • The equations of the state of a binary mixture weakly interacting Bose gases are established based on the Lagrangian of the system. • The formulae of the Casimir-type energy and the resulting Casimir-type force are produced from the quantum part of the grand canonical energy. • The results for the Casimir-type effect in the one-field case can be re-obtained by using the results of the two-field case. Within the framework of the quantum field theory (QFT) in the one-loop approximation, we investigate the Casimir-type effect in a mixture of two-component weakly interacting Bose-Einstein condensates (BECs) gases confined between two parallel plates at zero temperature. The Casimir-type forces in the miscible state of BECs for different boundary conditions (BCs) are discovered and compared to each other. The repulsive nature of the Casimir-type force in BECs is defined for the first time. Our calculations are finished in both grand canonical ensemble (GCE) and canonical ensemble (CE).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Solving scattering problems in the half-line using methods developed for scattering in the full line

We reduce the solution of the scattering problem defined on the half-line [0,∞) by a real or complex potential v(x) and a general homogeneous boundary condition at x=0 to that of the extension of v(x) to the full line that vanishes for x<0. We find an explicit expression for the reflection amplitude of the former problem in terms of the reflection and transmission amplitudes of the latter. We obtain a set of conditions on these amplitudes under which the potential in the half-line develops bound states, spectral singularities, and time-reversed spectral singularities where the potential acts as a perfect absorber. We examine the application of these results in the study of the scattering properties of a δ-function potential and a finite barrier potential defined in [0,∞), discuss optical systems modeled by these potentials, and explore the configurations in which these systems act as a laser or perfect absorber. In particular, we derive an explicit formula for the laser threshold condition for a slab laser with a single mirror and establish the surprising fact that a nearly perfect mirror gives rise to a lower threshold gain than a perfect mirror. We also offer a nonlinear extension of our approach which allows for utilizing a recently developed nonlinear transfer matrix method in the full line to deal with finite-range nonlinear scattering interactions defined in the half-line.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Improvement of Quantification in Non-Destructive Characterization of Radioactive Waste Packages -Introduction to the Ideas of a Currently Started Project - 20197

For the characterization of radioactive waste packages (e.g. 200-L drums) with inhomogeneous activity and density distributions preferably non-destructive methods like segmented gamma scanning and transmission measurements are applied. While the first is used for identification of the radioactive inventory, both are necessary for a reliable activity quantification. Although at nuclear facilities, segmented gamma scanner are state-of-the-art and are mostly equipped with a vertically moveable transmission source the combining of this data and further a-priori information improving accuracy in data estimation giving realistic, reliable and traceable uncertainty estimations are improvable. This challenge is tackled by a R and D project funded by the German Federal Ministry of Education and Research (BMBF). The boundary conditions for this project are (i) taking into account the typical instrumentation available at nuclear facilities, (ii) minimize the overall measurement time and (iii) combining data from emission and transmission measurements by Bayesian methods. For the R and D project, being still in its initial stage, the basic ideas, the data used for evaluation and the result of a preliminary investigation demonstrating the proof of principle are presented. Additionally, some first aspects on spatial resolution in emission and transmission measurements are addressed. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

A Formulation of the Many-Body Expansion (MBE) for Periodic Systems: Application to Several Ice Phases

A new formulation of the many-body expansion (MBE) for periodic systems, based on their inherent translational symmetry, is introduced and applied to 7 ice polymorphs. This new formulation is built via a hierarchical procedure connecting gas phase unit cells over finite supercells to infinite solids. For finite supercells, it is demonstrated that this method successfully recovers the energetics while reducing the scaling of the calculation of the many-body terms by a factor of ~N, where N is the system size. Furthermore, the proposed framework delineates a straightforward way to obtain properties of macroscopic systems in the limit of an infinite cell. For periodic systems, the success of this approach is demonstrated by showing that the lattice energies computed (up to the 4-th order in the MBE) reproduce the lattice energies obtained using periodic boundary conditions with an Ewald summation for 7 polymorphs of ice (Ih, II, VIII, IX, XIII, XIV, XV). This development makes it possible to quantify, for the first time, the many-body contributions to the lattice energy of various ice polymorphs. The many-body (three-body and higher) interactions were found to vary significantly among the 7 ice phases, amounting between 7-24% of the total lattice energies. So, this development opens the door for obtaining insights into solid-state properties, while leveraging the computational benefits of the MBE.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

DFT-based mesoscopic interactions of N adatoms on Cu(100)

The first-principles density functional theory (DFT) method is used to directly calculate the interactions of N atoms of various configurations on Cu(100) surfaces. This allows us to derive a mesoscopic substratemediated interaction model. We find that, for these Cu(100)-supported two-dimensional N clusters, oscillating substrate-mediated interactions obtained from an elastic eigenvector approach including a mesoscopic interacting range (called the mesoscopic model) have reasonable agreement with the DFT results. The DFT periodic boundary conditions are properly mapped into the mesoscopic model. The model is used to analyze the energies of c(2×2) square islands independent of their sizes. Previous experimental studies have shown that, in the low coverage region, N atoms agglomerate to isolated square islands of about 5 nm × 5 nm in size with a c(2×2) structure. With increasing coverage, those islands form regular patterns but do not touch below about 0.4 monolayers. From previous experiments and calculations, the square island pattern shows significant elastic effects which are now confirmed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Spectral properties of three-dimensional Anderson model

The three-dimensional Anderson model represents a paradigmatic model to understand the Anderson localization transition. In this work we first review some key results obtained for this model in the past 50 years, and then study its properties from the perspective of modern numerical approaches. Our main focus is on the quantitative comparison between the level sensitivity statistics and the level statistics. While the former studies the sensitivity of Hamiltonian eigenlevels upon inserting a magnetic flux, the latter studies the properties of unperturbed eigenlevels. We define two versions of dimensionless conductance, the first corresponding to the width of the level curvature distribution relative to the mean level spacing, and the second corresponding to the ratio of the Heisenberg time and the Thouless time obtained from the spectral form factor. We show that both conductances look remarkably similar around the localization transition, in particular, they predict a nearly identical critical point consistent with other well-established measures of the transition. We then study some further properties of those quantities: for level curvatures, we discuss particular similarities and differences between the width of the level curvature distribution and the characteristic energy studied by Edwards and Thouless (1972) in their pioneering work, in which the hopping at one lattice edge is changed from periodic to antiperiodic boundary conditions. In the context of the spectral form factor, we show that at the critical point it enters a broad time-independent regime, in which its value is consistent with the level compressibility obtained from the level variance. Finally, we test the scaling solution of the average level spacing ratio in the crossover regime using the cost function minimization approach introduced recently in Šuntajs et al. (2020). The latter approach seeks for the optimal scaling solution in the vicinity of the crossing point, while at the same time it allows for the crossing point to drift due to finite-size corrections. We find that the extracted transition point and the scaling coefficient agree with those from the literature to high numerical accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Constraint and size effects in confined layer plasticity

Plane strain finite element analyses are used to model the experiments for a thin metal layer confined between elastic solids. Here, the thin metal layer undergoes elastic-plastic deformations, here modelled by strain gradient plasticity, while the rest of the solid deforms only elastically. Plane strain is assumed and finite strains are accounted for. The microscopic boundary condition at the band interface is initially taken to be zero plastic strain increment, but plastic straining at the interface is allowed when the plastic strain gradient at the interface has exceeded a critical value. When the narrow metal band is inclined 45° relative to the direction of compression of the solid simple shear develops in the band and predictions are rather similar to those of a one-dimensional analysis for an infinite band. For a band at 90° the elastic parts give significant constraint on plastic flow, the stress fields are non-uniform along the band and the average nominal compressive stress shows strong sensitivity to the layer thickness and to the width of the specimen. Remarkably good qualitative, and in some cases quantitative, agreement is found with the experimental observations. Considering also other band angles, both the average nominal shear and compressive stress show variations with the angle.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Vacuum Neutral Transport Model in UEDGE for Tokamak Far Scrape‐Off Layer

A model for neutral transport in the far scrape-off layer (SOL) vacuum region (vacuum neutral model) has been developed and implemented in UEDGE. Free-streaming neutral trajectories between the outermost UEDGE boundary and the vessel wall are preprocessed using DEGAS2 to construct a tele-transport matrix that captures non-local neutral relocation through the vacuum. Here, this matrix is then used in UEDGE as a non-local boundary condition for the neutral equations, preserving the robustness and convergence of the implicit solver without introducing statistical noise. Simulations of a DIII-D lower single-null configuration show that the vacuum neutral model relocates neutrals from the divertor to the upstream region, increasing the outer midplane separatrix density required for detachment onset by about 30%. However, the characteristic target temperature at detachment (T e, osp ~ 3 - 4eV) and the radiation front behavior remain unchanged.

DEGAS2↗