Search NASA⌕ Search

SEARCH · Search NASA

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

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↗

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↗

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

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

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

Olivine Dissolution and Formation of Secondary phases in Ultramafic Soils

Introduction: Olivine has been proposed as an indicator for the duration of water-rock interaction within Martian rocks and sediments [1-3]. The use of olivine as a mineralogical indicator for past aqueous alteration on Mars requires interpretation of a complex combination of factors including pH, temperature, and composition [5,6]. Here, we examine the persistence of natural olivine within terrestrial ultramafic soils (Fe/Mg-rich, Al-poor) developing under different climatic conditions and the incipient dissolution of emplaced forsterite (Fo) and fayalite (Fa) surfaces to investigate environmental effects on incipient olivine dissolution, olivine persistence in soils, and formation of secondary phases. Methods: Field Sites. We examined olivine weathering and secondary material formation in ultramafic soils at 6 sites in the Klamath Mountains (KM) of northern California with a mean annual temperature of ~12.8℃ and precipitation of ~55.7-95.4 cm/year [7], and soil pH of ~6.5-7.3; 4 sites in the Tablelands (TB) of Newfoundland, Canada with a mean annual temperature of <3.9℃ and precipitation of ~120.0 cm/year [8], and soil pH of ~7.7; and at 3 sites at Pickhandle Gulch (PG), Nevada with a mean annual temperature of ~14.4℃ and precipitation of ~14.1 cm/year [7], and soil pH of ~8.5. Sampling sites span an age range of ~12.1-50+ kya in the Klamath Mountains [9,10] and ~13-30 kya in the Tablelands [11]. Pickhandle Gulch sites are undated. Parent Material and Soil Analyses. Polished thin sections of bulk soil prepared by Wagner Petrographic, Inc were carbon-coated and analyzed on a JEOL 2100 SEM in back-scattered electron mode in the EMIL lab at UNLV and at the 13-ID-E synchrotron beamline at Argonne National Laboratory using µXRF, µXRD, and XAS. Soil and parent material samples were powdered in a Fritsch pulverisette and analyzed by XRD and soil by VNIR. Soil preparation is further described in [12]. Disk Preparation, Burial, and Collection. Fo disks were cut from a column prepared via hot-pressing and Fa disks by sintering synthetic fayalite powder, see [13] for detail. Disks were polished to a 0.25-micron level with diamond grit. Disks were buried in 3 KM soils, 4 TB soils, and 3 PG soils, collected after exactly 365 days, and washed gently with 100% reagent grade ethanol to remove potential adhered soil material. Weathered disks and soil samples were stored in a -20℃ freezer until analysis. Unaltered control disks prepared identically to the buried disks were stored at -20℃ for the duration of the experiment. Disk Analyses. One Fo and Fa disk from each climate zone was analyzed on a variable pressure Zeiss Supra 40VP SEM at Northern Arizona University. A separate Fo and Fa disk from each climate zone was analyzed by XPS using a Physical Electronics VersaProbe II at the Penn State Univ. Materials Characterization Lab after a Na-dodecyl sulfate wash and ozonation to remove carbon contamination as in [14]. VNIR measurements were conducted at Johnson Space Center using an ASD FieldSpec3 under ambient lab conditions on a separate Fo and Fa disk from each climate zone. One separate Fo and Fa control sample was analyzed for each technique for comparison with weathered samples. XPS uncertainty was determined from 5 repeat measurements on controls. Results: Bedrock and Soil Results Olivine is present in the parent material in the KM and TB. Olivine is found in ~12.1 ka KM soils but is absent from all older soils, while persisting into the oldest (>20 ka) TB soil (Figure 1). In both locations, olivine is found as cores surrounded by a serpentine rind (Figure 2). VNIR spectra from the analyzed soils possess strong OH-associated spectral features at ~2.33 µm indicating the presence of Mg-rich phyllosilicates as well as ferric-oxide features at ~0.92 µm in the KM (Figure 3). Primary crystalline silicate grains mostly incorporate Fe2+, while poorly crystalline weathering rinds are best fit by ferric oxide XAS standards (Figure 4). µXRF also shows that Fe and Ni concentrate in weathering rinds and Cr remains within interior silicate grains (Figure 4). Buried Sample Results All Fo surfaces exhibited formation of dissolution features including shallow pitting not observed on controls. Dissolution features were most visually widespread on the KM disk (Figure 5). Leaching of Mg from KM and TB Fo disks was evident from <1.6 Mg/Si ratios measured by XPS (Figure 6). Fe-rich precipitates in SEM (Figure 5) and Fe presence in XPS scans (Figure 6) indicate Fe deposition onto KM and TB Fo disk surfaces. The appearance of a spectral feature at 0.55 µm in the VNIR spectra from the TB Fo suggests this Fe is ferric (Figure 7). The PG Fo appears least altered, with minimal formation of dissolution features in SEM (Figure 5), a Mg/Si ratio inconsistent with leaching (~2) (Figure 6), and VNIR spectra almost identical to the control sample. Analysis of Fa surfaces is ongoing. The higher temperatures and more acidic pH in the KM soils likely drive the faster dissolution of the Fo disks described above. While the TB soils experience greater precipitation than in the KM, the cooler temperatures and more basic soil pH facilitate observable but more limited alteration. The dry climate and basic soil pH at PG lead to minimal dissolution of the PG disk surfaces.

A D Feldman↗

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↗

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↗

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

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

Evaluation of Methods for In-Situ Calibration of Field-Deployable Microphone Phased Arrays

Current field-deployable microphone phased arrays for aeroacoustic flight testing require the placement of hundreds of individual sensors over a large area. Depending on the duration of the test campaign, the microphones may be required to stay deployed at the testing site for weeks or even months. This presents a challenge in regards to tracking the response (i.e., sensitivity) of the individual sensors as a function of time in order to evaluate the health of the array. To address this challenge, two different methods for in-situ tracking of microphone responses are described. The first relies on the use of an aerial sound source attached as a payload on a hovering small Unmanned Aerial System (sUAS) vehicle. The second relies on the use of individually excited ground-based sound sources strategically placed throughout the array pattern. Testing of the two methods was performed in microphone array deployments conducted at Fort A.P. Hill in 2015 and at Edwards Air Force Base in 2016. The results indicate that the drift in individual sensor responses can be tracked reasonably well using both methods. Thus, in-situ response tracking methods are useful as a diagnostic tool for monitoring the health of a phased array during long duration deployments.

Humphreys, William M.↗

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

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↗

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

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

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

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↗