Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Interface velocity”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

On the design of stable, consistent, and conservative high-order methods for multi-material hydrodynamics

Obtaining stable and high-order numerical solutions for multi-material hydrodynamics is an open challenge. Although slope limiters are widely used to maintain monotonicity near discontinuities, typical limiting procedures violate closure laws at the discrete level when applied to multi-material hydrodynamics equations. Due to this, the high-order expansions of quantities related by the closure laws are no longer consistent. The commonly observed symptom of this consistency-violation is that the numerical method fails to maintain constant pressure and velocity across material interfaces. This leads to sub-optimal convergence rates for smooth multi-material problems as well. Specialized limiting procedures that satisfy consistency while maintaining conservation need to be developed for such equations. A novel procedure that re-instates consistency into slope-limited high-order discretizations applied to the multi-material hydrodynamics equations is presented here. Using simple examples, it is demonstrated that the presented method satisfies closure laws at the discrete level, while maintaining conservative properties of the high-order method. Furthermore, this procedure involves a projection step which relies on the compact basis of the underlying spatial discretization, i.e. for discontinuous schemes (viz. DG and FV) the projection is local, and does not involve global matrix solves. Comparisons with conventional approaches emphasizes the necessity of the consistent closure-law preserving limiting approach, in order to maintain design order of accuracy for smooth multi-material problems.

36 MATERIALS SCIENCE↗
NASA NTRS1976

The design of digital-adaptive controllers for VTOL aircraft

Design procedures for VTOL automatic control systems have been developed and are presented. Using linear-optimal estimation and control techniques as a starting point, digital-adaptive control laws have been designed for the VALT Research Aircraft, a tandem-rotor helicopter which is equipped for fully automatic flight in terminal area operations. These control laws are designed to interface with velocity-command and attitude-command guidance logic, which could be used in short-haul VTOL operations. Developments reported here include new algorithms for designing non-zero-set-point digital regulators, design procedures for rate-limited systems, and algorithms for dynamic control trim setting.

Stengel, R. F.↗
NASA NTRS1976

Digital controllers for VTOL aircraft

Using linear-optimal estimation and control techniques, digital-adaptive control laws have been designed for a tandem-rotor helicopter which is equipped for fully automatic flight in terminal area operations. Two distinct discrete-time control laws are designed to interface with velocity-command and attitude-command guidance logic, and each incorporates proportional-integral compensation for non-zero-set-point regulation, as well as reduced-order Kalman filters for sensor blending and noise rejection. Adaptation to flight condition is achieved with a novel gain-scheduling method based on correlation and regression analysis. The linear-optimal design approach is found to be a valuable tool in the development of practical multivariable control laws for vehicles which evidence significant coupling and insufficient natural stability.

Stengel, R. F.↗
NASA NTRS1979

The violent interstellar medium

Observational evidence for high-velocity and high-temperature interstellar gas is reviewed. The physical processes that characterize this gas are described, including the ionization and emissivity of coronal gas, the behavior and appearance of high-velocity shocks, and interfaces between coronal gas and cooler interstellar gas. Hydrodynamical models for the action of supernova explosions and stellar winds on the interstellar medium are examined, and recent attempts to synthesize all the processes considered into a global model for the interstellar medium are discussed.

Mccray, R.↗
NASA NTRS1985

Numerical solution of the vertical structure equation in the normal mode method

In the present model of multilayered stability stratification, aimed at obtaining the analytic eigensolutions of the vertical structure equation, each layer is characterized by its own static stability value. By requiring continuity of pressure and vertical velocity across each interface level, and by imposing suitable upper and lower boundary conditions, matching eigensolutions are obtained in terms of the Bessel functions. Attention is given to an explicit example of a double-layered stratified atmosphere which demonstrates the mathematical manipulations involved; the resultant vertical structure functions are used to check the accuracy of the numerical solutions by the finite difference and finite element methods.

Sasaki, Y. K.↗
NASA NTRS1992

Modal element method for scattering of sound by absorbing bodies

The modal element method for acoustic scattering from 2-D body is presented. The body may be acoustically soft (absorbing) or hard (reflecting). The infinite computational region is divided into two subdomains - the bounded finite element domain, which is characterized by complicated geometry and/or variable material properties, and the surrounding unbounded homogeneous domain. The acoustic pressure field is represented approximately in the finite element domain by a finite element solution, and is represented analytically by an eigenfunction expansion in the homogeneous domain. The two solutions are coupled by the continuity of pressure and velocity across the interface between the two subdomains. Also, for hard bodies, a compact modal ring grid system is introduced for which computing requirements are drastically reduced. Analysis for 2-D scattering from solid and coated (acoustically treated) bodies is presented, and several simple numerical examples are discussed. In addition, criteria are presented for determining the number of modes to accurately resolve the scattered pressure field from a solid cylinder as a function of the frequency of the incoming wave and the radius of the cylinder.

Baumeister, Kenneth J.↗
NASA NTRS1993

Thermal effects in rapid directional solidification - Linear theory

We study the morphological instability of the planar solid/liquid interface for a unidirectionally-solidified dilute binary mixture. We use a model developed by Boettinger et al. (1985, 1986), Aziz (1982), and Jackson et al. (1980), which allows for nonequilibrium effects on the interface through velocity-dependent segregation and attachment kinetics. Two types of instabilities are found in the linear stability analysis: (1) a cellular instability, and (2) an oscillatory instability driven by disequilibrium effects. Merchant and Davis (1990) characterized these instabilities subject to the frozen-temperature approximation (FTA). The present work relaxes the FTA by including the effects of latent heat and the full temperature distribution. Thermal effects slightly postpone the onset of the cellular instability but dramatically postpone the onset of the oscillatory instability; however, the absolute-stability conditions, at which at high speed the cellular and oscillatory instabilities are suppressed, remain unchanged from the FTA.

Huntley, D. A.↗
NASA NTRS1993

Modal element method for scattering and absorbing of sound by two-dimensional bodies

The modal element method for acoustic scattering from a 2-D body is presented. The body may be acoustically soft (absorbing) or hard (reflecting). The infinite computational region is divided into two subdomains - the bounded finite element domain, which is characterized by complicated geometry and/or variable material properties, and the surrounding unbounded homogeneous domain. The acoustic pressure field is represented approximately in the finite element domain by a finite element solution, and is represented analytically by an eigenfunction expansion in the homogeneous domain. The two solutions are coupled by the continuity of pressure and velocity across the interface between the two subdomains. Also, for hard bodies, a compact modal ring grid system is introduced for which computing requirements are drastically reduced. Analysis for 2-D scattering from solid and coated (acoustically treated) bodies is presented, and several simple numerical examples are discussed. In addition, criteria are presented for determining the number of modes to accurately resolve the scattered pressure field from a solid cylinder as a function of the frequency of the incoming wave and the radius of the cylinder.

Baumeister, K. J.↗
NASA NTRS1998

3D Simulations of the Richtmyer-Meshkov Instability with Re-Shock

We present results of inviscid simulations, in three dimensions, of Richtmyer-Meshkov instability for high incident shock Mach number. The growth rate of a single harmonic perturbation is quantified and compared with the results of a 2D calculation. Upon re-shock, the perturbation amplitude undergoes a phase reversal while the mean velocity of the interface is zero. Before re-shock the normalized growth rate of a 2D and 3D interface are nearly the same, but the growth rate after re-shock is significantly larger for the 3D than the 2D case. We also examine the evolution of multiple harmonic perturbations. Computational and parallelization issues of the simulation code will also be briefly discussed. The computations were done on the T3E at Pittsburgh Supercomputing Center.

Meiron, Daniel I.↗
NASA NTRS2020

Composite Piezoelectric Resonator 1D Modeling with Loss

Abstract – Maximizing the benefit of applying piezoelectric materials in composite structures, which are a mix of piezoelectric and elastic elements, requires analytical modeling. In this paper we describe two main classes of resonators and how to model using network equivalent including material losses. These include extensional devices where the layer areas are parallel to the direction of wave travel as found in a variety of thickness or length mode transducers (BC. Stress T1 = T2 at boundary) and transverse mode resonators (BC. Strain S1 = S2 at boundary) where the major layer area is perpendicular to the wave direction as is found in radial mode or length thickness mode composite resonators. In order to illustrate these models we will use Mason’s equivalent circuits with elastic, dielectric and piezoelectric loss. It is noted that in order to maintain consistency with the linear equations of piezoelectricity and the wave equation care is required when applying complex loss coefficients to the models. Although these techniques are applicable to other models we use Mason’s network equivalent circuit because it is intuitive for composite modeling due to the frequency independent turns ratio and the independence of the equations describing the elastic layers from the adjoining layer properties which aids in understanding the boundary conditions to impose on the layer boundaries. We will present network models for a variety of applications with loss and show how to calculate their impedance spectra, interface stress, strain velocity and displacement spectra.

Sherrit, Stewart↗
NASA NTRS2006

A Magnetohydrodynamic Boost for Relativistic Jets

We have performed relativistic magnetohydrodynamic simulations of the hydrodynamic boosting mechanism for relativistic jets explored by Aloy & Rezzolla (2006) using the RAISHIN code. Simulation results show that the presence of a magnetic field may change the properties of the shock interface between the tenuous, overpressured jet (V(sub j) (sup z)) flowing tangentially to a dense external medium. Magnetic fields can lead to more efficient acceleration of the jet, in comparison to the pure-hydrodynamic case. A poloidal magnetic field (B(sup z)), tangent to the interface and parallel to the jet flow, produces both a stronger outward moving shock and inward moving rarefaction wave. This leads to a large velocity component normal to the interface in addition to acceleration tangent to the interface, and the jet is thus accelerated to a larger Lorentz factors than those obtained in the pure-hydrodynamic case. In contrast, a strong toroidal magnetic field (B(sup y)), tangent to the interface but perpendicular to the jet flow, also leads to stronger acceleration tangent to the shock interface relative to the pure-hydrodynamic case, but to a lesser extent than found for the poloidal case due to the fact that the velocity component normal to the shock interface is now much smaller. Overall, the acceleration efficiency in the toroidal case is less than that of the poloidal case but both geometries still result in higher Lorentz factors than the pure-hydrodynamic case. Thus, the presence and relative orientation of a magnetic field in relativistic jets can have a significant influence on the hydrodynamic boost mechanism studied by Aloy & Rezzolla (2006).

Mizuno, Yosuke↗
NASA NTRS1983

Developments toward an 18% efficient silicon solar cell

Limitations to increased open-circuit voltage were identified and experimentally verified for 0.1 ohm-cm solar cells with heavily doped emitters. After major reduction in the dark current contribution from the metal-silicon interface of the grid contacts, the surface recombination velocity of the oxide-silicon interface of shallow junction solar cells is the limiting factor. In deep junction solar cells, where the junction field does not aid surface collection, the emitter bulk is the limiting factor. Singly-diffused, shallow junction cells have been fabricated with open circuit voltages in excess of 645 mV. Double-diffusion shallow and deep junctions cells have displayed voltages above 650 mV. MIS solar cells formed on 0.1 ohm-cm substrates have exibited the lowest dark currents produced in the course of the contract work.

Meulenberg, A., Jr.↗
DOE OSTI2021

Multiphase phase-field approach for solid–solid phase transformations via propagating interfacial phase in HMX

In this study, a thermodynamically consistent multiphase phase-field theory has been formulated to describe temperature-induced solid–solid phase transformations via an interfacial phase. This theory satisfies all thermodynamic equilibrium and stability conditions. It allows us to formulate the analytical solutions of the interface profile, energy, width, and velocity for each of the propagating solid–melt and solid–solid interfaces. The formation of the intermediate melt during solid–solid phase transformation has been studied hundreds of degrees below melting temperature in the HMX energetic crystal. The effect of the penalizing term in the phase-field model and two nanoscale material parameters on the appearance and degree of disordering of intermediate melt has been analyzed for the non-equilibrium interface. The developed approach is applicable to phase transitions in geological, pharmaceutical, ferroelectric, colloidal, and superhard materials where phase transformations occur via the intermediate phase.

Roy, Arunabha M. (ORCID:000000034790726X)↗
DOE OSTI2023

Shock state distributions in porous tantalum and characterization with multipoint velocimetry

Heterogenous materials under shock compression can be expected to reach different shock states throughout the material according to local differences in microstructure and the history of wave propagation. Here, in this work, a compact, multiple-beam focusing optic assembly is used with high-speed velocimetry to interrogate the shock response of porous tantalum films prepared through thermal-spray deposition. The distribution of particle velocities across a shocked interface is compared to results obtained using a set of defocused interferometric beams that sampled the shock response over larger areas. The two methods produced velocity distributions along the shock plateau with the same mean, while a larger variance was measured with narrower beams. The finding was replicated using three-dimensional, mesoscopically resolved hydrodynamics simulations of solid tantalum with a pore structure mimicking statistical attributes of the material and accounting for radial divergence of the beams, with agreement across several impact velocities. Accounting for pore morphology in the simulations was found to be necessary for replicating the rise time of the shock plateau. The validated simulations were then used to show that while the average velocity along the shock plateau could be determined accurately with only a few interferometric beams, accurately determining the width of the velocity distribution, which here was approximately Gaussian, required a beam dimension much smaller than the spatial correlation lengthscale of the velocity field, here by a factor of ∼30×, with implications for the study of other porous materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Impact of Carrier Concentration and Carrier Lifetime on MgZnO/CdSeTe/CdTe Solar Cells

Here, the performance impact of multiple parameters related to the use of MgZnO (MZO) as the emitter for bilayer CdSeTe/CdTe solar cells has been investigated in detail through numerical simulations. Such a comprehensive study is particularly important, because while cell fabrication using MZO has been highly successful in some cases, it has been less so in others, and it has not been clear which combinations of parameter values are most effective. The parameters considered here include the recombination velocity at the emitter/absorber interface, bulk recombination lifetime, and the carrier concentrations of the emitter and absorber. The ranges chosen for the simulation parameters are those most likely to be found experimentally. The primary finding is that independent of the interfacial recombination velocity and bulk recombination lifetime, the MZO carrier density should be > 1017 cm-3 and in any case greater than that of the absorber to reduce interface recombination. At the same time, a shallow dopant the order of 50 meV or less should reduce SRH recombination in the bulk and enable VOC > 1V.

14 SOLAR ENERGY↗
DOE OSTI2022

Joint Inversion of Body Wave Arrival Times, Surface Wave Dispersion Data and Receiver Functions: Method and Application to South China

In this work, we have developed a new joint inversion method that incorporates body wave arrival times, surface wave dispersion and receiver functions to simultaneously update earthquake locations and constrain three-dimensional P-wave (Vp) and S-wave velocity (Vs) models. Due to complementary sensitivities of the three types of data, the proposed joint inversion algorithm can reduce the intrinsic non-uniqueness of inversions using fewer types of data and better determine smooth velocity variations and velocity discontinuities. Synthetic tests demonstrate the advantages of this new joint inversion algorithm in resolving velocity structures, especially in constraining velocity gradients across the Moho interface. We have applied the proposed joint inversion algorithm to image the lithosphere velocity structure of south China. The inverted Vp and Vs models fit body wave arrival times, surface wave dispersion and receiver functions well. We further analyzed the distribution of Vs gradients across the Moho interface in detail, which helps us better understand tectonics in south China.

58 GEOSCIENCES↗
NASA NTRS1988

On the prediction of the behavior of the velocity, temperature and humidity fields above a thermally stratified air-water interface

A turbulent kinetic energy-dissipation (k-epsilon) model is employed to numerically simulate the behavior of the mean velocity, temperature and humidity fields, along with such second-order statistics as Reynolds stresses, turbulent heat and mass transfer correlation, and bulk transfer coefficients in the boundary layer above wind-generated laboratory waves. The k and epsilon equations have been properly modified to account for thermal stratification, the presence of free surface, and favorable pressure gradient effects. Predictions are compared with available data.

Papadimitrakis, Yiannis A.↗
DOE OSTI2026

Optical transparency of LiF (100) shock compressed to ∼360 GPa

High-purity [100] lithium fluoride (LiF) is the most widely used optical window in dynamic compression experiments due to its wide bandgap and well-characterized mechanical response. Recent plate-impact experiments established the [100] LiF Hugoniot to ∼230 GPa and demonstrated shock-induced melting onset at 182 GPa with complete melting by 195 GPa; theoretical models predict LiF optical transparency to nearly 900 GPa. To experimentally examine the optical transparency of [100] LiF at higher pressures and in the liquid state, laser-driven shock experiments were performed at peak stresses ranging from 223 to 363 GPa. Optical response was examined by measuring the particle velocity histories at the Kapton/LiF interface using laser interferometry at 532 and 1550 nm wavelengths; in-material particle velocities were obtained using established refractive-index corrections. Continuous photonic Doppler velocimetry fringes were observed across the entire stress range, demonstrating that LiF remains transparent to 1550 nm light throughout the multi-megabar regime investigated. At 532 nm, fringe visibility depended on the reflector coating: aluminum mirrors provided signals to ∼235 GPa, while gold mirrors extended this limit to 270.5 GPa, indicating that the shorter-wavelength response is likely sensitive to experimental configuration rather than to the loss of LiF transparency. Continued optical transparency to at least 360 GPa indicates that shock-melted LiF does not display bandgap closure over the stress range explored. Furthermore, these results provide direct experimental constraints on the high-pressure optical response and establish LiF (100) as a robust optical window material for laser-driven dynamic compression experiments approaching 400 GPa.

Renganathan, P. [Argonne National Laboratory (ANL)↗