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 145 records · Page 8

Inhomogeneous diffusion-limited aggregation

It is demonstrated here that inhomogeneous diffusion-limited aggregation (DLA) model can be used to simulate viscous fingering in a medium with inhomogeneous permeability and homogeneous porosity. The medium consists of a pipe-pore square-lattice network in which all pores have equal volume and the pipes have negligible volume. It is shown that fluctuations in a DLA-based growth process may be tuned by noise reduction, and that fluctuations in the velocity of the moving interface are multiplicative in form.

Selinger, Robin Blumberg↗

Inverse heat transfer analysis of Bridgman crystal growth

The effects of asymmetry in furnace temperature profile and pulling velocity on the crystal interface shape are demonstrated while neglecting the latent heat of solidification. It is seen that the furnace temperature profile may be varied in order to influence the shape of the melt-crystal interface. An exact thermal analysis is then performed on the Bridgman technique by including the latent heat of solidification as a source term. The exact temperature field required for yielding a flat melt-crystal interface is obtained. The earlier observation regarding the influence of furnace temperature profile on the interface shape is confirmed and a criterion for achieving a flat interface is obtained. Various furnace temperature profiles are selected and their corresponding melt-crystal interface results are presented.

Taghavi, K.↗

Thermocapillary convection in floating zones under simulated reduced-gravity conditions

The present study demonstrated that calculated thermocapillary convection in a non-cylindrical floating zone can now be compared with measured one, by considering the lens effect of the floating zone. Flow visualization and computer simulation of thermocapillary convection in a silicone oil zone and a molten zone in an NaNO3 rod were conducted. The calculated results agree very well with the measured ones, including the free surface shapes, the solid/melt interface shapes and the velocity fields.

Tao, Y.↗

Dendritic Growth Velocity of Ice Crystals

The dendritic growth velocity of ice crystals is measured at undercooling levels between 5-8.5 K using a unique apparatus. The result shows that the measured velocities are smaller that the theoretical prediction based on the thermal diffusion theory with the scalling constant, Q= 0.025, which has been shown to agree with the experimental results at lower undercooling levels between 0.4 -2.5 K.

Dendritic Ice Growth velocity Undercooling Interfa↗

Convection and morphological stability during directional solidification

For growth of a vicinal face at constant velocity, the effect of anisotropic interface kinetics on morphological stability is calculated for a binary alloy. The dependence of the interface kinetic coefficient on crystallographic orientation is based on the motion and density of steps. Anisotropic kinetics give rise to traveling waves along the crystal-melt interface, and can lead to a significant enhancement of morphological stability. The stability enhancement increases as the orientation approaches a singular orientation and as the solidification velocity increases. Shear flows interact with the traveling waves and, depending on the direction of the flow, may either stabilize or destabilize the interface. Specific calculations are carried out for germanium-silicon alloys.

Coriell, Sam R.↗

SOFIA Observations of S106: Dynamics of the Warm Gas

Context The H II region/PDR/molecular cloud complex S106 is excited by a single O-star. The full extent of the warm and dense gas close to the star has not been mapped in spectrally resolved high-J CO or [C II] lines, so the kinematics of the warm. partially ionized gas, are unknown. Whether the prominent dark lane bisecting the hourglass-shaped nebula is due solely to the shadow cast by a small disk around the exciting star or also to extinction in high column foreground gas was an open question until now. Aims. To disentangle the morphology and kinematics of warm neutral and ionized gas close to the star, study their relation to the bulk of the molecular gas. and to investigate the nature of the dark lane. Methods. We use the heterodyne receiver GREAT on board SOFIA to observe velocity resolved spectral lines of [C II] and CO 11 yields 10 in comparison with so far unpublished submm continuum data at 350 micron (8HARC-Il) and complementary molecular line data. Results. The high angular and spectral resolution observations show a very complex morphology and kinematics of the inner S106 region, with many different components at different excitation conditions contributing to the observed emission. The [C II] lines are found to be bright and very broad. tracing high velocity gas close to the interface of molecular cloud and H II region. CO 11 yields 10 emission is more confined.. both spatially and in velocity, to the immediate surroundings of S 106 IR showing the presence of warm, high density (clumpy) gas. Our high angular resolution submm continuum observations rule out the scenario where the dark lane separating the two lobes is due solely to the shadow cast by a small disk close to the star. The lane is clearly seen also as warm, high column density gas at the boundary of the molecular cloud and H II region.

Simon, R.↗

Signal processing and display interface studies

Signal processing equipment specifications, operating and test procedures, and systems design and engineering are described. Five subdivisions of the overall circuitry are treated: (1) the spectrum analyzer; (2) the spectrum integrator; (3) the velocity discriminator; (4) the display interface; and (5) the formatter. They function in series: (1) first in analog form to provide frequency resolution, (2) then in digital form to achieve signal to noise improvement (video integration) and frequency discrimination, and (3) finally in analog form again for the purpose of real-time display of the significant velocity data. The formatter collects binary data from various points in the processor and provides a serial output for bi-phase recording. Block diagrams are used to illustrate the system.

Source record↗

Cross-Cutting Metrology Tools for In Operando Characterization of Carrier Dynamics in Photovoltaic Devices (Final Technical Report)

Understanding the nature of recombination and its dependence on defects and interfaces is essential for engineering materials and contacts for higher V oc and photovoltaic (PV) efficiency. Time-resolved photoluminescence (TRPL) has conventionally been used to evaluate recombination, but not all materials are strongly emissive or otherwise suitable. Time-resolved terahertz spectroscopy is a pump – probe method that presents valuable complementary information, wherein an optical pump pulse photoexcites carriers within the absorber and the transient photoconductivity is probed with a terahertz pulse. Until now, experimental constraints have prohibited the use of terahertz probes to interrogate full device stacks, and measurements were instead made on exfoliated films or films grown on unconventional substrates. However, interfaces are critical to the behavior of photoexcited carriers in solar cells, and the substrates themselves often influence the film growth and bulk properties. Therefore, it is important to probe the behavior of PV absorbers as close as possible to their normal operating conditions. Here we have developed cross-cutting metrology that enable in operando characterization of carrier dynamics and recombination mechanisms in working PV devices. Through a combination of complementary non-contact TRTS and TRPL experiments and modeling, we have obtained key parameters including photoexcited carrier bulk lifetime and interface and back surface recombination velocities with greater precision and accuracy than can be achieved with conventional TRPL alone. By varying pulsed photoexcitation conditions, the contributions of interface and bulk recombination mechanisms can be determined. We have developed and refined this characterization approach by investigating CdTe PVs, which is a well-established technology but with significant margin for further improvements in efficiency. Results were published in Journal of Applied Physics (DOI: 10.1063/5.0064730) and Proceedings of the IEEE PVSC (DOI: 10.1109/PVSC43889.2021.9518559). After validating the approach, we have applied it to correlate recombination rates to processing conditions and material properties, including composition and defects as well as interfaces and surface treatments. Such feedback can inform processing and design choices to enable higher PV efficiency, not only for CdTe, but also CdSe x Te 1-x , CIGS, perovskites, kesterites, and future technologies. A perspective on predicting solar cell performance from terahertz spectroscopy was published in Advanced Energy Materials (DOI: 10.1002/aenm.202102776). The cross-cutting metrology tools developed here will enable determination of the locus and mechanism of performance-limiting recombination in thin film PV devices. This feedback can guide the engineering of PV devices with higher V oc and efficiency, thereby leading to reductions in levelized cost of electricity to meet the SunShot 2030 target of $0.03/kWh for utility-scale PV. Lower costs will enable rapid expansion of renewable electricity generation that is nearly free of carbon emissions. Reducing carbon emissions of the electricity sector is one of the most important solutions to climate change. The project’s specific focus on CdTe PV can benefit First Solar, the US company that is the global leader in CdTe technology, through both technology and workforce development.

14 SOLAR ENERGY↗

Computationally Aided Design of a Small-Scale Waste-Glass Melter - 20452

To provide mission support to the WTP for the vitrification of legacy nuclear tank waste at the Hanford site, a reduced-scale vitrification pilot system is being designed to process simulated and actual radioactive tank wastes. The waste stored in underground tanks will be separated into HLW and LAW fractions, where the majority is LAW by mass (approximately 90%) and HLW by activity (greater than 95%). The first tank waste to be processed at the WTP will be LAW. Because Hanford tank waste and the resultant melter feed compositions are known to vary widely, pilot-scale operations are essential to identify potential problems, provide needed data, confirm assumptions, and determine impacts to full-scale melters and off-gas systems from these different feeds. In addition, the reduced-scale melter system must be capable of quickly providing results to operations without the typically high cost of full-scale radioactive operations and to provide an engineering platform in close proximity to local operations staff. The objectives of this computational modeling effort are to produce process conditions like those encountered in the WTP LAW melters and to minimize the volume of radioactive waste necessary for the evaluations. To this end, melter design features are being explored to increase production without significantly increasing surface area or melter-glass volumes or compromising data quality. In support of these objectives, a set of computational fluid dynamic models was developed to provide insight into the flow patterns within these small melters and evaluate design features to increase production without significantly increasing glass inventory. Computed velocities were evaluated at the interface between the cold cap and the molten-glass pool and within the bulk-glass pool. Since melt rate correlates with the velocity underneath the cold cap, this information can be used to select design features that produce the highest gains. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Magnetohydrodynamic Boost for Relativistic Jets

We 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 changes the properties of the shock interface between the tenuous, overpressured jet (V^z j) flowing tangentially to a dense external medium. We find that magnetic fields can lead to more efficient acceleration of the jet, in comparison to the pure-hydrodynamic case. A "poloidal" magnetic field (B^z), tangent to the interface and parallel to the jet flow, produces both a stronger outward moving shock and a stronger 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 larger Lorentz factors than those obtained in the pure-hydrodynamic case. Likewise, a strong "toroidal" magnetic field (B^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. Thus. the presence and relative orientation of a magnetic field in relativistic jets can significant modify the hydrodynamic boost mechanism studied by Aloy & Rezzolla (2006).

Mizuno, Yosuke↗

Experimental Studies of the Interaction Between a Parallel Shear Flow and a Directionally-Solidifying Front

It has long been recognized that flow in the melt can have a profound influence on the dynamics of a solidifying interface and hence the quality of the solid material. In particular, flow affects the heat and mass transfer, and causes spatial and temporal variations in the flow and melt composition. This results in a crystal with nonuniform physical properties. Flow can be generated by buoyancy, expansion or contraction upon phase change, and thermo-soluto capillary effects. In general, these flows can not be avoided and can have an adverse effect on the stability of the crystal structures. This motivates crystal growth experiments in a microgravity environment, where buoyancy-driven convection is significantly suppressed. However, transient accelerations (g-jitter) caused by the acceleration of the spacecraft can affect the melt, while convection generated from the effects other than buoyancy remain important. Rather than bemoan the presence of convection as a source of interfacial instability, Hurle in the 1960s suggested that flow in the melt, either forced or natural convection, might be used to stabilize the interface. Delves considered the imposition of both a parabolic velocity profile and a Blasius boundary layer flow over the interface. He concluded that fast stirring could stabilize the interface to perturbations whose wave vector is in the direction of the fluid velocity. Forth and Wheeler considered the effect of the asymptotic suction boundary layer profile. They showed that the effect of the shear flow was to generate travelling waves parallel to the flow with a speed proportional to the Reynolds number. There have been few quantitative, experimental works reporting on the coupling effect of fluid flow and morphological instabilities. Huang studied plane Couette flow over cells and dendrites. It was found that this flow could greatly enhance the planar stability and even induce the cell-planar transition. A rotating impeller was buried inside the sample cell, driven by an outside rotating magnet, in order to generate the flow. However, it appears that this was not a well-controlled flow and may also have been unsteady. In the present experimental study, we want to study how a forced parallel shear flow in a Hele-Shaw cell interacts with the directionally solidifying crystal interface. The comparison of experimental data show that the parallel shear flow in a Hele-Shaw cell has a strong stabilizing effect on the planar interface by damping the existing initial perturbations. The flow also shows a stabilizing effect on the cellular interface by slightly reducing the exponential growth rate of cells. The left-right symmetry of cells is broken by the flow with cells tilting toward the incoming flow direction. The tilting angle increases with the velocity ratio. The experimental results are explained through the parallel flow effect on lateral solute transport. The phenomenon of cells tilting against the flow is consistent with the numerical result of Dantzig and Chao.

Zhang, Meng↗

The growth and nanothermite reaction of 2Al/3NiO multilayer thin films

Nanothermite NiO–Al is a promising material system for low gas emission heat sources; still, its reactive properties are highly dependent on material processing conditions. In the current study, sputter deposition is used to fabricate highly controlled nanolaminates comprised of alternating NiO and Al layers. Films having an overall stoichiometry of 2Al to 3NiO were produced with different bilayer thicknesses to investigate how ignition and self-sustained, high temperature reactions vary with changes to nanometer-scale periodicity and preheat conditions. Ignition studies were carried out with both hot plate and laser irradiation and compared to slow heating studies in hot-stage x-ray diffraction. Ignition behavior has bilayer thickness and heating rate dependencies. The 2Al/3NiO with λ ≤ 300 nm ignited via solid/solid diffusion mixing (activation energy, E a = 49 ± 3 kJ/mole). Multilayers having λ ≥ 500 nm required a more favorable mixing kinetics of solid/liquid dissolution into molten Al (E a = 30 ± 4 kJ/mole). This solid/liquid dissolution E a is a factor of 5 lower than that of the previously reported powder compacts due to the elimination of a passivating Al oxide layer present on the powder. The reactant mixing mechanism between 300 and 500 nm bilayer thicknesses was dependent on the ignition source's heating rate. The self-propagating reaction velocities of 2Al/3NiO multilayers varied from 0.4 to 2.5 m/s. Pre-heating nanolaminates to temperatures below the onset reaction temperatures associated with forming intermediate nickel aluminides at multilayer interfaces led to increased propagation velocities, whereas pre-heating samples above the onset temperatures inhibited subsequent attempts at laser ignition.

36 MATERIALS SCIENCE↗

A sharp interface Lagrangian-Eulerian method for rigid-body fluid-structure interaction

Herein this paper introduces a sharp interface method to simulate fluid-structure interaction (FSI) involving rigid bodies immersed in viscous incompressible fluids. The capabilities of this methodology are benchmarked using a range of test cases and demonstrated using large-scale models of biomedical FSI. The numerical approach developed herein, which we refer to as an immersed Lagrangian-Eulerian (ILE) method, integrates aspects of partitioned and immersed FSI formulations by solving separate momentum equations for the fluid and solid subdomains, as in a partitioned formulation, while also using non-conforming discretizations of the dynamic fluid and structure regions, as in an immersed formulation. A simple Dirichlet-Neumann coupling scheme is used, in which the motion of the immersed solid is driven by fluid traction forces evaluated along the fluid-structure interface, and the motion of the fluid along that interface is constrained to match the solid velocity and thereby satisfy the no-slip condition. To develop a practical numerical method, we adopt a penalty approach that approximately imposes the no-slip condition along the fluid-structure interface. In the coupling strategy, a separate discretization of the fluid-structure interface is tethered to the volumetric solid mesh via stiff spring-like penalty forces. Our fluid-structure coupling scheme relies on an immersed interface method (IIM) for discrete geometries, which enables the accurate determination of both velocities and stresses along complex internal interfaces. Numerical methods for FSI can suffer from instabilities related to the added mass effect, but computational tests indicate that the methodology introduced here remains stable for selected test cases across a broad range of solid-fluid mass density ratios, including extremely small, nearly equal, equal, and large density ratios. Biomedical FSI demonstration cases include results obtained using this method to simulate the dynamics of a bileaflet mechanical heart valve in a pulse duplicator, and to model transport of blood clots in a patient-averaged anatomical model of the inferior vena cava.

97 MATHEMATICS AND COMPUTING↗

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

Large-scale atomistic studies of sliding friction in polycrystalline aluminum interfaces

In this work, we discuss large-scale non-equilibrium molecular dynamics (NEMD) simulations of ductile metal sliding comprising up to 1.8 × 10 9 atoms over time scales of 100 ns. The results of these simulations have identified a variety of physical mechanisms that are important in determining the steady-state frictional force for a wide range of velocities at compressed metal–metal interfaces. These include grain growth and refinement, the evolution of large plastic strains and strain rates, material mixing, and melting. These phenomena can be included in a strain, strain rate, and grain size model that gives good agreement with the NEMD simulations and can be applied to macroscopic continua.

74 ATOMIC AND MOLECULAR PHYSICS↗

Experimental study of the minority-carrier transport at the polysilicon-monosilicon interface

This paper presents the results of an experimental study designed to explore both qualitatively and quantitatively the mechanism of the improved current gain in bipolar transistors with polysilicon emitter contacts. Polysilicon contacts were deposited and heat treated at different conditions. The electrical properties were measured using p-n junction test structures that are much more sensitive to the contact properties than are bipolar transistors. A simple phenomenological model was used to correlate the structural properties with electrical measurements. Possible transport mechanisms are examined and estimates are made about upper bounds on transport parameters in the principal regions of the devices. The main conclusion of this study is that the minority-carrier transport in the polycrystalline silicon is dominated by a highly disordered layer at the polysilicon-monosilicon interface characterized by very low minority-carrier mobility. The effective recombination velocity at the n(+) polysilicon-n(+) monosilicon interface was found to be a strong function of fabrication conditions. The results indicate that the recombination velocity can be much smaller than 10,000 cm/s.

Neugroschel, A.↗

Understanding the Effects of Inhomogeneities at the Back Interface of CdTe‐Based Solar Cells Using 2D Modeling

One-dimensional modeling cannot capture lateral inhomogeneities in CdTe-based devices. Here, we use 2D modeling to investigate the role of varying energetics at the back interface. We consider improvements in the back interface layer (BIL) through either reducing back surface recombination velocity (BSRV) or decreasing the downward band bending near the back interface. We show that when the BSRV is reduced, but strong downward band bending remains, there is no change in the device performance until the BSRV of 90% of the back interface is improved by the BIL. On the other hand, any coverage with a BIL that improves band bending results in device improvements. We use band bending, back interface recombination current densities, and voltage dependent current flow through the device to understand these improvements. The modeling shows that lateral flow of carriers greatly affects device performance, which is not captured in parallel diode modeling, and demonstrates improved understanding with 2D modeling.

2D modelling↗

Inflectional instabilities in the wall region of bounded turbulent shear flows

The primary thrust of this research was to identify one or more mechanisms responsible for strong turbulence production events in the wall region of bounded turbulent shear flows. Based upon previous work in a transitional boundary layer, it seemed highly probable that the production events were preceded by an inflectional velocity profile which formed on the interface between the low-speed streak and the surrounding fluid. In bounded transitional flows, this unstable profile developed velocity fluctuations in the streamwise direction and in the direction perpendicular to the sheared surface. The rapid growth of these instabilities leads to a breakdown and production of turbulence. Since bounded turbulent flows have many of the same characteristics, they may also experience a similar type of breakdown and turbulence production mechanism.

Swearingen, Jerry D.↗