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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Inner Core Rotation from Geomagnetic Westward Drift and a Stationary Spherical Vortex in Earth's Core

The idea that geomagnetic westward drift indicates convective leveling of the planetary momentum gradient within Earth's core is pursued in search of a differentially rotating mean state, upon which various oscillations and secular effects might be superimposed. The desired state conforms to roughly spherical boundary conditions, minimizes dissipative interference with convective cooling in the bulk of the core, yet may aide core cooling by depositing heat in the uppermost core and lower mantle. The variational calculus of stationary dissipation applied to a spherical vortex within the core yields an interesting differential rotation profile akin to spherical Couette flow bounded by thin Hartmann layers. Four boundary conditions are required. To concentrate shear induced dissipation near the core-mantle boundary, these are taken to be: (i) no-slip at the core-mantle interface; (ii) geomagnetically estimated bulk westward flow at the base of the core-mantle boundary layer; (iii) no-slip at the inner-outer core interface; and, to describe magnetic locking of the inner core to the deep outer core, (iv) hydrodynamically stress-free at the inner-outer core boundary. By boldly assuming the axial core angular momentum anomaly to be zero, the super-rotation of the inner core is calculated to be at most 1.5 degrees per year.

Voorhies, C. V.↗
DOE OSTI2024

Parametric Instability of Alfvén Waves and Wave Packets in Periodic and Open Systems

The parametric decay instability of Alfvén waves has been widely studied, but few investigations have examined wave packets of finite size and the effect of different boundary conditions on the growth rate. In this paper, we perform a linear analysis of circular and arc-polarized wave trains and wave packets in periodic and open boundary systems in a low- β plasma. We find that both types of wave are 3–5 times more stable in open boundary conditions compared to periodic. Additionally, once the wave packet width ℓ becomes smaller than the system size L , the growth rate decreases nearly with a power law γ ∝ ℓ / L . This study demonstrates that the stability of a pump wave cannot be separated from the laboratory settings, and that the growth rate of daughter waves depends on the conditions downstream and upstream of the pump wave and on the fraction of volume it fills. Our results can explain simulations and experiments of localized Alfvén waves. They also suggest that Alfvénic fluctuations in the solar wind, including sharp impulses known as switchbacks, can be more stable than traditional theory suggests depending on wind conditions.

Alfven waves↗
NASA NTRS1981

The effect of zonal gradients of sea surface temperature on the Indian Ocean winter monsoon

Several global climate simulations by the 7-layer, 8 x 10 GISS climate model were designed to test the contributions of various surface boundary conditions to the global climate. The model was run with the sun fixed at a perpetual January. In a comparison of run #5, in which realistic January surface boundary conditions were used, with run #4, which was the same except that a zonally symmetric climatological January sea surface temperature (SST) field was used, one of the results was that run #5 provided a better simulation of the Indian Ocean monsoon. A further comparison of the wind fields over the Indian Ocean that were generated by these two model runs is presented.

Cohen, C.↗
NASA NTRS1988

Development of a thermal and structural analysis procedure for cooled radial turbines

A procedure for computing the rotor temperature and stress distributions in a cooled radial turbine are considered. Existing codes for modeling the external mainstream flow and the internal cooling flow are used to compute boundary conditions for the heat transfer and stress analysis. The inviscid, quasi three dimensional code computes the external free stream velocity. The external velocity is then used in a boundary layer analysis to compute the external heat transfer coefficients. Coolant temperatures are computed by a viscous three dimensional internal flow cade for the momentum and energy equation. These boundary conditions are input to a three dimensional heat conduction code for the calculation of rotor temperatures. The rotor stress distribution may be determined for the given thermal, pressure and centrifugal loading. The procedure is applied to a cooled radial turbine which will be tested at the NASA Lewis Research Center. Representative results are given.

Kumar, Ganesh N.↗
NASA NTRS1999

Electrohydrodynamic Stability of a Liquid Bridge: The "ALEX" Experiment

To provide insight into the roles of electrical forces, experiments on the stability of a liquid bridge were carried out during the 1996 Life And Microgravity Science Mission on the space shuttle Columbia. In terrestrial laboratories a Plateau configuration (where the bridge is surrounded by a matched density liquid) is necessary to avoid deformation due to buoyancy. This complicates the electrical boundary conditions, since charge is transported across the liquid-liquid interface. In the microgravity environment, a cylindrical bridge can be deployed in a gas which considerably simplifies the boundary condition. Nevertheless, to provide a tie-in to terrestrial experiments, two-phase experiments were carried out. The agreement with previous work was excellent. Then several experiments were conducted with a bridge deployed in a dielectric gas, SF6. In experiments with steady fields, it was found that the bridge was less stable than predicted by a linearized stability analysis using the Taylor-Melcher leaky dielectric model.

Burcham, C. L.↗
NASA NTRS2002

Discontinuous Galerkin Methods for Turbulence Simulation

A discontinuous Galerkin (DG) method is formulated, implemented, and tested for simulation of compressible turbulent flows. The method is applied to turbulent channel flow at low Reynolds number, where it is found to successfully predict low-order statistics with fewer degrees of freedom than traditional numerical methods. This reduction is achieved by utilizing local hp-refinement such that the computational grid is refined simultaneously in all three spatial coordinates with decreasing distance from the wall. Another advantage of DG is that Dirichlet boundary conditions can be enforced weakly through integrals of the numerical fluxes. Both for a model advection-diffusion problem and for turbulent channel flow, weak enforcement of wall boundaries is found to improve results at low resolution. Such weak boundary conditions may play a pivotal role in wall modeling for large-eddy simulation.

Collis, S. Scott↗
NASA NTRS2005

Source Term Model for Steady Micro Jets in a Navier-Stokes Computer Code

A source term model for steady micro jets was implemented into a non-proprietary Navier-Stokes computer code, OVERFLOW. The source term models the mass flow and momentum created by a steady blowing micro jet. The model is obtained by adding the momentum and mass flow created by the jet to the Navier-Stokes equations. The model was tested by comparing with data from numerical simulations of a single, steady micro jet on a flat plate in two and three dimensions. The source term model predicted the velocity distribution well compared to the two-dimensional plate using a steady mass flow boundary condition, which was used to simulate a steady micro jet. The model was also compared to two three-dimensional flat plate cases using a steady mass flow boundary condition to simulate a steady micro jet. The three-dimensional comparison included a case with a grid generated to capture the circular shape of the jet and a case without a grid generated for the micro jet. The case without the jet grid mimics the application of the source term. The source term model compared well with both of the three-dimensional cases. Comparisons of velocity distribution were made before and after the jet and Mach and vorticity contours were examined. The source term model allows a researcher to quickly investigate different locations of individual or several steady micro jets. The researcher is able to conduct a preliminary investigation with minimal grid generation and computational time.

Waithe, Kenrick A.↗
NASA NTRSDate not supplied

Beyond Guesswork: Code Verification for Acoustic Duct Mode Prediction

The code used in the work (Swirl) is a frequency-domain linearized Euler equation solver that is designed to predict the effect of acoustic liners on acoustic modes propagating in realistic sheared and swirling mean flows, guiding the design of more efficient liner configurations. The purpose of this study is to validate Swirl using the Method Of Manufactured Solutions (MMS). This study also investigated the effect of the integration and spatial differencing methods on the convergence for a given Manufactured Solution. In addition, the effect of boundary condition implementation was tested. Through rigorous validation of the Swirl code using the MMS, this research underscores the significance of integration and spatial differencing methods on convergence, evaluates boundary condition implementations, and systematically assesses the impact of various filters.

Code Verification↗
NASA NTRS1985

Three-dimensional elastic-plastic finite-element analysis of fatigue crack propagation

Fatigue cracks are a major problem in designing structures subjected to cyclic loading. Cracks frequently occur in structures such as aircraft and spacecraft. The inspection intervals of many aircraft structures are based on crack-propagation lives. Therefore, improved prediction of propagation lives under flight-load conditions (variable-amplitude loading) are needed to provide more realistic design criteria for these structures. The main thrust was to develop a three-dimensional, nonlinear, elastic-plastic, finite element program capable of extending a crack and changing boundary conditions for the model under consideration. The finite-element model is composed of 8-noded (linear-strain) isoparametric elements. In the analysis, the material is assumed to be elastic-perfectly plastic. The cycle stress-strain curve for the material is shown Zienkiewicz's initial-stress method, von Mises's yield criterion, and Drucker's normality condition under small-strain assumptions are used to account for plasticity. The three-dimensional analysis is capable of extending the crack and changing boundary conditions under cyclic loading.

Goglia, G. L.↗
NASA NTRS1985

Three-dimensional vibration analysis of a uniform beam with offset inertial masses at the ends

Analysis of a flexible beam with displaced end-located inertial masses is presented. The resulting three-dimensional mode shape is shown to consist of two one-plane bending modes and one torsional mode. These three components of the mode shapes are shown to be linear combinations of trigonometric and hyperbolic sine and cosine functions. Boundary conditions are derived to obtain nonlinear algebraic equations through kinematic coupling of the general solutions of the three governing partial differential equations. A method of solution which takes these boundary conditions into account is also presented. A computer program has been written to obtain unique solutions to the resulting nonlinear algebraic equations. This program, which calculates natural frequencies and three-dimensional mode shapes for any number of modes, is presented and discussed.

Robertson, D. K.↗
NASA NTRS1988

Development of a thermal and structural analysis procedure for cooled radial turbines

A procedure for computing the rotor temperature and stress distributions in a cooled radial turbine is considered. Existing codes for modeling the external mainstream flow and the internal cooling flow are used to compute boundary conditions for the heat transfer and stress analyses. An inviscid, quasi three-dimensional code computes the external free stream velocity. The external velocity is then used in a boundary layer analysis to compute the external heat transfer coefficients. Coolant temperatures are computed by a viscous one-dimensional internal flow code for the momentum and energy equation. These boundary conditions are input to a three-dimensional heat conduction code for calculation of rotor temperatures. The rotor stress distribution may be determined for the given thermal, pressure and centrifugal loading. The procedure is applied to a cooled radial turbine which will be tested at the NASA Lewis Research Center. Representative results from this case are included.

Kumar, Ganesh N.↗
NASA NTRS1993

Dispersion-relation-preserving finite difference schemes for computational acoustics

Time-marching dispersion-relation-preserving (DRP) schemes can be constructed by optimizing the finite difference approximations of the space and time derivatives in wave number and frequency space. A set of radiation and outflow boundary conditions compatible with the DRP schemes is constructed, and a sequence of numerical simulations is conducted to test the effectiveness of the DRP schemes and the radiation and outflow boundary conditions. Close agreement with the exact solutions is obtained.

Tam, Christopher K. W.↗
NASA NTRS2006

CDUCT-LaRC Status - Shear Layer Refraction and Noise Radiation

A proposed boundary condition accounting for shear layer effects within the Ffowcs Williams-Hawkings radiation module of the CDUCT-LaRC code is investigated. The development and numerical justification of the boundary condition formulation are reviewed. An initial assessment of the effectiveness of the shear layer correction is conducted through comparison with experimental data. Preliminary results indicate that the correction provides physically meaningful modifications of the baseline predicted directivity patterns. Trends of peak directivity steepening and shifting that appeared in predicted patterns were found to follow similar structures in measured data, particularly at higher radiation angles.

Nark, Douglas M.↗
NASA NTRS2016

Temporal Gain Correction for X-Ray Calorimeter Spectrometers

Calorimetric X-ray detectors are very sensitive to their environment. The boundary conditions can have a profound effect on the gain including heat sink temperature, the local radiation temperature, bias, and the temperature of the readout electronics. Any variation in the boundary conditions can cause temporal variations in the gain of the detector and compromise both the energy scale and the resolving power of the spectrometer. Most production X-ray calorimeter spectrometers, both on the ground and in space, have some means of tracking the gain as a function of time, often using a calibration spectral line. For small gain changes, a linear stretch correction is often sufficient. However, the detectors are intrinsically non-linear and often the event analysis, i.e., shaping, optimal filters etc., add additional non-linearity. Thus for large gain variations or when the best possible precision is required, a linear stretch correction is not sufficient. Here, we discuss a new correction technique based on non-linear interpolation of the energy-scale functions. Using Astro-HSXS calibration data, we demonstrate that the correction can recover the X-ray energy to better than 1 part in 104 over the entire spectral band to above 12 keV even for large-scale gain variations. This method will be used to correct any temporal drift of the on-orbit per-pixel gain using on-board calibration sources for the SXS instrument on the Astro-H observatory.

X-ray detector X-ray spectrometer X-ray calorimete↗
NASA NTRS2022

A Numerical Study of the Liquid Motion in Titan’s Subsurface Ocean

An ocean filled with liquid water lies beneath the icy surface of several Jovian and Saturnian moons. In such an ocean, the currents are driven by various phenomena such as the tidal forcing, the deformation of the ice shell lying at its top, the temperature gradient resulting from the surface and bottom heat fluxes…The flow induced by the first two forcings can be modelled by means of a 2D depth-averaged model, while the third one generates horizontal and vertical density variations whose effects can only be captured by a 3D baroclinic model. We study the tides of Titan’s subsurface ocean and the impact of the ice shell on the liquid motion by means of the Second-generation Louvain-la-Neuve Ice–ocean Model, SLIM (https://www.climate.be/slim). The impact of the ice shell lying at the top of the ocean is modelled by a surface friction term and surface pressure terms. The latter are a function of the difference between the ocean elevation and the vertical displacement of the shell and the time derivative of this difference. Because of Titan’s appreciable obliquity (0.306°), the tidal motion expected (and found) is similar to the Europa tidal scenario described by Tyler (2008): the surface elevation consists of two bulges rotating around Titan and the associated depth-averaged velocity field consists of two gyres, separated by an area of high speed flow, whose centre follows a sinusoidal path centred on the equator. The ice shell damps the surface motion, thus slowing down the flow, without significantly modifying the spatial patterns of these fields. The depth of the ocean and the mechanical characteristics of the ice shell being poorly constrained, a sensitivity analysis is conducted. The depth-averaged flow slows down when the depth is increased and a lag appears in the tidal phase but the tidal range remains similar. The ice shell mechanical characteristics influences both the elevation and depth-averaged velocity fields in terms of magnitude but does not modify the spatial patterns of these fields. The influence of the surface heat flux is studied by means of the 3D baroclinic version of SLIM. The heat flux derived from Titan’s topography by Kvorka et al. (2018) is used as surface boundary condition for the temperature equation while a uniform bottom heat flux is implemented. Its value is computed assuming that the heat budget of the ocean is at equilibrium. These boundary conditions cause density variations, which impact the hydrodynamics of the ocean. While the flow velocity induced by these variations is two orders of magnitude smaller than the tidal flow, its orientation is time-independent, hence impacting the orientation of the velocity field. Although the variations of ocean surface elevation and speed with respect to the shell mechanical properties can be larger than those induced by the surface heat flux, taking into account the latter results in large variations of the velocity field global patterns, which was not observed when modifying the shell mechanical properties. Future studies should therefore focus on modelling the surface and bottom heat fluxes while uncertainties about the mechanical characteristics of the shell can be tolerated.

Titan↗
NASA NTRS1972

Numerical calculation of boundary-induced interference in slotted or perforated wind tunnels including viscous effects in slots

A numerical method is presented for calculating the incompressible boundary-induced interference in wind tunnels of rectangular cross section with slotted or perforated walls. The method includes a wall representation which is capable of satisfying a generalized homogeneous boundary condition including the effects of viscosity within the slots. The effects of viscosity in the slots are found to be very significant. The method allows for a variation in the boundary conditions along the tunnel walls. The model can be any configuration and can be located anywhere in the test section. The interference can be computed at any point in the test section.

Keller, J. D.↗
NASA NTRS1973

A boundary layer approach to the analysis of atmospheric motion over a surface obstruction

A boundary layer approach for the solution of the flow field induced over a two-dimensional surface obstruction, such as a building or other man-made structure, is proposed. Adopting a specific geometry in the form of a semi-elliptical cylinder, the characteristics of atmospheric shear flow over a rough terrain are coupled with the well-known boundary layer equations. Two approaches are presented to incorporate the pressure field and boundary conditions which exist within the large viscous region over the obstruction. The first considers a region in the immediate vicinity of the body in which the pressure distribution and outer boundary condition on the velocity are computed from potential theory for flow over the elliptical cylinder. The second approach considers a much larger region of influence, extending from the surface to the undisturbed flow at large heights above the obstruction. Methods which appear to provide an improved theoretical model of the flow over the ellipse, such as a technique for simulating the effect of the seperation regions upstream and downstream of the body, are also presented.

Frost, W.↗
NASA NTRS1982

A full potential inverse method based on a density linearization scheme for wing design

A mixed analysis inverse procedure based on the full potential equation in conservation form was developed to recontour a given base wing to produce density linearization scheme in applying the pressure boundary condition in terms of the velocity potential. The FL030 finite volume analysis code was modified to include the inverse option. The new surface shape information, associated with the modified pressure boundary condition, is calculated at a constant span station based on a mass flux integration. The inverse method is shown to recover the original shape when the analysis pressure is not altered. Inverse calculations for weakening of a strong shock system and for a laminar flow control (LFC) pressure distribution are presented. Two methods for a trailing edge closure model are proposed for further study.

Shankar, V.↗