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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Optical Characterization of a Multipoint Lean Direct Injector for Gas Turbine Combustors: Velocity and Fuel Drop Size Measurements

Performance of a multipoint, lean direct injection (MP-LDI) strategy for low emission aero-propulsion systems has been tested in a Jet-A fueled, lean flame tube combustion rig. Operating conditions for the series of tests included inlet air temperatures between 672 and 828 K, pressures between 1034 and 1379 kPa and total equivalence ratios between 0.41 and 0.45, resulting in equilibrium flame temperatures approaching 1800 K. Ranges of operation were selected to represent the spectrum of subsonic and supersonic flight conditions projected for the next-generation of commercial aircraft. This document reports laser-based measurements of in situ fuel velocities and fuel drop sizes for the NASA 9-point LDI hardware arranged in a 3 3 square grid configuration. Data obtained represent a region of the flame tube combustor with optical access that extends 38.1-mm downstream of the fuel injection site. All data were obtained within reacting flows, without particle seeding. Two diagnostic methods were employed to evaluate the resulting flow path. Three-component velocity fields have been captured using phase Doppler interferometry (PDI), and two-component velocity distributions using planar particle image velocimetry (PIV). Data from these techniques have also offered insight into fuel drop size and distribution, fuel injector spray angle and pattern, turbulence intensity, degree of vaporization and extent of reaction. This research serves to characterize operation of the baseline NASA 9- point LDI strategy for potential use in future gas-turbine combustor applications. An additional motive is the compilation of a comprehensive database to facilitate understanding of combustor fuel injector aerodynamics and fuel vaporization processes, which in turn may be used to validate computational fluid dynamics codes, such as the National Combustor Code (NCC), among others.

Heath, Christopher M.↗
NASA NTRS1973

A conformal head-up display for the visual approach

The degree of conformity used in matching a superimposed display to its visual background is considered in relation to the information available for vertical guidance and control during a purely visual approach. The information may be represented by individual symbols or combined in a single symbol, and the relative merits of these methods are discussed. A fully conformal display format is developed for the purpose of showing both the position and direction of the flight path, with provision for the effects of disturbances, ILS compatibility, and control needs. The field of view needed for all conditions and phases of the visual approach with a fully conformal display is studied in relation to the limitations of conventional collimator systems. Methods are discussed which depend on deviation of the sight line, and on windshield reflection of the uncollimated image of a simple pointer. Limited flight tests show some promise for the uncollimated method.

Naish, J. M.↗
NASA NTRS2011

Analysis of the Radiative Transfer Equation with Highly Asymmetric Phase Function

This paper considers a scalar radiative transfer problem with high scattering anisotropy, Two computational methods are presented based on decomposition of the diffuse light field into a regular and anisotropic part. The first algorithm (DOMAS) singles out the anisotropic radiance in the forward scattering peak using the Small-Angle Modification of RTE. The second algorithm (DOM2+) separates the single scattering radiance as an anisotropic part, which largely defines the fine detail of the total radiance in the backscattering directions. In both cases, the anisotropic part is represented analytically. With anisotropy subtraction, the regular part of the signal. which requires a numerical solution, is essentially smoothed as a function of angles. Further, the transport equation is obtained for the regular part that contains an additional source function from the anisotropic part of the signal. This equation is solved with the discrete ordinates method. A conducted numerical analysis of this work showed that algorithm DOMAS has a strong advantage as compared to the standard discrete ordinates method for simulation of the radiance transmission, and DOM2 + is the best of the three for the reflection computations. Both algorithms offer at least a factor of three acceleration of convergence of the azimuthal series for highly anisotropic phase functions.

Korkin, Sergey V.↗
DOE OSTI2025

Quantum Mechanics of a Spherically Symmetric Causal Diamond in Minkowski Spacetime

We construct the phase space of a spherically symmetric causal diamond in ( d + 2 )-dimensional Minkowski spacetime. Utilizing the covariant phase space formalism, we identify the relevant degrees of freedom that localize to the d -dimensional bifurcate horizon and, upon canonical quantization, determine their commutators. On this phase space, we find two Iyer-Wald charges. The first of these charges, proportional to the area of the causal diamond, is responsible for shifting the null time along the horizon and has been well documented in the literature. The second charge is much less understood, being integrable for d ≥ 2 only if we allow for field-dependent diffeomorphisms and is responsible for changing the size of the causal diamond. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Field To Farm Aggregation For Agricultural Systems

The Fields to Farms methodology illustrates the generation of farm parcels from the Crop Data Layer (CDL), a raster dataset containing 133 categories representing various crop types and land uses. This methodology involves two primary steps: Field Delineation and Farm Aggregation. The code specifically addresses the aggregation of pre-delineated fields within a county to form farms, adhering to predefined criteria for farm size categories. It is assumed that the field delineation process, which involves creating vector polygons from CDL raster, has been completed beforehand, possibly through external tools or methods. Upon initialization, the script processes county-level fields, preparing them for farm aggregation. In the Farm Aggregation phase, the code iteratively combines delineated fields into farms based on specified criteria, continuing until the aggregated farm size meets predefined thresholds derived from data from the 2017 National Agricultural Statistics Service (NASS) census. Throughout this iterative process, the script dynamically adjusts the aggregation to ensure alignment with the desired distribution reported by NASS. The resulting output of the script is a GeoDataFrame containing classified farms, which are subsequently saved as GeoPackage files. These files enable further analysis and visualization, facilitating comprehensive exploration of the farm landscape generated through the methodology.

Paudel, Rajiv [Idaho National Laboratory (INL), Id↗
DOE OSTI2026

PyJMAK: An Open-Source Python Toolkit for Modeling Solid-State Metallurgical Phase Transformations

Accurate prediction of metallurgical phase transformations is an essential basis for autonomous optimization and rapid part qualification. Several methods can be used to estimate the evolution of phase fractions such as JMAK kinetics-based models, phase-field models, thermodynamic models, and data-driven machine learning models. Thermodynamic and phase-field-based methodologies solve multiphysics equations requiring numerous calibration parameters and significant computational resources. As a result, the computation domain is limited to a point or on order of micron-meters. The data-driven models rely on large datasets from experiments and simulations. While the JMAK model only provides information about phase fraction evolution, it can predict this evolution in near real-time using thermal history and thermodynamic data without restriction on the domain. JMAK models have been popularly used by researchers to model phase transformations occuring during additive manufacturing or over arbitrary temperature profiles. Commercial proprietary software such as Abaqus and Ansys or closed-source in-house implementations offer the ability to model JMAK based kinetics to predict phase transformation. However, these software packages are not open-source or freely available for use and development in conjunction with manufacturing machines, sensors, and machine learning algorithms. In addition, the use of the model is restricted by a license token. In contrast, given temperature profiles at multiple points in the domain, this Python-based PyJMAK model can compute phase evolution in parallel due to its stand-alone modular, voxel-based structure, and it can be executed on high-performance computing resources without any license restrictions.

Prabhune, Bhagya [Oak Ridge National Laboratory (O↗
NASA NTRS1993

Multilayer thin film design as far ultraviolet quarterwave retarders

At short wavelengths, such as FUV, transparent, optically active materials are scarce. Reflection phase retardation by a multilayer thin film can be a good alternative in this wavelength region. We design a multilayer quarterwave retarder by calculating the electric fields at each boundary in the multilayer thin film. Using this method, we achieve designs of FUV multilayers which provide high, matched reflectances for both s- and p-polarization states, and at the same time a phase difference between these two states of nearly 90 deg. For example, a quarterwave retarder designed at the Lyman-alpha line (121.6 nm) has 81.05 percent reflectance for the s-polarization and 81.04 percent for the p-polarization state. The phase difference between these two polarization states is 90.07 deg. For convenience the retarders are designed for 45 deg angle of incidence, but our design approach can be used for any other angle of incidence. Aluminum and MgF2 are used as film materials and an opaque thick film of aluminum as the substrate.

Kim, Jongmin↗
NASA NTRS1999

Effects of G-Jitter on Directional Solidification of a Binary Alloy

A study of directional solidification of a weak binary alloy (specifically, Bi - 1 at% Sn) based on the fixed grid single domain approach is being undertaken. The enthalpy method is used to solve for the temperature field over the computational domain including both the solid and liquid phases; latent heat evolution is treated with the aid of an effective specific heat coefficient. A source term accounting for the release of solute into the liquid during solidification has been incorporated into the solute transport equation. The vorticity-stream function formulation is used to describe thermo-solutal convection in the liquid region. In this paper we present a numerical simulation of g-jitter. A background gravity of 1 microgram has been assumed, and new results for the effects of periodic disturbances over a range of amplitudes and frequencies on solute field and segregation have been presented.

Benjapiyaporn, C.↗
NASA NTRS2000

Effects of Space Environment on Flow and Concentration During Directional Solidification

A study of directional solidification of a weak binary alloy (specifically, Bi - 1 at% Sn) based on the fixed grid single domain approach is being undertaken. The enthalpy method is used to solve for the temperature field over the computational domain including both the solid and liquid phases; latent heat evolution is treated with the aid of an effective specific heat coefficient. A source term accounting for the release of solute into the liquid during solidification has been incorporated into the solute transport equation. The vorticity-stream function formulation is used to describe thermosolutal convection in the liquid region. In this paper we numerically investigate the effects of g-jitter on directional solidification. A background gravity of 1 micro-g has been assumed, and new results for the effects of periodic disturbances over a range of amplitudes and frequencies on solute field and segregation have been presented.

Benjapiyaporn, C.↗
NASA NTRS2000

Effect of g-jitter on Directional Solidification of a Binary Alloy

A study of directional solidification of a weak binary alloy (specifically, Bi - 1 at% Sn) based on the fixed grid single domain approach is being undertaken. The enthalpy method is used to solve for the temperature field over the computational domain including both the solid and liquid phases; latent heat evolution is treated with the aid of an effective specific heat coefficient. A source term accounting for the release of solute into the liquid during solidification has been incorporated into the solute transport equation. The vorticity-stream function formulation is used to describe thermosolutal convection in the liquid region. In this paper we present a numerical simulation of g-jitter: the small, rapid fluctuations in gravitational acceleration which may be experienced in an orbiting space vehicle. A background gravity of 1 micro-g has been assumed, and new results for the effects of orientation angle of the periodic disturbances over a range of amplitudes and frequencies on solute field and segregation have been presented.

Santiviriyapanich, P.↗
NASA NTRS2003

Turbulent Mixing of Multiphase Flow

Thus we conduct numerical simulations of multiphase fluids stirred by two-dimensional turbulence to assess the possibility of self-similar drop size distribution in turbulence. In our turbulence simulations, we also explore the non-diffusive limit, where molecular mobility for the interface is vanishing. Special care is needed to transport the non-diffusive interface. Numerically, we use the particle level set method to evolve the interface. Instead of using the usual methods to calculate the surface tension force from the level set function, we reconstruct the interface based on phase- field modeling, and calculate the continuum surface tension forcing from the reconstructed interface.

Young, Y.-N.↗
NASA NTRS2021

Evidence from Galactic Cosmic Rays That the Sun Has Likely Entered A Secular Minimum in Solar Activity

Since the beginning of the space age, the Sun has been in a multi-cycle period of elevated activity (secular maximum). This secular maximum is the longest in the last 9300 years. Since the end of solar cycle 21 (SC21), however, the Sun has shown a decline in over-all activity, which has remarkably increased the fluxes of galactic cosmic rays (GCRs). Here, we investigate the correlation between the modulation of GCRs, the heliospheric magnetic field, and solar wind speed for the last 24 solar cycles to find trends that can potentially be used to predict future solar activity. Specifically, we develop a tool for predicting future magnetic field intensity, based on the hysteresis in the GCR variation, during the last phases of the current cycle. This method estimates that SC25 will be as weak as or weaker than SC24. This would mean that the Sun has likely entered a secular minimum, which, according to historical records, should last for another two cycles (SC2527and SC26).

Solar cycles↗
DOE OSTI2025

Searches for New Physics With Muon Conversion at Fermilab and Triboson Production at the LHC

We report on several efforts to search for physics beyond the standard model of particle physics at broad energy scales. The Mu2e experiment at Fermilab will search for charged lepton flavor violation via the muon to electron conversion process, which is suppressed in the Standard Model. Mu2e will be operated at a low energy, yet can probe New Physics at very high mass scales (O(1e3 - 1e4 ) TeV). At high energies, the CMS experiment at the CERN LHC continues to deliver an impressive suite of Standard Model measurements and limits on a variety of New Physics signatures. Mu2e is under construction and slated to collect its first physics data in the coming years. This thesis describes work done during the construction phase of Mu2e and focuses on two critical areas: magnetic field modeling and statistical analysis. We describe a novel method for field modeling which we validate using a simulated dataset representing the expected magnetic field in the Detector Solenoid. This method blends a standard least-squares fitting technique that utilizes physically motivated analytical model functions with a novel physics informed network that is constructed to obey Maxwell’s equations. We show the technique can model the field with an accuracy of 10−7 despite the presence of injected noise in the pseudo-measurements at the 10−5 level. We then present preliminary results of the calibration of 3D Hall probes at the sub-10−4 level. These probes will be used to directly measure the Mu2e Detector Solenoid magnetic field on a sparse grid; these measurements serve as the input to the field model fitting. Finally, we describe the first implementation of both an unbinned shape analysis and a Bayesian interpretation applied to Mu2e pseudo-data. Up to 20% tighter limits can be set by the shape analysis compared to a standard cut & count analysis. The AlCap experiment collected data at PSI in 2015 to measure several important quantities related to nuclear muon capture on an aluminum target, which is a significant background process for Mu2e. The neutron emission from muon capture can introduce background hits in the Mu2e detectors and can increase radiation damage in various elements of the apparatus. We present measurements of the neutron group fluence and mean neutron multiplicity for muon capture on aluminum nuclei. Finally, we discuss an analysis of triboson production at CMS using an Effective Field Theory framework. Standard Model triboson production, which was first observed at CMS in 2020, has a relatively small cross section and provides direct access to both anomalous triple gauge couplings and quartic gauge couplings. These couplings, interpreted in the Standard Model Effective Field Theory, are studied in the present work. We target the boosted regime where the background rate is low and yields are enhanced when dimension-6 and dimension-8 Wilson coefficients are non-zero. We do not observe an excess in the data and therefore set bounds on the Wilson coefficients. For dimension-6 coefficients the tightest observed (expected) bounds are set on cW /Λ2 where Λ is the mass scale of new physics; the bounds are [−0.13, 0.12] TeV−2 ([−0.12, 0.12] TeV−2 ) at 95% CL. The tightest bounds in dimension-8 are set on fT,0 / Λ4 ; the observed (expected) bounds at 95% CL are [−0.63, 0.69] TeV−4 ([−0.54, 0.62] TeV−4 ). Additional results are presented which include scenarios where multiple Wilson coefficients are non-zero, the application of signal model clipping to address unitarity violation in Effective Field Theories, and a novel template fit developed for easier reinterpretation of our results.

Kampa, Cole Erik [Northwestern U. (main)] (ORCID:0↗
NASA NTRS2007

Calibration of High Frequency MEMS Microphones

Understanding and controlling aircraft noise is one of the major research topics of the NASA Fundamental Aeronautics Program. One of the measurement technologies used to acquire noise data is the microphone directional array (DA). Traditional direction array hardware, consisting of commercially available condenser microphones and preamplifiers can be too expensive and their installation in hard-walled wind tunnel test sections too complicated. An emerging micro-machining technology coupled with the latest cutting edge technologies for smaller and faster systems have opened the way for development of MEMS microphones. The MEMS microphone devices are available in the market but suffer from certain important shortcomings. Based on early experiments with array prototypes, it has been found that both the bandwidth and the sound pressure level dynamic range of the microphones should be increased significantly to improve the performance and flexibility of the overall array. Thus, in collaboration with an outside MEMS design vendor, NASA Langley modified commercially available MEMS microphone as shown in Figure 1 to meet the new requirements. Coupled with the design of the enhanced MEMS microphones was the development of a new calibration method for simultaneously obtaining the sensitivity and phase response of the devices over their entire broadband frequency range. Over the years, several methods have been used for microphone calibration. Some of the common methods of microphone calibration are Coupler (Reciprocity, Substitution, and Simultaneous), Pistonphone, Electrostatic actuator, and Free-field calibration (Reciprocity, Substitution, and Simultaneous). Traditionally, electrostatic actuators (EA) have been used to characterize air-condenser microphones for wideband frequency ranges; however, MEMS microphones are not adaptable to the EA method due to their construction and very small diaphragm size. Hence a substitution-based, free-field method was developed to calibrate these microphones at frequencies up to 80 kHz. The technique relied on the use of a random, ultrasonic broadband centrifugal sound source located in a small anechoic chamber. Phase calibrations of the MEMS microphones were derived from cross spectral phase comparisons between the reference and test substitution microphones and an adjacent and invariant grazing-incidence 1/8-inch standard microphone.

Shams, Qamar A.↗
NASA NTRS1987

Electronic heterodyne moire deflectometry: A method for transient and three dimensional density fields measurements

Effects of diffraction and nonlinear photographic emulsion characteristics on the performance of deferred electronic heterodyne moire deflectometry are investigated. The deferred deflectometry is used for measurements of nonsteady phase objects where it is difficult to complete the analysis of the field in real time. The sensitivity, accuracy and resolution of the system are calculated and it is shown that they are weakly affected by diffraction and by nonlinear recording. The feactures of the system are significantly improved compared with the conventional deferred intensity moire technique, and are comparable with the online heterodyne moire. The system was evaluated experimentally by deferred measurements of the refractive index gradients of a weak phase object consisting of a large KD*P crystal. This was done by photographing the phase object through a Ronchi grating and analyzing the tranparency with the electronic heterodyne readout system. The results are compared with the measurements performed on the same phase object with online heterodyne moire deflectometry and with heterodyne holographic interferometry methods. Some practical considerations for system improvement are discussed.

Stricker, Josef↗
NASA NTRS1984

Natural convection in steady solidification - Finite element analysis of a two-phase Rayleigh-Benard problem

Galerkin finite-element approximations and Newton's method for solving free boundary problems are combined with computer-implemented techniques from nonlinear perturbation analysis to study solidification problems with natural convection in the melt. The Newton method gives rapid convergence to steady state velocity, temperature and pressure fields and melt-solid interface shapes, and forms the basis for algebraic methods for detecting multiple steady flows and assessing their stability. The power of this combination is demonstrated for a two-phase Rayleigh-Benard problem composed of melt and solid in a veritical cylinder with the thermal boundary conditions arranged so that a static melt with a flat melt-solid interface is always a solution. Multiple cellular flows bifurcating from the static state are detected and followed as Rayleigh number is varied. Changing the boundary conditions to approach those appropriate for the vertical Bridgman solidification system causes imperfections that eliminate the static state. The flow structure in the Bridgman system is related to those for the Rayleigh-Benard system by a continuous evolution of the boundary conditions.

Chang, C. J.↗
NASA NTRS1991

Neutral line chaos and phase space structure

Phase space structure and chaos near a neutral line are studied with numerical surface-of-section (SOS) techniques and analytic methods. Results are presented for a linear neutral line model with zero crosstail electric field. It was found that particle motion can be divided into three regimes dependening on the value of the conserved canonical momentum, Py, and the conserved Hamiltonian, h. The phase space structure, using Poincare SOS plots, is highly sensitive to bn = Bn/B0 variations, but not to h variations. It is verified that the slow motion preserves the action, Jz, as evaluated by Sonnerup (1971), when the period of the fast motion is smaller than the time scale of the slow motion. Results show that the phase space structure and particle chaos depend sensitively upon Py and bn, but are independent of h.

Burkhart, Grant R.↗
NASA NTRS1991

Development of N/P AlGaAs free-standing top solar cells for tandem applications

The combination of a free standing AlGaAs top solar cell and an existing bottom solar cell is the highest performance, lowest risk approach to implementing the tandem cell concept. The solar cell consists of an AlGaAs substrate layer, an AlGaAs base layer, an AlGaAs emitter, and an ultra-thin AlGaAs window layer. The window layer is compositionally graded which minimizes reflection at the window layer/emitter interface and creates a built-in electric field to improve quantum response in the blue region of the spectrum. Liquid phase epitaxy (LPE) is the only viable method to produce this free standing top solar cell. Small (0.125 sq cm), transparent p/n AlGaAs top solar cells were demonstrated with optimum bandgap for combination with a silicon bottom solar cell. The efficiency of an AlGaAs/Si stack using the free standing AlGaAs device upon an existing silicon bottom solar cell is 24 pct. (1X, Air Mass Zero (AM0). The n/p AlGaAs top solar cell is being developed in order to facilitate the wiring configuration. The two terminal tandem stack will retain fit, form, and function of existing silicon solar cells. Progress in the development of large area (8 and 16 sq cm), free standing AlGaAs top solar cells is discussed.

Negley, Gerald H.↗