Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constitutive models”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Computational Analysis of Different Sparging Systems and their Influence in the Fluid-Dynamic Behavior of Bubble Column Reactors

Bubble column bioreactors are being actively considered for gas fermentation applications, specifically for CO2 utilization, and sugars to fuels conversion. Their main advantages include good mass transfer without any moving parts and low-cost of operation and maintenance. However, the design and scale-up of such reactors is challenging specifically for carbon capture applications where a mixture of gases (e.g. CO2/CO/H2) with variable solubilities is used. The overall performance of scaled-up bioreactors (e.g., mass transfer rate) is largely affected by gas holdup, bubble size distribution (BSD), and multiphase hydrodynamics. We investigate the effect of gas sparger designs on the performance of these large-scale bioreactors using computational fluid dynamics simulations in this work, so as to improve CO2 conversion at scale. The gas distribution systems in bubble column reactors not only determines operational regime, but also affects the evolution of the BSD, which in turn influences interfacial mass transfer and ultimately the efficiency of the gas-liquid exchange process. In addition to the BSD, uniformity in gas sparging affects gas holdup and bubble residence time which constitute important metrics of performance in gas-liquid systems. In this work, we use computational models to simulate high fidelity representations of different sparger designs and their effect on the operation of a bubble column reactor. Four different types of spargers have been selected for the computational study (Fig. 1): ladder, multi-ring, single-ring and toroidal. Their effect on superficial velocity, gas holdup mixing efficiency, and BSD will be evaluated in this work. The model uses a multiphase Eulerian framework similar to [1] and include a composition of mixtures of H2/CO/CO2 gases, common in fermentation applications.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Protein cost minimization promotes the emergence of coenzyme redundancy

Metabolism relies on a small class of molecules (coenzymes) that serve as universal donors and acceptors of key chemical groups and electrons. Although metabolic networks crucially depend on structurally redundant coenzymes [e.g., NAD(H) and NADP(H)] associated with different enzymes, the criteria that led to the emergence of this redundancy remain poorly understood. Our combination of modeling and structural and sequence analysis indicates that coenzyme redundancy may not be essential for metabolism but could rather constitute an evolved strategy promoting efficient usage of enzymes when biochemical reactions are near equilibrium. Our work suggests that early metabolism may have operated with fewer coenzymes and that adaptation for metabolic efficiency may have driven the rise of coenzyme diversity in living systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Model Package Report: Composite Analysis Solid Waste Release Model (CASWR Model)

This document describes the implementation of a solid waste form release model in GoldSim for the Hanford Site Composite Analysis (CA) Update. This Composite Analysis Solid Waste Release model (CASWR model) was designed to generate deterministic radionuclide release rates for Hanford’s Central Plateau solid waste disposal sites using single realizations of release model coefficients. Five generalized waste form types are used for the conceptual model of waste release: surplus reactor block, cement, soil-debris, grouted residual waste, and ancillary equipment. The surplus reactor block waste form consists of radionuclide leaching from surplus production graphite reactor core blocks via unspecified processes (White et al., 1984 as cited in PNNL-15965). The cement waste form represents solidified wastes whose permeability is much lower than that of the surrounding soil. The soil-debris waste form type is defined as unconsolidated waste mixed with soil material. Tanks and canyon complexes comprise the grouted residual waste form such that their solid waste will be grouted and capped with a surface barrier at the completion of their cleanup. Finally, the ancillary equipment waste form constitutes contaminant releases from ancillary and auxiliary waste form residues associated with tank farms at closure. Individual sub-models are implemented to numerically represent a respective waste form within the CASWR Model: Surplus Reactor Block Sub-model, Cement Sub-model, SoilDebris Sub-model, Grouted Residual Waste Sub-model, and Ancillary Equipment Submodel. Advection is assumed to be the primary transport process governing the release of radionuclides in the Ancillary Equipment and Soil-Debris Sub-models. Diffusion is assumed to be the primary release process in the Grouted Residual Waste and Cement Sub-models. An unspecified zero-order release process is considered in the Surplus Reactor Block Sub-model due to the lack of information regarding actual processes involved in irradiated graphite leaching. The Surplus Reactor Block, Cement, and SoilDebris Sub-models were compared against analytical solutions (PNNL-11800, Composite Analysis for Low-Level Waste Disposal in the 200 Area Plateau of the Hanford Site). The agreement between results of these analytical solutions and the corresponding waste form models verified their correct implementation in GoldSim. A 1-D modeling abstraction approach for the Grouted Residual Waste and Ancillary Equipment Sub-models was adopted from existing Performance Assessment (PA) models. Despite the simplifications made in these sub-models, they were found to be appropriate representations of the waste forms, similar to what was used in the Waste Management Area C PA model (RPP-ENV-58782, Performance Assessment of Waste Management Area C, Hanford Site, Washington, Rev. 0). A sensitivity analysis was conducted to identify the most influential parameters in each waste form sub-model. Suggestions for considering pH-dependent and redox-dependent release mechanisms are formulated through the development of a conditional constant approach in GoldSim.

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

A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage

Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Shear Band Formation in Thin-Film Multilayer Columns Under Compressive Loading: A Mechanistic Study

Micro-pillar compression is a popular experimental technique used for characterizing the mechanical behavior of nano- and micro-laminates. The compressive stress–strain response of the column-shaped thin-film composite can be measured, and the deformation and damage features can be revealed by post-test cross-section microscopy. The development of plastic instability in the form of localized strain concentration (shear bands), leading to eventual failure, is frequently observed. In the present study, a computational approach is used to illustrate the commonality of shear band formation from a continuum standpoint. Systematic finite element analyses are conducted, showing that the strain field tends to become localized once plastic yielding commences. Distinct shear offsets of the layered structure can be revealed from the numerical model, which is similar to those observed in experiments. The actual appearance of shear bands depends on the materials’ constitutive behavior and precise geometries. Post-yield strain hardening reduces the propensity of shear band formation, while strain softening enhances it. Imperfections such as the undulated layer geometry, as well as the frictional characteristics between the specimen and test apparatus, can also influence the shear band morphology and overall stress–strain response.

finite element modeling↗
DOE OSTI2025

Thermal shock resistance of lightweight cements developed for geothermal conditions

Cements are a critical component in well construction, as they act to prevent well fluid and gas escape, prevent corrosion of the casing, and strengthen the wellbore to prevent deformation. Under the high temperature/pressure conditions common in geothermal systems, the injection of cold water for energy production is expected to induce cyclic damage to the borehole cement through the rapid temperature fluctuations. These “thermal shocks” are expected to cause casing shrinkage, annulus formation, and cement tensile stresses. To understand the effect of cold water injection on the wellbore environment, a set of rock-cement-steel samples were created to simulate the structure of a geothermal well. Lightweight thermally-insulating cement blends were tested under thermal shock conditions in this study. In each test, the samples were pressurized to an effective pressure of ~3.5 MPa and placed at high temperatures. Thermal shocks were performed by injecting cold water (~10-15 °C) through the samples at a constant rate while keeping the samples at high temperatures until the sample temperature stopped decreasing and deformation ceased. Eight thermal shock tests were conducted with each sample – two at 100 °C and six at 200 °C. Post-tests analysis was then conducted by cutting open each sample to examine the damage in each component of the simulated wellbore. Experimental results suggest that all samples experienced similar degrees of axial and lateral contraction during cold water injection, but for the most part this contraction is recoverable when injection halts. Post-test analysis revealed that fly ash cenosphere pre-treatment had the best effect on improving thermal shock resistance in the cement blends. Thermomechanical modeling of likely stress paths experienced by the cements during heating/cooling cycles shows that elasto-plastic cement constitutive behavior results in most plastic strain occurring during the initial heating steps, with mostly elastic strain occurring during the thermal shock cycles. In conclusion, this agrees with experimental evidence, suggesting that cement damage from shocking occurs via other mechanisms such as chemical alteration, corrosion, and fatigue.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Core formation in high-z massive haloes: heating by post-compaction satellites and response to AGN outflows

ABSTRACT Observed rotation curves in star-forming galaxies indicate a puzzling dearth of dark matter in extended flat cores within haloes of mass ≥1012M⊙ at z ∼ 2. This is not reproduced by current cosmological simulations, and supernova-driven outflows are not effective in such massive haloes. We address a hybrid scenario where post-compaction merging satellites heat up the dark-matter cusps by dynamical friction, allowing active galactic nucleus (AGN)-driven outflows to generate cores. Using analytic and semi-analytic models (SatGen), we estimate the dynamical friction heating as a function of satellite compactness for a cosmological sequence of mergers. Cosmological simulations (VELA) demonstrate that satellites of initial virial masses >1011.3M⊙, which undergo wet compactions, become sufficiently compact for significant heating. Constituting a major fraction of the accretion on to haloes ≥1012M⊙, these satellites heat up the cusps in half a virial time at z ∼ 2. Using a model for outflow-driven core formation (CuspCore), we demonstrate that the heated dark-matter cusps develop extended cores in response to removal of half the gas mass, while the more compact stellar systems remain intact. The mergers keep the dark matter hot, while the gas supply, fresh and recycled, is sufficient for the AGN outflows. AGNs indeed become effective in haloes ≥1012M⊙, where the black hole growth is no longer suppressed by supernovae and its compaction-driven rapid growth is maintained by a hot circumgalactic medium. For simulations to reproduce the dynamical friction effects, they should resolve the compaction of the massive satellites and avoid artificial tidal disruption. AGN feedback could be boosted by clumpy black hole accretion and clumpy response to AGN.

Dekel, Avishai↗
DOE OSTI2022

Exotic Higgs Decays

Exotic decays of the Standard Model (SM)-like Higgs boson into beyond-the-SM particles are predicted in a wide range of well-motivated theories. The enormous samples of Higgs bosons that have been and will be produced at the Large Hadron Collider thus constitute one of the key discovery opportunities at that facility, particularly in the upcoming high-statistics, high-luminosity run. Here we review recent theoretical work on models that predict or accommodate exotic Higgs decays, survey the status of current experimental searches, and look forward to future capabilities at dedicated Higgs factories and beyond.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Multi-body entanglement and information rearrangement in nuclear many-body systems: a study of the Lipkin–Meshkov–Glick model

Here, we examine how effective-model-space (EMS) calculations of nuclear many-body systems rearrange and converge multi-particle entanglement. The generalized Lipkin-Meshkov-Glick (LMG) model is used to motivate and provide insight for future developments of entanglement-driven descriptions of nuclei. The effective approach is based on a truncation of the Hilbert space together with a variational rotation of the qubits (spins), which constitute the relevant elementary degrees of freedom. The non-commutivity of the rotation and truncation allows for an exponential improvement of the energy convergence throughout much of the model space. Our analysis examines measures of correlations and entanglement, and quantifies their convergence with increasing cut-off. We focus on one- and two-spin entanglement entropies, mutual information, and $n$-tangles for $n=2,4$ to estimate multi-body entanglement. The effective description strongly suppresses entropies and mutual information of the rotated spins, while being able to recover the exact results to a large extent with low cut-offs. Naive truncations of the bare Hamiltonian, on the other hand, artificially underestimate these measures. The $n$-tangles in the present model provide a basis-independent measures of $n$-particle entanglement. While these are more difficult to capture with the EMS description, the improvement in convergence, compared to truncations of the bare Hamiltonian, is significantly more dramatic. We conclude that the low-energy EMS techniques, that successfully provide predictive capabilities for low-lying observables in many-body systems, exhibit analogous efficacy for quantum correlations and multi-body entanglement in the LMG model, motivating future studies in nuclear many-body systems and effective field theories relevant to high-energy physics and nuclear physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Validation of the Community Land Model Version 5 over the Contiguous United States (CONUS) using in situ and remote sensing data sets

The Community Land Model (CLM) is an effective tool to simulate the biophysical and biogeochemical processes and their interactions with the atmosphere. Although CLM Version 5 (CLM5) constitutes various updates in these processes, its performance in simulating energy, water and carbon cycles over the Contiguous United States (CONUS) at scales which land surface changes and hydrometeorological and hydroclimatological applications are more locally relevant is yet to be assessed. In this study, we conducted three simulations at 0.125? during 1979-2018 over the CONUS using different configurations of CLM, namely CLM5-biogeochemistry (CLM5BGC), CLM4.5BGC, and CLM5-satellite phenology (CLM5SP). We validated and compared their simulations against multiple remote-sensed and in-situ datasets. Overall, the parametric and structural updates (e.g., carbon cost for nitrogen uptake, variable soil thickness, dry surface layer) in CLM5 improve its ability in capturing terrestrial biogeochemical dynamics. The low evapotranspiration in CLM5BGC is associated with biases in simulating vegetation phenological characteristics rather than soil water limitations. The mismatch between CLM5BGC-simulated peak leaf area index and reference data can be attributed to CLM5BGC's inability in simulating phenology of trees and grasses. The differences between CLM-simulated irrigation and reference estimates can be attributed to differences between processes represented in models and in reality, and uncertainties in input and validation datasets. Evaluation against observations at small catchments suggest that hydrologic parameters needed to be calibrated to improve simulations of runoff, especially subsurface runoff. Additional efforts are needed to incorporate spatially-distributed plant phenology and physiology parameters and regional-specific agricultural management practices (e.g., planting, harvest).

Cheng, Yanyan↗
DOE OSTI2021

Dissolution of periodic arrays of grains: Upscaling of pore-scale simulations with fast reactions

Upscaling methods are frequently used to derive transport equations at the macroscopic scale from more fundamental equations formulated at the pore scale. These methods typically give a suitable structure for the macroscopic equations and can also provide explicit expressions for the constitutive parameters, such as permeability and dispersion coefficients. Introducing chemical reactions complicates upscaling in at least two important ways. First, the interplay between chemical reactions and transport processes introduces a new length scale, which can be much smaller than the convective or dispersive length scales. A small reactive length scale breaks one of the key assumptions in upscaling; that there is a significant separation in length between the pore-scale and macro-scale processes. The second complication is that if reactions take place at mineral surfaces (dissolution or precipitation) then the pore space itself is evolving in time. In this paper we suggest ways in which these difficulties can be approached, based on analysis of pore-scale simulation data. First, we noticed that the concentration field in successive unit cells has an almost identical spatial variation, with a single scaling factor for each unit cell that is proportional to the incoming reactant flux. Using pore-scale simulations to determine the mass transfer coefficient in a few unit cells, we can calculate the concentration field in the whole domain, even when dissolution is entirely transport limited. Second, we have noticed a time-dependent mapping of the grain shapes from different unit cells. From these observations, we can deduce constitutive relations where the only time-varying parameter is the porosity. Here, we show that a model based on these ideas can quantitatively account for the pore-scale simulation data.

58 GEOSCIENCES↗
DOE OSTI2024

Machine-learning-assisted deciphering of microstructural effects on ionic transport in composite materials: A case study of Li 7 La 3 Zr 2 O 12 -LiCoO 2

The effective diffusivity of ionic species in multiphase materials is critical for the design and function of composite materials for electrochemical energy storage. In practice, effective diffusivity depends sensitively not only on the intrinsic diffusivities of constituting materials but also on their topological arrangement; nevertheless, these coupled contributions are oversimplified in most analytical models. Here, we combine atomistically informed mesoscale modeling and machine learning (ML) analysis to unravel how such features affect effective diffusivity in two-phase composites. Using the Li 7 La 3 Zr 2 O 12 -LiCoO 2 composite solid-state battery cathode as a model system, we compute effective diffusivity for 600 distinct dense polycrystalline microstructures with different topological configurations of grains, grain boundaries, and heterointerfaces. We verify that in addition to atomic-scale variabilities, microstructural feature diversity can significantly impact effective transport properties. Across the ensemble of test microstructures, this often results in bimodal distributions of effective diffusivity that encompass two qualitatively distinct operating mechanisms, which we identify via flux analysis. An ML approach reveals that the most critical determining factors for effective diffusivity are the connectivity of bulk phases and their heterointerfaces. The role of ionic mobility at the heterointerfaces is also discussed. These insights highlight the combined importance of microstructure and interface engineering in tuning the transport properties of ionic species in composite materials. In conclusion, our framework can also be extended for understanding generic microstructure-property relationships in other complex multiphase materials.

25 ENERGY STORAGE↗
DOE OSTI2023

pvlib iotools—Open-source Python functions for seamless access to solar irradiance data

Access to accurate solar resource data is critical for numerous applications, including estimating the yield of solar energy systems, developing radiation models, and validating irradiance datasets. However, lack of standardization in data formats and access interfaces across providers constitutes a major barrier to entry for new users. pvlib python’s iotools subpackage aims to solve this issue by providing standardized Python functions for reading local files and retrieving data from external providers. All functions follow a uniform pattern and return convenient data outputs, allowing users to seamlessly switch between data providers and explore alternative datasets. The pvlib package is community-developed on GitHub: https://github.com/pvlib/pvlib-python. As of pvlib python version 0.9.5, the iotools subpackage supports 12 different datasets, including ground measurement, reanalysis, and satellite-derived irradiance data. The supported ground measurement networks include the Baseline Surface Radiation Network (BSRN), NREL MIDC, SRML, SOLRAD, SURFRAD, and the US Climate Reference Network (CRN). Additionally, satellite-derived and reanalysis irradiance data from the following sources are supported: PVGIS (SARAH & ERA5), NSRDB PSM3, and CAMS Radiation Service (including McClear clear-sky irradiance).

14 SOLAR ENERGY↗
DOE OSTI2022

Effects of laser-annealing on fixed-frequency superconducting qubits

As superconducting quantum processors increase in complexity, techniques to overcome constraints on frequency crowding are needed. The recently developed method of laser-annealing provides an effective post-fabrication method to adjust the frequency of superconducting qubits. In this work, we present an automated laser-annealing apparatus based on conventional microscopy components and demonstrate preservation of highly coherent transmons. In addition, we perform noise spectroscopy to investigate the change in defect features, in particular, two-level system defects, after laser-annealing. Finally, we present a local heating model as well as demonstrate aging stability for laser-annealing on the wafer scale. Our work constitutes an important step toward both understanding the underlying physical mechanism and scaling up laser-annealing of superconducting qubits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Galactomannan utilization by Cellvibrio japonicus relies on a single essential α‐galactosidase encoded by the aga27A gene

Plant mannans are a component of lignocellulose that can have diverse compositions in terms of its backbone and side-chain substitutions. Consequently, the degradation of mannan substrates requires a cadre of enzymes for complete reduction to substituent monosaccharides that can include mannose, galactose, and/or glucose. One bacterium that possesses this suite of enzymes is the Gram-negative saprophyte Cellvibrio japonicus, which has 10 predicted mannanases from the Glycoside Hydrolase (GH) families 5, 26, and 27. Here we describe a systems biology approach to identify and characterize the essential mannan-degrading components in this bacterium. The transcriptomic analysis uncovered significant changes in gene expression for most mannanases, as well as many genes that encode carbohydrate active enzymes (CAZymes) when mannan was actively being degraded. A comprehensive mutational analysis characterized 54 CAZyme-encoding genes in the context of mannan utilization. Growth analysis of the mutant strains found that the man26C, aga27A, and man5D genes, which encode a mannobiohydrolase, α-galactosidase, and mannosidase, respectively, were important for the deconstruction of galactomannan, with Aga27A being essential. Our updated model of mannan degradation in C. japonicus proposes that the removal of galactose sidechains from substituted mannans constitutes a crucial step for the complete degradation of this hemicellulose.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Design, tuning, and blackbox optimization of laser systems

Chirped pulse amplification (CPA) and subsequent nonlinear optical (NLO) systems constitute the backbone of myriad advancements in semiconductor manufacturing, communications, biology, defense, and beyond. Accurately and efficiently modeling CPA+NLO-based laser systems is challenging because of the complex coupled processes and diverse simulation frameworks. Our modular start-to-end model unlocks the potential for exciting new optimization and inverse design approaches reliant on data-driven machine learning methods, providing a means to create tailored CPA+NLO systems unattainable with current models. To demonstrate this new, to our knowledge, technical capability, we present a study on the LCLS-II photo-injector laser, representative of a high-power and spectro-temporally non-trivial CPA+NLO system.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

The Origin of Matter (Final Scientific Report)

The research of this project was concerned with the matter content of the universe. The first topic we studied was the amount of CP violation in kaon decays into two pions; with a new analysis we decreased the theoretical uncertainty by more than an order of magnitude. This constitutes a milestone in our understanding of the matter-antimatter asymmetry among the known standard-model particles. Furthermore, we studied the energy scale dependence of the self-interactions of spinless dark matter for the first time. This allows to study the consistency of theoretical models for dark matter, and is an important ingredient in the interpretation of various dark matter experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A Lamination Model for Pressure-Assisted Sintering of Multilayered Porous Structures

This work describes a lamination model for pressure-assisted sintering of thin, multilayered, and porous structures based on the linear viscous constitutive theory of sintering and the classical laminated plate theory of continuum mechanics. A constant out-of-plane normal stress is assumed in the constitutive relation. The lamination relations between the force/moment resultants and the strain/curvature rates are presented. Numerical simulations were performed for a symmetric tri-layer laminate consisting of a 10% gadolinia doped ceria (Ce 0.9 Gd 0.1 O 1.95-δ ) composite structure, where porous layers were adhered to the top and bottom of a denser layer under uniaxially-applied pressures and the sinter forging conditions. The numerical results show that, compared with free sintering, the applied pressure can significantly reduce the sintering time required to achieve given layer thicknesses and porosities. Unlike free sintering, which results in a monotonic decrease of the laminate in-plane dimension, pressure-assisted sintering may produce an in-plane dimension increase or decrease, depending on the applied pressure and sintering time. Finally, the individual layers in the laminate exhibit different stress characteristics under pressure-assisted sintering.

42 ENGINEERING↗