Engineering PapersSearch

SEARCH · Engineering Papers

Results for “phase field modeling”

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 19 records

Phase-field model for isothermal phase transitions in binary alloys

A new phase field model is described which models isothermal phase transitions between ideal binary alloy solution phases. Equations are developed for the temporal and spatial variation of the phase field, which describes the identity of the phase, and of the composition. An asymptotic analysis, as the gradient energy coefficient of the phase field becomes small, was conducted. From the analysis, it is shown that the model recovers classical sharp interface models of this situation when the interfacial layers are thin, and they show how to relate the parameters appearing in the phase field model to material and growth parameters in real systems. Further, three stages of temporal evolution are identified: the first corresponding to interfacial genesis which occurs very rapidly; the second to interfacial motion controlled by the local energy difference across the interface and diffusion; the last taking place on a long time scale in which curvature effects are important and which correspond to Ostwald ripening. The results of the numerical calculations are presented.

Wheeler, A. A.

Adaptive-Grid Methods for Phase Field Models of Microstructure Development

In this work the authors show how the phase field model can be solved in a computationally efficient manner that opens a new large-scale simulational window on solidification physics. Our method uses a finite element, adaptive-grid formulation, and exploits the fact that the phase and temperature fields vary significantly only near the interface. We illustrate how our method allows efficient simulation of phase-field models in very large systems, and verify the predictions of solvability theory at intermediate undercooling. We then present new results at low undercoolings that suggest that solvability theory may not give the correct tip speed in that regime. We model solidification using the phase-field model used by Karma and Rappel.

Provatas, Nikolas

Analytical prediction of the interior noise for cylindrical models of aircraft fuselages for prescribed exterior noise fields. Phase 2: Models for sidewall trim, stiffened structures and cabin acoustics with floor partition

An airplane interior noise prediction model is developed to determine the important parameters associated with sound transmission into the interiors of airplanes, and to identify apropriate noise control methods. Models for stiffened structures, and cabin acoustics with floor partition are developed. Validation studies are undertaken using three test articles: a ring stringer stiffened cylinder, an unstiffened cylinder with floor partition, and ring stringer stiffened cylinder with floor partition and sidewall trim. The noise reductions of the three test articles are computed using the heoretical models and compared to measured values. A statistical analysis of the comparison data indicates that there is no bias in the predictions although a substantial random error exists so that a discrepancy of more than five or six dB can be expected for about one out of three predictions.

Pope, L. D.

Use of Microgravity to Control the Microstructure of Eutectics

This grant began in June of 1996. Its long term goal is to be able to control the microstructure of directionally solidified eutectic alloys, through an improved understanding of the influence of convection. The primary objective of the present projects is to test hypotheses for the reported influence of microgravity on the microstructure of three fibrous eutectics (MnBi-Bi, InSb-NiSb, Al3Ni-Al). A secondary objective is to determine the influence of convection on the microstructure of other eutectic alloys. Two doctoral students and a masters student supported as a teaching assistant were recruited for this research. Techniques were developed for directional solidification of MnBi-Bi eutectics with periodic application of current pulses to produce an oscillatory freezing rate. Image analysis techniques were developed to obtain the variation in MnBi fiber spacing, which was found to be normally distributed. The mean and standard deviation of fiber spacing were obtained for several freezing conditions. Eighteen ampoules were prepared for use in the gradient freeze furnace QUELD developed at Queen's University for use in microgravity. Nine of these ampoules will be solidified soon at Queen's in a ground-based model. We hope to solidify the other nine in the QUELD that is mounted on the Canadian Microgravity Isolation Mount on MIR. Techniques are being developed for directional solidification of the Al-Si eutectic at different freezing rates, with and without application of accelerated crucible rotation to induce convection. For the first time, theoretical methods are being developed to analyze eutectic solidification with an oscillatory freezing rate. In a classical sharp-interface model, we found that an oscillatory freezing rate increases the deviation of the average interfacial composition from the eutectic, and increases the undercooling of the two phases by different amounts. This would be expected to change the volume fraction solidifying and the fiber spacing. Because of difficulties in tracking the freezing interfaces of the two solid phases, a phase-field model is also being developed. A paper demonstrating application of phase field methods to periodic structures has been submitted for publication.

Wilcox, William R.

Analytical prediction of the interior noise for cylindrical models of aircraft fuselages for prescribed exterior noise fields. Phase 1: Development and validation of preliminary analytical models

The basic theoretical work required to understand sound transmission into an enclosed space (that is, one closed by the transmitting structure) is developed for random pressure fields and for harmonic (tonal) excitation. The analysis is used to predict the noise reducton of unpressurized unstiffened cylinder, and also the interior response of the cylinder given a tonal (plane wave) excitation. Predictions and measurements are compared and the transmission is analyzed. In addition, results for tonal (harmonic) mechanical excitation are considered.

Pope, L. D.

Summary of nozzle-exhaust plume flowfield analyses performed in support of space shuttle applications

Results are summarized of various analyses of the exhaust plumes of the space shuttle. Specific topics discussed included: the development of the two-phase plume flow field model, including finite rate chemistry and free molecular effects; analyses of exhaust plume simulation studies; verification of the analytical two-phase plume flow field model; and complete documentation of the two-phase plume computer code.

Smith, S. D.

Advanced Multi-Phase Flow CFD Model Development for Solid Rocket Motor Flowfield Analysis

It is known that the simulations of solid rocket motor internal flow field with AL-based propellants require complex multi-phase turbulent flow model. The objective of this study is to develop an advanced particulate multi-phase flow model which includes the effects of particle dynamics, chemical reaction and hot gas flow turbulence. The inclusion of particle agglomeration, particle/gas reaction and mass transfer, particle collision, coalescence and breakup mechanisms in modeling the particle dynamics will allow the proposed model to realistically simulate the flowfield inside a solid rocket motor. The Finite Difference Navier-Stokes numerical code FDNS is used to simulate the steady-state multi-phase particulate flow field for a 3-zone 2-D axisymmetric ASRM model and a 6-zone 3-D ASRM model at launch conditions. The 2-D model includes aft-end cavity and submerged nozzle. The 3-D model represents the whole ASRM geometry, including additional grain port area in the gas cavity and two inhibitors. FDNS is a pressure based finite difference Navier-Stokes flow solver with time-accurate adaptive second-order upwind schemes, standard and extended k-epsilon models with compressibility corrections, multi zone body-fitted formulations, and turbulence particle interaction model. Eulerian/Lagrangian multi-phase solution method is applied for multi-zone mesh. To simulate the chemical reaction, penalty function corrected efficient finite-rate chemistry integration method is used in FDNS. For the AL particle combustion rate, the Hermsen correlation is employed. To simulate the turbulent dispersion of particles, the Gaussian probability distribution with standard deviation equal to (2k/3)(exp 1/2) is used for the random turbulent velocity components. The computational results reveal that the flow field near the juncture of aft-end cavity and the submerged nozzle is very complex. The effects of the turbulent particles affect the flow field significantly and provide better prediction of the ASRM performance. The multi-phase flow analysis using the FDNS code in the present research can be used as a design tool for solid rocket motor applications.

Liaw, Paul

Study of Nonclassical Fields in Phase-Sensitive Reservoirs

We show that the reservoir influence can be modeled by an infinite array of beam splitters. The superposition of the input fields in the beam splitter is discussed with the convolution laws for their quasiprobabilities. We derive the Fokker-Planck equation for the cavity field coupled with a phase-sensitive reservoir using the convolution law. We also analyze the amplification in the phase-sensitive reservoir with use of the modified beam splitter model. We show the similarities and differences between the dissipation and amplification models. We show that a super-Poissonian input field cannot become sub-Poissonian by the phase-sensitive amplification.

Kim, Myung Shik

Prompt particle acceleration around moving X-point magnetic field during impulsive phase of solar flares

We present a model for high-energy solar flares to explain prompt proton and electron acceleration, which occurs around moving X-point magnetic field during the implosion phase of the current sheet. We derive the electromagnetic fields during the strong implosion phase of the current sheets, which is driven by the converging flow derived from the magnetohydrodynamic equations. It is shown that both protons and electrons can be promptly (within 1 second) accelerated to approximately 70 MeV and approximately 200 MeV, respectively. This acceleration mechanism can be applicable for the impulsive phase of the gradual gamma ray and proton flares (gradual GR/P flare), which have been called two-ribbon flares.

Sakai, Jun-Ichi

Particle scattering and current sheet stability in the geomagnetic tail during the substorm growth phase

The particle scattering and current sheet stability features in the geomagnetic tail during the phase of substorm growth were investigated using Tsyganenko's (1989) magnetic field model. In a study of four substorm events which were observed both in the high-altitude nightside tail and in the auroral ionosphere, the model magnetic field was adjusted to each case so as to represent the global field development during the growth phase of the substorms. The model results suggest that the auroral brightenings are connected with processes taking place in the near-earth region inside about 15 earth radii. The results also suggest that there is a connection between the chaotization of the electrons and the auroral brightenings at substorm onset.

Pulkkinen, T. I.

Applicability of empirical data currently used in predicting solid propellant exhaust plumes

Theoretical and experimental approaches to exhaust plume analysis are compared. A two-phase model is extended to include treatment of reacting gas chemistry, and thermodynamical modeling of the gaseous phase of the flow field is considered. The applicability of empirical data currently available to define particle drag coefficients, heat transfer coefficients, mean particle size, and particle size distributions is investigated. Experimental and analytical comparisons are presented for subscale solid rocket motors operating at three altitudes with attention to pitot total pressure and stagnation point heating rate measurements. The mathematical treatment input requirements are explained. The two-phase flow field solution adequately predicts gasdynamic properties in the inviscid portion of two-phase exhaust plumes. It is found that prediction of exhaust plume gas pressures requires an adequate model of flow field dynamics.

Tevepaugh, J. A.

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record

Preliminary calibration plan for the Advanced Particles and Field Observatory (APAFO) magnetometer experiment

Prelaunch and postlaunch calibration plans for the APAFO magnetometer experiment are presented. A study of tradeoffs between boom length and spacecraft field is described; the results are summarized. The prelaunch plan includes: calibration of the Vector Fluxgate Magnetometer (VFM), Star Sensors, and Scalar Helium Magnetometer (SHM); optical bench integration; and acquisition of basic spacecraft field data. Postlaunch calibration has two phases. In phase one, SHM data are used to calibrate the VFM, total vector magnetic field data are used to calibrate a physical model of the spacecraft field, and both calibrations are refined by iteration. In phase two, corrected vector data are transformed into geocentric coordinates, previously undetected spacecraft fields are isolated, and initial geomagnetic field models are computed. Provided the SHM is accurate to the required 1.0 nT and can be used to calibrate the VFM to the required 3.0- nT accuracy, the tradeoff study indicates that a 12 m boom and a spacecraft field model uncertainty of 5 percent together allow the 1.0 nT spacecraft field error requirement to be met.

Voorhies, C. V.

Ordering and phase transitions in random-field Ising systems

An exact analysis of the Ising model with infinite-range interactions in a random field and a local mean-field theory in three dimensions is carried out leading to a phase diagram with several coexistence surfaces and lines of critical points. The results show that the phase diagram depends crucially on whether the distribution of random fields is symmetric or not. Thus, Ising-like phase transitions in a porous medium (the asymmetric case) are in a different universality class from the conventional random-field model (symmetric case).

Maritan, Amos

Novel local symmetries and chiral-symmetry-broken phases in S = 1/2 triangular-lattice Heisenberg model

Using a nonmean-field approach the triangular-lattice S = 1/2 Heisenberg antiferromagnet with nearest- and next-nearest-neighbor couplings is shown undergo an Ising-type phase transition into a chiral-symmetry-broken phase (Kalmeyer-Laughlin-like state) at small T. Removal of next-nearest-neighbor coupling introduces a local Z2 symmetry, thereby suppressing any finite-T chiral order.

Baskaran, G.

CDAW 9 analysis of magnetospheric events on May 3, 1986 - Event C

An intense geomagnetic substorm event on May 3, 1986, occurring toward the end of a strong storm period, is studied. The auroral electrojet indices and global imaging data from both the Northern and Southern Hemispheres clearly revealed the growth phase and expansion phase development for a substorm with an onset at 0111 UT. An ideally located constellation of four spacecraft allowed detailed observation of the substorm growth phase in the near-tail region. A realistic time-evolving magnetic field model provided a global representation of the field configuration throughout the growth and early expansion phase of the substorm. Evidence of a narrowly localized substorm onset region in the near-earth tail is found. This region spread rapidly eastward and poleward after the 0111 UT onset. The results are consistent with a model of late growth phase formation of a magnetic neutral line. This reconnection region caused plasma sheet current diversion before the substorm onset and eventually led to cross-tail current disruption at the time of the substorm onset.

Baker, D. N.

Assessment of analytical techniques for predicting solid propellant exhaust plumes

The calculation of solid propellant exhaust plume flow fields is addressed. Two major areas covered are: (1) the applicability of empirical data currently available to define particle drag coefficients, heat transfer coefficients, mean particle size and particle size distributions, and (2) thermochemical modeling of the gaseous phase of the flow field. Comparisons of experimentally measured and analytically predicted data are made. The experimental data were obtained for subscale solid propellant motors with aluminum loadings of 2, 10 and 15%. Analytical predictions were made using a fully coupled two-phase numerical solution. Data comparisons will be presented for radial distributions at plume axial stations of 5, 12, 16 and 20 diameters.

Tevepaugh, J. A.