Search NASA⌕ Search

SEARCH · Search NASA

Results for “Embedded boundaries”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Unraveling the dislocation–precipitate interactions in high-entropy alloys

The precipitates play a significant role in not only enhancing the strength, but also maintaining the high toughness in alloys. However, the interactions of the nanoscale precipitates with dislocations in the high entropy alloys (HEAs) are difficult to observe directly by in-situ TEM experiments due to the limits of the resolution and time. Here, using atomic simulations we report the synergistic strengthening of the coherent precipitate and atomic-scale lattice distortion in the HEAs at cryogenic/elevated temperatures. The effects of temperature, chemical disorder, precipitate spacing, precipitate size, elemental segregation, and dislocation-cutting number on the critical stress for the dislocation to overcome a row of precipitates are studied. A random stacking fault energy landscape along the slip plane, the lattice distortion at different temperatures, and the interface/surface energy at various crystallographic orientations are obtained. Compared with the traditional metals and alloys, HEAs have the severe atomic-scale lattice distortions to generate the local high tensile/compressive stress fields. This complex stress causes the dislocation line to bend, and thus improves the dislocation slip resistance, resulting in the strong solid-solution strengthening. The stacking fault strengthening induced by the obvious difference of the stacking fault energies between the HEA matrix and precipitate (within the inner of the HEA matrix), and the formation of the antiphase domain boundary contribute to the high strength. The precipitate embedded by the solute atoms produces the strong lattice distortion to enhance the dislocation slip resistance at high temperatures. Hence, the current results provide the mechanistic insight into the phenomenon that the coherent precipitate combined with the severe atomic-scale lattice distortion can enhance the strength at cryogenic/elevated temperatures to further broaden the scope of applications of advanced HEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A high-order Shifted Interface Method for Lagrangian shock hydrodynamics

Here, we present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that is, it uses the precise location of the material interface. The interface is represented by the zero level-set of a continuous high-order finite element function that moves with the material velocity. This strategy allows to evolve curved material interfaces inside curved elements. By reformulating the original interface problem over a surrogate (approximate) interface, located in proximity of the true interface, the SIM avoids cut cells and the associated problematic issues regarding implementation, numerical stability, and matrix conditioning. Accuracy is maintained by modifying the original interface conditions using Taylor expansions. We demonstrate the performance of the proposed algorithms on established numerical benchmarks in one, two and three dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Seismic Performance of Deeply Embedded Conceptual Advanced Reactors: Three-Dimensional Nonlinear Soil-Structure Interaction Analyses

The U.S. Nuclear Regulatory Commission’s functional containment concept provides advanced nuclear power plant designers with more flexibility in terms of the civil/structural design if the appropriate set of barriers for prevention of radioactive material release exist. Some of the conceptual advanced reactor structures, without the traditional pressure boundaries of large containment structures, are proposed to be deeply embedded or buried into soil. This approach is expected to provide (1) lesser seismic demands on the structures and safety-critical structures, (2) eased regulatory efforts and overall design against other external hazards such as aircraft impact, and (3) overall cost savings. Further, one of the important aspects of assessing the technical and economic viability of deeply embedding advanced reactor buildings is to assess the seismic performance with the understanding of effects with material and geometric nonlinearities. This study investigates the seismic response of deeply embedded or buried advanced reactors by conducting three-dimensional nonlinear soil-structure interaction analyses. Although the results indicate that there is a general trend of decreased seismic response with increased embedment depths, the change in the dynamic environment with different embedment depths and the nonlinear environment under high-intensity seismic inputs may result in increased peak response at increased embedment depths.

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

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Suppression of Collisionless Magnetic Reconnection in the High Ion β, Strong Guide Field Limit

Abstract In magnetic reconnection, the ion bulk outflow speed and ion heating have been shown to be set by the available reconnecting magnetic energy, i.e., the energy stored in the reconnecting magnetic field ( B r ). However, recent simulations, observations, and theoretical works have shown that the released magnetic energy is inhibited by upstream ion plasma beta β i —the relative ion thermal pressure normalized to magnetic pressure based on the reconnecting field—for antiparallel magnetic field configurations. Using kinetic theory and hybrid particle-in-cell simulations, we investigate the effects of β i on guide field reconnection. While previous works have suggested that guide field reconnection is uninfluenced by β i , we demonstrate that the reconnection process is modified and the outflow is reduced for sufficiently large β i > ( B r 2 + B g 2 ) / B r 2 . We develop a theoretical framework that shows that this reduction is consistent with an enhanced exhaust pressure gradient, which reduces the outflow speed as v out ∝ 1 / β i . These results apply to systems in which guide field reconnection is embedded in hot plasmas, such as reconnection at the boundary of eddies in fully developed turbulence like the solar wind or the magnetosheath as well as downstream of shocks such as the heliosheath or the mergers of galaxy clusters.

Giai, Carlos A.↗
DOE OSTI2026

Cooperative effect of local active stresses on the macroscopic contractility of elastic fiber networks

The collective action of actively contractile units embedded in elastic biopolymer networks plays a crucial role in regulating the network's macroscopic mechanical response. Here, in this study, we investigate how the macroscopic boundary stress in model elastic fiber networks depends on the number and nature of embedded contractile units, each exerting an isotropic force dipole, as well as on the bending stiffness of fibers. We find that the macroscopic stress increases nonlinearly with the number of dipoles due to mutual stiffening of initially soft, bending-dominated networks. Using effective medium theory, we relate this enhanced contractility to an increase in the effective average network coordination number due to constraints imposed by the force dipoles. By comparing three distinct force dipole models that differ in their local structures, we demonstrate that the specific manner in which an active unit constrains the network strongly influences the onset and nature of the stiffening transition. Our results highlight that not only the quantity but also the local geometry of force-generating units critically determines the macroscopic mechanical behavior. This framework provides a physical basis for understanding how biological systems—such as molecular motors in the cytoskeleton, or adherent cells in the extracellular matrix—can modulate network-scale nonlinear elastic properties through local tuning of active force-generating units.

Biological and medical sciences↗
DOE OSTI2021

Nonlinear convergence in contact mechanics: Immersed boundary finite volume

In this report we present an immersed boundary finite volume (IBM) method for simulating quasistatic contact mechanics of linearly elastic domains at small strains. In IBM, all external boundaries and internal contacts of an object are represented by embedded surfaces inside a Cartesian mesh, which need not conform to the grid lines. The contact constraints consist of the non-penetrability condition and Coulomb’s friction law, which are discretized using special interpolation stencils and enforced via penalty parameters. The resulting nonlinear system depends on displacement unknowns only. To solve it, we use the Newton method but find that it diverges frequently. To understand the divergence pattern, we analyze a simplified 2-cell problem and show that the global convergence of Newton cannot be ensured for any choice of penalty parameters. We thus propose a modified Newton solver, which guarantees convergence for the 2-cell problem and is numerically verified to converge for all the challenging simulations considered herein. While both 1 st - and 2 nd -order variants of IBM, in displacement unknowns, are proposed, the modified Newton solver applies only to the 1 st -order variant.

42 ENGINEERING↗
DOE OSTI2020

Helical spin structure in iron chains with hybridized boundaries

In this paper, we have compared the magnetic properties of well-controlled ultra-short (≤50 nm) atomic iron (Fe) chains embedded in Fe-phthalocyanine films with those in Fe–hydrogen (H2) phthalocyanine superlattices. Surprisingly, we found that the coercivity of the atomic chains with free boundary conditions is independent of the chain length, whereas the one subject to hybridization of the chain ends exhibits an unexpected length dependence. These findings suggest that ferromagnetism in the free-boundary condition system is caused by an intrinsic indirect exchange. On the other hand, controlled boundary conditions produce a helical spin structure due to an extrinsic indirect exchange, which arises from the interaction between iron atoms at the ends of the chain and the hydrogen in the H 2 phthalocyanine spacer. As a consequence, during magnetic reversal, ultra-short iron chains subject to boundary clamping develop a helical spin structure, leading to increased coercivity. These findings suggest unique insights and ideas for the design of atomic-scale ultra-dense magnetic storage nanodevices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A deep learning upscaling framework: Reactive transport and mineral precipitation in fracture-matrix systems

Pore-scale modeling has limited applicability at large scales due to its high computational cost. One common approach to upscale pore-scale models is the use of effective medium theories, which homogenize small-scale features in a porous structure and characterize the medium by macroscale properties (e.g., permeability) and equations (e.g., Darcy’s law). However, there are classes of physical processes for which effective medium approximations may become inaccurate, e.g., mineral precipitation and clogging during reactive transport. We have developed a deep learning upscaling framework, in which pore-scale modeling is directly employed in macroscale systems, without relying on effective medium approximations. The upscaling framework is first developed for general multiscale systems and then applied to modeling reactive transport with mineral precipitation in the altered layer in fracture-matrix structures. Solute transport from the fractures to the matrix is modeled as a wall boundary condition for the fractures, which, in turn, is predicted by recurrent neural networks using the concentration histories at the fracture-matrix boundary. Specifically, we consider a meter-scale fracture network embedded in sandstones, where the smallest feature is at the micron scale. Here the proposed framework allows us to span five orders of magnitude in length scales by capturing mineral precipitation in the altered layer of the rock matrix at the pore scale across the entire meter-scale fracture network.

42 ENGINEERING↗
DOE OSTI2022

Mapping out the internal space in AdS/BCFT with Wilson loops

We study Wilson loops in string theory realizations of AdS/BCFT and wedge holography. The field theories are based on 3d N = 4 long quiver gauge theories engineered by D3, D5 and NS5 branes, and on BCFTs involving 4d N = 4 SYM coupled to such 3d theories. The holographic duals have geometry AdS 4 × S 2 × S 2 warped over a strip. We identify the holographic representation of antisymmetric Wilson loops associated with individual 3d gauge nodes in terms of probe D5-branes. The expectation values obtained holographically are matched to supersymmetric localization computations. Our results yield an identification of regions in the internal space with individual 3d gauge nodes. Connecting to bottom-up braneworld models, this gives a concrete notion of which parts of the 10d solutions correspond to the end-of-the-world brane and which to the remaining bulk. We also construct supersymmetric Janus on the brane embeddings in AdS 5 ×S 5 , which describe surface defects with boundaries and interfaces in N = 4 SYM.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Liquid Phase Modeling in Porous Media: Adsorption of Methanol and Ethanol in H-MFI in Condensed Water

Zeolites are used in the chemical and separation industries for their exceptional selectivity, adsorption capacity, regenerability, and stability in gas and liquid phase processing. Here, we developed an explicit solvation method for predicting solvent/condensed phase effects on adsorption free energies in microporous media such as zeolites based on the hybrid quantum mechanical/molecular mechanical free energy perturbation (QM/MM-FEP) technique. Our explicit solvation method for zeolite systems, called eSZS, aims to capture site-specific interactions during the adsorption process at the Brønsted acid sites of H-MFI zeolite while still considering the diverse configuration space of the solvent molecules. This strategy is ideal for chemical reactions or adsorbates that interact with the microporous medium in few distinct adsorbate/transition state configurations, i.e., the harmonic or similar approximations are acceptable for the adsorbate/transition state while such approximations break down for the solvent molecules that require extensive configuration space sampling. In this way, our approach effectively overcomes the limitations of implicit solvation models and classical force field methods for describing solvation effects on chemical reactions within porous materials such as zeolites. Specifically, in this study, we investigated various aspects of our hybrid QM/MM approach, including QM cluster size dependencies in a periodic electrostatically embedded cluster model (PEECM), rules for link atoms at the QM/MM boundary, and functional and basis set considerations for converged and reasonably accurate gas and aqueous phase methanol and ethanol adsorption free energy predictions in H-MFI. For gas phase adsorption of methanol and ethanol in H-MFI at a Brønsted acid site in T12 position, we compute adsorption free energies at 298 K of −0.61 and −0.75 eV, respectively, using a PEECM containing 50 Si and 1 Al atom with ωB97x-D/def2-TZVP level of theory. For solvent effect calculations, we sample the aqueous phase using grand canonical Monte Carlo (GCMC) simulations to (1) obtain a mean field of electrostatic interactions in the reaction system and (2) perform a rigorous free energy perturbation calculation. Similar to the experimentally and computationally observed endergonic solvation effects observed for hydrocarbon adsorption on metal surfaces, we also observe that a condensed aqueous environment destabilizes methanol and ethanol at these acid sites in H-MFI at 298 K. Specifically, the computed solvation free energies of adsorption (ΔΔG solv ) for methanol and ethanol are +0.44 and +0.54 eV, respectively. From this study, it is evident that adsorbates (methanol and ethanol) are competing with water for adsorption space inside the H-MFI zeolite, leading to an endergonic solvation effect. Here, we expect that the endergonic, aqueous solvent effect during adsorption in microporous zeolites is highly tunable by changing the pore size and hydrophobicity of the microporous material as this will affect the water density inside the pore structure.

Adsorption↗
DOE OSTI2022

Towards Ab-Initio Simulations of Crystalline Defects at the Exascale Using Spectral Quadrature Density Functional Theory

Defects in crystalline solids play a crucial role in determining properties of materials at the nano, meso- and macroscales, such as the coalescence of vacancies at the nanoscale to form voids and prismatic dislocation loops or diffusion and segregation of solutes to nucleate precipitates, phase transitions in magnetic materials via disorder and doping. First principles Density Functional Theory (DFT) simulations can provide a detailed understanding of these phenomena. However, the number of atoms needed to correctly simulate these systems is often beyond the reach of many widely used DFT codes. The aim of this article is to discuss recent advances in first principles modeling of crystal defects using the spectral quadrature method. The spectral quadrature method is linear scaling with respect to the number of atoms, permits spatial coarse-graining, and is capable of simulating non-periodic systems embedded in a bulk environment, which allows the application of appropriate boundary conditions for simulations of crystalline defects. In this article, we discuss the state-of-the-art in ab-initio modeling of large metallic systems of the order of several thousand atoms that are suitable for utilizing exascale computing resourses.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Differences in chemical evolution between isolated and embedded prestellar cores

Models of prestellar cores often assume that the cores are isolated from their environment – material outside the core boundary plays no role in the subsequent evolution. This is unlikely to be the case in reality, where cores are located within hierarchically substructured molecular clouds. Here, we investigate the dynamical and chemical evolution of prestellar cores, modelled as Bonnor–Ebert spheres, and show that the density of the ambient medium has a large impact on the resulting chemical properties of the cores. Models embedded in high-density, low-temperature surroundings have greatly enhanced abundances of several molecules, such as CO and CS, compared to models with more diffuse surroundings, corresponding to relatively isolated cores. The predicted intensities and profile shapes of molecular lines are also affected. The density of the ambient medium has a stronger effect on the chemical evolution than whether the cores are initially in or out of equilibrium. This suggests that the impact of environment cannot be neglected when modelling chemistry in prestellar cores; the results of these models are highly sensitive to the assumptions made about the core surroundings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Data‐driven variational method for discrepancy modeling: Dynamics with small‐strain nonlinear elasticity and viscoelasticity

Abstract The effective inclusion of a priori knowledge when embedding known data in physics‐based models of dynamical systems can ensure that the reconstructed model respects physical principles, while simultaneously improving the accuracy of the solution in the previously unseen regions of state space. This paper presents a physics‐constrained data‐driven discrepancy modeling method that variationally embeds known data in the modeling framework. The hierarchical structure of the method yields fine scale variational equations that facilitate the derivation of residuals which are comprised of the first‐principles theory and sensor‐based data from the dynamical system. The embedding of the sensor data via residual terms leads to discrepancy‐informed closure models that yield a method which is driven not only by boundary and initial conditions, but also by measurements that are taken at only a few observation points in the target system. Specifically, the data‐embedding term serves as residual‐based least‐squares loss function, thus retaining variational consistency. Another important relation arises from the interpretation of the stabilization tensor as a kernel function, thereby incorporating a priori knowledge of the problem and adding computational intelligence to the modeling framework. Numerical test cases show that when known data is taken into account, the data driven variational (DDV) method can correctly predict the system response in the presence of several types of discrepancies. Specifically, the damped solution and correct energy time histories are recovered by including known data in the undamped situation. Morlet wavelet analyses reveal that the surrogate problem with embedded data recovers the fundamental frequency band of the target system. The enhanced stability and accuracy of the DDV method is manifested via reconstructed displacement and velocity fields that yield time histories of strain and kinetic energies which match the target systems. The proposed DDV method also serves as a procedure for restoring eigenvalues and eigenvectors of a deficient dynamical system when known data is taken into account, as shown in the numerical test cases presented here.

Masud, Arif↗
DOE OSTI2023

Rectified Rotational Dynamics of Mobile Inclusions in Two-Dimensional Active Nematics

We investigate the dynamics of mobile inclusions embedded in 2D active nematics. The interplay between the inclusion shape, boundary-induced nematic order, and autonomous flows powers the inclusion motion. Disks and achiral gears exhibit unbiased rotational motion, but with distinct dynamics. In comparison, chiral gear-shaped inclusions exhibit long-term rectified rotation, which is correlated with dynamics and polarization of nearby +1/2 topological defects. Furthermore, the chirality of defect polarities and the active nematic texture around the inclusion correlate with the inclusion’s instantaneous rotation rate. Inclusions provide a promising tool for probing the rheological properties of active nematics and extracting ordered motion from their inherently chaotic motion.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Reverse Osmosis Simulation Data

This dataset consists of computational fluid dynamics (CFD) output for various spacer configurations in a feed-water channel in reverse osmosis (RO) applications. Feed-water channels transport brine solution to the RO membrane surfaces. The spacers embedded in the channels help improve membrane performance by disrupting the concentration boundary layer growth on membrane surfaces. Refer to the "Related Work" resource below for more details. This dataset considers a feed-water channel of length 150mm. The inlet brine velocity and concentration are fixed at 0.1m/s and 100kg/m3 respectively. The diameter of the cylindrical spacers is fixed as 0.3mm and six varying inter-spacer distances of 0.75mm, 1mm, 1.5mm, 2mm, 2.5mm, and 3mm are simulated. The dataset comprising the steady, spatial fields of solute concentration, velocity, and density near each spacer is placed in the folder corresponding to the spacer configuration considered. We run two sets of CFD simulations and include the outputs from both sets for each configuration: (1) with a coarser mesh, producing low-resolution (LR) data of spatial resolution 20x20, and (2) with a finer mesh, producing high-resolution (HR) data of spatial resolution 100x100. These data points can be treated as images with the quantities of interest as their channels and can be used to train machine learning models to learn a mapping from the LR images as inputs to the HR images as outputs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Density functional theory calculation of helium segregation and decohesion effect in W110/Ni111 interphase boundary

Density functional theory (DFT) is employed to study the preferential distribution and decohesion effect of He in a W-NiFe composite consisting of W particles embedded in an Ni-based solid solution matrix. A slab containing {110}<100>W//{111}<110>Ni interface is employed as a surrogate model for the W-NiFe system. Firstly, the fracture energy of the W/Ni interphase boundary (IB) (4.37 J/m 2 ) is higher than the cleavage energy of Ni{111} (3.82 J/m^2) and lower than the cleavage energy of W{110} (6.60 J/m 2 ). The comparison shows that the cohesion of the IB is stronger than the Ni{111} planes that are away from the IB. However, the cohesion between the Ni{111} planes adjacent to the IB is found to be the weakest in this system, with a cleavage energy of 3.11 J/m 2 . Subsequently, formation energy of He is calculated in the Ni slab, W slab, and various interstitial sites in the IB. The DFT calculations show that He is significantly more stable in Ni than in W, by about 1.75 eV. Interestingly, He does not prefer to segregate at the IB as compared to bulk Ni. Nevertheless, it prefers to segregate to the region between the Ni{111} planes adjacent to the IB, and decreases the cohesion of the already weakest region. Based on an estimated amount of He gas production in 5 years under first wall neutron irradiation (neutron flux of 1.04x10 15 n/cm 2 /s), the He will decrease the cleavage energy of the weakest region by 21.2%, 15.4%, and 12.2% at 800, 1000, and 1200 C, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Embedding hard physical constraints in neural network coarse-graining of three-dimensional turbulence

In recent years, deep learning approaches have shown much promise in modeling complex systems in the physical sciences. A major challenge in deep learning of partial differential equations is enforcing physical constraints and boundary conditions. In this work, we propose a general framework to directly embed the notion of an incompressible fluid into convolutional neural networks, and apply this to coarse-graining of turbulent flow. These physics-embedded neural networks leverage interpretable strategies from numerical methods and computational fluid dynamics to enforce physical laws and boundary conditions by taking advantage the mathematical properties of the underlying equations. Here, we demonstrate results on three-dimensional fully developed turbulence, showing that this technique drastically improves local conservation of mass, without sacrificing performance according to several other metrics characterizing the fluid flow.

97 MATHEMATICS AND COMPUTING↗