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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Reduction of Convection in Closed Tube Vapor Transport Experiments

The primary objective of this effort was to develop a method for suppressing convective flows during the growth of mercurous chloride crystals by vapor transport in closed tubes to levels approaching those obtained in the microgravity environment. Mercurous chloride was chosen because it is a technologically interesting acoustical optical material whose optical properties are believed to be affected by convective flows. Since the Grashof number scales as the cube of the smallest dimension in the flow system, reduction of the size scale can be extremely effective in reducing unwanted convective flows. However, since materials of practical interest must be grown at least on the cm scale, reduction of the overall growth system is not feasible. But if the region just above the growing crystal could be restricted to a few mm, considerable reduction in flow velocity would result. By suspending an effusive barrier in the growth ampoule just above the growth interface, it should be possible to reduce the convective velocity in this vicinity to levels approaching flows in microgravity. If successful, this growth technique will offer a screening test for proposed space experiments that involve vapor transport to see if reduction of convection will result in improved material and will set a new standard against which the improvements obtained in microgravity may be judged. In addition, it may provide an improved method for preparing materials on Earth whose growth is affected adversely by convection. If the properties of this material can be improved there is a potential commercial interest from Brimrose Inc., who has agreed to fabricate and test devices from the crystals we have grown. This report describes the development of the growth facility, the purification processes developed for preparing the starting material, and the results from growth experiments with and without the effusive baffle. Mercurous chloride turned out to be a more difficult material to deal with than originally anticipated. At growth temperatures, it is extremely sensitive to practically any impurity which causes it to form oxychlorides and/or to decompose into elemental mercury and bichloride of mercury. We were unable to find a suitable method for protecting the magnetic material used to suspend the effusion barrier from the attack of mercurous chloride vapor. Although we were successful in growing single crystals of mercurous chloride without the effusion baffle, they exhibited severe microcracking which we attribute to wall-induced thermal stresses. This leads us to believe that uncontrolled convection may not be the most important problem in the development of this material and a new growth process was attempted that eliminates the wall-induced stress. Unfortunately, the grant ran out before this new method could be adequately tested.

Naumann, R. J.↗
NASA NTRS2003

Macrovoid Defect Growth during Evaporative Casting of Polymeric Membranes

Macrovoid (MV) formation is a significant problem in evaporatively cast polymeric membranes. MVs are large, elongated or teardrop-shaped pores (~10-50 micron) that can impair membrane structural integrity. Although MVs have been extensively studied, there is no general agreement on the mechanisms governing MV growth. Recently, our research group has formulated the solutocapillary convection (SC) hypothesis, which contends that MV growth involves three principal forces: a Marangoni force generated by surface tension gradients within the MV interface, a viscous drag force, and a gravitationally induced body force. Two sets of complementary experiments were conducted to test the SC hypothesis. Ground-based videomicroscopy flow-visualization (VMFV) was utilized to measure the flow velocities at the MV-casting solution interface and deep within the casting solution. The measurements were performed with casting solutions containing 10 wt% cellulose acetate (CA), 30 wt% H2O, 60 wt% acetone, and 200- ppm TiO2 particles for flow visualization, and the surface tension was controlled by surfactant addition. Qualitatively, the experiments indicated that MV growth occurs in three distinct phases: (1) a very rapid initial growth period, (2) a much slower growth phase, and (3) absorption of selected MVs into the expanding demixed region. The presence of tracer particles inside the MVs suggests the presence of a convective flow, which transfers the particles from the bulk solution to the MV interior. Although the VMFV experiments did not establish any surfactant effect on the interfacial velocities, a statistically significant effect on the MV number density was observed. In the second set of experiments, membranes were cast aboard a KC-135 aircraft under 0-g and 2-g conditions. Despite careful attention to the design and fabrication of the membrane casting apparatus (MCA), several problems were encountered, the most significant of which was the contamination of the casting solution by the activated carbon particles used for solvent absorption.

Greenberg, A. R.↗
NASA NTRS1990

Comparisons of one- and two- interface methods for tunnel wall interference calculation

Tunnel wall interference in two-dimensional subsonic flow can be analytically calculated on the basis of the one- or two-interface method. The application of numerical integration by means of matrix multiplication to both methods has been found to significantly reduce CPU time in on-line operation. An investigation of the sensitivity of these methods to random velocity-measurement disturbances shows the one-interface method to be less sensitive to measurement errors, and therefore to be preferable for practical applications.

Lo, C. F.↗
DOE OSTI2024

A high-order Shifted Interface Method for Lagrangian shock hydrodynamics

Here, we present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that is, it uses the precise location of the material interface. The interface is represented by the zero level-set of a continuous high-order finite element function that moves with the material velocity. This strategy allows to evolve curved material interfaces inside curved elements. By reformulating the original interface problem over a surrogate (approximate) interface, located in proximity of the true interface, the SIM avoids cut cells and the associated problematic issues regarding implementation, numerical stability, and matrix conditioning. Accuracy is maintained by modifying the original interface conditions using Taylor expansions. We demonstrate the performance of the proposed algorithms on established numerical benchmarks in one, two and three dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Multi-scale characterization and simulation of impact welding between immiscible Mg/steel alloys

Vaporizing foil actuator spot welding method is used in this paper to join magnesium alloy AZ31 and uncoated high-strength steel DP590, which are typically considered as un-weldable due to their high physical property disparities, low mutual solubility, and the lack of any intermetallic phases. Characterization results from scanning electron microscopy (SEM) and high-resolution transmission electron microscopy (HRTEM) of the weld interface indicate that the impact creates an Mg nanocrystalline interlayer with abundant Fe particles. The interlayer exhibits intact bonding with both DP590 and AZ31 substrates. To investigate the fundamental bond formation mechanisms at the interface, a finite element (FE)-based process simulation is first performed to calculate the local temperature and deformation at the interface under the given macroscopic experimental condition. Finally, taking the FE results at the interface as inputs, molecular dynamics (MD) simulations are conducted to study the interlayer formation at the Mg/Fe interface during the impact and cooling. The results found a high velocity shearing-induced mechanical mixing mechanism that mixes Mg/Fe atoms at the interface and creates the interlayer, leading to the metallurgical bond between Mg/steel alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A theoretical/computational framework to measure SiO2 and MgO viscosity at high pressure

The convection of the mantle of Earth and super-Earths is important for many terrestrial phenomena, from plate tectonics to outgassing. Rheological properties, such as viscosity, regulate the transport of thermal energy and mass. However, the viscosity of mantle-relevant materials, such as MgO, at relevant pressures (>120 GPa) are not well constrained. The objective of this work was to develop a computational platform to simulate novel experiments aiming to measure the viscosity of MgO at high pressures. Experiments performed by our collaborators at Johns Hopkins University and Lawrence Livermore National Laboratory use the OMEGA EP laser facility to shock a corrugated MgO interface to 170 GPa, with resulting velocity evolution governed by the viscous Richtmyer-Meshkov instability. We used an in-house hydrocode to simulate this process and thus provide a bound on viscosity by comparing our simulations results to these experiments, as well as to examine the physical processes at play. We simulated the experiments with different values of MgO viscosity (from inviscid to 10,000 Pa∙s) while taking into account the unsteady laser pulse, the material's equation of state, and the rate-dependent constitutive relation for MgO, and the materials' equations of state. Our results suggest that MgO at these conditions has a viscosity within an order of magnitude of 5000 Pa∙s.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Electrolyte-Assisted Hydrogen Cycling in Lithium and Sodium Alanates at Low Pressures and Temperatures

An investigation of electrolyte-assisted hydrogen storage reactions in complex aluminum hydrides (LiAlH4 and NaAlH4) reveals significantly reduced reaction times for hydrogen desorption and uptake in the presence of an electrolyte. LiAlH4 evolves ~7.8 wt% H2 over ~3 h in the presence of a Li-KBH4 eutectic at 130 °C compared to ~25 h for the same material without the electrolyte. Similarly, NaAlH4 exhibits 4.8 wt% H2 evolution over ~4 h in the presence of a diglyme electrolyte at 150 °C compared to 4.4 wt% in ~15 h for the same material without the electrolyte. These reduced reaction times are composed of two effects, an increase in reaction rates and a change in the reaction kinetics. While typical solid state dehydrogenation reactions exhibit kinetics with rates that continuously decrease with the extent of reaction, we find that the addition of an electrolyte results in rates that are relatively constant over the full desorption window. Fitting the kinetics to an Avrami-Erofe’ev model supports these observations. The desorption rate coefficients increase in the presence of an electrolyte, suggesting an increase in the velocities of the reactant-product interfaces. In addition, including an electrolyte increases the growth parameters, primarily for the second desorption steps, resulting in the observed relatively constant reaction rates. Similar effects occur upon hydrogen uptake in NaH/Al where the presence of an electrolyte enables hydrogenation under more practical low temperature (75 °C) and pressure (50 bar H2) conditions.

25 ENERGY STORAGE↗
NASA NTRS1985

Modeling of convection phenomena in Bridgman-Stockbarger crystal growth

Thermal convection phenomena in a vertically oriented Bridgman-Stockbarger apparatus were modeled by computer simulations for different gravity conditions, ranging from earth conditions to extremely low gravity, approximate space conditions. The modeling results were obtained by the application of a state-of-the art, transient, multi-dimensional, completely densimetrically coupled, discrete-element computational model which was specifically developed for the simulation of flow, temperature, and species concentration conditions in two-phase (solid-liquid) systems. The computational model was applied to the simulation of the flow and the thermal conditions associated with the convection phenomena in a modified Germanium-Silicon charge enclosed in a stationary fused-silica ampoule. The results clearly indicated that the gravitational field strength influences the characteristics of the coherent vortical flow patterns, interface shape and position, maximum melt velocity, and interfacial normal temperature gradient.

Carlson, F. M.↗
NASA NTRS1990

Analysis of the thermal performance of heat pipe radiators

A comprehensive mathematical model and computational methodology are presented to obtain numerical solutions for the transient behavior of a heat pipe radiator in a space environment. The modeling is focused on a typical radiator panel having a long heat pipe at the center and two extended surfaces attached to opposing sides of the heat pipe shell in the condenser section. In the set of governing equations developed for the model, each region of the heat pipe - shell, liquid, and vapor - is thermally lumped to the extent possible, while the fin is lumped only in the direction normal to its surface. Convection is considered to be the only significant heat transfer mode in the vapor, and the evaporation and condensation velocity at the liquid-vapor interface is calculated from kinetic theory. A finite-difference numerical technique is used to predict the transient behavior of the entire radiator in response to changing loads.

Boo, J. H.↗
NASA NTRS1991

Adaptive temperature profile control of a multizone crystal growth furnace

An intelligent measurement system is described which is used to assess the shape of a crystal while it is growing inside a multizone transparent furnace. A color video imaging system observes the crystal in real time, and determines the position and the shape of the interface. This information is used to evaluate the crystal growth rate, and to analyze the effects of translational velocity and temperature profiles on the shape of the interface. Creation of this knowledge base is the first step to incorporate image processing into furnace control.

Batur, C.↗
NASA NTRS1992

Identification and control of a multizone crystal growth furnace

This paper presents an intelligent adaptive control system for the control of a solid-liquid interface of a crystal while it is growing via directional solidification inside a multizone transparent furnace. The task of the process controller is to establish a user-specified axial temperature profile and to maintain a desirable interface shape. Both single-input-single-output and multi-input-multi-output adaptive pole placement algorithms have been used to control the temperature. Also described is an intelligent measurement system to assess the shape of the crystal while it is growing. A color video imaging system observes the crystal in real time and determines the position and the shape of the interface. This information is used to evaluate the crystal growth rate, and to analyze the effects of translational velocity and temperature profiles on the shape of the interface. Creation of this knowledge base is the first step to incorporate image processing into furnace control.

Batur, C.↗
NASA NTRS1996

Modeling of Bulk Evaporation and Condensation

This report describes the modeling and mathematical formulation of the bulk evaporation and condensation involved in liquid-vapor phase change processes. An internal energy formulation, for these phase change processes that occur under the constraint of constant volume, was studied. Compared to the enthalpy formulation, the internal energy formulation has a more concise and compact form. The velocity and time scales of the interface movement were obtained through scaling analysis and verified by performing detailed numerical experiments. The convection effect induced by the density change was analyzed and found to be negligible compared to the conduction effect. Two iterative methods for updating the value of the vapor phase fraction, the energy based (E-based) and temperature based (T-based) methods, were investigated. Numerical experiments revealed that for the evaporation and condensation problems the E-based method is superior to the T-based method in terms of computational efficiency. The internal energy formulation and the E-based method were used to compute the bulk evaporation and condensation processes under different conditions. The evolution of the phase change processes was investigated. This work provided a basis for the modeling of thermal performance of multi-phase nuclear fuel elements under variable gravity conditions, in which the buoyancy convection due to gravity effects and internal heating are involved.

Anghaie, S.↗
NASA NTRS1999

In-Situ Monitoring of a Bismuth-Tin Alloy Solidification Using MEPHISTO

Experiments were carried out to study the morphological stability of Bi- 1 atomic % Sn alloys using the MEPHISTO directional solidification apparatus aboard Space Shuttle Columbia (STS-87, launched Nov. 19, 1997) and in ground-based studies. The Seebeck signal and temperature measurements indicate that convection was significant for ground-based studies. In the space-based experiments, interface breakdown was observed at growth velocities of 6.7, 27, and 40 microns/sec, but not at 1.8 and 3.3 microns/sec.

Abbaschian, Reza↗
DOE OSTI2022

CFD simulations of electric motor end ring cooling for improved thermal management

Proper thermal management of an electric motor for vehicle applications extends its operating range. One cooling approach is to impinge Automatic Transmission Fluid (ATF) onto the rotor end ring. Increased ATF coverage correlates to enhanced heat transfer. Computational Fluid Dynamics (CFD) analytical tools provide a mechanism to assess motor thermal management prior to hardware fabrication. The complexity of the fluid flow (e.g., jet atomization, interface tracking, wall impingement) and heat transfer makes these simulations challenging. Computational costs are high when solving these flows on high-speed rotating meshes. Typically, a Volume-of Fluid (VOF) technique (i.e., two-fluid system) is used to resolve ATF dynamics within this rotating framework. Suitable numerical resolution of the relevant physics for thin films under strong inertial forces at high rotor speeds is computationally expensive, further increasing the run times. In this work, a numerical study of rotor-ring cooling by ATF is presented using a patent automated Cartesian cut-cell based method coupled with Automatic Mesh Refinement (AMR). This approach automatically creates the Cartesian mesh on-the-fly and can effectively handle complex rotating geometries by adaptively refining the mesh based on local gradients in the flow field which results in better resolution of the air-ATF interface. A Single non-inertial Reference Frame (SRF) approach is used to account for the rotating geometry and to further improve the overall computational efficiency. Quasi-steady state conditions are targeted in the analysis of the results. Important physics such as ATF jet structure, velocity detail near the air-jet interface, ATF coverage/accumulation on the ring surface, and cooling capacity are presented for a low-resolution Reynolds averaged Navier-Stokes (RANS), high-resolution RANS, and high-resolution Large-Eddy Simulation (LES) models. Computations are scaled over hundreds of cores on a supercomputer to maximize turnaround time. Each numerical approach is shown to capture the general trajectory of the oil jet prior to surface impingement. The high-resolution LES simulation, however, is superior in capturing small scale details and heat transfer between the free jet and surrounding air.

42 ENGINEERING↗
DOE OSTI2025

Phase-field model of alloy solidification far from chemical equilibrium at the solid-liquid interface

We further develop a recently introduced phase-field model of far-from-equilibrium alloy solidification under additive manufacturing conditions [K. Ji et al., Phys. Rev. Lett. 130, 026203 (2023)]. This model utilizes enhanced solute diffusivity within the spatially diffuse interface region to quantitatively capture solute trapping with a larger interface width, thereby making simulations on experimentally relevant length and timescales computationally feasible. The main developments presented here include testing the robustness of different variational formulations, extending the model to concentrated alloys by incorporating solid and liquid free energies from thermodynamic databases, as illustrated for hypoeutectic Al-Ag alloys with CALPHAD, extending convergence tests as a function of interface width to 3D, and carrying out simulations in both 2D and 3D to examine existing theories of microstructure development. Our results indicate that the simplest variational formulation that interpolates the bulk free-energy density between its solid and liquid forms is the most robust. Remarkably, for hypoeutectic Al-Ag alloys, this formulation yields a high-velocity nonequilibrium phase diagram that is independent of interface width, thereby demonstrating that the framework of enhanced solute diffusivity can be nontrivially extended to concentrated alloys. Other variational formulations have restricted ranges of materials or processing parameters that can be reliably modeled. We use 2D simulations to construct high-velocity microstructure selection maps for dilute Al-Cu alloys. The results validate the important role of latent heat rejection at the interface and extend the limited predictions of linear stability analysis [A. Karma and A. Sarkissian, Phys. Rev. E 47, 513 (1993)] and sharp-interface 1D simulations to fully nonlinear regimes. Furthermore, 3D simulations, carried out using a computationally tractable axisymmetric cellular/dendritic interface shape, demonstrate a good convergence similar to that observed in 2D as a function of interface width. Full 3D simulations, in turn, reveal that the standard theory of absolute stability is a good predictor of the upper critical velocity beyond which steady-state growth becomes unstable, despite the different morphological manifestations of this instability in 2D and 3D.

36 MATERIALS SCIENCE↗
NASA NTRS1990

The effect of an electric field on the morphological stability of the crystal-melt interface of a binary alloy. III - Weakly nonlinear theory

The effect of a constant electric current on the crystal-melt interface morphology during directional solidification at constant velocity of a binary alloy is considered. A linear temperature field is assumed, and thermoelectric effects and Joule heating are neglected; electromigration and differing electrical conductivities of crystal and melt are taken into account. A two-dimensional weakly nonlinear analysis is carried out to third order in the interface amplitude, resulting in a cubic amplitude equation that describes whether the bifurcation from the planar state is supercritical or subcritical. For wavelengths corresponding to the most dangerous mode of linear theory, the demarcation between supercritical and subcritical behavior is calculated as a function of processing conditions and material parameters. The bifurcation behavior is a sensitive function of the magnitude and direction of the electric current and of the electrical conductivity ratio.

Wheeler, A. A.↗
DOE OSTI2025

Sharp front tracking with geometric interface reconstruction

Here, this paper presents a novel sharp front-tracking method designed to address limitations in classical front-tracking approaches, specifically their reliance on smooth interpolation kernels and extended stencils for coupling the front and fluid mesh. In contrast, the proposed method employs exclusively sharp, localized interpolation and spreading kernels, restricting the coupling to the interfacial fluid cells–those containing the interface/front. This localized coupling is achieved by integrating a divergence-preserving velocity interpolation method with a piecewise parabolic interface calculation (PPIC) and a polyhedron intersection algorithm to compute the indicator function and local interface curvature. Surface tension is computed using the Continuum Surface Force (CSF) method, maintaining consistency with the sharp representation. Additionally, we propose an efficient local roughness smoothing implementation to account for surface mesh undulations, which is easily applicable to any triangulated surface mesh. Building on our previous work, the primary innovation of this study lies in the localization of the coupling for both the indicator function and surface tension calculations. By reducing the interface thickness on the fluid mesh to a single cell, as opposed to the 4–5 cell spans typical in classical methods, the proposed sharp front-tracking method achieves a highly localized and accurate representation of the interface. This sharper representation mitigates parasitic currents and improves force balancing, making it particularly suitable for scenarios where the interface plays a critical role, such as microfluidics, fluid-fluid interactions, and fluid-structure interactions. The proposed method is comprehensively validated and tested on canonical interfacial flow problems, including stationary and translating Laplace equilibria, oscillating droplets, and rising bubbles. The presented results demonstrate that the sharp front-tracking method significantly outperforms the classical approach in terms of accuracy, stability, and computational efficiency. Notably, parasitic currents are reduced by approximately two orders of magnitude and stable results are obtained for parameter ranges where classical front tracking fails to converge.

42 ENGINEERING↗
DOE OSTI2023

Human-in-the-Loop Motion Control of a Two-DOF Hydraulic Backhoe Powered by the Hybrid Hydraulic Electric Architecture (HHEA)

Abstract The Hybrid Hydraulic-Electric Architecture (HHEA) combines the respective power density and control advantages of hydraulic and electric actuation to save energy for off-road vehicles. It uses a set of selectable common pressure rails to transmit the majority of power and electric actuation to modulate that power. As it is critical that off-road vehicles can perform tasks dexterously and exactly as commanded by the operator, the switchings between discrete pressure rails pose a potential challenge for smooth and precise motion. A control strategy consisting of a backstepping nominal control and least norm transition control has previously been developed to address this issue. It has been tested on 1 degree of freedom (DOF) testbeds where known trajectories were able to be tracked precisely. This paper presents the implementation of the HHEA motion control strategy on a 2-DOF backhoe operated by a human operator via a 2-DOF joystick. Unlike previous studies, the duty cycle is unknown beforehand and the decision to change pressure rails is taken in real-time. The efficacy of the motion control strategy has been validated experimentally. Several strategies to improve the user interface: control in workspace coordinates, pressure feedback, and velocity field-based task specification, have also been implemented and demonstrated to make operating the multiple DOF, HHEA actuated machine more intuitive to novice operators.

24 POWER TRANSMISSION AND DISTRIBUTION↗