Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “characteristic equations”

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 253 records · Page 14

Multi-reward reinforcement learning based development of inter-atomic potential models for silica

Abstract Silica is an abundant and technologically attractive material. Due to the structural complexities of silica polymorphs coupled with subtle differences in Si–O bonding characteristics, the development of accurate models to predict the structure, energetics and properties of silica polymorphs remain challenging. Current models for silica range from computationally efficient Buckingham formalisms (BKS, CHIK, Soules) to reactive (ReaxFF) and more recent machine-learned potentials that are flexible but computationally costly. Here, we introduce an improved formalism and parameterization of BKS model via a multireward reinforcement learning (RL) using an experimental training dataset. Our model concurrently captures the structure, energetics, density, equation of state, and elastic constants of quartz (equilibrium) as well as 20 other metastable silica polymorphs. We also assess its ability in capturing amorphous properties and highlight the limitations of the BKS-type functional forms in simultaneously capturing crystal and amorphous properties. We demonstrate ways to improve model flexibility and introduce a flexible formalism, machine-learned ML-BKS, that outperforms existing empirical models and is on-par with the recently developed 50 to 100 times more expensive Gaussian approximation potential (GAP) in capturing the experimental structure and properties of silica polymorphs and amorphous silica.

36 MATERIALS SCIENCE↗

Theory of mean E × B shear in a stochastic magnetic field: ambipolarity breaking and radial current

The mean E $ \times $ B shear in a stochastic magnetic field is calculated, using the radial force balance relation and transport equations. This analysis is relevant to the L → H transition with resonant magnetic perturbations, and special focus is placed upon the physics of non-ambipolar transport and radial current. The key physical process is the flow of fluctuating currents along wandering magnetic fields. The increments in poloidal and toroidal rotation, density and ion pressure are calculated. The radial envelope of the magnetic perturbations inside the plasma defines a new scale ${\ell _{{\text{env}}}}$, which is the characteristic scale of the magnetic fluctuation intensity profile. Here, the net particle outflow due to stochastic magnetic fields is calculated and is determined by the net radial current through the separatrix. Implications for the L → H transition are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Seismic strain energy partitioning: estimating the strain energy of seismic body waves.

This report details a method to estimate the energy content of various types of seismic body waves. The method is based on the strain energy of an elastic wavefield and Hooke’s Law. We present a detailed derivation of a set of equations that explicitly partition the seismic strain energy into two parts: one for compressional (P) waves and one for shear (S) waves. We posit that the ratio of these two quantities can be used to determine the relative contribution of seismic P and S waves, possibly as a method to discriminate between earthquakes and buried explosions. We demonstrate the efficacy of our method by using it to compute the strain energy of synthetic seismograms with differing source characteristics. Specifically, we find that explosion-generated seismograms contain a preponderance of P wave strain energy when compared to earthquake-generated synthetic seismograms. Conversely, earthquake-generated synthetic seismograms contain a much greater degree of S wave strain energy when compared to explosion-generated seismograms.

58 GEOSCIENCES↗

Disentangling natural and anthropogenic effects on benthic macroinvertebrate assemblages in western US streams

Abstract Stream macroinvertebrate assemblages are shaped by natural and human‐related factors that operate through complex hierarchical pathways. Quantifying these relationships can provide additional insights into stream ecological assessment. We applied a structural equation modeling framework to evaluate hypothesized pathways by which watershed, riparian, and in‐stream factors affect benthic macroinvertebrate condition in the Western Mountains (WMT) and Xeric (XER) ecoregions in the United States. We developed a conceptual model grounded in theory, empirical evidence, and expert opinion to evaluate the following hypotheses: (1) macroinvertebrate assemblages are primarily driven by proximal, in‐stream factors (e.g., water quality and physical habitat); (2) anthropogenic land uses affect macroinvertebrates indirectly by altering in‐stream characteristics; and (3) riparian vegetation cover attenuates land use effects. We tested our model separately on three measures of benthic macroinvertebrate assemblage condition: ratio of observed‐to‐expected taxonomic richness (O/E); a multimetric index (MMI); and richness of Ephemeroptera, Plecoptera, and Trichoptera taxa (EPT). In the WMT, site‐level riparian cover, in‐stream physical habitat (relative bed stability), and water chemistry (total nitrogen) were the top three predictors of macroinvertebrate assemblages, each having over two times the magnitude of effect on macroinvertebrates compared with watershed‐level predictors. In the arid XER, annual precipitation and stream flow characteristics were top predictors of macroinvertebrate assemblages and had similar magnitudes of effect as in‐stream water chemistry. Path analyses revealed that land use activities in the watershed and at the stream site degraded macroinvertebrate assemblages indirectly by altering relative bed stability, water quality, and riparian cover/complexity. Increased riparian cover was associated with greater macroinvertebrate condition by reducing land use impacts on stream flow, streambed substrate, and water quality, but the pathways differed among ecoregions. In the WMT, site‐level riparian cover affected macroinvertebrate assemblages partly through indirect pathways associated with greater streambed stability and reduced total nitrogen concentrations. In contrast, in the XER, watershed‐level riparian cover affected macroinvertebrate assemblages through greater specific stream power. Identifying the relative effects of and pathways by which natural and anthropogenic factors affect macroinvertebrates can serve as a framework for prioritizing management and conservation efforts.

Environmental Sciences & Ecology↗

Rigid-flexible coupling multi-body dynamics modeling of a semi-submersible floating offshore wind turbine

A 14-degree-of-freedom (14-DOF) flexible multibody dynamics model is developed and verified for a semi-submersible floating offshore wind turbine (FOWT). The model considers the coupled dynamics of the platform, tower, nacelle, blades, and mooring subjected to external wind and wave loads. The platform is simplified as a rigid body connected to the seabed by the mooring system. The tower and blade are modeled as flexible cantilever beams. The energy method is used to derive the governing equations of motion, where the kinetic energy, potential energy, and work done by external wind-wave forces are all deduced in a global coordinate system. The 5-MW baseline semi-submersible FOWT is used to verify the derived model against the results simulated from FAST developed by the US National Renewable Energy Laboratory (NREL) at two scenarios: free decay state and different wind-wave load cases. Results show that the established FOWT model can well reflect the vibration characteristics of FAST model. The application of the model to the control of platform pitch with a tuned mass damper is studied. Furthermore, the simplified model could provide a low-order method for the structural dynamics analysis and advanced vibration control design for the multi-body components of the semi-submersible FOWT in the future.

17 WIND ENERGY↗

Protection and Restoration Solutions to Reliable and Resilient Integration of Grid-connected PV Installations and Distributed Energy Resources: Design, Testbed, Proof of Work and Impact Studies (Final Report)

To gain a better understanding of the complex transients in a utility grid with a large number of solar PV installations coupled by PWM inverters during the protection and restoration period of the power grid, in this project a kW-level experimental system with dominating inverter-based resources and a hardware-in-loop simulation system with transmission and distribution models, as well as several SEL protective relays have been built and used. The major research findings of the project are summarized below in two perspectives: Experimental research: A self-organized ultra-high frequency solitary waveform is discovered and demonstrated in the hardware experimental system. Despite the familiarities, such a solitary waveform is distinct from any harmonics, transient, or resonance waves in that, it is 1) non-dispersive over time and space, 2) not associated with any active source or the linear superposition of sources, 3) half-cycle asymmetric, 4) not responding to filters or change of system characteristic resonance frequency, 5) ubiquitous as it occurs simultaneously everywhere in the system from DC supply, power lines, and the utility grid, and 6) explosive through tripping the protection or damaging susceptible devices or circuits. Analytical study: The nature of such a new waveform and an analytical explanation of its formation are studied based on related physics and non-linear science principles, with the following highlights: 1) such a new waveform follows the solutions to the Non-linear Schrödinger equation, so the classic linear perturbation theory is unable to explain or predict such a unique waveform as confirmed with the research team of RTDS; 2) the critical condition of occurrence of such a waveform is derived, which indicates that the breakings of solitary wave is a system synchronous issue between the utility grid in the 60Hz phasor domain and duty-ratio modulation of DC sources in the inverter switching frequency domain; 3) such a waveform carries energy mass so it could be detrimental; 4) such a waveform is deceiving as it is not readily detectable in the energy propagation direction by the primary protection equipment, while it is detrimental in the perpendicular direction of energy propagation, i.e. voltage direction. Thus, it has more likely challenges to the second primary equipment and devices, particularly at the weakest link and point such as aging insulation and inappropriate setting of susceptible devices. Therefore, such a waveform can be easily ignored in the aftermath investigation. The intellectual merits: The experimental discovery reveals certain unfamiliar transient phenomena that could challenge the integration of large-scale grid-connected solar PV installation during the group ride-through period of solar PV installations, commissioning of large-solar farms, or dramatic change of solar radiation conditions. Particularly, the research findings suggest that the occurrence of the unfamiliar transient phenomena is uniquely associated with the inverter-based solar PV installation, which is less likely to happen for rotary energy systems. The physics and non-linear science-based research work laid down an analytical path to the challenging transient stability problems including those that have been observed currently and future calls. Both experimental and analytical research results explain the limitation of many current research effects including some DoE research undertakings as well as possible solutions. The broader impacts: The research outcome of the project advances the understanding of the possible transient problems for the integration of inverter-based solar PV installation. It also highlights the theoretical and technical barriers for applying the conventional theory and technology to design and implement countermeasures against their adverse impacts of the large-scale solar PV installation and operation on grid reliability and security. More broadly, the research outcomes have shed some light on the open, fundamental challenges in the integration of large-scale solar photovoltaic energy into current and next-generation power grids across the country, and worldwide. The advanced physics and non-linear science-based analysis open up a path to harmonization of the renewables for societal energy needs via building a resilient and sustainable electric power infrastructure.

14 SOLAR ENERGY↗

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. Furthermore, the thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 Wm -1 K -1 are estimated at ambient and extreme conditions (T = 4000K and P = 135 GPa), respectively.

36 MATERIALS SCIENCE↗

In-field critical current and pinning mechanisms at 4.2 K of Zr-added REBCO coated conductors

The critical current and pinning mechanisms at 4.2 K have been studied over a magnetic field range of 0–14 T for Zr-added (0, 5 and 15 mol.%) REBa 2 Cu 3 O 7-x (REBCO and RE = rare earth) coated conductors fabricated by advanced metal organic chemical vapor deposition (A-MOCVD). It is found that the (Ba + Zr)/Cu content in Zr-added (5 and 15 mol.%) REBCO affects the critical current at 77 K, 0 T as well as density, continuity and shape of BaZrO 3 (BZO) self-assembled nanocolumns and RE 2 O 3 in-plane precipitates that significantly enhance the pinning force density Fp(H) as well as isotropic pinning landscape at 4.2 K. In addition to bell-shape dependence of critical current density, J c , at 4.2 K with (Ba + Zr)/Cu content we observed an unusual Fp(H) behavior correlated to particular type of pinning centers, morphology and distribution that have been revealed by TEM microstructure analysis. By fitting the Dew–Hughes equation of the pinning force density Fp(H) at 4.2 K we extract the scaling behaviors of the Fp(H) associated with the competition of pinning mechanisms driven by vertically-aligned BZO nanorods and in-plane RE 2 O 3 pinning defects. This result sheds light on approaches towards interactive control of strong and isotropic pinning centers in Zr-added REBa 2 Cu 3 O 7-x (REBCO and RE = rare earth) coated conductors, and especially understanding the correlation between microstructural characteristics and vortex pinning mechanisms at 4.2 K in high magnetic fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Efficient Multigrid Reduction-in-Time for Method-of-Lines Discretizations of Linear Advection

Parallel-in-time methods for partial differential equations (PDEs) have been the subject of intense development over recent decades, particularly for diffusion-dominated problems. It has been widely reported in the literature, however, that many of these methods perform quite poorly for advection-dominated problems. In this report we analyze the particular iterative parallel-in-time algorithm of multigrid reduction-in-time (MGRIT) for discretizations of constant-wave-speed linear advection problems. We focus on common method-of-lines discretizations that employ upwind finite differences in space and Runge-Kutta methods in time. Using a convergence framework we developed in previous work, we prove for a subclass of these discretizations that, if using the standard approach of rediscretizing the fine-grid problem on the coarse grid, robust MGRIT convergence with respect to CFL number and coarsening factor is not possible. This poor convergence and non-robustness is caused, at least in part, by an inadequate coarse-grid correction for smooth Fourier modes in space-time known as characteristic components. We propose an alternative coarse-grid operator that provides a better correction of these modes. This coarse-grid operator is related to previous work and uses a semi-Lagrangian discretization combined with an implicitly treated truncation error correction. Theory and numerical experiments show the proposed coarse-grid operator yields fast MGRIT convergence for many of the method-of-lines discretizations considered, including for both implicit and explicit discretizations of high order. Parallel results demonstrate speed-up over sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗

A multiscale insight into the growth of h-BN: effect of the enclosure

Abstract There is a lack of knowledge on the fundamental growth mechanisms governing the characteristics of 2D materials synthesized by the chemical vapor deposition (CVD) technique and their correlation with experimentally controllable parameters, which hindered their wafer-scale synthesis. Here, we pursued an analytical and computational approach to access the system states that are not experimentally viable to address these critical needs. We developed a multiscale computational framework correlating the macroscale heat and mass flow with the mesoscale morphology of the as-grown 2D materials by solving the coupled system of heat/mass transfer and phase-field equations. We used hexagonal boron nitride (h-BN) as our model material and investigated the effect of substrate enclosure on its growth kinetics and final morphology. We revealed a lower concentration with a more uniform distribution on the substrate in an enclosed-growth than open-growth. It leads to a more uniform size distribution of the h-BN islands, consistent with existing experimental investigations.

Materials Science↗

Shock wave formation in radiative plasmas

The temporal evolution of weak shocks in radiative media is theoretically investigated in this work. The structure of radiative shocks has traditionally been studied in a stationary framework. Their systematic classification is complex because layers of optically thick and thin regions alternate to form a radiatively-driven precursor and a temperature-relaxation layer, between which the hydrodynamic shock is embedded. In this work, we analyze the formation of weak shocks when two radiative plasmas with different pressures are put in contact. Applying a reductive perturbative method yields a Burgers-type equation that governs the temporal evolution of the perturbed variables including the radiation field. The conditions upon which optically thick and thin solutions exist have been derived and expressed as a function of the shock strength and Boltzmann number. Below a certain Boltzmann number threshold, weak shocks always become optically thick asymptotically in time, while thin solutions appear as transitory structures. The existence of an optically thin regime is related to the presence of an overdense layer in the compressed material. Scaling laws for the characteristic formation time and shock width are provided for each regime. The theoretical analysis is supported by FLASH simulations, and a comprehensive testcase has been designed to benchmark radiative hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Study of Propane Combustion in a Spark-Ignited Cooperative Fuel Research (CFR) Engine

Liquefied petroleum gas (LPG), whose primary composition is propane, is a promising candidate for heavy-duty vehicle applications as a diesel fuel alternative due to its CO 2 reduction potential and high knock resistance. To realize diesel-like efficiencies, spark-ignited LPG engines are proposed to operate near knock-limit over a wide range of operating conditions, which necessitates an investigation of fuel-engine interactions that leads to end-gas autoignition with propane combustion. This work presents both experimental and numerical studies of stoichiometric propane combustion in a sparkignited (SI) cooperative fuel research (CFR) engine. Engine experiments are initially conducted at different compression ratio (CR) values, and the effects of CR on engine combustion are characterized. A three-pressure analysis (TPA) model based on the two-zone combustion concept is developed in GT-Power and validated using test results to estimate in-cylinder wall temperatures, residual gas fraction, etc. This model is further utilized to examine end-gas chemistry by enabling the SI turbulent flame combustion and unburned gas chemical kinetics modules. Finally, a three-dimensional (3D) computational fluid dynamic (CFD) model of the CFR engine is developed in CONVERGE, where the G-equation and SAGE detailed chemical kinetics models are implemented for combustion modeling. Here, a 153 species reduced chemical kinetics mechanism derived from the detailed NUIGMech1.1 mechanism based on the ignition delay and laminar flame speed (LFS) studies is used to generate an LFS lookup table and to describe end-gas autoignition chemistry. Multi-cycle Reynolds-averaged Navier-Stokes (RANS) simulations are then performed for the tested CRs, and the numerical model is shown to be capable of predicting the propane combustion characteristics, particularly the end-gas autoignition behavior.

20 FOSSIL-FUELED POWER PLANTS↗

A novel physics-based and data-supported microstructure model for part-scale simulation of laser powder bed fusion of Ti-6Al-4V

The elasto-plastic material behavior, material strength and failure modes of metals fabricated by additive manufacturing technologies are significantly determined by the underlying process-specific microstructure evolution. In this work a novel physics-based and data-supported phenomenological microstructure model for Ti-6Al-4V is proposed that is suitable for the part-scale simulation of laser powder bed fusion processes. The model predicts spatially homogenized phase fractions of the most relevant microstructural species, namely the stable β-phase, the stable α s -phase as well as the metastable Martensite α m -phase, in a physically consistent manner. In particular, the modeled microstructure evolution, in form of diffusion-based and non-diffusional transformations, is a pure consequence of energy and mobility competitions among the different species, without the need for heuristic transformation criteria as often applied in existing models. The mathematically consistent formulation of the evolution equations in rate form renders the model suitable for the practically relevant scenario of temperature- or time-dependent diffusion coefficients, arbitrary temperature profiles, and multiple coexisting phases. Due to its physically motivated foundation, the proposed model requires only a minimal number of free parameters, which are determined in an inverse identification process considering a broad experimental data basis in form of time-temperature transformation diagrams. Subsequently, the predictive ability of the model is demonstrated by means of continuous cooling transformation diagrams, showing that experimentally observed characteristics such as critical cooling rates emerge naturally from the proposed microstructure model, instead of being enforced as heuristic transformation criteria. Eventually, the proposed model is exploited to predict the microstructure evolution for a realistic selective laser melting application scenario and for the cooling/quenching process of a Ti-6Al-4V cube of practically relevant size. Numerical results confirm experimental observations that Martensite is the dominating microstructure species in regimes of high cooling rates, e.g., due to highly localized heat sources or in near-surface domains, while a proper manipulation of the temperature field, e.g., by preheating the base-plate in selective laser melting, can suppress the formation of this metastable phase.

Inverse parameter identification↗

MPACT Software Management Plan (V.4.3)

The MPACT code solves a discretized form of the Boltzmann transport equation on a wide variety of geometries and is distributed with a multigroup neutron cross section library. MPACT provides an advanced geometrically resolved neutral-particle transport capability to solve the flux distribution throughout the entire problem geometry, and it can model the isotopic depletion, decay, and activation of materials. The flux solution in MPACT is provided using a 2D/1D synthesis method within the framework of the 3D coarse mesh finite difference (CMFD) method for which axial and radial correction factors are obtained from 2D method of characteristics (MOC) and 1D nodal expansion method (NEM), PN, or SN. Other key characteristics of the MPACT code include the subgroup method and the embedded self-shielding method (ESSM) for resonance treatment, depletion capability based on the ORIGEN exponential matrix method, and a simplified thermal-hydraulics method for temperature/fluid feedback. The sole purpose of the simplified feedback model is to provide a mechanism for testing during code development and to provide a limited capability for educational applications. Work performed at the code level supports the VERA-QA-001, quality assurance program plan (QAPP) and VERA-QA-002, VERA Software Quality Assurance Plan.

97 MATHEMATICS AND COMPUTING↗

Carrier Dynamics of Polar, Semipolar, and Nonpolar InGaN/GaN LEDs Measured by Small-Signal Electroluminescence

The carrier dynamics in InGaN/GaN light-emitting diodes (LEDs) are directly tied to their efficiency and maximum modulation speed, which are important metrics for solid-state lighting, displays, and optical communication. In this work, we measure the carrier dynamics of a variety of InGaN/GaN LEDs using small-signal electroluminescence methods [1]. A rate equation approach and associated small-signal circuit are used to model carrier injection, recombination in the active region, recombination in the cladding regions, and carrier escape. The model is fit to the measured optical frequency response (S21) and input impedance (S11) of the LEDs to extract the various carrier lifetimes, the carrier density, and the radiative and non-radiative recombination rates. We specifically study planar nonpolar and semipolar LEDs, which show record-high modulation speeds for III-nitride LEDs and present the modulation characteristics of core-shell nanowire-based LEDs. The planar nonpolar m-plane ($101\bar{0}$) micro-LEDs achieve a record-high -3dB modulation bandwidth for a III-nitride LED of 1.5 GHz [2]. The -3dB response of an electrically injected nanowire-based micro-LED with nonpolar facets is also reported, showing a -3dB bandwidth of 1.2 GHz [3]. The high speed is attributed to the shorter carrier lifetime associated with the nonpolar orientation. We also study the carrier dynamics in semipolar ($20\bar{2}\bar{1}$) LEDs for various temperatures [4]. Finally, we present carrier dynamics measurements on commercial-grade c-plane epitaxy for various active region designs, including a wavelength series and a growth quality series. The wavelength series offers insight into the contributions of the quantum confined Stark effect (QCSE) and InGaN material quality on the green gap [5]. The growth quality series investigates the role of non-radiative centers on the LED performance. Extraction of the carrier dynamics using small-signal electroluminescence offers insight into the factors limiting the efficiency and high-speed performance of III-nitride emitters and can be leveraged to ultimately improve the devices.

LEDs, InGaN, electroluminescence, bandwidth, modul↗

Projected Cosmological Constraints from Strongly Lensed Supernovae with the Roman Space Telescope

One of the primary mission objectives of the Roman Space Telescope is to investigate the nature of dark energy with a variety of methods. Observations of Type I supernovae (SNe Ia) will be one of the principal anchors of the Roman cosmology program through traditional luminosity distance measurements. This SNe Ia cosmology program can provide another valuable cosmological probe, without altering the strategy of the mission: time delay cosmography with gravitationally lensed supernova (SN). In this work, we forecast lensed SN cosmology constraints with the Roman Space Telescope, while providing useful tools for future work. Using the anticipated characteristics of the Roman SNe Ia survey, we have constructed mock catalogs of expected resolved lensing systems, as well as strongly lensed Type Ia and core-collapse (CC) SN light curves, including microlensing effects. We predict Roman will find ~11 lensed SNe Ia and ~20 CCSNe, depending on the survey strategy. Next, we estimate the time delay precision obtainable with Roman (Ia: ~2 days, CC: ~3 days), and use a Fisher matrix analysis to derive projected constraints on H 0 ,Ω m , and the dark energy equation of state, w, for each SNe Ia survey strategy. Here, a strategy optimized for the discovery of high-redshift SNe Ia is preferred when considering the constraints possible from both SNe Ia and lensed SN cosmology, also delivering ~1.5 times more lensed SNe than other proposed survey strategies.

79 ASTRONOMY AND ASTROPHYSICS↗

Boundary-layer receptivity to oblique freestream vorticity waves for a high-enthalpy hypersonic flow

The receptivity of a Mach 15 straight-cone boundary layer to oblique freestream vorticity waves is investigated using direct numerical simulation (DNS) alongside linear stability theory and the linear parabolized stability equations. A thermochemical nonequilibrium gas model is used. Oblique freestream vorticity waves at frequencies of 400, 800, and 1200 kHz are considered, with incident angles ranging from 0° to 29.4° at 400 kHz, 0° to 15.7° at 800 kHz, and 0° to 20.6° at 1200 kHz. The 400 kHz case is of primary interest due to the strong second-mode amplification at this frequency. Although the underlying base flow is axisymmetric, the oblique vorticity waves lead to a fully three-dimensional boundary-layer disturbance whose characteristics vary depending on the azimuthal ray relative to the freestream wave. Moving downstream within the second-mode instability region, some clear trends emerge in terms of the boundary-layer disturbance amplitudes; that is, disturbance amplitudes are highest at the leeward ray (relative to the freestream wave), but weakest about halfway between the windward and leeward rays. Increasing the incident angle causes the amplitudes to increase on the leeward ray and decrease on the windward ray. Moreover, the boundary-layer disturbance throughout contains a wide spectrum of azimuthal wavenumbers in which the disturbance energy falls off at higher wavenumbers. Increasing the incident angle causes the azimuthal spectrum of the boundary-layer disturbance to broaden overall. Qualitatively similar results are found for the two higher frequencies leading up to the peak-amplitude locations corresponding to the second-mode instability.

42 ENGINEERING↗

Resonant x-ray emission across the L 3 edge of uranium compounds

A narrow bandwidth x-ray beamline and a multi-crystal von Hamos spectrometer were used to record x-ray absorption and x-ray emission (XES) across the uranium L 3 edge of UO 2 , UO 3 , Cs 2 UO 2 Cl 4 , and Cs 2 UCl 6 . Measurements were made over 17150–17250 eV with an instrumental resolution of ∼2 eV. This resolution allowed the use of Lorentzian peak fits to the L 3 N 4 and L 3 N 5 characteristic x-ray fluorescence lines that have ∼12–13 eV lifetime widths. The fluorescence yields of the four compounds display strong white lines near threshold and multiple scattering features at higher energies. The XES spectra were fit with three Lorentzians, two for the L 3 N 4 and L 3 N 5 lines and one for inelastic scattering by 2p–6d excitations and 4d–6d final states. The intensities of the white lines were largely due to the inelastic scattering component. To support the measurements, relativistic equation-of-motion coupled-cluster and restricted active space configuration interaction calculations were performed for uranium core-excitation energies in Cs 2 UO 2 Cl 4 and Cs 2 UCl 6 . These calculations with rigorous treatments of relativistic effects are shown to provide accurate uranium L 3 -edge binding energies, L 3 N 5 emission energies, as well as the energy losses in the inelastic scattering process. The measured and calculated results show variations between compounds with U(IV) oxidation states that contain 6d 0 5f 2 electrons in their ground configurations (UO 2 and Cs 2 UCl 6 ) compared with U(VI) compounds (UO 3 and Cs 2 UO 2 Cl 4 ) with empty 5f and 6d configurations.

74 ATOMIC AND MOLECULAR PHYSICS↗