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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Slow Crack Growth of Brittle Materials With Exponential Crack-Velocity Formulation: Analysis - Part 1

Extensive slow-crack-growth (SCG) analysis was made using a primary exponential crack-velocity formulation under three widely used load configurations: constant stress rate, constant stress, and cyclic stress. Although the use of the exponential formulation in determining SCG parameters of a material requires somewhat inconvenient numerical procedures, the resulting solutions presented gave almost the same degree of simplicity in both data analysis and experiments as did the power-law formulation. However, the fact that the inert strength of a material should be known in advance to determine the corresponding SCG parameters was a major drawback of the exponential formulation as compared with the power-law formulation.

Choi, Sung R.↗
NASA NTRS2002

Slow Crack Growth of Brittle Materials With Exponential Crack-Velocity Formulation: Constant Stress and Cyclic Stress Experiments - Part 3

The previously determined life prediction analysis based on an exponential crack-velocity formulation was examined using a variety of experimental data on advanced structural ceramics tested under constant stress and cyclic stress loading at ambient and elevated temperatures. The data fit to the relation between the time to failure and applied stress (or maximum applied stress in cyclic loading) was very reasonable for most of the materials studied. It was also found that life prediction for cyclic stress loading from data of constant stress loading in the exponential formulation was in good agreement with the experimental data, resulting in a similar degree of accuracy as compared with the power-law formulation. The major limitation in the exponential crack-velocity formulation, however, was that the inert strength of a material must be known a priori to evaluate the important slow-crack-growth (SCG) parameter n, a significant drawback as compared with the conventional power-law crack-velocity formulation.

Choi, Sung R.↗
NASA NTRS2005

Completed Beltrami-Michell Formulation in Polar Coordinates

A set of conditions had not been formulated on the boundary of an elastic continuum since the time of Saint-Venant. This limitation prevented the formulation of a direct stress calculation method in elasticity for a continuum with a displacement boundary condition. The missed condition, referred to as the boundary compatibility condition, is now formulated in polar coordinates. The augmentation of the new condition completes the Beltrami-Michell formulation in polar coordinates. The completed formulation that includes equilibrium equations and a compatibility condition in the field as well as the traction and boundary compatibility condition is derived from the stationary condition of the variational functional of the integrated force method. The new method is illustrated by solving an example of a mixed boundary value problem for mechanical as well as thermal loads.

Patnaik, Surya N.↗
NASA NTRS2012

A New and General Formulation of the Parametric HFGMC Micromechanical Method for Three-Dimensional Multi-Phase Composites

The recent two-dimensional (2-D) parametric formulation of the high fidelity generalized method of cells (HFGMC) reported by the authors is generalized for the micromechanical analysis of three-dimensional (3-D) multiphase composites with periodic microstructure. Arbitrary hexahedral subcell geometry is developed to discretize a triply periodic repeating unit-cell (RUC). Linear parametric-geometric mapping is employed to transform the arbitrary hexahedral subcell shapes from the physical space to an auxiliary orthogonal shape, where a complete quadratic displacement expansion is performed. Previously in the 2-D case, additional three equations are needed in the form of average moments of equilibrium as a result of the inclusion of the bilinear terms. However, the present 3-D parametric HFGMC formulation eliminates the need for such additional equations. This is achieved by expressing the coefficients of the full quadratic polynomial expansion of the subcell in terms of the side or face average-displacement vectors. The 2-D parametric and orthogonal HFGMC are special cases of the present 3-D formulation. The continuity of displacements and tractions, as well as the equilibrium equations, are imposed in the average (integral) sense as in the original HFGMC formulation. Each of the six sides (faces) of a subcell has an independent average displacement micro-variable vector which forms an energy-conjugate pair with the transformed average-traction vector. This allows generating symmetric stiffness matrices along with internal resisting vectors for the subcells which enhances the computational efficiency. The established new parametric 3-D HFGMC equations are formulated and solution implementations are addressed. Several applications for triply periodic 3-D composites are presented to demonstrate the general capability and varsity of the present parametric HFGMC method for refined micromechanical analysis by generating the spatial distributions of local stress fields. These applications include triply periodic composites with inclusions in the form of a cavity, spherical inclusion, ellipsoidal inclusion, discontinuous aligned short fiber. A 3-D repeating unit-cell for foam material composite is simulated.

Haj-Ali, Rami↗
NASA NTRS2015

Compact Assumption Applied to the Monopole Term of Farassat's Formulations

Farassat's formulations provide an acoustic prediction at an observer location provided a source surface, including motion and flow conditions. This paper presents compact forms for the monopole term of several of Farassat's formulations. When the physical surface is elongated, such as the case of a high aspect ratio rotorcraft blade, compact forms can be derived which are shown to be a function of the blade cross sectional area by reducing the computation from a surface integral to a line integral. The compact forms of all formulations are applied to two example cases: a short span wing with constant airfoil cross section moving at three forward flight Mach numbers and a rotor at two advance ratios. Acoustic pressure time histories and power spectral densities of monopole noise predicted from the compact forms of all the formulations at several observer positions are shown to compare very closely to the predictions from their non-compact counterparts. A study on the influence of rotorcraft blade shape on the high frequency portion of the power spectral density shows that there is a direct correlation between the aspect ratio of the airfoil and the error incurred by using the compact form. Finally, a prediction of pressure gradient from the non-compact and compact forms of the thickness term of Formulation G1A shows that using the compact forms results in a 99.6% improvement in computation time, which will be critical when noise is incorporated into a design environment.

Lopes, Leonard V.↗
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↗
DOE OSTI2020

In-Use Interfacial Stability of Monoclonal Antibody Formulations Diluted in Saline i.v. Bags

The use of monoclonal antibodies (mAbs) for the treatment of a variety of diseases is rapidly growing each year. Many mAbs are administered intravenously using i.v. bags containing 0.9% NaCl (normal saline). We studied the aggregation propensity of these antibody solutions in saline and compared it with a low ionic strength formulation buffer. Furthermore, the mAb studied in this work is prone to aggregate, and is known to form a viscoelastic network at the air-solution interface. We observed that this interfacial elasticity increased when formulated in saline. In the bulk, the mAbs exhibited a tendency to self-associate that was higher in saline. We also studied the aggregation of the mAbs in the presence of polysorbate-20, typically added to formulations to mitigate interfacial aggregation. We observed that with surfactants, the presence of salt in the buffer led to a greater mAb adsorption at the interface and resulted in the formation of more particulate aggregates. Our results show that the addition of salt to the buffer led to differences in the interfacial aggregation in mAb formulations, showing that stress studies used to screen for mAb aggregation intended for i.v. administration should be performed in conditions representative of their intended route of administration.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

On the cogent formulation of an elastomeric silicone ink material for direct ink write ( DIW ) 3D printing

Abstract Adhesives and sealants show fine rheology with good physical and mechanical properties as viscous pastes, a possible starting point for developing direct ink writing (DIW) 3D printing ink. However, many commercial adhesives and sealants take days or weeks to cure fully. DIW 3D‐printed parts made directly from these sealants are not designed for a scalable manufacturing process and high‐volume production. Moreover, most of these adhesives and sealants have volume shrinkage during cure. A systematic understanding of formulation methods and design principles for an elastomeric silicone DIW ink can overcome these issues. This study presents the cogent formulation development of a 3D printable thermoset elastomer silicone that gels and cures isotropically in minutes, reducing cycle time for rapid ink development with no shrinkage during cure. More specifically, we outline the principles of raw material selection of a formulation to achieve excellent rheology, printability, synchronized working, and gel time fitting requirements closer to scalable manufacturing. The reaction kinetics and their corresponding 3D‐printed structural properties are also described. Interest in future work is toward a rational DIW 3D printing ink material development protocol and use of machine learning (ML). Highlights Formulation method flexibility and design principle of DIW ink. Raw material selection principle to achieve optimal rheology for DIW printing. Ink gel kinetics for large‐scale DIW manufacturing. Hydrosilylation conversion over time at different ambient temperatures. Structural properties of DIW 3D printed parts.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Review of the Effects of Polymer Binder Properties on Microstructure and Irreversible Volume Growth of Plastic Bonded Explosives Formulations

The rational design of effective polymeric binders for the formulation of plastic-bonded explosives (PBX) is challenging due to their inherent compositional complexity. The composites comprise irregularly shaped energetic material (EM) powders coated with low weight fractions of polymer via non-equilibrium processes such as slurry coating. Defects can deleteriously affect PBX stability and performance: nano- to micrometer-scale voids can act as loci for hot spots, lowering deflagration and detonation temperatures in unpredictable ways. Furthermore, some nominally desirable polymer properties are at odds with each other: e. g. good flow characteristics are desirable for coating and adhesion, but mechanical stiffness is needed to prevent deformation and cracking of PBX under mechanical stress. Good binder adhesion is critical, but the best means to predict and measure adhesion in PBX is not obvious. Experimental methods of determining binder adhesion on model surfaces may not capture polymer structural configurations relevant to deposition during coating. Molecular dynamics-based computational models have predicted key observables in PBX formulation, suggesting that they may be powerful tools for binder selection. In this review, primarily recent (~2006 and later) literature on polymeric binders for insensitive HE (IHE) is surveyed. We focus on how binder properties influence observable PBX properties as resistance to irreversible volume growth and void formation in PBX formulations mainly (but not exclusively) featuring 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) as the EM (e. g. PBX 9502, LX-17 and others). Conclusions from these studies yield useful guidelines for choosing HE binder candidates, as well as for the general design of highly-filled polymer composite materials. Finally, studies describing challenges in PBX formulation with the newer and more energetically dense high explosive, LLM-105, will be discussed.

adhesion↗
DOE OSTI2021

Computational and numerical analysis of AC optimal power flow formulations on large-scale power grids

We report that alternating current optimal power flow (AC-OPF) is a fundamental tool in electric utilities to determine optimal operation of the various resources. Typically, the AC-OPF problem uses power balance formulation containing voltages and power equations. Yet, there is no comprehensive comparison of the different AC-OPF formulations, especially for large-scale networks. This paper presents a detailed comparative evaluation of different formulations of the AC-OPF problem on networks ranging from 9-bus to 25,000 buses. Three different formulations: 1) power balance with polar voltages, 2) power balance with Cartesian voltages, and 3) current balance with Cartesian voltages are discussed in detail by comparing their characteristics, and numerical and computational performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2013

Glass Property Models and Constraints for Estimating the Glass to be Produced at Hanford by Implementing Current Advanced Glass Formulation Efforts

Recent glass formulation and melter testing data have suggested that significant increases in waste loading in HLW and LAW glasses are possible over current system planning estimates. The data (although limited in some cases) were evaluated to determine a set of constraints and models that could be used to estimate the maximum loading of specific waste compositions in glass. It is recommended that these models and constraints be used to estimate the likely HLW and LAW glass volumes that would result if the current glass formulation studies are successfully completed. It is recognized that some of the models are preliminary in nature and will change in the coming years. Plus the models do not currently address the prediction uncertainties that would be needed before they could be used in plant operations. The models and constraints are only meant to give an indication of rough glass volumes and are not intended to be used in plant operation or waste form qualification activities. A current research program is in place to develop the data, models, and uncertainty descriptions for that purpose. A fundamental tenet underlying the research reported in this document is to try to be less conservative than previous studies when developing constraints for estimating the glass to be produced by implementing current advanced glass formulation efforts. The less conservative approach documented herein should allow for the estimate of glass masses that may be realized if the current efforts in advanced glass formulations are completed over the coming years and are as successful as early indications suggest they may be. Because of this approach there is an unquantifiable uncertainty in the ultimate glass volume projections due to model prediction uncertainties that has to be considered along with other system uncertainties such as waste compositions and amounts to be immobilized, split factors between LAW and HLW, etc.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Glass Property Models, Constraints, and Formulation Approaches for Vitrification of High-Level Nuclear Wastes at the US Hanford Site

The legacy nuclear wastes stored in underground tanks at the US Department of Energy’s Hanford site is planned to be separated into high-level waste and low-activity waste fractions and vitrified separately. Formulating optimized glass compositions that maximize the waste loading in glass is critical for successful and economical treatment and immobilization of nuclear wastes. Glass property-composition models have been developed and applied to formulate glass compositions for various objectives for the past several decades. The property models with associated uncertainties and combined with composition and property constraints have been used to develop preliminary glass formulation algorithms designed for vitrification process control and waste form qualification at the planned waste vitrification plant. This paper provides an overview of current status of glass property-composition models, constraints applicable to Hanford waste vitrification, and glass formulation approaches that have been developed for vitrification of hazardous and highly radioactive wastes stored at the Hanford site.

composition constraint↗
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 OSTI2011

Glass Formulation for Next Generation Cold Crucible Induction Melter

Transformational melter technologies are being considered to support mission acceleration within the U.S. Department of Energy (DOE) complex. New glass formulations are required to take full advantage of the next generation melters, for example, the cold crucible induction melter (CCIM). The key advantage of CCIM technology over current reference technologies is its capability to provide higher processing temperatures, which can lead to an increased waste throughput rate by achieving higher waste loadings and by increasing the feed processing rate. Various waste compositions within the DOE complex were evaluated to determine their potential for successfully demonstrating the unique advantages of the CCIM technology. Glass formulations that satisfy a set of constraints for product quality and assumed CCIM processing conditions were developed for two Hanford waste streams, AZ-101 high-level waste (HLW) and AN-105 low-activity waste (LAW). Three glasses selected for AZ-101 HLW have waste loadings of 40, 42.5, and 45 wt%. The 45-wt% waste loading corresponds to a 22% increase from 37 wt%, which is the maximum expected waste loading based on the current reference formulation. One glass selected for AN-105 LAW has a waste loading of 31.3 wt% at 24 wt% Na2O in glass, which is a 14% increase from the current reference formulation maximum of 21 wt% Na2O. These four glasses are planned for scaled melter tests for initial demonstration of the CCIM technologies for Hanford wastes.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

Synthesis and Evaluation of 3-Picrylamino-1,2,4-Triazole (PATO) Formulations

3–Picrylamino–1,2,4–triazole (PATO) is a thermally stable explosive invented at Los Alamos National Laboratory (LANL) almost half a century ago. Despite a rapid and high yielding synthesis, performance data for this promising explosive are scant. We prepared material using Coburn's synthesis, and discovered that the particle size distribution and morphology leads to difficulty in pressing formulations to high density. Three formulations were made using glycidyl azide polymer (GAP), FK–800, and Estane 5703/nitroplasticizer (NP) as binders. The maximum pressed density of these formulations was 89 % of theoretical maximum density (TMD). We determined detonation velocity and, for the FK–800 formulation, detonation pressure. Furthermore, the performance of PATO is similar to TATB at equivalent pressed density.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2018

Variational formulation of classical and quantum models for intense laser pulse propagation

Highlights: •We consider the theoretical description of intense laser pulses propagating through gases. •We derived reduced models using Lagrangian and Hamiltonian formulations. •We obtain models where the atomic electrons are treated classically or quantum-mechanically. •We derive a number of conserved quantities of the reduced models. -- Abstract: We consider the theoretical description of intense laser pulses propagating through gases. Starting from a first-principles description of both the electromagnetic field and the electron motion within the gas atoms, we derive a hierarchy of reduced models. We obtain a parallel set of models, where the atomic electrons are treated classically on the one hand, and quantum-mechanically on the other. By working consistently in either a Lagrangian formulation or a Hamiltonian formulation, we ensure that our reduced models preserve the variational structure of the parent models. Taking advantage of the Hamiltonian formulation, we deduce a number of conserved quantities of the reduced models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗