Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct coupling analysis”

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 883 records · Page 49

Uncovering hidden market opportunities for advanced nuclear reactors

Decarbonizing to meet aggressive climate change mitigation targets requires energy transition within all sectors. In the industrial sector, global emissions will need to decrease by 65–90% by 2050 to avert warming greater than 1.5°C (Pörtner et al., 2022). Recent U.S. laws, including the Inflation Reduction Act (IRA), Bipartisan Infrastructure Law (BIL), Defense Production Act, Creating Helpful Incentives to Produce Semiconductors (or CHIPS), state programs, and other recent laws, have clean energy requirements and provide financial incentives to accelerate the use of clean energy technologies in the industrial sector. These new laws include supporting mechanisms with direct financial support for nuclear power (e.g., advanced reactor development and hydrogen production). However, advanced nuclear could also gain these financial benefits by coupling with low-carbon industrial projects. For example, microreactors could supply low-carbon energy to producers of critical metals that (1) are eligible to receive investment and production tax credits and favorable loans, and (2) the low-carbon product could gain preference in emerging markets for green products.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Resolving Discrepancies in Spin-Torque Ferromagnetic Resonance Measurements: Lineshape versus Linewidth Analyses

When spin-orbit torques are measured using spin-torque ferromagnetic resonance, two alternative ways of analyzing the results to extract the torque efficiencies—lineshape analysis and analysis of the change in linewidth versus direct current—often give inconsistent results. We identify a source for these inconsistencies. We show that fits of spin-torque ferromagnetic resonance data to the standard analysis framework leave significant residuals that we identify as due to (i) current-induced excitations of a small volume of magnetic material with magnetic damping much larger than that of the bulk of the magnetic layer, which we speculate is associated with the heavy-metal–magnet interface and (ii) oscillations of the sample magnetization at the modulation frequency due to heating. Here, the dependence of the residual signals on direct current can interfere with an accurate extraction of spin-torque efficiencies by the linewidth method. We show that the discrepancies between the two types of analysis can be largely eliminated by extrapolation of the window of magnetic fields used in the linewidth fits to small values so as to minimize the influence of the residual signals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A formulation of rotor-airframe coupling for design analysis of vibrations of helicopter airframes

A linear formulation of rotor airframe coupling intended for vibration analysis in airframe structural design is presented. The airframe is represented by a finite element analysis model; the rotor is represented by a general set of linear differential equations with periodic coefficients; and the connections between the rotor and airframe are specified through general linear equations of constraint. Coupling equations are applied to the rotor and airframe equations to produce one set of linear differential equations governing vibrations of the combined rotor airframe system. These equations are solved by the harmonic balance method for the system steady state vibrations. A feature of the solution process is the representation of the airframe in terms of forced responses calculated at the rotor harmonics of interest. A method based on matrix partitioning is worked out for quick recalculations of vibrations in design studies when only relatively few airframe members are varied. All relations are presented in forms suitable for direct computer implementation.

Kvaternik, R. G.↗

The nearby 2-solar mass Bok globule LBN 11 - Sub-sonic molecular clumps in a magnetic environment

A multiwavelength study of the nearby small Bok globule LBN 11 is performed in order to investigate the relationship between embedded magnetic fields, cloud structure, and star formation. Optical polarimetry of background stars is used to characterize and trace the embedded magnetic field direction. High-dispersion CO isotopic mapping is employed to determine the radial density distribution of the gas, measure cloud rotation, and identify individual gas clumps. CS and SO mapping revealed dense cores within the (C-13)O clumps. Analysis of coadded IRAS images reveal the cloud to be free of current or recent star formation: there are no IR point sources of sufficient flux or proper colors to signify young stars associated with the cloud. The molecular line maps show the cloud to be extremely clumpy. Clumps are found to range in size from 0.2 pc for the largest CO clump to about 0.04 pc for the smallest CS feature. It is concluded that, at least in the envelope of the cloud, the magnetic field and the CO clumps are coupled.

Clemens, Dan P.↗

Time-Resolved OH-PLIF Assessment of Deflagrations Levels in a CH4-O2 Rotating Detonation Rocket Combustor

The parasitic loss incurred by deflagrative pre-burning is considered to be one of the key challenges to the effective implementation of rotating detonation engine (RDE) systems. Thus far, there have been relatively limited high fidelity, spatio-temporally resolved measurements of the pre-burning process as a function of the reactant conditions and composition within an RDE. In this work, simultaneous high-repetition-rate broadband OH* chemiluminescence and OH planar laser-induced fluorescence (PLIF) imaging are employed to investigate deflagrative burning dynamics in a fully optically accessible CH4-O2 RDE as compared to a H2-air systems under the same conditions. A custom-built optical parametric oscillator (OPO) is coupled with a high-repetition-rate burst-mode laser to generate the 284 nm source for the excitation of the Q1(9) transition of the OH radical. Significant deflagrative burning is observed throughout the chamber as a consequence of the oxygen-rich environment. Trailing wave systems that consume unburned reactants in the region immediately following the primary detonation wave are observed. The formation of a product gas recirculation zone that entrains and combusts incoming reactants is observed, and quantitative analysis is performed to gain valuable insight into deflagration characteristics. At an oxidizer mass flux of Gox ≈ 350kg/m2/s, it is found that the levels of pre-wave burning in CH4-O2 are a factor of ∼2-3 times higher than in H2-air in the region directly ahead of the detonation wave at a global equivalence ratio Φ ≈ 1.0. This study highlights the key roles of ignition delay and turbulent mixing of the combustion product and reactant fields on the levels of pre-burning and establishes a methodology for further investigation, such as for the effects of various RDE inlet configurations.

Propulsion↗

Acoustic properties, elasticity, and equation of state of glycerol under pressure

Herein we employed high-pressure Brillouin scattering to study the pressure dependencies of acoustic modes of glycerol up to 14 GPa at 300 K. We observed longitudinal acoustic velocities and transverse acoustic velocities for the first time from 5 to 14 GPa. The results allow the determination of a complete set of elastic properties and an accurate determination of the pressure–volume (P–V) equation of state (EOS). EOS parameters, K 0 = 14.9 ± 1.8 GPa and K' 0 = 5.6 ± 0.5, were determined from fits to the data from ambient pressure to 14 GPa. Direct volume measurements of the P–V EOS are consistent with those determined by Brillouin scattering. A deviation from a Cauchy-like relationship for elastic properties was observed, and the pressure dependencies of the photoelastic constants and relaxation times were documented from 5 to 14 GPa. These results have broad implications for glass-forming liquids, viscoelastic theory, and mode coupling theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gas-Phase Composition as a Predictive Metric for Calendar Life Behavior of Next-Generation Silicon Anodes

The expansion of renewable technologies and electrification of the transportation sector is driving increased demand for next-generation battery materials that provide higher power and energy density with superior cycling and calendar life stability. Silicon (Si) has a theoretical capacity nearly 10x that of graphite, and is therefore a promising anode material candidate to meet these rigorous performance demands. While leading Si anode battery demonstrations are approaching target metrics for cycle life, a series of complex and interrelated modes of reactivity lead to reduced calendar life and therefore challenge practical adoption of these materials. Deconvoluting the degradation processes that impact Si calendar life is critical to informing the rational and accelerated design of improved Si materials. In the present work, we employ novel sampling techniques and GC-MS-FID characterization to measure gas-phase composition during initial Si cycling, which we tie to selective mechanisms of Si passivation. We utilize a tiered analysis approach to identify and quantify the gas-phase reaction products associated with three advanced Si material candidates under practical operating conditions. Ex situ analysis of Si powders (pure chemical reactivity) is coupled with nondestructive in situ sampling of Si electrodes in a practical pouch-cell format (coupled chemical and electrochemical reactivity). We link the observed gas-phase species evolution to electrochemical behavior and measured calendar life of the three Si materials. Further, we evaluate the voltage-resolved evolution of gas-phase species for one such Si nanomaterial, where nonmonotonic gas generation implies competition between passivating reaction pathways. The measured gas-phase compositional data serves as a critical input for our advanced electrochemical SEI models to identify favorable vs unfavorable reaction pathways to stabilize Si. In addition to bolstering a fundamental understanding of Si reactivity, the present approach informs specific and quantifiable gas-phase metrics tied to calendar life improvements in Si, which can streamline and accelerate the process of next-generation material development.

DIRECT ENERGY CONVERSION,ENERGY STORAGE↗

Renewable Thermal Hybridization Framework for Industrial Process Heat Applications

Solar industrial process heat (SIPH) technologies, such as concentrating solar power collectors, could economically replace the steam or heat needs at many industrial sites by providing high-temperature heat transfer fluids (HTFs) such as pressurized water, synthetic-oil, or direct steam. Renewable thermal energy systems (RTES) could be hybridized with different renewable options e.g., flat plate collectors with parabolic trough collectors, or combined with existing heat supplies (e.g., fossil fuels), to give options for targeted SIPH applications, industrial decarbonization and the reduction of fuel consumption. Hybrid solutions and thermal energy storage will be important for the dispatch of heat at optimal times needed by the demand side of the buildings and industrial applications. At present, there is no integrated modeling tool for hybrid RTES, and this paper highlights the development of a renewable thermal hybridization framework for IPH use that is built from existing tools like System Advisor Model. The long-term vision for the framework (through significant further research) is to develop a coupled hybrid energy generation and cost analysis tool, where the tool could help the user in determining the most suitable and cost-effective technologies for their applications. Ongoing work will look to add costs for RTES options and further refinement on the selection of suitable technologies. This future tool could calculate the levelized cost of heat of various RTES hybrid options, by taking the user's solar resource, fuel costs, industrial heat demand profile, available land, and other factors into account to determine the applicability into their process.

concentrating solar power↗

Equilibration of the planar modes of ultracold two-dimensional ion crystals in a Penning trap

In this work, planar thermal equilibration is studied using direct numerical simulations of ultracold two-dimensional ion crystals in a Penning trap with a rotating wall. The large magnetic field of the trap splits the modes that describe in-plane motion of the ions into two branches: high-frequency cyclotron modes dominated by kinetic energy and low-frequency E × B modes dominated by potential energy associated with thermal position displacements. Using an eigenmode analysis we extract the equilibration rate between these two branches as a function of the ratio of the frequencies that characterize the two branches and observe this equilibration rate to be exponentially suppressed as the ratio increases. Under experimental conditions relevant for current work at NIST, the predicted equilibration time is orders of magnitude longer than any relevant experimental timescales. We also study the coupling rate dependence on the initial energy of the planar modes and the number of ions. In addition, we show how increasing the rotating wall strength improves crystal stability. These details of in-plane mode dynamics help set the stage for developing strategies to efficiently cool the in-plane modes and improve the performance of single plane ion crystals for quantum information processing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Control Algorithms Charge Batteries Faster

On March 29, 2011, NASA s Mercury Surface, Space Environment, Geochemistry and Ranging (MESSENGER) spacecraft beamed a milestone image to Earth: the first photo of Mercury taken from orbit around the solar system s innermost planet. (MESSENGER is also the first spacecraft to orbit Mercury.) Like most of NASA s deep space probes, MESSENGER is enabled by a complex power system that allows its science instruments and communications to function continuously as it travels millions of miles from Earth. "Typically, there isn't one particular power source that can support the entire mission," says Linda Taylor, electrical engineer in Glenn Research Center s Power Systems Analysis Branch. "If you have solar arrays and you are in orbit, at some point you re going to be in eclipse." Because of this, Taylor explains, spacecraft like MESSENGER feature hybrid power systems. MESSENGER is powered by a two-panel solar array coupled with a nickel hydrogen battery. The solar arrays provide energy to the probe and charge the battery; when the spacecraft s orbit carries it behind Mercury and out of the Sun s light, the spacecraft switches to battery power to continue operations. Typically, hybrid systems with multiple power inputs and a battery acting alternately as storage and a power source require multiple converters to handle the power flow between the devices, Taylor says. (Power converters change the qualities of electrical energy, such as from alternating current to direct current, or between different levels of voltage or frequency.) This contributes to a pair of major concerns for spacecraft design. "Weight and size are big drivers for any space application," Taylor says, noting that every pound added to a space vehicle incurs significant costs. For an innovative solution to managing power flows in a lightweight, cost-effective manner, NASA turned to a private industry partner.

Source record↗

The Reaction Rates of Amidogen and Ammonia with Nitrous Oxide: Implications for Combustion Mechanisms

Pulsed laser photolysis experiments with laser-induced fluorescence detection of NH 2 set an upper limit to the rate constant for reaction with N 2 O of k < 1 × 10 −15 cm 3 molecule −1 s −1 at 513 K. Computations were based on geometries and anharmonic frequency analysis (B2PLYP-D3/cc-pVTZ) followed by coupled cluster calculations extrapolated to the infinite basis set limit, with corrections for core−valence electron correlation, scalar relativistic effects, and correlation up to CCSDT(Q). Species that showed multireference character were quantified with MRCI(7,7)+Q/cc-pVTZ theory. Rate constants were obtained for the dominant product channel H 2 NN + NO, along with HN 3 + OH, H 2 NO + N 2 , NNH + HNO, ON(NH) 2 and HNNH + NO. The last channel is slow even at 2500 K, contrary to an early empirical estimate and confirming recent suggestions. Modeling of literature experiments on oxidation of NH 3 by N 2 O shows that all channels are too slow to make a significant impact on the loss of N 2 O in ammonia flames. Similarly, the direct NH 3 + N 2 O reaction is found to be negligibly slow.

Ammonia↗

Subsystem Estimation: A Modal Decomposition Approach

In structural analysis it is often the requirement that the response of a critical subsystem must be extracted from the response of the overall structure even when not directly excited or sometimes even measured. Based on this requirement, a number of signal processing approaches evolve that can provide a meaningful solution to this troubling problem; however, there is a variety of prior information (parameters or data) that must be available. For instance, knowledge of subsystem bandwidth may be enough to design a simple filter if there is little or no overlap with other component systems. If this is not satisfactory and a specific spectral response is required, then a model-based approach can be applied using an available free-free response of the subsystem from calibration or acceptance tests. However, when assembled and coupled with the overall system, the subsystem boundary conditions have been altered and therefore its response is altered as well.

42 ENGINEERING↗

Static aeroelastic analysis of wings using Euler/Navier-Stokes equations coupled with improved wing-box finite element structures

Strong interactions between flow about an aircraft wing and the wing structure can result in aeroelastic phenomena which significantly impact aircraft performance. Time-accurate methods for solving the unsteady Navier-Stokes equations have matured to the point where reliable results can be obtained with reasonable computational costs for complex non-linear flows with shock waves, vortices and separations. The ability to combine such a flow solver with a general finite element structural model is key to an aeroelastic analysis in these flows. Earlier work involved time-accurate integration of modal structural models based on plate elements. A finite element model was developed to handle three-dimensional wing boxes, and incorporated into the flow solver without the need for modal analysis. Static condensation is performed on the structural model to reduce the structural degrees of freedom for the aeroelastic analysis. Direct incorporation of the finite element wing-box structural model with the flow solver requires finding adequate methods for transferring aerodynamic pressures to the structural grid and returning deflections to the aerodynamic grid. Several schemes were explored for handling the grid-to-grid transfer of information. The complex, built-up nature of the wing-box complicated this transfer. Aeroelastic calculations for a sample wing in transonic flow comparing various simple transfer schemes are presented and discussed.

Guruswamy, Guru P.↗

Skin Temperature Analysis and Bias Correction in a Coupled Land-Atmosphere Data Assimilation System

In an initial investigation, remotely sensed surface temperature is assimilated into a coupled atmosphere/land global data assimilation system, with explicit accounting for biases in the model state. In this scheme, an incremental bias correction term is introduced in the model's surface energy budget. In its simplest form, the algorithm estimates and corrects a constant time mean bias for each gridpoint; additional benefits are attained with a refined version of the algorithm which allows for a correction of the mean diurnal cycle. The method is validated against the assimilated observations, as well as independent near-surface air temperature observations. In many regions, not accounting for the diurnal cycle of bias caused degradation of the diurnal amplitude of background model air temperature. Energy fluxes collected through the Coordinated Enhanced Observing Period (CEOP) are used to more closely inspect the surface energy budget. In general, sensible heat flux is improved with the surface temperature assimilation, and two stations show a reduction of bias by as much as 30 Wm(sup -2) Rondonia station in Amazonia, the Bowen ratio changes direction in an improvement related to the temperature assimilation. However, at many stations the monthly latent heat flux bias is slightly increased. These results show the impact of univariate assimilation of surface temperature observations on the surface energy budget, and suggest the need for multivariate land data assimilation. The results also show the need for independent validation data, especially flux stations in varied climate regimes.

Bosilovich, Michael G.↗

Direct observation of magnetic ordering induced via systematic lattice disorder in artificial rhombus spin ices

It is critical to understand the effect of lattice geometry on the order parameter of a condensed matter system, as it controls phase transitions in such systems. Artificial spin ices (ASIs) are two-dimensional lattices of Ising-like nanomagnets that provide an opportunity to explore such phenomena by lithographically controlling the lattice geometry to observe its influence on magnetic ordering and frustration effects. Here we report a systematic approach to studying the effects of disorder in rhombus ASIs generated from combinations of five vertex motifs. We investigate four geometries characterized by a geometric order parameter, with symmetries ranging from periodic to quasiperiodic to random. Lorentz transmission electron microscopy data indicates magnetic domain behavior depends on chains of strongly-coupled islands in the periodic and sixfold-twinned lattices, while the behavior of the disordered lattice is dominated by vertex motifs with large configurational degeneracy. Utilizing micromagnetic simulations, a quantitative analysis of the lattice energetics showed that the experimental rotationally-demagnetized state of the disordered ASI was closer in energy to the idealized ground state compared to other periodic and twinned ASIs. Finally, our work provides a unique pathway for using degeneracy, magnetic frustration, and order to control the magnetization behavior of designer disordered systems.

36 MATERIALS SCIENCE↗

On Boundary-Value Problems for RANS Equations and Two-Equation Turbulence Models

Currently, in engineering computations for high Reynolds number turbulent flows, turbulence modeling continues to be the most frequently used approach to represent the effects of turbulence. Such models generally rely on solving either one or two transport equations along with the Reynolds-Averaged Navier–Stokes (RANS) equations. The solution of the boundary-value problem of any system of partial differential equations requires the complete delineation of the equations and the boundary conditions, including any special restrictions and conditions. In the literature, such a description is often incomplete, neglecting important details related to the boundary conditions and possible restrictive conditions, such as how to ensure satisfying prescribed values of the dependent variables of the transport equations in the far field of a finite domain. In this article, we discuss the possible influence of boundary values, as well as near-field and far-field behavior, on the solution of the RANS equations coupled with transport equations for turbulence modeling. In so doing, we defne the concept of a welldefined boundary-value problem. Additionally, a three-dimensional, rather than a simpler one-dimensional analysis is performed to analyze the near-wall and far-field behavior of the turbulence model variables. This allows an assessment of the decay rate of these variables required to realize the boundary conditions in the far field. This paper also addresses the impact of various transformations of two-equation models (e.g., the model of Wilcox) to remove the singular behavior of the dissipation rate (ω) at the surface boundary. Finally, the issue of well-posedness regarding the governing equations is considered. A compelling argument (although not a proof) for ill-posedness is made for both direct and inverse problems.

Navier–Stokes equations↗

Atmospheric humidity deficits tell us how soil moisture deficits down-regulate ecosystem evaporation

Knowing how actual evaporation is down-regulated from potential evaporation during periods with soil moisture deficits is one of the greatest challenges towards computing evaporation everywhere on a regular basis. We propose the hypothesis that vegetated landscapes transmit information on soil moisture deficits through its down-regulation of evaporation, which in turn affects the humidification and growth of the planetary boundary layer. To test this hypothesis, here, we examined how the evaporative fraction, defined as the ratio between actual and potential evaporation, corresponds with the soil moisture stress index, defined as the relative humidity with a power law dependence on vapor pressure deficit times a parameter, β. We tested the parameterized soil moisture stress index with direct eddy covariance measurements of actual evaporation and computations of potential evaporation based on meteorological conditions, averaged on monthly time steps. The analysis was conducted on a dataset obtained from 144 FLUXNET sites. These sites spanned much of the world's climates and biomes and contained over 5900 months of observations. Observations of the evaporative fraction fit the model of the soil moisture stress index best for semi-arid and arid ecosystems, and least for the humid and wet ecosystems. Consequently, a significant relationship between measurements of the evaporative fraction and soil moisture stress index held for about one-half of the population of sites. Under this condition, the median value of the parameter, β, was 1.41. To investigate the mechanism of this empirical soil moisture stress index, we diagnosed it with a coupled evaporation-planetary boundary layer model. The soil moisture stress index is strongly related to surface resistance, as defined by inverting Penman-Monteith equation. Consequently, this index provides an independent estimate of surface resistance based on easy to measure mean monthly weather conditions like relative humidity and temperature. Thereby, this soil moisture stress index has potential to be applied to weather, climate and biogeochemical models and the interpretation of satellite derived evaporation products, like the one provided by the ECOSTRESS mission.

54 ENVIRONMENTAL SCIENCES↗

An Update on Cardinal: Toward Full Core Pebble Simulations

The new DOE-NE Center of Excellence for Thermal-fluids applications in Nuclear Energy inaugurated in April 2018 considers and researches novel new solution strategies for historically challenging flow issues that still plague the current fleet of deployed Light Water nuclear Reactors (LWRs) as well as predicting various fluid flow and fluid related issues with advanced reactor technologies. Our advanced thermal-fluids research and development approach synergistically combines three natural, though overlapping, length and time scales in a hierarchal multi-scale approach to avoid the temptation and pitfalls of attempting to develop a single solve-all algorithm for physical fluid flow problems that will span nine orders of magnitude in spatial and temporal scales. To demonstrate the multi-scale philosophy of the center we focus on Fluoride Cooled High Temperature Reactors (FHRs), and in particular on the Berkley’s PB-FHR Mark-I design. The Fluoride salt cooled High temperature Reactor (FHR) is a class of advanced nuclear reactors that combine the robust coated particle fuel form from high temperature gas cooled reactors, direct reactor auxiliary cooling system (DRACS) passive decay removal of liquid metal fast reactors, and the transparent, high volumetric heat capacitance liquid Fluoride salt working fluids - such as Flibe - from molten salt reactors. This combination of fuel and coolant enables FHRs to operate in a high-temperature low-pressure design space that has beneficial safety and economic implications. The PB-FHR reactor relies on a pebble bed approach and pebble bed reactors are in a sense the poster child for multi-scale analysis. The lower length-scale simulator for pebble reactor cores comprises three physics: neutronics (OpenMC), thermal-fluids (Nek5000/NekRS) and fuel performance (BISON). As part of Center of Excellence ongoing research efforts we have developed Cardinal, a new tool platform for lower length-scale simulation. Cardinal tightly couples all three physics and leverages advances in MOOSE such as the MultiApp system and the concept of MOOSE-wrapped Apps. The present report provides an update on the development of Cardinal with an extension of Cardinal to GPUs. We also perform a first-of-a-kind demonstration simulation on Summit representing a 10x capability increase in terms of pebble count for Cardinal.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗