Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupling method”

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 973 records · Page 54

Application of On-Line Monitoring and Real-Time Characterization of Low Level Waste Samples from Hanford Tanks

Waste and treatment stream characterization is an essential part of safe, efficient, and cost-effective processing of Hanford materials. The demand for materials characterization, including chemical, radionuclide, and physical attributes during waste management, transfer, and staging operations in the Hanford tank farms, is expected to significantly increase as the U.S. Department of Energy Office of River Protection (ORP) progresses toward Hanford waste treatment and immobilization. The application of on-line monitoring can greatly improve operation timeframes and reduce costs by providing real-time data and supporting the reduction of necessary grab sample collection. Recently, our team applied Raman online monitoring and chemometric modeling to the supernate of Hanford tank 241-AP-105, and demonstrated the quantitative measurement of nine analytes within this waste. Application of this method has been expanded to additional Hanford tank wastes, AW-102 and AP-107. The performance of the process monitoring instruments and the accuracy of the chemometric models will be assessed. On-line monitoring results will be compared to anion results determined from the current standard, ion-chromatographic and inductively coupled plasma mass spectrometry methods. Comparisons of both determined values and uncertainties of results will be presented. Results will be shared with personnel working on flowsheet interfaces for the Direct-Feed Low Activity Waste mission at Hanford.

Lines, Amanda M.↗

Application of On-Line Monitoring and Real-Time Characterization of Low Level Waste Samples from Hanford Tanks - 20380

Waste and treatment stream characterization is an essential part of safe, efficient, and cost-effective processing of Hanford materials. The demand for materials characterization, including chemical, radionuclide, and physical attributes during waste management, transfer, and staging operations in the Hanford tank farms, is expected to significantly increase as the U.S. Department of Energy Office of River Protection (ORP) progresses toward Hanford waste treatment and immobilization. The application of on-line monitoring can greatly improve operation time frames and reduce costs by providing real-time data and supporting the reduction of necessary grab sample collection. Recently, our team applied Raman online monitoring and chemometric modeling to the supernate of Hanford tank 241-AP-105, and demonstrated the quantitative measurement of nine analytes within this waste. Application of this method has been expanded to additional Hanford tank wastes, AW-102 and AP-107. The performance of the process monitoring instruments and the accuracy of the chemometric models will be assessed. On-line monitoring results will be compared to anion results determined from the current standard, ion-chromatographic and inductively coupled plasma mass spectrometry methods. Comparisons of both determined values and uncertainties of results will be presented. Results will be shared with personnel working on flowsheet interfaces for the Direct-Feed Low Activity Waste mission at Hanford. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Simulation scaling studies of reactor core two-phase flow using direct numerical simulation

Tremendous growth in supercomputing power in recent years has resulted in the emergence of high-resolution flow analysis methods as an advanced research tool to evaluate single and two-phase flow behavior. In particular, unstructured mesh-based methods have been applied to analyze flows in complex reactor core geometries, including those of light water reactors (LWR). The finite-element based code, PHASTA, is utilized to perform large-scale simulations of two-phase bubbly flows in LWR geometries. Given the large computational cost of direct numerical simulation (DNS) coupled with interface tracking methods (ITM), typical domains encompass a portion of a single subchannel. In the presented research, the state-of-the-art analysis of turbulent two-phase flows in complex LWR subchannel geometries are demonstrated at both prototypical reactor parameters as well as scaled low pressure conditions. Three different cases are studied, a high-pressure simulation in prototypical reactor subchannel geometry, a low-pressure case in prototypical geometry and a final low-pressure case in a geometry scaled up to conserve the ratio between the bubble size and the domain pitch. Utilizing advanced statistical processing tools, these simulation conditions are compared to shed light on the relevancy of two-phase flow characteristics given the significant differences between LWR and low-pressure conditions. These findings can lead to the generation of useful guiding principles when researchers need to scale the two-phase flow behavior captured at low pressure and temperature conditions to those at reactor operating conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Full potential integral solution for transonic flows with and without embedded Euler domains

Two methods are presented to solve for the transonic airfoil flow problems. The first method is based on the integral equation solution of the full-potential equation in terms of the velocity field, and a Shock Capturing-Shock Fitting (SCSF) scheme has been developed. The SCSF-scheme consists of a shock-capturing part and a shock-fitting part in which shock panels are introduced at the shock location. The sock panels are fitted and crossed by using the Rankine-Hugoniot relations. The second method is based on coupling the integral equation of the full-potential equation with the pseudo time integration of Euler equations in a small embedded region around the shock. The integral solution provides the initial and boundary conditions for the Euler domain. This scheme is named as the Integral Equation-Embedded Euler (IEEE) scheme. The two methods are applied to NACA 0012 and NACA 64A010A over a wide range of Mach numbers, and the reults are in good agreement with the experimental data and other computational results. The schemes converge within a number of iterations which is one-order of magnitude less than the finite-difference schemes.

Kandil, Osama A.↗

Hexagonal Geometries in MPACT

The MPACT code is a high-fidelity light-water reactor analysis code using whole-core pin-resolved neutron transport calculations on modern parallel-computing hardware. MPACT uses the 2D/1D method to solve 3D neutron transport problems by decomposing the problem into a stack of 2D slices, each of which is solved independently using the method of characteristics (MOC). The slices are then coupled axially using the P3 nodal expansion method (NEM-P3) for the 1D axial calculations. MPACT also employs the coarse mesh finite difference (CMFD) method to accelerate calculations. This manuscript details work supporting advanced reactor designs using hexagonal pins and hexagonal assemblies such as the VVER-1000. If performed correctly, MOC is geometry agnostic. However, MPACT previously had optimizations in place for Cartesian geometries, specifically in the modularization and current calculations. Sections 2 and 3 detail the changes made to MPACT to support MOC and CMFD calculations on hexagonal geometries. Section 4 reports results demonstrating solution consistency for problems run with and without CMFD acceleration, results demonstrating solution consistency when run in serial and parallel, and pincell results using the Monte Carlo code, McCard’s benchmark results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Islet: interpolation semi-Lagrangian element-based transport

Abstract. Advection of trace species, or tracers, also called tracer transport, in models of the atmosphere and other physical domains is an important and potentially computationally expensive part of a model's dynamical core. Semi-Lagrangian (SL) advection methods are efficient because they permit a time step much larger than the advective stability limit for explicit Eulerian methods without requiring the solution of a globally coupled system of equations as implicit Eulerian methods do. Thus, to reduce the computational expense of tracer transport, dynamical cores often use SL methods to advect tracers. The class of interpolation semi-Lagrangian (ISL) methods contains potentially extremely efficient SL methods. We describe a finite-element ISL transport method that we call the interpolation semi-Lagrangian element-based transport (Islet) method, such as for use with atmosphere models discretized using the spectral element method. The Islet method uses three grids that share an element grid: a dynamics grid supporting, for example, the Gauss–Legendre–Lobatto basis of degree three; a physics parameterizations grid with a configurable number of finite-volume subcells per element; and a tracer grid supporting use of Islet bases with particular basis again configurable. This method provides extremely accurate tracer transport and excellent diagnostic values in a number of verification problems.

58 GEOSCIENCES↗

Three-dimensional elastic stress and displacement analysis of finite circular geometry solids containing cracks

A seminumerical method is presented for solving a set of coupled partial differential equations subject to mixed and coupled boundary conditions. The use of this method is illustrated by obtaining solutions for two circular geometry and mixed boundary value problems in three-dimensional elasticity. Stress and displacement distributions are calculated in an axisymmetric, circular bar of finite dimensions containing a penny-shaped crack. Approximate results for an annular plate containing internal surface cracks are also presented.

Gyekenyesi, J. P.↗

Transonic CFD applications at Boeing

The use of computational methods for three dimensional transonic flow design and analysis at the Boeing Company is presented. A range of computational tools consisting of production tools for every day use by project engineers, expert user tools for special applications by computational researchers, and an emerging tool which may see considerable use in the near future are described. These methods include full potential and Euler solvers, some coupled to three dimensional boundary layer analysis methods, for transonic flow analysis about nacelle, wing-body, wing-body-strut-nacelle, and complete aircraft configurations. As the examples presented show, such a toolbox of codes is necessary for the variety of applications typical of an industrial environment. Such a toolbox of codes makes possible aerodynamic advances not previously achievable in a timely manner, if at all.

Tinoco, E. N.↗

Efficient Viscous Design of Realistic Aircraft Configurations

This paper addresses the use of the Constrained Direct Iterative Surface Curvature (CDISC) design method in the aircraft design process. A discussion of some of the requirements for practical use of CFD in the design process is followed by a description of different CFD design methods, along with their relative strengths and weaknesses. A detailed description of the CDISC design method highlights some of the aspects of the method that provide computational efficiency and portability, as well as the flow and geometry constraint capabilities. In addition, an efficient approach to multipoint design, the Weighted Averaging of Geometries (WAG) method, is described and illustrated using a couple of simple examples. The CDISC and WAG methods are then applied to a complex generic business jet geometry using an unstructured grid flow solver to demonstrate the multipoint and multicomponent design capabilities of these methods. Introduction

Campbell, Richard L.↗

Development and Performance Simulations of the MPACT-AGREE Code Coupling Interface for MAGNOX Reactors

MPACT is a whole-core 3D neutron transport code jointly developed by Oak Ridge National Laboratory and the University of Michigan to perform high-fidelity light water reactor (LWR) analysis. Its capability was recently extended for gas-cooled reactor analysis by coupling to the thermo-fluids code AGREE for thermal feedback. This paper presents the MPACT-AGREE code coupling for the simulation of MAGNOX-type gas-cooled graphite-moderated reactors. A coupling interface has been developed, along with methods to accurately simulate MAGNOX type reactors. The coupling mechanics were tested using two problems derived from the Calder Hall reactor. Simulations of these models demonstrated that the codes have been successfully coupled and can provide reasonable results in a tractable simulation time.

Gentry, Cole↗

Preliminary Monte Carlo and Thermal Hydraulic Analysis using a Hybrid ETF-Corrected-Diffusion Prediction Block

This paper builds upon previous work to accelerate the Picard iteration (PI) method typically applied for coupled Monte Carlo-Thermal hydraulic (MC-TH) solutions. Previously, the use of the generalized transfer functions (GTFs) to predict variation in macroscopic cross sections following a perturbation in TH properties was demonstrated for a subset of simple 3D problems. In addition, the reduced-order transport prediction block relied on the first order perturbation (FOP) method, which was shown to have computational overheads. Recent work replaced the FOP block with a 1-group nodal diffusion solver to eliminate these overheads. While the use of diffusion is desirable for large-scale problems, the new solver introduces significant homogenization error. This work aims to address this issue by using the Jacobian-Free Newton Krylov (JFNK) method to generate a set of super homogenization (SPH) factors to improve the accuracy of the diffusion solution. The SPH factors will be used in conjunction with an improved cross section prediction method – the expanded transfer function (ETF) method – to produce a highly accurate flux prediction for an axial 1D boiling water reactor (BWR) pincell following a large perturbation in moderator density. The ETF-corrected diffusion (ETF-CD) block is shown to be highly accurate for the 1D test case. Future work will investigate the accuracy of the method for a realistic 3D pressurized water reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A simple improved low temperature correction for the hierarchical equations of motion

The study of open system quantum dynamics has been transformed by the hierarchical equations of motion (HEOM) method, which gives the exact dynamics for a system coupled to a harmonic bath at arbitrary temperature and system–bath coupling strength. However, in its standard form, this method is only consistent with the weak-coupling quantum master equation at all temperatures when many auxiliary density operators are included in the hierarchy, even when low temperature corrections are included. Here, we propose a new low temperature correction scheme for the termination of the hierarchy based on Zwanzig projection, which alleviates this problem and restores consistency with the weak-coupling master equation with a minimal hierarchy. The utility of the new correction scheme is demonstrated on a range of model systems, including the Fenna–Matthews–Olson complex. The new closure is found to improve convergence of the HEOM even beyond the weak-coupling limit and is very straightforward to implement in existing HEOM codes.

Fay, Thomas P. (ORCID:000000030625731X)↗

Anderson acceleration stability in NDA-accelerated k-eigenvalue problems

Anderson acceleration (AA) has been used to improve the stability and convergence rate of multiphysics iterative methods for reactor analysis. Most applications studied assume a tightly converged solution for the different physics problems, and AA is usually applied to state variables like temperature, density, and heat generation rate. In this paper, we study the theoretical performance of AA in NDA-accelerated k-eigenvalue problems. The problems and algorithms studied are simplified from the coupled iteration scheme adopted by MPACT and many other high-fidelity whole-core reactor codes. Compared to previous analyses of AA for these iteration schemes, we study the case with a partially converged neutronics solution and possibly partially converged nonlinear diffusion acceleration (NDA)/coarse mesh finite difference (CMFD) solutions. We observe that the performance of the iteration scheme with AA is very sensitive to the initial guess and is affected by the partially converged CMFD solutions. When the NDA solution is fully converged, using AA cannot achieve the optimal convergence rate in large-sized problems. Conversely, if the NDA solution is partially converged, the iteration scheme with AA can diverge or converge extremely slowly. It is found that the loss of robustness for AA is due to the fact that it is applied to the iterative subspace of state variables rather than the fundamental unknowns of the governing equations. To improve the robustness, the scalar flux should also be considered in the implementation of AA. After considering the residuals of flux, we observe that the stability is regardless of the partial convergence of NDA solutions. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Technique for Radiometer and Antenna Array Calibration - TRAAC

Highly sensitive receivers are used to detect minute amounts of emitted electromagnetic energy. Calibration of these receivers is vital to the accuracy of the measurements. Traditional calibration techniques depend on calibration reference internal to the receivers as reference for the calibration of the observed electromagnetic energy. Such methods can only calibrate errors in measurement introduced by the receiver only. The disadvantage of these existing methods is that they cannot account for errors introduced by devices, such as antennas, used for capturing electromagnetic radiation. This severely limits the types of antennas that can be used to make measurements with a high degree of accuracy. Complex antenna systems, such as electronically steerable antennas (also known as phased arrays), while offering potentially significant advantages, suffer from a lack of a reliable and accurate calibration technique. The proximity of antenna elements in an array results in interaction between the electromagnetic fields radiated (or received) by the individual elements. This phenomenon is called mutual coupling. The new calibration method uses a known noise source as a calibration load to determine the instantaneous characteristics of the antenna. The noise source is emitted from one element of the antenna array and received by all the other elements due to mutual coupling. This received noise is used as a calibration standard to monitor the stability of the antenna electronics.

Meyer, Paul↗

NEAMS Milestone Report: M2MS–20OR030102 FW–CADIS PWR Ex-Core Analysis with Shift Through VERA

This report presents the work completed for the NEAMS milestone M2MS–20OR030102 titled "FW–CADIS PWR Ex-Core Analysis with Shift through VERA." The work completed for this milestone includes the implementation, integration, and optimization of memory and performance improvement methods in Shift for fully coupled ex-core calculations through Virtual Environment for Reactor Applications (VERA). Fully coupled in this context means the transfer of moderator boron concentration, pin-wise fission source, depleted compositions, temperatures, and moderator densities from MPACT (with COBRA-TF (CTF)) to Shift. The ability to run ex-core calculations with VERA has been enabled and used for several years by Consortium for Advanced Simulation of Light Water Reactors (CASL) partners. However, this implementation was limited and potentially computationally burdensome. This work has enabled the ability to run higher-fidelity ex-core calculations on moderate computing clusters by focusing on multithreading, domain decomposition, and Forward-Weighted CADIS (FW-CADIS) variance reduction. Tests performed on small cores, a small modular reactor (SMR), and CASL progression problems show very promising memory reduction and computational performance. Recommendations for settings when running fully coupled high-fidelity ex-core calculations with VERA are documented. Without these optimization methods, many processors on a compute node would be left unused for the entire ex-core calculation. Therefore, these methods enable the user to better use the resources available and reduce computation time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A robust direct-integration method for rotorcraft maneuver and periodic response

The Newmark-Beta method and the Newton-Raphson iteration scheme are combined to develop a direct-integration method for evaluating the maneuver and periodic-response expressions for rotorcraft. The method requires the generation of Jacobians and includes higher derivatives in the formulation of the geometric stiffness matrix to enhance the convergence of the system. The method leads to effective convergence with nonlinear structural dynamics and aerodynamic terms. Singularities in the matrices can be addressed with the method as they arise from a Lagrange multiplier approach for coupling equations with nonlinear constraints. The method is also shown to be general enough to handle singularities from quasisteady control-system models. The method is shown to be more general and robust than the similar 2GCHAS method for analyzing rotorcraft dynamics.

Panda, Brahmananda↗

Micromechanics-Based Progressive Failure Analysis of Composite Laminates Using Different Constituent Failure Theories

Predicting failure in a composite can be done with ply level mechanisms and/or micro level mechanisms. This paper uses the Generalized Method of Cells and High-Fidelity Generalized Method of Cells micromechanics theories, coupled with classical lamination theory, as implemented within NASA's Micromechanics Analysis Code with Generalized Method of Cells. The code is able to implement different failure theories on the level of both the fiber and the matrix constituents within a laminate. A comparison is made among maximum stress, maximum strain, Tsai-Hill, and Tsai-Wu failure theories. To verify the failure theories the Worldwide Failure Exercise (WWFE) experiments have been used. The WWFE is a comprehensive study that covers a wide range of polymer matrix composite laminates. The numerical results indicate good correlation with the experimental results for most of the composite layups, but also point to the need for more accurate resin damage progression models.

Moncada, Albert M.↗

Intercomparison of Evapotranspiration Measurement Methods for Vegetable Crops in California

Recent drought events in California and legislation passed with the goal of increasing the sustainability of groundwater supplies have led to increased interest in tools to optimize irrigation schedules and increase on-farm water used efficiency. With more than 400 different crops produced in California, evapotranspiration-based irrigation scheduling is a promising and well-established approach. However, there is a need for accurate methods to estimate crop evapotranspiration (ET(sub c)) across the diverse range of crops grown, coupled with cost-effective methods for quantifying the accuracy of these tools. In this study, we evaluated remotely sensed estimates of ET(sub c) and associated crop water requirements from NASA's Satellite Irrigation Support (SIMS) system for two vegetable crops and measured crop evapotranspiration ET(sub c) using multiple methods, including weighing lysimeters, eddy covariance towers (EC), and surface renewal stations. We compared ET(sub c) data from these measurements with remotely sensed basal crop evapotranspiration (ET(sub cb)) data from SIMS as well as ET(sub c) data from a standard FAO-56 crop coefficient approach. Studies were conducted for sugar beets in Five Points, CA from 2014 to 2015 and studies are ongoing for fresh market tomatoes in Firebaugh, CA. We present results from these intercomparison studies and describe implications for future studies to quantify the accuracy of remotely sensed measures of ET(sub c). Highlights from results to date include strong correlations between ET measured with both surface renewal instrumentation and eddy covariance calculations using a 3D sonic anemometer and ET(sub c) data measured with the weighing lysimeter, with respective R2 values of 0.7964 (surface renewal) and 0.8034 (eddy covariance). This study provides insights into agreement between different approaches for monitoring evapotranspiration and provides another reference point for the community working to develop accurate and cost-effective tools that support growers in optimizing irrigation management.

Measurement↗