Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Yield function”

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,063 records · Page 59

Combustion of solid carbon rods in zero and normal gravity

In order to investigate the mechanism of carbon combustion and to assess the importance of gravitational induced convection on the process, zero and normal gravity experiments were conducted in which spectroscopic carbon rods were resistance ignitied and burned in dry oxygen environments. In the zero-gravity drop tower tests, a blue flame surrounded the rod, showing that a gas phase reaction in which carbon monoxide was oxidized to carbon dioxide was taking place. The ratio of flame diameter to rod diameter was obtained as a function of time. It was found that this ratio was inversely proportional to both the oxygen pressure and the rod diameter. In the normal gravity tests, direct mass spectrometric sampling was used to measure gas phase concentrations. The gas sampling probe was positioned near the circumference of a horizontally mounted 0.615 cm diameter carbon rod, either at the top or at angles of 45 deg to 90 deg from the top, and yielded concentration profiles of CO2, CO, and O2 as a function of distance from the surface. The mechanism controlling the combustion process was found to change from chemical process control at the 90 deg and 45 deg probe positions to mass transfer control at the 0 deg probe position at the top of the rod. Under the experimental conditions used, carbon combustion was characterized by two surface reactions, 2C + O2 yields 2CO and CO2 + C yields 2CO, and a gas phase reaction, 2CO + O2 yields 2CO2.

Spuckler, C. M.↗

Thermal shock in a circumferentially cracked hollow cylinder with cladding

An theoretical analysis is presented which demonstrates the effect of cladding on the thermal resistance of a circumferentially cracked hollow cylinder. The cladding is assumed to be bonded to the inner wall of the hollow cylinder. The axisymmetric circumferential crack may be either embedded in the cylinder wall or may be an edge crack which passes through the clad and opens into the inner wall of the hollow cylinder. The problem is formulated mathematically and a solution is found which is in the form of a single integral equation. The integral equation is solved numerically and yields estimates of transient temperature distributions, thermal stresses in the uncracked cylinder, and stress intensity factors as a function of time for various cladding thickness to cylinder wall thickness ratios. It is shown that yielding of the clad under certain conditions can result in a reduction in the magnitude of the stress intensity factor for the crack tip in the elastic base material.

Nied, H. F.↗

Exchange-correlation energies of atoms from efficient density functionals: influence of the electron density

The primary goal of Kohn–Sham density functional theory is to evaluate the exchange-correlation contribution to electronic properties. However, the accuracy of a density functional can be affected by the electron density. Here we apply the nonempirical Tao–Mo (TM) semilocal functional to study the influence of the electron density on the exchange and correlation energies of atoms and ions, and compare the results with the commonly used nonempirical semilocal functionals local spin-density approximation (LSDA), Perdew–Burke–Ernzerhof (PBE), Tao–Perdew–Staroverov–Scuseria (TPSS), and hybrid functional PBE0. We find that the spin-restricted Hartree–Fock density yields the exchange and correlation energies in good agreement with the Optimized Effective Potential method, particularly for spherical atoms and ions. However, the errors of these semilocal and hybrid functionals become larger for self-consistent densities. We further find that the quality of the electron density have greater effect on the exchange-correlation energies of kinetic energy density-dependent meta-GGA functionals TPSS and TM than on those of the LSDA and GGA, and therefore, should have greater influence on the performance of meta-GGA functionals. Lastly, we show that the influence of the density quality on PBE0 is slightly reduced, compared to that of PBE, due to the exact mixing.

74 ATOMIC AND MOLECULAR PHYSICS↗

Simple Test Functions in Meshless Local Petrov-Galerkin Methods

Two meshless local Petrov-Galerkin (MLPG) methods based on two different trial functions but that use a simple linear test function were developed for beam and column problems. These methods used generalized moving least squares (GMLS) and radial basis (RB) interpolation functions as trial functions. These two methods were tested on various patch test problems. Both methods passed the patch tests successfully. Then the methods were applied to various beam vibration problems and problems involving Euler and Beck's columns. Both methods yielded accurate solutions for all problems studied. The simple linear test function offers considerable savings in computing efforts as the domain integrals involved in the weak form are avoided. The two methods based on this simple linear test function method produced accurate results for frequencies and buckling loads. Of the two methods studied, the method with radial basis trial functions is very attractive as the method is simple, accurate, and robust.

Raju, Ivatury S.↗

Exciting H 2 Molecules for Graphene Functionalization

Hydrogen functionalization of graphene by exposure to vibrationally excited H 2 molecules is investigated by combined scanning tunneling microscopy, high-resolution electron energy loss spectroscopy, X-ray photoelectron spectroscopy measurements, and density functional theory calculations. The measurements reveal that vibrationally excited H 2 molecules dissociatively adsorb on graphene on Ir(111) resulting in nanopatterned hydrogen functionalization structures. Calculations demonstrate that the presence of the Ir surface below the graphene lowers the H 2 dissociative adsorption barrier and allows for the adsorption reaction at energies well below the dissociation threshold of the H–H bond. The first reacting H 2 molecule must contain considerable vibrational energy to overcome the dissociative adsorption barrier. However, this initial adsorption further activates the surface resulting in reduced barriers for dissociative adsorption of subsequent H 2 molecules. This enables functionalization by H 2 molecules with lower vibrational energy, yielding an avalanche effect for the hydrogenation reaction. These results provide an example of a catalytically active graphene-coated surface and additionally set the stage for a re-interpretation of previous experimental work involving elevated H 2 background gas pressures in the presence of hot filaments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assembly and analysis of fragmentation data for liquid propellant vessels

Fragmentation data was assembled and analyzed for exploding liquid propellant vessels. These data were to be retrieved from reports of tests and accidents, including measurements or estimates of blast yield, etc. A significant amount of data was retrieved from a series of tests conducted for measurement of blast and fireball effects of liquid propellant explosions (Project PYRO), a few well-documented accident reports, and a series of tests to determine auto-ignition properties of mixing liquid propellants. The data were reduced and fitted to various statistical functions. Comparisons were made with methods of prediction for blast yield, initial fragment velocities, and fragment range. Reasonably good correlation was achieved. Methods presented in the report allow prediction of fragment patterns, given type and quantity of propellant, type of accident, and time of propellant mixing.

Baker, W. E.↗

Statistical measures of Planck scale signal correlations in interferometers

Here, a model-independent statistical framework is presented to interpret data from systems where the mean time derivative of positional cross correlation between world lines, a measure of spreading in a quantum geometrical wave function, is measured with a precision smaller than the Planck time. The framework provides a general way to constrain possible departures from perfect independence of classical world lines, associated with Planck scale bounds on positional information. A parametrized candidate set of possible correlation functions is shown to be consistent with the known causal structure of the classical geometry measured by an apparatus, and the holographic scaling of information suggested by gravity. Frequency-domain power spectra are derived that can be compared with interferometer data. As a result, simple projections of sensitivity for specific experimental set-ups suggests that measurements will directly yield constraints on a universal time derivative of the correlation function, and thereby confirm or rule out a class of Planck scale departures from classical geometry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Application of Hybrid Genetic Algorithm Routine in Optimizing Food and Bioengineering Processes

Optimization is a crucial step in the analysis of experimental results. Deterministic methods only converge on local optimums and require exponentially more time as dimensionality increases. Stochastic algorithms are capable of efficiently searching the domain space, however convergence is not guaranteed. This article demonstrates the novelty of hybrid genetic algorithm which combines both stochastic and deterministic routine for improved optimization results. A new hybrid genetic algorithm developed is applied to the Ackley benchmark function as well as case studies in foods, biofuels, and biotechnology processes. For each case study the hybrid genetic algorithm found a better optimum candidate than reported by the sources. In the case of food processing the hybrid genetic predicted 5% higher anthocyanin yield. Optimization of bio-oil production using HGA resulted in a 4.73 % higher yield. In enzyme production process HGA predicted 2.18 IU/gds of higher xylanase yield. In conclusion, hybridization of genetic algorithm with a deterministic algorithm resulted in maximum optimum compared to the statistical methods.

09 BIOMASS FUELS↗

Explosive Yield Estimation Using Regional Seismic Moment Tensors

Here, we use the Pasyanos and Chiang (2022) data set to calculate the seismic moment M 0 for each explosion and use the measured explosive yield W to validate the W~M 0 relationship in Denny and Johnson (1991; hereafter, DJ91). The M 0 is corrected by transforming to a potency tensor and applying more appropriate near-source geophysical parameter values in the moment estimate. The mean residual between observed and predicted yield is near zero; however, the standard deviation of the residuals results in an F-value (a 95% confidence factor) of about 5. We re-estimate the coefficients in the DJ91 model and find similar values and only a slight improvement in the F-value. Next, we embark on a similar model selection process as DJ91, allowing for non-cube-root yield scaling and other plausible near-source elastic moduli. As was found by DJ91, the yield dependence is not significantly different from unity, and a cube root assumption is valid. Therefore, we yield scale the seismic moment and test the significance of all plausible explanatory variables. Isotropic moment performs better in the response variable than total moment. The preference for isotropic moment could be due to its relationship to volume change, which would be more directly affected by explosive yield. Surprisingly, we find that the overburden pressure, which is a function of depth, is not a significant parameter in the model. We hypothesize that this is due to the competing depth effects on source asymmetry and the incorporation of depth in the Green’s functions used to calculate the seismic moment tensors. Importantly, this emphasizes that only seismic moment tensor-derived moments should be used in these models. After removing insignificant model parameters, we are left with a simple model to predict explosive yield $\widehat{W}$ in kt from isotropic moment M I in N·m, $\widehat{W}$=κ –1.4132 10 0.035626GP M I , in which κ and GP are the near-source bulk modulus and gas porosity in Pa and %, respectively. The F-value for this model is approximately 3.

58 GEOSCIENCES↗

Characteristics of Solar Energetic Ions as a Function of Longitude

Since the 2006 launch of STEREO , multi-spacecraft studies have yielded several surprising results regarding the spread of solar energetic particles (SEPs) within the inner heliosphere. We have investigated the role of energy and ridigity, using ACE and STEREO 10 MeV n{sup −1} oxygen data to identify 41 large SEP events observed by two or three spacecraft. We calculated fluence spectra from ∼0.1 to >10 MeV n{sup −1} for H, He, O, and Fe for each event at the observing spacecraft (including SOHO and GOES ). The particle fluences at 0.3, 1, and 10 MeV n{sup −1} were examined as a function of the distance between the associated solar flare longitude and the spacecraft magnetic footpoints at the Sun to determine the longitudinal spread of particles and study how the distribution centers and widths depend on energy and charge-to-mass (Q/M) for the first time. On average, the three-spacecraft event distributions were centered at 22 ± 4° west of the flare site and were 43 ± 1° wide, though there was substantial variability, while the fit to the aggregate of the two-spacecraft event fluences yielded significantly wider distributions at 0.3 and 1 MeV n{sup −1}. The widths derived from both the three- and two-spacecraft events show an energy dependence with distributions narrowing with increasing energy, consistent with lower energy ions experiencing more field line co-rotation, or being accelerated over a larger portion of the CME-driven shock or for longer times as the shock expands. Surprisingly, no clear evidence was found for a Q/M dependence to the widths or centers suggesting that rigidity-related processes are not the dominant means of spreading particles in longitude.

79 ASTRONOMY AND ASTROPHYSICS↗

CYP76BK1 orthologs catalyze furan and lactone ring formation in clerodane diterpenoids across the mint family

The Lamiaceae (mint family) is the largest known source of furanoclerodanes, a subset of clerodane diterpenoids with broad bioactivities including insect antifeedant properties. The Ajugoideae subfamily, in particular, accumulates significant numbers of structurally related furanoclerodanes. The biosynthetic capacity for formation of these diterpenoids is retained across most Lamiaceae subfamilies, including the early-diverging Callicarpoideae which forms a sister clade to the rest of Lamiaceae. VacCYP76BK1, a cytochrome P450 monooxygenase from Vitex agnus-castus, was previously found to catalyze the formation of the proposed precursor to furan and lactone-containing labdane diterpenoids. Through transcriptome-guided pathway exploration, we identified orthologs of VacCYP76BK1 in Ajuga reptans and Callicarpa americana. Functional characterization demonstrated that both could catalyze the oxidative cyclization of clerodane backbones to yield a furan ring. Subsequent investigation revealed a total of 10 CYP76BK1 orthologs across six Lamiaceae subfamilies. Through analysis of available chromosome-scale genomes, we identified four CYP76BK1 members as syntelogs within a conserved syntenic block across divergent subfamilies. This suggests an evolutionary lineage that predates the speciation of the Lamiaceae. Functional characterization of the CYP76BK1 orthologs affirmed conservation of function, as all catalyzed furan ring formation. Additionally, some orthologs yielded two novel lactone ring moieties. The presence of the CYP76BK1 orthologs across Lamiaceae subfamilies closely overlaps with the distribution of reported furanoclerodanes. Together, the activities and distribution of the CYP76BK1 orthologs identified here support their central role in furanoclerodane biosynthesis within the Lamiaceae family. Our findings lay the groundwork for biotechnological applications to harness the economic potential of this promising class of compounds.

59 BASIC BIOLOGICAL SCIENCES↗

Discontinuous Galerkin and Related Methods for ODE

A defining feature of the discontinuous Galerkin (DG) method for ODE is that the piecewise polynomial solution can have a jump discontinuity at the beginning of each step. Starting from the standard integral formulation, the DG method is derived here in differential form. The key ingredient is a polynomial called the correction function, which helps ‘correct’ the discontinuous solution by approximating the jump and yields a continuous one. Under the right Radau quadrature, this continuous solution is identical to the solutions by the right Radau collocation and the continuous Galerkin (CG) methods. Next, the correction function facilitates the construction of the associated implicit Runge-Kutta schemes (IRK-DG). Different quadratures for DG result in different IRK-DG methods: left Radau quadrature in Radau IA, right Radau quadrature in Radau IIA or right Radau collocation, and Gauss quadrature in a method called DG-Gauss. The construction of IRK-DG clarifies the meaning and facilitates the proofs of various 𝐵(𝑝), 𝐶(𝜂), and 𝐷(𝜁) conditions for accuracy. The two consequences of these conditions are that all 𝑠-stage IRK-DG methods are accurate to order 2𝑠 − 1, and the IRK-DG methods of Radau type are unique. Numerical examples showing the behavior of the DG solutions are provided. In all, the correction function plays a key role and helps establish the relations among the DG, IRK-DG, collocation, and CG schemes.

numerical methods↗

Discontinuous Galerkin and Related Methods for ODE

Starting from the standard integral formulation, the DG method is derived here in differential form. The key ingredient is a polynomial called the correction function, which helps ‘correct’ the discontinuous solution by approximating the jump and yields a continuous one. Under the right Radau quadrature, this continuous solution is identical to the solutions by the right Radau collocation and the continuous Galerkin (CG) methods. Next, the correction function facilitates the construction of the associated implicit Runge-Kutta schemes (IRK-DG). Different quadratures for DG result in different IRK-DG methods: left Radau quadrature in Radau IA, right Radau quadrature in Radau IIA or right Radau collocation, and Gauss quadrature in a method called DG-Gauss. The construction of IRK-DG clarifies the meaning and facilitates the proofs of various 𝐵(𝑝), 𝐶(𝜂), and 𝐷(𝜁) conditions for accuracy. The two consequences of these conditions are that all 𝑠-stage IRK-DG methods are accurate to order 2𝑠− 1, and the IRK-DG methods of Radau type are unique. Numerical examples showing the behavior of the DG solutions are provided. In all, the correction function plays a key role and helps establish the relations among the DG, IRK DG, collocation, and CG methods.

Numerical Methods for Ordinary Differential Equati↗

An Ab Initio Study of CuCO

Modified coupled-pair functional (MCPF) calculations and coupled cluster singles and doubles calculations, which include a perturbational estimate of the connected triples [CCSD(T)], yield a bent structure for CuCO, thus, supporting the prediction of a nonlinear structure based on density functional (DF) calculations. Our best estimate for the binding energy is 4.9 +/- 1.4 kcal/mol; this is in better agreement with experiment (6.0 +/- 1.2 kcal/mol) than the DF approach which yields a value (19.6 kcal/mol) significantly larger than experiment.

Bauschlicher, Charles W., Jr.↗

Communication: Improved ab initio molecular dynamics by minimally biasing with experimental data

Accounting for electrons and nuclei simultaneously is a powerful capability of ab initio molecular dynamics (AIMD). However, AIMD is often unable to accurately reproduce properties of systems such as water due to inaccuracies in the underlying electronic density functionals. This shortcoming is often addressed by added empirical corrections and/or increasing the simulation temperature. We present here a maximum-entropy approach to directly incorporate limited experimental data via a minimal bias. Biased AIMD simulations of water and an excess proton in water are shown to give significantly improved properties both for observables which were biased to match experimental data and for unbiased observables. Here, this approach also yields new physical insight into inaccuracies in the underlying density functional theory as utilized in the unbiased AIMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optical and radio variability in blazars

An analysis is presented of the optical (0.44 micron) and radio flux density variations (4.8, 8.0, and 14.5 GHz) of the five best observed variable AGN, of which three are BL Lac objects and two are optically violent variable (OVV) quasars. For each object the nature of the variability at individual wavelengths is determined by using a structure function analysis and the autocorrelation function as calculated from the discrete correlation technique. Both techniques yield similar results, showing that the optical variations are similar for all objects and are a combination of flicker and shot noise. The radio variability is distinct, being described by shot noise (less variation at short times), with the two OVVs exhibiting the most extreme behavior. The time lag between the variations at the three radio frequencies and the optical flux is studied with the discrete cross-correlation function and the mean variance methods. It is concluded that the optical and radio emitting regions are physically related, exist on distinct time scales, and excite their synchrotron emitting plasmas differently.

Hufnagel, Beth R.↗

Maize yield under a changing climate: The hidden role of vapor pressure deficit

Temperatures over the next century are expected to rise to levels detrimental to crop growth and yield. As the atmosphere warms without additional water vapor input, vapor pressure deficit (VPD) increases as well. Increased temperatures and accompanied elevated VPD levels can both lead to negative impacts on crop yield. The independent importance of VPD, however, is often neglected or conflated with that from temperature due to a tight correlation between the two climate factors. We used a coupled process-based crop (MAIZSIM) and soil (2DSOIL) model to gain a mechanistic understanding of the independent roles temperature and VPD play in crop yield projections, as well as their interactions with rising CO 2 levels and changing precipitation patterns. We found that by separating out the VPD effect from rising temperatures, VPD increases had a greater negative impact on yield (12.9 ± 1.8%, increase in VPD associated with 2 °C warming) compared to that from warming (8.5 ± 1.4%, the direct effect of 2 °C warming). The negative impact of these two factors varied with precipitation levels and influenced yield through separate mechanisms. Warmer temperatures caused yield loss mainly through shortening the growing season, while elevated VPD increased water loss and triggered several water stress responses such as reduced photosynthetic rates, lowered leaf area development, and shortened growing season length. Elevated CO 2 concentrations partially alleviated yield loss under warming or increased VPD conditions through water savings, but the impact level varied with precipitation levels and was most pronounced under drier conditions. These results demonstrate the key role VPD plays in crop growth and yield, displaying a magnitude of impact comparative to temperature and CO 2 . A mechanistic understanding of the function of VPD and its relation with other climate factors and management practices is critical to improving crop yield projections under a changing climate.

54 ENVIRONMENTAL SCIENCES↗

High-shade dryland agrivoltaic conditions enhanced carbon uptake and water-use efficiency in zucchini ( Cucurbita pepo )

Introduction: The increasing global demand for food and energy is intensifying land-use competition. Agrivoltaic systems are a multifunctional land-use approach that vertically integrates the production of agricultural crops and solar power on the same land area. Most food crops are adapted to full-sun conditions, and the physiological responses of these crops to the novel microclimate under solar panels remain poorly understood. We hypothesized that the microclimate beneath the high-density photovoltaic system would influence carbon uptake, water use, and yield outcomes of zucchini summer squash.Methods: We conducted a field experiment in a hot, semi-arid climate on zucchini (Cucurbita pepo). Plants were grown under an agrivoltaic system with a 75% ground cover ratio (GCR) and in a full-sun control plot, each with two irrigation regimes (100 and 50%). We measured leaf-level photosynthesis, microclimate variables, and fruit yield at plant maturity and throughout the growing season.Results: The agrivoltaic array reduced photosynthetically active radiation (PAR) by ~79%, resulting in a cooler (−1.1 °C), more humid environment with higher soil moisture. These microclimatic conditions enhanced midday photosynthesis and daily cumulative carbon uptake. However, fruit yield was consistently lower under the panels, indicating a shift in carbon allocation toward vegetative growth. Photosynthesis was primarily driven by PAR across treatments, while soil moisture significantly influenced photosynthesis only in the control plots, suggesting water limitation was alleviated under the panels.Discussion: These findings highlight a trade-off between improved physiological performance and reduced yield under high-density agrivoltaics. While the system buffered heat and drought stress and improved overall plant function, excessive shade reduced reproductive output. Optimizing panel density or selecting crops cultivated for non-fruit yields will be essential for balancing food production and energy generation in dryland agrivoltaic settings.

14 SOLAR ENERGY↗