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 577 records · Page 32

Microscopic Lagrangian description of warm plasmas. III - Nonlinear wave-particle interaction

The averaged-Lagrangian method is applied to nonlinear wave-particle interactions in an infinite, homogeneous, magnetic-field-free plasma. The specific example of Langmuir waves is considered, and the combined effects of four-wave interactions and wave-particle interactions are treated. It is demonstrated how the latter lead to diffusion in velocity space, and the quasilinear diffusion equation is derived. The analysis is generalized to the random phase approximation. The paper concludes with a summary of the method as applied in Parts 1-3 of the paper.

Galloway, J. J.↗

Liquid Phase Modeling in Porous Media: Adsorption of Methanol and Ethanol in H-MFI in Condensed Water

Zeolites are used in the chemical and separation industries for their exceptional selectivity, adsorption capacity, regenerability, and stability in gas and liquid phase processing. Here, we developed an explicit solvation method for predicting solvent/condensed phase effects on adsorption free energies in microporous media such as zeolites based on the hybrid quantum mechanical/molecular mechanical free energy perturbation (QM/MM-FEP) technique. Our explicit solvation method for zeolite systems, called eSZS, aims to capture site-specific interactions during the adsorption process at the Brønsted acid sites of H-MFI zeolite while still considering the diverse configuration space of the solvent molecules. This strategy is ideal for chemical reactions or adsorbates that interact with the microporous medium in few distinct adsorbate/transition state configurations, i.e., the harmonic or similar approximations are acceptable for the adsorbate/transition state while such approximations break down for the solvent molecules that require extensive configuration space sampling. In this way, our approach effectively overcomes the limitations of implicit solvation models and classical force field methods for describing solvation effects on chemical reactions within porous materials such as zeolites. Specifically, in this study, we investigated various aspects of our hybrid QM/MM approach, including QM cluster size dependencies in a periodic electrostatically embedded cluster model (PEECM), rules for link atoms at the QM/MM boundary, and functional and basis set considerations for converged and reasonably accurate gas and aqueous phase methanol and ethanol adsorption free energy predictions in H-MFI. For gas phase adsorption of methanol and ethanol in H-MFI at a Brønsted acid site in T12 position, we compute adsorption free energies at 298 K of −0.61 and −0.75 eV, respectively, using a PEECM containing 50 Si and 1 Al atom with ωB97x-D/def2-TZVP level of theory. For solvent effect calculations, we sample the aqueous phase using grand canonical Monte Carlo (GCMC) simulations to (1) obtain a mean field of electrostatic interactions in the reaction system and (2) perform a rigorous free energy perturbation calculation. Similar to the experimentally and computationally observed endergonic solvation effects observed for hydrocarbon adsorption on metal surfaces, we also observe that a condensed aqueous environment destabilizes methanol and ethanol at these acid sites in H-MFI at 298 K. Specifically, the computed solvation free energies of adsorption (ΔΔG solv ) for methanol and ethanol are +0.44 and +0.54 eV, respectively. From this study, it is evident that adsorbates (methanol and ethanol) are competing with water for adsorption space inside the H-MFI zeolite, leading to an endergonic solvation effect. Here, we expect that the endergonic, aqueous solvent effect during adsorption in microporous zeolites is highly tunable by changing the pore size and hydrophobicity of the microporous material as this will affect the water density inside the pore structure.

Adsorption↗

Numerical Predictions of Dust-Induced Heat Flux Augmentation in Hypersonic Blunt-Body Flows Using a Discontinuous Galerkin Multiphase Flow Solver

Recent interest in human-scale missions to Mars has motivated the need for high-fidelity simulations of reentry flows. During a dust storm, there can be high levels of suspended dust in the Martian atmosphere, which cannot only enhance erosion of thermal protection systems but also transfer energy and momentum to the shock layer, thereby significantly augmenting the surface heat flux. Second-order finite-volume schemes are typically employed for hypersonic flow simulations, but such schemes suffer from a number of disadvantages. An attractive alternative is discontinuous Galerkin methods, which benefit from arbitrarily high spatial order of accuracy, geometric flexibility, and other properties. To enable accurate computations of high-speed particle-laden flows, an Euler-Lagrange methodology was developed in which the Eulerian field of the carrier gas is calculated using a discontinuous Galerkin scheme while the disperse phase is treated with Lagrangian particle tracking. We discuss challenges associated with coupling these two formulations and how to handle them. Momentum and energy transfer between the carrier gas and the particle phase is considered, and the importance of accounting for interparticle collisions is assessed. In addition, we describe the physical model of the particle phase and examine effects of its uncertainties on the numerical solution. We demonstrate the performance of the Euler-Lagrange method in representative testcases, with focus on the accurate prediction of particle trajectories and heating augmentation. Quantitative comparisons with experiments are provided.

Ching, Eric J.↗

Modeling quasielastic interactions of monoenergetic kaon decay-at-rest neutrinos

Monoenergetic muon neutrinos at 236 MeV are readily produced in intense medium-energy proton facilities (≳2–3GeV) when a positive kaon decays at rest (KDAR; K + →μ + ν μ ). These neutrinos provide a unique opportunity to both study the neutrino interaction and probe the nucleus with a monoenergetic weak-interaction-only tool. We present cross-section calculations for quasielastic scattering of these 236-MeV neutrinos off 12 C and 40 Ar, paying special attention to low-energy aspects of the scattering process. Our model takes the description of the nucleus in a mean-field approach as the starting point, where we solve Hartree-Fock equations using a Skyrme type nucleon-nucleon interaction. Thereby, we introduce long-range nuclear correlations by means of a continuum random-phase approximation (CRPA) framework where we solve the CRPA equations using a Green's function method. The model successfully describes (e,e') data on 12 C and 40 Ca in the kinematic region that overlaps with the KDAR ν μ phase space. In addition to these results, we present future prospects for precision KDAR cross-section measurements and applications of our calculations in current and future experiments that will utilize these neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Machine learning the electric field response of condensed phase systems using perturbed neural network potentials

Abstract The interaction of condensed phase systems with external electric fields is of major importance in a myriad of processes in nature and technology, ranging from the field-directed motion of cells (galvanotaxis), to geochemistry and the formation of ice phases on planets, to field-directed chemical catalysis and energy storage and conversion systems including supercapacitors, batteries and solar cells. Molecular simulation in the presence of electric fields would give important atomistic insight into these processes but applications of the most accurate methods such as ab-initio molecular dynamics (AIMD) are limited in scope by their computational expense. Here we introduce Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) to push back the accessible time and length scales of such simulations. We demonstrate that important dielectric properties of liquid water including the field-induced relaxation dynamics, the dielectric constant and the field-dependent IR spectrum can be machine learned up to surprisingly high field strengths of about 0.2 V Å −1 without loss in accuracy when compared to ab-initio molecular dynamics. This is remarkable because, in contrast to most previous approaches, the two neural networks on which PNNP MD is based are exclusively trained on molecular configurations sampled from zero-field MD simulations, demonstrating that the networks not only interpolate but also reliably extrapolate the field response. PNNP MD is based on rigorous theory yet it is simple, general, modular, and systematically improvable allowing us to obtain atomistic insight into the interaction of a wide range of condensed phase systems with external electric fields.

Science & Technology - Other Topics↗

Comparison of the dielectric and magnetocaloric properties of bulk and film of GdFe 0.5 Cr 0.5 O 3

Here, we report a comparison of the magnetic, magnetocaloric, and dielectric properties of 50% iron substituted GdCrO 3 (GdFe 0.5 Cr 0.5 O 3 ) bulk pellet and 960 nm thick film of GdFe 0.5 Cr 0.5 O 3 (GFCO). The 960 nm film was synthesized on a platinized-silicon substrate by chem. soln. deposition and spin-coating methods. The X-ray diffraction scans of the bulk sample and the film as well as the morphology of the film as examined by the field-emission scanning electron microscope indicate phase-pure and polycrystalline nature of these samples. XPS was used to deermine the valence states of Gd, Fe, and Cr. The temp. dependence of the dielectric const. from 225 to 700 K shows peaks at T C = 525 K for the bulk and ≈450 K for the film due to ferroelectric to paraelectric transitions, since elec. polarization vs elec. field hysteresis loops are observed. at room temp. The dielectric studies in the bulk GFCO for T > T C indicate a relaxor-like behavior. The measurements of the magnetization (M) of the samples as a function of temp. (5-350 K) and magnetic field (H) up to 7 T (=70 kOe) depict hysteresis behavior at low temps. due to the canted antiferromagnetic order of Fe 3+ /Cr 3+ below the Néel temp. of ≈275 K. The M vs H isotherms at various temps. are used to det. and compare the magnetic entropy change (-ΔS) and relative cooling power (RCP) of the two samples, yielding (-ΔS) = 30.7 J/kg K (18.8 J/kg K) and RCP = 566.5 J/kg (375 J/kg) for the bulk (960 nm film) samples of GFCO at 7 K and 7 T, resp. The plot of RCP vs T shows that magnetic cooling for this system is most effective for T < 30 K. Comparatively smaller magnitudes of (-ΔS) and RCP for the film vis-a-vis the bulk sample of GFCO scale with its reduced magnetization. This suggests that further improvements in the quality of the films are needed to improve their magnetization and hence their magnetocaloric properties, possibly making them useful for on-chip cooling in miniaturized devices.

36 MATERIALS SCIENCE↗

Many routes to ferroelectric HfO 2 : A review of current deposition methods

Although 10 years have passed since the initial report of ferroelectricity in hafnia (HfO 2 ), researchers are still intensely fascinated by this material system and the promise it holds for future applications. A wide variety of deposition methods have been deployed to create ferroelectric HfO 2 thin films such as atomic layer deposition, chemical solution deposition, and physical vapor deposition methods such as sputtering and pulsed laser deposition. Process and design parameters such as deposition temperature, precursor choice, target source, vacuum level, reactive gases, substrate strain, and many others are often integral in stabilizing the polar orthorhombic phase and ferroelectricity. We examine processing parameters across four main different deposition methods and their effect on film microstructure, phase evolution, defect concentration, and resultant electrical properties. The goal of this review is to integrate the process knowledge collected over the past 10 years in the field of ferroelectric HfO 2 into a single comprehensive guide for the design of future HfO 2 -based ferroelectric materials and devices.

36 MATERIALS SCIENCE↗

A statistical formulation of one-dimensional electron fluid turbulence

A one-dimensional electron fluid model is investigated using the mathematical methods of modern fluid turbulence theory. Non-dissipative equilibrium canonical distributions are determined in a phase space whose co-ordinates are the real and imaginary parts of the Fourier coefficients for the field variables. Spectral densities are calculated, yielding a wavenumber electric field energy spectrum proportional to k to the negative second power for large wavenumbers. The equations of motion are numerically integrated and the resulting spectra are found to compare well with the theoretical predictions.

Fyfe, D.↗

A scatter model for vegetation up to Ku-band

A scatter model is developed based on the matrix doubling method for volume scattering and the Kirchhoff method in rough surface scattering. Scattering from vegetation is assumed to be dominated by leaves and a single leaf is modeled by a thin dielectric disk. In developing the phase matrix for the disk, field within the disk is taken to be constant over the disk thickness, but phase changes across the surface of the disk are accounted for. Comparisons of this scatter model with radar measurements indicate good agreements in polarization, angular trends, and frequency up to Ku-band. This represents a considerable improvement over low-frequency scatter models which are valid up to S-band.

Eom, H. J.↗

Remote sensing of the atmospheric aerosols and noctilucent clouds from space

The method for the determination of the optical density of tropospheric and mesospheric aerosols and for the estimation of the scattering phase function in the forward directions (0 to 30 deg) is proposed. The method is based on measurements of the brightness of the twilight horizon with the high resolution limb camera having five separated fields of view. In some detail the features of the camera are described and the mathematical aspects of remote sensing data inversion are discussed.

Avaste, O.↗

ATD-2 Benefits Mechanism

NASA has been developing and demonstrating a suite of decision support capabilities for integrated arrival, departure, and surface (IADS) operations in a metroplex environment. The effort is being made in three phases, under NASA’s Airspace Technology Demonstration 2 (ATD-2) sub-project, through a close partnership with the Federal Aviation Administration (FAA), air carriers, airport, and general aviation community. The Phase 1 Baseline IADS capabilities provide enhanced operational efficiency and predictability of flight operations through data exchange and integration, tactical surface metering, and automated coordination of release time of controlled flights for overhead stream insertion. The Phase 2 Fused IADS capabilities include the fusion of strategic and tactical surface metering, Atlanta Center airspace tactical scheduling, Electronic Flight Data (EFD) integration, Terminal Flight Data Manager (TFDM) Terminal Publication (TTP) prototype, and Mobile App for General Aviation (GA) community. In the Phase 2 field evaluation, strategic surface metering provides advance notice of metering and additional stability to the assigned gate holds. The users of the IADS system in Phases 1 and 2 include the personnel at Charlotte Douglas International Airport (CLT) air traffic control tower, American Airlines ramp tower, CLT terminal radar approach control (TRACON), and Washington and Atlanta Center. This document describes the ATD-2 benefits mechanism used to assess the Phases 1 and 2 IADS capabilities and field evaluation conducted at CLT since September 2017. The ATD-2 benefits mechanism mainly consists of surface metering and overhead stream insertion. This document provides detailed calculation methods of major benefit metrics, such as fuel savings, gas emissions savings, and engine runtime reduction, which can be obtained through surface metering, gate hold of Approval Request (APREQ) flights prior to pushback, and the renegotiation of release time while taxiing. As of March 31, 2020, it is estimated that 5,075,981 pounds of fuel savings and 15,634,022 pounds of CO2 emission reduction have been achieved so far, with a reduction of 3,832 hours in total engine runtime. The amount of CO2 savings is estimated to be equivalent to planting 116,254 urban trees. The pre- and post-metering comparison results using FAA’s Aviation System Performance Metrics (ASPM) data have also shown that the surface metering had no negative impact on the on-time arrival performance of both outbound and inbound flights at CLT.

ATD-2↗

ATD-2 Benefits Mechanism

NASA has been developing and demonstrating a suite of decision support capabilities for integrated arrival, departure, and surface (IADS) operations in a metroplex environment. The effort is being made in three phases, under NASA’s Airspace Technology Demonstration 2 (ATD-2) sub-project, through a close partnership with the Federal Aviation Administration (FAA), air carriers, airport, and general aviation community. The Phase 1 Baseline IADS capabilities provide enhanced operational efficiency and predictability of flight operations through data exchange and integration, tactical surface metering, and automated coordination of release time of controlled flights for overhead stream insertion. The Phase 2 Fused IADS capabilities include the fusion of strategic and tactical surface metering, Atlanta Center airspace tactical scheduling, Electronic Flight Data (EFD) integration, Terminal Flight Data Manager (TFDM) Terminal Publication (TTP) prototype, and Mobile App for General Aviation (GA) community. In the Phase 2 field evaluation, strategic surface metering provides advance notice of metering and additional stability to the assigned gate holds. The users of the IADS system in Phases 1 and 2 include the personnel at Charlotte Douglas International Airport (CLT) air traffic control tower, American Airlines ramp tower, CLT terminal radar approach control (TRACON), and Washington and Atlanta Center. This document describes the ATD-2 benefits mechanism used to assess the Phases 1 and 2 IADS capabilities and field evaluation conducted at CLT since September 2017. The ATD-2 benefits mechanism mainly consists of surface metering and overhead stream insertion. This document provides detailed calculation methods of major benefit metrics, such as fuel savings, gas emissions savings, and engine runtime reduction, which can be obtained through surface metering, gate hold of Approval Request (APREQ) flights prior to pushback, and the renegotiation of release time while taxiing. As of April 30, 2020, it is estimated that 5,097,173 pounds of fuel savings and 15,699,292 pounds of CO2 emission reduction have been achieved so far, with a reduction of 3,831 hours in total engine runtime. The amount of CO2 savings is estimated to be equivalent to planting 116,739 urban trees. The pre- and post-metering comparison results using FAA’s Aviation System Performance Metrics (ASPM) data have also shown that the surface metering had no negative impact on the on-time arrival performance of both outbound and inbound flights at CLT.

ATD-2↗

Phased Array Characterization of Slat Noise Radiation from a High-Lift Common Research Model

A test campaign was conducted in the open-jet test section of the NASA Langley 14- by 22-Foot Subsonic Tunnel on a 10%-scale, semispan version of a High-Lift Common Research Model (CRM-HL) incorporating a leading-edge slat, trailing-edge flap and removable high-fidelity main landing gear. Far field acoustic measurements were obtained on baseline and acoustically treated model configurations using a traversing 97-microphone phased array that viewed the pressure side of the airframe. A modified version of the DAMAS deconvolution method incorporating corrections for shear layer acoustic wave decorrelation was employed to determine the locations and strengths of relevant noise sources along the span of the slat and in the vicinity of the flap edges and landing gear. The main goal of the test campaign was to use the array data to evaluate the effectiveness of various slat noise reduction concepts. It was found that slat-gap fillers produced significant broadband noise reduction for all Mach numbers and airfoil angles of attack that were investigated. Mach number scaling revealed an approximate velocity to the sixth power relationship for the slat noise at higher frequencies. Directivity measurements obtained from integration of DAMAS pressure-squared values over defined geometric zones around the model showed that baseline slat noise radiation approximated a dipole as a function of polar emission angle, becoming approximately omnidirectional with the introduction of the slat-gap filler. The test campaign clearly demonstrated the noise reduction benefits that can be obtained by applying appropriate treatments to leading edge slats on commercial transport-class aircraft.

William M Humphreys, Jr.↗

Post-fire soil respiration in late growing season (2023 and 2024), Kougarok Fire Complex, Seward Peninsula, Alaska

Field soil respiration data collected in 2023 and 2024 from burned and unburned tussock tundra sites in the Kougarok Fire Complex, near Nome, on the Seward Peninsula of Alaska. Specifically, we measured soil properties and late-growing season CO2 fluxes in patches of unique plant functional types (forbs, shrubs, and graminoids) across two years in tundra recovering from repeated wildfires over the decade. The goal was to identify the main drivers of soil respiration in Arctic tundra underlain by discontinuous permafrost that is recovering from two recent, repeated wildfires that differed in fire age and number of times burned, thereby resulting in different levels of vegetation and subsurface property changes (i.e., successional trajectories). There are five files in *.csv format with one data file and four data description files including data, dictionary, methods, terminology, and file-level metadata. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), is a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic Phase 3 project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

Santos, Fernanda [ORNL] (ORCID:0000000191555623)↗

Improved Maneuver Reconstructions for the GRAIL Orbiters

Maneuver reconstructions for the Gravity Recovery and Interior Laboratory (GRAIL) A and B lunar orbiters were improved through updates to the orbit determination filter and dynamic models. Consistent reconstructions of the 27 GRAIL A and B maneuvers from the Trans-Lunar Cruise phase in the fall of 2011 through the Transition to Science Formation phase in February 2012 were performed. Standard methods of orbit determination were applied incorporating the latest dynamic models and filter strategies developed by the GRAIL Navigation and Science Teams, including a high resolution, 420 x 420 degree and order lunar spherical harmonic gravity field model. For Trans-Lunar Cruise for GRAIL-A (TLC-A), all maneuvers executed with delta V errors below 5.50 +/- 0.50 mm/s and pointing errors below 0.25 degrees. GRAIL-A lunar orbit maneuvers had delta V errors below 30.0 mm/s and pointing errors below 0.51 degrees. For TLC-B, all maneuvers executed with delta V errors below 8.60 +/- 1.41 mm/s and pointing errors below 0.300 degrees. GRAIL-B maneuvers in lunar orbit executed with maximum delta V errors of 25.0 mm/s and pointing error of 0.43 degrees. These maneuver reconstructions will enable the GRAIL Navigation Team to better characterize the main engine performance of each spacecraft. This will help the Navigation Team to navigate low (greater than 8 km) altitude orbits during the extended mission phase in the fall of 2012.

lunar mapping↗

Floating Block Method for Quantum Monte Carlo Simulations

Quantum Monte Carlo simulations are powerful and versatile tools for the quantum many-body problem. In addition to the usual calculations of energies and eigenstate observables, quantum Monte Carlo simulations can in principle be used to build fast and accurate many-body emulators using eigenvector continuation or design time-dependent Hamiltonians for adiabatic quantum computing. Furthermore, these new applications require something that is missing from the published literature, an efficient quantum Monte Carlo scheme for computing the inner product of ground state eigenvectors corresponding to different Hamiltonians. In this work, we introduce an algorithm called the floating block method, which solves the problem by performing Euclidean time evolution with two different Hamiltonians and interleaving the corresponding time blocks. We use the floating block method and nuclear lattice simulations to build eigenvector continuation emulators for energies of 4 He, 8 Be, 12 C, and 16 O nuclei over a range of local and nonlocal interaction couplings. From the emulator data, we identify the quantum phase transition line from a Bose gas of alpha particles to a nuclear liquid.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Electric-Field Manipulation of Magnetic Chirality in a Homo-Ferro-Rotational Helimagnet

Ferro-rotational (FR) materials, renowned for their distinctive material functionalities, present challenges in the growth of homo-FR crystals (i.e., single FR domain). This study explores a cost-effective approach to growing homo-FR helimagnetic RbFe(SO 4 ) 2 (RFSO) crystals by lowering the crystal growth temperature below the T FR threshold using the high-pressure hydrothermal method. Through polarized neutron diffraction experiments, it is observed that nearly 86% of RFSO crystals consist of a homo-FR domain. Notably, RFSO displays remarkable stability in the FR phase, with an exceptionally high T FR of ≈573 K. Furthermore, RFSO exhibits a chiral helical magnetic structure with switchable ferroelectric polarization below 4 K. Importantly, external electric fields can induce a single magnetic domain state and manipulate its magnetic chirality. The findings suggest that the search for new FR magnets with outstanding material properties should consider magnetic sulfates as promising candidates.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Impacts of Synoptic-Scale Dynamics on Clouds and Radiation in High Southern Latitudes

High-latitudinal mixed-phase clouds significantly affect Earth's radiative balance. Observations of cloud and radiative properties from two field campaigns in the Southern Ocean and Antarctica were compared with two global climate model simulations. A cyclone compositing method was used to quantify “dynamics-cloud-radiation” relationships relative to the extratropical cyclone centers. Observations show larger asymmetry in cloud and radiative properties between western and eastern sectors at McMurdo compared with Macquarie Island. Most observed quantities at McMurdo are higher in the western (i.e., post-frontal) than the eastern (frontal) sector, including cloud fraction, liquid water path (LWP), net surface shortwave and longwave radiation (SW and LW), except for ice water path (IWP) being higher in the eastern sector. Here, the two models were found to overestimate cloud fraction and LWP at Macquarie Island but underestimate them at McMurdo Station. IWP is consistently underestimated at both locations, both sectors, and in all seasons. Biases of cloud fraction, LWP, and IWP are negatively correlated with SW biases and positively correlated with LW biases. The persistent negative IWP biases may have become one of the leading causes of radiative biases over the high southern latitudes, after correcting the underestimation of supercooled liquid water in the older model versions. By examining multi-scale factors from cloud microphysics to synoptic dynamics, this work will help increase the fidelity of climate simulations in this remote region.

58 GEOSCIENCES↗