Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase field method”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Phase-Field Methods for Structure Evolution in Sheared Multiphase Systems

A homogeneous disordered phase separates into ordered structures when quenched into a broken-symmetry phase. The competition of broken-symmetry phases to select an equilibrium state may be studied in terms of coarse-grained order parameters described by a suitable Landau free-energy function. A network of equilibrium-phase domains develops on quenching and coarsens with time with a topology that may be controlled by shear. We use three-dimensional simulations, in which time-dependent models for conserved-order parameters coupled to Navier-Stokes fluid models are solved, to investigate the evolution of such domains, e.g. spinodal decompositions of polymeric materials under shear. The numerical problems are formidable because of the strong nonlinearities inherent in the coupled model, and these are amongst the first 3D calculations undertaken. In linear shear fields we find stable nanostrings, also recently seen in experiments. The affinity of the ordered phases to boundaries plays a role in the form of the structures that develop, with stacked plate-like phase distributions emerging under certain conditions. Such methods appear quite promising for design and analysis of multiphase and complex fluid formulations. The behavior of foams in such conditions is of particular interest in microgravity environments. Additional information can be found in the original extended abstract.

Badalassi, Vittorio↗
NASA NTRS1999

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↗
NASA NTRS2001

Adaptive-Grid Methods for Phase Field Models of Microstructure Development

Modeling solidification microstructures has become an area of intense study in recent years. The properties of large scale cast products, ranging from automobile engine blocks to aircraft components and other industrial applications, are strongly dependent on the physics that occur at the mesoscopic and microscopic length scales during solidification. The predominant morphology found in solidification microstructures is the dendrite, a tree-like pattern of solid around which solidification proceeds. The microscopic properties of cast products are determined by the length scales of these dendrites, and their associated segregation profiles. For this reason understanding the mechanisms for pattern selection in dendritic growth has attracted a great deal of interest from the experimental and theoretical communities. In particular, a great deal of research has been undertaken to understand such issues as dendrite morphology, shape and growth speed. Experiments on dendrite evolution in pure materials by Glicksman and coworkers on succinonitrile (SCN), and more recently pivalic acid (PVA), as well as other transparent analogs of metals, have provided tests of theories for dendritic growth, and have stimulated considerable theoretical progress. These experiments have clearly demonstrated that in certain parameter ranges the physics of the dendrite tip can be characterized by a steady value for the dendrite tip velocity, radius of curvature and shape. Away from the tip, the time-dependent dendrite exhibits a characteristic sidebranching as it propagates, which is not yet well understood. These experiments are performed by observing individual dendrites growing into an undercooled melt. The experiments are characterized by the dimensionless undercooling. Most experiments are performed at low undercooling.

Dantzig, Jonathan A.↗
NASA NTRS2004

Numerical Simulation of Nanostructure Growth

Nanoscale structures, such as nanowires and carbon nanotubes (CNTs), are often grown in gaseous or plasma environments. Successful growth of these structures is defined by achieving a specified crystallinity or chirality, size or diameter, alignment, etc., which in turn depend on gas mixture ratios. pressure, flow rate, substrate temperature, and other operating conditions. To date, there has not been a rigorous growth model that addresses the specific concerns of crystalline nanowire growth, while demonstrating the correct trends of the processing conditions on growth rates. Most crystal growth models are based on the Burton, Cabrera, and Frank (BCF) method, where adatoms are incorporated into a growing crystal at surface steps or spirals. When the supersaturation of the vapor is high, islands nucleate to form steps, and these steps subsequently spread (grow). The overall bulk growth rate is determined by solving for the evolving motion of the steps. Our approach is to use a phase field model to simulate the growth of finite sized nanowire crystals, linking the free energy equation with the diffusion equation of the adatoms. The phase field method solves for an order parameter that defines the evolving steps in a concentration field. This eliminates the need for explicit front tracking/location, or complicated shadowing routines, both of which can be computationally expensive, particularly in higher dimensions. We will present results demonstrating the effect of process conditions, such as substrate temperature, vapor supersaturation, etc. on the evolving morphologies and overall growth rates of the nanostructures.

Hwang, Helen H.↗
NASA NTRS1998

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.↗
NASA NTRS1996

Numerical Investigation of Two-Phase Flows With Charged Droplets in Electrostatic Field

A numerical method to solve two-phase turbulent flows with charged droplets in an electrostatic field is presented. The ensemble-averaged Navier-Stokes equations and the electrostatic potential equation are solved using a finite volume method. The transitional turbulence field is described using multiple-time-scale turbulence equations. The equations of motion of droplets are solved using a Lagrangian particle tracking scheme, and the inter-phase momentum exchange is described by the Particle-In-Cell scheme. The electrostatic force caused by an applied electrical potential is calculated using the electrostatic field obtained by solving a Laplacian equation and the force exerted by charged droplets is calculated using the Coulombic force equation. The method is applied to solve electro-hydrodynamic sprays. The calculated droplet velocity distributions for droplet dispersions occurring in a stagnant surrounding are in good agreement with the measured data. For droplet dispersions occurring in a two-phase flow, the droplet trajectories are influenced by aerodynamic forces, the Coulombic force, and the applied electrostatic potential field.

Kim, Sang-Wook↗
NASA NTRS1991

Convection of ion cyclotron waves to ion-heating regions

Results are presented on calculations of the convection, to lower altitudes, of low-frequency ion cyclotron waves generated in the central plasma sheet. In the calculations, ion cyclotron waves generated in the equatorial plane by a proton temperature anisotropy are considered, and the computed growth rates are used to create a model wave distribution. Spectral densities are then mapped along the magnetic field lines, using phase space methods of Ronnmark and Larsson (1988). It was found that, even in the absence of growth associated with a plasma instability, the electric fields of ion cyclotron emissions from the equatorial region will increase substantially as the waves propagate to higher latitudes. It is shown that ion cyclotron waves propagating down the field lines may contribute to the waves observed together with ion conics.

Ronnmark, Kjell↗
NASA NTRS1973

Study of high altitude plume impingement

The radiation intensities are determined in the base region of the space shuttle due to solid particle radiation emanating from the solid rocket motors of the shuttle. Results of an analysis of the Titan 3 and simulated solid rocket motor radiation intensities are presented. The gas particle flow fields of the Titan 3 nozzle and plume and a space shuttle solid rocket motor nozzle and plume are described. The gaseous Titan 3 flow fields are discussed utilizing the results of flow fields generated by a gaseous and two phase method-of-characteristics computer programs. A two phase computer flow field analysis program was developed. An outflow correction theory is developed which will be used to modify existing convection heat transfer methods for better heat transfer predictions on bodies immersed in rocket exhaust plumes.

Mcanally, J. V.↗
NASA NTRS1982

Superalloy composition modeling

Development of a predictive method for determination of the gamma/gamma prime phase fields, i.e., gamma prime volume fraction as a function of the multicomponent composition, is described. The cluster variation method used for binary alloys in which the precipitated phase is coherent with the matrix phase is extended for application to the multicomponent coherent gamma/gamma prime nickel-based superalloys. It is shown that the cluster variation method can accurately describe the equilibrium (incoherent) gamma/gamma prime phase fields in the binary Ni-Al phase diagram. The gamma/gamma prime phase field for the Ni-Cr-Al ternary phase diagram is computed as a function of temperature. A reasonable fit results between the calculated and the experimental diagrams. The modeling of the six-component Ni-Cr-Al-Co-Mo-Ti base superalloy and the effect of Ni substitution of Co are discussed.

Barefoot, J.↗
NASA NTRS1985

Methods for heat transfer and temperature field analysis of the insulated diesel phase 2 progress report

This report describes work done during Phase 2 of a 3 year program aimed at developing a comprehensive heat transfer and thermal analysis methodology for design analysis of insulated diesel engines. The overall program addresses all the key heat transfer issues: (1) spatially and time-resolved convective and radiative in-cylinder heat transfer, (2) steady-state conduction in the overall structure, and (3) cyclical and load/speed temperature transients in the engine structure. During Phase 2, radiation heat transfer model was developed, which accounts for soot formation and burn up. A methodology was developed for carrying out the multi-dimensional finite-element heat conduction calculations within the framework of thermodynamic cycle codes. Studies were carried out using the integrated methodology to address key issues in low heat rejection engines. A wide ranging design analysis matrix was covered, including a variety of insulation strategies, recovery devices and base engine configurations. A single cylinder Cummins engine was installed at Purdue University, and it was brought to a full operational status. The development of instrumentation was continued, concentrating on radiation heat flux detector, total heat flux probe, and accurate pressure-crank angle data acquisition.

Morel, T.↗
NASA NTRS1988

Methods for heat transfer and temperature field analysis of the insulated diesel, phase 3

Work during Phase 3 of a program aimed at developing a comprehensive heat transfer and thermal analysis methodology for design analysis of insulated diesel engines is described. The overall program addresses all the key heat transfer issues: (1) spatially and time-resolved convective and radiative in-cylinder heat transfer, (2) steady-state conduction in the overall structure, and (3) cyclical and load/speed temperature transients in the engine structure. These are all accounted for in a coupled way together with cycle thermodynamics. This methodology was developed during Phases 1 and 2. During Phase 3, an experimental program was carried out to obtain data on heat transfer under cooled and insulated engine conditions and also to generate a database to validate the developed methodology. A single cylinder Cummins diesel engine was instrumented for instantaneous total heat flux and heat radiation measurements. Data were acquired over a wide range of operating conditions in two engine configurations. One was a cooled baseline. The other included ceramic coated components (0.050 inches plasma sprayed zirconia)-piston, head and valves. The experiments showed that the insulated engine has a smaller heat flux than the cooled one. The model predictions were found to be in very good agreement with the data.

Morel, Thomas↗
NASA NTRS1981

Exploitation of SAR data for measurement of ocean currents and wave velocities

Methods of extracting information on ocean currents and wave orbital velocities from SAR data by an analysis of the Doppler frequency content of the data are discussed. The theory and data analysis methods are discussed, and results are presented for both aircraft and satellite (SEASAT) data sets. A method of measuring the phase velocity of a gravity wave field is also described. This method uses the shift in position of the wave crests on two images generated from the same data set using two separate Doppler bands. Results of the current measurements are pesented for 11 aircraft data sets and 4 SEASAT data sets.

Shuchman, R. A.↗
NASA NTRS2009

Phased-array optical whispering gallery mode modulation and method

A whispering gallery mode (WGM) resonator device and method capable of filtering sidebands of optical modulators are provided. The method includes providing an optical resonator adapted to support whispering gallery modes and forming a first field and a second field from a first location and a second location, respectively, at the circumference of the optical resonator and being separated by an arc angle, .alpha.. The method includes adjusting relative phase between the first field and the second field in accordance to a differential phase, .beta., and combining the first and the second fields into an output. Particular selection of the arc angle, .alpha., and the differential phase, .beta., can determine the function of the output.

Strekalov, Dmitry V.↗
NASA NTRS2004

Modeling Growth of Nanostructures in Plasmas

As semiconductor circuits shrink to CDs below 0.1 nm, it is becoming increasingly critical to replace and/or enhance existing technology with nanoscale structures, such as nanowires for interconnects. Nanowires grown in plasmas are strongly dependent on processing conditions, such as gas composition and substrate temperature. Growth occurs at specific sites, or step-edges, with the bulk growth rate of the nanowires determined from the equation of motion of the nucleating crystalline steps. Traditional front-tracking algorithms, such as string-based or level set methods, suffer either from numerical complications in higher spatial dimensions, or from difficulties in incorporating surface-intense physical and chemical phenomena. Phase field models have the robustness of the level set method, combined with the ability to implement surface-specific chemistry that is required to model crystal growth, although they do not necessarily directly solve for the advancing front location. We have adopted a phase field approach and will present results of the adatom density and step-growth location in time as a function of processing conditions, such as temperature and plasma gas composition.

Hwang, Helen H.↗
NASA NTRS1991

Detection of reflector surface from near field phase measurements

The deviation of a reflector antenna surface from a perfect parabolic shape causes degradation of the performance of the antenna. The problem of determining the shape of the reflector surface in a reflector antenna using near field phase measurements is not a new one. A recent issue of the IEEE tansactions on Antennas and Propagation (June 1988) contained numerous descriptions of the use of these measurements: holographic reconstruction or inverse Fourier transform. Holographic reconstruction makes use of measurement of the far field of the reflector and then applies the Fourier transform relationship between the far field and the current distribution on the reflector surface. Inverse Fourier transformation uses the phase measurements to determine the far field pattern using the method of Kerns. After the far field pattern is established, an inverse Fourier transform is used to determine the phases in a plane between the reflector surface and the plane in which the near field measurements were taken. These calculations are time consuming since they involve a relatively large number of operations. A much faster method can be used to determine the position of the reflector. This method makes use of simple geometric optics to determine the path length of the ray from the feed to the reflector and from the reflector to the measurement point. For small physical objects and low frequencies, diffraction effects have a major effect on the error, and the algorithm provides incorrect results. It is believed that the effect is less noticeable for large distortions such as antenna warping, and more noticeable for small, localized distortions such as bumps and depressions such as might be caused by impact damage.

Ida, Nathan↗
NASA NTRS1975

Generation of ULF waves by electric or magnetic dipoles

The generation of ULF waves by ground-based magnetic and electric dipoles is studied with a simplified model consisting of three adjoining homogeneous regions representing the groud, the vacuum (free space) region, and the ionosphere. The system is assumed to be immersed in a homogeneous magnetic field with an arbitrary tilt angle. By the use of Fourier techniques and the method of stationary phase, analytic expressions are obtained for the field strength of the compressional Alfven waves in the ionosphere. Expressions are also obtained for the strength of the torsional Alfven wave in the ionosphere and the ULF magnetic field at ground level. Numerical results are obtained for the compressional Alfven-wave field strength in the ionosphere with a nonvertical geomagnetic field and for the ULF magnetic field at ground level for a vertical geomagnetic field.

Harker, K. J.↗
NASA NTRS2008

Forecasting Lightning Threat using Cloud-Resolving Model Simulations

Two new approaches are proposed and developed for making time and space dependent, quantitative short-term forecasts of lightning threat, and a blend of these approaches is devised that capitalizes on the strengths of each. The new methods are distinctive in that they are based entirely on the ice-phase hydrometeor fields generated by regional cloud-resolving numerical simulations, such as those produced by the WRF model. These methods are justified by established observational evidence linking aspects of the precipitating ice hydrometeor fields to total flash rates. The methods are straightforward and easy to implement, and offer an effective near-term alternative to the incorporation of complex and costly cloud electrification schemes into numerical models. One method is based on upward fluxes of precipitating ice hydrometeors in the mixed phase region at the-15 C level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domain-wide statistics of the peak values of simulated flash rate proxy fields against domain-wide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. Our blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Exploratory tests for selected North Alabama cases show that, because WRF can distinguish the general character of most convective events, our methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because the models tend to have more difficulty in predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models,the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of forecasts become available.

McCaul, Eugene W., Jr.↗