Engineering Papers⌕ Search

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.

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

At least 109 records · Page 6

An outflow boundary condition for aeroacoustic computations

A formulation of boundary condition for flows with small disturbances is presented. The authors test their methodology in an axisymmetric jet flow calculation, using both the Navier-Stokes and Euler equations. Solutions in the far field are assumed to be oscillatory. If the oscillatory disturbances are small, the growth of the solution variables can be predicted by linear theory. Eigenfunctions of the linear theory are used explicitly in the formulation of the boundary conditions. This guarantees correct solutions at the boundary in the limit where the predictions of linear theory are valid.

Hayder, M. Ehtesham↗

Numerical boundary condition procedure for the transonic axisymmetric inverse problem

Two types of boundary condition procedures for the axisymmetric inverse problem are described. One is a Neumann type boundary condition (analogous to the analysis problem) and the other is a Dirichlet type boundary conditon, both requiring special treatments to make the inverse scheme numerically stable. The dummy point concept is utilized in implementing both. Results indicate the Dirichlet type inverse boundary condition is more robust and conceptually simpler to implement than the Neumann type procedure. A few results demonstrating the powerful capability of the newly developed inverse method that can handle both shocked as well as shockless body design are included.

Shankar, V.↗

Boundary conditions for subsonic compressible Navier-Stokes calculations

A systematic study of inflow and outflow boundary conditions for the numerical solution of the compressible Navier-Stokes equations is presented. Combinations of several representative inflow and outflow boundary conditions are applied in the solution of subsonic flow over a flat plate in a finite computational domain. These boundary conditions are evaluated in terms of their effect on the accuracy of the solution and the rate of convergence to a steady state. It is shown that errors in the data specified at the inflow boundary can produce significant errors in the computed flow field. It is also shown that a non-reflecting outflow boundary condition can significantly reduce the total computational time required.

Rudy, D. H.↗

Revisit boundary conditions for the self-adjoint angular flux formulation

We revisit the boundary conditions for SAAF. We derived the equivalent parity variational form ready for coding up. The more rigorous approach of evaluating odd parity should be solving the odd parity equation coupled with the even parity. We proposed a symmetric reflecting boundary condition although neither positive definiteness nor even-odd decoupling is achieved. A simple numerical test verifies the validity of these boundary conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Consistent boundary conditions for Reduced Navier-Stokes (RNS) scheme applied to 3-dimensional internal viscous flows

A consistent and efficient set of boundary conditions is developed for the multi-sweep space-marching pressure-elliptic Reduced Navier-Stokes (RNS) scheme as applied for 3-D internal viscous flow problems. No-slip boundary conditions are directly imposed on the solid walls. There is no iteration procedure required in the cross plane to ensure mass conservation across each marching plane. The finite difference equations forming the coefficient matrix are ordered such that the surface normal velocity is specified on all the solid walls; unlike external flows, a pressure boundary condition in the cross plane is not required. Since continuity is directly satisfied at all points in the flow domain, the first order momentum equations can be solved directly for the pressure without the need for a Poisson pressure correction equation. The procedure developed herein can also be applied with periodic boundary conditions. The analysis is given for general compressible flows. Incompressible flow solutions are obtained, for straight and curved ducts of square cross section, to validate the procedure. These solutions are used to demonstrate the applicability of the RNS scheme, with the improved boundary conditions for internal flows with strong interaction, as would be encountered in ducts and turbomachinery geometries.

Reddy, D. R.↗

Prediction of the Stress–Strain Behavior of Open-Cell Aluminum Foam under Compressive Loading and the Effects of Various RVE Boundary Conditions

The prediction of the mechanical behavior of metallic foams with realistic microstructure and the effects of various boundary conditions on the mechanical behavior is an important and challenging issue in modeling representative volume elements (RVEs). A numerical investigation is conducted to determine the effects of various boundary conditions and cell wall cross sections on the compressive mechanical properties of aluminum foam, including the stiffness, plateau stress and onset strain of densification. The open-cell AA6101-T6 aluminum foam Duocel is used in the analyses in this study. Geometrical characteristics including the cell size, foam relative density, and cross-sectional shape and thickness of the cell walls are extracted from images of the foam. Then, the obtained foam microstructure is analyzed as a 2D model. The ligaments are modeled as shear deformable beams with elastic-plastic material behavior. To prevent interpenetration of the nodes and walls inside the cells with large deformations, self-contact-type frictionless interaction is stipulated between the internal surfaces. Sensitivity analyses are performed using several boundary conditions and cells wall cross-sectional shapes. The predicted results from the finite element analyses are compared with the experimental results. Finally, the most appropriate boundary conditions, leading to more consistent results with the experimental data, are introduced.

36 MATERIALS SCIENCE↗

Boundary conditions for extremal black holes from 2d gravity

We devise new boundary conditions for the near-horizon geometries of extremal BTZ and Kerr black holes, as well as for the ultra-cold limit of the Kerr-de Sitter black hole. These boundary conditions are obtained as the higher-dimensional uplift of recently proposed boundary conditions in two-dimensional gravity. Their asymptotic symmetries consist in the semi-direct product of a Virasoro and a current algebra, of which we determine the central extensions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Multiple Boundary Condition Tests (MBCT) for verification of large space structures

The Multiple Boundary Condition Tests (MBCT) approach is evaluated and recent modifications are described. For the application of MBCT, 12 different boundary conditions are selected and the results of applying MBCT in conjunction with a nonlinear formulation are indicated schematically. It is concluded that the nonlinear formulation enhances the ability to implement the MBCT test approach on large space structures which cannot be ground tested without the artificial boundary conditions incorporated in MBCT. In addition, it leads to significant improvements in the convergence to the correct solution.

Wada, B. K.↗

On the Exit Boundary Condition for One-Dimensional Calculations of Pulsed Detonation Engine Performance

In one-dimensional calculations of pulsed detonation engine (PDE) performance, the exit boundary condition is frequently taken to be a constant static pressure. In reality, for an isolated detonation tube, after the detonation wave arrives at the exit plane, there will be a region of high pressure, which will gradually return to ambient pressure as an almost spherical shock wave expands away from the exit, and weakens. Initially, the flow is supersonic, unaffected by external pressure, but later becomes subsonic. Previous authors have accounted for this situation either by assuming the subsonic pressure decay to be a relaxation phenomenon, or by running a two-dimensional calculation first, including a domain external to the detonation tube, and using the resulting exit pressure temporal distribution as the boundary condition for one-dimensional calculations. These calculations show that the increased pressure does affect the PDE performance. In the present work, a simple model of the exit process is used to estimate the pressure decay time. The planar shock wave emerging from the tube is assumed to transform into a spherical shock wave. The initial strength of the spherical shock wave is determined from comparison with experimental results. Its subsequent propagation, and resulting pressure at the tube exit, is given by a numerical blast wave calculation. The model agrees reasonably well with other, limited, results. Finally, the model was used as the exit boundary condition for a one-dimensional calculation of PDE performance to obtain the thrust wall pressure for a hydrogen-air detonation in tubes of length to diameter ratio (L/D) of 4, and 10, as well as for the original, constant pressure boundary condition. The modified boundary condition had no performance impact for values of L/D > 10, and moderate impact for L/D = 4.

Wilson, Jack↗

Characteristic boundary conditions for three-dimensional transonic unsteady aerodynamics

Characteristic far-field boundary conditions for the three-dimensional unsteady transonic small disturbance potential equation have been developed. The boundary conditions were implemented in the XTRAN3S finite difference code and tested for a flat plate rectangular wing with a pulse in angle of attack; the freestream Mach number was 0.85. The calculated force response shows that the characteristic boundary conditions reduce disturbances that are reflected from the computational boundaries.

Whitlow, W., Jr.↗

Numerical Boundary Conditions for Computational Aeroacoustics Benchmark Problems

Category 1, Problems 1 and 2, Category 2, Problem 2, and Category 3, Problem 2 are solved computationally using the Dispersion-Relation-Preserving (DRP) scheme. All these problems are governed by the linearized Euler equations. The resolution requirements of the DRP scheme for maintaining low numerical dispersion and dissipation as well as accurate wave speeds in solving the linearized Euler equations are now well understood. As long as 8 or more mesh points per wavelength is employed in the numerical computation, high quality results are assured. For the first three categories of benchmark problems, therefore, the real challenge is to develop high quality numerical boundary conditions. For Category 1, Problems 1 and 2, it is the curved wall boundary conditions. For Category 2, Problem 2, it is the internal radiation boundary conditions inside the duct. For Category 3, Problem 2, they are the inflow and outflow boundary conditions upstream and downstream of the blade row. These are the foci of the present investigation. Special nonhomogeneous radiation boundary conditions that generate the incoming disturbances and at the same time allow the outgoing reflected or scattered acoustic disturbances to leave the computation domain without significant reflection are developed. Numerical results based on these boundary conditions are provided.

Tam, Chritsopher K. W.↗

A Streamlined Approach to Spatial Mapping of Complex 3D Thermal Boundary Condition Data

Spatial mapping of 1D, 2D, and particularly 3D thermal boundary condition data is essential for accurate thermal simulation of spacecraft and launch vehicles. This information typically is available in simple ASCII file format but can be problematic to use for thermal simulations. Difficulties that can be encountered include the ability to directly use the available file format, transformation of the data into the desired reference coordinate system in 3D space, and accurate interpolation of the data onto the thermal calculation domain. This paper presents a streamlined approach to spatial mapping of complex unsymmetrical 3D cold wall heat flux data using the fields’ capability in Simcenter3D. A recent example of 3D thermal plume simulation at MSFC will be employed to illustrate the use of text, csv, and Excel files, specifying coordinate systems for transformation of data in 3D space, and use of the various interpolation schemes available for spatial mapping. Spatial boundary condition verification is also very important, and the spacecraft thermal analyst needs effective visualization tools to develop confidence in the boundary condition definitions. Visualization of the spacecraft plume raw field data, the calculated field data, and the interpolated field data onto the thermal calculation domain will be discussed, with a particular focus on real-time viewing of the interpolated data for available interpolation schemes. The objective of this paper is to outline the modeling techniques used with the integrated thermal analysis of a terrestrial lunar lander demonstrator. Specifically, the techniques outlined herein focus on spatial mapping of convective and radiative exhaust plume results from a computational fluid dynamics (CFD) analysis. Two cases were investigated, however, the modeling techniques that were utilized translate to both cases seamlessly. Those two cases were as follows: • Lunar lander demonstrator firing on-pad (worst case). • Lunar lander demonstrator hovering at an elevated position above the launch pad. For the purposes of information control, this paper will focus on the modeling techniques only, using altered or dimensionless data where necessary. BACKGROUND The thermal model examined in this paper is of the lunar lander demonstrator known as XL-1T (terrestrial), born from a collaborative effort between Masten Space Systems (MSS) and NASA, under NASA’s Lunar CATALYST (also known as Lunar Cargo Transportation and Landing by Soft Touchdown) initiative. The lander is a reusable vertical takeoff/vertical landing (VTVL) test bed which is controlled by four throttleable main engines utilizing green hypergolic propellants.

Hawkins, Robert F.↗

Turbulance boundary conditions for shear flow analysis, using the DTNS flow solver

The effects of different turbulence boundary conditions were examined for two classical flows: a turbulent plane free shear layer and a flat plate turbulent boundary layer with zero pressure gradient. The flow solver used was DTNS, an incompressible Reynolds averaged Navier-Stokes solver with k-epsilon turbulence modeling, developed at the U.S. Navy David Taylor Research Center. Six different combinations of turbulence boundary conditions at the inflow boundary were investigated: In case 1, 'exact' k and epsilon profiles were used; in case 2, the 'exact' k profile was used, and epsilon was extrapolated upstream; in case 3, both k and epsilon were extrapolated; in case 4, the turbulence intensity (I) was 1 percent, and the turbulent viscosity (mu(sub t)) was equal to the laminar viscosity; in case 5, the 'exact' k profile was used and mu(sub t) was equal to the laminar viscosity; in case 6, the I was 1 percent, and epsilon was extrapolated. Comparisons were made with experimental data, direct numerical simulation results, or theoretical predictions as applicable. Results obtained with DTNS showed that turbulence boundary conditions can have significant impacts on the solutions, especially for the free shear layer.

Mizukami, M.↗

Boundary conditions at the walls with thermionic electron emission in two temperature modeling of “thermal” plasmas

In this paper, we propose new boundary conditions for the electric potential, the electron energy equation, and the energy equation for heavy particles (ions and neutrals) at the hot walls with thermionic electron emission for two-temperature thermal arc models. The derived boundary conditions assume that the walls are made from refractory metals and, consequently, the erosion of the wall is small and can be neglected. In these boundary conditions, the plasma sheath formed at the electrode is viewed as the interface between the plasma and the wall. The derived boundary conditions allow the calculation of the heat flux to the walls from the plasma. This allows the calculation of the thermionic electron current that makes the model of electrode-plasma interaction self-consistent.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonreflective boundary conditions for high-order methods

A different approach to nonreflective boundary conditions for the Euler equations is presented. This work is motivated by a need for in and outflow boundary conditions that do not limit the useful accuracy of high-order accurate methods. The primary interest is in the propagation and convection of continuous acoustic and convective waves. This new approach employs the exact solution to finite waves to relate interior values and ambient conditions to boundary values. The method is first presented in one dimension and then generalized to multidimensions. Grid refinement studies are used to demonstrate high-order convergence for both one-dimensional and two-dimensional flows.

Atkins, H. L.↗

Nonreflective boundary conditions for high-order methods

A different approach to nonreflective boundary conditions for the Euler equations is presented. This work is motivated by a need for inflow and outflow boundary conditions that do not limit the useful accuracy of high-order accurate methods. The primary interest is in the propagation and convection of continuous acoustic and convective waves. This new approach employs the exact solution to finite waves to relate interior values and ambient conditions to boundary values. The method is first presented in one dimension and then generalized to multidimensions. Grid refinement studies are used to demonstrate high-order convergence for both one-dimensional and two-dimensional flows.

Atkins, H.↗

On nonlocal problems with Neumann boundary conditions: scaling and convergence for nonlocal operators and solutions

Formulations of Neumann-type boundary conditions for boundary value problems in the nonlocal framework are beset with difficulties, some related to the choice of a proper scaling. Here we identify a space-dependent scaling for a nonlocal Neumann operator, for which we prove linear in δ (δ being the radius for the support for the kernel) convergence of the Neumann operator and $\mathcal{O}$(δ 2 ) convergence of solutions to their classical counterparts. The pointwise-like convergence of the nonlocal normal operator is cast as a new type of two-scale operator-point convergence, which we call condensated convergence . The results hold for general integrable kernels, a setting which is favored in numerical simulations. We support this analysis with numerical convergence studies using a piecewise linear discontinuous Galerkin discretization and show an $\mathcal{O}$(δ 2 ) rate of convergence of solutions, also exhibiting an $\mathcal{O}$(h 2 ) convergence, where h is the mesh size.

97 MATHEMATICS AND COMPUTING↗

Absorbing boundary conditions for relativistic quantum mechanics equations

This paper is devoted to the derivation of absorbing boundary conditions for the Klein–Gordon and Dirac equations modeling quantum and relativistic particles subject to classical electromagnetic fields. Microlocal analysis is the main ingredient in the derivation of these boundary conditions, which are obtained in the form of pseudo-differential equations. Basic numerical schemes are derived and analyzed to illustrate the accuracy of the derived boundary conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗