SEARCH · Engineering Papers

Results for “Reproducibility characterization”

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 811 records · Page 45

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Development of an Advanced Clothing Model that Considers Heat and Moisture Transport

Despite recent advancements in the development of human thermal physiology models, clothing is often considered in only a rudimentary manner, that is, as a boundary condition at the skin characterized by its bulk thermal and evaporative impedances. While such an approach may be sufficient for thermoregulation modeling of humans operating in relatively low thermal burden scenarios, e.g., within an office environment, a clothing thermal model that does not explicitly consider the layer-to-layer transport of vapor can yield incorrect results, especially under conditions corresponding to high thermal stress or extreme environments. For example, vapor diffusing through layers of clothing can reach its dew point, in which case a phase change will occur that can fundamentally affect the transport of moisture and heat within the clothing ensemble. Specifically, the evaporation and condensation of moisture within the fabric layers can affect liquid content and temperature and, ultimately, affect thermo-physiological response and thermal comfort perception. Development of an advanced clothing model that considers heat and moisture transport within the fabric layers can improve the accuracy of human thermal predictions in scenarios where clothing-based condensation and evaporation occur. This paper describes how a general heat transfer software modeling tool was modified to consider the transport of vapor, the storage of liquid, and the thermal effects of phase change within a clothing ensemble. The advanced clothing moisture model was tested by comparing its predictions to measurements derived from sweating hot plate experiments and human physiological experiments. The model was able to reproduce the thermal and moisture results from tests performed under conditions in which condensation within the clothing layers was known to occur. This study demonstrates the importance of considering layer-to-layer vapor transport in clothing models and its impact on human thermal modeling predictions.

Human thermal modeling↗

NASA NTRS1993

Models of the Cartwheel ring galaxy: Spokes and starbursts

Recent observations of this famous ring galaxy, including optical and near-infrared CCD surface photometry, and VLA radio continuum and 21 cm line mapping (Higdon 1992b, in prep.), have inspired a renewed modeling effort. Toomre's (1978, in The Large-scale Structure of the Universe, eds. Longair and Einasto) series of restricted three-body simulations demonstrated how the multiple rings could be produced in a nearly head-on galaxy collision. New models with a halo-dominated potential based on the 21 cm rotation curve are able to reproduce such details as the spacing between rings, ring widths, offset of the nucleus, and several kinematical features, thus providing strong support for the collisional theory. The new observations have shown there are little or no old stars in Cartwheel; it may consist almost entirely of gas and stars produced as a result of compression in the ring wave. To model this process Smooth Particle Hydrodynamics (SPH) simulations of the Cartwheel disk have been performed. Fixed gravitational potentials were used to represent the Cartwheel and a roughly 30 percent mass collision partner. The interaction dynamics was treated as in the usual restricted three-body approximation, and the effects of local self-gravity between disk particles were calculated. We are particularly interested in testing the theory that enhanced star formation in waves is the result of gravitational instability in the compressed region (see e.g. Kennicutt 1989, ApJ 344, 685). The gas surface density in a number of simulations was initialized to a value slightly below the threshold for local gravitational instability throughout most of the disk. The first ring wave produces relatively modest compressions (a factor of order a few), triggering instability in a narrow range of wavelengths. Self-gravity in the disk is calculated over a comparable range of scales. Simulations were run with isothermal, adiabatic, and adiabatic with radiative cooling characterized by a relatively short timescale. The isothermal approximation is good except in the vicinity of the strong second (inner) ring, and several snapshots from one case are shown in the figure below. Flocculent spiral segments are present before the collision, and these are compressed into dense knots in the ring wave. These knots are likely to be sites of vigorous star formation. In the strong rarefaction behind the outer ring most of the knots are radially stretched and sheared, giving rise to spoke-like features. A few dense knots are evidently very tightly bound, because they retain their coherence and are stretched relatively little through the rarefaction. This is in accord with evidence for continuing star formation in some spokes (Marcum, Appleton and Higdon 1992). The number and spacing of spokes is a direct function of the scale of the gravitational instability in the disk. Thus, the gravitational instability theory, together with the hypothesis that massive stars are only formed in dense knots of gas, can account for most of the distinct morphology of the Cartwheel.

Struck-Marcell, Curtis↗

DOE OSTI2024

Analysis of AP1000 Small-Break Loss-of-Coolant Accident Using Reactor Transient Simulator

The Westinghouse Electric Company’s Advanced Passive Reactor (AP1000) is characterized by the incorporation of passive safety systems (PSSs) designed to ensure core cooling during transient events. The assessment of PSSs requires evaluation of their performance through a combination of experiments and simulations employing various thermal-hydraulic codes. In addition, detailed evaluation of PSSs for a specific reactor system transient analysis such as loss-of-coolant-accident analysis supports understanding representative integral effects test facility development and the further evolution model development and assessment process. Developing a reactor system code is a complex and time-consuming process that requires significant engineering expertise and effort. It can take several months to even years to complete in the early stages of reactor system design and analysis. However, this process can be expedited through the use of transient simulator models for similar reactor systems, which can be used for lesson learning and training purposes. This study uses the Personal Computer Transient Analyzer (PCTRAN) code. The main advantage of PCTRAN is its ease of use and ability to run faster than real time. This study presents the results obtained for a small-break loss-of-coolant accident (SBLOCA) for two breaks using the full version (licensed) of PCTRAN. The purpose of this investigation is to evaluate the overall system behavior during the postulated SBLOCA event as well as assess the capability of the PCTRAN code to reproduce the system response during transient events. The obtained results were compared with the Westinghouse NOTRUMP system code. The PCTRAN code proved to be reliable in predicting the qualitative behavior of the system in both transient cases. As for the system response, it was found that it is contingent on the activation time of the PSSs. The differences in reactor coolant system pressure between the two codes were attributed to the critical flow model and simplification of mass and energy balance. Despite PCTRAN’s limitations, it can still provide a reasonable prediction of various reactor parameters such as pressure, mass flow rate, and void fraction during a SBLOCA scenario. It is worth noting that PCTRAN currently employs a bulk approach similar to that of the Modular Accident Analysis Program (MAAP) and MELCOR codes. However, the upcoming version of PCTRAN will include an artificial intelligence–based detection and accident prevention system, as well as different models for different reactor components. Consequently, PCTRAN has the potential to be upgraded to match the system thermal-hydraulic codes of the U.S. Nuclear Regulatory Commission and become more widely used in cybersecurity to safeguard nuclear power plants from cyberattacks.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Vortex Structure of Head Bubble in Convective Cloud Starting Plume

Abstract A developing cumulus cloud (Cu) was modeled, and dynamic, thermodynamic, and microphysical properties of an ascending head bubble reproducing the upper part of a developing Cu were investigated. The data for analysis are taken from 10-m-resolution LES of trade wind Cu under BOMEX conditions. The detection of a rising bubble is carried out using wavelet filtering of the velocity fields and microphysical fields, while a low-frequency signal of the filtering is associated with the convective-scale structure of cloud. We substantiate and discuss the representation of the bubble as a vortex ring, and estimate the parameters of this vortex ring. The simplest Hill’s vortex was chosen as a model of a vortex ring inside cloud. Analytical approximations of the radial profiles of the vertical velocity and of conservative quantities (such as total water mixing ratio and liquid water potential temperature inside and outside the bubble) are obtained. The spatial structure of these quantities is investigated using analytical expressions. Analytical models for spatial distributions of liquid water content (LWC) and adiabatic fraction (AF) are also designed and analyzed. The results demonstrate the existence of a cloud core with high values of LWC and AF up to the height of 1800 m. The horizontally averaged value of the adiabatic fraction, calculated analytically using the Hill’s vortex concept, is evaluated as 0.39, which is the typical AF value in the upper parts of such Cu. The vertical profiles of different important quantities characterizing cloud structure are presented. The analysis performed in this study allows us to conclude that a rising vortex ring plays the dominating role in formation of the thermodynamic and microphysical structure of developing Cu. Significance Statement 1) Dynamic and thermodynamic fields of a developing cumulus cloud simulated by high-resolution LES with spectral bin microphysics are separated into convective and turbulent components by means of the wavelet technique. 2) The analysis of convective component of the cloud revealed the existence of a vortex ring at the developing stage of the cloud and evaluate its parameters. 3) The analysis performed in this study allows us to conclude that a rising vortex ring plays the important role in formation of the thermodynamic and microphysical structure of developing Cu. 4) The study provides a novel insight into the cloud–environment interaction. 5) The approximating equations describing the vortex ring can be usefully applied for developing new schemes of convective parameterization.

Meteorology & Atmospheric Sciences↗

NASA NTRS1993

Planetary circulations in the presence of transient and self-induced heating

The research program focuses on large-scale circulations and their interaction with the global convective pattern. An 11-year record of global cloud imagery and contemporaneous fields of motion and temperature have been used to investigate organized convection and coherent variability of the tropical circulation operating on intraseasonal time scales. This study provides a detailed portrait of tropical variability associated with the so-called Madden-Julian Oscillation (MJO). It reveals the nature, geographical distribution, and seasonality of discrete convective signal, which is a measure of feedback between the circulation and the convective pattern. That discrete spectral behavior has been evaluated in light of natural variability of the ITCZ associated with climatological convection. A composite signature of the MJO, based on cross-covariance statistics of cloud cover, motion, and temperature, has been constructed to characterize the lifecycle of the disturbance in terms of these properties. The composite behavior has also been used to investigate the influence the MJO exerts on the zonal-mean circulation and the involvement of the MJO in transfers of momentum between the atmosphere and the solid Earth. The aforementioned observational studies have led to the production of two animations. One reveals the convective signal in band-pass filtered OLR and compares it to climatological convection. The other is a 3-dimensional visualization of the composite lifecycle of the MJO. With a clear picture of the MJO in hand, feedback between the circulation and the convective pattern can be diagnosed meaningfully in numerical simulations. This process is being explored in calculations with the linearized primitive equations on the sphere in the presence of realistic stability and shear. The numerical framework represents climatological convection as a space-time stochastic process and wave-induced convection in terms of the vertically-integrated moisture flux convergence. In these calculations, frictional convergence near the equator emerges as a key to feedback between the circulation and the convective pattern. At low latitudes, nearly geostrophic balance in the boundary layer gives way to frictional balance. This shifts the wave-induced convection into phase with the temperature anomaly and allows the attending heating to feed back positively onto the circulation. The calculations successfully reproduce the salient features of the MJO. They are being used to understand the growth and decay phases of the composite lifecycle and the conditions that favor amplification of the MJO.

Salby, Murry L.↗

NASA NTRS2005

The Far-Infrared Spectrum of Arp 220

ISO/LWS grating observations of the ultraluminous infrared galaxy Arp 220 shows absorption in molecular lines of OH, H(sub 2)O, CH, NH, and NH(sub 3), as well as in the [O I] 63 micron line and emission in the [C II] 158 micron line. We have modeled the continuum and the emission/absorption of all observed features by means of a non-local radiative transfer code. The continuum from 25 to 1300 microns is modeled as a warm (106 K) nuclear region that is optically thick in the far-infrared, attenuated by an extended region (size 2") that is heated mainly through absorption of nuclear infrared radiation. The molecular absorption in the nuclear region is characterized by high excitation due to the high infrared radiation density. The OH column densities are high toward the nucleus (2 - 6 x 10(exp 17) cm(exp -2)) and the extended region (approximately 2 x 10(exp 17) cm(exp -2)). The H(sub 2)O column density is also high toward the nucleus (2 - 10 x 10(exp 17) cm(exp -2)) and lower in the extended region. The column densities in a halo that accounts for the absorption by the lowest lying levels are similar to what are found in the diffuse clouds toward the star forming regions in the Sgr B2 molecular cloud complex near the Galactic Center. Most notable are the high column densities found for NH and NH(sub 3) toward the nucleus, with values of approximately 1.5 x 10(exp 16) cm(exp -2) and approximately 3 x 10(exp 16) cm(exp -2), respectively, whereas the NH(sub 2) column density is lower than approximately 2 x 10(exp 15) cm(exp -2). A combination of PDRs in the extended region and hot cores with enhanced H(sub 2)O photodissociation and a possible shock contribution in the nuclei may explain the relative column densities of OH and H(sub 2)O, whereas the nitrogen chemistry may be strongly affected by cosmic ray ionization. The [C II] 158 micron line is well reproduced by our models and its deficit relative to the CII/FIR ratio in normal and starburst galaxies is suggested to be mainly a consequence of the dominant non-PDR component of far-infrared radiation, although our models alone cannot rule out extinction effects in the nuclei.

Gonzalez-Alfonso, Eduardo↗

DOE OSTI2024

High performance power handling in the absence of an H-mode edge in negative triangularity DIII-D plasmas

Abstract Experiments performed during strongly-shaped high-power diverted negative triangularity (NT) experiments in DIII-D achieved detached divertor conditions and a transient-free edge, showcasing the potential for application of NT to a core-edge integrated reactor-like scenario and providing the first characterization of the parametric dependencies for detachment onset. Detached divertor conditions will be required in future devices to mitigate divertor heat fluxes. Access to dissipative divertor conditions was investigated via an increase in upstream density. Detachment onset at the outer strike point was achieved with H-mode level confinement H 98 − y 2 ∼ 1 and reactor-relevant normalized pressures β N ∼ 2 . Confinement degradation was observed with deeper detachment, associated with the loss of an electron temperature pedestal. Differences in geometry, radial transport, impact of cross field drifts are discussed to explain differences in access to detachment in NT discharges. Higher normalized densities, with respect to equivalent discharges in positive triangularity, were necessary to achieve detachment, partially explained by the shorter parallel connection length to the targets. The effect of cross-field particle drifts ( E × B , B × ∇ B ) on access to detachment was demonstrated by the lower upstream density needed to access detachment with ion B × ∇ B drift directed outside of the active divertor (Greenwald fraction f Gw ∼ 0.9–1.0 vs f Gw ∼ 1.3). The upstream density at detachment onset was observed to increase linearly with plasma current with ion B × ∇ B drift into the divertor, consistent with the observed narrowing of the scrape-off layer heat flux width λ q . Edge fluid simulations capture separatrix densities needed to achieve detachment in NT plasma and their dependence on drift direction. The ability to reproduce detachment dynamics in NT plasma increases the confidence in future design studies for NT divertors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗

DOE OSTI2025

Material Needs and Measurement Challenges for Advanced Semiconductor Packaging: Understanding the Soft Side of Science

This Perspective builds upon insights from the National Institute of Standards and Technology (NIST)-organized workshop, “Materials and Metrology Needs for Advanced Semiconductor Packaging Strategies,” held at the 35th annual Electronics Packaging Symposium in Binghamton, NY, on September 5, 2024. It outlines critical challenges and opportunities related to polymer-based “soft” materials in advanced semiconductor packaging, with emphasis on polymer science, measurement science (metrology), and the strategic development of Research-Grade Test Materials (RGTMs). These efforts, led by the NIST CHIPS team, aim to advance the fundamental understanding of structure-property-processing relationships, promote standardized guidelines and innovative methods for material characterization, and accelerate the development, qualification, and adoption of next-generation packaging materials. The Perspective also distills key insights from the panel discussion with industry experts, emphasizing the need for close collaboration among materials scientists, process engineers, and metrology experts to enable a holistic strategy, further highlighting the importance of cross-sector partnerships among industry, academia, and government to address pressing challenges in packaging materials and processes.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1992

Development of high T(sub c) (greater than 110 K) Bi, Tl, and Y-based materials as superconducting circuit elements

Experimental work has continued on the development and characterization of bulk and hot pressed powders and tapecast materials in the Bi-Sr-Ca-Cu-O and Tl-Ba-Ca-Cu-O systems. A process for producing warp-free, sintered, superconducting tapes of Bi composition Bi2Sr2Ca2Cu3O(x) with a mixed oxide process was established. This procedure required a triple calcination at 830 C for 24 hours and sintering at 845 C from 20 to 200 hours. Hot pressing the triple calcined powder at 845 C for 6 hours at 5000 psi yielded a dense material which on further heat treatment at 845 C for 24 hours exhibited a Tc of 108.2 K. A further improvement in the processing of the bismuth materials was achieved via a chemical coprecipitation process wherein the starting nitrate materials were coprecipitated with oxalic acid, thus yielding a more chemically homogeneous, more reactive powder. With the coprecipitated powders, only one calcine at 830 C for 12 hours and a final sinter at 845 C for 30 hours was sufficient to produce a bulk superconducting material with a Tc of 108.4 K. SAFIRE-type grounding links were successfully fabricated from sintered, tapecast, coprecipitated BSCCO 2223 powders. Compositional and processing investigations were continued on the Tl-based superconductors. Manganese and lithium additions and sintering temperature and time were examined to determine their influence on superconducting properties. It was found that lithium substitutions for copper enhance the transition temperatures while manganese additions produced deleterious effects on the superconducting properties. A suitable procedure for producing reproducible bulk and tapecast material of Tl composition Tl2Ba2Ca2Cu3O(x) was developed and used in fabricating uniform superconducting tapes. The highest transition temperature for Tl-based tapes was measured at 110.2 K. Thallium superconducting SAFIRE-type grounding links were fabricated from the tapes.

Haertling, Gene H.↗

DOE OSTI2024

Qualification and Quantification of Porosity at the Top of the Fuel Pins in Metallic Fuels Using Image Processing

Approximately 130,000 metal fuel pins were irradiated in the Experimental Breeder Reactor II (EBR-II) during its 30 years of operation to develop and characterize existing and prospective fuels. For many of the metal fuel irradiation experiments, neutron radiography imaging was performed to characterize fuel behavior, such as fuel axial expansion. While several fuel expansion results obtained from neutron radiography imaging have been published, the analysis of neutron radiography for the purpose of describing statistical properties of porous matter formed on top of the fuel pins, also referred to as fluff in previous publications, is significantly less represented in the literature with just a single paper so far. This study aims to validate and augment results reported in previous publications using automated image processing. The paper describes the statistical properties of the porous matter in terms of nine parameters derived from radiography images and correlates those parameters with such fuel properties as composition, expansion, temperature, and burnup. The reported results are based on 1097 fuel pins of eight different fuel compositions. For three major fuel types, U-10Zr, U-8Pu-10Zr, and U-19Pu-10Zr, a clear negative correlation is found between the Pu content and five parameters describing the amount of porous matter generated. The parameters describing granularity properties, however, showed either negative correlation or nonlinear dependency from fuel composition. The parameters describing the amount showed a positive correlation with fuel axial expansion, while granularity parameters showed a negative correlation with axial expansion. The dependency on cladding temperature was found to be weak. A positive correlation is demonstrated for volume parameters and fuel burnup. In general, reported results confirm and validate findings published in previous studies using a much larger number of pins and automated processing techniques, which easily lend themselves to reproducibility, thus avoiding subjective bias.

36 MATERIALS SCIENCE↗

NASA NTRS1993

Controls on the CO2 seasonal cycle

Surface pressure measurement performed by the Viking landers show substantial variations in pressure on seasonal timescales that are characterized by two local minima and two local maxima. These variations have widely been attributed to the seasonal condensation and sublimation of CO2 in the two polar regions. It has been somewhat of a surprise that the amplitude of the minimum and maximum that is dominated by the CO2 cycle in the north was much weaker than the corresponding amplitude of the south-dominated extrema. Another surprise was that the seasonal pressure cycle during years 2 and 3 of the Viking mission was so similar to that for year 1, despite the occurrence of two global dust storms during year 1 and none during years 2 and 3. An energy balance model that incorporates dynamical factors from general circulation model (GCM) runs in which the atmospheric dust opacity and seasonal date were systematically varied was used to model the observed seasonal pressure variations. The energy balance takes account of the following processes in determining the rates of CO2 condensation and sublimation at each longitudinal and latitudinal grid point: solar radiation, infrared radiation from the atmosphere and surface, subsurface heat conduction, and atmospheric heat advection. Condensation rates are calculated both at the surface and in the atmosphere. In addition, the energy balance model also incorporates information from the GCM runs on seasonal redistribution of surface pressure across the globe. Estimates of surface temperature of the seasonal CO2 caps were used to define the infrared radiative losses from the seasonal polar caps. The seasonal pressure variations measured at the Viking lander sites were closely reproduced.

Pollack, J. B.↗

NASA NTRS1999

Lidar Altimeter Measurements of Canopy Structure: Methods and Validation for Closed Canopy, Broadleaf Forests

Lidar altimeter observations of vegetated landscapes provide a time-resolved measure of laser pulse backscatter energy from canopy surfaces and the underlying ground. Airborne lidar altimeter data was acquired using the Scanning Lidar Imager of Canopies by Echo Recovery (SLICER) for a successional sequence of four, closed-canopy, deciduous forest stands in eastern Maryland. The four stands were selected so as to include a range of canopy structures of importance to forest ecosystem function, including variation in the height and roughness of the outer-most canopy surface and the vertical organization of canopy stories and gaps. The character of the SLICER backscatter signal is described and a method is developed that accounts for occlusion of the laser energy by canopy surfaces, transforming the backscatter signal to a canopy height profile (CHP) that quantitatively represents the relative vertical distribution of canopy surface area. The transformation applies an increased weighting to the backscatter amplitude as a function of closure through the canopy and assumes a horizontally random distribution of the canopy components. SLICER CHPs, averaged over areas of overlap where lidar ground tracks intersect, are shown to be highly reproducible. CHP transects across the four stands reveal spatial variations in vegetation, at the scale of the individual 10 m diameter laser footprints, within and between stands. Averaged SLICER CHPs are compared to analogous height profile results derived from ground-based sightings to plant intercepts measured on plots within the four stands. Tbe plots were located on the segments of the lidar ground tracks from which averaged SLICER CHPs were derived, and the ground observations were acquired within two weeks of the SLICER data acquisition to minimize temporal change. The differences in canopy structure between the four stands is similarly described by the SLICER and ground-based CHP results, however a Chi-square test of similarity documents differences that are statistically significant. The differences are discussed in terms of measurement properties that define the smoothness of the resulting CHPs and Lidar Altimeter Measurements of Canopy Structure - Harding et al. canopy properties that may vertically bias the CHP representations of canopy structure. The statistical differences are most likely due to the more noisy character of the ground-based CHPs, especially high in the canopy where ground-based sightings are rare resulting in an underestimate of canopy surface area and height, and to departures from the assumption of horizontal randomness which bias the CHPs toward the observer (upward for SLICER and downward for ground-based CHPs). The results demonstrate that the SLICER observations reliably provide a measure of canopy structure that reveals ecologically interesting structural variations such as those characterizing a successional sequence of closed-canopy, broadleaf forest stands.

Harding, D. J.↗

DOE OSTI2023

Plasma flows during the ablation stage of an over-massed pulsed-power-driven exploding planar wire array

We characterize the plasma flows generated during the ablation stage of an over-massed exploding planar wire array, fielded on the COBRA pulsed-power facility (1 MA peak current, 250 ns rise time). The planar wire array is designed to provide a driving magnetic field (80–100 T) and current per wire distribution (about 60 kA), similar to that in a 10 MA cylindrical exploding wire array fielded on the Z machine. Over-massing the arrays enables continuous plasma ablation over the duration of the experiment without implosion. The requirement to over-mass on the Z machine necessitates wires with diameters of 75–100μm, which are thicker than wires usually fielded on wire array experiments. To test ablation with thicker wires, we perform a parametric study by varying the initial wire diameter between 33 and 100 μm. The largest wire diameter (100 μm) array exhibits early closure of the cathode-wire gap, while the gap remains open over the duration of the experiment for wire diameters between 33 and 75 μm. Laser plasma interferometry and time-gated extreme-ultraviolet (XUV) imaging are used to probe the plasma flows ablating from the wires. The plasma flows from the wires converge to generate a pinch, which appears as a fast-moving (V≈100kms−1) column of increased plasma density (n¯e≈2×1018cm−3) and strong XUV emission. Finally, we compare the results with three-dimensional resistive-magnetohydrodynamic (MHD) simulations performed using the code GORGON, the results of which reproduce the dynamics of the experiment reasonably well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Automation of Nanoparticle Synthesis Processes in a Plasma Environment Using LabVIEW

This work presents an automated control system for the synthesis of nanomaterials by plasma-enhanced chemical vapor deposition (PECVD), implemented using the LabVIEW software environment. The main objective of the study is to develop an integrated hardware-software platform that enables sequential control of the key stages of the PECVD process, including vacuum chamber preparation, pressure monitoring, working gas supply, plasma ignition, power matching, cyclic nanomaterial growth, and optical monitoring of nanoparticles in the plasma environment. The use of LabVIEW made it possible to integrate actuator control, experimental parameter acquisition, and realtime process visualization within a single automated system. The automated cycle begins with evacuation of the reaction chamber to a predefined base pressure. Transition to the next stage is permitted only after the specified pressure threshold has been reached, ensuring reproducible initial conditions for each experiment. The program then controls the supply of the working gas through mass flow controllers (MFCs). In this work, two gas-flow control modes were considered: analog control using a 0-5 V voltage signal and digital communication via RS-232 interface. It was shown that the analog approach requires accurate scaling of the control voltage, since applying 5 V corresponds to full-scale opening of the controller and results in the maximum gas flow. In contrast, the RS232 interface enables the gas flow rate to be specified directly in sccm, improving the accuracy, flexibility, and convenience of gas-environment control. After pressure stabilization, LabVIEW initiates RF plasma ignition and executes the RF matching algorithm aimed at minimizing reflected power and improving the stability of the plasma process. A separate software module implements the cyclic nanomaterial growth mode, in which the plasma-on time, plasma duration, and total number of synthesis cycles are predefined. This approach makes it possible to control material accumulation on the substrate and to correlate the process parameters with the morphological characteristics of the resulting nanostructures. The final module of the system is designed for optical monitoring of the nanoparticle cloud density in dusty plasma. For this purpose, the change in the intensity of laser radiation passing through the plasma region is recorded using a photodetector and a Keithley 2401 measuring unit connected to LabVIEW via RS-232 interface. The difference between the initial and modified optical signal intensity is used as a diagnostic parameter characterizing the formation and temporal evolution of nanoparticles. The developed system demonstrates that LabVIEW can be effectively applied not only for the automation of individual instruments, but also for the implementation of a complete digital control cycle for PECVD-based nanomaterial synthesis.

PECVD↗

DOE OSTI2026

Findings of the African Combustion Aerosol Collaborative Intercomparison Analysis (ACACIA) Pilot Project to Understand the Optical Properties of Biomass Burning Smoke

Africa is a critical source of biomass burning (BB) aerosols, and its importance is increasing. The African Combustion Aerosol Collaborative Intercomparison Analysis (ACACIA) Pilot Project set to optically characterize BB aerosol generated from sub-Saharan African fuels. We used a photoacoustic spectrometer as a reference instrument to determine the multiple-scattering correction factor C λ for an AE33 aethalometer at three wavelengths, which produced weighted mean values of C 370 =3.69, C 470 =5.65, and C 520 =6.39. C λ increased with wavelength and C 370 was statistically independent of the others, suggesting a single C λ is insufficient, especially in BB scenarios. While a dependence of C λ on burning state was not found, C λ was shown to strongly relate to particle single scattering albedo (SSA, ω). When Cλ was plotted against SSA, values slowly rose at low SSA values, followed by a sharp rise around an SSA of ∼ 0.9; indicating a larger correction needed for less absorbing aerosol. A number of functions operating on either SSA or C λ were explored and the best function was -C λ /(1-C λ )=Aω+B. This is an important parametrization of C λ specifically geared towards BB aerosol from African fuels under different aging states, and is of particular importance for future field work in that continent. An Ångström matrix plot shows that African BB aerosol can have values more akin to dust, which demonstrates that these fuels are distinct in their wavelength dependence from more typical BB aerosol. Lastly, we examined the mass extinction and absorbance cross sections for BB aerosol generated for the same fuels with two different tube furnace setups. Not only is this combustion method flexible, it was found to be reproducible between labs.

47 OTHER INSTRUMENTATION↗

NASA NTRS2013

Validation of Proposed Metrics for Two-Body Abrasion Scratch Test Analysis Standards

Abrasion of mechanical components and fabrics by soil on Earth is typically minimized by the effects of atmosphere and water. Potentially abrasive particles lose sharp and pointed geometrical features through erosion. In environments where such erosion does not exist, such as the vacuum of the Moon, particles retain sharp geometries associated with fracturing of their parent particles by micrometeorite impacts. The relationship between hardness of the abrasive and that of the material being abraded is well understood, such that the abrasive ability of a material can be estimated as a function of the ratio of the hardness of the two interacting materials. Knowing the abrasive nature of an environment (abrasive)/construction material is crucial to designing durable equipment for use in such surroundings. The objective of this work was to evaluate a set of standardized metrics proposed for characterizing a surface that has been scratched from a two-body abrasion test. This is achieved by defining a new abrasion region termed Zone of Interaction (ZOI). The ZOI describes the full surface profile of all peaks and valleys, rather than just measuring a scratch width. The ZOI has been found to be at least twice the size of a standard width measurement; in some cases, considerably greater, indicating that at least half of the disturbed surface area would be neglected without this insight. The ZOI is used to calculate a more robust data set of volume measurements that can be used to computationally reconstruct a resultant profile for de tailed analysis. Documenting additional changes to various surface roughness par ameters also allows key material attributes of importance to ultimate design applications to be quantified, such as depth of penetration and final abraded surface roughness. Further - more, by investigating the use of custom scratch tips for specific needs, the usefulness of having an abrasion metric that can measure the displaced volume in this standardized manner, and not just by scratch width alone, is reinforced. This benefit is made apparent when a tip creates an intricate contour having multiple peaks and valleys within a single scratch. The current innovation consists of a software- driven method of quantitatively evaluating a scratch profile. The profile consists of measuring the topographical features of a scratch along the length of the scratch instead of the width at one location. The digitized profile data is then fed into software code, which evaluates enough metrics of the scratch to reproduce the scratch from the evaluated metrics. There are three key differences between the current art and this innovation. First, scratch width does not quantify how far from the center of the scratch damage occurs (ZOI). Second, scratch width does not discern between material displacement and material removal from the scratch. Finally, several scratches may have the same width but different zones of interactions, different displacements, and different material removals. The current innovation allows quantitative assessment of all three.

Street, Kenneth W., Jr.↗

NASA NTRS2020

Analysis of Sample Acquisition Dynamics Using Discrete Element Method

The analysis presented in this paper is conducted in the framework of the Ocean Worlds Autonomy Testbed for Exploration Research and Simulation (OceanWATERS) project, currently under development at NASA Ames Research Center. OceanWATERS aims at designing a simulation environment which allows for testing autonomy of scientific lander missions to the icy moons of our solar system. Mainly focused on reproducing the end effector interaction with the inherent terrain, this paper introduces a novel discrete element method (DEM)-based approach to determine forces and torques acting on the lander’s scoop during the sample acquisition process. An accurate force feedback from the terrain on the scoop is required by fault-detection and autonomous decision-making algorithms to identify when the requested torque on the robotic arm’s joints exceeds the maximum available torque. Knowledge of the terrain force feedback significantly helps evaluating the arm’s links structural properties and properly selecting actuators for the joints. Models available in literature constitute a partial representation of the dynamics of the interaction. As an example, Balovnev derived an analytical expression of the vertical and horizontal force acting on a bucket while collecting a sample as a function of its geometry and velocity, soil parameters and reached depth. Although the model represents an adequate approximation of the two force components, it ignores the direction orthogonal to the scoop motion and neglects the torque. This work relies on DEM analysis to compensate for analytical models’ deficiencies and inaccuracies, i. e. provide force and torque 3D vectors, defined in the moving reference (body) frame attached to the scoop, at each instant of the sample collection process. Results from the first presented analysis relate to the specific OceanWATERS sampling strategy, which consists of collecting the sample through five consecutive passes with increasing depth, each pass following the same circularlinear- circular trajectory. Data is collected given a specific scoop design interacting with two types of bulk materials, which may characterize the surface of icy planetary bodies: snow and ice. Although specifically concerned with the OceanWATERS design, this first analysis provides the expected force trends for similar sampling strategies and allows to deduce phenomenological information about the general scooping process. In order to further instruct the community on the use of DEM tools as a solution to the sampling collection problem, two more analyses have been carried out, mainly focused on reducing the DEM computation time, which increases with a decrease in particle size. After running a set of identical simulations, where the only changing parameter is the size of the spherical particle, it is observed that the resulting force trajectories, starting from a given particle size, converge to the true trend. It is deducible that a further decrease in size yields negligible improvements in the accuracy, while it sensibly increases computation time. A final analysis aims at discussing limitations of approximating bulk material particles having a complex shape, e. g. ice fragments, with spheres, by comparing force trends resulting in the two cases for the same simulation scenario.

Catanoso, Damiana↗