Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite density”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Predicting SLS Launch Environment using a Novel Multiphase Formulation

Powerful acoustic waves generated during ignition of launch vehicles may be dangerous to the vehicle, its payload, or the surrounding structures. The water-based Ignition Overpressure and Sound Suppression (IOP/SS) system at Kennedy Space Center’s (KSC) Launch Complex 39B (LC-39B) will be used to protect the Space Launch System (SLS) from the acoustic vibrations generated during launch. The IOP/SS system uses enormous amounts of water to dampen and attenuate these sound waves. To better understand the launch environment risks and to study the effectiveness of the IOP/SS system it is desirable to have time-accurate unsteady simulations of the vehicle ignition with water-based sound suppression. This paper presents results obtained with a novel, high-order accurate, and robust numerical method designed for simulating compressible multiphase flows. A positivity-preserving finite difference scheme is utilized which is formally high-order accurate and also provably robust. Robustness is critical due to the extreme nature of the flow which exhibits highly nonlinear shock and rarefaction waves interacting with liquid-gas interfaces with density ratios of the order of 1000:1. Furthermore, the high-order accuracy (and the high resolution property) is desirable for predicting wave phenomena like IOP waves since the signal can be resolved accurately and propagated long distances with fewer grid points. This finite-difference method was developed using NASA’s Launch, Ascent, and Vehicle Aerodynamics (LAVA) Cartesian immersed boundary framework. We present a validation case by applying our solver to the SLS Scale Model Acoustic Test (SMAT). The SLS SMAT is a well-instrumented 5% scale model test meant to represent the SLS at NASA KSC’s LC-39B pad. Scale IOP tests were performed with and without the sound suppression water and included many sensors which recorded the pressure waves produced during ignition. For this validation case we conduct two simulations, likewise with and without sound suppression water, and compare the SLS SMAT pressure sensor signals with our numerical signals at identical locations. Following this validation case we present a study of the SLS launch environment to examine engineering safety concerns about the mobile launch pad. Engineers at KSC redesigned the main flame deflector at LC-39B anticipating the increased loads from the SLS and to repair damage from prior Shuttle missions. This deflector redesign made use of surface pressure and temperature data from LAVA full-scale SLS simulations without the sound suppression system. The engineers were questioning the possibility of increased pressure loads on the underside of the mobile launcher due to the water in the flame trench. Based on the results established in our simulations of the SLS SMAT, we performed updated calculations for SLS at LC-39B with and without water systems active to assess the readiness of the launch pad for Artemis I launch. Our results show that the IOP/SS system is effective at reducing the overpressure signal and overall sound pressure levels felt by the vehicle and additionally that the pressure loads experienced by the mobile launcher (ML) during engine startup is not increased by the presence of water.

EGS↗
DOE OSTI2021

Quantum many-body simulations of the two-dimensional Fermi-Hubbard model in ultracold optical lattices

Understanding quantum many-body states of correlated electrons is one main theme in modern condensed-matter physics. Given that the Fermi-Hubbard model, the prototype of correlated electrons, was recently realized in ultracold optical lattices, it is highly desirable to have controlled numerical methodology to provide precise finite-temperature results upon doping to directly compare with experiments. Here, we demonstrate the exponential tensor renormalization group (XTRG) algorithm [Chen et al., Phys. Rev. X 8, 031082 (2018)], complemented by independent determinant quantum Monte Carlo, offers a powerful combination of tools for this purpose. XTRG provides full and accurate access to the density matrix and thus various spin and charge correlations, down to an unprecedented low temperature of a few percent of the tunneling energy. Finally, we observe excellent agreement with ultracold fermion measurements at both half filling and finite doping, including the sign-reversal behavior in spin correlations due to formation of magnetic polarons, and the attractive hole-doublon and repulsive hole-hole pairs that are responsible for the peculiar bunching and antibunching behaviors of the antimoments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Capacitance–voltage modeling of mid-wavelength infrared nBn detectors

Capacitance–voltage measurements are a powerful technique to determine doping profiles of semiconductor pn junctions and Schottky barrier diodes. The measurements were recently extended to III-V-based mid-wavelength nBn infrared detectors, and absorber doping densities have been extracted using the widely accepted Schottky approximation, where the potential drop across the device is assumed to be across the depleting absorber layer. However, this approach is limited to when the absorber region of the nBn is under high reverse bias and thus is only able to extract the absorber region doping profile. Here, we introduce a semi-analytical model that is capable of extracting barrier dopant polarity, doping concentration, and thickness, as well as contact and absorber layer doping concentrations, all from a capacitance–voltage measurement. Rather than solely considering the potential drop across the depleting layers, it considers the potential drop across the accumulating layer as well. This negative charge accumulation occurs for the contact and absorber layers in the case of reverse and forward biases, respectively. This allows for a single model to be applied to a capacitance–voltage curve at forward and reverse biases and it can provide regions of bias where the absorber transitions from depletion to accumulation. We compare the agreement of the semianalytical model with modeling results from commercially available finite element method software and experimental capacitance–voltage data. Finally, we show that the method is consistent with the Schottky approximation of extracting absorber doping densities at high reverse bias and discuss the model's limitations.

Physics↗
DOE OSTI2025

How are Heterogeneous Nucleation Rate Observations Influenced by Instrument Resolution?

Experimental measurements of the heterogeneous nucleation rate rely on counting the number of nuclei with time. However, the size of a thermodynamically stable nucleus is often a few nanometers in diameter and is below the resolution of most (in situ) measurement techniques that provide a statistically valid sample. Due to the finite resolution of the instruments and analysis methods, it is challenging to capture the incipient nuclei and the subsequent evolution of nuclei density over time. In this work, we demonstrate the impact of instrument resolution on observed nuclei densities by comparing numerical modeling with experimental results. Further, to achieve this, we implemented heterogeneous nucleation within the pore-scale reactive transport modeling framework using classical nucleation theory (CNT). We compared the modeling results with nucleation rates measured using X-ray nanotomography (XnT) and evaluated how these impact the apparent values of the prefactor and interfacial energy based on CNT and the crystal growth rate. Specifically, we applied a resolution threshold (artificial resolution limit) in the model during nuclei counting to resemble an experimental resolution, ranging from 15 to 500 nm. The findings reveal that the instrument resolution significantly impacts the apparent prefactor and interfacial energy. Both apparent prefactor and interfacial energy decrease with a decrease in the instrument resolution. While deviation in the prefactor due to resolution is anticipated, those in the interfacial energy are unexpected. The approach described here allows one to correct apparent nucleation rates that depend on the instrument’s resolution to derive “intrinsic” CNT parameters for the prefactor and interfacial energy.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Nanostructured $\mathrm{GaAs/(Al,Ga)As}$ Waveguide for Low-Density Polariton Condensation from a Bound State in the Continuum

Exciton-polaritons are hybrid light-matter states that arise from strong coupling between an exciton resonance and a photonic cavity mode. As bosonic excitations, they can undergo a phase transition to a condensed state that can emit coherent light without a population inversion. This aspect makes them good candidates for thresholdless lasers, yet short exciton-polariton lifetime has made it difficult to achieve condensation at very low power densities. In this sense, long-lived symmetry-protected states are excellent candidates to overcome the limitations that arise from the finite mirror reflectivity of monolithic microcavities. In this work we use a photonic symmetry-protected bound state in the continuum coupled to an excitonic resonance to achieve state-of-the-art polariton condensation threshold in a GaAs/(Al,Ga)As waveguide. Most important, we show the influence of fabrication control and how surface passivation via atomic layer deposition provides a way to reduce exciton quenching at the grating sidewalls.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Location and thermal evolution of the pseudogap due to spin fluctuations

We study pseudogap behavior in a metal near a spin density wave (SDW) instability due to thermal magnetic fluctuations. We consider the t–t' Hubbard model on a square lattice at a finite doping, at intermediate coupling strength, and analyze the thermal evolution of the electron spectral function between a SDW ordered state at low temperatures and a normal Fermi liquid at high temperatures. We argue that for proper description of the pseudogap one needs to sum up infinite series of diagrams for both the fermionic self-energy and the SDW order parameter in the SDW state or the magnetic correlation length in the paramagnetic state. We use the eikonal approach to sum up an infinite series of diagrammatic contributions from thermal fluctuations. Earlier studies found that in the SDW state the spectral function A k (ω) of a hot fermion at a finite T is exponentially small below the energy scale Δ(T), which scales with SDW order and vanishes at the ordering temperature T N , and has a hump at a larger frequency Δ PG , comparable to the zero-temperature SDW gap Δ(T = 0). We argue that the hump, which we associate with the pseudogap, survives in some T range above T N . We show that this range is split by regions of strong and weak pseudogap behavior. In the first region, Δ PG is weakly temperature dependent, despite the fact that it comes from thermal fluctuations. Such a behavior has been seen in numerical studies of the Hubbard model. We show that to obtain it one needs to go beyond the one-loop approximation and sum up the infinite series of diagrams. In the second regime, Δ PG decreases with increasing T and eventually vanishes. Here, we further argue that a magnetic pseudogap at a finite T emerges only if the ground state is magnetically ordered. We present the phase diagram and apply the results to high-T c cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Global description of $β$ – decay with the axially deformed Skyrme finite-amplitude method: Extension to odd-mass and odd-odd nuclei

Here we use the finite-amplitude method (FAM), an efficient implementation of the quasiparticle random phase approximation, to compute $β$-decay rates with Skyrme energy-density functionals for 3983 nuclei, essentially all the medium-mass and heavy isotopes on the neutron-rich side of stability. We employ an extension of the FAM that treats odd-mass and odd-odd nuclear ground states in the equal filling approximation. Our rates are in reasonable agreement both with experimental data where available and with rates from other global calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Multiscale Modeling of Nanoparticle Precipitation in Oxide Dispersion-Strengthened Steels Produced by Laser Powder Bed Fusion

Laser Powder Bed Fusion (LPBF) enables the efficient production of near-net-shape oxide dispersion-strengthened (ODS) alloys, which possess superior mechanical properties due to oxide nanoparticles (e.g., yttrium oxide, Y-O, and yttrium-titanium oxide, Y-Ti-O) embedded in the alloy matrix. To better understand the precipitation mechanisms of the oxide nanoparticles and predict their size distribution under LPBF conditions, we developed an innovative physics-based multiscale modeling strategy that incorporates multiple computational approaches. These include a finite volume method model (Flow3D) to analyze the temperature field and cooling rate of the melt pool during the LPBF process, a density functional theory model to calculate the binding energy of Y-O particles and the temperature-dependent diffusivities of Y and O in molten 316L stainless steel (SS), and a cluster dynamics model to evaluate the kinetic evolution and size distribution of Y-O nanoparticles in as-fabricated 316L SS ODS alloys. The model-predicted particle sizes exhibit good agreement with experimental measurements across various LPBF process parameters, i.e., laser power (110–220 W) and scanning speed (150–900 mm/s), demonstrating the reliability and predictive power of the modeling approach. The multiscale approach can be used to guide the future design of experimental process parameters to control oxide nanoparticle characteristics in LPBF-manufactured ODS alloys. Additionally, our approach introduces a novel strategy for understanding and modeling the thermodynamics and kinetics of precipitation in high-temperature systems, particularly molten alloys.

Wang, Zhengming (ORCID:0000000241627112)↗
NASA NTRS1986

An analysis of the role of drift waves in equatorial spread F

An account is given of results of rocket measurements of the wave number spectrum of equatorial spread F irregularities, with emphasis on wavelengths less than 100 m. The measurements were made from two sounding rockets launched from Peru as part of Project Condor. The Condor density fluctuation spectra display a break at a wavelength near 100 m, identical to that found in the PLUMEX experiment (Kelley et al., 1982). The Condor data also confirm a subrange in which the density and the wave potential obey the Boltzmann relation - a strong indication of the presence of low-frequency electrostatic waves with finite wavelength parallel to the magnetic field, perhaps low-frequency drift waves as proposed by Kelley et al. The Condor data are also consistent with the previous conjecture that drift waves only exist above 300 km altitude. To investigate the difference in spectra observed over two altitude ranges, the data must be fitted to a form for the power spectrum taken from Keskinen and Ossakow (1981). The fitted spectrum, along with empirically determined growth and dissipation rates, is used to calculate the energy pumped into the spectrum at long wavelengths as well as the energy dissipated at shorter wavelengths. It is found that the energy is balanced by classical collisional effects in the low-altitude case, but energy balance in the high-altitude case requires an enhanced dissipation of about 500 times that due to classical diffusion. The model is consistent with, but does not uniquely imply, an inverse cascade of drift wave turbulence in equatorial spread F.

Labelle, J.↗
NASA NTRS1990

Linear mixing rule in screened binary ionic mixtures

The validity of the linear mixing rule is examined for the following two cases (1) when the response of the electron gas is taken into account in the effective ionic interaction and (2) when finite-temperature effects are included in the dielectric response of the electrons, i.e., when the ions interact with both temperature- and density-dependent screened Coulomb potentials. It is found that the linear mixing rule remains valid when the electron response is taken into account in the interionic potential at any density, even though the departure from linearity can reach a few percent for the asymmetric mixtures in the region of weak degeneracy for the electron gas. A physical explanation of this behavior is proposed which is based on a simple additional length scale.

Chabrier, G.↗
NASA NTRS1990

Self-organized criticality in sandpiles - Nature of the critical phenomenon

The scaling behavior of sandpile models is investigated analytically. First, it is shown that sandpile models contain a set of domain walls, referred to as troughs, which bound regions that can experience avalanches. It is further shown that the dynamics of the troughs is governed by a simple set of rules involving birth, death, and coalescence events. A simple trough model is then introduced, and it is proved that the model has a phase transition with the density of the troughs as an order parameter and that, in the thermodynamic limit, the trough density goes to zero at the transition point. Finally, it is shown that the observed scaling behavior is a consequence of finite-size effects.

Carlson, J. M.↗
NASA NTRS1992

A six-degree-of-freedom guidance and control analysis of Mars aerocapture

A six-degree-of-freedom (6DOF) simulation is developed to investigate the control and guidance issues of a Mars aerobraking vehicle. The guidance algorithm used is a predictor-corrector guidance formulation designed to control the exit orbital apoapsis and wedge angle using bank-angle modulation. Major features of this predictor-corrector guidance algorithm include: (1) integration of the 3DOF equations of motion within an inner-loop simulation; (2) load-relief logic; (3) finite roll rates; and (4) an aerodynamic feedback multiplier. The algorithm is capable of successfully guiding the vehicle through combinations of atmospheric density dispersions, aerodynamic mispredictions, and off-nominal atmospheric interface conditions. This study demonstrated that the addition of vehicle dynamics to the Mars aerobraking simulation does not significantly impact mission feasibility. That is, a robust control system design coupled with an adaptive guidance algorithm can assure mission success in the presence of numerous off-nominal conditions.

Powell, Richard W.↗
NASA NTRS1992

Performance and state-space analyses of systems using Petri nets

The goal of any modeling methodology is to develop a mathematical description of a system that is accurate in its representation and also permits analysis of structural and/or performance properties. Inherently, trade-offs exist between the level detail in the model and the ease with which analysis can be performed. Petri nets (PN's), a highly graphical modeling methodology for Discrete Event Dynamic Systems, permit representation of shared resources, finite capacities, conflict, synchronization, concurrency, and timing between state changes. By restricting the state transition time delays to the family of exponential density functions, Markov chain analysis of performance problems is possible. One major drawback of PN's is the tendency for the state-space to grow rapidly (exponential complexity) compared to increases in the PN constructs. It is the state space, or the Markov chain obtained from it, that is needed in the solution of many problems. The theory of state-space size estimation for PN's is introduced. The problem of state-space size estimation is defined, its complexities are examined, and estimation algorithms are developed. Both top-down and bottom-up approaches are pursued, and the advantages and disadvantages of each are described. Additionally, the author's research in non-exponential transition modeling for PN's is discussed. An algorithm for approximating non-exponential transitions is developed. Since only basic PN constructs are used in the approximation, theory already developed for PN's remains applicable. Comparison to results from entropy theory show the transition performance is close to the theoretic optimum. Inclusion of non-exponential transition approximations improves performance results at the expense of increased state-space size. The state-space size estimation theory provides insight and algorithms for evaluating this trade-off.

Watson, James Francis, III↗
NASA NTRS2000

The Dynamics of Miscible Interfaces: A Space Flight Experiment

Experiments as well as accompanying simulations are described that serve in preparation of a space flight experiment to study the dynamics of miscible interfaces. The investigation specifically addresses the importance of both nonsolenoidal effects as well as nonconventional Korteweg stresses in flows that give rise to steep but finite concentration gradients. The investigation focuses on the flow in which a less viscous fluid displaces one of higher viscosity and different density within a narrow capillary tube. The fluids are miscible in all proportions. An intruding finger forms that occupies a fraction of the total tube diameter. Depending on the flow conditions, as expressed by the Peclet number, a dimensionless viscosity ratio, and a gravity parameter, this fraction can vary between approximately 0.9 and 0.2. For large Pe values, a quasi-steady finger forms, which persists for a time of O(Pe) before it starts to decay, and Poiseuille flow and Taylor dispersion are approached asymptotically. Depending on the specific flow conditions, we observe a variety of topologically different streamline patterns, among them some that leak fluid from the finger tip. For small Pe values, the flow decays from the start and asymptotically reaches Taylor dispersion after a time of O(Pe). Comparisons between experiments and numerical simulations based on the 'conventional' assumption of solenoidal velocity fields and without Korteweg stresses yield poor agreement as far as the Pe value is concerned that distinguishes these two regimes. As one possibility, we attribute this lack of agreement to the disregard of these terms. An attempt is made to use scaling arguments in order to evaluate the importance of the Korteweg stresses and of the assumption of solenoidality. While these effects should be strongest in absolute terms when steep concentration fronts exist, i.e., at large Pe, they may be relatively most important at lower values of Pe. We subsequently compare these conventional simulations to more complete simulations that account for nonvanishing divergence as well as Korteweg stresses. While the exact value of the relevant stress coefficients are not known, ballpark numbers do exist, and their use in the simulations indicates that these stresses may indeed be important. We plan to evaluate these issues in detail by means of comparing a space experiment with corresponding simulations, in order to extract more accurate Korteweg stress coefficients, and to confirm or deny the importance of such stresses.

Maxworthy, Tony↗
NASA NTRS2012

Compression After Impact on Honeycomb Core Sandwich Panels with Thin Facesheets, Part 2: Analysis

A two part research study has been completed on the topic of compression after impact (CAI) of thin facesheet honeycomb core sandwich panels. The research has focused on both experiments and analysis in an effort to establish and validate a new understanding of the damage tolerance of these materials. Part 2, the subject of the current paper, is focused on the analysis, which corresponds to the CAI testings described in Part 1. Of interest, are sandwich panels, with aerospace applications, which consist of very thin, woven S2-fiberglass (with MTM45-1 epoxy) facesheets adhered to a Nomex honeycomb core. Two sets of materials, which were identical with the exception of the density of the honeycomb core, were tested in Part 1. The results highlighted the need for analysis methods which taken into account multiple failure modes. A finite element model (FEM) is developed here, in Part 2. A commercial implementation of the Multicontinuum Failure Theory (MCT) for progressive failure analysis (PFA) in composite laminates, Helius:MCT, is included in this model. The inclusion of PFA in the present model provided a new, unique ability to account for multiple failure modes. In addition, significant impact damage detail is included in the model. A sensitivity study, used to assess the effect of each damage parameter on overall analysis results, is included in an appendix. Analysis results are compared to the experimental results for each of the 32 CAI sandwich panel specimens tested to failure. The failure of each specimen is predicted using the high-fidelity, physicsbased analysis model developed here, and the results highlight key improvements in the understanding of honeycomb core sandwich panel CAI failure. Finally, a parametric study highlights the strength benefits compared to mass penalty for various core densities.

Mcquigg, Thomas D.↗
DOE OSTI2023

Enhancement of charge density wave correlations in a Holstein model with an anharmonic phonon potential

The Holstein Hamiltonian describes itinerant electrons whose site density couples to local phonon degrees of freedom. In the single-site limit, at half filling, the electron-phonon coupling results in a double-well structure for the lattice displacement, favoring empty or doubly occupied sites. In two dimensions and on a bipartite lattice in d ≥ 2, intersite hopping causes these doubly occupied and empty sites to alternate in a charge density wave (CDW) pattern when the temperature is lowered. Because a discrete symmetry is broken, this occurs in a conventional second-order transition at finite T cdw . In this paper, we investigate the effect of changing the phonon potential energy to one with an intrinsic double-well structure even in the absence of electron-phonon coupling. While this aids in the initial process of pair formation, the implications for subsequent CDW order are nontrivial. One expects that, when the electron-phonon coupling is too strong, the double wells become deep and the polaron mass gets large, an effect which reduces T cdw . Here, we show here the existence of regions of parameter space where the double-well potential, while aiding local pair formation, does so in a way which also substantially enhances long-range CDW order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

On the mechanical behavior of a Nb 3 Sn superconducting coil during a quench: Two-dimensional finite element analysis of a quench heater protected magnet

New high-field accelerator magnets based on $Nb_3Sn$ superconductors are pushing the boundaries of magnet design and quench protection towards new limits. While their large stored energies and current densities result in a very challenging scenario for magnet protection, their great electromagnetic forces create also new requirements in terms of magnet design and stress management techniques. Furthermore, the strain sensitivity of $Nb_3Sn$ cables turns the electro-mechanical limits of the conductor into a parameter of the highest importance, where conductor degradation becomes a critical aspect in magnet operation. The coupling of all the above-mentioned considerations during quench is a case of special interest that adds further complexity to the design of $Nb_3Sn$ magnets. The objective of this paper is to provide a complete two-dimensional investigation of the coil and magnet structure mechanics during a quench event. Here, the analysis is performed using a combination of finite element codes that provide the necessary input for the mechanical study. The core of the modelling strategy relies on a thermal-electric model, whose results are directly used as loads in the mechanical simulation. In doing so, the stress evolution during and after a quench is obtained. We focus for this time in the analysis of a quench heater protected magnet, where dynamic effects are of less importance in contrast to other protection systems like CLIQ.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Functional Renormalization Group Study of Superconductivity in Rhombohedral Trilayer Graphene

We employ a functional renormalization group approach to ascertain the pairing mechanism and symmetry of the superconducting phase observed in rhombohedral trilayer graphene. Superconductivity in this system occurs in a regime of carrier density and displacement field with a weakly distorted annular Fermi sea. We find that repulsive Coulomb interactions can induce electron pairing on the Fermi surface by taking advantage of momentum-space structure associated with the finite width of the Fermi sea annulus. Furthermore, the degeneracy between spin-singlet and spin-triplet pairing is lifted by valley-exchange interactions that strengthen under the RG flow and develop nontrivial momentum-space structure. We find that the leading pairing instability is d-wave-like and spin-singlet, and that the theoretical phase diagram versus carrier density and displacement field agrees qualitatively with experiment.

36 MATERIALS SCIENCE↗