SEARCH · Search NASA

Results for “Variational 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2023

High-Order Mesh Morphing for Boundary and Interface Fitting to Implicit Geometries

Here, we propose a method that morphs high-order meshes such that their boundaries and interfaces coincide/align with implicitly defined geometries. Our focus is particularly on the case when the target surface is prescribed as the zero isocontour of a smooth discrete function. Common examples of this scenario include using level set functions to represent material interfaces in multimaterial configurations, and evolving geometries in shape and topology optimization. The proposed method formulates the mesh optimization problem as a variational minimization of the sum of a chosen mesh-quality metric using the Target-Matrix Optimization Paradigm (TMOP) and a penalty term that weakly forces the selected faces of the mesh to align with the target surface. The distinct features of the method are use of a source mesh to represent the level set function with sufficient accuracy, and adaptive strategies for setting the penalization weight and selecting the faces of the mesh to be fit to the target isocontour of the level set field. We demonstrate that the proposed method is robust for generating boundary- and interface-fitted meshes for curvilinear domains using different element types in 2D and 3D.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

A phase-field model of quasi-brittle fracture for pressurized cracks: Application to UO 2 high-burnup microstructure fragmentation

Here, in this paper, we present a phase-field model of quasi-brittle fracture with pressurized cracks, with dedicated applications for polycrystalline materials. The model is formulated as a minimization problem within the variational framework. The external work done by pressure on the crack surfaces is included in the objective function. Several careful modeling choices lead to a regularization-length-independent critical strength. The model is constructed to give a softening response with an underlying linear traction-separation law. The pressure-dependent softening response and the regularization of the prescribed pressure are demonstrated with a (quasi) one-dimensional numerical analysis. In a two-dimensional numerical analysis under plane strain assumptions, the critical stress corresponding to crack propagation (as predicted by our quasi-brittle fracture model) is compared with linear elastic fracture mechanics (LEFM) analytical solutions. Our model is further utilized to simulate fission-gas-induced fragmentation of the UO 2 high-burnup structure (HBS). Simulation results show that pressurized bubbles can cause crack nucleation and propagation, and that the bubble size and the surrounding external pressure affect the critical pressure corresponding to crack nucleation. Simulations of a partial HBS at different recrystallization stages show that different grain structures (due to recrystallization) also influence crack paths and fragmentation morphology.

36 MATERIALS SCIENCE↗

DOE OSTI2021

On a nonlocal Cahn–Hilliard model permitting sharp interfaces

A nonlocal Cahn–Hilliard model with a non-smooth potential of double-well obstacle type that promotes sharp interfaces in the solution is presented. To capture long-range interactions between particles, a nonlocal Ginzburg–Landau energy functional is defined which recovers the classical (local) model as the extent of nonlocal interactions vanish. In contrast to the local Cahn–Hilliard problem that always leads to diffuse interfaces, the proposed nonlocal model can lead to a strict separation into pure phases of the substance. In this work, the lack of smoothness of the potential is essential to guarantee the aforementioned sharp-interface property. Mathematically, this introduces additional inequality constraints that, in a weak formulation, lead to a coupled system of variational inequalities which at each time instance can be restated as a constrained optimization problem. We prove the well-posedness and regularity of the semi-discrete and continuous in time weak solutions, and derive the conditions under which pure phases are admitted. Moreover, we develop discretizations of the problem based on finite element methods and implicit–explicit time-stepping methods that can be realized efficiently. Finally, we illustrate our theoretical findings through several numerical experiments in one and two spatial dimensions that highlight the differences in features of local and nonlocal solutions and also the sharp interface properties of the nonlocal model.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Modeling of deformation of battery cells using thick shell element formulation

Here, we describe a new approach for modeling nonlinear deformation and stress distribution of battery cells using a new thick shell finite element formulation with a through-thickness calculation of stresses and strains that satisfy equilibrium conditions. Battery cells are transversely layered materials that contain numerous thin layers in a repeating sequence. The layers are made of materials with significantly different mechanical responses. Explicit discretization of the layers is computationally impractical except for very small domains, while homogenized material models cannot account for stress and strain variations and partition through the thickness. The new formulation allows for calculation of the transverse and interlayer stresses based on the material properties of the individual cell layers. The application to problems where the transverse stresses have strong influence is also possible.

25 ENERGY STORAGE↗

DOE OSTI2026

Variational neural network approach to QFT in the field basis

We present a variational neural network approach for solving quantum field theories in the field basis, focusing on the free Klein-Gordon model formulated in momentum space. While recent studies have explored neural-network-based variational methods for scalar field theory in position space, a systematic benchmark of the analytically solvable Klein-Gordon ground state—particularly in the momentum-space field basis—has been lacking. In this work, we represent the ground-state wavefunctional as a neural network defined on a discretized set of field configurations and train it by minimizing the Hamiltonian expectation value. This framework enables direct comparison to exact analytic results for a range of key observables, including the ground-state energy, two-point correlators, expectation value of the field, and the structure of the learned wavefunctional itself. Our results provide quantitative diagnostics of accuracy and establish a validated foundation for extending neural-network wavefunctional methods to interacting field theories and position-space formulations.

Klein-Gordon model↗

DOE OSTI2024

Lithium inventory tracking as a non-destructive battery evaluation and monitoring method

Tracking the active lithium (Li) inventory in an electrode shows the true state of a Li battery, akin to a fuel gauge for an engine. However, non-destructive Li inventory tracking is currently unavailable. Here, in this work, we used the theoretical capacity of a transition metal oxide to convert capacity into a Li inventory analysis. The Li inventory in electrodes was tracked reliably to show how battery formulations and test methods affect performance. Contrary to capacity, Li inventory tracking reveals stoichiometric variations near the electrode–electrolyte interface. Verifiable results rationalized differences in measurements, clarifying and reducing interferences from cell formulations and experimental manipulations. By tracing four variables from formation to end-of-life, we characterize electrode and cell performance with a thermodynamic framework. Accurate rationalization of subtle differences in Li inventory utilization promises precise battery engineering, evaluation, failure analysis and risk mitigation. The method could be applicable from cell design optimization and fabrication to battery management, improving battery performance and reliability.

25 ENERGY STORAGE↗

DOE OSTI2023

Biological toxicity and environmental hazards associated with PLGA nanoparticles

There have been tremendous advances in nanotechnology. More recently, the use of nanoparticles has expanded to applications including materials, packaging, energy, and medical uses such as drug delivery, diagnostics, and therapeutics. The American Society for Testing and Materials (ASTM) and the International Organization for Standardization (ISO) define nanoparticles as particles with at least one dimension measuring between 1 and 100 nm (2), but this view is generally too limited, and recommendations have been made to consider all materials with a dimension measuring under 1,000 nm as nanoparticles. The biomedical uses of nanoparticles are particularly promising because of their ability to reach and target various sites and organs. However, some nanoparticles can be composed of toxic materials or are limited by issues of biodistribution and bioaccumulation, which have hampered their use in biomedicine. The copolymer poly lactic-co-glycolic acid (PLGA) has gained use in biomedical applications as a delivery system because it is considered biocompatible and can be formulated with controlled degradation in physiological environments. PLGA has a history in biomedicine that dates to the 1970s when biodegradable sutures were developed using PLGA. In addition, the United States (US) Food and Drug Administration (FDA) and the European Medicine Agency (EMA) have approved various PLGA particle formulations as therapeutic delivery vehicles. Interest in biocompatible materials to deliver a range of therapeutic drugs, proteins, nucleic acids, and other molecules has risen recently. In addition to drug delivery, sustained drug release can be achieved by tuning the physical properties of PLGA, such as the molecular weight, ratio of lactic to glycolic acid, drug concentration, stabilizing molecules, and particle size. Given the rise of nanoparticles, human and environmental exposure to them is inevitable as more applications use free, unbound, and highly mobile particles. Although PLGA is generally considered safe, a detailed understanding of how PLGA and the other components used to help formulate this copolymer into a nanoparticle delivery vehicle is useful for assessing potential harmful effects. Despite recent advances in PLGA nanoparticle formulations, residual stabilizing molecules, inconsistent preparations, and batch-to-batch variations can lead to toxicity. Poor preparations of PLGA can lead to common mechanisms of toxicity observed with nanoparticles, such as inflammation and oxidative stress. Consistency in PLGA formulations and accurate methods to evaluate nanoparticle toxicity are needed to ensure safety. In order to better assess PLGA nanoparticles as non-toxic delivery systems, it is critical to understand their physical and chemical properties and the formulation protocols that may introduce toxic components into particles.

Biology, nanoparticle, toxicity↗

DOE OSTI2023

Variations in Creep Performance of P92 Steel Within Its Composition Range: Influence of N Solubility

Underspecified composition ranges often lead to alloys with unpredictable mechanical performance. To better understand the changes in microstructure and mechanical performance associated with variations of key elements, three versions of P92 are formulated within, or close to, the specified allowable for N, B, and C ranges. Chromium and Si are also varied to influence N solubility. Different service conditions (i.e., temperature and stress) are explored. It is observed that >80% decrease in creep life occurs at 625 °C and 155 MPa for the highest B and N containing alloy. Multiscale characterization reveals key changes due to the trace element variation. The high B and N containing alloy forms deleterious BN precipitates with morphology that promotes crack nucleation and damage accumulation, but this alloy additionally forms higher fractions of beneficial MX precipitates. The alloy with the lowest B and N concentrations but greater C content shows the best creep performance—a consequence of the refined M 23 C 6 carbide precipitate population and the absence of large-scale inclusions or BN precipitates. Calculations of creep activation energy reveal that the high B and N containing alloy is more prone to damage accumulation which causes an early onset of accelerating creep and greater minimum creep rate.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Effect of hydrogen on plasticity of $α$-Fe: A multi-scale assessment

A multi-scale study was carried out to quantify the effect of interstitial hydrogen concentration on plasticity in $α$-Fe. Here, in this work, the influence of hydrogen on the screw dislocation glide behavior was examined across several length-scales. The insights obtained were integrated to provide an accurate continuum description for the effect of hydrogen on the dislocation based plasticity in polycrystalline $α$-Fe. At the outset of this work, a new Fe—H interatomic potential was formulated that enhanced the atomistic estimation of the variation in dislocation glide behavior in presence of hydrogen. Next, the dislocation core reconstruction observed due to the addition of hydrogen using atomistic simulations was validated with the help of large-scale DFT calculations based on the DFT-FE framework. Several atomistic simulations were carried out to comprehensively quantify the effect of hydrogen on the non-Schmid behavior exhibited during the dislocation glide in $α$-Fe. Finally, crystal plasticity simulations were carried out to understand the effect of hydrogen on the meso-scale deformation behavior of polycrystalline $α$-Fe.

36 MATERIALS SCIENCE↗

DOE OSTI2022

A Bayesian Deep Learning Approach to Near-Term Climate Prediction

Since model bias and associated initialization shock are serious shortcomings that reduce prediction skills in state-of-the-art decadal climate prediction efforts, we pursue a complementary machine-learning-based approach to climate prediction. The example problem setting we consider consists of predicting natural variability of the North Atlantic sea surface temperature on the interannual timescale in the pre-industrial control simulation of the Community Earth System Model. While previous works have considered the use of recurrent networks such as convolutional LSTMs and reservoir computing networks in this and other similar problem settings, we currently focus on the use of feedforward convolutional networks. In particular, we find that a feedforward convolutional network with a Densenet architecture is able to outperform a convolutional LSTM in terms of predictive skill. Next, we go on to consider a probabilistic formulation of the same network based on Stein variational gradient descent and find that in addition to providing useful measures of predictive uncertainty, the probabilistic (Bayesian) version improves on its deterministic counterpart in terms of predictive skill. Finally, we characterize the reliability of the ensemble of machine learning models obtained in the probabilistic setting by using analysis tools developed in the context of ensemble numerical weather prediction.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Exposing and Reducing Biases of Simulating Mixed-Phase Clouds in the Convection-Permitting E3SM Atmosphere Model: Lessons From an Arctic Cold-Air Outbreak

Mixed-phase clouds modulate the water and energy cycles of high-latitude regions, yet their liquid-ice phase partitioning has long been poorly simulated in climate models. Here, simulations of Arctic mixed-phase clouds by the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) are assessed against large-eddy simulations, satellite data, and ground-based observations during the Cold-Air Outbreaks in the Marine Boundary Layer Experiment field campaign. SCREAM simulates nearly completely frozen clouds, which is attributed largely to the unreasonably strong Wegener–Bergeron–Findeisen (WBF) process that converts liquid to ice excessively and partly to the early over-abundant ice production at cold temperatures from a temperature-deterministic deposition ice nucleation scheme. Assuming no subgrid variation for the WBF process in the original formulation particularly conflicts with the instantaneous saturation adjustment assumption in the condensation scheme that assumes subgrid variability, leading to exaggerated WBF process rates. A proposed simple physically-based improvement on the treatment of subgrid cloud overlap substantially increases supercooled liquid water content and notably improves cloud-top phase partitioning, aligning better with observations. Improvement of supercooled liquid water content also converges with increasing horizontal resolution. The deposition ice nucleation scheme is found responsible for a falsely-produced ice cloud aloft that is not observed, biasing the simulated cloud radiative effects and top-of-atmosphere radiative fluxes. This study identifies key deficiencies in cloud parameterizations that continue to challenge convection-permitting models.

Geosciences↗

DOE OSTI2025

In‐situ Analysis of Paste Properties in Resonant Acoustic Mixers for Quality Monitoring

Formulation control is key to achieving consistent target properties of energetic materials, as feedstock variations and slight deviations in the ratios of different ingredients can have major effects on final product properties, particularly in dense pastes with high particle loading >65 vol.%. In large‐scale operations, it is imperative to either correct or remove batches of material that perform outside baseline property specifications as early as possible to avoid unnecessary processing of suboptimal material. Quality monitoring is the practice of measuring material properties during processing using process analytical technologies as opposed to only testing the properties of the final product; it is a key principle in the quality‐by‐design frameworks used for designing formulations and manufacturing processes. Herein, a process analytical technology method for correlating material properties of dense pastes directly after mixing in a Resonant Acoustic Mixer to motor data is developed and used to detect differences in the particle content of dense paste formulations. This method was also capable of detecting variations in powder feedstock properties, such as particle packing efficiency, and is sensitive enough to detect changes of 2 wt.% in the total solids content of the formulation. The techniques presented herein show excellent promise for use as a process analytical technology capable of quantifying formulation effects on material movement modes during resonant acoustic mixing.

Materials science↗

DOE OSTI2020

Laminar-to-turbulent flame transition and cycle-to-cycle variations in large eddy simulation of spark-ignition engines

This paper investigates the effect of laminar-to-turbulent flame transition modeling on the prediction of cycle-to-cycle variations (CCVs) in large eddy simulation (LES) of spark-ignition (SI) engines. A laminar-to-turbulent flame transition model that describes the non-equilibrium sub-filter flame speed evolution during an early stage of flame kernel growth is developed. In the present model, the flame transition is characterized by the flame kernel size at which the flame transition ends, defined here as the flame transition scale. The proposed model captures the effects that variations in a turbulent flow field have on the evolution of early-stage burning rates, through variations in the flame transition scale. The proposed flame transition model is combined with the front propagation formulation (FPF) method and a spark-ignition model to predict CCVs in a gasoline direct injection SI engine. It is found that multi-cycle LES with the proposed flame transition model reproduces experimentally-observed CCVs satisfactorily. When the transition model is not considered or when variations in the transition process are neglected, CCVs are significantly under-predicted for the case considered here. These results indicate the importance of modeling the laminar-to-turbulent flame transition and the effect of turbulence on the transition process, when predicting CCVs, under certain engine conditions. The LES results are also used to analyze sources for variations in the flame transition. It is found, for the present engine case, that the most important source is the cycle-to-cycle variation in the turbulence dissipation rate, which is used to measure the strength of turbulence in the proposed model, near a spark plug. The large-scale velocity field and the variations of the laminar flame speed due to the mixture composition and thermal stratification are also found to be important factors to contribute to the variations in the flame transition.

42 ENGINEERING↗

DOE OSTI2023

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

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

42 ENGINEERING↗

DOE OSTI2021

The quest for efficient oxygenated fuels: Examining interactions between lubricant components and oxygenates

This work empirically evaluated oxygenated blendstock-lubricant compatibility using a novel but simple method. The reactivity or compatibility was evaluated via gas chromatograph equipped with a flame ionization detector (GC-FID) before and after heating for 2h at 150 °C, which would roughly mimic conditions encountered in the combustion chamber of an Spark Ignition (SI) or Compression-Ignition (CI) engine. Five oxygenate blendstocks spanning several functional groups and chain lengths (prenol, isoprenol, butyl acetate, isohexanol, polyoxymethylene dimethyl ethers) were evaluated in a base fuel with nine variations of lubricant components of controlled concentrations, including fully formulated commercial lubricants. Baseline comparative examples were also generated, to decouple the effect of the oxygenate on the mixture, by heating only the base fuel with the lubricant component. Results demonstrate the relative compatibility between noted lubricants and oxygenates in base fuel upon heating. This was illustrated by minimal to no changes in the GC-FID traces. However, reactions of several of the chosen lubricants, especially Molybdenum di (2-ethylhexyl) phosphorodithioate (Molyvan L) and zinc dialkyldithiophosphate (ZDDP), with prenol and polyoxymethylene dimethyl ethers blendstocks showed substantial side reactions compared to corresponding baseline examples. Control experiments without oxygenates but comprising the same lubricant components showed no change in the mixture with heating. We concluded the incompatibility is due to the oxygenate and not the components of the base fuel. This preliminary work is valuable in establishing trends as well as provide useful information when a new component is introduced in the fuel.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Reliability-Informed Economic and Energy Evaluation for Design for Remanufacturing: A Case Study on a Hydraulic Manifold

Design for Remanufacturing (DfRem) is an attractive approach for sustainable product development. Evaluation of DfRem strategies, from both economic and environmental perspectives, at an early design stage can allow the designers to make informed decisions when choosing the best design option. Studying the long-term implications of a particular design scenario requires quantifying the benefits of remanufacturing for multiple life cycles while considering the reliability of the product. In addition to comparing designs on a one-to-one basis, we find that including reliability provides a different insight into comparing design strategies. We present a reliability-informed cost and energy analysis framework that accounts for product reliability for multiple remanufacturing cycles within a certain warranty policy. The variation of reuse rate over successive remanufacturing cycles is formulated using a branched power-law model which provides probabilistic scenarios of reusing or replacing with new units. To demonstrate the utility of this framework, we use the case study of a hydraulic manifold, which is a component of a transmission used in some agricultural equipment, and use real-world field reliability data to quantify the transmission’s reliability. Three design improvement changes are proposed for the manifold and we quantify the costs and energy consumption associated with each of the design changes for multiple remanufacturing cycles.

design for remanufacturing↗

DOE OSTI2021

Bayesian Spatially Varying Multi-Regularization Image Deblurring

Many scientific experiments such as those found in astronomy, geology, microbiology, and X-ray radiography require the use of high-energy instruments to capture images. Since blur and noise are inevitably present in any imaging system, the images must be \deblurred" to extract the full information content. Mathematically, image deblurring is an ill-posed inverse problem that requires regularization. The regularization, in turn, has a large effect on the deblurred image: different regularization strengths, and types, lead to drastically different reconstructions. Moreover, many images contain a mixture of smooth and sharp features which suggests the use of multi-regularization, i.e., varying the type of regularization (e.g. Tikhonov or total variation) across the image. We address these issues by formulating the image deblurring problem within a hierarchical Bayesian framework in which we spatially adapt the strength of the regularization and also vary the regularization type across the image. In this way, the image itself, along with corresponding regularization strength at each pixel, are described jointly by a posterior distribution which we can sample by Markov chain Monte Carlo (MCMC) methods. We illustrate our techniques on simplified test problems and apply them to high-energy X-ray images taken at the Nevada National Security Site. Numerical tests show that our new method is robustly applicable and increases the quality of the image reconstruction when compared to other (Bayesian) methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Solving Unit Commitment Problems with Demand Responsive Loads

This work focuses on using variations of the Frank-Wolfe (FW) algorithm for solving unit commitment problems with high volumes of demand responsive loads on the power grid. We present a formulation of the unit commitment problem with demand responsive loads. We then show through reformulation and relaxations of the problem that variations of the Frank-Wolfe algorithm can be used to determine the time series decisions for the demand responsive loads. We show through computational experiments on the IEEE Reliability Test System that the timeseries of demand responsive load decisions obtained through our approach are near optimal and describe how large-scale parallel implementations of our approach can be highly computationally efficient.

demand response↗