Search NASA⌕ Search

SEARCH · Search NASA

Results for “Formulation”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

A New Formulation of Equivalent Effective Stratospheric Chlorine (EESC)

Equivalent effective stratospheric chlorine (EESC) is a convenient parameter to quantify the effects of halogens (chlorine and bromine) on ozone depletion in the stratosphere. We show and discuss a new formulation of EESC that now includes the effects of age-of-air dependent fractional release values and an age-of-air spectrum. This new formulation provides quantitative estimates of EESC that can be directly related to inorganic chlorine and bromine throughout the stratosphere. Using this EESC formulation, we estimate that human-produced ozone depleting substances will recover to 1980 levels in 2041 in the midlatitudes, and 2067 over Antarctica. These recovery dates are based upon the assumption that the international agreements for regulating ozone-depleting substances are adhered to. In addition to recovery dates, we also estimate the uncertainties in the estimated time of recovery. The midlatitude recovery of 2041 has a 95% confidence uncertainty from 2028 to 2049, while the 2067 Antarctic recovery has a 95% confidence uncertainty from 2056 to 2078. The principal uncertainties are from the estimated mean age-of-air, and the assumption that the mean age-of-air and fractional release values are time independent. Using other model estimates of age decrease due to climate change, we estimate that midlatitude recovery may be accelerated from 2041 to 2031.

Newman, P. A.↗
NASA NTRS2006

Application of p-Multigrid to Discontinuous Galerkin Formulations of the Poisson Equation

We investigate p-multigrid as a solution method for several different discontinuous Galerkin (DG) formulations of the Poisson equation. Different combinations of relaxation schemes and basis sets have been combined with the DG formulations to find the best performing combination. The damping factors of the schemes have been determined using Fourier analysis for both one and two-dimensional problems. One important finding is that when using DG formulations, the standard approach of forming the coarse p matrices separately for each level of multigrid is often unstable. To ensure stability the coarse p matrices must be constructed from the fine grid matrices using algebraic multigrid techniques. Of the relaxation schemes, we find that the combination of Jacobi relaxation with the spectral element basis is fairly effective. The results using this combination are p sensitive in both one and two dimensions, but reasonable convergence rates can still be achieved for moderate values of p and isotropic meshes. A competitive alternative is a block Gauss-Seidel relaxation. This actually out performs a more expensive line relaxation when the mesh is isotropic. When the mesh becomes highly anisotropic, the implicit line method and the Gauss-Seidel implicit line method are the only effective schemes. Adding the Gauss-Seidel terms to the implicit line method gives a significant improvement over the line relaxation method.

Helenbrook, B. T.↗
NASA NTRS2008

An Improved Formulation for Calorimetric Emittance Testing of Spacecraft Thermal Control Coatings

Spacecraft often really heavily on passive thermal control to maintain operating temperature. An important parameter in the spacecraft heat balance equation is the emittance of thermal control coatings as a function of coating temperature. One method for determining the emittance of spacecraft thermal control from elevated temperature to cryogenic temperatures relies on a calorimetric technique. The fundamental equation governing this test method can be found in numerous places in the literature and although it generally provides reasonable results, its formulation is based on a conceptual flaw that only becomes apparent when the sample temperature approaches the wall temperature during testing. This paper investigates the cause for this error and develops the correct formulation for calorimetric emittance testing. Experimental data will also be presented that illustrates the difference between the two formulations and the resulting difference in the calculated emittance.

Kauder, Lonny R.↗
NASA NTRS2012

Broadband Noise Prediction When Turbulence Simulation Is Available - Derivation of Formulation 2B and Its Statistical Analysis

We show that a simple modification of Formulation 1 of Farassat results in a new analytic expression that is highly suitable for broadband noise prediction when extensive turbulence simulation is available. This result satisfies all the stringent requirements, such as permitting the use of the exact geometry and kinematics of the moving body, that we have set as our goal in the derivation of useful acoustic formulas for the prediction of rotating blade and airframe noise. We also derive a simple analytic expression for the autocorrelation of the acoustic pressure that is valid in the near and far fields. Our analysis is based on the time integral of the acoustic pressure that can easily be obtained at any resolution for any observer time interval and digitally analyzed for broadband noise prediction. We have named this result as Formulation 2B of Farassat. One significant consequence of Formulation 2B is the derivation of the acoustic velocity potential for the thickness and loading terms of the Ffowcs Williams-Hawkings (FW-H) equation. This will greatly enhance the usefulness of the Fast Scattering Code (FSC) by providing a high fidelity boundary condition input for scattering predictions.

Farassat, Fereidoun↗
NASA NTRS2016

JPL's Foundry Furnace: Web-Based Concurrent Engineering for Formulation

The Jet Propulsion Laboratory’s Innovation Foundry is an enterprise tasked with shepherding space mission concepts through the formulation lifecycle. It oversees a number of “virtual teams” for the various stages of formulation. Among these is Team X, which has had considerable success over its more than 20-year history. In a Team X study, domain experts (including engineers devoted to the various spacecraft subsystems) work concurrently and collaboratively over several days to arrive at a feasible point design with a reasonable cost estimate. They use a set of linked Excel workbooks, each developed and approved by a responsible “line organization” within JPL. This toolset has served Team X well over the years, and has evolved since its inception. At the same time, the Innovation Foundry’s portfolio of formulation teams has expanded, and so has the scope of the design challenges they face. The A-Team runs workshop-like architecture studies to focus science investigations, generate mission concepts, assess feasibility, and explore trade spaces. Team Xc performs rapid point design in the style of Team X, but for CubeSats and small spacecraft, using a different toolset. Proposal teams further mature concepts to the point where they can be proposed. Recognizing the importance of trade space modeling combined with new IT services for providing and integrating data, JPL is developing the Foundry Furnace web-based software infrastructure. It will support A-Team, Team Xc, and Team X, providing study management, a catalog of hardware components, a library of re-usable analyses, and a design environment. It is a modernization of JPL’s concurrent engineering infrastructure, embracing the core concepts of Model-Based Systems Engineering, and built with modern software design philosophies.

Murphy, Jonathan↗
NASA NTRS2022

A Cohesive Element Formulation for Simulation of General Mode I, II and III Delamination

Cohesive element formulations proposed thus far have assumed that Mode II and Mode III fracture processes, and hence traction-separation laws, are identical. However, experimental observations suggest that such an assumption may not hold in general. The approach described herein consists of a cohesive element formulation in which the assumption of identical Mode II and Mode III fracture processes is not required. This assumption is alleviated by estimating the normal to the delamination front based on the gradient of the displacement jumps determined within each element. The estimated normal is used to rotate the cohesive element coordinate system. Determining the displacement jumps in the rotated coordinate system enables the distinction between Mode II sliding, and Mode III tearing. This information is subsequently used to extend a recently proposed cohesive formulation, resulting in an approach that can simulate mixed-mode I/II/III fracture via prescribing piecewise-linear traction-separation laws combined with a 3D mixed-mode fracture criteria.

Delamination↗
NASA NTRSDate not supplied

A Physics-Based Work-Energy Formulation for Real-Time Trajectory Guidance of A Lunar Lander

Throughout the years, many researchers have calculated and optimized trajectory solutions for lunar landing systems by employing sophisticated mathematical methods, that include: Hamilton’s Principle of Variation, Pontryagin’s maximum principle, and well known convex-optimization techniques among others. Many of these approaches typically require expensive computational resources to achieve convergence in the solution. In an effort to reduce complexity and the computational load required to obtain real-time guidance commands, a simple physics-based work-energy approach has been formulated. This approach is based on the dissipation of the mechanical energy of the vehicle to its final desired energy state required to achieve a safe landing. The rocket engine(s) employed during landing (among other maneuvers) dissipates mechanical energy by both doing work against the velocity vector of the vehicle (thus defining the trajectory path), and by jettisoning mass. Therefore, by solving the energy dissipation problem at every step of the maneuver, a much simpler formulation that naturally and quickly attains convergence is obtained. This formulation is not limited to approach, landing, and divert maneuvers, but in principle it can be employed during de-orbiting, braking burn, ascent, as well as orbit insertion.

Guidance↗
NASA NTRSDate not supplied

A Physics-Based Work-Energy Formulation for Real-Time Trajectory Guidance of a Lunar Lander

Throughout the years, many researchers have calculated and optimized trajectory solutions for lunar landing systems by employing sophisticated mathematical methods, that include: Hamilton’s Principle of Variation, Pontryagin’s maximum principle, and well known convex-optimization techniques among others. Many of these approaches typically require expensive computational resources to achieve convergence in the solution. In an effort to reduce complexity and the computational load required to obtain real-time guidance commands, a simple physics-based work-energy approach has been formulated. This approach is based on the dissipation of the mechanical energy of the vehicle to its final desired energy state required to achieve a safe landing. The rocket engine(s) employed during landing (among other maneuvers) dissipates mechanical energy by both doing work against the velocity vector of the vehicle (thus defining the trajectory path), and by jettisoning mass. Therefore, by solving the energy dissipation problem at every step of the maneuver, a much simpler formulation that naturally and quickly attains convergence is obtained. This formulation is not limited to approach, landing, and divert maneuvers, but in principle it can be employed during de-orbiting, braking burn, ascent, as well as orbit insertion.

Guidance↗
DOE OSTI2025

A Domain-Decomposed A-ϕ Formulation Based on Lagrange Multipliers for Low-Frequency Problems

A domain-decomposed A-ϕ formulation based on Lagrange multipliers is proposed to simulate low-frequency elec- tromagnetic problems. This method partitions the computational domain into smaller subdomains, allowing each subdomain to be independently formulated using Lagrange multipliers as Dirichlet boundary conditions, while ensuring continuity of the fields across the interfaces. A mixed finite element method, utilizing both vector and scalar basis functions, is employed to discretize the formulation, resulting in a global system to be solved. The proposed method is validated using TEAM Problem 7 at 50 Hz, demonstrating its effectiveness in handling complex geometries and addressing the low-frequency breakdown issues commonly encountered in traditional finite element methods.

Hossain, Amzad↗
NASA NTRS2012

Early Formulation Model-centric Engineering on NASA's Europa Mission Concept Study

The proposed Jupiter Europa Orbiter and Jupiter Ganymede Orbiter missions were formulated using current state-of-the-art MBSE facilities: - JPL's TeamX, Rapid Mission Architecting - ESA's Concurrent Design Facility - APL's ACE Concurrent Engineering Facility. When JEO became an official "pre-project" in Sep 2010, we had already developed a strong partnership with JPL's Integrated Model Centric Engineering (IMCE) initiative; decided to apply Architecting and SysML-based MBSE from the beginning, begun laying these foundations to support work in Phase A. Release of Planetary Science Decadal Survey and FY12 President's Budget in March 2011 changed the landscape. JEO reverted to being a pre-phase A study. A conscious choice was made to continue application of MBSE on the Europa Study, refocused for early formulation. This presentation describes the approach, results, and lessons.

early formulation↗
DOE OSTI2024

Glass Design Using Machine Learning Property Models with Prediction Uncertainties: Nuclear Waste Glass Formulation

The United States Department of Energy is responsible for managing the legacy nuclear waste stored in underground tanks at the Hanford Site. The waste will be separately vitrified as low-activity waste and high-level waste fractions. Waste glass formulation algorithms have been traditionally developed using partial quadratic mixture property-composition models. Recently, machine learning (ML) techniques have been used to predict glass properties and discover new glass materials for nuclear waste vitrification, and these advancements can be utilized to improve waste glass composition design. In this proof-of-principle study, ML algorithms such as Gaussian process regression (GPR) were used to interpolate glass properties (e.g., viscosity, electrical conductivity, chemical durability). After selecting appropriate sets of GPR hyper-parameters for each property, an optimization program was developed to formulate glass compositions to maximize waste loading while simultaneously satisfying property within constraints. The results of the ML-based waste loadings and glass compositions were compared to those obtained using the traditional methods. Comparing to the previous glass design framework, the ML-based optimization methods offer improved glass designs and a streamlined approach to generation of optimally designed data and near real-time updates.

glass formulation, machine learning, constraints, ↗
DOE OSTI2025

Assessing the coating implications of slurry formulations in NCMA and LFMP blend cathodes for lithium-ion batteries

This paper investigates the implications of slurry formulations and electrode processing on the performance of lithium-ion battery (LIB) cathodes, focusing on Nickel Cobalt Manganese Aluminum (NCMA) and Lithium Ferro Manganese Phosphate (LFMP) blends. Through a comprehensive examination of electrode processing steps, from material selection to coating application, we elucidate the critical role of processing parameters in shaping electrode morphology and electrochemical behavior. Rheological studies reveal the influence of slurry composition on viscosity and flow behavior, highlighting the importance of achieving optimal rheological properties for uniform coating deposition. Electrochemical characterization, including cyclic voltammetry and rate capability tests, unveils the electrochemical behavior of NCMA, LFMP, and their blend, showcasing the synergistic effects of material blending on battery performance. In conclusion, our findings underscore the intricate relationship between slurry formulation, electrode processing, and LIB performance, offering valuable insights for the design and optimization of high-performance electrode materials for next-generation batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

A layered solid finite element formulation with interlaminar enhanced displacements for the modeling of laminated composite structures

Accurate modeling of layered composite structures often requires the use of detailed finite element models which can sufficiently resolve the kinematics and material behavior within each layer of the composite. However, individually discretizing each material layer into finite elements presents a prohibitive computational expensive given the large number of thin layers comprising some laminated composites. To address these challenges, an 8-node layered solid hexahedral finite element is formulated with the aim of striking an appropriate balance between efficiency and fidelity. The element is discretized into an arbitrary number of distinct material layers, and employs reduced in-plane integration within each layer. The chosen reduced integration scheme is supplemented by a novel physical stabilization approach which includes layerwise enhancements to mitigate various forms of locking phenomena. The proposed framework additionally supports the inclusion of interlaminar enhanced displacements to better represent the kinematics of general layered composite materials. Finally, the described element formulation has been implemented in the ParaDyn finite element code, and its efficacy for modeling laminated composite structures is demonstrated on a variety of verification problems.

42 ENGINEERING↗
DOE OSTI2025

Development of a Printable Prill Formulation Technique and Demonstration of Monomodal Prill Size on Compaction Density and Compressive Strength

Polymer-bonded explosive molding powder, or “prills,” are relied on for the fabrication of pressed high explosives since the 1950's. The wet granulation technique, also known as “slurry coating,” that is used to formulate prills, is a complex process that results in polydisperse and variable yields. This makes it difficult to study the mesoscale effect that prills have on the microstructure of a pressed article. The following study introduces a novel approach to energetic granulation that leverages techniques used in the additive manufacturing of paste-like energetic materials. This extrusion granulation, or prill printing technique, makes it possible to tailor the sizes and shapes of prills, allowing for their morphological influences to be studied in a controlled manner. The following work details the fabrication and characterization of four monomodal size lots of prills using an inert formulation (95 wt.% melamine, 5 wt.% polymer binder). Prills from each size lot were die-pressed using a fixed recipe to investigate how prill size impacts compaction density and therefore compressive strength. It was found that larger prills influence the pressing density by creating larger defects within the microstructure of a pressed article, resulting in a decrease in compressive strength.

direct ink write↗
DOE OSTI2025

Gauge loop-string-hadron formulation on general graphs and applications to fully gauge fixed Hamiltonian lattice gauge theory

We develop a gauge invariant, Loop-String-Hadron (LSH) based representation of SU(2) Yang-Mills theory defined on a general graph consisting of vertices and half-links. Inspired by weak coupling studies, we apply this technique to maximal tree gauge fixing. This allows us to develop a fully gauge-fixed representation of the theory in terms of LSH quantum numbers. We explicitly show how the quantum numbers in this formulation directly relate to the variables in the magnetic description. In doing so, we will also explain in detail how the Kogut-Susskind formulation, prepotentials, and point splitting work for general graphs. In the appendix of this work, we provide a self-contained exposition of the mathematical details of Hamiltonian pure gauge theories defined on general graphs.

Algorithms and Theoretical Developments↗
DOE OSTI2025

A Particle-in-cell Method for Plasmas with A Generalized Momentum Formulation, Part III: A family of Gauge Conserving Methods

In this paper, we introduce a new family of spatially co-located field solvers for particle-in-cell applications which evolve the potential formulation of Maxwell’s equations under the Lorenz gauge. Our recent work [2] introduced the concept of time-consistency, which connects charge conservation to the preservation of the gauge at the semi-discrete level. It will be shown that there exists a large family of time discretizations which satisfy this property. Additionally, it will be further shown that for large classes of time marching methods, the satisfaction of the gauge condition automatically implies the satisfaction of Gauss’s law for electricity, with the potential formulation ensuring that that Gauss’s law for magnetism is satisfied by definition. We focus on popular time marching methods including centered differences, backward differences, and diagonally-implicit Runge-Kutta methods, which are coupled to a spectral discretization in space. We demonstrate the theory by testing the methods on a relativistic Weibel instability and a drifting cloud of electrons.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Efficient Computation of Doppler-Broadened Elastic Scattering Kernel Moments Using Ladder-Operator Formulation

Anefficient routine for computing Legendre moments of the Doppler-broadened elastic scattering kernel, including resonance scattering effects, has been implemented in the ISOXML module of Griffin. Isotropic scattering in the center-of-mass system and the ideal gas model for target motion are assumed. A ladder-operator formulation is introduced to compute all Legendre moments from order 0 to N simultaneously, enabling near-linear scaling of computational cost with respect to the maximum Legendre order. A physics-based strategy for constructing outgoing energy grids has also been developed, in which a tailored base grid is combined with adaptive refinement to maintain accuracy while limiting the number of outgoing energy points. For energies between resonances, a constant cross-section model is employed to further reduce computational cost. In addition, a quantitative criterion is derived to determine isotope-wise cut-off incident energies based on a prescribed up-scattering probability coverage. For 238U, up to incident energies of approximately 75, 230, and 661 eV at 294, 900, and 2500 K (corresponding to a 2% up-scattering probability threshold), computation of P0 kernels requires 1–8 s and computation of P0–P5 kernels requires 0.4–4 min using a single thread, while maintaining 1–3% relative error in up-scattering probability. These results demonstrate that the proposed formulation enables accurate and computationally practical Doppler-broadened kernel generation for online multigroup cross-section production in Griffin.

Doppler-broadening↗
DOE OSTI2026

An immersed interface method for microstructure-scale electrochemical battery models: numerical formulation and performance portable implementation

We present the numerical formulation, verification, and performance portable implementation of an immersed interface method for microstructure scale electrochemical modeling of batteries. The innovation in this approach is the resolution of chemical species and electrostatic potential discontinuities at active interfaces without the use of interface conforming unstructured grids. A unified formulation on Cartesian grids for all domains (electrodes and electrolyte) is used with interfacial flux conditions applied using volume fraction or “color” function gradients. We have developed one dimensional and two dimensional test cases with analytic solutions for electrochemical modeling using which we verified the consistency and accuracy of our scheme. Our solver is also validated against solutions from a macroscale model and an unstructured multi-subdomain solver for a full lithium ion cell. We then demonstrated the utility of our solver on an image-based complex battery electrode microstructure at high charging rate. Our technique also exhibits good scalability on distributed memory architectures using central processing units (CPU), with problem sizes up to 1.8 billion degrees of freedom and with 5400 ranks. Initial performance studies of our open-source performance portable solver showed about 70 times speed up using a graphics processing unit (GPU) compared to single compute core for a problem with 4 million cells.

25 ENERGY STORAGE↗