Search NASA⌕ Search

SEARCH · Search NASA

Results for “Local sensitivity analysis”

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 379 records · Page 21

Probing the possibility of a C-12/C-13 galactic abundance gradient

High S/N (equal to or greater than 500) observations of interstellar CH+ with the 3.0 m telescope were performed at Lick Observatory and with the 4.0 m telescope at CTIO, of the reddened, early-type stars HD 183143, HD 24432, and HD 157038 in an effort to probe the existence of a C-12/C-13 abundance gradient in our Galaxy. Previous very high quality optical observations of interstellar CH+ toward five stars within 1 kpc of the Sun have yielded a precise weighted mean C-12/C-13 isotope ratio of 43 plus or minus 4 (1 sigma) (Hawkins, Jura, and Meyer 1985; Hawkins and Jura 1987). The isotope ratios derived toward four lines of sight in the local ISM are uniform within 12 percent. The similarity among these carbon isotope ratios determined in diffuse clouds possessing different physical conditions precludes the possibility that the CH+ molecule is suffering from isotope selective effects in these regions. The precise C-12/C-13 derivable from high quality observations of (12)CH+ and (13)CH+ provide the unique opportunity to probe the homogeneity of the ISM in a large scale and the history of nucleosynthesis in our Galaxy. Since CH+ seems to be the most sensitive probe of C-12/C-13 in the diffuse ISM, observations toward more distant stars located up to 2.5 kpc from the Sun are the best way to study the possibility of a Galactic (C-12/C-13) abundance gradient. The researcher obtained 4232 angstrom data toward all three of the stars mentioned above, 3957 angstrom data toward HD 183143 and HD 157038, and 3745 angstrom data toward HD 157038. Because of the poorer quality of the HD 24432 spectrum, and its weaker CH+ absorption lines, the satellite (13)CH+ line was not detected, and thus only a lower limit on the C-12/C-13 ratio toward this star was obtained. The results obtained from careful reduction and analysis of the data toward these stars in the Northern and Southern Hemispheres are presented.

Hawkins, Isabel↗

Simulated meteorological impacts of offshore wind turbines and sensitivity to the amount of added turbulence kinetic energy

Offshore wind energy projects are currently in development off the east coast of the United States and may influence the local meteorology of the region. Wind power production and other commercial uses in this area are related to atmospheric conditions, and so it is important to understand how future wind plants may change the local meteorology. In the absence of measurements of potential wind plant impacts on meteorology, simulations offer the next-best possible insight into wake effects on boundary layer height, temperature, fluxes, and wind speeds. However, simulation tools that capture these effects offer multiple options for representing the amount of turbine-added turbulence that may impact assessments of micrometeorological effects. To explore this sensitivity, we compare 1 year of simulations from the Weather Research and Forecasting (WRF) model with and without wind plants incorporated, focusing on the lease area south of Massachusetts and Rhode Island. The simulations with wind plants are repeated to include both the maximum and minimum amounts of added turbulence to provide bounds on the potential impacts. We assess changes in wind speeds, 2 m temperature, surface heat flux, turbulence kinetic energy (TKE), and boundary layer height during different stability classifications and ambient wind speeds over the entire year and compare results for the degree of added turbulence in the wind plant simulations. Because the wake behavior may be a function of boundary layer stability, in this paper, we also present a machine learning algorithm to quantify the area and distance of the wake generated by the wind plant. This analysis enables us to identify the relationship between wake extent and boundary layer height. We find that hub-height wind speed is reduced within and downwind of the wind plant, with the strongest impacts occurring during stable conditions and faster wind speeds in region 3 of the turbine power curve, although impacts lessen as wind speeds increase past 15 m s−1. In contrast, wind speeds near the surface decrease when no turbine-added turbulence is included but can increase for stably stratified conditions when 100 % of possible TKE is included in the simulations. TKE increases at hub height in the simulations with added TKE for all stability classes, suggesting that atmospheric stability does not immediately modify the TKE generated by turbines. Negligible changes in hub-height TKE manifest in the simulations without the added TKE. At the surface, TKE increases in the simulations with maximum added turbulence only for unstable conditions. In the no-added-turbulence simulations, surface TKE decreases slightly in neutral and unstable simulations. Differences in 2 m temperatures and surface heat fluxes are small but vary considerably with atmospheric stability and the amount of added TKE. Boundary layer heights increase within the wind plant when turbine-added turbulence is included and decrease slightly downwind during stable conditions. In contrast, with no added turbulence, the boundary layer height is in general reduced in stable conditions with wind speeds less than 15 m s −1 and slightly increased in neutral conditions. Finally, shallower upwind boundary layer heights tend to correlate with larger wake areas and distances, though other factors likely also play a role in determining the extent of the wind plant wake. These simulation-based results provide a bound for micrometeorological impacts of wind plant wakes: simulations that couple the atmosphere to the ocean may reduce these impacts, and we await observational verification.

17 WIND ENERGY↗

Adjoint Techniques for Topology Optimization of Structures Under Damage Conditions

The objective of this cooperative agreement was to seek computationally efficient ways to optimize aerospace structures subject to damage tolerance criteria. Optimization was to involve sizing as well as topology optimization. The work was done in collaboration with Steve Scotti, Chauncey Wu and Joanne Walsh at the NASA Langley Research Center. Computation of constraint sensitivity is normally the most time-consuming step of an optimization procedure. The cooperative work first focused on this issue and implemented the adjoint method of sensitivity computation (Haftka and Gurdal, 1992) in an optimization code (runstream) written in Engineering Analysis Language (EAL). The method was implemented both for bar and plate elements including buckling sensitivity for the latter. Lumping of constraints was investigated as a means to reduce the computational cost. Adjoint sensitivity computation was developed and implemented for lumped stress and buckling constraints. Cost of the direct method and the adjoint method was compared for various structures with and without lumping. The results were reported in two papers (Akgun et al., 1998a and 1999). It is desirable to optimize topology of an aerospace structure subject to a large number of damage scenarios so that a damage tolerant structure is obtained. Including damage scenarios in the design procedure is critical in order to avoid large mass penalties at later stages (Haftka et al., 1983). A common method for topology optimization is that of compliance minimization (Bendsoe, 1995) which has not been used for damage tolerant design. In the present work, topology optimization is treated as a conventional problem aiming to minimize the weight subject to stress constraints. Multiple damage configurations (scenarios) are considered. Each configuration has its own structural stiffness matrix and, normally, requires factoring of the matrix and solution of the system of equations. Damage that is expected to be tolerated is local and represents a small change in the stiffness matrix compared to the baseline (undamaged) structure. The exact solution to a slightly modified set of equations can be obtained from the baseline solution economically without actually solving the modified system.. Shennan-Morrison-Woodbury (SMW) formulas are matrix update formulas that allow this (Akgun et al., 1998b). SMW formulas were therefore used here to compute adjoint displacements for sensitivity computation and structural displacements in damaged configurations.

Akgun, Mehmet A.↗

Conceptual Designs for Irradiation Creep Testing of SiC in HFIR

Understanding irradiation creep of nuclear fuel cladding is important to properly size the initial fuel-cladding gap and understand when pellet-cladding contact is expected to occur due to a combination of fuel swelling and cladding creep-down. Irradiation creep also plays a role in relaxing stresses that develop in-pile. Silicon carbide fiber–reinforced silicon carbide matrix (SiC/SiC) composites are the leading long-term accident-tolerant fuel cladding concept for light-water reactors (LWRs). Although some limited data are available regarding irradiation creep of the individual constituents (fibers, matrix), data regarding irradiation creep of SiC/SiC composites are currently insufficient. Additional data regarding irradiation creep compliance and the rupture lifetime (combination of creep and slow crack growth) are needed to understand material limitations. This work describes the design and development of two irradiation vehicles that are being pursued for testing SiC/SiC concepts in the High Flux Isotope Reactor (HFIR). The first is a passive experiment, referred to as the PRECISE experiment, that leverages the constant coolant pressure of HFIR to compress a metallic bellows and provide a well-characterized load to drive creep in a SiC/SiC dog bone specimen. The total creep strain would be quantified post-irradiation by measuring dimensional changes of the specimen length as well as local dimensional changes within the gauge region. Non-stressed specimens would also be irradiated under the same conditions to provide an indication of dimensional changes due to radiation-induced swelling in the absence of creep. A second, more complex experiment, referred to as the INSITE experiment, is being designed in parallel that would use pneumatics to pressurize a metal bellows and linear variable differential transformers (LVDTs) to measure the specimen displacement in situ during irradiation. Such an experiment would provide significantly more data regarding the evolution of the creep compliance as a function of dose and applied stress within a single experiment but would require significantly more development time and cost to execute. The primary concern with the INSITE experiment is the accuracy, reliability, and expected lifetime of the LVDTs during irradiation at elevated temperatures. Efforts are being made to adjust the experiment design and operating procedure to limit LVDT temperatures and mitigate or otherwise compensate for uncertainties due to factors such as temperature fluctuations, creep in the surrounding structural materials, and drift of the LVDTs. This work describes the experiment designs, thermal and structural analysis that were performed to ensure that the desired temperature and stress conditions can be achieved, some initial sensitivity analyses to predict the evolution of the radiation-induced specimen displacements, and potential sources of uncertainty in the measurements. Out-of-pile testing is being performed in parallel to confirm that the test trains achieve the expected stress states in the specimens and do not result in prohibitive stress concentrators (e.g., in the grip regions) that might risk pre-mature failure. The PRECISE experiments are proceeding toward fabrication and assembly with HFIR insertion planned during fiscal year 2026. The INSITE experiment is progressing toward out-of-pile demonstrations, which will provide more conclusive evidence regarding the feasibility of executing these tests in HFIR or whether alternative displacement monitoring techniques may need to be considered.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Design and Analysis of Buckling-Critical Large-Scale Sandwich Composite Cylindrical Test Articles

It has long been established in the literature that the buckling response of thin-shell structures can be very sensitive to the presence of small geometric and loading imperfections. The Shell Buckling Knockdown Factor Project (SBKF) was established by the NASA Engineering and Safety Center (NESC) to develop analysis-based shell buckling design recommendations for stiffened-metallic and composite launch-vehicle shell structures. Large-scale buckling tests were used to validate the modeling and analysis methods applied in developing these analysis-based recommendations. Herein, the test article design methodology for 8-ft-diameter, honeycomb-core sandwich composite cylinder validation tests is discussed and cylinder designs are presented. In this methodology, first, the sandwich composite design space was defined using several nondimensional parameters, and the desired test article design space was determined by examining the designs of launch-vehicle cylinder structures. Essentially all test article designs within certain design parameters were generated and then downselected based on simple closed-form failure calculations and the nondimensional design-space parameters. Four of these designs that spanned a significant portion of the design space of interest and had global buckling as the first predicted failure mode were selected and subjected to higher-fidelity finite element analyses (FEAs): shell-element-based analyses, axisymmetric-element-based analyses, and global-local analyses. The analysis flow discussed in this report supported the design objective. As the analysis flow progressed, designs were downselected so the fidelity of the analysis methods, and consequently their computational cost and accuracy, was increased. The selection of the FEA types created an analysis framework where particular methods complemented each other and reduced the uncertainty of the predicted test article responses. The analysis results are illustrated using several designs when the computationally expeditious closed-form analysis stage is discussed. Once this stage is complete, the higher-fidelity FEA types are illustrated using one selected detailed test article design. Both perfect and imperfect test article geometries were considered.

Buckling↗

Modeling Microgravity Induced Fluid Redistribution Autoregulatory and Hydrostatic Enhancements

Space flight induces a marked cephalad (headward) redistribution of blood and interstitial fluid potentially resulting in a loss of venous tone and reduction in heart muscle efficiency upon introduction into the microgravity environment. Using various types of computational models, we are investigating how this fluid redistribution may induce intracranial pressure changes, relevant to reported reductions in astronaut visual acuity, part of the Visual Impairment and Intracranial Pressure (VIIP) syndrome. Methods: We utilize a lumped parameter cardiovascular system (CVS) model, augmented by compartments comprising the cerebral spinal fluid (CSF) space, as the primary tool to describe how microgravity, and the associated lack of hydrostatic gradient, impacts fluid redistribution. Models of ocular fluid pressures and biomechanics then accept the output of the above model as boundary condition input to allow more detailed, local analysis (see IWS Abstract by Ethier et al.). Recently, we enhanced the capabilities our previously reported CVS model through the implementation of robust autoregulatory mechanisms and a more fundamental approach to the implementation of hydrostatic mechanisms. Modifying the approach of Blanco et al., we implemented auto-regulation in a quasi-static manner, as an averaged effect across the span of one heartbeat. This approach reduced the higher frequency perturbations from the regulatory mechanism and was intended to allow longer simulation times (days) than models that implement within-beat regulatory mechanisms (minutes). A more fundamental approach to hydrostatics was implemented by a quasi-1D approach, in which compartment descriptions include compartment length, orientation and relative position, allowed for modeling of body orientation, relative body positioning and, in the future, alternative gravity environments. At this time the inclusion of hydrostatic mechanisms supplies additional capabilities to train and validate the CVS model with terrestrial data. Results and Conclusions: With the implementation of auto-regulation and hydrostatic modeling capabilities, the model performs as expected in the maintaining the CA (Central Artery) compartment pressure when simulating orientations ranging from supine to standing. The model appears to generally overpredict heart rate and thus cardiac output, possibly indicating sensitivity to the nominal heart rate, which is used as an initial set point of the regulation mechanisms. Despite this sensitivity, the model performs consistently for many hours of simulation time, indicating the success of our quasi-static implementation approach.

biomechanics↗

LiDAR-Based Map Relative Localization Performance Analysis for Landing on Europa

This paper presents preliminary simulations andanalyses done to assess the feasibility of performing Map RelativeLocalization (MRL) with the Europa Lander LiDAR beingdeveloped for the Europa Lander Pre-Phase A concept. MapRelative Localization is the process of determining the horizontalposition of a lander with respect to an onboard, a-priori map,by comparing the map to sensor observations of the terrain duringdeorbit, descent, and landing (DDL). Although kilometerscaleposition knowledge is commonly available during DDL,landing in hazard-rich environments requires position errors of100 m or less. Prior knowledge in the case of Europa Landerwill be visual and topographic maps collected by the upcomingEuropa Clipper mission. The Mars 2020 Lander Vision System(LVS) uses images from a camera to localize with respect tovisual maps. This technology, as well as a 3D imaging LiDAR indevelopment for hazard detection, is currently baselined for theEuropa Lander Pre-Phase A concept. This paper investigatesthe potential use of the hazard detection LiDAR to performMRL with respect to a 3D digital elevation model (DEM)provided by the Europa Clipper mission, as an alternative orbackup solution to passive optical MRL. Compared to passiveoptical MRL, one advantage of LiDAR-based localization isthat it is insensitive to lighting conditions, potentially relaxingrequirements on synchronizing map acquisition and landingtime of day. To analyze LiDAR based MRL performance,six representative terrains are synthetically up-sampled fromGalileo-derived maps of Europa to a resolution of 0.5 m/pxand covering an area of 4 km by 4 km. These maps are usedas ground-truth to generate simulated noisy a-priori onboardtopographic maps expected from Europa Clipper as well assimulated LiDAR DEMs generated at an altitude of 5 km duringEuropa Lander DDL. The simulated LiDAR DEM is matchedagainst the simulated map via 2D normalized cross-correlation,exploiting the accurately known spacecraft attitude to avoidthe need for more computationally intensive algorithms such asIterative Closest Point (ICP). Two sources of measurement errorare identified for analysis: 1) additive Gaussian noise in therange measurements from the Europa Lander LiDAR and theEuropa Clipper derived maps and 2) errors in the LiDAR DEMinduced by errors in the Europa Lander state estimate which isused to de-warp the LiDAR scan data into a DEM format. Weassess the effect of each of these types of errors independently onmatching performance as well as the overall performance whenall types of error are introduced. Additionally, we present theresult of a sensitivity study to terrain frequency content.

Trawny, Nikolas↗

Covalent Drug Binding in Live Cells Monitored by Mid-Infrared Quantum Cascade Laser Spectroscopy: Photoactive Yellow Protein as a Model System

The detection of drug-target interactions in live cells enables analysis of therapeutic compounds in a native cellular environment. Recent advances in spectroscopy and molecular biology have facilitated the development of genetically encoded vibrational probes like nitriles that can sensitively report on molecular interactions. Nitriles are powerful tools for measuring electrostatic environments within condensed media like proteins, but such measurements in live cells have been hindered by low signal-to-noise ratios. In this study, we design a spectrometer based on a double-beam quantum cascade laser (QCL)-based transmission infrared (IR) source with balanced detection that can significantly enhance sensitivity to nitrile vibrational probes embedded in proteins within cells compared to a conventional FTIR spectrometer. Here, using this approach, we detect small-molecule binding in Escherichia coli, with particular focus on the interaction between para-Coumaric acid (pCA) and nitrile-incorporated photoactive yellow protein (PYP). This system effectively serves as a model for investigating covalent drug binding in a cellular environment. Notably, we observe large spectral shifts of up to 15 cm –1 for nitriles embedded in PYP between the unbound and drug-bound states directly within bacteria, in agreement with observations for purified proteins. Such large spectral shifts are ascribed to the changes in the hydrogen-bonding environment around the local environment of nitriles, accurately modeled through high-level molecular dynamics simulations using the AMOEBA force field. Our findings underscore the QCL spectrometer’s ability to enhance sensitivity for monitoring drug–protein interactions, offering new opportunities for advanced methodologies in drug development and biochemical research.

chromophores↗

Comparison of The Dimensional Stability of Kel-F 81 and Neoflon CTFE M400H Polychlorotrifluoroethylenes Used in Valve Seat Application

The dimensional stability of polychlorotrifluoroethylene (PCTFE) valve seats used in oxygen regulator applications was determined by thermomechanical analysis (TMA). Two traceable grades of PCTFE were tested; Kel-F 81 and Neoflon CTFE M400H. For these particular resins, the effect of percent crystallinity, zero strength time (ZST) molecular weight, resin grade, process history (compression-molded versus extruded) on the dimensional stability and annealing behavior was determined. In addition to the traceable Kel-F 81 and Neoflon CTFE M400H grades, actual PCI'PH valve seats of differing geometry and design were tested by TMA. The PCTFE valve seats were of unspecified resin grade, although certain inferences about the grade could be drawn based on knowledge of the valve seat fabrication date. Results consistently revealed dimensional instability of varying magnitude at temperatures ranging from 40 to 70 degrees Celsius. Furthermore, some of the pre- 1 995 seats appeared to be more dimensionally stable than those fabricated after 1995. The TMA results are discussed in the context of several proposed ignition mechanisms; namely, particle impact, presence of contaminant oils and fibers, and localized heating by flow friction and/or resonance. The effect of metal constraint on the dimensional stability of PCTFE is also discussed. Finally, the effect of percent crystallinity, ZST molecular weight, resin grade, process history (compression-molded versus extruded) on the AIT, delta Hc and impact sensitivity of various types of Neoflon CTFE M400H was determined using Kel-F 81 as a control. Results show that the AIT, delta Hc and impact sensitivity were essentially independent of Neoflon CTFE process history and structure.

Waller, Jess M.↗

The Effects of Normal Metal Stripes on TES Performance

Exploring the effects of size and geometry of normal metal features on the transition shapes and performance of transition-edge sensor microcalorimeters. The spectral resolution of transition-edge sensor (TES) microcalorimeters is very sensitive to the specific dependencies of the resistance R in the superconducting transition on the current I, magnetic field B, and temperature T. In particular, it has been shown that transitions that are very steep in (R,T) space lead to a significant noise term, in excess of conventional expectations. This so-called unexplained noise is known to be reduced by the addition of normal metal stripes across the TES perpendicular to the direction of current flow. These normal metal stripes have been shown to drastically alter the oscillatory patterns seen in measurements of the critical current as a function of magnetic field. However, there are many remaining questions about the exact impact of the stripes on current distributions within the TES, the Fraunhofer pattern and, therefore, the shape of the R(I, B, T) surface. Through measurements of the resistance under DC bias of TES devices of various sizes, with different stripe patterns and dimensions, we will discuss how these stripes can affect the R(I, B, T) surface. In addition, using measurements and analysis of the noise spectra of various devices we will present how these changes to the stripe pattern may affect the performance of the TES. In particular, we will discuss strategies to reduce the presence of localized discontinuities in the derivative of R, associated with increased noise, while maintaining the globally low levels of unexplained noise currently achieved with conventional metal stripe patterns. Implementing these strategies is a path towards producing large arrays with highly uniform transitions and high spectral resolution. These large uniform arrays will be required for future x-ray astronomy applications, such as the X-IFU on ATHENA.

TES↗

Determining the Covering Factor of Compton-Thick Active Galactic Nuclei with NuSTAR

The covering factor of Compton-thick (CT) obscuring material associated with the torus in active galactic nuclei (AGNs) is at present best understood through the fraction of sources exhibiting CT absorption along the line of sight (N(sub H) greater than 1.5 x 10(exp 24) cm(exp −2)) in the X-ray band, which reveals the average covering factor. Determining this CT fraction is difficult, however, due to the extreme obscuration. With its spectral coverage at hard X-rays (greater than 10 keV), Nuclear Spectroscopic Telescope Array (NuSTAR) is sensitive to the AGNs covering factor since Compton scattering of X-rays off optically thick material dominates at these energies. We present a spectral analysis of 10 AGNs observed with NuSTAR where the obscuring medium is optically thick to Compton scattering, so-called CT AGNs. We use the torus models of Brightman and Nandra that predict the X-ray spectrum from reprocessing in a torus and include the torus opening angle as a free parameter and aim to determine the covering factor of the CT gas in these sources individually. Across the sample we find mild to heavy CT columns, with N(sub H) measured from 10(exp 24) to 10(exp 26) cm(exp −2), and a wide range of covering factors, where individual measurements range from 0.2 to 0.9. We find that the covering factor, f(sub c), is a strongly decreasing function of the intrinsic 2-10 keV luminosity, L(sub X), where f(sub c) = (−0.41 +/- 0.13)log(sub 10)(L(sub X)/erg s(exp −1))+18.31 +/- 5.33, across more than two orders of magnitude in L(sub X) (10(exp 41.5) - 10(exp 44) erg s(exp −1)). The covering factors measured here agree well with the obscured fraction as a function of LX as determined by studies of local AGNs with L(sub X) greater than 10(exp 42.5) erg s(exp −1).

NGC 4945↗

Evidence for Tropopause Layer Moistening by Convection During CRYSTAL-FACE

Measurements and analysis of the impact of deep convection on tropopause layer moisture are easily confounded by difficulties making precise observations with sufficient spatial coverage before and after convective events and difficulties distinguishing between changes due to local convection versus large-scale advection. The interactions between cloud microphysics and dynamics in the convective transport of moisture into the tropopause layer also result in a sufficiently complex and poorly characterized system to allow for considerable freedom in theoretical models of stratosphere-troposphere exchange. In this work we perform detailed large-eddy simulations with an explicit cloud microphysics model to study the impact of deep convection on tropopause layer moisture profiles observed over southern Florida during CRYSTALFACE. For four days during the campaign (July 11, 16, 28, and 29) we initialize a 100-km square domain with temperature and moisture profiles measured prior to convection at the PARSL ground site, and initiate convection with a warm bubble that produces an anvil at peak elevations in agreement with lidar and radar observations on that day. Comparing the moisture field after the anvils decay with the initial state, we find that convection predominantly moistens the tropopause layer (as defined by minimum temperature and minimum potential temperature lapse rate), although some drying is also predicted in localized layers. We will also present results of sensitivity tests designed to separate the roles of cloud microphysics and dynamics.

Ackerman, A.↗

Weather Research and Forecasting Model Sensitivity Comparisons for Warm Season Convective Initiation

Mesoscale weather conditions can significantly affect the space launch and landing operations at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS). During the summer months, land-sea interactions that occur across KSC and CCAFS lead to the formation of a sea breeze, which can then spawn deep convection. These convective processes often last 60 minutes or less and pose a significant challenge to the forecasters at the National Weather Service (NWS) Spaceflight Meteorology Group (SMG). The main challenge is that a "GO" forecast for thunderstorms and precipitation is required at the 90 minute deorbit decision for End Of Mission (EOM) and at the 30 minute Return To Launch Site (RTLS) decision at the Shuttle Landing Facility. Convective initiation, timing, and mode also present a forecast challenge for the NWS in Melbourne, FL (MLB). The NWS MLB issues such tactical forecast information as Terminal Aerodrome Forecasts (TAFs), Spot Forecasts for fire weather and hazardous materials incident support, and severe/hazardous weather Watches, Warnings, and Advisories. Lastly, these forecasting challenges can also affect the 45th Weather Squadron (45 WS), which provides comprehensive weather forecasts for shuttle launch, as well as ground operations, at KSC and CCAFS. The need for accurate mesoscale model forecasts to aid in their decision making is crucial. Both the SMG and the MLB are currently implementing the Weather Research and Forecasting Environmental Modeling System (WRF EMS) software into their operations. The WRF EMS software allows users to employ both dynamical cores - the Advanced Research WRF (ARW) and the Non-hydrostatic Mesoscale Model (NMM). There are also data assimilation analysis packages available for the initialization of the WRF model- the Local Analysis and Prediction System (LAPS) and the Advanced Regional Prediction System (ARPS) Data Analysis System (ADAS). Having a series of initialization options and WRF cores, as well as many options within each core, provides SMG and NWS MLB with a lot of flexibility. It also creates challenges, such as determining which configuration options are best to address specific forecast concerns. The goal of this project is to assess the different configurations available and to determine which configuration will best predict warm season convective initiation in East-Central Florida. Four different combinations of WRF initializations will be run (ADAS-ARW, ADAS-NMM, LAPS-ARW, and LAPS-NMM) at a 4-km resolution over the Florida peninsula and adjacent coastal waters. Five candidate convective initiation days using three different flow regimes over East-Central Florida will be examined, as well as two null cases (non-convection days). Each model run will be integrated 12 hours with three runs per day, at 0900, 1200, and 1500 UTe. ADAS analyses will be generated every 30 minutes using Level II Weather Surveillance Radar-1988 Doppler (WSR-88D) data from all Florida radars to verify the convection forecast. These analyses will be run on the same domain as the four model configurations. To quantify model performance, model output will be subjectively compared to the ADAS analyses of convection to determine forecast accuracy. In addition, a subjective comparison of the performance of the ARW using a high-resolution local grid with 2-way nesting, I-way nesting, and no nesting will be made for select convective initiation cases. The inner grid will cover the East-Central Florida region at a resolution of 1.33 km. The authors will summarize the relative skill of the various WRF configurations and how each configuration behaves relative to the others, as well as determine the best model configuration for predicting warm season convective initiation over East-Central Florida.

Watson, Leela R.↗

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING↗

Stellar Laboratories II. New Zn Iv and Zn v Oscillator Strengths and Their Validation in the Hot White Dwarfs G191-B2B and RE0503-289

Context. For the spectral analysis of high-resolution and high-signal-to-noise (SN) spectra of hot stars, state-of-the-art non-local thermodynamic equilibrium (NLTE) model atmospheres are mandatory. These are strongly dependent on the reliability of the atomic data that is used for their calculation. In a recent analysis of the ultraviolet (UV) spectrum of the DA-type white dwarf G191B2B,21 Zn iv lines were newly identified. Because of the lack of Zn iv data, transition probabilities of the isoelectronic Ge vi were adapted for a first, coarse determination of the photospheric Zn abundance.Aims. Reliable Zn iv and Zn v oscillator strengths are used to improve the Zn abundance determination and to identify more Zn lines in the spectra of G191B2B and the DO-type white dwarf RE 0503289. Methods. We performed new calculations of Zn iv and Zn v oscillator strengths to consider their radiative and collisional bound-bound transitions in detail in our NLTE stellar-atmosphere models for the analysis of the Zn iv v spectrum exhibited in high-resolution and high-SN UV observations of G191B2B and RE 0503289. Results. In the UV spectrum of G191B2B, we identify 31 Zn iv and 16 Zn v lines. Most of these are identified for the first time in any star. We can reproduce well almost all of them at log Zn 5.52 0.2 (mass fraction, about 1.7 times solar). In particular, the Zn iv Zn v ionization equilibrium, which is a very sensitive Teff indicator, is well reproduced with the previously determined Teff 60 000 2000 K and log g 7.60 0.05. In the spectrum of RE 0503289, we identified 128 Zn v lines for the first time and determined log Zn 3.57 0.2 (155 times solar). Conclusions. Reliable measurements and calculations of atomic data are a pre-requisite for stellar-atmosphere modeling. Observed Zn iv and Zn v line profiles in two white dwarf (G191B2B and RE 0503289) ultraviolet spectra were well reproduced with our newly calculated oscillator strengths. This allowed us to determine the photospheric Zn abundance of these two stars precisely.

strengths↗

Updated Human Mars Ascent Vehicle Concept in Support of NASA’s Strategic Analysis Cycle 2021

The NASA Artemis program has brought significant change to the agency’s human exploration strategy over the last several years. The human exploration of Mars remains an ultimate objective of the overall human exploration strategy. However, how the agency intends to execute initial human exploration of Mars has shifted to better align with current policy. These changes find their way into the design of key architecture elements, such as the Mars Ascent Vehicle (MAV). Of the numerous changes in the conceptual human Mars architecture, two in particular have had significant impact on the design of the human MAV: the desire for minimal surface infrastructure for initial human Mars missions, and technology investment timelines in support of initial human missions to Mars. The first leads to a surface mission architecture targeting a 30-day surface duration with two crewmembers. The second drives out in-situ resource utilization (ISRU) from the initial sortie. As a result, the reference MAV design for initial human Mars missions in the Human Exploration and Operations Mission Directorate’s Strategic Analysis Cycle 2021 (SAC21) is based around supporting the return of 2 crew from the surface of Mars, without surface ISRU. This has led to several design changes since 2019, when details of a reference human MAV were last published. The MAV concept for SAC21 consists of a two-stage vehicle capable of supporting two crew for 84 hours. This duration is derived from a baseline ascent trajectory targeting an apoapsis rendezvous with the habitation element in a nominal 5-sol Mars aggregation orbit. With the goal of minimizing technology investments for the initial surface sortie removing ISRU from the trade space, the MAV must now be landed either fully loaded with propellant, or with some portion of the propellant off-loaded. Landing a partially-fueled MAV will require additional surface assets to perform robotic propellant transfer on the surface of Mars from a pre-emplaced propellant depot. Furthermore, cryogenic fluid management and storage for extended duration with no losses represents additional technology investment that would have to be made to support cryogenic-based propellant concepts that have been the baseline in past architectures. As a result, the SAC21 architecture utilizes a storable propellant-based MAV propulsion system to minimize these potential technology investments incurred by cryogenic-based propulsion systems. Recent efforts also focused on evaluating the sensitivity of key design parameters, such as landing site latitude, elevation, and local atmospheric conditions. All play a large role in determining the predicted propellant requirements of the MAV. Understanding their potential impacts is important because of the rippling effects of changes to the MAV design on other architecture elements, such as entry, descent, landing, and transportation systems. Results indicate up to +2% to -4% wet mass variation due to landing site latitude, with higher latitudes greater than 70 degrees North resulting in greater than 6% wet mass increases. The MAV saw reduced impact due to atmospheric design parameters, which only accounting for -0.6% to +0.4% wet mass variation.

Douglas J. Trent↗

ALMA Observations of Molecular Clouds in Three Group-centered Elliptical Galaxies:NGC 5846, NGC 4636, and NGC 5044

We present new ALMA CO(2–1) observations of two well-studied group-centered elliptical galaxies: NGC 4636 and NGC 5846. In addition, we include a revised analysis of Cycle 0 ALMA observations of the central galaxy in the NGC 5044 group. We find evidence that molecular gas is a common presence in bright group-centered galaxies (BGG). CO line widths are broader than Galactic molecular clouds, and using the reference Milky Way XCO, the total molecular mass ranges from 2.6 x 10(exp 5)M(☉) in NGC 4636 to 6.1 x 10(exp 7) M(☉) in NGC 5044. Complementary observations using the ALMA Compact Array do not exhibit any detection of a CO diffuse component at the sensitivity level achieved by current exposures. The origin of the detected molecular features is still uncertain, but these ALMA observations suggest that they are the end product of the hot gas cooling process and not the result of merger events. Some of the molecular clouds are associated with dust features as revealed by HST dust extinction maps, suggesting that these clouds formed from dust-enhanced cooling. The global nonlinear condensation may be triggered via the chaotic turbulent field or buoyant uplift. The large virial parameter of the molecular structures and correlation with the warm (10(exp 3)–10(exp 5) K)/hot (≥ 10(exp 6)) phase velocity dispersion provide evidence that they are unbound giant molecular associations drifting in the turbulent field, consistent with numerical predictions of the chaotic cold accretion process. Alternatively, the observed large CO line widths may be generated by molecular gas flowing out from cloud surfaces due to heating by the local hot gas atmosphere.

Pasquale Temi↗

Using Whispering-Gallery-Mode Resonators for Refractometry

A method of determining the refractive and absorptive properties of optically transparent materials involves a combination of theoretical and experimental analysis of electromagnetic responses of whispering-gallery-mode (WGM) resonator disks made of those materials. The method was conceived especially for use in studying transparent photorefractive materials, for which purpose this method affords unprecedented levels of sensitivity and accuracy. The method is expected to be particularly useful for measuring temporally varying refractive and absorptive properties of photorefractive materials at infrared wavelengths. Still more particularly, the method is expected to be useful for measuring drifts in these properties that are so slow that, heretofore, the properties were assumed to be constant. The basic idea of the method is to attempt to infer values of the photorefractive properties of a material by seeking to match (1) theoretical predictions of the spectral responses (or selected features thereof) of a WGM of known dimensions made of the material with (2) the actual spectral responses (or selected features thereof). Spectral features that are useful for this purpose include resonance frequencies, free spectral ranges (differences between resonance frequencies of adjacently numbered modes), and resonance quality factors (Q values). The method has been demonstrated in several experiments, one of which was performed on a WGM resonator made from a disk of LiNbO3 doped with 5 percent of MgO. The free spectral range of the resonator was approximately equal to 3.42 GHz at wavelengths in the vicinity of 780 nm, the smallest full width at half maximum of a mode was approximately equal to 50 MHz, and the thickness of the resonator in the area of mode localization was 30 microns. In the experiment, laser power of 9 mW was coupled into the resonator with an efficiency of 75 percent, and the laser was scanned over a frequency band 9 GHz wide at a nominal wavelength of approximately equal to 780 nm. Resonance frequencies were measured as functions of time during several hours of exposure to the laser light. The results of these measurements, plotted in the figure, show a pronounced collective frequency drift of the resonator modes. The size of the drift has been estimated to correspond to a change of 8.5 x 10(exp -5) in the effective ordinary index of refraction of the resonator material.

Matsko, Andrey↗