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 37 records · Page 2

Feasibility of Formulating Ecosystem Biogeochemical Models From Established Physical Rules

Abstract To improve the predictive capability of ecosystem biogeochemical models (EBMs), we discuss the feasibility of formulating biogeochemical processes using physical rules that have underpinned the many successes in computational physics and chemistry. We argue that the currently popular empirically based approaches, such as multiplicative empirical response functions and the law of the minimum, will not lead to EBM formulations that can be continuously refined to incorporate improved mechanistic understanding and empirical observations of biogeochemical processes. Instead, we propose that EBM parameterizations, as a lossy data compression problem, can be better formulated using established physical rules widely used in computational physics and chemistry, and different biogeochemical processes can be more robustly integrated within a reactive‐transport framework. Through several examples, we demonstrate how mathematical representations derived from physical rules can improve understanding of relevant biogeochemical processes and enable more effective communication between modelers, observationalists, and experimentalists regarding essential questions, such as what measurements are needed to meaningfully inform models and how can models generate new process‐level hypotheses to test in empirical studies. Finally, while empirical models with more parameters are often less robust, physical rules‐based models can be more robust and show lower predictive equifinality, stemming from their enhanced consistency in representations of processes, interactions and spatial scaling.

54 ENVIRONMENTAL SCIENCES↗

Phase field dislocation dynamics formulation coupled with Fourier based micromechanics solver and its application to grain boundary–dislocation interactions

A new phase field dislocation dynamics (PFDD) formulation for homogeneous and heterogeneous materials is presented, which couples micromechanical solvers and the time-dependent Ginzburg–Landau equation. The strain fields are obtained from the micromechanical solver by solving the Lippmann–Schwinger equation and then used to define energy terms to model the evolution of the dislocations. Grain boundary (GB)–dislocation interactions are studied using the coupled PFDD formulation and by describing GBs as inclusions. GB energy and stiffness tensors are computed from molecular statics simulations, and a newly proposed lattice energy term that is dependent on the GB energy is considered in the calculations. Interaction of a screw dislocation with minimum energy and metastable states of low and high angle ⟨110⟩ symmetric tilt grain boundaries are studied. We show good agreement between predictions from our PFDD formulation and molecular dynamics simulations of grain boundary–dislocation interactions.

36 MATERIALS SCIENCE↗

Galerkin formulation of path integrals in lattice field theory

We present a mathematical framework for Galerkin formulations of path integrals in lattice field theory. The framework is based on using the degrees of freedom (DOFs) associated to a Galerkin discretization as the fundamental lattice variables. We formulate standard concepts in lattice field theory, such as the partition function and correlation functions, in terms of the DOFs. For example, using continuous finite element spaces, we show that the two-point spatial correlation function can be defined between any two points on the domain (as opposed to at just lattice sites) and furthermore, this two-point function satisfies a weak propagator (or Green’s function) identity, in analogy to the continuum case, as well as a convergence estimate obtained from the standard finite element techniques. Furthermore, this framework leads naturally to higher-order formulations of lattice field theories by considering higher-order finite element spaces for the Galerkin discretization. We consider analytical and numerical examples of scalar field theory to investigate how increasing the order of piecewise polynomial finite element spaces affect the approximation of lattice observables. Finally, we sketch an outline of this Galerkin framework in the context of gauge field theories.

97 MATHEMATICS AND COMPUTING↗

Pseudospectral particle-in-cell formulation with arbitrary charge and current-density time dependencies for the modeling of relativistic plasmas

This paper introduces a formulation of the particle-in-cell (PIC) method for the modeling of relativistic plasmas, that leverages the ability of the pseudospectral analytical time-domain solver (PSATD) to handle arbitrary time dependencies of the charge and current densities during one PIC cycle (applied to second-order polynomial dependencies here). Here, the formulation is applied to a modified set of Maxwell's equations that was proposed earlier in the context of divergence cleaning, and to recently proposed extensions of the PSATD-PIC algorithm. Detailed analysis and testings revealed that, under some condition, the formulation can expand the range of numerical parameters under which PIC simulations are stable and accurate when modeling relativistic plasmas such as, e.g., plasma-based particle accelerators.

43 PARTICLE ACCELERATORS↗

Curing inhibition and shelf stability in silicone formulations for direct ink writing

Various potential hydrosilylation inhibitors were incorporated into a silicone formulation for 3D printing. Different inhibitors resulted in peak cure temperatures from 42.8 to 180.3°C. The curing behavior of formulations with 1-ethynyl-1-cyclohexanol (ETCH), which exhibited a peak cure temperature of 98.4°C, was characterized for different ETCH concentrations and storage conditions. The curing was modeled using time temperature superposition, revealing an activation barrier of 89 kJ/mol. Overall, characterization resulted in a conservative recommended shelf life of 90 days. In conclusion, these studies offer insights into the stability of these formulations and the importance of inhibition for producing reliable, printable inks.

36 MATERIALS SCIENCE↗

Development of Magnesium Oxysulfate Formulation for SRPPF Aqueous Recovery Liquid Solidification

The liquid effluent from the Savannah River Plutonium Processing Facility (SRPPF) Aqueous Recovery Processes will be solidified into a stable form that is acceptable by Waste Isolation Pilot Plant (WIPP) for disposal. The current Aqueous Recovery flow sheet proposes to solidify the liquid effluent using a grout formula that was developed and tested for the former Waste Solidification Building process. This Portland cement based mixture results in a high pH (~13) leachate from the solidified waste form which is not acceptable to WIPP in the large quantities expected from production at SRPPF. Various cementitious materials were previously evaluated as alternative grout formulations to Portland cement and a MgO-based mix was identified as a promising alternative. A magnesium oxysulfate (MOS) cement formulation comprised of reactive magnesium oxide (MgO), anhydrous magnesium sulfate (MgSO 4 ), and sand, as a non-reactive heat sink provided good mixability, similar density to the original Portland-cement based mix, and a leachate pH of 9.4, within the assumed WIPP brine pH range. However, the MOS formulation exhibited an appreciable amount of heat generation, which resulted in premature setting of a large-scale test.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Optimizing high energy density sulfur cathodes: A multivariate approach to electrode formulation and processing

Lithium-sulfur (Li-S) batteries involve complex solid-liquid-solid phase transformations during both discharging and charging processes, where cathode materials, formulation, and structure play a crucial role. Here, a design of experiments (DoE) methodology and an empirical model are developed to systematically explore the interactions and trade-offs among cathode factors and process variables, and to obtain generalizable effects estimates for the multivariate system. Compared to the conventional one-factor-at-a-time (OFAT) approach, this work demonstrates advantages in both efficiency and accuracy by allowing the data to guide future research and decisions. Further, an optimized cathode formulation and processing parameters are predicted and validated experimentally, achieving over 1000 mAh g -1 in discharge capacity and improved cycling under practical lean electrolyte (4 µL mg -1 S) and high S-loading cathodes (>4 mg cm -2 ) conditions. The optimized cathode was scaled up and assembled into Li-S pouch cells, achieving 316 Wh kg -1 in cell-level energy, proving that the comprehensive and rigorous framework for optimizing complex systems with DoE leads to improved performance in a practical pouch cell system.

25 ENERGY STORAGE↗

Anti-symmetric and positivity preserving formulation of a spectral method for Vlasov-Poisson equations

We analyze the anti-symmetric properties of a spectral discretization for the one-dimensional Vlasov-Poisson equations. The discretization is based on a spectral expansion in velocity with the symmetrically weighted Hermite basis functions, central finite differencing in space, and an implicit Runge Kutta integrator in time. The proposed discretization preserves the anti-symmetric structure of the advection operator in the Vlasov equation, resulting in a stable numerical method. We apply such discretization to two formulations: the canonical Vlasov-Poisson equations and their continuously transformed square-root representation. The latter preserves the positivity of the particle distribution function. We derive analytically the conservation properties of both formulations, including particle number, momentum, and energy, which are verified numerically on the following benchmark problems: manufactured solution, linear and nonlinear Landau damping, two-stream instability, bump-on-tail instability, and ion-acoustic wave.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Minimized aging of isocyanurate-based rigid cellular foams for buildings through tailored barrier facers and optimized formulation

The thermal resistivity (h·ft 2 ·°F/Btu/in.) of closed-cell rigid foam insulation materials significantly decreases over time. Diffusion-tight facers are designed to significantly enhance initial thermal resistivity and long-term thermal performance by preventing gas diffusion into the foam cells and inhibiting the escape of low thermal conductivity blowing agents. In addition to the gas diffusion property of the facer film, adhesion between the foam and the facer is crucial for achieving diffusion-tight bonding. Here, this study addresses the challenge of thermal aging by investigating how facer film properties and foam formulation influence the durability of thermal performance. A systematic evaluation was conducted to understand the effects of polymeric barrier films, surface treatments, facer coverage, and foam matrix rigidity on thermal resistivity over time. Key findings reveal that diffusion-tight facers, particularly those with metallized layers and compatible heat seal layers, significantly reduce gas exchange and improve foam-facer adhesion. The optimized system, incorporating barrier facers and a tailored polyisocyanurate foam formulation, achieved initial and aged thermal resistivity values after 200 days of approximately 8.3 and 7.4 h ft 2 ·°F/Btu/in., respectively representing only a 10 % reduction compared to a 17 % reduction observed in control samples without facers. Notably, polyurethane spray foams with facers exhibited only a 4 % reduction in thermal resistivity, compared to a 23 % decrease in control samples, demonstrating nearly six times better retention of thermal performance. This innovative facer technology presents a promising solution for reducing energy costs and represents a significant advancement in optimizing energy management for building envelopes in future technologies. This technology can also be adapted for other applications that necessitate the preservation of long-term thermal performance.

Wanasinghe, Shiwanka Vidarshi [Oak Ridge National ↗

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗