Search NASASearch

SEARCH · Search NASA

Results for “formulation generation”

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 109 records · Page 6

Nonlinear aerodynamic modeling of flap oscillations in transonic flow - A numerical validation

The regime of validity of a nonlinear aerodynamic force and moment formulation, based on concepts from nonlinear functional analysis and applicable to a transonic airfoil with a deflecting flap, is investigated. A time-dependent finite difference technique is used to evaluate the aerodynamic data of the formulation in terms of specified, characteristic motions. Flap-motion histories are generated from the flap inertial equations of motion, with aerodynamic reactions specified by the moment formulation. The motion histories depicting the cases of decaying and growing flap oscillations are compared with histories generated through simultaneous, coupled solution of the fluid-dynamic equations and flap inertial equations of motion. The range of applicability of the formulation is discussed.

Chyu, W. J.

Retrofitting the AutoBayes Program Synthesis System with Concrete Syntax

AutoBayes is a fully automatic, schema-based program synthesis system for statistical data analysis applications. Its core component is a schema library. i.e., a collection of generic code templates with associated applicability constraints which are instantiated in a problem-specific way during synthesis. Currently, AutoBayes is implemented in Prolog; the schemas thus use abstract syntax (i.e., Prolog terms) to formulate the templates. However, the conceptual distance between this abstract representation and the concrete syntax of the generated programs makes the schemas hard to create and maintain. In this paper we describe how AutoBayes is retrofitted with concrete syntax. We show how it is integrated into Prolog and describe how the seamless interaction of concrete syntax fragments with AutoBayes's remaining legacy meta-programming kernel based on abstract syntax is achieved. We apply the approach to gradually mitigate individual schemas without forcing a disruptive migration of the entire system to a different First experiences show that a smooth migration can be achieved. Moreover, it can result in a considerable reduction of the code size and improved readability of the code. In particular, abstracting out fresh-variable generation and second-order term construction allows the formulation of larger continuous fragments.

Fischer, Bernd

Thermal Dynamics and Lithium Plating Detection in High‐Power Li‐Ion Batteries for eVTOL Applications

The rapid electrification of aerial transportation is driving the need for high‐performance Li‐ion batteries that can operate reliably under stringent thermal and safety constraints. The unique mission profile of electric Vertical Take‐off and Landing (eVTOL) aircraft necessitates a focused investigation into the thermal behavior and safety characteristics of these batteries. Here, in this study, operando isothermal microcalorimetry is employed to examine the thermal evolution of Li‐ion batteries under cycling conditions representative of eVTOL operations. These findings reveal that high‐power discharge events—such as those during take‐off and landing—shift the thermal response toward exothermic behavior, in contrast to the typically endothermic response expected under near‐equilibrium cycling conditions. Additionally, the results suggest that advanced electrolyte formulations may help suppress excess heat generation, thereby improving battery safety. Notably, the calorimetric results exhibit a distinct thermal signature associated with lithium plating, offering a potential diagnostic for detecting Li plating during eVTOL operation. Overall, this study demonstrates the utility of isothermal microcalorimetry as a valuable tool for assessing thermal risks in Li‐ion batteries for eVTOL applications, and highlights the importance of targeted design strategies to mitigate safety hazards during high‐power demand scenarios.

Li plating

Determining the profitability of energy storage over its life cycle using levelized cost of storage

Levelized cost of storage (LCOS) can be a simple, intuitive, and useful metric for determining whether a new energy storage plant would be profitable over its life cycle and to compare the cost of different energy storage technologies. However, researchers and industry decision makers still use conflicting definitions of LCOS. For example, some include charging cost, while others only include round trip efficiency (RTE) losses. Additionally, inputs to the existing formulations are not specific enough to generate repeatable results across studies, which reduces trust in the metric. To push for standardization in economic assessment of batteries and other energy storage devices, the authors review existing definitions of LCOS and identify the desired characteristics for a standard. They then propose a new definition and demonstrate that it fits these characteristics very well relative to other prominent options. Unit analysis is applied to this proposed definition to provide a deeper understanding of the equations and to demonstrate its effectiveness. Finally, the sensitivity of LCOS to different input parameters is investigated to help users understand how to compare analyses from literature to their own. The authors also provide a spreadsheet and a Python script to streamline adoption of the proposed definition.

25 ENERGY STORAGE

An efficient surrogate model of secondary electron formation and evolution

This work extends the adjoint-deep learning framework for runaway electron (RE) evolution, developed by McDevitt et al. [Phys. Plasmas 32, 042503 (2025)], to account for large-angle collisions. By incorporating large-angle collisions, the framework allows the avalanche of REs to be captured, an essential component of RE dynamics. This extension is accomplished by using a Rosenbluth–Putvinski approximation to estimate the distribution of secondary electrons generated by large-angle collisions. By evolving both the primary and multiple generations of secondary electrons, the present formulation can capture both the detailed temporal evolution of a RE population beginning from an arbitrary initial momentum space distribution, along with providing approximations to the saturated growth and decay rates of the RE population. Predictions of the adjoint-deep learning framework are verified against a traditional RE solver, with good agreement present across a broad range of parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Dynamic Modeling, Trajectory Optimization, and Linear Control of Cable-Driven Parallel Robots for Automated Panelized Building Retrofits

The construction industry faces a growing need for automation to reduce costs, improve accuracy and productivity, and address labor shortages. One area that stands to benefit significantly from automation is panelized prefabricated building envelope retrofits, which can improve a building’s energy efficiency in heating and cooling interior spaces. In this paper, we propose using cable-driven parallel robots (CDPRs), which can effectively lift and handle large objects, to install these panels. However, implementing CDPRs presents significant challenges because of their nonlinear dynamics, complex trajectory planning, and precise control requirements. To tackle these challenges, this work focuses on a new application of established control and trajectory optimization theories in a CDPR simulation of a building envelope retrofit under real-world conditions. We first model the dynamics of CDPRs, highlighting the critical role of damping in system behavior. Building on this dynamic model, we formulate a trajectory optimization problem to generate feasible and efficient motion plans for the robot under operational and environmental constraints. Given the high precision required in the construction industry, accurately tracking the optimized trajectory is essential. However, challenges such as partial observability and external vibrations complicate this task. To address these issues, a Linear Quadratic Gaussian control framework is applied, enabling the robot to track the optimized trajectories with precision. Simulation results show that the proposed controller enables precise end effector positioning with errors under 4 mm, even in the presence of external wind disturbances. Through comprehensive simulations, our approach allows for an in-depth exploration of the system’s nonlinear dynamics, trajectory optimization, and control strategies under controlled yet highly realistic conditions. The results demonstrate the feasibility of CDPRs for automating panel installation and provide insights into their practical deployment.

CDPR

Warm fog dispersal

The charged particle generator was further tested after some design modification. The generator performance was measured with additional instrumentation and found to confirm previous measurements. Plans for a field testing were than developed. The overall status of the program and the field test plans were presented to a group of atmospheric scientists and electrostatic experts at the NASA/MSFC sponsored USRA Workshop on Electrostatic Fog Dispersal at NCAR, Boulder, Colorado discussed in previous sections. The recommendations from this workshop are being evaluated as to whether NASA should proceed with the field test or whether further theoretical research on the phenomenon of electrostatic fog dispersal and additional development of the charged particle generator should be carried out. Information obtained from the USRA Workshop clearly identified three physical mechanisms that could possibly influence the fog dispersal process, which heretofore have not been considered, and which may provide additional insight to the direction of further fog dispersal work. These mechanisms are: the effect of corona discharge on the electric field strength at the surface, the influx of fog into the cleared volume by turbulent diffusion, and the increase in supersaturation as liquid water is removed, activating haze particles, and thus generating more fog. Plans are being formulated to investigate these mechanisms.

Frost, W.

The application of vortex theory to the optimum swept propeller

It is shown that an optimum propeller generating a swept wake must satisfy the Betz condition, at least to a first order. A numerical solution for swept propellers generating a rigid helicoidal wake is formulated and some results are presented. These results indicate that sweep has a significant effect on Goldstein's kappa factor, particularly at high advance ratios typical of those at which advanced turboprops operate.

Mccormick, B. W.

Rate dependent constitutive models for fiber reinforced polymer composites

A literature survey was conducted to assess the state-of-the-art in rate dependent constitutive models for continuous fiber reinforced polymer matrix composite (PMC) materials. Several recent models which include formulations for describing plasticity, viscoelasticity, viscoplasticity, and rate-dependent phenomenon such as creep and stress relaxation are outlined and compared. When appropriate, these comparisons include brief descriptions of the mathematical formulations, the test procedures required for generating material constants, and details of available data comparing test results to analytical predictions.

Gates, Thomas S.

Integrated analysis and design for controlled nonlinear multibody systems

A hybrid approach for integrated analysis and design of nonlinear multibody systems is proposed. Based on the approach, a general-purpose design tool is developed for nonlinear dynamic systems subjected to nonlinear design constraints. For analysis purposes, second-order nonlinear equations of motion of the system are automatically generated using general-purpose multibody formulations. Once the equations are solved, they are written in first-order form to take advantage of the first-order formulations of design sensitivity analysis for dynamic systems. A nonlinear programming technique is then used to optimize the nonlinear systems, such that a nonlinear cost function is minimized and performance constraints are satisfied. The approach proved to be very general and useful for design and analysis of controlled multibody systems. The approach is applied to the design of passive dynamic controllers for slewing control of multibody systems.

Wu, Shih-Chin

TE/TM Simulations of Interferometric Measurements

Interferometric synthetic aperture radar (IFSAR) measurements at X-, C-, L-, and P-band are used to derive ground topography at meter level resolution. Interpretation of the derived topography requires attention due to the complex interaction of the radar signal with ground cover. The presence of penetrable surfaces such as vegetation, and tree canopies poses a challenge since the depth of penetration depends on a number of parameters such as the operating radar frequency, polarization, incident angle, as well as terrain structure. The dependence of the reconstructed topography on polarization may lead to the characterization of the ground cover. Simulation of interferometric measurements is useful for interpretation of the derived topography (B. Houshmand, Proceedings of URSI, 314, 1997). In this talk , time domain simulations for interferometric measurement for TE- and TM- polarization are presented. Time domain simulation includes the effects of the surface material property as well geometry comparable the radar signal wavelength (B. Houshmand, Proceedings of the URSI, 25, 1998). The IFSAR simulation is carried out in two steps. First, the forward scattering data is generated based on full wave analysis. Next, the electromagnetic information is inverted to generate surface topography. This inversion is based on the well known IFSAR processing technique which is composed of signal compression, and formation of an interferogram. The full wave forward scattering data is generated by the scattered-field formulation of the FDTD algorithm. The simulation is carried out by exciting the computational domain by a radar signal. The scattered field is then computed and translated to the receiving interferometric antennas using the time-domain Huygen's principle. The inversion process starts by compressing the time-domain data. The range compressed data from both receivers are then coregistered to form an interferogram. The resulting interferogram is then related to the ground topography using the radar imaging geometry. In this talk, the simulation results are compared with the C-band TM IFSAR derived topography, and the TE/TM SAR images at L-Band.

Houshmand, Bijan

RainCube and its legacy for the next generation of spaceborne cloud and precipitation radars

The multidimensional challenge of observing cloud and precipitation processes from space has motivated and focused a number of technological developments in the last decade. They span from novel compact instrument architectures to modular millimeter wave phased array solutions. In this paper we summarize the main capabilities offered by three instruments from this new generation, and discuss how they affect formulation of future observing systems.

Rodriguez-Monje, Raquel

Unsteady heat transfer in turbine blade ducts: Focus on combustor sources

Thermal waves generated by either turbine rotor blades cutting through nonuniform combustor temperature fields or unsteady burning could lead to thermal fatigue cracking in the blades. To determine the magnitude of the thermal oscillation in blades with complex shapes and material compositions, a finite element Galerkin formulation has been developed to study combustor generated thermal wave propagation in a model two-dimensional duct with a uniform plug flow profile. The reflection and transmission of the thermal waves at the entrance and exit boundaries are determined by coupling the finite element solutions at the entrance and exit to the eigenfunctions of an infinitely long adiabatic duct. Example solutions are presented. In general, thermal wave propagation from an air passage into a metallic blade wall is small and not a problem. However, if a thermal barrier coating is applied to a metallic surface under conditions of a high heat transfer, a good impedance match is obtained and a significant portion of the thermal wave can pass into the blade material.

Baumeister, Kenneth J.

Unsteady heat transfer in turbine blade ducts - Focus on combustor sources

Thermal waves generated by either turbine rotor blades cutting through nonuniform combustor temperature fields or unsteady burning could lead to thermal fatigue cracking in the blades. To determine the magnitude of the thermal oscillation in blades with complex shapes and material compositions, a finite element Galerkin formulation has been developed to study combustor generated thermal wave propagation in a model two-dimensional duct with a uniform plug flow profile. The reflection and transmission of the thermal waves at the entrance and exit boundaries are determined by coupling the finite element solutions at the entrance and exit to the eigenfunctions of an infinitely long adiabatic duct. Example solutions are presented. In general, thermal wave propagation from an air passage into a metallic blade wall is small and not a problem. However, if a thermal barrier coating is applied to a metallic surface under conditions of a high heat transfer, a good impedance match is obtained and a significant portion of the thermal wave can pass into the blade material.

Baumeister, K. J.

The fragmentation of 510 MeV/nucleon iron-56 in polyethylene. II. Comparisons between data and a model

The results of a Monte Carlo model for calculating fragment fluences and LET spectra are compared to data taken with 600 MeV/nucleon iron ions incident on an accelerator beamline configured for irradiation of biological samples, with no target and with 2, 5 and 8 cm of polyethylene. The model uses a multi-generation nuclear fragmentation code, coupled with a formulation of ionization energy loss based on the Bethe-Bloch equation. In the region where the data are reliable and the experimental acceptance is well understood, many of the features of the experimental spectra are well replicated by the model. To obtain good agreement with the experimental data, the model must allow for at least two generations of fragment production in the target.

NASA Discipline Number 99-99

Review of Propulsion Technologies for N+3 Subsonic Vehicle Concepts

NASA has set aggressive fuel burn, noise, and emission reduction goals for a new generation (N+3) of aircraft targeting concepts that could be viable in the 2035 timeframe. Several N+3 concepts have been formulated, where the term "N+3" indicate aircraft three generations later than current state-of-the-art aircraft, "N". Dramatic improvements need to be made in the airframe, propulsion systems, mission design, and the air transportation system in order to meet these N+3 goals. The propulsion system is a key element to achieving these goals due to its major role with reducing emissions, fuel burn, and noise. This report provides an in-depth description and assessment of propulsion systems and technologies considered in the N+3 subsonic vehicle concepts. Recommendations for technologies that merit further research and development are presented based upon their impact on the N+3 goals and likelihood of being operational by 2035.

Ashcraft, Scott W.

A General Procedure to Formulate 3D Elements for Finite Element Applications

This paper presents a general procedure to formulate and implement 3D elements of arbitrary order in meshes with multiple element types. This procedure includes obtaining shape functions and integration quadrature and establishing an approach for checking the generated element’s compatibility with adjacent elements’ surfaces. This procedure was implemented in Matlab, using its symbolic and graphics toolbox, and complied as a GUI interface named ShapeGen3D to provide finite element users with a tool to tailor elements according to their analysis needs. ShapeGen3D also outputs files with the element formulation needed to enable users to implement the generated elements in other programming languages or through user elements in commercial finite element software. Currently, finite element (FE) users are limited to employing element formulation available in the literature, commercial software, or existing element libraries. Thus, the developed procedure implemented in ShapeGen3D offers FEM users the possibility to employ elements beyond those readily available. The procedure was tested by generating the formulation for a brick element, a brick transition element, and higher-order hexahedron and tetrahedron elements that can be used in a spectral finite element analysis. The formulation obtained for the 20-node element was in perfect agreement with the formulation available in the literature. In addition, the results showed that the interpolation condition was met for all the generated elements, which provides confidence in the implementation of the process. Researchers and educators can use this procedure to efficiently develop and illustrate three-dimensional elements.

higher-order elements

A Methodology to Evaluate the Grid Reliability Impact of Oscillations Induced by Large Loads

The rapid growth of hyperscale AI data centers is bringing renewed attention to the reliability risk that sustained forced oscillations pose to bulk power systems, with cyclic computational workloads emerging as a new forcing source. Unlike the broadband, stochastic disturbances from traditional industrial loads such as arc furnaces, AI training and inference facilities can inject large active power swings concentrated at specific frequencies over extended durations - characteristics that existing grid planning practices do not account for. While the North American Electric Reliability Corporation (NERC) has recognized this gap and called for system-level studies of large load interconnections, no standardized methodology exists to screen, simulate, and quantify these risks at the planning stage. This report presents the Risk Assessment Tool for Large Load-induced Events (RATLLE), a Python-based, publicly available script suite developed at the Pacific Northwest National Laboratory to evaluate bulk power system reliability risks from data center-induced oscillations. RATLLE implements a three-module workflow: a screening module that identifies vulnerable interconnection locations and excitable system modes; a simulation module that models cyclic data center load behavior using a commercial positive sequence simulation platform; and an analysis module that computes risk metrics and generates interactive visualization dashboards. The risk metrics, formulated around simulation observables, map oscillation impacts to a three-stage severity scale spanning latent equipment fatigue through imminent cascading failure. The methodology is demonstrated on two Western Electricity Coordinating Council (WECC) system models: a publicly available 240-bus reduced representation and a detailed 2031 Heavy Winter planning case. Case studies illustrate that even modest 50 MW forced oscillations at resonant frequencies can produce wide-area power swings, N-1 security constraint violations, and cascading generator trips through protection actions - outcomes that would not occur under normal operating conditions without oscillations present. The results underscore the need for standardized oscillation impact assessment in large load interconnection studies and provide a reproducible, extensible framework for utilities to adopt or customize within their existing planning workflows.

Biswas, Shuchismita