Modeling and Project Development
In order to reduce the costs of projects, it is necessary to make full use of today's modern simulation tools.
SEARCH · Engineering Papers
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.
In order to reduce the costs of projects, it is necessary to make full use of today's modern simulation tools.
The essential dynamical characteristics of a simple single spool turbojet engine were determined through simulation of low order system models on an analog computer. An accurate model was studied and system complexity was reduced through various linearizations and approximations. A derivation of a seventh order simplified simulation model is presented with a derivation of an even simpler third order model, and simulation results from each. The control problem studied is one of getting from zero fuel flow equilibrium to a high thrust equilibrium while taking into account surge margin and turbine inlet temperature constraints.
Large fission gas atoms such as Xe have very low solubility in the UO{sub 2} fuel matrix, which drives precipitation and growth of intra-granular fission gas bubbles, segregation of gas atoms to grain boundaries and evolution of intergranular bubbles. The inter-granular bubbles develop percolation networks leading to release of the fission gas. Fission gas bubbles induce fuel swelling and the fuel thermal conductivity is reduced by both fission gas bubbles and dispersed gas atoms. When fission gas is released from the fuel to the plenum, the pressure on the clad walls increases and the gap thermal conductivity decreases. All of these detrimentally affect fuel performance. The thermal conductivity of UO{sub 2} is low across a wide temperature range compared to iso-structural materials such as ThO{sub 2}, which is caused by the coupling between phonons carrying heat and uranium spins or magnetism providing a scattering mechanism with very low mean free path. For in-pile conditions the thermal conductivity is further reduced due to, for example, radiation damage, generation of fission products and fission gases and fuel microstructure evolution. This results in a complex burn-up dependent thermal conductivity. Thermal conductivity and fission gas evolution are two of the most important properties determining UO{sub 2} nuclear fuel performance. Both properties have significant uncertainty. In order to reduce this uncertainty, new models are being developed based on multi-scale simulations that attempt to improve the description of material properties and also couple them to the microstructure state of the material. In this paper we highlight atomistic calculations and simulations used to develop meso-scale models of the evolution of fission gases and thermal conductivity in UO{sub 2}. Literature examples of the application of these results and the corresponding models to fuel performance simulations will also be given. (author)
We investigate the effects of nonequilibrium dynamics in small colliding systems by comparing a nonequilibrium transport approach, Parton-Hadron-String Dynamics (PHSD), with a (2+1)-dimensional viscous hydrodynamic model, VISHNew. Focusing on p + Pb collisions at energies available at the Large Hadron Collider, we extract the initial conditions for the hydrodynamic model from PHSD, in order to reduce the impact of the early out-of-equilibrium dynamics and focus on the traces of nonequilibrium in the ensuing medium evolution. We find that in the transport approach quantities like energy density and bulk viscous pressure are highly inhomogeneous on the transverse plane during the whole evolution, whereas the hydrodynamic simulations dissolve more efficiently the initial spatial irregularities, still keeping a high degree of inhomogeneity due to the smaller size and lifetime of the medium produced in p + Pb collisions with respect to heavy-ion reactions. As a first step that will help to identify the impact of these nonequilibrium effects on final observables in proton-nucleus collisions, we perform an analysis of the transverse spherocity, an event-shape observable able to distinguish between jetty and isotropic configurations of transverse momenta. We found that the spherocity distribution in PHSD is slightly shifted towards the isotropic limit with respect to the hydrodynamic result. Even though this dissimilarity is partially due to the difference in the final charged particle production, it mainly comes from the different description within the two frameworks of the medium produced in small colliding systems. This finding supports the idea that multidifferential measurements, such as those based on event categorization according to multiplicity and spherocity, are useful to study final-state observables in ultrarelativistic proton-nucleus collisions.
Interconnected power systems have experienced a significant increase in its size and complexity, due to continuous topology and operating condition changes. It is computationally burdensome to represent entire system models in detail to conduct power system analysis. In the large-scale power system dynamic simulations, the models in the study area need to be retained in detail while the external area can be reduced without compromising accuracy, which can preserve the full system dynamic properties. This paper develops measurement-based continuous-time dynamic equivalents for power systems in order to both improve the model accuracy of the reduced system and increase the speed of dynamic simulation. Our method is based on a set of measurements at the boundary nodes between the study area and external area. Case studies demonstrate that the measurement-based model can capture the full system behavior sufficiently, and improve computational efficiency.Interconnected power systems experienced a significant increase in size and complexity. It is computationally burdensome to represent the entire system in detail to conduct power system analysis. Therefore, the model of the study system must be retained in detail while the external system can be reduced using system reduction techniques. This paper proposes a measurement-based dynamic equivalent in order to increase both model accuracy and simulation speed. The proposed method uses a set of measurements at the boundary nodes between the study area and external area for model parameter identification. Case studies demonstrate that the measurement-based technique can capture the main system behaviors accurately and improve computational efficiency.
Nuclear reactor digital twins (DTs) have been proposed for use as a safeguards technology to efficiently monitor new and novel reactors as they come online. A safeguards DT needs to be capable of detecting misuse and diversion as they occur, requiring physics models to be accurate and efficient. Mathematical surrogate models are capable of achieving the necessary efficiency and can largely maintain the accuracy of higher-order models given a quality training sample. The Multiphysics Object-Oriented Simulation Environment (MOOSE) code framework is specifically equipped to generate training samples and create surrogate models using full-order reactor physics models. Utilizing an operational AGN-201M reactor’s specifications, two surrogate types were trained on samples of variable size, and using Cartesian products, Latin hypercube sampling, and quadrature sampling, each was compared and evaluated on accuracy when compared to a full-order Monte Carlo model. Both surrogate types were able to capture reactivity changes within 0.05 $ of the Monte Carlo model while reducing the computation costs by eight orders of magnitude.
Comparison of hydrodynamic calculations with experimental data inevitably requires a model for converting the fluid to particles. For shear viscous fluids, several different conversion models have been proposed, such as the quadratic Grad corrections, the Strickland-Romatschke (SR) ansatz, self-consistent shear corrections from linearized kinetic theory, or corrections from the relaxation time approach. We show, via comparison to nonlinear 2 → 2 kinetic theory evolution, that even for low specific shear viscosities 0.03 ≲ $\textit{η/s}$ ≲ 0.2 , these four models have varying accuracy and still significant errors in reproducing particle phase space densities from hydrodynamic fields. Moreover, we demonstrate that the overall reconstruction error of additive shear viscous $f = f_{eq} + δf$ models dramatically reduces by up to one order of magnitude, if one ensures through exponentiation that $\textit{f}$ is always positive. Further, we also illustrate how even more accurate viscous $\textit{δf}$ models can be constructed by incorporating the first time derivatives of hydrodynamic fields to account for the past evolution of the system. Such time derivatives are readily available in hydrodynamic simulations, though usually not included in the output. Though our comparisons are limited to a one-component massless system undergoing a (0 + 1)-dimensional (0 + 1D) longitudinal boost-invariant expansion, we expect that the improved models will be useful in 2 + 1D and 3 + 1D hydrodynamic studies as well.
The uncertainties associated with neutrino-nucleus cross section modelling need to be reduced in order to measure neutrino oscillation parameters, the mass hierarchy and CP violation to a high degree of precision. SBND, a 112 ton liquid argon time projection chamber, will be the near detector of the short-baseline neutrino program at Fermilab. Once data taking begins in 2022, it will provide flux constraints for sterile neutrino searches and produce world leading neutrino-argon cross sections with up to seven million neutrino events in three years. This thesis describes the simulation, calibration and reconstruction tools developed to measure charged current muon neutrino interactions.Backgrounds from cosmic ray muons, neutrino interactions outside of the detector and different neutrino interaction channels were identified and removed as part of this process.The signal events were selected with purity of 90% and an efficiency of 67%, providing a strong starting point for meas urements of exclusive channels. The expected statistical and systematic uncertainties on rate predictions and cross section measurements at the full 6.6×10$^{20}$ proton on target exposure were evaluated.
The problem of constant gain output feedback regulator design for linear systems with ill-conditioned dynamics is considered in the context of singular perturbation theory. A design approach is developed in which gains can be separately calculated to stabilize reduced-order slow and fast subsystem models. By employing the notion of combined control and observation spillover suppression, conditions are derived assuring that these gains will stabilize the full-order system, assuming sufficient frequency separation between the slow and fast subsystems. An LQ design procedure is described in which the spillover suppression conditions are satisfied by adjoining penalty functions to the subsystem performance indices. The theory is demonstrated in a controller design for a flexible space structure.
Improved method for placement of eigenvalues and eigenvectors of closed-loop control system by use of either state or output feedback. Applied to reduced-order finite-element mathematical model of NASA's MAST truss beam structure. Model represents deployer/retractor assembly, inertial properties of Space Shuttle, and rigid platforms for allocation of sensors and actuators. Algorithm formulated in real arithmetic for efficient implementation. Choice of open-loop eigenvector matrix and its closest unitary matrix believed suitable for generating well-conditioned eigensystem with small control gains. Implication of this approach is that element of iterative search for "optimal" unitary matrix appears unnecessary in practice for many test problems.
This paper describes the formulation and validation of a high-order linearized mathematical model of helicopter flight mechanics, which includes rotor flap and lag degrees of freedom as well as inflow dynamics. The model is extracted numerically from an existing nonlinear, blade element, time simulation model. Extensive modifications in the formulation and solution process of the nonlinear model, required for a theoretically rigorous linearization, are described in detail. The validation results show that the linearized model successfully captures the coupled rotor-fuselage dynamics in the frequency band most critical for the design of advanced flight control systems. Additional results quantify the extent to which the order of the model can be reduced without loss of fidelity.
This paper demonstrates a new system identification approach of using Lanczos coordinates in place of modal coordinates. Identified experimental Lanczos vectors can be directly used in many structural dynamics analysis applications. A multi-input, multi-output frequency-domain technique was used to extract system matrices and an unsymmetric block Lanczos algorithm was used to reduce the order of the experimental model. A cantilever beam example showed promising results, indicating that a new system identification approach using Lanczos coordinates is worthy of further study.
This paper describes the formulation and validation of a high-order linearized mathematical model of helicopter flight mechanics, which includes rotor flap and lag degrees of freedom as well as inflow dynamics. The model is extracted numerically from an existing nonlinear, blade element, real-time simulation model. Extensive modifications in the formulation and solution process of the nonlinear model, required for a theoetically rigorous linearization, are described in detail. The validation results show that the linearized model successfully captures the coupled rotor-fuselage dynamics in the frequency band most critical for the design of advanced flight control systems. Additional results quantify the extent to which the order of the model can be reduced without loss of fidelity.
The Jet Propulsion Laboratory has learned that even innocuous errors in the spacecraft command process can have significantly detrimental effects on a space mission. Consequently, such Command File Errors (CFE), regardless of their effect on the spacecraft, are treated as significant events for which a root cause is identified and corrected. A CFE during space mission operations is often the symptom of imbalance or inadequacy within the system that encompasses the hardware and software used for command generation as well as the human experts and processes involved in this endeavor. As we move into an era of increased collaboration with other NASA centers and commercial partners, these systems become more and more complex. Consequently, the ability to thoroughly model and analyze CFEs formally in order to reduce the risk they pose is increasingly important. In this paper, we summarize the results of applying modeling techniques previously developed to the DAWN flight project. The original models were built with the input of subject matter experts from several flight projects. We have now customized these models to address specific questions for the DAWN flight project and formulating use cases to address their unique mission needs. The goal of this effort is to enhance the project's ability to meet commanding reliability requirements for operations and to assist them in managing their Command File Errors.
In this work, we study an effective spin model derived perturbatively from a random transverse-field Ising model on the pyrochlore lattice. The model consists of spin configurations on the pyrochlore lattice, restricted to the spin-ice subspace, with spins interacting with random Ising exchange couplings as well as ring exchanges along the hexagons of the lattice. This model is studied by exact diagonalization up to N=64 site systems. Furthermore, we calculate spin-glass correlation functions and local entanglement entropy S T between spins in a single tetrahedron and the rest of the system. We find that the model undergoes two phase transitions. At weak randomness the model is in a quantum spin-ice phase where S T =ln 6. Increasing randomness, at low transverse fields, first leads to a frozen phase, with long-range spin-glass order and S T =ln 2 corresponding to the Cat states associated with Ising order. Further increase in randomness leads to a random resonating-hexagon phase with a frozen backbone of spins and a broad distribution of entanglement entropies. The implications of these studies for non-Kramers rare-earth pyrochlores are discussed.
Abstract Radiation‐induced photocurrent in semiconductor devices can be simulated using complex physics‐based models, which are accurate, but computationally expensive. This presents a challenge for implementing device characteristics in high‐level circuit simulations where it is computationally infeasible to evaluate detailed models for multiple individual circuit elements. In this work we demonstrate a procedure for learning compact delayed photocurrent models that are efficient enough to implement in large‐scale circuit simulations, but remain faithful to the underlying physics. Our approach utilizes dynamic mode decomposition (DMD), a system identification technique for learning reduced‐order discrete‐time dynamical systems from time series data based on singular value decomposition. To obtain physics‐aware device models, we simulate the excess carrier density induced by radiation pulses by solving numerically the ambipolar diffusion equation, then use the simulated internal state as training data for the DMD algorithm. Our results show that the significantly reduced‐order delayed photocurrent models obtained via this method accurately approximate the dynamics of the internal excess carrier density—which can be used to calculate the induced current at the device boundaries—while remaining compact enough to incorporate into larger circuit simulations.
The problem of applying modular attitude control to a rigid body - flexible suspension model of a flexible spacecraft with some state variables inaccessible was addressed by developing a sequence of single axis models and generating a series of reduced state linear observers of minimum order to reconstruct those scalar variables that were inaccessible. The specific single axis models treated consisted of two, three and four rigid bodies, respectively, interconnected by a flexible shaft passing through the mass centers of the bodies. Reduced state linear observers of all orders up to one less than the total number of scalar state variables were generated for each of the three single axis models cited.
The problem of applying modular attitude control to a rigid body-flexible suspension model of a flexible spacecraft with some state variables inaccessible was addressed by developing a sequence of single axis models and generating a series of reduced state linear observers of minimum order to reconstruct those scalar state variables that were inaccessible. The specific single axis models treated consisted of two, three and four rigid bodies, respectively, interconnected by a flexible shaft passing through the mass centers of the bodies. Reduced state linear observers of all orders up to one less than the total number of scalar state variables were generated for each of the three single axis models cited.