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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Closed-form solutions for linear regulator design of mechanical systems including optimal weighting matrix selection

Vibration in modern structural and mechanical systems can be reduced in amplitude by increasing stiffness, redistributing stiffness and mass, and/or adding damping if design techniques are available to do so. Linear Quadratic Regulator (LQR) theory in modern multivariable control design, attacks the general dissipative elastic system design problem in a global formulation. The optimal design, however, allows electronic connections and phase relations which are not physically practical or possible in passive structural-mechanical devices. The restriction of LQR solutions (to the Algebraic Riccati Equation) to design spaces which can be implemented as passive structural members and/or dampers is addressed. A general closed-form solution to the optimal free-decay control problem is presented which is tailored for structural-mechanical system. The solution includes, as subsets, special cases such as the Rayleigh Dissipation Function and total energy. Weighting matrix selection is a constrained choice among several parameters to obtain desired physical relationships. The closed-form solution is also applicable to active control design for systems where perfect, collocated actuator-sensor pairs exist.

Hanks, Brantley R.↗
NASA NTRS1982

Determination of eigenvalues of dynamical systems by symbolic computation

A symbolic computation technique for determining the eigenvalues of dynamical systems is described wherein algebraic operations, symbolic differentiation, matrix formulation and inversion, etc., can be performed on a digital computer equipped with a formula-manipulation compiler. An example is included that demonstrates the facility with which the system dynamics matrix and the control distribution matrix from the state space formulation of the equations of motion can be processed to obtain eigenvalue loci as a function of a system parameter. The example chosen to demonstrate the technique is a fourth-order system representing the longitudinal response of a DC 8 aircraft to elevator inputs. This simplified system has two dominant modes, one of which is lightly damped and the other well damped. The loci may be used to determine the value of the controlling parameter that satisfied design requirements. The results were obtained using the MACSYMA symbolic manipulation system.

Howard, J. C.↗
NASA NTRS2015

An Efficient Solution Method for Multibody Systems with Loops Using Multiple Processors

This paper describes a multibody dynamics algorithm formulated for parallel implementation on multiprocessor computing platforms using the divide-and-conquer approach. The system of interest is a general topology of rigid and elastic articulated bodies with or without loops. The algorithm divides the multibody system into a number of smaller sets of bodies in chain or tree structures, called "branches" at convenient joints called "connection points", and uses an Order-N (O (N)) approach to formulate the dynamics of each branch in terms of the unknown spatial connection forces. The equations of motion for the branches, leaving the connection forces as unknowns, are implemented in separate processors in parallel for computational efficiency, and the equations for all the unknown connection forces are synthesized and solved in one or several processors. The performances of two implementations of this divide-and-conquer algorithm in multiple processors are compared with an existing method implemented on a single processor.

Multibody dynamics↗
DOE OSTI2015

Multiscale modeling and characterization for performance and safety of lithium-ion batteries

Lithium-ion batteries are highly complex electrochemical systems whose performance and safety are governed by coupled nonlinear electrochemical-electrical-thermal-mechanical processes over a range of spatiotemporal scales. Gaining an understanding of the role of these processes as well as development of predictive capabilities for design of better performing batteries requires synergy between theory, modeling, and simulation, and fundamental experimental work to support the models. This paper presents the overview of the work performed by the authors aligned with both experimental and computational efforts. In this paper, we describe a new, open source computational environment for battery simulations with an initial focus on lithium-ion systems but designed to support a variety of model types and formulations. This system has been used to create a three-dimensional cell and battery pack models that explicitly simulate all the battery components (current collectors, electrodes, and separator). The models are used to predict battery performance under normal operations and to study thermal and mechanical safety aspects under adverse conditions. This paper also provides an overview of the experimental techniques to obtain crucial validation data to benchmark the simulations at various scales for performance as well as abuse. We detail some initial validation using characterization experiments such as infrared and neutron imaging and micro-Raman mapping. In addition, we identify opportunities for future integration of theory, modeling, and experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1979

The computation of isentropic atmospheric trajectories using a 'discrete model' formulation

An explicit technique for calculating atmospheric trajectories is presented as an alternative method to the standard implicit scheme of Danielsen (1961). The technique uses the inviscid equations of motion and the discrete model formulation derived by Greenspan (1972, 1973) to compute trajectories on isentropic surfaces, assuming adiabatic flow. The discrete model formulation is designed specifically for a Lagrangian system and objectively accounts for the geostrophic departures, local psi-tendencies, and the subsequent accelerations along the entire length of the trajectory. Application of the discrete formulation to a diagnostic case study yielded favorable results.

Petersen, R. A.↗
NASA NTRSDate not supplied

Propulsion Airframe Aeroacoustics and Aircraft System Noise Flight Research Test: NASA Overview

Motivations and the formulation of the objectives are described for an ambitious flight research test conducted in collaboration between NASA and The Boeing Company. The Propulsion Airframe Aeroacoustics and Aircraft System Noise Flight Research Test was executed by the Boeing ecoDemonstrator Program with an Etihad Airways Boeing 787-10 aircraft. Five key technical approaches are described that were used to accomplish the more unconventional and challenging objectives of the research. In addition to the modern technology of the 787, these include evaluating multiple acoustic shielding and reflection effects using both the wing and fuselage in straight flight and banking flight paths, hardwalling the aft duct liner, and utilizing an extensive instrumentation array both on the ground and on the aircraft. The research level version of NASA’s aircraft system noise prediction capability including a new acoustic scattering method were used to confirm the technical approaches and provide guidance for the detailed design of instrumentation and flight test plans. Initial comparisons are shown of flight test data to predictions using the research level of NASA’s Aircraft Noise Prediction Program. Significant prediction challenges are revealed when compared to these high-quality data, while major progress is shown both in the measurement and the prediction of propulsion airframe aeroacoustic effects in flight. This comparison marks the first rigorous validation with modern flight test data and establishes a basis to further understand and quantify the state of NASA’s current capability and to develop improvements in the fidelity of aircraft system noise predictions. The context of multiple coordinated companion research papers is described as well as future plans for the continued analysis of this flight research data.

Propulsion Airframe Aeroacoustics↗
DOE OSTI2016

Tuning the driving force for exciton dissociation in single-walled carbon nanotube heterojunctions

Understanding the kinetics and energetics of interfacial electron transfer in molecular systems is crucial for the development of a broad array of technologies, including photovoltaics, solar fuel systems and energy storage. The Marcus formulation for electron transfer relates the thermodynamic driving force and reorganization energy for charge transfer between a given donor/acceptor pair to the kinetics and yield of electron transfer. Here we investigated the influence of the thermodynamic driving force for photoinduced electron transfer (PET) between single-walled carbon nanotubes (SWCNTs) and fullerene derivatives by employing time-resolved microwave conductivity as a sensitive probe of interfacial exciton dissociation. For the first time, we observed the Marcus inverted region (in which driving force exceeds reorganization energy) and quantified the reorganization energy for PET for a model SWCNT/acceptor system. The small reorganization energies (about 130 meV, most of which probably arises from the fullerene acceptors) are beneficial in minimizing energy loss in photoconversion schemes.

14 SOLAR ENERGY↗
DOE OSTI2013

Vapor condensation onto a non-volatile liquid drop

Molecular dynamics simulations of miscible and partially miscible binary Lennard–Jones mixtures are used to study the dynamics and thermodynamics of vapor condensation onto a non-volatile liquid drop in the canonical ensemble. When the system volume is large, the driving force for condensation is low and only a submonolayer of the solvent is adsorbed onto the liquid drop. A small degree of mixing of the solvent phase into the core of the particles occurs for the miscible system. At smaller volumes, complete film formation is observed and the dynamics of film growth are dominated by cluster-cluster coalescence. Mixing into the core of the droplet is also observed for partially miscible systems below an onset volume suggesting the presence of a solubility transition. We also develop a non-volatile liquid drop model, based on the capillarity approximations, that exhibits a solubility transition between small and large drops for partially miscible mixtures and has a hysteresis loop similar to the one observed in the deliquescence of small soluble salt particles. The properties of the model are compared to our simulation results and the model is used to study the formulation of classical nucleation theory for systems with low free energy barriers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

Parameter estimation in nonlinear distributed systems - Approximation theory and convergence results

An abstract approximation framework and convergence theory is described for Galerkin approximations applied to inverse problems involving nonlinear distributed parameter systems. Parameter estimation problems are considered and formulated as the minimization of a least-squares-like performance index over a compact admissible parameter set subject to state constraints given by an inhomogeneous nonlinear distributed system. The theory applies to systems whose dynamics can be described by either time-independent or nonstationary strongly maximal monotonic operators defined on a reflexive Banach space which is densely and continuously embedded in a Hilbert space. It is demonstrated that if readily verifiable conditions on the system's dependence on the unknown parameters are satisfied, and the usual Galerkin approximation assumption holds, then solutions to the approximating problems exist and approximate a solution to the original infinite-dimensional identification problem.

Banks, H. T.↗
DOE OSTI2022

Two-dimensional diamonds from sp 2 -to- sp 3 phase transitions

The ability to change material properties through phase engineering has long been sought, with the goal of ad hoc tunability of the physical and chemical properties of the transformed phases. The synthesis and study of graphene have made it possible to explore the mechanisms of 2D phase transformations, opening up paths towards the formation of 2D diamond and diamond thin films. This Review examines the state-of-the-art phase transformations in 2D graphitic systems and beyond. The theoretical models formulated to describe the sp2-to-sp3 transitions from graphene to 2D diamond and the experimental processes developed to induce the transition to 2D diamond are discussed, focusing on the transformations induced by chemical functionalization and pressure. The effects of different structural and environmental factors on the evolution of the phase transformations and on the properties of the transformed diamond phases are explored. Without comprehensively reviewing phase transitions in all 2D materials, we briefly mention hexagonal boron nitride, phosphorene, transition metal dichalcogenides and MXenes systems. Lastly, the Review delves into the technologies and applications of phase transformations in 2D materials and the opportunities for this field.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

An extension of Darcy’s law incorporating dynamic length scales

We propose a physics-based, macroscale formulation of multiphase porous-media flows that both honors the validity of Darcy’s law in steady or near-steady flows and accommodates the effects of heterogeneities and nonlinearities in unsteady flows. The new formulation recognizes that parameters characterizing the system operate at different length-scales. In particular, the use of Darcy’s law, predicated on the assumption of near-steady flows, requires dynamic length scales, owing to the possibility of rapid fluctuations in rock properties and fluid saturations attributable to heterogeneity and nonlinearity. We accommodate these dynamic length-scales through dynamic spatial averaging. The length and position of the averaging window are characterized by the length of the mixing zone and the direction of propagation of information in the transport process. Here, we validate the proposed formulation by comparing highly accurate,two dimensional numerical solutions against core-scale displacement experiments. The proposed paradigm is consistent with the classical multiphase Darcy formulation, in the sense that the latter also represents the results of an averaging approach.

42 ENGINEERING↗
NASA NTRS1980

MCC level C formulation requirements. Entry guidance and entry autopilot, optional TAEM targeting

The level C software formulations requirements for the entry guidance and the simplified autopilot to be used by the Mission Control Center (MCC) entry processor are presented. The modifications required to functionally simulate optional Terminal Area Energy Management (TAEM) targeting capability (OTT) are incorporated. Implementation of this logic in the MCC must be coordinated with flight software OTT implementation and MCC TAEM guidance OTT. The entry guidance logic is based on the orbiter avionics entry guidance software. Descriptions of the entry guidance detailed formulation requirements, the detailed autopilot formulation requirements, and the targeting routine are given. Also included are a definition of coordinate systems, a list of parameter definitions for the software formulations, and a set of formulation flow charts.

Harpold, J. C.↗
NASA NTRS1988

Characterization of aluminum/RP-1 gel propellant properties

Research efforts are being conducted by the NASA Lewis Research Center to formulate and characterize the properties of Al/RP-1 and RP-1 gelled propellants for rocket propulsion systems. Twenty four different compositions of gelled fuels were formulated with 5 and 16 micron, atomized aluminum powder in RP-1. The total solids concentration in the propellant varied from 5 to 60 wt percent. Tests were conducted to evaluate the stability and rheological characteristics of the fuels. Physical separation of the solids occurred in fuels with less than 50 wt percent solids concentration. The rheological characteristics of the Al/RP-1 fuels varied with solids concentration. Both thixotropic and rheopectic gel behavior were observed. The unmetallized RP-1 gels, which were formulated by a different technique than the Al/RP-1 gels, were highly viscoelastic. A history of research efforts which were conducted to formulate and characterize the properties of metallized propellants for various applications is also given.

Rapp, Douglas C.↗
NASA NTRS2006

NASA Langley Research Center Systems Analysis & Concepts Directorate Participation in the Exploration Systems Architecture Study

The NASA Langley Research Center (LaRC) Systems Analysis & Concepts Directorate (SACD) began studying human exploration missions beyond low Earth orbit (LEO) in the year 1999. This included participation in NASA s Decadal Planning Team (DPT), the NASA Exploration Team (NExT), Space Architect studies and Revolutionary Aerospace Systems Concepts (RASC) architecture studies that were used in formulating the new Vision for Space Exploration. In May of 2005, NASA initiated the Exploration Systems Architecture Study (ESAS). The primary outputs of the ESAS activity were concepts and functional requirements for the Crewed Exploration Vehicle (CEV), its supporting launch vehicle infrastructure and identification of supporting technology requirements and investments. An exploration systems analysis capability has evolved to support these functions in the past and continues to evolve to support anticipated future needs. SACD had significant roles in supporting the ESAS study team. SACD personnel performed the liaison function between the ESAS team and the Shuttle/Station Configuration Options Team (S/SCOT), an agency-wide team charged with using the Space Shuttle to complete the International Space Station (ISS) by the end of Fiscal Year (FY) 2010. The most significant of the identified issues involved the ability of the Space Shuttle system to achieve the desired number of flights in the proposed time frame. SACD with support from the Kennedy Space Center performed analysis showing that, without significant investments in improving the shuttle processing flow, that there was almost no possibility of completing the 28-flight sequence by the end of 2010. SACD performed numerous Lunar Surface Access Module (LSAM) trades to define top level element requirements and establish architecture propellant needs. Configuration trades were conducted to determine the impact of varying degrees of segmentation of the living capabilities of the combined descent stage, ascent stage, and other elements. The technology assessment process was developed and implemented by SACD as the ESAS architecture was refined. SACD implemented a rigorous and objective process which included (a) establishing architectural functional needs, (b) collection, synthesis and mapping of technology data, and (c) performing an objective decision analysis resulting in technology development investment recommendations. The investment recommendation provided budget, schedule, and center/program allocations to develop required technologies for the exploration architecture, as well as the identification of other investment opportunities to maximize performance and flexibility while minimizing cost and risk. A summary of the trades performed and methods utilized by SACD for the Exploration Systems Mission Directorate (ESAS) activity is presented along with how SACD is currently supporting the implementation of the Vision for Space Exploration.

Keyes, Jennifer↗
DOE OSTI2022

Model Predictive Control for a Grid-interactive Efficient Thermal Storage-integrated Heat Pump System

Building heating and cooling systems can be used to overcome the mismatch between the intermittent supply of renewable power and the fluctuating demand for electricity. A novel underground thermal energy storage integrated with a dual-source heat pump has been proposed to mitigate the mismatch while meeting the thermal demand of buildings efficiently. Conventional thermostat control with heuristic rules cannot provide intelligent decisions to maximize the thermal efficiency and flexibility of the proposed system. Advanced control strategies like model predictive control (MPC) have provided a new paradigm for grid-interactive efficient building operation with the advancement of computation and sensing. This study developed an MPC for the proposed system to provide grid service for Demand Side Management and minimize the operating cost of building owners. A control-oriented dynamic model of the proposed system has been developed. Given an objective function and proper constraints, an optimization problem is formulated to determine the optimal control strategy of the system. Dynamic Programming is adopted to solve the optimization problem. A rule-based control (RBC) is also developed to achieve similar goals. Short-term simulations are conducted to compare the system performance resulting from the two controls. The simulation results indicate that the MPC performs more intelligently than the RBC in charging thermal energy storage and selecting heat pump sources by taking advantage of the predicted cooling demands of the building and the performance of the integrated system. As a result, the MPC could save energy and reduce operating costs compared with the RBC. A case study shows that, for a 3-day operation, the MPC saves 36.9% energy and reduces 38.5% operating cost compared with the RBC.

Shi, Liang↗
NASA NTRS1974

A direct approach to the design of linear multivariable systems

Design of multivariable systems is considered and design procedures are formulated in the light of the most recent work on model matching. The word model matching is used exclusively to mean matching the input-output behavior of two systems. The term is used in the frequency domain to indicate the comparison of two transfer matrices containing transfer functions as elements. Design methods where non-interaction is not used as a criteria were studied. Two design methods are considered. The first method of design is based solely upon the specification of generalized error coefficients for each individual transfer function of the overall system transfer matrix. The second design method is called the pole fixing method because all the system poles are fixed at preassigned positions. The zeros of terms either above or below the diagonal are partially fixed via steady state error coefficients. The advantages and disadvantages of each method are discussed and an example is worked to demonstrate their uses. The special cases of triangular decoupling and minimum constraints are discussed.

Agrawal, B. L.↗