Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Moving interface”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Nonlinear Response of Layer Growth Dynamics in the Mixed Kinetics-Bulk-Transport Regime

In situ high-resolution interferometry on horizontal facets of the protein lysozyme reveal that the local growth rate R, vicinal slope p, and tangential (step) velocity v fluctuate by up to 80% of their average values. The time scale of these fluctuations, which occur under steady bulk transport conditions through the formation and decay of step bunches (macrosteps), is of the order of 10 min. The fluctuation amplitude of R increases with growth rate (supersaturation) and crystal size, while the amplitude of the v and p fluctuations changes relatively little. Based on a stability analysis for equidistant step trains in the mixed transport-interface-kinetics regime, we argue that the fluctuations originate from the coupling of bulk transport with nonlinear interface kinetics. Furthermore, step bunches moving across the interface in the direction of or opposite to the buoyancy-driven convective flow increase or decrease in height, respectively. This is in agreement with analytical treatments of the interaction of moving steps with solution flow. Major excursions in growth rate are associated with the formation of lattice defects (striations). We show that, in general, the system-dependent kinetic Peclet number, Pe(sub k) , i.e., the relative weight of bulk transport and interface kinetics in the control of the growth process, governs the step bunching dynamics. Since Pe(sub k) can be modified by either forced solution flow or suppression of buoyancy-driven convection under reduced gravity, this model provides a rationale for the choice of specific transport conditions to minimize the formation of compositional inhomogeneities under steady bulk nutrient crystallization conditions.

Vekilov, Peter G.↗
DOE OSTI2025

An end-to-end deep learning method for solving nonlocal Allen–Cahn and Cahn–Hilliard phase-field models

Here, we propose an efficient end-to-end deep learning method for solving nonlocal Allen–Cahn (AC) and Cahn–Hilliard (CH) phase-field models. One motivation for this effort emanates from the fact that discretized partial differential equation-based AC or CH phase-field models result in diffuse interfaces between phases, with the only recourse for remediation is to severely refine the spatial grids in the vicinity of the true moving sharp interface whose width is determined by a grid-independent parameter that is substantially larger than the local grid size. In this work, we introduce non-mass conserving nonlocal AC or CH phase-field models with regular, logarithmic, or obstacle double-well potentials. Because of non-locality, some of these models feature totally sharp interfaces separating phases. The discretization of such models can lead to a transition between phases whose width is only a single grid cell wide. Another motivation is to use deep learning approaches to ameliorate the otherwise high cost of solving discretized nonlocal phase-field models. To this end, loss functions of the customized neural networks are defined using the residual of the fully discrete approximations of the AC or CH models, which results from applying a Fourier collocation method and a temporal semi-implicit approximation. To address the long-range interactions in the models, we tailor the architecture of the neural network by incorporating a nonlocal kernel as an input channel to the neural network model. We then provide the results of extensive computational experiments to illustrate the accuracy, predictive capabilities, and cost reductions of the proposed method.

42 ENGINEERING↗
NASA NTRSDate not supplied

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment1. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗
NASA NTRS1998

Experiment 13: The Study of Dopant Segregation Behavior During the Growth of GaAs in Microgravity on USML-2

An investigation into the segregation behavior of selenium doped gallium arsenide (Se/GaAs) during directional solidification in the microgravity environment was conducted using the Crystal Growth Furnace (CGF) aboard the second United States Microgravity Laboratory (USML-2). Two crystals were successfully processed on USML-2, which lasted from October 20 to November 7, 1995. The first sample was processed for 67 hours, 45 minutes (MET 5/04:53:45-8/00:23:50) and included 19 hours of growth at 0.5 microns/sec which yielded 3.42 cm of sample length, and 5 hours of growth at 1.5 microns/sec which yielded 2.7 cm of sample. During the second experiment, the furnace temperature was adjusted to move the melt-solid interface position towards the hot end of the furnace. The second sample was processed for 50 hours, 10 minutes (MET 8/18:48:49-10/21:58:54) and included 11 hours of growth at 0.5 microns/sec which yielded 1.98 cm of sample, and 1 hour, 25 minutes of growth at 5.0 microns/sec which yielded 2.6 cm of sample. This sample provides an order of magnitude change in growth rate and reproduces one of the growth rates used during USML-1. In contrast to the results from USML-1, no voids were present in either crystal grown on USML-2. The absence of voids in either sample indicates that growth rate changes alone were not responsible for the formation of voids found in the crystals grown on USML-1. Sections of the ground-based and flight crystals grown on USML-2 were cut and polished. All of the interface demarcation lines expected from the current pulse interface demarcation (CPID) system have been identified. These measurements have been analyzed for interface positions, interface shapes, and growth rates. Using a newly developed technique, based on experimental and numerical results, the seeding interface reproducibility from run to run was <= 2.5 mm. The seeding interface position could be controllably moved, with respect to the furnace zones, by adjusting the control set points of the heating zones. The interface shapes flattened slightly as the interface position moved closer to the hot zone but was always an unfavorable concave into the solid shape. The growth rate was found to equal the furnace translation rate, after a 2 -hour transient, for growth rates <= 1.0 microns/sec. Segregation measurements for the ground-based crystals are indicative of complete mixing behavior, as expected. Segregation measurements of the flight crystals are still in progress.

Matthiesen, David H.↗
NASA NTRS1991

Qualitative investigation of cryogenic injected shock dissipation

The injection of fluid nitrogen into the flow discharged from a two-dimensional, ambient, static-temperature and static pressure Mach = 2.7 nitrogen flow tunnel was observed. A bow shock stands approximately one diameter off the end of the 3.2 mm o.d. x 1.6 mm i.d. injection tube placed directly into the tunnel exhaust. Cryogenic injection creates a high density region in the injected region and within the bow shock wake but the standoff distance remains unchanged. However as the temperature reaches a critical value the sharp shock interface begins to fade resembling the density gradient seen at the interface of a supercritical fluid and the shock interface begins to move into the supersonic stream. As the injection temperature decreases, the interface continues to move into the supersonic stream and becomes more diffuse; the shock interface can not be distinguished. The phenomena is completely reversible.

Hendricks, R. C.↗
NASA NTRS2005

Observation of an Aligned Gas - Solid "Eutectic" during Controlled Directional Solidification Aboard the International Space Station - Comparison with Ground-based Studies

Direct observation of the controlled melting and solidification of succinonitrile was conducted in the glovebox facility of the International Space Station (ISS). The experimental samples were prepared on ground by filling glass tubes, 1 cm ID and approximately 30 cm in length, with pure succinonitrile (SCN) in an atmosphere of nitrogen at 450 millibar pressure for eventual processing in the Pore Formation and Mobility Investigation (PFMI) apparatus in the glovebox facility (GBX) on board the ISS. Real time visualization during controlled directional melt back of the sample showed nitrogen bubbles emerging from the interface and moving through the liquid up the imposed temperature gradient. Over a period of time these bubbles disappear by dissolving into the melt. Translation is stopped after melting back of about 9 cm of the sample, with an equilibrium solid-liquid interface established. During controlled re-solidification, aligned tubes of gas were seen growing perpendicular to the planar solid/liquid interface, inferring that the nitrogen previously dissolved into the liquid SCN was now coming out at the solid/liquid interface and forming the little studied liquid = solid + gas eutectic-type reaction. The observed structure is evaluated in terms of spacing dimensions, interface undercooling, and mechanisms for spacing adjustments. Finally, the significance of processing in a microgravity environment is ascertained in view of ground-based results.

Grugel, R. N.↗
NASA NTRS1999

A Real Time Investigation of Morphological Evolution During Solidification of Different Alloy Systems

Solidification phenomenon which occur at the solid/liquid (s/I) interface play a major role in the determination of structure and hence the technologically important properties of a casting. However, metals being opaque, conclusions related to several important phenomenon such as boundary layer thickness, morphological evolution, and eutectic and cell spacing are deduced from quenching experiments and subsequent post solidification metallographic analysis. Consequently, limited information is obtained about the dynamics of the process. This paper will discuss the recent efforts at the Space Science Laboratory, NASA Marshall Space Flight Center, to view and quantify in-situ and in real time the dynamics of the solidification process and to measure interfacial undercooling. First, a high resolution x-ray transmission microscope (XTM) has been developed to monitor fundamental interfacial phenomena during directional solidification of metals and alloys. The XTM operates in the range of 10-100 KeV and through projection is capable of achieving magnification of up to 16OX. Secondly, an innovative collapsible furnace has been designed to quantify interfacial undercooling by measuring the temperature of a moving s/I interface in reference to a fixed s/l interface. This measurement technique is non-intrusive in nature and is based on the Seebeck principle. In this paper real time results obtained to characterize the dynamics of irregular eutectic spacing will be presented. As an example fiber to lamella or plate transition in the Al-Al2Au eutectic system will be discussed. Further, a resolution limit of 25 micron has permitted viewing in real time morphological instability and cellular growth in Al-Au and Al-Ag systems. Simultaneously, a systematic investigation has been carried out to measure interfacial undercooling for Pb-1 wt.% Sn at and near the marginal stability regime. In conjunction with the XTM observations this study attempts to validate existing relationships between undercooling and growth velocity during plane front growth, marginal stability regime, and stable cellular growth.

Sen, S.↗
DOE OSTI2025

Predictive dynamic wetting, fluid–structure interaction simulations for braze run-out

Brazing and soldering are metallurgical joining techniques that use a wetting molten metal to create a joint between two faying surfaces. Here, the quality of the brazing process depends strongly on the wetting properties of the molten filler metal, namely the surface tension and contact angle, and the resulting joint can be susceptible to various defects, such as run-out and underfill, if the material properties or joining conditions are not suitable. In this work, we implement a finite element simulation to predict the formation of such defects in braze processes. This model incorporates both fluid–structure interaction through an arbitrary Eulerian–Lagrangian technique and free surface wetting through conformal decomposition finite element modeling. Upon validating our numerical simulations against experimental run-out studies on a silver-Kovar system, we then use the model to predict run-out and underfill in systems with variable surface tension, contact angles, and applied pressure. Finally, we consider variable joint/surface geometries and show how different geometrical configurations can help to mitigate run-out. This work aims to understand how brazing defects arise and validate a coupled wetting and fluid–structure interaction simulation that can be used for other industrial problems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A fully-integrated lattice Boltzmann method for fluid–structure interaction

Here we present a fully-integrated lattice Boltzmann (LB) method for fluid–structure interaction (FSI) simulations that efficiently models deformable solids in complex suspensions and active systems. Our Eulerian method (LBRMT) couples finite-strain solids to the LB fluid on the same fixed computational grid with the reference map technique (RMT). An integral part of the LBRMT is a new LB boundary condition for moving deformable interfaces across different densities. With this fully Eulerian solid–fluid coupling, the LBRMT is well-suited for parallelization and simulating multi-body contact without remeshing or extra meshes. We validate its accuracy via a benchmark of a deformable solid in a lid-driven cavity, then showcase its versatility through examples of soft solids rotating and settling. The LBRMT achieves a spatial convergence rate between first-order and second-order for FSI simulations and is designed for low to intermediate Reynolds number flows with finite inertia at small Mach numbers. With simulations of complex suspensions mixing, we highlight the potential of the LBRMT for studying collective behavior in soft matter and biofluid dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Optical and electronic functionality arising from controlled defect formation in nanoscale complex oxide lateral epitaxy

Epitaxial crystallization of complex oxides provides the means to create materials with precisely selected composition, strain, and orientation, thereby controlling their functionalities. Extending this control to nanoscale three-dimensional geometries can be accomplished via a three-dimensional analog of oxide solid-phase epitaxy, lateral epitaxial crystallization. The orientation of crystals within laterally crystallized SrTiO 3 systematically changes from the orientation of the SrTiO 3 substrate. This evolution occurs as a function of lateral crystallization distance, with a rate of approximately 50° μm -1 . The mechanism of the rotation is consistent with a steady-state stress of tens of megapascal over a 100–nanometer scale region near the moving amorphous/crystalline interface arising from the amorphous-crystalline density difference. Second harmonic generation and piezoelectric force microscopy reveal that the laterally crystallized SrTiO 3 is noncentrosymmetric and develops a switchable piezoelectric response at room temperature, illustrating the potential to use lateral crystallization to control the functionality of complex oxides.

36 MATERIALS SCIENCE↗
NASA NTRS1974

Long Life Technology Work at Rockwell International Space Division

This paper presents highlights of long-life technology oriented work performed at the Space Division of Rockwell International Corporation under contract to NASA. This effort included evaluation of Saturn V launch vehicle mechanical and electromechanical components for potential extended life capabilities, endurance tests, and accelerated aging experiments. A major aspect was evaluation of the components at the subassembly level (i.e., at the interface between moving surfaces) through in-depth wear analyses and assessments. Although some of this work is still in progress, preliminary conclusions are drawn and presented, together with the rationale for each. The paper concludes with a summary of the effort still remaining.

D K Huzel↗
NASA NTRS1985

Finite element simulation of temperature dependent free surface flows

The method of Engelman and Sani (1984) for a finite-element simulation of incompressible surface flows with a free and/or moving fluid interface, such as encountered in crystal growth and coating and polymer technology, is extended to temperature-dependent flows, including the effect of temperature-dependent surface tension. The basic algorithm of Saito and Scriven (1981) and Ruschak (1980) has been generalized and implemented in a robust and versatile finite-element code that can be employed with relative ease for the simulation of free-surface problems in complex geometries. As a result, the costly dependence on the Newton-Raphson algorithm has been eliminated by replacing it with a quasi-Newton iterative method, which nearly retains the superior convergence properties of the Newton-Raphson method.

Engelman, M. S.↗
NASA NTRS1993

Numerical study of 3-D inducer and impeller for pump model development

Current design of high-performance turbopumps for rocket engines requires effective and robust analytical tools to provide design information in a productive manner. The main goal of this study is to develop a robust and effective CFD pump model for general turbopump design and analysis applications. A finite difference Navier-Stokes flow solver, FDNS, which includes an extended k-epsilon turbulence model and appropriate moving zonal interface boundary conditions, was developed to analyze turbulent flows in turbomachinery devices. In the present study, two key components of the turbopump, the inducer and impeller, were investigated by the proposed pump model, and the numerical results were benchmarked by the experimental data provided by Rocketdyne.

Cheng, G. C.↗
NASA NTRS1994

Inducer analysis/pump model development

Current design of high performance turbopumps for rocket engines requires effective and robust analytical tools to provide design information in a productive manner. The main goal of this study was to develop a robust and effective computational fluid dynamics (CFD) pump model for general turbopump design and analysis applications. A finite difference Navier-Stokes flow solver, FDNS, which includes an extended k-epsilon turbulence model and appropriate moving zonal interface boundary conditions, was developed to analyze turbulent flows in turbomachinery devices. In the present study, three key components of the turbopump, the inducer, impeller, and diffuser, were investigated by the proposed pump model, and the numerical results were benchmarked by the experimental data provided by Rocketdyne. For the numerical calculation of inducer flows with tip clearance, the turbulence model and grid spacing are very important. Meanwhile, the development of the cross-stream secondary flow, generated by curved blade passage and the flow through tip leakage, has a strong effect on the inducer flow. Hence, the prediction of the inducer performance critically depends on whether the numerical scheme of the pump model can simulate the secondary flow pattern accurately or not. The impeller and diffuser, however, are dominated by pressure-driven flows such that the effects of turbulence model and grid spacing (except near leading and trailing edges of blades) are less sensitive. The present CFD pump model has been proved to be an efficient and robust analytical tool for pump design due to its very compact numerical structure (requiring small memory), fast turnaround computing time, and versatility for different geometries.

Cheng, Gary C.↗
NASA NTRS1995

Analytical studies of hypersonic viscous dissociated flows

This project primarily dealt with integral boundary-layer solution techniques that are directly applicable to the problem of determining aerodynamic heating rates of hypersonic vehicles like X-33 in the vicinity of stagnation points, windward centerlines, and swept-wing leading edges. The analyses include effects of finite-rate gas chemistry across the boundary layer and finite-rate catalysis of atom recombination at the surface. A new approach for combining the insight afforded by integral boundary-layer analysis with comprehensive (and expensive) computational fluid dynamic (CFD) flowfield solutions of the thin-layer Navier-Stokes equations was developed. The approach extracts CFD derived quantities at the wall and at the boundary layer edge for inclusion in a post-processing boundary-layer analysis. The post-processed data base allows a designer at a workstation to ask and answer the following questions: (1) How much does the heating change if one uses a thermal protection system (TPS) with different catalytic properties than was used in the original CFD solution? (2) How does the heating change when one moves the interface of two different TPS materials with different catalytic efficiencies for the purpose of reducing vehicle weight and expense? The answer to the second question is particularly critical, because abrupt changes from low catalytic efficiency to high catalytic efficiency can lead to localized increase in heating which exceeds the usually conservative estimate provided by a fully catalytic wall assumption. A secondary issue that was addressed involves the prediction of heating levels in the vicinity of sharp corners that are transverse to or aligned with the flow. An example of the first case is heating at the edge of the COMET reentry module. An example of the second case is heating along the side edge of a deflected body flap on an SSV. The difficulty of putting grids in the vicinity of such corners with continuously varying metric coefficients causes problems in CFD predictions. A preliminary theory for prediction that says the heating at the corner is X percent of the heating N boundary-layer thicknesses inboard was developed. This will prove useful to analytically evaluate the possible benefits of rounding the edges of these configurations and defining how much rounding is sufficient.

Inger, George R.↗
NASA NTRS2000

Numerical Calculation of the Drag Force Acting on an Insoluble Particle Moving in Front of a Solidifying Interface

The distribution of insoluble particles in a metal casting depends primarily on the interaction of the particles with the solid/liquid interface (SLI) during the solidification process. Whether a particle will be engulfed or continuously pushed by SLI is essentially determined by the balance of forces acting on the particle. An important component of this force balance is the drag force generated by the particle motion in front of the SLI. Previously developed mathematical models for particle/SLI interaction made use of steady-state solutions of this force provided by the lubrication theory. However, our numerical model based on the SLI tracking approach shows that not only the steady-state approach is inappropriate to model the interaction process but also that even at steady-state the theoretical solution underestimates the drag force. A regression analysis of steady-state numerical solutions for cylindrical particles moving normal to a flat SLI gave a relationship of the form F(sub d) (exp num) = (square root of 3) pi eta V(sub p)((R(sub p) / d) (exp 10 gamma/3)) as compared to the theoretical solution F(sub d) (exp theor) = 3 X (square root of 2) pi eta V(sub p)((R(sub p) / d) (exp 3/2)) where F(sub d) is the drag force, eta is the dynamic viscosity of the fluid, V(sub p) is the particle velocity, R(sub p) is the particle radius, gamma is Euler's constant (gamma is approximately equal to 0.577), and d is the width of the gap between the particle and the SLI. Model validation as well as the influence of the interface shape on the value of F(sub d) will also be discussed.

Catalina, Adrian V.↗
NASA NTRS1995

Advanced Payload Transfer Measurement System (APT-MS) Mechanical Features

During the last two feet or so of transfer, for a large or heavy payload to its restraining fixture, the consequences of unplanned contact (or impact) between payload and support structure can range from merely annoying to something approaching disaster. Current transfer methods employ technicians with meter sticks stationed at the various hold-down locations to estimate the distances to contact. This information is communicated orally to the crane operator. It is understandable that this work proceeds carefully and therefore slowly. The objectives of this study are (1) to explore mechanical sources of measurement error; and (2) to develop an improved prototype design that is compact, inexpensive, and addresses the measurement error issues. APT is a relatively inexpensive electro-mechanical device that reduces both time and manpower required to make shuttle payload transfers. APT measurement system would provide a GUI for the 'move-conductor' (crane operator) so he could see the displacements of all hold-down interfaces as they move together and thus have a more accurate, comprehensive and 'real-time' picture of the engagement activity. An error model attempts to include all estimatable sources of mechanical error and design features were introduced to reduce or eliminate major sources of error.

Brewer, William V.↗