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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS1961

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗
NASA NTRS2010

Utilizing In Situ Directional Hyperspectral Measurements to Validate Bio-Indicator Simulations for a Corn Crop Canopy

Two radiative transfer canopy models, SAIL and the two-layer Markov-Chain Canopy Reflectance Model (MCRM), were coupled with in situ leaf optical properties to simulate canopy-level spectral band ratio vegetation indices with the focus on the photochemical reflectance index in a cornfield. In situ hyperspectral measurements were made at both leaf and canopy levels. Leaf optical properties were obtained from both sunlit and shaded leaves. Canopy reflectance was acquired for eight different relative azimuth angles (psi) at three different view zenith angles (Theta (sub v)), and later used to validate model outputs. Field observations of photochemical reflectance index (PRI) for sunlit leaves exhibited lower values than shaded leaves, indicating higher light stress. Canopy PRI expressed obvious sensitivity to viewing geometry, as a function of both Theta (sub v) and psi . Overall, simulations from MCRM exhibited better agreements with in situ values than SAIL. When using only sunlit leaves as input, the MCRM-simulated PRI values showed satisfactory correlation and RMSE, as compared to in situ values. However, the performance of the MCRM model was significantly improved after defining a lower canopy layer comprised of shaded leaves beneath the upper sunlit leaf layer. Four other widely used band ratio vegetation indices were also studied and compared with the PRI results. MCRM simulations were able to generate satisfactory simulations for these other four indices when using only sunlit leaves as input; but unlike PRI, adding shaded leaves did not improve the performance of MCRM. These results support the hypothesis that the PRI is sensitive to physiological dynamics while the others detect static factors related to canopy structure. Sensitivity analysis was performed on MCRM in order to better understand the effects of structure related parameters on the PRI simulations. Leaf area index (LAI) showed the most significant impact on MCRM-simulated PRI among the parameters studied. This research shows the importance of hyperspectral and narrow band sensor studies, and especially the necessity of including the green wavelengths (e.g., 531 nm) on satellites proposing to monitor carbon dynamics of terrestrial ecosystems.

Cheng, Yen-Ben↗
DOE OSTI2024

Embedded symmetric positive semi-definite machine-learned elements for reduced-order modeling in finite-element simulations with application to threaded fasteners

Here, we present a machine-learning strategy for finite element analysis of solid mechanics wherein we replace complex portions of a computational domain with a data-driven surrogate. In the proposed strategy, we decompose a computational domain into an “outer” coarse-scale domain that we resolve using a finite element method (FEM) and an “inner” fine-scale domain. We then develop a machine-learned (ML) model for the impact of the inner domain on the outer domain. In essence, for solid mechanics, our machine-learned surrogate performs static condensation of the inner domain degrees of freedom. This is achieved by learning the map from displacements on the inner-outer domain interface boundary to forces contributed by the inner domain to the outer domain on the same interface boundary. We consider two such mappings, one that directly maps from displacements to forces without constraints, and one that maps from displacements to forces by virtue of learning a symmetric positive semi-definite (SPSD) stiffness matrix. We demonstrate, in a simplified setting, that learning an SPSD stiffness matrix results in a coarse-scale problem that is well-posed with a unique solution. We present numerical experiments on several exemplars, ranging from finite deformations of a cube to finite deformations with contact of a fastener-bushing geometry. We demonstrate that enforcing an SPSD stiffness matrix drastically improves the robustness and accuracy of FEM–ML coupled simulations, and that the resulting methods can accurately characterize out-of-sample loading configurations with significant speedups over the standard FEM simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Renewable Thermal Energy Systems Designed for Industrial Process Solutions in Multiple Industries: Preprint

Industrial decarbonization is a key area that must be accelerated, to foster the removal of fossil fuels from the provision of heat, especially at low temperatures less than 300 degrees C. This paper looks at the results of two case studies for understanding the economics and potential for renewable thermal energy systems (RTES), particularly in hybrid configurations to provide industrial process heat (IPH). The first case study looks at heat pumps for district heating, and the second, the use of linear Fresnel collectors (LFCs) coupled with phase change material (PCM) thermal energy storage (TES) for direct steam generation (DSG). Using district heat as an input for the heat pump, three cases were run harvesting energy from ambient water (5 degrees C), sewage water (20 degrees C), and a solar collector (35 degrees C). Accounting for elevated costs of infrastructure for each heat source, the levelized cost of heat (LCOH) of the first case study ranged from $4-$10 per million British Thermal Units (MMBTU). For the second case study modeling LFCs with PCM and TES, the results show that a LCOH of $9-$15 per MMBTU is possible, depending on the direct normal irradiance.

concentrated solar thermal↗
DOE OSTI2024

Critical Role of Framework Flexibility and Disorder in Driving High Ionic Conductivity in LiNbOCl 4

Understanding Li-ion transport is key for the rational design of superionic solid electrolytes with exceptional ionic conductivities. LiNbOCl 4 is reported to be one of the most highly conducting materials in the recently realized new class of soft oxyhalide solid electrolytes, exhibiting an ionic conductivity of ~11 mS·cm -1 . Here, we apply X-ray/neutron diffraction and pair distribution function analysis - coupled with density functional theory/ab-initio molecular dynamics - to determine a structural model that provides a rationale for the high conductivity that we observe experimentally in this nanocrystalline solid. We show that it arises from unusually high framework flexibility at room temperature. This owes to isolated 1-D [NbOCl 4 ] - anionic chains which exhibit energetically favorable orientational disorder that is - in turn - correlated to multiple, disordered and equi-energetic Li + sites in the lattice. As the Li-ions sample the 3-D energy landscape with a fast predicted diffusion coefficient of 5.1 x 10 -7 cm 2 /s at room temperature (σ i calc = 17.4 mS·cm -1 ), the inorganic polymer chains can reorient or vice versa. The activation energy barrier for Li migration through the frustrated energy landscape is especially reduced by the elastic nature of the NbO 2 Cl 4 octahedra evident from very widely dispersed Cl-Nb-Cl bond angles in AIMD snapshots at 300 K. The phonon spectra are predominantly influenced by Cl vibrations in the low energy range, and there is strong overlap between the framework (Cl, Nb) and Li partial pDOS in the region between 1.2 - 4.0 THz. The framework flexibility is also reflected in a relatively low bulk modulus of 22 GPa. In conclusion, our findings pave the way for investigation of future “flex-ion” inorganic solids and open up a new direction for the design of high conductivity, soft solid electrolytes for all-solid-state batteries.

AIMD↗
DOE OSTI2020

Consequence Analysis of Residual Water in a Storage Canister - 20484

Recent observations from projects sponsored under the DOE-Nuclear Energy Spent Fuel and Waste Disposition campaign show that residual free water, well above the 0.4 gm-moles that correspond to a 3 torr rebound pressure after a 30-minute hold following prototypic drying, may remain within SNF canisters used for long-term dry storage. Specifically, experimental observations in an Integrated Research Project showed subzero temperatures and a local pocket of frozen water in a mockup fuel assembly subjected to an industry practiced drying process. In gas sampling and analysis from the High Burnup (HBU) Demonstration project, using a fully loaded and dried SNF canister, indicated presence of 100 g of residual water had remained in the canister. The effects of residual water in terms of potential adverse impacts on the condition of the SNF and the canister internals during long-term dry storage has been rigorously evaluated to address inadvertent residual water contents. The impact of residual water on the fuel and canister internals, with a focus on the oxidation of the cladding and of postulated exposed (breached-clad) fuel, were evaluated using an integrated set of time-dependent environmental conditions with semi-empirical materials oxidation models. An 'integration model,' with coupled sub-models, developed and previously used by the Center for Nuclear Waste Regulatory Analysis (CNWRA), was applied and expanded with improved sub-models for water radiolysis, and for cladding and fuel oxidations. Cladding oxidation models for various alloys were developed using the rate constants for various cladding alloys in the dry storage temperature range. A burnup-dependent fuel oxidation model was developed; it accounted for inhibition effect of fission products on fuel oxidation rates. The potential for extension of the initial breach flaw (1 mm hole) was evaluated by two separate mechanics criteria. In the first case, a critical cladding strain criterion was used. In the second case, a critical stress intensity was used to evaluate stability of a planar axial flaw postulate. A failure condition exists when a sufficient amount of fuel oxidation (at the postulated breached-clad location) causes fuel pellet swelling to load the cladding to meet a mechanics criterion for flaw extension. The integration model's results provide the following insights. The radiolysis-generated oxygen is consumed with a partitioning between the cladding and the fuel. The extent of cladding oxidation is no more than 2 μm additional consumption of the cladding metal even with 10 moles of residual water, indicating that changes in cladding conditions due to water are expected to be negligible. The conditions for cladding failure (breach opening size extension) directly correlates with the extent of fuel oxidation. Fuel oxidation significant to extend the breach could occur even with 5.5 moles of residual water. The results indicate that fuel oxidation vis-a-vis cladding oxidation is more likely under a higher radiation field, i.e., when radiolytic decomposition of the residual water occurs in few years compared to several decades, and that the conditions to cause failure are increased with fuel and cladding temperatures approaching the peak storage temperature limit of 400 deg. C. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1992

Dynamic interactions between hypersonic vehicle aerodynamics and propulsion system performance

Described here is the development of a flexible simulation model for scramjet hypersonic propulsion systems. The primary goal is determination of sensitivity of the thrust vector and other system parameters to angle of attack changes of the vehicle. Such information is crucial in design and analysis of control system performance for hypersonic vehicles. The code is also intended to be a key element in carrying out dynamic interaction studies involving the influence of vehicle vibrations on propulsion system/control system coupling and flight stability. Simple models are employed to represent the various processes comprising the propulsion system. A method of characteristics (MOC) approach is used to solve the forebody and external nozzle flow fields. This results in a very fast computational algorithm capable of carrying out the vast number of simulation computations needed in guidance, stability, and control studies. The three-dimensional fore- and aft body (nozzle) geometry is characterized by the centerline profiles as represented by a series of coordinate points and body cross-section curvature. The engine module geometry is represented by an adjustable vertical grid to accommodate variations of the field parameters throughout the inlet and combustor. The scramjet inlet is modeled as a two-dimensional supersonic flow containing adjustable sidewall wedges and multiple fuel injection struts. The inlet geometry including the sidewall wedge angles, the number of injection struts, their sweepback relative to the vehicle reference line, and strut cross-section are user selectable. Combustion is currently represented by a Rayleigh line calculation including corrections for variable gas properties; improved models are being developed for this important element of the propulsion flow field. The program generates (1) variation of thrust magnitude and direction with angle of attack, (2) pitching moment and line of action of the thrust vector, (3) pressure and temperature distributions throughout the system, and (4) performance parameters such as thrust coefficient, specific impulse, mass flow rates, and equivalence ratio. Preliminary results are in good agreement with available performance data for systems resembling the NASP vehicle configuration.

Flandro, G. A.↗
NASA NTRS2009

Detection of Cometary Amines in Samples Returned by the Stardust Spacecraft

The delivery of amino acids to the early Earth by comets and their fragments could have been a significant source of the early Earth's prebiotic organic inventory that led to the emergence of life (Chyba and Sagan, 1992). Over 20 organic molecules including methane, ethane, ammonia, cyanic acid, formaldehyde, formamide, acetaldehyde, acetonitrile, and methanol have been identified by radio spectroscopic observations of the comae of comets Hale-Bopp and Hyakutake (Crovisier et al. 2004). These simple molecules could have provided the organic reservoir to allow the formation of more complex prebiotic organic compounds such as amino acids. After a 7-year mission, the Stardust spacecraft returned to Earth samples from comet Wild 2 on January 15, 2006 providing the opportunity to analyze the organic composition and isotopic distribution of cometary material with state-of-the-art laboratory instrumentation. The Preliminary Examination Team analyses of organics in samples returned by Stardust were largely focused on particles that impacted the collector aerogel and aluminum foil (Sandford et al. 2006). However, it is also possible that Stardust returned a "diffuse" sample of gas-phase organic molecules that struck the aerogel directly or diffused away from the grains after impact. To test this possibility, samples of Stardust flight aerogel and foil were carried through a hot water extraction and acid hydrolysis procedure to see if primary amine compounds were present in excess of those seen in controls. Here we report highly sensitive liquid chromatography time-of-flight mass spectrometry measurements of amino acids and amines in samples returned from a comet (Glavin et al. 2008). A suite of amino acids and amines including glycine, L-alanine, methylamine (MA), and ethylamine (EA) were identified in the Stardust bulk aerogel. With the exception of MA and EA, all other primary amines detected in comet-exposed aerogels were also present in the aerogel witness tile that was not exposed to Wild 2, suggesting that most amines are terrestrial in origin. However, the enhanced abundances of MA, EA, and possibly glycine in comet-exposed aerogel compared to controls, coupled with MA to EA ratios (1 to 2) that are distinct from preflight aerogels (7 to 10), suggest that these amines were captured from Wild 2. It is possible that MA and EA were formed on energetically processed icy grains containing methane, ethane, and ammonia. The presence of cometary amines in Stardust material supports the hypothesis that comets were an important source of prebiotic organics on the early Earth. To better understand their origin, a systematic compound specific carbon isotopic analysis (C-CSIA) via gas chromatography quadrupole mass spectrometry in with parallel with combustion isotope ratio mass spectrometry (GCQMS/ IRMS) is being conducted. We will discuss our latest C-CSIA measurements and what they indicate about the origin of amino acids extracted from Stardust samples.

Glavin, Daniel P.↗
DOE OSTI2022

An economics-by-design approach to a radiant integrated thermophotovoltaic microreactor system

This paper presents an application of the economics-by-design approach to the Radiant Integrated Thermo-photovoltaic Microreactor System (RITMS). The RITMS design is unique in that it directly couples a critical fission reactor with thermo-photovoltaic (TPV) panels for high efficiency energy conversion. This significant shift from electric conversion using traditional dynamic heat pumps leads to a simpler and more reliable system without turbomachinery and high pressurization. In working towards wrapping up the early design work, the economics-by-design approach, which centers economic competitiveness as the optimization parameter, is well suited in making a final determination on viability. This paper describes the computational sequence that was developed to couple the radiative and conductive heat transfer and feed operational performance to cost estimation. The framework was applied to maximize the power of the system, while minimizing the fuel enrichment. This method is applied to a reference RITMS design as part of a parametric sensitivity study, which revealed that the system is under moderated and that single unit plants could produce power as low as 300 dollar/MWh. This price point supports the notion that early in the RITMS design implementation, adoption into niche markets as a first-of-a-kind technology is possible. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2010

Detection of Cometary Amines in Samples Returned by the Stardust Spacecraft

The delivery of amino acids to the early Earth by comets and their fragments could have been a significant source of the early Earth's prebiotic organic inventory that led to the emergence of life (Chyba and Sagan, 1992). Over 20 organic molecules including methane, ethane, ammonia, cyanic acid, formaldehyde, formamide, acetaldehyde, acetonitrile, and methanol have been identified by radio spectroscopic observations of the comae of comets Hale-Bopp and Hyakutake (Crovisier et al. 2004). These simple molecules could have provided the organic reservoir to allow the forn1ation of more complex prebiotic organic compounds such as amino acids. After a 7 -year mission, the Stardust spacecraft returned to Earth samples from comet Wild 2 on January 15, 2006 providing the opportunity to analyze the organic composition and isotopic distribution of cometary material with state-of-the-art laboratory instrun1entation. The Preliminary Examination Team analyses of organics in samples returned by Stardust were largely focused on particles that impacted the collector aerogel and aluminum foil (Sandford et al. 2006). However, it is also possible that Stardust returned a "diffuse" sample of gas-phase organic molecules that struck the aerogel directly or diffused away from the grains after impact. To test this possibility, samples of Stardust flight aerogel and foil were carried through a hot water extraction and acid hydrolysis procedure to see if primary amine compounds were present in excess of those seen in controls. Here we report highly sensitive liquid chromatography time-of-flight mass spectron1etry measurements of amino acids and amines in samples returned from a comet (Glavin et al. 2008). A suite of amino acids and amines including glycine, L-alanine, methylamine (MA), and ethylamine (EA) were identified in the Stardust bulk aerogel. With the exception of MA and EA, all other primary amines detected in comet-exposed aerogels were also present in the aerogel witness tile that was not exposed to Wild 2, suggesting that most amines are terrestrial in origin. However, the enhanced abundances of MA, EA, and possibly glycine in comet-exposed aerogel compared to controls, coupled with MA to EA ratios (1 to 2) that are distinct from preflight aerogels (7 to 10), suggest that these amines were captured from Wild 2. It is possible that MA and EA were formed on energetically processed icy grains containing methane, ethane, and ammonia. The presence of cometary amines in Stardust material supports the hypothesis that comets were an important source of pre biotic organics on the early Earth. To better understand their origin, a systematic compound specific carbon isotopic analysis (C-CSIA) via gas chromatography quadrupole mass spectrometry in with parallel with combustion isotope ratio mass spectrometry (GC-QMS/IRMS) is being conducted. We will discuss our latest C-CSIA measurements and what they indicate about the origin of amino acids extracted from Stardust samples.

Glavin, Daniel P.↗
DOE OSTI2026

General framework for quantifying dissipation pathways in open quantum systems. III. Off-diagonal subsystem–bath couplings

This paper extends the previously reported theory of dissipation pathways [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111 (2024)] to incorporate off-diagonal subsystem–bath coupling, which is often required to model molecular systems where the environment directly influences transitions and couplings between subsystem states. We systematically derive master equations for both population transfer and dissipation into individual bath components, for which we also rigorously prove energy conservation and detailed balance. The approach is based on second-order perturbation theory with respect to the subsystem–bath couplings, whose form is not limited to any specific model. The accuracy of the developed method is tested by applying it to diverse model Hamiltonians involving linearly coupled harmonic oscillator baths and comparing the outcomes against the hierarchical equations of motion (HEOM) method. Overall, our method accurately quantifies the contributions of specific bath components to the overall dissipation while significantly reducing the computational cost compared to numerically exact methods such as HEOM, thus offering a path to examine how vibronic interactions steer non-adiabatic processes in realistic chemical systems.

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

Effective field theory of Stückelberg vector bosons

We explore the effective field theory of a vector field $X^μ$ that has a Stückelberg mass. The absence of a gauge symmetry for $X^μ$ implies Lorentz-invariant operators are constructed directly from Xμ. Beyond the kinetic and mass terms, allowed interactions at the renormalizable level include $X_μX^μH^†H, (X_μX^μ)^2$, and $X_μj^μ$, where $j^μ$ is a global current of the SM or of a hidden sector. We show that all of these interactions lead to scattering amplitudes that grow with powers of $\sqrt{s}/m_X$, except for the case of $X_μj^μ$ where $j^μ$ is a nonanomalous global current. The latter is well-known when $\textit{X}$ is identified as a dark photon coupled to the electromagnetic current, often written equivalently as kinetic mixing between $\textit{X}$ and the photon. The power counting for the energy growth of the scattering amplitudes is facilitated by isolating the longitudinal enhancement. We examine in detail the interaction with an anomalous global vector current $X_μj^μ_{\text{anom}}$, carefully isolating the finite contribution to the fermion triangle diagram. We calculate the longitudinally-enhanced observables $Z → X_γ$ (when $m_X < m_Z), f\bar{f} → X_γ$, and $Z_γ → Z_γ$ when $\textit{X}$ couples to the baryon number current. Introducing a "fake" gauge-invariance by writing $X^μ = A^μ – ∂^μπ/m_X$, the would-be gauge anomaly associated with $A_μj^μ_{\text{anom}}$ is canceled by $j^μ_{\text{anom}}∂_μπ/m_X$; this is the four-dimensional Green-Schwarz anomaly-cancellation mechanism at work. Our analysis demonstrates there is a much larger set of possible interactions that an EFT with a Stückelberg vector field can have, revealing scattering amplitudes that grow with energy. The growth of these amplitudes can be tamed by a dark Higgs sector, but this requires dark Higgs boson interactions (and reintroduces fine-tuning in the dark Higgs sector) that can be separated from X interactions only in the limit $g \ll 1$ .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRSDate not supplied

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗
DOE OSTI2024

Rapid dissolution without elemental fractionation by laser driven hydrothermal processing

Traditional dissolution of geologic samples often requires a significant time investment. Here, we present an alternative method for the dissolution of geologic materials using laser-driven hydrothermal processing (LDHP). LDHP uses laser energy directed onto a submerged sample, which increases the temperature and pressure at the liquid–sample interface and drives the hydrothermal dissolution coupled with photomechanical spallation, an ablative process. This uses focused 527 nm laser energy at 40 W average power, 1 kHz pulse repetition rate, and 115 ns pulse duration. Importantly, when LDHP is performed on basalt geostandards (BCR-2 and BHVO-2) using the conditions outlined, we show that LDHP does not produce significant elemental fractionation and, thus, can be considered an alternative processing method to traditional mechanical crushing and acid digestion. Additionally, it is possible using LDHP to utilize the spatially confined beam to target and selectively isolate individual phases in a rock, potentially alleviating the need for mechanical separation of inclusions that are difficult to physically isolate. Furthermore, using this outlined method of LDHP, we demonstrate full dissolution of 120 mg of obsidian in 85 minu, meaning that LDHP is a potentially very useful method when sample processing is time sensitive.

58 GEOSCIENCES↗
DOE OSTI2022

A Realistic Full-Scale 3D Modeling of Turning Using Coupled Smoothed Particle Hydrodynamics and Finite Element Method for Predicting Cutting Forces

Computational modelling is an effective technique for understanding the complex physics of machining. Large deformations, material separation, and high computational requirements are the key challenges faced while simulating machining. This work introduces a full-scale three-dimensional model of turning operations using a combined approach based on the Smoothed Particle Hydrodynamics (SPH) and Finite Element (FE) methods. By exploiting the advantages of each method, this approach leads to high-fidelity coupled SPH-FE machining models. Cutting forces and chip morphology are the primary results of interest. The machining models are validated with the results of turning experiments. Two-dimensional machining model underpredicts the cutting force and feed force by approximately 49% and 70%, respectively. Moreover, passive force cannot be predicted using the two-dimensional model. On the other hand, with the three-dimensional models developed in this manuscript, the difference between the total simulated force and experimentally measured force is ∼17%. The chip morphologies correlate with experiments in terms of the direction of the chip movement and the “long” continuous chips observed while turning Al 6061. This work expands the realm of machining simulations from two-dimensional orthogonal machining or sectional three-dimensional model to a full-scale realistic simulation. The encouraging simulation results show the potential to study more complex phenomena, such as machining stability and tool path modulation.

42 ENGINEERING↗
DOE OSTI2022

Charting C–C coupling pathways in electrochemical CO 2 reduction on Cu(111) using embedded correlated wavefunction theory

The electrochemical CO 2 reduction reaction (CO 2 RR) powered by excess zero-carbon-emission electricity to produce especially multicarbon (C 2+ ) products could contribute to a carbon-neutral to carbon-negative economy. Foundational to the rational design of efficient, selective CO 2 RR electrocatalysts is mechanistic analysis of the best metal catalyst thus far identified, namely, copper (Cu), via quantum mechanical computations to complement experiments. Here, we apply embedded correlated wavefunction (ECW) theory, which regionally corrects the electron exchange-correlation error in density functional theory (DFT) approximations, to examine multiple C–C coupling steps involving adsorbed CO (*CO) and its hydrogenated derivatives on the most ubiquitous facet, Cu(111). We predict that two adsorbed hydrogenated CO species, either *COH or *CHO, are necessary precursors for C–C bond formation. The three kinetically feasible pathways involving these species yield all three possible products: *COH–CHO, *COH–*COH, and *OCH–*OCH. The most kinetically favorable path forms *COH–CHO. In contrast, standard DFT approximations arrive at qualitatively different conclusions, namely, that only *CO and *COH will prevail on the surface and their C–C coupling paths produce only *COH–*COH and *CO–*CO, with a preference for the first product. This work demonstrates the importance of applying qualitatively and quantitatively accurate quantum mechanical method to simulate electrochemistry in order ultimately to shed light on ways to enhance selectivity toward C 2+ product formation via CO 2 RR electrocatalysts.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Evidence of cryptic methane cycling and non-methanogenic methylamine consumption in the sulfate-reducing zone of sediment in the Santa Barbara Basin, California

The recently discovered cryptic methane cycle in the sulfate-reducing zone of marine and wetland sediment couples methylotrophic methanogenesis to anaerobic oxidation of methane (AOM). Here we present evidence of cryptic methane cycling activity within the upper regions of the sulfate-reducing zone, along a depth transect within the Santa Barbara Basin, off the coast of California, USA. The top 0–20 cm of sediment from each station was subjected to geochemical analyses and radiotracer incubations using 35 S–SO$^{2-}_{4}$, 14 C–mono-methylamine, and 14 C–CH 4 to find evidence of cryptic methane cycling. Methane concentrations were consistently low (3 to 16 µM) across the depth transect, despite AOM rates increasing with decreasing water depth (from max 0.05 nmol cm -3 d -1 at the deepest station to max 1.8 nmol cm -3 d -1 at the shallowest station). Porewater sulfate concentrations remained high (23 to 29 mM), despite the detection of sulfate reduction activity from 35 S–SO$^{2-}_{4}$ incubations with rates up to 134 nmol cm -3 d -1 . Metabolomic analysis showed that substrates for methanogenesis (i.e., acetate, methanol and methylamines) were mostly below the detection limit in the porewater, but some samples from the 1–2 cm depth section showed non-quantifiable evidence of these substrates, indicating their rapid turnover. Estimated methanogenesis from mono-methylamine ranged from 0.2 to 0.5 nmol cm -3 d -1 . Discrepancies between the rate constants (k) of methanogenesis (from 14 C–mono-methylamine) and AOM (from either 14 C–mono-methylamine-derived 14 C–CH 4 or from directly injected 14 C–CH 4 ) suggest the activity of a separate, concurrent metabolic process directly metabolizing mono-methylamine to inorganic carbon. We conclude that the results presented in this work show strong evidence of cryptic methane cycling occurring within the top 20 cm of sediment in the Santa Barbara Basin. The rapid cycling of carbon between methanogenesis and methanotropy likely prevents major build-up of methane in the sulfate-reducing zone. Furthermore, our data suggest that methylamine is utilized by both methanogenic archaea capable of methylotrophic methanogenesis and non-methanogenic microbial groups. We hypothesize that sulfate reduction is responsible for the additional methylamine turnover, but further investigation is needed to elucidate this metabolic activity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Review—In Situ X-ray and Infrared Spectroscopic Studies of Electrochemical Systems

Despite of intense research and a wealth of data, the phenomena occurring during electrocatalysis are still a major obstacle in many chemical processes. Molecular analysis of the electrode/electrolyte interface is needed to correctly describe the reaction through identifying the species involved, their interaction with the environment and kinetics in situ, i.e. while the reaction is taking place. That can be done by coupling the electrochemical system with complementary non-electrochemical techniques. Particularly revealing are in situ X-ray spectroscopic techniques to analyze the electrode itself, providing the information on the changes in the catalyst during the reaction. The synergy of the traditional electrochemical techniques with the complementary spectroscopic methodologies offer understanding of the electrode/electrolyte interface above and beyond traditional experimental mainframe. Here we demonstrate how in situ X-ray absorption spectroscopy (XAS), in situ infrared reflection/absorption spectroscopy (IRRAS), and traditional voltammetric studies can increase our understanding of electrochemical processes during oxidation of ethanol. The results show the pronounced role of electrode surface in determining reaction kineticks and revealed the selectivity of the catalyst to complete oxidation pathway. They further provide understanding of the parameters that enhance its oxidation for future designing catalysts for alcohol oxidation fuel cells.

30 DIRECT ENERGY CONVERSION↗