SEARCH · Engineering Papers
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.
An entry from the Cambridge Structural Database, the world’s repository for small molecule crystal structures. The entry contains experimental data from a crystal diffraction study. The deposited dataset for this entry is freely available from the CCDC and typically includes 3D coordinates, cell parameters, space group, experimental conditions and quality measures.
An entry from the Cambridge Structural Database, the world’s repository for small molecule crystal structures. The entry contains experimental data from a crystal diffraction study. The deposited dataset for this entry is freely available from the CCDC and typically includes 3D coordinates, cell parameters, space group, experimental conditions and quality measures.
The transport of molecules in microporous material is a significant and active area of research in separation applications. The movement of vapor/liquid molecules in a high-loading condition in microporous spaces is especially challenging to interrogate. This project aims to understand and control the transport of complex water and organic solvent mixtures in varying structures of microporous carbon molecular sieve (CMS) and activated carbon membranes. Modeling of transport behavior in such carbonaceous samples using computational modeling is difficult due to their amorphous structure. Therefore, this research intends to study the transport mechanism of water-organic mixtures by experimentally revealing fundamental transport properties, such as guest sorption amounts as well as diffusion and permeation rates. Considering the different size and guest-host affinity of water and organic solvent molecules, the objective of this research is to understand the structural conditions within the ultramicropores and micropores that generate different types of molecular transport and selection mechanisms within such complex systems. Carbonaceous materials are developed using tailored pyrolysis techniques to pyrolyze polymeric precursors, including polyvinylidene chloride (PVDC),polyvinylidene fluoride (PVDF), polymer of intrinsic microporosity (PIM) -1, and fully aromatic polyamide. The initial stage of the research will focus on the microscopic diffusion and sorption studies of pure component water, xylene isomers, and n,n- dimethylformamide within various carbons. Additional structural investigation on CMS, such as gas physisorption and neutron scattering studies, will be executed to gain deeper insight into these structure-transport relationships. Scanning electron microscopy and X-ray photoelectron microscopy are also used to further characterize these materials. The research will progress to probe mutual diffusion of complex water-organic mixtures within various CMS microstructures. Competitive sorption and permeation studies of the binary mixtures in different microstructural CMS will be followed. The obtained binary mixture transport parameters were used to predict binary water-organic mixtures transport.
Perovskite photovoltaics are attractive for both terrestrial and space applications. Although terrestrial conditions require durability against stressors such as moisture and partial shading, space poses different challenges: radiation, atomic oxygen, vacuum and high-temperature operation. Here we demonstrate a silicon oxide layer that hardens perovskite photovoltaics to critical space stressors. A 1-um-thick silicon oxide layer evaporated atop the device contacts blocks 0.05 MeV protons at fluences of 10 15 cm -2 without a loss in power conversion efficiency, which results in a device lifetime increase in low Earth orbit by x20 and in highly elliptical orbit by x30. Silicon-oxide-protected Cs 0.05 (MA 0.17 FA 0.83 ) 0.95 Pb(I 0.83 Br 0.17 ) 3 ) (MA, methylammonium; FA, formamidinium cation) and CsPbI 2 Br cells survive submergence in water and N,N-dimethylformamide. Furthermore, moisture tolerance of Sn-Pb and CsPbI 2 Br devices is boosted. Devices are also found to retain power conversion efficiencies on exposure to alpha irradiation and atomic oxygen. So, this barrier technology is a step towards lightweight packaging designs for both space and terrestrial applications.
The thickness uniformity of chemical vapor deposited (CVD) diamond coatings on non-planar, non-stationary substrates depends on both the intrinsic instantaneous coating thickness distribution (ICTD) of the coating conditions used and, if applicable, on the frequency of substrate reorientation. While important for many CVD diamond applications, the relative impact of the ICTD and substrate reorientation on the coating thickness uniformity has not been studied. In this work, we systematically investigate the effect of these factors for microwave-plasma chemical vapor deposition (MPCVD) of diamond (referred to as high density carbon (HDC) in the inertial confinement fusion (ICF) community) coatings on spherical, rolling substrates. This coating technique is used to fabricate capsules for ICF experiments, which require extreme coating uniformity with <0.3 % thickness variation (so-called Mode 1 or M1) to ensure symmetric compression of imploding targets. To extract the otherwise unobservable reorientation timescale (Δt), Monte Carlo simulations were performed using experimental ICTD data as input. This combined approach confirms scaling relationships between the substrate reorientation timescale as well as coating thickness and coating uniformity, as expected from a 3D random walk. Simulations confirm that M1 is Rayleigh-distributed and scales as (Δt) 1/2 , consistent with the randomization of two angles that determine orientation of a sphere. We also demonstrate that, under the conditions studied, Δt is the dominant factor in determining thickness uniformity while the intrinsic ICTD has minimal impact. Finally, experiments show that Δt can be affected by total batch size under constant agitation conditions due to space constraints that limit the capsule reorientation kinetics. In conclusion, this study highlights the utility of a combined experiment-simulation approach as a general methodology for understanding and improving coating uniformity on non-planar, non-stationary substrates.
Significant progress has been made in designing magnetic fields that provide excellent confinement of the guiding-centre trajectories of alpha particles using quasisymmetry (QS). Given the reduction in this transport channel, we assess the impact of resonant Alfvén eigenmodes (AEs) on the guiding-centre motion. The AE amplitudes are chosen to be consistent with experimental measurements and large-scale simulations. We evaluate the drift resonance condition, phase-space island width and island overlap criterion for quasisymmetric configurations. Kinetic Poincaré plots elucidate features of the transport, including stiff transport above a critical perturbation amplitude. Our analysis highlights key departures from the AE-driven transport in tokamaks, such as the avoidance of phase-space island overlap in quasihelical configurations and the enhanced transport due to wide phase-space islands in low magnetic shear configurations. In configurations that are closer to QS, with QS deviations $\delta B/B_0 \lesssim 10^{-3}$ , the transport is primarily driven by the AE, while configurations that are further from QS, $\delta B/B_0 \sim 10^{-2}$ , experience significant transport due to the QS-breaking fields in addition to the AE.
Forward osmosis (FO) may serve as a near-zero energy approach for wastewater management and selective chemical enrichment with an appropriate FO membrane. There is a need for novel FO membranes with better physical stability, chemical compatibility, and enhanced reverse flux selectivity. In this work, we have systematically studied a variety of factors for fabricating thermally reduced graphene oxide (RGO) FO filtration membranes and have demonstrated that with a comparable water flux, the RGO membranes have an exceptionally large reverse flux selectivity with 1.5 M sodium sulfate draw solution, about 7 times larger than that of cellulose triacetate membranes under the same testing condition. The interlayer spacing of the membranes can be fine-tuned by varying temperatures, so that the free interlayer spacing is less than 0.7 nm after exposure to water. The RGO-based membranes offer a few advantages with high reverse flux selectivity, mechanical robustness with strong adhesion to nylong support membranes without damage after FO tests and tape peel tests, and enhanced chromate and chlorine resistance. Finally, we have further demonstrated that in the FO processes, one can concentrate waste brine while simultaneously harvesting the osmotic energy as electricity, offering potential applications of forward osmosis-based osmotic batteries for powering electronic devices.
Abstract Two-dimensional carbides and nitrides, known as MXenes, are promising for water-processable coatings due to their excellent electrical, thermal, and optical properties. However, depositing hydrophilic MXene nanosheets onto inert or hydrophobic polymer surfaces requires plasma treatment or chemical modification. This study demonstrates a universal salt-assisted assembly method that produces ultra-thin, uniform MXene coatings with exceptional mechanical stability and washability on various polymers, including high-performance polymers for extreme temperatures. The salt in the Ti 3 C 2 T x colloidal suspension reduces surface charges, enabling electrostatically hydrophobized MXene deposition on polymers. A library of salts was used to optimize assembly kinetics and coating morphology. A 170 nm MXene coating can reduce radiation temperature by ~200 °C on a 300 °C PEEK substrate, while the coating on Kevlar fabric provides comfort in extreme conditions, including outer space and polar regions.
In magnetic confinement fusion devices, the ratio of the plasma pressure to the magnetic field energy, β, can become sufficiently large that electromagnetic microinstabilities become unstable, driving turbulence that distorts or reconnects the equilibrium magnetic field. In this paper, a theory is proposed for electromagnetic, electron-driven linear instabilities that have current layers localised to mode-rational surfaces and binormal wavelengths comparable to the ion gyroradius. The model retains axisymmetric toroidal geometry with arbitrary shaping, and consists of orbit-averaged equations for the mode-rational surface layer, with a ballooning space kinetic matching condition for passing electrons. The matching condition connects the current layer to the large scale electromagnetic fluctuations, and is derived in the limit that β is comparable to the square root of the electron-to-ion-mass ratio. Electromagnetic fluctuations only enter through the matching condition, allowing for the identification of an effective β that includes the effects of equilibrium flux surface shaping. The scaling predictions made by the asymptotic theory are tested with comparisons to results from linear simulations of micro-tearing and electrostatic microinstabilities in MAST discharge #6252, showing excellent agreement. In particular, it is demonstrated that the effective β can explain the dependence of the local micro-tearing mode (MTM) growth rate on the ballooning parameter θ 0 – possibly providing a route to optimise local flux surfaces for reduced MTM-driven transport.
Here we report neutron diffraction results for the body-centered-tetragonal series Eu(Co 1-x Ni x ) 2-y As 2 , x=0.10, 0.20, 0.42, and 0.82, y≤0.10, that detail changes to the magnetic ordering with nominal hole doping. We report the antiferromagnetic (AFM) propagation vectors, magnetic transition temperatures, and the ordered magnetic moments. We find a nonmonotonic change of the AFM propagation vector with x, with a minimum occurring at the tetragonal to collapsed-tetragonal phase crossover. For x=0.10 and 0.82 we find c-axis helix ordering of the Eu magnetic moments (spins) similar to x=0 and 1, with the spins oriented within the ab-plane. For x=0.20 and 0.42 we find higher-temperature c-axis ferromagnetic (FM) order and lower-temperature c-axis cone order. Using the extinction conditions for the space group, we discovered that the Eu spins are ordered in the higher-temperature c-axis FM phase for intermediate values of x, contrary to a previous report suggesting only Co/Ni spin ordering. Although we cannot directly confirm that the Co/Ni spins are also ordered, we suggest that c-axis itinerant-FM ordering of the Co/Ni spins could provide a molecular field that drives FM ordering of the Eu spins, which in turn provides the anisotropy for the lower-temperature c-axis cone order.
Hamiltonian formulation of lattice gauge theories (LGTs) is the most natural framework for the purpose of quantum simulation, an area of research that is growing with advances in quantum-computing algorithms and hardware. It, therefore, remains an important task to identify the most accurate, while computationally economic, Hamiltonian formulation(s) in such theories, considering the necessary truncation imposed on the Hilbert space of gauge bosons with any finite computing resources. This paper is a first step toward addressing this question in the case of non-Abelian LGTs, which further require the imposition of non-Abelian Gauss’s laws on the Hilbert space, introducing additional computational complexity. Focusing on the case of SU(2) LGT in 1+1 dimensions coupled to one flavor of fermionic matter, a number of different formulations of the original Kogut-Susskind framework are analyzed with regard to the dependence of the dimension of the physical Hilbert space on boundary conditions, system’s size, and the cutoff on the excitations of gauge bosons. The impact of such dependencies on the accuracy of the spectrum and dynamics obtained from a Hamiltonian computation is examined, and the (classical) computational-resource requirements given these considerations are studied. Besides the well-known angular-momentum formulation of the theory, the cases of purely fermionic and purely bosonic formulations (with open boundary conditions), and the Loop-String-Hadron formulation are analyzed, along with a brief discussion of a Quantum-Link-Model formulation of the same theory. Clear advantages are found in working with the Loop-String-Hadron framework which implements non-Abelian Gauss’s laws a priori using a complete set of gauge-invariant operators. Although small lattices are studied in the numerical analysis of this work, and only the simplest algorithms are considered, a range of conclusions will be applicable to larger systems and potentially to higher dimensions. Future studies will extend this investigation to the analysis of resource requirements for quantum simulations of non-Abelian LGT, with the goal of shedding light on the most efficient Hamiltonian formulation of gauge theories of relevance in nature
Increasing concern over higher frequency extreme weather events is driving a push towards a more resilient built environment. In recent years there has been growing interest in understanding how to evaluate, measure, and improve building energy resilience, i.e., the ability of a building to provide energy-related services in the event of a local or regional power outage. In addition to human health and safety, many stakeholders are keenly interested in the ability of a building to allow continuity of operations and minimize business disruption. Office buildings are subject to significant economic losses when building operations are disrupted due to a power outage. We propose “occupant hours lost” (OHL) as a means to measure the business productivity lost as the result of a power outage in office buildings. OHL is determined based on indoor conditions in each space for each hour during a power outage, and then aggregated spatially and temporally to determine the whole building OHL. We used quasi-Monte Carlo parametric energy simulations to demonstrate how the OHL metric varies due to different building characteristics across different climate zones and seasons. The simulation dataset was then used to develop simple regression models for assessing the impact of ten key building characteristics on OHL. The most impactful were window-to-wall ratio and window characteristics. The regression models show promise as a simple means to assess and screen for resilience using basic building characteristics, especially for non-critical facilities where it may not be viable to conduct detailed engineering analysis.
Energy efficient heat pump water heaters (HPWHs) reduce customer energy costs; grid-connected controls further increase HPWH value by enabling energy storage and shifting demand to cheaper, off-peak periods. However, a HPWH’s first cost premium can put it out of reach for low-income customers. This paper will present results of an ongoing study exploring the benefits of load shifting HPWHs for 24 low-income households, aiming to support increased product deployment in this demographic. The study uses EcoPort modules to shift HPWH load for a study sample consisting primarily of low-income adults and seniors in North Carolina. Building upon prior studies, the study’s HPWH operating schedules are designed to maximize shifted energy and minimize participant electricity costs based on local time-of-use rates. The study documents load profiles of grid-connected HPWHs in a less-studied demographic group and explores how certain groups may have particularly flexible loads due to their unique usage profiles. The results are relevant for HPWH product performance in the Southeast. The study will quantify HPWH load shifting performance during North Carolina’s hot summer, temperate shoulder, and occasionally sub-freezing winter seasons, and the study includes installations in conditioned, semi-conditioned, and unconditioned spaces. Finally, study results will also explore how controlled HPWHs can provide value to the regional grid by reducing seasonal peaks via demand response. Lessons learned include HPWH acceptance for low-income and senior users, best practices for maintaining HPWH connectivity among users without prior product experience, and strategies to leverage remote monitoring to ensure optimal operation and enhance HPWH reliability.
Over the last several decades, designers have used digital screens to view images of real and simulated spaces and make critical design decisions. Screen technology has improved during this time, as technologies like OLED have replaced legacy displays (CRT, plasma, and LCD). These new screens provide a higher pixel resolution, luminous output and contrast ratio. Immersive head-mounted displays now allow designers to view immersive images, and recent developments in real-time rendering have encouraged the uptake of virtual reality (VR) head-mounted displays in mainstream practice and design education. This paper presents an experiment on lighting perception using a series of LED lighting conditions in a real space and a virtual representation of those conditions captured using a 360° high-dynamic-range camera and presented on an HTC Vive Pro HMD. Fifty-three participants were asked to rate each lighting condition by viewing it in a real space (n = 30) or via immersive HDR photographs displayed in a VR HMD (n = 23). The results show that ratings of visual comfort, pleasantness, evenness, contrast and glare are similar between the HTC Vive Pro HMD and our real space when evaluating well-lit scenes, but significant differences emerge in dim and highly contrasted scenes for a number of rating scales.
We review some aspects of the cutting and gluing law in local quantum field theory (QFT) and study it from a new point of view. In particular, we emphasize the description of gluing by a path integral over a space of polarized boundary conditions, which are given by leaves of some Lagrangian foliation in the phase space. We think of this path integral as a non-local (d – 1)-dimensional gluing theory associated to the parent local d-dimensional QFT. This is a novel point of view paving the way for applications of the standard QFT techniques (that do not rely on locality) to the gluing theory. We describe various properties of this procedure and spell out conditions under which symmetries of the parent theory lead to symmetries of the gluing theory. The purpose of this paper is to set up a playground for the companion paper where these techniques are applied to obtain new results in supersymmetric theories.
A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.
Abstract Both experimental and computational methods for the exploration of structure, functionality, and properties of materials often necessitate the search across broad parameter spaces to discover optimal experimental conditions and regions of interest in the image space or parameter space of computational models. The direct grid search of the parameter space tends to be extremely time-consuming, leading to the development of strategies balancing exploration of unknown parameter spaces and exploitation towards required performance metrics. However, classical Bayesian optimization (BO) strategies based on the Gaussian process (GP) do not readily allow for the incorporation of the known physical behaviors or past knowledge. Here we explore a hybrid optimization/exploration algorithm created by augmenting the standard GP with a structured probabilistic model of the expected system’s behavior. This approach balances the flexibility of the non-parametric GP approach with a rigid structure of physical knowledge encoded into the parametric model. The fully Bayesian treatment of the latter allows additional control over the optimization via the selection of priors for the model parameters. The method is demonstrated for a noisy version of a standard univariate test function used to evaluate optimization algorithms and further extended to physical lattice models. This methodology is expected to be universally suitable for injecting prior knowledge in the form of physical models and past data in the BO framework.
A challenge to simulating turbulent flow in multiscale atmospheric applications is the efficient generation of resolved turbulence motions over an area of interest. One approach is to apply small perturbations to flow variables near the inflow planes of turbulence-resolving simulation domains nested within larger mesoscale domains. While this approach has been examined in numerous idealized and simple terrain cases, its efficacy in complex terrain environments has not yet been fully explored. Here, we examine the benefits of the stochastic cell perturbation method (CPM) over real complex terrain using data from the 2017 Perdigão field campaign, conducted in an approximately 2-km wide valley situated between two nearly parallel ridges. Following a typical configuration for multiscale simulation using nested domains within the Weather Research and Forecasting (WRF) model to downscale from the mesoscale to a large-eddy simulation (LES), we apply the CPM on a domain with horizontal grid spacing of 150 m. At this resolution, spurious coherent structures are often observed under unstable atmospheric conditions with moderate mean wind speeds. Results from such an intermediate resolution grid are often nested down for finer, more detailed LES, where these spurious structures adversely affect the development of turbulence on the subsequent finer grid nest. We therefore examine the impacts of the CPM on the representation of turbulence within the nested LES domain under moderate mean flow conditions in three different stability regimes: weakly convective, strongly convective, and weakly stable. In addition, two different resolutions of the underlying terrain are used to explore the role of the complex topography itself in generating turbulent structures. We demonstrate that the CPM improves the representation of turbulence within the LES domain, relative to the use of high-resolution complex terrain alone. During the convective conditions, the CPM improves the rate at which smaller-scales of turbulence form, while also accelerating the attenuation of the spurious numerically generated roll structures near the inflow boundary. During stable conditions, the coarse mesh spacing of the intermediate LES domain used herein was insufficient to maintain resolved turbulence using CPM as the flow develops downstream, highlighting the need for yet higher resolution under even weakly stable conditions, and the importance of accurate representation of flow on intermediate LES grids.