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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Non-equilibrium simulations of hydraulic permeation: Role of mechanical boundary conditions in dense membranes

Understanding the mechanisms of water transport in reverse osmosis membranes is critical for improving membrane performance and guiding material design. While classical models describe transport as either solution-diffusion (SD)—involving concentration-driven diffusion through a homogeneous medium—or pore-flow (PF)—involving pressure-driven convection through percolated water channels—their applicability to crosslinked polyamide membranes remains debated. Here, using non-equilibrium molecular dynamics simulations, we investigate the impact of mechanical support conditions on pressure-driven water transport in polyamide membranes across varying crosslink densities and pressure differentials (1000–5000 bar). Two support conditions are considered: graphene-restrained, representing experimentally relevant supported membranes, and freeze-restrained, mimicking a self-supported structure. In graphene-restrained systems, water concentration gradients and constant pressure profiles emerge, consistent with SD theory and incompatible with PF assumptions due to the absence of percolated pores and sub-nanometer voids. In contrast, freeze-restrained systems display uniform water concentration and linearly decreasing pressure at 1000 bar, and exhibit compressibility-induced water gradients and partial percolation at 5000 bar, resembling PF-like behavior. However, the underlying assumptions of PF theory—continuous solvent pathways and pressure transmission through water-filled pores—are not met under most conditions. Our results demonstrate that accurate modeling of reverse osmosis membranes must incorporate realistic mechanical boundary conditions to distinguish between transport mechanisms. For dense polyamide membranes supported by porous substrates, graphene-restrained simulations best reflect experimental setups and support the SD model as the dominant mechanism of water permeation.

molecular dynamics↗
DOE OSTI2022

Slip velocity boundary conditions for the lattice Boltzmann modeling of microchannel flows

Abstract Slip flows in ducts are important in numerous engineering applications, most notably in microchannel flows. Compared to the standard no‐slip Dirichlet condition, the case of slip formulates as a Robin‐type condition for the fluid tangential velocity. Such an increase in mathematical complexity is accompanied by a more challenging numerical transcription. The present work concerns with this topic, addressing the modeling of the slip velocity boundary condition in the lattice Boltzmann method (LBM) applied to steady slow viscous flows inside ducts of nontrivial shapes. As novelty, we extend the newly revised local second‐order boundary (LSOB) Dirichlet fluid flow method [ Philos. Trans. R. Soc. A 378, 20190404 (2020)] to implement the slip velocity condition within the two‐relaxation‐time (TRT) framework. The LSOB follows an in‐node philosophy where its operation principle seeks to explicitly reconstruct the unknown boundary populations in the form of a third‐order accurate Chapman–Enskog expansion, where the wall slip condition is built‐in as a normal Taylor‐type condition. The key point of this approach is that the required first‐ and second‐order momentum derivatives, rather than computed through nonlocal finite difference approximations, are locally determined through a simple local linear algebra procedure, whose formulation is particularly aided by the TRT symmetry argument. To express the obtained derivatives, two approaches are considered, called and , which operate with node and wall variables, respectively. These two formulations are developed to prescribe the physical slip condition over plane and curved walls, including the corners. Their consistency and accuracy characteristics are examined against alternative linkwise strategies to impose the wall slip velocity, such as the kinetic‐based diffusive bounce‐back scheme, the central linear interpolation slip scheme, and the multireflection slip scheme. The several slip schemes are tested over different 3D microchannel configurations, with walls not conforming with the LBM uniform mesh. Numerical tests confirm the advanced accuracy characteristics of the proposed LSOB slip boundary scheme, revealing the added challenge of the wall slip modeling, and that parabolic accuracy is a necessary requirement to reach second‐order accuracy within this problem class.

Silva, Goncalo↗
DOE OSTI2024

Small circle expansion for adjoint QCD 2 with periodic boundary conditions

We study 1 + 1-dimensional SU(N) gauge theory coupled to one adjoint multiplet of Majorana fermions on a small spatial circle of circumference L. Using periodic boundary conditions, we derive the effective action for the quantum mechanics of the holonomy and the fermion zero modes in perturbation theory up to order (gL) 3 . When the adjoint fermion mass-squared is tuned to g 2 N/(2π), the effective action is found to be an example of supersymmetric quantum mechanics with a nontrivial superpotential. We separate the states into the ℤN center symmetry sectors (universes) labeled by p = 0, . . . , N – 1 and show that in one of the sectors the supersymmetry is unbroken, while in the others it is broken spontaneously. These results give us new insights into the (1, 1) supersymmetry of adjoint QCD 2 , which has previously been established using light-cone quantization. When the adjoint mass is set to zero, our effective Hamiltonian does not depend on the fermions at all, so that there are 2 N−1 degenerate sectors of the Hilbert space. This construction appears to provide an explicit realization of the extended symmetry of the massless model, where there are 2 2N−2 operators that commute with the Hamiltonian. We also generalize our results to other gauge groups G, for which supersymmetry is found at the adjoint mass-squared g 2 h ∨ /(2π), where h ∨ is the dual Coxeter number of G.

effective field theories↗
DOE OSTI2023

Isospin 0 and 2 two-pion scattering at physical pion mass using all-to-all propagators with periodic boundary conditions in lattice QCD

A study of two-pion scattering for the isospin channels, 𝐼 = 0 and 𝐼 = 2, using lattice QCD is presented. Möbius domain-wall fermions, on top of the Iwasaki-DSDR gauge action for gluons with periodic boundary conditions, are used for the lattice computations, which are carried out on two ensembles of gauge field configurations generated by the RBC and UKQCD Collaborations with physical masses, inverse lattice spacings of 1.023 and 1.378 GeV, and spatial extents of 𝐿 = 4.63 and 4.58 fm, respectively. The all-to-all propagator method is employed to compute a matrix of correlation functions of two-pion operators. The generalized eigenvalue problem (GEVP) is solved for a matrix of correlation functions to extract phase shifts with multiple states—two pions with a nonzero relative momentum, as well as two pions at rest. Our results for phase shifts for both the 𝐼 = 0 and 𝐼 = 2 channels are consistent with the Roy equation and chiral perturbation theory, though at this preliminary stage our errors for 𝐼 = 0 are large. An important outcome of this work is that we are successful in extracting two-pion excited states, which are useful for studying 𝐾 → 𝜋⁢𝜋 decay, on physical-mass ensembles using the GEVP.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Impact of cement composition, brine concentration, diffusion rate, reaction rate and boundary condition on self-sealing predictions for cement-CO 2 systems

Geological CO 2 storage (GCS) plays an important role in curbing CO 2 emissions by reducing the carbon footprint of difficult to decarbonize operations and achieve negative CO 2 emissions through activities like Bioenergy with Carbon Capture and Storage (BECCS) and Direct Air Carbon Capture and Storage (DACCS). Leakage of CO 2 through wells is an important concern when it comes to deployment of large-scale GCS. Existing wells in sites that are otherwise suitable for GCS can act as conduits for stored CO 2 to escape the reservoir. There is broad consensus that the main risk of leakage through wellbores is via fractures/damaged pathways. The results from several studies evaluating the permeability evolution of cement fractures in wells upon leakage of CO 2 agree that smaller fracture apertures, slower brine velocities and higher brine residence times promote self-sealing of fractures by mineral precipitation. Quantitatively, however, the differences in sealing conditions are significant and are typically attributed to differences in experimental conditions or model assumptions. Here we examine the sensitivity of our model, describing CO 2 leakage through wellbores, to cement composition, brine concentration, diffusion rates, and reaction rates. We also evaluate the impact of the boundary condition to allow comparisons between observations from experiments performed at constant flow rate and model predictions made at constant pressure conditions. Our results show that diffusion and reactions rates have the most impact on the self-sealing criteria for cement-CO 2 systems. In addition, conditions associated with self-sealing of fractures at constant flow rate require longer fractures, smaller fracture apertures, and slower velocities than under constant pressure.

58 GEOSCIENCES↗
DOE OSTI2021

Effects of as-received defects on ceramic matrix composites properties using high-fidelity microstructures with periodic boundary conditions

The presence of microstructural defects in as-received specimens of ceramic matrix composites (CMCs) significantly influences their constitutive response and damage, highlighting the importance of characterization and quantification of these defects for accurate assessment of damage and failure in the service environment. In a recent effort, the authors developed an algorithm to generate stochastic representative volume elements (SRVEs) of Carbon fiber Silicon-Carbide-Nitride matrix (C/SiNC) CMCs based on extensive multiscale material and defect characterization data. This paper implements this algorithm within a commercial finite element solver with periodic boundary conditions (PBCs) for high-fidelity micromechanics analysis and investigation of macroscopic material behavior of C/SiNC composites. Different loading directions are used to predict the global mechanical properties, and the results are in excellent agreement with theoretical (rule of mixture) predictions. Subsequently, the effects of as-received defects on the global and local responses are investigated. The results show that intratow porosity has pronounced degradation effects on the global elastic properties and results in complex stress localization patterns, which can be attributed to potential damage initiation sites.

Khafagy, Khaled H.↗
DOE OSTI2023

Computational Study on the Effect of Thermal Boundary Conditions and Axial Aspect Ratio on Catalytic Oxidative Coupling of Methane

Packed-bed reactor models with a large reactor-diameter-to-particle-diameter ratio (N > 10) are commonly used in catalysis research. Their use is underpinned by strong simplifying assumptions including unidirectional plug flow, negligible radial gradients, small axial diffusion, and isothermal operation. However, the complex exothermic–endothermic nature of reactions like oxidative coupling of methane (OCM) coupled with heat transfer from the catalyst surface and convective heat losses from the reactor wall provides an interesting case study to evaluate packed-bed assumptions and understand axial and radial temperature variations that may impact conversion and selectivity. This computational study investigates the effect of thermal boundary conditions on catalyst performance in laboratory OCM reactors (N = ~13) through computational fluid dynamics (CFD) simulation. Further, the 1-D model typically used for modeling plug flow reactors assumes negligible radial temperature and species concentration gradients. This assumption is evaluated for OCM reactors over a range of geometric axial aspect ratios, defined as the reactor length divided by the diameter, comparing the 1-D with 3-D CFD results. 1-D model predictions of catalyst performance start to deviate from the 3-D model predictions as the axial aspect ratio decreases, specifically as diameter increases at constant length. In the work, the 1-D model assumption is only found to be valid for axial aspect ratios higher than 1.6, a result considered to be relevant for all exothermic catalytic reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A Filon-like integration strategy for calculating exact exchange in periodic boundary conditions: a plane-wave DFT implementation

An efficient and accurate approach for calculating exact exchange and other two-electron integrals has been developed for periodic electronic structure methods. Traditional approaches used for integrating over the Brillouin zone in band structure calculations, e.g. trapezoidal or Monkhorst-Pack, are not accurate enough for two-electron integrals. This is because their integrands contain multiple singularities over the double integration of the Brillouin zone, which with simple integration methods lead to very inaccurate results. A common approach to this problem has been to replace the Coulomb interaction with a screened Coulomb interaction that removes singularities from the integrands in the two-electron integrals, albeit at the inelegance of having to introduce a screening factor which must precomputed or guessed. Instead of introducing screened Coulomb interactions in an ad hoc way, the method developed in this work derives an effective screened potential using a Filon-like integration approach that is based only on the lattice parameters. This approach overcomes the limitations of traditionally defined screened Coulomb interactions for calculating two-electron integrals, and makes chemistry many-body calculations tractable in periodic boundary conditions. This method has been applied to several systems for which conventional DFT methods do not work well, including the reaction pathways for the addition of H 2 to phenol and Au$^{–}_{20}$ nanoparticle, and the electron transfer of a charge trapped state in the Fe(II) containing mica, annite.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Impact of Backing Plate and Thermal Boundary Conditions for High-Speed Friction Stir Welding of 25-mm Thick Aluminum Alloy 7175-T79

Here this study demonstrates high-speed (150 mm/min) single pass friction stir butt joining of 25 mm thick aluminum alloy 7175-T79. To understand the impact of quenching and cooling rate on process responses, joint strength, and grain size distributions across the weld thickness, we performed a series of friction stir welding (FSW) trials in air and with a trailing water spray using steel and composite backing plates (BP). Welds made in air and trailing water spray exhibit significantly different hardness distributions across the nugget, heat affected zone (HAZ), and HAZ minimum hardness as evidenced by detailed microhardness mapping. The effect of trailing water spray (TWS) on joint efficiency overshadows the influence of BP composition due to the vastly different contribution to quenching. TWS also resulted in a multifaceted effect on FSW such as lowering processing temperature, increasing X and Z forces while lowering Y force, and narrow heat affecting zone. Finally, digital image correlation (DIC)-based fracture mode analysis and grain size measurements correlate with the hardness distribution.

42 ENGINEERING↗
DOE OSTI2022

Modelling of edge plasma dynamics with active wall boundary conditions

A self-consistent 2D model is presented for transport in boundary plasma and plasma-facing material walls. Plasma dynamics in the domain is represented by a 2D collisional plasma fluid model in the edge-plasma code UEDGE, and transport of hydrogen and heat in the wall is represented by a system of reaction–diffusion equations in the 1D wall code FACE. To account for variation of parameters along the wall, in the coupled model multiple instances of the FACE code run in parallel. Here, the coupled model provides a tool for investigating a range of dynamic plasma–material interactions phenomena in 2D. For demonstration of its capability, one application of particular interest is the role of active wall in tokamak strike point sweeping proposed for mitigation of divertor heat loads. In the present study, the coupled calculations are applied to investigation of the impact of heat and hydrogen transport in the material wall on the divertor plasma and target heat load during sweeping of the target strike point for parameters of a high-power tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Consistent boundary conditions for PDE filter regularization in topology optimization

Abstract Design variables in density-based topology optimization are typically regularized using filtering techniques. In many cases, such as stress optimization, where details at the boundaries are crucially important, the filtering in the vicinity of the design domain boundary needs special attention. One well-known technique, often referred to as “padding,” is to extend the design domain with extra layers of elements to mitigate artificial boundary effects. We discuss an alternative to the padding procedure in the context of PDE filtering. To motivate this augmented PDE filter, we make use of the potential form of the PDE filter which allows us to add penalty terms with a clear physical interpretation. The major advantages of the proposed augmentation compared with the conventional padding is the simplicity of the implementation and the possibility to tune the boundary properties using a scalar parameter. Analytical results in 1D and numerical results in 2D and 3D confirm the suitability of this approach for large-scale topology optimization.

42 ENGINEERING↗
DOE OSTI2023

Computational Analysis of the Impact of Boundary Conditions on a Particle-Laden Flow: A Case Study in a Pressurized Oxy-Coal Combustor

Designing an effective burner is vital for the development of coal combustion technologies. Because of high pressure, the volumetric fraction of the coal particles in the injected fuel in a pressurized oxy-combustion (POC) burner approaches or even exceeds the limitations allowed by the commercial computational fluid dynamics codes (e.g., Ansys Fluent). Consequently, for such high particle volumetric fractions, the interplay between the particles, the fluid flow, and the burner wall needs to be re-evaluated. The present computational work is a first step in a systematic analysis of the roles of various characteristics involved in the POC process, such as the method of particle release, its location, and the particle size. Specifically, pulverized coal is burned under an elevated pressure of 15 bar in an O 2 /CO 2 environment. A 100 kW, a POC combustor, is modeled with Ansys Fluent using the Reynolds-averaged Navier−Stokes approach. It is revealed that for this pilot-scale, pressurized burner, the gas phase flow velocity in the near-wall region exhibits anomalies. With the major focus on POC, this work aims to eliminate/reduce the impact of high particle loading on the gas-phase flow. To scrutinize the role of particle loading in the near-wall region and eliminate the impact of this velocity on POC downstream, the particle−gas interplay in the boundary layer is investigated by means of the computational simulations incorporating the coupling between the turbulent flow and the particles. It is found that the tuning of the particle release location makes the gas-phase flow velocity in the presence of particles consistent with the pure gas flow velocity profile. The particles size is also found to have a significant impact on the particle trajectory.

CFD↗
DOE OSTI2025

Nucleation and Growth of Blue Phase 3D Cubic Structure Under Continuously Changing Curved Boundary Conditions

The blue phase of liquid crystals (BPLCs) with a cubic lattice of disclination lines and 3D nanostructure enables the modulation of photonic bandgap thus casting them in the category of photonic crystals. Its unique nature promises applications in display technologies, electro‐optics, and sensors. To integrate these ordered materials into wearable devices a fundamental understanding of curvature, and spatial confinement is necessary. Although continuous confinement in flat geometries have been studied, confining curvature has shown to induce strong destabilization effects on the cubic structure and formation of topological defects, thereby deteriorating their optical performance. Moreover, limitations in controlling the curvature of droplets further hinder studies of nucleation and growth of BPLCs. Here, to address these challenges, micro‐scale patterned surfaces of concentric cylinders are exploited to systematically control curvatures. The impact of curvature on the confined BPLCs is revealed in terms of phase transition temperature, nucleation and growth, morphology, as well as phase transformation. This research offers valuable insights into the stability, and structural characteristics of BPLCs in adaptive photonic devices, paving the way for future advancements in flexible displays, sensors, and other technologies leveraging liquid crystal (LC) materials.

3D Cubic Structure↗
DOE OSTI2021

Mutual Diffusion Coefficients of Na 2 SO 4 at C = 0.02989 mol·dm -3 and 298.15 K by Using Rayleigh Interferometry with Free Diffusion Boundary Conditions: Experimental Test of the Effect of the C 3/2 Concentration Dependence of Refractive Index at Low Concentration

Rayleigh interferometry is one of the most precise and accurate methods for the experimental determination of the mutual diffusion coefficients of both electrolytes and non-electrolytes in liquid solutions. For binary solutions at moderate and high concentrations where the concentration dependences of both diffusion coefficient and refractive index can be assumed to be linear or almost linear and the concentration difference ΔC between a pair of solutions undergoing free diffusion is not large compared to their average concentration, and by using proper combinations of the interference fringe positions as they vary with time, under these conditions diffusion coefficients are obtained that are independent of the size of ΔC. However, this situation becomes more complicated at low concentrations for electrolyte solutions where the concentration dependences of both the diffusion coefficient and refractive index are non-linear, with C 1/2 dependence for the diffusion coefficient and C 3/2 dependence for the refractive index. As a test of these effects on calculated mutual diffusion of electrolyte solutions at low concentrations as measured by Rayleigh interferometry, apparent mutual diffusion coefficients D v,a (volume-fixed reference frame) were measured for fixed average concentrations C=(0.02989 ± 0.00001) mol·dm –3 of Na 2 SO 4 (aq) at 298.15 K while varying ΔC/C from 0.66673 to 2, a threefold variation. The dependence of the diffusion coefficient on ΔC/C was found to be linear, yielding D v = 1.041 6 × 10 –9 m 2 ·s –1 for an infinitely small concentration difference between the diffusing solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗