Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase modeling”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Improved pore network models to simulate single-phase flow in porous media by coupling with lattice Boltzmann method

In this paper, different pore network models to simulate single-phase flow in porous media are built and their accuracy is evaluated. In addition to the conventional pore network model (CPNM) which consists of regular pore bodies and throat bonds, three improved pore network models (IPNMs) are developed allowing to better describing the real pore and throat geometry. The first improved pore network model (IPNM1) replaces the regular throat bond with a throat bond showing the real throat cross section. The second improvement (IPNM2) uses a series of sub-throat bonds with varying cross sections to better describe the real throat geometry, which is firstly proposed in this paper. The third model (IPNM3) extracts the real pore-throat-pore geometry without simplification. The conductance of fluid flow through these more realistic throat bonds is calculated by the lattice Boltzmann method (LBM). The accuracy and computational efficiency of the different pore network models are evaluated taking the LBM simulation over the whole porous medium as reference solution. The global permeability and detailed pressure distributions in the pores for the different pore network models are validated. The results show that the accuracy of the pore network model increases from CPNM to IPNM3, but at the expense of increasing computational cost. This study suggests that IPNM3 can replace a whole-domain LBM simulation with similar accuracy but much lower computational cost. As a first-order approximation the newly proposed IPNM2 is suggested as good compromise between accuracy and computational cost.

42 ENGINEERING↗
DOE OSTI2022

Modeling of high intensity orbital angular momentum beams for laser–plasma interactions

In this work, we explore the field of high intensity orbital angular momentum (OAM) beams, their generation with spiral phase mirrors, and the theory behind modeling both ideal and realistic beams. We explore OAM beam asymmetries introduced by aberrations in the beam, manufacturing defects, and bandwidth. A full three-dimensional description of the paraxial Laguerre–Gaussian (LG) modes suitable for modeling beams down to f/2 focal geometries is derived. A perturbative approach to modeling asymmetric OAM beams is introduced showing that only three LG modes are sufficient to model a wide variety of OAM asymmetries. The models are compared with experimental results followed by a discussion on the future of high intensity OAM beams in plasma physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Progress Toward Simulating Departure from Nucleate Boiling at High-Pressure Applications with Selected Wall Boiling Closures

Recently, a Eulerian-based two-fluid computational fluid dynamics (CFD) framework with a wall heat flux partitioning approach has been intensively investigated for departure from nucleate boiling (DNB) simulation under the U.S. Department of Energy–funded Consortium for Advanced Simulation of Light Water Reactors (CASL) program. Understanding of the DNB characteristics over a range of pressurized water reactor–like operating conditions and accurate prediction of boiling crisis in the nuclear power system have been grand challenges because of the large impact of DNB on reactor safety and operational economics. The ultimate goal of this task in the CASL program is to introduce a robust multiphase CFD–based DNB modeling framework that is capable of characterizing an entire boiling history in which the wall boiling mode experiences the following through multiple stages of heat transfer mode: (1) single-phase convective heat transfer, (2) nucleate boiling heat transfer, and (3) identification of the departure of nucleate boiling. To validate the CASL boiling model, we have benchmarked simulated DNB over three different flow channel configurations (pipe flow, 5 × 5 fuel bundle with mixing vane tests, and 5 × 5 fuel bundle without mixing vane tests) against experimental measurements, and the validation result with open literature is reported. The DNB detection criteria in the simulation are checked by monitoring the peak wall temperature, wall dryout factor, and net energy balance. In addition to the DNB performance test, some preliminary sensitivity results on closure model selection are reported to address the prediction capability of local void profile against measurements. The boiling simulation tested in this study exhibits a maximum deviation of 24% from the measured DNB value in a high-pressure (i.e., 138 bars) subcooled pipe flow test. The ranges of operating conditions are as follows: 1650 to 2650 kg/m 2 ·s for mass flux and 8.5 to 96 K for subcooled inlet temperature. The deviation is even reduced to 7% when the subcooled temperature is less than 40 K. Besides accuracy, base practice guidelines for DNB detection criteria are tested by monitoring three simulation variables: (1) maximum wall temperature, (2) wall dryout factor (i.e., K-value), and (3) energy balance. Numerical robustness of DNB simulation is largely achieved in most of the validation test except for a few high subcooled test cases.

42 ENGINEERING↗
DOE OSTI2020

Elliptic flow coalescence to identify the $f_0$(980) content

Here we use a simple coalescence model to generate $f_0$(980) particles for three configurations: a $s\bar{s}$ meson, a $u\bar{u}s\bar{s}$ tetraquark, and a $K^+K^–$ molecule. The phase-space information of the coalescing constituents is taken from a multiphase transport (AMPT) simulation of heavy-ion collisions. It is shown that the number of constituent quarks scaling of the elliptic flow anisotropy can be used to discern $s\bar{s}$ from $u\bar{u}s\bar{s}$ and $K^+K^–$ configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Spin dynamics in natural multiferroic pyroxene NaFeSi 2 ⁢O 6

Spin dynamics in the natural mineral aegirine, NaFeSi 2⁢ O 6 , a member of the pyroxene family, was studied by elastic and inelastic neutron scattering. Magnetization and specific-heat measurements as well as single-crystal neutron diffraction maps, taken in the temperature range 2–20 K, confirm two successive magnetic transitions at 8.8 and 5.8 K, consistent with previous studies. The observed spin-wave excitations emerge from the incommensurate magnetic Bragg peaks corresponding to the propagation vector 𝑘 ICM = (0,0.77,0), and extend up to energies of about 1.5 meV. In the low-temperature helical phase, the spin dynamics of the Fe3+ ions is well described by a simple linear spin-wave model. The observed excitations can be modeled using a spin Hamiltonian that includes three primary exchange interactions—intrachain coupling 𝐽 = 0.142⁢(2) ⁢meV, and interchain couplings 𝐽 1 = 0.083⁢(1) ⁢meV and 𝐽 2 = 0.186⁢(1) meV—and an easy-plane anisotropy 𝐷 = 0.020⁢(6)⁢ meV. Furthermore, our results show that no single exchange interaction dominates the spin dynamics. The similar strengths of the intrachain and interchain couplings point to the fact that the magnetic interactions in aegirine are three dimensional rather than confined along one direction. As a result, the system cannot be considered quasi one dimensional, as previously suggested, and calls for further investigations.

Magnetic coupling↗
DOE OSTI2020

Energy Saving Quantification on Ductless Heat Pump (DHP) in Existing Homes

In residential retrofit applications, ductless mini-split heat pumps (DHP) are often reported to have high-energy savings potential, depending on the system they are supplementing or replacing. However, recently, there have been a number of utility studies and analyses indicating these energy savings are not being achieved when the existing system is left in place for backup heating or air conditioning. The Pacific Northwest National Laboratory conducted a three-phase project to help determine which control strategies would have the most energy savings impact in various climate zones around the United States. The first phase developed the standalone simulation model for the PNNL Lab Homes and investigated the energy-saving potential from different control strategies and HVAC system configurations. The second phase focused on conducting experiments in the PNNL Lab Homes, which tested the most promising solutions that were modeled in the first phase. The third phase used the field data to calibrate the simulation model and then extrapolated the results to different climate locations and different building sizes. This report focuses on the third phase of the study, including five major parts: model calibration, parametric model setup, results, a sensitivity analysis of air leakage rate, and conclusions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Assessment of Carbonated Brine Injection as Low-risk Strategy for Geologic Carbon Storage

The main objective of this early-phase research was to evaluate the techno-economic feasibility and risk associated with combined brine and CO 2 storage in SWD wells using brine dissolution in the North Dakota portion of the Williston Basin. Three simulation studies were conducted to investigate: (1) CO 2 phases at different conditions, (2) wellbore compatibility, and (3) long-term storage fate in reservoir. (1) A simple reservoir model and injection simulations were created using data to represent the BEST (brine extraction and storage test) site, an operational SWD facility located near Watford City, North Dakota. The pressure evolution caused by CO 2 comingled in produced water injectate in a layer cake reservoir was then modeled while tracking aqueous CO 2 throughout the project. The salinity of the injection water, the salinity of the reservoir brine, and the amount of dissolved CO 2 comingled in the injection water were varied. (2) A wellbore corrosion model was performed using the CO 2 concentrations selected based on the reservoir modeling to examine the carbonated produced water impact on wellbore. (3) Reactive transport modeling was conducted with the optimal CO 2 concentration for this injection site to study the rock-fluid interactions and CO 2 fate in the reservoir. Results suggest that CO 2 dissolved in produced water can be injected without appreciably increasing subsurface pressure or leakage risks. Pressure buildup was found to vary with salinity but not with CO 2 mass fraction. Simulation results show that lower CO 2 percent mass fraction leads to a higher amount of CO 2 that can be dissolved at a higher injection salinity. Furthermore, the long-term goal of dissolution trapping in a traditional carbon storage project is accomplished from the start, mitigating risks associated with potential migration of buoyant CO 2 , so long as the reservoir pressure and temperature are used to determine the maximum mass fraction of the dissolved CO 2 .

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1974

HiMAT - A new approach to the design of highly maneuverable aircraft

Needed improvements in the maneuvering performance of combat aircraft appear to be possible through the simultaneous application of advances in various disciplines in such a way that they complement one another and magnify the benefits derived. The highly maneuverable aircraft technology (HiMAT) program is being conducted to investigate such multidisciplinary concepts. The program has three phases: preliminary studies, conceptual design studies, and the final design and construction of a test airplane. Work is now in the second phase. The test airplane will be a scaled model flown by a remotely piloted research vehicle technique. This paper outlines the HiMAT program and indicates the types of concepts being considered.

Bellman, D. R.↗
NASA NTRS1990

Differential memory in the trilinear model magnetotail

The previously proposed concept of 'differential memory' is quantitatively demonstrated using an idealized analytical model of particle dynamics in the magnetotail geometry. In this model (the 'trilinear' tail model) the magnetotail is divided into three regions. The particle orbits are solved exactly in each region, thus reducing the orbit integration to an analytical mapping. It is shown that the trilinear model reproduces the essential phase space features of the earlier model (Chen and Palmadesso, 1986), possessing well-defined entry and exit regions, and stochastic, integrable (regular), and transient orbits, occupying disjoint phase space regions. Different regions have widely separated characteristic time scales corresponding to different types of particle motion. Using the analytical model, the evolution of single-particle distribution functions is calculated.

Chen, James↗
NASA NTRS1991

The first three minutes - 1990 version

The present state of understanding of what occurred in the universe's first three minutes is reviewed. Emphasis is on the events that lead to potentially observable consequences and that are model-independent or at least generic to broad classes of models. Inflation, phase transitions, dark matter, and nucleosynthesis are summarized.

Schramm, David N.↗
NASA NTRS1993

Remote sensing estimates of cirrus particle size for tropical and midlatitude cirrus: Hexagonal crystals and ice spheres

A large discrepancy exists in current estimates of a mean cirrus particle size appropriate for calculations of the effects of these ice clouds on solar and thermal infrared radiative fluxes. For spheres with large size parameter (x = (2 pi r / lambda) is greater than 30, where r is particle radius), and moderate absorption (n(sup i) x less than 1, where n(sup i) is imaginary index of refraction for ice), the optimal effective particle radius is given by: r(sub e) = integral of r(exp 3)n(r)dr / integral of r(exp 2)n(r)dr. For the remote sensing of cirrus particle size at wavelengths of 0.83, 1.65, and 2.21 mu m, a 50 mu m ice sphere would have a size parameter of about 200, and values of n(sup i) x of 0, 0.045, and 0.06, satisfying the above conditions. However, while r(sub e) is a well-defined parameter for spheres, this cross-section area-weighted particle radius can only be extended to non-spherical particles by defining some equivalent sphere, typically an equivalent volume or equivalent cross-section area sphere. Using equivalent volume spheres, values of r(sub e) obtained over Lake Michigan on October 28, 1986, during FIRE phase I varied from 200 mu m (King Air 2D Imaging probes) to 60 mu m (Landsat reflectances at 0.83, 1.65, and 2.2 mu m), to 25 mu m (HIS spectrometer thermal emission between 8 and 12 mu m). Three major uncertainties were identified in this comparison: small ice particles missed by the 2D-C aircraft probes, uncertain ice refractive index, and uncertainties in the single scatter albedos and scattering phase functions used in the radiative calculations. Since the first FIRE cirrus results, advances have been made in all three areas. The present paper reports on improvements in the radiative modeling of ice particles at 0.83, 1.65, and 2.21 mu m wavelengths appropriate for comparisons to Landsat Thematic Mapper data. The paper also includes new results for Landsat observations of ice clouds in the eastern and western tropical Pacific.

Wielicki, Bruce A.↗
NASA NTRS2002

The Long-term Light Curves of X-ray Binaries Contain Simultaneous Periodic and Random Components

LMC X-3 and Cyg X-2 show large amplitude X-ray fluctuations that have been attributed to a warped accretion disk. Cyg X-3 displays high amplitude, apparently non-periodic oscillations. We reanalyze these systems using RXTE ASM data and time-frequency decomposition techniques. We find that the long-term variations in Cyg X-2 can be completely characterized by excursions whose durations are integer multiples of the orbital period, including one essentially identical to the reported "period" of 78 days. Cyg X-3 can be characterized in terms of integer multiples of a 71-day fundamental period unrelated to the 4.8 day orbital period, but suggestively close to the approximately equal to greater than 60 day reported precession period of the relativistic jet inferred from recent radio observations. The long-term excursions of LMC X-3 are related to each other by rational fractions, suggesting the characteristic time scale is 10.594 days, shorter than any observed excursion to date. We explore the phase space evolution of the light curves using a natural embedding and find that all three systems possess two rotation centers that organize the phase space trajectories, one of low luminosity and the other of high luminosity. The implications of this repeatable behavior on generic models of accretion disk dynamics and mass transfer variability are explored.

White, Nicholas E.↗
NASA NTRSDate not supplied

Growth Curve Parameterization of Metabolic Activity of Yeast Cells for BioSentinel

The goal of the BioSentinel small satellite payload is to measure the effect of deep space radiation on the growth and metabolic activity of yeast cells. Raw test data is generated by fluidics cards containing yeast cells rehydrated at different periods, with metabolic activity measured by the reduction of alamarBlue. Each card well has a sensor array that measures the amount of red, green, and infrared light transmitted through the yeast culture. This illumination data is then converted to absorbance values, which are further converted into concentrations. The ultimate objective is to convert these concentrations into biologically-relevant metrics that can be compared against one another to determine changes due to differential radiation exposure. Beginning with IR absorbance data (corresponding to cell density) from ground studies, three parameters from a sigmoidal growth curve were extracted and analyzed: 𝜆 (lag phase), 𝜇 (max growth rate), and A (max cell growth). The data was fit to the Gompertz model of microbial growth using non-linear regression (Minitab), as the fit error was reduced compared to the simpler logistic growth curve. Graphs showed that the data contained a discrepancy (drift) in the lag phase that is attributable to a slow, constant loss of moisture. Correcting this discrepancy by fitting the first 25 hours of the data to a power function and subtracting these values from the absorbance readings obtained a better statistical fit to the growth curve in the lag phase. A power fit was selected over a linear fit because it reflected the effects of constant volume loss. This correction to the BioSentinel data analysis pipeline will enable quantitative statistical analysis of the effect of different levels of deep space radiation on yeast cells. Future work includes automation of drift correction and curve modeling to extract these parameters directly from data.

Growth Curve↗
NASA NTRS2022

Modeling the Lunar Wake Response to a CME Using a Hybrid PIC Model

In the solar wind, a low-density wake region forms downstream of the nightside lunar surface. In this study, we use a series of 3D hybrid particle-in-cell simulations to model the response of the lunar wake to a passing coronal mass ejection (CME). Average plasma parameters are derived from the Wind spacecraft located at 1 au during three distinct phases of a passing halo (Earth-directed) CME on 2015 June 22. Each set of plasma parameters, representing the shock/plasma sheath, a magnetic cloud, and plasma conditions we call the mid-CME phase, are used as the time-static upstream boundary conditions for three separate simulations. These simulation results are then compared with results that use nominal solar wind conditions. Results show a shortened plasma void compared to nominal conditions and a distinctive rarefaction cone originating from the terminator during the CME's plasma sheath phase, while a highly elongated plasma void reforms during the magnetic cloud and mid-CME phases. Developments of electric and magnetic field intensification are also observed during the plasma sheath phase along the central wake, while electrostatic turbulence dominates along the plasma void boundaries and 2–3 lunar radii RM downstream in the central wake during the magnetic cloud and mid-CME phases. The simulations demonstrate that the lunar wake responds in a dynamic way with the changes in the upstream solar wind during a CME.

Anthony P. Rasca↗
DOE OSTI2020

X-Ray Diffraction and Electron Microscopy Studies of the Size Effects on Pressure-Induced Phase Transitions in CdS Nanocrystals

In recent years, investigations of the phase transition behavior of semiconducting nanoparticles under high pressure has attracted increasing attention due to their potential applications in sensors, electronics, and optics. However, current understanding of how the size of nanoparticles influences this pressure-dependent property is somewhat lacking. In particular, phase behaviors of semiconducting CdS nanoparticles under high pressure have not been extensively reported. Therefore, in this work, CdS nanoparticles of different sizes are used as a model system to investigate particle size effects on high-pressure-induced phase transition behaviors. In particular, 7.5, 10.6, and 39.7 nm spherical CdS nanoparticles are synthesized and subjected to controlled high pressures up to 15 GPa in a diamond anvil cell. Analysis of all three nanoparticles using in-situ synchrotron wide-angle X-ray scattering (WAXS) data shows that phase transitions from wurtzite to rocksalt occur at higher pressures than for bulk material. Bulk modulus calculations not only show that the wurtzite CdS nanomaterial is more compressible than rocksalt, but also that the compressibility of CdS nanoparticles depends on their particle size. Furthermore, sintering of spherical nanoparticles into nanorods was observed for the 7.5 nm CdS nanoparticles. Our results provide new insights into the fundamental properties of nanoparticles under high pressure that will inform designs of new nanomaterial structures for emerging applications.

Meng, Lingyao↗
NASA NTRS2013

Swept-Wing Ice Accretion Characterization and Aerodynamics

NASA, FAA, ONERA, the University of Illinois and Boeing have embarked on a significant, collaborative research effort to address the technical challenges associated with icing on large-scale, three-dimensional swept wings. The overall goal is to improve the fidelity of experimental and computational simulation methods for swept-wing ice accretion formation and resulting aerodynamic effect. A seven-phase research effort has been designed that incorporates ice-accretion and aerodynamic experiments and computational simulations. As the baseline, full-scale, swept-wing-reference geometry, this research will utilize the 65% scale Common Research Model configuration. Ice-accretion testing will be conducted in the NASA Icing Research Tunnel for three hybrid swept-wing models representing the 20%, 64% and 83% semispan stations of the baseline-reference wing. Three-dimensional measurement techniques are being developed and validated to document the experimental ice-accretion geometries. Artificial ice shapes of varying geometric fidelity will be developed for aerodynamic testing over a large Reynolds number range in the ONERA F1 pressurized wind tunnel and in a smaller-scale atmospheric wind tunnel. Concurrent research will be conducted to explore and further develop the use of computational simulation tools for ice accretion and aerodynamics on swept wings. The combined results of this research effort will result in an improved understanding of the ice formation and aerodynamic effects on swept wings. The purpose of this paper is to describe this research effort in more detail and report on the current results and status to date. 1

Broeren, Andy P.↗
DOE OSTI2023

Lattice dynamics and thermodynamics for δ -plutonium from density functional theory

Here, we present results from density functional theory (DFT) calculations of the lattice dynamics (phonons) and thermodynamics for δ-phase plutonium. The fully relativistic electronic structure is calculated assuming a three-dimensional noncollinear magnetic structure in conjunction with DFT and the general gradient approximation for the electron exchange and correlation interactions. The electronic-structure model is further enhanced by addressing strong orbital-orbital coupling via the conventional orbital-polarization (OP) scheme as has been successfully done for plutonium. The temperature dependence of the phonons is calculated within the self-consistent ab initio lattice dynamics approach. The obtained phonons compare very well with measurements although a modest overestimation of the transverse L-point [ξξξ] phonon is acknowledged. Calculated thermal vibration amplitudes and the associated Debye-Waller temperatures are close to experiments. Lattice, electronic, and magnetic contributions to the heat capacity are predicted and consistent to a few percent with that deduced from experimental data. Good agreement is only achieved when a magnetic contribution to the specific heat is recognized. The parameter-free DFT+OP electronic model is thus capable of predicting phonon properties and thermodynamic behavior of δ-phase plutonium rather accurately.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Strain-tuned quantum criticality in electronic Potts-nematic systems

Motivated by recent observations of threefold rotational symmetry breaking in twisted moiré systems, cold-atom optical lattices, quantum Hall systems, and triangular antiferromagnets, we phenomenologically investigate the strain-temperature phase diagram of the electronic 3-state Potts-nematic order. While in the absence of strain the quantum Potts-nematic transition is first-order, quantum critical points (QCP) emerge when uniaxial strain is applied, whose nature depends on whether the strain is compressive or tensile. In one case, the nematic amplitude jumps between two non-zero values while the nematic director remains pinned, leading to a symmetry-preserving metanematic transition that terminates at a quantum critical end-point. For the other type of strain, the nematic director unlocks from the strain direction and spontaneously breaks an in-plane twofold rotational symmetry, which in twisted moiré superlattices triggers an electric polarization. Such a piezoelectric transition changes from first to second-order upon increasing strain, resulting in a quantum tricritical point. Using a Hertz-Millis approach, we show that these QCPs share interesting similarities with the widely studied Ising-nematic QCP. The existence of three minima in the nematic action also leaves fingerprints in the strain-nematic hysteresis curves, which display multiple loops. At non-zero temperatures, because the upper critical dimension of the 3-state Potts model is smaller than three, the Potts-nematic transition is expected to remain first-order in 3D, but to change to second-order in 2D. As a result, the 2D strain-temperature phase diagram displays two first-order transition wings bounded by lines of critical end-points or tricritical points, reminiscent of the phase diagram of metallic ferromagnets. Furthermore, we discuss how our results can be used to unambiguously identify spontaneous Potts-nematic order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗