SEARCH · Engineering Papers

Results for “linear response”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Impacts of year-to-year weather variability and inter-panel spacing on agrivoltaic crop yields in Massachusetts

Deployment of utility-scale solar power plants could lead to agricultural land-use changes. Agrivoltaics, a dual land use combining solar and agriculture on the same land, can provide multiple environmental benefits, including improving soil quality and water use efficiency. The body of agrivoltaic field data is still growing, and crop responses to different solar configurations under different local climates are highly varied. We investigate the impact of adding spacing between adjacent solar panels in a fixed-tilt system to improve light diffusion to crops. For four crops (broccoli, peppers, kale, Swiss chard) grown across 3 years in an agrivoltaic system in Massachusetts, we found that only kale had a linearly increasing trend as the inter-panel spacing increased from 0.6 m to 1.5 m (2 ft to 5 ft). However, there were significant year-to-year differences in the yield of agrivoltaic versus control fields. Agrivoltaic and full sun fields produced equivalent yields in a hot, dry year, whereas the full-sun control beds produced more salable yield for all four crops in a warm, wet year. This demonstrates variability of agricultural outcomes and the need for more multi-year studies to ensure agrivoltaic impacts are not under- or overestimated.

14 SOLAR ENERGY↗

DOE OSTI2020

Linking Hydrology and Dissolved Organic Matter Characteristics in a Subtropical Wetland: A Long-Term Study of the Florida Everglades

Dissolved organic matter (DOM) acts as an important biogeochemical component of aquatic ecosystems that controls nutrient cycling, influences water quality, and links terrestrial and oceanic carbon pools, yet long-term studies of how changing environmental drivers alter its abundance and composition are rare. Using a ten-year, spatially explicit dataset from Everglades National Park (ENP), a globally significant wetland, we investigated the relationships between DOM quality/quantity and hydrologic/climatic drivers along two contrasting marsh-estuarine transects based on generalized linear modeling and a cumulative sums analysis. Analyses revealed distinct spatial, seasonal, and interannual patterns in variability of DOC and optical properties. Landscape-scale seasonal patterns showed an enrichment in microbial-like and protein-like DOM during the dry season relative to the wet season. We found that, while some compositional constituents varied with the solar calendar, responsive to temperature and photoperiod, others varied with the hydrologic calendar. Independent water level and discharge effects indicated strong hydrologic control on DOM quality that differed between the two transects, evidencing differences in their connectivity to areas of high agricultural activity. Across all sites, a significant long-term increasing trend in the fluorescence index was observed, associated with a positive correlation with precipitation and also potential changes in agricultural inputs, with other features associated with drought and hurricanes. Lastly, the cumulative sums analysis revealed differences between the two transects in the sensitivity of DOM composition to decreased water levels associated with 30-year climate scenarios, with the less hydrologically dynamic transect exhibiting greater potential sensitivity.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Third integer resonant extraction transit time simulation studies

In this work, we present the investigation of transit time of particles in the non-linear third-integer resonant extraction process. Transit time is defined as the number of turns a particle takes to get extracted once it is in the unstable region in the phase space, i.e., outside the triangular separatrix in case of third-integer resonance. The study of transit time is important because transit time directly contributes to the beam response time during resonant extraction and thus knowing it apriori would be practically useful in designing of the extraction system. In this work, we shall investigate the analytical derivation of the transit time of particles (to the first order Kobayashi Hamiltonian) in different parts of the phase space distribution and compare against the analytical results. We also compare the simulation result of the transit time of particles (with higher statistics) for the static as well as dynamic extraction conditions cases, particularly in the context of resonant extraction parameters for Mu2e experiment at Fermilab.

Narayanan, Aakaash [Fermilab]↗

DOE OSTI2025

Third Integer Resonant Extraction Transit Time Simulation Studies

Narayanan, Aakaash [Fermilab] (ORCID:0000000157944↗

DOE OSTI2026

Impacts of Year-to-Year Weather Variability and Inter-Panel Spacing on Crop Yields in a Massachusetts Agrivoltaics System

This presentation summarizes recently published work on agrivoltaic field data and irradiance modeling. The body of agrivoltaic field data is still growing, and crop responses to different solar configurations under different local climates are highly varied. We investigate the impact of adding spacing between adjacent solar panels in a fixed-tilt system to improve light diffusion to crops. For four crops (broccoli, peppers, kale, Swiss chard) grown across 3 years in an agrivoltaic system in Massachusetts, we found that only kale had a linearly increasing trend as the inter-panel spacing increased from 0.6 m to 1.5 m (2 ft to 5 ft). However, there were significant year-to-year differences in the yield of agrivoltaic versus control fields. Agrivoltaic and full sun fields produced equivalent yields in a hot, dry year, whereas the full-sun control beds produced more salable yield for all four crops in a warm, wet year. This demonstrates variability of agricultural outcomes and the need for more multi-year studies to ensure agrivoltaic impacts are not under- or overestimated.

14 SOLAR ENERGY↗

DOE OSTI2022

Intrinsic Junction Thermocouples For Surface Temperature Measurement

The material selection, attachment, and calibration of intrinsic junction thermocouples is discussed. Intrinsic junction thermocouples are utilized for their fast response times in surface measurements during high thermal transients. The as-is measurement and calibrated corrections are presented for both the out-of-pile pulse power system and the critical heat flux experiments in the Idaho National Laboratory Transient Reactor Test facility. The fin effect from temperature measurements on the surface of fuel cladding is corrected using a linear coefficient, m, and temperature corrections are sometimes 100 °C or more—representing 18% of the original temperature measurement.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Comparative Analysis of HEATNETS for Geothermal Network Performance: Preprint

Thermal energy networks (TENs), also known as 5th generation district energy systems, or more specifically geothermal networks when exchanging heat with geothermal boreholes, are an important technology for decarbonization. In these networks an ambient loop connects buildings and thermal sources, such as a borehole field, to exchange energy and maintain a desired loop temperature. Water-source heat pumps are used at the buildings to connect to the ambient or thermal loop to meet to the building heating and cooling loads and maintain comfort. A semi-transient, reduced-order technical model and techno-economic model, called HEATNETS, has been developed at NREL that captures the flow of energy around a TEN. In this work, a comparison of the HEATNETS technical model and a well-known coding platform used for modeling geothermal networks, TRNSYS, has been completed for a proposed geothermal network as a verification process. Hourly data provided from the TRNSYS simulation included building loads, pumping power, heat pump power, temperature entering and leaving the borehole field, and mass flow rates. The hourly borehole temperatures were used to create a linear regression model utilized in HEATNETS to estimate the borehole field heat exchange. The building loads and mass flow rates were direct inputs to HEATNETS while the pumping power, heat pump power, borehole temperatures, and coefficients of performance were all simulated and calculated by HEATNETS, allowing for direct comparison of the thermal energy transfer, rather than also comparing control systems responses. HEATNETS considers the full process from design inputs to economic outputs and can provide modeling options for high-level initial system design and operational optimization. This study shows that HEATNETS, while not intended to replace other modeling tools, can be a unique modeling tool for the performance of a full geothermal network system.

15 GEOTHERMAL ENERGY↗

DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗

DOE OSTI2020

Dispersive corrections in elastic electron-nucleus scattering: an investigation in the intermediate energy regime and their impact on the nuclear matter

Measurements of elastic electron scattering data within the past decade have highlighted two-photon exchange contributions as a necessary ingredient in theoretical calculations to precisely evaluate hydrogen elastic scattering cross sections. This correction can modify the cross section at the few percent level. In contrast, dispersive effects can cause significantly larger changes from the Born approximation. The purpose of this experiment is to extract the carbon-12 elastic cross section around the first diffraction minimum, where the Born term contributions to the cross section are small to maximize the sensitivity to dispersive effects. The analysis uses the LEDEX data from the high resolution Jefferson Lab Hall A spectrometers to extract the cross sections near the first diffraction minimum of 12 C at beam energies of 362 MeV and 685 MeV. The results are in very good agreement with previous world data, although with less precision. The average deviation from a static nuclear charge distribution expected from linear and quadratic fits indicate a 30.6% contribution of dispersive effects to the cross section at 1 GeV. The magnitude of the dispersive effects near the first diffraction minimum of 12 C has been confirmed to be large with a strong energy dependence and could account for a large fraction of the magnitude for the observed quenching of the longitudinal nuclear response. These effects could also be important for nuclei radii extracted from parity-violating asymmetries measured near a diffraction minimum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Center for Integrated Simulation of Fusion Relevant RF Actuators

This project was part of the “Center for Integrated Simulation of Fusion Relevant RF Actuators” SciDAC-4 project, led by Dr. Paul Bonoli (MIT). Rather than use an acronym (CISFRRFA), the project will be referred to in this document as the “RF-SciDAC4”. The larger SciDAC-4 project goals were to: 1. Develop an integrated simulation of the antenna + sheath + scrape-off-layer + core plasma system which fully utilizes leadership class computing resources to move towards a quantitative predictive capability for the response to RF power. 2. Work closely with the SciDAC-4 Whole Device Modeling (WDM) community to make both our new code development efforts, as well as the established hierarchy of RF tools, available within their environment, and to utilize WDM technologies to implement the couplings below. 3. Validate this predictive capability on appropriately diagnosed experiments including dedicated RF test stands, linear devices, and existing tokamaks. 4. Use these tools to inform design of robust, impurity-mitigating RF heating and current drive sources for future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Characterization of Pinhole Collimators for High-Resolution Gamma Imaging of Irradiated Fuel

Post-irradiation examination (PIE) of nuclear fuels requires imaging tools capable of resolving isotopic and spatial features with high throughput. This project contributes to a proof-of-concept effort aimed at advancing gamma emission tomography (GET) by evaluating novel fine-aperture pinhole collimators. Two Rose’s metal collimators, 100 µm (20° acceptance angle) and 350 µm (30° acceptance angle), were prototyped and characterized for their effectiveness in transporting gamma rays through the pinhole aperture. To support data collection, a Python-based data acquisition system was developed to coordinate a rotation stage, linear stage, and CZT detector, reducing latency in high-rate gamma event logging to one second per acquisition. Queue-based file writing enabled seamless real-time data capture for count rates up to 35,000 counts per second (cps). List-mode parsing algorithms were implemented to differentiate single and simultaneous gamma interactions for future tomographic reconstruction. Detector response was evaluated in both spectroscopy and list mode acquisition methods across varying source distances to confirm absolute and collimator efficiencies. Preliminary efficiency figures suggest effective collimation of gamma-rays with energies below 700 keV, with ~4% residual intensity through the aperture for Cs-137. The impact of collimator geometry on image quality is currently being evaluated. This groundwork supports the ongoing development of a sub mm resolution cone-beam CT system for imaging fuel phantoms, an essential step toward improving GET efficiency and accelerating nuclear fuel qualification efforts.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

DOE OSTI2024

Local carbon reserves are insufficient for phloem terpene induction during drought in Pinus edulis in response to bark beetle–associated fungi

Stomatal closure during drought inhibits carbon uptake and may reduce a tree's defensive capacity. Limited carbon availability during drought may increase a tree's mortality risk, particularly if drought constrains trees' capacity to rapidly produce defenses during biotic attack. We parameterized a new model of conifer defense using physiological data on carbon reserves and chemical defenses before and after a simulated bark beetle attack in mature Pinus edulis under experimental drought. Attack was simulated using inoculations with a consistent bluestain fungus (Ophiostoma sp.) of Ips confusus, the main bark beetle colonizing this tree, to induce a defensive response. Trees with more carbon reserves produced more defenses but measured phloem carbon reserves only accounted for c. 23% of the induced defensive response. Our model predicted universal mortality if local reserves alone supported defense production, suggesting substantial remobilization and transport of stored resin or carbon reserves to the inoculation site. Our results show that de novo terpene synthesis represents only a fraction of the total measured phloem terpenes in P. edulis following fungal inoculation. Without direct attribution of phloem terpene concentrations to available carbon, many studies may be overestimating the scale and importance of de novo terpene synthesis in a tree's induced defense response.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Shock compression of additively manufactured high solids loaded polymer composites with aligned porosity

The effects of aligned porosity on the shock-compression response of additively manufactured (AM) composite structures with collinear and log-cabin geometries are investigated. The composites contain two inorganic and two organic particle populations and a UV-curable polymer binder. X-ray phase contrast imaging (X-PCI) is used as an in-material diagnostic in plate-impact experiments performed at the Dynamic Compression Sector of the Advanced Photon Source at the Argonne National Laboratory. Macroscale shock and particle velocities and microscale pore collapse velocities are measured from X-PCI images obtained from samples impacted in various orientations relative to the AM build directions. The linear shock and particle velocity relation for the structures with collinear build geometry shows directional differences, with the build direction exhibiting the highest slope and the filament print direction having the shallowest slope. Composites with log-cabin geometry, due to their slightly lower porosity and smaller sized pores, show higher shock velocities than the collinear structures at the same particle velocity. Microscale directionality effects are also observed, with the pore collapse velocity dependent on the pore configuration. Specifically, the pore collapse velocity is highest for pores elongated along the impact direction which is the case for the colinear structures impacted in the filament print direction, which is consistent with the directionality trend observed in the data for the shock and particle velocity equation of state. The pore collapse velocity approaches ∼80 % of the shock velocity and is several times the average bulk particle velocity, suggesting a highly hydrodynamic mechanism of pore collapse. Overall, both the macroscale response as measured by the shock and particle velocity relationship, and the microscale response as measured by the pore collapse velocity, are found to be anisotropic and dependent on the AM composite structure.

Wagner, K. B. [Georgia Institute of Technology, At↗

DOE OSTI2021

Considering computational speed vs. accuracy: Choosing appropriate mesoscale RVE boundary conditions

Modeling a material’s microstructure using continuum theories allows for inspection of the relationship between coarse scale and fine scale behaviors. Computational limits generally require selection of a sub-volume from a bulk sample in order to directly model the microstructure. Boundary conditions are applied to the sub-volume to mimic the excluded bulk material. Appropriate selection of boundary conditions helps effectively determine the appropriate spatial scale required of the sub-volume. Applicable boundary conditions include direct displacement, periodic, and uniform traction. While direct displacement and periodic boundary conditions are commonly used, uniform traction boundary conditions have seen limited use due to rigid body stability issues in simulations of compression or shear deformation. A new application of uniform traction boundary conditions was developed through linear constraint equations, similar to approaches employed by direct displacement and periodic boundary conditions, to quench rigid body motions with minimal interference of the relative deformation of the model. These boundary conditions were tested by compressing several synthetically generated periodic microstructures using the finite element method. Evaluating the effective stiffness along the compression axis, the direct displacement boundary condition produced the stiffest response, whereas the uniform traction boundary condition produced the most compliant. Periodic boundary conditions produced the same response for all volumes analyzed and both the direct displacement and uniform traction boundary conditions trended toward the periodic response as the domain volume increased. Computational performance was also evaluated for each boundary condition using implicit and explicit solvers. Direct displacement boundary conditions presented the lowest computational cost of all of the boundary conditions followed by periodic then uniform traction. The computational expense of periodic and uniform traction boundary conditions limited the viable spatial scale and mesh resolutions able to be simulated. Selection of appropriate boundary conditions for specific uses need to be a balance between allowable computational expense and accuracy of the method. Techniques for evaluating which boundary conditions to use are discussed.

42 ENGINEERING↗

DOE OSTI2025

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following Du Toit et al. (2018), the scheme transforms mixed-derivative diffusion fluxes (which are responsible for the lack of a discrete maximum principle) into nonlinear advective fluxes, amenable to nonlinear-solver-friendly monotonicity-preserving limiters. The scheme enables accurate multi-dimensional heat transport simulations with up to seven orders of magnitude of heat-transport-coefficient anisotropies with low cross-field numerical error pollution and excellent algorithmic performance, with the number of linear iterations scaling very weakly with grid resolution and grid anisotropy, and scaling with the square-root of the implicit timestep. We propose a multigrid preconditioning strategy based on a lower-order approximation that renders the scheme efficient and scalable under grid refinement. Several numerical tests are presented that display the expected spatial convergence rates and strong algorithmic performance, including fully nonlinear magnetohydrodynamics simulations of kink instabilities in a Bennett pinch in 2D helical geometry and of ITER in 3D toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Bathymetric Influences on Antarctic Ice-Shelf Melt Rates

Ocean bathymetry exerts a strong control on ice sheet-ocean interactions within Antarctic ice-shelf cavities, where it can limit the access of warm, dense water at depth to the underside of floating ice shelves. However, ocean bathymetry is challenging to measure within or close to ice-shelf cavities. It remains unclear how uncertainty in existing bathymetry datasets affect simulated sub-ice-shelf melt rates. In this work, we infer linear sensitivities of ice-shelf melt rates to bathymetric shape with grid-scale detail by means of the adjoint of an ocean general circulation model. Both idealized and realistic-geometry experiments of sub-ice-shelf cavities in West Antarctica reveal that bathymetry has a strong impact on melt in localized regions such as topographic obstacles to flow. Moreover, response of melt to bathymetric perturbation is found to be non-monotonic, with deepening leading to either increased or decreased melt depending on location. Our computational approach provides a comprehensive way of identifying regions, where refined knowledge of bathymetry is most impactful, and also where bathymetric errors have relatively little effect on modeled ice sheet-ocean interactions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

High-throughput platform for yeast morphological profiling predicts the targets of bioactive compounds

Abstract Morphological profiling is an omics-based approach for predicting intracellular targets of chemical compounds in which the dose-dependent morphological changes induced by the compound are systematically compared to the morphological changes in gene-deleted cells. In this study, we developed a reliable high-throughput (HT) platform for yeast morphological profiling using drug-hypersensitive strains to minimize compound use, HT microscopy to speed up data generation and analysis, and a generalized linear model to predict targets with high reliability. We first conducted a proof-of-concept study using six compounds with known targets: bortezomib, hydroxyurea, methyl methanesulfonate, benomyl, tunicamycin, and echinocandin B. Then we applied our platform to predict the mechanism of action of a novel diferulate-derived compound, poacidiene. Morphological profiling of poacidiene implied that it affects the DNA damage response, which genetic analysis confirmed. Furthermore, we found that poacidiene inhibits the growth of phytopathogenic fungi, implying applications as an effective antifungal agent. Thus, our platform is a new whole-cell target prediction tool for drug discovery.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Local electric field measurement in GaN diodes by exciton Franz–Keldysh photocurrent spectroscopy

The eXciton Franz–Keldysh (XFK) effect is observed in GaN p–n junction diodes via the spectral variation of photocurrent responsivity data that redshift and broaden with increasing reverse bias. Photocurrent spectra are quantitatively fit over a broad photon energy range to an XFK model using only a single fit parameter that determines the line shape and the local bias (V l ), uniquely determining the local electric field maximum and depletion widths. As expected, the spectrally determined values of V l vary linearly with the applied bias (V) and reveal a large reduction in the local electric field due to electrostatic non-uniformity. The built-in bias (V bi ) is estimated by extrapolating V l at V=0, which, when compared with independent C-V measurements, indicates an overall ±0.31 V accuracy of V l . This demonstrates sub-bandgap photocurrent spectroscopy as a local probe of electric field in wide bandgap diodes that can be used to map out regions of device breakdown (hot spots) for improving electrostatic design of high-voltage devices.

42 ENGINEERING↗