Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear Parameter Varying”

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.

484 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects of internal swelling on residual elasticity of a quasi-brittle material through a composite sphere model

This work describes the development of a micromechanical-based constitutive model accounting for the effect of internal expansion on the residual elasticity of a Hashin composite material. This material is made of spherical inclusions that are subjected to gradual swelling within a quasi-brittle matrix. The main focus of this work is to describe and analyze the material mechanical response, with an additional focus on the internal swelling’s effect on the stress–strain response and residual elasticity. Microstructural features and parameters of major importance for the mechanical responses were identified. The innovative characteristics of the proposed approach are summarized as follows: (1) a full determination of the physics of a complete-damage problem throughout the whole process of inclusion swelling with upscaling techniques, which transfers the microcrack-related properties from the lower scale to upper scale; and (2) an evolution of the mechanical fields and the corresponding residual elasticity for various inclusion swelling levels. Concerning the matrix–inclusion composite, it was hypothesized that only the matrix was susceptible to cracking, with varied degrees of damage, whereas the inclusions behave elastically and the elastic modulus of the expanding inclusions remains constant. It is to emphasize that the gradual swelling of inclusions is modeled by an increasing strain in the current micromechanical-based constitutive model, and the microcracks are represented by a set of randomly oriented penny-shaped microcracks with identical radii (namely, closed cracks). The main contribution of the current research is to establish the exact mathematical solutions for the mechanical fields (stress, strain) caused by the swelling of the inclusions (mechanical loading) and derive the effective residual elasticity of a composite (structural response) subject to internal expansion and quasi-brittle damage. Based on the assumption of closed cracks that could be extended to open cracks in the upcoming work, the results proposed by the present paper help to understand the non-linear mechanical behavior of quasi-brittle materials subject to microcracking and provide a theoretical framework to be used as academical benchmark for numerical simulations.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Chronic mitral regurgitation and Doppler estimation of left ventricular filling pressures in patients with heart failure

Previous studies relating Doppler parameters and pulmonary capillary wedge pressures (PCWP) typically exclude patients with severe mitral regurgitation (MR). We evaluated the effects of varying degrees of chronic MR on the Doppler estimation of PCWP. PCWP and mitral Doppler profiles were obtained in 88 patients (mean age 55 +/- 8 years) with severe left ventricular (LV) dysfunction (mean ejection fraction 23% +/- 5%). Patients were classified by severity of MR. Patients with severe MR had greater left atrial areas, LV end-diastolic volumes, and mean PCWPs and lower ejection fractions (each P <.01). In patients with mild MR, multiple echocardiographic parameters correlated with PCWP; however, with worsening MR, only deceleration time strongly related to PCWP. From stepwise multivariate analysis, deceleration time was the best independent predictor of PCWP overall, and it was the only predictor in patients with moderate or severe MR. Doppler-derived early mitral deceleration time reliably predicts PCWP in patients with severe LV dysfunction irrespective of degree of MR.

NASA Discipline Cardiopulmonary↗
DOE OSTI2023

Temperature effects in plasma-based positron acceleration schemes using electron filaments

Preserving the quality of a positron beam in a plasma-based accelerator, where a wakefield suitable for positron transport and acceleration is generated by means of an electron filament, is challenging. This is due to the nature of the wakefields, characterized by focusing fields that vary nonlinearly in the transverse direction, and by accelerating fields that are non-uniform. These fields also change slice-by-slice along the beam. Maintaining a high beam quality is pivotal for application of positron beams in a plasma-based collider. In this paper, we show that an initial background plasma temperature can help mitigate the positron beam quality degradation in plasma-based accelerators that rely on electron filaments. We show that temperature effects broaden the electron filament and smooth radially both the non-linear transverse and the non-uniform longitudinal wakefields. Using warm plasmas opens up new possibilities to improve beam quality in several plasma-based positron acceleration concepts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1999

Combustion of a Polymer (PMMA) Sphere in Microgravity

A series of low gravity, aircraft-based, experiments was conducted to investigate the combustion of supported thermoplastic polymer spheres under varying ambient conditions. The three types of thermoplastic investigated were polymethylmethacrylate (PMMA), polypropylene (PP). and polystyrene (PS). Spheres with diameters ranging from 2 mm to 6.35 mm were tested. The total initial pressure varied from 0.05 MPa to 0. 15 MPa whereas the ambient oxygen concentration varied from 19 % to 30 % (by volume). The ignition system consisted of a pair of retractable energized coils. Two CCD cameras recorded the burning histories of the spheres. The video sequences revealed a number of dynamic events including bubbling and sputtering, as well as soot shell formation and break-up during combustion of the spheres at reduced gravity. The ejection of combusting material from the burning spheres represents a fire hazard that must be considered at reduced gravity. The ejection process was found to be sensitive to polymer type. All average burning rates were measured to increase with initial sphere diameter and oxygen concentration, whereas the initial pressure had little effect. The three thermoplastic types exhibited different burning characteristics. For the same initial conditions, the burning rate of PP was slower than PMMA, whereas the burning rate of PS was comparable to PMMA. The transient diameter of the burning thermoplastic exhibited two distinct periods: an initial period (enduring approximately half of the total burn duration) when the diameter remained approximately constant, and a final period when the square of the diameter linearly decreased with time. A simple homogeneous two-phase model was developed to understand the changing diameter of the burning sphere. Its value is based on a competition between diameter reduction due to mass loss from burning and sputtering, and diameter expansion due to the processes of swelling (density decrease with heating) and bubble growth. The model relies on empirical parameters for input, such as the burning rate and the duration of the initial and final burning periods.

Yang, Jiann C.↗
NASA NTRS1959

Investigation at Mach Numbers of 0.20 to 3.50 of a Blended Diamond Wing and Body Combination of Sonic Design but with Low Wave-Drag Increase with Increasing Mach Number

A diamond wing and body combination was designed to have an area distribution which would result in near optimum zero-lift wave-drag coefficients at a Mach number of 1.00, and decreasing wave-drag coefficient with increasing Mach number up to near sonic leading-edge conditions for the wing. The airfoil section were computed by varying their shape along with the body radii (blending process) to match the selected area distribution and the given plan form. The exposed wing section had an average maximum thickness of about 3 percent of the local chords, and the maximum thickness of the center-line chord was 5.49 percent. The wing had an aspect ratio of 2 and a leading-edge sweep of 45 deg. Test data were obtained throughout the Mach number range from 0.20 to 3.50 at Reynolds numbers based on the mean aerodynamic chord of roughly 6,000,000 to 9,000,000. The zero-lift wave-drag coefficients of the diamond model satisfied the design objectives and were equal to the low values for the Mach number 1.00 equivalent body up to the limit of the transonic tests. From the peak drag coefficient near M = 1.00 there was a gradual decrease in wave-drag coefficient up to M = 1.20. Above sonic leading-edge conditions of the wing there was a rise in the wave-drag coefficient which was attributed in part to the body contouring as well as to the wing geometry. The diamond model had good lift characteristics, in spite of the prediction from low-aspect-ratio theory that the rear half of the diamond wing would carry little lift. The experimental lift-curve slope obtained at supersonic speeds were equal to or greater than the values predicted by linear theory. Similarly the other basic aerodynamic parameters, aerodynamic center position, and maximum lift-drag ratios were satisfactorily predicted at supersonic speeds.

Holdaway, George H.↗
DOE OSTI2020

Study of a linear surface wave plasma source for tin removal in an extreme ultraviolet source

Tin deposition mitigation employs hydrogen radicals and ions, formed in a hydrogen plasma, to interact with tin to form tin hydride (SnH 4 ) in the gaseous state, which is then pumped away. Surface wave plasma (SWP) technology developed at Illinois generates hydrogen radicals and ions, resulting in tin etch rates that are high enough to keep extreme ultraviolet (EUV) lithographic tools clean. An advantage of an SWP antenna is the ability to generate a high density of hydrogen radicals and hydrogen ions directly at the desired etching location. In situ etching of tin enables high availability EUV tools by maintaining high reflectivity of the multilayer mirror of the collector. Additionally, the SWP is characterized with low ion energies and low electron temperature, such that the multilayer mirror does not suffer any damage from sputtering or implantation of hydrogen ions during operation. Here, experiments elucidating the fundamental processes of tin removal are conducted by varying pressure, power, surface temperature, and gas flow rate in order to observe the etch rate behavior. Our results have shown that the presence of hydrogen ions increases etch rates because ion bombardment weakens Sn–Sn bonds, which, in turn, allows for a higher rate of chemical etching by the radicals. The ion bombardment reduces the number of radicals needed to etch a single tin atom to the range of 10 2 –10 3 . The linear SWP antenna yields plasma densities on the order of 10 17 to 10 18 m -3 and radical densities on the order of 10 18 to 10 19 m -3 , allowing for greater utilization of ion etch enhancement. Etch rates of up to 200 nm/min have been achieved. The surface temperature of the samples is an important factor in the etching process such that the decrease in the surface temperature increases the etch rates and decreases the hydrogen desorption rates. In addition, a kinetic etch model is developed to explain the behavior of etch rates as a function of surface temperature. Furthermore, results from experiments performed in an Illinois NXE:3100 chamber will be discussed.

42 ENGINEERING↗
NASA NTRS1994

Reionization in a cold dark matter universe: The feedback of galaxy formation on the intergalactic medium

We study the coupled evolution of the intergalactic medium (IGM) and the emerging structure in the universe in the context of the cold dark matter (CDM) model, with a special focus on the consequences of imposing reionization and the Gunn-Peterson constraint as a boundary condition on the model. We have calculated the time-varying density of the IGM by coupling our detailed, numerical calculations of the thermal and ionization balance and radiative transfer in a uniform, spatially averaged IGM of H and He, including the mean opacity of an evolving distribution of gas clumps which correspond to quasar absorption line clouds, to the linearized equations for the growth of density fluctuations in both the gaseous and dark matter components in a CDM universe. We use the linear growth equations to identify the fraction of the gas which must have collapsed out at each epoch, an approach similar in spirit to the so-called Press-Schechter formalism. We identify the IGM density with the uncollapsed baryon fraction. The collapsed fraction is postulated to be a source of energy injection into the IGM, by radiation or bulk hydrodynamical heating (e.g., via shocks) or both, at a rate which is marginally enough to satisfy the Gunn-Peterson constraint at z less than 5. Our results include the following: (1) We find that the IGM in a CDM model must have contained a substantial fraction of the total baryon density of the universe both during and after its reionization epoch. (2) As a result, our previous conclusion that the observed Quasi-Stellar Objects (QSOs) at high redshift are not sufficient to ionize the IGM enough to satisfy the Gunn-Peterson constraint is confirmed. (3) We predict a detectable He II Gunn-Peterson effect at 304(1 + z) A in the spectra of quasars at a range of redshift z greater than or approx. 3, depending on the nature of the sources of IGM reionization. (4) We find, moreover, that a CDM model with high bias parameter b (i.e., b greater than or approx. 2) cannot account for the baryon content of the universe at z approximately 3 observed in quasar absorption line gas unless Omega (sub B) significantly exceeds the maximum value allowed by big bang nucleocynthesis. (5) For a CDM model with bias parameter within the allowed range of (lower) values, the lower limit to Omega(sub B) imposed by big bang nucleosynthesis (Omega(sub B) h(sup 2) greater than or equal to 0.01) combines with our results to yield the minimum IGM density for the CDM fodel. For CDM with b = 1 (Cosmic Background Explorer (COBE) normalization), we find Omega(sub IGM)(sup min) (z approximately 4) approx. equal 0.02-0.03, and Omega(sub IGM)(sup min)(z approximately 0) approx. equal 0.005-0.03, depending upon the nature of the sources of IGM reionization. (6) In general, we find that self-consistent reionization of the IGM by the collapsed baryon fraction has a strong effect on the rate of collapse. (7) As a further example, we show that the feedback effect on the IGM of energy release by the collapsed baryon fraction may explain the slow evolution of the observed comoving QSO number density between z = 5 and z = 2, followed by the sharp decline after z = 2.

Shapiro, Paul R.↗
DOE OSTI2022

Methods for Computing Physically Realistic Estimates of Electric Water Heater Demand Response Resource Suitable for Bulk Power System Planning Models

Demand response is commonly called on to reduce load during system peak times or to respond to contingency events. In future power systems with higher shares of wind and solar generation (which we describe together as variable generation [VG]), demand response could have more opportunities to provide energy shifting or operating reserve services. This report evaluates the ability of residential electric water heaters, both electric resistance water heaters (ERWHs) and heat pump water heaters (HPWHs), to provide such services starting from detailed whole-building energy models that realistically represent New England single family home stock. We use a parsimonious surrogate model to represent operational flexibility in a form suitable for linear and mixed integer programming. This enables relatively fast determination of aggregate contingency reserve resource, price-taking energy shifting outcomes, and in some cases the determination of aggregate models at the megawatt (MW) scale that can be directly included in large-scale grid models. After selecting modeling methods and parameters through various computational experiments, we find interquartile ranges of contingency reserve resource in ISO-NE for about 603,400 ERWHs of 45 MW - 69 MW for Claim10 (50 minute responses provided with 10 minutes of advanced notification) and 65 MW - 102 MW for Claim30 (30 minute responses provided with 30 minutes of advanced notification), and for about 619,000 HPWHs of 48 MW - 88 MW for Claim10 and 52 MW - 90 MW for Claim30. The overall reserve resource is up to 32% of total load for ERWHs providing Claim10 service, 47% for ERWHs providing Claim30 service, 93% for HPWHs providing Claim10 service, and 97% for HPWHs providing Claim30 service. More work is required to determine if HPWHs are inherently more suitable than ERWHs for providing contingency reserve or if these results reflect idiosyncrasies of the single family home stock model used in this study. The value of this contingency resource in a Near-term VG model of ISO-NE is $\$ 0.40$ to $\$1.20$ per water heater-year, and significantly larger, $\$ 3.80$ to $\$ 5.30$ per water heater-year in a Mid-term VG model of ISONE. Aggregating surrogate models to the MW-scale for energy shifting service is more challenging than for contingency service and we only present such results for ERWHs, because we were unable to determine satisfactory ways to deal with HPWHs' time-varying and path dependent operational characteristics. Individual surrogate models suitable for evaluating the energy shifting resource from both ERWHs and HPWHs are created, however, and dispatched against day-ahead prices from the Near-Term VG and Mid-Term VG models of ISO-NE. The individual surrogate models are able to access and potentially shift all 640 GWh of HPWH load and 1,547 GWh of ERWH load we modeled in two different single family home stock models. In contrast, the most effective model of aggregate ERWH shifting resource we created only captured 34.7% of the total ERWH load. Energy shifting affected by price-taking dispatch against modeled day-ahead energy prices produces per water heater year profits of $\$19.44$ - $\$22.93$ for individual HPWHs, $\$39.11$ - $\$40.54$ for individual ERWHs, and up to $\$4.00$ - $\$4.24$ for aggregated ERWHs, with the variations mainly due to grid conditions (more or less VG). When the supply-side response to these changes is accounted for, the per water heater year production cost savings for ISO-NE are $\$7.50$ to $\$17.70$ for the most effective set of endogenously dispatched aggregate ERWHs, $\$15.60$ to $\$15.70$ for individual ERWHs dispatched against the DA prices, and $\$10.70$ to $\$11.20$ for individual HPWHs dispatched against DA prices. Those ranges primarily represent the difference between Near-Term VG and Mid-Term VG grid conditions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Analysis of heat transfer and AuNPs-mediated photo-thermal inactivation of E. coli at varying laser powers using single-phase CFD modeling

In the wake of the COVID-19 pandemics, the demand for innovative and effective methods of bacterial inactivation has become a critical area of research, providing the impetus for this study. The purpose of this research is to analyze the AuNPs-mediated photothermal inactivation of E. coli. Gold nanoparticles irradiated by laser represent a promising technique for combating bacterial infection that combines high-tech and scientific progress. The intermediate aim of the work was to present the calibration of the model with respect to the gold nanorods experiment. The purpose of this work is to study the effect of initial concentration of E. coli bacteria, the design of the chamber and the laser power on heat transfer and inactivation of E. coli bacteria. Using the CFD simulation, the work combines three main concepts. 1. The conversion of laser light to heat has been described by a combination of three distinctive approximations: a- Discrete particle integration to take into account every nanoparticle within the system, b- Rayleigh-Drude approximation to determine the scattering and extinction coefficients and c- Lambert–Beer–Bourger law to describe the decrease in laser intensity across the AuNPs. 2. The contribution of the presence of E. coli bacteria to the thermal and fluid-dynamic fields in the microdevice was modeled by single-phase approach by determining the effective thermophysical properties of the water-bacteria mixture. 3. An approach based on a temperature threshold attained at which bacteria will be inactivated, has been used to predict bacterial response to temperature increases. The comparison of the thermal fields and temporal temperature changes obtained by the CFD simulation with those obtained experimentally confirms the accuracy of the light-heat conversion model derived from the aforementioned approximations. The results show a linear relationship between maximum temperature and variation in laser power over the range studied, which is in line with previous experimental results. It was also found that the temperature inside the microchamber can exceed 55 °C only when a laser power higher than 0.8 W is used, so bacterial inactivation begins. The experimental data allows to determinate the concentration of nanoparticles. This parameter is introduced into the mathematical model obtaining the same number of AuNPs. However, this assumption introduces a certain simplification, as in the mathematical model the distribution of nanoparticles is uniform. This work is directly connected to the use of gold nanoparticles for energy conversion, as well as the field of bacterial inactivation in microfluidic systems such as lab-on-a-chip. Presented mathematical and numerical models can be extended to the entire spectrum of wavelengths with particular use of white light in the inactivation of bacteria. This work represents a significant advancement in the field, as to the best of the authors’ knowledge, it is the first to employ a single-phase computational fluid dynamics (CFD) approach specifically combined with the thermal inactivation of bacteria. Moreover, this research pioneers the use of a numerical simulation to analyze the temperature threshold of photothermal inactivation of E. coli mediated by gold nanorods (AuNRs). The integration of these methodologies offers a new perspective on optimizing bacterial inactivation techniques, making this study a valuable contribution to both computational modeling and biomedical applications.

36 MATERIALS SCIENCE↗
NASA NTRS2016

Feasibility of Varying Geo-Fence Around an Unmanned Aircraft Operation Based on Vehicle Performance and Wind

Managing trajectory separation is critical to ensuring accessibility, efficiency, and safety in the unmanned airspace. The notion of geo-fences is an emerging concept, where distance buffers enclose individual trajectories and areas of operation in order to manage the airspace. Currently, the Air Traffic Management system for commercial travel defines static distance buffers around the aircraft; however, commercial UASs are envisioned to operate in significantly closer proximity to other UAS requiring a geo-fence for spacing operations. The geo-fence size can be determined based on vehicle performance characteristics, state of the airspace, weather, and other unforeseen events such as emergency or disaster response. Calculation of the geo-fence size could be determined as part of pre-flight planning and during real-time operations. A largely non-homogeneous fleet of UASs will be operating in low altitude and will likely be commercially developed. Due to intellectual property concerns, the operators may not provide detailed specifications of the control system to UTM. In addition, the huge variety of UAS makes modeling each control system prohibitive and flight data for these vehicles may not exist. Therefore, a generalized, simple geo-fence sizing algorithm must be developed such that it does not rely on detailed knowledge of the vehicle control system, accounts for the presence of urban winds, and is sufficiently accurate. In this work, two simple models are investigated to determine its feasibility as an adequate means for calculating the geo-fence size. The vehicle data used in this work are provided by UAS manufactures who have partnered with NASA's UTM project and some publicly available websites. The first model utilizes wind data processed from the NOAA HRRR (Hourly Rapid Refresh) product and Sonar Annemometer data provided by San Jose State. The second model utilizes OpenFOAM which is a CFD code used to generate a wind field for flow around a single building. The key vehicle performance parameters can include UAS response time to disturbances, command to actuation latency, control system rate limits, time to recovery to desired path, and aerodynamics. It was found that the first model provides an initial understanding of geo-fence sizing, but does not provide enough accuracy to provide UTM with an efficient means of scheduling vehicles. The results of the second model reveal that modeling UAS controls systems with a linearized plant and gain scheduled PID controller does not allow capture the UAS flight dynamics within a significant envelope of the wind disturbances.

uncertainty↗
NASA NTRS2024

Global Sensitivity Analysis of Simulated Remote Sensing Polarimetric Observations Over Snow

This study presents a detailed theoretical assessment of the information content of passive polarimetric observations over snow scenes, using a global sensitivity analysis (GSA) method. Conventional sensitivity studies focus on varying a single parameter while keeping all other parameters fixed. In contrast, the GSA correctly addresses the covariance of state parameters across their entire parameter space, hence favoring a more correct interpretation of inversion algorithms and the optimal design of their state vectors. The forward simulations exploit a vector radiative transfer model to obtain the Stokes vector emerging at the top of the atmosphere for different solar zenith angles, when the bottom boundary consists of a vertically resolved snowpack of non-spherical grains. The presence of light-absorbing particulates (LAPs), either embedded in the snow or aloft in the atmosphere above in the form of aerosols, is also considered. The results are presented for a set of wavelengths spanning the visible (VIS), near-infrared (NIR), and shortwave infrared (SWIR) region of the spectrum. The GSA correctly captures the expected, high sensitivity of the reflectance to LAPs in the VIS–NIR and to grain size at different depths in the snowpack in the NIR–SWIR. With adequate viewing geometries, mono-angle measurements of total reflectance in the VIS–SWIR (akin to those of the Moderate Resolution Imaging Spectroradiometer, MODIS) resolve grain size in the top layer of the snowpack sufficiently well. The addition of multi-angle polarimetric observations in the VIS–NIR provides information on grain shape and microscale roughness. The simultaneous sensitivity in the VIS–NIR to both aerosols and snow-embedded impurities can be disentangled by extending the spectral range to the SWIR, which contains information on aerosol optical depth while remaining essentially unaffected when the same particulates are mixed with the snow. Multi-angle polarimetric observations can therefore (i) effectively partition LAPs between the atmosphere and the surface, which represents a notorious challenge for snow remote sensing based on measurements of total reflectance only and (ii) lead to better estimates of grain shape and roughness and, in turn, the asymmetry parameter, which is critical for the determination of albedo. The retrieval uncertainties are minimized when the degree of linear polarization is used in place of the polarized reflectance. The Sobol indices, which are the main metric for the GSA, were used to select the state parameters in retrievals performed on data simulated for multiple instrument configurations. Improvements in retrieval quality with the addition of measurements of polarization, multi-angle views, and different spectral channels reflect the information content, identified by the Sobol indices, relative to each configuration. The results encourage the development of new remote sensing algorithms that fully leverage multi-angle and polarimetric capabilities of modern remote sensors. They can also aid flight planning activities, since the optimal exploitation of the information content of multi-angle measurements depends on the viewing geometry. The better characterization of surface and atmospheric parameters in snow-covered regions advances research opportunities for scientists of the cryosphere and ultimately benefits albedo estimates in climate models.

remote sensing↗
DOE OSTI2025

Efficient Measurement of Length Distribution of 1D Nanoparticles in Solution via Optical Polarimetry

The efficient measurement of the length distribution of nanotubes, nanowires, and other one-dimensional (1D) nanoparticles in solution is important to enable their incorporation into materials and devices and to optimize their processing for properties of interest, such as thermal/electrical conductivity or mechanical strength, in suspensions and composites. We report an electric-field (E-field)-assisted optical-polarimetry technique to measure the length distribution of ensembles of high-aspect-ratio particles in dilute suspension. The degree of alignment of polarizable 1D particles suspended in a fluid under Brownian motion explicitly depends on the E-field strength and the particle length. We show that it is possible to extract the length distribution of 1D nanoparticles suspended in an insulating fluid by applying a range of E-fields and using optical polarimetry to measure the corresponding alignment order parameter. Notably, the method is relatively insensitive to the diameter of the 1D particles, which can be poorly known or vary within a sample. The technique is validated with silver nanowires and carbon nanotubes of known lengths, as well as polymer-depletion-length-separated single-wall carbon nanotube samples with length distributions independently measured with analytical ultracentrifugation. Finally, we demonstrate the ability of the optical-polarimetry technique to quantify changes in the length distribution of ultranarrow, sub-nanometer-diameter single-wall carbon nanotubes under different types and durations of ultrasonication. Within its range of applicability (polarizable 1D nanoparticles in the 0.5 to 15 μm length range, constrained by the voltage stability of the media and the suspended particles), the E-field-assisted optical-polarimetry method is a particularly efficient and accurate method to measure the length distribution of nanowires and nanotubes in suspension.

1D nanoparticles↗
NASA NTRS1993

Relationships between pigment composition variation and reflectance for plant species from a coastal savannah in California

Advances in imaging spectroscopy have indicated that remotely sensed reflectance measurements of the plant canopy may be used to identify and qualify some classes of canopy biochemicals; however, the manner in which differences in biochemical compositions translate into differences is not well understood. Most frequently, multiple linear regression routines have been used to correlate narrow band reflectance values with measured biochemical concentrations. Although some success has been achieved with such methods for given data sets, the bands selected by multiple regression are not consistent between data sets, nor is it always clear what physical or biological basis underlies the correlation. To examine the relationship between biochemical concentration and leaf reflectance signal we chose to focus on the visible spectrum where the primary biochemical absorbances are due to photosynthetic pigments. Pigments provide a range of absorbance features, occur over a range of concentrations in natural samples, and are ecophysiologically important. Concentrations of chlorophyll, for example, have been strongly correlated to foliar nitrogen levels within a species and to photosynthetic capacity across many species. In addition pigments effectively absorb most of the photosynthetically active radiation between 400-700 nm, a spectral region for which silicon detectors have good signal/noise characteristics. Our strategy has been to sample a variety of naturally occurring species to measure leaf reflectance and pigment compositions. We hope to extend our understanding of pigment reflectance effects to interpret small overlapping absorbances of other biochemicals in the infrared region. For this reason, selected samples were also tested to determine total nitrogen, crude protein, cellulose, and lignin levels. Leaf reflectance spectra measured with AVIRIS bandwidths and wavelengths were compared between species and within species and for differences between seasons, for changes in the the shape of the spectra. We attempt to statistically correlate these shape changes with differences in pigment compositions. In parallel with our comparisons of pigment composition and leaf reflectance, we have modified the PROSPECT leaf reflectance model to test the contributions of pigments or pigment group concentrations. PROSPECT considers a leaf as a multi-layer dielectric plane with an uneven surface. Jacquemoud adapted the basic analysis of Allen for surface effects, a leaf thickness factor, and the absorption of water and chlorophyll (actually all pigments) and the plant matrix. Our modifications to PROSPECT in the forward direction include breaking out the pigment concentration parameter into separate components for chlorophyll a and b and a number of xanthophylls and carotenes, and introducing a shift and convolution function to model the spread and shift from their in vitro measurements to their in vivo state. Further, we have considered how the matrix elements (i.e., all biochemicals and structural effects not modeled explicity) vary with species.

Ustin, Susan L.↗
NASA NTRS2018

Implementation of J-A Methodology Elastic-Plastic Crack Instability Analysis Capability into the WARP-3D Code

Characterization of the near crack-tip stress/strain fields is the foundation of fracture mechanics. The description of the near tip stress field and the prediction of when fracture occurs is well established for brittle materials that exhibit linear elastic behavior. However, in ductile materials or conditions that violate linear elastic assumptions (Aluminum alloys, Al 2024-T3, Al 2024- T351 etc.), the elastic-plastic crack-tip stress fields are characterized by the Hutchison-Rice-Rosengren (HRR) field. The J-integral is commonly used to characterize amplitude of the HRR field under elastic-plastic conditions. The J-integral has been demonstrated for crack-tip fields that are under high constraint conditions (i.e., small-scale plasticity where the J-dominance is maintained). However, as the external load increases, yielding changes from small- to largescale plasticity and usually a loss of constraint (i.e., reduction in the triaxial stress field along the crack front). The loss of constraint leads to the deviation of the crack-tip stress fields from that given by the HRR field. Hence, the J-dominance will be gradually lost and additional parameter(s) are required to quantify the crack-tip stress fields and predict fracture behavior. The assessment objectives were to: 1) implement a two-parameter (i.e., J-A) fracture criterion into an elastic-plastic three-dimensional (3D) finite element analysis (FEA), 2) validate the implementation by comparison with the A parameter from literature data, 3) conduct material characterization tests to quantify the material behavior and provide fracture data for validation of the J-A fracture criteria, and (4) perform evaluations to establish if the J-A criteria can be used to predict fracture in a ductile metallic material (e.g., aluminum alloys). The A parameter in these criteria is the second parameter in a three-term elastic-plastic asymptotic expansion of the neartip stress behavior. A series of extensive FEAs were performed using WARP3D software package to obtain solutions for the A parameter for different specimen configurations. The methodology needed for the estimation of the A parameter in the asymptotic expansion was developed and implemented using Matlab®. A user material (UMAT) routine was used to model the material stress-strain response using a Ramberg-Osgood power law with a hardening exponent (n) and a material coefficient (alpha). This UMAT routine was successfully implemented in WARP3D software and validated through comparison with the experimental data. Three configurations were extracted from published results: 1) center cracked plate (CCP), 2) single edge-cracked plate (SECP), and 3) double edge-cracked plate (DECP). These configurations and four other configurations (three-hole tension (THT)), three-point bend (3PTB), three-hole compact tension (3PCT), and compact tension (CT)) were analyzed to verify the methodology that was developed and implemented into WARP3D. Solutions of the A parameter were obtained for remote tension loading conditions that started with small-scale yielding and continued into the large-scale plasticity regime. The results indicate that the methodology developed can be used to calculate the elastic-plastic J-A parameters for test specimens with a range of crack geometries, material strain hardening behaviors, and loading conditions. The J-A parameters were implemented as fracture criteria and used to predict the test results. For comparison, other fracture criteria were used to predict the same test results. Major findings include: The A constraint parameter A varies with specimen type and applied load thus accurate determination is crucial in predicting the failure load, and the A parameter is asymptotic as the failure load is approached, making an accurate determination difficult (i.e., small differences in the A parameter can cause large variations in failure load) for materials exhibiting elastic-plastic behavior. The failure predictions from J-A methodology were more accurate than the traditionally used KC and J methods, and have comparable scatter to that observed when using the crack-tip opening angle (CTOA) method. However, the J-A methodology requires considerable effort (expertise level and labor) to implement and to evaluate the A parameter for different specimen types and materials, or to apply this methodology to part-through crack (e.g., 3D problems) structural applications.

Hamm, Kenneth R., Jr.↗
DOE OSTI2025

A comparison of effective field theory models of redshift space galaxy power spectra for DESI 2024 and future surveys

In preparation for the next generation of galaxy redshift surveys, and in particular the year-one data release from the Dark Energy Spectroscopic Instrument (DESI), we investigate the consistency of a variety of effective field theory models that describe the galaxy-galaxy power spectra in redshift space into the quasi-linear regime using 1-loop perturbation theory. These models are employed in the pipelines velocileptors, PyBird, and Folpsν. While these models have been validated independently, a detailed comparison with consistent choices has not been attempted. After briefly discussing the theoretical differences between the models we describe how to provide a more apples-to-apples comparison between them. We present the results of fitting mock spectra from the AbacusSummit suite of N-body simulations provided in three redshift bins to mimic the types of dark time tracers targeted by the DESI survey. We show that the theories behave similarly and give consistent constraints in both the forward-modeling and ShapeFit compressed fitting approaches. We additionally generate (noiseless) synthetic data from each pipeline to be fit by the others, varying the scale cuts in order to show that the models agree within the range of scales for which we expect 1-loop perturbation theory to be applicable. Finally, this work lays the foundation of Full-Shape analysis with DESI Y1 galaxy samples where in the tests we performed, we found no systematic error associated with the modeling of the galaxy redshift space power spectrum for this volume.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Dark Energy Survey Year 3 results: Constraints on extensions to Λ CDM with weak lensing and galaxy clustering

We constrain six possible extensions to the Λ cold dark matter (CDM) model using measurements from the Dark Energy Survey’s first three years of observations, alone and in combination with external cosmological probes. The DES data are the two-point correlation functions of weak gravitational lensing, galaxy clustering, and their cross-correlation. We use simulated data vectors and blind analyses of real data to validate the robustness of our results to astrophysical and modeling systematic errors. In many cases, constraining power is limited by the absence of theoretical predictions beyond the linear regime that are reliable at our required precision. The Λ CDM extensions are dark energy with a time-dependent equation of state, nonzero spatial curvature, additional relativistic degrees of freedom, sterile neutrinos with eV-scale mass, modifications of gravitational physics, and a binned σ 8 ( z ) model which serves as a phenomenological probe of structure growth. For the time-varying dark energy equation of state evaluated at the pivot redshift we find ( w p , w a ) = ( - 0.9 9 - 0.17 + 0.28 , - 0.9 ± 1.2 ) at 68% confidence with z p = 0.24 from the DES measurements alone, and ( w p , w a ) = ( - 1.0 3 - 0.03 + 0.04 , - 0. 4 - 0.3 + 0.4 ) with z p = 0.21 for the combination of all data considered. Curvature constraints of Ω k = 0.0009 ± 0.0017 and effective relativistic species N eff = 3.1 0 - 0.16 + 0.15 are dominated by external data, though adding DES information to external low-redshift probes tightens the Ω k constraints that can be made without cosmic microwave background observables by 20%. For massive sterile neutrinos, DES combined with external data improves the upper bound on the mass m eff by a factor of 3 compared to previous analyses, giving 95% limits of ( Δ N eff , m eff ) ≤ ( 0.28 , 0.20 eV ) when using priors matching a comparable Planck analysis. For modified gravity, we constrain changes to the lensing and Poisson equations controlled by functions Σ ( k , z ) = Σ 0 Ω Λ ( z ) / Ω Λ , 0 and μ ( k , z ) = μ 0 Ω Λ ( z ) / Ω Λ , 0 , respectively, to Σ 0 = 0. 6 - 0.5 + 0.4 from DES alone and ( Σ 0 , μ 0 ) = ( 0.04 ± 0.05 , 0.0 8 - 0.19 + 0.21 ) for the combination of all data, both at 68% confidence. Overall, we find no significant evidence for physics beyond Λ CDM .

79 ASTRONOMY AND ASTROPHYSICS↗