Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary condition”

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 271 records · Page 15

Mathematical analysis of the Navier-Stokes equations with non standard boundary conditions

One of the major applications of the domain decomposition time marching algorithm is the coupling of the Navier-Stokes systems with Boltzmann equations in order to compute transitional flows. Another important application is the coupling of a global Navier-Stokes problem with a local one in order to use different modelizations and/or discretizations. Both of these applications involve a global Navier-Stokes system with nonstandard boundary conditions. The purpose of this work is to prove, using the classical Leray-Schauder theory, that these boundary conditions are admissible and lead to a well posed problem.

Tidriri, M. D.↗

Inhomogeneous Radiation Boundary Conditions Simulating Incoming Acoustic Waves for Computational Aeroacoustics

A set of nonhomogeneous radiation and outflow conditions which automatically generate prescribed incoming acoustic or vorticity waves and, at the same time, are transparent to outgoing sound waves produced internally in a finite computation domain is proposed. This type of boundary condition is needed for the numerical solution of many exterior aeroacoustics problems. In computational aeroacoustics, the computation scheme must be as nondispersive ans nondissipative as possible. It must also support waves with wave speeds which are nearly the same as those of the original linearized Euler equations. To meet these requirements, a high-order/large-stencil scheme is necessary The proposed nonhomogeneous radiation and outflow boundary conditions are designed primarily for use in conjunction with such high-order/large-stencil finite difference schemes.

Tam, Christopher K. W.↗

Investigating the Numerical Stability of Using an Impedance Boundary Condition to Model Broadband Noise Scattering With Acoustic Liners

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for achieving aircraft noise reduction and are characterized by a frequency-dependent impedance value. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. A Broadband Impedance Model will be discussed in which the liner impedance is specified along a wide range of frequencies. The liner impedance boundary condition will be derived and coupled with a time-domain boundary integral equation to model acoustic scattering by a flat plate consisting of both un-lined and lined surfaces. The scattering solution will be obtained iteratively using both spatial and temporal basis functions and the stability will be demonstrated through eigenvalue analysis. Stability will be assessed for its dependence on time step, spatial discretization, as well as temporal basis function order. Both second- and third-order Lagrange temporal basis functions are considered.

acoustics↗

Investigating the Numerical Stability of Using an Impedance Boundary Condition to Model Broadband Noise Scattering With Acoustic Liners

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for achieving aircraft noise reduction and are characterized by a frequency-dependent impedance value. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. A Broadband Impedance Model will be discussed in which the liner impedance is specified along a wide range of frequencies. The liner impedance boundary condition will be derived and coupled with a time-domain boundary integral equation to model acoustic scattering by a flat plate consisting of both un-lined and lined surfaces. The scattering solution will be obtained iteratively using both spatial and temporal basis functions and the stability will be demonstrated through eigenvalue analysis. Stability will be assessed for its dependence on time step, spatial discretization, as well as temporal basis function order. Both second- and third-order Lagrange temporal basis functions are considered.

acoustics↗

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↗

A method for developing K/S boundary conditions

A method for developing the missing general K/S (Kustaanheimo/Stiefel) boundary conditions is presented, with use of the formalism of optimal control theory. As illustrative examples, the method is applied to the transfer between two position and velocity vectors and to the K/S Lambert problem to derive the missing terminal conditions. The necessary equations for a solution are then developed to the K/S Lambert problem with both the fictitious time, s, and the generalized eccentric anomaly, E, as the independent variables. The latter formulation, requiring the solution of only one nonlinear, well-behaved equation in one unknown, E, results in considerable simplification of the problem. This simplification is possible because the energy equation, in the E-formulation, is separable.

Jezewski, D. J.↗

A parametric experimental study of the slotted-wall boundary condition

An experimental study of slotted upper and lower walls in a two dimensional transonic wind tunnel with solid sidewalls is reported. Results are presented for several slot spacings and slot openness ratios. The experimental data were pressure measurements which were made on an airfoil model and on a sidewall near one of the slotted walls. The slotted-wall boundary condition coefficient, which related the pressure and streamline curvature near the wall, was determined from the wall pressure measurements. The measured wall-induced interference was correlated with the experimental values for the boundary condition coefficient. This correlation was compared with theory.

Everhart, J. L.↗

Boundary conditions for the Navier-Stokes equations

Using ideas from kinetic theory, the Navier-Stokes equations are modified in such a way that they can be cast as a set of first order hyperbolic equations. This is achieved by incorporating time dependent terms into the definition of the stress tensor and the heat flux vectors. The boundary conditions are then determined from the theory of characteristics. Because the resulting equations reduce to the traditional Navier-Stokes equations when the steady state is reached, the present approach provides a straightforward scheme for the determination of inflow and outflow boundary conditions. The method is validated by comparing its predictions with exact and numerical solutions of the steady Navier-Stokes equations.

Rudisill, Edgar N., Jr.↗

A Numerical Combination of Extended Boundary Condition Method and Invariant Imbedding Method Applied to Light Scattering by Large Spheroids and Cylinders

The extended boundary condition method (EBCM) and invariant imbedding method (IIM) are two fundamentally different T-matrix methods for the solution of light scattering by nonspherical particles. The standard EBCM is very efficient but encounters a loss of precision when the particle size is large, the maximum size being sensitive to the particle aspect ratio. The IIM can be applied to particles in a relatively large size parameter range but requires extensive computational time due to the number of spherical layers in the particle volume discretization. A numerical combination of the EBCM and the IIM (hereafter, the EBCM+IIM) is proposed to overcome the aforementioned disadvantages of each method. Even though the EBCM can fail to obtain the T-matrix of a considered particle, it is valuable for decreasing the computational domain (i.e., the number of spherical layers) of the IIM by providing the initial T-matrix associated with an iterative procedure in the IIM. The EBCM+IIM is demonstrated to be more efficient than the IIM in obtaining the optical properties of large size parameter particles beyond the convergence limit of the EBCM. The numerical performance of the EBCM+IIM is illustrated through representative calculations in spheroidal and cylindrical particle cases.

invariant imbeddings↗

A parametric experimental study of the interference effects and the boundary-condition coefficient of slotted wind-tunnel walls

An experimental study of slotted upper and lower walls in a two-dimensional transonic wind tunnel with solid sidewalls is reported. Results are presented for several slot spacings and slot openness ratios. The experimental data are pressure measurements which were made on an airfoil model and on a sidewall near one of the slotted walls. The slotted-wall boundary-condition coefficient, which relates the pressure and streamline curvature near the wall, is determined from the wall pressure measurements. The measured wall-induced interference is correlated with the experimental values for the boundary-condition coefficient. This correlation is compared with theory.

Everhart, J. L.↗

Influence of Spanwise Boundary Conditions on Slat Noise Simulations

The slat noise from the 30P/30N high-lift system is being investigated through computational fluid dynamics simulations with the OVERFLOW code in conjunction with a Ffowcs Williams-Hawkings acoustics solver. In the present study, two different spanwise grids are being used to investigate the effect of the spanwise extent and periodicity on the near-field unsteady structures and radiated noise. The baseline grid with periodic boundary conditions has a short span equal to 1/9th of the stowed chord, whereas the other, longer span grid adds stretched grids on both sides of the core, baseline grid to allow inviscid surface boundary conditions at both ends. The results indicate that the near-field mean statistics obtained using the two grids are similar to each other, as are the directivity and spectral shapes of the radiated noise. However, periodicity forces all acoustic waves with less than one wavelength across the span to be two-dimensional, without any variation in the span. The spanwise coherence of the acoustic waves is what is needed to make estimates of the noise that would be radiated from realistic span lengths. Simulations with periodic conditions need spans of at least six slat chords to allow spanwise variation in the low-frequencies associated with the peak of broadband slat noise. Even then, the full influence of the periodicity is unclear, so employing grids with a fine, central region and highly stretched meshes that go to slip walls may be a more efficient means of capturing the spanwise decorrelation of low-frequency acoustic phenomena.

Lockard, David P.↗

Boundary Condition Controls on the High-Sand-Flux Regions of Mars

Wind has been an enduring geologic agent throughout the history of Mars, but it is often unclear where and why sediment is mobile in the current epoch. We investigated whether eolian bed-form (dune and ripple) transport rates are depressed or enhanced in some areas by local or regional boundary conditions (e.g., topography, sand supply/availability). Bedform heights, migration rates, and sand fluxes all span two to three orders of magnitude across Mars, but we found that areas with the highest sand fluxes are concentrated in three regions: Syrtis Major, Hellespontus Montes, and the north polar erg. All regions are located near prominent transition zones of topography (e.g., basins, polar caps) and thermophysical properties (e.g., albedo variations); these are not known to be critical terrestrial boundary conditions. The two regions adjacent to major impact basins (Hellas and Isidis Planitia) showed radially outward upslope winds driving sand movement, although seasonally reversing wind regimes were also observed. The northern polar dunes yielded the highest known fluxes on the planet, driven by summer katabatic winds modulated by the seasonal CO2 cap retreat—processes not known to affect terrestrial dunes. In contrast, southern dune fields (<45°S) were less mobile, likely as a result of seasonal frost and ground ice suppressing sand availability. Results suggest that, unlike on Earth, large-scale topographic and thermophysical variabilities play a leading role in driving sand fluxes on Mars.

Matthew Chojnacki↗

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↗

The Influence of Prescribed Boundary Conditions on Near-Surface Temperature over the Arctic in the MERRA-2 Atmospheric Model

An accurate historical record of evolving Arctic conditions is integral to furthering our understanding of climate processes and to providing a foundation for predicting future climate scenarios in northern high latitudes. Atmospheric reanalyses are seen as an important source of information on the recent past for the data-sparse Arctic region. An assessment of near-surface Arctic air temperatures finds significant discrepancies among the various modern reanalyses. An important point is the treatment of surface boundary conditions: specifically, the sea ice cover and sea surface temperatures (SSTs) over the Arctic Ocean. Reanalyses use different methodologies and data sources for SSTs and sea ice concentration boundary forcing. Notably, the Modern Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) and the European Centre for Medium-Range Weather Forecasts Interim Re-Analysis (ERA-Interim) both use boundary forcing derived from the Operational Sea Surface Temperature and Sea Ice Analysis (OSTIA) over an extended, overlapping period of time. This allows for an examination of differences between the two systems while both concurrently employ the same fractional sea ice coverage. To further understand these differences, an ensemble of AMIP-style simulations using the MERRA-2 atmospheric model - but without data assimilation - shows considerable differences in Arctic temperatures as compared to reanalyses, particularly in autumn and winter months. Results from the AMIP simulations suggest that the surface representation over sea ice used in the MERRA-2 model provides an intrinsic warm bias and obfuscates Arctic Amplification, an established feature present in observations and reanalyses. An additional ensemble of AMIP-style simulations using the MERRA-2 atmospheric model was performed using boundary conditions derived from the ERA-Interim reanalysis. An in-depth comparison of surface temperatures over the Arctic from the two reanalyses and two AMIP-style ensembles will be presented, along with an assessment of the effects of the varying Arctic temperature time series on the atmospheric general circulation and energy budget.

Collow, Allison B. Marquardt↗

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.↗

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↗

Well-posedness of one-way wave equations and absorbing boundary conditions

A one-way wave equation is a partial differential which, in some approximate sense, behaves like the wave equation in one direction but permits no propagation in the opposite one. The construction of such equations can be reduced to the approximation of the square root of (1-s sup 2) on -1, 1 by a rational function r(s) = p sub m (s)/q sub n(s). Those rational functions r for which the corresponding one-way wave equation is well-posed are characterized both as a partial differential equation and as an absorbing boundary condition for the wave equation. We find that if r(s) interpolates the square root of (1-s sup 2) at sufficiently many points in (-1,1), then well-posedness is assured. It follows that absorbing boundary conditions based on Pade approximation are well-posed if and only if (m, n) lies in one of two distinct diagonals in the Pade table, the two proposed by Engquist and Majda. Analogous results also hold for one-way wave equations derived from Chebyshev or least-squares approximation.

Trefethen, L. N.↗