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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Particle Interaction Physics Model Formulation for Plume-Surface Interaction Erosion and Cratering

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The Predictive Simulation Capability (PSC) Element is focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the predictive simulation capability, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper will present the status of implementing and maturing particle-particle interaction constituent physics models essential in simulating the landing surface granular material flow under PSI effects. This gas-particle multi-phase interaction modeling of plume impingement flow on the extra-terrestrial soil material is performed with the Gas-Granular Flow Solver (GGFS) addressed in a companion paper. The response of regolith particle flow induced by lander PSI requires accurate representation of the regolith granular material fluidic behavior and gas-granular interactions. The lunar regolith, as the extreme example, is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking particle shapes for the very jagged particles. The combination of particle shape and size distribution has been identified as major drivers in the complex particle flow response and resulting crater shape characteristics of extraterrestrial granular material. Constituent models for spherical particles can be formulated directly from particle kinetics theory. Complex particle shapes can be modeled by gluing together elemental spherical shapes into composite particles, requiring a Discrete Element Model (DEM) particle kinetics modeling approach to extract data and formulate constituent models. Mixture constituent models for poly-disperse mixtures (i.e, containing distribution of particle sizes) have recently been developed. The required non-spherical particle mixture granular material response closure models are then obtained through small-scale unit physics DEM simulations for the range of particle shapes, mixtures and packing densities. The granular material response closure models are then implemented in the Eulerian granular flow formulation. This DEM-based constituent model extraction process and formulation of poly-disperse particle mixtures has been successfully developed by small business and academic partners in the development of the Gas-Granular Flow Solver (GGFS) simulation program simulation framework. The currently implemented capabilities have reached the capability level of modeling bi-disperse, non-spherical particle mixtures is being continuously extended towards computational modeling of full range irregular particle mixtures. Under the GCD project, this technology is being further developed, transferred to NASA analysts, and matured towards application readiness. The predictive simulation capability team under the GCD project has acquired the modeling tools and processes of the DEM based constituent model formulation from the GGFS development team and is developing the capability to replicate the existing process. This is the first important step towards the ability of the NASA team to independently perform such model development in a production setting. Further efforts are underway to migrate the DEM based model simulation process performed with the academic based tools to more capable Open Source, highly parallelized simulation tools for efficient operation on NASA HPC assets. Evaluation of the currently implemented (such as mono-disperse and bi-disperse spherical and irregular shape particle constituent model applications) and continuously evolving full-range particle physics models in the GGFS tool is performed by the NASA team to advance application readiness of the simulations. Application testing for complex PSI erosions and cratering scenarios such as the Apollo LM is performed for axi-symmetric and full 3D simulations to aid the tool developers in achieving practical application readiness for NASA projects. Important validation and application testing will further be performed against experimental data generated under the GCD PSI project experimental component.

Peter A Liever↗
DOE OSTI2025

Time-domain all-frequency stable formulation for low-frequency electromagnetic simulation with Newmark-β time integration

An implicitly Coulomb-gauged A-ϕ formulation has previously been proposed and validated for finite ele- ment simulations of low-frequency and multiscale electromag- netic problems in the frequency domain. This formulation has demonstrated numerical stability across all frequencies, with its accuracy, efficiency, and iterative convergence established in various frequency-domain scenarios. However, direct time- domain computation is often preferable for wideband electro- magnetic problems and is typically indispensable in nonlinear and multiphysics simulations. In this work, the A-ϕ formulation is extended to the time domain. By incorporating the well-known Newmark-β time integration scheme, the proposed formulation is validated through capacitive and inductive test cases. The results confirm the solution’s accuracy and demonstrate the formulation’s stability in the time domain.

Mekonnen, Minyechil↗
NASA NTRS2006

Nonlinear consider covariance analysis using a sigma-point filter formulation

The research reported here extends the mathematical formulation of nonlinear, sigma-point estimators to enable consider covariance analysis for dynamical systems. This paper presents a novel sigma-point consider filter algorithm, for consider-parameterized nonlinear estimation, following the unscented Kalman filter (UKF) variation on the sigma-point filter formulation, which requires no partial derivatives of dynamics models or measurement models with respect to the parameter list. It is shown that, consistent with the attributes of sigma-point estimators, a consider-parameterized sigma-point estimator can be developed entirely without requiring the derivation of any partial-derivative matrices related to the dynamical system, the measurements, or the considered parameters, which appears to be an advantage over the formulation of a linear-theory sequential consider estimator. It is also demonstrated that a consider covariance analysis performed with this 'partial-derivative-free' formulation yields equivalent results to the linear-theory consider filter, for purely linear problems.

covariance↗
DOE OSTI2024

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗
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 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 OSTI2024

Fragmentation in Collisions of Snow with Graupel/Hail: New Formulation from Field Observations

Abstract Secondary ice production (SIP) has been attributed to the generation of most ice particles observed in precipitating clouds with cloud tops warmer than −36°C, from various aircraft- and ground-based field observations across the globe. One of the known SIP mechanisms is fragmentation during collisions among ice particles. It has been studied with our theoretical formulation, which has been applied in microphysical schemes of atmospheric models in a few studies. These have predicted an extensive impact on cloud glaciation and radiative properties. However, there has been a lack of experimental field studies, especially involving naturally falling snowflakes, to better understand this particular mechanism of SIP. This study reports the first field measurements with modern technology for fragmentation during collisions between naturally falling snowflakes and graupel/hail particles. This was observed with an innovatively designed portable chamber that was deployed outdoors in northern Sweden. Applying the observations from this field-based study, we optimized the existing formulation for predicting numbers of fragments from collisions of snow with graupel/hail. The observations show the average numbers of fragments per collision for dendritic (3–12 mm) and nondendritic (1–3 mm) snow were about 12 and 1, respectively. This represents a boost of predicted fragment numbers relative to our original formulation published in 2017. The updated formulation for breakup in ice–ice collisions can be implemented in the microphysical schemes of atmospheric models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1973

Boundary lubrication of formulated C-ethers in air to 300 C. 2: Organic acid additives

Friction and wear measurements were made on CVM M-50 steel lubricated with three C-ether (modified polyphenyl ether) formulations in dry and moist air. Results were compared to those obtained with a formulated Type 2 ester and the C-ether base fluid. A ball-on-disk sliding friction apparatus was used. Experimental conditions were a 1-kilogram load, a 17-meter/minute surface speed, and a 25 to 300 C (77 to 572 F) disk temperature range. The three C-ether formulations yielded better boundary lubricating characteristics than the Type 2 ester under most test conditions. All C-ether formulations exhibited higher friction coefficients than the ester from 150 to 300 C (302 to 572 F) and similar or lower values from 25 to 150 C (77 to 302 F).

Jones, W. R., Jr.↗
NASA NTRS1974

Finite element formulations for problems of large elastic-plastic deformation

An Eulerian finite element formulation is presented for problems of large elastic-plastic flow. The method is based on Hill's variational principle for incremental deformations, and is suited to isotropically hardening Prandtl-Reuss materials. The formulation is given in a manner which allows any conventional finite element program, for "small strain" elasticplastic analysis, to be simply and rigorously adapted to problems involving arbitrary amounts of deformation and arbitrary levels of stress in comparison to plastic deformation moduli. The method is applied to a necking bifurcation analysis of a bar in plane-strain tension. A unified general formulation of finite element equations, both Lagrangian and Eulerian, for large deformations, with arbitrary choice of the conjugate stress and strain measures, and a discussion is given of other proposed formulations for elastic-plastic finite element analysis at large strain.

Mcmeeking, R. M.↗
NASA NTRS1975

Boundary lubrication of formulated C-ethers in air to 300 C

Friction and wear measurements were made in dry and moist air on CVM M-50 steel lubricated with six C-ether formulations containing phosphorus ester and organic acid additives. Results were compared to those obtained with a formulated Type 2 ester and the C-ether base fluid. A ball-on-disk sliding friction apparatus was used. Experimental conditions were a 1-kilogram load, 17 meter-per-minute (100 rpm) surface speed, and a 25 to 300 C disk temperature range. The C-ether base fluid and the C-ether formulations yielded lower wear than the ester under most test conditions. The C-ether formulations exhibited higher friction coefficients than the ester from 150 to 300 C and similar or lower values from 25 to 150 C.

Jones, W. R., Jr.↗
NASA NTRS1975

Finite-element formulations for problems of large elastic-plastic deformation

An Eulerian finite element formulation is presented for problems of large elastic-plastic flow. The method is based on Hill's variational principle for incremental deformations, and is ideally suited to isotropically hardening Prandtl-Reuss materials. Further, the formulation is given in a manner which allows any conventional finite element program, for 'small strain' elastic-plastic analysis, to be simply and rigorously adapted to problems involving arbitrary amounts of deformation and arbitrary levels of stress in comparison to plastic deformation moduli. The method is applied to a necking bifurcation analysis of a bar in plane-strain tension. The paper closes with a unified general formulation of finite element equations, both Lagrangian and Eulerian, for large deformations, with arbitrary choice of the conjugate stress and strain measures. Further, a discussion is given of other proposed formulations for elastic-plastic finite element analysis at large strain, and the inadequacies of some of these are commented upon.

Mcmeeking, R. M.↗
NASA NTRS1975

Boundary lubrication of formulated C-ethers in air to 300 C

Six C-ether formulations (containing organic acid and phosphorus ester additives) were tested in dry (less than 100 ppm H2O) and moist air (relative humidity 50 percent at 25 C) and temperatures from 25 to 300 C, in order to determine the boundary lubricating characteristics and to compare these results to those obtained with a fully formulated Type II ester and the unformulated C-ether base fluid. The apparatus, procedures and detailed results of the test are described. Friction and wear measurements were made on CVM M-50 steel under 1 kg load and 17 meter per minute (100 rpm) surface speed conditions, for a duration of 25 minutes. In general, the C-ether formulations exhibited higher friction coefficients than the ester from 150 to 300 C and similar or lower values from 25 to 150 C, furthermore, in most test conditions the C-ether base fluid and the six C-ether formulations yielded lower wear than the ester. Possible applications in the aerospace industry are mentioned.

Jones, W. R., Jr.↗
NASA NTRS1977

Formulation of electrically conductive, thermal-control coatings

Formulation of electrically conductive, thermal-control coatings was undertaken for use on the International Sun Earth Explorer (ISEE) spacecraft. Unsuccessful formulation efforts as well as the successful use of oxide pigments fired at 1175 C are described. Problems attributed to reactivity of specific coating vehicles exposed to high humidity are discussed. Measurement and testing methods, including resulting data are mentioned, but the emphasis, in this report, is placed on coating formulation and application techniques. Methods of varying, as desired, optical properties are also described as well as formulations of white, low-absorptance coatings.

Shai, M. C.↗
NASA NTRS1980

Corotational velocity strain formulations for nonlinear analysis of beams and axisymmetric shells

Finite element formulations for large strain, large displacement problems are formulated using a kinematic description based on the corotational components of the velocity strain. The corotational components are defined in terms of a system that rotates with each element and approximates the rotation of the material. To account for rotations of the material relative to this element system, extra terms are introduced in the velocity strain equations. Although this formulation is incremental, in explicitly integrated transient problems it compares very well with formulations that are not.

Belytschko, T.↗
NASA NTRS1982

Improved boundary lubrication with formulated C-ethers

A comparison of five recently developed C-ether-formulated fluids with an advanced formulated MIL-L-27502 candidate ester is described. Steady state wear and friction measurements were made with a sliding pin on disk friction apparatus. Conditions included disk temperatures up to 260 C, dry air test atmosphere, 1 kilogram load, 50 rpm disk speed, and test times to 130 minutes. Based on wear rates and coefficients of friction, three of the C-ether formulations as well as the C-ether base fluid gave better boundary lubrication than the ester fluid under all test conditions. The susceptibility of C-ethers to selective additive treatment (phosphinic esters or acids and other antiwear additives) was demonstrated when two of the formulations gave somewhat improved lubrication over the base fluid. The increased operating potential for this fluid was shown in relationship to bulk oil temperature limits for MIL-L-23699 and MIL-L-27502 type esters.

Loomis, W. R.↗
NASA NTRS1982

Thermal and oxidative degradation studies of formulated C-ethers by gel-permeation chromatography

Gel-permeation chromatography was used to analyze C-ether lubricant formulations from high-temperature bearing tests and from micro-oxidation tests. Three mu-styragel columns (one 500 and two 100 A) and a tetrahydrofuran mobile phase were found to adequately separate the C-ether degradation products. The micro-oxidation tests yielded degradation results qualitatively similar to those observed from the bearing tests. Micro-oxidation tests conducted in air yielded more degradation than did tests in nitrogen. No great differences were observed between the thermal-oxidative stabilities of the two C-ether formulations or between the catalytic degradation activities of silver and M-50 steel. C-ether formulation I did yield more degradation than did formulation II in 111- and 25-hour bearing tests, respectively.

Jones, W. R., Jr.↗
NASA NTRS1982

A finite element formulation of Euler equations for the solution of steady transonic flows

The main objective of the considered investigation is related to the development of a relaxation scheme for the analysis of inviscid, rotational, transonic flow problems. To formulate the equations of motion for inviscid flows in a fixed coordinate system, an Eulerian type variational principle is required. The derivation of an Eulerian variational principle which is employed in the finite element formulation is discussed. The presented numerical method describes the mathematical formulation and the application of a numerical process for the direct solution of steady Euler equations. The development of the procedure as an extension of existing potential flow formulations provides the applicability of previous procedures, e.g., proper application of the artificial viscosity for supersonic elements, and the accurate modeling of the shock.

Ecer, A.↗