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.

SourceResearch connectionFrom yearThrough year

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

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Evaluation of advanced materials through experimental mechanics and modelling

Composite materials have been frequently used in aerospace vehicles. Very often defects are inherited during the manufacture and damages are inherited during the construction and services. It becomes critical to understand the mechanical behavior of such composite structure before it can be further used. One good example of these composite structures is the cylindrical bottle of solid rocket motor case with accidental impact damages. Since the replacement of this cylindrical bottle is expensive, it is valuable to know how the damages affects the material, and how it can be repaired. To reach this goal, the damage must be characterized and the stress/strain field must be carefully analyzed. First the damage area, due to impact, is surveyed and identified with a shearography technique which uses the principle of speckle shearing interferometry to measure displacement gradient. Within the damage area of a composite laminate, such as the bottle of solid rocket motor case, all layers are considered to be degraded. Once a lamina being degraded the stiffness as well as strength will be drastically decreased. It becomes a critical area of failure to the whole bottle. And hence the stress/strain field within and around a damage should be accurately evaluated for failure prediction. To investigate the stress/strain field around damages a Hybrid-Numerical method which combines experimental measurement and finite element analysis is used. It is known the stress or strain at the singular point can not be accurately measured by an experimental technique. Nevertheless, if the location is far away from the singular spot, the displacement can be found accurately. Since it reflects the true displacement field locally regardless of the boundary conditions, it is an excellent input data for a finite element analysis to replace the usually assumed boundary conditions. Therefore, the Hybrid-Numerical method is chosen to avoid the difficulty and to take advantage of both experimental technique and finite element analysis. Experimentally, the digital image correlation technique is employed to measure the displacement field. It is done by comparing two digitized images, before and after loading. Numerically, the finite element program, ABAQUS (version 5.2), is used to analyze the stress and strain field. It takes advantage of the high speed and huge memory size of modern supercomputer, CRAY Y-MP, at NASA Marshall Space Flight Center.

Yang, Yii-Ching↗

NASA NTRS1993

Model development for Ulysses and SOHO

The purpose of this research is to provide scientific expertise in solar physics and in the development and use of magnetohydrodynamic (MHD) models of coronal structures for the computation of Lyman alpha scattered radiation in these structures. The specific objectives will be to run MHD models with new boundary conditions and compute resulting scattered solar Lyman alpha intensities, guided by results from the first series of boundary conditions.

Wu, S. T.↗

DOE OSTI2025

Explicit block encodings of boundary value problems for many-body elliptic operators

Simulation of physical systems is one of the most promising use cases of future digital quantum computers. In this work we systematically analyze the quantum circuit complexities of block encoding the discretized elliptic operators that arise extensively in numerical simulations for partial differential equations, including high-dimensional instances for many-body simulations. When restricted to rectangular domains with separable boundary conditions, we provide explicit circuits to block encode the many-body Laplacian with separable periodic, Dirichlet, Neumann, and Robin boundary conditions, using standard discretization techniques from low-order finite difference methods. To obtain high-precision, we introduce a scheme based on periodic extensions to solve Dirichlet and Neumann boundary value problems using a high-order finite difference method, with only a constant increase in total circuit depth and subnormalization factor. We then present a scheme to implement block encodings of differential operators acting on more arbitrary domains, inspired by Cartesian immersed boundary methods. We then block encode the many-body convective operator, which describes interacting particles experiencing a force generated by a pair-wise potential given as an inverse power law of the interparticle distance. This work provides concrete recipes that are readily translated into quantum circuits, with depth logarithmic in the total Hilbert space dimension, that block encode operators arising broadly in applications involving the quantum simulation of quantum and classical many-body mechanics.

Kharazi, Tyler [University of California, Berkeley↗

DOE OSTI2024

A weighted shifted boundary method for immersed moving boundary simulations of Stokes' flow

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions: hence the name of the method, that shifts the location and values of the boundary conditions. Here, in this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1994

Separated and Recovering Turbulent Boundary Layer Flow Behind a Backward Facing Step For Different Reynolds Numbers

Experimental results for a two-dimensional separated turbulent boundary layer behind a backward facing step for five different Reynolds numbers are reported. Results are presented in the form of tables, graphs and a floppy disk for an easy access of the data. Reynolds number based on the step height was varied by changing the reference velocity upstream of the step, U(sub o), and the step height, h. Hot-wire measurement techniques were used to measure three Reynolds stresses and four triple-velocity correlations. In addition, surface pressure and skin friction coefficients were measured. All hot-wire measurements were acquired in a measuring domain which excluded recirculating flow region due to the directional insensitivity of hot-wires. The downstream extent of the domain from the step was 51 h for the largest and I 14h for the smallest step height. This significant downstream length permitted extensive study of the flow recovery. Prediction of perturbed flows and their recovery is particularly attractive for popular turbulence models since variations of turbulence length and time scales and flow interactions in different regions are generally inadequately predicted. The data indicate that the flow in the free shear layer region behaves like the plane mixing layer up to about 2/3 of the mean reattachment length when the flow interaction with the wall commences the flow recovery to that of an ordinary turbulent boundary layer structure. These changes of the flow do not occur abruptly with the change of boundary conditions. A reattachment region represents a transitional region where the flow undergoes the most dramatic adjustments to the new boundary conditions. Large eddies, created in the upstream free-shear layer region, are being torn, recirculated, reentrained back into the main stream interacting with the incoming flow structure. It is foreseeable that it is quite difficult to describe the physics of this region in a rational and quantitative manner other than statistical. Downstream of the reattachment point the flow recovers at different rates near the wall, in the newly developing internal boundary layer, and in the outer part of the flow. It appears that Reynolds stresses do not fully recover up to the longest recovery length of 114 h.

Jovic, Srba↗

NASA NTRS1974

Unsteady Subsonic and Transonic Potential Flow over Helicopter Rotor Blades

Differential equations and boundary conditions for a rotor blade in forward flight, with subsonic or transonic tip Mach number, are derived. A variety of limiting flow regimes determined by different limits involving blade thickness ratio, aspect ratio, advance ratio and maximum tip Mach number is discussed. The transonic problem is discussed in some detail, and in particular the conditions that make this problem quasi-steady or essentially unsteady are determined. Asymptotic forms of equations and boundary conditions that are valid in an appropriately scaled region of the tip and an azimuthal sector on the advancing side are derived. The equations are then put in a form that is valid from the blade tip inboard through the strip theory region.

Isom, M. P.↗

NASA NTRS1976

Graded-porosity heat-pipe wicks

To maximize the capacity of a nonarterial heat pipe, a wick is considered whose porosity is allowed to vary axially along its length. At every axial location the porosity is set no lower than required to maintain the wick in a nearly saturated state under the maximum heat-transport rate. The result is a wick whose permeability is everywhere as high as possible. The differential equation that governs the optimum porosity variation is solved numerically between a condenser-end boundary condition that just prevents a liquid slug or puddle in the vapor spaces and an evaporator-end boundary condition that just prevents circumferential groove dry-up. Experimental performance measurements for an ammonia heat pipe are presented.

Eninger, J. E.↗

NASA NTRS1983

Higher modes of the Orr-Sommerfeld problem for boundary layer flows

The discrete spectrum of the Orr-Sommerfeld problem of hydrodynamic stability for boundary layer flows in semi-infinite regions is examined. Related questions concerning the continuous spectrum are also addressed. Emphasis is placed on the stability problem for the Blasius boundary layer profile. A general theoretical result is given which proves that the discrete spectrum of the Orr-Sommerfeld problem for boundary layer profiles (U(y), 0,0) has only a finite number of discrete modes when U(y) has derivatives of all orders. Details are given of a highly accurate numerical technique based on collocation with splines for the calculation of stability characteristics. The technique includes replacement of 'outer' boundary conditions by asymptotic forms based on the proper large parameter in the stability problem. Implementation of the asymptotic boundary conditions is such that there is no need to make apriori distinctions between subcases of the discrete spectrum or between the discrete and continuous spectrums. Typical calculations for the usual Blasius problem are presented.

Lakin, W. D.↗

NASA NTRS1989

Euler procedure for three-dimensional transonic wall interference

Based on an optimization formulation, a procedure has been developed to evaluate Mach number and angle-of-attack corrections. The Euler equations are assumed to be the flow governing equations. To obtain efficient solutions for the optimization problem, the iterative solutions for the flow variables and the design parameters are simultaneously updated. In addition to the model lift and geometry, the procedure requires pressure measurements near the tunnel walls. The tunnel boundary conditions are based on the introduction of Reimann invariants for a one-dimensional flow normal to the boundary. Computations are performed to verify that the errors introduced by this approximate boundary-condition formulation are acceptably small. The correction scheme is applied to an aircraft configuration in an open jet. The results indicate that the optimization scheme is highly efficient with the rate of convergence of the flow solution nearly equal to the corresponding rate of a regular analysis problem.

Rizk, Magdi H.↗

NASA NTRS1992

Long time behavior of unsteady flow computations

This paper addresses a specific issue of time accuracy in the calculation of external aerodynamic problems. The class of problems discussed consists of inviscid compressible subsonic flows. These problems are governed by a convective equation. A key issue that is not understood is the long time behavior of the solution. This is important if one desires transient calculations of problems governed by the Euler equations or its derivatives such as the small disturbance equations or the potential formulations for the gust problem. Difficulties arise for two dimensional problems where the time rate decay solutions of the wave equation is slow. In concert with the above mentioned problem, exterior flows require proper modeling of the boundary conditions. In particular, this requires the truncation of infinite regions into finite regions with the aid of artificial boundaries. These boundary conditions must be consistent with the physics of the unbounded problem as well as consistent in time and space. Our treatment of the problem is discussed in detail and examples are given to verify the results.

Hariharan, S. I.↗

NASA NTRS1994

The resonance lines of MG2 as diagnostics of the upper solar chromosphere

The resonance lines of singly ionized magnesium, the MgII h&k lines at about 280 nm, are two of the small number of lines in the solar spectrum that are optically thick in the chromospheric part of the solar atmosphere. Potentially these lines contain information on the initial temperature rise that occurs at the top of the photosphere. Unfortunately, few good observations of the lines exist due to their wavelength near 280 nm the ultraviolet. However, a fair number of observations (on the order of 200) are available from the data base of the UltraViolet Polarimeter and Spectrometer (UVSP) instrument that flew on board of NASA's Solar Maximum Mission (SMM) satellite. In addition, this data base contains a number of spectra that include the Mg I resonance line at (lambda)285.2nm, just longward of the h&k lines. The neutral magnesium line is not as strong as its ionic counterparts and samples slightly lower parts of the atmosphere. Its width is a sensitive diagnostic of the ionization balance between neutral and singly ionized magnesium, which determines the opacity scale (and formation height) of other diagnostically important MgI lines like the 457.1 nm intercombination line, the magnesium b lines and the infrared MgI emission lines near 12 microns. Analysis of the observed line profiles shows that it is necessary to include the effects of partial frequency redistribution (PRD) in the formation of the line as in the case of the h&k lines. This implies that the core of the line is very sensitive to the way scattering is treated in the modeling of the line, and in turn this allows us to separate the uncertain effects in the atomic data (viz. the Van der Waals broadening) from the uncertainties in the underlying atmospheric model. The main objective of this research was to compare observed spectra of the magnesium resonance lines against theoretical line profiles calculated from recent models of the solar atmosphere by Fontenla et al., hereafter called FAL. These models extend earlier solar models by Avrett et al. and provide yet a better match between observed and calculated spectra, in particular of the hydrogen Lyman line, by including the effects of ambipolar diffusion of hydrogen and helium in the mass and energy transport equations for the chromosphere-corona transition region. Although the transition region models can be built self-consistently (without additional observational input) once the proper boundary conditions are given at the top of the chromosphere, there are uncertainties in the models precisely because these boundary conditions are not well known. With the analysis of UVSP magnesium observations we strive to constrain these uncertainties.

Avrett, Eugene H.↗

NASA NTRS2012

Analytical Round Robin for Elastic-Plastic Analysis of Surface Cracked Plates: Phase I Results

An analytical round robin for the elastic-plastic analysis of surface cracks in flat plates was conducted with 15 participants. Experimental results from a surface crack tension test in 2219-T8 aluminum plate provided the basis for the inter-laboratory study (ILS). The study proceeded in a blind fashion given that the analysis methodology was not specified to the participants, and key experimental results were withheld. This approach allowed the ILS to serve as a current measure of the state of the art for elastic-plastic fracture mechanics analysis. The analytical results and the associated methodologies were collected for comparison, and sources of variability were studied and isolated. The results of the study revealed that the J-integral analysis methodology using the domain integral method is robust, providing reliable J-integral values without being overly sensitive to modeling details. General modeling choices such as analysis code, model size (mesh density), crack tip meshing, or boundary conditions, were not found to be sources of significant variability. For analyses controlled only by far-field boundary conditions, the greatest source of variability in the J-integral assessment is introduced through the constitutive model. This variability can be substantially reduced by using crack mouth opening displacements to anchor the assessment. Conclusions provide recommendations for analysis standardization.

Wells, D. N.↗

DOE OSTI2021

AFM Slip Length Measurements for Water at Selected Phyllosilicate Surfaces

Most reported slip length measurements have been made at the surfaces of synthetic materials and modified synthetic materials. In contrast, few slip length measurements at the surface of unmodified natural mineral surfaces have been reported. In this regard, flow at the silica face surfaces of the phyllosilicate minerals, talc and mica, was considered. A slip boundary condition was expected at the nonpolar hydrophobic silica surface of talc leading to enhanced flow, and a no-slip boundary condition was expected at the hydrophilic silica surface of mica. Atomic force microscopy (AFM) slip length measurements were made at the talc and mica surfaces. The slip length results for the hydrophobic silica surface of talc were contrasted to the results for the hydrophilic silica surface of mica (no-slip flow). The results are discussed based on molecular dynamics simulations (MDS), as reported in the literature, and AFM images of surface nanobubbles. For nonpolar hydrophobic surfaces (such as talc), it is doubtful that the MDS interfacial water structure and the water exclusion zone (3.2 Å) account for the AFM slip flow with slip lengths as great as 95 nm. Rather, a better explanation for the AFM slip flow condition is based on reduced interfacial viscosity due to the presence of dissolved gas and the accommodation of pancake nanobubbles at the talc surface having a height dimension of magnitude similar to the slip length.

42 ENGINEERING↗

DOE OSTI2025

Exact spectral gaps of random one-dimensional quantum circuits

The spectral gap of local random quantum circuits is a fundamental property that determines how close the moments of the circuit's unitaries match those of a Haar random distribution. When studying spectral gaps, it is common to bound these quantities using tools from statistical mechanics or via quantum information-based inequalities. Here, by focusing on the second moment of one-dimensional unitary circuits where nearest-neighboring gates act on sets of qudits (with open and closed boundary conditions), we show that one can exactly compute the associated spectral gaps. Indeed, having access to their functional form allows us to prove several important results, such as the fact that the spectral gap for closed boundary condition is exactly the square of the gap for open boundaries, as well as improve on previously known bounds for approximate design convergence. Finally, we verify our theoretical results by numerically computing the spectral gap for systems of up to 70 qubits, as well as comparing them to gaps of random orthogonal and symplectic circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1981

Measurements in ducted flows by laser-Doppler anemometry

The merits of small-scale, water-flow experiments with measurements by laser-Doppler anemometry are indicated in relation to previously reported results. New measurements are reported in a square-to-round transition duct and show that, for example, the pressure forces overwhelm any small normal-stress driven secondary flows in the inlet plane to given maximum values of around 0.07 of the bulk velocity in the exit plane. The precision of manufacture of the ducts means that the geometric boundary conditions are accurately known and the results are therefore suitable for calculation purposes. Similarly, the instrumentation allows the flow boundary conditions to be accurately determined.

Taylor, A. M. K. P.↗

NASA NTRS1989

Time domain numerical calculations of unsteady vortical flows about a flat plate airfoil

A time domain numerical scheme is developed to solve for the unsteady flow about a flat plate airfoil due to imposed upstream, small amplitude, transverse velocity perturbations. The governing equation for the resulting unsteady potential is a homogeneous, constant coefficient, convective wave equation. Accurate solution of the problem requires the development of approximate boundary conditions which correctly model the physics of the unsteady flow in the far field. A uniformly valid far field boundary condition is developed, and numerical results are presented using this condition. The stability of the scheme is discussed, and the stability restriction for the scheme is established as a function of the Mach number. Finally, comparisons are made with the frequency domain calculation by Scott and Atassi, and the relative strengths and weaknesses of each approach are assessed.

Hariharan, S. I.↗

NASA NTRS1991

Sensing and actuation in a liquid propellant rocket engine

Equations describing the conservation of mass, momentum, and energy, and multicomponent transport and chemical reaction in a liquid propellant rocket (LPR) engine are given, along with appropriate boundary conditions. The formulation includes a detailed description of the two-phase (gas/liquid) flow in terms of a moving boundary problem, with droplets and vapor occupying time-dependent subdomains of the LPR chamber. Several sensor and actuator types are considered in detail, and their mathematical descriptions are incorporated into the differential equations and boundary conditions.

Bentsman, Joseph↗

NASA NTRS1992

Time domain numerical calculations of unsteady vortical flows about a flat plate airfoil

Hariharan, S. I.↗