Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic boundaries”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Flutter study of an advanced composite wing with external stores

A flutter test using a scaled model of an advanced composite wing for a Navy attack aircraft has been conducted in the NASA Langley Research Center Transonic Dynamics Tunnel. The model was a wall-mounted half-span wing with a semi-span of 6.63 ft. The wing had an aspect ratio of 5.31, taper ratio of 0.312, and quarter-chord sweep of 25 degrees. The model was supported in a manner that simulated the load path in the carry-through structure of the aircraft and the symmetric boundary condition at the fuselage centerline. The model was capable of carrying external stores from three pylon locations on the wing. Flutter tests were conducted for the wing with and without external stores. No flutter was encountered for the clean wing at test conditions which simulated the scaled airplane operating envelope. Flutter boundaries were obtained for several external store configurations. The flutter boundaries for the fuel tanks were nearly Mach number independent (occurring at constant dynamic pressure). To study the aerodynamic effect of the fuel tank stores, pencil stores (slender cylindrical rods) which had the same mass and pitch and yaw inertia as the fuel tanks were tested on the model. These pencil store configurations exhibited a transonic dip in the flutter dynamic pressure, indicating that the aerodynamic effect of the actual fuel tanks on flutter was significant. Several flutter analyses methods were used in an attempt to predict the flutter phenomenon exhibited during the wind-tunnel test. The analysis gave satisfactory predictions of flutter for the pencil store configurations, but unsatisfactory correlation for the actual fuel tank configurations.

Cole, Stanley R.↗
NASA NTRS1995

A numerical model for dynamic wave rotor analysis

A numerical model has been developed which can predict the dynamic (and steady state) performance of a wave rotor, given the geometry and time dependent boundary conditions. The one-dimensional, perfect gas, CFD based code tracks the gasdynamics in each of the wave rotor passages as they rotate past the various ducts. The model can operate both on and off-design, allowing dynamic behavior to be studied throughout the operating range of the wave rotor. The model accounts for several major loss mechanisms including finite passage opening time, fluid friction, heat transfer to and from the passage walls, and leakage to and from the passage ends. In addition, it can calculate the amount of work transferred to and from the fluid when the flow in the ducts is not aligned with the passages such as occurs in off-design operation. Since it is one-dimensional, the model runs reasonably fast on a typical workstation. This paper will describe the model and present the results of some transient calculations for a conceptual four port wave rotor designed as a topping cycle for a small gas turbine engine.

Paxson, D. E.↗
NASA NTRS2008

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

Interaction of solar wind with the magnetopause-boundary layer and generation of magnetic impulse events

The transport of mass, momentum, energy and waves from the solar wind to the Earth's magnetosphere takes place in the magnetopause-boundary layer region. Various plasma processes that may occur in this region have been proposed and studied. In this paper, we present a brief review of the plasma processes in the dayside magnetopause-boundary layer. These processes include (1) flux transfer events at the dayside magnetopause, (2) formation of plasma vortices in the low-latitude boundary layer by the Kelvin-Helmholtz instability and coupling to the polar ionosphere, (3) the response of the magnetopause to the solar wind dynamic pressure pulses, and (4) the impulsive penetration of solar wind plasma filaments through the dayside magnetopause into the magnetospheric boundary layer. Through the coupling of the magnetopause-boundary layer to the polar ionosphere, those above processes may lead to occurrence of magnetic impulse events observed in the high-latitude stations.

Lee, L. C.↗
NASA NTRSDate not supplied

Validation of Cryogenic Propellant Tank Self-Pressurization by Leveraging Reduced Order Modeling within Computational Fluid Dynamics Simulation

Validation of cryogenic propellant tank self-pressurization was performed using a hybrid Computational Fluid Dynamics (CFD) and reduced order modeling methodology. Data from a liquid hydrogen ground test conducted at the K-site facility at the National Aeronautics and Space Administration (NASA) Glenn Research Center was used for the validation effort. Liquid phase dynamics were explicitly resolved with a CFD tool. Vapor phase dynamics were modeled as a point mass that communicated heat from the tank wall to the liquid phase via a boundary condition used at the gas-liquid interface. The method proved to be more accurate, robust, and efficient than explicit resolution of the dynamics using a standard Volume of Fluid (VOF) methodology. The subject pressurization process was found to be heavily dependent upon both the relatively high liquid temperature gradient near the gas-liquid interface and the natural convection flow path. Modeling the gas-liquid interface as an immovable surface eliminated temperature gradient destroying gas-liquid interface velocities observed in VOF simulations, and correspondingly enabled more rapid simulation since interface advection was not allowed. The single phase computational domain also facilitated the ability to demonstrate spatial resolution convergence of natural convection cells within the liquid which significantly impacted the tank pressurization rate. This work was used to demonstrate the critical physics for tank self-pressurization and numerical methodologies that may be used to best resolve those physics. The findings informed development and operation of production level CFD tools used in the Fluid Dynamics Branch at NASA Marshall Space Flight Center.

J. M. Brodnick↗
NASA NTRS1995

Eliminating Flow Separation and Reducing Viscous Drag Through Boundary Layer Analysis and Manipulation

As both computers and flow-analyzing equations have increased in sophistication, Computational Fluid Dynamics (CFD) has evolved into a fixture for advanced aircraft design. While CFD codes have improved in accuracy and efficiency, their ability to encompass viscous effects is lacking in certain areas. For example, current CFD codes cannot accurately predict or correct for the increased drag due to these viscous effects at some flow conditions. However, by analyzing an airfoil's turbulent boundary layer, one can predict not only flow separation via the shape factor parameter, but also viscous drag via the momentum thickness. Various codes have been written which can calculate turbulent boundary layer parameters. The goal of my research is to develop procedures for modifying an airfoil (via its local pressure distribution) to eliminate boundary layer separation and/or to reduce viscous drag. The modifications to the local pressure distribution necessary to achieve these objectives will be determined using a direct-iterative method installed into a turbulent boundary layer analyzer. Furthermore, the modifications should preserve the basic characteristics of the original airfoil.

Oser, Matthew D.↗
DOE OSTI2023

California's harvested wood products: A time-dependent assessment of life cycle greenhouse gas emissions

Following life-cycle assessment (LCA) methodology, this study presents a state-level estimation of embodied carbon of wood products harvested in 2019 from California and subsequently processed, manufactured, transported, used, and disposed at the end-of-life (EoL). In a conventional static approach to LCA, all GHG emissions were aggregated and considered to occur at year 0 of the given time horizon (500 years in this study) and used a static characterization factor (CF). In dynamic LCA, GHG emissions occurring in different years were considered, and their global warming impact (GWI) was determined using a time-dependent CF over the selected time horizon of 500 years. Four scenarios were developed to examine the impact of EoL choices on GWI. It was found that dynamic GWI for all scenarios ranged from 0.27 to 0.93 million tonne CO₂e, which were 45–73 % lower than those estimated with static LCA approach, indicating that the static LCA approach could lead to an underestimation of the benefits of substituting wood for non-wood products, compared to those based on dynamic LCA approach. This analysis also demonstrated that the choice of EoL treatment option is a key factor affecting the estimated GWI as it directly determines the annual emission of GHGs released into atmosphere and subsequently their warming effect depending on the time harvested wood products (HWPs) spend in the horizon of assessment. Altogether, the dynamic LCA performed in this study enabled more robust interpretations of embodied carbon by including temporal boundaries associated with the HWPs life cycle.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1996

Artificial Boundary Conditions Based on the Difference Potentials Method

While numerically solving a problem initially formulated on an unbounded domain, one typically truncates this domain, which necessitates setting the artificial boundary conditions (ABC's) at the newly formed external boundary. The issue of setting the ABC's appears to be most significant in many areas of scientific computing, for example, in problems originating from acoustics, electrodynamics, solid mechanics, and fluid dynamics. In particular, in computational fluid dynamics (where external problems present a wide class of practically important formulations) the proper treatment of external boundaries may have a profound impact on the overall quality and performance of numerical algorithms. Most of the currently used techniques for setting the ABC's can basically be classified into two groups. The methods from the first group (global ABC's) usually provide high accuracy and robustness of the numerical procedure but often appear to be fairly cumbersome and (computationally) expensive. The methods from the second group (local ABC's) are, as a rule, algorithmically simple, numerically cheap, and geometrically universal; however, they usually lack accuracy of computations. In this paper we first present a survey and provide a comparative assessment of different existing methods for constructing the ABC's. Then, we describe a relatively new ABC's technique of ours and review the corresponding results. This new technique, in our opinion, is currently one of the most promising in the field. It enables one to construct such ABC's that combine the advantages relevant to the two aforementioned classes of existing methods. Our approach is based on application of the difference potentials method attributable to V. S. Ryaben'kii. This approach allows us to obtain highly accurate ABC's in the form of certain (nonlocal) boundary operator equations. The operators involved are analogous to the pseudodifferential boundary projections first introduced by A. P. Calderon and then also studied by R. T. Seeley. The apparatus of the boundary pseudodifferential equations, which has formerly been used mostly in the qualitative theory of integral equations and PDE'S, is now effectively employed for developing numerical methods in the different fields of scientific computing.

Tsynkov, Semyon V.↗
DOE OSTI2020

< c+a > dislocations in {101 - 2} twins in Mg: A kinematic and energetic requirement

Consideration of published TEM and molecular dynamics studies of interactions between matrix dislocations and {101¯2} twin boundaries (TB) in hexagonal close packed metals inspires a hypothesis regarding the formation of dislocations with < c+a > Burgers vectors inside the twins. It was previously reported that when basal < a > dislocations from the matrix are swept by an advancing TB, they are converted into partial dislocations bounding basal I 1 stacking faults (SFs) within the twin. In the present work, we show that this dislocation configuration near TB renders it energetically favorable for a pair of partial dislocations to constrict into a dislocation within the twin. New in situ TEM observations confirm this hypothesis by revealing long < c+a > dislocation dipoles terminated by the TB at one end and by I1 SFs within the twin at the other end. Finally, the resulting near-screw < c+a > dislocation dipoles are observed to be glissile, enabling self-annihilation in some circumstances. Thus, rather than requiring a high energy dislocation pile-up to press two < a > dislocations through a TB into a single < c+a > dislocation, it is demonstrated that < c+a > dislocations naturally form, within {101¯2} twins in Mg, as TB migrates into a matrix, which inevitably contains dislocations with < a > Burgers vectors. Furthermore, these observations help to explain how the twins themselves may subsequently deform and why they may be particularly vulnerable to plastic, shear localization.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Towards understanding particle rigid-body motion during solid-state sintering

A quantitative understanding of particle rigid body (RB) motion that inherently accompanies grain boundary (GB) diffusion is highly desirable to understand and control the dynamic interplay between coarsening and densification during solid state sintering. In this work, by computer simulation using a multi-phase-field approach, we analyze systematically the roles played by each of these processes at different stages of the shrinkage of the internal pore in a three-particle green body as a function of particle size as well as thermodynamic and kinetic factors of interfaces. We demonstrate that particle RB translation promotes both neck growth, and pore rounding and shrinkage. Moreover, the forces acting at GBs and pulling neighboring particles towards one another dynamically evolve as particles fuse. In contrast, particle RB rotation has no contribution to pore shrinkage. The translational force acting on an individual particle varies with not only its size, but also the number and sizes of its neighboring particles.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Improving Stratocumulus Cloud Amounts in a 200‐m Resolution Multi‐Scale Modeling Framework Through Tuning of Its Interior Physics

Abstract High‐Resolution Multi‐scale Modeling Frameworks (HR)—global climate models that embed separate, convection‐resolving models with high enough resolution to resolve boundary layer eddies—have exciting potential for investigating low cloud feedback dynamics due to reduced parameterization and ability for multidecadal throughput on modern computing hardware. However low clouds in past HR have suffered a stubborn problem of over‐entrainment due to an uncontrolled source of mixing across the marine subtropical inversion manifesting as stratocumulus dim biases in present‐day climate, limiting their scientific utility. We report new results showing that this over‐entrainment can be partly offset by using hyperviscosity and cloud droplet sedimentation. Hyperviscosity damps small‐scale momentum fluctuations associated with the formulation of the momentum solver of the embedded large eddy simulation. By considering the sedimentation process adjacent to default one‐moment microphysics in HR, condensed phase particles can be removed from the entrainment zone, which further reduces entrainment efficiency. The result is an HR that can produce more low clouds with a higher liquid water path and a reduced stratocumulus dim bias. Associated improvements in the explicitly simulated sub‐cloud eddy spectrum are observed. We report these sensitivities in multi‐week tests and then explore their operational potential alongside microphysical retuning in decadal simulations at operational 1.5° exterior resolution. The result is a new HR having desired improvements in the baseline present‐day low cloud climatology, and a reduced global mean bias and root mean squared error of absorbed shortwave radiation. We suggest it should be promising for examining low cloud feedbacks with minimal approximation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Preface for frontiers of magnetic reconnection research in heliophysical, astrophysical, and laboratory plasmas

Magnetic reconnection—the topological rearrangement of magnetic field—underlies many explosive phenomena across a wide range of natural and laboratory plasmas.3 It plays a pivotal role in electron and ion heating, particle acceleration to high energies, energy transport, and self-organization. Reconnection can have a complex relationship with turbulence at both large and small scales, leading to various effects that are only beginning to be understood. In heliophysics, magnetic reconnection plays a key role in solar flares, coronal mass ejections, coronal heating, solar wind dissipation, the interaction of interplanetary plasma with magnetospheres, dynamics of planetary magnetospheres such as magnetic substorms, and the heliospheric boundary with the interstellar medium. Here, the magnetic reconnection is integral to the solar and planetary dynamo processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Vortex splitting in two-dimensional fluids and non-neutral electron plasmas with smooth vorticity profiles

Initially, elliptical, quasi-two-dimensional (2D) fluid vortices can split into multiple pieces if the aspect ratio is sufficiently large due to the growth and saturation of perturbations known as Love modes on the vortex edge. Presented here are experiments and numerical simulations, showing that the aspect ratio threshold for vortex splitting is significantly higher for vortices with realistic, smooth edges than that predicted by a simple “vortex patch” model, where the vorticity is treated as piecewise constant inside a deformable boundary. The experiments are conducted by exploiting the E × B drift dynamics of collisionless, pure electron plasmas in a Penning–Malmberg trap, which closely model 2D vortex dynamics due to an isomorphism between the Drift–Poisson equations describing the plasmas and the Euler equations describing ideal fluids. The simulations use a particle-in-cell method to model the evolution of a set of point vortices. The aspect ratio splitting threshold ranges up to about twice as large as the vortex patch prediction and depends on the edge vorticity gradient. This is thought to be due to spatial Landau damping, which decreases the vortex aspect ratio over time and, thus, stabilizes the Love modes. Near the threshold, asymmetric splitting events are observed in which one of the split products contains much less circulation than the other. These results are relevant to a wide range of quasi-2D fluid systems, including geophysical fluids, astrophysical disks, and drift-wave eddies in tokamak plasmas.

Physics↗
DOE OSTI2021

Phase coexistence at the first-order Mott transition revealed by pressure-dependent dielectric spectroscopy of κ - (BEDT-TTF) 2 - Cu 2 (CN) 3

The dimer Mott insulator kappa - (BEDT-TTF) 2 - Cu 2 (CN) 3 can be tuned into metallic and superconducting states on applying pressure of 1.5 kbar and more. We have performed dielectric measurements (7.5 kHz to 5 MHz) on kappa - (BEDT-TTF) 2 - Cu 2 (CN) 3 single crystals as a function of temperature (down to T = 8 K) and pressure (up to p = 4.3 kbar). In addition to the relaxor-like dielectric behavior seen below 50 K at p = 0, that moves toward lower temperatures with pressure, a second peak emerges in ε 1 (T) around T = 15 K. When approaching the insulator-metal boundary, this peak diverges rapidly reaching ε 1 ≈ 10 5 . Our dynamical mean-field theory calculations substantiate that the dielectric catastrophe at the Mott transition is not caused by closing the energy gap, but due to the spatial coexistence of correlated metallic and insulating regions. We discuss the percolative nature of the first-order Mott insulator-to-metal transition in all details.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Origin of micrometer-scale dislocation motion during hydrogen desorption

Hydrogen, while being a potential energy solution, creates arguably the most important embrittlement problem in high-strength metals. However, the underlying hydrogen-defect interactions leading to embrittlement are challenging to unravel. Here, we investigate an intriguing hydrogen effect to shed more light on these interactions. By designing an in situ electron channeling contrast imaging experiment of samples under no external stresses, we show that dislocations (atomic-scale line defects) can move distances reaching 1.5 μm during hydrogen desorption. Combining molecular dynamics and grand canonical Monte Carlo simulations, we reveal that grain boundary hydrogen segregation can cause the required long-range resolved shear stresses, as well as short-range atomic stress fluctuations. Thus, such segregation effects should be considered widely in hydrogen research.

36 MATERIALS SCIENCE↗
DOE OSTI2020

( Invited ) Ruggedness of SiC and GaN Power Transistors in Switching Based Tests

Currently, a spectrum of reliability tests are being performed by both device manufacturers and end users. Many of these qualification tests are operated within the device safe-operating-area (SOA), however, the devices often undergo dynamic events in many converter applications that can exceed the SOA boundaries. This paper presents our recent work on quantifying the out-of-SOA robustness of GaN power high-electron-mobility transistors (HEMTs) and SiC power MOSFETs by switching based tests. The unclamped inductive switching stresses are applied to p-gate GaN HEMTs to evaluate their surge-energy and overvoltage robustness. The continuous, hard-switching, turn-off stresses are applied to SiC MOSFETs to evaluate their overvoltage switching robustness. As a result, both experiments demonstrate the robustness of SiC and GaN devices under out-of-SOA switching events, and the physics of degradation and failure are understood.

42 ENGINEERING↗
DOE OSTI2023

The saturation mechanism of thermal instability

The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing the recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates; whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change the sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of condensation. For isobaric evolution, this “pressure reversal” occurs nearly simultaneously for every fluid element in condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies the gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave mediated. We show that both the contraction rebound event and subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal non-equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1972

Space shuttle booster multi-engine base flow analysis

A comprehensive review of currently available techniques pertinent to several prominent aspects of the base thermal problem of the space shuttle booster is given along with a brief review of experimental results. A tractable engineering analysis, capable of predicting the power-on base pressure, base heating, and other base thermal environmental conditions, such as base gas temperature, is presented and used for an analysis of various space shuttle booster configurations. The analysis consists of a rational combination of theoretical treatments of the prominent flow interaction phenomena in the base region. These theories consider jet mixing, plume flow, axisymmetric flow effects, base injection, recirculating flow dynamics, and various modes of heat transfer. Such effects as initial boundary layer expansion at the nozzle lip, reattachment, recompression, choked vent flow, and nonisoenergetic mixing processes are included in the analysis. A unified method was developed and programmed to numerically obtain compatible solutions for the various flow field components in both flight and ground test conditions. Preliminary prediction for a 12-engine space shuttle booster base thermal environment was obtained for a typical trajectory history. Theoretical predictions were also obtained for some clustered-engine experimental conditions. Results indicate good agreement between the data and theoretical predicitons.

Tang, H. H.↗