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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Time-stepping DPG formulations for the heat equation

For a wide range of PDEs, the discontinuous Petrov–Galerkin (DPG) methodology of Demkowicz and Gopalakrishnan provides discrete stability starting from a coarse mesh and minimization of the residual in a user-controlled norm, among other appealing features. Research on DPG for transient problems has mainly focused on spacetime discretizations, which has theoretical advantages, but practical costs for computations and software implementations. The sole examination of time-stepping DPG formulations was performed by Führer, Heuer, and Gupta, who applied Rothe’s method to an ultraweak formulation of the heat equation to develop an implicit time-stepping scheme; their work emphasized theoretical results, including error estimates in time and space. Here, we follow Führer, Heuer, and Gupta in examining the heat equation; our focus is on numerical experiments, examining the stability and accuracy of several formulations, including primal as well as ultraweak, and explicit as well as implicit and Crank–Nicolson time-stepping schemes. We are additionally interested in communication-avoiding algorithms, and we therefore include a highly experimental formulation that places all the trace terms on the right-hand side of the equation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The impacts of convex piecewise linear cost formulations on AC optimal power flow

Despite strong connections through shared application areas, research efforts on power market optimization (e.g., unit commitment) and power network optimization (e.g., optimal power flow) remain largely independent. A notable illustration of this is the treatment of power generation cost functions, where nonlinear network optimization has largely used polynomial representations and market optimization has adopted piecewise linear encodings. This work combines state-of-the-art results from both lines of research to understand the best mathematical formulations of the nonlinear AC optimal power flow problem with piecewise linear generation cost functions. An extensive numerical analysis of non-convex models, linear approximations, and convex relaxations across fifty-four realistic test cases illustrates that nonlinear optimization methods are surprisingly sensitive to the mathematical formulation of piecewise linear functions. The results indicate that a poor formulation choice can slow down algorithm performance by a factor of ten, increasing the runtime from seconds to minutes. Furthermore, these results provide valuable insights into the best formulations of nonlinear optimal power flow problems with piecewise linear cost functions, an important step towards building a new generation of energy markets that incorporate the nonlinear AC power flow model.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2016

A new efficient formulation of the HLLEM Riemann solver for general conservative and non-conservative hyperbolic systems

In this paper a new, simple and universal formulation of the HLLEM Riemann solver (RS) is proposed that works for general conservative and non-conservative systems of hyperbolic equations. For non-conservative PDE, a path-conservative formulation of the HLLEM RS is presented for the first time in this paper. The HLLEM Riemann solver is built on top of a novel and very robust path-conservative HLL method. It thus naturally inherits the positivity properties and the entropy enforcement of the underlying HLL scheme. However, with just the slight additional cost of evaluating eigenvectors and eigenvalues of intermediate characteristic fields, we can represent linearly degenerate intermediate waves with a minimum of smearing. For conservative systems, our paper provides the easiest and most seamless path for taking a pre-existing HLL RS and quickly and effortlessly converting it to a RS that provides improved results, comparable with those of an HLLC, HLLD, Osher or Roe-type RS. This is done with minimal additional computational complexity, making our variant of the HLLEM RS also a very fast RS that can accurately represent linearly degenerate discontinuities. Our present HLLEM RS also transparently extends these advantages to non-conservative systems. For shallow water-type systems, the resulting method is proven to be well-balanced. Several test problems are presented for shallow water-type equations and two-phase flow models, as well as for gas dynamics with real equation of state, magnetohydrodynamics (MHD & RMHD), and nonlinear elasticity. Since our new formulation accommodates multiple intermediate waves and has a broader applicability than the original HLLEM method, it could alternatively be called the HLLI Riemann solver, where the “I” stands for the intermediate characteristic fields that can be accounted for. -- Highlights: •New simple and general path-conservative formulation of the HLLEM Riemann solver. •Application to general conservative and non-conservative hyperbolic systems. •Inclusion of sub-structure and resolution of intermediate characteristic fields. •Well-balanced for single- and two-layer shallow water equations and multi-phase flows. •Euler equations with real equation of state, MHD equations, nonlinear elasticity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A regularized-interface method as a unified formulation for simulations of high-pressure multiphase flows

The injection of multi-species fluids into high-pressure and high-temperature environments beyond the species' critical points is commonly found in engineering applications. At these conditions, for immiscible species, both subcritical interfacial dynamics and supercritical mixing can coexist due to variations in temperature around the mixture critical point. The modeling of these complex transcritical phenomena for large-scale configurations is so far not possible. To address this issue, we propose the Regularized-Interface Method (RIM) as a unified formulation that can describe both sub- and supercritical processes as well as the transition between them. The proposed method is derived via filtering of the nanoscale interface-resolving formulation based on van der Waals' linear gradient theory. Thus, this approach allows for the consistent modeling of interfacial dynamics that vanishes at supercritical conditions, while significantly reducing the temporal and spatial resolution constraints of the original nanoscale formulation. The resulting RIM formulation is examined in interface-capturing simulations of sub-, trans-, and supercritical fuel injection processes, involving droplets and jets. Furthermore, these results highlight the importance of resolving spatio-temporal transitions from subcritical interfacial dynamics to supercritical mixing in high-pressure multiphase simulations, in contrast to commonly employed diffused-interface methods, where interfacial dynamics are often neglected.

Interface capturing↗
DOE OSTI2026

Synergistic cellulase–xylanase formulations for enhanced dewatering and fiber bonding toward energy-efficient and sustainable paper and packaging production

A mechanistic understanding of the synergistic effects of enzymes on cellulosic fibers dewatering and fiber bonding is essential for advancing energy-efficiency and lightweight production of paper and packaging materials. This study investigates the impact of varying cellulase and xylanase formulations on equilibrium moisture content (EMC) after pressing and tensile strength of cellulosic fiber webs using a factorial experimental design. Nine custom enzyme formulations were evaluated at controlled dosages ranging from ∼20 to 76 FPU/mL cellulase and ∼500–1130 IU/mL xylanase. Response surface modeling revealed a significant synergy, particularly at 20–40 FPU/mL cellulase combined with ≥ 1000 IU/mL xylanase. Under these conditions, EMC decreased by up to 3.6% compared with the untreated refined control, while tensile index gains of up to 16% were statistically significant for optimized blends (p < 0.05). The interaction between cellulase and xylanase was also significant for the tensile response (p = 0.010). Protein efficiency analysis showed that optimized formulations containing 30–40% less protein outperformed the commercial benchmark. The nonlinear synergy between cellulase and xylanase is attributed to their complementary substrate specificities. Endoglucanase- and β‑glucosidase‑rich cellulases hydrolyze internal β‑1,4‑glycosidic bonds in amorphous cellulose, loosening fiber walls and increasing flexibility, while xylanases target hemicellulose, primarily xylan-rich domains, enhancing porosity and improve cellulase accessibility. Tailoring enzyme formulations at low loadings overcomes traditional trade-offs between strength and dewatering, enabling cost-effective, energy-efficient, low-carbon solutions for sustainable packaging and hygiene products.

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

New large-strain FFT-based formulation and its application to model strain localization in nano-metallic laminates and other strongly anisotropic crystalline materials

This paper presents a new robust large-strain (LS) elasto-viscoplastic (EVP) formulation based on Fast Fourier Transforms (FFTs) for the prediction of the micro-mechanical response and microstructure evolution of polycrystalline and multiphase materials, with emphasis on the effect of strong crystallographic and/or morphologic anisotropy on localization of plastic deformation. In this work, the novel LS-EVPFFT formulation allows treatment of complex initial geometries and large deformations considering three grids of material points: a regular grid in the reference configuration, where FFTs can be performed; an irregular grid in the initial configuration, created by applying a stress-free displacement field to the reference regular grid; and an irregular grid in the current configuration, undergoing large strains and rotations as the material is loaded. Further numerical stability of the new formulation also required the use of a novel expression for the discrete modified Green’s operator, which reduces spurious field oscillations. After presenting and validating the new formulation by comparison with preexisting implementations and analytical solutions, LS-EVPFFT is applied to the prediction of slip and kink bands formation in polycrystalline columnar ice, and kink bands in single crystal zinc wires, showing good agreement with classic experiments. Finally, the model is used to study kink band formation during compression of Cu–Nb nano-metallic laminates (NMLs), in which accurate treatment of the complex geometry associated with the tortuosity of interfaces and large deformations become critical, showing consistency with corresponding micropillar experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Piecewise polyhedral formulations for a multilinear term

Herein, we present a mixed-integer linear programming (MILP) formulation of a piecewise, polyhedral relaxation (PPR) of a multilinear term using its convex-hull representation. Based on the PPR’s solution, we also present a MILP formulation whose solutions are feasible for nonconvex, multilinear equations. We then present computational results showing the effectiveness of proposed formulations on standard benchmark nonlinear programs (NLPs) with multilinear terms and compare with a traditional formulation that is built using recursive bilinear groupings of multilinear terms.

42 ENGINEERING↗
DOE OSTI2023

On formulations for modeling pressurized cracks within phase-field methods for fracture

Over the past few decades, the phase-field method for fracture has seen widespread appeal due to the many benefits associated with its ability to regularize a sharp crack geometry. Along the way, several different models for including the effects of pressure loads on the crack faces have been developed. Further, this work investigates the performance of these models and compares them to a relatively new formulation for incorporating crack-face pressure loads. It is shown how the new formulation can be obtained either by modifying the trial space in the traditional variational principle or by postulating a new functional that is dependent on the rates of the primary variables. The key differences between the new formulation and existing models for pressurized cracks in a phase-field setting are highlighted. Model-based simulations developed with discretized versions of the new formulation and existing models are then used to illustrate the advantages and differences. In order to analyze the results, a domain form of the J-integral is developed for diffuse cracks subjected to pressure loads. Results are presented for a one-dimensional cohesive crack, steady crack growth, and crack nucleation from a pressurized enclosure.

42 ENGINEERING↗
DOE OSTI2025

Data-Efficient Methods for Determining Flory–Huggins χ Parameters in Multicomponent Polymer Formulations

Polymer formulations are essential in diverse applications including personal care products, coatings, paints, adhesives, and plastic materials. Designing these formulations requires navigating large, complex design spaces, where phase and self-assembly behavior critically impact performance. The Flory–Huggins χ parameter, which quantifies segmental miscibility, is widely used to parametrize the excess free energy of mixing in formulation models. In this work, we introduce two data-efficient, top-down methods for estimating χ parameters using the Random Phase Approximation (RPA): (i) Boundary Nonlinear Regression (Boundary-NLR), which fits theoretical spinodal boundaries to experimental phase boundaries, and (ii) Surrogate Model Inverse Parameter Estimation (SMIPE), which uses a Gaussian Process Classifier to fit sparse phase maps via a surrogate model. Both methods allow rapid parametrization of polymer field-theoretic models without the need for additional experiments. We evaluate these approaches on data sets involving polymer–solvent–nonsolvent ternary mixtures and block copolymer–solvent systems, demonstrating their robustness to experimental noise and their relevance for real-world formulation design.

copolymers↗
DOE OSTI2021

QUBO formulations for training machine learning models

Abstract Training machine learning models on classical computers is usually a time and compute intensive process. With Moore’s law nearing its inevitable end and an ever-increasing demand for large-scale data analysis using machine learning, we must leverage non-conventional computing paradigms like quantum computing to train machine learning models efficiently. Adiabatic quantum computers can approximately solve NP-hard problems, such as the quadratic unconstrained binary optimization (QUBO), faster than classical computers. Since many machine learning problems are also NP-hard, we believe adiabatic quantum computers might be instrumental in training machine learning models efficiently in the post Moore’s law era. In order to solve problems on adiabatic quantum computers, they must be formulated as QUBO problems, which is very challenging. In this paper, we formulate the training problems of three machine learning models—linear regression, support vector machine (SVM) and balanced k-means clustering—as QUBO problems, making them conducive to be trained on adiabatic quantum computers. We also analyze the computational complexities of our formulations and compare them to corresponding state-of-the-art classical approaches. We show that the time and space complexities of our formulations are better (in case of SVM and balanced k-means clustering) or equivalent (in case of linear regression) to their classical counterparts.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Importance of microstructural features in mechanical response of cast-cured HMX formulations

We report that cast-cure formulations of plastic-bonded explosives (PBX) present some advantages in manufacturing and application compared to traditional pressed slurry formulated materials. For example, these formulations can be cast to specific shapes and then cured in place, avoiding the need for machining precision parts. However, the microstructure of these materials can be greatly affected by the specific manufacturing process. Here, we evaluate the effect of minor changes to the formulation and manufacturing process on several cyclotetramethylene-tetranitramine (HMX) PBXs. The binders were based on hydroxyl-terminated polybutadiene (HTPB) and cast-cured using diphenylmethane diisocyanate (isonate) as a curing agent. We examined the materials using X-ray computed tomography (CT) imaging and uniaxial compression testing. The isonate content was found to significantly affect the modulus and strength of the binder. The presence of significant void content could be controlled by adding a centrifuging step during the curing process, but the resulting effect on mechanical properties was relatively minor. Finally, mesoscale simulations showed that differences in the mechanical strength of the binder were not sufficient to describe the differences observed in mechanical testing, indicating that the HMX-binder adhesion was also being changed by the manufacturing process.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Novel PBX formulations containing thermally-expandable microspheres for on-demand control of explosive behavior

In this work, we present the formulation and analysis of inert Plastic-Bonded eXplosive (PBX)-surrogates loaded with thermally expandable microspheres (TEMs). TEMs consist of a thermoplastic acrylonitrile shell (10-50 µm diameter) encapsulating an inert low boiling hydrocarbon. Upon heating, the TEMs expand as the shell softens while the hydrocarbon gasifies, increasing the internal pressure and expanding the particle by as much as 120 vol%. We hypothesize that TEM expansion within a PBX will introduce tunable changes in local density and porosity, and ultimately the shock sensitivity. This paper focuses on microstructural analysis of surrogate PBX-TEM formulations, both before and after thermally- induced TEM expansion. The formulations were analyzed by scanning electron microscopy (SEM), ultra-small angle neutron scattering (USANS), and X-ray Computed Tomography (x-ray CT) under various thermal conditions. Evidence suggests that TEM expansion within a PBX is a viable method to impart microstructural changes that are predicted to have a measurable, on-demand effect on the initiation sensitivity. Future work includes development of PBX-TEM formulations with HMX and shock initiation experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Real-space formulation of the stress tensor for O (N) density functional theory: Application to high temperature calculations

We present an accurate and efficient real-space formulation of the Hellmann–Feynman stress tensor for O (N) Kohn–Sham density functional theory (DFT). While applicable at any temperature, the formulation is most efficient at high temperature where the Fermi–Dirac distribution becomes smoother and the density matrix becomes correspondingly more localized. We first rewrite the orbital-dependent stress tensor for real-space DFT in terms of the density matrix, thereby making it amenable to O (N) methods. We then describe its evaluation within the O (N) infinite-cell Clenshaw–Curtis Spectral Quadrature (SQ) method, a technique that is applicable to metallic and insulating systems, is highly parallelizable, becomes increasingly efficient with increasing temperature, and provides results corresponding to the infinite crystal without the need of Brillouin zone integration. We demonstrate systematic convergence of the resulting formulation with respect to SQ parameters to exact diagonalization results and show convergence with respect to mesh size to the established plane wave results. Here, we employ the new formulation to compute the viscosity of hydrogen at 10 6 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

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

Benchmarking NIMROD continuum kinetic formulations through the steady-state poloidal flow

In this work, continuum kinetic formulations are employed as a mechanism to include closure physics in an extended magnetohydrodynamics model. Two continuum kinetic approaches have been implemented in the plasma fluid code NIMROD [Sovinec et al., “Nonlinear magnetohydrodynamics with high-order finite elements,” J. Comput. Phys. 195, 355 (2004)] including a Chapman–Enskog-like (CEL) formulation and a more conventional df approach. Ion kinetic closure schemes are employed to describe the neoclassical flow properties in axisymmetric toroidal geometry. In particular, predictions for steady-state values of poloidal flow profiles in tokamak geometry are provided using both the df formulation and two different solution techniques for the CEL approach. These results are benchmarked against analytic theory predictions as well as results from the drift kinetic code DK4D. The continuum kinetic formulations employed here show agreement with both the analytic theory and DK4D results, and offer a novel velocity space representation involving higher-order finite elements in pitch angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

THE EFFECT OF MAGNETIC TOPOLOGY ON THERMALLY DRIVEN WIND: TOWARD A GENERAL FORMULATION OF THE BRAKING LAW

Stellar wind is thought to be the main process responsible for the spin down of main-sequence stars. The extraction of angular momentum by a magnetized wind has been studied for decades, leading to several formulations for the resulting torque. However, previous studies generally consider simple dipole or split monopole stellar magnetic topologies. Here we consider, in addition to a dipolar stellar magnetic field, both quadrupolar and octupolar configurations, while also varying the rotation rate and the magnetic field strength. Sixty simulations made with a 2.5D cylindrical and axisymmetric set-up, and computed with the PLUTO code, were used to find torque formulations for each topology. We further succeed to give a unique law that fits the data for every topology by formulating the torque in terms of the amount of open magnetic flux in the wind. We also show that our formulation can be applied to even more realistic magnetic topologies, with examples of the Sun in its minimum and maximum phases as observed at the Wilcox Solar Observatory, and of a young K-star (TYC-0486-4943-1) whose topology has been obtained by Zeeman-Doppler Imaging.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Benefits of Sparse Tableau Over Nodal Admittance Formulation for Power-Flow Studies

Assembly of power-flow equations has traditionally begun from a nodal analysis formulation of the underlying transmission circuit's behavior. Most power-flow formulations encapsulate network constraints in the bus admittance matrix, Y bus . From a circuit perspective, this admittance representation restricts the network elements to be voltage controlled; resulting drawbacks treating zero-impedance branches in such applications as state estimation have long been recognized. This paper explores the advantages of the alternatives to Y bus -based formulations in power flow. It proposes a sparse tableau formulation (STF), and it demonstrates its computational efficiency, robustness, and generality in detailed comparison to traditional Y bus -based solution algorithms. In the examples of power networks ranging from 1888 to 82 000 buses, computational case studies indicate that STF provides comparable computational speed, while allowing simple treatment of zero-impedance branches and more reliably converging to solutions in many cases for which Y bus -based Newton algorithms diverge.

42 ENGINEERING↗
DOE OSTI2023

Modeling and Rapid Prototyping of Integrated Transmission-Distribution OPF Formulations with PowerModelsITD.jl

Conventional electric power systems are composed of different unidirectional power flow stages of generation, transmission, and distribution, managed independently by transmission system and distribution system operators. However, as distribution systems increase in complexity due to the integration of distributed energy resources, coordination between transmission and distribution networks will be imperative for the optimal operation of the power grid. However, coupling models and formulations between transmission and distribution is non-trivial, in particular due to the common practice of modeling transmission systems as single-phase, and distribution systems as multi-conductor phase-unbalanced. To enable the rapid prototyping of power flow formulations, in particular in the modeling of the boundary conditions between these two seemingly incompatible data models, we introduce PowerModelsITD.jl, a free, open-source toolkit written in Julia for integrated transmission-distribution (ITD) optimization that leverages mature optimization libraries from the InfrastructureModels.jl-ecosystem. The primary objective of the proposed framework is to provide baseline implementations of steady-state ITD optimization problems, while providing a common platform for the evaluation of emerging formulations and optimization problems. In this work, we introduce the nonlinear formulations currently supported in PowerModelsITD.jl, which include AC-polar, AC-rectangular, current-voltage, and a linear network transportation model. Results are validated using combinations of IEEE transmission and distribution networks.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2015

SU-E-T-515: Investigating the Linear Energy Transfer Dependency of Different PRESAGE Formulations in a Proton Beam

Purpose Previous studies have reported an under-response of PRESAGE in a proton beam as a Result of the extremely high LET in the distal end of the spread out Bragg peak (SOBP). This work is a preliminary investigation to quantify the effect of the formulation, specifically the concentration of halocarbon radical initiator relative to leuco dye, on radical recombination resulting in LET dependence. Methods The traditional PRESAGE formulation developed by Heuris Pharma was altered to constitute radical initiator concentrations of 5, 15, and 30% (low, medium, and high) by weight with all other components balanced to maintain proportionality. Chloroform was specifically examined in this study and all dosimeters were made in-house. Cylindrical PRESAGE dosimeters (3.5cm diameter and 6cm length) were made for each formulation and irradiated by a 200-MeV proton beam to 500 cGy across a 2cm SOBP. Dosimeters were read out using the DMOS optical-CT scanner. The dose distributions were analyzed and dose profiles were used to compare the relative dose response to find the stability across the high-LET region of the SOBP. LET dependence was measured by the variation to ion chamber measurements for the final 25% of the SOBP (∼0.5cm) prior to the distal-90 of each profile. Results Relative to ion chamber data, all PRESAGE dosimeters showed an under-response at the distal end of the SOBP. The medium concentration formulation matched most closely with an average 8.3% under-response closely followed by the low concentration at 12.2% and then the high concentration at 22.8%. In all three cases, the highest points of discrepancy were in the distal most regions. Conclusion The radical initiator concentration in PRESAGE can be tailored to reduce the LET dependence in a proton beam. This warrants further study to quantify comprehensively the effect of concentration of different halocarbon radical initiators on LET dependency. Grant number 5RO1CA100835.

60 APPLIED LIFE SCIENCES↗