Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural model decomposition”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Predictive Engineering Tools for Injection-Molded Long-Carbon-Fiber Thermoplastic Composites (FY 2015 Third Quarterly Report)

During the third quarter of FY 2015, the following technical progress has been made toward project milestones: 1) Magna oversaw the tool build and prepared the molding plan for the complex part of Phase II. 2) PlastiComp hosted a visit by Magna and Toyota on April 23rd to finalize the molding scope and schedule. The plan for molding trials including selection of molding parameters for both LFT and D-LFT for the U-shape complex part was established. 3) Toyota shipped the U-shape complex part tool to Magna on May 28th, 2015. 4) Plasticomp provided 30wt% LCF/PP and 30wt% LCF/PA66 compounded pellets to Magna for molding the complex part. 5) Magna performed preliminary molding trials on June 2nd, 2015 to validate wall thickness, fill profile, tool temperature and shot size requirements for the complex part. 6) Magna performed the first complex part run on June 16th and 17th, 2015 at Magna’s Composite Centre of Excellence in Concord, ON, Canada. Dale Roland of Plasticomp, and Umesh Gandhi of Toyota also attended the molding. 7) Magna discussed and finalized the plan with PNNL and the team for cutting samples from molded parts at selected locations for fiber orientation and length measurements. 8) Magna provided the computer-aided design (CAD) files of the complex parts with and without ribs to PNNL and Autodesk to build the corresponding ASMI models for injection molding simulations. Magna also provided the actual parameters used. 9) Plasticomp’s provided knowledge and experience of molding LCF materials essential to the successful molding of the parts including optimization of fill speed, tool temperatures, and plasticizing conditions for the 30wt% LCF/PP and 30wt% LCF/PA66 materials in both rib and non-rib versions. 10) Magna molded additional parts for evaluation of mechanical property testing including torsional stiffness on June 29th and 30th, 2015 at Magna’s Composite Center of Excellence. 11) Toyota began preparation for the torsion test of the specimens. Preparation of a computer-aided engineering (CAE) model to predict the performance is in progress. 12) Autodesk fixed an error in the implementation of the proper orthogonal decomposition (POD) calculation of fiber length that had caused the ASMI solution to crash and provided an updated build of ASMI containing the fix. 13)Autodesk reviewed and provided feedback for the complex part molding and measurement locations. 14) Autodesk provided support to set up the workflow for ASMI-ABAQUS® analysis, and provided a fix and workaround for a bug in the ASMI-ABAQUS® output command. 15) Autodesk helped build ASMI analysis models for the complex parts with and without ribs. 16) Autodesk worked on improving the orientation prediction accuracy in the shearing layer for 3D meshes based on comparison to measured data of the plaque moldings. 17) PNNL installed a new ASMI version received from Autodesk and performed comparative analyses to assess mid-plane versus 3D fiber length predictions using the full fiber length model and the reduced-order model (ROM) using POD. 18) PNNL presented the project scope, accomplishments, significant results and future plans to DOE and the USCAR Materials Tech Team on June 3rd, 2015. 19) PNNL discussed the cutting of samples from molded parts and finalized a plan with Magna and the team suggesting the sample size, locations and number of samples per location. 20) PNNL and Autodesk built ASMI models for the complex parts with and without ribs, and preliminary analyses of the part with ribs were conducted using the actual molding parameters received from Magna. 21) PNNL worked on a procedure to extract fiber orientation and length results from a 3D ASMI analysis to a 3D ABAQUS model. This procedure is essential to import ASMI fiber orientation and length to a 3D ABAQUS model of the part allowing future part structural analysis for weight reduction study.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Pyrolysis Pathways of the Furanic Ether 2-Methoxyfuran

Substituted furans, including furanic ethers, derived from nonedible biomass have been proposed as second-generation biofuels. In order to use these molecules as fuels, it is important to understand how they break apart thermally. In this work, a series of experiments were conducted to study the unimolecular and low-pressure bimolecular decomposition mechanisms of the smallest furanic ether, 2-methoxyfuran. Electronic structure (CBS-QB3) calculations indicate this substituted furan has an unusually weak O - CH 3 bond, approximately 190 kJ mol -1 (45 kcal mol -1 ); thus, the primary decomposition pathway is through bond scission resulting in CH 3 and 2-furanyloxy (O - C 4 H 3 O) radicals. Final products from the ring opening of the furanyloxy radical include 2 CO, HC$\equiv$CH, and H. The decomposition of methoxyfuran is studied over a range of concentrations (0.0025-0.1%) in helium or argon in a heated silicon carbide (SiC) microtubular flow reactor (0.66-1 mm i.d., 2.5-3.5 cm long) with reactor wall temperatures from 300 to 1300 K. Inlet pressures to the reactor are 150-1500 Torr, and the gas mixture emerges as a skimmed molecular beam at a pressure of approximately 10 μTorr. Products formed at early pyrolysis times (100 μs) are detected by 118.2 nm (10.487 eV) photoionization mass spectrometry (PIMS), tunable synchrotron VUV PIMS, and matrix infrared absorption spectroscopy. Secondary products resulting from H or CH 3 addition to the parent and reaction with 2-furanyloxy were also observed and include CH 2 =CH-CHO, CH 3 -CH=CH-CHO, CH 3 =CO-CH=CH 2 , and furanones; under the conditions in the reactor, we estimate these reactions contribute to at most 1-3% of total methoxyfuran decomposition. This work also includes observation and characterization of an allylic lactone radical, 2-furanyloxy (O-C 4 H 3 O), with the assignment of several intense vibrational bands in an Ar matrix, an estimate of the ionization threshold, and photoionization efficiency. A pressure-dependent kinetic mechanism is also developed to model the decomposition behavior of methoxyfuran and provide pathways for the minor bimolecular reaction channels that are observed experimentally.

09 BIOMASS FUELS↗
DOE OSTI2019

Support for Next-Generation Ecosystem Experiments (NGEE Arctic) Field Campaign Report

An important challenge for Earth System Models (ESMs) is to represent land surface and subsurface processes and their complex interactions in a warming climate. This is especially important for arctic ecosystems where permafrost extent, topography, hydrology, vegetation, disturbance, and biogeochemistry are inextricably linked. The implications of such linkages include permafrost thaw and deepening of the active layer, microbial decomposition of vulnerable soil organic matter, altered productivity and migration of tall woody shrubs, and watershed-scale changes in surface and groundwater transport and storage. Although ESMs describe some of these interactions for high-latitude ecosystems, their representation requires extensive confrontation with field and laboratory observations to test and improve models, and to use those models to inspire new observations and experiments. The Next-Generation Ecosystem Experiments (NGEE Arctic) is a 10-year project (2012 to 2022) to improve our predictive understanding of carbon-rich arctic system processes and feedbacks to climate. This is achieved through experiments, observations, and synthesis of existing data sets that strategically inform model process representation and parameterization, and that enhance the knowledge base required for model initialization, calibration, and evaluation. One question of special interest to the NGEE Arctic project addresses how above- and below-ground plant functional traits might change across environmental gradients, and what are the consequences for arctic ecosystem C, water, and nutrient fluxes? Arctic plant traits, and their variation in response to changing environmental conditions, will play a key role in the response of tundra ecosystems to warming, permafrost thaw, and the wetter or drier conditions expected in the future. The appropriate representation of these functional traits in models is necessary to accurately represent ecosystem C, water, and nutrient cycling in tundra ecosystems, now and in the future. We characterized the variation in above- and below-ground plant traits in response to varying edaphic and environmental conditions in Utqiaġvik (formerly Barrow), Alaska. These new data, combined with remote-sensing and synthesis activities across the Arctic and the globe, have been used to inform model structure and parameterization of key processes such as photosynthesis, nutrient dynamics, and trait variation across the landscape.

North Slope of Alaska, arctic processes (meteorolo↗
DOE OSTI2017

The fate of the tert-butyl radical in low-temperature autoignition reactions

Alkyl combustion models depend on kinetic parameters derived from reliable experimental or theoretical energetics that are often unavailable for larger species. Here, we have performed a comprehensive investigation of the tert -butyl radical ( R • in this paper) autoignition pathways. CCSD(T)/ANO0 geometries and harmonic vibrational frequencies were obtained for key stationary points for the R • + O 2 and QOOH + O 2 mechanisms. Relative energies were computed to chemical accuracy (±1 kcal mol -1 ) via extrapolation of RCCSD(T) energies to the complete basis-set limit, or usage of RCCSD(T)-F12 methods. At 0 K, the minimum energy R • + O 2 pathway involves direct elimination of HO $^{•}_{2}$ (30.3 kcal mol -1 barrier) from the tert -butyl peroxy radical ( ROO • ) to give isobutene. This pathway lies well below the competing QOOH -forming intramolecular hydrogen abstraction pathway (36.2 kcal mol -1 barrier) and ROO • dissociation (35.9 kcal mol -1 barrier). The most favorable decomposition channel for QOOH radicals leads to isobutene oxide (12.0 kcal mol -1 barrier) over isobutene (18.6 kcal mol -1 barrier). For the QOOH + O 2 pathways, we studied the transition states and initial products along three pathways: (1) α -hydrogen abstraction (42.0 kcal mol -1 barrier), (2) γ -hydrogen abstraction (27.0 kcal mol -1 barrier), and (3) hydrogen transfer to the peroxy moiety (24.4 kcal mol -1 barrier). The barrier is an extensive modification to the previous 18.7 kcal mol -1 value and warrants further study. However, it is still likely that the lowest energy QOOH + O 2 pathway corresponds to pathway (3). We found significant spin contamination and/or multireference character in multiple stationary points, especially for transition states stemming from QOOH . Lastly, we provide evidence for an $\tilde{A} - \tilde{X}$ surface crossing at a C s -symmetric, intramolecular hydrogen abstraction structure.

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

The Cementitious Barriers Partnership Experimental Programs and Software Advancing DOE’s Waste Disposal/Tank Closure Efforts – 15436

The U.S. Department of Energy Environmental Management (DOE-EM) Office of Tank Waste Management-sponsored Cementitious Barriers Partnership (CBP) is chartered with providing the technical basis for implementing cement-based waste forms and radioactive waste containment structures for long-term disposal. DOE needs in this area include the following to support progress in final treatment and disposal of legacy waste and closure of High-Level Waste (HLW) tanks in the DOE complex: long-term performance predictions, flow sheet development and flow sheet enhancements, and conceptual designs for new disposal facilities. The DOE-EM Cementitious Barriers Partnership is producing software and experimental programs resulting in new methods and data needed for end-users involved with environmental cleanup and waste disposal. Both the modeling tools and the experimental data have already benefited the DOE sites in the areas of performance assessments by increasing confidence backed up with modeling support, leaching methods, and transport properties developed for actual DOE materials. In 2014, the CBP Partnership released the CBP Software Toolbox –“Version 2.0” which provides concrete degradation models for 1) sulfate attack, 2) carbonation, and 3) chloride initiated rebar corrosion, and includes constituent leaching. These models are applicable and can be used by both DOE and the Nuclear Regulatory Commission (NRC) for service life and long-term performance evaluations and predictions of nuclear and radioactive waste containment structures across the DOE complex, including future SRS Saltstone and HLW tank performance assessments and special analyses, Hanford site HLW tank closure projects and other projects in which cementitious barriers are required, the Advanced Simulation Capability for Environmental Management (ASCEM) project which requires source terms from cementitious containment structures as input to their flow simulations, regulatory reviews of DOE performance assessments, and Nuclear Regulatory Commission reviews of commercial nuclear power plant (NPP) structures which are part of the overall US Energy Security program to extend the service life of NPPs. In addition, the CBP experimental programs have had a significant impact on the DOE complex by providing specific data unique to DOE sodium salt wastes at Hanford and SRS which are not readily available in the literature. Two recent experimental programs on cementitious phase characterization and on technetium (Tc) mobility have provided significant conclusions as follows: recent mineralogy characterization discussed in this paper illustrates that sodium salt waste form matrices are somewhat similar to but not the same as those found in blended cement matrices which to date have been used in long-term thermodynamic modeling and contaminant sequestration as a first approximation. Utilizing the CBP generated data in long-term performance predictions provides for a more defensible technical basis in performance evaluations. In addition, recent experimental studies related to technetium mobility indicate that conventional leaching protocols may not be conservative for direct disposal of Tc-containing waste forms in vadose zone environments. These results have the potential to influence the current Hanford supplemental waste treatment flow sheet and disposal conceptual design.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Towards a self-consistent analysis of the anisotropic galaxy two- and three-point correlation functions on large scales: application to mock galaxy catalogues

We establish a practical method for the joint analysis of anisotropic galaxy two- and three-point correlation functions (2PCF and 3PCF, respectively) on the basis of the decomposition formalism of the 3PCF using tripolar spherical harmonics. We perform such an analysis with MultiDark-Patchy mock catalogues to demonstrate and understand the benefit of the anisotropic 3PCF. We focus on scales above $80\, h^{-1}\, {\rm Mpc}$, and use information from the shape and the baryon acoustic oscillation (BAO) signals of the 2PCF and 3PCF. We also apply density field reconstruction to increase the signal-to-noise ratio of BAO in the 2PCF measurement, but not in the 3PCF measurement. In particular, we study in detail the constraints on the angular diameter distance and the Hubble parameter. We build a model of the bispectrum or 3PCF that includes the non-linear damping of the BAO signal in redshift space. We carefully account for various uncertainties in our analysis including theoretical models of the 3PCF, window function corrections, biases in estimated parameters from the fiducial values, the number of mock realizations to estimate the covariance matrix, and bin size. The joint analysis of the 2PCF and 3PCF monopole and quadrupole components shows a $30{{\ \rm per\ cent}}$ and $20{{\ \rm per\ cent}}$ improvement in Hubble parameter constraints before and after reconstruction of the 2PCF measurements, respectively, compared to the 2PCF analysis alone. Furthermore, this study clearly shows that the anisotropic 3PCF increases cosmological information from galaxy surveys and encourages further development of the modelling of the 3PCF on smaller scales than we consider.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Primordial black hole dark matter: A quantitative parameter sensitivity comparison across formation mechanisms and particle candidates

Primordial black holes (PBHs) in the asteroid-mass window ( 10 17 – 10 22 g ) can account for all of the dark matter without violating any observational constraint, yet are routinely dismissed as fine-tuned. I put that dismissal to the test by applying three complementary sensitivity measures uniformly across a broad landscape: three noninflationary PBH production mechanisms, six classes of inflationary PBH models, and seven particle dark matter benchmarks, all evaluated against the same observable target. Three distinct naturalness universality classes emerge, determined entirely by the analytic structure of the abundance map rather than by the nature of the dark matter candidate. Biased-domain-wall PBHs, in their least model-dependent (free- V b ) form, have the same low sensitivity, Δ = 4.5 , as off-resonance weakly interacting massive particles and freeze-in particles ( Δ = 2 ), a sensitivity that, because it is constant over the entire parameter space of the construction, also coincides trivially with its own Wilson-normalized average within that parameter space (Section Definition and conventions), an equivalence that concerns only the space over which Δ is computed and is not a naturalness statement about the construction as a whole; a further reduction to Δ = 2 is possible only under the additional, independently motivated but not required, assumption that the domain-wall bias is generated by Planck-suppressed operators; early matter-domination PBHs occupy an intermediate tier alongside coannihilating weakly interacting massive particles (WIMPs), unified by a structural identity in which the sensitivity measure equals the logarithm of the ratio of the formation scale to the matter–radiation equality scale; first-order phase transition PBHs, once the more accurate super-exponential collapse probability is used in place of the single-exponential approximation, instead belong to the same highly sensitive tier as resonant WIMP annihilation and single-field inflationary collapse, for a structurally distinct reason; single-field ultraslow-roll inflationary collapse is severely tuned for a distinct reason: a double exponential in which the power spectrum amplitude is itself exponentially sensitive to the inflaton potential coefficients, on top of the exponential collapse sensitivity of the abundance map. My main conclusion is that the claim that PBH dark matter is generically fine-tuned conflates the worst case with a landscape spanning every naturalness tier. The Barbieri-Giudice sensitivity computed here and the Wilson-normalized measure of Iovino and Riotto answer distinct and mutually consistent questions about the same construction, a distinction I clarify within the two-layer decomposition.

Profumo, Stefano [University of California, Santa ↗
DOE OSTI2025

Fused filament fabrication of thermoplastic polyurethane composites with microencapsulated phase-change material

Here, the present study examines the thermal energy storage (TES) effectiveness and printability of microencapsulated phase-change material (MEPCM) combined with thermoplastic polyurethane (TPU) for fused filament fabrication (FFF). Two formulations were assessed: 24D MEPCM, which changes phase at 24 ° C, compounded with TPU pellets and 43D MEPCM, which changes phase at 43 ° C, integrated with TPU powder. These combinations are designed to evaluate the effectiveness of the form of the TPU (pellets versus powder) in the FFF process. The investigation includes a comprehensive analysis of thermal characteristics, encompassing phase-change temperature, latent heat of fusion, thermal conductivity, and thermal decomposition, which are assessed through differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA). Additionally, mechanical properties, including stress-strain behavior, are examined to evaluate material suitability for TES applications, while microstructural visualization is used to provide deeper insights into material performance, structural integrity, and the quality of printed components. The 24D MEPCM and TPU pellets formulation experienced a significant loss of approximately 39.6% of PCM during filament extrusion and printing, resulting in a reduced effective latent heat. Therefore, further characterization of the pellet formulation was discontinued due to excessive leakage. In contrast, the 43D MEPCM and TPU powder formulation demonstrated minimal PCM loss, with the 60 wt.% composition achieving an effective latent heat of 132 J/g. This value represents the highest effective latent heat currently documented in the literature for PCM-polymer-composite materials produced using an FFF-based additive manufacturing process.

25 ENERGY STORAGE↗
NASA NTRS2011

NASA Tech Briefs, July 2011

Topics covered include: 1) Collaborative Clustering for Sensor Networks; 2) Teleoperated Marsupial Mobile Sensor Platform Pair for Telepresence Insertion Into Challenging Structures; 3) Automated Verification of Spatial Resolution in Remotely Sensed Imagery; 4) Electrical Connector Mechanical Seating Sensor; 5) In Situ Aerosol Detector; 6) Multi-Parameter Aerosol Scattering Sensor; 7) MOSFET Switching Circuit Protects Shape Memory Alloy Actuators; 8) Optimized FPGA Implementation of Multi-Rate FIR Filters Through Thread Decomposition; 9) Circuit for Communication Over Power Lines; 10) High-Efficiency Ka-Band Waveguide Two-Way Asymmetric Power Combiner; 11) 10-100 Gbps Offload NIC for WAN, NLR, and Grid Computing; 12) Pulsed Laser System to Simulate Effects of Cosmic Rays in Semiconductor Devices; 13) Flight Planning in the Cloud; 14) MPS Editor; 15) Object-Oriented Multi Disciplinary Design, Analysis, and Optimization Tool; 16) Cryogenic-Compatible Winchester Connector Mount and Retaining System for Composite Tubes; 17) Development of Position-Sensitive Magnetic Calorimeters for X-Ray Astronomy; 18) Planar Rotary Piezoelectric Motor Using Ultrasonic Horns; 19) Self-Rupturing Hermetic Valve; 20) Explosive Bolt Dual-Initiated from One Side; 21) Dampers for Stationary Labyrinth Seals; 22) Two-Arm Flexible Thermal Strap; 23) Carbon Dioxide Removal via Passive Thermal Approaches; 24) Polymer Electrolyte-Based Ambient Temperature Oxygen Microsensors for Environmental Monitoring; 25) Pressure Shell Approach to Integrated Environmental Protection; 26) Image Quality Indicator for Infrared Inspections; 27) Micro-Slit Collimators for X-Ray/Gamma-Ray Imaging; 28) Scatterometer-Calibrated Stability Verification Method; 29) Test Port for Fiber-Optic-Coupled Laser Altimeter; 30) Phase Retrieval System for Assessing Diamond Turning and Optical Surface Defects; 31) Laser Oscillator Incorporating a Wedged Polarization Rotator and a Porro Prism as Cavity Mirror; 32) Generic, Extensible, Configurable Push-Pull Framework for Large-Scale Science Missions; 33) Dynamic Loads Generation for Multi-Point Vibration Excitation Problems; 34) Optimal Control via Self-Generated Stochasticity; 35) Space-Time Localization of Plasma Turbulence Using Multiple Spacecraft Radio Links; 36) Surface Contact Model for Comets and Asteroids; 37) Dust Mitigation Vehicle; 38) Optical Coating Performance for Heat Reflectors of the JWST-ISIM Electronic Component; 39) SpaceCube Demonstration Platform; 40) Aperture Mask for Unambiguous Parity Determination in Long Wavelength Imagers; 41) Spaceflight Ka-Band High-Rate Radiation-Hard Modulator; 42) Enabling Disabled Persons to Gain Access to Digital Media; 43) Cytometer on a Chip; 44) Principles, Techniques, and Applications of Tissue Microfluidics; and 45) Two-Stage Winch for Kites and Tethered Balloons or Blimps.

Source record↗
DOE OSTI2023

Non-singlet quark helicity $\mathrm{PDF}$s of the nucleon from pseudo-distributions

The non-singlet helicity quark parton distribution functions (PDFs) of the nucleon are determined from lattice QCD, by jointly leveraging pseudo-distributions and the distillation spatial smearing paradigm. A Lorentz decomposition of appropriately isolated space-like matrix elements reveals pseudo-distributions that contain information on the leading-twist helicity PDFs, as well as an invariant amplitude that induces an additional z 2 contamination of the leading-twist signal. An analysis of the short-distance behavior of the space-like matrix elements using matching coefficients computed to next-to-leading order (NLO) exposes the desired PDF up to this additional z 2 contamination. Due to the non-conservation of the axial current, we elect to isolate the helicity PDFs normalized by the nucleon axial charge at the same scale μ 2 . The leading-twist helicity PDFs as well as several sources of systematic error, including higher-twist effects, discretization errors, and the aforementioned z 2 contaminating amplitude are jointly determined by characterizing the computed pseudo-distribution in a basis of Jacobi polynomials. The Akaike Information Criterion is exploited to effectively average over distinct model parameterizations and cuts on the pseudo-distribution. Encouraging agreement is observed with recent global analyses of each non-singlet quark helicity PDF, notably a rather small non-singlet anti-quark helicity PDF for all quark momentum fractions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Ecohydrologic processes and soil thickness feedbacks control limestone-weathering rates in a karst landscape

Here, chemical weathering of bedrock plays an essential role in the formation and evolution of Earth's critical zone. Over geologic time, the negative feedback between temperature and chemical weathering rates contributes to the regulation of Earth climate. The challenge of understanding weathering rates and the resulting evolution of critical zone structures lies in complicated interactions and feedbacks among environmental variables, local ecohydrologic processes, and soil thickness, the relative importance of which remains unresolved. We investigate these interactions using a reactive-transport kinetics model, focusing on a low-relief, wetland-dominated karst landscape (Big Cypress National Preserve, South Florida, USA) as a case study. Across a broad range of environmental variables, model simulations highlight primary controls of climate and soil biological respiration, where soil thickness both supplies and limits transport of biologically derived acidity. Consequently, the weathering rate maximum occurs at intermediate soil thickness. The value of the maximum weathering rate and the precise soil thickness at which it occurs depend on several environmental variables, including precipitation regime, soil inundation, vegetation characteristics, and rate of groundwater drainage. Simulations for environmental conditions specific to Big Cypress suggest that wetland depressions in this landscape began to form around beginning of the Holocene with gradual dissolution of limestone bedrock and attendant soil development, highlighting large influence of age-varying soil thickness on weathering rates and consequent landscape development. While climatic variables are often considered most important for chemical weathering, our results indicate that soil thickness and biotic activity are equally important. Weathering rates reflect complex interactions among soil thickness, climate, and local hydrologic and biotic processes, which jointly shape the supply and delivery of chemical reactants, and the resulting trajectories of critical zone and karst landscape development.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1998

Seasonal Cycles of Meridional Overturning and Heat Transport of the Indian Ocean

A general circulation model of the Indian Ocean is fitted to monthly averaged climatological temperatures, salinities, and surface fluxes using the adjoint method. Interannual variability is minimized by penalizing the temporal drift from one seasonal cycle to another during a two-year integration. The resultant meridional overturning and heat transport display large seasonal variations, with maximum amplitudes of 18 and 22 (x 10(exp 6) cubic m/s) for the overturning and 1.8 and 1.4 (x 10(exp 15) W) for heat transport near 10 S and 10 N, respectively. A dynamical decomposition of the overturning and heat transport shows that the time-varying Ekman How plus its barotropic compensation can explain a large part of the seasonal variations in overturning and heat transport. The maximum variations at 10 deg N and 10 deg S are associated with monsoon reversal over the northern Indian Ocean and changes of the easterlies over the southern Indian Ocean. An external mode with variable topography has a moderate contribution where the Somali Current and the corresponding gyre reverse direction seasonally. Contribution front vertical shear (thermal wind and ageostrophic shear) is dominant near the southern boundary and large near the Somali Current latitudes. The dominant balance in the zonally integrated heat budget is between heat storage change and heat transport convergence except south of 15 S. Optimization with seasonal forcings improves estimates of sea surface temperatures, but the annual average overturning and heat transport are very similar to previous results with annual mean forcings. The annual average heat transport consists of roughly equal contributions from time-mean and time-varying fields of meridional velocities and temperatures in the northern Indian Ocean. indicating a significant rectification to the heat transport due to the time-varying fields. The time-mean and time-varying contributions are primarily due to the overturning and horizontal gyre, respectively. Inclusion of TOPEX data enhances the seasonal cycles of the estimated overturning and heat transport in the central Indian Ocean significantly and improves the estimated equatorial zonal flows but leads 10 unrealistic estimates of the velocity structure near the Indonesian Throughflow region, most likely owing to the deficiencies in the lateral boundary conditions.

Lee, Tong↗
DOE OSTI2025

Sign Problem in Tensor-Network Contraction

We investigate how the computational difficulty of contracting tensor networks depends on the sign structure of the tensor entries. Using results from computational complexity, we observe that the approximate contraction of tensor networks with only positive entries has lower computational complexity as compared to tensor networks with general real or complex entries. This raises the question of how this transition in computational complexity manifests itself in the hardness of different tensor-network-contraction schemes. We pursue this question by studying random tensor networks with varying bias toward positive entries. First, we consider contraction via Monte Carlo sampling and find that the transition from hard to easy occurs when the tensor entries become predominantly positive; this can be understood as a tensor-network manifestation of the well-known negative-sign problem in quantum Monte Carlo. Second, we analyze the commonly used contraction based on boundary tensor networks. The performance of this scheme is governed by the number of correlations in contiguous parts of the tensor network (which by analogy can be thought of as entanglement). Remarkably, we find that the transition from hard to easy—i.e., from a volume-law to a boundary-law scaling of entanglement—already occurs for a slight bias of the tensor entries toward a positive mean, scaling inversely with the bond dimension D , and thus the problem becomes easy the earlier the larger D occurs. This is in contrast both to expectations and to the behavior found in Monte Carlo contraction, where the hardness at fixed bias increases with the bond dimension. To provide insight into this early breakdown of computational hardness and the accompanying entanglement transition, we construct an effective classical statistical-mechanical model that predicts a transition at a bias of the tensor entries of 1 / D , confirming our observations. We conclude by investigating the computational difficulty of computing expectation values of tensor-network wave functions (projected entangled-pair states, PEPSs) and find that in this setting, the complexity of entanglement-based contraction always remains low. We explain this by providing a local transformation that maps PEPS expectation values to a positive-valued tensor network. This not only provides insight into the origin of the observed boundary-law entanglement scaling but also suggests new approaches toward PEPS contraction based on positive decompositions. Published by the American Physical Society 2025

Chen, Jielun (ORCID:0000000178411545)↗
DOE OSTI2025

Calcium is associated with specific soil organic carbon decomposition products at Blodgett Forest Research Center, Georgetown, California as analysed with scanning transmission X-ray microscopy carbon near-edge X-ray absorption fine structure spectroscopy

This data is from the paper calcium is associated with specific soil organic carbon decomposition products, published in SOIL. DOI: https://doi.org/10.5194/soil-11-381-2025, 2025.This file contains CSVs with spectral data and bulk soil data and there is no specific program required to open this data. The data includes Scanning transmission X-ray microscopy carbon near-edge X-ray absorption fine structure spectroscopy. data from the measurement of samples from the Whole-soil Warming project, run by the Belowground Biogeochemistry team at Blodgett Forest Research Center, Georgetown, California run by the University of California, Berkeley. It also includes bulk soil chemical properties. The University of California's Blodgett Forest Research Station (Forest) is situated in the Sierra Nevada foothills (1370 m a.s.l.) near Georgetown, California. The samples were collected from here: 38.912013, -120.661469, https://maps.app.goo.gl/291bCJ1zVqUhgktz6. The Forest soils were characterised as Alfisols, which are equivalent to Dystric Cambisols (IUSS Working Group WRB, 2015), and formed in granitic parent materials, in a temperate climate, under thinned, mixed-coniferous forest (Fig. S3; Gaudinski et al., 2009). With these analyses we aimed to answer the question, is calcium associated with a specific type of organic matter enriched in aromatic and phenolic carbon at the microscale in samples from Blodgett Forest Research Center? and how does this specific type of carbon respond to experiments targetted at removing and adding calcium to the soils, specifically cation exchange and incubation after calcium addition? Abstract from the paper can be found below: Calcium (Ca) may contribute to the preservation of soil organic carbon (SOC) in more ecosystems than previously thought. Here we provide evidence that Ca is co-located with SOC compounds that are enriched in aromatic and phenolic groups, across different acidic soil-types and locations with different ecosystem properties, differing in terms of climate, parent material, soil type, and vegetation. In turn, this co-localised fraction of Ca-SOC is removed through cation-exchange, and the association is then only re-established during decomposition in the presence of Ca (Ca addition incubation). Thus, highlighting a causative link between decomposition and the co-location of Ca with a characteristic fraction of SOC. Decomposition increases the relative proportion of negatively charged functional groups, which can increase the propensity for the association between SOC and Ca, and in turn, this association inhibits dissolved organic carbon export or further decomposition. We propose that this mechanism could be driven by Ca hotspots on the microscale shifting local decomposition processes and thereby explaining the colocation of Ca with SOC of a specific composition across different acidic soil environments. Incorporating this biogeochemical process into Earth System Models could improve our understanding, predictions, and management of carbon dynamics in soils, and account for their response to Ca-rich amendments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Degradation of Adenine on the Martian Surface in the Presence of Perchlorates and Ionizing Radiation: A Reflectron Time-of-flight Mass Spectrometric Study

The aim of the present work is to unravel the radiolytic decomposition of adenine (C{sub 5}H{sub 5}N{sub 5}) under conditions relevant to the Martian surface. Being the fundamental building block of (deoxy)ribonucleic acids, the possibility of survival of this biomolecule on the Martian surface is of primary importance to the astrobiology community. Here, neat adenine and adenine–magnesium perchlorate mixtures were prepared and irradiated with energetic electrons that simulate the secondary electrons originating from the interaction of the galactic cosmic rays with the Martian surface. Perchlorates were added to the samples since they are abundant—and therefore relevant oxidizers on the surface of Mars—and they have been previously shown to facilitate the radiolysis of organics such as glycine. The degradation of the samples were monitored in situ via Fourier transformation infrared spectroscopy and the electron ionization quadruple mass spectrometric method; temperature-programmed desorption profiles were then collected by means of the state-of-the-art single photon photoionization reflectron time-of-flight mass spectrometry (PI-ReTOF-MS), allowing for the detection of the species subliming from the sample. The results showed that perchlorates do increase the destruction rate of adenine by opening alternative reaction channels, including the concurrent radiolysis/oxidation of the sample. This new pathway provides a plethora of different radiolysis products that were identified for the first time. These are carbon dioxide (CO{sub 2}), isocyanic acid (HNCO), isocyanate (OCN{sup −}), carbon monoxide (CO), and nitrogen monoxide (NO); an oxidation product containing carbonyl groups (R{sub 1}R{sub 2}–C=O) with a constrained five-membered cyclic structure could also be observed. Cyanamide (H{sub 2}N–C≡N) was detected in both irradiated samples as well.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Structure and formation of cD galaxies: NGC 6166 in Abell 2199

Hobby–Eberly Telescope (HET) spectroscopy is used to measure the velocity dispersion profile of the nearest prototypical cD galaxy, NGC 6166 in the cluster Abell 2199. We also present composite surface photometry from many telescopes. We confirm the defining feature of a cD galaxy; i.e., (we suggest), a halo of stars that fills the cluster center and that is controlled dynamically by cluster gravity, not by the central galaxy. Our HET spectroscopy shows that the velocity dispersion of NGC 6166 rises from σ≃300 km s{sup −1} in the inner r∼10′′ to σ=865±58 km s{sup −1} at r ∼ 100″ in the cD halo. This extends published observations of an outward σ increase and shows for the first time that σ rises all the way to the cluster velocity dispersion of 819 ± 32 km s{sup −1}. We also observe that the main body of NGC 6166 moves at +206 ± 39 km s{sup −1} with respect to the cluster mean velocity, but the velocity of the inner cD halo is ∼70 km s{sup −1} closer to the cluster velocity. These results support our picture that cD halos consist of stars that were stripped from individual cluster galaxies by fast tidal encounters. However, our photometry does not confirm the widespread view that cD halos are identifiable as an extra, low-surface-brightness component that is photometrically distinct from the inner, steep-Sérsic-function main body of an otherwise-normal giant elliptical galaxy. Instead, all of the brightness profile of NGC 6166 outside its core is described to ±0.037 V mag arcsec{sup −2} by a single Sérsic function with index n≃8.3. The cD halo is not recognizable from photometry alone. This blurs the distinction between cluster-dominated cD halos and the similarly-large-Sérsic-index halos of giant, core-boxy-nonrotating ellipticals. These halos are believed to be accreted onto compact, high-redshift progenitors (“red nuggets”) by large numbers of minor mergers. They belong dynamically to their central galaxies. Still, cDs and core-boxy-nonrotating Es may be more similar than we think: both may have outer halos made largely via minor mergers and the accumulation of tidal debris. We construct a main-body+cD-halo decomposition that fits both the brightness and dispersion profiles. To fit σ(r), we need to force the component Sérsic indices to be smaller than a minimum-χ{sup 2} photometric decomposition would suggest. The main body has M{sub V}≃−22.8≃30% of the total galaxy light. The cD halo has M{sub V}≃−23.7, ∼1/2 mag brighter than the brightest galaxy in the Virgo cluster. A mass model based on published cluster dynamics and X-ray observations fits our observations if the tangential dispersion is larger than the radial dispersion at r≃20′′ –60′′. The cD halo is as enhanced in α element abundances as the main body of NGC 6166. Quenching of star formation in ≲1 Gyr suggests that the center of Abell 2199 has been special for a long time during which dynamical evolution has liberated a large mass of now-intracluster stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Quantifying Rates of Erosion Using Historical Aerial Photos - 19478

Erosion can play an important role in the fate and transport of waste at sites where disposal of long-lived waste is anticipated. The statistical characterization of key processes related to erosion is essential to the understanding of site stability through time. In this work multi-temporal historical aerial images are analyzed in conjunction with orthophotography and Lidar data to develop probability distributions for variables representing important erosion processes. Scans of photographic negatives were obtained for multiple photographs and several negatives were scanned with a calibrated photogrammetric scanner by Geomni (formerly Blue Skies Consulting) in Belen, New Mexico. Other negatives were provided by USDA or scanned on a high-resolution, large format media scanner in the United States National Archives Cartographic Research Room, College Park, Maryland (1939). An additional set of aerial photographs was obtained from the USGS EROS Data Center web site. For other years, photographic prints were available for aerial photography. Several prints were scanned with a large-format media scanner at 1200 dpi. Ground control points for selected digital images used in the stereoscopic analyses were identified in ArcMap. Ground control points were collected by matching features which could be identified in the aerial photography as well as the orthophotography or shaded relief imagery from the Lidar dataset. Such features included fence line intersections, road intersections, railroad bridges, and structures, depending on the age of the photography. Many old fence lines and abandoned roads remained visible in the shaded relief imagery while they were completely obscured by tree canopy in modern photography. Ground control elevations were sampled from the Lidar-based DEM or directly from point clouds. The image coordinates and matching three-dimensional projected coordinates of the ground control points were exported to IDL where parameters for a three-dimensional direct linear transformation (DLT) were obtained. Built-in singular value decomposition routines were used to solve an overdetermined set of equations for the eleven parameters of the DLT for a particular historical photograph. A solution to the DLT parameters was obtained for each historical aerial photograph in which features of interest were identified using the stereoscope. Ultimately, this process resulted in the identification of locations for features like gully heads using both the Lidar dataset and projection of the estimated location from the historical aerial photo under consideration back to the Lidar dataset. The distance between the estimate of the location from the historical image and that of the Lidar dataset was the estimated distance the feature has moved. This distance was then divided by the number of years between the Lidar dataset and the year of the historical aerial image of interest to get a rate of movement through time. This was done for several dozen points on each historical aerial image. There was inevitable uncertainty associated with the estimation of points from the historical aerial imagery. This includes estimation of changes in locations of features like gully heads through time. In this work a method to estimate the uncertainty associated with measurements of locations of features using the stereoscopic analysis is developed and applied. For features with locations in space that do not change (e.g., road intersections, etc.) the difference in measurements of locations from the Lidar dataset and that of the estimate from the historical aerial image provides an observation of the measurement error. These points used to estimate the measurement error are referred to as control points. For each of the historical aerial images, a different measurement error distribution is estimated using control points. These measurement errors are combined with the estimated distances from the stereoscopic analyses to produce distributions of gully retreat rates and hillslope erosion rates used along with output from Erosion Working Group (EWG) Landscape Evolution Models (LEMs) to provide statistical probability distributions for input into the PPA Model. Specifically, gully head retreat rate and gully widening rate were characterized using statistical probability distributions. The central limit theorem was then applied to develop rates that represent long-term averages that are used to inform a PPA Model built in GoldSim. The importance of these rates with respect to characterization of future risk can then assessed using sensitivity analyses. This provides an approach to determine if subsequent reduction of the uncertainty associated with the probability distributions of erosion rates will correspond to meaningful reduction in the estimation of future risk from the site. Ultimately this is a general method that utilizes historical aerial imagery to estimate rates of erosion with uncertainty. This information can be used to inform PPA models to project future risks from a given site. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

Roofline Analysis in the Intel® Advisor to Deliver Optimized Performance for applications on Intel® Xeon Phi™ Processor

In this session we show, in two case studies, how the roofline feature of Intel Advisor has been utilized to optimize the performance of kernels of the XGC1 and PICSAR codes in preparation for Intel Knights Landing architecture. The impact of the implemented optimizations and the benefits of using the automatic roofline feature of Intel Advisor to study performance of large applications will be presented. This demonstrates an effective optimization strategy that has enabled these science applications to achieve up to 4.6 times speed-up and prepare for future exascale architectures. # Goal/Relevance of Session The roofline model [1,2] is a powerful tool for analyzing the performance of applications with respect to the theoretical peak achievable on a given computer architecture. It allows one to graphically represent the performance of an application in terms of operational intensity, i.e. the ratio of flops performed and bytes moved from memory in order to guide optimization efforts. Given the scale and complexity of modern science applications, it can often be a tedious task for the user to perform the analysis on the level of functions or loops to identify where performance gains can be made. With new Intel tools, it is now possible to automate this task, as well as base the estimates of peak performance on measurements rather than vendor specifications. The goal of this session is to demonstrate how the roofline feature of Intel Advisor can be used to balance memory vs. computation related optimization efforts and effectively identify performance bottlenecks. A series of typical optimization techniques: cache blocking, structure refactoring, data alignment, and vectorization illustrated by the kernel cases will be addressed. # Description of the codes ## XGC1 The XGC1 code [3] is a magnetic fusion Particle-In-Cell code that uses an unstructured mesh for its Poisson solver that allows it to accurately resolve the edge plasma of a magnetic fusion device. After recent optimizations to its collision kernel [4], most of the computing time is spent in the electron push (pushe) kernel, where these optimization efforts have been focused. The kernel code scaled well with MPI+OpenMP but had almost no automatic compiler vectorization, in part due to indirect memory addresses and in part due to low trip counts of low-level loops that would be candidates for vectorization. Particle blocking and sorting have been implemented to increase trip counts of low-level loops and improve memory locality, and OpenMP directives have been added to vectorize compute-intensive loops that were identified by Advisor. The optimizations have improved the performance of the pushe kernel 2x on Haswell processors and 1.7x on KNL. The KNL node-for-node performance has been brought to within 30% of a NERSC Cori phase I Haswell node and we expect to bridge this gap by reducing the memory footprint of compute intensive routines to improve cache reuse. ## PICSAR is a Fortran/Python high-performance Particle-In-Cell library targeting at MIC architectures first designed to be coupled with the PIC code WARP for the simulation of laser-matter interaction and particle accelerators. PICSAR also contains a FORTRAN stand-alone kernel for performance studies and benchmarks. A MPI domain decomposition is used between NUMA domains and a tile decomposition (cache-blocking) handled by OpenMP has been added for shared-memory parallelism and better cache management. The so-called current deposition and field gathering steps that compose the PIC time loop constitute major hotspots that have been rewritten to enable more efficient vectorization. Particle communications between tiles and MPI domain has been merged and parallelized. All considered, these improvements provide speedups of 3.1 for order 1 and 4.6 for order 3 interpolation shape factors on KNL configured in SNC4 quadrant flat mode. Performance is similar between a node of cori phase 1 and KNL at order 1 and better on KNL by a factor 1.6 at order 3 with the considered test case (homogeneous thermal plasma).

97 MATHEMATICS AND COMPUTING↗