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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2024

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

Anisotropic finite hyper-elastoplasticity of geomaterials with Drucker-Prager/Cap type constitutive model formulation

The formulation of large strain anisotropic hyper-elastoplasticity of geomaterials is examined. Attention is given to the role of structure tensors (also called fabric tensors), especially in context of the Eshelby-Mandel stress and large inelastic volume changes attributable to porosity. Both (hyper-)elastic and inelastic orthotropic symmetry, reducing to the particular case of transverse isotropy, are considered. Specific material assumptions and constitutive choices are identified for the development of a novel Anisotropic Drucker-Prager/Cap (ADPC) model formulated within the intermediate configuration consistent with multiplicative split of the deformation gradient. The model is calibrated to existing experimental measurements, including high pressure large strain triaxial compression of lithographic (Solnhofen) limestone and triaxial compression measurements on Tournemire shale assessing elastoplastic anisotropy. Manifest implications of constitutive theory are investigated, including consequences of recognizing (or not) the Eshelby-Mandel stress as energy conjugate to the plastic velocity gradient and including (or not) contribution from the skew-symmetric parts of the Mandel stress to the plastic anisotropy. Numerical simple shear experiments and large deformation simulated indentation experiments are provided in order to investigate model predictions and demonstrate the overall robustness in finite element modeling.

42 ENGINEERING↗
DOE OSTI2012

Glioblastoma Treatment: Bypassing the Toxicity of Platinum Compounds by Using Liposomal Formulation and Increasing Treatment Efficiency With Concomitant Radiotherapy

Purpose: Treatments of glioblastoma with cisplatin or oxaliplatin only marginally improve the overall survival of patients and cause important side effects. To prevent adverse effects, improve delivery, and optimize the tumor response to treatment in combination with radiotherapy, a potential approach consists of incorporating the platinum agent in a liposome. Methods and Materials: In this study, cisplatin, oxaliplatin, carboplatin, Lipoplatin (the liposomal formulation of cisplatin), and Lipoxal (the liposomal formulation of oxaliplatin) were tested on F98 glioma orthotopically implanted in Fischer rats. The platinum compounds were administered by intracarotid infusion and were assessed for the ability to reduce toxicity, improve cancer cell uptake, and increase survival of animals when combined or not combined with radiotherapy. Results: The tumor uptake was 2.4-fold more important for Lipoxal than the liposome-free oxaliplatin. Lipoxal also improved the specificity of oxaliplatin as shown by a higher ratio of tumor to right hemisphere uptake. Surprisingly, Lipoplatin led to lower tumor uptake compared with cisplatin. However, Lipoplatin had the advantage of largely reducing the toxicity of cisplatin and allowed us to capitalize on the anticancer activity of this agent. Conclusion: Among the five platinum compounds tested, carboplatin showed the best increase in survival when combined with radiation for treatment of glioma implanted in Fischer rats.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Augmented Lagrangian formulation of orbital-free density functional theory

We present an Augmented Lagrangian formulation and its real-space implementation for non-periodic Orbital-Free Density Functional Theory (OF-DFT) calculations. In particular, we rewrite the constrained minimization problem of OF-DFT as a sequence of minimization problems without any constraint, thereby making it amenable to powerful unconstrained optimization algorithms. Further, we develop a parallel implementation of this approach for the Thomas–Fermi–von Weizsacker (TFW) kinetic energy functional in the framework of higher-order finite-differences and the conjugate gradient method. With this implementation, we establish that the Augmented Lagrangian approach is highly competitive compared to the penalty and Lagrange multiplier methods. Additionally, we show that higher-order finite-differences represent a computationally efficient discretization for performing OF-DFT simulations. Overall, we demonstrate that the proposed formulation and implementation are both efficient and robust by studying selected examples, including systems consisting of thousands of atoms. We validate the accuracy of the computed energies and forces by comparing them with those obtained by existing plane-wave methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Physical formulation and numerical algorithm for simulating N immiscible incompressible fluids involving general order parameters

We present a family of physical formulations, and a numerical algorithm, based on a class of general order parameters for simulating the motion of a mixture of N (N⩾2) immiscible incompressible fluids with given densities, dynamic viscosities, and pairwise surface tensions. The N-phase formulations stem from a phase field model we developed in a recent work based on the conservations of mass/momentum, and the second law of thermodynamics. The introduction of general order parameters leads to an extremely strongly-coupled system of (N−1) phase field equations. On the other hand, the general form enables one to compute the N-phase mixing energy density coefficients in an explicit fashion in terms of the pairwise surface tensions. We show that the increased complexity in the form of the phase field equations associated with general order parameters in actuality does not cause essential computational difficulties. Our numerical algorithm reformulates the (N−1) strongly-coupled phase field equations for general order parameters into 2(N−1) Helmholtz-type equations that are completely de-coupled from one another. This leads to a computational complexity comparable to that for the simplified phase field equations associated with certain special choice of the order parameters. We demonstrate the capabilities of the method developed herein using several test problems involving multiple fluid phases and large contrasts in densities and viscosities among the multitude of fluids. In particular, by comparing simulation results with the Langmuir–de Gennes theory of floating liquid lenses we show that the method using general order parameters produces physically accurate results for multiple fluid phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A signed particle formulation of non-relativistic quantum mechanics

A formulation of non-relativistic quantum mechanics in terms of Newtonian particles is presented in the shape of a set of three postulates. In this new theory, quantum systems are described by ensembles of signed particles which behave as field-less classical objects which carry a negative or positive sign and interact with an external potential by means of creation and annihilation events only. This approach is shown to be a generalization of the signed particle Wigner Monte Carlo method which reconstructs the time-dependent Wigner quasi-distribution function of a system and, therefore, the corresponding Schrödinger time-dependent wave-function. Its classical limit is discussed and a physical interpretation, based on experimental evidences coming from quantum tomography, is suggested. Moreover, in order to show the advantages brought by this novel formulation, a straightforward extension to relativistic effects is discussed. To conclude, quantum tunnelling numerical experiments are performed to show the validity of the suggested approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Gaussian-windowed frame based method of moments formulation of surface-integral-equation for extended apertures

Scattering or coupling of electromagnetic beam-field at a surface discontinuity separating two homogeneous or inhomogeneous media with different propagation characteristics is formulated using surface integral equation, which are solved by the Method of Moments with the aid of the Gabor-based Gaussian window frame set of basis and testing functions. The application of the Gaussian window frame provides (i) a mathematically exact and robust tool for spatial-spectral phase-space formulation and analysis of the problem; (ii) a system of linear equations in a transmission-line like form relating mode-like wave objects of one medium with mode-like wave objects of the second medium; (iii) furthermore, an appropriate setting of the frame parameters yields mode-like wave objects that blend plane wave properties (as if solving in the spectral domain) with Green's function properties (as if solving in the spatial domain); and (iv) a representation of the scattered field with Gaussian-beam propagators that may be used in many large (in terms of wavelengths) systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

High order ADER schemes for a unified first order hyperbolic formulation of continuum mechanics: Viscous heat-conducting fluids and elastic solids

Highlights: • High order schemes for a unified first order hyperbolic formulation of continuum mechanics. • The mathematical model applies simultaneously to fluid mechanics and solid mechanics. • Viscous fluids are treated in the frame of hyper-elasticity as generalized visco-plastic solids. • Formal asymptotic analysis reveals the connection with the Navier–Stokes equations. • The distortion tensor A in the model appears to be well-suited for flow visualization. - Abstract: This paper is concerned with the numerical solution of the unified first order hyperbolic formulation of continuum mechanics recently proposed by Peshkov and Romenski [110], further denoted as HPR model. In that framework, the viscous stresses are computed from the so-called distortion tensor A, which is one of the primary state variables in the proposed first order system. A very important key feature of the HPR model is its ability to describe at the same time the behavior of inviscid and viscous compressible Newtonian and non-Newtonian fluids with heat conduction, as well as the behavior of elastic and visco-plastic solids. Actually, the model treats viscous and inviscid fluids as generalized visco-plastic solids. This is achieved via a stiff source term that accounts for strain relaxation in the evolution equations of A. Also heat conduction is included via a first order hyperbolic system for the thermal impulse, from which the heat flux is computed. The governing PDE system is hyperbolic and fully consistent with the first and the second principle of thermodynamics. It is also fundamentally different from first order Maxwell–Cattaneo-type relaxation models based on extended irreversible thermodynamics. The HPR model represents therefore a novel and unified description of continuum mechanics, which applies at the same time to fluid mechanics and solid mechanics. In this paper, the direct connection between the HPR model and the classical hyperbolic–parabolic Navier–Stokes–Fourier theory is established for the first time via a formal asymptotic analysis in the stiff relaxation limit. From a numerical point of view, the governing partial differential equations are very challenging, since they form a large nonlinear hyperbolic PDE system that includes stiff source terms and non-conservative products. We apply the successful family of one-step ADER–WENO finite volume (FV) and ADER discontinuous Galerkin (DG) finite element schemes to the HPR model in the stiff relaxation limit, and compare the numerical results with exact or numerical reference solutions obtained for the Euler and Navier–Stokes equations. Numerical convergence results are also provided. To show the universality of the HPR model, the paper is rounded-off with an application to wave propagation in elastic solids, for which one only needs to switch off the strain relaxation source term in the governing PDE system. We provide various examples showing that for the purpose of flow visualization, the distortion tensor A seems to be particularly useful.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

A generalized transport-velocity formulation for smoothed particle hydrodynamics

The standard smoothed particle hydrodynamics (SPH) method suffers from tensile instability. In fluid-dynamics simulations this instability leads to particle clumping and void regions when negative pressure occurs. In solid-dynamics simulations, it results in unphysical structure fragmentation. In this work the transport-velocity formulation of Adami et al. (2013) is generalized for providing a solution of this long-standing problem. Other than imposing a global background pressure, a variable background pressure is used to modify the particle transport velocity and eliminate the tensile instability completely. Furthermore, such a modification is localized by defining a shortened smoothing length. The generalized formulation is suitable for fluid and solid materials with and without free surfaces. The results of extensive numerical tests on both fluid and solid dynamics problems indicate that the new method provides a unified approach for multi-physics SPH simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

High-order gradients with the shifted boundary method: An embedded enriched mixed formulation for elliptic PDEs

Here, we propose an extension of the embedded boundary method known as “shifted boundary method” to elliptic diffusion equations in mixed form (e.g., Darcy flow, heat diffusion problems with rough coefficients, etc.). Our aim is to obtain an improved formulation that, for linear finite elements, is at least second-order accurate for both flux and primary variable, when either Dirichlet or Neumann boundary conditions are applied. Following previous work of Nishikawa and Mazaheri in the context of residual distribution methods, we consider the mixed form of the diffusion equation (i.e., with Darcy-type operators), and introduce an enrichment of the primary variable. This enrichment is obtained exploiting the relation between the primary variable and the flux variable, which is explicitly available at nodes in the mixed formulation. The proposed enrichment mimics a formally quadratic pressure approximation, although only nodal unknowns are stored, similar to a linear finite element approximation. We consider both continuous and discontinuous finite element approximations and present two approaches: a non-symmetric enrichment, which, as in the original references, only improves the consistency of the overall method; and a symmetric enrichment, which enables a full error analysis in the classical finite element context. Combined with the shifted boundary method, these two approaches are extended to high-order embedded computations, and enable the approximation of both primary and flux (gradient) variables with second-order accuracy, independently on the type of boundary conditions applied. We also show that the primary variable is third-order accurate, when pure Dirichlet boundary conditions are embedded.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Formulations and analysis of Compton scattering for deterministic thermal radiation transport

In high-energy-density physics at low densities and high temperatures, Compton scattering can be the dominant interaction mechanism between thermal radiation and the material. Boson statistics leads to induced scattering, where the rate of scattering is proportional to the square of the photon intensity. To explore the importance of this nonlinear effect in the context of thermal radiative transfer, in this work we formulate several linear approximations to the scattering rate and compare them to the proper nonlinear rate in an infinite-medium setting. With the goal of implementing Compton scattering in a radiation transport code, we study heating-rate accuracy and iterative convergence properties for the linear methods as we vary the fidelity of the discretization in photon frequency. We find the equilibrium out-scattering formulation to be the most accurate linear method, with heating rate errors around 15% when the radiation is heating the material and up to 50% when the radiation is cooling the material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

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

From coherent to incoherent mismatched interfaces. A generalized continuum formulation of surface stresses

The equilibrium of coherent and incoherent mismatched interfaces is reformulated in the context of continuum mechanics based on the Gibbs dividing surface concept. Two surface stresses are introduced: a coherent surface stress and an incoherent surface stress, as well as a transverse excess strain. Additionally, the coherent surface stress and the transverse excess strain represent the thermodynamic driving forces of stretching the interface while the incoherent surface stress represents the driving force of stretching one crystal while holding the other fixed and thereby altering the structure of the interface. These three quantities fully characterize the elastic behavior of coherent and incoherent interfaces as a function of the in-plane strain, the transverse stress and the mismatch strain. The isotropic case is developed in detail and particular attention is paid to the case of interfacial thermo-elasticity. This exercise provides an insight on the physical significance of the interfacial elastic constants introduced in the formulation and illustrates the obvious coupling between the interface structure and its associated thermodynamics quantities. Finally, an example based on atomistic simulations of Cu/Cu 2 O interfaces is given to demonstrate the relevance of the generalized interfacial formulation and to emphasize the dependence of the interfacial thermodynamic quantities on the incoherency strain with an actual material system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Non-local large-strain FFT-based formulation and its application to interface-dominated plasticity of nano-metallic laminates

Here, this paper presents a novel formulation and its robust numerical implementation of strain-gradient (SG) crystal plasticity within a large-strain (LS) elasto-viscoplastic (EVP) fast Fourier transform (FFT)-based micromechanical model. The resulting non-local SG-LS-EVPFFT formulation is used to model and understand the process of kink band formation during layer-parallel compression of nano-metallic laminates (NMLs). NMLs are layered composites with nanoscale thicknesses, thus requiring consideration of the interaction between dislocations and interfaces within the micromechanical model. The length-scale parameter of the SG model is calibrated by simulating a double pile-up and comparing predictions to analytical solution. This required new expressions for the defect energy, resulting in more accurate double pile-up predictions. The calibrated SG-LS-EVPFFT model is then used to simulate layer-parallel compression of copper-niobium NML. Formation of kink bands is predicted, and the model is used to rationalize the microscopic mechanisms enabling the formation process. It is found that accumulation of dislocations at interfaces leads to activation of layer-parallel slip, which in turn leads to kink band formation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Formulation and calibration of two-dimensional constitutive models for composite structures based on panel tests

Steel Plate Concrete (SC) composite members have been widely adopted because of its cost-efficiency and enhanced structural behavior. While researchers have attempted to study its in-plane shear behavior in the past twenty years, very limited number of large-scale pure shear tests were performed due to the challenge of experimental set-up and the availability of facilities. In this paper, a series of uniaxial loading tests and two full-scale pure shear panel tests of SC members were reported, on which the “mechanics-based Membrane Model of SC elements (MM-SC)” is developed. The MM-SC model is based on the fixed-angle crack formulation and the smeared-crack formulation, in which the experimental-based uniaxial constitutive models are implemented, considering the local buckling of faceplate, the tension stiffening of steel plate, the strength degradation of cracked concrete and the confinement effect of concrete. The proposed MM-SC model is subsequently incorporated into the object-oriented software OpenSEES. Finally, the simulation results of proposed model well predict the SC test observations in terms of critical branch points and structural behaviors, including initial stiffness, cracking strength, post-crack stiffness, yield strength, maximum strength, and failure modes.

42 ENGINEERING↗
DOE OSTI2026

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 ↗