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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Aerothermal modeling program. Phase 2, element A: Improved numerical methods for turbulent viscous recirculating flows

The objective of this effort is to develop improved numerical schemes for predicting combustor flow fields. Various candidate numerical schemes were evaluated, and promising schemes were selected for detailed assessment. The criteria for evaluation included accuracy, computational efficiency, stability, and ease of extension to multidimensions. The candidate schemes were assessed against a variety of simple one- and two-dimensional problems. These results led to the selection of the following schemes for further evaluation: flux spline schemes (linear and cubic) and controlled numerical diffusion with internal feedback (CONDIF). The incorporation of the flux spline scheme and direct solution strategy in a computer program for three-dimensional flows is in progress.

Karki, K. C.↗

NASA NTRS2002

High-Fidelity Generalization Method of Cells for Inelastic Periodic Multiphase Materials

An extension of a recently-developed linear thermoelastic theory for multiphase periodic materials is presented which admits inelastic behavior of the constituent phases. The extended theory is capable of accurately estimating both the effective inelastic response of a periodic multiphase composite and the local stress and strain fields in the individual phases. The model is presently limited to materials characterized by constituent phases that are continuous in one direction, but arbitrarily distributed within the repeating unit cell which characterizes the material's periodic microstructure. The model's analytical framework is based on the homogenization technique for periodic media, but the method of solution for the local displacement and stress fields borrows concepts previously employed by the authors in constructing the higher-order theory for functionally graded materials, in contrast with the standard finite-element solution method typically used in conjunction with the homogenization technique. The present approach produces a closed-form macroscopic constitutive equation for a periodic multiphase material valid for both uniaxial and multiaxial loading. The model's predictive accuracy in generating both the effective inelastic stress-strain response and the local stress said inelastic strain fields is demonstrated by comparison with the results of an analytical inelastic solution for the axisymmetric and axial shear response of a unidirectional composite based on the concentric cylinder model, and with finite-element results for transverse loading.

Aboudi, Jacob↗

DOE OSTI2020

Covariant phase space with boundaries

The covariant phase space method of Iyer, Lee, Wald, and Zoupas gives an elegant way to understand the Hamiltonian dynamics of Lagrangian field theories without breaking covariance. The original literature however does not systematically treat total derivatives and boundary terms, which has led to some confusion about how exactly to apply the formalism in the presence of boundaries. In particular the original construction of the canonical Hamiltonian relies on the assumed existence of a certain boundary quantity “B”, whose physical interpretation has not been clear. We here give an algorithmic procedure for applying the covariant phase space formalism to field theories with spatial boundaries, from which the term in the Hamiltonian involving B emerges naturally. Our procedure also produces an additional boundary term, which was not present in the original literature and which so far has only appeared implicitly in specific examples, and which is already nonvanishing even in general relativity with sufficiently permissive boundary conditions. The only requirement we impose is that at solutions of the equations of motion the action is stationary modulo future/past boundary terms under arbitrary variations obeying the spatial boundary conditions; from this the symplectic structure and the Hamiltonian for any diffeomorphism that preserves the theory are unambiguously constructed. We show in examples that the Hamiltonian so constructed agrees with previous results. We also show that the Poisson bracket on covariant phase space directly coincides with the Peierls bracket, without any need for non-covariant intermediate steps, and we discuss possible implications for the entropy of dynamical black hole horizons.

36 MATERIALS SCIENCE↗

DOE OSTI2022

A review of the metastable omega phase in beta titanium alloys: the phase transformation mechanisms and its effect on mechanical properties

Since its discovery in 1954, the omega (ω) phase in titanium and its alloys has attracted substantial attention from researchers. The β-to-ω and ω-to-α phase transformations are central to β-titanium alloy design, but the transformation mechanisms have been a subject of debate. With new generations of aberration-corrected transmission electron microscopy and atom probe tomography, both the spatial resolution and compositional sensitivity of phase transformation analysis have been rapidly improving. Furthermore, this review provides a detailed assessment of the new understanding gained and related debates in this field enabled by advanced characterization methods. Specifically, new insights into the possibility of a coupled diffusional-displacive component in the β-to-ω transformation and key nucleation driving forces for the ω-assisted α phase formation are discussed. Additionally, the influence of ω phase on the mechanical properties of β-titanium alloys is also reviewed. Lastly, a perspective on open questions and future direction for research is discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Influence of Wetting on Viscous Fingering Via 2D Lattice Boltzmann Simulations

We present simulations of two-phase flow using the Rothman and Keller colour gradient Lattice Boltzmann method to study viscous fingering when a “red fluid” invades a porous model initially filled with a “blue” fluid with different viscosity. We conducted eleven suites of 81 numerical experiments totalling 891 simulations, where each suite had a different random realization of the porous model and spanned viscosity ratios in the range $$M\in [0.01,100]$$ M ϵ [ 0.01 , 100 ] and wetting angles in the range $$\theta _w\in [180^\circ ,0^\circ ]$$ θ w ϵ [ 180 ° , 0 ° ] to allow us to study the effect of these parameters on the fluid-displacement morphology and saturation at breakthrough (sweep). Although sweep often increased with wettability, this was not always so and the sweep phase space landscape, defined as the difference in saturation at a given wetting angle relative to saturation for the non-wetting case, had hills, ridges and valleys. At low viscosity ratios, flow at breakthrough is localized through narrow fingers that span the model. After breakthrough, the flow field continues to evolve and the saturation continues to increase albeit at a reduced rate, and eventually exceeds 90% for both non-wetting and wetting cases. The existence of a complicated sweep phase space at breakthrough, and continued post-breakthrough evolution suggests the hydrodynamics and sweep is a complicated function of wetting angle, viscosity ratio and time, which has major potential implications to Enhanced Oil Recovery by water flooding, and hence, on estimates of global oil reserves. Validation of these results via experiments is required to ensure they translate to field studies.

Engineering↗

DOE OSTI2022

Deformations, relaxation, and broken symmetries in liquids, solids, and glasses: A unified topological field theory

Here, we combine hydrodynamic and field theoretic methods to develop a general theory of phonons as Goldstone bosons in crystals, glasses, and liquids based on nonaffine displacements and the consequent Goldstone phase relaxation. We relate the conservation, or lack thereof, of specific higher-form currents with properties of the underlying deformation field—nonaffinity—which dictates how molecules move under an applied stress or deformation. In particular, the single-valuedness of the deformation field is associated with conservation of higher-form charges that count the number of topological defects. Our formalism predicts, from first principles, the presence of propagating shear waves above a critical wave vector in liquids, thus giving a formal derivation of the phenomenon in terms of fundamental symmetries. The same picture provides also a theoretical explanation of the corresponding “positive sound dispersion” phenomenon for longitudinal sound. Importantly, accordingly to our theory, the main collective relaxation timescale of a liquid or a glass (known as the α relaxation for the latter) is given by the phase relaxation time, which is not necessarily related to the Maxwell time. Finally, we build a nonequilibrium effective action using the in-in formalism defined on the Schwinger-Keldysh contour, that further supports the emerging picture. In summary, our work suggests that the fundamental difference between solids, fluids, and glasses has to be identified with the associated generalized higher-form global symmetries and their topological structure, and that the Burgers vector for the displacement fields serves as a suitable topological order parameter distinguishing the solid (ordered) phase and the amorphous ones (fluids, glasses).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Recent advances in describing and driving crystal nucleation using machine learning and artificial intelligence

With the advent of faster computer processors and especially graphics processing units (GPUs) over the last few decades, the use of data-intensive machine learning (ML) and artificial intelligence (AI) has increased greatly, and the study of crystal nucleation has been one of the beneficiaries. In this study, we outline how ML and AI have been applied to address four outstanding difficulties of crystal nucleation: how to discover better reaction coordinates (RCs) for describing accurately non-classical nucleation situations; the development of more accurate force fields for describing the nucleation of multiple polymorphs or phases for a single system; more robust identification methods for determining crystal phases and structures; and as a method to yield improved course-grained models for studying nucleation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Sub-tesla On-Chip Nanomagnetic Metamaterial Platform for Angle-Resolved Photoemission Spectroscopy

Magnetically controlled states in quantum materials are central to their unique electronic and magnetic properties. However, direct momentum-resolved visualization of these states via angle-resolved photoemission spectroscopy (ARPES) has been hindered by the disruptive effect of magnetic fields on photoelectron trajectories. Here, we introduce an in situ method that is, in principle, capable of applying magnetic fields up to 1 T. This method uses substrates composed of nanomagnetic metamaterial arrays with alternating polarity. Such substrates can generate strong, homogeneous, and spatially confined fields applicable to samples with thicknesses up to the micron scale, enabling ARPES measurements under magnetic fields with minimal photoelectron trajectory distortion. Here, we demonstrate this minimal distortion with ARPES data taken on monolayer graphene. Our method paves the way for probing magnetic field-dependent electronic structures and studying field-tunable quantum phases with state-of-the-art energy-momentum resolutions.

Magnetic properties↗

DOE OSTI2023

Statistical Generation of Ocean Forcing With Spatiotemporal Variability for Ice Sheet Models

Melting of ice at the base of floating ice shelves that fringe the Antarctic ice sheet has been identified as a significant source of uncertainty in sea level rise projections. Part of this uncertainty derives from chaotic internal variability of the coupled ocean-atmosphere system. For numerical ice sheet model projections, this uncertainty has not previously been quantified because of the prohibitive computational expense of running large climate model ensembles. Here, we develop and demonstrate a technique that generates independent realizations of internal climate variability from a single climate model simulation. Building on previous developments in model emulation, this technique uses empirical orthogonal function decomposition and Fourier-phase randomization to generate statistically consistent realizations of spatiotemporal variability fields for the target climate variable. The method facilitates efficient sampling of a wide range of climate trajectories, which can also be incorporated within ice sheet or other physical models to represent feedback processes.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1984

Indirect Phase Height Measurements in Central and Eastern Europe for Monitoring D Region Plasma

Low-frequency propagation experiments for the investigation of the lower part of the ionospheric D region were at first used by BRACEWELL et al. (1951) in the early fifties. Among these was the method of indirect phase height measurements for continuous monitoring of the lower ionosphere. It is based upon field strength measurements of commercial radio transmitters in the frequency range between 50 and 200 kHz at distances from 500 to 1500 km. The field strength records show characteristic diurnal variations with maxima and minima, produced by interference between the ground wave and the ionospherically reflected sky wave, the phase difference between varies in correspondence to the diurnal variation of the reflection height. In order to check the validity of interpretations of indirect phase height data, comparisons were made with simultaneous rocket soundings. Results are summarized.

Cossart, G. V.↗

NASA NTRS2001

Restraint of Liquid Jets by Surface Tension in Microgravity Modeled

Tension in Microgravity Modeled Microgravity poses many challenges to the designer of spacecraft tanks. Chief among these are the lack of phase separation and the need to supply vapor-free liquid or liquidfree vapor to the spacecraft processes that require fluid. One of the principal problems of phase separation is the creation of liquid jets. A jet can be created by liquid filling, settling of the fluid to one end of the tank, or even closing a valve to stop the liquid flow. Anyone who has seen a fountain knows that jets occur in normal gravity also. However, in normal gravity, the gravity controls and restricts the jet flow. In microgravity, with gravity largely absent, jets must be contained by surface tension forces. Recent NASA experiments in microgravity (Tank Pressure Control Experiment, TPCE, and Vented Tank Pressure Experiment, VTRE) resulted in a wealth of data about jet behavior in microgravity. VTRE was surprising in that, although it contained a complex geometry of baffles and vanes, the limit on liquid inflow was the emergence of a liquid jet from the top of the vane structure. Clearly understanding the restraint of liquid jets by surface tension is key to managing fluids in low gravity. To model this phenomenon, we need a numerical method that can track the fluid motion and the surface tension forces. The fluid motion is modeled with the Navier-Stokes equation formulated for low-speed incompressible flows. The quantities of velocity and pressure are placed on a staggered grid, with velocity being tracked at cell faces and pressure at cell centers. The free surface is tracked via the introduction of a color function that tracks liquid as 1/2 and gas as -1/2. A phase model developed by Jacqmin is used. This model converts the discrete surface tension force into a barrier function that peaks at the free surface and decays rapidly. Previous attempts at this formulation have been criticized for smearing the interface. However, by sharpening the phase function, double gridding the fluid function, and using a higher order solution for the fluid function, interface smearing is avoided. These equations can be rewritten as two coupled Poisson equations that also include the velocity. The method of solution is as follows: first, the phase equations are solved from this solution, a velocity field is generated, then a successive overrelaxation scheme is used to solve for a pressure field consistent with the velocity solution. After the code was implemented in axisymmetric form and verified by several test cases, the drop tower runs of Aydelott were modeled. The model handed the free-surface deformation quite nicely, even to the point of modeling geyser growth in the regime where the free surface was no longer restrained. A representative run is shown.

Chato, David J.↗

DOE OSTI2022

Modulating the Pseudoelastic Response of Nitinol Using Ion Implantation

This work explores whether ion beam modification can be used to modulate the austenite to martensite phase transformation in Nickel-Titanium (NiTi), thereby achieving novel or localized transformation properties in near-surface regions. We report this could provide alternatives to laser shot peening or other surface treatment methods and possibly expand applications in biomedical, aerospace, and other fields. Irradiation induces defects and internal stress that can serve as nucleation and/or pinning sites for the phase transformation. Thus, it can augment more convention- al approaches, including alloying, severe mechanical work, grain size reduction, and precipitation of coherent precipitates. A range of outcomes is possible in principle, including a shift of the critical stress or temperature for onset of the transformation, linearization, reduction of hysteresis, stabilization, and extent of transformation strain.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Alignment Method for Synchronized Phase Angle Measurement With Presence of Practical Time Shift

Synchronized Phase Angle Measurements (SPAMs) are widely used in power systems in the applications of situational awareness. However, the practical time shifts can lead to unexpected angle differences among Phasor Measurement Units (PMUs) manufactured from various vendors. Moreover, since most of PMUs calculate phase angle via Discrete Fourier Transform based approaches, this issue becomes even worse under the off-nominal frequency condition. To mitigate the impact of practical timing shift, this letter presents a fast method to rectify the SPAM for appropriate alignment. To verify the performance of the proposed method, laboratory and field tests have been conducted by implementing the method in PMUs and phasor data concentrators, respectively. The results demonstrate that the standard deviations of the phase angle differences have significantly decreased to 0.1$^{\circ }$ order with the adoption of the proposed method.

42 ENGINEERING↗

NASA NTRS2000

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Tryggvason, G.↗

DOE OSTI2024

Thermodynamic stability of a spin microemulsion in Rashba spin-orbit-coupled bosons

Recent finite-temperature numerical simulations have unveiled a quantum “spin” microemulsion analog, found by raising the temperature of a stripe supersolid phase in a Rashba spin-orbit coupled Bose gas. This microemulsion state is a highly correlated, isotropic normal fluid where atoms self-arrange based on their internal pseudospin into patterns that resemble bicontinuous microemulsions. This finding leaves several open questions regarding the broader accessibility of this phase in experiments. Here, we use equilibrium finite-temperature numerical simulations based on a coherent-state path integral representation to perform a computational investigation into the thermodynamic stability of the spin microemulsion state. Numerical simulations emphasize the requirement of a nearly, but not perfectly, isotropic spin-orbit coupling in order to achieve the microemulsion phase in cold-atom experiments. Moreover, the microemulsion state exists independent of miscibility of the pseudospin components and for a wide range of pseudospin population imbalance, suggesting a high degree of flexibility in choosing the atom and hyperfine states in an experimental realization. Lastly, we demonstrate this feasibility by mimicking a Rashba spin-orbit-coupled 87 Rb experiment in an isotropic harmonic trap, where we confirm the microemulsion's existence via its density profile and equilibrium quasimomentum distribution.

Complex Langevin dynamics↗

DOE OSTI2024

Segregation of Al and its effect on coercivity in Nd-Fe-B

Doping a small amount of Al can effectively enhance coercivity in Nd-Fe-B magnets. Here we investigated the partitioning behavior of Al and its effect on coercivity in Nd-Fe-B using first principles DFT (density functional theory) calculation and micromagnetic simulation. The calculated substitution energies of Fe by Al are negative at the crystallographic sites of 4c and 8j 2 while they are positive values at the other sites in Nd 2 Fe 14 B (2:14:1), implying a small solubility of Al in 2:14:1. Further, Al prefers to segregate at grain boundary (GB) and stabilize the Nd-Fe-Al phase with a Nd 6 Fe 13 Si-type tetragonal structure (6:13:1). The formation of the antiferromagnetic or weak ferrimagnetic 6:13:1-like phase depletes Fe and reduces the amount of ferromagnetic Nd–Fe type grain boundary phase (GBP), which weakens the inter-grain magnetic interaction. Micromagnetic simulations indicate that the 6:13:1-like GBP increases the pinning field of magnetic domain wall at GB and suppresses the nucleation of reversal magnetic domain on the grain surface of 2:14:1 during demagnetization process. The formation of Al-rich shell on 2:14:1 grain surface can further moderately increase the domain pinning field at GB and the nucleation field of 2:14:1 grain. Developing novel processing method to tailor Al segregation and promote formation of 6:13:1-like phase at GB can be a promising approach to improve coercivity in Nd-Fe-B magnet.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Coded-mask-based X-ray phase-contrast and dark-field imaging

Phase contrast and dark-field X-ray imaging enable imaging of objects that absorb or reflect very little X-ray light. Disclosed is a method and systems for performing coded-mask-based multi-contrast imaging (CMMI). The method includes providing radiation to a coded mask that has a known phase and absorption profile according to a pre-determined pattern. The radiation is then impingent upon a sample, and the radiation is detected to perform phase-reconstruction and image processing. The method and associated systems allow for the use of maximum-likelihood and machine learning methods for reconstruction images of the sample from the detected radiation.

Shi, Xianbo↗

DOE OSTI2021

Understanding Reentrance in Frustrated Magnets: The Case of the Er 2 Sn 2 O 7 Pyrochlore

Reentrance, the return of a system from an ordered phase to a previously encountered less-ordered one as a controlled parameter is continuously varied, is a recurring theme found in disparate physical systems, yet its microscopic cause is often not investigated thoroughly. Here, through detailed characterization and theoretical modeling, we uncover the microscopic mechanism behind reentrance in the strongly frustrated pyrochlore antiferromagnet Er 2 Sn 2 O 7 . We use single crystal heat capacity measurements to expose that Er 2 Sn 2 O 7 exhibits multiple instances of reentrance in its magnetic field B vs temperature T phase diagram for magnetic fields along three cubic high symmetry directions. Through classical Monte Carlo simulations, mean field theory, and classical linear spin-wave expansions, we argue that the origins of the multiple occurrences of reentrance observed in Er 2 Sn 2 O 7 are linked to soft modes. These soft modes arise from phase competition and enhance thermal fluctuations that entropically stabilize a specific ordered phase, resulting in an increased transition temperature for certain field values and thus the reentrant behavior. Our work represents a detailed examination into the mechanisms responsible for reentrance in a frustrated magnet and may serve as a template for the interpretation of reentrant phenomena in other physical systems.

36 MATERIALS SCIENCE↗