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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Impact mechanics of the Cretaceous-Tertiary extinction bolide

An examination of the mechanics of asteroidal, cometary, and meteor swarm impact on the earth determined if the enrichment of projectile material in the K-T layer is consistent with melts and impact breccias on the earth and moon, the size of the impacters, the distribution of the kinetic energy, and the sequence of impacts that could give rise to observed extinction phenomena. Flows resulting from spherical projectile impacts onto layers of air, water, and silicates were modeled and Eulerian finite difference algorithms were employed to solve conservation equations and equations of state. A range of speeds and impacter densities were considered, along with sizes from 0.17 km, which would be consumed in the atmosphere, to a 10 km object, which would have had a diameter greater than a reference 7.1 km atmosphere depth. It is concluded that an impact of the K-T bolide could result in global biotic extinction and worldwide material deposition.

Okeefe, J. D.↗
NASA NTRS2001

Impact and Collisional Processes in the Solar System

A series of impact experiments on anhydrite CaSO4, in which vaporized sample accelerates an element in a velocity interferometer, generate velocity data that we have recently reanalyzed using an explicit entropy generating finite difference code. The shock pressure required from the onset, and complete vaporization of 30% porous and 70% crystal density anhydrite is 52 +/- 3 and 122 +/- 13 GPa. Using observed acid leaching in non-marine K/T ejecta in North America, and the sharp global increase in Sr-87/Sr-86 ratios recorded at 65 Ma in marine rocks, we demonstrated that global acidification is primarily due to the SO2 released by anhydrite volatilization, and not HNO3 formed from bolide-induced air pyrolysis. Shock temperatures for crystal CaCO3 are measured from 3000 to 7000 K in the 90 to 160 GPa pressure range. These temperatures are much lower than calculated theoretically indicating that possibly bond breakdown at the shock front is occurring. This is the first mineral in which this effect has ever been seen. New data defining the ion species which are produced upon impact of volatilization of metals and minerals using a pulse ultraviolet laser to simulate intense shock heating from a projectile impact indicate that in shock experiments we can for the first time study the speciation of neutrals using a moderate resolution time-of-flight mass spectrometer. Measurements of the gas species from a series of proposed impact experiments appear to be quite feasible. We will attempt these experiments in the next year. Measurements of the impact induced shock wave decay in SiO2 and GeO2 glass are underway to measure these pure oxide properties. Predictive calculations indicate that the pressure will decay as r(exp -2.7) in the phase transition regime, versus a much lower rate of r(exp -1.18), if a phase transition does not occur.

Ahrens, Thomas J.↗
NASA NTRS1999

Transverse Mode Dynamics of VCSELs Through Space-Time Domain Simulation

Modeling and simulation are important to understand laser operation and to optimize and design device functions. Numerical simulation of VCSEL (Vertical Cavity Surface Emitting Lasers) has been largely based on solving time-independent Helmholtz equation or time dependent coupled mode equations. There are various advantages for choosing these approaches. However, the disadvantages are also apparent. The former cannot handle dynamical mode competition seen in VCSELs, while the latter assumes a given type and number of modes a priori. Furthermore, the microscopic physics of heterstructures and electron-hole plasma is very often represented by a few parameters such as linear gain coefficients and the linewidth enhancement factor. These are over simplification of space and frequency (wavelength) dependent gain and refractive index functions. When the space-time dynamical operation of VCSELs becomes important, these simple approximations become questionable. In this paper, we apply a recently developed model for edge-emitting lasers to a gain guided VCSEL for space-time domain simulation. This model takes into account the actual nonlinear dependence of gain and refractive index on frequency and carrier density within the frame work of the effective Bloch equations. The corresponding partial differential equations are solved directly by finite difference methods. Laser behavior with increasing pumping current is investigated in detail. Special attention is paid to the dynamical competition of the transverse modes.

Goorjian, Peter M.↗
DOE OSTI2024

Thermal quasiparticle theory

The widely used thermal Hartree–Fock (HF) theory is generalized to include the effect of electron correlation while maintaining its quasi-independent-particle framework. An electron-correlated internal energy (or grand potential) is postulated in consultation with the second-order finite-temperature many-body perturbation theory (MBPT), which then dictates the corresponding thermal orbital (quasiparticle) energies in such a way that all fundamental thermodynamic relations are obeyed. The associated density matrix is of a one-electron type, whose diagonal elements take the form of the Fermi–Dirac distribution functions, when the grand potential is minimized. The formulas for the entropy and chemical potential are unchanged from those of Fermi–Dirac or thermal HF theory. The theory thus stipulates a finite-temperature extension of the second-order Dyson self-energy of one-particle many-body Green’s function theory and can be viewed as a second-order, diagonal, frequency-independent, thermal inverse Dyson equation. At low temperatures, the theory approaches finite-temperature MBPT of the same order, but it may outperform the latter at intermediate temperatures by including additional electron-correlation effects through orbital energies. Here, a physical meaning of these thermal orbital energies is proposed (encompassing that of thermal HF orbital energies, which has been elusive) as a finite-temperature version of Janak’s theorem.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRSDate not supplied

Finite Element Modeling for Compression Strength Prediction of Honeycomb Cores with Geometric Imperfections Measured using X-ray CT Imaging

Aluminum honeycomb cores have been used extensively in composite sandwich panels due to having high bending rigidity while maintaining low density. During manufacturing of the aluminum honeycomb, the thin metal walls are susceptible to imperfections that deviate from an ideal honeycomb shape. Computational tools to quantify imperfections and investigate their effects on quasi-static compression and impact response can inform design criteria to construct safer and lighter launch vehicle structures. This paper presents the results of a study on finite element modeling of metallic honeycomb cores (HCC) with geometric imperfections measured using X-ray Computer Tomography (CT), to predict the compression response. Finite element (FE) models are constructed with measured imperfections with appropriate boundary conditions. It is shown that the average behavior of larger 10x11 cell models can be predicted by sampling single cell models from the larger domain. Amplification of imperfections for cells at the center of sampled specimens are mode switched from the larger out-of-plane imperfection to the mode shape of the adjacent cells resulting in strengthening of the center specimens.

aluminum honeycomb↗
NASA NTRS2013

A Patch Density Recommendation based on Convergence Studies for Vehicle Panel Vibration Response resulting from Excitation by a Diffuse Acoustic Field

Producing fluid structural interaction estimates of panel vibration from an applied pressure field excitation are quite dependent on the spatial correlation of the pressure field. There is a danger of either over estimating a low frequency response or under predicting broad band panel response in the more modally dense bands if the pressure field spatial correlation is not accounted for adequately. It is a useful practice to simulate the spatial correlation of the applied pressure field over a 2d surface using a matrix of small patch area regions on a finite element model (FEM). Use of a fitted function for the spatial correlation between patch centers can result in an error if the choice of patch density is not fine enough to represent the more continuous spatial correlation function throughout the intended frequency range of interest. Several patch density assumptions to approximate the fitted spatial correlation function are first evaluated using both qualitative and quantitative illustrations. The actual response of a typical vehicle panel system FEM is then examined in a convergence study where the patch density assumptions are varied over the same model. The convergence study results illustrate the impacts possible from a poor choice of patch density on the analytical response estimate. The fitted correlation function used in this study represents a diffuse acoustic field (DAF) excitation of the panel to produce vibration response.

Smith, Andrew↗
DOE OSTI2020

Data-Driven Approach for Uncertainty Propagation and Reachability Analysis in Dynamical Systems

In this paper, we propose a data-driven approach for uncertainty propagation and reachability analysis in a dynamical system. The proposed approach relies on the linear lifting of a nonlinear system using linear Perron-Frobenius (P-F) and Koopman operators. Here, the uncertainty can be characterized in terms of the moments of a probability density function. We demonstrate how the P-F and Koopman operators are used for propagating the moments. Time-series data is used for the finite-dimensional approximation of the linear operators, thereby enabling data-driven approach for moment propagation. Simulation results are presented to demonstrate the effectiveness of the proposed method.

42 ENGINEERING↗
DOE OSTI2021

A REDUCED ORDER MODELING APPROACH TO PROBABILISTIC CREEP-DAMAGE PREDICTIONS IN FINITE ELEMENT ANALYSIS

This paper introduces a computationally efficient Reduced Order Modeling (ROM) approach for the probabilistic prediction of creep-damage failure. Component-level probabilistic simulations are needed to assess the reliability and safety of high-temperature components. Full-scale probabilistic creep-damage modeling in finite element (FE) approach is computationally expensive requiring many hundreds of simulations to replicate the uncertainty of component failure. To that end, ROM is proposed to minimize the elevated computational cost while controlling the loss of accuracy. It is proposed that full-scale probabilistic simulations can be completed in 1D at a reduced cost, the extremum conditions extracted, and those conditions applied for lower-cost 2D/3D probabilistic simulations of components that capture the mean and uncertainty of failure. The probabilistic Sine-hyperbolic (Sinh) model is selected which in previous work was calibrated to alloy 304 stainless steel. The Sinh model includes probability density functions (pdfs) for test condition (stress and temperature), initial damage (i.e. microstructure), and material properties uncertainty. The Sinh model is programmed into ANSYS finite element software using the USERCREEP.F material subroutine. First, the Sinh model and FE code are subject to verification and validation to ensure the accuracy of the simulations. Numerous Monte Carlo simulations are executed in a 1D model to generate probabilistic creep deformation, damage, and rupture data. This data is analyzed and the probabilistic parameters corresponding to extreme creep response are extracted. The ROM concept is applied where only the extreme conditions are applied in the 2D probabilistic prediction of a component. The probabilistic predictions between the 1D and 2D geometry is compared to assess ROM for creep. The accuracy of the probabilistic prediction employing the ROM approach will potentially reduce the time and cost of simulating complex engineering systems. Future studies will introduce multi-stage Sinh, stochasticity, and spatial uncertainty for improved prediction.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mechanism for fluctuating pair density wave

Abstract In weakly coupled BCS superconductors, only electrons within a tiny energy window around the Fermi energy, E F , form Cooper pairs. This may not be the case in strong coupling superconductors such as cuprates, FeSe, SrTiO 3 or cold atom condensates where the pairing scale, E B , becomes comparable or even larger than E F . In cuprates, for example, a plausible candidate for the pseudogap state at low doping is a fluctuating pair density wave, but no microscopic model has yet been found which supports such a state. In this work, we write an analytically solvable model to examine pairing phases in the strongly coupled regime and in the presence of anisotropic interactions. Already for moderate coupling we find an unusual finite temperature phase, below an instability temperature T i , where local pair correlations have non-zero center-of-mass momentum but lack long-range order. At low temperature, this fluctuating pair density wave can condense either to a uniform d -wave superconductor or the widely postulated pair-density wave phase depending on the interaction strength. Our minimal model offers a unified framework to understand the emergence of both fluctuating and long range pair density waves in realistic systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1996

Assessment of thermochemical nonequilibrium and slip effects for Orbital Reentry Experiment (OREX)

Results are provided from a viscous shock layer (VSL) analysis of the reentry flowfield around the forebody of the Japanese Orbital Reentry Experiment (OREX) vehicle. This vehicle is a 50 deg. spherically blunted cone with a nose radius of 1.35 m and a base diameter of 3.4 m. Calculations are done for the OREX trajectory from 105 to 48.4 km altitude range. A 7-species chemical model is found adequate for the flowfield analysis. However, for altitudes greater than 84 km, the low density effects (such as thermal nonequilibrium and slip) are to be implemented for good agreement between the predictions and flight inferred heat-transfer rate data. Further, at altitudes lower than 84 km, a finite surface recombination probability is to be employed in place of a non-catalytic surface for better comparison between the calculations and data. VSL results are also compared with the direct simulation Monte Carlo (DSMC) predictions at high altitudes (greater than 80 km) and the electron number density data for three altitudes in the OREX trajectory. Overall, there is a good comparison between the flight data and calculated results. With the ongoing refinements in data extraction procedures, the OREX data should prove valuable for validating theoretical models employed in flowfield codes for calculation of reacting-gas flowfields.

Gupta, Roop N.↗
DOE OSTI2020

Quantifying leakage fields at ionic grain boundaries using off-axis electron holography

The electrical properties of interfaces in semiconductors and ionic conductors are immensely important in a wide range of applications. Electron holography is ideally suited for the direct measurement of the electrostatic potential of such interfaces. A key challenge with this approach is the contribution of the leakage field from the sample to the observed electron phase shift. This leakage field cannot be a priori independently determined and can cause an overestimation of the phase shift. In this work, we use finite element simulations to compute the three-dimensional electrostatic potential in the vicinity of an interface associated with a given interfacial charge density distribution. We then evaluate the predicted phase shift and demonstrate that the leakage field strongly affects the recovery of the projected interface potential. From the difference between the true potential and uncorrected, recovered potential, we propose a method to correct for this effect. We then demonstrate the application of this methodology to the analysis of experimental off-axis electron holography data acquired from the grain boundaries in lightly doped ceria.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Locally Engineering and Interrogating the Photoelectrochemical Behavior of Defects in Transition Metal Dichalcogenides

Transition metal dichalcogenides (TMDs) are attractive materials for a variety of applications in solar energy conversion and electrocatalysis, due to their favorable optical and electrical properties and their unique two-dimensional structures which facilitate the fabrication of wide-area, ultrathin layers. Unfortunately, the basal planes which make up the majority of these materials are relatively inert, and thus a great deal of effort has been directed to engineering favorable, catalytically active defects into these materials. In this paper, we demonstrate how probe-based electrochemical techniques can be employed as multifunctional tools for locally modifying TMD materials and probing the electrochemical behavior of the resulting defects. Scanning Electrochemical Cell Microscopy (SECCM) was employed to locally anodize exfoliated p-type WSe 2 nanosheets, creating hole-like defects within individual basal planes in a highly controllable fashion. Photoelectrochemical SECCM imaging was then employed to characterize the chemical behavior of these engineered defects, revealing significantly enhanced activity toward the Hydrogen Evolution Reaction (HER). Atomic force microscopy studies are presented which suggest these enhancements result from an increased density of monolayer-high step features within the anodized defects. Analysis of the SECCM data in the context of finite element simulations revealed that these enhancements increased with increasing anodization time, with local kinetic rates over 2 orders of magnitude higher than unaltered basal planes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Advection algorithms for quantum neutrino moment transport

Neutrino transport in compact objects is an inherently challenging multidimensional problem. Here, this difficulty is compounded if one includes flavor transformation—an intrinsically quantum phenomenon requiring one to follow the coherence between flavors and thus necessitating the introduction of complex numbers. To reduce the computational burden, simulations of compact objects that include neutrino transport often make use of momentum-angle-integrated moments (the lowest order ones being commonly referred to as the energy density and flux) and these quantities can be generalized to include neutrino flavor, i.e., they become quantum moments. Numerous finite-volume approaches to solving the moment evolution equations for classical neutrino transport have been developed based on solving a Riemann problem at cell interfaces. In this paper we describe our generalization of a Riemann solver for quantum moments, specifically decomposing complex numbers in terms of a (signed) magnitude and phase instead of real and imaginary parts. We then test our new algorithm in numerous cases showing a neutrino fast flavor instability, varying from toy models with analytic solutions to snapshots from neutron star merger simulations. Compared to previous algorithms for neutrino transport with flavor mixing, we find uniformly smaller growth rates of the flavor transformation along with concomitantly larger length-scales, and that the results are a better match with the growth rates seen from multiangle codes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Thermal Design and Optimization of High- Power Wireless Charging System

In this paper, the thermal design and optimization of a high-power wireless charging system (WCS) is proposed. An integrated electromagnetic and thermal co-optimization is essential to design a high power-density WCS pad. This study presents the thermal analysis and the potential optimization scopes for a polyphase WCS pad. The coil and core causes most of the power losses of a WCS pad and causes thermal hotspot in the pad. In this paper, the thermally conductive epoxy is introduced to balance the effective volumetric loss density in the pad and mitigate the thermal hotspots in the coil and core. The proposed design is simulated through finite element analysis (FEA) and tested experimentally for a 50 kW three-phase WCS. Simulation and experimental results show that the conventional design has localized hotspots in the coil and core. The proposed thermal design mitigates the thermal hotspots without increasing the coil or core volume.

Mohammad, Mostak↗
DOE OSTI2023

Imaging the breaking of electrostatic dams in graphene for ballistic and viscous fluids

The charge carriers in a material can, under special circumstances, behave as a viscous fluid. Here, in this work, we investigated such behavior by using scanning tunneling potentiometry to probe the nanometer-scale flow of electron fluids in graphene as they pass through channels defined by smooth and tunable in-plane p-n junction barriers. We observed that as the sample temperature and channel widths are increased, the electron fluid flow undergoes a Knudsen-to-Gurzhi transition from the ballistic to the viscous regime characterized by a channel conductance that exceeds the ballistic limit, as well as suppressed charge accumulation against the barriers. Our results are well modeled by finite element simulations of two-dimensional viscous current flow, and they illustrate how Fermi liquid flow evolves with carrier density, channel width, and temperature.

Science & Technology - Other Topics↗
NASA NTRS1987

Nonlinear susceptibilities of finite conjugated organic polymers

Tight-binding calculations of the length dependence of the third-order molecular hyperpolarizability for polyenes and polyynes are reported. The pi-electron wave functions were determined by exploiting the limited translational symmetry of the molecules. Perturbation theory was used to calculate the longitudinal component of the electronic nonresonant hyperpolarizability. This is the first two-'band' calculation of third-order hyperpolarizabilities on finite pi-electron systems of varying length. In contrast to the results of the one-'band' models, the hyperpolarizability densities increase rapidly and then, after about 10-15 repeating units, approach an asymptotic value.

Beratan, David N.↗
NASA NTRS1990

The nonlinear evolution of magnetized solar filaments

Thermal instability driven by optically thin radiation is believed to initiate the formation of plasma filaments in the solar corona. The fact that filaments are observed generally to separate regions of opposite, line-of-sight, magnetic polarity in the underlying photosphere suggests that filament formation requires the presence of a highly sheared, local magnetic field. Two-dimensional, nonlinear, magnetohydrodynamic simulations of the local genesis and growth of solar filaments in a force-free, sheared, magnetic field were performed, and the evolution of generic perturbations possessing broad spatial profiles was traced. It was found that simulations of the evolution of initial random-noise perturbations produce filamentary plasma structures that exhibit densities and temperatures characteristic of observed solar filaments. Furthermore, in each of these simulations, the filament axis lies at a finite angle with respect to the local magnetic field, consistent with solar observations.

Sparks, L.↗
NASA NTRS1990

The pdf approach to turbulent flow

This paper provides a detailed discussion of the theory and application of probability density function (pdf) methods, which provide a complete statistical description of turbulent flow fields at a single point or a finite number of points. The basic laws governing the flow of Newtonian fluids are set up in the Eulerian and the Lagrangian frame, and the exact and linear equations for the characteristic functionals in those frames are discussed. Pdf equations in both frames are derived as Fourier transforms of the equations of the characteristic functions. Possible formulations for the nonclosed terms in the pdf equation are discussed, their properties are assessed, and closure modes for the molecular-transport and the fluctuating pressure-gradient terms are reviewed. The application of pdf methods to turbulent combustion flows, supersonic flows, and the interaction of turbulence with shock waves is discussed.

Kollmann, W.↗