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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Numerical simulation of a relativistic magnetron using a fluid electron model

An approach to numerically modeling relativistic magnetrons, in which the electrons are represented with a relativistic fluid, is described. A principal effect in the operation of a magnetron is space-charge-limited (SCL) emission of electrons from the cathode. We have developed an approximate SCL emission boundary condition for the fluid electron model. This boundary condition prescribes the flux of electrons as a function of the normal component of the electric field on the boundary. Here we show the results of a benchmarking activity that applies the fluid SCL boundary condition to the one-dimensional Child–Langmuir diode problem and a canonical two-dimensional diode problem. Simulation results for a two-dimensional A6 magnetron are then presented. Computed bunching of the electron cloud occurs and coincides with significant microwave power generation. Numerical convergence of the solution is considered. Sharp gradients in the solution quantities at the diocotron resonance, spanning an interval of three to four grid cells in the most well-resolved case, are present and likely affect convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2007

Validation of Blockage Interference Corrections in the National Transonic Facility

A validation test has recently been constructed for wall interference methods as applied to the National Transonic Facility (NTF). The goal of this study was to begin to address the uncertainty of wall-induced-blockage interference corrections, which will make it possible to address the overall quality of data generated by the facility. The validation test itself is not specific to any particular modeling. For this present effort, the Transonic Wall Interference Correction System (TWICS) as implemented at the NTF is the mathematical model being tested. TWICS uses linear, potential boundary conditions that must first be calibrated. These boundary conditions include three different classical, linear. homogeneous forms that have been historically used to approximate the physical behavior of longitudinally slotted test section walls. Results of the application of the calibrated wall boundary conditions are discussed in the context of the validation test.

Walker, Eric L.↗
NASA NTRS1981

Spectral multigrid methods for elliptic equations

An alternative approach which employs multigrid concepts in the iterative solution of spectral equations was examined. Spectral multigrid methods are described for self adjoint elliptic equations with either periodic or Dirichlet boundary conditions. For realistic fluid calculations the relevant boundary conditions are periodic in at least one (angular) coordinate and Dirichlet (or Neumann) in the remaining coordinates. Spectral methods are always effective for flows in strictly rectangular geometries since corners generally introduce singularities into the solution. If the boundary is smooth, then mapping techniques are used to transform the problem into one with a combination of periodic and Dirichlet boundary conditions. It is suggested that spectral multigrid methods in these geometries can be devised by combining the techniques.

Zang, T. A.↗
NASA NTRS1992

Surface tension and buoyancy-driven flow in a non-isothermal liquid bridge

The Navier-Stokes-Boussinesq equations governing the transport of momentum, mass and heat in a nonisothermal liquid bridge with a temperature-dependent surface tension are solved using a vorticity-stream-function formulation together with a nonorthogonal coordinate transformation. The equations are discretized using a pseudo-unsteady semi-implicit finite difference scheme and are solved by the ADI method. A Picard-type iteration is adopted which consists of inner and outer iterative processes. The outer iteration is used to update the shape of the free surface. Two schemes have been used for the outer iteration; both use the force balance normal to the free surface as the distinguished boundary condition. The first scheme involves successive approximation by the direct solution of the distinguished boundary condition. The second scheme uses the artificial force imbalance between the fluid pressure, viscous and capillary forces at the free surface which arises when the boundary condition for force balance normal to the surface is not satisfied. This artificial imbalance is then used to change the surface shape until the distinguished boundary condition is satisfied. These schemes have been used to examine a variety of model liquid bridge situations including purely thermocapillary-driven flow situations and mixed thermocapillary- and bouyancy-driven flow.

Zhang, Yiqiang↗
DOE OSTI2022

Detection of Damage Inducing Impact Events in Vibration Test Fixtures Using Holder Exponents

The Shock and Vibration team within the Test Engineering group at Los Alamos National Laboratory (LANL) aids in the development and qualification of weapon systems for the United States Department of Energy (DOE). Performing tests that replicate shock and vibration environments representative of a weapon’s lifetime is critical in ensuring safety, survivability, and reliability. A key factor in shock and vibration testing is accurate representation of the weapon’s boundary conditions in its end use configuration. The boundary conditions representative of a reentry body in flight are particularly difficult to replicate on an electrodynamic shaker. When executing shock and vibration tests on a reentry body, two system axes are tested with the system in a cantilevered configuration. This cantilevered testing configuration results in a moment load being applied to the system while undergoing vibration environments. A special fixture, known as the Flex-Web, was developed to decouple this moment load from the reentry body in order to more accurately represent real world boundary conditions while replicating flight environments.

42 ENGINEERING↗
NASA NTRS2001

A {1,2}-Order Plate Theory Accounting for Three-Dimensional Thermoelastic Deformations in Thick Composite and Sandwich Laminates

A {1,2}-order theory for laminated composite and sandwich plates is extended to include thermoelastic effects. The theory incorporates all three-dimensional strains and stresses. Mixed-field assumptions are introduced which include linear in-plane displacements, parabolic transverse displacement and shear strains, and a cubic distribution of the transverse normal stress. Least squares strain compatibility conditions and exact traction boundary conditions are enforced to yield higher polynomial degree distributions for the transverse shear strains and transverse normal stress through the plate thickness. The principle of virtual work is used to derive a 10th-order system of equilibrium equations and associated Poisson boundary conditions. The predictive capability of the theory is demonstrated using a closed-form analytic solution for a simply-supported rectangular plate subjected to a linearly varying temperature field across the thickness. Several thin and moderately thick laminated composite and sandwich plates are analyzed. Numerical comparisons are made with corresponding solutions of the first-order shear deformation theory and three-dimensional elasticity theory. These results, which closely approximate the three-dimensional elasticity solutions, demonstrate that through - the - thickness deformations even in relatively thin and, especially in thick. composite and sandwich laminates can be significant under severe thermal gradients. The {1,2}-order kinematic assumptions insure an overall accurate theory that is in general superior and, in some cases, equivalent to the first-order theory.

Tessler, A.↗
NASA NTRS1997

Low-Dispersion Scheme for Nonlinear Acoustic Waves in Nonuniform Flow

The linear dispersion-relation-preserving scheme and its boundary conditions have been extended to the nonlinear Euler equations. This allowed computing, a nonuniform flowfield and a nonlinear acoustic wave propagation in such a medium, by the same scheme. By casting all the equations, boundary conditions, and the solution scheme in generalized curvilinear coordinates, the solutions were made possible for non-Cartesian domains and, for the better deployment of the grid points, nonuniform grid step sizes could be used. It has been tested for a number of simple initial-value and periodic-source problems. A simple demonstration of the difference between a linear and nonlinear propagation was conducted. The wall boundary condition, derived from the momentum equations and implemented through a pressure at a ghost point, and the radiation boundary condition, derived from the asymptotic solution to the Euler equations, have proven to be effective for the nonlinear equations and nonuniform flows. The nonreflective characteristic boundary conditions also have shown success but limited to the nonlinear waves in no mean flow, and failed for nonlinear waves in nonuniform flow.

Baysal, Oktay↗
NASA NTRS1995

Wave propagation and scattering in computational aeroacoustics

This paper describes predictions of model problems in computational aeroacoustics. Two problem classes are considered: The first (Workshop Category 3, Problem 1 and 2) considers two-dimensional wave propagation and non-reflecting boundary conditions in the presence of a mean flow. The second (Workshop Category 4, Problem 1) examines wall boundary conditions. For the last problem we introduce the Impedance Mismatch Method (IMM) to treat the solid wall boundaries. In this method the solid wall is simulated using a wall region in which the characteristic impedance is set to a different value to that in the fluid region. This method has advantages over traditional solid wall boundary conditions including simplicity of coding, speed of computations, and the ability to treat curved boundaries efficiently. Several numerical examples are given in addition to the Workshop Problems. The discretization of the Euler equations is performed in all cases with a Dispersion-Relation-Preserving (DRP) algorithm. The numerical results are compared with either analytical solutions or solutions obtained using traditional solid wall boundary conditions.

Chung, Cathy↗
DOE OSTI2023

The Influence of Carbon Fiber Composite Specimen Design Parameters on Artificial Lightning Strike Current Dissipation and Material Thermal Damage

Previous artificial lightning strike direct effect research has examined a broad range of specimen design parameters. No works have studied how such specimen design parameters and electrical boundary conditions impact the dissipation of electric current flow through individual plies. This article assesses the influence of carbon fiber composite specimen design parameters (design parameters = specimen size, shape, and stacking sequence) and electrical boundary conditions on the dissipation of current and the spread of damage resulting from Joule heating. Thermal-electric finite element (FE) modelling is used and laboratory scale (<1 m long) and aircraft scale (>1 m long) models are generated in which laminated ply current dissipation is predicted, considering a fixed artificial lightning current waveform. The simulation results establish a positive correlation between the current exiting the specimen from a given ply and the amount of thermal damage in that ply. The results also establish that the distance to ground, from the strike location to the zero potential boundary conditions (ground), is the controlling factor which dictates the electric current dissipation in each ply. Significantly, this distance to ground is dependent on each of the specimen shape, dimensions, stacking sequence, and location of ground boundary conditions. Therefore, it is not possible to decouple current dissipation and damage from specimen design and boundary condition setup. However, it is possible to define a specimen size for a given specimen shape, stacking sequence, and waveform which limit the influence of specimen dimensions on the resulting current distribution and damage. For a rectangular specimen design which appears in literature multiple times, as 100 × 150 mm and with a stacking sequence of [45/0/-45/90] 4s , a specimen design of greater than 300 × 200 mm is required to limit the influence of specimen dimensions on current distribution and damage.

36 MATERIALS SCIENCE↗
DOE OSTI2019

What Is the Fractional Laplacian? A Comparative Review with New Results

The fractional Laplacian in $\mathbb{R}^{d}$, which we write as (–Δ) α/2 with, α ϵ (0, 2) has multiple equivalent characterizations. Furthermore, in bounded domains, boundary conditions must be incorporated in these characterizations in mathematically distinct ways, and there is currently no consensus in the literature as to which definition of the fractional Laplacian in bounded domains is most appropriate for a given application. The Riesz (or integral) definition, for example, admits a nonlocal boundary condition, where the value of a function must be prescribed on the entire exterior of the domain in order to compute its fractional Laplacian. Yet, the spectral definition requires only the standard local boundary condition. These differences, among others, lead us to ask the question: “What is the fractional Laplacian?” Beginning from first principles, we compare several commonly used definitions of the fractional Laplacian theoretically, through their stochastic interpretations as well as their analytical properties. Next, we present quantitative comparisons using a sample of state-of-the-art methods. We discuss recent advances on nonzero boundary conditions and present new methods to discretize such boundary value problems: radial basis function collocation (for the Riesz fractional Laplacian) and nonharmonic lifting (for the spectral fractional Laplacian). In our numerical studies, we aim to compare different definitions on bounded domains using a collection of benchmark problems. We consider the fractional Poisson equation with both zero and nonzero boundary conditions, where the fractional Laplacian is defined according to the Riesz definition, the spectral definition, the directional definition, and the horizon-based nonlocal definition. We verify the accuracy of the numerical methods used in the approximations for each operator, and we focus on identifying differences in the boundary behaviors of solutions to equations posed with these different definitions. Through our efforts, we aim to further engage the research community in open problems and assist practitioners in identifying the most appropriate definition and computational approach to use for their mathematical models in addressing anomalous transport in diverse applications.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

A physical approach to the numerical treatment of boundaries in gas dynamics

Two types of boundaries are considered: rigid walls, and artificial (open) boundaries which were arbitrarily drawn somewhere across a wider flow field. A set of partial differential equations (typically, the Euler equations) has an infinite number of solutions, each one defined by a set of initial and boundary conditions. The initial conditions remaining the same, any change in the boundary conditions will produce a new solution. To pose the problem well, a necessary and sufficient number of boundary conditions are prescribed.

Moretti, G.↗
NASA NTRS1998

Application of Aeroelastic Solvers Based on Navier-Stokes Equations

A pre-release version of the Navier-Stokes solver (TURBO) was obtained from MSU. Along with Dr. Milind Bakhle of the University of Toledo, subroutines for aeroelastic analysis were developed and added to the TURBO code to develop versions 1 and 2 of the TURBO-AE code. For specified mode shape, frequency and inter-blade phase angle the code calculates the work done by the fluid on the rotor for a prescribed sinusoidal motion. Positive work on the rotor indicates instability of the rotor. The version 1 of the code calculates the work for in-phase blade motions only. In version 2 of the code, the capability for analyzing all possible inter-blade phase angles, was added. The version 2 of TURBO-AE code was validated and delivered to NASA and the industry partners of the AST project. The capabilities and the features of the code are summarized in Refs. [1] & [2]. To release the version 2 of TURBO-AE, a workshop was organized at NASA Lewis, by Dr. Srivastava and Dr. M. A. Bakhle, both of the University of Toledo, in October of 1996 for the industry partners of NASA Lewis. The workshop provided the potential users of TURBO-AE, all the relevant information required in preparing the input data, executing the code, interpreting the results and bench marking the code on their computer systems. After the code was delivered to the industry partners, user support was also provided. A new version of the Navier-Stokes solver (TURBO) was later released by MSU. This version had significant changes and upgrades over the previous version. This new version was merged with the TURBO-AE code. Also, new boundary conditions for 3-D unsteady non-reflecting boundaries, were developed by researchers from UTRC, Ref. [3]. Time was spent on understanding, familiarizing, executing and implementing the new boundary conditions into the TURBO-AE code. Work was started on the phase lagged (time-shifted) boundary condition version (version 4) of the code. This will allow the users to calculate non-zero interblade phase angles using, only one blade passage for analysis.

Keith, Theo G., Jr.↗
DOE OSTI2022

Multicomponent gas mixture parametric CFD study of condensation heat transfer in small modular reactor system safety

Safety is always the primary concern for designing and analyzing nuclear reactor systems. The requirements for the safety margin for advanced small modular reactor (SMR) systems are targeted even higher than the conventional commercial large-scale nuclear reactors incorporating the passive and inherent safety systems. The SMR systems are designed with the condensation passive containment cooling system (PCCS), which plays a critical role in removing reactor heat during a steam release accident case. However, the presence of non-condensable gas (NCG), like air, reduces the heat transfer performance. This physics phenomenon becomes multifactorial for nuclear reactor containment during a fuel failure accident case that releases hydrogen gas. Besides, the mixture component of steam-air-hydrogen varies in reactor accident cases, which need simulation and validation keeping parameters of importance. Reviews showed that previous studies for SMR’s PCCCS did not cover the condensation heat transfer (CHT) in the presence of multicomponent NCG mixture parametric computational fluid dynamics (CFD) simulation and validation, making a research gap in the SMR design safety. A comprehensive CHT parametric CFD study was performed for SMR PCCS to fill this research gap. This study used experimental data as simulation 3D physics domain inlet and outlet boundary conditions. However, the wall boundary conditions were constant temperature, curve-fit, and annular coolant for verifying the turbulence models. Parametric simulations were performed, verified, and optimized for steam-NCGs mixtures. The multicomponent gases, multiphase mixtures, and fluid film condensation models were applied with associated turbulence models. The results of the parametric study were evaluated for realistic reactor conditions. Results showed that parametric study provided critical insight about the dependency of multicomponent gas mixture parameters that supports reactor safety design, analysis, and licensing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Fuel-air mixing in motored CFR engine at research octane number (RON) relevant condition

This paper presents a three-dimensional (3-D) computational fluid dynamics (CFD) study of a motored cooperative fuel research (CFR) engine at research octane number (RON) relevant condition. The boundary conditions for 3-D simulations were generated with a one-dimensional GT-Power model. For the first time in literature, a carburetor was added to a virtual CFR engine model with 3-D CFD. Therefore, the proposed setup can simulate the fuel and thermal stratifications inside the engine cylinder with realistic detail. The transient simulations in this work were performed within the Reynolds-averaged Navier-Stokes (RANS) framework with a Realizable k-ε turbulence model. Major conclusions from the present work are: (1) The in-cylinder flow of the CFR engine is swirl-dominated due to the existence of the intake valve shroud. (2) There is a significant amount of liquid droplets entering the cylinder during the intake stroke. The maximum instantaneous amount of liquid for 50% PRF 87 (containing 87% iso-octane and 13% n-heptane (v/v)) and 50% ethanol mixture is indicated to be around 26% of total injected fuel mass. (3) The heat of vaporization (HoV) of the fuel is responsible for creating both temperature and charge stratification inside the cylinder.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Spatiotemporal refraction of light in an epsilon-near-zero indium tin oxide layer: frequency shifting effects arising from interfaces

When light travels through a medium in which the refractive index is rapidly changing with time, the light will undergo a shift in its frequency. Significant frequency shifting effects have recently been reported for transparent conductive oxides. These observations have been interpreted as emerging from temporal changes to the propagation phase in a bulk medium resulting from temporal variations in the refractive index. It is an effect referred to as temporal refraction. Here, we show that the frequency shift in an epsilon-near-zero layer made of indium tin oxide originates not only from this bulk response but includes a significant effect resulting from temporal changes to the spatial boundary conditions. This boundary effect can lead to a dominant, opposing shift to the bulk effect for certain angles. Hence, this process gives rise to a frequency shift that can be tailored through the angle, decoupling the amplitude and phase modulation.

Bohn, Justus (ORCID:0000000214376450)↗
NASA NTRS1978

Large-eddy simulation of a turbulent mixing layer

The three dimensional, time dependent (incompressible) vorticity equations were used to simulate numerically the decay of isotropic box turbulence and time developing mixing layers. The vorticity equations were spatially filtered to define the large scale turbulence field, and the subgrid scale turbulence was modeled. A general method was developed to show numerical conservation of momentum, vorticity, and energy. The terms that arise from filtering the equations were treated (for both periodic boundary conditions and no stress boundary conditions) in a fast and accurate way by using fast Fourier transforms. Use of vorticity as the principal variable is shown to produce results equivalent to those obtained by use of the primitive variable equations.

Mansour, N. N.↗
NASA NTRS1980

Computer-aided analysis of nonlinear problems in transport phenomena

The paper describes algorithms for equilibrium and steady-state problems with coefficients in the expansions derived by the Galerkin weighted residual method and calculated from the resulting sets of nonlinear algebraic equations by the Newton-Raphson method. Initial approximations are obtained from nearby solutions by continuation techniques as parameters are varied. The Newton-Raphson technique is preferred because the Jacobian of the solution is useful for continuation, for analyzing the stability of solutions, for detecting bifurcation of solution families, and for computing asymptotic estimates of the effects on any solution of small changes in parameters, boundary conditions, and boundary shape.

Brown, R. A.↗
NASA NTRS1980

Effect of grid system on finite element calculation

Detailed parametric studies of the effect of grid system on finite element calculation for potential flows were made. These studies led to the formulation of a design criteria for optimum mesh system and the development of two methods to generate the optimum mesh system. The guidelines for optimum mesh system are: (1) the mesh structure should be regular; (2) the element should be as regular and equilateral as possible; (3) the distribution of size of element should be consistent with that of flow variables to insure maximum uniformity in error distribution; (4) for non-Dirichlet boundary conditions, smaller boundary elements or higher order interpolation functions should be used; and (5) the mesh should accommodate the boundary geometry as accurately as possible. The results of the parametric studies are presented.

Lee, K. D.↗