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 613 records · Page 34

Computational techniques for flows with finite-rate condensation

A computational method to simulate the inviscid two-dimensional flow of a two-phase fluid was developed. This computational technique treats the gas phase and each of a prescribed number of particle sizes as separate fluids which are allowed to interact with one another. Thus, each particle-size class is allowed to move through the fluid at its own velocity at each point in the flow field. Mass, momentum, and energy are exchanged between each particle class and the gas phase. It is assumed that the particles do not collide with one another, so that there is no inter-particle exchange of momentum and energy. However, the particles are allowed to grow, and therefore, they may change from one size class to another. Appropriate rates of mass, momentum, and energy exchange between the gas and particle phases and between the different particle classes were developed. A numerical method was developed for use with this equation set. Several test cases were computed and show qualitative agreement with previous calculations.

Candler, Graham V.↗

Exploration of a Novel Technique for Waste Heat Recovery Through Molecular Dynamics: Influence of Wettability and Electric Field on Water and Water-Based Nanofluids

Most of the energy produced globally comes by way of a heat engine. The Carnot principle places a limit as to how thermodynamically efficient a heat engine can be. There is no heat engine that can be 100% thermodynamically efficient and as such a substantial proportion of all heat supplied to a heat engine is lost as waste heat. Waste heat therefore is a large energy source ready to be properly utilized. Herein, a novel approach for converting waste heat to electricity is discussed. It involves the use of the liquid to vapor phase change of a material dielectric (water) or electrolyte (nanofluid) in the embodiment of a capacitor for direct thermal to electrostatic energy conversion. While this method of waste heat recovery could potentially be added to the ever expanding portfolio of energy conversion techniques, a number of aspects must be addressed before it can be brought into practice. Water was seen as an ideal dielectric phase change material given its high relative permittivity ratio when in the liquid form as compared to its vapor form. However, given its short voltage holdoff time the phase change of water would need to occur rapidly. This brings up concerns of explosive boiling. Herein, molecular dynamics analysis into the explosive boiling behavior of thin water films gave more insight into how the interaction between the surface and liquid affected explosive boiling onset time. A Lennard-Jones potential with one interaction site and a Morse potential with three interaction sites between water and solid substrate were used. It was found generally that a stronger interaction between water film and substrate led to faster explosive boiling onset times but an increase in the number of interaction sites delayed explosive boiling, even at the same wettability (contact angle). Understanding changes in the density and enthalpy of vaporization of a liquid dielectric such as water in the presence of an electric field is of importance due to the electrostatic nature of the waste heat conversion method under consideration. Specifically, if both density and enthalpy of vaporization are increased, the thermodynamic efficiency of the waste heat conversion method under consideration is decreased. Electric field effects are explored herein via molecular dynamics using two water models, the TIP4P-Ew and SWM4-NDP. The SWM4-NDP model is polarizable while the TIP4P-Ew model is not, which allows for a determination of the importance of model polarizability (i.e. variation in water model dipole moment) on these two properties of water when subjected to an electric field. Herein it was found that both water models respond similarly in terms of density and vaporization enthalpy variance upon the introduction of an electric field. Comparison was also made to the pressure induced by the electric field (electrostriction pressure) by way of a density comparison and it was found that the predicted electrostriction pressure overestimates the pressure experienced by water. Water by itself has a high enthalpy of vaporization, which limits the efficiency of the newly proposed conversion method. Research both experimental and through simulation has shown that the vaporization enthalpy of nanofluids can be engineered via nanoparticle size and material selection. An avenue less explored is manipulating the enthalpy of vaporization by altering the interaction strength between the nanoparticles and the base fluid. In practice this could be achieved through the addition of coatings to the nanoparticles to alter their wettability to the base fluid. This was explored by using a Lennard-Jones potential and Morse potential to model the interaction between base fluid (water) and the nanoparticle. For nanoparticles 2nm in diameter and at weight percentages up to 6%, the change in vaporization enthalpy due to alterations of the interaction strength between nanoparticle and base fluid was not significant (less than a 1% difference) when compared to the effect of altering the weight percentage of nanoparticles in the nanofluid or introducing an electric field. However, the effect of wettability may still become important at other nanoparticle concentrations and sizes. In all, the studies presented here further the understanding of phase change and thermodynamic properties of water and water based nanofluids under an electrostatic field which will help inform the development of a novel approach to waste heat conversion. The reduction of waste heat will improve energy sustainability outlooks.

30 DIRECT ENERGY CONVERSION↗

Multipath analysis diffraction calculations

This report describes extensions of the Kirchhoff diffraction equation to higher edge terms and discusses their suitability to model diffraction multipath effects of a small satellite structure. When receiving signals, at a satellite, from the Global Positioning System (GPS), reflected signals from the satellite structure result in multipath errors in the determination of the satellite position. Multipath error can be caused by diffraction of the reflected signals and a method of calculating this diffraction is required when using a facet model of the satellite. Several aspects of the Kirchhoff equation are discussed and numerical examples, in the near and far fields, are shown. The vector form of the extended Kirchhoff equation, by adding the Larmor-Tedone and Kottler edge terms, is given as a mathematical model in an appendix. The Kirchhoff equation was investigated as being easily implemented and of good accuracy in the basic form, especially in phase determination. The basic Kirchhoff can be extended for higher accuracy if desired. A brief discussion of the method of moments and the geometric theory of diffraction is included, but seems to offer no clear advantage in implementation over the Kirchhoff for facet models.

Statham, Richard B.↗

Anisotropic index of refraction and structural properties of hexagonal boron nitride epilayers probed by spectroscopic ellipsometry

The anisotropic index of refraction of 200 μm thick boron-10 enriched hexagonal boron nitride (h-BN) freestanding epilayers grown by metalorganic chemical vapor deposition has been measured using spectroscopic ellipsometry in the UV (4.0–5.1 eV) spectral range. It was found that the index of refraction for the polarization mode with an electric field perpendicular to the c-axis (ordinary, n o ) is much higher than that with an electric field parallel to the c-axis (extraordinary, n e ). By inclusion of turbostratic- (t-) phase layers within h-BN having an average inclination angle (θ) with respect to the ideal c-plane, a simple method for quantifying θ has been deduced. Our results revealed that the presence of t-phase layers decreases the optical anisotropy of h-BN and that a signature of improved crystalline quality is an increase in the ordinary index of refraction (n o ) as a result of the average incline angle θ approaching 0° and predicted that n o = 2.7 and n e = 1.5 at 280 nm for single crystalline h-BN epilayers. More importantly, our results demonstrated that spectroscopic ellipsometry is an effective technique for characterizing the crystalline quality of h-BN epilayers with the advantages of being noninvasive and highly sensitive.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Neutral beam prompt loss in LTX-β

Prompt loss of beam injected fast ions approaches 100% in lithium tokamak experiment-beta (LTX-β) discharges, though significantly improved confinement is expected for the higher current plasmas made available by a recent upgrade to the Ohmic heating power supply. Modeling of fast ions using TRANSP/NUBEAM finds a maximum coupled beam fraction of 76% at the near-term limits of the LTX-β operating space. The full ion orbit code POET is employed to validate NUBEAM results against possible non-adiabatic effects on fast ion orbits, but corrections to the prompt loss fraction due to collisionless transport are found to be small. The graphical method code CONBEAM is used to investigate the topology of fast ion phase space as it relates to neutral beam deposition, and counter-injected NBI is considered as a way to access a region of high field side beam deposition. A metric is developed within the CONBEAM using a beam filament model to estimate the prompt loss fraction and shown to agree well with both POET and NUBEAM, enabling near real-time analysis and potential feedback to operators between plasma discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Universal Tripartite Entanglement in One-Dimensional Many-Body Systems

Motivated by conjectures in holography relating the entanglement of purification and reflected entropy to the entanglement wedge cross section, we introduce two related non-negative measures of tripartite entanglement g and h. We prove structure theorems which show that states with nonzero g or h have nontrivial tripartite entanglement. We then establish that in one dimension these tripartite entanglement measures are universal quantities that depend only on the emergent low-energy theory. For a gapped system, we argue that either g≠0 and h=0 or g=h=0, depending on whether the ground state has long-range order. For a critical system, we develop a numerical algorithm for computing g and h from a lattice model. We compute g and h for various CFTs and show that h depends only on the central charge whereas g depends on the whole operator content.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Real-Time Lossless Compression for Ultra-High-Density Synchrophasor and Point on Wave Data

Modern advanced Phasor Measurement Units (PMUs) are developed with ultra-high reporting rates to meet the demand for monitoring the power systems dynamics in detail. Due to the large volume of data, the communication and storage systems are seriously challenged with the presence of Ultra-High-Density (UHD) synchrophasor and Point on Wave (POW) data. Therefore, it is an urgent task to compress the UHD data for more efficient communication and data storage. This paper proposes several methods to compress the synchrophasor and POW data in a lossless manner. First, an Improved-Time-Series-Special Compression (ITSSC) method is proposed to compress the UHD frequency data. Second, a Delta-difference Huffman method is combined with the TSSC algorithm to compress the UHD phase angle data. Finally, a cyclical high-order delta modulation method is proposed to compress the UHD POW data. The proposed models are extensively tested and compared with different existing lossless compression algorithms using the field-collected synchrophasor and POW data at different reporting rates. The results indicate that the proposed algorithms are efficient in performing lossless compression for the UHD synchrophasor and POW data in real time.

42 ENGINEERING↗

Multigrid Optimization for Large-Scale Ptychographic Phase Retrieval

Ptychography is a popular imaging technique that combines diffractive imaging with scanning microscopy. The technique consists of a coherent beam that is scanned across an object in a series of overlapping positions, leading to reliable and improved reconstructions. Ptychographic microscopes allow for large fields to be imaged at high resolution at the cost of additional computational expense. Herein, we propose a multigrid-based optimization framework to reduce the computational burdens of large-scale ptychographic phase retrieval. Our proposed method exploits the inherent hierarchical structures in ptychography through tailored restriction and prolongation operators for the object and data domains. Our numerical results show that our proposed scheme accelerates the convergence of its underlying solver and outperforms the ptychographical iterative engine, a workhorse in the optics community.

97 MATHEMATICS AND COMPUTING↗

Applications of ultrafast nano-spectroscopy and nano-imaging with tip-based microscopy

Innovation in microscopy has often been critical in advancing both fundamental science and technological progress. Notably, the evolution of ultrafast near-field optical nano-spectroscopy and nano-imaging has unlocked the ability to image at spatial scales from nanometers to ångströms and temporal scales from nanoseconds to femtoseconds. This approach revealed a plethora of fascinating light-matter states and quantum phenomena, including various species of polaritons, quantum phases, and complex many-body effects. This review focuses on the working principles and state-of-the-art development of ultrafast tip-enhanced and near-field microscopy, integrating diverse optical pump-probe methods across the terahertz (THz) to ultraviolet (UV) spectral ranges. It highlights their utility in examining a broad range of materials, including two-dimensional (2D), organic molecular, and hybrid materials. The review concludes with a spatio-spectral-temporal comparison of ultrafast nano-imaging techniques, both within already well-defined domains, and offering an outlook on future developments of ultrafast tip-based microscopy and their potential to address a wider range of materials.

Zhao, Zhichen [Tongji University, Shanghai (China)↗

A method for unique identification of relativistic /greater than 0.5/ magnetic monopoles with a fast film Cerenkov detector

A magnetic monopole traversing a dielectric medium at a velocity greater than the phase velocity of light in that medium, will give rise to Cerenkov radiation with the electric field tangent to the cone generated by the photon wave propagation vector, and the magnetic field normal to that surface. This is the opposite polarization to that encountered with an electric charge. It is proposed that either by inserting a linearly polarizing layer between the radiator and the photographic emulsion, or by selecting a linearly polarizing material as the radiator, one could directly observe the field polarization by examining the photographic image and thus uniquely identify a magnetic monopole. The ability of the detector is further enhanced by the index of refraction dependence of the Cerenkov output from a magnetic monopole.

Pinsky, L. S.↗

Theory and application of calibration techniques for an NDBC directional wave measurements buoy

The National Data Buoy Center (NDBC) of the National Oceanic and Atmospheric Administration (NOAA) deployed a 10-m-diameter discus-type hull in the Pacific Ocean some 185 km southwest of Los Angeles, CA, in April 1984. Aboard this hull was an electronic system capable of acquiring, processing, and transmitting to shore directional wave measurements. For this system to produce accurate data, a number of factors had to be taken into account. These factors included noise, amplitude and phase alterations due to mechanical and electrical components, and magnetic fields arising from the hull. Comprehensive calibration and verification techniques were developed and applied to ensure data quality. The system configuration is described with emphasis on the methods used in the data processing to correct for the various factors. Examples of the resulting corrected data are given.

Steele, K. E.↗

High mass resolution, high angular acceptance time-of-flight mass spectroscopy for planetary missions under the Planetary Instrument Definition and Development Program (PIDDP)

This final report covers three years and several phases of work in which instrumentation for the Planetary Instrument Definition and Development Program (PIDDP) were successfully developed. There were two main thrusts to this research: (1) to develop and test methods for electrostatically scanning detector field-of-views, and (2) to improve the mass resolution of plasma mass spectrometers to M/delta M approximately 25, their field-of-view (FOV) to 360 degrees, and their E-range to cover approximately 1 eV to 50 keV. Prototypes of two different approaches to electrostatic scanning were built and tested. The Isochronous time-of-flight (TOF) and the linear electric field 3D TOF devices were examined.

Young, David T.↗

Subgrid Combustion Modeling for the Next Generation National Combustion Code

In the first year of this research, a subgrid turbulent mixing and combustion methodology developed earlier at Georgia Tech has been provided to researchers at NASA/GRC for incorporation into the next generation National Combustion Code (called NCCLES hereafter). A key feature of this approach is that scalar mixing and combustion processes are simulated within the LES grid using a stochastic 1D model. The subgrid simulation approach recovers locally molecular diffusion and reaction kinetics exactly without requiring closure and thus, provides an attractive feature to simulate complex, highly turbulent reacting flows of interest. Data acquisition algorithms and statistical analysis strategies and routines to analyze NCCLES results have also been provided to NASA/GRC. The overall goal of this research is to systematically develop and implement LES capability into the current NCC. For this purpose, issues regarding initialization and running LES are also addressed in the collaborative effort. In parallel to this technology transfer effort (that is continuously on going), research has also been underway at Georgia Tech to enhance the LES capability to tackle more complex flows. In particular, subgrid scalar mixing and combustion method has been evaluated in three distinctly different flow field in order to demonstrate its generality: (a) Flame-Turbulence Interactions using premixed combustion, (b) Spatially evolving supersonic mixing layers, and (c) Temporal single and two-phase mixing layers. The configurations chosen are such that they can be implemented in NCCLES and used to evaluate the ability of the new code. Future development and validation will be in spray combustion in gas turbine engine and supersonic scalar mixing.

Menon, Suresh↗

Dynamic Characterization of the Crew Module Uprighting System for NASA’s Orion Crew Module

The Orion Crew Module Uprighting System is a set of five airbag that are responsible for the uprighting of the crew module in the case of an inverted splashdown. A series of tests during the Underway Recovery Test 7 (URT-7) were run in preparation for the Artemis I mission, where the dynamic characterization of the CMUS in an ocean wave environment was performed. A Datawell Waverider DWR-G4 wave buoy was deployed to the characterize the wave environment during these tests. The heave measurements from this buoy were projected to the Orion Crew Module Buoyancy Test Article location by two different methods: (1) directly time-shifting the data, and (2) performing a frequency-domain, phase-shifting operation. Results demonstrate that the phase-shifting operation led to better correlation with the true crew module response to wave excitation as compared with the purely time-shifted method. Additionally, a novel approach to localize an object in a bidirectional wave field based on its heave response is presented and validated with URT- 7 data. Given a wave measurement device at a known location, one can estimate the relative distance to another object based solely off its heave response. Results show that if signals have sufficiently good correlation, this method can be used to estimate the relative separation between two objects in the same wave field.

Vertaska, Ivan Bertaska↗

Measurement of Phased Array Point Spread Functions for Use with Beamforming

Microphone arrays can be used to localize and estimate the strengths of acoustic sources present in a region of interest. However, the array measurement of a region, or beam map, is not an accurate representation of the acoustic field in that region. The true acoustic field is convolved with the array s sampling response, or point spread function (PSF). Many techniques exist to remove the PSF's effect on the beam map via deconvolution. Currently these methods use a theoretical estimate of the array point spread function and perhaps account for installation offsets via determination of the microphone locations. This methodology fails to account for any reflections or scattering in the measurement setup and still requires both microphone magnitude and phase calibration, as well as a separate shear layer correction in an open-jet facility. The research presented seeks to investigate direct measurement of the array's PSF using a non-intrusive acoustic point source generated by a pulsed laser system. Experimental PSFs of the array are computed for different conditions to evaluate features such as shift-invariance, shear layers and model presence. Results show that experimental measurements trend with theory with regard to source offset. The source shows expected behavior due to shear layer refraction when observed in a flow, and application of a measured PSF to NACA 0012 aeroacoustic trailing-edge noise data shows a promising alternative to a classic shear layer correction method.

Bahr, Chris↗

Method of Distributions for Two‐Phase Flow in Heterogeneous Porous Media

Abstract Multiscale heterogeneity and insufficient characterization data for a specific subsurface formation of interest render predictions of multi‐phase fluid flow in geologic formations highly uncertain. Quantification of the uncertainty propagation from the geomodel to the fluid‐flow response is typically done within a probabilistic framework. This task is computationally demanding due to, for example, the slow convergence of Monte Carlo simulations (MCS), especially when computing the tails of a distribution that are necessary for risk assessment and decision‐making under uncertainty. The frozen streamlines method (FROST) accelerates probabilistic predictions of immiscible two‐phase fluid flow problems; however, FROST relies on MCS to compute the travel‐time distribution, which is then used to perform the transport (phase saturation) computations. To alleviate this computational bottleneck, we replace MCS with a deterministic equation for the cumulative distribution function (CDF) of travel time. The resulting CDF‐FROST approach yields the CDF of the saturation field without resorting to sampling‐based strategies. Our numerical experiments demonstrate the high accuracy of CDF‐FROST in computing the CDFs of both saturation and travel time. For the same accuracy, it is about 5 and 10 times faster than FROST and MCS, respectively.

Yang, Hyung Jun↗

The magnetic, electronic, and light-induced topological properties in two-dimensional hexagonal FeX 2 (X = Cl, Br, I) monolayers

Using Floquet–Bloch theory, we propose to realize chiral topological phases in two-dimensional (2D) hexagonal FeX 2 (X = Cl, Br, I) monolayers under irradiation of circularly polarized light. Such 2D FeX 2 monolayers are predicted to be dynamically stable and exhibit both ferromagnetic and semiconducting properties. To capture the full topological physics of the magnetic semiconductor under periodic driving, we adopt ab initio Wannier-based tight-binding methods for the Floquet–Bloch bands, with the light-induced bandgap closings and openings being obtained as the light field strength increases. The calculations of slabs with open boundaries show the existence of chiral edge states. Interestingly, the topological transitions with branches of chiral edge states changing from zero to one and from one to two by tuning the light amplitude are obtained, showing that the topological Floquet phase of high Chern number can be induced in the present Floquet–Bloch systems.

36 MATERIALS SCIENCE↗

Atomic-Scale Structural Mapping of Active Sites in Monolayer PGM-Free Catalysts by Low-Voltage 4D-STEM

Two-dimensional (2D) materials have attracted a large amount of attention in both basic and applied fields, and scanning transmission electron microscopy (STEM) is often uniquely well-suited for characterizing the atomic-scale structure of these materials [1-4]. As a result, STEM is poised to significantly impact progress on platinum group metal (PGM)-free catalysts, which are currently under intense development to enable low-cost, commercially viable hydrogen fuel cells [5]. While recent advancements have resulted in fuel cell performance comparable to Pt catalysts by some measures [6], cell durability remains a significant challenge, limiting practical applications [7]. Catalytically active sites in PGM-free materials are proposed to be FeN4 complexes embedded in a graphene lattice (Fig. 1b) within layered or other larger materials, but this is still under debate largely due to the range of potential actives sites predicted by computational methods and lack of methods for directly validating these models [5]. Fundamental insights into the atomic structure and resulting degradation pathways of proposed active sites are therefore needed to fully understand and control cell performance and durability [6].2D materials typically make ideal samples for STEM, but those within PGM-free catalysts present additional challenges since these materials are often defect-rich, with a high density of edges, dopant atoms, etc., which significantly increase susceptibility to beam damage at standard operating voltages. This makes analysis of potential FeN4 active sites particularly challenging, since a large proportion of Fe exists at edge sites where beam-induced atomic displacements can prohibit high-resolution structural characterization [6]. Conventional dark-field imaging compounds this problem by producing less signal for a given dose and being less sensitive to light elements than dose-efficient phase contrast imaging techniques such as those enabled by four-dimensional (4D)-STEM [8-10] (Fig. 1a). Consequently, active site structural analysis is often left to methods such as low-resolution imaging combined with quantum chemical calculations [6], which hinders accurate determination of reaction and degradation mechanisms.Here, we demonstrate direct atomic-scale structural mapping of FeN4 sites by performing low-voltage 4D-STEM on a model PGM-free catalyst system with many exposed monolayer regions. To accomplish this, we pair a 30 keV aberration-corrected probe with a fast pixelated detector that has optimal performance at low beam voltages [11]. This enables us to simultaneously image light and heavy elements with high signal-to-noise by center-of-mass analysis (Fig. 1c) while minimizing beam-induced atomic displacements at sensitive sites. The monolayer nature of these materials additionally allows for experimental validation by direct comparison with multislice simulations [12] of model structures (Fig. 1d-e). This work demonstrates how low-voltage 4D-STEM will provide new insights into the atomic-scale structure and degradation mechanisms of active sites in PGM-free catalysts, facilitating the development of low-cost hydrogen fuel cells and other energy conversion technologies in the future [13].

Zachman, Michael↗