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 1,081 records · Page 60

The Prediction of Noise Scattered by a Wing/Ducted Fan Configuration

In this proof of concept research, a computational method was developed for predicting the sound field created by the scattering of ducted fan engine noise by a blended wing-body (BWB). It was assumed that all acoustic processes were linear and time harmonic with excitation frequency co. Inflow effects were neglected and no penetration boundary conditions were applied to the engine nacelle and BWB surfaces. A scattering approach was adopted in which the total acoustic field is written as the sum of known incident (from the engine duct) and unknown scattered parts. We further assume that the incident field is independent of the scattered field. Application of the above conditions to the equations of linearized acoustics yields the Helmholtz equation (reduced wave equation) for the scattered pressure with Neumann boundary conditions.

Tweed, John↗

Reynolds-Averaged Navier-Stokes Solutions to Flat Plate Film Cooling Scenarios

The predictions of several Reynolds-Averaged Navier-Stokes solutions for a baseline film cooling geometry are analyzed and compared with experimental data. The Fluent finite volume code was used to perform the computations with the realizable k-epsilon turbulence model. The film hole was angled at 35 to the crossflow with a Reynolds number of 17,400. Multiple length-to-diameter ratios (1.75 and 3.5) as well as momentum flux ratios (0.125 and 0.5) were simulated with various domains, boundary conditions, and grid refinements. The coolant to mainstream density ratio was maintained at 2.0 for all scenarios. Computational domain and boundary condition variations show the ability to reduce the computational cost as compared to previous studies. A number of grid refinement and coarsening variations are compared for further insights into the reduction of computational cost. Liberal refinement in the near hole region is valuable, especially for higher momentum jets that tend to lift-off and create a recirculating flow. A lack of proper refinement in the near hole region can severely diminish the accuracy of the solution, even in the far region. The effects of momentum ratio and hole length-to-diameter ratio are also discussed.

Johnson, Perry L.↗

PAN AIR: A Computer Program for Predicting Subsonic or Supersonic Linear Potential Flows About Arbitrary Configurations Using a Higher Order Panel Method: Theory Document (Version 1.1) - Volume 1

An outline of the derivation of the differential equation governing linear subsonic and supersonic potential flow is given. The use of Green's Theorem to obtain an integral equation over the boundary surface is discussed. The engineering techniques incorporated in the PAN AIR (Panel Aerodynamics) program (a discretization method which solves the integral equation for arbitrary first order boundary conditions) are then discussed in detail. Items discussed include the construction of the compressibility transformations, splining techniques, imposition of the boundary conditions, influence coefficient computation (including the concept of the finite part of an integral), computation of pressure coefficients, and computation of forces and moments.

Magnus, Alfred E.↗

Nonlinear 3-D evolution of bounded kinetic Alfven waves due to shear flow and collisionless tearing instability

The three-dimensional nonlinear simulation of the evolution of a kinetic Alfven wave is presented for simple boundary conditions that resemble those of a discrete auroral breakup arc. The model used is a variant of the reduced megnetohydrodynamic description that includes, in addition, the important dispersive effect of electron inertia. It is found that with a 'generator' boundary condition corresponding to a current generator and a fully conducting 'load' boundary, the wave evolves three-dimensionally through a combination of shear flow and collisionless tearing mode instability.

Seyler, Charles E., Jr.↗

TPSAS-NF1676L-36025-DND

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.

Michelle E Rodio↗

Covariant phase space with boundaries

The covariant phase space method of Iyer, Lee, Wald, and Zoupas gives an elegant way to understand the Hamiltonian dynamics of Lagrangian field theories without breaking covariance. The original literature however does not systematically treat total derivatives and boundary terms, which has led to some confusion about how exactly to apply the formalism in the presence of boundaries. In particular the original construction of the canonical Hamiltonian relies on the assumed existence of a certain boundary quantity “B”, whose physical interpretation has not been clear. We here give an algorithmic procedure for applying the covariant phase space formalism to field theories with spatial boundaries, from which the term in the Hamiltonian involving B emerges naturally. Our procedure also produces an additional boundary term, which was not present in the original literature and which so far has only appeared implicitly in specific examples, and which is already nonvanishing even in general relativity with sufficiently permissive boundary conditions. The only requirement we impose is that at solutions of the equations of motion the action is stationary modulo future/past boundary terms under arbitrary variations obeying the spatial boundary conditions; from this the symplectic structure and the Hamiltonian for any diffeomorphism that preserves the theory are unambiguously constructed. We show in examples that the Hamiltonian so constructed agrees with previous results. We also show that the Poisson bracket on covariant phase space directly coincides with the Peierls bracket, without any need for non-covariant intermediate steps, and we discuss possible implications for the entropy of dynamical black hole horizons.

36 MATERIALS SCIENCE↗

Use-Inspired, Process-Oriented GCM Selection: Prioritizing Models for Regional Dynamical Downscaling

Dynamical downscaling is a crucial process for providing regional climate information for broad uses, using coarser-resolution global models to drive higher-resolution regional climate simulations. The pool of global climate models (GCMs) providing the fields needed for dynamical downscaling has increased from the previous generations of the Coupled Model Intercomparison Project (CMIP). However, with limited computational resources, the need for prioritizing the GCMs for subsequent downscaling studies remains. GCM selection for dynamical downscaling should focus on evaluating processes relevant for providing boundary conditions to the regional models and be inspired by regional uses such as the response of extremes to changes in the boundary conditions. This leads to the need for metrics representing processes of relevance to diverse stakeholders and subregions of a domain. Procedures to account for metric redundancy and the statistical distinguishability of GCM rankings are required. Further, procedures for selecting realizations from ensembles of top-performing GCM simulations can be used to span the range of climate change signals in multiple ways. As a result, distinct weighting of metrics and prioritization of particular realizations may depend on user needs. We provide high-level guidelines for such region-specific evaluations and address how CMIP7 might enable dynamical downscaling of a representative sample of high-quality models across representative shared socioeconomic pathways (SSPs).

54 ENVIRONMENTAL SCIENCES↗

Finite Volume Discretization of the Euler Equations in Pronghorn

Modeling flow and heat transfer in high temperature gas reactors (HTGR) requires the ability to model a wide range of flow speeds from slow (natural convection), to intermediate (forced-flow conditions), to supersonic regimes (depressurization) for a wide range of geometries including the pebble bed, upper and lower plenum, and risers. In previous work, Pronghorn has effectively modeled low-to-medium speed flows in scenarios such as the one described in the two-dimensional PBMR-400 benchmark, using its finite-element-based streamline-upwind Petrov-Galerkin (SUPG) stabilized implementation of the Euler equations. However, limitations of this method become apparent when dealing with more complicated geometries (e.g. imposing slip boundary conditions at nodes belonging to two different boundaries) and when gas speeds are fast enough for shocks and supersonic flow to occur. For these problems, the finite-element-based solver lacks robustness and is plagued by slow iterative convergence or even divergence. In order to address these challenges, the Pronghorn code at INL has been updated with new, modified versions of its original equations. The new Pronghorn models are built on the finite volume method with a Harten-Lax-van Leer-Contact (HLLC) Riemann solver based numerical flux method, which (1) allows imposing slip boundary conditions much more robustly and (2) performs well for a wide range of flow speeds. The finite-volume-based flow solver will form the basis for a robust coarse-mesh thermal-hydraulics capability in Pronghorn.

97 MATHEMATICS AND COMPUTING↗

Variations in density and chemical composition at 120 km from chemical and dynamical processes.

Atmospheric parameters show both systematic and random patterns of behavior at 120 km altitude. Variations in density, temperature, pressure, and especially atomic oxygen concentration are important. The variations are particularly significant because constant boundary conditions at this altitude have commonly been assumed in atmospheric model making. It is surprising how well the assumption of constant boundary conditions has served over the past decade. However, their use has probably introduced erroneous concepts into the field of atmospheric structure, or at least held off the recognition and introduction of important new concepts. The variations at 120 km are caused by changes in energy input into the upper atmosphere, including internal gravity waves and tidal energy from below, and changes in the transport processes within the atmosphere.

Johnson, F. S.↗

A Theoretical and Experimental Study of Acoustic Propagation in Multisectioned Circular Ducts

The propagation of plane waves and higher order acoustic modes in a circular multisectioned duct was studied. A unique source array consisting of two concentric rings of sources, providing phase and amplitude control in the radial, as well as circumferential direction, was developed to generate plane waves and both spinning and nonspinning higher order modes. Measurements of attenuation and radial mode shapes were taken with finite length liners inserted between the hard wall sections of an anechoically terminated duct. Materials tested as liners included a glass fiber material and both sintered fiber metals and perforated sheet metals with a honeycomb backing. The fundamental acoustic properties of these materials were studied with emphasis on the attenuation of sound by the liners and the determination of local versus extended reaction behavior for the boundary condition. A search technique was developed to find the complex eigenvalues for a liner under the assumption of a locally reacting boundary condition.

Wyerman, B. R.↗

Off-axis impact of unidirectional composites with cracks: Dynamic stress intensification

The dynamic response of unidirectional composites under off axis (angle loading) impact is analyzed by assuming that the composite contains an initial flaw in the matrix material. The analytical method utilizes Fourier transform for the space variable and Laplace transform for the time variable. The off axis impact is separated into two parts, one being symmetric and the other skew-symmetric with reference to the crack plane. Transient boundary conditions of normal and shear tractions are applied to a crack embedded in the matrix of the unidirectional composite. The two boundary conditions are solved independently and the results superimposed. Mathematically, these conditions reduce the problem to a system of dual integral equations which are solved in the Laplace transform plane for the transformation of the dynamic stress intensity factor. The time inversion is carried out numerically for various combinations of the material properties of the composite and the results are displayed graphically.

Sih, G. C.↗

Fuel optimal maneuvers of spacecraft about a circular orbit

Fuel optimal maneuvers of spacecraft relative to a body in circular orbit are investigated using a point mass model in which the magnitude of the thrust vector is bounded. All nonsingular optimal maneuvers consist of intervals of full thrust and coast and are found to contain at most seven such intervals in one period. Only four boundary conditions where singular solutions occur are possible. Computer simulation of optimal flight path shapes and switching functions are found for various boundary conditions. Emphasis is placed on the problem of soft rendezvous with a body in circular orbit.

Carter, T. E.↗

Investigations of flow field perturbations induced on slotted transonic-tunnel walls

The free-stream interference caused by the flow through the slotted walls of the test sections of transonic wind tunnels has continuously a problem in transonic tunnel testing. The adaptive-wall transonic tunnel is designed to actively control the near-wall boundary conditions by sucking or blowing through the wall. In order to make the adaptive-wall concept work, parameters for computational boundary conditions must be known. These parameters must be measured with sufficient accuracy to allow numerical convergence of the flow field computations and must be measured in an inviscid region away from the model that is placed inside the wind tunnel. The near-wall flow field was mapped in detail using a five-port cone probe that was traversed in a plane transverse to the free-stream flow. The initial experiments were made using a single slot and recent measurements used multiple slots, all with the tunnel empty. The projection of the flow field velocity vectors on the transverse plane revealed the presence of a vortex-like flow with vorticity in the free stream. The current research involves the measurement of the flow field above a multislotted system with segmented plenums behind it, in which the flow is controlled through several plenums simultaneously. This system would be used to control a three-dimensional flow field.

Wu, J. M.↗

Fuel-optimal maneuvers of a spacecraft relative to a point in circular orbit

Fuel optimal maneuvers of spacecraft relative to a body in circular orbit are investigated using a point mass model in which the magnitude of the thrust vector is bounded. All nonsingular optimal maneuvers consist of intervals of full thrust and coast and are found to contain at most seven such intervals in one period. Only four boundary conditions where singular solutions occur are possible. Computer simulation of optimal flight path shapes and switching functions are found for various boundary conditions. Emphasis is placed on the problem of soft rendezvous with a body in circular orbit.

Carter, T. E.↗

A computational study of multiple jet and wall interaction

In the fuel preburner of the Space Shuttle Main Engine, face plate, injector, and baffle erosion have been observed. The observed patterns of erosion suggest that flame attachment to the walls is a contributing factor. To better understand the physical phenomena involved, a portion of the preburner was modeled computationally. The simulated 'preburner' had three two-dimensional jets entering a cavity adjacent to a baffle. The computational model employed the Patankar Spalding algorithm with upwind differencing. The turbulence model was a standard k-epsilon model with wall functions. The effects of incoming boundary conditions on turbulent kinetic energy and dissipation, k and epsilon, was studied. The results indicate a very strong sensitivity to these boundary conditions over certain ranges of values.

Campbell, W.↗

An improved three-dimensional aerodynamics model for helicopter airloads prediction

The importance of unsteady aerodynamic effects due to rotor blade motion is demonstrated for a helicopter in high-speed forward flight. These unsteady effects are modeled as surface boundary conditions to a three-dimensional, unsteady Computational Fluid Dynamics (CFD) code called the Full-Potential Rotor Code (FPR). These boundary conditions cause significant changes in the computed lift and pitching moment at the front and rear of the rotor disk. Airloads from the modified FPR code are then iteratively coupled with the helicopter comprehensive code, CAMRAD/JA. Computed airloads show good agreement with flight-test data when lift values from the FPR code are used in the coupled calculation. However, the computed airloads from CAMRAD/JA along also show good agreement with the experimental data. Thus for this case one cannot demonstrate a significant improvement in computed airloads with the hybrid coupled scheme. The addition of the pitching moment values from the FPR code into the CAMRAD/JA calculation slows down the overall iterative convergence and does not yield any improvement in the final results.

Strawn, Roger C.↗

On the assumption of vanishing temperature fluctuations at the wall for heat transfer modeling

Boundary conditions for fluctuating wall temperature are required for near-wall heat transfer modeling. However, their correct specifications for arbitrary thermal boundary conditions are not clear. The conventional approach is to assume zero fluctuating wall temperature or zero gradient for the temperature variance at the wall. These are idealized specifications and the latter condition could lead to an ill posed problem for fully-developed pipe and channel flows. In this paper, the validity and extent of the zero fluctuating wall temperature condition for heat transfer calculations is examined. The approach taken is to assume a Taylor expansion in the wall normal coordinate for the fluctuating temperature that is general enough to account for both zero and non-zero value at the wall. Turbulent conductivity is calculated from the temperature variance and its dissipation rate. Heat transfer calculations assuming both zero and non-zero fluctuating wall temperature reveal that the zero fluctuating wall temperature assumption is in general valid. The effects of non-zero fluctuating wall temperature are limited only to a very small region near the wall.

Sommer, T. P.↗

New Developments in the Embedded Statistical Coupling Method: Atomistic/Continuum Crack Propagation

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain has been enhanced. The concurrent MD-FEM coupling methodology uses statistical averaging of the deformation of the atomistic MD domain to provide interface displacement boundary conditions to the surrounding continuum FEM region, which, in turn, generates interface reaction forces that are applied as piecewise constant traction boundary conditions to the MD domain. The enhancement is based on the addition of molecular dynamics-based cohesive zone model (CZM) elements near the MD-FEM interface. The CZM elements are a continuum interpretation of the traction-displacement relationships taken from MD simulations using Cohesive Zone Volume Elements (CZVE). The addition of CZM elements to the concurrent MD-FEM analysis provides a consistent set of atomistically-based cohesive properties within the finite element region near the growing crack. Another set of CZVEs are then used to extract revised CZM relationships from the enhanced embedded statistical coupling method (ESCM) simulation of an edge crack under uniaxial loading.

Saether, E.↗