SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Implementation of wall boundary conditions for transpiration in F3D thin-layer Navier-Stokes code

Numerical boundary conditions for mass injection/suction at the wall are incorporated in the thin-layer Navier-Stokes code, F3D. The accuracy of the boundary conditions and the code is assessed by a detailed comparison of the predictions of velocity distributions and skin-friction coefficients with exact similarity solutions for laminar flow over a flat plate with variable blowing/suction, and measurements for turbulent flow past a flat plate with uniform blowing. In laminar flow, F3D predictions for friction coefficient compare well with exact similarity solution with and without suction, but produces large errors at moderate-to-large values of blowing. A slight Mach number dependence of skin-friction coefficient due to blowing in turbulent flow is computed by F3D code. Predicted surface pressures for turbulent flow past an airfoil with mass injection are in qualitative agreement with measurements for a flat plate.

Kandula, M.↗

NASA NTRS1986

Low density gas dynamic wall boundary conditions

Low density nozzles or large expansion ratio nozzles used in space experience rarefaction effects near their exit in the form of velocity slip and temperature jump at the walls. In addition, the boundary layers become very thick and there is a very strong viscous/inviscid interaction. For these reasons no existing design technique has been found to accurately predict the nozzle flow properties up to the nozzle exit. The objective of this investigation was to examine the slip boundary conditions and formulate them in a form appropriate for use with a full Navier-Stokes numerical code. The viscous/inviscid interaction would automatically be accounted for by using a compressible Navier-Stokes code. Through examination of the interaction of molecules with solid surfaces, a model for the distribution function of the reflected molecules has been determined and this distribution function has been used to develop a new slip boundary condition that can be shown to yield more realistic surface boundary conditions.

Collins, F. G.↗

NASA NTRS1964

Space vehicle guidance-a boundary value formulation. part ii- boundary conditions with parameters

Boundary value formulation of space vehicle guidance problem

GUIDANCE↗

NASA NTRS1961

THE FREE-SURFACE BOUNDARY CONDITIONS WITH SURFACE TENSION

Free-surface boundary conditions with surface tension

FREE BOUNDARY↗

DOE OSTI2021

Verification of Boundary Conditions in MCNP

To verify in MCNP 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 OSTI2015

Soft-spring wall based non-periodic boundary conditions for non-equilibrium molecular dynamics of dense fluids

Non-equilibrium molecular dynamics (MD) simulations require imposition of non-periodic boundary conditions (NPBCs) that seamlessly account for the effect of the truncated bulk region on the simulated MD region. Standard implementation of specular boundary conditions in such simulations results in spurious density and force fluctuations near the domain boundary and is therefore inappropriate for coupled atomistic-continuum calculations. In this work, we present a novel NPBC model that relies on boundary atoms attached to a simple cubic lattice with soft springs to account for interactions from particles which would have been present in an untruncated full domain treatment. We show that the proposed model suppresses the unphysical fluctuations in the density to less than 1% of the mean while simultaneously eliminating spurious oscillations in both mean and boundary forces. The model allows for an effective coupling of atomistic and continuum solvers as demonstrated through multiscale simulation of boundary driven singular flow in a cavity. The geometric flexibility of the model enables straightforward extension to nonplanar complex domains without any adverse effects on dynamic properties such as the diffusion coefficient.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1995

Grid resolution and turbulent inflow boundary condition recommendations for NPARC calculations

The effects of grid resolution and specification of turbulent inflow boundary conditions were examined using the NPARC code with the Baldwin-Lomax and Chien k-e turbulence models. Three benchmark turbulent test cases were calculated: two were wall bounded flows and the third was a compressible mixing layer. The wall bounded flows were essentially insensitive to axial grid density; however, the location of the first point off the wall had a substantial effect on flow solutions. It was determined that the first point off the wall must be in the laminar sublayer (y+ less than or equal to 5) for the entire boundary layer. For the compressible mixing layer cases, the axial grid density affected the capturing of oblique shock waves in the mixing region, but the overall mixing rate was not strongly dependent on grid resolution. In specifying the inflow turbulent boundary conditions, it was very important to match the boundary layer and momentum thicknesses of the two flows entering the mixing region; calculations obtained with smaller or no boundary layers resulted in substantially reduced mixing. The solutions were relatively insensitive to freestream turbulence level.

Georgiadis, Nicholas J.↗

DOE OSTI2017

A fast parallel 3D Poisson solver with longitudinal periodic and transverse open boundary conditions for space-charge simulations

A three-dimensional (3D) Poisson solver with longitudinal periodic and transverse open boundary conditions can have important applications in beam physics of particle accelerators. Here, we present a fast efficient method to solve the Poisson equation using a spectral finite-difference method. This method uses a computational domain that contains the charged particle beam only and has a computational complexity of $O(N_u (logN_{mode}))$, where $N_u$ is the total number of unknowns and $N_{mode}$ is the maximum number of longitudinal or azimuthal modes. This saves both the computational time and the memory usage of using an artificial boundary condition in a large extended computational domain. The new 3D Poisson solver is parallelized using a message passing interface (MPI) on multi-processor computers and shows a reasonable parallel performance up to hundreds of processor cores.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Effect of geometric and magnetic boundary conditions on magnetic islands in 3D force-free ideal MHD equilibria

Three-dimensional boundary features and magnetic perturbations affect the magnetic topology of the force-free ideal MHD equilibria formed within a given volume. The use of the PSI-Tet 3D finite-element equilibrium code reveals the existence of magnetic islands inside such equilibria within boundaries with three-dimensional features in configurations that directly depend on the natures of those features. Here, the nature and limits of the influence of such features are explored, as is the influence of flux boundary conditions to qualitatively identify the structure of magnetic islands in Taylor state equilibria with three-dimensional boundary conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Singularities, boundary conditions and gauge link in the light cone gauge

In this work, we first review the issues on the singularities and the boundary conditions in light cone gauge and how to regularize them properly. Then we will further review how these singularities and the boundary conditions can result in the gauge link at the infinity in the light cone direction in the Drell-Yan process. Except for reviewing, we also have verified that the gauge link at the light cone infinity has no dependence on the path not only for the Abelian field but also for non-Abelian gauge field.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1994

Boundary Conditions for Unsteady Compressible Flows

This paper explores solutions to the spherically symmetric Euler equations. Motivated by the work of Hagstrom and Hariharan and Geer and Pope, we modeled the effect of a pulsating sphere in a compressible medium. The literature available on this suggests that an accurate numerical solution requires artificial boundary conditions which simulate the propagation of nonlinear waves in open domains. Until recently, the boundary conditions available were in general linear and based on nonreflection. Exceptions to this are the nonlinear nonreflective conditions of Thompson, and the nonlinear reflective conditions of Hagstrom and Hariharan. The former are based on the rate of change of the incoming characteristics; the latter rely on asymptotic analysis and the method of characteristics and account for the coupling of incoming and outgoing characteristics. Furthermore, Hagstrom and Hariharan have shown that, in a test situation in which the flow would reach a steady state over a long time, Thompson's method could lead to an incorrect steady state. The current study considers periodic flows and includes all possible types and techniques of boundary conditions. The technique recommended by Hagstrom and Hariharan proved superior to all others considered and matched the results of asymptotic methods that are valid for low subsonic Mach numbers.

Hariharan, S. I.↗

NASA NTRS1993

Nonreflecting boundary conditions for linearized unsteady aerodynamic calculations

The present method for the implementation of nonreflecting boundary conditions in 2D and 3D linearized unsteady flow computations is applied to cases of unsteady flows in turbomachine blade rows. The eigenmodes of a discrete representation of the governing equations are computed and used to construct nonreflecting boundary conditions. In 3D, a mixed numerical method is used; in 2D, the discrete representation of the governing equations is obtained from the discretized equations used by the flow solver itself. Wavenumbers and radial mode shapes are computed.

Hall, Kenneth C.↗

DOE OSTI2017

Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes

The present contribution focuses on the accuracy of reflection-type boundary conditions in the Stokes–Brinkman–Darcy modeling of porous flows solved with the lattice Boltzmann method (LBM), which we operate with the two-relaxation-time (TRT) collision and the Brinkman-force based scheme (BF), called BF-TRT scheme. In parallel, we compare it with the Stokes–Brinkman–Darcy linear finite element method (FEM) where the Dirichlet boundary conditions are enforced on grid vertices. In bulk, both BF-TRT and FEM share the same defect: in their discretization a correction to the modeled Brinkman equation appears, given by the discrete Laplacian of the velocity-proportional resistance force. This correction modifies the effective Brinkman viscosity, playing a crucial role in the triggering of spurious oscillations in the bulk solution. While the exact form of this defect is available in lattice-aligned, straight or diagonal, flows; in arbitrary flow/lattice orientations its approximation is constructed. At boundaries, we verify that such a Brinkman viscosity correction has an even more harmful impact. Already at the first order, it shifts the location of the no-slip wall condition supported by traditional LBM boundary schemes, such as the bounce-back rule. For that reason, this work develops a new class of boundary schemes to prescribe the Dirichlet velocity condition at an arbitrary wall/boundary-node distance and that supports a higher order accuracy in the accommodation of the TRT-Brinkman solutions. For their modeling, we consider the standard BF scheme and its improved version, called IBF; this latter is generalized in this work to suppress or to reduce the viscosity correction in arbitrarily oriented flows. Our framework extends the one- and two-point families of linear and parabolic link-wise boundary schemes, respectively called B-LI and B-MLI, which avoid the interference of the Brinkman viscosity correction in their closure relations. The performance of LBM and FEM is thoroughly evaluated in three benchmark tests, which are run throughout three distinctive permeability regimes. The first configuration is a horizontal porous channel, studied with a symbolic approach, where we construct the exact solutions of FEM and BF/IBF with different boundary schemes. The second problem refers to an inclined porous channel flow, which brings in as new challenge the formation of spurious boundary layers in LBM; that is, numerical artefacts that arise due to a deficient accommodation of the bulk solution by the low-accurate boundary scheme. The third problem considers a porous flow past a periodic square array of solid cylinders, which intensifies the previous two tests with the simulation of a more complex flow pattern. The ensemble of numerical tests provides guidelines on the effect of grid resolution and the TRT free collision parameter over the accuracy and the quality of the velocity field, spanning from Stokes to Darcy permeability regimes. It is shown that, with the use of the high-order accurate boundary schemes, the simple, uniform-mesh-based TRT-LBM formulation can even surpass the accuracy of FEM employing hardworking body-fitted meshes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2014

Global and blowup solutions of a mixed problem with nonlinear boundary conditions for a one-dimensional semilinear wave equation

A mixed problem for a one-dimensional semilinear wave equation with nonlinear boundary conditions is considered. Conditions of this type occur, for example, in the description of the longitudinal oscillations of a spring fastened elastically at one end, but not in accordance with Hooke's linear law. Uniqueness and existence questions are investigated for global and blowup solutions to this problem, in particular how they depend on the nature of the nonlinearities involved in the equation and the boundary conditions. Bibliography: 14 titles. (paper)

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1981

Unconditional instability of inflow dependent boundary conditions in difference approximations to hyperbolic systems

The stability of finite difference approximations to initial boundary hyperbolic systems is studied. As is well known, a proper specification of boundary conditions for such systems is essential for their solutions to be well defined. A discrete analogue of the above is proved - if the numerical boundary conditions are consistent with an inflow part of the problem, they render the overall computation unstable. An example of the inviscid gasdynamics equations is considered.

Tadmor, E.↗

NASA NTRS1983

The unconditional instability of inflow-dependent boundary conditions in difference approximations to hyperbolic systems

Tadmor, E.↗

NASA NTRS1998

Global Artificial Boundary Conditions for Computation of External Flow Problems with Propulsive Jets

We propose new global artificial boundary conditions (ABC's) for computation of flows with propulsive jets. The algorithm is based on application of the difference potentials method (DPM). Previously, similar boundary conditions have been implemented for calculation of external compressible viscous flows around finite bodies. The proposed modification substantially extends the applicability range of the DPM-based algorithm. In the paper, we present the general formulation of the problem, describe our numerical methodology, and discuss the corresponding computational results. The particular configuration that we analyze is a slender three-dimensional body with boat-tail geometry and supersonic jet exhaust in a subsonic external flow under zero angle of attack. Similarly to the results obtained earlier for the flows around airfoils and wings, current results for the jet flow case corroborate the superiority of the DPM-based ABC's over standard local methodologies from the standpoints of accuracy, overall numerical performance, and robustness.

Tsynkov, Semyon↗

DOE OSTI2024

Boundary condition effects on runaway electron mitigation coil modeling for the SPARC and DIII-D tokamaks

Abstract Extended-MHD modeling of planned Runaway Electron Mitigation Coils (REMC) for SPARC and DIII-D is performed with the NIMROD code. A coil has been designed for each machine, with the two differing in shape and location, but both having n = 1 symmetry (with n the toroidal mode number). Compared to previous modeling efforts, three improvements are made to the simulations boundary conditions. First a resistive wall model is used in place of an ideal wall. Second, the ThinCurr code is used to compute the time-dependent 3D fields used as magnetic boundary conditions for the simulations. Third, the simulation boundary is moved from the first-wall location to the Vacuum Vessel (VV), which extends the boundary past the location of the internal REMC. To remove the 3D coil from the simulation domain, an equivalent set of 3D fields is calculated at the VV boundary that produce approximately the same field distribution at the last closed flux surface assuming vacuum between the two. Each of these three boundary condition improvements leads to an improvement in the predicted performance of the REMC for both machines. The resistive wall alone primarily effects the resonance of the coil with the plasma after the TQ, affecting the q-profile evolution in the SPARC modeling, and allowing the applied spectrum to be modified in response to the plasma in the DIII-D modeling. The movement of the simulation boundary has the most significant effect on the RE confinement overall, including in the early stages, particularly for a DIII-D inner wall limited equilibrium, where the RE loss fraction increases from 90% to > 99%, with SPARC RE losses also occurring much earlier when the boundary is placed at the VV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗