Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,009 records · Page 56

Lieb-Schultz-Mattis, Luttinger, and 't Hooft - anomaly matching in lattice systems

We analyze lattice Hamiltonian systems whose global symmetries have ’t Hooft anomalies. As is common in the study of anomalies, they are probed by coupling the system to classical background gauge fields. For flat fields (vanishing field strength), the nonzero spatial components of the gauge fields can be thought of as twisted boundary conditions, or equivalently, as topological defects. The symmetries of the twisted Hilbert space and their representations capture the anomalies. We demonstrate this approach with a number of examples. In some of them, the anomalous symmetries are internal symmetries of the lattice system, but they do not act on-site. (We clarify the notion of “on-site action.”) In other cases, the anomalous symmetries involve lattice translations. Using this approach we frame many known and new results in a unified fashion. In this work, we limit ourselves to 1+1d systems with a spatial lattice. In particular, we present a lattice system that flows to the c=1 c = 1 compact boson system with any radius (no BKT transition) with the full internal symmetry of the continuum theory, with its anomalies and its T-duality. As another application, we analyze various spin chain models and phrase their Lieb-Shultz-Mattis theorem as an ’t Hooft anomaly matching condition. We also show in what sense filling constraints like Luttinger theorem can and cannot be viewed as reflecting an anomaly. As a by-product, our understanding allows us to use information from the continuum theory to derive some exact results in lattice model of interest, such as the lattice momenta of the low-energy states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Apparent fine tunings for field theories with broken space-time symmetries

We exhibit a class of effective field theories that have hierarchically small Wilson coefficients for operators that are not protected by symmetries but are not finely tuned. These theories possess bounded target spaces and vacua that break space-time symmetries. We give a physical interpretation of these theories as generalized solids with open boundary conditions. We show that these theories realize unusual RG flows where higher dimensional (seemingly irrelevant) operators become relevant even at weak coupling. Finally, we present an example of a field theory whose vacuum energy relaxes to a hierarchically small value compared to the UV cut-off.

Physics↗

Lieb-Schultz-Mattis anomalies as obstructions to gauging (non-on-site) symmetries

We study ’t Hooft anomalies of global symmetries in 1+1d lattice Hamiltonian systems. We consider anomalies in internal and lattice translation symmetries. We derive a microscopic formula for the “anomaly cocycle” using topological defects implementing twisted boundary conditions. The anomaly takes value in the cohomology group H^3(G,U(1)) × H^2(G,U(1)) H 3 ( G , U ( 1 ) ) × H 2 ( G , U ( 1 ) ) . The first factor captures the anomaly in the internal symmetry group G, and the second factor corresponds to a generalized Lieb-Schultz-Mattis anomaly involving G and lattice translation. We present a systematic procedure to gauge internal symmetries (that may not act on-site) on the lattice. We show that the anomaly cocycle is the obstruction to gauging the internal symmetry while preserving the lattice translation symmetry. As an application, we construct anomaly-free chiral lattice gauge theories. We demonstrate a one-to-one correspondence between (locality-preserving) symmetry operators and topological defects, which is essential for the results we prove. We also discuss the generalization to fermionic theories. Finally, we construct non-invertible lattice translation symmetries by gauging internal symmetries with a Lieb-Schultz-Mattis anomaly.

Physics↗

Lattice phase shifts and mixing angles for an arbitrary number of coupled channels

We present a lattice method for determining scattering phase shifts and mixing angles for the case of an arbitrary number of coupled channels. The proposed method combines a spherical wall boundary condition and a channel-mixing auxiliary potential to extract the full-rank S-matrix from the radial wave functions. We consider the scattering problem of two spin-1 bosons interacting with a test potential involving up to four coupled channels. For this benchmark system, the phase shifts and mixing angles are shown to agree on the lattice and in the continuum. Our method should allow to extend previous two-channel nuclear lattice EFT simulations to mixing of more than two partial waves.

Bovermann, Lukas↗

A Guideline for the Assessment of Uniaxial Creep and Creep-Fatigue Data and Models

Of primary concern to FE practitioners is a determination of which constitutive models are the “best”, capable of reproducing the mechanisms expected in an intended design accurately; as well as what experimental datasets are proper or “best” to use for fitting the constitutive parameters needed for the model(s) of interest. In addition, the inherent uncertainty in experiments and the repeatability and stability of extrapolations (for long-lived FE components) is a major concern. The goal of this project is to develop a guideline for the assessment of uniaxial creep and creep-fatigue data and models. The research objectives are RO1) develop an aggregated database of creep and creep-fatigue data for P91 steel and 316 stainless steel and perform a systematic assessment of the uncertainty and integrity of the datasets. RO2) employ the database to benchmark existing creep and creep-fatigue models to evaluate the model’s ability to predict standard and nonstandard boundary conditions.

20 FOSSIL-FUELED POWER PLANTS↗

Porflow Implementation of Vadose Zone Conceptual Model for Slit and Engineered Trenches in the E-Area Low Level Waste Facility Performance Assessment

A summary of the conceptual models to be used in representing slit (ST) and engineered trenches (ET) in PORFLOW (ACRi 2018) simulations as part of the performance assessment (PA) of the E-Area Low Level Waste Facility (ELLWF) is presented. Key details that are discussed include: model geometries, spatially dependent hydro-stratigraphic representations, model dimensionality, and boundary conditions (i.e., infiltration rates, cover overhangs, subsidence). Accounting for differences in the percent of non-crushable materials, eighteen unique models, defined by seven hydro-stratigraphic groupings, will be used to represent slit and engineered trenches.

42 ENGINEERING↗

Characterizing Dynamic Test Fixtures Through the Modal Projection Error

Across many industries and engineering disciplines, physical components and systems of components are designed and deployed into their environment of intended use. It is the desire of the design agency to be able to predict whether their component or system will survive its physical environment or if it will fail due to mechanical stresses. One method to determine if the component will survive the environment is to expose the component to a simulation of the environment in a laboratory. One difficulty in doing this is that the component may not have the same boundary condition in the laboratory as is in the field configuration. This paper presents a novel method of quantifying the error in the modal domain that occurs from the impedance difference between the laboratory test fixture and the next level of assembly in the field configuration. The error is calculated from the projection from one mode shape space to the other, and the error is in terms of each mode of the field configuration. This provides insight into the effectiveness of the test fixture with respect to the ability to recreate the mode shapes of the field configuration. A case study is presented to show that the error in the modal projection between two configurations is a lower limit for the error that can be achieved by a laboratory test.

42 ENGINEERING↗

Shocks in the Scrape-off Layer: Hydrodynamical Aspects (Summary) [Slides]

Hydrodynamical model provides a robust description of mass, momentum, and energy transfer. It can be used for benchmarking hydrodynamical modules of SOL codes. It clarifies the role of boundary conditions at the end-plates. It provides a “geometrical” model for partition of power flows between two divertor legs, both for the steady-state and pulsed flows. It indicates possible presence of shocks in the divertor legs. A framework for finding the shock location and strength has been developed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

ZND Verification Tests for Reactive Burn Models in FLAG

The ZND theory, named after Zeldovich, von Neumann, and Doering, provides a simple model for one-dimensional ideal steady-state detonation. It assumes that the detonation wave front starts with a shock that is a discontinuous jump and is followed by a finite-length reaction zone. Reactive burn (also called reactive flow) models are based on ZND theory, as they model the shock initiation and detonation process with a finite reaction rate. The ZND wave propagation test is essentially the only available test case where an analytic solution exists for verification of reactive burn models in numerical codes. However, there are extensions and variants of the ZND test that have been devised for verification of multidimensional flows. The objective of this work is to provide verification of the reactive burn models currently implemented in the Lagrangian hydrocode FLAG and investigate the influence of mesh resolution, artificial viscosity models, and the Arbitrary Lagrangian-Eulerian (ALE) Euler relaxer on the simulation results. The burn models of interest are the Wescott-Stewart-Davis (WSD) model, the Scaled Uniform Reactive Flow (SURF) model (specifically with the SURFplus model extension), and the Arrhenius shock temperature state dependent WSD (AWSD) model. Previously, Ralph Menikoff has used ZND tests for verification of the SURF and SURFplus models in the Eulerian hydrocode xRAGE. The ZND tests here are somewhat different than the approach by Menikoff. In particular, we use a piston-driven ZND detonation wave (via a prescribed constant velocity boundary condition) in a Lagrangian framework whereas Menikoff had a ZND wave followed by a invariant rarefaction wave. The xRAGE simulations were carried out on uniform grids and adaptive mesh refinement (AMR) grids. Although AMR was recently implemented in FLAG and now fully functional for 2D simulations, it will not be evaluated with ZND tests at this time. Some work has been done previously for validating the reactive burn models in FLAG. For example, the validation studies of SURF with shock-to-detonation (SDT) tests, cylinder tests, and gap-stick tests. Further validation of the AWSD, WSD, and SURF models is described in for SDT, multi-shock, cylinder, and corner-turning tests. Recently, a large parameter study with approximately twenty SDT tests was performed to validate the AWSD, WSD, SURF, and SURFplus models while also investigating mesh resolution and artificial viscosity settings. To the best of our knowledge, the current work represents the first documented verification of these burn models in FLAG. The reader should note that many of the tables and figures in this report use units of cm/g/µs, which are the default units for FLAG. However, some lengths and velocities are expressed in µm and mm/µs, respectively, which are typical units for detonation analyses.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

nPINNs: nonlocal Physics-Informed Neural Networks for a parametrized nonlocal universal Laplacian operator. Algorithms and Applications

Physics-informed neural networks (PINNs) are effective in solving inverse problems based on differential and integro-differential equations with sparse, noisy, unstructured, and multifidelity data. PINNs incorporate all available information, including governing equations (reflecting physical laws), initial-boundary conditions, and observations of quantities of interest, into a loss function to be minimized, thus recasting the original problem into an optimization problem. In this paper, we extend PINNs to parameter and function inference for integral equations such as nonlocal Poisson and nonlocal turbulence models, and we refer to them as nonlocal PINNs (nPINNs). The contribution of the paper is three-fold. First, we propose a unified nonlocal Laplace operator, which converges to the classical Laplacian as one of the operator parameters, the nonlocal interaction radius δ goes to zero, and to the fractional Laplacian as δ goes to infinity. This universal operator forms a super-set of classical Laplacian and fractional Laplacian operators and, thus, has the potential to fit a broad spectrum of data sets. We provide theoretical convergence rates with respect to δ and verify them via numerical experiments. Second, we use nPINNs to estimate the two parameters, δ and α, characterizing the kernel of the unified operator. The strong non-convexity of the loss function yielding multiple (good) local minima reveals the occurrence of the operator mimicking phenomenon, that is, different pairs of estimated parameters could produce multiple solutions of comparable accuracy. Third, we propose another nonlocal operator with spatially variable order α(γ), which is more suitable for modeling turbulent Couette flow. Our results show that nPINNs can jointly infer this function as well as δ. More importantly, these parameters exhibit a universal behavior with respect to the Reynolds number, a finding that contributes to our understanding of nonlocal interactions in wall-bounded turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Structure and Dynamics of Domains in Ferroelectric Nanostructures – Phase-Field Modeling (Final Technical Report)

Domain pattern formation is one of the most common phenomena in nature and is a topic of immense interest in many fields ranging from materials science and physics to chemistry and biology. This DOE sponsored research project explored the basic science concerning the thermodynamic stability of mesoscale polarization domain patterns and their temporal evolution mechanisms during formation and subsequent switching in ferroelectric nanostructures and heterostructures. The project employed the computational phase-field method in combination of microelasticity and electrostatic theories. The research was carried out in close collaboration with a number of experimental groups who used High Resolution Transmission Electron Microscopy (HRTEM), in situ TEM with Scanning Probe Microscopy (SPM), or Piezoresponse Force Microscopy (PFM) to characterize the domain structures and dynamics of ferroelectric thin films and heterostructures and who grow high-quality ferroelectric and multiferroic thin films using advanced growth techniques such as Molecular Beam Epitaxy (MBE), Pulsed Laser Deposition (PLLD), and sputtering. The research efforts of the project help establish the phase-field method as the most powerful method for understanding and predicting domain structures in ferroelectric thin films and nanostructures. The findings of the project led to the basic understanding of stability and switching mechanisms of ferroelectric domains under different mechanical boundary conditions and under either homogeneous capacitor configurations or local fields using metallic probes as electrodes and the emergence of charged domain walls during domain switching. The project predicted the spatial length scales, temperature ranges, and electromechanical conditions for different polar states in heterostructures and guided the discovery of both transient and stable novel polarization states containing vortex lattices in oxide superlattices. The project resulted in 134 journal publications and 6 PhD theses with all the PhD graduates currently working in either academia or industry within the United States. The basic understanding on the stability of mesoscale polar states and pattern evolution achieved by the project improved our ability to control and engineer properties of ferroelectric thin films and heterostructures for potential applications in nanoscale electronic devices.

36 MATERIALS SCIENCE↗

Direct E&M Field Solver for a Plane Wave without Sources

A three dimensional direct field solver for the electric and magnetic fields using the Finite Difference Time Domain (FDTD) method is described. Example simulations using periodic and absorbing boundary conditions in a vacuum medium for a plane wave propagating in three dimensions are explored and are compared to theory.

97 MATHEMATICS AND COMPUTING↗

Enduring Transuranic Liquid Waste Volume Forecast and Constituent Evaluation (FY20)

The “FY11 Radioactive Liquid Waste (RLW) Influent Boundary Condition Assessment and Forecast” report provided influent volume forecasts that were critical in establishing the design capacity for the new Radioactive Liquid Waste Transuranic (RLW-TRU) capability. The design for the new RLW-TRU capability was paused as the final design was being completed and prior to the transition of the new management and operations contract for Los Alamos National Laboratory (LANL). Recently, NNSA has directed the project to restart and complete the design effort. As part of the restart process, LANL was asked to validate or document any material changes in the influent volumes or constituents that formed the design basis for the new capability. This analysis updates the transuranic liquid waste influents volumes, incorporating the changes in the LANL program of record since the FY11 assessment, and includes a cursory validation of influent constituents. The results of the FY19 update to influent volume forecasts shows: 1) the influent volume of 29,000 L should continue to be used as the project design-basis; 2) There have been no material changes in the operational processes that would change the constituents of concern from either a treatment or material compatibility perspective.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Advanced Modular Sub-Atmospheric Hybrid Heat Engine (Phase 1 Final Report)

The Phase 1 Final Report describes the results of the work completed during the “Advanced Modular Sub-Atmospheric Hybrid Heat Engine” project. The objectives of Phase 1 work were to develop a conceptual design of the turbine-based modular hybrid heat engine (MHHE), including thermodynamic cycle analysis, define the nominal engine component boundary conditions, identify technology gaps, and develop a detailed test plan to address these gaps through bench scale testing in a potential subsequent Phase 2 project. To achieve the project objective, the following activities were performed in Phase 1 work: (1) Market analysis and primary fuel selection, (2) Conceptual design and integrated layout development, (3) Design support and key components characterization, (4) Planning to move technology to pre-commercial demonstration.

42 ENGINEERING↗

Estimation of Respirable Aerosol Release Fractions through Stress Corrosion Crack-Like Geometries

The formation of a stress corrosion crack (SCC) in the canister wall of a dry cask storage system (DCSS) has been identified as a potential issue for the long-term storage of spent nuclear fuel. The presence of an SCC in a storage system could represent a through-wall flow path from the canister interior to the environment. Modern, vertical DCSSs are of particular interest due to the significant backfill pressurization of the canister, up to approximately 800 kPa. This pressure differential offers a relatively high driving potential for blowdown of any particulates that might be present in the canister. In this study, the carrier gas flow rates and aerosol transmission properties were evaluated for an engineered microchannel with characteristic dimensions similar to those of an SCC. The microchannel was formed by mating two gage blocks with a slot orifice measuring 28.9 μm (0.0011 in.) tall by 12.7 mm (0.500 in.) wide by 8.86 mm (0.349 in.) long (flow length). Surrogate aerosols of cerium oxide, Ce0 2 , were seeded and mixed inside a pressurized tank. The aerosol characteristics were measured immediately upstream and downstream of the simulated SCC at elevated and ambient pressures, respectively. These data sets are intended to demonstrate a new capability to characterize SCCs under well-controlled boundary conditions. Separate modeling efforts are also underway that will be validated using these data. The test apparatus and procedures developed in this study can be easily modified for the evaluation of more complex SCC-like geometries including laboratory-grown SCC samples.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Stability and Control of Burning Tokamak Plasmas with Resistive Walls (Second Year Report/ Final Report)

This research will focus on quantitative prediction of the stability, control, and equilibrium state evolution in toroidal burning plasmas. The stability of long pulse burning plasmas is highly sensitive to the physics of resonant layers in the plasma, sources of momentum and ow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal ow equilibrium state, sustained primarily by energetic alpha particles from fusion reactions, will require the consideration of all of these key elements to predict quantitatively the stability and evolution. The principal investigators on this proposal are leading experts in the relevant theoretical and computational areas, and aim to perform computations guided by analytic modeling, to address this physics in realistic configurations. The overall goal is to understand the key physics mechanisms that describe resistive toroidal burning plasmas, surrounded by a resistive wall, under active feedback control. With the physics of the energetic ions, resonant layers, resistive wall, and toroidal momentum transport included, this study will extend from recent publications in theory and simulation of individual effects and move toward predictive modeling for burning plasmas.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Incorporating physical constraints into Gaussian process surrogate models (LDRD Project Summary)

This report summarizes work done under the Laboratory Directed Research and Development (LDRD) project titled "Incorporating physical constraints into Gaussian process surrogate models?' In this project, we explored a variety of strategies for constraint implementations. We considered bound constraints, monotonicity and related convexity constraints, Gaussian processes which are constrained to satisfy linear operator constraints which represent physical laws expressed as partial differential equations, and intrinsic boundary condition constraints. We wrote three papers and are currently finishing two others. We developed initial software implementations for some approaches. This report summarizes the work done under this LDRD.

97 MATHEMATICS AND COMPUTING↗

Blind Modeling Validation Exercises Using the Horizontal Dry Cask Simulator

The U.S. Department of Energy (DOE) established a need to understand the thermal-hydraulic properties of dry storage systems for commercial spent nuclear fuel (SNF) in response to a shift towards the storage of high-burnup (HBU) fuel (> 45 gigawatt days per metric ton of uranium, or GWd/MTU). This shift raises concerns regarding cladding integrity, which faces increased risk at the higher temperatures within spent fuel assemblies present within HBU fuel compared to low-burnup fuel (≤ 45 GWd/MTU). A dry cask simulator (DCS) was built at Sandia National Laboratories (SNL) in Albuquerque, New Mexico to produce validation-quality data that can be used to test the accuracy of the modeling used to predict cladding temperatures. These temperatures are critical to evaluating cladding integrity throughout the storage cycle of commercial spent nuclear fuel. A model validation exercise was previously carried out for the DCS in a vertical configuration. Lessons learned during the previous validation exercise have been applied to a new, blind study using a horizontal dry cask simulator (HDCS). Three modeling institutions – the Nuclear Regulatory Commission (NRC), Pacific Northwest National Laboratory (PNNL), and Empresa Nacional del Uranio, S.A., S.M.E. (ENUSA) – were granted access to the input parameters from the DCS Handbook, SAND2017-13058R, and results from a limited data set from the horizontal BWR dry cask simulator tests reported in the HDCS update report, SAND2019-11688R. With this information, each institution was tasked to calculate peak cladding temperatures and air mass flow rates for ten HDCS test cases. Axial as well as vertical and horizontal transverse temperature profiles were also calculated. These calculations were done using modeling codes (ANSYS/Fluent, STAR-CCM+, or COBRA-SFS), each with their own unique combination of modeling assumptions and boundary conditions. For this validation study, the ten test cases of the horizontal dry cask simulator were defined by three independent variables – fuel assembly decay heat (0.5 kW, 1 kW, 2.5 W, and 5 kW), internal backfill pressure (100 kPa and 800 kPa), and backfill gas (helium and air). The plots provided in Chapter 3 of this report show the axial, vertical, and horizontal temperature profiles obtained from the dry cask simulator experiments in the horizontal configuration and the corresponding models used to describe the thermal-hydraulic behavior of this system. The tables provided in Chapter 3 illustrate the closeness of fit of the model data to the experiment data through root mean square (RMS) calculations of the error in peak cladding temperatures (PCTs), PCT axial locations, axial temperature profiles, vertical and horizontal temperature profiles at two different axial locations, and air mass flow rates for the ten test cases, normalized by the experimental results. The model results are assigned arbitrary model numbers to retain anonymity. Due to the relatively flat axial temperature profiles, small temperature gradients resulted in large deviations of all models’ PCT axial location from the experimental PCT axial location. When the PCT axial location error is excluded in the calculation of the combined RMS of the normalized errors that considers PCT, the temperature profiles, and the air mass flow rates, the model data fits the experimental data to within 5%. When the vault information is excluded, the model data fits the experimental data to within 2.5%. An error analysis was developed further for one model, using the model and experimental uncertainties in each validation parameter to calculate validation uncertainties. The uncertainties for each parameter were used to define quantifiable validation criteria. For this analysis, the model was considered validated for a given comparison metric if the normalized error in that metric divided by the validation uncertainty was less than or equal to 1. When considering the combined RMS of the normalized errors of all metrics divided by their validation uncertainties, the model was found to have satisfied the criterion for model validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗