Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “formulation systems”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Assessing Component Effects in Formulation Systems

The purpose of Snee (2011) was to discuss and illustrate how using Six Sigma concepts and approaches can lead to a better understanding of formulation systems. Central to the recommended approach is estimating and interpreting the effects of the components in the formulation on the response variable. Knowledge of components that have no effect, small or large effects (positive or negative), and similar effects is very helpful in understanding how the formulation system works. This understanding helps one decide 1) which components having large effects should be increased or decreased, and 2) what to do about components with no effects or similar effects. Correspondence between the authors following Snee (2011), along with further review of the literature, have pointed out that the approaches discussed in Snee (2011), and the interpretation of component effects in general, need further clarification. That is the focus of this article.

component effect directions↗
DOE OSTI2018

Hamiltonian formulation of systems with balanced loss–gain and exactly solvable models

A Hamiltonian formulation of generic many-body systems with balanced loss and gain is presented. It is shown that a Hamiltonian formulation is possible only if the balancing of loss and gain terms occurs in a pairwise fashion. It is also shown that with the choice of a suitable co-ordinate, the Hamiltonian can always be reformulated in the background of a pseudo-Euclidean metric. If the equations of motion of some of the well-known many-body systems like Calogero models are generalized to include balanced loss and gain, it appears that the same may not be amenable to a Hamiltonian formulation. A few exactly solvable systems with balanced loss and gain, along with a set of integrals of motion are constructed. The examples include a coupled chain of nonlinear oscillators and a many-particle Calogero-type model with four-body inverse square plus two-body pair-wise harmonic interactions. For the case of nonlinear oscillators, stable solution exists even if the loss and gain parameter has unbounded upper range. Further, the range of the parameter for which the stable solutions are obtained is independent of the total number of the oscillators. The set of coupled nonlinear equations are solved exactly for the case when the values of all the constants of motions except the Hamiltonian are equal to zero. Exact, analytical classical solutions are presented for all the examples considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Demonstration of emulator-based Bayesian calibration of safety analysis codes: Theory and formulation

System codes for simulation of safety performance of nuclear plants may contain parameters whose values are not known very accurately. New information from tests or operating experience is incorporated into safety codes by a process known as calibration, which reduces uncertainty in the output of the code and thereby improves its support for decision-making. The work reported here implements several improvements on classic calibration techniques afforded by modern analysis techniques. The key innovation has come from development of code surrogate model (or code emulator) construction and prediction algorithms. Use of a fast emulator makes the calibration processes used here with Markov Chain Monte Carlo (MCMC) sampling feasible. This study uses Gaussian Process (GP) based emulators, which have been used previously to emulate computer codes in the nuclear field. The present work describes the formulation of an emulator that incorporates GPs into a factor analysis-type or pattern recognition-type model. This “function factorization” Gaussian Process (FFGP) model allows overcoming limitations present in standard GP emulators, thereby improving both accuracy and speed of the emulator-based calibration process. Calibration of a friction-factor example using a Method of Manufactured Solution is performed to illustrate key properties of the FFGP based process.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

The integrated Earth system model version 1: formulation and functionality

The integrated Earth system model (iESM) has been developed as a new tool for projecting the joint human/climate system. The iESM is based upon coupling an integrated assessment model (IAM) and an Earth system model (ESM) into a common modeling infrastructure. IAMs are the primary tool for describing the human–Earth system, including the sources of global greenhouse gases (GHGs) and short-lived species (SLS), land use and land cover change (LULCC), and other resource-related drivers of anthropogenic climate change. ESMs are the primary scientific tools for examining the physical, chemical, and biogeochemical impacts of human-induced changes to the climate system. The iESM project integrates the economic and human-dimension modeling of an IAM and a fully coupled ESM within a single simulation system while maintaining the separability of each model if needed. Both IAM and ESM codes are developed and used by large communities and have been extensively applied in recent national and international climate assessments. By introducing heretofore-omitted feedbacks between natural and societal drivers, we can improve scientific understanding of the human–Earth system dynamics. Potential applications include studies of the interactions and feedbacks leading to the timing, scale, and geographic distribution of emissions trajectories and other human influences, corresponding climate effects, and the subsequent impacts of a changing climate on human and natural systems. This paper describes the formulation, requirements, implementation, testing, and resulting functionality of the first version of the iESM released to the global climate community.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Dirac open-quantum-system dynamics: Formulations and simulations

In this work, we present an open-system interaction formalism for the Dirac equation. Overcoming a complexity bottleneck of alternative formulations, our framework enables efficient numerical simulations (utilizing a typical desktop) of relativistic dynamics within the von Neumann density matrix and Wigner phase-space descriptions. Employing these instruments, we gain important insights into the effect of quantum dephasing for relativistic systems in many branches of physics. In particular, the conditions for robustness of Majorana spinors against dephasing are established. Using the Klein paradox and tunneling as examples, we show that quantum dephasing does not suppress negative energy particle generation. Hence, the Klein dynamics is also robust to dephasing.

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

Computation of Direct Sensitivities of Spatial Multibody Systems With Joint Friction

Abstract Friction exists in most mechanical systems and may have a major influence on the dynamic performance of the system. The incorporation of friction in dynamic systems has been a subject of active research for several years owing to its high nonlinearity and its dependence on several parameters. Consequently, optimization of dynamic systems with friction becomes a challenging task. Gradient-based optimization of dynamical systems is a prominent technique for optimal design and requires the computation of model sensitivities with respect to the design parameters. The novel contribution of this paper is the derivation of the analytical methodology for the computation of direct sensitivities for smooth multibody systems with joint friction using the Lagrangian index-1 formulation. System dynamics have been computed using two different friction models; the Brown and McPhee, and the Gonthier et al. model. The methodology proposed to obtain model sensitivities has also been validated using the complex finite difference method. A case study has been conducted on a spatial multibody system to observe the effect of friction on the dynamics and model sensitivities, compare sensitivities with respect to different parameters and demonstrate the numerical and validation aspects. Since design parameters can have very different magnitudes and units, the sensitivities have been scaled with the parameters for comparison. Finally, a discussion has been presented on the interpretation of the case study results. Due to the incorporation of joint friction, ‘jumps’ or discontinuities are observed in the model sensitivities akin to those observed for hybrid dynamical systems.

Engineering↗
DOE OSTI2019

Effect of chemical environment on the radiation chemistry of N , N -di-(2-ethylhexyl)butyramide (DEHBA) and plutonium retention

N,N-di-(2-ethylhexyl)butyramide (DEHBA) has been proposed as part of a hydro-reprocessing solvent extraction system for the co-extraction of uranium and plutonium from spent nuclear fuel, owing to its selectivity for hexavalent uranium and tetravalent plutonium. However, there is a critical lack of quantitative understanding regarding the impact of chemical environment on the radiation chemistry of DEHBA, and how this would affect process performance. Here we present a systematic investigation into the radiolytic degradation of DEHBA in a range of n-dodecane solvent system formulations, where we subject DEHBA to gamma irradiation, measure reaction kinetics, ligand integrity, degradation product formation, and investigate solvent system performance through uranium and plutonium extraction and strip distribution ratios. The rate of DEHBA degradation in n-dodecane was found to be slow (G = -0.31 ± 0.02 μmol J -1 ) but enhanced upon contact with the oxidizing conditions of the investigated solvent systems (organic-only, or in contact with either 0.1 or 3.0 M aqueous nitric acid). Two major degradation products were identified in the organic phase, bis-2-ethylhexylamine (b2EHA) and N-(2-ethylhexyl)butyramide (MEHBA), resulting from the cleavage of C–N bonds, and could account for the total loss of DEHBA up to ~300 kGy for organic-only conditions. Both b2EHA and MEHBA were also found to be susceptible to radiolytic degradation, having G-values of -0.12 ± 0.01 and -0.08 ± 0.01 μmol J -1 , respectively. Solvent extraction studies showed: (i) negligible change in uranium extraction and stripping with increasing absorbed dose; and (ii) plutonium extraction and retention exhibits complex dependencies on absorbed dose and chemical environment. Organic-only conditions afforded enhanced plutonium extraction and retention attributed to b2EHA, while acid contacts inhibited this effect and promoted significant plutonium retention for the highest acidity. Overall it has been demonstrated that chemical environment during irradiation has a significant influence on the extent of DEHBA degradation and plutonium retention.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Radiolytic transformation of AHA under single-cycle conditions

Historically, the radiation robustness of extractant ligands and the impact of their radiolysis products have been used to evaluate the radiolytic feasibility of solvent system formulations for used nuclear fuel (UNF) reprocessing. However, other key additives (e.g., reducing and hold-back reagents) are employed by these solvent systems. These chemical species have typically been overlooked from a radiolytic perspective, and consequently little is known on how they impact process performance. Therefore, to holistically support the development of advanced, simplified, single-cycle flowsheets for UNF reprocessing, the radiation robustness of two key proposed additives – (i) the complexing reductant acetohydroxamic acid (AHA) and (ii) the hold-back reagent 1,2-cyclohexylenedinitrilotetraacetic acid (CDTA) – have been evaluated under envisioned process conditions using a combination of steady-state gamma and time-resolved pulsed electron irradiation techniques.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2019

Open Data, Models, and Codes for Vanadium Redox Batch Cell Systems: A Systems Approach Using Zero-Dimensional Models

Abstract In this paper, we study, analyze, and validate some important zero-dimensional physics-based models for vanadium redox batch cell (VRBC) systems and formulate an adequate physics-based model that can predict the battery performance accurately. In the model formulation process, a systems approach to multiple parameters estimation has been conducted using VRBC systems at low C-rates (~C/30). In this batch cell system, the effect of ions' crossover through the membrane is dominant, and therefore, the capacity loss phenomena can be explicitly observed. Paradoxically, this means that using the batch system might be a better approach for identifying a more suitable model describing the effect of ions transport. Next, we propose an efficient systems approach, which enables to help understand the battery performance quickly by estimating all parameters of the battery system. Finally, open source codes, executable files, and experimental data are provided to enable people's access to robust and accurate models and optimizers. In battery simulations, different models and optimizers describing the same systems produce different values of the estimated parameters. Providing an open access platform can accelerate the process to arrive at robust models and optimizers by continuous modification from the users' side.

25 ENERGY STORAGE↗
DOE OSTI2023

A Robust Parallel Distributed State Estimation for Large Scale Distribution Systems

The growing need and interest in real-time monitoring of large distribution networks motivated by the rapid population of renewable sources, EVs and etc. demand a computationally efficient state estimation framework. Furthermore, this paper presents an improved computational framework for implementing a robust state estimator using a multi-core processor. The main contribution of the paper is the proposed computational framework along with two partitioning strategies which enable fast and robust state estimation for large scale radial and/or meshed distribution systems. Formulation of the proposed method and its implementation are described in detail. Performance of the estimator is tested by simulations first using a small 84-bus radial distribution system. Then the method’s scalability is demonstrated by simulations on two very large scale distribution networks one configured radially and the other meshed each containing over 12,500 buses.

42 ENGINEERING↗
DOE OSTI2023

Dioctyl ether radiolysis under used nuclear fuel reprocessing conditions: Foundational knowledge for the development of sacrificial ligand grafts

Radiation-induced ligand destruction and concomitant degradation product formation is unavoidable under envisioned used nuclear fuel reprocessing conditions, and ultimately limits the recyclability of a given solvent system formulation. With this in mind, a novel approach to ligand innovation for advanced separations processes is in the design of tailored “protomolecules,” which when exposed to ionizing radiation undergo desired chemical transformations, thereby allowing for a reprocessing solvent system to advantageously evolve with absorbed radiation dose. Here we provide foundational knowledge for the time-resolved (electron pulse) and steady-state (cobalt-60 gamma) radiolytic behavior of dioctyl ether, a protomolecule grafting surrogate. Gamma irradiation of single-phase solutions of dioctyl ether (5–100 vol.%) in n-dodecane resulted in the loss of dioctyl ether (G = -0.72 µmol J –1 ) and the formation of octane (G = 0.11 µmol J –1 ) and octanol (G = 0.15 µmol J –1 ) degradation products, the latter of which is a known reprocessing phase modifier and radioprotectant. Further, these radiation-induced changes were attributed to direct radiation effects, for more concentrated dioctyl ether solutions, and indirect radiation effects, predominantly driven by the reaction of dioctyl ether with the dodecane radical cation, for which we report a second-order rate coefficient of k = (1.53 ± 0.05) × 10 10 L mol –1 s –1 . Under typical biphasic extraction (4.0 mol L –1 HNO 3 ) and strip (0.1 mol L –1 HNO 3 ) reprocessing conditions, gamma irradiation of 5 vol.% dioctyl ether solvent systems afforded negligible change in the rate of parent molecule destruction but promoted significant differences in the radiolytic behavior of its degradation products. These differences are attributed to their respective interactions with [dioctyl ether/HNO 3 /H 2 O] and [octanol/HNO 3 /H 2 O] adducts extracted into the organic phase. Overall, these results support the grafting of ether linkages to advanced separations ligands (e.g., modified diglycolamides).

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2016

On dynamical systems approaches and methods in f ( R ) cosmology

We discuss dynamical systems approaches and methods applied to flat Robertson-Walker models in f ( R )-gravity. We argue that a complete description of the solution space of a model requires a global state space analysis that motivates globally covering state space adapted variables. This is shown explicitly by an illustrative example, f ( R ) = R + α R {sup 2}, α > 0, for which we introduce new regular dynamical systems on global compactly extended state spaces for the Jordan and Einstein frames. This example also allows us to illustrate several local and global dynamical systems techniques involving, e.g., blow ups of nilpotent fixed points, center manifold analysis, averaging, and use of monotone functions. As a result of applying dynamical systems methods to globally state space adapted dynamical systems formulations, we obtain pictures of the entire solution spaces in both the Jordan and the Einstein frames. This shows, e.g., that due to the domain of the conformal transformation between the Jordan and Einstein frames, not all the solutions in the Jordan frame are completely contained in the Einstein frame. We also make comparisons with previous dynamical systems approaches to f ( R ) cosmology and discuss their advantages and disadvantages.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

A new perspective on the integrability of Inozemtsev’s elliptic spin chain

The aim of this paper is studying from an alternative point of view the integrability of the spin chain with long-range elliptic interactions introduced by Inozemtsev. Our analysis relies on some well-established conjectures characterizing the chaotic vs. integrable behavior of a quantum system, formulated in terms of statistical properties of its spectrum. More precisely, we study the distribution of consecutive levels of the (unfolded) spectrum, the power spectrum of the spectral fluctuations, the average degeneracy, and the equivalence to a classical vertex model. Our results are consistent with the general consensus that this model is integrable, and that it is closer in this respect to the Heisenberg chain than to its trigonometric limit (the Haldane–Shastry chain). On the other hand, we present some numerical and analytical evidence showing that the level density of Inozemtsev’s chain is asymptotically Gaussian as the number of spins tends to infinity, as is the case with the Haldane–Shastry chain. We are also able to compute analytically the mean and the standard deviation of the spectrum, showing that their asymptotic behavior coincides with that of the Haldane–Shastry chain. - Highlights: • Construction of Inozemtsev’s elliptic spin chain using Polychronakos’s freezing trick. • Numerical evidence of the Gaussian character of the level density. • Exact computation and asymptotics of the mean and standard deviation of the spectrum. • Evidence of the chain’s integrability from key statistical properties of its spectrum. • Exact evaluation of finite sums of powers of Weierstrass’s elliptic function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Evaluation of the Index and Singular Points of Linear Differential-Algebraic Equations of Higher Order

We study linear systems of ordinary differential equations of higher order with an identically singular matrix at the higher derivative of the required vector-valued function in the domain of definition. We give definitions of the index and singular points for these systems, formulate the conditions of solvability, and deduce a formula for the general solution. The algorithms for finding the index and singular points are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A special generalized HSS method for discrete ill-posed problems

In this paper, we use a new special generalized Hermitian and skew-Hermitian splitting (SGHSS) method for solving ill-posed inverse problems. Based on an augmented system formulation, we apply a new splitting for the Hermitian part of the coefficient matrix and discuss on the convergence of this iterative method. Moreover, the optimal parameter minimizing the spectral radius of iteration matrix is derived. Finally, we present some numerical examples to demonstrate the performance of the presented method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

2016 High Performance Computing Annual Report

Sandia National Laboratories’ California site is celebrating its 60th anniversary (1956 to 2016), and high performance computing has been a key enabler for its scientists and engineers throughout much of its history. Since its founding, Sandia California has helped pioneer the use of HPC platforms including hosting Sandia’s first Cray-1 supercomputer in the 1970s and supporting development of scalable cluster computing platforms to create a new paradigm for cost-effective supercomputing in the 1990s. Recent decades of investment in creation of scalable application frameworks for scientific computing have also enabled new generations of modeling and simulation codes. These resources have facilitated computational analysis of complex phenomena in diverse applications spanning national defense, energy, and homeland security. Today, Sandia California researchers work with partners in academia, industry, and national labs to evolve the state-of-the-art in HPC, modeling, and data analysis (including foundational capabilities for exascale computing platforms) and apply them in transformational ways. Research efforts include mitigating the effects of silent hardware failures that can jeopardize the results of large-scale computations, developing exascale-capable asynchronous task-parallel programming models and runtime systems, formulating new techniques to better explore and analyze extreme-scale data sets, and increasing our understanding of materials and chemical sciences which has applications spanning nuclear weapons stockpile stewardship to more efficient automobile engines. The following section highlights some of these research and applications projects and further illustrates the breadth of our HPC capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Resilient Computing with Dynamical Systems

We reformulate fundamental numerical problems to run on novel hardware inspired by the brain. Such "neuromorphie hardware consumes less energy per computation, promising a means to augment next-generation exascale computers. However, their programming model is radically different from floating-point machines, with fewer guarantees about precision and communication. The approach is to pass each given problem through a sequence of transformations (algorithmic "reductions") which change it from conventional form into a dynamical system, then ultimately into a spiking neural network. Results for the eigenvalue problem are presented, showing that the dynamical system formulation is feasible.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Autoregressive Neural Network for Simulating Open Quantum Systems via a Probabilistic Formulation

The theory of open quantum systems lays the foundation for a substantial part of modern research in quantum science and engineering. Rooted in the dimensionality of their extended Hilbert spaces, the high computational complexity of simulating open quantum systems calls for the development of strategies to approximate their dynamics. In this Letter, we present an approach for tackling open quantum system dynamics. Using an exact probabilistic formulation of quantum physics based on positive operator-valued measure, we compactly represent quantum states with autoregressive neural networks; such networks bring significant algorithmic flexibility due to efficient exact sampling and tractable density. We further introduce the concept of string states to partially restore the symmetry of the autoregressive neural network and improve the description of local correlations. Efficient algorithms have been developed to simulate the dynamics of the Liouvillian superoperator using a forward-backward trapezoid method and find the steady state via a variational formulation. Our approach is benchmarked on prototypical one-dimensional and two-dimensional systems, finding results which closely track the exact solution and achieve higher accuracy than alternative approaches based on using Markov chain Monte Carlo method to sample restricted Boltzmann machines. Our Letter provides general methods for understanding quantum dynamics in various contexts, as well as techniques for solving high-dimensional probabilistic differential equations in classical setups.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗