Search NASA⌕ Search

SEARCH · Search NASA

Results for “BOUNDARY CONDITIONS”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

An optimization-based strategy for peridynamic-FEM coupling and for the prescription of nonlocal boundary conditions

We develop and analyze an optimization-based method for the coupling of a static peridynamic (PD) model and a static classical elasticity model. The approach formulates the coupling as a control problem in which the states are the solutions of the PD and classical equations, the objective is to minimize their mismatch on an overlap of the PD and classical domains, and the controls are virtual volume constraints and boundary conditions applied at the local-nonlocal interface. Our numerical tests performed on three-dimensional geometries illustrate the consistency and accuracy of our method, its numerical convergence, and its applicability to realistic engineering geometries. We demonstrate the coupling strategy as a means to reduce computational expense by confining the nonlocal model to a subdomain of interest, and as a means to transmit local (e.g., traction) boundary conditions applied at a surface to a nonlocal model in the bulk of the domain.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Implementing contact angle boundary conditions for second-order Phase-Field models of wall-bounded multiphase flows

In the present work, a general formulation is proposed to implement the contact angle boundary conditions for the second-order Phase-Field models, which is applicable to N-phase (N ≥ 2) moving contact line problems. To remedy the issue of mass change due to the contact angle boundary condition, a source term or Lagrange multiplier is added to the original second-order Phase-Field models, which is determined by the consistent and conservative volume distribution algorithm so that the summation of the order parameters and the consistency of reduction are not influenced. To physically couple the proposed formulation to the hydrodynamics, especially for large-density-ratio problems, the consistent formulation is employed. The reduction-consistent conservative Allen-Cahn models are chosen as examples to illustrate the application of the proposed formulation. The numerical scheme that preserves the consistency and conservation of the proposed formulation is employed to demonstrate its effectiveness. Results produced by the proposed formulation are in good agreement with the exact and/or asymptotic solutions. The proposed method captures complex dynamics of moving contact line problems having large density ratios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Verification of Boundary Conditions in MCNP

MCNP (Monte Carlo N-Particle) is code that is used to model the transport of different particles. To verify in MCNP6.1 that arrays formed with reflective and periodic boundary conditions produce a neutron multiplication factor (keff) that is equal to or higher than an array of the same size with no boundary conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Calibration of hyperelastic constitutive models: the role of boundary conditions, search algorithms, and experimental variability

The calibration of hyperelastic constitutive models of soft tissue and tissue surrogates is often treated as an exercise in curve-fitting to the average experimental response, and many of the complicating factors such as experimental boundary conditions and data variability are ignored. Here, we focus on three questions that arise in this area: the ramifications of ignoring the experimental boundary conditions, the use of local optimizers, and the role of data variability. Using data from a uniaxial extension experiment on a tissue surrogate, we study how these three factors affect the calibration of isotropic hyperelastic constitutive models. Our results show that even with the simplest of constitutive models, it is necessary to look beyond a "good fit" to the average.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

On the prescription of boundary conditions for nonlocal Poisson's and peridynamics models

We introduce a technique to automatically convert local boundary conditions into nonlocal volume constraints for nonlocal Poisson’s and peridynamic models. The proposed strategy is based on the approximation of nonlocal Dirichlet or Neumann data with a local solution obtained by using available boundary, local data. The corresponding nonlocal solution converges quadratically to the local solution as the nonlocal horizon vanishes, making the proposed technique asymptotically compatible. The proposed conversion method does not have any geometry or dimensionality constraints and its computational cost is negligible, compared to the numerical solution of the nonlocal equation. The consistency of the method and its quadratic convergence with respect to the horizon is illustrated by several two-dimensional numerical experiments conducted by meshfree discretization for both the Poisson’s problem and the linear peridynamic solid model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A boundary value “reservoir problem” and boundary conditions for multi-moment multifluid simulations of sheaths

Multifluid simulations of plasma sheaths are increasingly used to model a wide variety of problems in plasma physics ranging from global magnetospheric flows around celestial bodies to plasma–wall interactions in thrusters and fusion devices. For multifluid problems, accurate boundary conditions to model an absorbing wall that resolves a classical sheath remain an open research area. This work justifies the use of vacuum boundary conditions for absorbing walls to show comparable accuracy between a multifluid sheath and lower moments of a continuum-kinetic sheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A System-Level Model for Estimating Residual Strain and Life of Nuclear Reactor Coolant System Components Under Connected-System-Thermal–Mechanical Boundary Conditions

Background: Environmental-assisted fatigue (EAF) is a major issue for the long-term survival of nuclear power plant fleets in the U.S. and worldwide. Multi-material welded regions (e.g., nozzles) and other high-stress regions of reactor coolant system (RCS) components are prone to EAF-related damage. Objective: The discussed work describes a system-level finite element (FE) model of RCS components of a pressurized water reactor (PWR). This is with the goal of predicting the stress hotspots, strain residuals, strain amplitudes and the resulting fatigue lives. Methods: Herein, the FE model was developed considering system-level loading conditions (under connected system thermal-mechanical boundary conditions). Thermal-mechanical stress analysis was performed considering thermal stratification and a design-basis reactor loading cycle. Based on the FE model results, the strain residuals, strain amplitudes and resulting fatigue lives of RCS components were predicted. Results: The results show that some of the RCS components can have significantly different strain amplitudes, residual strain, and fatigue lives, despite having similar geometry and material. Higher residual strain can lead to accelerated cyclic hardening of material and the associated effect of EAF. The simulated component-level strain profile (under realistic multi-axial-multi-physics loading cycle) can guide the selection of appropriate test inputs for conducting laboratory-scale EAF tests, which is a focus of future works. Conclusions: Despite similar geometry and material the RCS component can have significantly different strain profiles and resulting fatigue lives. The discussed approach can help to identify and prioritize the RCS components for conducting expensive nondestructive evaluation (NDE) inspections.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Role of Inlet Boundary Conditions on Fuel-Air Mixing at Supercritical Conditions

Advanced engine design and alternative fuels present the possibility of fuel injection at purely supercritical conditions in diesel engines and gas turbines. The complex interactions that govern this phenomenon still need significant research, particularly the boundary conditions for fuel injection are critical for accurate simulation. However, the flow inside the injector itself is often omitted to reduce the computational efforts, and thus, velocity, mass flux, or total pressure is specified at the injector exit (or domain inlet), often with simplified velocity profiles and turbulence levels. This simplified inlet boundary treatment has minimal effects on results for conventional fuel injection conditions, however, the validity of this approach at supercritical conditions has not been assessed. Comprehensive real-gas and binary fluid mixing models have been implemented for computational fluid dynamics (CFD) analysis of fuel-air mixing at supercritical conditions. The model is verified using prior CFD results from the literature. The model is used to investigate the effects of the shape of axial velocity and mass fraction profiles at the inlet boundary with the goal to improve the comparison of predictions to experimental data. Results show that the boundary conditions have a significant effect on the predictions, and none of the cases match precisely with experimental data. Furthermore, the study reveals that the physical location of the inlet boundary might be difficult to infer correctly from the experiments and highlights the need for high-quality, repeatable measurements at supercritical conditions to support the development of relevant high-fidelity models for fuel-air mixing.

02 PETROLEUM↗
DOE OSTI2022

Closed-form existence conditions for bandgap resonances in a finite periodic chain under general boundary conditions

Bragg scattering in periodic media generates bandgaps, frequency bands where waves attenuate rather than propagate. Yet, a finite periodic structure may exhibit resonance frequencies within these bandgaps. This is caused by boundary effects introduced by the truncation of the nominal infinite medium. Previous studies of discrete systems determined existence conditions for bandgap resonances, although the focus has been limited to mainly periodic chains with free–free boundaries. In this paper, we present closed-form existence conditions for bandgap resonances in discrete diatomic chains with general boundary conditions (free–free, free–fixed, fixed–free, or fixed–fixed), odd or even chain parity (contrasting or identical masses at the ends), and the possibility of attaching a unique component (mass and/or spring) at one or both ends. The derived conditions are consistent with those theoretically presented or experimentally observed in prior studies of structures that can be modeled as linear discrete diatomic chains with free–free boundary conditions. An intriguing case is a free–free chain with even parity and an arbitrary additional mass at one end of the chain. Introducing such an arbitrary mass underscores a transition among a set of distinct existence conditions, depending on the type of chain boundaries and parity. The proposed analysis is applicable to linear periodic chains in the form of lumped-parameter models, examined across the frequency spectrum, as well as continuous granular media models, or similar configurations, examined in the low-frequency regime.

Bastawrous, Mary V.↗
DOE OSTI2020

Application of the Recharge Estimation Tool (RET) to Prepare Spatially and Temporally Variable Recharge Boundary Conditions for Hanford Site Composite Analysis Vadose Zone Models

This environmental calculation file (ECF) describes the development of a tool for translating recharge estimates into readable input for STOMP© (Subsurface Transport Over Multiple Phases) vadose zone models primarily supporting the vadose zone (VZ) facets of the updated Hanford Site Composite Analysis (CA) and the Hanford Site Cumulative Impact Evaluation (CIE). The recharge estimates are spatiotemporally variable and are produced by the Recharge Evolution Tool (RET) described in Hanford Site-wide Natural Recharge Boundary Condition for Groundwater Models (ECF-HANFORD-15-0019). Outputs from the RET are given in Esri’s™ feature class format with yearly estimates and associated metadata encapsulated in file geodatabase objects. For STOMP models, the translated output is a text file in the format of an input boundary condition card, consistent with STOMP software requirements. The text file contains assimilated spatiotemporal recharge estimates produced by the RET. The tool discussed in this document will be referred to as the “RET2STOMP” tool.

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

Peridynamic elastic waves in two-dimensional unbounded domains: Construction of nonlocal Dirichlet-type absorbing boundary conditions

The focus of this paper is on application of peridynamics (PD) to propagation of elastic waves in unbounded domains. We construct absorbing boundary conditions (ABCs) derived from a semi-analytical solution of the PD governing equation at the exterior region. This solution is made up of a finite series of plane waves, as fundamental solutions (modes), which satisfy the PD dispersion relations. The modes are adjusted to transmit the energy from the interior region (near field) to the exterior region (far field). The corresponding unknown coefficients of the series are found in terms of the displacement field at a layer of points adjacent to the absorbing boundary. This is accomplished through a collocation procedure at subregions (clouds) around each absorbing point. The proposed ABCs offer appealing advantages, which facilitate their application to PD. They are of Dirichlet-type, hence their implementation is relatively simple as no derivatives of the field variables are required. They are constructed in the time and space domains and thus application of Fourier and Laplace transforms, cumbersome for nonlocal models, is not required. At the discrete level, the modes satisfy the same numerical dispersion relations of the near field, which makes the far-field solution compatible with that of the near field. We scrutinize the performance of the proposed ABCs through several examples. So our investigation shows that the proposed ABCs perform stably in time with an appropriate level of accuracy even in problems characterized by highly-dispersive propagating waves, including crack propagation in semi-unbounded brittle solids.

42 ENGINEERING↗
DOE OSTI2020

Dirichlet absorbing boundary conditions for classical and peridynamic diffusion-type models

Diffusion-type problems in (nearly) unbounded domains play important roles in various fields of fluid dynamics, biology, and materials science. The aim of this paper is to construct accurate absorbing boundary conditions (ABCs) suitable for classical (local) as well as nonlocal peridynamic (PD) diffusion models. The main focus of the present study is on the PD diffusion formulation. The majority of the PD diffusion models proposed so far are applied to bounded domains only. In this study, we propose an effective way to handle unbounded domains both with PD and classical diffusion models. For the former, we employ a meshfree discretization, whereas for the latter the finite element method (FEM) is employed. The proposed ABCs are time-dependent and Dirichlet-type, making the approach easy to implement in the available models. The performance of the approach, in terms of accuracy and stability, is illustrated by numerical examples in 1D, 2D, and 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Conservative discontinuous Galerkin interpolation: Sheared boundary conditions

Local studies of accretion disks and laboratory magnetized plasmas employ analytical coordinate mappings that introduce sheared boundary conditions (BCs). Here, we present a discontinuous Galerkin (DG) algorithm to apply such BCs based on projections and quadrature-free integration. The procedure is high-order accurate, preserves moments exactly and works in multiple dimensions. Tests of increasing complexity are provided, beginning with translations of one and two dimensional fields, followed by 3D and 5D simulations with sheared (twist-shift) BCs. Results show that the algorithm is (p + 1)-order accurate in the DG representation and (p + 2) -order accurate in the cell averages, with p being the order of the polynomial basis. Quantification of the algorithm's hyperdiffusion and discussion of aliasing errors are given. This technique enables conservative local simulations of plasma turbulence with DG, not possible until now.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Numerical and boundary condition effects on the prediction of detonation engine behavior using detailed numerical simulations

High-fidelity numerical simulations of an experimental rotating detonation engine with discrete fuel/air injection were conducted. A series of configurations with different feed-plenum pressures but with constant equivalence ratio were studied. Detailed chemical kinetics for the hydrogen/air system is used. A resolution study for the full rotating detonation engine (RDE) system simulation is also conducted. Two kinds of boundary conditions, a total pressure boundary and a constant mass flow rate boundary, are used to assess the effects of the inlet boundary. As mass flow rate is increased, the total pressure boundary causes more error in the axial pressure distribution while the constant mass flow rate gives a better solution for all cases ran. The simulations confirm experimental findings, and reproduce qualitative as well as some of the quantitative trends. These results demonstrate that a) fuel-air mixing is highly non-uniform within the detonation chamber, leading to variations in local equivalence ratio, b) the fuel and oxidizer injectors experience significant backflow as the detonation wave passes over, but recover at different rates which further augments the inefficiencies in mixing, and c) parasitic combustion in the mixing region makes the detonation wave weak by extending the reaction zone across the wave.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Dislocation-graphene interactions in Cu/graphene composites and the effect of boundary conditions: a molecular dynamics study

Here, in this paper, systematic molecular dynamics simulations were employed to study dislocation-graphene interactions in Cu/graphene systems using two carefully selected configurations: (i) edge dislocations in the pileup interacting with graphene during nanoindentation and (ii) complex dislocations interacting with graphene in compression of layered nanopillars. The intrinsic and extrinsic size effects were investigated with respect to varying the Cu lamella thickness and pillar sizes, revealing an anomalous extrinsic size effect: the smaller, the weaker. To understand which boundary conditions are physically consistent, both free and periodic graphene were considered. It was seen that edge dislocations in the pileup continuously transmitted across the free graphene through Cu/graphene interfacial sliding and graphene reorientation, whereas transmission was very difficult for periodic graphene as its in-plane deformation was required. In compression of Cu/graphene nanopillars, the strengthening effect of free graphene was found to be less obvious than that of the periodic case, which could be attributed to free graphene edges acting as dislocation sources after the compressed Cu overflowed the graphene sheet from all directions. We therefore conclude that the strengthening effect of periodic graphene inclusions appears to be overestimated, such that free graphene should be more appropriate.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effective data sampling strategies and boundary condition constraints of physics-informed neural networks for identifying material properties in solid mechanics

Material identification is critical for understanding the relationship between mechanical properties and the associated mechanical functions. However, material identification is a challenging task, especially when the characteristic of the material is highly nonlinear in nature, as is common in biological tissue. In this work, we identify unknown material properties in continuum solid mechanics via physics-informed neural networks (PINNs). To improve the accuracy and efficiency of PINNs, we develop efficient strategies to nonuniformly sample observational data. We also investigate different approaches to enforce Dirichlet-type boundary conditions (BCs) as soft or hard constraints. Finally, we apply the proposed methods to a diverse set of time-dependent and time-independent solid mechanic examples that span linear elastic and hyperelastic material space. The estimated material parameters achieve relative errors of less than 1%. As such, this work is relevant to diverse applications, including optimizing structural integrity and developing novel materials.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Absorbing boundary conditions in material point method adopting perfectly matched layer theory

This study focuses on solving the numerical challenges of imposing absorbing boundary conditions for dynamic simulations in the material point method (MPM). To attenuate elastic waves leaving the computational domain, the current work integrates the Perfectly Matched Layer (PML) theory into the implicit MPM framework. The proposed approach introduces absorbing particles surrounding the computational domain that efficiently absorb outgoing waves and reduce reflections, allowing for accurate modeling of wave propagation and its further impact on geotechnical slope stability analysis. The study also includes several benchmark tests to validate the effectiveness of the proposed method, such as several types of impulse loading and symmetric and asymmetric base shaking. The conducted numerical tests also demonstrate the ability to handle large deformation problems, including the failure of elasto-plastic soils under gravity and dynamic excitations. The findings extend the capability of MPM in simulating continuous analysis of earthquake-induced landslides, from shaking to failure.

58 GEOSCIENCES↗
DOE OSTI2020

Analysis of the performance of different culvert boundary conditions in 2D shallow flow models

Culverts allow roads to safely traverse small streams or drainage ditches, and their proper design is critical to ensure a safe and reliable transportation network. A correct modelization of these hydraulic structures becomes crucial in the assessment of flood footprints or discharge peak estimation in a risk evaluation plan. The question of how to include culverts comes up frequently when assembling a hydraulic model that requires the presence of as many singular elements as possible. In this work, three different culvert integrations with the surface domain are studied and compared in the context of a 2D shallow water (SW) model. All of them are based on the Federal Highway Administration (FHWA) formulation for the culvert discharge estimation but differ in complexity and in the interaction with the numerical model for surface flow, some of them as internal boundary conditions. Several steady and unsteady validation test cases are presented and the numerical results are compared with the predictions from HEC-RAS 1D and HY-8 software. Furthermore, the culvert area, shape and their sensitivity to the 2D computational mesh is also analyzed.

54 ENVIRONMENTAL SCIENCES↗