Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary condition”

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

Effect of Soil Properties and Input Motion on Site Amplification Using Validated Nonlinear Soil Model

Here, three-dimensional (3-D) nonlinear site response analyses are conducted using finite element models of actual soil profiles from ten nuclear power plant (NPP) sites in the United States to investigate the effects of soil properties and input motions on site amplification. The modeling approach developed in this study combines several novel elements, such as 3-D analysis (including vertical motions), nonlinear inelastic behavior of soil (strain-dependent shear modulus reduction and hysteretic damping), formulation of nonreflecting boundary conditions at the base, and generation of realistic outcrop ground motions for specific sites. All these elements of the modeling approach are first validated using actual data from five earthquakes at three downhole array stations recorded in the Kiban-Kyoshin network (KiK-net), Japan. The same approach is then used to develop site models of ten NPP sites in the United States and corresponding ground motions that are spectrally matched to the site hazard spectra. Eight sets of three-component input motions are used in the study and are categorized on the basis of presence or absence of a near-field pulse in the seed ground motions used for spectral matching. It is found that all sites retain a definite site amplification function regardless of the input motion, provided that the seed motion is spectrally matched to the site hazard spectra. The magnitude of site amplification and frequencies at which they occur depend upon soil properties, particularly the shear wave velocity profile and the constitutive relationship (strain-dependent shear modulus reduction and hysteretic damping) of soil. Amplification of spectral acceleration in the vertical direction (up-down motion) is found to be just as much as, if not more than, the amplification in the horizontal direction. Peak shear strain is found to be about 20% larger for near-field motions compared to far-field motions whereas maximum horizontal site amplification for far-field motions is found to be consistently larger than that of near-field motions, even though the differences between the two remain within the scatter resulting from individual ground motions.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Coupling of CTF and RELAP5-3D Within an Enhanced Fidelity Nuclear Power Plant Simulator

A robust and accurate multiphysics engineering simulator is being developed to model the core behavior and system response of pressurized water reactors. This simulator relies on the NESTLE and CTF computer codes to model the neutronics and thermal hydraulics (TH), respectively, inside the core on a nodal scale and on the Reactor Excursion and Leak Analysis Program—Three Dimensional (RELAP5-3D) to model the entire nuclear steam supply system. The RELAP5-3D model includes highly detailed nodalization and multidimensional flow modeling throughout the vessel. Previously, pin-resolved data generated via the Virtual Environment for Reactor Analysis core simulator were used to improve the accuracy of the NESTLE core predictions. The engineering simulator being developed as part of this work uses the 3KEYMASTER platform to couple the enhanced NESTLE model to a nodal-fidelity CTF model to balance run time with accuracy; NESTLE provides node-dependent powers to CTF, and CTF provides node-dependent coolant densities and fuel temperatures to NESTLE.An overlapping domain approach is used for the core TH in which RELAP5-3D provides core boundary conditions based on the system response and CTF provides a node-dependent coolant heating rate to the RELAP5-3D core solution. In the preliminary TH demonstration discussed in this paper, CTF and RELAP5-3D provided similar steady-state core predictions, indicating the hydraulic compatibility between the codes, as well as reasonable and expected behavior under hypothetical transient conditions. This provides an initial step in ongoing efforts toward a robust, multiscale TH/neutronics engineering simulator capability.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Investigation of Abnormal Level Control Oscillations in a BWR Feedwater System

In the long-term operation of nuclear power plants, the aging of systems, structures, and components can lead to maintenance issues that must be dealt with to maintain cost-effective plant operations. One common issue affecting the currently operated boiling water reactors is the onset of unexpected level oscillations in feedwater heaters. This phenomenon can cause excessive cycling of drain valves and lead to premature failures. In this work, we develop a dynamic model of a set of feedwater heaters to determine the root cause of oscillations observed in an operating plant. Simulation results of various transient scenarios were used to investigate the effects of the controller parameters, boundary conditions, and possible valve and instrument issues. The analysis led to the conclusion that the most likely causes of the observed self-sustained oscillations in the system are the nonlinear behaviors of the drain valve and the level transmitter induced by degraded equipment condition. In conclusion, a partial plug of the pressure line used for level sensing in the system can account for a significant deadtime in the level transmitter, a nonlinear effect shown to induce self-sustained oscillatory behaviors.

Boiling water reactors↗
DOE OSTI2023

Code Benchmark of the HTTF Pressurized Conduction Cooldown Test Using SAM

The High Temperature Test Facility (HTTF) at Oregon State University is an integral system test facility to simulate postulated reactor transients of prismatic high-temperature gas-cooled reactors(HTGRs). A series of test campaigns was launched, providing abundant test data that could be used to benchmark reactor system analysis codes like the System Analysis Module (SAM). In this study, a SAM model of the facility is developed based on the two-dimensional (2D) ring model approach. All components including the ceramic matrix, graphite heaters, helium coolant channels, core barrel, upcomer, pressure vessel, and reactor cavity cooling system are modeled as concentric cylindrical rings. The model is used to simulate one of the benchmark problems-Pressurized Conduction Cooldown (PCC)-within the scope of the Organisation for Economic Co-operation and Development Nuclear Energy Agency International HTTF Benchmark. The simulations consist of two parts. In the first part, operating and boundary conditions as well as thermophysical properties of materials are specified for the benchmark problem. In this work, results from the first part will be used in code-to-code comparison. In the second part, the SAM model is used to simulate Test PG-27, which is the first PCC test carried out in the HTTF, with only two of the ten heater banks activated. The results in the second part are used for code-to-data comparison. Because the helium coolant flow rate is not measured in this facility, it is estimated using the input power and inlet/outlet coolant temperatures. Additionally, radial heat flow in the ceramic blocks is complicated by hundreds of cylindrical coolant channels and heater rods embedded in them. As such, it is necessary to deduce an effective thermal conductivity for the ceramic to analyze the core thermal behavior. SAM predictions of the helium coolant and ceramic temperatures are compared with test data measured in three equivalent sectors. Overall, the SAM results agree reasonably well with test data within the variation of data among the three sectors, which demonstrates SAM's capability in capturing transient effects in HTGR using the simplified 2D ring model.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Droplet and Particle Methods to Investigate Turbulent Particle Laden Jets

The SARS-CoV-2 pandemic has heightened the interest in particle-laden turbulent jets generated by breathing, talking, coughing and sneezing, and how these can contribute to disease transmission. We present quantitative measurement methods for such flows, while exploring and offering improvements for common shortcomings. We generate jets consisting of either liquid droplets or solid particles in an isothermal, quiescent and electrically isopotential experimental chamber that was constructed to control the effects of ambient forcing on jet behavior. For liquid droplets, we find promise in surface deposition analysis based on fluorescent tracer use. For particles, we explore the performance of commercially available adhesive sampling strips and develop conductive grounded carbon tape based sampling strips. We explore ways in which the smallest of thermal gradients or electrostatic charge issues can affect particle dispersion, and suggest practical methods to address these issues. The developed methods are applied to study the simultaneous deposition of (Formula presented.) 25, 50 and 200 μm solid particles from a particle laden turbulent jet with a mean velocity of 33.2 m/s. The deposition location as a function of particle size was compared to results from a simple numerical RANS model, and illustrates ways in which imprecise initial or boundary conditions can lead to a notable deviation from experimental results. The differences in deposition pattern seen in experimental and numerical results despite a carefully controlled environment and characterized particle ejection indicate the need for a more stringent numerical model validation, especially when studying fate and transport of mid-range (neither purely aerosol or ballistic) sized particles.

42 ENGINEERING↗
DOE OSTI2021

Modelling Airborne Transmission and Ventilation Impacts of a COVID-19 Outbreak in a Restaurant in Guangzhou, China

Computational fluid dynamics (CFD) modelling was performed to simulate spatial and temporal airborne pathogen concentrations during an observed COVID-19 outbreak in a restaurant in Guangzhou, China. The reported seating configuration, overlap durations, room ventilation, layout, and dimensions were modelled in the CFD simulations to determine relative exposures and probabilities of infection. Results showed that the trends in the simulated probabilities of infection were consistent with the observed rates of infection at each of the tables surrounding the index patient. Alternative configurations that investigated different boundary conditions and ventilation conditions were also simulated. Increasing the fresh-air percentage to 10%, 50%, and 100% of the supply air reduced the accumulated pathogen mass in the room by an average of ~30%, ~70%, and ~80%, respectively, over 73 min. Overall, the probability of infection was reduced by ~10%, 40%, and 50%, respectively.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Mathematical Models and Numerical Methods for High-Fidelity Simulation of Ignition of Reactive Mixtures by Nanosecond Plasma Discharges in Realistic Configurations

We present a newly developed framework for the numerical simulation of ignition of reactive mixtures using single or repeated nanosecond discharge pulses. The framework builds upon the AMReX library, using the existing compressible solver PeleC and low-Mach solver PeleLMeX and allowing for adaptive mesh refinement, complex geometries, and execution on next-generation high-performance computing (HPC) systems. High-fidelity elementary models are adopted for weakly-ionised plasma discharges with significant energy deposition, consistent with nanosecond discharge pulses, and then implemented in the solver. The treatment of non-thermal electrons and charged species, thermodynamics of non-equilbrium species, plasma kinetics, limiting time scales, and boundary conditions for charged species are discussed and addressed for computational efficiency. The framework is demonstrated for three relevant applications: single and multi-pulse discharges in air, single pulse ignition of an ethylene/air mixture, and a three-dimensional plasma discharge in air with temperature stratification. The successful application of the framework demonstrates the feasibility of high-fidelity simulation of ignition of air/hydrocarbon mixtures in three-dimensions with multiple discharge pulses.

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

Quantification of air leakage paths: a comparison of airflow and acoustic measurements

Today, the fan pressurization method is the most frequently used method to evaluate a building's airtightness. However, the localization and quantification of leaks remain difficult. In this paper, an acoustic method is introduced to estimate the leakage size of single leaks. Acoustic and airflow measurements were conducted and compared in laboratory tests within the same boundary conditions. This work aims to investigate if various leak sizes can be predicted using acoustic measurement methods. The test apparatus consists of two chambers, separated by a test wall. This wall represents a single characteristic air leakage path. Various types of wall structures with different slit geometries, wall thicknesses, and insulation materials were investigated. The acoustic measurements were performed with a sound source placed in one chamber and ultrasonic microphones located in both chambers. These results were compared to measured airflows through the test wall to provide estimates of uncertainty in the acoustic approach, which indicate a linear trend. Finally, these laboratory measurements were compared to the same measurements at a real office building. Finally, although the acoustic measurement uncertainty is still significant (greater than +/- 50%), the acoustic method has the potential to give an order of magnitude of single leak sizes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Structural Analysis of the US ITER ECH Transmission Line System

The electron cyclotron heating (ECH) and current drive system is one of the main plasma heating systems for ITER. It uses high-power microwave beams with the power deposition location steerable across the plasma cross section. Microwave power is conveyed via transmission lines (TLs) that run from the gyrotrons in the radio frequency building through the assembly hall and tokamak building to the ECH launchers within the tokamak vacuum vessel. The ECH system includes a vast array of interconnected TL waveguides, in-line components, and support structures.Finite element (FE) modeling provides an essential means of simulating the system, applying loads and determining deflections, rotations, forces, moments, and stresses in order to evaluate various structural and microwave transmission performance metrics. A representative FE model of the overall ECH TL system is developed in ANSYS®. This top-level model defines the centerline of the waveguide system. Waveguide segments are represented by line elements (beams and pipes) with equivalent section properties, and support structures are represented by boundary conditions. A systematic approach is used to model each ECH component with lumped masses and structurally equivalent stiffness matrices or ANSYS superelements.The top-level TL FE model is used to evaluate the various loads (thermal, vacuum, seismic, etc.) and operating scenarios. The top-level model directly calculates stresses in the straight aluminum waveguide segments. The model provides the forces and moments acting on the in-line components for detailed submodel assessments. Displacement results from the top-level analysis feed into a separate microwave performance model to help determine operational efficiency. All TL performance and thermal-structural requirements are met, as specified by the applicable codes and standards, and successfully documented in numerous technical reports and demonstrated at the final design review.

Shanmugasundaram, Aravind↗
DOE OSTI2022

Thermo-Mechanical Distortion of Tungsten-Coated Steel During High Heat Flux Testing Using Plasma Arc Lamps

An experimental setup and a test section were designed and fabricated for high heat flux testing (HHFT) of neutron-irradiated specimens using water-wall plasma arc lamps. Because of the radiological considerations and limitations of reactor irradiation, the size of the test articles was limited to disks less than 10 mm in diameter. The specimen was clamped onto an actively cooled block, and clamping allowed the insertion of several thermocouples on the back surface of the specimen through a copper (Cu) block. Five vacuum plasma sprayed tungsten (W)–coated F82H steel specimens were subjected to HHFT. Surface profilometry measurements, which were conducted after HHFT, revealed central bowing of the top W surface. This type of residual distortion occurred for all of the specimens, and the larger the specimens were, the larger was the distortion.In an attempt to understand specimen distortion and address the science questions related to the testing of subsize specimens during HHFT, a simplified thermo-mechanical model was developed. By using a measured temperature in the Cu as an isothermal boundary condition, the model eliminated the need for coupling cooling fluid flow models with stress models, greatly simplifying the analysis. The main variable in the proposed model is hC, i.e., the thermal contact conductance between the F82H and the Cu washer. Inelastic properties, including hardening properties, were considered for F82H steel and Cu. Numerical simulation results demonstrated a buildup of residual deformation during HHFT and a very complex state of stress and deformation during typical heat flux (HF) cycling. Additionally, hoop stress evolution during a high heat flux cycle reveals that F82H at an interface with W would be mainly in compression during HF application and experienced a transition to a tension state during cooldown. Also, specimen distortion evolves during each HF cycle, as the specimen bows downward during HF application and upward during the cooldown period between HF cycles. The final specimen distortion, i.e., upward bowing of the specimen center, was qualitatively predicted for hC values of 4000 to 5000 W/(m 2 ·K). This hC range of values, for which bulging is obtained, is at the lower spectrum of the range of values for hC, consistent with the low thermal contact conductance expected from the unpolished F82H surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Design and Analysis of the Open-Surface Slow Li Flow Divertor and Comparison to the Fast Li Flow Divertor

In the ongoing U.S. project, “Liquid Metal Plasma Facing Components,” sponsored by the U.S. Department of Energy, efforts have been taken to develop two open-surface divertor designs for the Fusion Nuclear Science Facility using liquid lithium (Li) as a heat and particle flux removal media. The main focus of this study is the design and analysis of a slow (~1 mm/s) and thin (<1 mm) open-surface Li flow divertor with a Li-cooled substrate, which is then compared with an earlier design of a fast (up to 10 m/s) and thick (~0.5 cm) Li flow divertor with the substrate cooled with helium. The slow Li flow divertor design is based on the original LiWall concept developed at the Princeton Plasma Physics Laboratory. Such a thin and slow Li layer can remove the particle flux by reducing the recycling flux, while the heat flux is removed mainly through the heat sink located beneath. In the present study, the heat sink is provided through a Li cooling flow inside the substrate of reduced activation ferritic/martensitic steel. By performing a multiphysics analysis with COMSOL that included liquid-metal magnetohydrodynamics (MHD), heat transfer, and structural mechanics, the impact of various factors on the divertor heat removal capability, such as Li flow velocity, MHD effects, and inlet velocity boundary condition, were examined. Based on comparisons of the two divertor designs, it was shown that the fast-flow divertor significantly outperformed the slow-flow design, whose heat removal capability was limited to ~1 to 2 MW/m2.

Jiang, Yuchen [ORNL] (ORCID:0000000283280284)↗
DOE OSTI2022

Alpert multi-wavelets for functional inverse problems: direct optimization and deep learning

Computational engineering models often contain unknown entities (e.g. parameters, initial and boundary conditions) that require estimation from other measured observable data. Estimating such unknown entities is challenging when they involve spatio-temporal fields because such functional variables often require an infinite-dimensional representation. Here, we address this problem by transforming an unknown functional field using Alpert wavelet bases and truncating the resulting spectrum. Hence the problem reduces to the estimation of few coefficients that can be performed using common optimization methods. We apply this method on a one-dimensional heat transfer problem where we estimate the heat source field varying in both time and space. The observable data is comprised of temperature measured at several thermocouples in the domain. This latter is composed of either copper or stainless steel. The optimization using our method based on wavelets is able to estimate the heat source with an error between 5% and 7%. We analyze the effect of the domain material and number of thermocouples as well as the sensitivity to the initial guess of the heat source. Finally, we estimate the unknown heat source using a different approach based on deep learning techniques where we consider the input and output of a multi-layer perceptron in wavelet form. We find that this deep learning approach is more accurate than the optimization approach with errors below 4%.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Estimating the value of jointly optimized electric power generation and end use: a study of ISO-scale load shaping applied to the residential building stock

A generation-to-load simulation estimated the impact, in terms of production costs and CO2 emissions, attributable to the joint optimization of electric power generation and flexible end uses to support increasing penetrations of renewable energy. Newly conceived, evaluated, and foundational in developing a U.S. National Standard was a transaction-less yet continuous demand response system based on a day-ahead optimum load shape (OLS) designed to encourage Internet-connected devices to autonomously and voluntarily explore options to favour lowest cost generators - without requiring two-way communications, personally identifiable information, or customer opt-in. Boundary conditions used for model calibration included historical weather, residential building stock construction attributes, home appliance and device empirical operating schedules, prototypical power distribution feeder models, thermal generator heat rates, startup and ramping constraints, and fuel costs. Results of an hourly-based annual case study of Texas indicate a 1/3 reduction in production costs and a 1/5 reduction in CO2 emissions are possible.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

A world without pythons would be so simple

We show that bulk operators lying between the outermost extremal surface and the asymptotic boundary admit a simple boundary reconstruction in the classical limit. This is the converse of the Python's lunch conjecture, which proposes that operators with support between the minimal and outermost (quantum) extremal surfaces—e.g. the interior Hawking partners—are highly complex. Our procedure for reconstructing this 'simple wedge' is based on the HKLL construction, but uses causal bulk propagation of perturbed boundary conditions on Lorentzian timefolds to expand the causal wedge as far as the outermost extremal surface. As a corollary, we establish the Simple Entropy proposal for the holographic dual of the area of a marginally trapped surface as well as a similar holographic dual for the outermost extremal surface. Here, we find that the simple wedge is dual to a particular coarse-grained CFT state, obtained via averaging over all possible Python's lunches. An efficient quantum circuit converts this coarse-grained state into a 'simple state' that is indistinguishable in finite time from a state with a local modular Hamiltonian. Under certain circumstances, the simple state modular Hamiltonian generates an exactly local flow; we interpret this result as a holographic dual of black hole uniqueness.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

AdS 3 Einstein gravity and boundary description: pedagogical review

We review the various aspects of the 3D Einstein gravity theory with a negative cosmological constant and its boundary description. We also explore its connections to conformal field theories (CFTs), modular symmetry, and holography. It is worth noting that this particular theory is topological in nature, which means that all the physical degrees of freedom are located on the boundary. Additionally, we can derive the boundary description on a torus, which takes the form of a 2D Schwarzian theory. This observation suggests that the relevant degrees of freedom for the theory can be described using this 2D theory. Because of the renormalizability of the 3D gravity theory, one can probe the quantum regime. This suggests that it is possible to investigate quantum phenomena. Unlike the conventional CFTs, when considering the AdS 3 background, the boundary theory loses modular symmetry. This represents a departure from the usual behavior of CFT and is quite intriguing. The Weyl transformation induces anomaly in CFTs, and we indicate that applying this transformation to the 2D Schwarzian theory leads to similar results. Summing over all geometries with the asymptotic AdS 3 boundary condition is equivalent to summing over a modular group. The partition function is one-loop exact and therefore an analytical expression from the summation. This theory holds potential applications in Quantum Information and is a recurring theme in the study of holography, where gravitational theories are connected with CFTs.

2D Schwarzian theory↗
DOE OSTI2020

Double differential distributions of e-emission in ionization of N 2 by 3,4 and 5 keV electron impact

In this paper, we report the measurement of the absolute double differential cross sections (DDCS) of secondary electrons emitted due to the ionization of N2 molecule in collisions with fast electrons having energies between 3 and 5 keV. The emitted electrons with energies from 1-500 eV have been measured for different forward and backward emission angles. The measured DDCS have been compared with the state-of-the-art first Born approximation with correct boundary condition (CB1) model calculations as well as with the classical trajectory Monte Carlo (CTMC) method. From the measured DDCS, the single differential cross sections (SDCS) as a function of the emission energies have been computed and eventually the total ionization cross sections (TCS) have been derived. The TCS values are also compared with a semi-empirical calculation, namely, the CSP-ic (complex scattering potential-ionization contribution) model.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Electron emission from fluorene (C 13 H 10 ) upon 3.5 MeV/u Si 8+ ion impact: double differential distributions

In this work, we have studied the electron emission from one of the polycyclic aromatic hydrocarbon (PAH) molecules namely, fluorene (C13H10), upon 3.5 MeV/u Si8+ ion impact. The experimentally measured absolute double differential cross sections (DDCS) are compared with the continuum distorted wave-eikonal initial state (CDW-EIS) model and the first Born approximation including correct boundary conditions (CB1). The measurements are carried out in the ejected e⁻-energy range of 1 eV–400 eV and in the angular range of 20°–160°. We have obtained the single differential and the total cross sections (TCSs) of e⁻-emission as well. The CB1 calculation largely underestimates the data. The CDW-EIS model, which is applied for the PAH molecule for the first time, provides an overall better agreement with the double differential, single differential and TCS data. The DDCS data for fluorene has also been compared with that for CH4 molecule, at a few angles. The forward–backward angular asymmetry shows a relatively flatter distribution compared to the theoretical predictions. The contribution due to the giant plasmon resonance could not be clearly observed except a mild indication in the asymmetry parameter. The angular distribution of the carbon KLL Auger electron cross section shows certain variations. The study of the KLL hyper-satellite component indicates the double K-ionization cross section is about 8.6% of the single K-ionization one.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

2D photofragmentation LIF imaging of H 2 O 2 and HO 2 in the effluent of an atmospheric-pressure plasma jet: effects of solid and liquid interfaces

Two-dimensional (2D) absolute measurements of hydrogen peroxide (H 2 O 2 ) and approximations of the hydroperoxyl radical (HO 2 ) in the effluent of a COST Reference Microplasma Jet operated with a He/H 2 O feed gas are presented. Gas-phase densities are mapped using photofragmentation laser-induced fluorescence (PF-LIF) under three boundary conditions: open effluent, a solid target, and a liquid target. A novel method is presented for separating PF-LIF signals from H 2 O 2 and HO 2 using comparative measurements in oxygen-rich and oxygen-free environments to exploit the preferential formation of HO 2 in the presence of molecular oxygen. This separation strategy is supported by results from a plug-flow plasma chemistry model. Measured densities agree closely with model predictions in both magnitude and trend, while the 2D experimental distributions provide additional insight into the spatial dependencies of these species. In particular, the results show distinct differences in species transport depending on the target type: solid surfaces induce lateral deflection and reduced centerline densities, whereas liquid interfaces promote axial accumulation and higher near-axis concentrations.

atmospheric-pressure plasma jet (APPJ)↗