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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Unravelling donor–acceptor film morphology formation for environmentally-friendly OPV ink formulations

The challenge of coating organic photovoltaics (OPV) from green solvents is to achieve the required nanostructured interpenetrating network of donor and acceptor domains based on a rational choice of solvent approach as opposed to the usual trial-and-error methods. We demonstrate here that we can achieve a bicontinuous interpenetrating network with nanoscale phase separation for the chosen donor–acceptor material system poly[2,3-bis-(3-octyloxyphenyl)quinoxaline-5,8-diyl-alt-thiophene-2,5-diyl]:phenyl-C 61 butyric acid methyl ester (TQ1:PC 61 BM) when processing from green solvent ink formulations. This structure is achieved by first calculating the Hansen solubility parameters (HSP) of the donor and acceptor materials, followed by careful choice of solvents with selective relative solubilities for the two materials based on the desired order of precipitation necessary for forming a nanostructured interdigitated network morphology. We found that the relative distances in Hansen space (R a ) between TQ1 and the primary solvent, on the one hand, and PC 61 BM and the primary solvent, on the other hand, could be correlated to the donor–acceptor morphology for the formulations based on the solvents d-limonene, anisole, and 2-methyl anisole, as well as the halogenated reference solvent o-dichlorobenzene. This nanostructured blend film morphology was characterised with scanning transmission X-ray microscopy (STXM) and transmission electron microscopy (TEM), and the film surface composition was analysed by near edge X-ray absorption fine structure (NEXAFS) spectroscopy. Hansen solubility theory, based on solution thermodynamics, has been used and we propose an HSP-based method that is a general platform for the rational design of ink formulations for solution-based organic electronics, in particular facilitating the green solvent transition of organic photovoltaics. Our results show that the bulk heterojunction morphology for a donor–acceptor system processed from customised solvent mixtures can be predicted by the HSP-based method with good reliability.

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

A full-dimensional quantum dynamical study of H{sub 2}+H{sub 2} collisions: Coupled-states versus close-coupling formulation

Collision-induced energy transfer involving H{sub 2} molecules plays an important role in many areas of physics. Kinetic models often require a complete set of state-to-state rate coefficients for H{sub 2}+H{sub 2} collisions in order to interpret results from spectroscopic observations or to make quantitative predictions. Recent progress in full-dimensional quantum dynamics using the numerically exact close-coupling (CC) formulation has provided good agreement with existing experimental data for low-lying states of H{sub 2} and increased the number of state-to-state cross sections that may be reliably determined over a broad range of energies. Nevertheless, there exist many possible initial states (e.g., states with high rotational excitation) that still remain elusive from a computational standpoint even at relatively low collision energies. In these cases, the coupled-states (CS) approximation offers an alternative full-dimensional formulation. We assess the accuracy of the CS approximation for H{sub 2}+H{sub 2} collisions by comparison with benchmark results obtained using the CC formulation. The results are used to provide insight into the orientation effects of the various internal energy transfer mechanisms. A statistical CS approximation is also investigated and cross sections are reported for transitions which would otherwise be impractical to compute.

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

Novel method to control explosive shock sensitivity: Formulations, experiments and modeling

The microstructure of a heterogeneous high explosive (HE) affects its shock initiation sensitivity and detonation performance. Alteration of the void content and/or void structure (i.e., bulk heating or mechanical damage) therefore changes the shock initiation behavior. Controlling and predicting the change in shock sensitivity after an HE has undergone microstructural changes addresses an important and challenging goal for the design and understanding of novel energetic material formulations. In this study, we aimed to develop HE systems to precisely tune shock sensitivity by a thermal treatment prior to use. Specifically, we incorporated a small fraction (1 wt% or less) of thermally expandable microspheres (TEMs) during the formulation process of various plastic-bonded explosives (PBX). TEMs typically consist of a thermoplastic acrylonitrile shell (10-50 µm diameter), which encapsulates 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% (irreversibly). Here, we present our progress on shock sensitivity comparisons of HE formulations doped with TEMs after heating and expansion. Experiments, using Composition C-4 thermally cycled to 120°C, showed an increase in shock sensitivity. Additionally, mesoscale modeling revealed that the TEM itself does not act as a hotspot, but instead has a secondary effect on the response of nearby voids. We conclude that TEMs indeed do provide a method to tune shock sensitivity. Further, our results emphasize the non-uniformity of shock waves at the mesoscale, and that upstream defects influence the reactivity of downstream defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Hamiltonian formulations for perturbed dissipationless plasma equations

The Hamiltonian formulations for the perturbed Vlasov–Maxwell equations and the perturbed ideal magnetohydrodynamics (MHD) equations are expressed in terms of the perturbation derivative ∂ F / ∂ ϵ ≡ [ F , S ] of an arbitrary functional F [ ψ ] of the Vlasov–Maxwell fields ψ = ( f , E , B ) or the ideal MHD fields ψ = ( ρ , u , s , B ) , which are assumed to depend continuously on the (dimensionless) perturbation parameter ϵ . In this study, [ , ] denotes the functional Poisson bracket for each set of plasma equations and the perturbation action functional S is said to generate dynamically accessible perturbations of the plasma fields. The new Hamiltonian perturbation formulation introduces a framework for functional perturbation methods in plasma physics and highlights the crucial roles played by polarization and magnetization in Vlasov–Maxwell and ideal MHD perturbation theories. One application considered in this paper is a formulation of plasma stability that guarantees dynamical accessibility and leads to a natural generalization to higher-order perturbation theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Alternative formulation of weak magnetohydrodynamic turbulence theory

In a recent paper, the weak turbulence theory for incompressible magnetohydrodynamics is formulated by employing the method customarily applied in the context of kinetic weak plasma turbulence theory. Such an approach simplified certain mathematical procedures including achieving the closure relationship. The formulation in the above-cited paper starts from the equations of incompressible magnetohydrodynamic (MHD) theory expressed via Elsasser variables. The derivation of nonlinear wave kinetic equation therein is obtained via a truncated solution at the second-order of iteration following the standard practice. In the present paper, the weak MHD turbulence theory is alternatively formulated by employing the pristine form of incompressible MHD equation rather than that expressed in terms of Elsasser fields. The perturbative expansion of the nonlinear momentum equation is carried out up to the third-order iteration rather than imposing the truncation at the second order. Furthermore, it is found that while the resulting wave kinetic equation is identical to that obtained in the previous paper cited above, the third-order nonlinear correction plays an essential role for properly calculating derived quantities such as the total and residual energies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

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↗
DOE OSTI2025

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↗
DOE OSTI2010

Nonadiabatic couplings from the Kohn-Sham derivative matrix: Formulation by time-dependent density-functional theory and evaluation in the pseudopotential framework

We study the time-dependent density-functional theory formulation of nonadiabatic couplings (NAC's) to settle problems regarding practical calculations. NAC's have so far been rigorously formulated on the basis of the density response scheme and expressed using the nuclear derivative of the Hamiltonian, {partial_derivative}H/{partial_derivative}R, whereby causing the pseudopotential problem. When rewritten using the nuclear derivative operator, {partial_derivative}/{partial_derivative}R, or the d operator, the formula is found free of the problem and thus provides a working numerical scheme. The d-operator-based formulation also allows us to lay a foundation on the empirical Slater transition-state method and to show an improved way of using the auxiliary excited-state wave-function ansatz, both of which have been utilized in previous works. Evaluation of NAC near either the Jahn-Teller or the Renner-Teller intersection in various molecular systems shows that the values of NAC are much improved over previous calculations when the d-operator formula is implemented in the pseudopotential framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Majorization formulation of uncertainty in quantum mechanics

Heisenberg's uncertainty principle is formulated for a set of generalized measurements within the framework of majorization theory, resulting in a partial uncertainty order on probability vectors that is stronger than those based on quasientropic measures. The theorem that emerges from this formulation guarantees that the uncertainty of the results of a set of generalized measurements without a common eigenstate has an inviolable lower bound which depends on the measurement set but not the state. A corollary to this theorem yields a parallel formulation of the uncertainty principle for generalized measurements corresponding to the entire class of quasientropic measures. Optimal majorization bounds for two and three mutually unbiased bases in two dimensions are calculated. Similarly, the leading term of the majorization bound for position and momentum measurements is calculated which provides a strong statement of Heisenberg's uncertainty principle in direct operational terms. Another theorem provides a majorization condition for the least-uncertain generalized measurement of a given state with interesting physical implications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Theoretical formulation of phase space microbunching instability in the presence of intrabeam scattering for single-pass or recirculation accelerators

Microbunching instability (MBI) has been one of the most challenging issues in designs of high-brightness beam transport lines for single-pass or recirculating accelerators. Although the intrabeam scattering (IBS) has long been studied in lepton or hadron storage rings as a slow diffusion process or in high-intensity proton linear accelerators as one mechanism for the beam halo, the effects of IBS on single-pass or recirculating electron accelerators have drawn attention only in the recent two decades due to emergence of linac-based or energy-recovery-linac-based fourth-generation light sources, which require high-quality electron beams during the beam transport. In this paper we develop a theoretical formulation of microbunching instability in the presence of IBS for single-pass or recirculation accelerators. To quantify MBI with inclusion of IBS, we start from the Vlasov-Fokker-Planck (VFP) equation, combining both collective interactions and incoherent IBS effects. The linearized VFP equation and the corresponding friction and diffusion coefficients are derived. The evolutions of the resultant density and energy modulations are formulated as a set of coupled integral equations. The theoretical formulation is then applied to a recirculating beamline design. The results from the semianalytical calculation are compared and show good agreement with massive particle tracking simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Euclidean formulation of relativistic quantum mechanics of N particles

A Euclidean formulation of relativistic quantum mechanics for systems of a finite number of degrees of freedom is discussed. Relativistic treatments of quantum theory are needed to study hadronic systems at subhadronic distance scales. While direct interaction approaches to relativistic quantum mechanics have proved to be useful, they have two disadvantages. One is that cluster properties are difficult to realize for systems of more than two particles. The second is that the relation to quantum field theories is indirect. Euclidean formulations of relativistic quantum mechanics provide an alternative representation that does not have these difficulties. More surprising, the theory can be formulated entirely in the Euclidean representation without the need for analytic continuation. In this work a Euclidean representation of a relativistic N-particle system is discussed. Kernels for systems of N free particles of any spin are given and shown to be reflection positive. Explicit formulas for generators of the Poincaré group for any spin are constructed and shown to be self-adjoint on the Euclidean representation of the Hilbert space. The structure of correlations that preserve both the Euclidean covariance and reflection positivity is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

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↗
DOE OSTI2022

A General State Estimation Formulation for Three-Phase Unbalanced Power Systems

Almost all of the three-phase state estimation algorithms assume existence of a reference bus whose phase angles are perfectly balanced. This assumption is quite realistic for transmission systems, and also for most distribution systems that are connected to a strong transmission system where transmission side can modeled by a balanced reference bus. However, for distribution systems having high penetration of renewable sources or for microgrids operating in islanded mode, the assumption of a balanced reference bus will not be realistic. While there are recent publications focusing on this problem, formulation of the three-phase unbalanced state estimation problem with proper treatment of the reference bus remains unaddressed. In this paper, a new formulation will be described where an accurate state estimation solution can be obtained for any unbalanced threephase system irrespective of its operating conditions (balanced or highly unbalanced), configuration (isolated microgrid, connected to transmission system, etc.) and whether or not it contains any synchronous generators. Validation of the proposed formulation will be carried out via simulations.

State Estimation, Distribution System, Unbalanced ↗
DOE OSTI2016

An Improved Stochastic Unit Commitment Formulation to Accommodate Wind Uncertainty

The United States targets to supply 20% of its electricity generation using wind energy by 2030. The expansion of renewable resources, especially weather-based resources such as wind, creates more uncertainty and variability in the operation of the power grid. New methods and approaches in electricity market operations are needed to efficiently manage the continuing increase in variability and uncertainty caused by expanding intermittent wind. This paper proposes an improved stochastic programming approach for incorporating wind uncertainty into energy markets. The proposed formulation improves the two-stage stochastic unit commitment problem by introducing a dynamic decision making approach similar to a multi-stage formulation in the presence of wind power scenarios which are not well represented by a scenario tree. The numerical results present up to 1%-2% decrease in operational costs compared to the two-stage stochastic unit commitment formulation.

Electricity markets↗
DOE OSTI2012

Comparison of signal to noise ratios from spatial and frequency domain formulations of nonprewhitening model observers in digital mammography

Purpose: Image quality indices based upon model observers are promising alternatives to laborious human readings of contrast-detail images. This is especially appealing in digital mammography as limiting values for contrast thresholds determine, according to some international protocols, the acceptability of these systems in the radiological practice. The objective of the present study was to compare the signal to noise ratios (SNR) obtained with two nonprewhitening matched filter model observer approaches, one in the spatial domain and the other in the frequency domain, and with both of them worked out for disks as present in the CDMAM phantom. Methods: The analysis was performed using images acquired with the Siemens Novation and Inspiration digital mammography systems. The spatial domain formulation uses a series of high dose CDMAM images as the signal and a routine exposure of two flood images to calculate the covariance matrix. The frequency domain approach uses the mathematical description of a disk and modulation transfer function (MTF) and noise power spectrum (NPS) calculated from images. Results: For both systems most of the SNR values calculated in the frequency domain were in very good agreement with the SNR values calculated in the spatial domain. Both the formulations in the frequency domain and in the spatial domain show a linear relationship between SNR and the diameter of the CDMAM discs. Conclusions: The results suggest that both formulations of the model observer lead to very similar figures of merit. This is a step forward in the adoption of figures of merit based on NPS and MTF for the acceptance testing of mammography systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Modeling Single and Two-Phase Transport in Thin Porous Layers Using a Composite Continuum-Pore Network Formulation

In this work, a composite continuum-pore network formulation is presented to model single and two-phase transport thin porous layers, such as gas diffusion layers in polymer electrolyte fuel cells (PEFCs) and active electrodes in redox flow batteries (RFBs). The formulation can be integrated into CFD codes, thus combining the ease of implementation of continuum-based modeling and the computational power of pore-network modeling. The composite model includes a control volume (CV) mesh at the layer scale, which embeds a cubic pore network [1,2]. The pore-network model is used to determine analytically local anisotropic effective transport properties (local effective diffusivity and permeability), which are mapped onto the CV mesh to simulate transport in the porous transport layer [3,4]. Good agreement is found between the predicted global effective transport properties (global effective diffusivity and permeability) under dry and wet conditions and previous experimental data reported in the literature for Toray TGP-H series carbon paper. Water saturation distributions are also compared with results obtained using X-ray computed tomography [4]. [1] P.A. García-Salaberri, I.V. Zenyuk, J.T. Gostick, A.Z. Weber, Modeling Gas Diffusion Layer in Polymer Electrolyte Fuel Cells Using a Continuum-Based Pore-Network Formulation, ECS Trans. 97 (2020) 615. [2] P.A. García-Salaberri, Modeling diffusion and convection in thin porous transport layers using a composite continuum-network model: Application to gas diffusion layers in polymer electrolyte fuel cells, Int J. Heat Mass Transf. (2020), submitted. [3] P.A. García-Salaberri, J.T. Gostick, G. Hwang, A.Z. Weber, M. Vera, Effective diffusivity in partially-saturated carbon-fiber gas diffusion layers: Effect of local saturation and application to macroscopic continuum models, J. Power Sources 296 (2015) 440–453. [4] P.A. García-Salaberri, G. Hwang, M. Vera, A.Z. Weber, J.T. Gostick, Effective diffusivity in partially-saturated carbon-fiber gas diffusion layers: Effect of through-plane saturation distribution, Int. J. Heat Mass Transf. 86 (2015) 319–333.

Garcia-Salaberri, Pablo↗
DOE OSTI2022

A deformation-dependent coupled Lagrangian/semi-Lagrangian meshfree hydromechanical formulation for landslide modeling

Abstract The numerical modelling of natural disasters such as landslides presents several challenges for conventional mesh-based methods such as the finite element method (FEM) due to the presence of numerically challenging phenomena such as severe material deformation and fragmentation. In contrast, meshfree methods such as the reproducing kernel particle method (RKPM) possess unique features conducive to modelling extreme events such as the absence of a structured mesh and the ease of adaptive refinement, among others. While the semi-Lagrangian reproducing kernel (SL-RK) shape functions of RKPM defined in the current configuration have proven to be effective in extreme event modelling, the computational cost for the re-evaluation of the shape functions at every time step is costly. In this work, a deformation-dependent coupling of the Lagrangian reproducing kernel (L-RK) and SL-RK approximations is proposed for the solution of a hydro-mechanical formulation for effective simulations of landslides. The ramp function is constructed based on an equivalent plastic strain as a deformation-dependent transition from L-RK shape functions to SL-RK ones as the deformation progresses. The particular focus of the paper will be on modelling seepage-induced landslides with a mixed $$u$$ u – $$p$$ p formulation to couple the solid and fluid phases. Examples are presented to examine the effectiveness of this coupled Lagrangian/semi-Lagrangian reproducing kernel (L–SL RK) formulation and to highlight its performance in landslide modelling.

Baek, Jonghyuk↗
DOE OSTI2024

Decomposable Formulation of Transmission Constraints for Decentralized Power Systems Optimization

One of the most complicating factors in decentralized solution methods for a broad range of power system optimization problems is the modeling of power flow equations. Existing formulations for direct current power flows either have limited scalability or are very dense and unstructured, making them unsuitable for large-scale decentralized studies. Here, in this work, we present a novel sparsified variant of the injection shift factors formulation, which has a decomposable block-diagonal structure and scales well for large systems. We also propose a decentralized solution method, based on the alternating direction multiplier method, that efficiently handles transmission line outages in N-1 security requirements. Benchmarks on multizonal security-constrained unit commitment problems show that the proposed formulation and algorithm can reliably and efficiently solve interconnection-level test systems with up to 6,515 buses with no convergence or numerical issues.

Alternating-method of multipliers (ADMM)↗