Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “embedded boundary”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Embedded function methods for compressible high speed turbulent flow

Fundamental issues relating to compressible turbulent flow are addressed. The focus has been on developing methods and testing concepts for attached flows rather than trying to force a conventional law of the wall into a zone of backflow. Although the dynamics of the near-wall flow in an attached turbulent boundary layer are relatively well documented, the dynamical features of a zone of reversed turbulent flow are not, nor are they well understood. Incompressibility introduces effects and issues that have been dealt with only marginally in the literature, therefore, the present work has been focussed on attached high-speed flows. The wall function method has been extended up through the supersonic to hypersonic speeds. Algorithms have been successfully introduced into the code that calculates the flow all the way to the wall, and testing is being carried out for progressively more complex flow situations.

Walker, J. D. A.↗
NASA NTRS1991

Kinematic domain decomposition for boundary-motion-induced flow simulations

A method is developed to solve the unsteady Navier-Stokes equations on a composite grid, which consists of subdomain grids moving with respect to each other. These subdomains are structured grids with different topologies. This method eliminates assuming the moving components to be instantaneously stationary, where deciding on the particular frozen instants is difficult and affects the solution adversely. Moreover, this method captures the boundary-motion-induced flow component. The method is demonstrated through a transonic flow past an airfoil, which experiences a combined motion of pitching and plunging. An O-grid around the airfoil is overlapped on a fine Cartesian grid, which is zonally embedded in a coarse Cartesian grid. The coarse grid is stationary but the other two grids are plunging. Only the O-grid is also sinusoidally pitching. The results are compared successfully with the experimental data.

Baysal, Oktay↗
NASA NTRS1994

Coupled electromechanical response of composite beams with embedded piezoelectric sensors and actuators

Unified mechanics are developed with the capability to model both sensory and active composite laminates with embedded piezoelectric layers. A discrete-layer formulation enables analysis of both global and local electromechanical response. The mechanics include the contributions from elastic, piezoelectric, and dielectric components. The incorporation of electric potential into the state variables permits representation of general electromechanical boundary conditions. Approximate finite element solutions for the static and free-vibration analysis of beams are presented. Applications on composite beams demonstrate the capability to represent either sensory or active structures and to model the complicated stress-strain fields, the interactions between passive/active layers, interfacial phenomena between sensors and composite plies, and critical damage modes in the material. The capability to predict the dynamic characteristics under various electrical boundary conditions is also demonstrated.

Saravanos, D. A.↗
NASA NTRS1995

The decay of longitudinal vortices shed from airfoil vortex generators

An experimental study is conducted to examine the crossplane structure and streamwise decay of vortices shed from airfoil-type vortex generators. The vortex generators are set in a counter-rotating array spanning the full circumference of a straight pipe. The span of the vortex generators above the duct surface, h, is approximately equal to the local turbulent boundary layer thickness, delta. Measurement of three-component mean flow velocity in downstream crossplanes are used to characterize the structure of the shed vortices. Measurements in adjacent crossplanes (closely spaced along the streamwise coordinate) characterize the interaction and decay of the embedded vortices. A model constructed by the superposition of Oseen vortices is compared to the data for one test case.

Wendt, Bruce J.↗
NASA NTRS1994

Auroral activity associated with Kelvin-Helmholtz instability at the inner edge of the low-latitude boundary layer

Auroral activity occurred in the late afternoon sector (approx. 16 MLT) in the northern hemisphere during the passage at Earth of an interplanetary magnetic cloud on January 14, 1988. The auroral activity consisted of a very dynamic display which was preceded and followed by quiet auroral displays. During the quiet displays, discrete rayed arcs aligned along the geomagnetic L shells were observed. In the active stage, rapidly evolving spiral forms centered on magnetic zenith were evident. The activity persisted for many minutes and was characterized by the absence of directed motion. They were strongly suggestive of intense filaments of upward field-aligned currents embedded in the large-scale region 1 current system. Distortions of the flux ropes as they connect from the equatorial magnetosphere to the ionosphere were witnessed. We assess as possible generating mechanisms three nonlocal sources known to be associated with field-aligned currents. Of these, partial compressions of the magnetosphere due to variations of solar wind dynamic pressure seem an unlikely source. The possibility that the auroral forms are due to reconnection is investigated but is excluded because the active aurora were observed on the closed field line region just equatorward of the convection reversal boundary. To support this conclusion further, we apply recent results on the mapping of ionospheric regions to the equatorial plane based on the Tsyganenko 1989 model (Kaufmann et al., 1993). We find that for comparable magnetic activity the aurora map to the equatorial plane at X(sub GSM) = approx. 3 R(sub E) and approx. 2 R(sub E) inward of the magnetopause, that is, the inner edge of the boundary layer close to dusk. Since the auroral forms are manifestly associated with magnetic field shear, a vortical motion at the equatorial end of the flux rope is indicated, making the Kelvin-Helmholtz instability acting at the inner edge of the low-latitude boundary layer the most probable generating source.

Farrugia, C. J.↗
NASA NTRS1993

Dynamic characteristics of specialty composite structures with embedded damping layers

Damping mechanics for simulating the damped dynamic characteristics in specialty composite structures with compliant interlaminar damping layers are presented. Finite-element based mechanics incorporating a discrete layer (or layer-wise) laminate damping theory are utilized to represent general laminate configurations in terms of lay-up and fiber orientation angles, cross-sectional thickness, shape, and boundary conditions. Evaluations of the method with exact solutions and experimental data illustrate the accuracy of the method. Additional applications investigate the potential for significant damping enhancement in angle-ply composite laminates with cocured interlaminar damping layers.

Saravanos, D. A.↗
DOE OSTI2022

Creep anisotropy modeling and uncertainty quantification of an additively manufactured Ni-based superalloy

The advantages offered by additive manufacturing over traditional processes has driven a great deal of industrial and academic interest in recent years. However, the process is relatively new and requires additional investigation to become sufficiently mature for wide scale industrial adoption. Electron beam melting powder bed fusion is one technology that has shown promise for fabricating high temperature resistant materials such as nickel based superalloys. The resulting microstructures typically exhibit a strong fiber texture in the build direction giving rise to anisotropic time-dependent deformation behavior. In order to accelerate the qualification of these materials for industrial adoption accurate numerical models are needed for simulating their behavior. In this work a crystal plasticity model including non-Schmid effects is presented for capturing creep anisotropy observed in additively manufactured IN738LC. The model is calibrated via a probabilistic framework where model parameters are treated as random variables. An iterative sequential design strategy is utilized to efficiently identify the probability density of the unknown model parameters. As a case study the model is utilized to investigate the behavior of randomly oriented equiaxed grain clusters sometimes observed embedded in the additively manufactured columnar structure. A synthetic realization is simulated and uncertainty is propagated through to the full-field response. Results indicate that these features are the source of significant creep relaxation and strain accumulation which partially explains observed grain boundary decohesion at these locations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Single- and two-particle finite size effects in interacting lattice systems

Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with nonlocal correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles.

36 MATERIALS SCIENCE↗
NASA NTRS1985

A 3-D chimera grid embedding technique

A three-dimensional (3-D) chimera grid-embedding technique is described. The technique simplifies the construction of computational grids about complex geometries. The method subdivides the physical domain into regions which can accommodate easily generated grids. Communication among the grids is accomplished by interpolation of the dependent variables at grid boundaries. The procedures for constructing the composite mesh and the associated data structures are described. The method is demonstrated by solution of the Euler equations for the transonic flow about a wing/body, wing/body/tail, and a configuration of three ellipsoidal bodies.

Benek, J. A.↗
DOE OSTI2022

Design of a Combined Proportional Integral Derivative Controller to Regulate the Temperature Inside a High-Temperature Tubular Solar Reactor

Solar fuels are proven to be promising candidates for thermochemical energy storage. However, the transient nature of solar radiation is an obstacle to maintaining a stable operational temperature inside a solar reactor. To overcome this challenge, the temperature of a solar reactor can be regulated by controlling the incoming solar radiation or the feedstock flowrate inside the reactor. In this work, a combined proportional integral derivative (PID) controller is implemented to regulate the temperature inside a high-temperature tubular solar reactor with counter-current flowing gas/particles. The control model incorporates two control systems to regulate incoming solar radiation and gas flow simultaneously. The design of the controller is based on a reduced-order numerical model of a high-temperature tubular solar reactor that is vertically oriented with an upward gas flow and downward particle flow. The reactor receives heat circumferentially through its wall over a finite segment of its length. Formulation of the heat transfer model is presented by applying the energy balance for the reactor tube and considering heat and mass transfer inside. A set of governing differential equations are solved numerically by using the finite volume method to obtain reactor wall, particles, and gas temperatures along the reactor length with various boundary conditions. Simulation results are used to tune the PID controller parameters by utilizing the Ziegler–Nichols tuning method. Both the simulation results and the controller performance are visualized on the LABVIEW platform. The controller is challenged to track different temperature setpoints with different scenarios of transient solar radiation. Here, the performance of the PID controller was compared to experimental results obtained from an industrial PID controller embedded in a 7 kW electric furnace. Results show that the combined PID controller is successful in maintaining a stable temperature inside the reactor by regulating the incoming solar radiation and the flowrate via small steady-state error and reasonable settling time and overshoot.

14 SOLAR ENERGY↗
NASA NTRS1987

Observations of large-scale Birkeland currents with Viking

Data from the Viking Magnetic Field Experiment and Hot Plasma Experiment acquired on March 25, 1986 are analyzed. The large-scale region 1 and region 2 Birkeland current systems (with densities of about 0.1 microA/sq m) were readily identified in the morning sector at altitudes of 12,500 to 13,000 km. Embedded in the large-scale region 1 system are small-scale intense currents that are not necessarily aligned in the east-west geomagnetic direction. The major carriers of the large-scale region 1 Birkeland currents appear to be upward-flowing electrons with energies below 10 eV. Embedded within this flux can be spatially limited regions of 'hotter' electron fluxes that carry small-scale currents. The correlation of magnetic and particle observations obtained here supports previous suggestions that the region 1 currents in the morning sector have their source in the boundary layer and that the region 2 currents are associated with the plasma sheet.

Potemra, T. A.↗
NASA NTRS1993

A Novel Optical/digital Processing System for Pattern Recognition

This paper describes two processing algorithms that can be implemented optically: the Radon transform and angular correlation. These two algorithms can be combined in one optical processor to extract all the basic geometric and amplitude features from objects embedded in video imagery. We show that the internal amplitude structure of objects is recovered by the Radon transform, which is a well-known result, but, in addition, we show simulation results that calculate angular correlation, a simple but unique algorithm that extracts object boundaries from suitably threshold images from which length, width, area, aspect ratio, and orientation can be derived. In addition to circumventing scale and rotation distortions, these simulations indicate that the features derived from the angular correlation algorithm are relatively insensitive to tracking shifts and image noise. Some optical architecture concepts, including one based on micro-optical lenslet arrays, have been developed to implement these algorithms. Simulation test and evaluation using simple synthetic object data will be described, including results of a study that uses object boundaries (derivable from angular correlation) to classify simple objects using a neural network.

Boone, Bradley G.↗
NASA NTRS2017

Design of an Axisymmetric Afterbody Test Case for CFD Validation

As identified in the CFD Vision 2030 Study commissioned by NASA, validation of advanced RANS models and scale-resolving methods for computing turbulent flow fields must be supported by continuous improvements in fundamental, high-fidelity experiments designed specifically for CFD implementation. In accordance with this effort, the underpinnings of a new test platform referred to herein as the NASA Axisymmetric Afterbody are presented. The devised body-of-revolution is a modular platform consisting of a forebody section and afterbody section, allowing for a range of flow behaviors to be studied on interchangeable afterbody geometries. A body-of-revolution offers advantages in shape definition and fabrication, in avoiding direct contact with wind tunnel sidewalls, and in tail-sting integration to facilitate access to higher Reynolds number tunnels. The current work is focused on validation of smooth-body turbulent flow separation, for which a six-parameter body has been developed. A priori RANS computations are reported for a risk-reduction test configuration in order to demonstrate critical variation among turbulence model results for a given afterbody, ranging from barely-attached to mild separated flow. RANS studies of the effects of forebody nose (with/without) and wind tunnel boundary (slip/no-slip) on the selected afterbody are presented. Representative modeling issues that can be explored with this configuration are the effect of higher Reynolds number on separation behavior, flow physics of the progression from attached to increasingly-separated afterbody flows, and the effect of embedded longitudinal vortices on turbulence structure.

Disotell, Kevin J.↗
NASA NTRS1981

Trailing edge noise prediction from measured surface pressures

Trailing edge (TE) noise is investigated for the case of a two-dimensional airfoil embedded in a uniform low Mach number flow, and the usefulness of several TE noise theories is examined by applying them to the measured data. The TE noise spectra and directivity are quantitatively determined for the case of a high Reynolds number and a fully turbulent boundary layer. Parameters include angle of attack, flow velocity and TE bluntness. Evanescent wave theories by Chase (1975) and Chandiramani (1974) are compared to the present results and show good agreement. Agreement of the near field pressure scatter phenomenon analysis with measurements implies that the basic assumptions used in the analysis are correct, i.e., the turbulent boundary layer (TBL) flow passes the trailing edge into the wake region. No hydrodynamic wake shedding activity is confirmed for the two-sided TBL flow, and a method incorporating the principles of the coherent output power method is used to determine the sound field. The near field edge scatter model is found to not only establish optimum sizing of edge treatment for noise control, but also to separate and identify the scattered field from the incident hydrodynamic field.

Brooks, T. F.↗
DOE OSTI2025

Temperature effects of ion irradiation on the nanostructural features in ductile-phase-toughened tungsten composites

Ductile-phase toughened tungsten (DPT W) composites have emerged as promising candidates for load-bearing components behind the plasma-facing tungsten armor in fusion reactors due to their enhanced thermomechanical properties. This study focuses on a composite consisting of W particles embedded in a ductile NiFeW solution matrix, hot-rolled to a thickness reduction of 87% (87R DPT W). Sequential irradiations with Ni2+ and He+ ions were performed to identical doses and helium concentrations at room temperature (RT) and 1273 K. Irradiation at RT produced no discernible nanostructural features due to the immobility of mono-vacancies, whereas cavity formation was observed at 973 K. At 1273 K, the W phase exhibited larger cavities, reduced cavity number density, and lower volumetric swelling compared to 973 K. Notably, nanosized NiFeW precipitates formed within the W phase at 1273 K, a phenomenon absent at 973 K. A new phase of cubic (NiFe)6W6C was also observed at the interphase boundary. In contrast, the NiFeW matrix showed no nanostructural changes at 1273 K, likely due to cavity dissociation. Separate irradiations at 1273 K indicated that Ni2+ ions induced precipitate formation in the W phase, while He+ ions exclusively caused cavity formation. The microstructure of 87R DPT W irradiated at RT and subsequently annealed at 1273 K closely resembled that of material irradiated directly at 1273 K. Like oxide-dispersion-strengthened steels, the observed nanoparticle-embedded W can inhibit dislocation propagation, potentially delaying the ductile-to-brittle transition temperature. These findings highlight the potential of NiFeW nanoparticle-reinforced W composites as irradiation-resistant materials for fusion reactors.

Jiang, Weilin (ORCID:0000000183028313)↗
DOE OSTI2022

Incorporating two-dimensional conduction modeling techniques into an energy simulation program: The EnergyPlus radiant system example

Due to the inherent complexity of multi-dimensional solutions for transient conduction in a multilayered building element, most heat balance-based algorithms assume that conduction occurs only in one dimension—through the building element between the inside and outside environments. Usually, since the length and width of the elements are much greater than the thickness of the element, this is a valid assumption. However, there are some cases such as hydronic radiant slabs where this assumption may not lead to an accurate estimate of the energy consumption of the system. In a hydronic radiant slab, water tubes are embedded in a slab, and hot or chilled water is circulated to meet any heating or cooling loads. The energy consumption of such a system will depend highly on the temperature of the water being sent to the slab (i.e., produced by the boiler or chiller). The temperature of the water that must be sent to the slab will depend on the spacing of the tubes, the multilayer construction of the building element, and the boundary conditions to which the construction is exposed. Thus, to accurately predict the energy consumption of such a system, a transient two-dimensional heat conduction model must be integrated with a heat balance-based simulation. This paper demonstrates the potential loss of accuracy associated with the one-dimensional radiant system models by providing a comparison of one- and two-dimensional solutions using EnergyPlus.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1983

The stability of coronal loops - Finite-length and pressure-profile limits

Results are described from a quickly converging, necessary-and-sufficient, MHD-stability test for coronal-loop models. The primary stabilizing influence arises from magnetic line tying at the photosphere, and this end conditions requires a series expansion of possible loop excitations. The stability boundary is shown to quickly approach a limit as the number of terms increases, providing a critical length for the loop in proportion to its transverse magnetic scale. Several models of force-free-field profiles are tested and the stability behavior of a localized current channel, embedded in an external current-free region, is shown to be superior to that of other, broader, current profiles. Pressure-gradient effects, leading to increased or decreased stability, are shown to be amplified by line tying. Long loops must either conduct low net current, or exhibit an axial-field reversal coexisting with a low-pressure core. The limits on stability depend on the magnetic aspect ratio, the plasma-to-magnetic pressure ratio, and the field orientation at the loop edge. Applications of these results to the structure of coronal loops are described.

Einaudi, G.↗
NASA NTRS1989

Low level remote sensing: Orographic winds

The orographic flow data set was obtained from a flight program to measure the influence of orographic features on turbulence momentum, heat, and moisture fluxes. The NASA B-57 aircraft instrumented with probes for measuring the three fluctuating wind speed components, temperature, and humidity was the primary measuring vehicle. Ancillary measurements were made with several ground-based sensors. These include the NOAA radar wind profilers, the Boulder wind network, the PROFS mesoscale surface network, the Boulder Atmospheric Observatory 300 m tower, special rawinsonde observations, and the NOAA/WPL Doppler lidar. The major objective of the flight program was to provide planetary boundary layer parameter information for new and current general circulation computer models. A numerical code, WINDER, based on a discrete element technique was run to numerically model the water tunnel simulated flow. Comparison of the analytical model with the experimental results is very good. Physical fluid dynamic principles embedded in the computational model and visual and hot wire anemometer measurements from the simulation are being rationalized to develop a physical understanding of the vortex flow. The results will be interpreted as they pertain to full scale atmospheric flows.

Frost, Walter↗