SEARCH · Engineering Papers

Results for “Distributed model predictive control”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Using Strong Magnetic Fields to Control Solutal Convection

An important component in biotechnology, particularly in the area of protein engineering and rational drug design is the knowledge of the precise three-dimensional molecular structure of proteins. The quality of structural information obtained from X-ray diffraction methods is directly dependent on the degree of perfection of the protein crystals. As a consequence, the growth of high quality macromolecular crystals for diffraction analyses has been the central focus for biochemists, biologists, and bioengineers. Macromolecular crystals are obtained from solutions that contain the crystallizing species in equilibrium with higher aggregates, ions, precipitants, other possible phases of the protein, foreign particles, the walls of the container, and a likely host of other impurities. By changing transport modes in general, i.e., reduction of convection and sedimentation, as is achieved in microgravity , we have been able to dramatically affect the movement and distribution of macromolecules in the fluid, and thus their transport, formation of crystal nuclei, and adsorption to the crystal surface. While a limited number of high quality crystals from space flights have been obtained, as the recent National Research Council (NRC) review of the NASA microgravity crystallization program pointed out, the scientific approach and research in crystallization of proteins has been mainly empirical yielding inconclusive results. We postulate that we can reduce convection in ground-based experiments and we can understand the different aspects of convection control through the use of strong magnetic fields and field gradients. We postulate that limited convection in a magnetic field will provide the environment for the growth of high quality crystals. The approach exploits the variation of fluid magnetic susceptibility with concentration for this purpose and the convective damping is realized by appropriately positioning the crystal growth cell so that the magnetic susceptibility force counteracts terrestrial gravity. The general objective is to test the hypothesis of convective control using a strong magnetic field and magnetic field gradient and to understand the nature of the various forces that come into play. Specifically we aim to delineate causative factors and to quantify them through experiments, analysis and numerical modeling. The paper will report on the experimental results using paramagnetic salts and solutions in magnetic fields and compare them to analytical predictions.

Ramachandran, N.↗

NASA NTRS2018

Towards an Enhanced Droplet Activation Scheme for Multi-Moment Bulk Microphysics Schemes

Initial droplet spectra produced upon activation impact the ensuing chain of microphysical processes andtherefore play a crucial role in cloud evolution. This work re-examines dependencies of newly formed clouddroplet size distribution (CDSD) characteristics on environmental and aerosol properties via parcel model simulationsthat serve as the basis for a multi-moment bulk microphysics droplet activation scheme suitable for acloud-resolving model (CRM). It is found that applying a fixed size threshold to define activated droplets versusemploying physical considerations can lead to erroneous activation and overly broad CDSDs for high aerosolconcentration and weak updraft conditions. Aerosol distributions characterized by larger median sizes and/orincreased solubility can result in greater activated droplet numbers, whereas impacts of these parameters onCDSD spectral width depend on both aerosol number concentration and updraft velocity. An expansion of theactivation scheme to include CDSD spectral width is proposed to aid efforts to extend high-order momentprediction to cloud droplet categories in CRMs as well as better represent variability in the activation process onthe cloud scale.simulations to investigate the regime dependence of the relative dispersion(d)1 of newly activated CDSDs, where d is the ratio of dropletradius standard deviation (σ) to the mean radius (r ). C16 demonstratedthat increasing Na resulted in increasing (decreasing) d values via reducedcondensational narrowing (spectral broadening) rates within theAL (UL) regime, with d values peaking in the TR regime. Their findingssuggest a similar regime dependence for d as R09 noted for Nc and helpexplain reportedly conflicting relationships between Na and CDSDspectral characteristics (cf. Hudson and Noble, 2014; Liu et al., 2014),although the applicability of these results within bulk microphysicalschemes was not addressed.Simulating aerosol-cloud interactions with CRMs employing bulkmicrophysics requires that the scheme minimally predict two CDSDparameters, namely mass and number concentrations, and represent thedroplet activation process. Various activation schemes aim to determineNc from aerosol and environmental properties and include analyticalexpressions (e.g., Abdul-Razzak et al., 1998; Morrison et al., 2005) aswell as lookup tables (LUTs) based on detailed parcel model calculations(e.g., Saleeby and Cotton, 2004, hereafter SC04; Segal and Khain,2006; Thompson and Eidhammer, 2014). Expressions to diagnose CDSDspectral width from Nc (Grabowski, 1998; Liu et al., 2006; Morrison andGrabowski, 2007) or cloud water content (Geoffroy et al., 2010) havealso been developed, although more robust methods to obtain CDSDspectral width upon activation are presently lacking. This latter point isrelevant for triple-moment (3 M) bulk microphysics that aim to predictdistribution spectral width alongside number and mass concentrations(e.g., Loftus et al., 2014; Milbrandt and Yau, 2005).The current work extends the findings of C16 to the current LUTbasedaerosol activation scheme used in the Regional AtmosphericModeling System (RAMS) (Cotton et al., 2003; SC04; Saleeby and vanden Heever, 2013, hereafter SvdH13) and additionally examinesaerosol size and solubility impacts on newly activated CDSD properties.Because early cloud development processes such as condensationalgrowth, evaporation, and droplet self-collection depend on and impactCDSD spectral width (Hudson and Yum, 1997; Seifert and Beheng 2001;Lu and Seinfeld, 2006; Igel and van den Heever, 2017), an expansion ofthe activation LUTs to include CDSD spectral width is proposed as apreliminary step for extending 3M prediction to CDSDs in CRMs forimproved simulations of aerosol-cloud interactions.2. MethodologyThe current RAMS two-moment microphysics module determinesthe fractional number of aerosol particles that activate to cloud dropletsfrom five-dimensional LUTs based on model predicted air temperature(T), w, Na, and the geometric median radius (rg) and soluble fraction (ε)of the aerosol size distribution (SvdH13). These LUTs are created offlineusing a one-dimensional Lagrangian adiabatic parcel model (Feingoldand Heymsfield, 1992; Heymsfield and Sabin, 1989; SC04) to simulateexplicit droplet activation and initial CDSD growth for a range of ambientatmospheric conditions [T, w] and binned lognormal aerosol sizedistributions given by= ⎡⎣ ⎢− ⎤⎦ ⎥N r Nr π σr rσ( )2 lnexp[ln( / )]2(ln )aggg22 (1)where r is the dry aerosol particle bin radius and σg is the geometricstandard deviation of the distribution. As the parcel model simulationsfocus on the activation process, other processes such as coalescence,sedimentation, and mixing are not considered. Details of the parcelmodel can be found in SC04 and SvdH13, and only a brief description isprovided here. At the onset of parcel model calculations, the initiallydry aerosol particles in all bins first deliquesce and reach theirequilibrium diameters in a sub-saturated environment based on theKöhler equation for solution droplets. The parcel is then lifted at a fixedupward velocity w, and particle growth by vapor diffusion, along withconcurrent changes in the ambient environment, are iteratively computedusing the Variable-coefficient Ordinary Differential Equation(VODE) solver (Brown et al., 1989). The time resolution of these calculationsis determined within the VODE solver, and the frequency atwhich the solver is called is controlled by a longer model time stepbased on fixed upward parcel displacement increments (Δz) at thespecified w (Δt=Δz/w). Model calculations proceed until the parcelreaches a height 50m beyond the level of maximum saturation ratio(Smax) or total parcel displacement exceeds 2 km. Upon model termination,Smax and the fractional number of aerosols (factv) resulting innewly formed cloud droplets, defined as particles having diameters of atleast 2 μm, are cataloged in the LUTs according to the specified T, w, Na,rg, and ε parameter values.A critical point regarding the creation of these LUTs is the use of afixed minimum diameter (Dmin) to define cloud droplets in the parcelmodel, which can produce erroneous CDSD characteristics, particularlywithin the UL regime. For aerosol distributions with large rg valuesunder low SS conditions, for example, deliquesced aerosols within thelarge tail of the distribution can exceed 2 μm in diameter yet remainunactivated as ‘haze’ particles (Levin and Cotton, 2009; McFigganset al., 2006). For this study, aerosol particles activate to cloud dropletsbased on the critical diameter Dcrit as a function of parcel supersaturationratio (Sr) as in R09:D = σ MS RTρ83 ln( ) critsol wr w (2)where σsol is the surface tension of a solution droplet, Mw and ρw are themolar mass and density of liquid water, respectively, and R is theuniversal gas constant. Additionally, at relatively large w values withinthe AL regime, Nc stabilizes shortly after reaching supersaturation.However, parcel ascent and condensational growth continue beyondthe level of Smax, potentially causing additional narrowing of the CDSD.In the current work, model calculations terminate upon reaching Smaxas changes in Nc are negligible with continued ascent (Peng et al., 2007;R09).Parcel model simulations are performed to examine the sensitivitiesof CDSD characteristics to w, Na, rg, and ε, with the ranges for theseparameters listed in Table 1. Aerosols are assumed to be a mix of solubleand insoluble material of equal density, specified by ε, where fullysoluble aerosols correspond to ammonium sulfate with hygroscopicityparameter κ=0.61 (Petters and Kreidenweis, 2007). FollowingSvdH13, aerosol geometric standard deviation is fixed at σg=1.8, andaerosol distributions (Eq. 1) are partitioned into 100 logarithmicallyspacedbins spanning a size range specific to each rg value. For all simulations,Δz=1 m, and initial values of relative humidity, air temperatureand pressure are set to RH=0.99, T=10 °C and p=900 hPa,respectively.

Loftus, Adrian M.↗

DOE OSTI2021

Topographical Controls on Hillslope-Scale Hydrology Drive Shrub Distributions on the Seward Peninsula, Alaska

Observations indicate shrubs are expanding across the Arctic tundra, mainly on hillslopes and primarily in response to climate warming. However, the impact topography exerts on hydrology, nutrient dynamics, and plant growth can make untangling the mechanisms behind shrub expansion difficult. We examined the role topography plays in determining shrub expansion by applying a coupled transect version of a mechanistic ecosystem model (ecosys) in a tundra hillslope site in the Seward Peninsula, Alaska. Modeled biomass of the dominant plant functional types agreed well with field measurements (R 2 = 0.89) and accurately represented shrub expansion over the past 30 years inferred from satellite observations. In the well-drained crest position, canopy water potential and plant nitrogen (N) uptake was modeled to be low from plant and microbial water stress. Intermediate soil water content in the mid-slope position enhanced mineralization and plant N uptake, increasing shrub biomass. The deciduous shrub growth in the mid-slope position was further enhanced by symbiotic N 2 fixation primed by increased root carbon allocation. The gentle slope in the poorly drained lower-slope position resulted in saturated soil conditions that reduced soil O 2 concentrations, leading to lower root O 2 uptake and lower nutrient uptake and plant biomass. Additionally, a simulation that removed topographical interconnectivity between grid cells resulted in (1) a 28% underestimate of mean shrub biomass and (2) over or underestimated shrub productivity at the various hillslope positions. Our results indicate that land models need to account for hillslope-scale coupled surface and subsurface hydrology to accurately predict current plant distributions and future trajectories in Arctic ecosystems.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2000

The Use of Ferrofluids to Model Materials Processing (MSFC Center Director's Discretionary Fund)

Many crystals grown in space have structural flaws believed to result from convective motions during the growth phase. The character of these instabilities is not well understood but is associated with thermal and solutal density variations near the solidification interface in the presence of residual gravity and g-jitter. To study these instabilities in a separate, controlled space experiment, a concentration gradient would first have to be artificially established in a timely manner as an initial condition. This is generally difficult to accomplish in a microgravity environment because the momentum of the fluid injected into a test cell tends to swirl around and mix in the absence of a restoring force. The use of magnetic fields to control the motion and position of liquids has received recent, growing interest. The possibility of using the force exerted by a non-uniform magnetic field on a ferrofluid to not only achieve fluid manipulation but also to actively control fluid motion makes it an attractive candidate for space applications. This paper describes a technique for quickly establishing a linear or exponential fluid concentration gradient using a magnetic field in place of gravity to stabilize the deployment. Also discussed is a photometric technique for measuring the concentration profile using light attenuation. Although any range of concentrations can be realized, photometric constraints impose some limitations on measurements. Results of the ground-based experiments indicate that the species distribution is within 3 percent of the predicted value.

Leslie, F.↗

NASA NTRS1995

Parasitic current collection by PASP Plus solar arrays

Solar cells at potentials positive with respect to a surrounding plasma collect electrons. Current is collected by the exposed high voltage surfaces: the interconnects and the sides of the solar cells. This current is a drain on the array power that can be significant for high-power arrays. In addition, this current influences the current balance that determines the floating potential of the spacecraft. One of the objectives of the Air Force (PL/GPS) PASP Plus (Photovoltaic Array Space Power Plus Diagnostics) experiment is an improved understanding fo parasitic current collection. We have done computer modeling of parasitic current collection and have examined current collection flight data from the first year of operations. Prior to the flight we did computer modeling to improve our understanding of the physical processes that control parasitic current collection. At high potentials, the current rapidly rises due to a phenomenon called snapover. Under snapover conditions, the equilibrium potential distribution across the dielectric surface is such that part of the area is at potentials greater than the first crossover of the secondary yield curve. Therefore, each incident electron generates more than one secondary electron. The net effect is that the high potential area and the collecting area increase. We did two-dimensional calculations for the various geometries to be flown. The calculations span the space of anticipated plasma conditions, applied potential, and material parameters. We used the calculations and early flight data to develop an analytic formula for the dependence of the current on the primary problem variables. The analytic formula was incorporated into the EPSAT computer code. EPSAT allows us to easily extend the results to other conditions. PASP Plus is the principal experiment integrated onto the Advanced Photovoltaic and Electronics Experiments (APEX) satellite bus. The experiment is testing twelve different solar array designs. Parasitic current collection is being measured for eight of the designs under various operational and environment conditions. We examined the current collected as a function of the various parameters for the six non-concentrator designs. The results are similar to those obtained in previous experiments and predicted by the calculations. We are using the flight data to validate the analytic formula developed. The formula can be used to quantify the parasitic current collected. Anticipating the parasitic current value allows the spacecraft designer to include this interaction when developing the design.

Davis, Victoria Ann↗

NASA NTRS1995

Combining Speed Information Across Space

We used speed discrimination tasks to measure the ability of observers to combine speed information from multiple stimuli distributed across space. We compared speed discrimination thresholds in a classical discrimination paradigm to those in an uncertainty/search paradigm. Thresholds were measured using a temporal two-interval forced-choice design. In the discrimination paradigm, the n gratings in each interval all moved at the same speed and observers were asked to choose the interval with the faster gratings. Discrimination thresholds for this paradigm decreased as the number of gratings increased. This decrease was not due to increasing the effective stimulus area as a control experiment that increased the area of a single grating did not show a similar improvement in thresholds. Adding independent speed noise to each of the n gratings caused thresholds to decrease at a rate similar to the original no-noise case, consistent with observers combining an independent sample of speed from each grating in both the added- and no-noise cases. In the search paradigm, observers were asked to choose the interval in which one of the n gratings moved faster. Thresholds in this case increased with the number of gratings, behavior traditionally attributed to an input bottleneck. However, results from the discrimination paradigm showed that the increase was not due to observers' inability to process these gratings. We have also shown that the opposite trends of the data in the two paradigms can be predicted by a decision theory model that combines independent samples of speed information across space. This demonstrates that models typically used in classical detection and discrimination paradigms are also applicable to search paradigms. As our model does not distinguish between samples in space and time, it predicts that discrimination performance should be the same regardless of whether the gratings are presented in two spatial intervals or two temporal intervals. Our last experiment largely confirmed this prediction.

Verghese, Preeti↗

NASA NTRS1995

Combining Speed Information Across Space

Verghese, Preeti↗

DOE OSTI2016

A spatially distributed model for assessment of the effects of changing land use and climate on urban stream quality: Development of a Spatially Distributed Urban Water Quality Model

While the effects of land use change in urban areas have been widely examined, the combined effects of climate and land use change on the quality of urban and urbanizing streams have received much less attention. We describe a modeling framework that is applicable to the evaluation of potential changes in urban water quality and associated hydrologic changes in response to ongoing climate and landscape alteration. The grid-based spatially distributed model, DHSVM-WQ, is an outgrowth of the Distributed Hydrology-Soil-Vegetation Model (DHSVM) that incorporates modules for assessing hydrology and water quality in urbanized watersheds at a high spatial and temporal resolution. DHSVM-WQ simulates surface runoff quality and in-stream processes that control the transport of nonpoint-source (NPS) pollutants into urban streams. We configure DHSVM-WQ for three partially urbanized catchments in the Puget Sound region to evaluate the water quality responses to current conditions and projected changes in climate and/or land use over the next century. Here we focus on total suspended solids (TSS) and total phosphorus (TP) from nonpoint sources (runoff), as well as stream temperature. The projection of future land use is characterized by a combination of densification in existing urban or partially urban areas, and expansion of the urban footprint. The climate change scenarios consist of individual and concurrent changes in temperature and precipitation. Future precipitation is projected to increase in winter and decrease in summer, while future temperature is projected to increase throughout the year. Our results show that urbanization has a much greater effect than climate change on both the magnitude and seasonal variability of streamflow, TSS and TP loads largely due to substantially increased streamflow, and particularly winter flow peaks. Water temperature is more sensitive to climate warming scenarios than to urbanization and precipitation changes. Future urbanization and climate change together are predicted to significantly increase annual mean streamflow (up to 55%), water temperature (up to 1.9 ºC), TSS load (up to 182%), and TP load (up to 74%).

Sun, Ning↗

DOE OSTI2015

Novel Sensor for the In Situ Measurement of Uranium Fluxes

The goal of this project was to develop a sensor that incorporates the field-tested concepts of the passive flux meter to provide direct in situ measures of flux for uranium and groundwater in porous media. Measurable contaminant fluxes [J] are essentially the product of concentration [C] and groundwater flux or specific discharge [q ]. The sensor measures [J] and [q] by changes in contaminant and tracer amounts respectively on a sorbent. By using measurement rather than inference from static parameters, the sensor can directly advance conceptual and computational models for field scale simulations. The sensor was deployed in conjunction with DOE in obtaining field-scale quantification of subsurface processes affecting uranium transport (e.g., advection) and transformation (e.g., uranium attenuation) at the Rifle IFRC Site in Rifle, Colorado. Project results have expanded our current understanding of how field-scale spatial variations in fluxes of uranium, groundwater and salient electron donor/acceptors are coupled to spatial variations in measured microbial biomass/community composition, effective field-scale uranium mass balances, attenuation, and stability. The coupling between uranium, various nutrients and micro flora can be used to estimate field-scale rates of uranium attenuation and field-scale transitions in microbial communities. This research focuses on uranium (VI), but the sensor principles and design are applicable to field-scale fate and transport of other radionuclides. Laboratory studies focused on sorbent selection and calibration, along with sensor development and validation under controlled conditions. Field studies were conducted at the Rifle IFRC Site in Rifle, Colorado. These studies were closely coordinated with existing SBR (formerly ERSP) projects to complement data collection. Small field tests were conducted during the first two years that focused on evaluating field-scale deployment procedures and validating sensor performance under controlled field conditions. In the third and fourth year a suite of larger field studies were conducted. For these studies, the uranium flux sensor was used with uranium speciation measurements and molecular-biological tools to characterize microbial community and active biomass at synonymous wells distributed in a large grid. These field efforts quantified spatial changes in uranium flux and field-scale rates of uranium attenuation (ambient and stimulated), uranium stability, and quantitatively assessed how fluxes and effective reaction rates were coupled to spatial variations in microbial community and active biomass. Analyses of data from these field experiments were used to generate estimates of Monod kinetic parameters that are ‘effective’ in nature and optimal for modeling uranium fate and transport at the field-scale. This project provided the opportunity to develop the first sensor that provides direct measures of both uranium (VI) and groundwater flux. A multidisciplinary team was assembled to include two geochemists, a microbiologist, and two quantitative contaminant hydrologists. Now that the project is complete, the sensor can be deployed at DOE sites to evaluate field-scale uranium attenuation, source behavior, the efficacy of remediation, and off-site risk. Because the sensor requires no power, it can be deployed at remote sites for periods of days to months. The fundamental science derived from this project can be used to advance the development of predictive models for various transport and attenuation processes in aquifers. Proper development of these models is critical for long-term stewardship of contaminated sites in the context of predicting uranium source behavior, remediation performance, and off-site risk.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2012

NEAMS Software Licensing, Release, and Distribution: Implications for FY2013 Work Package Planning

The vision of the NEAMS program is to bring truly predictive modeling and simulation (M&S) capabilities to the nuclear engineering community in order to enable a new approach to the analysis of nuclear systems. NEAMS anticipates issuing in FY 2018 a full release of its computational 'Fermi Toolkit' aimed at advanced reactor and fuel cycles. The NEAMS toolkit involves extensive software development activities, some of which have already been underway for several years, however, the Advanced Modeling and Simulation Office (AMSO), which sponsors the NEAMS program, has not yet issued any official guidance regarding software licensing, release, and distribution policies. This motivated an FY12 task in the Capability Transfer work package to develop and recommend an appropriate set of policies. The current preliminary report is intended to provide awareness of issues with implications for work package planning for FY13. We anticipate a small amount of effort associated with putting into place formal licenses and contributor agreements for NEAMS software which doesn't already have them. We do not anticipate any additional effort or costs associated with software release procedures or schedules beyond those dictated by the quality expectations for the software. The largest potential costs we anticipate would be associated with the setup and maintenance of shared code repositories for development and early access to NEAMS software products. We also anticipate an opportunity, with modest associated costs, to work with the Radiation Safety Information Computational Center (RSICC) to clarify export control assessment policies for software under development.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS2014

Analysis and Prediction of Weather Impacted Ground Stop Operations

When the air traffic demand is expected to exceed the available airport's capacity for a short period of time, Ground Stop (GS) operations are implemented by Federal Aviation Administration (FAA) Traffic Flow Management (TFM). The GS requires departing aircraft meeting specific criteria to remain on the ground to achieve reduced demands at the constrained destination airport until the end of the GS. This paper provides a high-level overview of the statistical distributions as well as causal factors for the GSs at the major airports in the United States. The GS's character, the weather impact on GSs, GS variations with delays, and the interaction between GSs and Ground Delay Programs (GDPs) at Newark Liberty International Airport (EWR) are investigated. The machine learning methods are used to generate classification models that map the historical airport weather forecast, schedule traffic, and other airport conditions to implemented GS/GDP operations and the models are evaluated using the cross-validations. This modeling approach produced promising results as it yielded an 85% overall classification accuracy to distinguish the implemented GS days from the normal days without GS and GDP operations and a 71% accuracy to differentiate the GS and GDP implemented days from the GDP only days.

Analysis↗

DOE OSTI2017

Adaptive RAS/SPS System Settings for Improving Grid Reliability and Asset Utilization through Predictive Simulation and Controls

Remedial Action Schemes (RAS) or Special Protection Schemes (SPS) are used throughout the bulk transmission system as a non-wires method of increasing transmission transfer capability. However, as the population grows within the Western Interconnection, more power needs to be transferred with less ability to build lines. This leads to a greater number of RAS, and the subsequent interaction between these disparate schemes can have unintended consequences. While the large number of RAS within the Western Electricity Coordinating Council (WECC) protect for different things, which can range from the violation of thermal limitations to transient stability issues, they can be described by common methods of operation. In order to understand these methods, a survey was sent out to members of the WECC Modeling & Validation Working Group (MVWG). The response received provides a necessary base in which to understand and build models that can be used within WECC. This report provides information on the survey results regarding RAS/SPS modeling practices for planning and operational studies within the WECC footprint, which is foundational for on-going and future work on RAS. Pacific Northwest National Laboratory (PNNL) and its industry partners are developing innovative methods for adaptively setting RAS parameters based on realistic and near real-time operational conditions to improve power grid reliability and grid asset utilization.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2013

DARK MATTER HALO PROFILES OF MASSIVE CLUSTERS: THEORY VERSUS OBSERVATIONS

Dark-matter-dominated cluster-scale halos act as an important cosmological probe and provide a key testing ground for structure formation theory. Focusing on their mass profiles, we have carried out (gravity-only) simulations of the concordance {Lambda}CDM cosmology, covering a mass range of 2 Multiplication-Sign 10{sup 12} to 2 Multiplication-Sign 10{sup 15} h {sup -1} M{sub Sun} and a redshift range of z = 0-2, while satisfying the associated requirements of resolution and statistical control. When fitting to the Navarro-Frenk-White profile, our concentration-mass (c-M) relation differs in normalization and shape in comparison to previous studies that have limited statistics in the upper end of the mass range. We show that the flattening of the c-M relation with redshift is naturally expressed if c is viewed as a function of the peak height parameter, {nu}. Unlike the c-M relation, the slope of the c-{nu} relation is effectively constant over the redshift range z = 0-2, while the amplitude varies by {approx}30% for massive clusters. This relation is, however, not universal: using a simulation suite covering the allowed wCDM parameter space, we show that the c-{nu} relation varies by about {+-}20% as cosmological parameters are varied. At fixed mass, the c(M) distribution is well fit by a Gaussian with {sigma}{sub c}/(c) {approx_equal} 1/3, independent of the radius at which the concentration is defined, the halo dynamical state, and the underlying cosmology. We compare the {Lambda}CDM predictions with observations of halo concentrations from strong lensing, weak lensing, galaxy kinematics, and X-ray data, finding good agreement for massive clusters (M{sub vir} > 4 Multiplication-Sign 10{sup 14} h {sup -1} M{sub Sun }), but with some disagreements at lower masses. Because of uncertainty in observational systematics and modeling of baryonic physics, the significance of these discrepancies remains unclear.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Geoelectrical Measurement of Multi-Scale Mass Transfer Parameters

Mass transfer affects contaminant transport and is thought to control the efficiency of aquifer remediation at a number of sites within the Department of Energy (DOE) complex. An improved understanding of mass transfer is critical to meeting the enormous scientific and engineering challenges currently facing DOE. Informed design of site remedies and long-term stewardship of radionuclide-contaminated sites will require new cost-effective laboratory and field techniques to measure the parameters controlling mass transfer spatially and across a range of scales. In this project, we sought to capitalize on the geophysical signatures of mass transfer. Previous numerical modeling and pilot-scale field experiments suggested that mass transfer produces a geoelectrical signature—a hysteretic relation between sampled (mobile-domain) fluid conductivity and bulk (mobile + immobile) conductivity—over a range of scales relevant to aquifer remediation. In this work, we investigated the geoelectrical signature of mass transfer during tracer transport in a series of controlled experiments to determine the operation of controlling parameters, and also investigated the use of complex-resistivity (CR) as a means of quantifying mass transfer parameters in situ without tracer experiments. In an add-on component to our grant, we additionally considered nuclear magnetic resonance (NMR) to help parse mobile from immobile porosities. Including the NMR component, our revised study objectives were to: 1. Develop and demonstrate geophysical approaches to measure mass-transfer parameters spatially and over a range of scales, including the combination of electrical resistivity monitoring, tracer tests, complex resistivity, nuclear magnetic resonance, and materials characterization; and 2. Provide mass-transfer estimates for improved understanding of contaminant fate and transport at DOE sites, such as uranium transport at the Hanford 300 Area. To achieve our objectives, we implemented a 3-part research plan involving (1) development of computer codes and techniques to estimate mass-transfer parameters from time-lapse electrical data; (2) bench-scale experiments on synthetic materials and materials from cores from the Hanford 300 Area; and (3) field demonstration experiments at the DOE’s Hanford 300 Area. In a synergistic add-on to our workplan, we analyzed data from field experiments performed at the DOE Naturita Site under a separate DOE SBR grant, on which PI Day-Lewis served as co-PI. Techniques developed for application to Hanford datasets also were applied to data from Naturita. 1. Introduction The Department of Energy (DOE) faces enormous scientific and engineering challenges associated with the remediation of legacy contamination at former nuclear weapons production facilities. Selection, design and optimization of appropriate site remedies (e.g., pump-and-treat, biostimulation, or monitored natural attenuation) requires reliable predictive models of radionuclide fate and transport; however, our current modeling capabilities are limited by an incomplete understanding of multi-scale mass transfer—its rates, scales, and the heterogeneity of controlling parameters. At many DOE sites, long “tailing” behavior, concentration rebound, and slower-than-expected cleanup are observed; these observations are all consistent with multi-scale mass transfer [Haggerty and Gorelick, 1995; Haggerty et al., 2000; 2004], which renders pump-and-treat remediation and biotransformation inefficient and slow [Haggerty and Gorelick, 1994; Harvey et al., 1994; Wilson, 1997]. Despite the importance of mass transfer, there are significant uncertainties associated with controlling parameters, and the prevalence of mass transfer remains a point of debate [e.g., Hill et al., 2006; Molz et al., 2006] for lack of experimental methods to verify and measure it in situ or independently of tracer breakthrough. There is a critical need for new field-experimental techniques to measure mass transfer in-situ and estimate multi-scale and spatially variable mass-transfer parameters. The current lack of such techniques results in large parameter uncertainty, which in turn translates into enormous prediction uncertainty and cost to DOE. In this project, we considered three hydrogeophysical approaches for providing information about mass-transfer parameters: (1) the combination of electrical-resistivity tomography (ERT) and ionic tracer experiments to explore rates of exchange and relative mobile and immobile porosities; (2) complex resistivity (CR) measurements to infer the distribution of diffusive length scales active in a porous medium; and (3) nuclear magnetic resonance (NMR) to estimate mobile and immobile porosity.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2022

Random Forest Regressor-Based Approach for Detecting Fault Location and Duration in Power Systems

Power system failures or outages due to short-circuits or “faults” can result in long service interruptions leading to significant socio-economic consequences. It is critical for electrical utilities to quickly ascertain fault characteristics, including location, type, and duration, to reduce the service time of an outage. Existing fault detection mechanisms (relays and digital fault recorders) are slow to communicate the fault characteristics upstream to the substations and control centers for action to be taken quickly. Fortunately, due to availability of high-resolution phasor measurement units (PMUs), more event-driven solutions can be captured in real time. In this paper, we propose a data-driven approach for determining fault characteristics using samples of fault trajectories. A random forest regressor (RFR)-based model is used to detect real-time fault location and its duration simultaneously. This model is based on combining multiple uncorrelated trees with state-of-the-art boosting and aggregating techniques in order to obtain robust generalizations and greater accuracy without overfitting or underfitting. Four cases were studied to evaluate the performance of RFR: 1. Detecting fault location (case 1), 2. Predicting fault duration (case 2), 3. Handling missing data (case 3), and 4. Identifying fault location and length in a real-time streaming environment (case 4). A comparative analysis was conducted between the RFR algorithm and state-of-the-art models, including deep neural network, Hoeffding tree, neural network, support vector machine, decision tree, naive Bayesian, and K-nearest neighborhood. Experiments revealed that RFR consistently outperformed the other models in detection accuracy, prediction error, and processing time.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

Measurements of TRACER pre-convective conditions and mesoscale circulations using small unmanned aircraft systems (sUAS)

Improved comprehension of the physical processes governing convective cloud formation and lifecycle are of critical importance for understanding and predicting future climate states. The influence of these clouds on the planetary energy budget, including on precipitation, is significant. Things are particularly complex in coastal regimes, where gradients in aerosol particle properties, localized circulations such as sea breezes, and large population centers are found. To date, numerical models struggle to accurately represent these critical clouds and are therefore challenged to provide a realistic view on the planetary energy budget. Through the proposed research, we deployed two uncrewed aircraft systems (UAS) equipped with a variety of instruments alongside sensors deployed by the US Department of Energy Atmospheric Radiation Measurement (ARM) program for the TRACER (Tracking Aerosol Convection Interactions Experiment) field campaign. The two small UAS platforms consisted of a CU RAAVEN fixed-wing airplane and an OU CopterSonde system, with the copter collecting frequent vertical profiles of thermodynamic and kinematic variables such as temperature, pressure, wind and humidity. At the same time, the fixed-wing captured horizontal gradients of these quantities and aerosol size distribution. These systems were deployed south of the Houston metro area, in an area that is impacted by the Gulf of Mexico sea breeze on a daily basis. These observations offer enhanced and complementary perspectives to those provided by the DOE ARM Mobile Facility (AMF) which is was deployed in southeast Houston, and an ancillary site in a more rural location west of the urban Houston area. Quality-controlled versions of the UAS data were collected and posted on the DOE ARM data archive after the conclusion of the campaign where they are accessible by the research community and general public. The UAS perspective offers revolutionary insight into key spatial and temporal effects that have not been evaluated previously.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRSDate not supplied

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand [1] while simultaneously maintaining one of the safest forms of transportation [2], [3]. One of the reasons for this success is the ability of the system and the operators to adapt and accommodate to situations that routinely disrupt optimal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators ability to control. These factors can lead to states where automation is unable to properly handle these issues and therefore air traffic controllers and pilots have to intervene, ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions this can be viewed as an increase in complexity. The reasoning for this is because humans are now required to make tactical decisions in response to external factors, resulting in a departure from the strategic plan where operations would be more efficiently managed. Human operators control airspace complexity under rigid regulations that are constantly changing. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. There has been past work that devised airspace complexity metrics in commercial aviation and related these metrics to controller workload (e.g., [4],[5]). The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic, including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly [6]. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to our proposed effort that identifies such contributing factors or precursor patterns. To define the scope we are proposing to measure complexity from the viewpoint of the Terminal Radar Approach Control Facilities (TRACON) controller’s perspective. In particular we are analyzing arrivals into KSFO. With safety as the top concern for airspace operators, it is important to recognize that as density and heterogeneity grow, the focus of the system will change. Times of the day when the airspace has low density and heterogeneity, the flights will follow more efficient paths where the aircraft move on established routes that are more or less directly to the destination. However, when density and heterogeneity increases, the system will begin changing focus to avoiding conflicts and collisions and route the flights in a more flexible way. Higher flexibility requires more communication and coordination between controllers and pilots which the current automation is unable to handle. This paper proposes a novel approach that monitors airspace complexity at multiple scales, uses a Machine Learning-based tool that predicts when operations will transition to a regime of greater complexity, and identifies actions that can reduce the complexity while still maintaining efficient and safe operations. We demonstrate our proposed approach using data from multiple complementary sources. This includes, but is not limited to: historical aircraft surveillance data from NASA’s Sherlock Data Warehouse [7], METAR weather data, and airport configuration data from Aviation System Performance Metrics (ASPM). The surveillance data flight paths are sampled at a variable sample rate — increasing as the aircraft approaches the airport. This is due to how Sherlock manages flight track stitching between different radar facilities which have different sampling rates. The weather and performance data are logged at defined intervals throughout the day at a courser refresh rate. In addition to the logged data and metrics, we leverage pre-defined Standard Terminal Arrival Routes (STARs) procedures to characterize the path of each flight. Each flight files for one of these routes in the flight plan well before entering the terminal airspace, and approximately follows the route until it leaves the STAR, typically on the final fix of a runway transition. However, most flights do not always fly the full STAR procedure to completion [8], but the majority do adhere to the fixes within the common route of the procedure. Our approach leverages fixes in the common route of each of the STARs to build a reference path to the airport. This allows us to characterize the flight paths in what we are defining as the “maneuvering area” (the airspace between the STAR and before the flight is lined up on the runway’s final approach) to determine how off nominal the flights are to calculate its complexity score. Determining airspace complexity is a concept that does not have a concrete answer. In designing this metric, we consider what increases the workload for the air traffic controllers. Consequently more specialized vectoring maneuvers results in higher workload. Accordingly, we start with a theory: each flight has a direct path it takes from the STAR’s common route to the final approach’s outer marker fix for the flight’s landing runway. It is important to note that the direct path is only used as a reference. If the majority of the flights have a large consistent offset as compared to other routes it does not necessarily mean that those flights have higher complexity. We are merely building a distribution based on this direct path for that particular STAR and runway pair to determine the normal mode of operations for that route. Flights that are in the upper tail of these distributions will result in higher complexity scores and flights that fly in the median will represent the normal mode of operations and therefore will have lower complexity scores. Since flights following each STAR route take different paths to the airport, we have a different distribution for each STAR route and therefore can model these distributions to compute a complexity score from their respective normalized distributions. To evaluate the effectiveness of our proposed airspace complexity metric we will compare against an established approach based on trajectory clustering [9]. This unsupervised learning technique consists of the following steps: (1) identify the general maneuvering areas (waypoints) by performing $\kappa$-means or DBSCAN clustering on locations where aircraft frequently turn based on the surveillance radar track data, (2) map flight trajectories onto sequences of waypoints, and (3) cluster the sequences based on their common subsequences. From a high-level perspective, this baseline model learns nominal operations in the airspace through the sequence of waypoints that are representative of where aircraft change direction and defines deviations from the nominal operations as “complex.” Therefore, more deviations from the nominal operations correspond to higher complexity values. For our validation, we re-implemented this technique and tune model hyper-parameters to correctly detect waypoints for the arrival traffic into the San Francisco bay area. We will compute the complexity measure over a one-year period using our proposed technique as well as the baseline. Our validation will be based on each technique’s ability to detect a set of undesirable outcomes (e.g., go-arounds, holding patterns, average time in the airspace, etc.). Since our current complexity metric is derived from the offset from the direct reference path, it’s important to understand what causes these offsets. In many of the flights with high offset distance, flights performing holding patterns and S turns can be observed. These maneuvering tactics are utilized to add distance between the aircraft and the destination runway to prevent multiple flights from having conflicting arrival times. In order to predict a rise in complexity (or the precursor to complexity), it’s necessary to be able to identify these potential conflicts (which in turn, result in higher offsets). To do this, we define a “representative flight” for each STAR route and runway pair. This flight is approximately the path the flight would take if there was a clear path with no other flights in the airspace — including the time remaining to the airport. We first identify the flights for a given STAR runway pair using the offset to the reference path distributions that fall between the 44-55 percentiles. This yields the flights that conform to the most normal mode of operation. Each of these flights is partitioned based on the percent complete from the entry point into the maneuvering areas from 0\% – 100\% complete. Then for each percent “bin”, we take the median value of the flight’s latitude/longitude coordinates, airspeed, and (non causal) time remaining to the airport to construct a lookup table for each percent complete bin on a given route. As a flight enters the maneuvering area, we can find the estimated arrival time of a flight to the airport by finding the closest point to the representative path’s percent complete bin (relative to the flight’s current position at any snapshot in the airspace) and therefore retrieve the corresponding remaining time left on the “representative path”. We assume that the flight will follow the representative path to completion when deriving these estimates. We can then compare these estimated arrival times against other flights for the same snapshot in time to identify potential conflicts. If more flights are estimated to arrive within a tolerance window than there are runways available, then we have a potential conflict. We can use this derived measure along with other factors expected to add disruption to the operation such as weather and runway configuration changes as an input to machine learning tools to detect precursors that increases in our complexity measure. This novel method will assist in uncovering insights into the contributing factors that lead to increased complexity that may allow for in-time responses to avoid reaching a high complexity state in the airspace.

complexity↗

NASA NTRS2021

Adaptive Aerostructures for Revolutionary Civil Supersonic Transportation

To enable commercially viable civil supersonic transport (SST) aircraft, innovative solutions must be developed to meet noise and efficiency requirements for overland flight. This research effort consists of a multidisciplinary team of academic and industrial experts exploring for the first time the potential of small real-time geometric outer mold line (OML) reconfigurations to minimize sonic boom signatures and aircraft drag in response to changing ambient conditions, thereby enabling noise-compliant overland supersonic flight. The team utilizes recent advances in supersonic computational fluid dynamic (CFD) methods, new noise prediction tools, and new design approaches to consider embedded highly energy-dense shape memory alloy (SMA) actuators for local shape modifications to an SST aircraft leading to optimal low boom signature and low drag in different environments. This university-led program will provide strategic leadership toward technology convergence that advances NASA's Aerospace Research Mission Directorate's (ARMD) research objectives with regard to Thrust 2: “Innovation in Commercial Supersonic Aircraft” by exploring for the first time enabling low-boom operation across a range of flight conditions via structural adaptivity, and will promote education of the next generation of engineers. The overall research strategy is to pursue three critical areas: the design of configurations for reducing boom, SMA material development and modeling, and technology feasibility demonstration in a relevant environment. Initially, the team will identify potential applications where structure or geometry adaptivity provides a benefit in noise or drag across the entire flight envelope. For selected applications/structural locations, required OML geometry changes will be determined based on analysis of sonic boom ground signature and drag reduction using new design tools, trade studies, and atmospheric sensing techniques. Designs will be developed and evaluated against requirements on loading, stroke length, and operational temperature. New alloy formulations will be developed tailored for both autonomous and controlled actuation modes. As the SMA material development matures, integrated system-level factors will be investigated. Optimized designs for small-scale distributed adaptivity applications of maximum benefit will then be matured and tested, moving toward demonstration of the innovative technology approaches at a TRL 4-5 and showing that sonic booms can be reduced by reconfiguration on demand.

Koushik Datta↗