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,081 records · Page 60

Atomistic theories and Simulations of Multiferroics

The goal of this award was to build a research program aimed at investigating puzzling complex phenomena in multiferroics, that are systems possessing electric and magnetic dipoles. We first emphasized the study of multiferroic nanostructures, that are 0D, 1D, 2D systems as well as superlattices and then also study multiferroics as a whole, i.e., bulks as well. In order to fully understand such complex systems, we also decided to study magnetic systems and ferroelectrics, both in their nanostructure and bulk forms. For that and in addition to use standard Density Functional Theory codes, we developed and used ab-initio tools that are able to model properties of large supercells (about a million atoms nowadays) at finite temperature and for any mechanical and boundary conditions. We also collaborated with experimentalists to fully understand the systems under study. The whole program resulted in 108 publications that are listed at the end of this report. In the following, I will put an emphasis on some publications before the last award period and also indicate a description of the last award period.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

PipeSight: A High-Performance Computing Platform for Pipeline Integrity Management

The Phase I feasibility study completed as part of this project has led to a number of innovative technologies being developed and has laid the foundation for a successful Phase II effort to commercialize a platform for managing the integrity of pipelines for the damage mechanisms of the new, hybrid-energy based economy. To ground the development efforts and direction of the project, an extensive market research and customer discovery effort was undertaken early in Phase I. Through this effort, a number of pipeline owners and operators were interviewed, and the following key findings were discovered about the pipeline industry: • Small pipeline operators do not have the central engineering groups necessary to perform their own independent analysis of inspection data, but instead rely on summarized tally sheets provided to them by inspection service providers. • The time it takes to go from an inspection to a completed engineering assessment, even for small segments of pipeline, can take anywhere from 30-120 days. During this delay, critical threats can (and have been known to) cause failures. • Uncertainty is often not accounted for in the assessment of pipeline integrity. The tally sheets provided by third-party service providers are almost always deterministic in nature, identifying threats that present a concern only to the current (not the future) integrity of the pipeline. • It is uncommon to apply the latest technologies to perform advanced assessments of damaged pipelines. There is a desire to use more advanced analysis capabilities to assess threats. Many pipeline operators indicated that they would often excavate a pipeline to perform an inspection and find that the damage was not as bad as they anticipated, thus using limited resources unnecessarily. Companies are not consistent in their use of inspection data to determine corrosion rates, and those that do only calculate deterministic corrosion rates. • The industry has prominently relied on time-based inspections but has recently started to transition to risk-based inspections. However, there appears to be no uniform guidance on how to do so while properly accounting for all sources of uncertainty. • Companies are not storing inspection data in a manner that allows for the ready determination of temporal trends. • Predictive maintenance principles and practices are beginning to be used by early adopters • Some pipelines are being re-purposed to transport different process fluids than they were designed for, e.g., H 2 and CO 2 rich process streams to serve the new hybrid-energy based economy, which are presenting new integrity concerns for the existing pipeline network that crisscrosses the United States. As a result of these discoveries, we were able to target the development efforts in Phase I to best serve the needs of the industry. In Phase I, we developed a way to correlate multiple large-scale scans of the pipeline to determine a probabilistic corrosion rate that accounts for all sources of error and uncertainty in the inspection process. This probabilistic corrosion rate can be used to predict the future thickness distribution of the pipe wall. We demonstrate how this analysis may be performed in an analytical fashion and has been implemented in such a manner that it can be readily distributed using GPU computing through integration of the Kokkos programming model. We also make a very novel extension of the analytical corrosion rate model to Bayesian Networks (an explainable AI technique) that can account for non-parametric distributions of corrosion rates. With the predictions made above for the probabilistic corrosion rate and corresponding future distribution of the pipe wall thickness, we can assess the integrity of the pipeline through the use of a probabilistic engineering assessment. We developed a novel screening data analysis approach that can rapidly identify ‘hotspots’ (local thin areas) where the integrity of the pipeline is a concern. Once more, we implemented this screening approach in C++ to leverage GPU computing via the Kokkos programming model. After the critical hotspots are identified, we developed a program that can automatically generate an advanced finite element model of the damaged regions. Since the number of damaged regions that require advanced analysis can number in the thousands, we integrated an open-source container-native workflow engine for orchestrating parallel jobs on the cloud. Initially, these advanced numerical models were only designed to account for loading due to internal pressure. However, in a slight pivot from the initial Phase I proposal, we developed a complete pipe stress analysis program (called Simflex) which can simulate the complete pipeline and its response to thermal expansion, pressure, thermal bowing, weight, wind, earthquake, support displacement, support friction and external forces. This pipe stress analysis program was written generically, to handle any piping system, but contains the features needed to model long pipelines (i.e., it incorporates a model for soil mechanics and can account for the nonlinear boundary conditions necessary to simulate long underground pipelines). This pipe stress analysis program can simulate any segment of the pipeline (simple or complex) under any set of conditions and loads, to determine the supplemental loads (axial forces and bending moments) at the location of damage. This enables the most accurate state of stress to be accounted for in the pipeline, which can prove critical when evaluating the integrity of a damaged region. In the process of developing the technologies to perform the integrity assessment of the pipeline, we also extended one of the industry standard approaches for performing the assessment of local thin areas that extend more in the circumferential direction than the longitudinal direction of the pipeline. This approach was presented to the API 579-1/AS ME FFS-1 steering committee in November 2021 for consideration in the next edition of the industry standard for Fitness-For-Service (expected to be released in 2023). To help pipeline operators make decisions with the results on any integrity assessment, we developed a new approach to the life-cycle management of pipelines which uses a Bayesian Decision Network. The network is designed to help pipeline operators plan and prioritize inspection activities and ultimately make smarter, more cost-effective decisions. The Bayesian approach accounts for all sources of uncertainty and carries them through to the final optimal decisions, providing a probabilistic framework for optimizing inspection intervals. The proof-of-concept networks developed in the feasibility study are complete, verified, and are focused on a subset of the pipeline. To expand this novel approach to the scale necessary for an entire network of pipelines in Phase II, we will leverage the DOE-funded Bengi solver for industrial-scale decision making with Bayesian Networks [22]. Once implemented, we will be able to provide the pipeline industry with a much-needed tool for optimal inspection planning using truly explainable artificial intelligence (XAI). To handle all of these advanced capabilities into a cloud-based platform, the architecture of the Equity Engineering Cloud (EEC) was extended to include Argo Workflows, a framework capable of distributing and managing a massive number of jobs that consume their own resources, such that thousands of serial finite element simulations can be run in parallel. As part of this substantial undertaking, we also integrated Argo Continuous Delivery (CD) into the EEC, to aid with the rapid prototyping and iterations that will be imperative to the success of the PipeSight platform’s Agile development process in Phase II. As part of the pipe stress analysis program, we also developed a custom visualizer that leverages the DOE-funded VTK visualization library. We added custom contouring capabilities and a means for interacting visually with both the inputs and outputs of the pipe stress analysis program. We also developed routines for automating the post-processing of the finite element simulations to determine if any failure criteria are met and to visualize the deformations, stresses and strains in ParaView using the exodus II file format (a subset of netCDF).

24 POWER TRANSMISSION AND DISTRIBUTION↗

An Infinite Domain 3D Poisson Solver Based on the Barnes-Hut Algorithm

We present a domain decomposition method for the solution of the 3D Poisson equation with infinite domain boundary conditions. The method is based on an application of the Barnes-Hut tree particle scheme adapted to gridded data. Long range interactions are computed using the first two terms in the Cartesian multipole expansion of Green’s function convoluted with the charge while short range computations are performed using Hockney’s domain doubling algorithm. A standard domain decomposition strategy requires O(N 2 ) applications of Hockney’s algorithm, where N is the number of subdomains that intersect that charge support, while in the present approach only O(Nlog 2 N) such computations suffice. The discretization scheme employed is a sixth order Mehrstellen approximation of the 3D Laplace opera tor. The method exhibits satisfactory accuracy at a substantially reduced computational cost compared to the full domain decomposition Hockney’s algorithm.

97 MATHEMATICS AND COMPUTING↗

Earth Sciences Applications for Energy and Global Security: A Few Vignettes [Slides]

Understanding the physics of flow and transport in the subsurface environment is crucial to several National Security applications. In the Energy Security realm, some examples include characterizing radioactive waste disposal, exploring natural gas extraction and sequestering CO2. An example in the Global Security arena is to be able to distinguish between anthropogenic and naturally occurring seismic disturbances, the former being of particular concern as underground nuclear testing activity has seen an increase in recent years. A common theme across these applications is the ability to accurately model fluid flow in porous and/or fractured subsurface media. For decades, researchers have grappled with how to include the effects of uncertainties in these models, including but not limited to heterogeneity, unknown initial and boundary conditions and issues with scalability. More recently advances in computing and machine learning methodologies have allowed a broader exploration of the uncertainty space as we strive for real time decision making. I will present an overview of the field and discuss a few examples of how my own research in these topics has evolved over nearly two decades.

58 GEOSCIENCES↗

A simulation study of test fixture attachment points for vibration test optimization (MEEN 5440 Final Project)

The aerospace industry uses vibration shaker tables to perform component durability testing. These dynamic environments tests often have a hard time replicating the service environment of a component due to a mismatch in the applied force or the boundary condition. Figure 1 shows the desired outcome of a dynamic environments test. The measured responses in a lab test, ideally, should match the measured response of the service environment field test. One way to address the mismatch between field environment and lab environment is with test fixture design. The test fixture attaches the component under test to the shaker table. Single-axis shakers are the industry standard type of shaker table, and, therefore, test fixtures are often rigid in order to minimize cross-axis responses.

42 ENGINEERING↗

FAST-1.0 Software Release Document

The purpose of this document is to detail the release of FAST-1.0 and the documents that accompany the development and testing of FAST-1.0. Design requirements were for FAST-1.0 developed in the NRC Statement of Work (SOW) for NRC Agreement Number NRC-HQ-25-14-D-001 and are transcribed to this document as part of the Software Quality Assurance Plan for the FAST Code System, PNNL-28767. FAST 1.0 represents the merger of the FRAPCON and FRAPTRAN codes and as such FAST-1.0 performs steady state and transient fuel performance calculations described below. FAST-1.0 calculates the steady state response of light-water reactor fuel rods during long-term burnup. The code calculates temperature, pressure, and deformation of a fuel rod as functions of time-dependent fuel rod power and coolant boundary conditions. The phenomena modeled by the code include: 1) heat conduction through fuel and cladding to the coolant; 2) cladding elastic and plastic deformation; 3) fuel-cladding mechanical interaction; 4) fission gas release from the fuel and rod internal pressure; and 5) cladding oxidation. FAST is used to perform independent calculations for regulatory evaluations of fuel performance under normal operation and anticipated operational occurrences (AOOs).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Imaging Nanoscale Energy Transport and Conversion with Ultrafast Electron Microscopy (Final Technical Report)

Light-matter interactions are ubiquitous in nature and reside at the heart of innumerable technologies. The cascade of processes that occur when a material absorbs a photon of light are exceedingly complex and are interwoven in both space and time, rendering precise determination of the atomic-scale and ultrafast mechanisms immensely challenging. The advent of methods for generating short pulses of light several decades ago led to major advances in understanding the initial moments of light absorption and the resultant effects, though directly interrogating the response of the atoms within the material continued to prove challenging. More recently, methods for generating ultrashort pulses of X-rays and fast electrons have opened the way to probing photoinduced structural dynamics of a wide range of matter in multiple phases. An especially promising laboratory-scale method is ultrafast electron microscopy (UEM), wherein the modalities of conventional transmission electron microscopes (imaging, diffraction, spectroscopy) are extended into the femtosecond temporal regime. Here, the ultrafast imaging and diffraction modalities of UEM were used to study the transient structural aspects of photoexcitation of semiconducting materials (e.g., spatially-resolved electron-phonon coupling, excitation and emission of acoustic phonons, and discrete nanoscale scattering processes). The project had three main objectives: (1) determination of the excitation mechanisms of dense, hypersonic charge-carrier waves and the spatially mediated means by which they couple to the lattice via coherent phonon emission, especially with UEM imaging, wherein effects of nanoscale structural and morphological features on the coupling and relaxation dynamics are expected to dictate nucleation sites and preferred wave vectors, (2) elucidation of photoinduced acoustic-phonon seeding, emergence, propagation, and decay over nanoscale crystal regions and especially with respect to local strain fields and atomic-scale disruptions in lattice order, and (3) interwoven with the first two were efforts aimed at realizing combined angstrom-femtosecond spatiotemporal imaging with UEM. The outcomes and impacts of this project were the generation of new knowledge with respect to fundamental light-matter interactions and, especially, the spatially-mediated excitation and evolution of the structural response of materials following coherent photoexcitation. Importantly, the spatial and temporal resolutions of the UEM imaging modalities used are well-suited for such studies and enable spatially-resolved mechanisms to be determined as a function of atomic order, structural features, and morphology. In addition, ultrafast crystallographic measurements were used to correlate real- and reciprocal-space dynamics in order to determine atomic-scale preferential wave vectors and ultrafast scattering mechanisms, especially as dictated by specimen boundary conditions. The obtained results, as detailed in peer-reviewed publications and presentations, illustrate the importance of ultrafast, angstrom-scale real-space imaging for developing a comprehensive understanding of energy transport and conversion in materials.

47 OTHER INSTRUMENTATION↗

Time-average holographic estimation of Poisson's ratio by Cornu's method from mode shapes

More than 150 years ago, Marie Alfred Cornu first demonstrated his method of using interferometry to directly measure Poisson’s ratio from a radially bent beam. It is an elegant approach to measuring Poisson’s ratio of a beam requiring only two length measurements. Cornu’s method combines the high precision of interferometry with a four-point bending apparatus. This research simulates an interferometry technique called time-average scanning digital holography to make direct estimates of Poisson’s ratio from the antin odes of a mode shape. We show it is possible to make direct estimates of Poisson’s ratio from the antinode of a mode shape although the estimates diverge from the true value of Poisson’s ratio. We find that the di vergence is proportional to the distance between the nodes of a mode shape. We then develop an expression for converting our direct estimate back to within two percent of the true value of Poisson’s ratio for a thin, homogeneous, and hinged-hinged beam, or plate. The remainder of the research seeks to extend this method to combinations of clamped, free, and hinged boundary conditions.

42 ENGINEERING↗

Effects of spin-orbit coupling and very large supercells on the description of acceptors in CdTe

Predicting accurate shallow donor and acceptor levels in semiconductors has been quite challenging using periodic boundary conditions as implemented in current density functional theory codes. The reason is that the wave functions associated with the shallow centers are quite extended and are not fully contained in the supercell (typically with a few hundred atoms) so impurities in the periodically repeated image cells interact with each other. Errors of ~0.1 eV are expected, and these are of the same order of magnitude as the ionization energies themselves. In the case of acceptors in CdTe, this problem is exacerbated by the strong spin-orbit coupling that split the Te-related states at the top of the valence band, making the calculations at least 8 times more expensive. CdTe is an important solar-cell material with record high efficiency of 22%. One of the main limiting factors to increasing the efficiency towards the theoretical limit of ~30% is the often reported very low hole concentration. Pushing to the limit of computational capability by using very large supercells with spin-orbit coupling we report the results of hybrid functional calculations of group-V acceptors in CdTe. We show that extrapolation to the dilute limit leads to an interpretation of the experimental data that is qualitatively different from previous DFT and hybrid functional calculation reports. We find that the group-V impurities indeed behave as shallow acceptors and that the corresponding compensating AX-centers are unstable and do not limit p-type doping. We address the differences between our results and previous theoretical predictions and show that our calculated ionization energies predict hole concentrations that are in excellent agreement with recent temperature-dependent Hall measurements on high-quality single-crystal samples.

14 SOLAR ENERGY↗

Nuclear Data Evaluations in the Resolved Resonance Region and Proposed Updates to R­matrix Modeling in the SAMMY Code [Slides]

This detailed presentation starts off with an evaluation workflow scheme along with observables definition. Provided next are typical cases of explicit experimental effects and the basic quantities for reproducibility. Lastly the presentation mentions the Inverse channel, Complex radius, and the level­ dependent boundary condition. A visual workflow scheme is provided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Beam Induced Heating Analysis Update for the EIC Vacuum Chamber Components

One of the challenges of designing the Electron-Ion Col lider (EIC) is to mitigate beam-induced heating due to the intense electron and hadron beams. Heating of the ESR vac uum chamber components is mainly due to beam-induced resistive wall (RW) losses and synchrotron radiation. For the HSR heating is mainly due to large radial offsets and heat conduction from room temperature to cryogenic tem perature for cryo-components. In this paper, we provide an update on the beam-induced heating and thermal analysis of some critical ESR vacuum chamber components, such as the ESR Large Angle Bremhustraulung Monitor (LABM). We also provide a similar update for critical HSR vacuum com ponents, including the cryo-cooled BPM button assembly, beam screen, abort kicker, and polarimeter. To perform the thermal analyses, the resistive wall loss on individual compo nents is calculated with CST and the synchrotron radiation (if exists) is evaluated using SynRad. These losses, along with realistic boundary conditions are then fed to ANSYS to determine the temperature distribution.

43 PARTICLE ACCELERATORS↗

Preliminary Development of Heat Transfer Model-Based Control Algorithms of Liquid Sodium Purification System: Advanced Sensors and Instrumentation Advanced Controls

Monitoring the operation of sodium purification system is essential for efficient operation of sodium fast reactors. In this work, a heat transfer model has been developed for monitoring the plugging meter and cold trap systems at the Mechanisms Engineering Test Loop (METL) liquid sodium facility at Argonne National Laboratory. The model of the purification system was developed by treating the respective aspects of the cold trap purification loop and plugging meter diagnostic loop as two separate control volumes using information from the METL piping and instrumentation diagram (P&ID). A model predictive controller was designed using first order differential equations with the specified boundary conditions. The system behavior was studied with a tuned optimized procedure using the internal cold trap temperature and plugging meter outlet temperature as control variables, and the air blower temperature as an independent variable respectively. Results of computer simulations obtained in this study compared favorably with experimental data showing very good reference tracking response with negligible overshoot as both plugging meter and cold trap physical models approach the setpoint.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Coupling hydrodynamics to a rigid-body motion solver for fluid-structure interaction [Slides]

Quinoa is a massively parallel computational fluid dynamics with multi-material and programmed burn capabilities developed from Programmatic and LDRD funds. Overset is a well-established method of using an “inset” mesh that communicates with a background mesh. Solutions transfer freely from one to another and operate as boundary conditions on the opposite mesh. Either mesh can move at specified velocity! This project couples the existing Mesh-to-mesh transfer with Quinoa and demonstrates use of resulting Overset solver in computational simulation of blast effects on a re-entry body for survivability assessments.

42 ENGINEERING↗

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

This research has focused 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 flow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal flow 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 aimed 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.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SAM Assessment Report

The System Analysis Module (SAM) is a modern system analysis tool being developed at Argonne National Laboratory (Argonne) for advanced non-LWR safety analysis. A series of assessment cases are modeled using SAM and, when applicable, compared against analytical solutions or engineering judgement. This report documents the assessment of SAM for eleven different categories. Each category contains a series of tests that examine different boundary conditions, system configurations, or modeling options. Although relatively simple, these tests cover the basic equation models, basic component models, and basic system level processes and phenomena that may be utilized during advanced reactor safety analyses.

1-D fluid flow↗

NEAMS-Multiphysics MOOSE End Year Framework Activities FY20

The Multiphysics Object Oriented Simulation Environment (MOOSE) is a general finite element package meant for high performance solution of multiphysics problems in science and engineering. In this report, we present additions and enhancements to MOOSE funded by the Nuclear Engineering Advanced Modeling and Simulation (NEAMS) program. NEAMS-funded MOOSE framework library improvements include: expansion of support for multi-level multi-application restart; enhancement of coupling between native MOOSE applications and external libraries; addition of a sparse automatic differentiation (AD) container enabling non-local degree of freedom coupling; block-specific quadrature rules; further development of an eigenvalue executioner; faster setup of periodic boundary conditions; and creation of a mesh meta-data system streamlining simulation startup. Besides developing the framework library, NEAMS funds were used to overhaul a swath of important MOOSE infrastructure in order to substantially improve user experience. These critical infrastructure changes included adapting MOOSE to python 3, improving the test harness to catch race conditions, and most importantly transitioning MOOSE to use Conda, a globally known packaging system, that greatly eases software adoption.

97 MATHEMATICS AND COMPUTING↗

Plasma-assisted surface modification of biopolymers

The interaction of non-thermal plasmas with surfaces is an important problem of interest in various applications. This work focuses on two main aspects of the problem - simulation and experiments. The development of a low-temperature plasma modeling is described along with experiments focusing on the interaction of a low-temperature plasma with liquid water. The simulation framework is capable of handling arbitrary geometries, chemistry and a fully-coupled boundary condition for plasma-dielectric interface. Representative examples for the verification and validation of the framework are provided. The experimental work deals with the study of low-temperature plasmas with water with the nitrate injection rate compared for three different configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tokamak Disruption Simulation (Final Technical Report)

This was a multi-institutional project concerned with simulation of tokamak disruption – rapid termination of the plasma and confining magnetic field in a tokamak fusion reactor – with the aim of quantifying the physics associated with disruption mitigation. The PI’s role is as one of a group of computational scientists concerned with development of effective solution algorithms for the underlying models of magnetohydrodynamics and transport that represent the physical processes, and to develop techniques of uncertainty quantification for treating unknown aspects of the model such as boundary conditions and effects of impurities in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗