Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
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 OSTI2017

Finite element formulation with embedded weak discontinuities for strain localization under dynamic conditions

Here, we present an explicit finite element formulation designed for the treatment of strain localization under highly dynamic conditions. We also used a material stability analysis to detect the onset of localization behavior. Finite elements with embedded weak discontinuities are employed with the aim of representing subsequent localized deformation accurately. The formulation and its algorithmic implementation are described in detail. Numerical results are presented to illustrate the usefulness of this computational framework in the treatment of strain localization under highly dynamic conditions, and to examine its performance characteristics in the context of two-dimensional plane-strain problems.

36 MATERIALS SCIENCE↗
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 OSTI2019

Constrained Mixture Design to Optimize Formulation and Performance of Foams Based on Cassava Starch and Peanut Skin

Foams composed of cassava starch, peanut skin, and glycerol were developed by thermal expansion process. A constrained mixture design was applied to evaluate the effect of components on flexural mechanical properties and water absorption capacity (WAC). The interaction between cassava starch and peanut skin exhibited a significant antagonistic effect in WAC. Effects of the concentration of cassava starch, glycerol and the interactions between cassava starch–glycerol and peanut skin–glycerol were significant to mechanical properties. Desirability function was applied to optimize the formulation in order to obtain less hydrophilic, more flexible and more resistant foams. Optimized formulation consisted of 23.71% (w/w) peanut skin and 12.68% (w/w) glycerol, based on solids content (cassava starch + peanut skin). Optimized foam showed experimental values similar to those predicted, validating the fitted models for defined constraints. This foam showed an important improvement in hydrophilicity indicated by a decrease in WAC, compared to cassava starch-based foams without peanut skin.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Solubility Comparison of Cyclotetramethylenetetranitramine (HMX) and Potential Mock-HMX Candidates in Solvents Relevant to Slurry Formulating Polymer-Bonded Materials

Understanding solubility behavior of organic crystalline materials is critical in the development of slurry-formulated plastic bonded composites. In the field of high explosives (HE), a common method for obtaining high density, greater than 95 % theoretical maximum density, is to encapsulate HE in a polymer matrix resulting in the formation of plastic bonded explosive (PBX) prills. Prill formation is highly nuanced and requires understanding of particle size distributions, solvent, antisolvent, and binder system as well as interactions thereof. In this paper, in an effort to create mock-HE that mimics HE in density and mechanical strength, solubility of candidate molecules was compared to cyclotetramethylenetetranitramine (HMX) in binder solvents relevant to prill formation. The results indicate that the mock candidates chosen for this study have different solubilities compared to HMX in water and both organic solvents, 2-butanone and ethyl acetate, investigated. With these solubility differences known, prills to match PBXs based on HMX can be formulated, pressed and tested.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Formulation and Assessment of a Wash-Primer Containing Lanthanum “Tannate” for Steel Temporary Protection

Tannins are polyphenols synthesized by plants and useful for the coating industry as corrosion inhibitors. In addition, lanthanum salts have a great inhibitory effect on steel corrosion. The aim of this study was to obtain lanthanum “tannate” with adequate solubility to be incorporated as the corrosion inhibitor in a wash-primer. The “tannate” was obtained from commercial “Quebracho” tannin and 0.1 M La(NO{sub 3}){sub 3}. The soluble tannin was determined by the Folin–Denis reagent, while the concentration of Lanthanum was obtained by a gravimetric procedure. The protective action of “tannate” on SAE 1010 steel was evaluated by linear polarization curves and corrosion potential measurements. Lanthanum “tannate” was incorporated in a wash-primer formulation and tested by corrosion potential and ionic resistance measurements. The corrosion rate was also determined by the polarization resistance technique. Besides, the primer was incorporated in an alkyd paint system and its anticorrosion performance assessed in the salt spray cabinet and by electrochemical impedance spectroscopy. Results showed that lanthanum “tannate” primer inhibits the development of deleterious iron oxyhydroxides on the steel substrate and incorporated into a paint system had a similar behavior to the primer formulated with zinc tetroxychromate.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A Bi-Level Formulation And Solution Method For The Integration Of Process Design And Scheduling

Design and operation decisions of manufacturing plants can often be expressed within a hierarchical structure, where optimal decisions at lower operating level provide constraints for the decision raking at a higher decision level. In this work, we are addressing the case of a hierarchical design and scheduling optimization problem. The integration of design and scheduling decisions can play a big role in designing economically profitable plants and improving their operational performance. This problem can be expressed as a bi-level problem, where design related decisions occur at the upper level and operational scheduling decisions at the lower level. Since discrete decisions are involved in both optimization levels, the resulting formulation typically corresponds to bi-level mixed-integer programming problems (B-MIP). The solution of B-MIP problems is very challenging, and typically requires the use of global optimization techniques, with many algorithms not able to guarantee feasibility. To overcome the challenges in solving this class of problems we propose the use of a multi-parametric based algorithm for the solution of bi-level mixed-integer linear and quadratic programming problems, capable of providing the exact and global solution. The main idea of this approach is to treat the lower scheduling problem as a multi-parametric programming problem in which the design decisions are considered as parameters. Here, the resulting parametric solutions can then be substituted into the upper level design problem, which can be solved as a set of single-level deterministic programming problems. Through the developed formulation and algorithm we are able to supply the decision makers with the exact solution of the integrated design and scheduling problem.

42 ENGINEERING↗
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 OSTI2014

Affine group formulation of the Standard Model coupled to gravity

In this work we apply the affine group formalism for four dimensional gravity of Lorentzian signature, which is based on Klauder’s affine algebraic program, to the formulation of the Hamiltonian constraint of the interaction of matter and all forces, including gravity with non-vanishing cosmological constant Λ, as an affine Lie algebra. We use the hermitian action of fermions coupled to gravitation and Yang–Mills theory to find the density weight one fermionic super-Hamiltonian constraint. This term, combined with the Yang–Mills and Higgs energy densities, are composed with York’s integrated time functional. The result, when combined with the imaginary part of the Chern–Simons functional Q, forms the affine commutation relation with the volume element V(x). Affine algebraic quantization of gravitation and matter on equal footing implies a fundamental uncertainty relation which is predicated upon a non-vanishing cosmological constant. -- Highlights: •Wheeler–DeWitt equation (WDW) quantized as affine algebra, realizing Klauder’s program. •WDW formulated for interaction of matter and all forces, including gravity, as affine algebra. •WDW features Hermitian generators in spite of fermionic content: Standard Model addressed. •Constructed a family of physical states for the full, coupled theory via affine coherent states. •Fundamental uncertainty relation, predicated on non-vanishing cosmological constant.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

A polymeric nanoparticle formulation of curcumin in combination with sorafenib synergistically inhibits tumor growth and metastasis in an orthotopic model of human hepatocellular carcinoma

Curcumin, a yellow polyphenol extracted from the rhizome of turmeric root (Curcuma longa) has potent anti-cancer properties in many types of tumors with ability to reverse multidrug resistance of cancer cells. However, widespread clinical application of this agent in cancer and other diseases has been limited due to its poor aqueous solubility. The recent findings of polymeric nanoparticle formulation of curcumin (NFC) have shown the potential for circumventing the problem of poor solubility, however evidences for NFC's anti-cancer and reverse multidrug resistance properties are lacking. Here we provide models of human hepatocellular carcinoma (HCC), the most common form of primary liver cancer, in vitro and in vivo to evaluate the efficacy of NFC alone and in combination with sorafenib, a kinase inhibitor approved for treatment of HCC. Results showed that NFC not only inhibited the proliferation and invasion of HCC cell lines in vitro, but also drastically suppressed primary tumor growth and lung metastases in vivo. Moreover, in combination with sorafenib, NFC induced HCC cell apoptosis and cell cycle arrest. Mechanistically, NFC and sorafenib synergistically down-regulated the expression of MMP9 via NF-κB/p65 signaling pathway. Furthermore, the combination therapy significantly decreased the population of CD133-positive HCC cells, which have been reported as cancer initiating cells in HCC. Taken together, NanoCurcumin provides an opportunity to expand the clinical repertoire of this agent. Additional studies utilizing a combination of NanoCurcumin and sorafenib in HCC are needed for further clinical development. - Highlights: • Polymeric nanoparticle formulation of curcumin not only inhibited the proliferation and invasion of HCC cell lines in vitro, but also drastically suppressed primary tumor growth and lung metastases in vivo. • In combination with sorafenib, NanoCurcumin induced HCC cell apoptosis and cell cycle arrest. • NanoCurcumin and sorafenib synergistically down-regulated the expression of MMP9 via NF-κB/p65 signaling pathway. • NanoCurcumin and sorafenib significantly decreased the population of CD133-positive HCC cells.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Formulation and numerical analysis of a fully-coupled dynamically deforming electromagnetic wire

Here, an electromagnetic beam model is developed for the simulation of actuated electronic textiles. The beam is solved using a nonlinear director-based kinematic description with additional temperature and electric potential fields along its length. The three fields are fully coupled by mutual dependences on the deformation, Lorenz force, back electromotive force, temperature dependent constitutive responses, and the Seebeck effect. Instead of solving Maxwell's equations in full detail, a quasistatic approximation is used to solve the electric potential in the presence of a moving material medium. The current-carrying beam approximation is used to further simplify the solution space for the potential. While this formulation alleviates the spatial and temporal discretization restrictions, the coupled problem is an index-1 semi-explicit Differential Algebraic Equation requiring special treatment. The time dependent problem is solved using different Runge-Kutta methods. Diagonally implicit Runge-Kutta methods and explicit Runge-Kutta methods using implicit solution of the electric potential problem are explored. The finite element model is implemented using the open source package FEniCS, which is able to automatically generate the linearizations of the multiphysics equations required for the implicit solutions. A model problem is constructed with which to test and analyze the physical formulation and numerical solution techniques. The time stepping methods are verified using the convergence orders of the higher-order Runge-Kutta methods. Runtime comparisons show that the explicit methods are generally more computationally efficient than the implicit schemes used for this problem. For the implicit schemes, a staggered solution is significantly faster than a monolithic solution at most time step sizes. However, at very large time steps, such as those that would be used for dynamic relaxation, the monolithic solution can be more efficient than the staggered solution.

36 MATERIALS SCIENCE↗
DOE OSTI2019

An isogeometric finite element formulation for phase transitions on deforming surfaces

This paper presents a general theory and isogeometric finite element implementation for studying mass conserving phase transitions on deforming surfaces. The mathematical problem is governed by two coupled fourth-order nonlinear partial differential equations (PDEs) that live on an evolving two-dimensional manifold. For the phase transitions, the PDE is the Cahn–Hilliard equation for curved surfaces, which can be derived from surface mass balance in the framework of irreversible thermodynamics. For the surface deformation, the PDE is the (vector-valued) Kirchhoff–Love thin shell equation. Both PDEs can be efficiently discretized using C 1 -continuous interpolations without derivative degrees-of-freedom (dofs). Structured NURBS and unstructured spline spaces with pointwise C 1 -continuity are utilized for these interpolations. The resulting finite element formulation is discretized in time by the generalized-α scheme with adaptive time-stepping, and it is fully linearized within a monolithic Newton–Raphson approach. A curvilinear surface parameterization is used throughout the formulation to admit general surface shapes and deformations. The behavior of the coupled system is illustrated by several numerical examples exhibiting phase transitions on deforming spheres, tori and double-tori.

42 ENGINEERING↗
DOE OSTI2019

Three-dimensional explicit finite element formulation for shear localization with global tracking of embedded weak discontinuities

Shear localization is a damage mechanism characterized by the concentration of plastic deformation within thin bands of materials, often observed as a precursor to failure in metals. In this study, we present a three-dimensional explicit finite element formulation with embedded weak discontinuities for the treatment of shear localization under dynamic loading conditions. In this formulation, the onset of localization is detected via a material stability analysis suitable for rate-sensitive materials, and the continuity of propagating shear bands is ensured using a global tracking strategy involving a heat-conduction type boundary value problem, solved for a scalar level set function over the global domain of the problem. In addition, we propose a modified quadrature rule that is used to compute the contribution of an individual element to the global finite element arrays, taking into account the position of the embedded shear band within that element.

42 ENGINEERING↗
DOE OSTI2021

Galerkin formulations of isogeometric shell analysis: Alleviating locking with Greville quadratures and higher-order elements

Here we propose new quadrature schemes that asymptotically require only four in-plane points for Reissner–Mindlin shell elements and nine in-plane points for Kirchhoff–Love shell elements in B-spline and NURBS-based isogeometric shell analysis, independent of the polynomial degree p of the elements. The quadrature points are Greville abscissae associated with pth-order B-spline basis functions whose continuities depend on the specific Galerkin formulations, and the quadrature weights are calculated by solving a linear moment fitting problem in each parametric direction. The proposed shell element formulations are shown through numerical studies to be rank sufficient and to be free of spurious modes. The studies reveal comparable accuracy, in terms of both displacement and stress, compared with fully integrated spline-based shell elements, while at the same time reducing storage and computational cost associated with forming element stiffness and mass matrices and force vectors. The high accuracy with low computational cost makes the proposed quadratures along with higher-order spline bases, in particular polynomial orders, p = 5 and 6, good choices for alleviating membrane and shear locking in shells.

42 ENGINEERING↗
DOE OSTI2021

A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems

This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

On a fully-implicit VMS-stabilized FE formulation for low Mach number compressible resistive MHD with application to MCF

This study presents the development and evaluation of a fully-implicit variational multiscale (VMS) stabilized unstructured finite element (FE) formulation for compressible magnetohydrodynamics (MHD) model, at low Mach number regime. The model describes the dynamics of a compressible conducting fluid in the low Mach number limit in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics phenomena. The specific applications that motivate this study are macroscopic simulations of the longer time-scale stability and disruptions of magnetic confinement fusion (MCF) devices, specifically the ITER tokamak. The discussion considers the development of the VMS FE representation, the structure of the stabilizing terms that deal with significant convective flows, the stabilization of the nearly incompressible response of the fluid flow, and the stabilization of the constraint that enforces the solenoidal involution on the magnetic field. The nonlinear discretized system is solved with scalable preconditioned Newton–Krylov iterative methods, which employs a multiphysics block preconditioning method based on approximate block factorizations and Schur complements. The study presents an evaluation of the VMS method on a 2D cartesian tearing mode instability, and illustrates the scalability of the solvers on MCF relevant problems. A set of results are also presented for longer time-scale stability and disruptions for the ITER tokamak. These include a vertical displacement event (VDE), and a (1,1) internal kink mode. Here, the formulation is demonstrated to be scalable and also reasonably robust with respect to the Lundquist number scaling.

42 ENGINEERING↗
DOE OSTI2024

Pathways to formulate lightweight and ultra-lightweight 3D printable cementitious composites

Here, this paper studies the pathways to formulate lightweight and ultra-lightweight 3D printable cementitious composites. A hybrid approach was proposed by combining the advantages of traditional chemical-induced foaming (foaming approach) and lightweight particulate inclusions (synthetic foam approach). A comprehensive experimental program was conducted to evaluate the effects of foaming agents and fly ash cenosphere (FAC) on the printability, microstructure, mechanical and thermal properties of 3D printed samples. The results showed that the hybrid approach could produce a mixture with a density as low as 470 kg/m 3 while ensuring good flowability and buildability owing to the lubricating effect of foaming and supporting skeleton formed by FAC. In addition, a three-step homogenization procedure was also developed to predict the effective elastic modulus and thermal conductivity of 3D printable cementitious composites and cementitious foam. The findings of the study highlighted the effectiveness of the hybrid approach in formulating 3D printable ultra-lightweight cementitious composites in thermal insulation and acoustic applications.

36 MATERIALS SCIENCE↗
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↗