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 1,279 records · Page 71

Robust Fixed-Structure Controller Synthesis

The ability to develop an integrated control system design methodology for robust high performance controllers satisfying multiple design criteria and real world hardware constraints constitutes a challenging task. The increasingly stringent performance specifications required for controlling such systems necessitates a trade-off between controller complexity and robustness. The principle challenge of the minimal complexity robust control design is to arrive at a tractable control design formulation in spite of the extreme complexity of such systems. Hence, design of minimal complexitY robust controllers for systems in the face of modeling errors has been a major preoccupation of system and control theorists and practitioners for the past several decades.

Corrado, Joseph R.↗

Ultra reliability at NASA

Ultra reliable systems are critical to NASA particularly as consideration is being given to extended lunar missions and manned missions to Mars. NASA has formulated a program designed to improve the reliability of NASA systems. The long term goal for the NASA ultra reliability is to ultimately improve NASA systems by an order of magnitude. The approach outlined in this presentation involves the steps used in developing a strategic plan to achieve the long term objective of ultra reliability. Consideration is given to: complex systems, hardware (including aircraft, aerospace craft and launch vehicles), software, human interactions, long life missions, infrastructure development, and cross cutting technologies. Several NASA-wide workshops have been held, identifying issues for reliability improvement and providing mitigation strategies for these issues. In addition to representation from all of the NASA centers, experts from government (NASA and non-NASA), universities and industry participated. Highlights of a strategic plan, which is being developed using the results from these workshops, will be presented.

risk↗

A Moving Discontinuous Galerkin Method with Interface Condition Enforcement for Reacting Hypersonic Flows

The necessity of enforcing conservation in computational elements or cells (element conservation) for discontinuous solutions is well understood and respected for solving conservation laws in computational fluid dynamics (CFD). In contrast, interface conservation, where the conservation across cell interfaces is enforced, is long ignored, and yet is also ruled and required by the underlying physics just like element conservation. Violation of the interface conservation across discontinuities is the root cause why an exact discontinuous solution can never be achieved in shock capturing methods. The interface conservation is examined and explored in this talk. Moving discontinuous Galerkin (MDG) finite element method with interface conservation enforcement (MDG-ICE) [1],[2] are then presented for solving compressible flow problems with discontinuities based on the observation that the interface conservation can only be satisfied, only when mesh interfaces are aligned with discontinuities. In the MDG-ICE formulation, both conservative quantities and grid geometry are considered as independent variables. A space-time DG formulation is used to solve the multi-material compressible Euler equations in the standard discontinuous solution space and the discrete grid geometry is solved using a variational formulation in a continuous space. A self-adaptive Levenberg-Marquardt method is utilized to solve the resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation. A number of numerical experiments for a variety of flow problems are conducted to assess the accuracy and performance of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to deliver the designed order of both h- and p-convergence even for discontinuous solutions, and detect all types of interfaces, via interface condition enforcement and satisfy, via grid movement, the compressible Euler equations and the associated interface condition.

Discontinuous Galekin↗

Robust Distributed Volt/VAR Control of Power Distribution Systems

Volt/VAR control (VVC) problem is important for reliable and efficient operation of modern power distribution systems. In this paper, we present a novel formulation of the VVC problem in the presence of uncertainties adopting the robust control perspective. To handle the structural constraints imposed on the controller, a set of sufficient conditions transforming a bilinear matrix inequality into a linear matrix inequality are derived. Simulation results using a real world power distribution systems are presented to illustrate the contribution of the proposed solution.

Chen, Hua↗

Extension of the Kohn-Sham formulation of density functional theory to finite temperature

Based on Mermin's extension of the Hohenberg and Kohn theorems to non-zero temperature, the Kohn-Sham formulation of density functional theory (KS-DFT) is generalized to finite temperature. Here, we show that present formulations are inconsistent with Mermin's functional containing expressions, in particular describing the Coulomb energy, that defy derivation and are even in violation of rules of logical inference. More; current methodology is in violation of fundamental laws of both quantum and classical mechanics. Based on this feature, we demonstrate the impossibility of extending the KS formalism to finite temperature through the self-consistent solutions of the single-particle Schrödinger equation of T>0. Guided by the form of Mermin's functional that depends on the eigenstates of a Hamiltonian, determined at T>0 we base our extension of KS-DFT on the determination of the excited states of a non-interacting system at the zero of temperature. The resulting formulation is consistent with that of Mermin constructing the free energy at T>0 in terms of the excited states of a non-interacting Hamiltonian (system) that, within the KS formalism, are described by Slater determinants. To determine the excited states at T=0 use is made of the extension of the Hohenberg and Kohn theorems to excited states presented in previous work applied here to a non-interacting collection of replicas of a non-interacting N-particle system, whose ground state density is taken to match that of K non-interacting replicas of an interacting N-particle system at T>0. The formalism allows for an ever denser population of the excitation spectrum of a Hamiltonian, within the KS approximation. The form of the auxiliary potential, (Kohn-Sham potential), is formally identical to that in the ground state formalism with the contribution of the Coulomb energy provided by the derivative of the Coulomb energy in all excited states taken into account. Once the excited states are determined, the minimum of the free energy within the KS formalism follows immediately in the form of Mermin's functional, but with the exact excited states in that functional represented by Slater determinants obtained through self-consistency conditions at the zero of temperature. Lastly, it is emphasized that, in departure from all existing formulations, no self-consistency conditions are implemented at finite T; as we show, in fact, such formulations are rigorously blocked.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Efficient discontinuous Galerkin implementations and preconditioners for implicit unsteady compressible flow simulations

This work presents and compares efficient implementations of high-order discontinuous Galerkin methods: a modal matrix-free discontinuous Galerkin (DG) method, a hybridizable discontinuous Galerkin (HDG) method, and a primal formulation of HDG, applied to the implicit solution of unsteady compressible flows. The matrix-free implementation allows for a reduction of the memory footprint of the solver when dealing with implicit time-accurate discretizations. HDG reduces the number of globally-coupled degrees of freedom relative to DG, at high order, by statically condensing element-interior degrees of freedom from the system in favor of face unknowns. The primal formulation further reduces the element-interior degrees of freedom by eliminating the gradient as a separate unknown. This paper introduces a p-multigrid preconditioner implementation for these discretizations and presents results for various flow problems. Benefits of the p-multigrid strategy relative to simpler, less expensive, preconditioners are observed for stiff systems, such as those arising from low-Mach number flows at high-order approximation. The p-multigrid preconditioner also shows excellent scalability for parallel computations. Finally, additional savings in both speed and memory occur with a matrix-free/reduced version of the preconditioner.

97 MATHEMATICS AND COMPUTING↗

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗

Development of ADOCS controllers and control laws. Volume 3: Simulation results and recommendations

The Advanced Cockpit Controls/Advanced Flight Control System (ACC/AFCS) study was conducted by the Boeing Vertol Company as part of the Army's Advanced Digital/Optical Control System (ADOCS) program. Specifically, the ACC/AFCS investigation was aimed at developing the flight control laws for the ADOCS demonstator aircraft which will provide satisfactory handling qualities for an attack helicopter mission. The three major elements of design considered are as follows: Pilot's integrated Side-Stick Controller (SSC) -- Number of axes controlled; force/displacement characteristics; ergonomic design. Stability and Control Augmentation System (SCAS)--Digital flight control laws for the various mission phases; SCAS mode switching logic. Pilot's Displays--For night/adverse weather conditions, the dynamics of the superimposed symbology presented to the pilot in a format similar to the Advanced Attack Helicopter (AAH) Pilot Night Vision System (PNVS) for each mission phase is a function of SCAS characteristics; display mode switching logic. Results of the five piloted simulations conducted at the Boeing Vertol and NASA-Ames simulation facilities are presented in Volume 3. Conclusions drawn from analysis of pilot rating data and commentary were used to formulate recommendations for the ADOCS demonstrator flight control system design. The ACC/AFCS simulation data also provide an extensive data base to aid the development of advanced flight control system design for future V/STOL aircraft.

Landis, Kenneth H.↗

On the manipulability of dual cooperative robots

The definition of manipulability ellipsoids for dual robot systems is given. A suitable kineto-static formulation for dual cooperative robots is adopted which allows for a global task space description of external and internal forces, and relative velocities. The well known concepts of force and velocity manipulability ellipsoids for a single robot are formally extended and the contributions of the two single robots to the cooperative system ellipsoids are illustrated. Duality properties are discussed. A practical case study is developed.

Chiacchio, P.↗

Novel method to control explosive shock sensitivity: Formulations, experiments and modeling

The microstructure of a heterogeneous high explosive (HE) affects its shock initiation sensitivity and detonation performance. Alteration of the void content and/or void structure (i.e., bulk heating or mechanical damage) therefore changes the shock initiation behavior. Controlling and predicting the change in shock sensitivity after an HE has undergone microstructural changes addresses an important and challenging goal for the design and understanding of novel energetic material formulations. In this study, we aimed to develop HE systems to precisely tune shock sensitivity by a thermal treatment prior to use. Specifically, we incorporated a small fraction (1 wt% or less) of thermally expandable microspheres (TEMs) during the formulation process of various plastic-bonded explosives (PBX). TEMs typically consist of a thermoplastic acrylonitrile shell (10-50 µm diameter), which encapsulates an inert low boiling hydrocarbon. Upon heating, the TEMs expand as the shell softens while the hydrocarbon gasifies, increasing the internal pressure and expanding the particle by as much as 120 vol% (irreversibly). Here, we present our progress on shock sensitivity comparisons of HE formulations doped with TEMs after heating and expansion. Experiments, using Composition C-4 thermally cycled to 120°C, showed an increase in shock sensitivity. Additionally, mesoscale modeling revealed that the TEM itself does not act as a hotspot, but instead has a secondary effect on the response of nearby voids. We conclude that TEMs indeed do provide a method to tune shock sensitivity. Further, our results emphasize the non-uniformity of shock waves at the mesoscale, and that upstream defects influence the reactivity of downstream defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Voltage Probability Density Function Shaping Control Strategy Considering Grid Operational Uncertainties

It is well-known that power systems operation always affected by various uncertainties which make the bus voltage a random process that can be characterized by its probability density function (PDF) at any time instant. In this context, this paper presents a novel PDF-based voltage control framework for power systems. By modeling voltage as a stochastic process, we formulate a stochastic differential equationthat captures grid uncertainties. The associated Fokker-Planck-Kolmogorov equation is derived to describe the evolution of the voltage PDF, which enables the formulation of a PDF-shaping control strategy. To simplify the PDF control formulation, a B-spline neural network is introduced for real-time estimation and regulation of the voltage distribution. The proposed PDF control law updates voltage references for energy storage systems and synchronous generators using real-time PDF measurements and feedback signals. The proposed method is validated on a modified Kundur’s two-area system. Simulation results demonstrate that the controller can significantly improve the voltage stability under stochastic conditions, highlighting its effectiveness in modern inverter-rich grids.

Gui, Yonghao [ORNL] (ORCID:0000000250435534)↗

Future applications of solid-state lasers for space remote sensing

This paper describes the Laser Atmospheric Sounder and Altimetry (LASA) and the Lidar In-Space Technology Experiment (LITE) projects, both of which will use solid-state lasers for remote sensing. LASA will provide long-term global measurements of aerosols and clouds in the earth's atmosphere. The LITE project is designed to demonstrate the first use of a lidar system in space, proving the technology required for future space-based lidar facilities such as LASA. A Shuttle flight experiment is currently being developed to enhance the technology and to formulate the design of the solid-state lidar system. Design diagrams of LITE and LASA facilities and their subsystems are presented.

Cox, John W.↗

On the bulk viscosity of relativistic matter

It is shown that the Kubo formulation leads to an expression for bulk viscosity involving the differences of the trace equilibrium and nonequilibrium hydrodynamic tensors just as in the Boltzmann formulation. It is argued that, if a dense system is represented by a quantum field, the field-theoretic energy-momentum tensor and, therefore, the hydrodynamic tensor should always have a vanishing trace in the limit of high energies. This bulk viscosity should vanish in the same limit. Finally, it is noted that the explicit use of well-known statistical mechanical formulation can be considered as a justification of Weinberg's intuitive argument from a fundamental point of view.

Canuto, V.↗

A multiscale anisotropic polymer network model coupled with phase field fracture

The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.

42 ENGINEERING↗

Steady-state radiochemical transport model of the molten salt reactor experiment

A radiochemical transport model of the Molten Salt Reactor Experiment (MSRE) was formulated based on the material derivative featuring system level 1-D flow providing spatial dependence accounting for multiple phenomena including radioactive decay, chemical removal, and fuel depletion across custom volume regions. Further, the steady-state version, which is based on a system of ODE's, can be solved efficiently and is useful for sensitivity studies utilizing Monte Carlo sampling. Accounting for multiple phenomena in a fully coupled model is necessary to understand the complexities of transport and distribution of radionuclides throughout the primary loop of a flowing fuel reactor.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Gauge-invariant Green function dynamics: A unified approach

We present a gauge-invariant description of Green function dynamics introduced by means of a generalized Peirels phase involving an arbitrary differentiable path in space–time. Two other approaches to formulating a gauge-invariant description of systems, the Green function treatment of Levanda and Fleurov [M. Levanda, V. Fleurov, J. Phys.: Condens. Matter 6 (1994) 7889] and the usual multipolar expansion for an atom, are shown to arise as special cases of our formalism. We argue that the consideration of paths in the generalized Peirels phase that do not lead to introduction of an effective gauge-invariant Hamiltonian with polarization and magnetization fields may prove useful for the treatment of the response of materials with short electron correlation lengths. -- Highlights: •Peirels phase for an arbitrary path in space–time established. •Gauge-invariant Green functions and the Power–Zienau–Wooley transformation connected. •Limitations on possible polarization and magnetization fields established.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Basis-independent quantum coherence and its distribution

We analyze a basis-independent definition of quantum coherence. The maximally mixed state is used as the reference state, which allows for a way of defining coherence that is invariant under arbitrary unitary transformations. Such a basis independent formulation is related to the purity of system rather than the superposition in a particular basis. The basis-independent approach is applied to finding the distribution of the coherence within a multipartite system, where the contributions due to correlations between the subsystems and within each subsystem are isolated. The use of the square root of the Jensen–Shannon divergence allows for inequality relations to be derived between these quantities, giving a geometrical picture within the Hilbert space of the system. We describe the relationship between the basis-independent and the basis-dependent approaches, and argue that many advantages exist for the former method. The formalism is illustrated with several numerical examples which show that the states can be characterized in a simple and effective manner.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗