Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “embedded boundary”

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 91 records · Page 5

Materials Genomics Search for Possible Helium‐Absorbing Nano‐Phases in Fusion Structural Materials

Abstract Civilian fusion demands structural materials that can withstand the harsh environments imposed inside fusion plasma reactors. The structural materials often transmute under 14.1 MeV fast neutrons, producing helium (He), which embrittles the grain boundary (GB) network. Here, it is shown that neutron‐friendly and mechanically strong nano‐phases with atomic‐scale free volume can have low He‐embedding energy and >10 at.% He‐absorbing capacity, and can be especially advantageous for soaking up He on top of resisting radiation damage and creep, provided they have thermodynamic compatibility with the matrix phase, satisfactory equilibrium wetting angle, as well as a high enough melting point. The preliminary experimental demonstration proves that is a good ab initio predictor of He shielding potency in nano‐heterophase materials, and thus, is used as a key feature for computational screening. In this context, a list of viable compounds expected to be good He‐absorbing nano‐phases is presented, taking into account , the neutron absorption and activation cross‐sections, the elastic moduli, melting temperature, the thermodynamic compatibility, and the equilbrium wetting angle of the nano‐phases with the Fe matrix as an example.

36 MATERIALS SCIENCE↗

Soft, malleable double diamond twin

A twin boundary (TB) is a common low energy planar defect in crystals including those with the atomic diamond structure (C, Si, Ge, etc.). We study twins in a self-assembled soft matter block copolymer (BCP) supramolecular crystal having the double diamond (DD) structure, consisting of two translationally shifted, interpenetrating diamond networks of the minority polydimethyl siloxane block embedded in a polystyrene block matrix. The coherent, low energy, mirror-symmetric double tubular network twin has one minority block network with its nodes offset from the (222) TB plane, while nodes of the second network lie in the plane of the boundary. The offset network, although at a scale about a factor of 10 3 larger, has precisely the same geometry and symmetry as a (111) twin in atomic single diamond where the tetrahedral units spanning the TB retain nearly the same strut (bond) lengths and strut (bond) angles as in the normal unit cell. In DD, the second network undergoes a dramatic restructuring—the tetrahedral nodes transform into two new types of mirror-symmetric nodes (pentahedral and trihedral) which alternate and link to form a hexagonal mesh in the plane of the TB. The collective reorganization of the supramolecular packing highlights the hierarchical structure of ordered BCP phases and emphasizes the remarkable malleability of soft matter.

36 MATERIALS SCIENCE↗

Dynamic Mode Decomposition of Random Pressure Fields over Bluff Bodies

Fluctuating surface pressures on a bluff body exposed to a boundary layer flow generally are characterized as a spatiotemporally varying random field. In this paper, a dynamic mode decomposition (DMD) was applied to extract dominant features embedded in these random pressure fields. Utilizing an unsupervised machine learning algorithm, spatial modes and their temporal variations were grouped into different clusters at scales, e.g., macro, meso, and micro. A proper orthogonal decomposition (POD) of the experimental data was carried out to observe commonalities and distinctive perspectives each decomposition offers. Here, a comprehensive examination of the DMD/POD for their convergence criteria, data sufficiency, and modal components analysis was conducted. The physical interpretation of the spatiotemporal pressure field based on these decomposition schemes was discussed. At different scales, the DMD modes can capture the evolution of aerodynamic features, e.g., convection of vortices (or vortex tubes) and other structures. The distribution of energy among these three broad scales also reflects an energy cascade in pressure fluctuations akin to turbulence.

97 MATHEMATICS AND COMPUTING↗

A Perspective on Quantum Computing Applications in Quantum Chemistry Using 25-100 Logical Qubits

The intersection of quantum computing and quantum chemistry represents a promising frontier for achieving quantum utility in domains of both scientific and societal relevance. Owing to the exponential growth of classical resource requirements for simulating quantum systems, quantum chemistry has long been recognized as a natural candidate for quantum computation. This perspective focuses on identifying scientifically meaningful use cases where early fault-tolerant quantum computers, which are considered to be equipped with approximately 25-100 logical qubits, could deliver tangible impact. While recent advances in classical computing have pushed the boundaries of tractable simulations to unprecedented scales, this logical-qubit regime represents the first window where quantum devices can pursue qualitatively distinct strategies, such as polynomial-scaling phase estimation, direct simulation of quantum dynamics, and active-space embedding, that remain challenging for classical solvers, such as multireference charge-transfer and conical-intersection states central to photochemistry and materials design. We highlight near-term opportunities in algorithm and software design, discuss representative chemical problems suited for quantum acceleration, and propose strategic roadmaps and collaborative pathways for advancing practical quantum utility in quantum chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Revealing the corrosion mechanism of an Al 0.1 CoCrFeNi high entropy alloy in high temperature carbon dioxide environment

The corrosion behavior of an Al 0.1 CoCrFeNi HEA in high temperature CO 2 at 850 °C for different exposure periods was investigated using different characterization methods. Further, weight gain increased with extended exposure periods, yet it remained lower than most conventional alloys. A double-layer oxide film consisting of Cr 2 O 3 and Al 2 O 3 formed after 100 and 500 h of exposure, while a single layer of Cr 2 O 3 developed after 1000 h of exposure. CoCrFeNi particles with an FCC structure were embedded in the oxide films irrespective of the exposure periods. Moreover, a transition zone containing numerous striped-like Al 2 O 3 oxides along the grain boundaries was developed underneath the oxide films.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Isothermal and Cyclic Oxidation of Haynes 282 Processed by Electron Beam Melting (EBM) and Laser Powder Bed Fusion (LPBF) in Dry Air at 800 and 950°C

Here, the isothermal and cyclic oxidation behavior of Haynes 282 (H282) processed by laser powder bed fusion and electron beam melting was compared to its cast counterpart during exposures in air at 800 and 950°C. The specific microstructure of the AM alloys compared to coarse-grain cast 282 resulted in differences in oxidation rates, internal oxidation and spallation behavior. At 800°C, faster Ti diffusion to the surface due to the smaller grain size of the AM alloys led to Ti doping of the Cr 2 O 3 scale and faster oxidation rates. The higher density of grain boundaries also resulted in more pronounced internal oxidation. At 950°C, a duplex Cr 2 O 3 scale was observed with a dense inner layer and porous cracked outer layer with embedded TiO 2 particles. The impact of the external and internal oxide compositions and oxidation-induced elemental depletions on the oxidation lifetime, creep properties and spallation behavior is discussed.

36 MATERIALS SCIENCE↗

Predicting interface structure using the minima hopping method

Here, we adapt the minima hopping method (MHM) to the problem of interfacial structure prediction and apply it to study a canonical problem, the tilt grain boundaries in SrTiO 3 . Our method employs a hybrid approach by first exploring the potential energy surface (PES) of different grain boundary samplings with an empirical force field, among which the fifteen candidates with lower energies are then refined using ab initio density functional theory (DFT) calculations. During the exploratory stage, we bias the search using a local order parameter to primarily sample various reconstructions in the vicinity of the interface, while preserving the crystallinity of the bulk regions. We further enhance the search by incorporating initial structures with rigid body displacements to account for translational variations between bulk phases, enabling the MHM to effectively generate both stoichiometric and nonstoichiometric SrTiO 3 Σ⁢3(111)[110] and Σ⁢3(112)[110] grain boundaries. From an algorithmic standpoint, MHM outperforms earlier studies based on genetic algorithms (GA) by identifying more stable interfacial structures of several SrTiO 3 grain boundaries. The performance of the present implementation of the MHM approach is primarily limited by exploring an approximate description of the PES with a rather simple Buckingham potential. This limitation leads to variations in performance when compared to approaches utilizing more advanced surrogate PES models, such as direct DFT-PES sampling or GA with the embedded atom method (EAM). Despite the present limitations, the MHM approach is able to yield interfacial structures with comparable or lower interfacial energies in specific cases, such as Σ⁢3(111)[110] Γ=1, ±0.5 and Σ⁢3(112)[110] Γ= ±1, −2, underscoring the robustness of the MHM approach even with a simple approximation of the DFT PES. The MHM interfacial structure prediction method thus offers an efficient approach to understanding the grain boundaries and heterointerfaces at the atomic scale, providing an important prerequisite for effective materials design.

density functional theory↗

Templates for Risk Informed Assurance with Curvature Embeddings (TRACE)

We investigate recovery of geometric structure from networks embedded in manifolds with spatially varying curvature, extending the constant-curvature framework of Lubold et al. (2023). Our work supports cascade risk assessment in critical infrastructure through the Templates for Risk-informed Assurance with Curvature Embeddings (TRACE) framework. Simulations on a bi-modal Gaussian surface show that constant-curvature methods yield weighted averages shaped by clique patterns, while hierarchical clustering identifies distinct regimes. Localized estimation, however, reveals boundary contamination in transitional regions. To address heterogeneity, we develop distance metrics for graphs with edge and node features, proving their metric validity, and validate them via deterministic graph generation from canonical tilings. We further propose a diffusion-based anomaly detection approach that treats networks as glued manifolds, using curvature discontinuities to detect structural anomalies. Employing the carré-du-champ operator and scalar curvature, we achieve robust anomaly discrimination, demonstrated on the Singapore Water Treatment (SWaT) dataset with joint network-traffic and sensor features. Integration with TRACE reveals how curvature shapes cascade dynamics: positive curvature impedes, while negative curvature accelerates propagation. This geometric perspective provides interpretable risk metrics and visualization tools for critical infrastructure managers. While full validation remains ongoing, our contributions establish a rigorous foundation for geometric analysis of network resilience and cascade vulnerability.

97 MATHEMATICS AND COMPUTING↗

TRUST Sensors in Environments: Thermocouples (SE-TC), Release FY25

The Delivery Environments Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) project is a broad project intended to analyze simplified problems experimentally and with modeling and simulation. The purpose of analyzing these simplified problems is to extend solution methods to more complex problems, as well as understand deficiencies and gaps in knowledge of methods currently used in more complex analyses. The TRUST project encompasses several smaller testbeds intended to isolate individual phenomena. The testbed under consideration in this report is the Sensors in Environments: Thermocouples testbed. In previous years, the purpose of this testbed was to quantify uncertainty of thermocouple sensors. To accomplish this, an aluminum plate was placed in a thermal chamber and subject to various types of thermal loading. Thermocouples were placed in various locations on the aluminum plate in various configurations (e.g., embedded in the plate, placed under Kapton tape), and an effort was made to quantify uncertainty in these measurements. Finite element simulations were performed to investigate how sensitive these measurements were to parameters such as the boundary conditions on the plate and material properties. However, a fundamental source of uncertainty in this analysis was the convective heat transfer from the plate. Convective heat transfer is a complex physical phenomenon comprised of a number of interacting sub-processes and is difficult to predict accurately a priori. As such, the main purpose of this testbed in FY25 was to better understand, both experimentally and numerically, the convective heat transfer from the plate. This is a highly applicable problem to several more complex problems, as convective heat transfer occurs in nearly all problems where a body is moving through air. Numerically, this required a two-step approach. First, the air flow in the thermal chamber was in vestigated using computational fluid dynamics. The commercial solver Fluent was used to perform these simulations. From these simulations, a heat transfer coefficient over the surface of the plate was calculated. This heat transfer was then used as boundary conditions for finite element heat transfer simulations within the plate, which were performed using Abaqus. Significant effort was devoted to automating the handoff between these two solvers. Experimentally, previous thermocouple results in the plate were used to validate the time-dependent thermal profiles produced from Abaqus. Further experimental efforts were performed both to help validate the Fluent simulations and to inform its boundary conditions. For example, hot-wire anemometers were used to measure the velocity in the chamber, which would be particularly useful in understanding the chamber inlet velocity. Thermocouple measurements were also taken in the chamber, instead of only on the plate, to serve as validation evidence for the Fluent simulations. Numerical results showed that the Fluent to Abaqus workflow matched previous plate thermocouple measurements well. This type of handoff is useful for more complex experiments, or those that are not able to be examined in as great of detail as this testbed, as it was performed without any experimental input. Experimental results, however, were more mixed. The anemometers proved unreliable, with inconsistent measurements across all anemometers, even at locations that were nearly identical. On the other hand, the thermocouples provided a relatively rich view of the temperature field in the chamber.

42 ENGINEERING↗

Impact of Lateral Flow on Surface Water and Energy Budgets over the Southern Great Plains – A Modeling Study

As we see the horizontal grid spacing decrease, treatment of hydrologic processes in land surface models, such as the lateral flow of surface and subsurface flow, need to be explicitly represented. Unlike previous studies that mainly focused on the mountainous regions, in this study the offline WRF-Hydro model is employed to study the impact of lateral flow on soil moisture and energy fluxes over the relatively flat southern Great Plains (SGP). The vast amount of measurements over the SGP provide an unique opportunity to assess the model behavior. In addition, newly developed land surface properties and input forcing are ingested into the model, in an attempt to reduce uncertainties associated with the initial and boundary forcing and help to identify model deficiencies. Our results show that the more realistic inputs (parameters, soil types, forcing) lead to larger underestimation of latent heat flux and dry bias, indicating the existence of model structural uncertainty (embedded errors) in WRF-Hydro that need to be characterized to inform future model development efforts. Including lateral flow processes partly mitigates the model deficiencies in representing hydrologic processes and alleviates the dry bias. In particular, both surface and subsurface lateral flow increase soil moisture mainly over the lower elevations, except that subsurface flow also affects soil moisture over steeper terrains. Additional simulations are performed to assess the effect of routing resolution on model results. When LSM resolution is high, noticeable differences in soil moisture are produced between different routing resolutions especially over steep terrain. Whereas when LSM resolution is coarse, differences between routing resolutions become negligible, especially over flat terrain.

54 ENVIRONMENTAL SCIENCES↗

Cell invasion during competitive growth of polycrystalline solidification patterns

Abstract Spatially extended cellular and dendritic array structures forming during solidification processes such as casting, welding, or additive manufacturing are generally polycrystalline. Both the array structure within each grain and the larger scale grain structure determine the performance of many structural alloys. How those two structures coevolve during solidification remains poorly understood. By in situ observations of microgravity alloy solidification experiments onboard the International Space Station, we have discovered that individual cells from one grain can unexpectedly invade a nearby grain of different misorientation, either as a solitary cell or as rows of cells. This invasion process causes grains to interpenetrate each other and hence grain boundaries to adopt highly convoluted shapes. Those observations are reproduced by phase-field simulations further demonstrating that invasion occurs for a wide range of misorientations. Those results fundamentally change the traditional conceptualization of grains as distinct regions embedded in three-dimensional space.

36 MATERIALS SCIENCE↗

Multiscale analysis in solids with unseparated scales: fine-scale recovery, error estimation, and coarse-scale adaptivity

There are several engineering applications in which the assumptions of homogenization and scale separation may be violated, in particular, for metallic structures constructed through additive manufacturing. Instead of resorting to direct numerical simulation of the macroscale system with an embedded fine scale, an alternative approach is to use an approximate macroscale constitutive model, but then estimate the model-form error using a posteriori error estimation techniques and subsequently adapt the macroscale model to reduce the error for a given boundary value problem and quantity of interest. Here, we investigate this approach to multiscale analysis in solids with unseparated scales using the example of an additively manufactured metallic structure consisting of a polycrystalline microstructure that is neither periodic nor statistically homogeneous. As a first step to the general nonlinear case, we focus here on linear elasticity in which each grain within the polycrystal is linear elastic but anisotropic.

42 ENGINEERING↗

Embedded Ferroelectric Nanoclusters can drive Polarization Reversal in a Non-Ferroelectric Polar Film via the Proximity Effect

Heterogeneous nucleation from defects dominates the electric field required for polarization switching of ferroelectrics. Here, we consider the switching of a nominally non-switchable polar thin film of AlN due to the proximity effect arising from embedded ferroelectric nanoclusters of Al 1-x Sc x N. Using a Landau-Ginzburg-Devonshire thermodynamic approach and finite element modeling, we study the influence of nanocluster shape on polarization switching and domain nucleation emerging in AlN. The ferroelectric nanocluster boundary is modeled as a thin layer transitioning from Al 1-x Sc x N to AlN. We analyze the conditions under which polarization switching in the AlN film occurs at coercive fields significantly lower than its dielectric breakdown field. In the presence of spike-like Al 1-x Sc x N nanoclusters, the proximity effect enables switching of the spontaneous polarization in AlN and significantly reduces the corresponding coercive field. The internal field, which is depolarizing inside the AlN (due to its larger spontaneous polarization) and polarizing within the ferroelectric Al 1-x Sc x N nanoclusters (due to its smaller spontaneous polarization), lowers the potential barrier in the clusters and nucleates nanodomains at the Al 1-x Sc x N-AlN interface, forming localized regions of reversed polarization. Proximity effect can thus provide a pathway towards "thawing" previously "frozen" ferroelectrics through engineered nucleation for memory, actuation and optical technologies.

36 MATERIALS SCIENCE↗

Characterization of mechanical discontinuities based on data-driven classification of compressional-wave travel times

Wave propagation and diffusive transport phenomena are influenced by the mechanical discontinuities in material. Here this study shows that certain bulk properties of the network of low-velocity mechanical discontinuities (e.g. air-filled cracks) in a material can be characterized by processing compressional-wave travel times using traditional data-driven classification techniques. To that end, we perform three tasks in chronological order: (1) use the discrete fracture network (DFN) method to create two-dimensional (2D) numerical models of crack-bearing material embedded with various types of low-velocity mechanical discontinuities, (2) use the fast marching method (FMM) to simulate the propagation of the wave/diffusion front from a single source through the 2D crack-bearing material to multiple receivers placed on the boundary of the material, and (3) train 9 data-driven classifiers to characterize the crack-bearing materials (i.e. bulk properties of the network of mechanical discontinuities in the crack-bearing material) by learning from the simulations of travel times detected by multiple receivers placed around the crack-bearing material. The classifiers identified the orientation, spatial distribution, and dispersion of the low-velocity mechanical discontinuities. Voting classifier performs the best among the 9 classifiers. For the characterization of bulk dispersion and distribution of discontinuities, the sensors located on the adjacent boundaries are more important; whereas for the characterization of bulk orientation of discontinuities, the sensors located on the opposite side are more important.

42 ENGINEERING↗

Creep anisotropy modeling and uncertainty quantification of an additively manufactured Ni-based superalloy

The advantages offered by additive manufacturing over traditional processes has driven a great deal of industrial and academic interest in recent years. However, the process is relatively new and requires additional investigation to become sufficiently mature for wide scale industrial adoption. Electron beam melting powder bed fusion is one technology that has shown promise for fabricating high temperature resistant materials such as nickel based superalloys. The resulting microstructures typically exhibit a strong fiber texture in the build direction giving rise to anisotropic time-dependent deformation behavior. In order to accelerate the qualification of these materials for industrial adoption accurate numerical models are needed for simulating their behavior. In this work a crystal plasticity model including non-Schmid effects is presented for capturing creep anisotropy observed in additively manufactured IN738LC. The model is calibrated via a probabilistic framework where model parameters are treated as random variables. An iterative sequential design strategy is utilized to efficiently identify the probability density of the unknown model parameters. As a case study the model is utilized to investigate the behavior of randomly oriented equiaxed grain clusters sometimes observed embedded in the additively manufactured columnar structure. A synthetic realization is simulated and uncertainty is propagated through to the full-field response. Results indicate that these features are the source of significant creep relaxation and strain accumulation which partially explains observed grain boundary decohesion at these locations.

36 MATERIALS SCIENCE↗

Single- and two-particle finite size effects in interacting lattice systems

Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with nonlocal correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles.

36 MATERIALS SCIENCE↗

Design of a Combined Proportional Integral Derivative Controller to Regulate the Temperature Inside a High-Temperature Tubular Solar Reactor

Solar fuels are proven to be promising candidates for thermochemical energy storage. However, the transient nature of solar radiation is an obstacle to maintaining a stable operational temperature inside a solar reactor. To overcome this challenge, the temperature of a solar reactor can be regulated by controlling the incoming solar radiation or the feedstock flowrate inside the reactor. In this work, a combined proportional integral derivative (PID) controller is implemented to regulate the temperature inside a high-temperature tubular solar reactor with counter-current flowing gas/particles. The control model incorporates two control systems to regulate incoming solar radiation and gas flow simultaneously. The design of the controller is based on a reduced-order numerical model of a high-temperature tubular solar reactor that is vertically oriented with an upward gas flow and downward particle flow. The reactor receives heat circumferentially through its wall over a finite segment of its length. Formulation of the heat transfer model is presented by applying the energy balance for the reactor tube and considering heat and mass transfer inside. A set of governing differential equations are solved numerically by using the finite volume method to obtain reactor wall, particles, and gas temperatures along the reactor length with various boundary conditions. Simulation results are used to tune the PID controller parameters by utilizing the Ziegler–Nichols tuning method. Both the simulation results and the controller performance are visualized on the LABVIEW platform. The controller is challenged to track different temperature setpoints with different scenarios of transient solar radiation. Here, the performance of the PID controller was compared to experimental results obtained from an industrial PID controller embedded in a 7 kW electric furnace. Results show that the combined PID controller is successful in maintaining a stable temperature inside the reactor by regulating the incoming solar radiation and the flowrate via small steady-state error and reasonable settling time and overshoot.

14 SOLAR ENERGY↗

Temperature effects of ion irradiation on the nanostructural features in ductile-phase-toughened tungsten composites

Ductile-phase toughened tungsten (DPT W) composites have emerged as promising candidates for load-bearing components behind the plasma-facing tungsten armor in fusion reactors due to their enhanced thermomechanical properties. This study focuses on a composite consisting of W particles embedded in a ductile NiFeW solution matrix, hot-rolled to a thickness reduction of 87% (87R DPT W). Sequential irradiations with Ni2+ and He+ ions were performed to identical doses and helium concentrations at room temperature (RT) and 1273 K. Irradiation at RT produced no discernible nanostructural features due to the immobility of mono-vacancies, whereas cavity formation was observed at 973 K. At 1273 K, the W phase exhibited larger cavities, reduced cavity number density, and lower volumetric swelling compared to 973 K. Notably, nanosized NiFeW precipitates formed within the W phase at 1273 K, a phenomenon absent at 973 K. A new phase of cubic (NiFe)6W6C was also observed at the interphase boundary. In contrast, the NiFeW matrix showed no nanostructural changes at 1273 K, likely due to cavity dissociation. Separate irradiations at 1273 K indicated that Ni2+ ions induced precipitate formation in the W phase, while He+ ions exclusively caused cavity formation. The microstructure of 87R DPT W irradiated at RT and subsequently annealed at 1273 K closely resembled that of material irradiated directly at 1273 K. Like oxide-dispersion-strengthened steels, the observed nanoparticle-embedded W can inhibit dislocation propagation, potentially delaying the ductile-to-brittle transition temperature. These findings highlight the potential of NiFeW nanoparticle-reinforced W composites as irradiation-resistant materials for fusion reactors.

Jiang, Weilin (ORCID:0000000183028313)↗